Abstract

Within the framework of the effective mass approximation, the ground-state binding energy of a hydrogenic impurity is investigated in cylindrical wurtzite GaN/ strained quantum ring (QR) by means of a variational approach, considering the influence of the applied hydrostatic pressure along the QR growth direction and the strong built-in electric field (BEF) due to the piezoelectricity and spontaneous polarization. Numerical results show that the donor binding energy for a central impurity increases inchmeal firstly as the QR radial thickness decreases gradually and then begins to drop quickly. In addition, the donor binding energy is an increasing (a decreasing) function of the inner radius (height). It is also found that the donor binding energy increases almost linearly with the increment of the applied hydrostatic pressure. Moreover, we also found that impurity positions have an important influence on the donor binding energy. The physical reasons have been analyzed in detail.

1. Introduction

In recent years, nitride-based wurtzite semiconductor materials AlN, GaN, and InN and their alloys have been widely studied due to their promising applications in optoelectronic devices such as light-emitting devices and laser diodes [13]. The electronic and optical properties of various confined systems such as quantum wells [48], quantum well wires [9, 10], and (double) quantum dots [1115] can be changed by doped impurities. Therefore, many theoretical investigations concerning the impurity states have been done in semiconducting nanostructures by well-known approaches such as variational approach [415], strong confinement approach [16, 17], and perturbation approach [18]. As expected, the hydrostatic pressure applied on a bulk material can not only modify the parameters, such as the band gaps, the potential barriers, the conduction effective masses, the static dielectric constants, and the lattice constants, but also change the dimension of the low-dimensional systems, which is associated with the fractional change in the volume. Moreover, the strong built-in electric field also affects obviously the electric and optical properties of the wurtzite nitride-based quantum heterostructures. According to the above characteristics, the hydrogenic donor impurity states in QWs (QWWs) [410] have been extensively studied under the influence of the electric and magnetic fields. Various studies [11, 12, 1922] related to the effects of hydrostatic pressure and built-in electric field on the donor binding energy in QDs have also been reported.

Since the first GaAs/ quantum ring (QR) of nanoscopic sizes was fabricated [22, 23] and exhibited fascinating behaviors such as interesting electronic and optical properties. The electronic structure and optical properties of arsenic based QR under external influences have been extensively studied. Li and Xia [24, 25] studied the electronic states of InAs/GaAs QR (GaAs/ DQR) in the frame of effective mass envelope function theory by employing the matrix diagonalization method. The papers show that the electron energy levels are sensitively dependent on the radial thickness, the height of the QR (DQR), and the magnetic field strength. Except for the energy levels of the carriers, other effects can be obtained in the QR, such as the Aharonov-Bohm oscillation induced by the magnetic field [26]. On the other hand, a theoretical study of shallow donor states in GaAs-(Ga,Al)As QRs, within the effective mass approximation and using a variational method, which is presented in [27]. The explicit dependencies of the impurity binding energy on the magnetic field strength, the structure parameters of the QR, and the impurity positions are obtained. In addition, considering the effects of hydrostatic pressure, temperature, aluminum concentration, and impurity position, Baghramyan et al. [28] investigated the binding energy of hydrogenic donor impurity in GaAs/ concentric double quantum rings. And Monozon and Schmelcher [29] have studied the problem of an impurity in a QR in the presence of a radially directed strong external electric field by means of the analytical and numerical approaches. On the other hand, Barseghyan and coworkers [30, 31] have investigated the behavior of the binding energy and photoionization cross section of a hydrogenic-like donor impurity in InAs quantum ring with Pöschl-Teller confinement potential along the axial direction by using the variational method. In the above papers, the combined effects of hydrostatic pressure and electric and magnetic fields applied in the growth direction have been taken into account. Moreover, we [32] also studied the hydrogenic impurity binding energy in GaAs/ QR under the external electric field in finite potential barrier by means of a variational approach.

As far as we know, there are no reports on the hydrostatic pressure and built-in electric field (BEF) effects on the donor binding energy in wurtzite GaN/ strained QR to date. Therefore, it is very necessary to investigate the donor impurity states in wurtzite GaN/ strained QR. In our work, the effects of hydrostatic pressure, impurity positions, and sizes of the structure on the donor binding energy in wurtzite GaN/ strained QR are investigated by means of a variational approach. The paper is organized as follows. In Section 2 we describe the theoretical framework. Section 3 is dedicated to the results and discussions, and our conclusions are given in Section 4.

2. Theoretical Framework

2.1. Hamiltonian and Impurity States

In Figure 1, the schematic view of a cylindrical wurtzite GaN/ QR is depicted, with a detailed description of the different dimensions of the QR (inner radius , outer radius , height , and the QR radial thickness ). Additionally, the GaN ring is embedded in a host matrix material, and the -axis is defined to be the QR growth direction. Within the frame of the effective mass approximation, the Hamiltonian for a hydrogenic donor impurity in cylindrical wurtzite GaN/ QR under the influence of hydrostatic pressure can be written as where denotes the position vector of the electron (impurity ion), is the absolute value of the electron charge, is the permittivity of free space, and is the pressure-dependent effective mean relative dielectric constant of GaN and materials. The Hamiltonian is given by [11] where and are the pressure-dependent effective mass of the electron along and perpendicular to the -direction. In (2), is the pressure-dependent BEF in the finitely thick barrier layer. The values of the BEF along the QR growth direction in the ring () and the barrier , which results from the difference in the total electric polarizations in each region, are given by simple formulas [14] as follows: where is the QR height and is the thickness of the barrier layer along the QR growth direction. is the dielectric constant for GaN () material; , and , are the spontaneous and piezoelectric polarizations for GaN and materials. And is the electron confinement potential due to the band offset () and is given by The wave function of an electron in the wurtzite GaN/ QR can be written as [24] where the constants (, , and ) are determined by the continuity of the derivative of the radial wave function at the QR boundary and is the electron -component angular momentum quantum number. The radial wave function of the electron can be obtained using the Bessel function , and the modified Bessel function , . The wave function can be expressed by means of the Airy functions Ai and Bi where . The coefficients and (, and 5) can also be determined by the transfer matrix methods [20].

In order to calculate the donor binding energy, the trial wave function can be chosen as [6] where is the normalization constant. With the adiabatic approximation, the donor binding energy of a hydrogenic impunity is defined as the difference between the ground-state energy of the system without impurity and the ground-state energy of the system with impurity [26]; that is,

2.2. Pressure and Strain Dependence of Physical Parameters

In this model, we take the strains induced by the biaxial lattice mismatch into account. The components of biaxial stress tensors of GaN and AlGaN materials are given under the hydrostatic pressure [33] where the equilibrium lattice constant for the strained layer under the hydrostatic pressure depends on the lattice constants of the component materials, and it is weighted by their relative thicknesses The lattice constant of the material GaN(AlN) dependence of hydrostatic pressure [34] satisfies where the actual lattice constant of the material can be obtained by the linear interpolation method from the corresponding values of GaN and AlN. The strain induced by the biaxial lattice mismatch along the -direction in the heterostructure is [35] The coefficient of the GaN (AlN) material is given by [35] where is the pressure-dependent elastic stiffness constant of material and is given by [36] . In addition, the coefficient and the pressure-dependent elastic stiffness constant of the material can be obtained by the linear interpolation method from the corresponding values of GaN and AlN. The strain-dependent energy gaps of GaN and AlN are given by means of the following expression [37]: where , , , and ( = GaN and AlN) are the deformation potentials. The dependence of the energy gap on hydrostatic pressure is considered by the following equation [38]: The energy gap of the alloy can be calculated from energy gaps of GaN and AlN with the simplified coherent potential approximation (SCPA) [39]: The biaxial strain and hydrostatic pressure dependence of the electron effective masses in the radial direction and the -direction can be calculated by [40] where is a fixed value for a given material and can be derived from the values of and . The electron effective mass of the material can be obtained by the linear interpolation method from the corresponding values of GaN and AlN. In (5), considering the biaxial and uniaxial strain and hydrostatic pressure effects, the dielectric constants and the phonon eigenfrequencies will be changed. The static dielectric constant in material can be derived from the generalized Lyddane-Sachs-Teller relation [35] where the LO- and TO-phonon frequencies influenced by biaxial strain and hydrostatic pressure can be written as [35] Furthermore, the hydrostatic pressure dependence of can be determined by the given mode-Grüneisen parameters [41] where the subscript represents LO- or TO-phonon, ( = or ) denotes plane or -direction, is the zone-center phonon frequency of material , is Grüneisen parameter of phonon mode given in [35], is bulk modulus, and and are the strain coefficients of zone-center phonon modes.

Following [14], considering the influence of hydrostatic pressure, the high frequency dielectric constants in (19) can be written as where ( = GaN and AlN) is the Phillips iconicity parameter of material and is the static bulk modulus under hydrostatic pressure and is given by [41] where , , , and are the elastic constants of material .

The effective mean relative dielectric constant in (1) is defined as [39] The piezoelectric polarization along the -oriented wurtzite GaN/ QR can be calculated as [14] where and are the pressure-dependent piezoelectric constants of material and are given by [14] Here, is the Born effective charge, is the pressure-dependent equilibrium lattice constant, and is the anion-cation bond length along the -direction in units of . Based on the definition of the strain tensor and the bulk modulus, we can derive the following formula for the pressure dependence of the QR inner (outer) radius , the height of the QR, and the thickness of barrier along the QR growth direction: where , , and are the QR inner (outer) radius, the QR height, and the thickness of the barrier layer along the QR growth direction at zero pressure, respectively. and are the compliance constants of GaN and materials, respectively [42].

3. Results and Discussions

Considering the effects of the hydrostatic pressure and strong built-in electric field due to the spontaneous and piezoelectric polarizations, we calculated the donor binding energy of a hydrogenic impurity as a function of structural parameters, such as inner radius , outer radius , height , and impurity position along the growth direction -axis in cylindrical wurtzite GaN/ QR. All material parameters used in our calculations are listed in Tables 14. The material parameters can be calculated from the corresponding parameters of GaN and AlN by using a linear interpolation (Vegard’ slaw). In this paper, the Al concentration is 0.15.

Figure 2 displays the built-in electric field (BEF) FGaN (FAlGaN) as a function of the hydrostatic pressure for different QR heights with barrier thickness ( nm) in wurtzite GaN/ strained QR. Numerical results show that the BEF () increases linearly with increasing the hydrostatic pressure . For example, when the hydrostatic pressure increases from 0 to 8 GPa, the BEF () increases to 0.23 (0.024) MV/cm. For a given pressure and barrier thickness, Figure 2 also shows that the BEF () increases as the height decreases and the BEF remains insensitive to the variation of the hydrostatic pressure. This is caused by the change of the pressure-dependent piezoelectric constants, the biaxial strains, the dielectric constants, and the structural parameters of wurtzite GaN/ strained QR. According to (11), the equilibrium lattice constant of the strained layer increases with decreasing the QR height . Thus, the absolute value of the strain tensor of the ring layer along the QR growth direction increases, which induces the fact that the strength of the built-in electric field increases with the decrement of the height ; this can be understood based on (3).

In Figure 3, the ground-state donor binding energy in cylindrical wurtzite GaN/ QR is shown as a function of the radial thickness with the parameters  nm,  nm, and  nm for different hydrostatic pressures ( Gpa, 2 Gpa, 4 Gpa). The impurity is placed at the center of the QR. Set and . As shown in Figure 3, the donor binding energy increases with decreasing the radial thickness () in all cases, reaches a maximum value, and then decreases sharply. The behavior is related to the variation of the electron confinement in QR; the electron wave function is more firmly localized inside the QR with the decrease of the QR radial thickness. Thus, the Coulomb interaction between the electron and the impurity ion is enhanced and the donor binding energy increases correspondingly. However, when the radial thickness decreases continuously to a certain value, the kinetic energy of the confined electron raises greatly, which increases greatly the probability of the electron penetrating into the potential barrier by the uncertainty principle. Therefore, the donor binding energy starts decreasing quickly due to the weakness of the geometrical confinement in the radial direction. Moreover, Figure 3 also displays that the stronger the applied hydrostatic pressures is, the bigger the donor binding energy is. It can be explained that the QR sizes become small gradually with the larger hydrostatic pressure, and much stronger quantum confinement effect results in an increase of the donor binding energy.

In Figure 4, we present the ground-state donor binding energy in cylindrical wurtzite GaN/ strained QR as a function of the inner radius with the parameters ( nm,  nm, and  nm) for different hydrostatic pressures ( Gpa, 1 Gpa, 2 Gpa). The impurity is placed at the center of the QR. Set and . For a fixed value of the outer radius, the QR radial thickness () decreases as the inner radius augments; hence, the size quantization becomes stronger which leads to an increase in the donor binding energy. It can be understood that the expectation value of the electron-impurity distance along the radial direction decreases with the strengthening of the size quantization, and the probability of the electron localized near the impurity ion increases greatly; therefore, the donor binding energy becomes larger. The curves in Figure 6 also show that the donor binding energy increases quickly with the inner radius  nm. It is because the localization effect of the electron around the impurity ion is caused mainly by the radial confining potential and the influences of the hydrostatic pressure on the Coulomb interaction between the electron and the impurity ion are weaker; therefore, the donor binding energy increases gradually as the radial thickness becomes narrow. In addition, for a given structure parameters of the QR the stronger the applied hydrostatic pressure is, the bigger the donor binding energy is. This can be explained by the fact that a growth in the hydrostatic pressure leads to the increase in the quantum localization effect of the electron wave function.

To clarify the effect of the applied hydrostatic pressure on the ground-state donor binding energy, we investigated the ground-state donor binding energy in cylindrical wurtzite GaN/ strained QR with the parameters ( nm,  nm,  nm, and , ) and several values of the QR height . From Figure 5, one can observe that the donor binding energy increases almost linearly as the hydrostatic pressure increases. The behavior can be explained as follows. Firstly, the changes of the pressure-dependent piezoelectric constants, biaxial strains, and dielectric constants can cause the modification of the polarization in the ring layer, which leads to a significant increase of the BEF . Secondly, the bigger hydrostatic pressure leads to an increase in the electron effective masses, which can cause the increment of the electron confinement potential in -direction. At last, the bigger hydrostatic pressure can also lead to the decrement of the expected value of the distance between the electron and the impurity ion. Therefore, the larger the applied hydrostatic pressure is, the stronger the Coulomb interaction between the electron and the impurity is, and the donor binding energy increases correspondingly. Taking the solid curve for example, for the QR height  nm, the donor binding energy increases 24 meV approximately if the hydrostatic pressure increases from 0 to 8 GPa. Thus, the applied hydrostatic pressure has an important influence on the donor binding energy.

Figure 6 displays the ground-state donor binding energy as a function of cylindrical wurtzite GaN/ strained QR height with the parameters ( nm,  nm,  nm, and , ) and several values of the applied hydrostatic pressure. We can see from Figure 6 that the donor binding energy increases monotonically as the height decreases from 7 nm to 1 nm in all case. As expected, the electron wave function is much strongly compressed in the QR for the less height of the QR, and the electron-impurity interaction becomes stronger, which leads to the enhancement of the binding energy correspondingly. Figure 6 also shows that the donor binding energy increases with the increment of the hydrostatic pressure . This is because the quantum localization effect of the electron wave function in wurtzite GaN/ strained QR is strengthened with the increase of the applied hydrostatic pressure. Therefore, the decrement of the mean relative electron-impurity distance leads to the increase in the donor binding energy with the same spatial confinement.

In Figure 7, the ground-state donor binding energy is investigated as a function of the impurity position along the QR growth direction with the parameters ( nm,  nm,  nm, and ) and several values of the applied hydrostatic pressure ( Gpa). It is seen that the curves in Figure 7 are absolutely asymmetric, and the donor binding energy firstly increases, reaches a maximum value, and then decreases gradually when the impurity is moved from the plane to the plane along the QR growth direction. The maximum value of the donor binding energy is not located at the point and shifted towards positive -direction. It is because the topology structure of the QR and the strong BEF modify the spread of the electron wave function in the QR. The direction of the built-in electric field in the ring layer is opposite to the QR growth direction. Thus, the built-in electric field pushes the electron towards positive -direction. This behavior is in agreement with the result of [19]. In addition, the larger the applied hydrostatic pressure is, the stronger the localization effect of the electron wave function is with the same parameters (, , and ), so that the peak value of the donor binding energy increases accordingly. Therefore, the distribution of the electron wave function is not central symmetrical about the QR in presence of the strong BEF.

4. Conclusions

With the framework of the effective mass approximation, the ground-state donor binding energies in cylindrical wurtzite GaN/ strained QR are investigated theoretically in the presence of built-in electric field and hydrostatic pressure by means of a variational approach. The ground-state donor binding energy depends strongly on ring geometry, applied hydrostatic pressure, and impurity position in the finite confinement potential. Numerical results show that the donor binding energy increases firstly, reaches a maximum value, and then drops quickly as the radial thickness of the QR decreases. The donor binding energy is an increasing (a decreasing) function of inner radius (height). In addition, the donor binding energy has a maximum when the impurity ion is moved along the -axis of the QR from the bottom of the QR to the top. Moreover, the stronger the applied hydrostatic pressure is, the bigger the peak value of the impurity binding energy is, with the same spatial confinement. And the position of the peak value of the binding energy is also shifted towards positive -direction. The electronic wave function distribution in the QR is also obviously modified by the applied hydrostatic pressure. These theoretical results obtained in this paper are useful for design of some photoelectric devices constructed based on GaN/ QR structures.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported by the Scientific and Technological Department Foundation of Hebei Province (no. 12210617) and the Natural Science Foundation of Hebei Province (no. A201420308).