The Emission Properties of Light Emitting Diodes using InGaN/AlGaN/GaN Multiple Quantum Wells

Abstract

Luminescence spectra of Light Emitting Diodes (LEDs) with Multiple Quantum Wells (MQWs) were studied at currents J = 0.15 μA - 150 mA. A high quantum efficiency at low J is caused by a low probability of the tunnel current J (which is maximum at J_m ≈ 0.5-1.0 mA). J(V) curves were measured in the range J= 10^−12-10⁻¹ A; at J > 10⁻³A they may be approximated by a sum of four parts: V= φ_k+ mkT·[ln(J/J₀)+(J/J₁)^0.5] + J·R_s. The part V ~ (J/J₁)^0.5is the evidence of a double-injection into i-layers near MQWs. Their presence is confirmed by capacitance measurements. An overflow of carriers through the MQW causes a lower quantum efficiency at high J. A model of a 2D-density of states with exponential tails fits the spectra. The value of T in the active layer was estimated. A new band was detected at high J; it can be caused by non-uniformity of In content in MQWs.

1. Introduction

Recombination mechanisms in InGaN/AlGaN/GaN heterostructures are not fully understood in spite of the great progress in the development of GaN-based light-emitting diodes (LEDs). A model of radiative recombination in 2D-structures with band tails caused by potential fluctuations was successfully applied to luminescence spectra of LEDs with single quantum wells (SQW) [1] [2] [3] [4] [5] and radiative recombination by tunneling was detected at low currents [3] [4] [5] [6] [7] [8].

There is evidence of the fact that phase separation can take place during the growth of InGaN active thin layers (quantum wells). Clusters (quantum dots) with a higher In content may be formed [9]. But no attempts have been reported to describe the spontaneous emission spectra of LEDs using a model of recombination in such clusters.

It was interesting to study details of luminescence spectra of working LEDs recently developed with multiple quantum wells (MQW) [10] [11] and to compare their properties with those of LEDs with SQWs.

In this work samples of blue and green LEDs with MQW InGaN/GaN active layers [10] [11] were studied at a wide range of currents. A model of radiative recombination in 2D-structures with band tails is applied to describe the luminescence spectra. Charge and electric field distributions for LEDs with SQWs and MQWs are compared. The mechanisms of recombination in GaN-based MQWs are discussed.

2. Experimental

Blue and green LEDs based on In_xGa_1-xN/Al_yGa_1-y/GaN heterostructures were studied [10] [11]. Structures were grown by MOCVD on sapphire substrates with an AlN buffer layer (30 nm) followed by a base n-GaN: Si layer (4-5 μm). An In_xGa_1-xN/GaN MQW structure was grown on the base. The number of periods in the MQW varied; samples with 5 periods were chosen for the study; the thickness of each period was less then 8 nm. The upper layer of Al_yGa_1-yN (50 nm) and a cap layer GaN (0.5 μm) were Mg-doped. The indium content in the wells varied, with x = 0.2-0.4. This value determined the spectral range of the luminescence, blue (x ≈ 0.2) or green (x ≈ 0.4). In order to look at details of the spectra, a wide interval of forward currents J was used (0.1 μA-200 mA); pulsed measurements were used at J>10 mA (50 Hz, 5 μs).

3. Experimental results

3.1 Luminescence spectra of LEDs

Spectra of 10 blue and 10 green LEDs were studied. The room temperature spectral maxima of the blue LEDs at J = 10 mA were hω_max= 2.64-2.67 eV, (λ_max = 465-467 nm), and the spectral width was Δ(hω)_1/2= 0.21 eV (Δ(λ)_1/2= 36-37 nm). The maxima for green LEDs were hω_max= 2.35-2.37 eV, (λ_max = 465-467 nm), spectral width Δ(hω)_1/2= 0.21 eV (Δ(λ)_1/2= 36-37 nm).

Spectra of blue and green LEDs at currents in the range J = 10⁻⁷−10⁻¹ A are shown in Figure 1 and Figure 2. The lower currents at which spectra are shown is ≈0.15 μA for blue and ≈0.5 mA for green LEDs. We have not seen room-temperature spectra of GaN-based LEDs at such low currents in the literature. The maxima of the spectra of the blue LEDs move with the current in the range hω_max= 2.57-2.67 eV, in contrast to blue SQW LEDs in which the blue maximum does not shift with the current [3] [4] [5]. There is no additional band in the yellow-green region moving with the voltage at low currents. Such a band was described as a tunnel band in blue SQW LEDs [5] [6] [7]. The maxima of the spectra of the green LEDs move in the range hω_max= 2.2-2.45 eV, a wider range than in green SQW LEDs [3] [4] [5].

Figure 1. Luminescence spectra of a blue LED N17 (room temperature) at different currents, numbers - J, mA; points - approximations by equation 2.

Figure 2. Luminescence spectra of a green diode N18 (room temperature) at different currents, numbers - current J; points - approximations by equation 2.

The low-energy sides of the spectra have an exponential form I ~ exp(hω/E₀). The parameter E₀had the value E₀ ≈ 50-60 meV, and changed only slightly with the current, as occurred in the spectra of SQW LEDs [3] [4] [5]. The high-energy sides also have an exponential form, I ~ exp~. The value of E₁ was about 40-50 meV, not equal to kT. A new band could be detected (as shoulders, hω=2.7-2.8 eV) on the high-energy tails of the spectra of green LEDs at higher currents (see Figure 2). The value of E₁in the high-energy tails of the spectra of blue LEDs was proportional to T in the range T = 220-290 K, E₁ = m·kT, m = 1.3-1.6.

3.2 Spectral shift with current and voltage

Spectra of blue LEDs at higher currents are shown in Figure 3 (J = 20-150 mA). The maxima of the spectra at constant (dc) current move to lower energies for J > 40 mA (see Figure 3a). The parameter E₁ in the high-energy exponential tails grew with this shift. The maxima of the spectra at pulsed currents (50 Hz, 5 μs) moved to higher energies; the parameter E₁ remained unchanged (see Figure 3b). Heating of LEDs at high dc currents may explain these facts. A dependence of hω_max of the energy eV (V-voltage) is shown in Figure 4. In a comparatively wide range of voltage this function is linear, but the slope of the line is << 1, (in contrast with the tunnel band reported in [3] [4] [5] [6]). Filling of the tail states in the active layer causes this shift.

Figure 3a. Spectra of a blue LED N13 at high currents in dc conditions; points - approximation by equation 2.

Figure 3b. Spectra of a blue LED N13 at high currents in pulse conditions (50 Hz, 5 μs); points - approximation by equation 2.

Figure 4a. Dependence of spectral maxima of a blue LED B17 versus voltage in dc conditions.

Figure 4b. Dependence of spectral maxima of a green LED G18 versus voltage in dc conditions.

3.3 Current-voltage characteristics

Current-voltage characteristics J(V) of blue and green LEDs are shown in Figure 5. There is an exponential part at low currents, J < 10⁻⁷ A at 300 K, a steep exponential growth in the range V = 2.3-2.7 V, a linear part at higher currents, J > 20 mA. Low currents can be understood as a tunnel component; tunnel currents in these LEDs play some role at J 3-4 orders of magnitude lower than that for SQW-based LEDs [3] [4] [5] [6] [7]; J(V) curves of SQW-based LEDs are shown in Figure 5 for comparison. The difference can be explained as a consequence of a wider active layer of MQW structures.

Figure 5. Current-voltage characteristics of blue (B1, B2) and green (G1, G2) LEDs with single (B1, G1) and multiple (B2, G2) QW at room temperature (solid curves) and at 80 K (dashed).

A good approximation of the J(V) curves for MQW LEDs was made when not only a series resistance R_s at the linear part at higher J was taken into account, but also the quadratic part: J ~ (V-V₁)². The fit of the curve J(V) at J > 0.1 mA by the equation:

(1)

is shown in Figure 5. The fitting parameters are φ_k (contact potential), E_J(E_J=c·kT, c = 1-2), J₀ (saturation current), and J₁, R_s. One part ~(J/J₁)^0.5 is sufficient between an exponential (injection) and a linear parts, in the usual working current range J=2-30 mA..

3.4 Quantum efficiency

Dependencies of the integrated intensity Φ(J) and external quantum efficiency η_e(J) = eΦ/J versus J are shown in Figure 6. Measurements of η_e were done by a method described in [12]. The efficiency η_e(J) has a maximum at low currents J≈0.5-1.0 mA, at the start of the steep exponential growth of J(V). The value of η_e goes down logarithmically with J at high currents (linearly with V).

Figure 6a. Dependence of integrated intensity and quantum efficiency for a blue B17 diode on the current.

Figure 6b. Dependence of integrated intensity and quantum efficiency for a green G18 diode on the current.

3.5 Distribution of charged centers

The distributions of charged centers in p- regions of MQW and SQW InGaN/AlGaN/GaN p-n- heterostructures are shown in Figure 7 (see the measurement method in [13]). The MQW LEDs space charge is wider than that of SQW LEDs [3] [4] [5] [6]; in both cases the width for green LEDs is wider than for blue ones. This fact corresponds to a low probability of tunneling in the MQW LEDs.

Figure 7. Distributions of charged centers in p-regions of p-n heterojunctions; points - values of N_A ^- at V = 0; begining of abscissa is at the n-interface. Blue (1,3) and green (2,4-6) LEDs with single (1,2) and multiple (3-6) quantum wells.

It seems that high Mg-doping of p-AlGaN and GaN layers is more difficult for higher In concentration in InGaN active layers.

4. Discussion

4.1 Spectral fit by the model of 2D-density of states with low-energy tails

We describe the spectra with a model previously applied for fitting the spectra of SQW LEDs [1] [2] [3] [6]. The model implies that an effective radiative recombination takes place when carriers of both signs are injected into the active layer at voltages on the layer U < V. The value of U is close to φ_k. Optical transitions at hω are going between states E_(c) and E_(v) in the tails of the 2D-structure caused by potential fluctuations. A model 2D joint density of states is

(2)

an effective energy gap E_g ^eff is E_g ^eff= E*_c- E*_v . The parameter E₀ is determined by potential fluctuations. A discussion of possible sources of these fluctuations (well and barrier inhomogeneities, fields due to impurities, or piezoelectric effects) will be published elsewhere.

The spectral intensity I( hω) is proportional to the Fermi-functions of electrons and holes with quasi-Fermi levels F_n, F_p as parameters (details in Ref. [5]):

(3)

Examples of the fit are shown in Figure 1 and Figure 2; parameters of the fit are summarized in Table 1. It is possible to describe a change of hω_max in a certain range of J by changes of the parameter F_n - the parameters E_g ^eff, E₀ and E₁= m·kT may be unchanged. This is evidence of the fact that the mechanism of recombination in the 2D tail-states is not changed.

Table 1 Fitting parameters for approximating blue LED spectra

4.2 Parameters of the approximation

This description is valid only in some range of J. The parameter E₁= m·kT changes at higher J. This is caused first of all by heating at J > 10 mA. Curves of approximation are shown in Figures 3a, 3b. It is possible to describe shifts of pulse spectra without changing the parameter E₁, and shifts of the dc spectra - without changing the parameter m, supposing a change of temperature T. The low energy shift of spectral maxima at higher J corresponds to the empirical equation of Varshni:

(4)

The parameters in this equation are E(0)=3.07 eV; α=12.8·10⁻⁴ eV/K; β=1190 K (see [14]).

4.3 Possible origin of the new spectral band

The short wavelength tail changes not only by heating, but also with the current. It depends also on the new spectral band that is clearly seen in the logarithmic scale on the spectra (see Figure 2). We suppose that this band is caused by large-scale inhomogeneities - separation of phases with different content of indium in In_xGa_1-xN. Models of recombination either in the band tails or quantum dots were examined as an alternative in the discussion at the Tokushima Conference [10]. It seems that our results confirm that both possibilities are realized. The proof of our supposition may be obtained by studying the luminescence of MQWs with microstructures revealed by electron microscopy and SIMS.

4.4 Maximum quantum efficiency

The problem of the maximum η_e versus J is a very important one. It is connected to the number of QWs and to the properties of p-AlGaN layers made by various technologies.

This maximum can be understood as follows. Nonradiative channels of recombination (for example, tunnel recombination) take place at low J. Electrons are filling the active MQW layer by injection. This is the region of maximum η_e. At higher J electrons overflow the active layer and are pulled by an electric field into the i-layer and p-AlGaN (see an analogous model in [15]). The quadratic part of the J(V) and the linear dependence η_e(V) show the role of electric field and of a drift component of the current.

5. Conclusions

1. 1. Luminescence spectra of LEDs based on InGaN/AlGaN/GaN heterostructures were studied in a wide range of currents, down to J = 10⁻⁷ A at room temperature. Filling the tail states in MQWs causes shifts of the spectral maxima with the voltage.

2. 2. The model of recombination in a 2D-structure with exponential band tails describes the spectra with good accuracy. A new spectral band was detected; it is supposed that this band can be caused by phases of higher indium concentrations in InGaN QWs.

3. 3. The tunnel component of the current is 3-4 orders of magnitude lower than in analogous SQW LEDs. The LEDs have a current component described by a model of double injection of carriers into the i-layers adjacent to the active MQW layer.

4. 4. The quantum efficiency has a maximum depending on the current. Overflow of electrons through the active layer can cause lower quantum efficiency at higher currents; an electric field is pulling the carriers into compensated i-layers.