Abstract

Within the framework of the effective mass approximation, barrier thickness and hydrostatic pressure effects on the ground-state binding energy of hydrogenic impurity are investigated in wurtzite (WZ) GaN/AlxGa1−xN strained quantum dots (QDs) by means of a variational approach. The hydrostatic pressure dependence of physical parameters such as electron effective mass, energy band gaps, lattice constants, and dielectric constants is considered in the calculations. Numerical results show that the donor binding energy for any impurity position increases when the hydrostatic pressure increases. The donor binding energy for the impurity located at the central of the QD increases firstly and then begins to drop quickly with the decrease of QD radius (height) in strong built-in electric fields. Moreover, the influence of barrier thickness along the QD growth direction and Al concentration on donor binding energy is also investigated. In addition, we also found that impurity positions have great influence on the donor binding energy.

1. Introduction

Recently, much attention has been paid to wide-band gap wurtzite (WZ) GaN/AlGaN quantum heterostructures due to their promising applications in optoelectronic devices such as light-emitting devices (LEDs) and laser diodes (LDs) [13]. Doping impurities in GaN-based confined systems is an effective method for controlling the electronic and optical properties of optoelectronic devices. It is well known that the hydrostatic pressure applied on a GaN-based semiconductor 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 induced by the spontaneous and piezoelectric polarizations also affects obviously the distribution of the carrier wave function in WZ nitride-based quantum heterostructures. According to the above characteristics, various studies concerning impurity states are reported in GaN-based nanostructures such as quantum wells [48], quantum well wires [9, 10], and (double) quantum dots [1117]. In these references, considering the strong built-in electric field and/or hydrostatic pressure effects, the hydrogenic donor impurity states are mainly discussed in a radial infinite confinement potential barrier. Their results demonstrate that quantum structure size, strong built-in electric field, and hydrostatic pressure have a significant influence on the donor impurity binding energies, but there are few reports involved in barrier thickness effects on the donor binding energy in WZ GaN/AlxGa1−xN strained QD in finite potential barrier to date.

To further demonstrate the barrier thickness effect on impurity states in a WZ GaN/AlGaN strained QD, we calculate the donor binding energy of hydrogenic impurity in a WZ GaN/AlGaN QD under a strong built-in electric field by means of a variational approach in the finite confinement potential. In our calculation, effective mass of the electron, dielectric constants, phonon frequencies, energy gaps, sizes (radius and height) of QD, and piezoelectric polarizations are considered as a function of hydrostatic pressure. The paper is organized as follows. In Section 2, we describe the theoretical framework. Then, the pressure and the strain coefficients of GaN and AlxGa1−xN are discussed in Section 3. In Section 4, the numerical results are discussed. Finally, conclusions are drawn from the present study in Section 5.

2. Theoretical Framework

In Figure 1, the schematic view of a cylindrical WZ GaN/AlxGa1−xN QD is depicted, with a detailed description of the different dimensions of the QD (dot radius , height , and barrier thickness ). Additionally, the GaN quantum dot is embedded in an AlxGa1−xN host matrix material, and the -axis is defined to be the growth direction. Within the frame of the effective mass approximation, the Hamiltonian for a hydrogenic donor impurity in cylindrical WZ GaN/AlxGa1−xN QD under the influence of hydrostatic pressure can be written [4] aswhere 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 AlxGa1−xN materials. The Hamiltonian is given by [11]where and are the pressure- and strain-dependent effective masses of electron along and perpendicular to the -direction. And is the electron confinement potential due to the band offset () and is given by [13]In (2), F(p) is the pressure-dependent built-in electric field (BEF) in the finitely thick barrier layer for WZ GaN/AlxGa1−xN QD. The values of the BEF along the growth direction in the well () and barrier () that result from the difference in the total electric polarizations in each region are given by simple formulas [14]:where , and , are the thicknesses and electronic dielectric constants for GaN and AlxGa1−xN materials, , and , are the spontaneous and piezoelectric polarizations for GaN and AlxGa1−xN materials, and the bowing parameter is chosen as −0.019 C/m2 [14]:The wave function of the electron confined in the WZ GaN/AlxGa1−xN QD can be written as where the constants ( and ) are determined by the continuity of the derivative of the radial wave function at the QD 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 . Wave function can be expressed by means of the Airy functions Ai and Bi [16]: Here . The coefficients and can also be obtained by the transfer matrix methods [16].

In order to calculate the donor binding energy, the trial wave function can be chosen as [17]where denotes 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 impunity and the ground-state energy of the system with impurity; that is,

3. 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 [18]:where the equilibrium lattice constant for the strained layer under hydrostatic pressure depends on the lattice constants of the component materials and is weighted by their relative thicknesses [19]:The lattice constant of the GaN (AlN) material dependence of hydrostatic pressure satisfies [20] Here the actual lattice constant of the AlxGa1−xN material can be obtained by the linear interpolation method from the corresponding values of GaN and AlN. is the thickness of the dot (barrier) layer. The strain induced by the biaxial lattice mismatch along the z-direction in the heterostructure is [21]The coefficient of the GaN (AlN) material is given by [21]where is the pressure-dependent elastic stiffness constant of material and is given by [21] and the coefficient and the pressure-dependent elastic stiffness constant of AlxGa1−xN 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 [22]where , , , and ( and AlN) are the deformation potentials. The energy gap dependence of hydrostatic pressure is considered by the following equation [22]:The energy gap of the AlxGa1−xN alloy can be calculated by the following formula [14]:The biaxial strain and hydrostatic pressure dependences of the electron effective masses in the plane and the -direction can be calculated by [23]where is a fixed value for a given material and can be derived from the values of and . In (5), under consideration of the biaxial strain in a wurtzite structure modified by hydrostatic pressure, the dielectric constants and the phonon eigen-frequencies will also change. The static dielectric constant is influenced by biaxial strain and hydrostatic pressure. The tensor components of for the WZ structure are derived from the generalized Lyddane-Sachs-Teller relation [24]:where the LO- and TO-phonon frequencies influenced by biaxial strain and hydrostatic pressure can be written asFurthermore, the hydrostatic pressure dependence of can be given by [25]where the subscript represents LO or TO phonon, denotes plane or -direction, is the zone-center phonon frequency of material , is Grüneisen parameter of phonon mode given in [25], is bulk modulus, and and are the strain coefficients of zone-center phonon modes. Following [14], the influence of hydrostatic pressure on the high frequency dielectric constants can be written asHere ( GaN and AlN) is Phillips ionicity parameter of material; the static bulk modulus under hydrostatic pressure is given by [25]where , , , and are the elastic constants of material . The effective mean relative dielectric constant in (1) is defined as [26] Then the hydrostatic-pressure-modified biaxial and uniaxial strain dependence of the static dielectric constant is fully considered, whereas the dielectric constant of AlxGa1−xN can be calculated by the SCPA [13]. The piezoelectric polarization along the oriented WZ GaN/AlxGa1−xN QD can be calculated as [27]where and are the pressure-dependent piezoelectric constants of material .

4. Results and Discussions

Under the strong built-in electric field induced by the spontaneous and piezoelectric polarizations, barrier thickness and hydrostatic pressure effects on the donor binding energy of hydrogenic impurity are investigated in WZ GaN/AlxGa1−xN strained QD. All material parameters used in our calculations are listed in Tables 14. Material parameters of AlxGa1−xN are estimated using a linear interpolation between the values of the corresponding parameters of GaN and AlN.

Figure 2 presents that the built-in electric field (BEF) () in the well (barrier) layer along the growth direction -axis as a function of the barrier thickness in WZ GaN/AlxGa1−xN QD for different hydrostatic pressures . Numerical results show that the BEF () is an increasing (a decreasing) function of the barrier thickness . This is because the fact that the equilibrium lattice constant of the strained layer and the components of the strain tensor of the dot layer decrease with the increase of the barrier thickness according to (11), which induces the decrement of the piezoelectric polarization in the dot layer and the increment of the piezoelectric polarization in the barrier layer along the QD growth direction. Therefore, the BEF () gradually increases (decreases) (see (5)). Moreover, Figure 2 also shows that the built-in electric fields ( and ) increase correspondingly as the hydrostatic pressure increases; this is caused by the change of the pressure-dependent piezoelectric constants, the biaxial strains, and the dielectric constants of WZ GaN/AlxGa1−xN QD when the hydrostatic pressure increases.

In Figure 3, the built-in electric field (BEF) () in the well (barrier) layer along the QD growth direction -axis is displayed as a function of the barrier thickness for different Al compositions in WZ GaN/AlxGa1−xN QD. Numerical results show that the BEF () increases (decreases) gradually when the barrier thickness increases. This is caused by the change of the pressure-dependent piezoelectric constants, the biaxial strains, and the dielectric constants of WZ GaN/AlxGa1−xN strained QD. In addition, Figure 3 also shows that the BEF () increases as Al composition increases, but the BEF remains insensitive to the bigger barrier thickness . The reason can be given as follows. When Al composition increases, the equilibrium lattice constant of the strained layer decreases according to (11); therefore, the absolute values of strain tensor and of the barrier layer along the QD growth direction increase, which induces that the piezoelectric polarization in the barrier layer increases, and then the BEF () increases correspondingly in (4).

In Figure 4, the ground-state donor binding energy is shown as a function of the QD radius in WZ GaN/AlxGa1−xN QD with the parameters ( nm,  nm, and ) for different hydrostatic pressures ( Gpa, 4 Gpa, and 8 Gpa). The impurity ion is placed at the center of the QD. As shown in Figure 4, the donor binding energy increases with decreasing the radius in all cases, reaches a maximum value, and then decreases rapidly. The behavior is related to the variation of the electron confinement in QD; the electron wave function is firmly localized inside the QD with decreasing the QD radius. Therefore, the Coulomb interaction between the electron and the impurity ion is enhanced, and the donor binding energy increases correspondingly. Moreover, when the radial thickness decreases continuously to a certain value, the kinetic energy of the confined electron rises greatly, which increases greatly the probability of the electron penetrating into the potential barrier by the uncertainty principle and therefore the donor binding energy starts decreasing quickly. Moreover, Figure 4 also displays that the larger the hydrostatic pressure is, the bigger the donor binding energy is. The main reasons can be given as follows. With the increase of the hydrostatic pressure, the dielectric constants, the electron effective mass, and the finite confinement potential barrier will increase, which will result in bigger donor binding energy.

Figure 5 demonstrates that the ground-state donor binding energy in cylindrical WZ GaN/AlxGa1−xN strained QD as a function of the height with the parameters ( nm,  nm, and ) for different hydrostatic pressures ( Gpa, 4 Gpa, and 8 Gpa). The impurity ion is placed at the center of the QD, setting and . As expected, in all cases of hydrostatic pressures, the donor binding energy increases with a decrease of the QD height, reaches a maximum value, and then decreases quickly in finite potential barrier. For a fixed value of the barrier thickness , the size quantization confinement of the electron wave function goes stronger with the decrease of the QD height, the electron-impurity Coulomb interaction becomes larger, and when the QD height becomes small enough, the probability of the electron leaking into the potential barrier increases greatly by the uncertainly principle. Therefore, the donor binding energy decreases correspondingly. In addition, the curves in Figure 5 also show that the stronger the applied hydrostatic pressure is, the bigger the donor binding energy is. Take the QD height  nm for example; a change of the hydrostatic pressure from 0 to 8 Gpa results in an increase of the impurity binding energy from 34.03 meV to 48.35 meV. As expected, with the increase of the hydrostatic pressure , electron effective masses and dielectric constants of GaN and AlxGa1−xN materials become lager, and a growth of the hydrostatic pressure leads to the increase of the finite potential barrier, which will lead to the electron wave function being firmly squeezed around the impurity ion, and consequently the donor binding energy increases correspondingly.

The ground-state donor binding energy as a function of Al composition in WZ GaN/AlxGa1−xN QD is displayed in Figure 6 for different hydrostatic pressures . Numerical results show that the donor binding energy for the central impurity in WZ GaN/AlxGa1−xN strained QD increases with the increase of Al composition. This is because that the competition effects between the built-in electric field and the potential barrier confinement will change the strength of the electron-impurity interaction. For the small QD height  nm, as the Al concentration increases, conductor band offset of WZ GaN/AlxGa1−xN QD increases, and the potential barriers on the surfaces of QD play the main role in the distribution of the electron wave function. In addition, the bigger the Al composition is, the larger the potential barrier is, which results in the fact that the probability of the electron leaking into the barrier region becomes small. Accordingly, increasing the Coulomb effect between the electron and the impurity ion leads to the enhancement of the binding energy. Figure 6 also shows that the donor binding energy increases with the increment of the hydrostatic pressure . This is due to the fact that the electron wave function is strongly compressed in the QD as hydrostatic pressure increases, and the strength of the electron-impurity interaction becomes larger, leading to the enhancement of the binding energy correspondingly.

To clarify the effect of the hydrostatic pressure on the ground-state donor binding energy, we investigated the donor binding energy in cylindrical WZ GaN/AlxGa1−xN strained QD with the parameters ( nm, ,  nm, and  nm) for several values of QD height and barrier thickness. From Figure 7, one can observe that the donor binding energy increases almost linearly as the hydrostatic pressure increases. The pressure behavior in WZ GaN/AlxGa1−xN QD can be explained by the modification of the polarization in the dot layer by the pressure induced strain, which leads to a significant increase of the BEF in the QD. This behavior is in agreement with the result of [13]. In addition, the stronger the applied hydrostatic pressures is, the bigger the electron effective masses and dielectric constants of GaN and AlxGa1−xN materials are, and the finite confinement potential at the boundary of GaN QD also becomes large under the bigger hydrostatic pressure; hence the expected value of the distance between the electron and the impurity ion reduces and the strength of the electron-impurity interaction becomes larger, which will lead to the increase of the donor binding energy correspondingly. Taking the solid curve for example, the donor binding energy increases by 21.15 meV approximately if the hydrostatic pressure increases from 0 to 8 GPa. Thus, the hydrostatic pressure has an important influence on the donor binding energy.

Figure 8 displays the ground-state donor binding energy as a function of the barrier thickness along the QD growth direction with the parameters ( nm,  nm,  nm, and  nm) and different values of the hydrostatic pressure . The impurity ion is located at the centre of the QD. We can see from Figure 8 that the donor binding energy increases, reaching a maximum value, and then reduces gradually with the increase of the barrier thickness in all cases. As expected, the cures in Figure 8 also show that the donor binding energy has a maximum value. This is because the enhancement of the barrier thickness leads to the change of the finite confinement potential along the QD growth direction. As the barrier thickness increases, the finite confinement potential at the bottom of the QD becomes small, and the one at the top of the QD becomes big due to the strong built-in electric field effects. When the QD barrier thickness increases to about 11 nm, the finite confinement potentials at two sides of the QD along the growth direction are approximately equal, which leads to the fact that the electron wave function is strongly compressed around the central impurity ion. Therefore, the Coulomb action between the electron and the impurity ion reaches a maximum.

In Figure 9, the ground-state donor binding energy is investigated as a function of the impurity position along the QD growth direction with the parameters ( nm,  nm, and  nm) and different values of the hydrostatic pressure . As shown in Figure 9 the curves in all cases are absolutely asymmetry, and the donor binding energy demonstrates a maximum value when the impurity is located from the plane to the symmetry plane along the growth direction of the QD, and the maximum value of the donor binding energy is not located at the point . This is because the fact that the strong BEF modifies the spread of the electron wave function in the QD, and the direction of the built-in electric field in the dot layer is opposite to the growth direction of the QD. Thus, the built-in electric field pushes the electron toward the right side of the QD. This behavior is in agreement with the result of [17]. In addition, the curves also show that the stronger the hydrostatic pressure is, the larger the peak value of the donor binding energy is with the same parameters (, , and ). Moreover, the position of the peak value of the donor binding energy is also shifted to positive -direction. This is because the fact that the electronic wave function is obviously modified, and the bigger concentration of the electron wave function is squeezed strongly around the impurity ion. In addition, the stronger the hydrostatic pressure is, the bigger the localization effect of the electron wave function is, so that the peak value of the binding energy increases accordingly. Therefore, the distribution of the electron wave function is not central symmetrical about the QD in the presence of the strong BEF.

5. Conclusions

With the framework of the effective mass approximation, the ground-state donor binding energy in a cylindrical WZ GaN/AlxGa1−xN strained QD is investigated theoretically in the presence of built-in electric field and hydrostatic pressure by using a variational approach. The ground-state donor binding energy depends strongly on dot radius, hydrostatic pressure, impurity position, and barrier thickness in the finite confinement potential. Numerical results show that the donor binding energy increases firstly, reaches a maximum value, and then drops slowly as the QD radius (height) decreases. And the donor binding energy is an increasing function of Al composition and/or hydrostatic pressure. In addition, the donor binding energy has a maximum value, when the impurity position moves along the symmetry axis of the QD from the bottom of the QD to the top, and the position of the peak value of the donor binding energy is also shifted towards positive -direction. Moreover, the stronger the hydrostatic pressure is, the larger the peak value of the donor binding energy is with the same spatial confinement. The electronic wave function distribution in the QD is also obviously modified by the hydrostatic pressure. We hope that our results would stimulate further researches and lead to some potential applications on group-III nitride materials.

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).