Time-domain ab initio modeling of excitation dynamics in quantum dots

Nanocrystal Assemblies: Current Advances and Open Problems

We explore the potential of nanocrystals (a term used equivalently to nanoparticles) as building blocks for nanomaterials, and the current advances and open challenges for fundamental science developments and applications. Nanocrystal assemblies are inherently multiscale, and the generation of revolutionary material properties requires a precise understanding of the relationship between structure and function, the former being determined by classical effects and the latter often by quantum effects. With an emphasis on theory and computation, we discuss challenges that hamper current assembly strategies and to what extent nanocrystal assemblies represent thermodynamic equilibrium or kinetically trapped metastable states. We also examine dynamic effects and optimization of assembly protocols. Finally, we discuss promising material functions and examples of their realization with nanocrystal assemblies.

Ultrafast Charge Transfer and Recombination Dynamics in Monolayer-Multilayer WSe2 Junctions Revealed by Time-Resolved Photoemission Electron Microscopy

The ultrafast carrier dynamics of junctions between two chemically identical, but electronically distinct, transition metal dichalcogenides (TMDs) remains largely unknown. Here, we employ time-resolved photoemission electron microscopy (TR-PEEM) to probe the ultrafast carrier dynamics of a monolayer-to-multilayer (1L-ML) WSe2 junction. The TR-PEEM signals recorded for the individual components of the junction reveal the sub-ps carrier cooling dynamics of 1L- and 7L-WSe2, as well as few-ps exciton-exciton annihilation occurring on 1L-WSe2. We observe ultrafast interfacial hole (h) transfer from 1L- to 7L-WSe2 on an ∼0.2 ps time scale. The resultant excess h density in 7L-WSe2 decays by carrier recombination across the junction interface on an ∼100 ps time scale. Reminiscent of the behavior at a depletion region, the TR-PEEM image reveals the h density accumulation on the 7L-WSe2 interface, with a decay length ∼0.60 ± 0.17 μm. These charge transfer and recombination dynamics are in agreement with ab initio quantum dynamics. The computed orbital densities reveal that charge transfer occurs from the basal plane, which extends over both 1L and ML regions, to the upper plane localized on the ML region. This mode of charge transfer is distinctive to chemically homogeneous junctions of layered materials and constitutes an additional carrier deactivation pathway that should be considered in studies of 1L-TMDs found alongside their ML, a common occurrence in exfoliated samples.

Control of Polaronic Behavior and Carrier Lifetimes via Metal and Anion Alloying in Chalcogenide Perovskites

Transition-metal perovskite chalcogenides (TMPCs) have emerged as lead-free alternatives to lead-halide perovskites and have been currently of increasing interest for optoelectronic applications because of their suitable band gaps, high carrier mobility, strong light absorption, and high stability. Here, we systematically report a study on the effects of Ti- and Se-alloying strategies on polaron behavior and carrier lifetimes in nonradiative recombination. Although such alloying can effectively tune the band gap of BaZrS₃, we observe localized small polaron formation upon Ti alloying and large polarons generating in Se alloying. Ti-alloying strengthens the electron-phonon coupling, leading to a reduced carrier lifetime. Remarkably, Se-alloying weakens the electron-phonon coupling and prolongs the nonradiative electron-hole recombination lifetime by up to 60% compared to that in pristine BaZrS₃ material. The simulations rationalize the difference in carrier lifetimes in TMPC alloys and provide guidelines for further improvements in TMPC-based photoelectronic devices.

Ultra-narrow room-temperature emission from single CsPbBr3 perovskite quantum dots

Semiconductor quantum dots have long been considered artificial atoms, but despite the overarching analogies in the strong energy-level quantization and the single-photon emission capability, their emission spectrum is far broader than typical atomic emission lines. Here, by using ab-initio molecular dynamics for simulating exciton-surface-phonon interactions in structurally dynamic CsPbBr₃ quantum dots, followed by single quantum dot optical spectroscopy, we demonstrate that emission line-broadening in these quantum dots is primarily governed by the coupling of excitons to low-energy surface phonons. Mild adjustments of the surface chemical composition allow for attaining much smaller emission linewidths of 35−65 meV (vs. initial values of 70–120 meV), which are on par with the best values known for structurally rigid, colloidal II-VI quantum dots (20−60 meV). Ultra-narrow emission at room-temperature is desired for conventional light-emitting devices and paramount for emerging quantum light sources.

Narrow emission is desired for light-emitting devices. Here, Kovalenko et al. demonstrate that the emission line-broadening in perovskite quantum dots is dominated by the coupling between excitons and surface phonon modes which can be controlled by minimal surface modifications.

Increasing Efficiency of Nonadiabatic Molecular Dynamics by Hamiltonian Interpolation with Kernel Ridge Regression

Nonadiabatic (NA) molecular dynamics (MD) goes beyond the adiabatic Born-Oppenheimer approximation to account for transitions between electronic states. Such processes are common in molecules and materials used in solar energy, optoelectronics, sensing, and many other fields. NA-MD simulations are much more expensive compared to adiabatic MD due to the need to compute excited state properties and NA couplings (NACs). Similarly, application of machine learning (ML) to NA-MD is more challenging compared with adiabatic MD. We develop an NA-MD simulation strategy in which an adiabatic MD trajectory, which can be generated with a ML force field, is used to sample excitation energies and NACs for a small fraction of geometries, while the properties for the remaining geometries are interpolated with kernel ridge regression (KRR). This ML strategy allows for one to perform NA-MD under the classical path approximation, increasing the computational efficiency by over an order of magnitude. Compared to neural networks, KRR requires little parameter tuning, saving efforts on model building. The developed strategy is demonstrated with two metal halide perovskites that exhibit complicated MD and are actively studied for various applications.

Engineering of Exciton Spatial Distribution in CdS Nanoplatelets

Zinc-blende CdS nanoplatelets with atomically flat and very large {100} basal planes terminated solely by one type of element (either Cd or S atoms) are synthesized. Optical spectroscopy, X-ray diffraction, X-ray absorption, and transmission electron microscopy confirm that the surface structures of newly developed S-terminated CdS nanoplatelets are at least as well-defined as the original Cd-terminated nanoplatelets. Band gaps of the nanoplatelets are found to depend on not only the quantum-confined dimension (thickness) but also the nature of the surface termination. The facet structure dictates the packing of the ligands (carboxylate for Cd-terminated nanoplatelets and alkyl for S-terminated nanoplatelets), which causes a difference in the lattice strain and significantly affects the optical spectral width. Experimental and theoretical results reveal that engineering the exciton spatial distribution by the tailored synthesis of semiconductor nanocrystals with a precisely controlled surface structure is fully possible, which should open a new door for delivering the long-promised potential of semiconductor nanocrystals.

Dynamic Formation of Metal-Based Traps in Photoexcited Colloidal Quantum Dots and Their Relevance for Photoluminescence

Trap states play a crucial role in the design of colloidal quantum dot (QD)-based technologies. The presence of these in-gap states can either significantly limit the efficiency of devices (e.g., in solar cells or LEDs) or play a pivotal role in the functioning of the technology (e.g., in catalysis). Understanding the atomistic nature of traps is therefore of the highest importance. Although the mechanism through which undercoordinated chalcogenide atoms can lead to trap states in II-VI QDs is generally well understood, the nature of metal-based traps remains more elusive. Previous research has shown that reduction of metal sites in negatively charged QDs can lead to in-gap states. Here, we use density functional theory to show that metal-based traps are also formed in charge-neutral but photoexcited CdSe QDs. It is found that Cd-Cd dimers and the concomitant trap states are transient in nature and appear and disappear on the picosecond time scale. Subsequent nonradiative recombination from the trap is shown to be much faster than radiative recombination, indicating that dimer-related trap states can quench the photoluminescence. These results are expected to be transferable to other II-VI materials and highlight the importance of surface redox reactions for the optical properties of QDs. Moreover, they show that photoexcitation can lead to atomic rearrangements on the surface and thus create transient in-gap states.

Effect of Ligand Functionalization on the Rate of Charge Carrier Recombination in Metal-Organic Frameworks: A Case Study of MIL-125

Ligand functionalization is a powerful approach for modifying the electronic structure of metal-organic frameworks when targeting the optimal electronic properties for photocatalysis and photovoltaics. However, its effect on the charge carrier lifetimes and recombination pathways remains unexplored. In this work, first-principles simulations, including nonadiabatic molecular dynamics, are performed for the representative TiO₂-based metal-organic framework systems MIL-125-X to unravel the impact of ligand functionalization on the nonradiative electron-hole recombination process, decoherence rates, and phonon modes giving the largest contribution to the nonradiative decay. Nonradiative recombination rates, simulated using the PBE0 density functional, are in excellent agreement with experiment. The ligand functionalization in MIL-125-X influences the recombination rates, unraveling the trend opposite to the evolution of the band gap and affecting the nonadiabatic coupling coefficients. Ligand modification impacts the phonon modes, which contribute most to the recombination process, altering the distribution between soft phonon modes and vibrational modes associated with specific structural motifs.

Ab initio quantum dynamics of charge carriers in graphitic carbon nitride nanosheets

Graphitic carbon nitride (g-C₃N₄), a metal-free and visible light responsive photocatalyst, has garnered much attention due to its wide range of applications. In order to elucidate the role of dimensionality on the properties of photo-generated charge carriers, we apply nonadiabatic (NA) molecular dynamics combined with time-domain density functional theory to investigate nonradiative relaxation of hot electrons and holes, and electron-hole recombination in monolayer and bulk g-C₃N₄. The nonradiative charge recombination occurs on a nanosecond timescale and is faster in bulk than the nanosheet, in agreement with the experiment. The difference arises due to the smaller energy gap and participation of additional vibrations in the bulk system. The long carrier lifetimes are favored by small NA coupling and rapid phonon-induced loss of quantum coherence between the excited and ground electronic states. Decoherence is fast because g-C₃N₄ is soft and undergoes large scale vibrations. The NA coupling is small since electrons and holes are localized on different atoms, and the electron-hole overlap is relatively small. Phonon-driven relaxation of hot electrons and holes takes 100-200 fs and is slightly slower at higher initial energies due to participation of fewer vibrational modes. This feature of two-dimensional g-C₃N₄ contrasts traditional three-dimensional semiconductors, which exhibit faster relaxation at higher energies due to larger density of states, and can be used to extract hot carriers to perform useful functions. The ab initio quantum dynamics simulations present a comprehensive picture of the photo-induced charge carrier dynamics in g-C₃N₄, guiding design of photovoltaic and photocatalytic devices.

Optoelectronic Properties of Two-Dimensional Bromide Perovskites: Influences of Spacer Cations

Two-dimensional (2D) halide perovskites have displayed unique emission properties, making them potential candidates for next-generation light-emitting devices. Here, we combine nonadiabatic molecular dynamics and time-domain density functional theory to investigate the fundamental mechanisms of carrier recombination processes. Considering monolayer bromide perovskites with dissimilar organic spacer molecules, n-butylammonium (BA) and phenylethylammonium (PEA) cations, we find a strong correlation between temperature-induced structural fluctuations and nonradiative carrier recombination rates in these materials. The more flexible geometry of (BA)₂PbBr₄ compared to that of (PEA)₂PbBr₄, results in faster electron-hole recombination and shorter carrier lifetime, diminishing the photoluminescence quantum yield for softer 2D perovskites. Reduced structural fluctuations in relatively rigid (PEA)₂PbBr₄ not only indicate of a longer carrier lifetime but also suggest a narrower emission line width, implying a higher purity of the emitted light. Our ab initio modeling of excited state properties in 2D perovskites conveys material designing strategies to fine-tune perovskite emissions for solid-state lighting applications.

Phonon-Mediated and Weakly Size-Dependent Electron and Hole Cooling in CsPbBr3 Nanocrystals Revealed by Atomistic Simulations and Ultrafast Spectroscopy

We combine state-of-the-art ultrafast photoluminescence and absorption spectroscopy and nonadiabatic molecular dynamics simulations to investigate charge-carrier cooling in CsPbBr3 nanocrystals over a very broad size regime, from 0.8 to 12 nm. Contrary to the prevailing notion that polaron formation slows down charge-carrier cooling in lead-halide perovskites, no suppression of carrier cooling is observed in CsPbBr3 nanocrystals except for a slow cooling (over ∼10 ps) of "warm" electrons in the vicinity (within ∼0.1 eV) of the conduction band edge. At higher excess energies, electrons and holes cool with similar rates, on the order of 1 eV ps-1 carrier-1, increasing weakly with size. Our ab initio simulations suggest that cooling proceeds via fast phonon-mediated intraband transitions driven by strong and size-dependent electron-phonon coupling. The presented experimental and computational methods yield the spectrum of involved phonons and may guide the development of devices utilizing hot charge carriers.

Photoexcited carrier dynamics in colloidal quantum dot solar cells: insights into individual quantum dots, quantum dot solid films and devices

The certified power conversion efficiency (PCE) record of colloidal quantum dot solar cells (QDSCs) has considerably improved from below 4% to 16.6% in the last few years. However, the record PCE value of QDSCs is still substantially lower than the theoretical efficiency. So far, there have been several reviews on recent and significant achievements in QDSCs, but reviews on photoexcited carrier dynamics in QDSCs are scarce. The photovoltaic performances of QDSCs are still limited by the photovoltage, photocurrent and fill factor that are mainly determined by the photoexcited carrier dynamics, including carrier (or exciton) generation, carrier extraction or transfer, and the carrier recombination process, in the devices. In this review, the photoexcited carrier dynamics in the whole QDSCs, originating from individual quantum dots (QDs) to the entire device as well as the characterization methods used for analyzing the photoexcited carrier dynamics are summarized and discussed. The recent research including photoexcited multiple exciton generation (MEG), hot electron extraction, and carrier transfer between adjacent QDs, as well as carrier injection and recombination at each interface of QDSCs are discussed in detail herein. The influence of photoexcited carrier dynamics on the physiochemical properties of QDs and photovoltaic performances of QDSC devices is also discussed.

Investigating exciton structure and dynamics in colloidal CdSe quantum dots with two-dimensional electronic spectroscopy

Two-Dimensional Electronic Spectroscopy (2DES) is performed on CdSe colloidal quantum dots. These experiments reveal new observations on exciton structure and dynamics in quantum dots, expanding upon prior transient absorption measurements of excitonics in these systems. The 2DES method enables the separation of line broadening mechanisms, thereby better revealing the excitonic lineshapes and biexcitonic interactions. 2DES enables more information rich spectral probing of coherent phonons and their coupling to excitons. The data show spectral modulations and drifts, with differences based upon whether one monitors the excitation energy (E₁) or emission energy (E₃). These measurements reveal both homogeneous and inhomogeneous broadenings, as well as static and dynamic line broadening. The longitudinal optical phonon modulates the dynamic absorption spectrum both in energy and linewidth. These experiments enable measurement of hot exciton cooling with improved resolution in energy and time. These 2DES results are consistent with prior excitonic state-resolved transient absorption measurements, albeit with the addition of contributions due to coherent phonons. Finally these 2DES experiments enable disentangling of coupling versus relaxation contributions to the signals, further offering a test of electronic structure theory.

Decoherence Allows Model Reduction in Nonadiabatic Dynamics Simulations

A nonadiabatic (NA) molecular dynamics (MD) simulation requires calculation of NA coupling matrix elements, the number of which scales as a square of the number of basis states. The basis size can be huge in studies of nanoscale materials, and calculation of the NA couplings can present a significant bottleneck. A quantum-classical approximation, NAMD overestimates coherence in the quantum, electronic subsystem, requiring decoherence correction. Generally, decoherence times decrease with increasing energy separation between pairs of states forming coherent superpositions. Since rapid decoherence stops quantum dynamics, one expects that decoherence-corrected NAMD can eliminate the need for calculation of NA couplings between energetically distant states, notably reducing the computational cost. Considering several types of dynamics in a semiconductor quantum dot, we demonstrate that indeed, decoherence allows one to reduce the number of needed NA coupling matrix elements. If the energy levels are spaced closer than 0.1 eV, one obtains good results while including only three nearest-neighbor couplings, and in some cases even with just the first nearest-neighbor coupling scheme. If the energy levels are spaced by about 0.4 eV, the nearest-neighbor model fails, while three or more nearest-neighbor schemes also provide good results. In comparison, the results of NAMD simulation without decoherence vary continuously with changes in the number of NA couplings. Thus, decoherence effects induced by coupling to a quantum-mechanical environment not only provide the physical mechanism for NAMD trajectory branding and improve the accuracy of NAMD simulations, but also afford significant computational savings.

Two-photon absorption in CdSe colloidal quantum dots compared to organic molecules

We discuss fundamental differences in electronic structure as reflected in one- and two-photon absorption spectra of semiconductor quantum dots and organic molecules by performing systematic experimental and theoretical studies of the size-dependent spectra of colloidal quantum dots. Quantum-chemical and effective-mass calculations are used to model the one- and two-photon absorption spectra and compare them with the experimental results. Currently, quantum-chemical calculations are limited to only small-sized quantum dots (nanoclusters) but allow one to study various environmental effects on the optical spectra such as solvation and various surface functionalizations. The effective-mass calculations, on the other hand, are applicable to the larger-sized quantum dots and can, in general, explain the observed trends but are insensitive to solvent and ligand effects. Careful comparison of the experimental and theoretical results allows for quantifying the range of applicability of theoretical methods used in this work. Our study shows that the small clusters can be in principle described in a manner similar to that used for organic molecules. In addition, there are several important factors (quality of passivation, nature of the ligands, and intraband/interband transitions) affecting optical properties of the nanoclusters. The larger-size quantum dots, on the other hand, behave similarly to bulk semiconductors, and can be well described in terms of the effective-mass models.

Second-quantized surface hopping

The trajectory surface hopping method for quantum dynamics is reformulated in the space of many-particle states to include entanglement and correlation of trajectories. Used to describe many-body correlation effects in electronic structure theories, second quantization is applied to semiclassical trajectories. The new method allows coupling between individual trajectories via energy flow and common phase evolution. It captures the properties of a wave packet, such as branching, Heisenberg uncertainty, and decoherence. Applied to a superexchange process, the method shows very accurate results, comparable to exact quantum data and improving greatly on the standard approach.

Charge Separation and Exciton Dynamics at Polymer/ZnO Interface from First-Principles Simulations

Charge separation and exciton dynamics play a crucial role in determining the performance of excitonic photovoltaics. Using time-dependent density functional theory with a range-separated exchange-correlation functional as well as nonadiabatic ab initio molecular dynamics, we have studied the formation and dynamics of charge-transfer (CT) excitons at polymer/ZnO interface. The interfacial atomic structure, exciton density of states and conversions between exciton species are examined from first-principles. The exciton dynamics exhibits both adiabatic and nonadiabatic characters. While the adiabatic transitions are facilitated by C═C vibrations along the polymer (P3HT) backbone, the nonadiabatic transitions are realized by exciton hopping between the excited states. We find that the localized ZnO surface states lead to localized low-energy CT states and poor charge separation. In contrast, the surface states of crystalline C60 are indistinguishable from the bulk states, resulting in delocalized CT states and efficient charge separation in polymer/fullerene (P3HT/PCBM) heterojunctions. The hot CT states are found to cool down in an ultrafast time scale and may not play a major role in charge separation of P3HT/ZnO. Finally we suggest that the dimensions of nanostructured acceptors can be tuned to obtain both efficient charge separation and high open circuit voltages.