Introduction
Wide-bandgap semiconductors such as GaN, AlN, InN and their ternary compounds have been studied so as to realize the semiconductor blue-green laser. As a matter of fact, the demonstration of high-brightness blue-green light-emitting diodes from the wurtzite GaN semiconductor has stimulated a lot of research on the fabrication of blue-green laser diodes using nitride-based semiconductors [Reference Nakamura1-Reference Mack, Abare, Hansen, Kozodoy, Keller, Mishra, Coldren and Denbaars3]. InGaN/GaN single-quantum-well structures have been used to achieve high-power blue and green light-emitting diodes. Also, room temperature pulsed lasing of a strained InGaN/GaN multiquantum-well laser diode has recently been demonstrated.
The nature crystal structure of most of the group III nitrides is wurtzite. Strain is also present in most of the group III nitrides because of the lack of substrate materials with a matched lattice constant and a thermal expansion coefficient. Indeed, fundamental studies of the effects of strain on wurtzite band structure play an important role in understanding the electronic and optical properties of GaN-based optoelectronic materials and devices.
In order to tune the interesting wavelength, the concept of band-gap engineering is a useful tool to accept the particular devices operation [Reference Li4]. Basically, intermixing process is one of easy way to achieve the modification of bandstructure. The intermixing process involves the interdiffusion of the constituent atoms of the heterostructure, the process temperature, and time. During the process the as-grown square-QW compositional profile is modified to a graded profile thereby altering the confinement profile and subband structure in the QW. The interdiffusion of InGaN/GaN QW system is more complex compared with AlGaAs/GaAs and InGaAs/GaAs because of the immiscibility of InN and GaN [Reference McCluskey, Romano, Krusor, Johnson and Suski5].
Theoretical Model
This paper presents a theoretical study of the effects of interdiffusion on a 25Å well/70Å barrier of an In0.2Ga0.8N/GaN MQW structure, whose properties leads to the shift in spontaneous emission rate. The effects of interdiffusion on the strains, the splitting of the heavy-hole (HH), the light-hole (LH), the crystal-field split-hole (CH) subbands, the conduction band (CB) and their dispersions, are calculated theoretically, based on the Hamiltonian derived using the multiband k•p model by Chuang and Chang [Reference Chuang and Chang6].
The QW composition intermixing refer to the thermal induced interdiffusion of the constituent atoms through the hetero-interface. In fact, the interdiffused QW (DFQW) structures created by both impurity induced and impurity-free vacancy promoted processes have recently attracted much attention [Reference Yuan, Kim, Jagadish, Burke, Gal, Cai, Zou, Cai, Cockayne and Cohen7-Reference Krauss, Marsh and De La Rue9]. The interdiffusion mechanism can be a single-phase diffusion for two constituent atoms or it can be a two-phase or multiple-phase one and/or for multiple species. In InGaN/GaN DFQW structures, only the interdiffusion of group-III atoms, i.e. In and Ga atoms, occurs as there is no N concentration gradient across the interface. The diffusion of group III atoms in the QW structure is usually described by the Fick’s law with constant diffusion coefficients in both the well and the barrier layers. The composition profile after the interdiffusion is characterized by a diffusion length Ld defined as Ld=√(Dt), where D is the diffusion coefficient and t is the annealing time of thermal processing. The In mole fraction across the InGaN/GaN DFQW structure after the interdiffusion is given by:
where z denotes the co-ordinate along the crystal growth direction, Lw is the as-grown well width and x0 is the as-grown In mole fraction.
To calculate the electron and hole wave functions in QW’s, we use the multiband effective mass theory. It is a good approximation that the conduction and valence bands are decoupled. A parabolic band model and the Luttinger-Kohn Hamiltonian with strain components are used for the conduction and valence bands respectively. The envelope function scheme is adopted to describe the slowly varying part of the wavefunction. The effects of the quantum well confinement potential on the energies and envelope functions of the electron and hole subband edge at the zone centre of the Brillouin zone can be calculated separately, according to the Ben-Daniel and Duke model, using the one-dimensional Schrodinger-like equation that follows:
where ψrl(z) is the envelope function of the lth subband for electrons(r=c) or holes (r=hh, lh for heavy hole and light holes respectively); m⊥r *(z) is the corresponding carrier effective mass in the z direction; Erl is the subband-edge energy. Ur(z) is the confinement potential of the quantum well. Equation (2) is solved numerically using a finite difference method with the corresponding confinement profile. The strain of InGaN/GaN DFQW is compressive in nature, pushes the valence subbands downwards and further increases the bandgap.
The in-plane valence band dispersion of the wurtzite-InGaN/GaN QW is calculated using the parameters of GaN, InN as listed in Table I. Meanwhile, alloy properties of InxGa1-xN are obtained by linear interpolation, except for the bandgap energy, in which a bowing parameter of b=1.8eV for InGaN/GaN is used [Reference Lambrecht10]. The partition ratio for the band edge discontinuity at the heterojunction for the valence and conduction bands is assumed to be 33:67 for both QW systems.
Parameters | GaN | InN |
a (Å) | 3.190 | 3.544 |
Eg(eV) at 300K | 3.44 | 1.89 |
Δ1=Δcr (meV) | 10 | 27.8 |
Δ2=Δ3 (meV) | 5.17 | 3.9 |
mez (mo) | 0.2 | 0.11 |
met (mo) | 0.18 | 0.1 |
A1 | −6.56 | −9.28 |
A2 | −0.91 | −0.6 |
A3 | 5.65 | 8.68 |
A4 | −2.83 | −4.34 |
A5 | −3.13 | −4.32 |
A6 | −4.86 | −6.08 |
ac (eV) | −11.8 | −1.2 |
D1 (eV) | 0.7 | 0.35 |
D2 (eV) | 2.1 | 1.05 |
D3 (eV) | 1.4 | 1.32 |
D4 (eV) | −0.7 | −0.74 |
C13 (1011 dyn cm−2) | 10 | 9.4 |
C33 (1011 dyn cm−2) | 39.2 | 20 |
We assume that the injected carrier density is 1 × 1013cm−2 in the InGaN-GaN DFQW. The spontaneous emission rate is calculated by integrating the emission rate over all the electrons and holes recombination processes. The spontaneous emission is given by the following equation:
where q is the electric charge, nr is the refractive index, ε0 is the dielectric constant of the vacuum, c is the speed of light, Lw is the width of the quantum well, E is the photon energy, Ppq is the optical matrix element, and fe and fh are the Fermi distribution functions for electrons in the conduction subband and holes in the valence subband respectively. To include the spectral broadening of each transition, the total spontaneous emission spectrum in a single DFQW structure is obtained using some spectral lineshape function over all transition energies E’. Such function is given by:
where L(E-E’) is a Lorentzian line-broadening function given by
where τin is the intraband relaxation time. Then, the spontaneous emission spectrum can be obtained.
The spontaneous emission spectrum of the InGaN/GaN DFQW is shown in Fig.1. The light emission can shift the photon energy from 2.96eV to 3.02eV within the range of Ld=0-4Å. The results show that a magnitude of spontaneous emission enhancement of more than 5%. The InGaN/GaN DFQW can be used in the blue LEDs, which can adjust the wavelength of the emitted photons according to the extent of intermixing.
Experiment
The InGaN/GaN MQW in this study were grown on c-plane sapphire substrates in an atmospheric pressure MOCVD reactor. The MQW structure consist of a 570Å GaN layer, a 10 period superlattice of 25Å In0.2Ga0.8N well/70Å In0.05Ga0.95N barrier, and a 2300nm GaN:Si layer on a sapphire substrate. The SiO2 film was deposited on the sample surface using plasma enhanced chemical vapor deposition. Thermal annealing treatment was performed in a flowing N2 ambient using the diffusion furnace at 1060°C for 60minutes. Photoluminescence spectra were taken with the 325nm line of a 15mW He-Cd laser at 10K. The PL spectra of the as-grown and intermixed MQW are shown in Fig. 2. The spectrum of the as-grown MQW is dominated by the intense localized exciton (2.950eV) related to the very strong blue luminescence. After the thermal annealing, there are six weak luminescence peaks in the region of 2.955-3.45eV. The transition line (at 2.955eV) for the InGaN/GaN MQW is due to the localized exciton transition from the InGaN well region and the first quantized level of exciton localization located at 3.197eV. Two emission peaks, at approximately 3.162 and 3.071eV, correspond to the LO phonon replicas (1 and 2LO) of the main emission line at 3.251eV from the GaN layer. The energy positions are separated by 90meV which is close to the longitudinal optical (LO) phonon energy of bulk InN (86meV) and GaN (91meV) [Reference Smith, Lin, Jiang, Khan, Chen, Salvador, Botchkarev, Kim and Morkoc11]. Also, the GaN layer is observed at 3.45eV. The energy shift of localized excitons for the as-grown and the intermixed InGaN/GaN MQW is only observed to be 5meV. The effect of impurities enhanced intermixing is required to increase the rate of diffusion in the system.
Conclusions
In this paper, the InGaN/GaN MQW intermixed structure is analyzed, and used to tune the operation wavelength of light emission in the device. Theoretically, the results show that the photon energy can be shifted from 2.96 to 3.02eV within the range of Ld = 0-4Å diffusion length, resulting in an enhancement of the spontaneous light emission. Experimentally, the PL of the intermixed MQW obtained the energy shifts from 2.95eV to 2.955eV (5meV blue-shift).
Acknowledgments
This work is supported by HKU-CRCG Grants. The authors would like to thank Dr. K.W. Cheah, Hong Kong Baptist University, for helping with the PL measurements, and Dr. C. Jagadish and Dr. H. H. Tan, Australian National University, for helping with the deposition of SiO2 film.