Abstract
GaN-based semiconductors are promising materials for solid-state optoelectronic applications. However, the strong internal electrostatic field (IEF) along the [0001] direction is a serious problem that harms the efficiency of lighting devices based on GaN-based semiconductors due to the quantum confined Stark effect. Here we theoretically predict a method, reducing the dimensions from bulk to two-dimensional (2D) structures, to fundamentally remove the IEF. After thinning the materials to several nanometers, the wurtzite configuration (with strong IEF) spontaneously transform to the haeckelite (4 | 8) configuration (without IEF) due to the more stable neutral surface in the 4 | 8 configuration. Meanwhile, the 4 | 8 configuration maintain optoelectronic properties comparable to or even better than those of the wurtzite configuration. By carefully analyzing the interaction between 2D GaN and different types of substrates (SiC and graphene), we not only provide clear physical insights for experimental results but also address a “thickness-controlled” vdW epitaxy scheme to experimentally realize the 4 | 8 configuration. We believe that the 4 | 8 configuration without IEF is a prospective material for diverse optoelectronic applications. In addition, we propose a point of view in engineering the properties of GaN-based semiconductors.
Similar content being viewed by others
Introduction
GaN-based semiconductors are highlighted for their excellent properties, such as tunable and direct band gap, high thermal conductivity, and good chemical stability1. The materials and devices have been extensively investigated in recent decades, which has led to significant progress in light-emitting diodes2,3,4,5,6, photodetectors7,8,9, and laser diodes10,11,12. However, the great success in lighting devices is mainly limited by the polar material growth along the c-axis [0001]. GaN-based semiconductors crystallize into the wurtzite structure with a c : a ratio lower than that for perfectly hexagonally close-packed atoms (c : a = 1.633 for perfect hcp, 1.601, 1.627, and 1.623 for AlN, GaN, and InN, respectively). This, together with the lack of inversion symmetry and the strong ionic nature of III-N bonds, produces a spontaneous dipole moment along the [0001] direction. The dipole moment causes opposite bound polarization charges on the metal-polar [0001] and N-polar [000–1] surfaces, which leads to a strong internal electrostatic field (IEF) along the [0001] direction13,14,15. The IEF separates electrons and holes spatially on either side and induces an energy band slope16, thus reducing the emission efficiency and leading to a redshift in lighting devices (quantum confined Stark effect)17,18. Therefore, eliminating the IEF in GaN-based materials is highly expected to improve the efficiency of lighting devices. Many efforts have been devoted to this. Doping, which may be the earliest attempt to screen the electrostatic field, fails, because the doping concentrations required to achieve the magnitude of the field are as high as 1019 cm−3 19,20. A promising way to overcome the fundamental challenges is to realize semi- or nonpolar materials that can sharply reduce or even eliminate the IEF. In 2000, Ploog and colleagues21 reported the molecular beam epitaxy (MBE) deposition technique by using LiAlO2 as a substrate to obtain nonpolar GaN. Atomically flat surfaces have been achieved by the classic two-step approach or a high-temperature AlN buffer method22. However, the high density of stacking faults makes the quality of the materials far from the requirements for device structures. Another intuitive alternative is to replace the wurtzite phase by the unpolarized zinc-blende phase of nitride semiconductors. However, the intrinsic thermodynamic metastability and the lack of ideal substrates hinder the realization of device-quality material even after extensive efforts23. The fundamental solution to this problem is to find a stable alternative material without IEF but that maintains optoelectronic performance comparable to that of wurtzite structures.
Recently, diverse two-dimensional (2D) GaN-based semiconductors were explored by both theoretical and experimental scientists24,25,26,27,28,29,30,31,32,33,34,35,36,37. Interestingly, many reported structures (including hexagonal24,25,26, haeckelite (4 | 8)27,28, and tetragonal29,34) were nonpolarized. Tsipas et al.32 obtained hexagonal AlN nanosheets by MBE on Ag(111) and confirmed the reduced band gap for hexagonal AlN in 2013. One year later, Malin et al.33 realized hexagonal AlN on the semiconductor substrate Si(111). However, the huge cost and uncontrollable quality of the sample obtained by MBE remarkably hindered the large-scale commercial application of 2D III-nitride materials. In the past 3 years, scientists from the Pennsylvania State University35 and South China University of Technology36,37 successfully synthesized 2D GaN and AlN by metal organic chemical vapor deposition with the assistance of graphene encapsulation. The exciting achievements significantly accelerated the steps of industrialization for 2D GaN-based semiconductors. Kolobov et al.28 predicted that graphitic and wurtzite few-layer GaN was unstable and spontaneously reconstructed into the haeckelite structure, in which the interlayer bonds along the c-axis exhibited a 4 | 8 motif (highlighted by yellow shadow), as shown in Fig. 1a. Notably, the [0001] and [000–1] surfaces with opposite charges in the wurtzite configuration are replaced by neutral surfaces in 4 | 8 configurations. The IEF should be reduced or even eliminated in these configurations. This is highly expected to increase the quantum efficiency of lighting devices based on GaN-based semiconductors.
First-principles calculations have obtained great success in the research of GaN-based semiconductors24,28,29,38,39. In this work, we theoretically demonstrated by first-principles calculations that in GaN-based semiconductors, the 4 | 8 configuration is more energetically favorable than the wurtzite configuration after reducing the dimensions from bulk to 2D structures. This supports a convenient method to obtain GaN-based semiconductor (represented by AlN, GaN, and InN in this work) samples without IEF. The structural phase transition of wurtzite (w) → haeckelite (4 | 8) → hexagonal (h) induced by surface energy is revealed, as shown in Fig. 1a. Compared with the wurtzite configuration, the 4 | 8 configuration possesses comparable optoelectronic performance but no IEF. Meanwhile, we also provide a feasible approach, “thickness-controlled” vdW epitaxy, to promote the experimental realization of the 4 | 8 configuration. We find that the surface atomic structure of the substrate and the substrate/epilayer interaction are two decisive factors determining the configuration of materials.
Results and discussions
Structural stability
In the bulk phase, the total free energies of the wurtzite configurations are smaller than those of the 4 | 8 configurations by 78, 110, and 48 meV for AlN, GaN, and InN, respectively, as the bond lengths and angles of the wurtzite configurations are closer to perfect tetrahedral coordination. The bond angles are all ~109.5 ± 1° in the wurtzite configurations, which are close to those of the ideal tetrahedral configuration. Meanwhile, one bond angle decreases to 90 ± 2° and the others increase to 112.5 ± 1° in the 4 | 8 configurations, which significantly deviates from the ideal tetrahedral configuration. Therefore, the 4 | 8 configurations endure serious distortions and the distortion will accumulate with increasing thickness. The lattice parameters of the wurtzite and 4 | 8 configurations are listed in Table 1. For 2D structures, the stability of the materials is contributed by both the bulk and surface parts. The calculated surface energies (defined in the methodology section) are as follows: \(\gamma _{{\it{\mathrm {AlN}}}}^w\) is 2.02 eV f.u.−1, \(\gamma _{{\it{\mathrm {GaN}}}}^w\) is 1.72 eV f.u.−1, \(\gamma _{{\it{\mathrm {InN}}}}^w\) is 0.91 eV f.u.−1, \(\gamma _{{\it{\mathrm {AlN}}}}^{{\it{4/8}}}\) is 0.93 eV f.u.−1, \(\gamma _{{\it{\mathrm {GaN}}}}^{{\it{4/8}}}\) is 0.74 eV f.u.−1, and \(\gamma _{{\it{\mathrm {InN}}}}^{{\it{4/8}}}\) is 0.46 eV f.u.−1. The \(\gamma _{\it{m}}^{{\it{4/8}}}\) values are generally much lower than the \(\gamma _{\it{m}}^w\) values regardless of the thickness (as shown in Supplementary Fig. 2), indicating that the surfaces of the 4 | 8 configurations are obviously more stable than those of the wurtzite configurations. The surfaces of the wurtzite configurations were classified as type III by Tasker and are intrinsically unstable due to the divergence of the surface energy (\(\gamma _{\it{m}}^c\)), as well as numerous dangling bonds40. The strong surface dipoles associated with the presence of oppositely charged cations or anions on each surface introduce instability into the wurtzite structures. These issues are considerably optimized in the 4 | 8 structures. As shown in Fig. 2b, the outermost atomic planes in 4 | 8 slabs are nearly flat, thus minimizing the surface dipoles and neutralizing the surface charge. However, with increasing thickness, the rise in energy induced by accumulated distortion becomes greater than the energy compensation of the surface energy, leading to the wurtzite configurations being more stable than the 4 | 8 configurations in thick films or bulk phases. According to our results, the 4 | 8 slabs show lower energy than wurtzite slabs when 4 < n ≤ 28 (n is the number of layers) for AlN, 2 < n ≤ 18 for GaN, and 4 < n ≤ 14 for InN, as shown in Fig. 1b, c, although the wurtzite structures are the most stable configurations when n > 28 for AlN, n > 18 for GaN, and n > 14 for InN. The wurtzite → 4 | 8 transition occurs when the thickness is reduced to several nanometers. In fact, the dangling bonds on the surface are quite unstable; thus, we consider the hydrogen-passivated surface structures following the same method and similar trends are observed as shown in Supplementary Fig. 3.
As the thickness is continuously reduced to a few atomic layers, the 4 | 8 configurations transform to hexagonal configurations. First, we extract buckled monolayers from bulk wurtzite or 4 | 8 structures as initial structures. After relaxation, we obtain planar hexagonal structures such as graphene or h-BN, as shown in Fig. 1a. To remove the dangling bonds and stabilize the surface, the tetrahedral coordination is reconstructed into triangle-planar coordination. However, the sp2 hybridization is still not stable enough. If two monolayers are stacked together to form a bilayer, the sp2 hybridization is slightly distorted and the interlayer distances are 2.14, 2.65, and 2.47 Å, which are much shorter than that in graphene (3.30 Å). Hence, the interlayer coupling is much stronger than the normal vdW interaction. The calculated \(\gamma _{{\it{\mathrm {AlN}}}}^h\) is 1.04 eV f.u.−1, \(\gamma _{{\it{\mathrm {GaN}}}}^h\) is 0.48 eV f.u.−1, and \(\gamma _{{\it{\mathrm {InN}}}}^h\) is 0.86 eV f.u.−1, indicating the quasi-bond character between the layers of hexagonal structures. According to our calculations, the planar hexagonal structure can only be stable for at most two layers of GaN and InN, and four layers of AlN.
Optoelectronic properties
First, we employ the ten-layer slabs of 2D GaN-based semiconductors to evaluate the IEF of the 4 | 8, hexagonal, and wurtzite configurations, as shown in Fig. 2a. We find that the IEF pointing from the metal toward the N surface is enhanced by decreasing the radius of the metal atoms. The magnitudes are calculated to be 22, 16, and 7 MV cm−1 for ten-layer w-AlN, w-GaN, and w-InN, respectively, as determined by the slope of the internal macroscopic-average electrostatic potential. We also repeat the results of the bilayer w-GaN reported in ref. 39 to confirm the reliability of the method. In contrast, the IEFs of the 4 | 8 and hexagonal configurations are almost zero, demonstrating that the IEF is almost eliminated in 2D structures. This can be explained by the different charge distributions along the c-axis displayed in Fig. 2b. The 4 | 8 and hexagonal configurations are formed by stacked neutral atomic layers, so the electrostatic potentials are equal for each layer, whereas positive (metal atoms) and negative (N atoms) charged atomic layers are alternately stacked in the wurtzite configuration. This causes an electrostatic potential difference, which accumulates with increasing layer number. The results for the ten-layer slabs without hydrogen passivation are displayed in Supplementary Fig. 4. The same trend can be seen and the values are 16, 10, and 6 MV cm−1 for the wurtzite structures.
The ideal materials should not only be IEF free but also maintain optoelectronic properties comparable to those of the wurtzite configurations. The HSE06-level electronic band structures for the different configurations are displayed in Fig. 3 and the band gaps are listed in Table 1. First, the band structures of the bulk wurtzite structures are calculated as references and the band gaps are 5.98, 3.24, and 0.69 eV for w-AlN, w-GaN, and w-InN, respectively. The calculated results agree with the experimental observations. The band structures of the bulk 4 | 8 structures are almost the same as those for the bulk wurtzite structures and the band gaps are 5.57, 3.24, and 0.76 eV, because the contribution of the band edge states has not been affected by the phase transition. The spatial wavefunction distributions of the valence band maximum (VBM) and conduction band minimum (CBM) are displayed in Fig. 4. We find that the CBM is always contributed by the s orbital of N atoms for both the wurtzite and 4 | 8 structures; however, the contributions to the VBM are slightly different. Overall, the VBM is contributed by the p orbital of the N atoms, but the component changes with the wurtzite-to-4 | 8 transition. The VBM is contributed by the pz orbital for w-AlN, by the px orbital for 4 | 8-AlN, by the px orbital for w-GaN, by the pz orbital for 4 | 8-GaN, by the px orbital for w-InN, and by the py orbital for 4 | 8-InN. In the hexagonal structures, the VBM is between the Y and G points (defined as the Q point and marked by blue dots) and contributed by the pz orbital of N atoms, and the CBM is contributed by the hybridized state including N-s and Al/Ga-sp for AlN and GaN; thus, the hexagonal structures show indirect band gaps. h-InN is still a direct band gap material, because the VBM and CBM are contributed by the bonding and antibonding states of the N-s and In-sp orbitals. The band gaps are 4.32, 3.53, and 1.72 eV for h-AlN, h-GaN, and h-InN, respectively. Therefore, by reducing the thickness, the direct → indirect band gap transition occurs for AlN and GaN but not for InN. The Perdew–Burke–Ernzerhof parametrization (PBE) level band structure evolutions for wurtzite and 4 | 8 2D III-nitrides with increasing thickness are exhibited in Supplementary Figs. 5–7.
The absorption coefficient and carrier effective mass for the hexagonal, 4 | 8, and wurtzite structures are calculated and displayed in Fig. 5 and Table 2, which are key parameters for optoelectronic devices. First, the wurtzite (blue line) and 4 | 8 (red line) configurations exhibit stronger band edge absorption intensities than the hexagonal (black line) configuration, as shown in Fig. 5; hence, the wurtzite and 4 | 8 configurations are more suitable for application in lighting devices than the hexagonal configuration. Second, for the hexagonal configuration, the edge absorption intensities are at least one order of magnitude weaker than the second absorption. By analyzing the band structures, as shown in the insets of Fig. 3 (bottom row), we find that the band gaps at the G point (green) are 4.65, 3.76, and 1.71 eV, which match the band edge absorption. On the other hand, the band gaps at the Q point (blue) are 6.37, 6.01, and 4.79 eV, which agree with the second absorption peak in Fig. 5. The VBMs of h-AlN and h-GaN are located at the Q points, so the absorption intensities are stronger than those at the G point. The electron effective mass is hardly changed for the different configurations based on a comparison of the results with those of wurtzite structures. We can find that the in-plane hole effective masses for the 4 | 8 and hexagonal configurations are much smaller than those for the wurtzite configurations. As the hole is the majority carrier in nitride materials, the in-plane mobility of the 4 | 8 and hexagonal structures is much better than that of the wurtzite structures. This will improve the efficiency of horizontal structure optoelectronic devices.
Substrates effect
In experiments, wurtzite and hexagonal structures have already been realized, but the 4 | 8 structure has not yet been observed, even though the 4 | 8 structure is energetically favorable for nanothickness structures. This observation should be attributed to the substrate effect. Here we employ a four-layer GaN slab and two substrates, graphene (Gra) and silicon carbide (SiC), as examples to investigate the substrate effect. First, we cover the wurtzite or 4 | 8 slabs on graphene, acting as vdW epitaxy41,42. After structural relaxation, the 4 | 8 structure remains in the initial configuration, as shown in Fig. 6a, indicating the good stability of the 4 | 8 structures for the four-layer GaN. However, the wurtzite configuration is too unstable to be maintained and changes into the 4 | 8 configuration, as shown in Fig. 6b. However, the situation for the SiC substrate is entirely different. Both the wurtzite and 4 | 8 configurations maintain their initial states after structural relaxation. The interaction between the four-layer GaN slab and SiC resulted in a different behavior. For freestanding four-layer GaN, the 4 | 8 configuration is more stable than the wurtzite configuration. After stacking on graphene, the interaction between the four-layer GaN slab and graphene is weak and cannot restrict the phase transition. On the other hand, the four-layer GaN slab bonded on the SiC substrate has a much stronger interaction than that on graphene, so the interaction is hard to overcome. It is observed from Fig. 6 (middle row) that charge transfer is almost negligible at the graphene/GaN interface but is quite obvious at the SiC/GaN interface. Moreover, the structure of the SiC surface matches that of w-GaN, so w-GaN is more energetically favorable than 4 | 8-GaN on the SiC substrate by 63.9 meV Å−2.
To consider the temperature effect, ab initio molecular dynamics (AIMD) is employed. The temperature is set to 1000 K, which refers to the experimental environment35. The total energy of the system converged after a 4 ps AIMD simulation, as shown in Supplementary Fig. 8, and the structures at the end of the simulations are displayed in the bottom row of Fig. 6. The graphene and materials are separated due to the destruction of the weak interaction. The atoms of 4 | 8-GaN just vibrate around the equilibrium position and the w-GaN transformed to 4 | 8-GaN, indicating that 1000 K is high enough to overcome the wurtzite-to-4 | 8 phase transition barriers (0.30, 0.31, and 0.28 eV f.u.−1 for AlN, GaN, and InN, respectively, as shown in Supplementary Fig. 9). On the SiC substrate, w-GaN is still stable and 4 | 8-GaN on SiC transforms to the wurtzite configuration, demonstrating that w-GaN is more stable than 4 | 8-GaN. Therefore, 4 | 8 configuration materials are expected to be experimentally achieved if we use a weakly bonded substrate such as graphene and control the thickness of the materials to within several nanometers by controlling the growth time. As far as we know, the reported experimental work for the growth of 2D GaN and AlN used SiC or Si as substrates35,36,37. The interface interaction is strong and the surface structures match w-GaN. The growth temperature is usually as high as 1000 K; thus, the 4 | 8 configuration prefers to transform to wurtzite under this condition.
In summary, we have theoretically provided an effective approach, reducing the thickness, to overcome the problem of strong IEF in wurtzite GaN-based semiconductors. The 4 | 8 configuration is predicted to be the most stable configuration when the thickness is reduced to several nanometers. We reveal the wurtzite → 4 | 8 → hexagonal phase transition as the thickness decreases. Interestingly, the 4 | 8 configurations exhibit almost zero IEF, whereas the optoelectronic properties are comparable to those of the wurtzite configurations. The in-plane hole mobility for the hexagonal and 4 | 8 structures is much better than that for the wurtzite structures. The substrate plays an important role in determining the configurations of the materials. The 4 | 8 configuration is expected to be achieved by using the so-called “thickness-controlled” vdW substrates such as graphene and the wurtzite configuration is energetically favorable on SiC or Si substrates. Our results demonstrate that 2D GaN-based semiconductors exhibit many desirable properties and show great potential for application in novel optoelectronic devices.
Methods
General calculation details
First-principles calculations were performed by employing the generalized gradient approximation (PBE) and projector-augmented wave method, as implemented in the Vienna Ab initio Simulation Package (VASP)43. A dispersion correction (DFT-D3) was used to incorporate the long-range vdW interaction44,45. The unit cell for the 4 | 8 structure was rectangular and contained eight atoms. To be comparable, we used the same supercell for the hexagonal and wurtzite structures in our calculations. The lattice stability of the bulk 4 | 8 configurations for group-III nitrides was confirmed by phonon spectrum calculations, as displayed in Supplementary Fig. 1. The simulation cells were built with vacuum slabs thicker than 20 Å to avoid artificial interaction of periodic images. The plane-wave basis set with a kinetic energy of 500 eV was used. We employed 5 × 8 × 1 Monkhorst-Pack k-grids in self-consistent calculations and 51 k-points between high symmetry points of the Brillouin zone in the band structure calculations. In addition, the force and total energy convergence thresholds were set to 0.01 eV Å−1 and 10−6 eV, respectively. The lattice parameter differences were <1% compared with the experimental results of wurtzite structures, as shown in Table 1.
HSE band structures and absorptions
To avoid the underestimation of the band gap from the PBE function, the hybrid exchange and correlation functional of Heyd–Scuseria–Ernzerhof (HSE06)46 was used to obtain an accurate band structure. Here we used a mixing parameter (α) of 33% for AlN and 31% for GaN and InN so that band gaps of 5.98, 3.24, and 0.69 eV, respectively, were obtained for the wurtzite structures. Absorption spectra were calculated based on HSE06 including the random phase approximation. The HSE06 and absorption spectra were obtained by the PWmat code47.
Nudged elastic band and phonon
The climbing image nudged elastic band method48,49 using eight images to connect the initial and final configurations (force convergence of 0.01 eV Å−1) was employed to obtain the structural phase transition barrier. Force constants for phonon calculations were found by using density functional perturbation theory50, as implemented in VASP. On the basis of the calculated force constants, phonon dispersion curves were obtained with the PHONOPY package51.
Ab initio molecular dynamics
Finally, 3 × 3 supercells (more than 300 atoms) were used for AIMD simulations with the canonical ensemble, which were performed to consider the temperature effect on the structural transition at 1000 K (Nosé–Hoover thermostat, total simulation time of 3 ps, time step of 1 fs). AIMD was carried out using the CP2K package52 with the hybrid Gaussian and plane wave method implemented in the Quickstep module53. Three-dimensional periodic boundary conductions were applied for the system under NVT conditions. Temperature control was implemented for the ionic degrees of freedom by using Nosé–Hoover thermostats. The PBE exchange and correlation functional was used in combination with Godecker–Teter–Hutter (GTH)54 pseudopotentials and the DZVP-MOLOPT-SR-GTH basis sets55. The auxiliary plane-wave basis was set to an energy cutoff of 300 Ry and the time step was set to 1 fs.
Surface energy
The surface energy (\(\gamma _{\it{m}}^c\)) for different materials (m = AlN, GaN, or InN) in different atomic configurations (c = w, 4 | 8, or h) per molecular formula unit (f.u.) was employed to estimate the stability of the 2D group-III nitride slabs, which is defined as:
Eslab and \(E_{\mathrm {bulk}}\) represent the total energy of slabs (contain two dangling surfaces) and bulk materials (the number of atoms is the same as in the slab), respectively. N is the number of molecular formula units on the surface. A high \(\gamma _{\it{m}}^c\) implied poor stability.
Carrier effective mass
The carrier effective mass (\(m_{e/h}^ \ast\)) obtained by fitting the band edge was calculated according to the following equation:
where \(\hbar\) is the reduced ∂ Planck constant, ε(k) ε(k) is the band edge eigenvalue, and k is the wavevector. A small \({\it{m}}_{{\it{e/h}}}^ \ast\) implied good carrier mobility. In this work, we studied \({\it{m}}_{{\it{e/h}}}^ \ast\) along Γ (0, 0, 0) → X (0.5, 0, 0), Γ → Y (0, 0.5, 0), and Γ → Z (0, 0, 0.5).
Data availability
All data generated or analyzed during this study are included in this published article (and its Supplementary Information files), and are available from the authors.
Code availability
The calculations were implemented using the VASP, PWmat, CP2K, and Phonopy packages.
References
Mohammad, S. N. & Morkoç, H. Progress and prospects of group-III nitride semiconductors. Prog. Quant. Electron. 20, 361–525 (1996).
Nakamura, S., Mukai, T. & Senoh, M. Candela-class high-brightness InGaN/AlGaN double-heterostructure blue-light-emitting diodes. Appl. Phys. Lett. 64, 1687–1689 (1994).
Nakamura, S. & Krames, M. R. History of gallium–nitride-based light-emitting diodes for illumination. Proc. IEEE 101, 2211–2220 (2013).
Li, D. et al. Realization of a high-performance GaN UV detector by nanoplasmonic enhancement. Adv. Mater. 24, 845–849 (2012).
Sun, Y. et al. Room-temperature continuous-wave electrically pumped InGaN/GaN quantum well blue laser diode directly grown on Si. Light Sci. Appl. 7, 1–7 (2018).
Lim, S.-H. et al. Electrically driven, phosphor-free, white light-emitting diodes using gallium nitride-based double concentric truncated pyramid structures. Light Sci. Appl. 5, e16030 (2016).
Morkoc, H. et al. Large-band-gap SiC, III-V nitride, and II-VI ZnSe-based semiconductor device technologies. J. Appl. Phys. 76, 1363–1398 (1994).
Li, D. et al. Influence of threading dislocations on GaN-based metal-semiconductor-metal ultraviolet photodetectors. Appl. Phys. Lett. 98, 011108 (2011).
Li, D., Jiang, K., Sun, X. & Guo, C. AlGaN photonics: recent advances in materials and ultraviolet devices. Adv. Opt. Photonics 10, 43–110 (2018).
Nakamura, S. et al. InGaN/GaN/AlGaN-based laser diodes with modulation-doped strained-layer superlattices grown on an epitaxially laterally overgrown GaN substrate. Appl. Phys. Lett. 72, 211–213 (1998).
Li, Z. et al. Suppression of thermal degradation of InGaN/GaN quantum wells in green laser diode structures during the epitaxial growth. Appl. Phys. Lett. 103, 152109 (2013).
Zhang, Y. et al. Effects of a step-graded AlxGa1−xN electron blocking layer in InGaN-based laser diodes. J. Appl. Phys. 109, 083115 (2011).
Bernardini, F., Fiorentini, V. & Vanderbilt, D. Spontaneous polarization and piezoelectric constants of III-V nitrides. Phys. Rev. B 56, R10024 (1997).
Fiorentini, V. et al. Effects of macroscopic polarization in III-V nitride multiple quantum wells. Phys. Rev. B 60, 8849 (1999).
Bernardini, F. & Fiorentini, V. Macroscopic polarization and band offsets at nitride heterojunctions. Phys. Rev. B 57, R9427 (1998).
Langer, R. et al. Giant electric fields in unstrained GaN single quantum wells. Appl. Phys. Lett. 74, 3827–3829 (1999).
Morkoc, H. & Mohammad, S. N. High-luminosity blue and blue-green gallium nitride light-emitting diodes. Science 267, 51–55 (1995).
Deguchi, T. et al. Quantum-confined stark effect in an AlGaN/GaN/AlGaN single quantum well structure. Jpn J. Appl. Phys. 38, L914 (1999).
Thamm, A. et al. Optical properties of heavily doped GaN/(Al, Ga)N multiple quantum wells grown on 6H−SiC (0001) by reactive molecular-beam epitaxy. Phys. Rev. B 61, 16025 (2000).
Jiang, H. et al. Photoluminescence and photoluminescence excitation spectra of In0.2Ga0.8N-GaN quantum wells: comparison between experimental and theoretical studies. IEEE J. Quantum Electron. 35, 1483–1490 (1999).
Waltereit, P. et al. Nitride semiconductors free of electrostatic fields for efficient white light-emitting diodes. Nature 406, 865–868 (2000).
Ranalli, F. et al. Non-polar AlN and GaN/AlN on r-plane sapphire. Phys. Stat. Solid. C 6, S780–S783 (2009).
Gil, B. Group III Nitride Semiconductor Compounds: Physics and Applicationss (Clarendon Press, 1998).
Freeman, C. L., Claeyssens, F., Allan, N. L. & Harding, J. H. Graphitic nanofilms as precursors to wurtzite films: theory. Phys. Rev. Lett. 96, 066102 (2006).
Prete, M. S., Pulci, O. & Bechstedt, F. Strong in-and out-of-plane excitons in two-dimensional InN nanosheets. Phys. Rev. B 98, 235431 (2018).
Pashartis, C. & Rubel, O. Alloying strategy for two-dimensional GaN optical emitters. Phys. Rev. B 96, 155209 (2017).
Gürbüz, E., Cahangirov, S., Durgun, E. & Ciraci, S. Single layers and multilayers of GaN and AlN in square-octagon structure: stability, electronic properties, and functionalization. Phys. Rev. B 96, 205427 (2017).
Kolobov, A. et al. Instability and spontaneous reconstruction of few-monolayer thick GaN graphitic structures. Nano Lett. 16, 4849–4856 (2016).
Singh, A. K., Zhuang, H. L. & Hennig, R. G. Ab initio synthesis of single-layer III-V materials. Phys. Rev. B 89, 245431 (2014).
Zhuang, H. L., Singh, A. K. & Hennig, R. G. Computational discovery of single-layer III-V materials. Phys. Rev. B 87, 165415 (2013).
Kecik, D. et al. Fundamentals, progress, and future directions of nitride-based semiconductors and their composites in two-dimensional limit: a first-principles perspective to recent synthesis. Appl. Phys. Rev. 5, 011105 (2018).
Tsipas, P. et al. Evidence for graphite-like hexagonal AlN nanosheets epitaxially grown on single crystal Ag (111). Appl. Phys. Lett. 103, 251605 (2013).
Malin, T., Mansurov, V., Galitsyn, Y. G. & Zhuravlev, K. In 2014 Int. Conference on Manipulation, Manufacturing and Measurement on the Nanoscale (3M-NANO) 308–311 (IEEE, 2014).
Singh, A. K. & Hennig, R. G. Computational synthesis of single-layer GaN on refractory materials. Appl. Phys. Lett. 105, 051604 (2014).
Al Balushi, Z. Y. et al. Two-dimensional gallium nitride realized via graphene encapsulation. Nat. Mater. 15, 1166 (2016).
Wang, W. et al. Lattice structure and bandgap control of 2D GaN grown on graphene/Si heterostructures. Small 15, 1802995 (2019).
Wang, W. et al. 2D AlN layers sandwiched between graphene and Si substrates. Adv. Mater. 31, 1803448 (2019).
Xiao, Y. et al. Anharmonic multi-phonon nonradiative transition: an ab initio calculation approach. Sci. China Phys. Mech. 63, 277312 (2020).
Sanders, N. et al. Electronic and optical properties of two-dimensional GaN from first-principles. Nano Lett. 17, 7345–7349 (2017).
Tasker, P. The stability of ionic crystal surfaces. J. Phys. C Solid State Phys. 12, 4977 (1979).
Kim, J. et al. Principle of direct van der Waals epitaxy of single-crystalline films on epitaxial graphene. Nat. Commun. 5, 4836 (2014).
Shi, Z. et al. Construction of van der Waals substrates for largely mismatched heteroepitaxy systems using first principles. Sci. China Phys. Mech. 62, 127311 (2019).
Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 59, 1758 (1999).
Grimme, S. Semiempirical GGA-type density functional constructed with a long-range dispersion correction. J. Comput. Chem. 27, 1787–1799 (2006).
Barone, V. et al. Role and effective treatment of dispersive forces in materials: Polyethylene and graphite crystals as test cases. J. Comput. Chem. 30, 934–939 (2009).
Heyd, J., Scuseria, G. E. & Ernzerhof, M. Hybrid functionals based on a screened Coulomb potential. J. Chem. Phys. 118, 8207–8215 (2003).
Jia, W. et al. Fast plane wave density functional theory molecular dynamics calculations on multi-GPU machines. J. Comput. Phys. 251, 102–115 (2013).
Sheppard, D., Terrell, R. & Henkelman, G. Optimization methods for finding minimum energy paths. J. Chem. Phys. 128, 134106 (2008).
Sheppard, D. et al. A generalized solid-state nudged elastic band method. J. Chem. Phys. 136, 074103 (2012).
Baroni, S., De Gironcoli, S., Dal Corso, A. & Giannozzi, P. Phonons and related crystal properties from density-functional perturbation theory. Rev. Mod. Phys. 73, 515 (2001).
Chaput, L., Togo, A., Tanaka, I. & Hug, G. Phonon-phonon interactions in transition metals. Phys. Rev. B 84, 094302 (2011).
Hutter, J., Iannuzzi, M., Schiffmann, F. & VandeVondele, J. cp2k: atomistic simulations of condensed matter systems. Wiley Interdiscip. Rev. Comput. Mol. Sci. 4, 15–25 (2014).
VandeVondele, J. et al. Quickstep: fast and accurate density functional calculations using a mixed Gaussian and plane waves approach. Comput. Phys. Commun. 167, 103–128 (2005).
Goedecker, S., Teter, M. & Hutter, J. Separable dual-space Gaussian pseudopotentials. Phys. Rev. B 54, 1703 (1996).
VandeVondele, J. & Hutter, J. Gaussian basis sets for accurate calculations on molecular systems in gas and condensed phases. J. Chem. Phys. 127, 114105 (2007).
Acknowledgements
The research reported in this publication was supported by the National Science Fund for Distinguished Young Scholars (61725403), the National Natural Science Foundation of China (61804152, 61874118, 61834008, and 61922078), the Special Fund for Research on National Major Research Instruments (61827813) Key Research Program of Frontier Sciences, CAS (Grant Number ZDBS-LY-JSC026), the Jilin Provincial Science and Technology Department (20180201026GX), and the Youth Innovation Promotion Association of CAS. It was also supported by the CAS Talents Program.
Author information
Authors and Affiliations
Contributions
D. L. and Z. S. conceived the idea and provided advice on this work. Y. J., Z. H., W. H., Z. Q. and Z. S. carried out the calculations. Y. J. and Z. S. wrote the manuscript. K. J., Y. C., S. Z. and X. S. provided very useful comments. All authors participated in the discussions on the results.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Jia, Y., Shi, Z., Hou, W. et al. Elimination of the internal electrostatic field in two-dimensional GaN-based semiconductors. npj 2D Mater Appl 4, 31 (2020). https://doi.org/10.1038/s41699-020-00165-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41699-020-00165-1