We present a semiclassical approach for nonadiabatic molecular dynamics based on the Ehrenfest method with corrections for decoherence and detailed balance. Decoherence is described via a coherence penalty functional that drives dynamics away from regions in Hilbert space characterized by large values of coherences. Detailed balance is incorporated by modification of the off-diagonal matrix elements with a quantum correction factor used in semiclassical approximations to quantum time-correlation functions. Both decoherence and detailed balance corrections introduce nonlinear terms to the Schrödinger equation. At the same time, the simplicity of fully deterministic dynamics and a single trajectory for each initial condition is preserved. In contrast, surface hopping is stochastic and requires averaging over multiple realization of the stochastic process for each initial condition. The Ehrenfest-decoherence-detailed-balance (Ehrenfest-DDB) method is adapted to the classical path approximation and ab initio time-dependent density functional theory and applied to an experimentally studied nanoscale system consisting of a fluorophore molecule and an scanning tunneling microscopy tip and undergoing current-induced charge injection, cooling, and recombination. Ehrenfest-DDB produces time scales that are similar to those obtained with decoherence induced surface hopping, which is a popular nonadiabatic molecular dynamics technique applied to condensed matter. At long times, Ehrenfest-DDB dynamics slows down considerably because the detailed balance correction makes off-diagonal elements go to zero on approach to Boltzmann equilibrium. The Ehrenfest-DDB technique provides efficient means to study quantum dynamics in large systems.

Recent experiments on CdSe nanoplatelets synthesized with precisely controlled thickness that eliminates ensemble disorder have allowed accurate measurement of quantum coherence at room temperature. Matching exactly the CdSe cores of the experimentally studied particles and considering several defects, we establish the atomistic origins of the loss of coherence between heavy and light hole excitations in two-dimensional CdSe and CdSe/CdZnS core/shell structures. The coherence times obtained using molecular dynamics based on tight-binding density functional theory are in excellent agreement with the measured values. We show that a long coherence time is a consequence of both small fluctuations in the energy gap between the excited state pair, which is much less than thermal energy, and a slow decay of correlation between the energies of the two states. Anionic defects at the core/shell interface have little effect on the coherence lifetime, while cationic defects strongly perturb the electronic structure, destroying the experimentally observed coherence. By coupling to the same phonon modes, the heavy and light holes synchronize their energy fluctuations, facilitating long-lived coherence. We further demonstrate that the electronic excitations are localized close to the surface of these narrow nanoscale systems, and therefore, they couple most strongly to surface acoustic phonons. The established features of electron-phonon coupling and the influence of defects, surfaces, and core/shell interfaces provide important insights into quantum coherence in nanoscale materials in general.

The microscopic origin and timescale of the fluctuations of the energies of electronic states has a significant impact on the properties of interest of electronic materials, with implication in fields ranging from photovoltaic devices to quantum information processing. Spectroscopic investigations of coherent dynamics provide a direct measurement of electronic fluctuations. Modern multidimensional spectroscopy techniques allow the mapping of coherent processes along multiple time or frequency axes and thus allow unprecedented discrimination between different sources of electronic dephasing. Exploiting modern abilities in coherence mapping in both amplitude and phase, we unravel dissipative processes of electronic coherences in the model system of CdSe quantum dots (QDs). The method allows the assignment of the nature of the observed coherence as vibrational or electronic. The expected coherence maps are obtained for the coherent longitudinal optical (LO) phonon, which serves as an internal standard and confirms the sensitivity of the technique. Fast dephasing is observed between the first two exciton states, despite their shared electron state and common environment. This result is contrary to predictions of the standard effective mass model for these materials, in which the exciton levels are strongly correlated through a common size dependence. In contrast, the experiment is in agreement with ab initio molecular dynamics of a single QD. Electronic dephasing in these materials is thus dominated by the realistic electronic structure arising from fluctuations at the atomic level rather than static size distribution. The analysis of electronic dephasing thereby uniquely enables the study of electronic fluctuations in complex materials.

Motivated by the need to study nonequilibrium evolutions of many-electron systems at the atomistic ab initio level, as they occur in modern devices and applications, we developed a quantum dynamics approach bridging master equations and surface hopping (SH). The Lindblad master equation (LME) allows us to propagate efficiently ensembles of particles, while SH provides nonperturbative evaluation of transition rates that evolve in time and depend explicitly on nuclear geometry. We implemented the LME-SH technique within real-time timedependent density functional theory using global flux SH, and we demonstrated its efficiency and utility by modeling metallic films, in which charge-phonon dynamics was studied experimentally and showed an unexpectedly strong dependence on adhesion layers. The LME-SH approach provides a general framework for modeling efficiently quantum dynamics in a broad range of complex many-electron condensed-matter and nanoscale systems.

Simulation of electronic dynamics in realistically large molecular systems is a demanding task that has not yet achieved the same level of quantitative prediction already realized for its static counterpart. This is particularly true for processes occurring beyond the Born–Oppenheimer regime. Non-adiabatic molecular dynamics (NAMD) simulations suffer from two convoluted sources of error: numerical algorithms for dynamics and electronic structure calculations. While the former has gained increasing attention, particularly addressing the validity of ad hoc methodologies, the effect of the latter remains relatively unexplored. Indeed, the required accuracy for electronic structure calculations to reach quantitative agreement with experiment in dynamics may be even more strict than that required for static simulations. Here, we address this issue by modeling the electronic energy transfer in a donor–acceptor–donor (D–A–D) molecular light harvesting system using fewest switches surface hopping NAMD simulations. In the studied system, time-resolved experimental measurements deliver complete information on spectra and energy transfer rates. Subsequent modeling shows that the calculated electronic transition energies are “sufficiently good” to reproduce experimental spectra but produce over an order of magnitude error in simulated dynamical rates. We further perform simulations using artificially shifted energy gaps to investigate the complex relationship between transition energies and modeled dynamics to understand factors affecting non-radiative relaxation and energy transfer rates.