This article is part of the series Low dimensional electric and photonic materials and devices.

Open Access Highly Accessed Nano Express

High Mg effective incorporation in Al-rich Al x Ga1 - x N by periodic repetition of ultimate V/III ratio conditions

Tongchang Zheng, Wei Lin*, Duanjun Cai, Weihuang Yang, Wei Jiang, Hangyang Chen, Jinchai Li, Shuping Li and Junyong Kang*

Author Affiliations

Department of Physics, Fujian Key Laboratory of Semiconductor Materials and Applications, Xiamen University, Xiamen 361005, People’s Republic of China

For all author emails, please log on.

Nanoscale Research Letters 2014, 9:40  doi:10.1186/1556-276X-9-40


The electronic version of this article is the complete one and can be found online at: http://www.nanoscalereslett.com/content/9/1/40


Received:23 December 2013
Accepted:13 January 2014
Published:21 January 2014

© 2014 Zheng et al.; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

According to first-principles calculations, the solubility of Mg as a substitute for Ga or Al in AlxGa1 – xN bulk is limited by large, positive formation enthalpies. In contrast to the bulk case, the formation enthalpies become negative on AlxGa1 – xN surface. In addition, the N-rich growth atmosphere can also be favorable to Mg incorporation on the surface by changing the chemical potentials. On the basis of these special features, we proposed a modified surface engineering technique that applies periodical interruptions under an ultimate V/III ratio condition (extremely N-rich), to enhance Mg effective incorporation. By optimizing the interruption conditions (2 nm interruption interval with 2 s interruption time), the enhancement ratio can be up to about 5 in the Al0.99Ga0.01N epilayer.

Keywords:
First-principles calculations; MOVPE; AlGaN epilayer; High Mg incorporation; Ultimate V/III ratio

Background

AlxGa1 – xN alloys have attracted considerable attention in recent years because of their great potential for applications in UV and deep UV optoelectronic devices with spectral lengths as short as 200 nm [1]. Both high-quality p-type and n-type AlGaN epilayers are strongly demanded for electrical injection in constructing these short wavelength devices. However, similar to most wide bandgap semiconductors, AlGaN suffers from the ‘asymmetric doping’ limitation [2,3], i.e., doping AlGaN to form n-type layer is easy, but achieving p-type doping is difficult [4,5]. Although Mg is the most widely adopted p-type dopant for AlGaN, its doping efficiency is extremely low, particularly for high Al content AlxGa1 – xN [6]. The low doping efficiency of Mg is mainly attributed to its limited solubility, high activation energy, and compensation effect with impurities or native donor defects [2,7]. In spite of the extensive efforts to improve the Mg activation efficiency [5,6,8,9], the bottleneck of low Mg solubility in GaN [10] and AlN [11] materials strongly restricts the overall p-type doping in AlGaN.

Regarding the dopant solubility issue, an extremely high carbon dopant concentration was shown to exist on the epitaxial surface of Si system [12]. This high concentration can be attributed to the surface enhancement effect caused by the partial release of atom mismatch strain. As the epitaxy continues, part of this high concentration dopant segregates to the new surface, and the residual components freezes into the host matrix [12] which corresponds to the final dopant concentration. In other words, the growing surface plays a critical role in determining dopant solubility. If the transient solubility on the growing front surface can be effectively enhanced, high dopant incorporation can be achieved. Theoretical simulations have recently predicted that a N-rich condition is beneficial for Mg incorporation in GaN and AlN [10,11]. However, high V/III ratio was determined to be unfavorable for high-quality AlxGa1 – xN crystal growth [13-16]. Thus, the dilemma between maintaining high V/III ratio to promote Mg incorporation and maintaining low V/III ratio to ensure high crystal quality presents a long-standing challenge for deep UV optoelectronic devices.

In this work, we proposed a method to solve this V/III ratio dilemma by periodically interrupting the AlGaN growth (using usual V/III ratio as the AlGaN growth) and by shortly producing an ultimate V/III ratio condition (extremely N-rich). First-principles simulations were utilized to analyze the behavior of substituting Mg for Al and Ga in the bulk and on the surface of AlxGa1 – xN under different growth atmospheres and to demonstrate the mechanism for the preferred Mg incorporation. On the basis of the analysis results, a modified surface engineering (MSE) technique that utilizes periodical interruptions under an extremely N-rich atmosphere was applied to enhance Mg effective incorporation by metalorganic vapor phase epitaxy (MOVPE). Significant Mg incorporation improvements in Al-rich AlxGa1 – xN epilayer were achieved.

Methods

The first-principles total energy calculations based on density functional theory were performed by using the Vienna ab initio simulation package [17]. Pseudopotentials were specified by the projector augmented wave [18,19] and by generalized gradient approximation [20]. Ga 3d electrons were treated as part of the valence band, and the plane wave cutoff energy was set at 520 eV. Geometry optimizations were performed until the total energy converged to 1 meV. For the bulk calculations, a 2 × 2 × 4 supercell containing 64 atoms [7] and a 5 × 5 × 3 Monkhorst-Pack grid [21] of k-points were used. All atoms were allowed to relax fully for energy minimization. For the surface calculations, we employed a 2 × 2 supercell with six AlxGa1 – xN bilayers separated by a 13-Å wide vacuum region [22] and a 4 × 4 × 1 k-point mesh. The back side of the slab was saturated with hydrogen atoms of fractional charge. The three bottom AlxGa1 – xN bilayers were fixed in the appropriate bulk-optimized configuration to simulate the growth surface, in which all the other layers was relaxed fully.

The Mg-doped AlxGa1 – xN samples were grown on (0001) sapphire substrates via MOVPE. Trimethylgallium (TMGa), trimethylaluminum (TMAl), bis-cyclopentadienylmagnesium (Cp2Mg), and ammonia (NH3) were used as precursors, and H2 was used as carrier gas. Buffer layers with a 20-nm low temperature AlN nucleation layer, a 1-μm high temperature AlN layer, and a graded composition AlGaN layer have been used for initial growth on sapphire. The conventional method for fabricating Mg-doped AlxGa1 – xN was conducted by retaining all the precursors during Mg-doped growth with usual V/III ratio. The MSE technique was implemented by periodically interrupting the conventional growth mode with closing the metal flows (TMAl, TMGa, and Cp2Mg) and continuously maintaining the NH3 flow to shortly produce an ultimate V/III ratio. The Mg and H concentrations were measured by using the Quad PHI 6600 secondary ion mass spectrometer (SIMS) system with depth resolution of approximately 2 nm, and Cs+ ion beams were used as primary ion sources.

Results and discussion

Considering that MOVPE growth is usually characterized by N-rich growth, we first discuss the formation enthalpies of neutral charge state Mg substituting for Al (MgAl) and Ga (MgGa) in AlxGa1 – xN bulk as a function of Al content under N-rich condition. The calculated results are shown in Figure 1a, wherein both the MgAl and MgGa formation enthalpies are positive and large, thus indicating limited Mg solubility. The formation enthalpies of MgAl in AlN and MgGa in GaN are comparable with previous results [10,11]. As the Al content in AlxGa1 – xN increases, both the MgAl and MgGa formation enthalpies monotonically increase. The formation enthalpy ΔHf is closely related to the equilibrium Mg solubility C, which is given by [10]:

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/40/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/40/mathml/M1">View MathML</a>

(1)

where Nsites is the number of sites on which the dopant can be incorporated, kB is the Boltzmann constant, and T denotes the temperature. Large formation enthalpy yields low dopant solubility. At the growth temperature (T = 1,000°C), the Mg solubility in bulk GaN is approximately 1.65 × 1017 cm-3. Considering that ΔHf increases with increasing Al content, AlxGa1 – xN experiences an aggravating Mg solubility limit. The Mg solubility limit may even decrease to approximately 2.32 × 1016 cm-3 in AlN (for T = 1,200°C). On the basis of this tendency, incorporating Mg becomes more difficult in Al-rich AlxGa1 – xN. Notably, the formation enthalpy for MgAl is larger than that for MgGa over the entire Al content range. This characteristic demonstrates that substituting Mg for Al is more energetically unfavorable than substituting Mg for Ga, which also explains the low Mg incorporation in Al-rich AlxGa1 - xN. Such behavior of Mg is partly attributable to its larger covalent radius (1.36 Å) compared with those of Al (1.18 Å) and Ga (1.26 Å), as well as the compressive strain after Mg substitution [23,24]. As shown in the inset of Figure 1a, the AlxGa1 – xN lattice constants a and c decrease as the Al content increases, thus making the mismatch strain caused by substituting Mg for Al or Ga atoms with smaller radii becomes more considerable.

thumbnailFigure 1. Formation enthalpies of MgGa/MgAl and normalized CMg cprofile of AlGaN films. (a) In the bulk and (b) on the surface of AlxGa1 -xN as a function of Al content under N-rich condition. (c) Normalized CMg of AlxGa1 -xN (x = 0.33, 0.54) epilayers from the surface to bulk. The inset in (a) shows the calculated AlxGa1 -xN lattice constants a and c as a function of Al content. The inset in (c) illustrates the source supply sequence of the conventional method.

Reducing the formation enthalpy is believed to be the key issue in solving the problem of Mg incorporation. The formation enthalpy is governed by two important factors, as given by [11]:

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/40/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/40/mathml/M2">View MathML</a>

(2)

Here, ΔE = EMg- Ehost, where EMg and Ehost are the total energies of the supercell with and without Mg substitution; Δμ = μAl/Ga– μMg, where μi (i = Al, Ga, Mg) is the chemical potential. ΔE is mainly induced by the strain caused by the atom size mismatch. Consequently, larger atom size mismatch results in large ΔE, thus resulting in larger ΔHf as mentioned above. The strain induced by the C-dopant into the Si matrix becomes smaller on the surface than in the bulk [12]. The question of whether the surface also plays a similar role in the Mg incorporation in AlxGa1 – xN arises. To address this issue, we further investigated the formation enthalpies of MgAl and MgGa on AlxGa1 – xN surface, and the results are shown in Figure 1b. In contrast to the bulk case, both of the formation enthalpies in the surface are negative, suggesting favorable Mg substitution. The value of MgAl becomes lower than that of MgGa and decreases as the Al content increases. These interesting reversed tendencies provide us a possible way to promote the Mg incorporation in Al-rich AlxGa1 – xN by utilizing the surface effect. An epitaxy growth, e.g., MOVPE and molecular beam epitaxy systems, is conducted under an inherently non-equilibrium process with the surface transforming into a bulk [12]. Therefore, enhancing the Mg incorporation by using the surface effect should be practically feasible.

Two Mg-doped AlxGa1 – xN (x = 0.33, 0.54) films were grown by MOVPE using the conventional method (the inset of Figure 1c) to validate the prediction of the surface effect on Mg incorporation. As shown in Figure 1c, the Mg concentration (CMg) on the surface is about one order higher than that of in the bulk for both samples. Although CMg rapidly falls beneath the topmost surface (about 10 nm), CMg is still several orders higher than the theoretical prediction by Equation 1. This phenomenon can be understood in terms of the competition between the Mg incorporation enhancement on the growing surface due to the surface effect and the Mg segregation as the epitaxy continues. Simply, when the surface with the enhanced Mg solubility is covered by newly added layers during further growth, most of these Mg segregates to the new surface to regain equilibrium because the surface transforms into a bulk. Meanwhile, considerable part of these Mg is frozen in because of solidification. As a result, the CMg in the bulk is lower than that of in the final epilayer surface but is much higher than the equilibrium value of the ideal bulk. Considering this competition, Mg incorporation can be modified by balancing the surface time and solidification time.

As shown in Equation 2, the factor Δμ also affects Mg incorporation. In principle, the chemical potential μ sensitively depends on the growth atmosphere (N-rich or Al/Ga-rich). The μ of a given species under equilibrium conditions is equal in all phases that are in contact [22]. Therefore, we can obtain

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/40/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/40/mathml/M3">View MathML</a>

(3)

In addition, CMg is limited by the formation of Mg3N2 to substitute Mg for Ga or Al as an acceptor [10]. This limitation meets the relation

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/40/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/40/mathml/M4">View MathML</a>

(4)

By substituting Equations 3 and 4 into Equation 2, we can obtain

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/40/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/40/mathml/M5">View MathML</a>

(5)

which, aside from ΔE, depends only on μN, since the μAlN/GaN and <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/40/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/40/mathml/M6">View MathML</a>are constants [25]. μN should be limited between μN (Al/Ga-rich) ≤ μN≤ μN (N-rich) [11], namely, <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/40/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/40/mathml/M7">View MathML</a>, to drive the source materials to form AlxGa1 – xN alloys instead of the undesirable phases (bulk Ga, Al, and N2). Our calculated ΔHGaN value of -1.01 eV is higher than the ΔHAlN value of -2.97 eV, which are consistent with the experimental values of -1.08 and -3.13 eV [25]. Therefore, as the growth condition varies from Ga-rich to N-rich conditions, μN changes from <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/40/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/40/mathml/M8">View MathML</a> to <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/40/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/40/mathml/M9">View MathML</a>. Thus, ΔHf varies over a range corresponding to 1/3ΔHGaN of 0.337 eV, as shown in Figure 2a. This variation indicates that the N-rich growth atmosphere favor the Mg incorporation effectively in AlGaN. Generally, the N-rich condition is modulated by increasing the V/III ratio. However, for a fixed III flow, the AlxGa1 – xN growth has an optimal V/III ratio for the best crystal quality [13-16]. Nonetheless, the max flow limitation of the MOVPE system does not allow the V flow to be increased infinitely. Accounting for these limitations, an inspiration can be obtained from Figure 1c, in which the protecting atmosphere with NH3 flow just provides an ultimate V/III ratio condition (extremely N-rich) for CMg enhancement when the epitaxy ends with the III flow becoming zero. Simultaneously, the stopped growth avoids the formation of low-quality AlxGa1 – xN crystal. If this special condition is introduced as an intentional interruption during the continuous p-AlxGa1 – xN growth, then the overall Mg incorporation could be improved.

thumbnailFigure 2. Formation enthalpy difference of MgGa/MgAl and CMg profile of Al0.49Ga0.51N film. (a) Formation enthalpy difference of MgGa and MgAl between Ga-rich and N-rich condition. (b)CMg profile of Al0.49Ga0.51N film with three different Cp2Mg flows grown by the MSE technique. The inset in (b) illustrates the source supply sequence of the MSE technique, an ultimate V/III ratio condition is shortly produced during the interruption.

To validate this hypothesis, a growth interruption experiment was designed, as shown schematically in the inset of Figure 2b. We closed the metal flows (TMAl, TMGa, and Cp2Mg flows) three times. In these three periods (35 nm thick), different Cp2Mg flows (0.45, 0.81, and 0.99 nmol/min) were applied to investigate the interruption effect systematically. Figure 2b shows the SIMS CMg profile of Al0.49Ga0.51N film across three periods. Obviously, CMg exhibits three distinct peaks at the interruptions, and the profile of each interruption is similar to that of the final epilayer surface (Figure 1c). Since there are always some Mg floating on the surface during growth because of segregation [26], the interruption will drive the floating Mg to incorporate into the AlxGa1 – xN crystal, thus greatly enhancing Mg solubility. This result confirms that the Mg incorporation on the growing surface can be transiently enhanced further by an extremely N-rich condition interruption, thereby increasing the CMg that would reside at the interrupting region.

However, the CMg enhancement at the interruption region is much smaller than that on the final epilayer surface (Figure 1c), and the CMg far from the interruption region remains low. This result is caused by the wide interval between consecutive interruptions, considerably decreasing the CMg at the interruption regions and resulting in the non-uniformity of the CMg distribution by Mg segregation and diffusion after interruption (Figure 3a). Therefore, the interruption interval, interruption time, and growth rate should play critical roles in affecting the CMg overlap. As illustrated in Figure 3b, we further proposed the MSE technique, optimizing the interruption conditions, to incorporate surface Mg atoms before they can re-segregate to the surface, thus further increasing the average Mg incorporation and approaching a uniform Mg distribution over the entire AlGaN epilayer instead of being distributed locally.

thumbnailFigure 3. Schematic diagram of the Mg incorporation behavior in the AlGaN grown by the MSE technique. As the interruption interval is long, only some peaks distribute locally at the interruptions after Mg segregation and diffusion (a), optimizing the interruption interval, a high and uniform Mg distribution over the entire AlGaN epilayer could be achieved (b).

Three Mg-doped AlxGa1 – xN (x = 0.54, 0.76, 0.99) samples were grown by using the MSE technique (the inset of Figure 2b). An optimized 2-nm interruption interval combining with 2-s interruption time were used for all samples, with Cp2Mg flux of 0.81 nmol/min. As shown in Figure 4a, the samples with different Al contents exhibit high CMg range from 4 × 1019 cm-3 to 5 × 1019 cm-3 and homogeneous distribution at a wide region as expected, whereas the CMg of the samples grown via conventional method decrease with increasing Al content, which is consistent with the theoretical prediction. By comparison, the average CMg in the samples with different Al contents increase several times, and the enhancement ratios increase as the Al content increases, as shown in Figure 4b. Particularly, the enhancement ratio is approximately up to 5 in the Al0.99Ga0.01N. These results indicate that a high CMg can be easily achieved in Al-rich AlGaN by combining the surface effect with the N-rich growth atmosphere modulation.

thumbnailFigure 4. Bulk CMg of the samples and enhancement ratios of Mg/H concentrations. (a) Bulk CMg of the samples with different Al contents grown by the MSE technique and the conventional method. (b) Enhancement ratios of Mg and H concentrations by the MSE technique as a function of Al content compared with that of the conventional method.

High Mg doping was reported to result in Mg-rich precipitates. The primary Mg-rich precipitates were presumed to be Mg3N2[27,28], which can be formed when Mg do not incorporate as acceptors in the desired substitutional sites. The substitutional Mg was suggested to be usually passivated by H during growth, and the corresponding Mg acceptor can be activated by postgrowth thermal annealing to dissociate the Mg - H complex [29]. The correlation between the substitutional Mg and H was verified by previous theoretical and experimental investigations [30,31]. Thus, the H concentration is most likely associated with CMg if Mg is effectively incorporated in the desired substitutional sites. The enhancement ratios of H concentration for the MSE technique increase from 1.2 to 10 with increasing Al content, compared with that of the conventional method, as shown in Figure 4b. This simultaneous enhancement in H concentration demonstrates that the Mg was effectively incorporated in the desired substitutional sites by the MSE technique. In this work, the high CMg is the important basis for improving the hole concentration in p-type AlGaN epilayer. Besides the solubility limit, the high activation energy of Mg acceptors is another contribution for the low p-type doping of AlxGa1 – xN, leading to a low acceptor activation probability [5,8]. In order to increase the overall p-type doping, more efforts on activating the obtained high CMg will be included in future progress.

Conclusions

The MSE technique, which utilizes periodical interruptions under an extremely N-rich atmosphere, was proposed to enhance Mg incorporation, base on the first-principles total energy calculations. During the interruption, metal flows were closed to produce an ultimate V/III ratio condition without affecting the AlGaN growth. By optimizing the interruption conditions, we obtained a high concentration and uniform distribution Mg in the AlGaN epilayer. The CMg enhancements increase with increasing Al content through this method. Particularly, for the Al0.99Ga0.01N, the enhancement ratio can be achieved up to about 5 and the final Mg concentration was determined to be 5 × 1019 cm–3. Meanwhile, the simultaneous increase of the H concentration confirms the Mg effective incorporation in the desired substitutional sites instead of forming Mg3N2. The proposed approach, which is convenient as well as effective, could be used as a general strategy to promote dopant incorporation in wide bandgap semiconductors with stringent dopant solubility limits.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

TCZ carried out the experiments and drafted the manuscript. WHY, WJ and HYC helped in the preparation and characterization of the samples. JCL and SPL took part in the data analysis. WL, DJC and JYK conceived the study and participated in the data analysis and the critical review of the manuscript. All authors read and approved the final manuscript.

Acknowledgments

This work was supported by the ‘973’ program (2012CB619301 and 2011CB925600), the ‘863’ program (2011AA03A111 and 2014AA032608), the National Natural Science Foundation of China (61106008, 60827004, and 90921002), the Natural Science Foundations of Fujian Province (2010 J01343 and 2012 J01024), and the fundamental research funds for the central universities (2011121042).

References

  1. Nakarmi ML, Nepal N, Lin JY, Jiang HX: Photoluminescence studies of impurity transitions in Mg-doped AlGaN alloys.

    Appl Phys Lett 2009, 94:9. OpenURL

  2. Yan Y, Li J, Wei SH, Al-Jassim MMA: Possible approach to overcome the doping asymmetry in wideband gap semiconductors.

    Phys Rev Lett 2007, 98(13):135506. PubMed Abstract | Publisher Full Text OpenURL

  3. Yan Y, Zhang SB, Pantelides ST: Control of doping by impurity chemical potentials: predictions for p-type ZnO.

    Phys Rev Lett 2001, 86(25):5723-5726. PubMed Abstract | Publisher Full Text OpenURL

  4. Nam KB, Nakarmi ML, Li J, Lin JY, Jiang HX: Mg acceptor level in AlN probed by deep ultraviolet photoluminescence.

    Appl Phys Lett 2003, 83(5):878-880. Publisher Full Text OpenURL

  5. Li JC, Yang W, Li S, Chen H, Liu D, Kang J: Enhancement of p-type conductivity by modifying the internal electric field in Mg- and Si-delta-codoped AlxGa1-xN/AlyGa1-yN superlattices.

    Appl Phys Lett 2009, 95:15. OpenURL

  6. Szabo A, Son NT, Janzen E, Gail A: Group-II acceptors in wurtzite AlN: a screened hybrid density functional study.

    Appl Phys Lett 2010, 96:19. OpenURL

  7. Wei S–H, Zhang SB: Chemical trends of defect formation and doping limit in II-VI semiconductors: the case of CdTe.

    Phys Rev B 2002, 66(15):155211. OpenURL

  8. Simon J, Protasenko V, Lian C, Xing H, Jena D: Polarization-induced hole doping in wide-band-gap uniaxial semiconductor heterostructures.

    Science 2010, 327(5961):60-64. PubMed Abstract | Publisher Full Text OpenURL

  9. Schubert EF, Grieshaber W, Goepfert ID: Enhancement of deep acceptor activation in semiconductors by superlattice doping.

    Appl Phys Lett 1996, 69(24):3737-3739. Publisher Full Text OpenURL

  10. Neugebauer J, VandeWalle CG: Role of hydrogen in doping of GaN.

    Appl Phys Lett 1996, 68(13):1829-1831. Publisher Full Text OpenURL

  11. Stampfl C, Van de Walle CG: Theoretical investigation of native defects, impurities, and complexes in aluminum nitride.

    Phys Rev B 2002, 65(15):155212. OpenURL

  12. Tersoff J: Enhanced solubility of impurities and enhanced diffusion near crystal surfaces.

    Phys Rev Lett 1995, 74(25):5080-5083. PubMed Abstract | Publisher Full Text OpenURL

  13. Keller S, Parish G, Fini PT, Heikman S, Chen CH, Zhang N, DenBaars SP, Mishra UK, Wu YF: Metalorganic chemical vapor deposition of high mobility AlGaN/GaN heterostructures.

    J Appl Phys 1999, 86(10):5850-5857. Publisher Full Text OpenURL

  14. Allerman AA, Crawford MH, Fischer AJ, Bogart KHA, Lee SR, Follstaedt DM, Provencio PP, Koleske DD: Growth and design of deep-UV (240–290 nm) light emitting diodes using AlGaN alloys.

    J Cryst Growth 2004, 272(1–4):227-241. OpenURL

  15. Imura M, Fujimoto N, Okada N, Balakrishnan K, Iwaya M, Kamiyama S, Amano H, Akasaki I, Noro T, Takagi T, Bandoh A: Annihilation mechanism of threading dislocations in AlN grown by growth form modification, method using V/III ratio.

    J Cryst Growth 2007, 300(1):136-140. Publisher Full Text OpenURL

  16. Banal RG, Funato M, Kawakami Y: Growth characteristics of AlN on sapphire substrates by modified migration-enhanced epitaxy.

    J Cryst Growth 2009, 311(10):2834-2836. Publisher Full Text OpenURL

  17. Kresse G, Furthmüller J: Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set.

    Phys Rev B 1996, 54(16):11169-11186. Publisher Full Text OpenURL

  18. Blöchl PE: Projector augmented-wave method.

    Phys Rev B 1994, 50(24):17953-17979. Publisher Full Text OpenURL

  19. Kresse G, Joubert D: From ultrasoft pseudopotentials to the projector augmented-wave method.

    Phys Rev B 1999, 59(3):1758-1775. OpenURL

  20. Perdew JP, Burke K, Ernzerhof M: Generalized gradient approximation made simple.

    Phys Rev Lett 1996, 77(18):3865-3868. PubMed Abstract | Publisher Full Text OpenURL

  21. Monkhorst HJ, Pack JD: Special points for Brillouin-zone integrations.

    Phys Rev B 1976, 13(12):5188-5192. Publisher Full Text OpenURL

  22. Timon V, Brand S, Clark SJ, gibson MC, Abram RA: First-principles calculations of 2 × 2 reconstructions of GaN(0001) surfaces involving N, Al, Ga, In, and As atoms.

    Phys Rev B 2005, 72(3):035327. OpenURL

  23. Sadigh B, Lenosky TJ, Caturla MJ, Quong AA, Benedict LX, de la Rubia TZ, Giles MM, Foad M, Spataru CD, Louie SG: Large enhancement of boron solubility in silicon due to biaxial stress.

    Appl Phys Lett 2002, 80(25):4738-4740. Publisher Full Text OpenURL

  24. Zhu J, Liu F, Stringfellow GB, Wei SH: Strain-enhanced doping in semiconductors: effects of dopant size and charge state.

    Phys Rev Lett 2010, 105(19):195503. PubMed Abstract | Publisher Full Text OpenURL

  25. Zoroddu A, Bernardini F, Ruggerone P: First-principles prediction of structure, energetics, formation enthalpy, elastic constants, polarization, and piezoelectric constants of AlN, GaN, and InN: comparison of local and gradient-corrected density-functional theory.

    Phys Rev B 2001, 64(4):045208. OpenURL

  26. Bungaro C, Rapcewicz K, Bernholc J: Surface sensitivity of impurity incorporation: Mg at GaN (0001) surfaces.

    Phys Rev B 1999, 59(15):9771-9774. Publisher Full Text OpenURL

  27. Hansen M, Chen LF, Lim SH, DenBaars SP, Speck JS: Mg-rich precipitates in the p-type doping of InGaN-based laser diodes.

    Appl Phys Lett 2002, 80(14):2469-2471. Publisher Full Text OpenURL

  28. Vennéguès P, Leroux M, Dalmasso S, Benaiisa M, De Mierry P, Lorenzini P, Damilano B, Beaumont B, Massies J, Gibart P: Atomic structure of pyramidal defects in Mg-doped GaN.

    Phys Rev B 2003, 68(23):235214. OpenURL

  29. Nakamura S, Iwasa N, Senoh M, Mukai T: Hole compensation mechanism of p-type GaN films.

    Japanese Journal of Applied Physics Part 1-Regular Papers Short Notes & Review Papers 1992, 31(5A):1258-1266. OpenURL

  30. Clerjaud B, Côte D, Lebkiri A, Naud C: Infrared spectroscopy of Mg-H local vibrational mode in GaN with polarized light.

    Phys Rev B 2000, 61(12):8238-8241. Publisher Full Text OpenURL

  31. Limpijumnong S, Northrup JE, Van de Walle CG: Entropy-driven stabilization of a novel configuration for acceptor-hydrogen complexes in GaN.

    Phys Rev Lett 2001, 87(20):205505. PubMed Abstract | Publisher Full Text OpenURL