Using photoluminescence spectroscopy, we construct a recombination model for state-of-the-art CdTe solar cells doped with Cu. We observe that Cu on Cd sites form a dominant acceptor state about 150 meV from the valence band. Although it is intuitive that this state can increase hole density, we also find that this relatively shallow dopant can also limit lifetime. Consequently, CdTe solar cells doped with Cu could have a lifetime limitation inversely proportional to the hole concentration.

II–VI semiconductors are used in lasers, light-emitting diodes, detectors, and solar cells.1 The most common approach to dope II–VI semiconductors is to substitute type IB or type V elements on the cation or anion sites, respectively. Because of the relatively large effective hole mass in CdTe, these states are approximately 50–200 meV above the valence band.2 Type IB dopants such as Ag, Au, and Cu are more mobile and easier to insert at lower temperatures than type V dopants.3–7 However, we find that these relatively shallow acceptors can influence recombination.

This has significant implications particularly for photovoltaics. CdTe solar cells are challenging the dominance of polycrystalline silicon. Record CdTe solar cell efficiency is higher than that for multicrystalline Si (21.5% and 20.8%, respectively),8 and costs at current manufacturing efficiencies are already similar or lower.9 Yet, CdTe technology has room to improve. The open-circuit voltage (Voc), fill factor, and efficiency depend on hole concentration, which for CdTe is currently on the order of 1014 cm−3. If the hole concentration can be increased by two orders of magnitude without decreasing lifetime, VOC can be significantly increased and efficiencies of 25% are possible.

Throughout the history of CdTe solar technology, Cu has been critical to attain hole concentration and to form the back contact.10–13 However, the role of type IB dopants is complex. First-principles calculations indicate that increasing Cu levels can lead to formation of Cui donors that compensate CuCd acceptors.14 Lifetime is affected by Cu levels, and Cu is mobile and can segregate to grain boundaries or adjacent material layers.10,13,15–17 Cu concentration is often 1017–1018 cm−3, yet hole concentration is several orders of magnitude lower.15,17 Here, we construct a recombination model for state-of-the-art CdTe devices and demonstrate that the same CuCd site that increases hole density can also reduce lifetime.

Five polycrystalline films with large grains (6–10 μm) were chosen for photoluminescence (PL) analysis, because for smaller grains inhomogeneous broadening could obscure some features in the PL spectra.6,11,12,18,19 The films were grown in a superstrate configuration by vapor-transport deposition, annealed in a CdCl2 containing atmosphere, and doped with Cu from a wet-applied Cu source followed by a thermal bake. For polycrystalline CdTe, the resulting devices have an extremely high Voc of 890–900 mV. A reference sample without intentional Cu doping had a much lower Voc of 777 mV. The high Voc of these Cu-doped solar cells can be attributed largely to an increase in the hole concentration and long lifetimes in the CdTe absorber layer.20,21 Capacitance-voltage measurements indicate a net acceptor concentration of 1 × 1015 cm−3, which is approximately an order of magnitude greater than typical CdTe solar cells.10,21 PL and electrical measurements indicate that CuCd sites cause this hole concentration.4,11,22,23

For PL emission spectroscopy, excitation was at 632.8 nm with a continuous wave (cw) HeNe laser. Time-resolved photoluminescence (TRPL) was measured using time-correlated single-photon counting and laser excitation at 640 nm with 0.3-ps pulses fired with a frequency of 1.1 MHz.24 The photons were detected with a Si avalanche photodiode and appropriate bandpass filters, or with a monochromator and photomultiplier tube.25 The resulting decay curves were fit with a bi-exponential function using deconvolution26 

(1)

The initial component τ1 is attributed to the combination of interface recombination, drift, and diffusion, whereas τ2 better reflects the bulk lifetime.27–30 For the samples studied here, τ2 ranges from 13 to 24 ns, which is state of the art for polycrystalline CdTe.21 

To understand how hole concentration and lifetime may be further increased, we use energy- and time-resolved spectroscopy to form a detailed recombination model. Electron-hole pairs were generated in a Beer's law distribution with a penetration depth of ∼500 nm (2.3/α640 nm, where α is absorption coefficient31). By measuring time-resolved emission spectra starting from such well-defined electron-hole pair distribution, we can analyze emission signatures and carrier dynamics even when spectral lines are heterogeneously broadened due to the S interdiffusion at the CdS/CdTe interface.19 Fig. 1 illustrates the time-resolved emission spectra at 5 K for a sample doped with Cu. The PL decay is much faster for the high-energy (exciton) region of the spectrum than for the lower-energy (defect) region. This leads to a time-dependent PL spectral shift attributed to faster decay of the excitonic PL and donor-acceptor pair (DAP) recombination. The red shift at >17 ns is largely attributed to DAP recombination. DAP recombination can be analyzed from32 

(2)

where EZPL is zero phonon energy, Eg is bandgap, Ea/Ed are acceptor/donor energies, e is electron charge, ε is dielectric constant, and ε0 is permittivity of vacuum. As carriers recombine, the average distance between donors and acceptors, r, increases so that the spectral shift toward the red is observed. To determine all terms in Eq. (2), we varied temperature and excitation power.

FIG. 1.

Time-resolved PL emission spectra measured at different times (see legend) after excitation on the junction side of the device. T = 5 K and the average excitation power 1 mW.

FIG. 1.

Time-resolved PL emission spectra measured at different times (see legend) after excitation on the junction side of the device. T = 5 K and the average excitation power 1 mW.

Close modal

Figure 2 illustrates the temperature-dependent shift of the exciton band. Data show that the device has different Eg values at the junction and back-contact regions. The difference ΔEg ≈ 40 meV is attributed to formation of a CdSTe alloy at the junction.19 Ellipsometry data indicate this energy shift corresponds to a CdSxTe1−x composition with x ≈ 0.05.33 The measurements described below were repeated for both sides of the photovoltaic device: the CdS/CdTe junction where CdSTe alloy properties were evaluated, and the CdTe film side (see inset in Fig. 2). The absolute Eg values were estimated from the data in Fig. 2 by adding the free exciton binding energy (≈10 meV)34 and exciton-defect binding energy (10–20 meV).35,36 From 50 to 250 K, the Eg shift is measured to be approximately 0.3 meV/K. A similar result was reported for CdTe single crystals.32 

FIG. 2.

Exciton PL emission temperature dependence for measurements from the film (CdTe, green) and junction (CdSTe, red) sides as illustrated in the inset. When lower excitation power is used (black), the exciton PL energy is essentially identical.

FIG. 2.

Exciton PL emission temperature dependence for measurements from the film (CdTe, green) and junction (CdSTe, red) sides as illustrated in the inset. When lower excitation power is used (black), the exciton PL energy is essentially identical.

Close modal

Next, we determined EZPL from the analysis of well-resolved phonon replicas in Fig. 3. Peaks from 1.35 to 1.45 eV can be fit with five Gaussians spaced by 21.3 meV, which is the LO phonon energy in CdTe.35,36 Therefore, the sub-bands are phonon replicas for a defect with EZPL = 1.450 eV. The intensity distribution shown in the inset is well described by the Poisson distribution with the Huang-Rhys factor S = 1.58 ± 0.02. EZPL and S values are in agreement with Stadler et al., who derived EZPL = 1.45 eV and S = 1.60 for single-crystal CdTe with Cu.32 The data also agree with other reports for the Cu on Cd site.4,22,23,37 The clear PL signal from CuCd indicates that this dopant is important for the recombination kinetics.

FIG. 3.

Low-temperature (5 K) PL emission spectrum from the film side of the device. Excitation power 10 μW. Gaussian sub-bands (blue) produce an excellent fit (red) for DAP emission. The inset: symbols show the intensity of Gaussian sub-bands, I(n). Solid line shows fit to Huang-Rhys model I(n)=I0eSSnn!.

FIG. 3.

Low-temperature (5 K) PL emission spectrum from the film side of the device. Excitation power 10 μW. Gaussian sub-bands (blue) produce an excellent fit (red) for DAP emission. The inset: symbols show the intensity of Gaussian sub-bands, I(n). Solid line shows fit to Huang-Rhys model I(n)=I0eSSnn!.

Close modal

Having determined the EZPL, we examined the DAP binding energy by varying the electron-hole injection level. Fig. 4 shows PL emission spectra from the junction side when excitation power was varied from 3 nW to 11 μW. Heterogeneous broadening does not allow observation of the phonon sub-bands,6,11,12,18,19 but the DAP band is similar to that in Fig. 3. As shown in the inset, the Coulombic interaction for the DAP leads to the blue shift of the PL emission. The dependence of the broad DAP peak on the excitation power, IEXC, can be analytically described as38 

(3)

where D is a coefficient, hνmax is the experimental DAP energy, hν = Eg − Ea − Ed is the limiting photon energy of infinitely distant pairs, and hνB = EB + hν with EB = e2/4πεε0RB (RB is Bohr radius of the shallow impurity). Data in the inset indicate hν = 1.336 eV. Eq. (3) predicts that the blue shift due to Coulombic interaction is at most the binding energy of the shallow impurity. For CdTe, shallow donors form with Ed ≈ 10 meV below the conduction band.32 Consistent with this estimate, when excitation power is increased by almost four orders of magnitude, the DAP band shifts by ≈10 meV. Fitting with Eq. (3) gives hνB = 1.36 eV and RB = 66 Å. The large Bohr radius is consistent with the small donor ionization energy. Taken together, the data indicate Ea + Ed = 0.16 eV, and most of this energy is attributed to Ea because Ed is small.

FIG. 4.

PL emission spectra at 5 K from the junction side of the device as a function of excitation power (see legend). Arrows indicate the DAP band maxima. The inset shows the DAP PL shift with excitation power.

FIG. 4.

PL emission spectra at 5 K from the junction side of the device as a function of excitation power (see legend). Arrows indicate the DAP band maxima. The inset shows the DAP PL shift with excitation power.

Close modal

Next, we examine the concentration of DAPs from the TRPL data measured at 5 K (Fig. 5). DAP recombination can be described as consisting of two steps—electron (minority carrier) capture by donors, followed by tunneling recombination.39 These processes occur on very different time scales. By examining the DAP PL rise, we find that the donor sites fill at a rate of (400 ± 30 ps)−1 independent of injection (Fig. 5, inset). Consequently, it appears that the donor density is high and the donor states are not saturated for the injection levels used in the experiments. DAP recombination is much slower (>700 ns time scale). We determine the acceptor concentration, Nt, by fitting the data to the DAP recombination model39,40

(4)

where ⟨A(t)⟩ is the ensemble average for the probability that the electron is on the donor, and the recombination rate W(r) exponentially depends on the DAP separation, r, where the maximal rate is Wmax39,40

(5)

When fitting data to the model of Eq. (4), the best fit occurs for Nt ≈ 1016 cm−3. Experimental TRPL decays for the DAP band were essentially unchanged from 5 to 150 K, and the PL emission intensity was also approximately constant in this temperature range.41 These data are consistent with DAP tunneling recombination at low temperatures because tunneling rates are expected to have weak temperature dependence.

FIG. 5.

Donor-acceptor pair (DAP) PL kinetics (junction side) at 5 K. The average excitation power was 0.5 mW (black), 1 mW (blue), or 3 mW (green). The solid red line shows the fit to Eq. (4). The inset shows the rise of the DAP PL on a much faster time scale. The instrumental response function used in the fitting is shown in grey. All four decays (excitation power indicated in the inset) are fit with a 400 ± 30 ps rise time.

FIG. 5.

Donor-acceptor pair (DAP) PL kinetics (junction side) at 5 K. The average excitation power was 0.5 mW (black), 1 mW (blue), or 3 mW (green). The solid red line shows the fit to Eq. (4). The inset shows the rise of the DAP PL on a much faster time scale. The instrumental response function used in the fitting is shown in grey. All four decays (excitation power indicated in the inset) are fit with a 400 ± 30 ps rise time.

Close modal

Shockley-Read-Hall (SRH) recombination at higher temperatures was analyzed by comparing TRPL lifetimes for undoped single-crystalline CdTe, a polycrystalline CdTe device without intentional Cu, and a high-voltage device with Cu-doping (Fig. 6). For the single-crystal reference sample, the lifetime has relatively weak temperature dependence and changes from 360 to 383 ns at 275–200 K.42 

FIG. 6.

Recombination rates (TRPL lifetime)−1, as a function of temperature for single-crystal CdTe (red circles), undoped CdS/CdTe (blue triangles), and doped CdS/CdTe (open black squares). The solid line shows an Arrhenius fit for the Cu-doped sample. Filled squares show recombination rates from TCAD simulations assuming NA = 1017 cm−3, ND = 9 × 1016 cm−3, electron capture cross-section 5 × 10−15 cm2, and temperature-independent electron mobility 50 cm2/(V s). The inset shows the recombination model that includes DAP tunneling and SRH recombination.

FIG. 6.

Recombination rates (TRPL lifetime)−1, as a function of temperature for single-crystal CdTe (red circles), undoped CdS/CdTe (blue triangles), and doped CdS/CdTe (open black squares). The solid line shows an Arrhenius fit for the Cu-doped sample. Filled squares show recombination rates from TCAD simulations assuming NA = 1017 cm−3, ND = 9 × 1016 cm−3, electron capture cross-section 5 × 10−15 cm2, and temperature-independent electron mobility 50 cm2/(V s). The inset shows the recombination model that includes DAP tunneling and SRH recombination.

Close modal

For the Cu-doped device, the lifetime at 170 K is 134 ns, which is on the order of single crystal lifetimes.20,24,42 An apparent Arrhenius trend with the activation energy EArrhenius = 86 ± 8 meV is observed in τ(T) data for the Cu-doped sample in Fig. 6. Arrhenius-like temperature dependence for capture cross-sections has been identified for non-radiative multiphonon transitions due to deep levels in some semiconductors, such as GaAs and GaP.43 Because the lifetime temperature dependence is very different for the single-crystal and Cu doped/undoped polycrystalline CdTe, it does not appear that a common deep defect can explain our experimental data, unless deep levels are created by Cu.

It is plausible that the lifetime decrease from 134 ns (170 K) to 11 ns (295 K) occurs due to relatively shallow recombination centers. The PL emission data indicate the presence of acceptors ≈150 meV above the valence band and shallow donors. Radiative lifetime at this doping is ≈10 μs,44 while experimental lifetimes range from 500 to 11 ns. Therefore, SRH recombination is dominant. For a single recombination center in p-type material, the Shockley-Read-Hall equation can be written as

(6)

where σn and vth,n are the electron capture cross-section and thermal velocity, respectively, kb is Boltzmann's constant, and T is temperature. The shift of the Fermi level Ef(T) with temperature will depend on the distribution of defects, and the PL spectra indicate the presence of both donors and acceptors. Time-resolved photoluminescence simulations were performed with code written in Sentaurus Device to analyze the role of compensation on lifetime temperature dependence.28,30 Without a large concentration of compensating donors, the Fermi level can move both above and below the Ea = 150 meV defect and produce trends that are not observed here. In the case of strong compensation, the Fermi energy temperature variation is reduced, and the CuCd recombination site can generate results consistent with the experimental data (see Fig. 6). Both donors and acceptors have been observed in the PL spectra here, and compensation has been widely reported in CdTe. So the PL and TRPL data can be consistently explained with the recombination model in Fig. 6. At low temperatures, recombination happens by DAP tunneling as described by Eq. (4). At higher temperatures, nonradiative SRH recombination from CuCd states describes the data.

The estimated capture cross-sections and acceptor densities needed to fit these data are on the order of 10−16 cm−2 and 1015−1017 cm−3, respectively, which is typical for neutral traps.45 Consequently, it is reasonable that the relatively shallow states can influence recombination. PL spectra for II–VI semiconductors are often dominated by emission from acceptors such as CuCd in CdTe, and our data indicate that these defects can limit lifetime.

This research was funded by NREL-First Solar Cooperative Research and Development Agreement (CRADA). The work at NREL was supported by the U.S. Department of Energy, Energy Efficiency and Renewable Energy, under Contract No. DE-AC36-08GO28308.

1.
Properties of Wide Bandgap II-VI Semiconductors
, edited by
R.
Bhargava
(
INSPEC
,
London
,
1997
).
2.
J. L.
Paurat
,
J. M.
Francou
,
N.
Magnea
,
E.
Molva
, and
K.
Saminadayar
,
J. Cryst. Growth
72
,
194
(
1985
).
3.
J.
Ma
,
J.
Yang
,
S. H.
Wei
, and
J. L. F.
Silva
,
Phys. Rev. B
90
,
155208
(
2014
).
4.
J. P.
Chamonal
,
E.
Molva
, and
J. L.
Pautrat
,
Solid State Commun.
43
,
801
(
1982
).
5.
H. H.
Woodbury
and
M.
Aven
,
J. Appl. Phys.
39
,
5485
(
1968
).
6.
D.
Grecu
and
A. D.
Compaan
,
Appl. Phys. Lett.
75
,
361
(
1999
).
7.
G.
Zoth
,
F. G.
Riedel
, and
W.
Schröter
,
Phys. Status Solidi B
172
,
187
(
1992
).
8.
See http://www.nrel.gov/ncpv/images/efficiency_chart.jpg for Best Research-Cell Efficiencies.
9.
M.
Woodhouse
,
K.
Horowitz
,
R.
Fu
,
D.
Chung
,
D.
Feldman
,
T.
Remo
, and
R.
Margolis
, “
On the path to sunshot: The role of advancements in photovoltaics efficiency, reliability, and costs
” (unpublished).
10.
L.
Kranz
,
C.
Gretener
,
J.
Perrenoud
,
D.
Jaeger
,
S. S. A.
Gerstl
,
R.
Schmitt
,
S.
Buecheler
, and
A. N.
Tiwari
,
Adv. Energy Mater.
4
,
1301400
(
2014
).
11.
C. R.
Corwine
,
A. O.
Pudov
,
M.
Gloeckler
,
S. H.
Demtsu
, and
J. R.
Sites
,
Sol. Energy Mater. Sol. Cells
82
,
481
(
2004
).
12.
S. H.
Demtsu
,
D. S.
Albin
,
J. R.
Sites
,
W. K.
Metzger
, and
A.
Duda
,
Thin Solid Films
516
,
2251
(
2007
).
13.
X.
Wu
,
J.
Zhou
,
A.
Duda
,
Y.
Yan
,
G.
Teeter
,
S.
Asher
,
W. K.
Metzger
,
S.
Demtsu
,
S. H.
Wei
, and
R.
Noufi
,
Thin Solid Films
515
,
5798
(
2007
).
14.
J.
Ma
,
S. H.
Wei
,
T. A.
Gessert
, and
K. K
Chin
,
Phys. Rev. B
83
,
245207
(
2011
).
15.
T. A.
Gessert
,
W. K.
Metzger
,
P.
Dippo
,
S. E.
Asher
,
R. G.
Dhere
, and
M. R.
Young
,
Thin Solid Films
517
,
2370
(
2009
).
16.
G.
Teeter
and
S.
Asher
, in
IEEE 33th Photovoltaic Specialist Conference, San Diego, CA
(
2008
), pp.
1
6
.
17.
J.
Perrenoud
,
L.
Kranz
,
C.
Gretener
,
F.
Pianezzi
,
S.
Nishiwaki
,
S.
Buecheler
, and
A. N.
Tiwari
,
J. Appl. Phys.
114
,
174505
(
2013
).
18.
S.
Vatavu
,
H.
Zhao
,
I.
Caraman
,
P.
Gaşin
, and
C.
Ferekides
,
Thin Solid Films
517
,
2195
(
2009
).
19.
D.
Kuciauskas
,
A.
Kanevce
,
J. N.
Duenow
,
P.
Dippo
,
M.
Young
,
J. V.
Li
,
D. H.
Levi
, and
T. A.
Gessert
,
Appl. Phys. Lett.
102
,
173902
(
2013
).
20.
J. N.
Duenow
,
J. M.
Burst
,
D. S.
Albin
,
D.
Kuciauskas
,
S. W.
Johnston
,
R. C.
Reedy
, and
W. K.
Metzger
,
Appl. Phys. Lett.
105
,
053903
(
2014
).
21.
M.
Gloeckler
,
I.
Sankin
, and
Z.
Zhao
,
IEEE J. Photovoltaics
3
,
1389
(
2013
).
22.
D.
Grecu
,
A. D.
Compaan
,
D.
Young
,
U.
Jayamaha
, and
D. H.
Rose
,
J. Appl. Phys.
88
,
2490
(
2000
).
23.
J. P.
Laurenti
,
G.
Bastide
,
M.
Rouzeyre
, and
R.
Triboulet
,
Solid State Commun.
67
,
1127
(
1988
).
24.
D.
Kuciauskas
,
A.
Kanevce
,
J. M.
Burst
,
J. D.
Duenow
,
R.
Dhere
,
D.
Albin
,
D. H.
Levi
, and
R.
Ahrenkiel
,
IEEE J. Photovoltaics
3
,
1319
(
2013
).
25.
D.
Kuciauskas
,
J. N.
Duenow
,
A.
Kanevce
,
J. V.
Li
,
M. R.
Young
,
P.
Dippo
, and
D. H.
Levi
, in
38th IEEE Photovoltaic Specialists Conference
(
2012
), pp.
1721
1726
.
26.
D. V.
O'Connor
and
D.
Phillips
,
Time Correlated Single Photon Counting
(
Academic Press
,
San Diego
,
1984
).
27.
W. K.
Metzger
,
R. K.
Ahrenkiel
,
J.
Dashdorj
, and
D. J.
Friedman
,
Phys. Rev. B
71
,
035301
(
2005
).
28.
A.
Kanevce
,
D.
Kuciauskas
,
T. A.
Gessert
,
D. H.
Levi
, and
D. S.
Albin
, in
38th IEEE Photovoltaic Specialists Conference, Austin, USA
(
2012
), pp.
848
853
.
29.
W. K.
Metzger
,
D.
Albin
,
D.
Levi
,
P.
Sheldon
,
X.
Li
,
B. M.
Keyes
, and
R. K.
Ahrenkiel
,
J. Appl. Phys.
94
,
3549
(
2003
).
30.
A.
Kanevce
,
D. H.
Levi
, and
D.
Kuciauskas
,
Prog. Photovoltaics: Res. Appl.
22
,
1138
(
2014
).
31.
D. E.
Aspnes
and
H.
Arwin
,
J. Vac. Sci. Technol., A
2
,
1309
1315
(
1984
).
32.
W.
Stadler
,
D. M.
Hofmann
,
H. C.
Alt
,
T.
Muschik
,
B. K.
Meyer
,
E.
Weigel
,
G.
Müller-Vogt
,
M.
Salk
,
E.
Rupp
, and
K. W.
Benz
,
Phys. Rev. B
51
,
10619
(
1995
).
33.
K.
Wei
,
F. H.
Pollak
,
J. L.
Freeouf
,
D.
Shvydka
, and
A. D.
Compaan
,
J. Appl. Phys.
85
,
7418
(
1999
).
34.
D. T. F.
Marple
,
Phys. Rev.
129
,
2466
(
1963
).
35.
C.
Kraft
,
H.
Metzner
,
M.
Hädrich
,
U.
Reislöhner
,
P.
Schley
,
G.
Gobsch
, and
R.
Goldhahn
,
J. Appl. Phys.
108
,
124503
(
2010
).
36.
D. S.
Albin
,
D.
Kuciauskas
,
J.
Ma
,
W. K.
Metzger
,
J. M.
Burst
,
H. R.
Moutinho
, and
P. C.
Dippo
,
Appl. Phys. Lett.
104
,
092109
(
2014
).
37.
E.
Molva
,
J. P.
Chamonal
,
G.
Milchberg
,
K.
Saminadayar
,
B.
Pajot
, and
G.
Neu
,
Solid State Commun.
44
,
351
(
1982
).
38.
E.
Zacks
and
A.
Halperin
,
Phys. Rev. B
6
,
3072
(
1972
).
39.
D. G.
Thomas
,
J. J.
Hopfield
, and
W. M.
Augustyniak
,
Phys. Rev.
140
,
A202
(
1965
).
40.
P. J.
Roland
,
N. R.
Paudel
,
Y.
Yan
, and
R. J.
Ellingson
, in
2014 IEEE 40th Photovoltaic Specialists Conference
(
IEEE
,
2014
), pp.
3266
3271
.
41.
D.
Kuciauskas
,
P.
Dippo
,
Z.
Zhao
,
L.
Cheng
,
A.
Kanevce
,
W. K.
Metzger
, and
M.
Gloeckler
, “
Recombination Analysis in Cadmium Telluride Photovoltaic Solar Cells With Photoluminescence Spectroscopy
,”
IEEE J. Photovoltaics
(to be published).
42.
D.
Kuciauskas
,
A.
Kanevce
,
P.
Dippo
,
S.
Seyedmohammadi
, and
R.
Malik
,
IEEE J. Photovoltaics
5
,
366
(
2015
).
43.
C. H.
Henry
and
D. V.
Lang
,
Phys. Rev. B
15
,
989
(
1977
).
44.
C. H.
Swartz
,
M.
Edirisooriya
,
E. G.
LeBlanc
,
O. C.
Noriega
,
P. A. R. D.
Jayathilaka
,
O. S.
Ogedengbe
,
B. L.
Hancock
,
M.
Holtz
,
T. H.
Myers
, and
K. N.
Zaunbrecher
,
Appl. Phys. Lett.
105
,
222107
(
2014
).
45.
Collaboration: Authors and editors of the volumes III/44B
,
CdTe: Impurity Complexes
, Springer Materials—The Landolt-Börnstein Database, edited by
O.
Madelung
,
U.
Rössler
, and
M.
Schulz
(
Springer
,
2008
).