The incorporation of electronic states in the bandgap, with high radiative efficiency, provides a means for optical up/down conversion and for tuning the carrier lifetime of a material. Such intermediate states in highly mismatched alloys have recently received much attention due to their potential for an efficient solar energy conversion. Understanding the carrier dynamics and pathways for the charge transfer of sub-bandgap transitions is critical for understanding energy conversion processes. In this work, time-resolved photoluminescence of ZnTeO is reported, revealing electron relaxation from the conduction band to the intermediate band and the carrier transfer between intermediate states and the conduction band utilizing two time-delayed optical excitations. This work demonstrates the utility of time-resolved techniques for characterizing energy conversion mechanisms in intermediate band materials, and the ability to use intermediate states to intentionally alter carrier lifetimes in materials for applications such as ultrafast scintillation.

The dilute incorporation of atoms with a large electronegativity difference, often referred to as highly mismatched alloys (HMAs), can result in a dramatic alteration of the electronic structure near the conduction band (CB) or the valence band (VB) edge. The modified band structure can lead to severe bandgap bowing or the formation of an intermediate band (IB) within the bandgap.1 The intermediate states will alter the carrier generation-recombination mechanisms, often undesirably in the case of nonradiative Shockley–Read–Hall recombination. However, the introduction of intermediate states with strong radiative strength can be highly desirable for energy conversion processes such as optical up- or down conversion or combined with an electronic response to realize lifetime-controlled fast detectors or intermediate band solar cells (IBSCs). The intermediate band solar cell (IBSC) concept utilizes sequential photogeneration across intermediate sub-bandgap states, providing a means to exceed the Shockley-Queisser efficiency limit by recovering losses due to bandgap transparency while maintaining an operating voltage that can track the bandgap of the material.2 The two-photon current generation process has been confirmed in several material systems in the steady state.3–7 

The success of HMAs for an efficient energy conversion requires that the dilute addition of impurities does not introduce nonradiative recombination pathways that counteract any advantage provided by sub-bandgap states. Increasing the density of states to induce an insulator-to-metal transition is predicted to result in “lifetime recovery” associated with delocalization of the electron wavefunction for impurity concentrations beyond the Mott transition.8 The concept of lifetime recovery has been debated, based on the argument that electron-phonon interactions will localize extended states, preventing lifetime recovery. In the absence of lifetime recovery, the efficiency may still be improved if the carrier lifetime in the conduction/valence band is much longer than the carrier transit time for collection.9 Carrier lifetime is clearly a critical parameter, along with carrier transport, that will determine the ultimate conversion efficiency.

There have been several studies reporting on carrier lifetimes in HMA systems: ZnTeO has shown promise for IBSCs, where the introduction of oxygen into ZnTe (EG = 2.3 eV at 300 K) results in electronic states approximately 0.4 eV below the conduction band edge.1,10,11 Intermediate band formation is proposed at sufficient oxygen introduction into ZnTe, described by a band anticrossing behavior.12,13 Time-resolved photoluminescence (TRPL) experiments have reported long lifetimes for carriers in the oxygen states of ZnTeO but with short lifetimes in the ZnTeO conduction band.11,14,15 Transient absorption spectroscopy on an analogous intermediate band material, GaPyAs1−x−yNx, similarly shows fast decay of carriers from the conduction band to the intermediate band with a time constant of approximately 23 ps.16 These prior experiments provide an initial indication of carrier recombination mechanisms but do not provide a full picture of carrier dynamics related to the sub-bandgap carrier generation in intermediate band materials. In this work, TRPL experiments are performed on ZnTeO using two sub-bandgap pump sources to directly probe the temporal response of the carrier transfer between the intermediate band and the conduction band, providing a direct confirmation of this process. Complementary transient absorption (TA) and steady state photoluminescence (PL) measurements indicate that stronger optical excitation can counteract some of the radiative recombination losses.

ZnTeO samples were grown by molecular beam epitaxy with solid sources for Zn, Te, and a plasma source for oxygen. Samples under study included growth on both GaAs (001) and GaSb (001) substrates. TRPL experiments were performed with a Ti-sapphire ultrafast laser system with two optical parametric amplifiers, allowing for wavelength tuning across the visible to near IR. The sample is first excited with a 400 nm (3.1 eV) pulse (VIS Pump 1) to promote electrons from VB to CB, followed by a second pulse centered at 1270 nm to provide sufficient energy for IB-CB absorption (absorption coefficient approximately 700 cm−1; Ref. 17), but less energy than required for VB-IB absorption, targeted at the CB-IB transition. For some experiments, involving excitations from the valence band to the intermediate band, a second optical pulse is applied (VIS Pump 2) with a wavelength of 590 nm (2.10 eV). A motorized delay stage spanning 100's of picoseconds controls the delay between the two pump pulses. A streak camera then records the TRPL signal. Transient absorption experiments are performed by an asynchronous optical sampling (ASOPS) system that offers four-color separation between the pump and the probe, a sensitivity of <10−7, and a dynamic range from femtosecond to nanosecond timescales,18 while preserving spatial overlap and spot size through this dynamic range. The energetic positions of the CB, IB, VB, and excitation sources are shown in Fig. 1.

FIG. 1.

Schematic illustration of energy levels in ZnTeO, where an IB due to OTe lies approximately 0.4 eV below the conduction band edge and 1.9 eV above the valence band edge. Also shown are various visible (VIS) and near-infrared (NIR) pumping schemes used in this study. VIS Pump 1 is 400 nm (3.1 eV); VIS Pump 2 is 590 nm (2.1 eV); and NIR is 1270 nm (0.98 eV).

FIG. 1.

Schematic illustration of energy levels in ZnTeO, where an IB due to OTe lies approximately 0.4 eV below the conduction band edge and 1.9 eV above the valence band edge. Also shown are various visible (VIS) and near-infrared (NIR) pumping schemes used in this study. VIS Pump 1 is 400 nm (3.1 eV); VIS Pump 2 is 590 nm (2.1 eV); and NIR is 1270 nm (0.98 eV).

Close modal

A comparison of the room-temperature steady-state photoluminescence (PL) spectra of ZnTe and ZnTe:O is illustrated in Fig. 2. A dramatic increase in the PL intensity is evident in the oxygen doped sample in the region of 1.8 eV. The steady-state PL intensity in the vicinity of this peak, widely held to be associated with the IB, is shown in Fig. 3(a) for three different excitation energies. No response is observed for excitation in the near-infrared (NIR) at 1270 nm, where absorption is not sufficient to promote carriers from VB to either IB or CB. A broad peak centered near 680 nm is observed for excitation by VIS Pump 2 attributed to emission from IB to VB. When the NIR source is added to the VIS Pump 2 VB-CB excitation, the PL intensity is reduced, indicating a reduction in the IB electron population. Optical generation of carriers from IB to CB may explain the reduced intensity and IB electron population. However, the steady state confirmation of PL from CB-VB as a result of the CB-IB carrier transfer is difficult in the case of short lifetimes for radiative recombination since the steady-state measurement duty cycle may dramatically reduce the observed PL amplitude. To circumvent this problem, we measured time-resolved PL (TRPL) spectra, which demonstrate emission from both CB-VB and IB-VB transitions, as shown in Fig. 3(b). The temporal resolution of these measurements using a streak camera is ∼30 ps. Next, we focus on the behavior of these two transitions under time-delayed excitation with a subsequent sub-bandgap optical pulse excitation.

FIG. 2.

Comparison of ambient temperature steady-state PL spectra from ZnTe (blue) and ZnTe : O (red).

FIG. 2.

Comparison of ambient temperature steady-state PL spectra from ZnTe (blue) and ZnTe : O (red).

Close modal
FIG. 3.

(a) Steady state photoluminescence spectra demonstrating IB-VB emission under visible pump 2 and depletion by introducing an additional near NIR pump, implying that IB electrons are radiatively excited into the CB. (b) TRPL spectra showing a fast decay for the CB-VB emission and a slow decay for the IB-VB emission.

FIG. 3.

(a) Steady state photoluminescence spectra demonstrating IB-VB emission under visible pump 2 and depletion by introducing an additional near NIR pump, implying that IB electrons are radiatively excited into the CB. (b) TRPL spectra showing a fast decay for the CB-VB emission and a slow decay for the IB-VB emission.

Close modal

The nature of the sub-bandgap emission, attributed to IB-VB transitions, can be further studied based on the spectral dependence of the TRPL quenching due to a subsequent NIR pump pulse. The spectral dependence of the quenching can help to determine the origin of the electronic states, and whether they are bandlike states related to oxygen, or discrete states related to native defects in ZnTe where emission at 1.8 eV has been reported.19 Following an initial VIS Pump 1, the PL spectra for time delays before and after a subsequent NIR pulse (time delay of 290 ps) are shown in Fig. 4. The PL spectra following the NIR pulse exhibits a clear quenching behavior on a faster timescale than the TRPL decay for IB-VB shown in Fig. 3(b), suggesting that quenching is attributable to the delayed NIR pulse. The difference between the TRPL spectra at the time instants before and after the NIR pulse is also shown in Fig. 4 to reveal the spectral shape of the PL quenching behavior. This illustrates how the PL spectrum associated with the IB evolves with a delay of the subgap NIR pulse, in particular, showing quenching of the TRPL signal. The line shape and the spectral position follow the emission spectrum of OTe substitutional impurities.10 Under the Franck-Condon description, the line shape and associated phonon replicas describe the overlap between the wavefunction of the ground vibrational state of the IB electronic state and the successive vibrational wavefunctions of the CB electronic state. As the IB vibrational wavefunction decreases its amplitude as a result of the PL depletion, the quenching is maximum at the peak and decreases at energies above or below the peak. Similarly, spectral ranges with stronger overlap are expected to exhibit stronger quenching, where the PL depletion is expected to qualitatively follow the steady state PL spectrum. This is also in agreement with the dominance of IB behavior, and in contrast to multiple discrete defect states that would show an energy dispersion difference in comparison to the original spectrum due to an interplay between the local density of states and their radiative recombination efficiency.

FIG. 4.

PL spectra recorded at two delay times following an initial VIS Pump 1 and a subsequent NIR pulse with a time delay of 290 ps. The spectra at the two time instants demonstrate quenching as a result of the NIR pulse. The blue line represents adjacent averaging of 5 time delay points binned together and represents the decrease of the IB emission after illumination by the NIR pulse. The line shape of the PL quench follows the emission spectrum of OTe, confirming that the oxygen-related energy level is involved in the sub-bandgap carrier generation.

FIG. 4.

PL spectra recorded at two delay times following an initial VIS Pump 1 and a subsequent NIR pulse with a time delay of 290 ps. The spectra at the two time instants demonstrate quenching as a result of the NIR pulse. The blue line represents adjacent averaging of 5 time delay points binned together and represents the decrease of the IB emission after illumination by the NIR pulse. The line shape of the PL quench follows the emission spectrum of OTe, confirming that the oxygen-related energy level is involved in the sub-bandgap carrier generation.

Close modal

While the broad spectral line shape in Fig. 3(a) is consistent with excitation from bandlike states, and spectral narrowing indicative of discrete states is not observed at low temperatures,19 steady-state PL data such as those shown in Figs. 2 and 3(a) are insufficient for a definitive identification of IB behavior. A deeper insight into the nature of these in-gap states can be obtained by performing time-resolved measurements that probe the transitions between individual states as depicted in Fig. 1 and their relaxation. TRPL experiments commonly involve initial photo-excitation of carriers into the conduction band using a single pump pulse. On the other hand, sequential two-photon excitation can provide a means to probe transitions that would occur in an operational IBSC. TRPL of the CB-VB transition is shown in Fig. 5(a) for a sequence of sub-bandgap illumination pulses. The sample is initially excited with below-bandgap visible light (VIS Pump 2), where a weak CB-VB PL is observed, likely due to nonlinear two-photon absorption present at high excitation intensities. After a time delay of ∼150 ps, an NIR pulse results in a significant enhancement in CB-VB emission. The enhanced emission after the time-delayed NIR pulse can again be explained by carrier excitation from IB to CB due to the NIR pulse. Carrier transfer from IB to VB is further observed by probing the TRPL for IB-VB as shown in Fig. 5(b). The data in Fig. 5(b) are for the case of a visible pump above the bandgap (VIS Pump 1 VB-CB) at time t = 0, followed by an NIR pulse at time t ∼ 300 ps. The initial visible pump results in a fast PL decay, corresponding to a combination of CB-VB PL emission and carrier transfer from CB to IB. PL emission for the IB-VB transition shows a response from the VIS Pump 1 VB-CB pump with a time constant of approximately 100 ps. At time t ∼ 300 ps, a rapid decrease of the IB-VB PL emission is observed, coinciding with a peak in the CB-VB emission. This behavior may be understood by the depletion of the IB population by the NIR pump resulting from the IB-CB transition. This demonstrates that the IB band can be filled and further depleted using a below-gap NIR pulse.

FIG. 5.

(a) Time-resolved photoluminescence of CB-VB emission due to sub-bandgap excitation (VIS pump 1) of the IB-VB transition followed by a delayed NIR (1270 nm) sub-bandgap pulse at approximately 150 ps later, demonstrating the carrier transfer to CB via the intermediate states; (b) time-resolved photoluminescence of CB-VB and IB-VB emission after initial VIS excitation (VIS pump 1) demonstrating a simultaneous enhancement of the CB PL and depletion of the IB PL upon introduction of a time-delayed NIR pump with approximately 300 ps delay.

FIG. 5.

(a) Time-resolved photoluminescence of CB-VB emission due to sub-bandgap excitation (VIS pump 1) of the IB-VB transition followed by a delayed NIR (1270 nm) sub-bandgap pulse at approximately 150 ps later, demonstrating the carrier transfer to CB via the intermediate states; (b) time-resolved photoluminescence of CB-VB and IB-VB emission after initial VIS excitation (VIS pump 1) demonstrating a simultaneous enhancement of the CB PL and depletion of the IB PL upon introduction of a time-delayed NIR pump with approximately 300 ps delay.

Close modal

Carrier transitions at faster time resolution were investigated using the ASOPS technique to investigate the femtosecond TA dynamics of this system. Figure 6 shows TA data for ZnTe and ZnTeO conduction band carriers. A VIS Pump 1 VB-CB pulse injects electrons into the conduction band, followed by a time-delayed IR probe (1550 nm). The transient electron population in the conduction band may be inferred from the transmittance of the IR probe, as shown in Fig. 6 for ZnTe and ZnTeO. The changes in transmittance T of the IR probe are quantified by the ratio dT/T with absorption indicated by a decrease of this ratio. At early time delays, a much shorter (few picoseconds) lifetime is observed for ZnTeO in comparison to ZnTe (>100 ps), indicating that the introduction of sub-bandgap electronic states related to oxygen results in ultrafast carrier relaxation. At longer time delays (>100 ps), the carrier relaxation in ZnTeO slows down in comparison to ZnTe and is attributed to the carriers present in the IB of ZnTeO. A simple exponential fit was performed on the relaxation data in two time regimes, short and long (>100 ps); for example, fitting the short-time regime data in Fig. 7 we see that the relaxation time follows an increasing trend with increasing illumination power in the range of 3–10 μJ/cm2 and corresponding relaxation times increasing from 5–10 ps. Moreover, at longer times (>100 ps) where free-carrier absorption is likely to dominate, ZnTeO shows a marked increase in relaxation times (few hundred picoseconds) compared to the undoped ZnTe samples. This establishes a correlation between the free carrier absorption and the TA relaxation. It should be noted that the sub-bandgap states in ZnTeO are initially unpopulated prior to optical excitation, due to the background p-type nature of the material.20 Altering the electron population in the intermediate states would, therefore, change the carrier relaxation time. Doping, electrical injection, or optical excitation may be used to increase electron population in the intermediate states. A portion of the initial pump will inject some fraction of charge carriers into the intermediate states and provides a means to alter the intermediate state occupation. Transient absorption for variable pump power is shown in Fig. 7(a), where an increase in time constant is observed ranging from approximately 3.1 ps to 8.3 ps. This carrier lifetime increase is attributed to an increased occupation of intermediate states and a corresponding suppression of carrier relaxation from the conduction band to intermediate states due to IB filling that can reduce the CB-IB recombination. Similarly, an increased radiative efficiency is observed in the steady state due to carrier filling, as shown in Fig. 7(b). The PL peak position related to intermediate state emission demonstrates a Burstein-Moss effect blue shift with increasing pump power. While these experiments do not investigate a large range of carrier occupation or provide extensive quantitative information, there is a demonstrable dependence of a conduction band carrier lifetime on occupation. Although not fully quantified here, capturing the multiband charge transfer and recombination dynamics directly using our TRPL and TA approaches is a substantial step forward in enabling future studies of multiband semiconductors that undergo nonlinear carrier dynamics largely dependent on carrier densities and band filling effects, to assist device prospects such as solar cells.21 Using feedback from multipump excitation, TRPL, and TA, the fast relaxation pathway associated with oxygen incorporation may, therefore, be tuned by doping, electrical injection, or selective optical excitation.

FIG. 6.

Transient absorption for VIS excitation above bandgap, followed by an NIR probe of the CB population, where a substantial reduction in the CB carrier relaxation lifetime is observed at early time delays for ZnTeO in comparison to ZnTe [varying time segments for (a) and (b)]. Conversely, an increase in the carrier lifetime is observed in ZnTeO on longer timescales (>100 ps), reflecting the important role of free carrier CB absorption.

FIG. 6.

Transient absorption for VIS excitation above bandgap, followed by an NIR probe of the CB population, where a substantial reduction in the CB carrier relaxation lifetime is observed at early time delays for ZnTeO in comparison to ZnTe [varying time segments for (a) and (b)]. Conversely, an increase in the carrier lifetime is observed in ZnTeO on longer timescales (>100 ps), reflecting the important role of free carrier CB absorption.

Close modal
FIG. 7.

Excitation power dependent (a) transient absorption and corresponding single exponential fit demonstrating an increase in decay time constant and (b) TRPL demonstrating an increase in PL intensity around 1.9 eV as well as a blue shift with increasing carrier occupation. Illumination power increases over the range 3–10 μJ/cm2.

FIG. 7.

Excitation power dependent (a) transient absorption and corresponding single exponential fit demonstrating an increase in decay time constant and (b) TRPL demonstrating an increase in PL intensity around 1.9 eV as well as a blue shift with increasing carrier occupation. Illumination power increases over the range 3–10 μJ/cm2.

Close modal

We investigated the carrier dynamics of sub-bandgap electronic states induced by oxygen incorporation in ZnTe with time-resolved photoluminescence and transient absorption measurements. From these measurements, it is clear that systematic time-resolved studies with sub-picosecond resolution provide deep insights into the charge transfer mechanisms. Moreover, the results serve to reveal the bandlike characteristics of the in-gap states and to distinguish them from discrete impurity states. The time-resolved response of carrier transfer was explored with high temporal resolution for the first time between IB and CB for both photogeneration from IB to CB and relaxation from IB to CB. These transitions demonstrate the critical mechanism for realizing optical up conversion and an intermediate band solar cell. Oxygen incorporation results in a fast decay from CB to IB, which may be altered by varying the carrier population in the IB. The fast carrier relaxation from CB to IB reveals a weakness for the intermediate band solar energy conversion for an empty IB in equilibrium and underscores the importance of controlling the carrier occupation in the IB through doping and/or a device design.21–23 The tunable carrier lifetime and ultrafast carrier relaxation time with doping may also provide a new degree of freedom for achieving devices with a high radiative efficiency and a fast time-response, such as an ultrafast scintillator.

This work was supported in part by the Center for Solar and Thermal Energy Conversion and an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences under Award No. DE-SC0000957. This work was performed, in part, at the Center for Nanoscale Materials, a U.S. Department of Energy Office of Science User Facility and supported by the U.S. Department of Energy, Office of Science, under Contract No. DE-AC02-06CH11357.

1.
J.
Wu
,
W.
Walukiewicz
, and
E. E.
Haller
,
Phys. Rev. B
65
(
23
),
233210
(
2002
).
2.
A.
Luque
and
A.
Martí
,
Phys. Rev. Lett.
78
(
26
),
5014
5017
(
1997
).
3.
A.
Martí
,
E.
Antolín
,
E.
Cánovas
,
N.
López
,
P. G.
Linares
,
A.
Luque
,
C. R.
Stanley
, and
C. D.
Farmer
,
Thin Solid Films
516
(
20
),
6716
6722
(
2008
).
4.
Y.
Okada
,
T.
Morioka
,
K.
Yoshida
,
R.
Oshima
,
Y.
Shoji
,
T.
Inoue
, and
T.
Kita
,
J. Appl. Phys.
109
(
2
),
024301
(
2011
).
5.
A.
Scaccabarozzi
,
S.
Adorno
,
S.
Bietti
,
M.
Acciarri
, and
S.
Sanguinetti
,
Phys. Status Solidi Rapid Res. Lett.
7
(
3
),
173
176
(
2013
).
6.
N.
Ahsan
,
N.
Miyashita
,
M. M.
Islam
,
K. M.
Yu
,
W.
Walukiewicz
, and
Y.
Okada
,
Appl. Phys. Lett.
100
(
17
),
172111
(
2012
).
7.
H.
Jussila
,
P.
Kivisaari
,
J.
Lemettinen
,
T.
Tanaka
, and
M.
Sopanen
,
Phys. Rev. Appl.
3
(
5
),
054007
(
2015
).
8.
A.
Luque
,
A.
Martí
,
E.
Antolín
, and
C.
Tablero
,
Physica B
382
(
1–2
),
320
327
(
2006
).
9.
J. J.
Krich
,
B. I.
Halperin
, and
A.
Aspuru-Guzik
,
J. Appl. Phys.
112
(
1
),
013707
(
2012
).
10.
J. L.
Merz
,
Phys. Rev.
176
(
3
),
961
968
(
1968
).
11.
W.
Wang
,
A.
Lin
,
J. D.
Phillips
, and
W.
Metzger
,
Appl. Phys. Lett.
95
,
261107
(
2009
).
12.
T.
Tanaka
,
T.
Mochinaga
,
K.
Saito
,
Q.
Guo
,
M.
Nishio
,
K. M.
Yu
, and
W.
Walukiewicz
, paper presented at the 2012 38th IEEE Photovoltaic Specialists Conference, Austin, Texas, June 3–8, 2012.
13.
T.
Tanaka
,
K.
Matsuo
,
K.
Saito
,
Q.
Guo
,
T.
Tayagaki
,
K. M.
Yu
, and
W.
Walukiewicz
,
J. Appl. Phys.
125
(
24
),
243109
(
2019
).
14.
Y.-C.
Lin
,
M.-J.
Tasi
,
W.-C.
Chou
,
W.-H.
Chang
,
W.-K.
Chen
,
T.
Tanaka
,
Q.
Guo
, and
M.
Nishio
,
Appl. Phys. Lett.
103
(
26
),
261905
(
2013
).
15.
J.
Li
,
J.
Ye
,
F.
Ren
,
D.
Tang
,
Y.
Yang
,
K.
Tang
,
S.
Gu
,
R.
Zhang
, and
Y.
Zheng
,
Sci. Rep.
7
,
44399
(
2017
).
16.
J. N.
Heyman
,
A. M.
Schwartzberg
,
K. M.
Yu
,
A. V.
Luce
,
O. D.
Dubon
,
Y. J.
Kuang
,
C. W.
Tu
, and
W.
Walukiewicz
,
Phys. Rev. Appl.
7
(
1
),
014016
(
2017
).
17.
E.
Antolín
,
C.
Chen
,
I.
Ramiro
,
J.
Foley
,
E.
López
,
I.
Artacho
,
J.
Hwang
,
A.
Teran
,
E.
Hernández
,
C.
Tablero
,
A.
Martí
,
J. D.
Phillips
, and
A.
Luque
,
IEEE J. Photovoltaics
4
(
4
),
1091
1094
(
2014
).
18.
V. A.
Stoica
,
Y.-M.
Sheu
,
D. A.
Reis
, and
R.
Clarke
,
Opt. Express
16
(
4
),
2322
2335
(
2008
).
19.
C.
Chen
,
J.
Zheng
,
K.
Nguy
,
F.
Naab
, and
J.
Phillips
,
J. Electron. Mater.
43
(
4
),
879
883
(
2014
).
20.
J. D.
Dow
,
R.-D.
Hong
,
S.
Klemm
,
S. Y.
Ren
,
M. H.
Tsai
,
O. F.
Sankey
, and
R. V.
Kasowski
,
Phys. Rev. B
43
(
5
),
4396
4407
(
1991
).
21.
A. S.
Lin
,
W.
Wang
, and
J. D.
Phillips
,
J. Appl. Phys.
105
,
064512
(
2009
).
22.
A.
Lin
and
J.
Phillips
,
IEEE Trans. Electron Devices
56
(
12
),
3168
3174
(
2009
).
23.
A.
Lin
and
J.
Phillips
,
Prog. Photovolt.
22
(
10
),
1062
1069
(
2014
).