Exact-Factorization-Based Surface Hopping without Velocity Adjustment

While surface hopping has emerged as a powerful method for simulating non-adiabatic dynamics in large molecules, the ad hoc nature of the necessary velocity adjustments and decoherence corrections in the algorithm somewhat reduces its reliability. Here we propose a new scheme that eliminates these aspects by combining the nuclear equation from the quantum-trajectory surface-hopping approach with the electronic equation derived from the exact-factorization approach. The resulting method, denoted QTSH-XF, yields a surface-hopping method on firmer ground than previous and is shown to successfully capture dynamics in Tully models and in a linear vibronic coupling model of the photoexcited uracil cation.

Nonadiabatic dynamics with classical trajectories: The problem of an initial coherent superposition of electronic states

Advances in coherent light sources and development of pump-probe techniques in recent decades have opened the way to study electronic motion in its natural time scale. When an ultrashort laser pulse interacts with a molecular target, a coherent superposition of electronic states is created and the triggered electron dynamics is coupled to the nuclear motion. A natural and computationally efficient choice to simulate this correlated dynamics is a trajectory-based method where the quantum-mechanical electronic evolution is coupled to a classical-like nuclear dynamics. These methods must approximate the initial correlated electron-nuclear state by associating an initial electronic wavefunction to each classical trajectory in the ensemble. Different possibilities exist that reproduce the initial populations of the exact molecular wavefunction when represented in a basis. We show that different choices yield different dynamics and explore the effect of this choice in Ehrenfest, surface hopping, and exact-factorization-based coupled-trajectory schemes in a one-dimensional two-electronic-state model system that can be solved numerically exactly. This work aims to clarify the problems that standard trajectory-based techniques might have when a coherent superposition of electronic states is created to initialize the dynamics, to discuss what properties and observables are affected by different choices of electronic initial conditions and to point out the importance of quantum-momentum-induced electronic transitions in coupled-trajectory schemes.

Oscillator strengths and excited-state couplings for double excitations in time-dependent density functional theory

Although useful to extract excitation energies of states of double-excitation character in time-dependent density functional theory that are missing in the adiabatic approximation, the frequency-dependent kernel derived earlier [Maitra et al., J. Chem. Phys. 120, 5932 (2004)] was not designed to yield oscillator strengths. These are required to fully determine linear absorption spectra, and they also impact excited-to-excited-state couplings that appear in dynamics simulations and other quadratic response properties. Here, we derive a modified non-adiabatic kernel that yields both accurate excitation energies and oscillator strengths for these states. We demonstrate its performance on a model two-electron system, the Be atom, and on excited-state transition dipoles in the LiH molecule at stretched bond-lengths, in all cases producing significant improvements over the traditional approximations.

Significance of Energy Conservation in Coupled-Trajectory Approaches to Nonadiabatic Dynamics

Through approximating electron-nuclear correlation terms in the exact factorization approach, trajectory-based methods have been derived and successfully applied to the dynamics of a variety of light-induced molecular processes, capturing quantum (de)coherence effects rigorously. These terms account for the coupling among the trajectories, recovering the nonlocal nature of quantum nuclear dynamics that is completely overlooked in traditional independent-trajectory algorithms. Nevertheless, some of the approximations introduced in the derivation of some of these methods do not conserve the total energy. We analyze energy conservation in the coupled-trajectory mixed quantum-classical (CTMQC) algorithm and explore the performance of a modified algorithm, CTMQC-E, where some of the terms are redefined to restore energy conservation. A set of molecular models is used as a test, namely, 2-cis-penta-2,4-dienimium cation, bis(methylene) adamantyl radical cation, butatriene cation, uracil radical cation, and neutral pyrazine.

Different flavors of exact-factorization-based mixed quantum-classical methods for multistate dynamics

The exact factorization approach has led to the development of new mixed quantum-classical methods for simulating coupled electron-ion dynamics. We compare their performance for dynamics when more than two electronic states are occupied at a given time, and analyze: (1) the use of coupled versus auxiliary trajectories in evaluating the electron-nuclear correlation terms, (2) the approximation of using these terms within surface-hopping and Ehrenfest frameworks, and (3) the relevance of the exact conditions of zero population transfer away from nonadiabatic coupling regions and total energy conservation. Dynamics through the three-state conical intersection in the uracil radical cation as well as polaritonic models in one dimension are studied.

Energy-conserving coupled trajectory mixed quantum-classical dynamics

The coupled-trajectory mixed quantum-classical method (CTMQC), derived from the exact factorization approach, has successfully predicted photo-chemical dynamics in a number of interesting molecules, capturing population transfer and decoherence from first principles. However, due to the approximations made, CTMQC does not guarantee energy conservation. We propose a modified algorithm, CTMQC-E, which redefines the integrated force in the coupled-trajectory term so to restore energy conservation, and demonstrate its accuracy on scattering in Tully's extended coupling region model and photoisomerization in a retinal chromophore model.

Curing the Divergence in Time-Dependent Density Functional Quadratic Response Theory

The adiabatic approximation in time-dependent density functional theory is known to give an incorrect pole structure in the quadratic response function, leading to unphysical divergences in excited state-to-state transition probabilities and hyperpolarizabilties. We find the form of the exact quadratic response kernel and derive a practical and accurate approximation that cures the divergence. We demonstrate our results on excited state-to-state transition probabilities of a model system and of the LiH molecule.

DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science

In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.

Exact Factorization Adventures: A Promising Approach for Non-Bound States

Modeling the dynamics of non-bound states in molecules requires an accurate description of how electronic motion affects nuclear motion and vice-versa. The exact factorization (XF) approach offers a unique perspective, in that it provides potentials that act on the nuclear subsystem or electronic subsystem, which contain the effects of the coupling to the other subsystem in an exact way. We briefly review the various applications of the XF idea in different realms, and how features of these potentials aid in the interpretation of two different laser-driven dissociation mechanisms. We present a detailed study of the different ways the coupling terms in recently-developed XF-based mixed quantum-classical approximations are evaluated, where either truly coupled trajectories, or auxiliary trajectories that mimic the coupling are used, and discuss their effect in both a surface-hopping framework as well as the rigorously-derived coupled-trajectory mixed quantum-classical approach.

Exact-Factorization-Based Surface Hopping for Multistate Dynamics

A surface-hopping algorithm recently derived from the exact factorization approach, SHXF [Ha et al. J. Phys. Chem. Lett. 2018, 9, 1097], introduces an additional term in the electronic equation of surface hopping that couples electronic states through the quantum momentum. This term not only provides a first-principles description of decoherence, but here we show it is crucial to accurately capture nonadiabatic dynamics when more than two states are occupied at any given time. Using a vibronic coupling model of the uracil cation, we show that the lack of this term in traditional surface-hopping methods, including those with decoherence corrections, leads to failure to predict the dynamics through a three-state intersection, while SHXF performs similarly to the multiconfiguration time-dependent Hartree quantum dynamics benchmark.

Double and Charge-Transfer Excitations in Time-Dependent Density Functional Theory

Time-dependent density functional theory has emerged as a method of choice for calculations of spectra and response properties in physics, chemistry, and biology, with its system-size scaling enabling computations on systems much larger than otherwise possible. While increasingly complex and interesting systems have been successfully tackled with relatively simple functional approximations, there has also been increasing awareness that these functionals tend to fail for certain classes of approximations. Here I review the fundamental challenges the approximate functionals have in describing double excitations and charge-transfer excitations, which are two of the most common impediments for the theory to be applied in a black-box way. At the same time, I describe the progress made in recent decades in developing functional approximations that give useful predictions for these excitations. Expected final online publication date for the Annual Review of Physical Chemistry, Volume 73 is April 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Minimizing the Time-Dependent Density Functional Error in Ehrenfest Dynamics

Simulating electron-ion dynamics using time-dependent density functional theory within an Ehrenfest dynamics scheme can be done in two ways that are in principle exact and identical: propagating time-dependent electronic Kohn-Sham equations or propagating electronic coefficients on surfaces obtained from linear-response. We show here that using an approximate functional leads to qualitatively different dynamics in the two approaches. We argue that the latter is more accurate because the functionals are evaluated on domains close to the ground state where currently used approximations perform better. We demonstrate this on an exactly solvable model of charge transfer and discuss implications for time-resolved spectroscopy.

Study of the Decoherence Correction Derived from the Exact Factorization Approach for Nonadiabatic Dynamics

We present a detailed study of the decoherence correction to surface hopping that was recently derived from the exact factorization approach. Ab initio multiple spawning calculations that use the same initial conditions and the same electronic structure method are used as a reference for three molecules: ethylene, the methaniminium cation, and fulvene, for which nonadiabatic dynamics follows a photoexcitation. A comparison with the Granucci–Persico energy-based decoherence correction and the augmented fewest-switches surface-hopping scheme shows that the three decoherence-corrected methods operate on individual trajectories in a qualitatively different way, but the results averaged over trajectories are similar for these systems.

Case studies of the time-dependent potential energy surface for dynamics in cavities

The exact time-dependent potential energy surface driving the nuclear dynamics was recently shown to be a useful tool to understand and interpret the coupling of nuclei, electrons, and photons in cavity settings. Here, we provide a detailed analysis of its structure for exactly solvable systems that model two phenomena: cavity-induced suppression of proton-coupled electron-transfer and its dependence on the initial state, and cavity-induced electronic excitation. We demonstrate the inadequacy of simply using a weighted average of polaritonic surfaces to determine the dynamics. Such a weighted average misses a crucial term that redistributes energy between the nuclear and the polaritonic systems, and this term can in fact become a predominant term in determining the nuclear dynamics when several polaritonic surfaces are involved. Evolving an ensemble of classical trajectories on the exact potential energy surface reproduces the nuclear wavepacket quite accurately, while evolving on the weighted polaritonic surface fails after a short period of time. The implications and prospects for application of mixed quantum-classical methods based on this surface are discussed.

Effect of many modes on self-polarization and photochemical suppression in cavities

The standard description of cavity-modified molecular reactions typically involves a single (resonant) mode, while in reality, the quantum cavity supports a range of photon modes. Here, we demonstrate that as more photon modes are accounted for, physicochemical phenomena can dramatically change, as illustrated by the cavity-induced suppression of the important and ubiquitous process of proton-coupled electron-transfer. Using a multi-trajectory Ehrenfest treatment for the photon-modes, we find that self-polarization effects become essential, and we introduce the concept of self-polarization-modified Born-Oppenheimer surfaces as a new construct to analyze dynamics. As the number of cavity photon modes increases, the increasing deviation of these surfaces from the cavity-free Born-Oppenheimer surfaces, together with the interplay between photon emission and absorption inside the widening bands of these surfaces, leads to enhanced suppression. The present findings are general and will have implications for the description and control of cavity-driven physical processes of molecules, nanostructures, and solids embedded in cavities.

Embedding via the Exact Factorization Approach

We present a quantum electronic embedding method derived from the exact factorization approach to calculate static properties of a many-electron system. The method is exact in principle but the practical power lies in utilizing input from a low-level calculation on the entire system in a high-level method computed on a small fragment, as in other embedding methods. Here, the exact factorization approach defines an embedding Hamiltonian on the fragment. Various Hubbard models demonstrate that remarkably accurate ground-state energies are obtained over the full range of weak to strongly correlated systems.

Developing new and understanding old approximations in TDDFT

When a system has evolved far from a ground-state, the adiabatic approximations commonly used in time-dependent density functional theory calculations completely fail in some applications, while giving qualitatively good predictions in others, and sometimes even quantitative predictions. It is not clearly understood why this is so, and developing practical approximations going beyond the adiabatic approximation remains a challenge. This paper explores three different lines of investigation. First, an expression for the exact time-dependent exchange-correlation potential suggests that the accuracy of an adiabatic approximation is intimately related to the deviation between the natural orbital occupation numbers of the physical system and those of the Kohn-Sham system, and we explore this on some exactly-solvable model systems. The exact expression further suggests a path to go beyond the adiabatic approximations, and in the second part we discuss a newly proposed class of memory-dependent approximations developed in this way. Finally, we derive a new expression for the exact exchange-correlation potential from a coupling-constant path integration.

Exact Potential Energy Surface for Molecules in Cavities

We find and analyze the exact time-dependent potential energy surface driving the proton motion for a model of cavity-induced suppression of proton-coupled electron transfer. We show how, in contrast to the polaritonic surfaces, its features directly correlate to the proton dynamics and we discuss cavity modifications of its structure responsible for the suppression. The results highlight the interplay between nonadiabatic effects from coupling to photons and coupling to electrons and suggest caution is needed when applying traditional dynamics methods based on polaritonic surfaces.

On the numerical solution of the exact factorization equations

The exact factorization (EF) approach to coupled electron-ion dynamics recasts the time-dependent molecular Schrödinger equation as two coupled equations, one for the nuclear wavefunction and one for the conditional electronic wavefunction. The potentials appearing in these equations have provided insight into non-adiabatic processes, and new practical non-adiabatic dynamics methods have been formulated starting from these equations. Here, we provide a first demonstration of a self-consistent solution of the exact equations, with a preliminary analysis of their stability and convergence properties. The equations have an unprecedented mathematical form, involving a Hamiltonian outside the class of Hermitian Hamiltonians usually encountered in time-propagation, and so the usual numerical methods for time-dependent Schrödinger fail when applied in a straightforward way to the EF equations. We find an approach that enables stable propagation long enough to witness non-adiabatic behavior in a model system before non-trivial instabilities take over. Implications for the development and analysis of EF-based methods are discussed.

Exploring non-adiabatic approximations to the exchange–correlation functional of TDDFT

Decomposition of the exact time-dependent exchange–correlation potential offers a new starting point to build approximations with memory.

Exact Time-Dependent Exchange-Correlation Potential in Electron Scattering Processes

We identify peak and valley structures in the exact exchange-correlation potential of time-dependent density functional theory that are crucial for time-resolved electron scattering in a model one-dimensional system. These structures are completely missed by adiabatic approximations that, consequently, significantly underestimate the scattering probability. A recently proposed nonadiabatic approximation is shown to correctly capture the approach of the electron to the target when the initial Kohn-Sham state is chosen judiciously, and it is more accurate than standard adiabatic functionals but ultimately fails to accurately capture reflection. These results may explain the underestimation of scattering probabilities in some recent studies on molecules and surfaces.

Charge transfer in time-dependent density functional theory

Charge transfer plays a crucial role in many processes of interest in physics, chemistry, and bio-chemistry. In many applications the size of the systems involved calls for time-dependent density functional theory (TDDFT) to be used in their computational modeling, due to its unprecedented balance between accuracy and efficiency. However, although exact in principle, in practise approximations must be made for the exchange-correlation functional in this theory, and the standard functional approximations perform poorly for excitations which have a long-range charge-transfer component. Intense progress has been made in developing more sophisticated functionals for this problem, which we review. We point out an essential difference between the properties of the exchange-correlation kernel needed for an accurate description of charge-transfer between open-shell fragments and between closed-shell fragments. We then turn to charge-transfer dynamics, which, in contrast to the excitation problem, is a highly non-equilibrium, non-perturbative, process involving a transfer of one full electron in space. This turns out to be a much more challenging problem for TDDFT functionals. We describe dynamical step and peak features in the exact functional evolving over time, that are missing in the functionals currently used. The latter underestimate the amount of charge transferred and manifest a spurious shift in the charge transfer resonance position. We discuss some explicit examples.

Electronic non-adiabatic dynamics in enhanced ionization of isotopologues of hydrogen molecular ions from the exact factorization perspective

An exact-factorization perspective of enhanced ionization in isotopologues of H₂⁺ demonstrates the concept of the exact potential driving the electrons in non-adiabatic motion of molecules in strong fields, and sets a new platform for introducing various approximations.

Time-dependent density functional theory beyond Kohn–Sham Slater determinants

Different choices of initial Kohn Sham wavefunction shape the time-dependent exchange–correlation potential.

Exact Potential Driving the Electron Dynamics in Enhanced Ionization of H(2)(+)

It was recently shown that the exact factorization of the electron-nuclear wave function allows the construction of a Schrödinger equation for the electronic system, in which the potential contains exactly the effect of coupling to the nuclear degrees of freedom and any external fields. Here we study the exact potential acting on the electron in charge-resonance enhanced ionization in a model one-dimensional H(2)(+) molecule. We show there can be significant differences between the exact potential and that used in the traditional quasistatic analyses, arising from nonadiabatic coupling to the nuclear system, and that these are crucial to include for accurate simulations of time-resolved ionization dynamics and predictions of the ionization yield.

The exact forces on classical nuclei in non-adiabatic charge transfer

The decomposition of electronic and nuclear motion presented in Abedi et al. [Phys. Rev. Lett. 105, 123002 (2010)] yields a time-dependent potential that drives the nuclear motion and fully accounts for the coupling to the electronic subsystem. Here, we show that propagation of an ensemble of independent classical nuclear trajectories on this exact potential yields dynamics that are essentially indistinguishable from the exact quantum dynamics for a model non-adiabatic charge transfer problem. We point out the importance of step and bump features in the exact potential that are critical in obtaining the correct splitting of the quasiclassical nuclear wave packet in space after it passes through an avoided crossing between two Born-Oppenheimer surfaces and analyze their structure. Finally, an analysis of the exact potentials in the context of trajectory surface hopping is presented, including preliminary investigations of velocity-adjustment and the force-induced decoherence effect.

Time-resolved spectroscopy in time-dependent density functional theory: an exact condition

A fundamental property of a quantum system driven by an external field is that when the field is turned off the positions of its response frequencies are independent of the time at which the field is turned off. We show that this leads to an exact condition for the exchange-correlation potential of time-dependent density functional theory. The Kohn-Sham potential typically continues to evolve after the field is turned off, which leads to time dependence in the response frequencies of the Kohn-Sham response function. The exchange-correlation kernel must cancel out this time dependence. The condition is typically violated by approximations currently in use, as we demonstrate by several examples, which has severe consequences for their predictions of time-resolved spectroscopy.

Kinetic and interaction components of the exact time-dependent correlation potential

The exact exchange-correlation (xc) potential of time-dependent density functional theory has been shown to have striking features. For example, step and peak features are generically found when the system is far from its ground-state, and these depend nonlocally on the density in space and time. We analyze the xc potential by decomposing it into kinetic and interaction components and comparing each with their exact-adiabatic counterparts, for a range of dynamical situations in model one-dimensional two-electron systems. We find that often, but not always, the kinetic contribution is largely responsible for these features that are missed by the adiabatic approximation. The adiabatic approximation often makes a smaller error for the interaction component, which we write in two parts, one being the Coulomb potential due to the time-dependent xc hole. Non-adiabatic features of the kinetic component were also larger than those of the interaction component in cases that we studied when there is negligible step structure. In ground-state situations, step and peak structures arise in cases of static correlation, when more than one determinant is essential to describe the interacting state. We investigate the time-dependent natural orbital occupation numbers and find the corresponding relation between these and the dynamical step is more complex than for the ground-state case.

Laser-induced electron localization in H2+: mixed quantum-classical dynamics based on the exact time-dependent potential energy surface

The exact nuclear time-dependent potential energy surface for laser-induced electron localization is studied with a view to developing a mixed quantum-classical dynamics method for strong-field processes.

Challenging adiabatic time-dependent density functional theory with a Hubbard dimer: the case of time-resolved long-range charge transfer

Adiabatic TDDFT fails for time-resolved charge-transfer in a simple model, even when linear response spectra is accurate.

Dynamics of Charge-Transfer Processes with Time-Dependent Density Functional Theory

We show that whenever an electron transfers between closed-shell molecular fragments, the exact correlation potential of time-dependent density functional theory develops a step and peak structure in the bonding region. This structure has a density dependence that is nonlocal both in space and in time that even the exact adiabatic ground-state exchange-correlation functional fails to capture it. For charge-transfer between open-shell fragments, an initial step and peak vanish as the charge-transfer state is reached. The inability of usual approximations to develop these structures leads to inaccurate charge-transfer dynamics. This is illustrated by the complete lack of Rabi oscillations in the dipole moment under conditions of resonant charge transfer for an exactly solvable model system. The results transcend the model and are applicable to more realistic molecular complexes.

Correlated electron-nuclear dynamics: exact factorization of the molecular wavefunction

It was recently shown [A. Abedi, N. T. Maitra, and E. K. U. Gross, Phys. Rev. Lett. 105, 123002 (2010)] that the complete wavefunction for a system of electrons and nuclei evolving in a time-dependent external potential can be exactly factorized into an electronic wavefunction and a nuclear wavefunction. The concepts of an exact time-dependent potential energy surface (TDPES) and exact time-dependent vector potential emerge naturally from the formalism. Here, we present a detailed description of the formalism, including a full derivation of the equations that the electronic and nuclear wavefunctions satisfy. We demonstrate the relationship of this exact factorization to the traditional Born-Oppenheimer expansion. A one-dimensional model of the H(2)(+) molecule in a laser field shows the usefulness of the exact TDPES in interpreting coupled electron-nuclear dynamics: we show how features of its structure indicate the mechanism of dissociation. We compare the exact TDPES with potential energy surfaces from the time-dependent Hartree-approach, and also compare traditional Ehrenfest dynamics with Ehrenfest dynamics on the exact TDPES.

Universal dynamical steps in the exact time-dependent exchange-correlation potential

We show that the exact exchange-correlation potential of time-dependent density-functional theory displays dynamical step structures that have a spatially nonlocal and time nonlocal dependence on the density. Using one-dimensional two-electron model systems, we illustrate these steps for a range of nonequilibrium dynamical situations relevant for modeling of photochemical or physical processes: field-free evolution of a nonstationary state, resonant local excitation, resonant complete charge transfer, and evolution under an arbitrary field. A lack of these steps in the usual approximations yields inaccurate dynamics, for example, predicting faster dynamics and incomplete charge transfer.

Comparison of some dispersion-corrected and traditional functionals with CCSD(T) and MP2 ab initio methods: dispersion, induction, and basis set superposition error

We compare dispersion and induction interactions for noble gas dimers and for Ne, methane, and 2-butyne with HF and LiF using a variety of functionals (including some specifically parameterized to evaluate dispersion interactions) with ab initio methods including CCSD(T) and MP2. We see that inductive interactions tend to enhance dispersion and may be accompanied by charge-transfer. We show that the functionals do not generally follow the expected trends in interaction energies, basis set superposition errors (BSSE), and interaction distances as a function of basis set size. The functionals parameterized to treat dispersion interactions often overestimate these interactions, sometimes by quite a lot, when compared to higher level calculations. Which functionals work best depends upon the examples chosen. The B3LYP and X3LYP functionals, which do not describe pure dispersion interactions, appear to describe dispersion mixed with induction about as accurately as those parametrized to treat dispersion. We observed significant differences in high-level wavefunction calculations in a basis set larger than those used to generate the structures in many of the databases. We discuss the implications for highly parameterized functionals based on these databases, as well as the use of simple potential energy for fitting the parameters rather than experimentally determinable thermodynamic state functions that involve consideration of vibrational states.

Effect of cusps in time-dependent quantum mechanics

Within nonrelativistic quantum mechanics, spatial cusps in initial wave functions can lead to nonanalytic behavior in time. We suggest a method for calculating the short-time behavior in such situations. For these cases, the density does not match its Taylor expansion in time, but the Runge-Gross proof of time-dependent density functional theory still holds, as it requires only the potential to be time analytic.

Electron correlation via frozen Gaussian dynamics

We investigate the accuracy and efficiency of the semiclassical frozen Gaussian method in describing electron dynamics in real time. Model systems of two soft-Coulomb-interacting electrons are used to study correlated dynamics under non-perturbative electric fields, as well as the excitation spectrum. The results show that a recently proposed method that combines exact-exchange with semiclassical correlation to propagate the one-body density-matrix holds promise for electron dynamics in many situations that either wavefunction or density-functional methods have difficulty describing. The results also however point out challenges in such a method that need to be addressed before it can become widely applicable.

Exact factorization of the time-dependent electron-nuclear wave function

We present an exact decomposition of the complete wave function for a system of nuclei and electrons evolving in a time-dependent external potential. We derive formally exact equations for the nuclear and electronic wave functions that lead to rigorous definitions of a time-dependent potential energy surface (TDPES) and a time-dependent geometric phase. For the H(2)(+) molecular ion exposed to a laser field, the TDPES proves to be a useful interpretive tool to identify different mechanisms of dissociation.

Semiclassical electron correlation in density-matrix time propagation

Lack of memory (locality in time) is a major limitation of almost all present time-dependent density-functional approximations. By using semiclassical dynamics to compute correlation effects in the time propagation of the density matrix, we incorporate memory, including initial-state dependence, as well as changing occupation numbers, and predict more observables in strong-field applications.