Gang
Qiang‡
*a,
Aleksandr A.
Golovatenko‡
*b,
Elena V.
Shornikova
a,
Dmitri R.
Yakovlev
*ab,
Anna V.
Rodina
b,
Evgeny A.
Zhukov
ab,
Ina V.
Kalitukha
b,
Victor F.
Sapega
b,
Vadim Kh.
Kaibyshev
b,
Mikhail A.
Prosnikov
c,
Peter C. M.
Christianen
c,
Aleksei A.
Onushchenko
d and
Manfred
Bayer
ab
aExperimentelle Physik 2, Technische Universität Dortmund, 44227 Dortmund, Germany. E-mail: gang.qiang@tu-dortmund.de; dmitri.yakovlev@tu-dortmund.de
bIoffe Institute, Russian Academy of Sciences, 194021 St. Petersburg, Russia. E-mail: sasha.pti@mail.ioffe.ru
cHigh Field Magnet Laboratory (HFML-EMFL), Radboud University, 6525 ED Nijmegen, The Netherlands
dITMO University, 199034 St Petersburg, Russia
First published on 7th December 2020
The recombination dynamics and spin polarization of excitons in CdSe nanocrystals synthesized in a glass matrix are investigated using polarized photoluminescence in high magnetic fields up to 30 Tesla. The dynamics are accelerated by increasing temperature and magnetic field, confirming the dark exciton nature of low-temperature photoluminescence (PL). The circularly polarized PL in magnetic fields reveals several unusual appearances: (i) a spectral dependence of the polarization degree, (ii) its low saturation value, and (iii) a stronger intensity of the Zeeman component which is higher in energy. The latter feature is the most surprising being in contradiction with the thermal population of the exciton spin sublevels. The same contradiction was previously observed in the ensemble of wet-chemically synthesized CdSe nanocrystals but was not understood. We present a theory which explains all the observed features and shows that the inverted ordering of the circularly polarized PL maxima from the ensemble of nanocrystals is a result of competition between the zero phonon (ZPL) and one optical phonon-assisted (1PL) emission of the dark excitons. The essential aspects of the theoretical model are different polarization properties of the dark exciton emission via ZPL and 1PL recombination channels and the inhomogeneous broadening of the PL spectrum from the ensemble of nanocrystals exceeding the optical phonon energy.
For direct bandgap semiconductor NCs, e.g. made up of CdSe, CdTe or InP, most of the optical properties can be explained within the exciton fine structure model.17–20 In ideal spherical NCs, the band-edge exciton state 1S3/21Se is eight-fold degenerate. 1S3/2 is the lowest quantum size level of the holes and is four-fold degenerate with respect to its total angular momentum components. 1Se is the lowest quantum size level of the electrons, which is doubly degenerate with respect to its spin components. However, due to the NC shape asymmetry, the intrinsic crystal field, and the electron–hole exchange interaction, the degeneracy is lifted resulting in the formation of five exciton states. In the case of wurtzite CdSe NCs, the two lowest states are a typically optically forbidden (dark) exciton state |F〉 (angular momentum projection along the NC quantization axis ±2) and a higher-lying optically allowed (bright) exciton state |A〉 (momentum projection ±1).19 The presence of the bright and dark excitons allows one to explain the recombination dynamics, spin dynamics, and magneto-optical properties, taking into account also the Zeeman effect.17,21–27
CdSe NCs are a test bed for investigation of colloidal NCs; therefore, their optical properties have been extensively studied. External magnetic fields have been used as a powerful tool not only to address magneto-optical properties and spin-dependent phenomena, but also to identify the exciton energy, spin structure and mechanisms responsible for recombination of the dark excitons. In particular, the following magneto-optical properties in external magnetic fields of wet-chemically synthesized colloidal CdSe NCs were reported: shortening of the dark exciton lifetime,17,28 circular polarization of photoluminescence (PL),21,22,29,30 exciton fine structure splitting including the Zeeman effect in single NCs,25,31 anisotropic exchange interaction,32,33 and electron spin coherence.34,35 All these make the colloidal NCs promising for spintronics and quantum information applications based on the spin-dependent phenomena. No magneto-optical study has been performed so far on CdSe in the glass matrix.
Studies on the circularly polarized PL of CdSe-based NCs in the magnetic field22,29 revealed a puzzling behavior: with increasing field the σ− polarized PL, which is stronger in intensity compared to the σ+ component, shifts to higher energy than the σ+ polarized emission. This contradicts with the expectation for the thermal population of the exciton spin sublevels split in the magnetic field. Later, the same contradiction was reported for colloidal CdTe NCs.24 This contradiction, however, is absent in fluorescence line narrowing (FLN) experiments and in measurements of single NCs, i.e. for the experimental conditions when the inhomogeneous broadening of the PL due to NC size dispersion in the ensemble is suppressed.25,29 In the FLN experiments, pronounced optical phonon-assisted recombination of the dark exciton was observed.17,28,29,36–38 Under nonresonant excitation, the low-temperature PL of the NC ensemble is contributed by zero phonon and optical phonon-assisted exciton emission, which form overlapping bands due to NCs of different sizes.17,37 In ref. 22 and 29, the contradiction with respect to the energy shift was tentatively assigned to the interplay between different emission channels of the dark exciton. However, a corresponding theoretical model has not been proposed so far.
For the present magneto-optical study, we choose CdSe NCs which are synthesized in a glass matrix. One of the motivations is to compare their magneto-optical properties with those of the wet-chemically synthesized CdSe NCs. Further reasons are based on the simpler experimental situation, namely the simpler surface conditions due to the absence of organic ligand passivation, long-term stability due to the matrix encapsulation, and the ideally random orientation of the NC quantization axis in the glass-based NC ensemble which allow us to make an easier model description. The wet-chemically grown and glass-based NCs differ by their surface conditions. In the first case, a NC is passivated using organic ligands, in the second case is surrounded by a glass matrix. Our study shows that the magneto-optical properties of the dark excitons in both cases are similar and not affected by the difference in the NC surrounding media.
In this paper, we study CdSe NCs in glass with diameters varying from 3.3 nm up to 6.1 nm and investigate the recombination dynamics and spin polarization properties of excitons in strong magnetic fields up to 30 T. We find experimentally the same contradiction between intensities and energies of the circularly polarized PL as reported for wet-chemically prepared CdSe NCs, and a strong spectral dependence of the circular polarization degree. We develop a theoretical model accounting for the contributions of both zero phonon emission (ZPL) and first optical phonon-assisted emission (1PL) of the dark excitons, which allows us to describe all experimental signatures.
Fig. 1 Photoluminescence and recombination dynamics. (a) PL spectra of the studied CdSe NCs measured at T = 4.2 K. (b) Recombination dynamics measured for sample D6.1 at various temperatures, B = 0 T. (c) Recombination dynamics measured for sample D6.1 at B = 0 and 17 T, T = 2.2 K. (d) Scheme of exciton energy levels. |A〉, |F〉, and |G〉 denote the bright and dark exciton states, and crystal ground state, respectively. ΔEAF is the bright–dark splitting. γ0 is the zero-temperature relaxation rate from the bright to the dark exciton state. γ0NB is the thermal activation rate of the reverse process. ΓA and ΓF are the recombination rates of the bright and dark excitons. (e) Temperature dependence of the long component of the PL decay, τLong, for all samples at B = 0 T. The lines are fits with eqn (S5†). (f) Magnetic field dependence of τLong at T = 4.2 K for all samples. The lines show the averaged lifetime of the dark exciton calculated with eqn (S23†) in a randomly oriented ensemble of NCs as described in ESI section S4.2.† The color codes in panels (e and f) are same as in panel (a). |
Sample | D3.3 | D4.1 | D4.9 | D6.1 |
---|---|---|---|---|
NC diameter (nm) | 3.30 ± 0.17 | 4.10 ± 0.21 | 4.9 ± 0.25 | 6.1 ± 0.31 |
PL FWHM (meV) | 120 | 99 | 96 | 101 |
PL peak energy (eV) | 2.366 | 2.210 | 2.111 | 2.017 |
Absorption peak energy (eV) | 2.455 | 2.267 | 2.160 | 2.039 |
Stokes shift (meV) | 89 | 57 | 49 | 22 |
Bright–dark splitting from PL decay (meV) | 3.9 | 4.3 | 2.6 | 0.9 |
Bright–dark splitting from FLN (meV) | — | 8.4 | 7.5 | 3.5 |
The low-temperature PL in the studied samples is dominated by exciton recombination, which has a characteristic decay variation with temperature and magnetic field. This allows one to distinguish neutral and charged excitons in the NCs.40–42 Taking sample D6.1 as a representative example, we show in Fig. 1b the recombination dynamics at various temperatures. With increasing temperature from 2.2 to 20 K, the decay time of the long component, τLong, shortens from 212 ns down to 58 ns. The temperature dependence of τLong for all the studied samples is given in Fig. 1e. One can see that the decay times shorten with increasing temperature as a result of thermal activation from the dark exciton to the bright exciton state. A similar behavior is observed when external magnetic field is applied; in Fig. 1c, the long component measured at T = 2.2 K shortens from 212 ns down to 119 ns with increasing field up to B = 17 T. The magnetic field dependence of τLong for all studied samples is given in Fig. 1f.
This is the typical behavior for exciton emission in CdSe NCs and can be well explained with the diagram as shown in Fig. 1d (a detailed scheme of the exciton levels is given in Fig. S13a†). The optically forbidden dark state |F〉 with angular momentum projection ±2 is the exciton ground state. The closest optically allowed bright state |A〉 with angular momentum projection ±1 is shifted to higher energy by the bright–dark splitting ΔEAF.17 The fast component of the PL decay is provided by the recombination of bright excitons with rate ΓA and their rapid thermalization to the dark state with rate γ0(1 + NB), where γ0 is the spin–flip rate, NB = 1/[exp(ΔEAF/kBT) − 1] is the Bose–Einstein phonon occupation at temperature T. The long component corresponds to the dark exciton recombination with rate ΓF. With an increase in temperature, the dark excitons are activated to the bright state. In the case of a one-phonon process, the activation rate is given by γ0NB. The shortening of the long component in a magnetic field is due to the mixing of bright and dark exciton states by the field component perpendicular to the NC quantization axis.17,40 These properties allow us to uniquely assign the long decay component τLong to the recombination of the dark excitons. More information about the PL decay fitting is given in ESI section S2.†
The analysis of the τLong temperature dependence allows us to evaluate the bright–dark splitting ΔEAF (see ESI section S2†). The corresponding fits are shown by the solid lines in Fig. 1e, and the evaluated ΔEAF values are given in Table 1. In Fig. S6c,† we compare the ΔEAF values with the literature data17,37–39 and a good consistency is found. However, ΔEAF in D3.3 is deviating from the general trend, which may point toward a slightly prolate shape of the CdSe NCs.43 The ΔEAF values determined from fluorescence line narrowing (see ESI section S2.3†) are found to be larger compared to the values determined from the τLong temperature dependence. Two explanations for this discrepancy were proposed: (i) The ΔEAF values determined from the τLong temperature dependence correspond to the true values of the bright–dark splitting, while the ΔEAF values from FLN can show an additional energy shift due to polaron formation,44 (ii) The ΔEAF values determined from the τLong temperature dependence correspond to the energy of confined acoustic phonons, which provide thermal activation to the bright excitons.45
Fig. 2 Polarized photoluminescence in magnetic field. (a) (left axis) σ− (blue) and σ+ (red) circularly polarized emission, measured at B = 30 T for sample D6.1. The black spectrum at B = 0 T is shown for comparison. (right axis) Green and red-dashed curves show the experimental and calculated spectral dependence of the time-integrated DCP, respectively. The orange curve is the function fZPL(E) and the cyan curve is the calculated PL spectrum at B = 0 T, which accounts for the ZPL and 1PL contributions. (b) Magnetic field dependences of Pintc (B) measured at the PL maximum in all samples. Lines are fits with eqn (1) and (3). |
Fig. 2b shows the magnetic field dependence of the time-integrated degree of circular polarization (DCP), Pintc, measured at the PL maximum. It is expected that in the ensemble of NCs with randomly oriented c-axis (the case for NCs in glass), the DCP saturation value should reach −0.75.21 However, as one can see in Fig. 2b, this is obviously not the case, Pintc does not reach −0.75 and has a clear size-dependence. At B = 30 T, Pintc = −0.56 in D6.1 and −0.66 in D3.3.
Fig. 3a–d show the magnetic field dependence of the polarized PL intensity. The σ− polarized component increases monotonously with the magnetic field and the σ+ polarized component remains almost constant after an initial decrease. Most surprising are the spectral shifts of the σ+ and σ− polarized components (Fig. 3e and f). For neutral excitons, it is expected that the PL maximum of the stronger component, i.e. σ−, shifts in the magnetic field to lower energy, as the lower energy spin state of the exciton has the higher thermal population. This expectation is shown by the dashed lines for the sample D6.1 in Fig. 3e (more details will be given below), but the experimental shift of the PL maximum demonstrates a very different behavior. First, for the D6.1, D4.9 and D4.1 samples, the σ+ component shifts to lower energy and second, the splitting between σ+ and σ− is several times larger than expected (Fig. 3e–g). In the sample D3.3, both components shift initially to higher energy and the σ+ component lowers its energy for fields above 15 T (Fig. 3h).
Fig. 3 Photoluminescence intensity and spectral shifts in magnetic fields. (a–d) PL intensity of the σ+ (red) and σ− (blue) polarized PL as functions of the magnetic field in the CdSe NCs. (e–h) Magnetic field dependences of the corresponding PL peak energies. For all panels, the symbols correspond to the experimental data, and the lines show the modeling results. Dashed lines in panel (e) show the Zeeman splitting of the dark exciton spin sublevels −2 (blue) and +2 (red) in a nanocrystal with the c-axis parallel to the magnetic field direction. We note that the energy of ±2 states at B = 0 T is taken the same as the PL maxima position for the sake of demonstrativeness. Relative positions of ±2 energy state distribution maxima and polarized PL maxima in the applied magnetic field are shown in Fig. S15.† |
The modeling results presented in Fig. 2 and 3 are based on the consideration of the ZPL and one optical phonon-assisted line (1PL) contributions to the PL. The contribution of the linearly polarized recombination of the dark exciton with assistance of acoustic phonons is considered in ESI section S4.† Our task is to calculate the σ+ and σ− polarized PL spectra and the spectral dependence of Pintc. For this purpose, we sum the ZPL and 1PL contributions, which have different spectral distributions, and average them over the random orientation of the hexagonal c-axis in an ensemble of NCs:
(1) |
Here, E is the spectral energy, δE±2(x, B) = ±gFμBBx/2 are the Zeeman shifts of the ±2 dark exciton states, x = cosθ with θ being the angle between the c-axis of the nanocrystal and the magnetic field applied in the Faraday geometry. The explicit forms of Ii,ZPL±(x) and Ii,1PL±(x) are given in ESI section S4.† Functions fZPL(E) and f1PL(E) describing the inhomogeneous broadening of ZPL and 1PL emission due to size dispersion of NCs are determined as follows:
(2) |
Here, E0ZPL corresponds to the energy of the ZPL in QDs at the maximum of the size distribution, ELO = 26 meV is the energy of the optical phonon in CdSe, and w is the standard deviation. The function fZPL(E) used for the fitting of the experimental data for the D6.1 sample is shown in Fig. 2a by the orange line. The cyan line in Fig. 2a models the PL spectrum at B = 0 T with accounting for the 1PL contribution.
Using eqn (1), we calculate the DCP spectral dependence as:
(3) |
The fit parameters are: the dark exciton g-factor gF, the characteristic energy of interaction ε which results in the admixture of the dark and ±1L bright exciton states in the zero magnetic field, the ratio of the 1PL and ZPL recombination rates at zero magnetic field χ0, and a phenomenological parameter c1PL, which determines an increase in the 1PL recombination rate in the magnetic field, while an increase in the ZPL recombination rate in the magnetic field is described within second order perturbation theory.17,48 The parameter χ0 can be estimated independently from the ratio of the 1PL to ZPL intensities measured in fluorescence line narrowing experiments. Typically, this ratio is of the order of unity.17,28,29,36–38 The FLN spectra for samples D4.1, D4.9 and D6.1 are shown in Fig. S7.† Fitting of the magnetic field dependence of the long decay time (Fig. 1f) imposes a restriction on the χ0, ε and c1PL values. Finally, the linked sets of parameters χ0, ε and c1PL together with gF are used for joint fitting of the magnetic field dependences of Pintc, the intensities and maxima positions of the σ− and σ+ polarized PL.
Besides the saturation of the DCP in high magnetic fields, the 1PL emission causes the DCP spectral dependence (calculation is shown by the orange dashed curve in Fig. 2a). The reason is the weaker relative intensity of the 1PL emission at the high-energy side of the PL spectrum due to the small number of NCs that can provide this contribution. The pronounced spectral dependence of the DCP in large NCs correlates with its low saturation level, i.e. with the strong 1PL emission (Fig. S18†).
We turn now to the modeling of the intensities and maxima positions of the σ+ and σ− polarized PL. The calculated magnetic field dependences of the PL intensities (Fig. 3a–d) are in good agreement with the experimental data for all the studied samples. In a randomly oriented ensemble of NCs, the emission of σ+ polarized photons from the predominantly populated −2 dark exciton state results in the saturation of the σ+ PL intensity at a level of 0.4 of the zero-field intensity. This is observed for all the studied samples (Fig. 3a–d).
The ability of the dark exciton states to emit both σ− and σ+ polarized light also explains the energy shifts of the polarized PL maxima in a magnetic field. The spectral maximum of σ− polarized PL (solid line in Fig. 4e) is determined by the relationship between the σ− polarized emission from the −2 (dashed line) and +2 (dash-dotted line) states. In a magnetic field, the ZPL emission from the −2 state dominates, and the PL maximum shifts toward higher energy. This behavior is observed in all the studied samples (Fig. 3e and h). Similarly (Fig. 4f), the spectral maximum of σ+ polarized PL is determined by the emission from the −2 and +2 states. In high magnetic fields, the 1PL emission from the −2 state dominates, and the PL maximum shifts to lower energy. Thus, positions of polarized PL maxima to a large extent are controlled by the emission of the −2 state: ZPL in the case of σ−, and 1PL in the case of σ+.
For comparison, the Zeeman splitting of the dark exciton spin sublevels −2 (blue) and +2 (red) in a nanocrystal with the c-axis parallel to the magnetic field direction is presented in Fig. 3e by the dash lines. One can see that it is much smaller than the difference of PL maxima energies and also results in reversed order (see also Fig. S15†).
The developed model allows us to describe all unusual features of the polarized PL in CdSe NCs by accounting for the superposition of the ZPL and 1PL emission of the dark excitons. It is important to note that the model is universal and can be readily applied to describe the properties of the wet-chemically grown NCs. Their magneto-optical properties are very similar to what we have found in CdSe NCs in glass.22,29 This allows us to conclude that the spectral shifts of the polarized emission in magnetic fields are controlled by the mechanism suggested in this paper and are weakly dependent on the surface conditions. However, the surface can be strongly modified by the growth conditions and post-growth treatment in the wet-chemically grown NCs. Furthermore, in this case, the splitting of the PL lines in magnetic fields can be additionally contributed by the surface states.
The dark exciton g-factor comprises the electron, ge, and hole, gh, g-factors as17gF = ge − 3gh. The size dependence of ge is well known for CdSe NCs34,35,51 (Fig. S26†). This allows us to evaluate the hole g-factor as gh = (ge − gF)/3. Using the gF values determined from the fitting with the linearly polarized contribution to the ZPL, we find gh ≈ −0.45 and without it gh ≈ −0.15 for NCs with diameters of 4–6 nm (Fig. S27†). It is smaller than the theoretically estimated value of gh ≈ −1 from ref. 17 and 52. However, the first value is close to gh ≈ −0.56 measured in giant shell CdSe/CdS colloidal nanocrystals from the polarized PL of negative trions.40
The values of the dark exciton g-factor used for the modeling of experimental data in Fig. 2 and 3 are obtained in the assumption of thermal equilibrium between the ±2 states. However, the spin relaxation between these states requires a change in the exciton total spin projection by ΔF = 4. If the rate of this process is much slower than the dark exciton lifetime, the populations of the ±2 states are determined by the relaxation from the bright ±1 exciton states. As shown in ESI section S4.7–8,† the fitting of the experimental dependences in this case results in larger values of the dark exciton g-factor, which are close to the theoretical estimations. The spin relaxation time, or spin dynamics, is often evaluated from the rise of the degree of circular polarization in the magnetic field.21,22 However, such evaluation might be complicated in the case of an ensemble with a nonmonoexponential PL decay.53 We have presented the analysis of the DCP dynamics in ESI section S3,† showing that in any case it is faster than the dark exciton recombination and there is not much of a difference between the time-integrated and saturated DCP values. However, as this time is of the order of the relaxation time from the bright to the dark exciton state (the first time component in the PL decay), it is not possible to conclude whether the DCP rise and the population of the dark exciton sublevels is determined by the spin relaxation between ±1 and ±2 states. Further clarification of this question is beyond the scope of this paper.
In previous publications, a saturation value of the DCP in high magnetic fields lower than −0.75 was associated with nonradiative recombination of the dark excitons,21 partial activation of the dark exciton recombination via an admixture of the 0U bright exciton resulting in linearly polarized ZPL emission,54 or the presence of prolate and oblate nanocrystals in the same ensemble.24 Accounting for the linearly polarized dark exciton emission assisted by optical phonons is in line with the first two explanations. Here, we show that consideration of the 1PL emission of the dark excitons besides the saturation of the DCP in high magnetic fields determines the DCP spectral dependence. Similar to the samples studied in the present paper, the DCP increase towards the high energy part of the PL spectrum was previously observed for CdSe/CdS dot-in-rods54 and CdTe NCs.24 Recently, it was also reported for InP NCs,26 where also the strong phonon-assisted emission of dark excitons was observed. We suggest that the developed theory of the circular polarized PL with account of the phonon-assisted recombination of dark excitons can explain these and other experimental observations for different ensembles of colloidal nanocrystals.
We introduce a phenomenological magnetic field-induced increase of the 1PL recombination rate. This increase is required for the modeling of the experimentally observed σ− and σ+ polarized PL maxima positions. Note that within our model, the σ− and σ+ polarized PL maxima from the ensemble have the inverted ordering with respect to the exciton Zeeman splitting; however, both shift to higher energies if the 1PL recombination rate is constant in the external magnetic field. Within the second-order perturbation theory, the external magnetic field results only in an increase in the ZPL recombination rate.17,48 An increase in the 1PL recombination rate in the external magnetic field, as well as in the recombination rate with the assistance of two or three optical phonons, requires consideration of higher-order corrections to the recombination rate and will be presented elsewhere. It should be noted that some additional influence on the shift of the polarized PL maxima can come from the Förster energy transfer between NCs in dense ensembles.22,24 It results in spectral diffusion towards lower energy due to the exciton transfer from smaller to larger NCs. These remarks indicate that the analysis of the polarized PL of dark excitons in ensembles of CdSe nanocrystals is a complicated task due to multiple recombination channels with different energy and polarization properties of the emitted light.
The prepared glass samples with CdSe NCs were examined by small angle X-ray scattering (SAXS) with an infinitely high primary beam using Ni-filtered CuKα radiation. We processed the scattering curves, namely the angular dependence of the scattered X-ray intensity, using the Guinier plot56 to obtain the radius of gyration, Rg, of the CdSe crystallites and evaluate the NC diameters with D = 2.58Rg.
(4) |
Here, I+(t) and I−(t) are the intensity of the σ+ and σ− polarized PL, respectively, measured at time delay t after pulsed excitation. For cw laser excitation, the measured polarization degree corresponds to the time-integrated DCP, Pintc, which is calculated as
(5) |
Footnotes |
† Electronic supplementary information (ESI) available: Absorption spectra, bright–dark splitting of excitons in CdSe NCs in glass, degree of circular polarization (DCP) dynamics, model description. See DOI: 10.1039/D0NR07117J |
‡ These authors contributed equally to this work. |
This journal is © The Royal Society of Chemistry 2021 |