Measurement of the non-radiative minority recombination lifetime and the effective radiative recombination coefficient in GaAs

The combination of time-resolved (TR) and power-dependent relative (PDR) photoluminescence (PL) measurements reveals the possibility of separating the radiative and non-radiative minority carrier lifetimes and measuring the sample-dependent effective radiative recombination coefﬁcient in direct bandgap semiconductors. To demonstrate the method, measurements on 2 μm thick p-type GaAs double-hetero structures were conducted for various doping concentrations in the range of 5x10 16 and 1x10 18 cm -3 . With a photon recycling factor of 0.76 ± 0.04 the radiative recombination coefﬁcient was determined to be (3.3 ± 0.6) × 10 -10 cm 3 s -1 for the structures with a doping concentration below 1 ∗ 10 18 cm -3 , whereas the effective radiative recombination parameter for an absorber thickness of 2 μm was directly measured to be (0.78 ± 0.07) × 10 -10 cm 3 s -1 . For a doping concentration of 1 × 10 18 cm -3 , the radiative recombination coefﬁcient decreases signiﬁcantly probably due to the degeneracy of the semiconductor.


I. INTRODUCTION
It has recently been shown that photon recycling must be considered for the modelling of highly-efficient optoelectronic devices. [1][2][3] Photon recycling within the active layer is effectively increasing the minority carrier lifetime associated with the radiative recombination process τ eff rad . 4 This is often taken into account by introducing a photon-recycling factor f. 4 Other authors 5 have included this factor in defining an effective radiative recombination coefficient as shown in equation (1), where N denotes the activated dopant density and Δn denotes the excess carrier density in this context.
where B eff rad depends primarily on the radiative recombination parameter of the active material, but also depends strongly on the geometry of the sample, especially the thickness of the active layer and the presence of any mirrors on the rear-side (after substrate removal). 6 It is possible to calculate both the photon recycling factor and B eff rad . 1,5 However, these calculations are strongly dependent on accurate knowledge of the sample geometry and absorption data close to the bandgap. As ternary and quaternary alloys are becoming more adopted, such absorption data may not be known to sufficient accuracy. Thus, a direct measurement approach is of great interest.
Furthermore, from an experimental material development point of view, the quality of optoelectronic semiconductor material is often rated in terms of the minority carrier recombination lifetime. This lifetime is a combination of several recombination channels: radiative and non-radiative recombination, such as Shockley-Read-Hall (SRH) recombination 7-10 including interface recombination via traps, and finally Auger recombination. While radiative recombination, i.e. the spontaneous emission of photons, is unavoidable, AIP Advances 9, 045034 (2019); doi: 10.1063/1.5051709 9, 045034-1 ARTICLE scitation.org/journal/adv non-radiative recombination can be reduced by advanced preparation, thus creating a material with a low defect density. For this optimization, a direct measurement of the lifetime associated with the non-radiative recombination channel is very helpful. There are other methods to disentangle these two recombination mechanisms, such as fitting of nonlinear high-injection time-resolved photoluminescence (TR-PL) measurements, 11 analyzing the temperature dependence of TR-PL 12 or both. 13 In order to extract the lowinjection bulk-lifetime via monoexponential decay, sufficiently high initial carrier-density must be established in order to achieve quasiequilibrium with metastable trap states, 11 and the carrier density must also be sufficiently low to not enter the non-monoexponential high-injection regime.
There are several methods to determine the radiative efficiency at low injection. 14,15,12 Recently, a simple photoluminescence (PL) method was reported based on power-dependent relative PL (PDR-PL) which measures the effective radiative efficiency from low to high injection. 5 However, all steady-state measurement methods for radiative efficiency face the challenge of calibration as they are not absolute. This is solved by the assumption of pure radiative recombination at high injection 5 or low temperature, 14 which is only valid for certain samples.
The PDR-PL measurement routine as described has two limitations: firstly, it is necessary to calculate the photon recycling factor from known material properties (i.e. absorption and emission, see for example Refs. 1 and 5) and secondly, purely radiative recombination has to be achieved for self-calibration. For the latter Auger and defect recombination have to be negligible compared to radiative recombination at high injection, which may or may not be the case depending on the material and its quality. Furthermore, band filling of the indirect bands as well as carrier escape out of the doubleheterostructure (DH) via thermionic emission must both be negligible; fortunately for the AlGaAs/GaAs material system, the barrier height is sufficient. 16 In other words, these constraints are not necessarily given for every sample investigated. Thus, two components are required to further validate the PDR-PL method: first a direct method of extracting B eff rad and second, a calibration of the PDR-PL signal. For this purpose, TR-PL can be used to quantify the effective radiative efficiency under low injection conditions when combined with the knowledge of the effective radiative recombination coefficient B eff rad . The result is a calibrated low-injection radiative efficiency used to calibrate the PDR-PL signal, which serves to quantify the high injection radiative efficiency of the system.
In the approach suggested, it is possible to determine the sample-dependent B eff rad of a DH, and its excess carrier density dependent effective radiative efficiency. However, the Auger coefficient of the material must be known, and the TR-PL lifetime has to be limited by defect recombination, i.e. not trap filling. In the case of known photon recycling properties (i.e. absorption and emission), the material-specific B rad can be extracted.
As an example, measurements are performed on p-GaAs DHs with different doping concentrations in the test layer. These samples are described in section II, as well as the measurement methods of TR-PL and PDR-PL. In section III, the calculation of B eff rad is given, followed by the calibration of the PDR-PL measurement in section IV. As both calculations depend on each other, as a last step, the calculations are performed iteratively. Results from the investigated samples are discussed and the benefits of combining PDR-PL with TR-PL are demonstrated. Finally, the results are summarized in section VI.

A. Sample description
Four Al 0.70 Ga 0.30 As/GaAs/Al 0.70 Ga 0.30 As DH were grown in an AIX2800-G4-TM metal organic vapor phase (MOVPE) reactor with a p-doping concentration of 5×10 16 cm -3 , 1×10 17 cm -3 , 2.5×10 17 cm -3 , and 1×10 18 cm -3 in the test layer, respectively. The layer stack used for each sample is shown in Figure 1. Nine periods of strain-compensated quantum wells were included to remove any luminescence coupling between the substrate and the test structure. No increase of threading dislocation density was observed by plain view cathodoluminescence measurements due to the quantum wells. The four samples were deposited on semi-insulating GaAs substrates at a growth temperature of 590 ○ C and a V/III ratio of 32.
The doping concentrations of the test layers were measured by electrochemical capacitance-voltage; a measurement uncertainty of 10% was assumed for the doping concentration. The effective lifetime was determined by TR-PL and the unscaled effective radiative efficiency by PDR-PL. 5 Both photoluminescence measurement setups are described below.

Time-resolved photoluminescence
The effective minority carrier lifetime was measured by TR-PL. The DHs were excited by a pulsed laser light employing a 640 nm diode laser (PicoQuant). The decay of the radiative band-to-band recombination was measured by a single photon counting avalanche photodiode (id Quantique) combined with a time-correlated single photon counting unit (PicoHarp 300). The excess carrier density per pulse (cm -3 ) was calculated by the laser spot size, the test layer thickness, the surface reflection, the repetition rate of the laser and the cw laser power, which was measured by a calibrated photodiode at the position of the sample. The overall time resolution given by the width of the system response function was 150 ps. The effective lifetime was extracted from the measured photoluminescence signal by fitting a single exponential decay function.

ARTICLE scitation.org/journal/adv
In order to extract the true bulk recombination lifetime from TR-PL measurements, it is essential that the lifetime does not change with excess carrier density when fitting the decay of the PL signal by a monoexponential decay function; otherwise, carrier trapping effects are impacting the signal as a function of excess carrier density. In the samples discussed, the bulk minority lifetime is extracted for low excess carrier densities Δn, where Δn is negligible compared to the doping density N. Additionally, TR-PL is not a steady state measurement. Therefore carriers trapped by shallow metastable states can have a strong impact on the TR-PL signal. 11,13,17 The mono-exponential behavior of the PL decay combined with our estimated initial carrier density suggest that the excess carrier density is significant compared to the shallow metastable trap states; this is confirmed by the agreement between the extracted and calculated values of B eff rad . This renders the decay lifetime as approximately representative of the effective bulk lifetime.

Power-dependent relative photoluminescence
The PDR-PL measurement setup and routine are described in detail elsewhere 5 and will therefore only be summarized concisely. The sample was illuminated by a 532 nm 0.5 W laser in cw operation. The resulting wavelength-resolved PL signal of the test layer was measured by a charged coupled device (Andor Technology, DU401A-BR-DD). The laser intensity on the sample was decreased stepwise over seven orders of magnitude by a neutral density filter wheel and the resulting PL signal was measured for each step. The PL signal of each intensity step was integrated over the photoluminescence peak and divided by the laser power to get the integrated so-called relative PL intensity. The relative PL signal was normalized to the maximum relative PL signal which is now defined as I rel . If there is a stable plateau at the maximum of I rel , recombination can be assumed to be purely radiative at this plateau, and the measured normalized relative PL signal is equal to the effective radiative efficiency η eff . However, 100% radiative recombination cannot be reached in every sample.

A. Effective lifetime at low injection measured by TR-PL
The decay of the photoluminescence signal was measured by TR-PL for different laser intensities, which are equivalent to different excess carrier densities Δn. The normalized TR-PL measurement signal for the lowest laser intensity of each sample is shown in Figure 2. The signals measured at low injection are all monoexponential, whereas the analysis of non-monoexponential decays would need special attention. 11,13,17 At these low excitation densities, the monoexponential time constant is assumed to be equivalent to the effective bulk minority carrier lifetime; this assumption will be tested later when comparing the extracted effective radiative recombination coefficient to that calculated from theory.
The resulting time constants from the monoexponential fits depicted in Figure 2 show an increase in effective lifetime with increasing doping concentration. Note that these lifetimes are significantly shorter than the expected radiative lifetimes assuming the often quoted value of B rad for GaAs of 10 -10 cm 3 /s. 6,18,19 This indicates that nonradiative recombination is critical in these samples

FIG. 2.
Normalized time-resolved photoluminescence decay measurements (solid lines) and fits (dashed) for the lowest laser pulse energy applied. The initial carrier density for these measurements is in the range of 10 13 −10 14 cm −3 , as shown in Figure 3.
for low injection, which can explain the increase in effective lifetime with increasing doping concentration.
The extracted monoexponential time constants from TR-PL measurements with varying laser pulse intensity are shown in Figure 3 as a function of excess carrier concentration. The lifetime is constant when Δn ≪ N is fulfilled. The monoexponential time constants at low injection are therefore equivalent to the effective bulk minority carrier lifetime at these excitation densities.

B. PDR-PL
The relative photoluminescence signal I rel normalized to 1 and measured by PDR-PL is shown in Figure 4. For all samples a constant I rel was measured at low laser power, concluding that the samples were probed in low injection conditions and a constant effective lifetime is given in this regime. Additionally signals from all samples reach a plateau at high laser power before they drop again at the highest laser power. Further analysis is not possible at this point as each PDR-PL signal has to be scaled individually. Furthermore laser power and excess carrier density are linked by the lifetime which depends strongly on doping concentration, which varies between samples. Therefore jumping to conclusions based purely on the PDR-PL measurements can be misleading.

IV. CALCULATION OF THE EFFECTIVE RADIATIVE RECOMBINATION COEFFICIENT
As mentioned, the defect lifetime of minority carriers is a crucial measure of semiconductor material quality. However, TR-PL measures the effective lifetime τ eff which depends on the effective radiative τ eff rad (i.e. the spontaneous emission lifetime weighted by the probability of photon escape out of the DH) and non-radiative τNR lifetime as given by: Typically τ NR is linked to Auger recombination (τAug) and defect recombination (τSRH) including interface recombination. The effective radiative efficiency η eff is the probability of radiative recombination which results in emission out of the front surface of the DH, and can be expressed as: Knowing τ eff and η eff of a sample thus allows for the separation of the radiative and non-radiative components to the effective lifetime.
Therefore, equations (2) and (3) can be combined, resulting in the effective radiative lifetime and the non-radiative lifetime A calculation of the effective radiative recombination coefficient B eff rad for a given structure is possible as τ eff rad can be approximated by 1 if N is orders of magnitude higher than the intrinsic free carrier concentration, which is impacted by the doping concentration in degenerate semiconductors, 20 but is in the range of 10 6 cm -3 for GaAs at room temperature. 21 Combining equations (4) and (6) reveals the possibility to calculate B eff rad . Note that B eff rad depends on the absorption and emission properties of the material, as well as its thickness since this strongly impacts photon recycling in contrast to the Einstein coefficient B rad . It can be shown that B eff rad can be calculated as follows: where B eff rad depends on the low injection effective radiative efficiency η eff ,LI . If an effective radiative recombination of 1 is reached within the measurement, η eff is equivalent to the normalized PDR-PL signal. Otherwise a calibration of the PDR-PL signal is necessary, which is discussed in the next section.

V. CALIBRATION OF THE PDR-PL SIGNAL
We will now take a look closer look at steady state conditions, under which PDR-PL measurements are performed. At steady state, the generation rate G is balanced with the recombination rate U via the continuity equation: where Δn is the excess carrier concentration in the test layer (with thickness d and associated relative absorption A) depending on illumination with a cw laser (photon energy E ph , intensity I laser ), which is therefore defined by where the effective lifetime (defined by equation (2)) is related to defect recombination τ srh , radiative recombination (equation (1)), and Auger recombination τaug = (C aug (N + Δn) 2 ) −1 , where C aug is the Auger recombination coefficient. The Auger lifetime can be calculated if the Auger recombination coefficient Caug is known, as it is for GaAs: Caug = (7 ± 4) * 10 −30 . 22 Therefore the effective radiative efficiency η eff can be calculated using equations (2) and (3) assuming a known τ srh . Vice versa: for a known η eff the excess carrier density dependent τ srh can be calculated and fitted for m different types of defects by the minority (τmin) and majority carrier defect lifetimes (τmaj) of each defect type: 7 Note that τ srh is dominated by the defect type with the shortest lifetime, which could be surface recombination. It is possible that the defect type dominating the SRH lifetime changes with Δn. The measured PDR-PL signal represents the effective radiative efficiency only if a radiative efficiency of unity is reached in the measurement. In this case, the measurement is self-calibrated. Most samples, however, do not reach complete effective radiative efficiency due to material quality, or a low bandgap (which results in significant Auger recombination) for example, and a calibration of the PDR-PL measurement is necessary. This calibration is done using a known B eff rad and a calculation of the excess carrier density dependent τ eff . The calibration factor x rad for the PDR-PL signal is defined as: x rad can be calculated iteratively with the excess carrier density dependent lifetimes for Auger, radiative and defect recombination (fit from equation (11) to measured data) at the maximum of the PDR-PL signal: Fitted SRH lifetime with two defect types (type one: black, type two: red, resulting: blue solid line), as defined in equation (11), which results in a minority carrier and majority carrier lifetime for both defects.
x rad = τ eff (Δn(I rel,max )) τ rad (Δn(I rel,max )) The process is iterative and converges for all samples investigated. A fit of the SRH-lifetime with m = 2 defect types to the data obtained for the sample with N = 2.5 × 10 17 cm −3 is shown in Figure 5.

VI. RESULTS AND DISCUSSION
Calculating B eff rad and x rad is an iterative process since both depend on each other; in other words, a decreased radiative efficiency η eff leads to a corrected B eff rad (see equation (7)) which itself leads to a corrected η eff . For all samples investigated, the iterative process converges independently of the starting parameters for τ eff and x rad .
The excess carrier density dependent SRH lifetimes (equation (6)) are calculated and shown in Figure 7. The SRH lifetimes of all samples saturate within the low injection regime and increase with excess carrier density. At the highest radiative efficiency the measured lifetime does not reach a plateau, highlighting the inadequate assumption of a purely radiative system, and afterwards subsequently drops for most samples. This could be attributed to a number of factors, including measurement uncertainty at high injection, the uncertainty of the Auger coefficient of GaAs, and lastly, carrier spilling over the barriers via thermionic emission. In contrast to low injection, the SRH lifetime of the samples decreases at high injection for increasing doping concentration.
Using the measurement results shown in Figure 3 and Figure 4, the radiative efficiencies η eff were calculated for the four samples and are shown in Figure 6 together with the calculated radiative efficiency, assuming two types of defects are limiting the SRH lifetime.

ARTICLE scitation.org/journal/adv
Comparing the onset of the increase in τ eff measured by TR-PL ( Figure 3) and τSRH measured by PDR-PL ( Figure 7) reveals a difference in excess carrier density of this point with a factor of about one order of magnitude. One reason for that is the change in excess carrier density of several orders of magnitude over one TR-PL measurement. A monoexponential fit is not adequate in the case of a non-constant lifetime. A second reason is that both measurement methods suffer from the fact that the excess carrier density is assumed constant in the probing area. Since the laser beam cross section is not a rectangular function but rather Gaussian in nature, this assumption can lead to uncertainties in the calculated excess carrier density in both measurement methods.
Lifetimes at low injection are shown in Figure 8 for all samples. The non-radiative lifetime includes defect recombination in the bulk as well as at the interfaces. A separation of both is possible by a thickness variation of the test layer. As shown in previous measurements, 5 the defect lifetime of the samples investigated is dominated by bulk recombination. The measured non-radiative lifetimes are comparable to some published data (2 μm thick p-type sample: 1.0×10 18 cm -3 , τ eff : 8.5 ns 4,23 ) but orders of magnitude smaller than others reported for samples showing a biexponential decay (2 μm thick p-type sample: 2.4×10 17 cm -3 , τ eff : 350 ns 17 ). The low growth temperature of 590 ○ C could be a reason for this limited lifetime. 24 Usually the opposite behavior is expected as the dopant itself, or impurities in the dopant source, may decrease the material quality with increasing concentration. Here, a decline of the recombination rate with doping is observed. This behavior can be explained by the uncertainty of the measurement, but a similar behavior has been reported for the EL2 defect in n-GaAs samples. 25 Calculations suggest that the formation energy of a Ga antisite or a Ga vacancy is increased with p-doping concentration for As-rich growth, 26 which would hint at these defects limiting the samples defect lifetime. Deep-level transient spectroscopy conducted on specially designed samples could clarify this behavior. Even though the SRH lifetime at low injection increases with doping concentration, the trends are completely reversed for increasing injection. This example shows that a measure of the lifetime at low injection is not sufficient for device optimization, as the operation point of minority carrier devices, such as concentrator solar cells operating under up to 2000 times concentrated solar radiation, is mostly at higher carrier densities.
Furthermore, the effective radiative recombination coefficient B eff rad of the DH-structures were deduced from equation (7) and are shown in Figure 9. For the same samples the effective radiative recombination coefficients were also calculated in a similar fashion to Lumb et al. 1 using doping dependent refractive indices of GaAs 18 to compute the relative absorption within and outside the critical angle of emission inside the DH for a particular thickness of 2 μm. Note that the model assumes that the photon recycling factor is not a function of bias (or light intensity), which is a limitation of this model; accounting for the dependence of bias or excess carrier density results in significant complications in the calculations. The bandgap is extracted optically from measured PL spectra at room temperature to account for bandgap narrowing effects. 27 The results of these calculations are also shown in Figure 9.
The agreement between the extracted values of the B eff rad and the calculated B eff rad are rather poor. However, large uncertainties of the calculated values of B eff rad arise from the impact of doping concentration on the intrinsic carrier concentration, which has been shown to be significant. 28,29 Unfortunately, this impact is not clearly visible analyzing the TR-PL lifetimes of p-GaAs double-heterostructures as a function of doping, as done in by Steiner et al., 6 since the range of doping concentration was not fully explored and there is significant scatter in the data. In fact, there is subtle evidence that the radiative lifetime is increasing beyond the expected trend for a constant value of B rad . This highlights the need for more thorough experimental validation of the impact of high carrier concentrations on both radiative recombination as well as the radiative lifetime.
The B rad of GaAs was then calculated using equation (1), since the photon recycling factor in GaAs was previously determined to be 0.74 5 or 0.80 1 for 2 μm GaAs absorber structures. For the low-doped samples, a radiative recombination coefficient B rad of (3.3±0.6)×10 -10 cm 3 /s was measured. This is around a factor of two lower than the value calculated by Varshni in 1976, 30 but in agreement with values calculated by Casey et al. 19 For a doping concentration of 10 18 cm -3 , B rad decreases to (2.1±0.4)×10 -10 cm 3 /s. For highly doped samples (N > 1×10 18 cm -3 ), a decrease of B rad is expected due to the degenerate state of the semiconductor, leading to a decrease of the slope of the absorption coefficient close to the bandgap. A similar effect was also measured for n-doped GaAs. 31

VII. CONCLUSION
The combination of time-resolved photoluminescence and power-dependent relative photoluminescence reveals a method to determine the radiative and non-radiative lifetimes of a double heterostructure. Furthermore, the combination of those two measurement methods enables the measurement of the radiative efficiency and the effective radiative recombination coefficient.
It was shown that the low injection non-radiative lifetime of the investigated MOVPE grown Zn-doped GaAs samples increased with doping concentration for the growth conditions chosen. The radiative recombination coefficient for bulk GaAs was determined to be (3.3±0.6)×10 -10 cm 3 s -1 for the structure with a doping concentration below 1 * 10 18 cm -3 , whereas the effective radiative recombination parameter for the sample with an absorber thickness of 2 μm was directly measured to be (0.78±0.07) ×10 -10 cm 3 s -1 . For a doping concentration of 1×10 18 cm -3 , the radiative recombination coefficient decreases significantly probably due to the degeneracy of the semiconductor.