Introduction
Advances in nitride materials have led to the development of state-of-the-art LEDs in the visible range and lasers in the near uv [Reference Nakamura and Fasol1]. While the success of these nitride devices has been remarkable, there is a need for an improved understanding of the recombination process. Interest at the moment is focused on the influence of the spontaneous polarization, the piezoelectric effect arising mismatch-induced strain and the potential fluctuations associated with the tendency in InGaN alloys toward phase separation. Numerous studies of these effects have been published. Among the optical characterizations are studies of photoluminescence, photoluminescence excitation and photoluminescence decay [Reference Nakamura and Fasol1-Reference Schmidt, Yong-Hoon, Gainer, Song, Keller, Mishra and DenBaars5]. These studies reveal the effects of charge separation of the electron and hole in the growth direction and localization of the carriers or excitons by potential fluctuations in the growth plane.
Optically detected magnetic resonance (ODMR) has recently been observed in green and extrablue InGaN single-quantum well diode structures [Reference Glaser, Kennedy, Carlos, Ruden and Nakamura6]. The distinct electron and hole resonances denote significant charge separation in the lightly photo-excited structures. Analysis of the frequency dependence of the ODMR for magnetic-dipole microwave transitions gives long lifetimes [Reference Carlos, Glaser, Kennedy and Nakamura7]. The ODMR technique is particularly valuable since the g-factors directly reflect the symmetry of the recombining electron and hole. This symmetry reveals the detailed structure determining the electron- and hole-states.
This paper explores the link between the symmetry as revealed by ODMR and the vertical and lateral structure of the InGaN LED. Theories for the g-factor of electrons and localized holes are presented, extended and compared with experiment. For the electron, the present extensions of the theory are still inadequate to achieving a good understanding of the data. For the hole, new spectral data is presented and discussed using a recent theory for acceptors in GaN [Reference Malyshev, Merkulov and Rodina8]. The analysis shows that the hole states are affected by strong in-plane strain.
Symmetry of the Electron
The g-factor for a conduction electron in a direct-gap semiconductor can be calculated using k • p theory [Reference Hermann, Weisbuch, Meier and Zakharchenya9]. The g-values for GaN were calculated within a cubic approximation [Reference Carlos, Freitas, Khan, Olson and Kuznia10]. It is necessary to include both the valence bands and the upper conduction bands in order to achieve agreement with experiment. The experiments in InGaN LEDs show changes in the average g-factor and in the anisotropy of the g-factor[Reference Glaser, Kennedy, Carlos, Ruden and Nakamura6]. Hence it is necessary to extend the calculations to include the wurtzitic crystal fields, the changes in band structure in moving from GaN to the InGaN alloy, and effects due to the quantum well and differing barrier layers. The following paragraphs address these issues.
First, the anisotropic g-tensor for electrons in wurtzitic GaN is calculated. Two sets of band parameters were chosen [Reference Kim, Lambrecht, Segall and Schilfgaarde11]. The first is a “best” set from a mix of theory and experiment within a quasi-cubic approximation. The second is a full set of calculated values. The energy differences , crystal-field and spin-orbit splittings and interband matrix elements are shown in Fig. 1 and collected in Table 1. For coupling between the lowest conduction band and the uppermost valence band, the expressions for the g-factors are
where
E0 | Δ1 | Δ2 | Δ3 | P2 | E0’ | Δ1’ | Δ2’ | Δ3’ | P’2 | |
---|---|---|---|---|---|---|---|---|---|---|
Best | 3.5eV | 10 | 6.3 | (6.3 | 22.5 | 9.015 | −0.406 | 0.112 | (0.112 | 9 eV |
GaN | meV | meV | meV) | eV | eV | eV | eV | eV) | ||
Calcul | 36 | 5.0 | 5.9 | |||||||
ated | meV | meV | meV | |||||||
GaN | ||||||||||
InN | 1.89 | 43 | 0.3 | 0.3 | 22.5 | 7.91 | 0.284 | 0.284 | 9 eV | |
eV | meV | meV | meV | eV | eV | meV | meV |
Similarly, the contributions from the upper conduction bands were calculated using the same equations with the upper conduction band parameters replacing the valence band contributions. The contributions from both bands add and are of the same magnitude. The resulting total g-factors are given in Table 2 for both the “best” parameters and the calculated parameters. The average shift from the free-electron g-factor is predicted well. The anisotropy is of the right order but with differing signs: disagreeing with experiment for the “best” parameters but with the right sense for the calculated parameters. Thus the anisotropy in wurtzitic GaN depends in a delicate way on the band parameters.
g∥ | g⊥ | |
---|---|---|
Experiment | 1.9510 | 1.9483 |
Best | 1.952 | 1.953 |
Calculated | 1.955 | 1.954 |
In previous work the g-values have been calculated for In-mole fractions up to 0.4 [Reference Glaser, Kennedy, Carlos, Ruden and Nakamura6,Reference Carlos, Glaser, Kennedy and Nakamura7]. The calculations were done for unstrained alloys assumed to be cubic. Here these calculations are improved by incorporating the recent results for the energy-gap dependence on alloy composition [Reference Wetzel, Takeuchi, Yamaguchi, Katoh, Amano and Akasaki12]. The new relationship allows a comparison of the g-factors with calculations for both pseudomorphic and fully relaxed InGaN.
The new calculations were done by combining the gap equation with a k• p calculation within a quasi-cubic approximation. The equation for the energy gap is
where the bowing parameter, b, is 2.6 eV for a pseudomorphic layer and 3.8 eV for a fully relaxed layer. Since the original estimates for x for the green and extra-blue diodes were made with a bowing parameter of 1 eV, these need to be revised. The green changes from 0.43 to 0.22 and the extra-blue from 0.3 to 0.16. The band parameters for InN are shown in Table 1. The interband matrix elements are taken from GaN. E0’ and Δ2’ are scaled from GaN using the ratio of the values of atomic transitions. Δ1 and Δ2 for InN are taken from first-principles calculations [Reference Wei and Zunger13]. For the pseudomorphic layers, strain introduces a Δ1 by a deformation potential taken also from the first-principles calculations. Δ1, Δ2, Δ2’ and E0’ are then scaled to different In mole fractions by a virtual crystal approach. The g-shifts versus mole-fraction obtained in this way needed to be scaled by 15 % to match at GaN.
The calculated g-factors versus mole fraction for both relaxed and pseudomorphic layers are plotted along with the data in Fig. 2. The experimental values for the green diode, expected to be relaxed, fall close to the calculated values. However, the g-values for the extra-blue diode lie substantially above the calculated lines. From the reversal in anisotropy for this structure and its lower mole-fraction, the extra-blue diode is expected to be pseudomorphic. Thus while the calculation approximates the g-factors for the green diode, other interactions must be important for the extra-blue diode.
The next step in this analysis is to consider the effect of the quantum well potential. Theories have been developed for other heterostructure systems [Reference Ivchenko and Kisilev14]. The strain, confinement and piezoelectric effect present in the InGaN LED structure certainly provide enough interactions to produce the reversal in symmetry observed experimentally. Hence, it does seem likely that a clear understanding of the symmetry of the electronic state can be deduced from the experimental data.
Symmetry of the Hole
The g-factors for an acceptor in a wide-bandgap semiconductor have been calculated [Reference Malyshev, Merkulov and Rodina8] for an acceptor including both Coulomb and central-cell potentials. Formulation of the problem in k-space avoids a divergence at the origin. The calculation is first done in a spherical approximation with only the Coulomb and central-cell potentials. The Coulomb interaction gives a binding energy of 200 meV. The crystal field is introduced as a strain in the [111] direction which corresponds to a crystal-field splitting of 37 meV. Spin-orbit interaction is then added to produce a splitting of 12 meV. This reproduces the conditions in GaN with perfect crystal structure and leads to highly anisotropic g-values. Since previous ODMR studies of acceptors in GaN revealed g-tensors which are only weakly anisotropic [Reference Kunzer, Baur, Kaufmann, Schneider, Amano and Akasaki15,Reference Glaser, Kennedy, Doverspike, Rowland, Gaskill, Freitas, Asif Khan, Olson, Kuznia and Wickenden16], the authors then consider an additional strain which lowers the symmetry at the acceptor from axial to rhombic and produces a reduction in the g-anisotropy. The authors attribute the reduction in symmetry to the Jahn-Teller effect but the symmetry lowering could arise from other sources.
The results of this theory can be summarized by considering the g-factors as a function of the symmetry-lowering strain energy (See Fig. 3). For an acceptor in purely axial symmetry g∥ is −3.3 and g⊥ is 0. In-plane strain becomes important as its associated energy becomes comparable to the spin-orbit splitting of 12 meV. For large in-plane strain, the g-tensor becomes less anisotropic and reverts to the free-electron value of 2.00.
When all the luminescence is collected from the green LED, the g-tensor for the hole is g∥ =2.08 (1) and g⊥ = 1.996 (3)[Reference Glaser, Kennedy, Carlos, Ruden and Nakamura6]. The anisotropy is strongly reduced from that expected for holes in purely axial symmetry. Using filters, a spectral study has been performed which reveals that the g∥ changes for emission at different energies (See Fig. 4). In the study, long-pass filters were used to detect all the emission across the green band and the emission from the lower-energy half of the band. A subtraction reveals the part characteristic of the upper energy-range with g∥ = 2.20 and the one filter gives the lower range with g∥ = 2.09. There is a direct correlation between decreasing energy and decreasing g-value. This may arise in part from a changing hole g-factor with changing In-concentration. However, the character of the hole may also be changing. Since a deepening in energy is expected with increasing in-plane strain, the reduction in anisotropy with energy may also arise from increased in-plane strain.
In summary, the g-tensor for holes in lightly photo-excited InGaN reveals states which differ from those of the band edge and those for bound holes in purely axial symmetry. The binding could arise from either an acceptor or a fluctuation in composition of the alloy. The theory which describes the changes in g-tensor for different values of in-plane strain provides a solid framework to interpret the data. Spectral studies of the g-factor show that the anisotropy decreases with decreasing emission-energy. Compositional change and greater non-axial strain produce the reduction in anisotropy. These results are consistent with the behavior of acceptors in other semiconductors [Reference Morgan17].
Summary
The symmetry of the g-factors of recombining electrons and holes in InGaN is related to the band structure and local perturbations of that structure. For GaN, the average g-factor for electrons can be deduced from known band parameters but the anisotropy is difficult to calculate. Improved calculations are made for both pseudomorphic and relaxed InGaN. Again the anisotropy is beyond the current calculations. The hole g-anisotropy is strongly reduced from that for the band edge. A recent theory for acceptors in GaN demonstrates the experimental features and reconciles them with strong in-plane strain.
Acknowledgments
The work at NRL was supported by ONR. We thank I. Merkulov for helpful conversations. PPR thanks W.R.L. Lambrecht for unpublished GaN band structure data.