Abstract
The fundamental bandgap Eg of a semiconductor—often determined by means of optical spectroscopy—represents its characteristic fingerprint and changes distinctively with temperature. Here, we demonstrate that in magic sized II-VI clusters containing only 26 atoms, a pronounced weakening of the bonds occurs upon optical excitation, which results in a strong exciton-driven shift of the phonon spectrum. As a consequence, a drastic increase of dEg/dT (up to a factor of 2) with respect to bulk material or nanocrystals of typical size is found. We are able to describe our experimental data with excellent quantitative agreement from first principles deriving the bandgap shift with temperature as the vibrational entropy contribution to the free energy difference between the ground and optically excited states. Our work demonstrates how in small nanoparticles, photons as the probe medium affect the bandgap—a fundamental semiconductor property.
Similar content being viewed by others
Introduction
The tunability of the fundamental bandgap, Eg, through size control on the nanometer scale is a key attribute of semiconductor nanocrystals, enabling applications1,2 spanning from biomedical imaging3 to optoelectronic devices such as light emitting diodes4,5 and photodetectors6. The bandgap of semiconductors with reduced dimensions can usually be derived from the bulk value, correcting for electronic perturbations7,8, such as charge carrier quantization and Coulomb interaction. Although the electronic structure of semiconductor nanocrystals is changed by quantum confinement, the temperature dependence of the bandgap usually appears to be dominated by bulk mechanisms in both epitaxial9,10 as well as most colloidal quantum dots (QDs)11,12,13,14,15,16. However, there are few examples17,18,19 that reported conflicting results on the temperature dependence of the bandgap in quantum dots, opening room for debate. For bulk, it is widely accepted that the temperature dependence of the bandgap results from a superposition of electron–phonon interactions and thermal lattice expansion20,21. Besides widely used empirical equations20,22, the AHC approach23,24,25 (after Allen, Heine and Cardona), which is based on second-order perturbation theory, approaches the temperature-induced bandgap shift as the change in energy levels when phonons are thermally populated. This has been used to obtain the temperature dependent bandgap computationally for a number of bulk semiconductors26. Following semi-empirical approaches, AHC was subsequently combined with density functional theory (DFT) or density functional perturbation theory (DFPT)27,28,29,30. For nanocrystals, the temperature dependence of the bandgap has been derived either (i) as time average of the bandgap with first principle molecular dynamic (MD) calculations31,32,33, or (ii) via the frozen-phonon approach34,35,36, summing over the Bose–Einstein weighted changes in eigenenergies due to atomic displacements along the different phonon modes.
Here, we demonstrate that the impact of temperature on the optical bandgap is fundamentally altered in materials at the border between solids and molecules. In magic sized clusters37,38,39,40 (MSC) containing a defined number of 26 atoms, a pronounced exciton-induced weakening of the crystal bonds occurs, resulting in a giant shift of the phonon spectrum driven by the exciton. As a consequence, a drastic increase of dEg/dT (up to a factor of 2) with respect to bulk material or nanocrystals of typical size is found. Using a first principles approach, we computed the exciton-induced change of the phonon density of states from molecular dynamics simulations on the excited state potential energy surface. This enabled us to quantitatively describe our experimental data by deriving the bandgap shift with temperature as the vibrational entropy contribution to the free energy difference between the ground and optically excited states.
Results
Our study benefits from the spectrally narrow absorption and emission resonances of high quality II–VI MSCs that consist of a well-defined number of atoms and vanishing size distribution. This allows us to trace the bandgap energy shift between 5 and 300 K with various optical techniques, as summarized in Fig. 1 for (CdSe)13 MSCs. The absorption of (CdSe)13 MSCs, which is dominated by multiple resonances affiliated to the fine structure of the lowest excitonic transition41, exhibits a redshift with increasing temperature for the whole collectivity of peaks of about 160 meV (Fig. 1a) between cryogenic and room temperature. The short-time photoluminescence (PL) from the bandgap states (Fig. 1b), which is obtained by integrating the emission within the first 20 ps after laser excitation, is found to coincide with the low energy absorption feature of the band edge fine structure, with virtually no Stokes shift within our resolution (compare Supplementary Fig. 1b). This is in distinct contrast to what is observed in typical QDs with Stokes shifts of up to 100 meV for diameters below 2 nm42, and enables us to trace the temperature dependence of an individual state (Fig. 1b). In agreement with absorption, an energy shift of 160 meV is observed between cryogenic and room temperature.
a Absorption and b short-time photoluminescence (PL), measured between 5 K and 300 K. c, d Temperature dependence of the bandgap in Mn2+-doped (CdSe)13 MSCs using optical signatures introduced by the magnetic dopants. c Magnetic circular dichroism (MCD) spectra of Mn:(CdSe)13 MSCs at 1.6 T between 5.2 and 300 K. The signals allow tracking of the bandgap shift with temperature by monitoring the energetic position of the first zero-crossing. d, Photoluminescence excitation (PLE) spectra probed at the 4T1 → 6A1 Mn2+ internal emission between 5.2 and 300 K. The inset depicts the time integrated PL signal of Mn2+-doped (CdSe)13 MSCs revealing the specific orange Mn2+ luminescence, used for probing bandgap states in PLE experiments.
Transition metal doping with manganese (Mn2+)41 offers additional possibilities to track the bandgap shift for specific fine structure states. The sp-d exchange interactions, which persist up to room temperature in these clusters43, generate a pronounced magneto-optical response at the bandgap, representing a powerful tool to distinguish between magneto-optically active and inactive states via magnetic circular dichroism (MCD) spectroscopy (Fig. 1c)41. For MSCs with xMn = 2% (average Mn2+ content among the cations in all clusters of the measured ensemble), the shift of the bandgap slightly exceeds that of the undoped clusters (Fig. 1c), while the mean energy shift between cryogenic and room temperature among a concentration series (xMn = 2–10%, average out of 5 data sets) accounts for 164 ± 16 meV (Supplementary Table 1).
In addition, Mn2+ ions modify the emission of the MSCs, leading to a characteristic orange luminescence originated in the internal Mn2+ 4T1 → 6A1 transition. Monitoring this emission in photoluminescence excitation (PLE) spectroscopy reveals a shift of about 200 meV of the bandgap absorption between 5 and 300 K (see Fig. 1d, 174 ± 23 meV in average among three doped samples). Mn2+-doped (ZnSe)13 MSCs44 prepared for comparison exhibit a similarly enhanced temperature dependence as the CdSe based MSCs (Supplementary Fig. 2). As a reference, we synthesized CdSe QDs covered with amine ligands or a ZnSe shell, which exhibit significantly smaller bandgap shifts with temperature in this regime (average shift of both samples is 93 ± 1 meV in agreement with literature, see Fig. 2 and Supplementary Fig. 3). This demonstrates that our findings are related to the small size of the clusters, irrespective of the material system.
a, b Comparison of the shift in transition energies from cryogenic to room temperature (ΔEg) for undoped (a) and Mn2+-doped (CdSe)13 MSCs (b) as extracted from different experimental approaches. Data for conventional CdSe (light grey) and CdSe/ZnS (dark grey) quantum dots (QDs) are displayed as dots for comparison (see Supplementary Fig. 3 for the spectra). The peak positions for absorption, photoluminescence (PL) and photoluminescence excitation spectroscopy (PLE) signals are extracted via Gaussian fits from the spectra. The peak positions for the energetically lowest magneto-optically active transition are taken from MCD zero-crossings (Fig. 1c). Solid lines represent guides to the eye. c Comparison of the temperature dependence of the band edge transitions for the entity of (CdSe)13 (blue area) and (ZnSe)13 MSCs (black symbols) analyzed within this study to some of the most important group IV, III–V, and II–VI semiconductors (data sets are plotted using parametrization following ref. 21). The bandgap shifts for Mn2+-doped (ZnSe)13 MSCs are extracted from MCD zero-crossings (black circles) and absorption (black stars) (spectra shown in Supplementary Fig. 2).
The bandgap shift with temperature as derived by different techniques for various doped and undoped MSCs are summed up in Fig. 2. The entity of evaluated data in the MSCs exhibits drastically enhanced values of dEg/dT in the high temperature regime as compared to typical CdSe QDs or the most common bulk semiconductors (ref. 21 and references therein). For (CdSe)13 MSCs, our data reveals a mean increase of 89 % compared to bulk CdSe (wurtzite) for the high temperature limit of dEg/dT (between 150 and 300 K) among 14 data sets (Supplementary Table 1).
To reconstruct the impact of electron–phonon coupling on the temperature dependent bandgap shift, it is instructive to interpret the bandgap energy as the standard Gibbs free energy for the formation of an electron-hole pair45. As first suggested in the context of Brook’s theorem46, the bandgap change with temperature due to electron–phonon interactions can be expressed as the change in the vibronic free energy caused by the excitation of an electron from a bonding state in the valence band into a non-bonding state in the conduction band24,47. The limited number of atoms in our material allows us to directly calculate the free energy difference for (CdSe)13 MSCs between the ground and the optically excited states as a function of temperature using constrained DFT. We performed the same set of calculations for the two most prominent structural models for the (CdSe)13 MSCs48: the so-called sliced-wurtzite MSC without ligands (see Supplementary Fig. 4) and a core/cage MSC with methylamine ligands (Fig. 3a). In harmonic approximation, the free energy of the phonons is given by
where g(ω, T) is the phonon density of states (PDOS). The first optically bright transition in the (CdSe)13 MSC is the HOMO→LUMO+1 transition (highest occupied/lowest unoccupied molecular orbital, note that the LUMO describes a mid-gap state, which is related to the surface and frequently observed in CdSe nanoclusters49,50,51,52, see Supplementary Fig. 6).
a Relaxed structural model of a core/cage (CdSe)13 magic sized cluster (MSC) with methylamine (MA) ligands in the ground state electronic configuration. b Normalized PDOS for core/cage (CdSe)13 MSC as shown in a in the ground and HOMO→LUMO+1 excited state electronic configurations, respectively, used to calculate ΔEg(T, gcorecage(ω, T = 0 K)) (blue straight line in c). Note, that beside the modes inherent to the CdSe cluster, the PDOS also contains vibrational and rotational modes stemming from the MA ligands. c Temperature-dependent bandgap shift of (CdSe)13 MSCs calculated using phonon frequencies obtained in harmonic approximation g(ω, T = 0 K)) (blue), including anharmonicities g(ω, T = 300 K) (red) and with full temperature dependence g(ω, T) (green), in comparison to experiment (grey area). Solid and dashed lines denote the core/cage MSC with MA ligands and the sliced-wurtzite MSC without ligands, respectively.
To obtain the excited state PDOS, we constrained the electron occupations of the orbitals, so that one electron from the HOMO is transferred to the LUMO+1 state. We calculated the phonon frequencies, ωi, for the cluster in the ground and the excited state from the force constants. Compared to the MSC in the ground state electronic configuration, the presence of the excited electron weakens the bonds and causes a red-shift of the average phonon frequency of 4.2% for the sliced-wurtzite structure (116.0 cm−1→111.1 cm−1) (individual phonon frequencies in Supplementary Tables 2 and 3). This shift in the mean frequency is almost two orders of magnitude larger than calculated for bulk Si (<0.1% at the melting temperature)24. Inserting the obtained phonon frequencies, the term for the free energy is reduced to a summation. Neglecting the explicit temperature and volume dependence of the PDOS (g(ω, T) ≈ g(ω, T = 0 K)), we obtain the bandgap shift with temperature
as the free energy difference
between the ground and excited state electronic configurations (blue lines in Fig. 3c). For T < 150 K, the measured and calculated temperature-dependence of ΔEg agree very closely. At higher temperatures, the experimental data deviate slightly from a perfectly linear behavior, indicating anharmonic contributions to the free energy.
In order to explicitly account for the temperature-dependence of the phonon frequencies, we used the same constrained DFT approach to perform ab initio MD calculations to obtain the PDOS at 300 K (see Methods). Inserting the ground and excited state PDOSs for T = 300 K, we obtained a slight overestimation for the temperature-dependence of the bandgap (red lines in Fig. 3c). In the case of the sliced-wurtzite (CdSe)13 MSCs, we then approximated the temperature-dependent PDOS g(ω, T) using Gaussian functions as a basis set and interpolated between the PDOSs at T = 0 K and T = 300 K (see “Methods”). As shown in Fig. 3c, the bandgap calculated ab initio from the fully temperature-dependent g(ω, T) (green dashed line) is in excellent agreement with the experiments within the whole measured temperature range.
Discussion
We now discuss our results in connection to current literature. Previous DFT calculations for nanoclusters determined the temperature dependent bandgap either (i) as the mean HOMO-LUMO energy difference in a MD31,32,33, which distributes the atomic velocities according to classical Maxwell-Boltzmann statistics instead of Bose–Einstein occupation of the vibrational modes, and thus inherently overestimates dEg/dT at low temperatures, or (ii) by summing over Bose–Einstein weighted eigenenergy changes induced by frozen-phonon modes34,35,36. The latter is known to severely overestimate the temperature dependence of the bandgap in nanocrystals35,36, often attributed to overrated contributions of surface and ligand modes. Indeed, we obtain an overestimation by almost an order of magnitude, calculating the bandgap shift with temperature based on a frozen-phonon approach (compare Supplementary Fig. 5). The failure of the commonly used frozen-phonon approach, which is very successful in the bulk, indicates that there is something fundamentally different about the temperature-dependence of the band gap in NCs. The approach presented in this work, which explicitly incorporates the effect of photon absorption on phonon modes, is able to successfully describe dEg/dT over the entire temperature range between 5 and 300 K with quantitative accuracy and without any free parameters. Brook’s theorem46 states that the temperature-dependence of the bandgap can be obtained equally from either the changes in the eigenenergies induced by the phonons or the changes in the vibrational frequencies induced by the optical excitations. Clearly, our results demonstrate that Brook’s theorem is invalid for small nanoclusters and they highlight the importance of taking excitonic effects into account explicitly.
For high temperatures (ℏω ≪ kBT), the slope of the temperature-induced band gap shift can be estimated as
disclosing two factors which might enhance dEg/dT in nanoclusters compared to bulk: (i) Softer phonon modes in general compared to bulk, e.g. due to surface related acoustic modes, and (ii) an enhanced change Δω in the phonon energies due to the absorption process. The former has been predicted for PbS53 and Si35, while for CdSe nanoparticles theoretically36,54,55 and experimentally54 a blueshift of the acoustic surface modes has been observed. Our calculations reveal a distinct redshift in the acoustic modes for the core-cage nanocluster model compared to bulk, while the data exhibit no clear shift for the sliced-wurtzite structure (see Supplementary Figure 4c). Both structures reveal enhanced temperature-induced bandgap shifts (red lines in Fig. 3c), indicating that the main origin lies in the significant change in the PDOS of clusters in the presence of an exciton. Due to the small size of the clusters, the exciton is distributed over a minimal number of 40–50 bonds. Note that one exciton per cluster corresponds to a photo-generated charge carrier density in the range of 1021 cm−3, representing the regime of particularly high optical pumping in conventional semiconductors.
Our results not only highlight how the absorption of one photon as the probe medium influences the observed bandgap, but also visibly demonstrate how the temperature dependence of the bandgap immediately depends on the phonon dispersion of a material. We expect these findings to be relevant in materials that are either quantum-confined, strongly correlated or exhibit strong electron–phonon coupling and thus optical excitations significantly alter the vibrational free energies and in turn derived properties.
Methods
Sample preparation
Nanocrystals and MSCs used in this study were synthesized using standard Schlenk techniques under an argon atmosphere. The MSCs were obtained from a Lewis acid-base reaction between metal-ammine halide complexes (CdCl2(octylamine)2, ZnCl2(octylamine)2, and MnCl2(octylamine)2) and octylammonium selenocarbamate (0.67 M) in n-octylamine (Aldrich, 99%). For Mn2+:(CdSe)13 MSCs, 1.0 mmol of CdCl2 and 0.1 mmol of MnCl2 were heated in 7 ml of n-octylamine at 120 °C for 2 h. The average doping concentration is changed by adjusting the initial precursor ratio (MnCl2/CdCl2). Octylammonium selenocarbamate was prepared by bubbling CO gas into 3 ml of n-octylamine containing 2.0 mmol Se powder for 1 h at room temperature. This solution was injected into the as-prepared precursor mixture containing metal-ammine halide complexes and kept at room temperature for more than 24 h. The MSCs were precipitated by adding ethanol containing trioctylphosphine followed by washing several times with ethanol1,2,3. For the synthesis of Mn2+:(ZnSe)13 MSCs, ZnCl2 was used instead of CdCl2. The average doping concentration of the MSCs was measured by inductively coupled plasma-atomic emission spectrometry (ICP-AES, Shimadzu ICPS-7500). To enhance the dispersibility of MSCs, which results in minimization of light scattering, the initial octylamine ligands were exchanged with long-chain oleylamine ligands1. CdSe QDs and CdSe/ZnS QDs were synthesized by conventional hot injection method using Cd(oleate)2 and trioctylphosphine selenide (TOPSe) as precursors4,56. In all, 1.2 mmol CdO were added to a mixture of 1.5 mL of oleic acid and 20 mL of 1-octadecene and heated under vacuum for 2 h at 120 °C to form Cd(oleate)2 complexes. The mixture was then heated to 300 °C under Ar and TOPSe (1 M) was rapidly injected to obtain CdSe QDs. For the growth of ZnS shells, 4.8 mmol of Zn(oleate)2 and 4.8 mmol of tributylphosphine sulfide were additionally introduced to the as-prepared solution containing CdSe QDs and the reaction mixture was kept at 300 °C for 20 min. For CdSe QDs and CdSe/ZnS QDs, we also performed a surface modification with oleylamine to neglect the effect of the surface ligands.
Optical characterization
For absorption, photoluminescence excitation spectroscopy (PLE) and magnetic circular dichroism (MCD) measurements the MSCs passivated by oleylamine were prepared as thin films between two quartz glass substrates, and for PL measurements between a silicon and a quartz glass substrate. MCD spectroscopy was conducted on a homemade setup consisting of a 75 W Xenon lamp (Lot-Oriel) equipped with a monochromator (omni-λ 150, Lot-Oriel) and a photomultiplier (R928, Hamamatsu). The excitation light was modulated using a photoelastic modulator (PEM-90, Hinds Instruments). The sample was placed in a helium vapor cryostat (ST-300, Janis) between two poles of an electromagnet (EM4-HVA, Lake Shore) in Faraday geometry. Absorption measurements were either extracted from the DC signal during MCD measurements or collected with an UV-VIS spectrophotometer (Shimadzu 2550) equipped with a helium vapor cryostat (ST-300, Janis). PLE measurements were conducted with a helium vapor cryostat (ST-300, Janis) placed in a Fluorolog-3 (FL3-22, Horiba Scientific). For time resolved analysis of the PL signal the samples placed in a helium vapor cryostat (ST-300, Janis) were excited with a frequency tripled Ti:sapphire laser (Mira 900, Coherent, 270 nm) pumped with an Nd:YAG laser (Verdi-V10, Coherent) with 100 fs pulses at a repetition rate of 76 MHz. The signal was detected with a streak camera system consisting of a monochromator (250IS, Bruker Optics) and a Synchroscan Streak camera (C5680-24 c, Hamamatsu), providing a temporal resolution of 4 ps.
DFT calculations
We carried out first principles calculations with density functional theory and plane-wave basis sets, within the Perdew-Burke-Ernzerhof (PBE)57,58 approximation as implemented in the Quantum Espresso package59. We used norm-conserving pseudopotentials and a wavefunction energy cutoff of 45 Ry. The molecular dynamics (MD) calculations were performed at 300 K, using a time step of 2 fs and the Berendsen60 thermostat with a rise time of 10 time steps. The clusters were equilibrated for 25 ps and subsequently data was sampled for 100 ps. Excited states were described by pair excitations61 within constrained DFT, where the Kohn-Sham energy is minimized under the constraint that the occupation of the highest occupied molecular orbital (HOMO) is reduced by 1 e−, compared to the ground state, and the occupation of the lowest energy state directly above the lowest unoccupied molecular orbital (LUMO + 1) is increased by 1 e−. To prevent convergence issues within the electronic loop due to level crossings and concurrent discontinuous changes in the charge density, we implemented a smearing scheme into Quantum Espresso, where a Gaussian smearing is applied individually to the electron and the hole with a smearing of σ = 0.001 Ry. At T = 0 K, phonon eigenfrequencies and eigenvectors are calculated from finite differences. At T = 300 K, phonon densities of states (PDOS) are computed from the MD trajectories using the maximum entropy method (MEM)62,63, the same numerical setup as above and 4096 poles in the frequency interval from 0 cm−1 to 300 cm−1. The spectra are renormalized so that the integrals ∫dωg(ω) are equal to the number of phonon modes. To calculate g(ω, T), we approximate the set of discrete phonon frequencies for T = 0 K by Gaussians with a very small but finite width σ = 2 cm−1 according to:
Subsequently. g(ω, T = 300 K) is expanded into the same number of Gaussians, albeit with modified positions and widths. All resulting PDOS within the Gaussian basis set have been confirmed to yield the same free energies as before. We then linearly interpolate the widths and positions of the Gaussians as a function of temperature to obtain g(ω, T).
Code availability
The computer code used for the DFT and MD calculations during the study is available from S.W. on request.
References
Kovalenko, M. V. et al. Prospects of Nanoscience with Nanocrystals. ACS Nano 9, 1012–1057 (2015).
Akkerman, Q. A., Rainò, G., Kovalenko, M. V. & Manna, L. Genesis, challenges and opportunities for colloidal lead halide perovskite nanocrystals. Nat. Mater. 17, 394–405 (2018).
Medintz, I. L., Uyeda, H. T., Goldman, E. R. & Mattoussi, H. Quantum dot bioconjugates for imaging, labelling and sensing. Nat. Mater. 4, 435–446 (2005).
Choi, M. K., Yang, J., Hyeon, T. & Kim, D.-H. Flexible quantum dot light-emitting diodes for next-generation displays. npj Flex. Electron. 2, 10 (2018).
Shirasaki, Y., Supran, G. J., Bawendi, M. G. & Bulović, V. Emergence of colloidal quantum-dot light-emitting technologies. Nat. Photonics 7, 13–23 (2013).
García de Arquer, F. P., Armin, A., Meredith, P. & Sargent, E. H. Solution-processed semiconductors for next-generation photodetectors. Nat. Rev. Mater. 2, 16100 (2017).
Efros, A. et al. Band-edge exciton in quantum dots of semiconductors with a degenerate valence band: dark and bright exciton states. Phys. Rev. B 54, 4843–4856 (1996).
Norris, D., Efros, A., Rosen, M. & Bawendi, M. Size dependence of exciton fine structure in CdSe quantum dots. Phys. Rev. B. Condens. Matter 53, 16347–16354 (1996).
Ortner, G. et al. Temperature dependence of the excitonic band gap in InxGa1−xAs/GaAs self-assembled quantum dots. Phys. Rev. B 72, 085328 (2005).
Sanguinetti, S. et al. Electron-phonon interaction in individual strain-free GaAs∕Al0.3Ga0.7As quantum dots. Phys. Rev. B 73, 1–7 (2006).
Liptay, T. J. & Ram, R. J. Temperature dependence of the exciton transition in semiconductor quantum dots. Appl. Phys. Lett. 89, 223132 (2006).
Joshi, A., Narsingi, K. Y., Manasreh, M. O., Davis, E. A. & Weaver, B. D. Temperature dependence of the band gap of colloidal CdSe/ZnS core/shell nanocrystals embedded into an ultraviolet curable resin. Appl. Phys. Lett. 89, 89–92 (2006).
Valerini, D. et al. Temperature dependence of the photoluminescence properties of colloidal CdSe∕ZnS core/shell quantum dots embedded in a polystyrene matrix. Phys. Rev. B 71, 1–6 (2005).
Jing, P. et al. Temperature-dependent photoluminescence of CdSe-Core CdS/CdZnS/ZnS-multishell quantum dots. J. Phys. Chem. C. 113, 13545–13550 (2009).
Nomura, S. & Kobayashi, T. Exciton–LO-phonon couplings in spherical semiconductor microcrystallites. Phys. Rev. B 45, 1305–1316 (1992).
Korsunska, N. E., Dybiec, M., Zhukov, L., Ostapenko, S. & Zhukov, T. Reversible and non-reversible photo-enhanced luminescence in CdSe/ZnS quantum dots. Semicond. Sci. Technol. 20, 876–881 (2005).
Olkhovets, A., Hsu, R.-C., Lipovskii, A. & Wise, F. Size-dependent temperature variation of the energy gap in lead-salt quantum dots. Phys. Rev. Lett. 81, 3539–3542 (1998).
Wise, F. W. Lead salt quantum dots: the limit of strong quantum confinement. Acc. Chem. Res. 33, 773–780 (2000).
Vossmeyer, T. et al. CdS nanoclusters: synthesis, characterization, size dependent oscillator strength, temperature shift of the excitonic transition energy, and reversible absorbance shift. J. Phys. Chem. 98, 7665–7673 (1994).
O’Donnell, K. P. & Chen, X. Temperature dependence of semiconductor band gaps. Appl. Phys. Lett. 58, 2924–2926 (1991).
Pässler, R. Dispersion-related description of temperature dependencies of band gaps in semiconductors. Phys. Rev. B 66, 085201 (2002).
Varshni, Y. P. Temperature dependence of the energy gap in semiconductors. Physica 34, 149–154 (1967).
Allen, P. B. & Heine, V. Theory of the temperature dependence of electronic band structures. J. Phys. C. Solid State Phys. 9, 2305–2312 (1976).
Heine, V. & Van Vechten, J. A. Effect of electron-hole pairs on phonon frequencies in Si related to temperature dependence of band gaps. Phys. Rev. B 13, 1622–1626 (1976).
Allen, P. B. & Cardona, M. Temperature dependence of the direct gap of Si and Ge. Phys. Rev. B 27, 4760–4769 (1983).
Ridley, B. K. Quantum Process in Semiconductors. (Oxford University Press, 1993).
Giustino, F., Louie, S. G. & Cohen, M. L. Electron-phonon renormalization of the direct band gap of diamond. Phys. Rev. Lett. 105, 1–4 (2010).
Antonius, G., Poncé, S., Boulanger, P., Côté, M. & Gonze, X. Many-body effects on the zero-point renormalization of the band structure. Phys. Rev. Lett. 112, 1–5 (2014).
Saidi, W. A., Poncé, S. & Monserrat, B. Temperature dependence of the energy levels of methylammonium lead iodide perovskite from first principles. J. Phys. Chem. Lett. 7, 5247–5252 (2016).
Poncé, S. et al. Temperature dependence of the electronic structure of semiconductors and insulators. J. Chem. Phys. 143, 102813 (2015).
Franceschetti, A. First-principles calculations of the temperature dependence of the band gap of Si nanocrystals. Phys. Rev. B - Condens. Matter Mater. Phys. 76, 16–18 (2007).
Kamisaka, H., Kilina, S. V., Yamashita, K. & Prezhdo, O. V. Ab initio study of temperature and pressure dependence of energy and phonon-induced dephasing of electronic excitations in CdSe and PbSe quantum dots †. J. Phys. Chem. C. 112, 7800–7808 (2008).
Ibrahim, Z. A. et al. Temperature dependence of the optical response: application to bulk GaAs using first-principles molecular dynamics simulations. Phys. Rev. B 77, 125218 (2008).
Capaz, R. B., Spataru, C. D., Tangney, P., Cohen, M. L. & Louie, S. G. Temperature dependence of the band gap of semiconducting carbon nanotubes. Phys. Rev. Lett. 94, 1–4 (2005).
Han, P. & Bester, G. Large nuclear zero-point motion effect in semiconductor nanoclusters. Phys. Rev. B 88, 165311 (2013).
Han, P. & Bester, G. Fundamental difference between measured and calculated exciton-phonon coupling in nanostructures. Phys. Rev. B 99, 100302 (2019).
Kasuya, A. et al. Ultra-stable nanoparticles of CdSe revealed from mass spectrometry. Nat. Mater. 3, 99–102 (2004).
Yu, J. H. et al. Giant Zeeman splitting in nucleation-controlled doped CdSe:Mn2+ quantum nanoribbons. Nat. Mater. 9, 47–53 (2010).
Cossairt, B. M. & Owen, J. S. CdSe clusters: at the interface of small molecules and quantum dots. Chem. Mater. 23, 3114–3119 (2011).
Bowers, M. J., McBride, J. R. & Rosenthal, S. J. White-light emission from magic-sized cadmium selenide nanocrystals. J. Am. Chem. Soc. 127, 15378–15379 (2005).
Yang, J. et al. Route to the smallest doped semiconductor: Mn+2-Doped (CdSe)13 clusters. J. Am. Chem. Soc. 137, 12776–12779 (2015).
Kuno, M., Lee, J. K., Dabbousi, B. O., Mikulec, F. V. & Bawendi, M. G. The band edge luminescence of surface modified CdSe nanocrystallites: probing the luminescing state. J. Chem. Phys. 106, 9869–9882 (1997).
Muckel, F. et al. Digital doping in magic-sized CdSe clusters. ACS Nano 10, 7135–7141 (2016).
Yang, J. et al. Chemical synthesis, doping, and transformation of magic-sized semiconductor alloy nanoclusters. J. Am. Chem. Soc. 139, 6761–6770 (2017).
Thurmond, C. D. The standard thermodynamic functions for the formation of electrons and holes in Ge, Si, GaAs, and GaP. J. Electrochem. Soc. 122, 1133 (1975).
Brooks, H. Theory of the electrical properties of germanium and silicon. Adv. Electr. 7, 85–182 (1955).
King-Smith, R. D., Needs, R. J., Heine, V. & Hodgson, M. J. A first-principle calculation of the temperature dependence of the indirect band gap of silicon. Epl 10, 569–574 (1989).
Azpiroz, J. M., Matxain, J. M., Infante, I., Lopez, X. & Ugalde, J. M. A DFT/TDDFT study on the optoelectronic properties of the amine-capped magic (CdSe)13 nanocluster. Phys. Chem. Chem. Phys. 15, 10996 (2013).
Bawendi, M. G., Carroll, P. J., Wilson, W. L. & Brus, L. E. Luminescence properties of CdSe quantum crystallites: resonance between interior and surface localized states. J. Chem. Phys. 96, 946–954 (1992).
Landes, C. F., Braun, M. & El-Sayed, M. A. On the nanoparticle to molecular size transition: fluorescence quenching studies. J. Phys. Chem. B 105, 10554–10558 (2001).
Puzder, A., Williamson, A. J., Gygi, F. & Galli, G. Self-healing of CdSe nanocrystals: first-principles calculations. Phys. Rev. Lett. 92, 1–4 (2004).
Vörös, M., Galli, G. & Zimanyi, G. T. Colloidal nanoparticles for intermediate band solar cells. ACS Nano 9, 6882–6890 (2015).
Bozyigit, D. et al. Soft surfaces of nanomaterials enable strong phonon interactions. Nature 531, 618–622 (2016).
Shi, C. et al. Size-dependent lattice dynamics of atomically precise cadmium selenide quantum dots. Phys. Rev. Lett. 122, 26101 (2019).
Han, P. & Bester, G. Insights about the surface of colloidal nanoclusters from their vibrational and thermodynamic properties. J. Phys. Chem. C. 116, 10790–10795 (2012).
Choi, M. K. et al. Extremely vivid, highly transparent, and ultrathin quantum dot light-emitting diodes. Adv. Mater. 30, 1703279 (2018).
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865–3868 (1996).
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865 (1996); erratum 78, 1396 (1997).
Giannozzi, P. et al. QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials. J. Phys. Condens. Matter 21, 395502 (2009).
Berendsen, H. J. C., Postma, J. P. M., van Gunsteren, W. F., DiNola, A. & Haak, J. R. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 81, 3684–3690 (1984).
Weissker, H.-C., Furthmüller, J. & Bechstedt, F. Structure- and spin-dependent excitation energies and lifetimes of Si and Ge nanocrystals from ab initio calculations. Phys. Rev. B 69, 115310 (2004).
Jaynes, E. T. Information theory and statistical mechanics. Phys. Rev. 106, 620–630 (1957).
Press, W., Teukolsky, S., Vetterling, W. & Flannery, B. Numerical Recipes in C. (Cambridge University Press, 2002)
Acknowledgements
F.M., S.L., and G.B. gratefully acknowledge financial support from the German Research Foundation DFG under contract Ba 1422/13-2. F.M. in addition was supported by the German Academic Exchange Service (DAAD) with funds from the German Federal Ministry of Education and Research (BMBF) and the European Union (FP7-PEOPLE-2013-COFUND - grant agreement n° 605728). T.A.N., E.S., and S.W. were supported by the German Federal Ministry of Education and Research (BMBF) within the NanoMatFutur programme, grant no. 13N12972. Supercomputer time provided by NERSC (project no. 35687) and the Max-Planck Computing and Data Facility, Garching, is gratefully acknowledged. T.H. acknowledges financial support by the Research Center Program of the Institute for Basic Science (IBS) in Korea (IBS-R006-D1). J.Y. acknowledges the financial support by the DGIST Start-up Fund Program of the Ministry of Science, ICT and Future Planning (2019070014).
Author information
Authors and Affiliations
Contributions
F.M. and G.B. conceived and designed the experiments, led the research and wrote the paper. J.Y. and T.H. designed and prepared the materials. S.L. contributed to the optical measurements, data analysis and interpretation. S.W. conceived and designed the DFT calculations, T.A.N. and E.S. conducted the DFT calculations. All authors discussed the results and assisted in paper preparation.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Peer review information Nature Communications thanks the anonymous reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Source data
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Muckel, F., Lorenz, S., Yang, J. et al. Exciton-driven change of phonon modes causes strong temperature dependent bandgap shift in nanoclusters. Nat Commun 11, 4127 (2020). https://doi.org/10.1038/s41467-020-17563-0
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-020-17563-0