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.