The Effect of Dislocation and Interface-Roughness Scattering on Electron Mobility in the MgZnO/ZnO Heterostructure

Liu, DongFeng

doi:https://doi.org/10.1155/2022/2770047

Advances in Materials Science and Engineering

On this page

Abstract Introduction Results and Discussion Conclusion Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2022 | Article ID 2770047 | https://doi.org/10.1155/2022/2770047

The Effect of Dislocation and Interface-Roughness Scattering on Electron Mobility in the MgZnO/ZnO Heterostructure

DongFeng Liu¹

Academic Editor: Samia A. Gad

Received24 May 2022

Revised18 Jun 2022

Accepted30 Jun 2022

Published02 Aug 2022

Abstract

In this work, the electron mobility in the MgZnO/ZnO heterostructure at room temperature is theoretically studied by considering interface roughness (IFR), dislocation (DIS), and polar optical phonon (POP) scattering. Analytical formulae are introduced to calculate the critical thickness and dislocation density in the barrier layer of MgZnO/ZnO heterostructures. The calculated critical thickness for the MgZnO/ZnO heterostructure is much smaller than that for the AlGaN/GaN heterostructure system. At room temperatures, POP scattering is found to be the most important scattering mechanism. On the other hand, the change of electron mobility limited by IFR as a function of the barrier thickness in the MgZnO layer is found to be quite different to that limited by DIS. High-density (>) 2DEG can be obtained in the MgZnO/ZnO interface by increasing the thickness and Mg composition in the MgZnO layer.

1. Introduction

The II-VI semiconductor ZnO, a wide band-gap metal oxide material, has been proposed as new wide band-gap semiconductors for short-wavelength optoelectronic applications. ZnO and MgZnO can form heterostructures with polarization at the interface [1–4], which causes very high two-dimensional electron gas (2DEG). High 2DEG is very crucial for the realization of MgZnO/ZnO high electron mobility transistor (HEMT) applications [5]. The MgZnO/ZnO heterosystem is an effective candidate to confine carriers in the ZnO layer [6]. Recently, high electron mobility exceeding in MgZnO/ZnO heterostructures was reported [7].

The electron mobility in heterostructures is limited by a combination of various scattering mechanisms related to ionized impurities (II), dislocations (DIS), interface roughness (IFR), acoustic phonons (AP), alloy disorder (ADO), and polar optical phonon (POP) scatterings. Begum et al. [8] investigated the role of AP and POP scatterings in the electron mobility at the MgZnO/ZnO heterojunction. Recently, Zan and Ban theoretically analyzed the electron mobility limited by interface optical (IF) and confined optical (CO) phonon modes in a symmetric MgZnO/ZnO quantum well structure [9]. Thongnum et al. reported the role of IFR in electron transport in MgZnO/ZnO [2], finding the scattering potential in the in-plane direction which can affect the electron mobility. GaN can have a dislocation density up to , which may dominate the electron mobility in AlGaN/GaN heterostructure systems [10, 11].

But in MgZnO/ZnO heterostructure systems, IFR and DIS scattering limitation to electron mobility is not fully investigated. In this work, we first deduced expressions to analytically calculate the critical thickness and dislocation density of in-barrier layer with a different Mg composition x so that the electron mobility limitation due to DIS and IFR scatterings is analytically investigated.

2. Calculation Methods

In the effective mass approximation, the energy levels and the corresponding wave functions for electrons in the heterostructure are obtained by solving the Schrodinger and Poisson equations self-consistently. For more details of the procedure used in the self-consistent numerical solutions, we can see [12, 13] and references therein.

In the case of unintentionally doping, the 2DEG in heterostructures is assumed to originate from the surface donor distribution [14]. The positive surface donor sheet density is determined by the following [15]:

2DEG density is equal to due to the charge neutrality condition. Here, t is the barrier thickness, is the critical barrier thickness, and is polarization-induced charges at the surface . Piezoelectric polarization is expressed by [16]:where is the lattice mismatch, which is determined bywhere and are the lattice constants of ZnO and MgZnO, respectively. The critical thickness is calculated according to the expression [17]:where is the Poisson ratio. B is the magnitude of the Burgers vector. The strain-related dislocation density can be expressed as [18]where is the material constant. is the strain, which can be estimated by [19]where , , and are Pauli covalent radii of Zn and O atoms, and x is the Mg composition.

Considering that the material system of MgZnO/ZnO is analogous to that of AlGaN/GaN, it is helpful to prepare the DIS scattering calculation based on the data from AlGaN/GaN. The real AlGaN/GaN HEMTs contain a high density of dislocations in the epitaxial layers of GaN, which forms because of the absence of a suitable lattice-matched substrate for growth. State-of-the-art AlGaN/GaN HEMTs have a defect density of the order of 10⁸–10¹⁰/cm². At the same time, when the thickness of the Al_xGa_1-xN alloy layer exceeds the critical limit, the interface between the AlGaN barrier and the GaN channel suffers from crystal defects, resulting in misfit dislocation [20]. The effect of dislocation on electron transport can be analyzed by introduction of dislocation scattering. The corresponding scattering rate is obtained as [21].where is the Fermi wave vector, which depends on the 2DEG carrier concentration . is the dislocation density. is the 2D Thomas Fermi wave vector, and is the effective Bohr radius in the material. is the line charge density, where is the lattice constant in the (0001) direction of ZnO. The 2-DEG mobility inhibited by dislocation scattering can be given by .

The scattering rate for electrons being scattered from interface roughness (IFR) is obtained as [22].where is the average height of the roughness, and is the correlation length of the roughness.

, is the scattering angle, and is the screening constant,where is the form factor defined by the envelope wave function

The 2-DEG mobility inhibited by interface roughness scattering can be given by .

The polar optical phonon scattering between subband i to subband j is calculated as follows [23]:where are initial and final wave vectors, respectively. is the high frequency dielectric constant, and is the static dielectric constant of the system. The electron mobility inhibited by polar phonon scattering is calculated using the expressionwhere is the Fermi–Dirac equilibrium distribution function, is the electron kinetic energy, and N is the number of the subbands.

The total mobility then is given by Matthiessen’s rule as follows:

3. Results and Discussion

The band diagram, eigenvalues, and subband wave functions are calculated using the material parameters adopted from [9] for MgO and ZnO, listed in Table 1. Mg_xZn_1-xO material constants are determined by the interpolation calculation of MgO and ZnO material constants. The main parameters used for calculating DIS, POP, and IFR scattering rates are listed in Table 2. Because of the quantum confinement, electrons fly and are scattered mainly in the channel, so we only present the material parameters of ZnO in Table 2.

Figure 1 shows the conduction band edge profiles of a MgZnO/ZnO heterostructure which consists of an undoped MgZnO layer and a thick undoped ZnO buffer layer. The envelope wave functions and subband structures are calculated by self-consistently solving the one-dimensional Poisson–Schrodinger equations. We use a Schottky barrier height of 0.7 eV for the metal-semiconductor interface [24]. Here, the Mg composition of the MgZnO layer is fixed as 0.3. The thickness of the MgZnO and ZnO layers is 30 nm and 400 nm, respectively. The energy difference between the ground subband and the first excited subband is approximately 0.17 eV; therefore, the electrons mainly occupy the ground subband.

(a)

(b)

The critical thickness as a function of Mg composition x is calculated according to Equation (4). Poisson ratio of Mg_xZn_1-xO is fixed as 0.351 according to the work [27]. Here, the magnitude of the Burgers vector [28]. is the lattice constant of MgZnO, which is calculated from the corresponding values of MgO and ZnO for the given Mg mole fraction x using linear relation. Figure 2 shows the calculated critical thickness as a function of Mg composition x. For a specific x, the critical thickness for the MgZnO/ZnO heterostructure is found much smaller than that for the AlGaN/GaN heterostructure.

Various design constraints exist to change the 2DEG density. For example, one may change the barrier thickness . Figure 3 shows the MgZnO layer thickness dependence of carrier sheet density. We first choose x to be 0.30 so that all layer thicknesses ranging from 10 nm to 30 nm are larger than the corresponding critical thickness. We find that the 2DEG density increases with increasing the barrier thickness, which is consistent with the experimental data of AlGaN/GaN heterostructure systems [15, 17] and theoretical data for MgZnO/ZnO heterostructure systems [29]. We can see that electron sheet density can reach a value of at barrier thickness 10 nm, at 14 nm, at 30 nm, respectively. We then change x to be 0.35. The corresponding results are shown by the dashed line in Figure 3. It is found that for a specific thickness, higher Mg mole fraction in the MgZnO barrier layer can result in higher 2DEG density. Hence, high-density 2DEG can be obtained in MgZnO/ZnO heterostructures by changing the barrier thickness and the Mg composition.

Dislocation (DIS) scattering plays important role in electron transport and the resulted mobility. We know from Equation (7) that the DIS scattering rate is closely related with the dislocation density. Referring to the experimental data of dislocation density in [30], , we determined the material constant in Equation (5), . Based on Equation (5), we draw in Figure 4 the calculated results of dislocation densities in the MgZnO barrier layer as a function of Mg composition x.

Since the electron transport in HEMT is closely related with the band offset of heterostructures, we first investigate the electron mobility as a function of the Mg composition x in the Mg_xZn_1-xO layer. Correlation length 6.0 nm and averaged height = 0.22 nm are used for the calculation of IFR scattering rates. The MgZnO and ZnO layer thicknesses are fixed as 40 nm and 60 nm, respectively. We only consider the mobility limitations by IFR, POP, and DIS scatterings, so the total electron mobility is calculated by Equation (14). The dislocation densities used for the calculation of DIS scattering rates are chosen from Figure 4. The corresponding results are displayed in Figure 5. The solid line in Figure 5(a) is for the IFR mobility , and the dashed line is for the total mobility. For comparison, the electron mobility limited by DIS scattering is displayed in the inset. This gradual decrease for the mobility with increasing Mg composition is consistent with those in [24]. On the other hand, we can see from Figure 5(a) that is close to , which means that IFR scattering is dominant over DIS and can describe the electron transport of MgZnO/ZnO heterostructures in the low Mg composition (x < 0.4). Because the dislocation scattering rate [18], and increases with x, so we see from the inset in Figure 5(a) that the DIS mobility increases at low Mg composition, and then it decreases with increasing x. POP scattering is the most important one for polar wide band-gap semiconductor materials in room temperature. For comparison, we present the POP-limited electron mobilities in Figure 5(b). We find that POP scattering is the dominant mechanism in limiting electron mobility at room temperature.

(a)

(b)

(c)

(d)

We reduce that were used in Figure 5(a), setting = 0.11 nm, and 6.0 nm, 8.0 nm, 10.0 nm, respectively. The corresponding calculated IFR electron mobilities are shown in Figure 5(c). We find that IFR mobility strongly depends on the average height of roughness. Here, the value used for average height is very small (0.11 nm), compared with the typical value of roughness (0.286∼3.38 nm) [31]. In Figure 5(d), we fix correlation length 6.0 nm and slightly vary the averaged height: 0.535 nm, 0.684 nm, and 0.832 nm. In these cases, the obtained IFR electron mobilities are still much smaller than DIS electron mobility, which further confirms IFR scattering as the domain mechanism for the electron mobility limitation in MgZnO/ZnO heterostructures.

Figure 6 shows the dependence of electron mobility on the Mg_0.35Zn_0.65O layer thickness. We find from Figure 6(a) that electron mobility is limited by IFR scattering drops with the increase of barrier thickness. This phenomenon is due to the fact that the electron sheet density increases with increasing MgZnO barrier thickness, as shown in Figure 3 and that the IFR scattering rate is proportional to , as presented in Equation (8). Figure 6(b) shows that the electron mobility limited by DIS scattering increases with increasing MgZnO barrier thickness, which displays a reverse behavior in contrast to IFR electron mobility shown in Figure 6(a). The result in Figure 6(a) is due to the decrease of the DIS scattering rate with increasing , as pointed out by Equation (7).

(a)

(b)

4. Summary and Conclusion

In this study, we theoretically investigated electron mobility limitations by the interface roughness, dislocation, and polar optical phonon scatterings in MgZnO/ZnO heterostructures at room temperature. Analytical formulae are introduced to calculate dislocation density in the barrier layer of heterostructures. The results confirmed that the strength of the electron-optical-phonon interaction in MgZnO/ZnO heterostructures is high enough as to constitute a strong mobility limiter. The 2DEG mobility as a function of the MgZnO barrier thickness in the MgZnO layer which is limited by interface roughness is found to be quite different to that limited by dislocation scattering; the former increases and the latter decreases with increasing the barrier thickness. Higher Mg mole fraction and a thicker MgZnO barrier layer can generate higher 2DEG density in MgZnO/ZnO heterostructures. The mobility limited by IFR is an order of magnitude higher than that of DIS, suggesting that IFR scattering, compared to DIS scattering, is the major limitation to electron mobility at room temperature. The work presented here could help in the design of future devices, such as high electron mobility transistors (HEMT).

Data Availability

Data available on request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the Science and Technology Program of Guangzhou, China (Grant no. 201804010444).

References

Bo-S. Wang, Y.-S. Li, and I. C. Cheng, “Mobility enhancement in RF-sputtered MgZnO/ZnO heterostructure thin-film transistors,” IEEE Transactions on Electron Devices, vol. 63, no. 4, pp. 1545–1549, 2016.
View at: Publisher Site | Google Scholar
A. Thongnum, V. Sa-yakanit, and U. Pinsook, “Two-dimensional electron transport in MgZnO/ZnO heterostructures: role of interface roughness,” Journal of Physics D: Applied Physics, vol. 44, no. 32, Article ID 325109, 2011.
View at: Publisher Site | Google Scholar
H.-A. Chin, I. C. Cheng, C. I. Huang et al., “Two dimensional electron gases in polycrystalline MgZnO/ZnO heterostructures grown by rf-sputtering process,” Journal of Applied Physics, vol. 108, no. 5, Article ID 054503, 2010.
View at: Publisher Site | Google Scholar
A. Gold, “Transport properties of the electron gas in ZnO/MgZnO heterostructures,” Applied Physics Letters, vol. 96, no. 24, Article ID 242111, 2010.
View at: Publisher Site | Google Scholar
Y. K. Verma, V. Mishra, and S. K. Gupta, “A physics-based analytical model for MgZnO/ZnO HEMT,” Journal of Circuits, Systems, and Computers, vol. 29, no. 1, Article ID 2050009, 2020.
View at: Publisher Site | Google Scholar
S. H. Park, W. P. Hong, J.-J. Kim, and D. Ahn, “Theoretical study of a two-dimensional electron gas in wurtzite ZnO/MgZnO heterostructures and comparison with experiment,” Journal of the Korean Physical Society, vol. 67, no. 10, pp. 1844–1847, 2015.
View at: Publisher Site | Google Scholar
J. Falson, Y. Kozuka, M. Uchida et al., “MgZnO/ZnO heterostructures with electron mobility exceeding 1 × 106 cm2/Vs,” Scientific Reports, vol. 6, no. 1, Article ID 26598, 2016.
View at: Publisher Site | Google Scholar
K. Begum, S. Kolvekar, and N. Sankeshwar, “Electron mobility in MgZnO/ZnO heterojunctions,” Proceedings of the International Conference on Trends in Applied Physics and Material Science, vol. 1536, pp. 447-448, 2013.
View at: Publisher Site | Google Scholar
Y. H. Zan and S. L. Ban, “Electronic mobility limited by optical phonons in symmetric MgxZn1-xO/ZnO quantum wells with mixed phases,” Superlattices and Microstructures, vol. 150, Article ID 106782, 2021.
View at: Publisher Site | Google Scholar
X. Xu, X. Liu, X. Han et al., “Dislocation scattering in AlxGa1−xN∕GaN heterostructures,” Applied Physics Letters, vol. 93, no. 18, Article ID 182111, 2008.
View at: Publisher Site | Google Scholar
F. Carosella and J.-L. Farvacque, “Effect of threading dislocations on carrier mobility in AlGaN/GaN quantum wells,” Journal of Physics: Condensed Matter, vol. 20, no. 32, Article ID 325210, 2008.
View at: Publisher Site | Google Scholar
D. Liu, Y. Cheng, and J. F. He, “Nanogroove control of intersubband absorption in AlGaAs/GaAs step quantum wells,” Semiconductor Science and Technology, vol. 31, no. 3, Article ID 035006, 2016.
View at: Publisher Site | Google Scholar
L. Ramdas Ram-Mohan, Finite Element and Boundary Element Application in Quantum Mechanics, Oxford University Press, Oxford, UK, 2002.
G. Koley and M. G. Spencer, “On the origin of the two-dimensional electron gas at the AlGaN∕GaN heterostructure interface,” Applied Physics Letters, vol. 86, no. 4, Article ID 042107, 2005.
View at: Publisher Site | Google Scholar
J. P. Ibbetson, P. T. Fini, K. D. Ness, S. P. DenBaars, J. S. Speck, and U. K. Mishra, “Polarization effects, surface states, and the source of electrons in AlGaN/GaN heterostructure field effect transistors,” Applied Physics Letters, vol. 77, no. 2, pp. 250–252, 2000.
View at: Publisher Site | Google Scholar
D. Liu, D. Lin, Z. Li, and K. Guo, “Electron mobility in abrupt-interface and step-graded AlGaN/GaN Heterostructures,” Semiconductor Science and Technology, vol. 32, no. 10, Article ID 105008, 2017.
View at: Publisher Site | Google Scholar
J. Qu, J. Li, and G. Zhang, “Heterostructure grown by metalorganic vapor phase epitaxy,” Solid State Communications, vol. 107, no. 9, pp. 467–470, 1998.
View at: Publisher Site | Google Scholar
G. K. Williamson and R. E. Smallman, “III. Dislocation densities in some annealed and cold-worked metals from measurements on the X-ray debye-scherrer spectrum,” Philosophical Magazine, vol. 1, no. 1, pp. 34–46, 1956.
View at: Publisher Site | Google Scholar
A. B. M. A. Ashrafi and Y. Segawa, “Determination of Mg composition in Mg[sub x]Zn[sub 1−x]O alloy: validity of Vegard's law,” Journal of Vacuum Science & Technology B: Microelectronics and Nanometer Structures, vol. 23, no. 5, p. 2030, 2005.
View at: Publisher Site | Google Scholar
J. A. Floro, D. M. Follstaedt, P. Provencio, S. J. Hearne, and S. R. Lee, “Misfit dislocation formation in the AlGaN∕GaN heterointerface,” Journal of Applied Physics, vol. 96, no. 12, pp. 7087–7094, 2004.
View at: Publisher Site | Google Scholar
D. Jena, A. C. Gossard, and U. K. Mishra, “Dislocation scattering in a two-dimensional electron gas,” Applied Physics Letters, vol. 76, no. 13, pp. 1707–1709, 2000.
View at: Publisher Site | Google Scholar
S. Gökden, R. Baran, N. Balkan, and S. Mazzucato, “The effect of interface roughness scattering on low field mobility of 2D electron gas in GaN/AlGaN heterostructure,” Physica E: Low-Dimensional Systems and Nanostructures, vol. 24, no. 3-4, pp. 249–256, 2004.
View at: Publisher Site | Google Scholar
A. Mátyás, “Modeling of Quantum cascade Lasers by Ensemble Monte-Carlo Methods,” Technische Universitat Munchen, München, Germany, 2011, Ph.D Dissertation.
View at: Google Scholar
S. H. Park, “Properties of the two-dimensional electron-gas of a hybrid MgZnO/InGaN/ZnO heterostructure with an InGaN channel layer,” Journal of the Korean Physical Society, vol. 75, no. 4, pp. 326–330, 2019.
View at: Publisher Site | Google Scholar
D. C. Look, D. C. Reynolds, J. R. Sizelove et al., “Electrical properties of bulk ZnO,” Solid State Communications, vol. 105, no. 6, pp. 399–401, 1998.
View at: Publisher Site | Google Scholar
G. Milano, L. D’Ortenzi, K. Bejtka et al., “Metal-insulator transition in single crystalline ZnO nanowires,” Nanotechnology, vol. 32, no. 18, Article ID 185202, 2021.
View at: Publisher Site | Google Scholar
Z. Liu, X. He, Z. Mei et al., “Polarity-manipulation based on nanoscale structural transformation on strained 2D MgO,” Journal of Physics D: Applied Physics, vol. 47, no. 10, Article ID 105303, 2014.
View at: Publisher Site | Google Scholar
S. Murugesan, P. Kuppusami, E. Mohandas, and M. Vijayalakshmi, “X-ray diffraction Rietveld analysis of cold worked austenitic stainless steel,” Materials Letters, vol. 67, no. 1, pp. 173–176, 2012.
View at: Publisher Site | Google Scholar
S. Ait Adi, A. Kadri, and K. Zitouni, “Combined ab-initio and Schrödinger-Poisson simulation study of spontaneous and piezoelectric polarizations effects on I-V characteristics of ZnO/MgxZn1-xO heterostructures,” Computational Condensed Matter, vol. 21, Article ID e00430, 2019.
View at: Publisher Site | Google Scholar
Z. Vashaei, T. Minegishi, H. Suzuki, M. W. Cho, and T. Yao, “Defect and interface studies of ZnO/MgxZn1−xO heterostructures,” Journal of Physics and Chemistry of Solids, vol. 69, no. 2-3, pp. 497–500, 2008.
View at: Publisher Site | Google Scholar
D. Zanato, S. Gokden, N. Balkan, B. K. Ridley, and W. J. Schaff, “The effect of interface-roughness and dislocation scattering on low temperature mobility of 2D electron gas in GaN/AlGaN,” Semiconductor Science and Technology, vol. 19, no. 3, pp. 427–432, 2004.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 DongFeng Liu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies