Decoherence ensures convergence of non-adiabatic molecular dynamics with number of states

Non-adiabatic (NA) molecular dynamics (MD) is a powerful approach for studying far-from-equilibrium quantum dynamics in photophysical and photochemical systems. Most NA-MD methods are developed and tested with few-state models, and their validity with complex systems involving many states is not well studied. By modeling intraband equilibration and interband recombination of charge carriers in MoS2, we investigate the convergence of three popular NA-MD algorithms, fewest switches surface hopping (FSSH), global flux surface hopping (GFSH), and decoherence induced surface hopping (DISH) with the number of states. Only the standard DISH algorithm converges with the number of states and produces Boltzmann equilibrium. Unitary propagation of the wave function in FSSH and GFSH violates the Boltzmann distribution, leads to internal inconsistency between time-dependent Schrödinger equation state populations and trajectory counts, and produces non-convergent results. Introducing decoherence in FSSH and GFSH by collapsing the wave function fixes these problems. The simplified version of DISH that omits projecting out the occupied state and is applicable to few-state systems also causes problems when the number of states is increased. We discuss the algorithmic application of wave function collapse and Boltzmann detailed balance and provide detailed FSSH, GFSH, and DISH flow charts. The use of convergent NA-MD methods is highly important for modeling complicated quantum processes involving multiple states. Our findings provide the basis for investigating quantum dynamics in realistic complex systems.

Quantifying Ultrafast Energy Transfer from Plasmonic Hot Carriers for Pulsed Photocatalysis on Nanostructures

Photocatalysis with plasmonic nanostructures has lately emerged as a transformative paradigm to drive and alter chemical reactions using light. At the surface of metallic nanoparticles, photoexcitation results in strong near fields, short-lived high-energy "hot" carriers, and light-induced heating, thus creating a local environment where reactions can occur with enhanced efficiencies. In this context, it is critical to understand how to manipulate the nonequilibrium processes triggered by light, as their ultrafast (femto- to picoseconds) relaxation dynamics compete with the process of energy transfer toward the reactants. Accurate predictions of the plasmon photocatalytic activity can lead to optimized nanophotonic architectures with enhanced selectivity and rates, operating beyond the intrinsic limitations of the steady state. Here, we report on an original modeling approach to quantify, with space, time, and energy resolution, the ultrafast energy exchange from plasmonic hot carriers (HCs) to molecular systems adsorbed on the metal nanoparticle surface while consistently accounting for photothermal bond activation. Our analysis, illustrated for a few typical cases, reveals that the most energetic nonequilibrium carriers (i.e., with energies well far from the Fermi level) may introduce a wavelength-dependence of the reaction rates, and it elucidates on the role of the carriers closer to the Fermi energy and the photothermally heated lattice, suggesting ways to enhance and optimize each contribution. We show that the overall reaction rates can benefit strongly from using pulsed illumination with the optimal pulse width determined by the properties of the system. Taken together, these results contribute to the rational design of nanoreactors for pulsed catalysis, which calls for predictive modeling of the ultrafast HC-hot adsorbate energy transfer.

Detailed Complementary Consistency: Wave Function Tells Particle How to Hop, Particle Tells Wave Function How to Collapse

In mixed quantum-classical dynamics, the quantum subsystem can have both wave function and particle-like descriptions. However, they may yield inconsistent results for the expectation value of the same physical quantity. We here propose a novel detailed complementary consistency (DCC) method based on the principle of detailed internal consistency. Namely, the wave function along each trajectory tells the particle how to hop, while the particle tells the wave function how to collapse based on active states in the trajectory ensemble. As benchmarked in a diverse array of representative models with localized nonadiabatic couplings, DCC not only achieves fully consistent results (i.e., identical populations calculated based on wave functions and active states) but also closely reproduces the exact quantum results. Due to the high performance, our new DCC method has great potential to give a consistent and accurate mixed quantum-classical description of general nonadiabatic dynamics after further development.

General framework for quantifying dissipation pathways in open quantum systems. I. Theoretical formulation

We present a general and practical theoretical framework to investigate how energy is dissipated in open quantum system dynamics. This is performed by quantifying the contributions of individual bath components to the overall dissipation of the system. The framework is based on the Nakajima-Zwanzig projection operator technique, which allows us to express the rate of energy dissipation into a specific bath degree of freedom by using traces of operator products. The approach captures system-bath interactions to all orders, but is based on second-order perturbation theory on the off-diagonal subsystem's couplings and a Markovian description of the bath. The usefulness of our theory is demonstrated by applying it to various models of open quantum systems involving harmonic oscillators or spin baths and connecting the outcomes to existing results such as our previously reported formula derived for locally coupled harmonic baths [Kim and Franco, J. Chem. Phys. 154, 084109 (2021)]. We also prove that the dissipation calculated by our theory rigorously satisfies thermodynamic principles such as energy conservation and detailed balance. Overall, the strategy can be used to develop the theory and simulation of dissipation pathways to interpret and engineer the dynamics of open quantum systems.

Nuclear Quantum Effects on Nonadiabatic Dynamics of a Green Fluorescent Protein Chromophore Analogue: Ring-Polymer Surface-Hopping Simulation

Herein, we have used the "on-the-fly" ring-polymer surface-hopping simulation method with the centroid approximation (RPSH-CA), in combination with the multireference OM2/MRCI electronic structure calculations to study the photoinduced dynamics of a green fluorescent protein (GFP) chromophore analogue in the gas phase, i.e., o-HBI, at 50, 100, and 300 K with 1, 5, 10, and 15 beads (3600 1 ps trajectories). The electronic structure calculations identified five new minimum-energy conical intersection (MECI) structures, which, together with the previous one, play crucial roles in the excited-state decay dynamics of o-HBI. It is also found that the excited-state intramolecular proton transfer (ESIPT) occurs in an ultrafast manner and is completed within 20 fs in all the simulation conditions because there is no barrier associated with this ESIPT process in the S1 state. However, the other excited-state dynamical results are strongly related to the number of beads. At 50 and 100 K, the nuclear quantum effects (NQEs) are very important; therefore, the excited-state dynamical results change significantly with the bead number. For example, the S1 decay time deduced from time-dependent state populations becomes longer as the bead number increases. Nevertheless, an essentially convergent trend is observed when the bead number is close to 10. In contrast, at 300 K, the NQEs become weaker and the above dynamical results converge very quickly even with 1 bead. Most importantly, the NQEs seriously affect the excited-state decay mechanism of o-HBI. At 50 and 100 K, most trajectories decay to the S0 state via perpendicular keto MECIs, whereas, at 300 K, only twisted keto MECIs are responsible for the excited-state decay. The present work not only comprehensively explores the temperature-dependent photoinduced dynamics of o-HBI, but also demonstrates the importance and necessity of NQEs in nonadiabatic dynamics simulations, especially at relatively low temperatures.

Prediction challenge: First principles simulation of the ultrafast electron diffraction spectrum of cyclobutanone

Computer simulation has long been an essential partner of ultrafast experiments, allowing the assignment of microscopic mechanistic detail to low-dimensional spectroscopic data. However, the ability of theory to make a priori predictions of ultrafast experimental results is relatively untested. Herein, as a part of a community challenge, we attempt to predict the signal of an upcoming ultrafast photochemical experiment using state-of-the-art theory in the context of preexisting experimental data. Specifically, we employ ab initio Ehrenfest with collapse to a block mixed quantum-classical simulations to describe the real-time evolution of the electrons and nuclei of cyclobutanone following excitation to the 3s Rydberg state. The gas-phase ultrafast electron diffraction (GUED) signal is simulated for direct comparison to an upcoming experiment at the Stanford Linear Accelerator Laboratory. Following initial ring-opening, dissociation via two distinct channels is observed: the C3 dissociation channel, producing cyclopropane and CO, and the C2 channel, producing CH2CO and C2H4. Direct calculations of the GUED signal indicate how the ring-opened intermediate, the C2 products, and the C3 products can be discriminated in the GUED signal. We also report an a priori analysis of anticipated errors in our predictions: without knowledge of the experimental result, which features of the spectrum do we feel confident we have predicted correctly, and which might we have wrong?

Nonadiabatic Dynamics Simulations for Photoinduced Processes in Molecules and Semiconductors: Methodologies and Applications

Nonadiabatic dynamics (NAMD) simulations have become powerful tools for elucidating complicated photoinduced processes in various systems from molecules to semiconductor materials. In this review, we present an overview of our recent research on photophysics of molecular systems and periodic semiconductor materials with the aid of ab initio NAMD simulation methods implemented in the generalized trajectory surface-hopping (GTSH) package. Both theoretical backgrounds and applications of the developed NAMD methods are presented in detail. For molecular systems, the linear-response time-dependent density functional theory (LR-TDDFT) method is primarily used to model electronic structures in NAMD simulations owing to its balanced efficiency and accuracy. Moreover, the efficient algorithms for calculating nonadiabatic coupling terms (NACTs) and spin-orbit couplings (SOCs) have been coded into the package to increase the simulation efficiency. In combination with various analysis techniques, we can explore the mechanistic details of the photoinduced dynamics of a range of molecular systems, including charge separation and energy transfer processes in organic donor-acceptor structures, ultrafast intersystem crossing (ISC) processes in transition metal complexes (TMCs), and exciton dynamics in molecular aggregates. For semiconductor materials, we developed the NAMD methods for simulating the photoinduced carrier dynamics within the framework of the Kohn-Sham density functional theory (KS-DFT), in which SOC effects are explicitly accounted for using the two-component, noncollinear DFT method. Using this method, we have investigated the photoinduced carrier dynamics at the interface of a variety of van der Waals (vdW) heterojunctions, such as two-dimensional transition metal dichalcogenides (TMDs), carbon nanotubes (CNTs), and perovskites-related systems. Recently, we extended the LR-TDDFT-based NAMD method for semiconductor materials, allowing us to study the excitonic effects in the photoinduced energy transfer process. These results demonstrate that the NAMD simulations are powerful tools for exploring the photodynamics of molecular systems and semiconductor materials. In future studies, the NAMD simulation methods can be employed to elucidate experimental phenomena and reveal microscopic details as well as rationally design novel photofunctional materials with desired properties.

Learning Decoherence Time Formulas for Surface Hopping from Quantum Dynamics

Surface hopping simulations have achieved great success in many different fields, but their reliability has long been limited by the overcoherence problem. We here present a general machine learning assisted approach to identify optimal decoherence time formulas for surface hopping using exact quantum dynamics as references. In order to avoid computationally expensive force calculations, we use the nuclear kinetic energy and the adiabatic energy difference to iteratively generate the descriptor space. Through multilayer screening of the candidate descriptors and discrete optimization of the relevant parameters, we obtain new energy-based decoherence time formulas. As benchmarked in thousands of diverse multilevel systems and six standard scattering models, surface hopping with our new decoherence time formulas nearly reproduces the exact quantum dynamics while maintaining high efficiency. Thereby, our approach provides a promising avenue for systematically improving the accuracy of surface hopping simulations in complex systems from quantum dynamics data.

The landscape of computational approaches for artificial photosynthesis

Artificial photosynthesis is an attractive strategy for converting solar energy into fuels, largely because the Earth receives enough solar energy in one hour to meet humanity's energy needs for an entire year. However, developing devices for artificial photosynthesis remains difficult and requires computational approaches to guide and assist the interpretation of experiments. In this Perspective, we discuss current and future computational approaches, as well as the challenges of designing and characterizing molecular assemblies that absorb solar light, transfer electrons between interfaces, and catalyze water-splitting and fuel-forming reactions.

Two-Dimensional Octuple-Atomic-Layer M2Si2N4 (M = Al, Ga and In) with Long Carrier Lifetime

Bulk III-nitride materials MN (M = Al, Ga and In) and their alloys have been widely used in high-power electronic and optoelectronic devices, but stable two-dimensional (2D) III-nitride materials, except h-BN, have not been realized yet. A new kind of 2D III-nitride material M2Si2N4 (M = Al, Ga and In) is predicted by choosing Si as the appropriate passivation element. The stability, electronic and optical properties of 2D M2Si2N4 materials are studied systematically based on first-principles calculations. The results show that Al2Si2N4 and Ga2Si2N4 are found to be indirect bandgap semiconductors, while In2Si2N4 is a direct bandgap semiconductor. Moreover, Al2Si2N4 and In2Si2N4 have good absorption ability in the visible light region, while Ga2Si2N4 is an ultraviolet-light-absorbing material. Furthermore, the carrier lifetimes of Ga2Si2N4 and In2Si2N4 are as large as 157.89 and 103.99 ns, respectively. All these desirable properties of M2Si2N4 materials make them attractive for applications in electronics and photoelectronics.

Non-negligible roles of charge transfer excitons in ultrafast excitation energy transfer dynamics of a double-walled carbon nanotube

Herein, we employed a developed linear response time dependent density functional theory-based nonadiabatic dynamics simulation method that explicitly takes into account the excitonic effects to investigate photoinduced excitation energy transfer dynamics of a double-walled carbon nanotube (CNT) model with different excitation energies. The E₁₁ excitation of the outer CNT will generate a local excitation (LE) |out*〉 exciton due to its low energy, which does not induce any charge separation. In contrast, the E₁₁ excitation of the inner CNT can generate four kinds of excitons with the LE exciton |in*〉 dominates. In the 500-fs dynamics simulation, the LE exciton |in*〉 and charge transfer (CT) excitons |out^-in⁺〉 and |out⁺in^-〉 are all gradually converted to the |out*〉 exciton, corresponding to a photoinduced excitation energy transfer, which is consistent with experimental studies. Finally, when the excitation energy is close to the E₂₂ state of the outer CNT (∼1.05 eV), a mixed population of different excitons, with the |out*〉 exciton dominated, is generated. Then, photoinduced energy transfer from the outer to inner CNTs occurs in the first 50 fs, which is followed by an inner to outer excitation energy transfer that is completed in 400 fs. The present work not only sheds important light on the mechanistic details of wavelength-dependent excitation energy transfer of a double-walled CNT model but also demonstrates the roles and importance of CT excitons in photoinduced excitation energy transfer. It also emphasized that explicitly including the excitonic effects in electronic structure calculations and nonadiabatic dynamics simulations is significant for correct understanding/rational design of optoelectronic properties of periodically extended systems.

Exciton Dissociation in a Model Organic Interface: Excitonic State-Based Surface Hopping versus Multiconfigurational Time-Dependent Hartree

Quantum dynamical simulations are essential for a molecular-level understanding of light-induced processes in optoelectronic materials, but they tend to be computationally demanding. We introduce an efficient mixed quantum-classical nonadiabatic molecular dynamics method termed eXcitonic state-based Surface Hopping (X-SH), which propagates the electronic Schrödinger equation in the space of local excitonic and charge-transfer electronic states, coupled to the thermal motion of the nuclear degrees of freedom. The method is applied to exciton decay in a 1D model of a fullerene-oligothiophene junction, and the results are compared to the ones from a fully quantum dynamical treatment at the level of the Multilayer Multiconfigurational Time-Dependent Hartree (ML-MCTDH) approach. Both methods predict that charge-separated states are formed on the 10-100 fs time scale via multiple "hot-exciton dissociation" pathways. The results demonstrate that X-SH is a promising tool advancing the simulation of photoexcited processes from the molecular to the true nanomaterials scale.

Quantum dynamical study of inter-chain exciton transport in a regioregular P3HT model system at finite temperature: HJ vs. H-aggregate models

We report on quantum dynamical simulations of inter-chain exciton transport in a model of regioregular poly(3-hexylthiophene), rr-P3HT, at finite temperature, using the Multi-Layer Multi-Configuration Time-Dependent Hartree (ML-MCTDH) method for a system of up to 63 electronic states and 180 vibrational modes. A Frenkel Hamiltonian of HJ aggregate type is used, along with a reduced H-aggregate representation; electron-phonon coupling includes local high-frequency modes as well as anharmonic intermolecular modes. The latter are operative in mediating inter-chain transport, by a mechanism of transient localization type. Strikingly, this mechanism is found to be of quantum coherent character and involves non-adiabatic effects. Using periodic boundary conditions, a normal diffusion regime is identified from the exciton mean-squared displacement, apart from early-time transients. Diffusion coefficients are found to be of the order of 3 x 10-3 cm2/s, showing a non-monotonous increase with temperature.

Simulations of molecular photodynamics in long timescales

Nonadiabatic dynamics simulations in the long timescale (much longer than 10 ps) are the next challenge in computational photochemistry. This paper delimits the scope of what we expect from methods to run such simulations: they should work in full nuclear dimensionality, be general enough to tackle any type of molecule and not require unrealistic computational resources. We examine the main methodological challenges we should venture to advance the field, including the computational costs of the electronic structure calculations, stability of the integration methods, accuracy of the nonadiabatic dynamics algorithms and software optimization. Based on simulations designed to shed light on each of these issues, we show how machine learning may be a crucial element for long time-scale dynamics, either as a surrogate for electronic structure calculations or aiding the parameterization of model Hamiltonians. We show that conventional methods for integrating classical equations should be adequate to extended simulations up to 1 ns and that surface hopping agrees semiquantitatively with wave packet propagation in the weak-coupling regime. We also describe our optimization of the Newton-X program to reduce computational overheads in data processing and storage. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.

Multilayer Subsystem Surface Hopping Method for Large-Scale Nonadiabatic Dynamics Simulation with Hundreds of Thousands of States

We present a multilayer subsystem surface hopping (MSSH) method to deal with nonadiabatic dynamics in large-scale systems. A small subsystem instead of the full system is adopted for surface hopping and is updated on-the-fly to achieve a reliable description of important adiabatic states and the wave function evolution. Additional subsystems for molecular dynamics and statistical description are introduced to further improve the simulation reliability. The global flux hopping probabilities with optimal state assignments are utilized to treat the complex surface crossings. As demonstrated in a series of one- and two-dimensional Holstein models with up to hundreds of thousands of states, MSSH shows weak parameter dependence in all investigated systems. Especially, the computational costs are reduced by 2-6 orders of magnitude compared to traditional surface hopping simulations in full systems, and size-independent results are achieved with a large time-step size of 2-5 fs. The new method is compatible with different decoherence correction strategies and achieves a much better balance between efficiency and reliability, thus promising for applications in general charge and exciton dynamics simulations.

Generalized Ab Initio Nonadiabatic Dynamics Simulation Methods from Molecular to Extended Systems

Nonadiabatic dynamics simulation has become a powerful tool to describe nonadiabatic effects involved in photophysical processes and photochemical reactions. In the past decade, our group has developed generalized trajectory-based ab initio surface-hopping (GTSH) dynamics simulation methods, which can be used to describe a series of nonadiabatic processes, such as internal conversion, intersystem crossing, excitation energy transfer and charge transfer of molecular systems, and photoinduced nonadiabatic carrier dynamics of extended systems with and without spin-orbit couplings. In this contribution, we will first give a brief introduction to our recently developed methods and related numerical implementations at different computational levels. Later, we will present some of our latest applications in realistic systems, which cover organic molecules, biological proteins, organometallic compounds, periodic organic and inorganic materials, etc. Final discussion is given to challenges and outlooks of ab initio nonadiabatic dynamics simulations.

Effects of oxygen vacancies on the photoexcited carrier lifetime in rutile TiO2

The photoexcited carrier lifetime in semiconductors plays a crucial role in solar energy conversion processes. The defects or impurities in semiconductors are usually proposed to introduce electron-hole (e-h) recombination centers and consequently reduce the photoexcited carrier lifetime. In this report, we investigate the effects of oxygen vacancies (O_V) on the carrier lifetime in rutile TiO₂, which has important applications in photocatalysis and photovoltaics. It is found that an O_V introduces two excess electrons which form two defect states in the band gap. The lower state is localized on one Ti atom and behaves as a small polaron, and the higher one is a hybrid state contributed by three Ti atoms around the O_V. Both the polaron and hybrid states exhibit strong electron-phonon (e-ph) coupling and their charge distributions become more and more delocalized when the temperature increases from 100 to 700 K. Such strong e-ph coupling and charge delocalization enhance the nonadibatic coupling between the electronic states along the hole relaxation path, where the defect states behave as intermediate states, leading to a distinct acceleration of e-h recombination. Our study provides valuable insights to understand the role of defects on photoexcited carrier lifetime in semiconductors.

Excited State Dynamics in Dual-Defects Modified Graphitic Carbon Nitride

Significant efforts are focused on defect-engineering of metal-free graphitic carbon nitride (g-C₃N₄) to amplify its efficacy. A conceptually new multidefect-modified g-C₃N₄ having simultaneously two or more defects has attracted strong attention for its enhanced photocatalytic properties. We model and compare the excited state dynamics in g-C₃N₄ with (i) nitrogen defects (N vacancy and CN group) and (ii) dual defects (N vacancy, CN group, and O doping) and show that the nonradiative recombination of charge carriers in these systems follows the Shockley-Read-Hall mechanism. The nitrogen defects create three midgap states that trap charges and act as recombination centers. The dual-defect modified systems exhibit superior properties compared with pristine g-C₃N₄ because the defects facilitate rapid charge separation and extend the spectrum of absorbed light. The system doped with O shows better performance due to enhanced carrier lifetime and higher oxidation potential caused by a downshifted valence band. The study provides guidance for rational design of stable and efficient photocatalytic materials.

All-Atom Nonadiabatic Dynamics Simulation of Hybrid Graphene Nanoribbons Based on Wannier Analysis and Machine Learning

Trajectory surface hopping combined with ab initio electronic structure calculations is a popular and powerful approach for on-the-fly nonadiabatic dynamics simulations. For large systems, however, this remains a significant challenge because of the unaffordable computational cost of large-scale electronic structure calculations. Here, we present an efficient divide-and-conquer approach to construct the system Hamiltonian based on Wannier analysis and machine learning. In detail, the large system under investigation is first decomposed into small building blocks, and then all possible segments formed by building blocks within a cutoff distance are found out. Ab initio molecular dynamics is carried out to generate a sequence of geometries for each equivalent segment with periodicity. The Hamiltonian matrices in the maximum localized Wannier function (MLWF) basis are obtained for all geometries and utilized to train artificial neural networks (ANNs) for the structure-dependent Hamiltonian elements. Taking advantage of the orthogonality and spatial locality of MLWFs, the one-electron Hamiltonian of a large system at arbitrary geometry can be directly constructed by the trained ANNs. As demonstrations, we study charge transport in a zigzag graphene nanoribbon (GNR), a coved GNR, and a series of hybrid GNRs with a state-of-the-art surface hopping method. The interplay between delocalized and localized states is found to determine the electron dynamics in hybrid GNRs. Our approach has successfully studied GNRs with >10 000 atoms, paving the way for efficient and reliable all-atom nonadiabatic dynamics simulation of general systems.

Extending the Time Scales of Nonadiabatic Molecular Dynamics via Machine Learning in the Time Domain

A novel methodology for direct modeling of long-time scale nonadiabatic dynamics in extended nanoscale and solid-state systems is developed. The presented approach enables forecasting the vibronic Hamiltonians as a direct function of time via machine-learning models trained directly in the time domain. The use of periodic and aperiodic functions that transform time into effective input modes of the artificial neural network is demonstrated to be essential for such an approach to work for both abstract and atomistic models. The best strategies and possible limitations pertaining to the new methodology are explored and discussed. An exemplary direct simulation of unprecedentedly long 20 picosecond trajectories is conducted for a divacancy-containing monolayer black phosphorus system, and the importance of conducting such extended simulations is demonstrated. New insights into the excited states photophysics in this system are presented, including the role of decoherence and model definition.

Conditional Wave Function Theory: A Unified Treatment of Molecular Structure and Nonadiabatic Dynamics

We demonstrate that a conditional wave function theory enables a unified and efficient treatment of the equilibrium structure and nonadiabatic dynamics of correlated electron-ion systems. The conditional decomposition of the many-body wave function formally recasts the full interacting wave function of a closed system as a set of lower-dimensional (conditional) coupled "slices". We formulate a variational wave function ansatz based on a set of conditional wave function slices and demonstrate its accuracy by determining the structural and time-dependent response properties of the hydrogen molecule. We then extend this approach to include time-dependent conditional wave functions and address paradigmatic nonequilibrium processes including strong-field molecular ionization, laser-driven proton transfer, and nuclear quantum effects induced by a conical intersection. This work paves the road for the application of conditional wave function theory in equilibrium and out-of-equilibrium ab initio molecular simulations of finite and extended systems.

An accurate, non-empirical method for incorporating decoherence into Ehrenfest dynamics

In mixed quantum-classical nonadiabatic molecular dynamics methods, the anchoring of the electronic wave function to a single nuclear geometry results in both quantitative and qualitative errors in the dynamics. In the context of both Ehrenfest and trajectory surface hopping methods, methods for incorporating decoherence are widely used to eliminate these errors. However, the accuracy of these methods often depends strongly on the parameterization of the decoherence time and/or other related quantities. Here, we present a refinement of the recently introduced collapse to a block (TAB) scheme for incorporating decoherence into Ehrenfest dynamics. The proposed approach incorporates an approximation to the history of the population dynamics and treats the coherence decay as Gaussian, rather than exponential. This method uses parameters that can be obtained from first principles, rather than empirical fitting. Application to one-dimensional models indicates excellent agreement with numerically exact simulations. We also introduce a second refinement to the TAB method: a robust linear least-squares algorithm for determining collapse probabilities.

An efficient implementation of spin-orbit coupling within the framework of semiempirical orthogonalization-corrected methods for ultrafast intersystem crossing dynamics

We implement spin-orbit coupling (SOC) within the framework of semiempirical orthogonalization-corrected methods (OMx). The excited-state wavefunction is generated from configuration interaction with single excitations (CIS). The SOC Hamiltonian in terms of the one-electron Breit-Pauli operator with effective nuclear charges is adopted in this work. Benchmark calculations show that SOCs evaluated using the OMx/CIS method agree very well with those obtained from time-dependent density functional theory. As a particularly attractive application, we incorporate SOCs between singlet and triplet states into Tully's fewest switches surface hopping algorithm to enable excited-state nonadiabatic dynamics simulations, treating internal conversion and intersystem crossing on an equal footing. This semiempirical dynamics simulation approach is applied to investigate ultrafast intersystem crossing processes in core-substituted naphthalenediimides.

Dependence between Structural and Electronic Properties of CsPbI3: Unsupervised Machine Learning of Nonadiabatic Molecular Dynamics

Using unsupervised machine learning on the trajectories from a nonadiabatic molecular dynamics simulation with time-dependent Kohn-Sham density functional theory, we elucidated the structural parameters with the largest influence on nonradiative recombination of charge carriers in CsPbI₃, which forms the basis for solar energy and optoelectronic applications. The I-I-I angles between PbI₆ octahedra, followed by the Cs-I distance, have the strongest impact on the bandgap and the nonadiabatic coupling. The importance of the Cs-I distance is unexpected, because Cs does not contribute to electron and hole wave functions. The nonadiabatic coupling is most influenced by static properties, which is also surprising, given its explicit dependence on atomic velocities.

Interfacial photoinduced carrier dynamics tuned by polymerization of coronene molecules encapsulated in carbon nanotubes: bridging type-I and type-II heterojunctions

Carbon nanomaterials play important roles in modern scientific research. Integrating different carbon-based building blocks into nano-hybrid architectures not only takes full advantage of each component, but also brings in novel interfacial properties. Herein, we have employed density functional theory (DFT) calculations to investigate the effects of polymerization degree of coronene molecules encapsulated in single-walled carbon nanotubes (SWNTs) (19,0) on their interfacial properties. The present results reveal that the interfacial properties of the formed heterojunctions are remarkably regulated by the polymerization degree. For example, monomer- and dimer-encapsulated SWNTs are type-I heterojunctions in which interfacial excitation energy transfer is preferred, whereas interfacial charge carrier transfer is favorable in trimer- and polymer-encapsulated SWNTs because they are type-II heterojunctions. On the other hand, we have employed the time-domain nonadiabatic dynamics simulation approach to explore the interfacial carrier dynamics in type-II polymer-encapsulated SWNT heterojunctions. It is found that the electron and hole transfer processes are asymmetric and occur in opposite directions and at different rates. The former takes place from polymers to SWNTs in an ultrafast way (ca. 370 fs), whereas the latter occurs slowly from SWNTs to polymers (ca. 24 ps). A closer analysis uncovers the fact that the different carrier transfer rates mainly originate from the different densities of the acceptor states, energy differences and inter-state couplings between the donor and acceptor states. Finally, the present work demonstrates that the polymerization degree could act as a new regulating strategy to tune the interfacial properties of molecule-encapsulated SWNT heterojunctions.

Ab initio nonadiabatic molecular dynamics of charge carriers in metal halide perovskites

Photoinduced nonequilibrium processes in nanoscale materials play key roles in photovoltaic and photocatalytic applications. This review summarizes recent theoretical investigations of excited state dynamics in metal halide perovskites (MHPs), carried out using a state-of-the-art methodology combining nonadiabatic molecular dynamics with real-time time-dependent density functional theory. The simulations allow one to study evolution of charge carriers at the ab initio level and in the time-domain, in direct connection with time-resolved spectroscopy experiments. Eliminating the need for the common approximations, such as harmonic phonons, a choice of the reaction coordinate, weak electron-phonon coupling, a particular kinetic mechanism, and perturbative calculation of rate constants, we model full-dimensional quantum dynamics of electrons coupled to semiclassical vibrations. We study realistic aspects of material composition and structure and their influence on various nonequilibrium processes, including nonradiative trapping and relaxation of charge carriers, hot carrier cooling and luminescence, Auger-type charge-charge scattering, multiple excitons generation and recombination, charge and energy transfer between donor and acceptor materials, and charge recombination inside individual materials and across donor/acceptor interfaces. These phenomena are illustrated with representative materials and interfaces. Focus is placed on response to external perturbations, formation of point defects and their passivation, mixed stoichiometries, dopants, grain boundaries, and interfaces of MHPs with charge transport layers, and quantum confinement. In addition to bulk materials, perovskite quantum dots and 2D perovskites with different layer and spacer cation structures, edge passivation, and dielectric screening are discussed. The atomistic insights into excited state dynamics under realistic conditions provide the fundamental understanding needed for design of advanced solar energy and optoelectronic devices.

Multiscale Models for Light-Driven Processes

Multiscale models combining quantum mechanical and classical descriptions are a very popular strategy to simulate properties and processes of complex systems. Many alternative formulations have been developed, and they are now available in all of the most widely used quantum chemistry packages. Their application to the study of light-driven processes, however, is more recent, and some methodological and numerical problems have yet to be solved. This is especially the case for the polarizable formulation of these models, the recent advances in which we review here. Specifically, we identify and describe the most important specificities that the polarizable formulation introduces into both the simulation of excited-state dynamics and the modeling of excitation energy and electron transfer processes.

Quantum Dynamics of Exciton Transport and Dissociation in Multichromophoric Systems

Due to the subtle interplay of site-to-site electronic couplings, exciton delocalization, nonadiabatic effects, and vibronic couplings, quantum dynamical studies are needed to elucidate the details of ultrafast photoinduced energy and charge transfer events in organic multichromophoric systems. In this vein, we review an approach that combines first-principles parameterized lattice Hamiltonians with accurate quantum dynamical simulations using advanced multiconfigurational methods. Focusing on the elementary transfer steps in organic functional materials, we address coherent exciton migration and creation of charge transfer excitons in homopolymers, notably representative of the poly(3-hexylthiophene) material, as well as exciton dissociation at polymer:fullerene heterojunctions. We emphasize the role of coherent transfer, trapping effects due to high-frequency phonon modes, and thermal activation due to low-frequency soft modes that drive a diffusive dynamics.

Quantum dynamical simulations of intra-chain exciton diffusion in an oligo (para-phenylene vinylene) chain at finite temperature

We report on quantum dynamical simulations of exciton diffusion in an oligo(para-phenylene vinylene) chain segment with 20 repeat units (OPV-20) at finite temperature, complementary to our recent study of the same system at T = 0 K [R. Binder and I. Burghardt, J. Chem. Phys. 152, 204120 (2020)]. Accurate quantum dynamical simulations are performed using the multi-layer multi-configuration time-dependent Hartree method as applied to a site-based Hamiltonian comprising 20 electronic states of Frenkel type and 460 vibrational modes, including site-local quinoid-distortion modes along with site-correlated bond-length alternation (BLA) modes, ring torsional modes, and an explicit harmonic-oscillator bath. A first-principles parameterized Frenkel-Holstein type Hamiltonian is employed, which accounts for correlations between the ring torsional modes and the anharmonically coupled BLA coordinates located at the same junction. Thermally induced fluctuations of the torsional modes are described by a stochastic mean-field approach, and their impact on the excitonic motion is characterized in terms of the exciton mean-squared displacement. A normal diffusion regime is observed under periodic boundary conditions, apart from transient localization features. Even though the polaronic exciton species are comparatively weakly bound, exciton diffusion is found to be a coherent-rather than hopping type-process, driven by the fluctuations of the soft torsional modes. Similar to the previous observations for oligothiophenes, the evolution for the most part exhibits a near-adiabatic dynamics of local exciton ground states (LEGSs) that adjust to the local conformational dynamics. However, a second mechanism, involving resonant transitions between neighboring LEGSs, gains importance at higher temperatures.

Branching Corrected Mean Field Method for Nonadiabatic Dynamics

When nonadiabatic dynamics are described on the basis of trajectories, severe trajectory branching occurs when the nuclear wave packets on some potential energy surfaces are reflected while those on the remaining surfaces are not. As a result, the traditional Ehrenfest mean field (EMF) approximation breaks down. In this study, two versions of the branching corrected mean field (BCMF) method are proposed. Namely, when trajectory branching is identified, BCMF stochastically selects either the reflected or the nonreflected group to build the new mean field trajectory or splits the mean field trajectory into two new trajectories with the corresponding weights. As benchmarked in six standard model systems and an extensive model base with two hundred diverse scattering models, BCMF significantly improves the accuracy while retaining the high efficiency of the traditional EMF. In fact, BCMF closely reproduces the exact quantum dynamics in all investigated systems, thus highlighting the essential role of branching correction in nonadiabatic dynamics simulations of general systems.

First-Principles Quantum and Quantum-Classical Simulations of Exciton Diffusion in Semiconducting Polymer Chains at Finite Temperature

We report on first-principles quantum-dynamical and quantum-classical simulations of photoinduced exciton dynamics in oligothiophene chain segments, representative of intrachain exciton migration in the poly(3-hexylthiophene) (P3HT) polymer. Following up on our recent study (Binder R.; Burghardt, I. Faraday Discuss. 2020, 221, 406), multilayer multiconfiguration time-dependent Hartree calculations for a short oligothiophene segment comprising 20 monomer units (OT-20) are carried out to obtain full quantum-dynamical simulations at finite temperature. These are employed to benchmark mean-field Ehrenfest calculations, which are shown to give qualitatively correct results for the present system. Periodic boundary conditions turn out to significantly improve earlier estimates of diffusion coefficients. Using the Ehrenfest approach, a series of calculations are subsequently carried out for larger lattices (OT-40 to OT-80), leading to estimates for temperature-dependent mean-squared displacements, which are found to exhibit a near-linear dependence as a function of time. The resulting diffusion coefficient estimates are an increasing function of temperature, whose detailed functional form depends on the degree of static disorder. With a realistic static disorder parameter (σ_s ≃ 0.06 eV), the diffusion coefficients decrease from D ∼ 1 × 10^-2 cm² s^-1 to D ∼ 1 × 10^-3 cm² s^-1, in qualitative agreement with experimental data for P3HT. The dynamical scenario obtained from our simulations shows that exciton migration in P3HT-type chains is a largely adiabatic process throughout the temperature regime we investigated (i.e., T = 50-300 K). The resulting picture of exciton migration is a coherent, but not bandlike, motion of an exciton-polaron driven by fluctuations induced by low-frequency modes. This process acquires partial hopping character if static disorder becomes prominent and Anderson localization sets in.

Long-Lived Hot Electron in a Metallic Particle for Plasmonics and Catalysis: Ab Initio Nonadiabatic Molecular Dynamics with Machine Learning

Multiple experiments provide evidence for photovoltaic, catalytic, optoelectronic, and plasmonic processes involving hot, i.e., high energy, electrons in nanoscale materials. However, the mechanisms of such processes remain elusive, because electrons rapidly lose energy by relaxation through dense manifolds of states. We demonstrate a long-lived hot electron state in a Pt nanocluster adsorbed on the MoS₂ substrate. For this purpose, we develop a simulation technique, combining classical molecular dynamics based on machine learning potentials with ab initio nonadiabatic molecular dynamics and real-time time-dependent density functional theory. Choosing Pt₂₀/MoS₂ as a prototypical system, we find frequent shifting of a top atom in the Pt particle occurring on a 50 ps time scale. The distortion breaks particle symmetry and creates unsaturated chemical bonds. The lifetime of the localized state associated with the broken bonds is enhanced by a factor of 3. Hot electrons aggregate near the shifted atom and form a catalytic reaction center. Our findings prove that distortion of even a single atom can have important implications for nanoscale catalysis and plasmonics and provide insights for utilizing machine learning potentials to accelerate ab initio investigations of excited state dynamics in condensed matter systems.

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.

State-pairwise decoherence times for nonadiabatic dynamics on more than two electronic states

Independent trajectory (IT) nonadiabatic molecular dynamics simulation methods are powerful tools for modeling processes involving transitions between electronic states. Incorporation and refinement of decoherence corrections into popular IT methods, e.g., Ehrenfest dynamics and trajectory surface hopping, is an important means of improving their accuracies. In this work, we identify a new challenge in the development of such decoherence corrections; when a system exists in a coherent superposition of three or more electronic states, coherences may decay unphysically when the decoherence correction is based on decoherence times assigned on a state-wise basis. As a solution, we introduce decoherence corrected Ehrenfest schemes based on decoherence times assigned on a state-pairwise basis. By application of these methods to a set of very simple one-dimensional model problems, we show that one of these state-pairwise methods ("collapse to a block") correctly describes the loss of coherence between all pairs of states in our multistate model problems, whereas a method based on a state-wise description of coherence loss does not. The new one-dimensional models introduced here can serve as useful tests for other decoherence correction schemes.

MAI Termination Favors Efficient Hole Extraction and Slow Charge Recombination at the MAPbI3/CuSCN Heterojunction

Photoinduced charge separation is the key step determining the efficiency of photon-to-electron conversion in solar cells, while charge carrier lifetimes govern the overall solar cell performance. Experiments report that copper(I) thiocyanate (CuSCN) is a very promising hole extraction layer for perovskite solar cells. Using nonadiabatic molecular dynamics combined with ab initio time-domain density functional theory, we show that termination of CH₃NH₃PbI₃ (MAPbI₃) at MAPbI₃/CuSCN heterojunctions has a strong influence on both charge separation and recombination. Both processes are favored by MAI termination, compared to PbI₂ termination. Because the MAPbI₃ valence band originates from iodine orbitals while the conduction band arises from Pb orbitals, MAI termination places holes close to CuSCN, favoring extraction, and creates an MAI barrier for recombination of electrons in MAPbI₃ and holes in CuSCN. The opposite is true for PbI₂ termination. The origin of these effects is attributed solely to the properties of the MAPbI₃ surfaces, and therefore, the conclusions should apply to other hole-transporting materials and can be generalized to other perovskites. Importantly, the simulations show that the injected hole remains hot for several hundreds of femtoseconds, allowing it to escape the interfacial region and prevent formation of bound excitons. This study suggests that metal halide perovskites should be treated with an organic precursor, such as MAI, prior to the formation of their interfaces with hole-transporting materials. The reported results advance the fundamental understanding of the highly unusual properties of metal halide perovskites and provide specific guidelines for optimizing the performance of perovskite solar cells and other devices.

First-principles description of intra-chain exciton migration in an oligo(para-phenylene vinylene) chain. I. Generalized Frenkel-Holstein Hamiltonian

A generalized Frenkel-Holstein Hamiltonian is constructed to describe exciton migration in oligo(para-phenylene vinylene) chains, based on excited state electronic structure data for an oligomer comprising 20 monomer units (OPV-20). Time-dependent density functional theory calculations using the ωB97XD hybrid functional are employed in conjunction with a transition density analysis to study the low-lying singlet excitations and demonstrate that these can be characterized to a good approximation as a Frenkel exciton manifold. Based on these findings, we employ the analytic mapping procedure of Binder et al. [J. Chem. Phys. 141, 014101 (2014)] to translate one-dimensional (1D) and two-dimensional (2D) potential energy surface (PES) scans to a fully anharmonic, generalized Frenkel-Holstein (FH) Hamiltonian. A 1D PES scan is carried out for intra-ring quinoid distortion modes, while 2D PES scans are performed for the anharmonically coupled inter-monomer torsional and vinylene bridge bond length alternation modes. The kinetic energy is constructed in curvilinear coordinates by an exact numerical procedure, using the TNUM Fortran code. As a result, a fully molecular-based, generalized FH Hamiltonian is obtained, which is subsequently employed for quantum exciton dynamics simulations, as shown in Paper II [R. Binder and I. Burghardt, J. Chem. Phys. 152, 204120 (2020)].

Hot Electron Cooling in Silicon Nanoclusters via Landau-Zener Nonadiabatic Molecular Dynamics: Size Dependence and Role of Surface Termination

We develop a new express methodology for modeling excited-state dynamics occurring in dense manifolds of electronic states in atomistic systems. The approach leverages a modified Landau-Zener formula, the neglect of a back-reaction approximation, and the highly efficient density functional tight-binding method. We study the hot electron dynamics in a series of H- and F-terminated silicon nanocrystals (NCs) containing up to several hundred atoms. We explain the slower electron cooling dynamics in F-terminated NCs by the larger energy gaps between the adjacent electronic states in these systems as well as their slower fluctuations. We conclude that both the mass and chemical identity of the surface termination groups equally influence the electron dynamics, on average. However, the mass effect becomes dominant for higher-energy excitations. We find that the electron decay dynamics in F-terminated NCs has a greater sensitivity to the mass of the surface ligands than do the H-terminated NCs and explain this observation by the details of the electron-phonon coupling in the systems. We find that in the H-terminated NCs, electronic transitions in the cooling process occur predominantly between the surface states, whereas in F-terminated Si NCs, both surface and NC core states are coupled to the nuclear vibrations. We find that electron energy relaxation is accelerated in larger NCs and attribute this effect to the higher densities of states and smaller energy gaps in these systems.

Reorganization energies and spectral densities for electron transfer problems in charge transport materials

Various contributions to the outer reorganization energy of an electron transfer system and their theoretical and computational aspects have been discussed.