From time- and spatially resolved optical measurements, we show that extended defects can have a large effect on the charge-carrier recombination in II–VI semiconductors. In CdTe double heterostructures grown by molecular beam epitaxy on the InSb (100)-orientation substrates, we characterized the extended defects and found that near stacking faults the space-charge field extends by 2–5 μm. Charge carriers drift (with the space-charge field strength of 730–1,360 V cm−1) and diffuse (with the mobility of 260 ± 30 cm2 V−1 s−1) toward the extended defects, where the minority-carrier lifetime is reduced from 560 ns to 0.25 ns. Therefore, the extended defects are nonradiative recombination sinks that affect areas significantly larger than the typical crystalline grains in II–VI solar cells. From the correlative time-resolved photoluminescence and second-harmonic generation microscopy data, we developed a band-diagram model that can be used to analyze the impact of extended defects on solar cells and other electronic devices.
To investigate the fundamental performance limits of II–VI solar cells, single crystals and epitaxial heterostructures have recently been used in solar cells that reached an open-circuit voltage of >1 V.1,2 Epitaxial materials are also used in other optoelectronic devices, such as sensors and detectors. All applications require an understanding of defects and their impact on electronic properties such as charge-carrier transport and recombination. Photoluminescence (PL) microscopy and imaging were applied to study the transport in III–V3–6 and II–VI heterostructures7–9 and in polycrystalline thin films.10 We extended this approach by using correlative second-harmonic generation (SHG) measurements. Interaction between the incident light wave and a DC electric field in the material generate EFISH (electric-field-induced second-harmonic) signals, and the optical EFISH measurements allow the analysis of the space-charge fields (SCFs) near the extended defects.
We report the electronic properties of CdTe heterostructures grown by molecular beam epitaxy (MBE) on InSb (100) substrates. Polar (100) surfaces have high surface energy and can nucleate the extended defects, such as twins and stacking faults.11 We found that in the high-quality epitaxial heterostructures, a single extended defect can affect carrier recombination in an area larger than 100 μm2. We also found that drift due to the SCF can increase the impact of recombination.
Double heterostructures (DHs) were grown by MBE and characterized as described in Zaunbrecher et al.12 Confocal PL (c-PL) was used to identify the extended defects and to determine their density.12 The incorporation of 0.8% Se in the CdTe absorber and the buffer provided better lattice matching to InSb and allowed the growth of thicker absorbers without strain relaxation. For the DH studied here, the CdTe epilayer thickness was 5 μm, the extended defect density estimated from c-PL was ≈104 cm−2, and the minority-carrier lifetime was 560 ns. Passivation with 30-nm-thick epitaxial MgxCd1-xTe (x = 0.25–30) reduced the interface recombination velocity to 160 cm s−1.12 The DHs were not intentionally doped; the background doping was p-type at ≈3 × 1014 cm−3.
The nonlinear optical microscopy employed here uses the local generation/local signal collection geometry (Figure 1(a)). Such measurements have two-times-higher spatial resolution than microscopy in the local generation/global collection configuration.6 We used 1030 nm excitation with 300-fs pulses at 1.1 MHz repetition rate. For the PL/time-resolved PL (TRPL) measurements on semiconductors with a bandgap of 1.5 eV, this wavelength allows two-photon excitation (2PE).13 The excitation energy was 0.05–0.20 pJ pulse−1, and the excitation beam was focused with a Zeiss A-Plan 40× NA0.65 objective.14 As shown below, such excitation corresponds to low-injection conditions. Photon counting detectors (PDM from Micro Photon Devices) with 10-nm bandpass interference filters were used in PL (840 nm) and SHG (515 nm) measurements. Time-correlated single-photon counting (TCSPC) with a PHR800 router and Picoharp 300 (Picoquant) afforded ≈90 ps time resolution that was mostly limited by the time response of the single-photon detectors. The sample was translated with mechanical xyz stages, and at each coordinate the TRPL/SHG data were integrated for 0.1–1 s. Because deconvolution was not used in data analysis, the PL and SHG data were time-integrated in 0.1-ns steps. Overall, the microscope provides an estimated 0.5-μm lateral (x and y) and 0.1-ns temporal resolution for the PL and SHG measurements. The axial (z) resolution for 2PE PL microscopy in a semiconductor with a refractive index n = 2.85 is ≈3 μm.14,15 The axial resolution for SHG microscopy is ≈0.5 μm because the SHG signal at 515 nm is strongly reabsorbed in CdTe.
(a) Experimental setup and DH sample structure. (b) and (c): Integrated PL and SHG images for the same area. Letters (A–H) indicate pixel locations for which TRPL decays are shown in (f). Dashed lines V, L, and S in (b) mark directions where charge-carrier transport was analyzed in Figure 3 and where linear intensity profiles are shown for PL (d) and SHG (e). Dashed lines in (d) and (e) show the Gaussian fits to the data, and legends indicate the full-width-at-half-maximum values for the Gaussian fits.
(a) Experimental setup and DH sample structure. (b) and (c): Integrated PL and SHG images for the same area. Letters (A–H) indicate pixel locations for which TRPL decays are shown in (f). Dashed lines V, L, and S in (b) mark directions where charge-carrier transport was analyzed in Figure 3 and where linear intensity profiles are shown for PL (d) and SHG (e). Dashed lines in (d) and (e) show the Gaussian fits to the data, and legends indicate the full-width-at-half-maximum values for the Gaussian fits.
To characterize the influence of the extended defects on the electronic characteristics of II–VI heterostructures, Figure 1 shows PL (b) and SHG (c) images near one extended defect. The images were obtained by time-integrating the PL and SHG data measured in an area 20 × 25 μm2 with a pixel size of 0.25 μm. Several TRPL decays at individual pixels are shown in (f). At the center of the extended defect (pixel E), the PL decay can be fit with a 0.25 ns lifetime. This short lifetime indicates fast nonradiative recombination due to the Shockley-Read-Hall (SRH) recombination centers associated with the extended defect. When the distance to the center of the image increases, the decays have longer lifetimes and non-exponential shapes. For example, the decays at pixels A and H (same distance to the center) are nearly identical. We have previously shown that when the carrier diffusion length is >1 μm, 2PE TRPL microscopy data in the local excitation/local collection geometry indicate not only recombination but also charge-carrier transport.16 These characteristics were also analyzed with numerical simulations17 and analytical modeling.18 Therefore, the decays in Figure 1(f) are dominated by charge-carrier transport. Unlike the case analyzed earlier,16–18 transport near the extended defect is not isotropic due to the recombination sink (extended defect). Therefore, we apply time-resolved microscopic analysis.
Figure 1(c) shows an SHG image measured simultaneously with the PL data. The weaker SHG background is present everywhere in the image, but SHG is approximately three times stronger near the extended defect (see the SHG intensity profiles in Figure 1(e)). This enhancement is attributed to the EFISH mechanism, where the second harmonic of the incident light is generated due to the SCF near the extended defect in the heterostructure. EFISH microscopy has been applied to organic bulk heterojunctions,19 organic transistors,20 and metal-silicon junctions.21 Our data show that EFISH microscopy could be applied for the SCF analysis near the extended defects. Because the SCF affects the carrier dynamics, the field strength can be estimated from the TRPL data (see Figure 3). This is an important advantage of correlative microscopy, where the SHG/EFISH data indicate the spatial distribution of the SCF, and analysis of the carrier dynamics provides the strength of the SCF.
In contrast, when similar measurements were applied near the threading-dislocation extended defects in DHs grown on (211)B-orientation CdTe substrates, the second-harmonic enhancement was not observed.9 Comparison suggests that space-charge fields are absent (or are weaker) near threading dislocations. Transport near dislocations in the GaAs/GaInP DH could also be described by the diffusion model.4–6
When Gaussian fits are used to describe the intensity profiles in PL (Figure 1(d)) and SHG (Figure 1(e)) images, the full width at half maximum (FWHM) is larger for the PL data (e.g., in the S direction, FWHMPL = 10.9 μm and FWHMSHG = 2.2 μm). Because the reduced PL intensity at the center is attributed to recombination, the data indicate that carriers can reach nonradiative recombination sites over a distance of >10 μm. Such data suggest transport to the SRH recombination centers.3–7,9,12 Drift and diffusion processes in the dynamics could be resolved from the time-resolved microscopy data described below.
The time-dependent PL and SHG images in Figure 2(a) were obtained by normalizing the intensity in each frame. The size of the dark area in the PL images (corresponding to low PL intensity) is increasing over time. Because the charge-carrier lifetime at the center is very short (0.25 ns) due to the fast nonradiative recombination, carriers generated in other areas can diffuse (and drift; see below) toward the extended defect and nonradiatively recombine there. As shown in Figure 2(b), the Gaussian fits to PL intensity profiles were used to numerically evaluate the time-dependent FWHM for the region affected by recombination. For this orientation (vertical line in the middle of the frame), the FWHM increases from 5.2 ± 0.7 μm at 0.1 ns to 9.2 ± 0.2 μm at 1.9 ns (1.8 × increase) and 16.9 ± 1.2 μm at 20.9 ns (additional 1.8 × increase). At > 20 ns, almost all of the 20 × 25 μm2 region is affected by recombination due to the extended defect.
(a) Time-resolved PL and SHG images at times indicated in the insets. The data in each image are normalized to illustrate the time-dependent increase in the area affected by recombination. Vertical line in the frame “PL 20.3 ns” indicates the location for which several time-dependent linear intensity profiles are shown for PL (b) and SHG (c). Gaussian fits in (b) and (c) are shown as dashed lines, and FWHMs are given in the legends.
(a) Time-resolved PL and SHG images at times indicated in the insets. The data in each image are normalized to illustrate the time-dependent increase in the area affected by recombination. Vertical line in the frame “PL 20.3 ns” indicates the location for which several time-dependent linear intensity profiles are shown for PL (b) and SHG (c). Gaussian fits in (b) and (c) are shown as dashed lines, and FWHMs are given in the legends.
In contrast to the PL data, the size of the bright feature in the SHG images in Figure 2 does not increase (FWHM = 2.9 ± 0.1 μm at 0.1 ns and 3.0 ± 0.1 μm at 1 ns), and the SHG time dependence is limited by the temporal response of the detector. Therefore, SHG/EFISH probes the interaction between the space-charge field in the material and the optical field,19–21 whereas PL probes carrier dynamics.3,7,9
To describe the charge-carrier transport characteristics, in Figure 3 we analyze the time-dependent FWHM for the region where the PL intensity is reduced due to nonradiative recombination. At 0–2 ns, the size of the defect feature increases more rapidly and then continues to expand more slowly. The second phase of the increase is linear on the (FWHM)2 vs. time scale and independent of the direction in the epilayer plane. By fitting these data to the two-dimensional diffusive transport model, we determine the charge-carrier diffusion coefficient D by using3,7
where t is time and FWHM0 is the initial FWHM for the carrier distribution. The fits in Figure 3(a) yield D = 6.7 ± 0.3 cm2 s−1 (V direction), D = 6.1 ± 0.2 cm2 s−1 (S direction), and D = 7.7 ± 0.2 cm2 s−1 (L direction). The averaged charge-carrier mobility calculated using Einstein's relationship is μ = 260 ± 30 cm2 V−1 s−1. Figure 3(b) shows that when the excitation power is reduced 4 times (due to 2PE, injection is reduced 16 times), the time dependence for the FWHM is essentially unchanged. An intensity-independent carrier mobility indicates that the analysis is for electrons (minority carriers). Low injection is also suggested from the two-photon absorption coefficient (β ≈ 25–30 cm GW−1 at 1000–1400 nm)22 and the excitation peak power (I ≈ 0.8 GW cm−2). The product β × I ≈ 20 cm−1 is equivalent to the linear absorption coefficient. Because the excitation photon flux is 1.3 × 1015 photons cm−2 and the absorption coefficient is small, injection is estimated to be lower than the background doping.
(a) Time-dependent (FWHM)2 and FWHM (inset) for the three directions indicated in Figure 1(b). (b) Comparison of FWHM data when injection differs by 16×. Inset in (b) shows band diagram that summarizes the results of this study.
(a) Time-dependent (FWHM)2 and FWHM (inset) for the three directions indicated in Figure 1(b). (b) Comparison of FWHM data when injection differs by 16×. Inset in (b) shows band diagram that summarizes the results of this study.
The initial growth phase in the FWHM data in Figure 3 can be related to drift in the SCF. If the extended defects have positive charge, the minority carriers (electrons) will drift toward this area and recombination will affect a larger area than it would in the absence of SCF. For example, in the V direction, the size of the recombination region approximately doubles in size due to drift (see Figure 3(b)).
The inset in Figure 3(a) shows that the drift velocity (estimated as slope of FWHM vs. time) is the highest in the V direction. We estimate the SCF strength ESCF from
If we assume that drift and diffusion mobilities are the same, the fits indicate ESCF,V = 1360 ± 80 V cm−1, ESCF,L = 730 ± 35 V cm−1, and ESCF,S = 820 ± 47 V cm−1.
The SCF near the extended defects could be created by charged impurities, strain,23 or charge separation. The charge separation at pairs of partial dislocations was predicted from first-principles calculations.24 Stacking faults are usually bordered by partial dislocations, and such dislocations could be recombination centers.25,26 The charge-carrier lifetimes could be reduced due to recombination at dislocations in II–VI heterostructures.7,8 Calculated band-bending due to the charge separation24 is similar to our experimental results, and we also observe fast nonradiative recombination at the extended defects. It was suggested that the extended defects in CdTe could include wurtzite layers formed by closely spaced lamellar twins.27 According to the calculated band alignment at the zinc-blende/wurtzite interface,27 the wurtzite region would repel electrons, but we found that the SCF leads to electron drift toward the extended defects. Therefore, our data are in agreement with the model of recombination at partial dislocations/stacking faults, but not with the model of recombination due to the defects at the zinc-blende/wurtzite interface.
Data from our study can be summarized in the band diagram shown in the inset of Figure 3(b). The charge-carrier lifetime in the DH (τB = 560 ns) is limited by the SRH recombination due to point defects in the bulk and at the interfaces.12 The lifetime near the extended defects is reduced to 0.25 ns; therefore, the extended defects are recombination sinks. The carrier diffusion length L = sqrt(D τB) = 20 μm, and the extent of recombination due to a single defect is consistent with this estimate. The microscopic SHG/EFISH images indicate the spatial distribution of the space-charge field. Depending on the orientation in the epilayer plane, the SCF can be described by the FWHM = 2–5 μm. The free carriers in the DH would screen the fixed charges. The screening length (depletion width) with the background doping of 3 × 1014 cm−3 is about 1 μm. Because the SCF extends up to 5 μm, our data suggest that recombination could occur at paired extended defects identified in an earlier study.12
The band bending (ranging from 0.2 to 0.4 eV) is calculated from the dimensions and strength of the SCF. Because the minority carriers (electrons) are attracted toward the SCF, we infer the downward direction of the band bending for the conduction band. Because the PL emission energy was the same near the extended defect, the same band bending is present for the valence band.
The band diagram is quite similar to that of the grain boundaries in polycrystalline CdTe.28 It is generally accepted that grain-boundary recombination can limit the open-circuit voltage and efficiency of II–VI solar cells.29,30 Our data show that recombination at other extended defects could also be important—we found that in the high-quality heterostructures (τB = 560 ns) a single extended defect could cause recombination in an area much larger than the typical crystalline grain size in CdTe, Cu(In,Ga)Se2, kesterite, perovskite, and other thin-film solar cells. The recent research efforts have increased the charge-carrier lifetimes in polycrystalline materials (e.g., τB = 10–100 ns in CdTe31 and τB ≈ 400 ns in Cu(In,Ga)Se232), and the passivation of other extended defects—not only grain boundaries—could be increasingly important for thin-film solar cells. A balance of the crystalline grain size vs. the extended defect density is carefully considered for multicrystalline silicon solar cells,33 and thin film solar cell technologies could follow similar approach.
This work was supported by the U.S. Department of Energy under Contract No. DE-AC36-08-GO28308 with the National Renewable Energy Laboratory. The U.S. Government retains and the publisher, by accepting the article for publication, acknowledges that the U.S. Government retains a nonexclusive, paid up, irrevocable, worldwide license to publish or reproduce the published form of this work, or allow others to do so, for the U.S. Government purposes.