Exact time-dependent density-functional theory (DFT) separates the exchange-correlation potential into interaction and kinetic-correlation components, but the structural relation between them remains unclear. We establish a representability constraint based on the off-diagonal expansion of the one-electron reduced density matrix (1RDM) relative to a chosen reference dynamics. A non-zero linear off-diagonal term of the 1RDM generates a current density beyond the representable space of the reference, while a nontrivial quadratic off-diagonal structure activates the kinetic-correlation component. The resulting hierarchy provides a criterion for assessing, from the spatial locality of electronic correlation, whether a given exchange-correlation functional or a more general electron-correlation approximation remains adequate. In a one-dimensional two-electron correlation quench, retaining the exact interaction component while setting the kinetic component to zero suppresses the density broadening of the exact evolution. These results establish a spatial-order hierarchy among the force components governing the density dynamics and provide an exact constraint for the construction and assessment of nonadiabatic approximations.
Exact time-dependent density-functional theory (DFT) separates the exchange-correlation potential into interaction and kinetic-correlation components, but the structural relation between them remains unclear. We establish a representability constraint based on the off-diagonal expansion of the one-electron reduced density matrix (1RDM) relative to a chosen reference dynamics. A non-zero linear off-diagonal term of the 1RDM generates a current density beyond the representable space of the reference, while a nontrivial quadratic off-diagonal structure activates the kinetic-correlation component. The resulting hierarchy provides a criterion for assessing, from the spatial locality of electronic correlation, whether a given exchange-correlation functional or a more general electron-correlation approximation remains adequate. In a one-dimensional two-electron correlation quench, retaining the exact interaction component while setting the kinetic component to zero suppresses the density broadening of the exact evolution. These results establish a spatial-order hierarchy among the force components governing the density dynamics and provide an exact constraint for the construction and assessment of nonadiabatic approximations.
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.
A decomposition of the exact exchange-correlation potential of time-dependent density functional theory into an interaction component and a kinetic component offers a new starting point for non-adiabatic approximations. The components are expressed in terms of the exchange-correlation hole and the difference between the one-body density matrix of the interacting and Kohn-Sham systems, which must be approximated in terms of quantities accessible from the Kohn-Sham evolution. We explore several preliminary approximations, evaluate their fulfillment of known exact conditions, and test their performance on simple model systems for which available exact solutions indicate the significance of going beyond the adiabatic approximation.
Time-dependent quantal density-functional theory (Q-DFT) is a description of the s-system of noninteracting fermions with electronic density equivalent to that of Schr\"odinger theory, in terms of fields whose sources are quantum-mechanical expectations of Hermitian operators. The theory delineates and defines the contribution of each type of electron correlation to the local electron-interaction potential ${\ensuremath{\nu}}_{\mathrm{ee}}(\mathbf{r},t)$ of the s system. These correlations are due to the Pauli exclusion principle, Coulomb repulsion, correlation-kinetic, and correlation-current-density effects, the latter two resulting, respectively, from the difference in kinetic energy and current density between the interacting Schr\"odinger and noninteracting systems. We employ Q-DFT to prove the following sum rules and properties of the s system: (i) the components of the potential due to these correlations separately exert no net force on the system; (ii) the torque of the potential is finite and due solely to correlation-current-density effects; (iii) two sum rules involving the curl of the dynamic electron-interaction kernel defined as the functional derivative of ${\ensuremath{\nu}}_{\mathrm{ee}}(\mathbf{r},t)$ are derived and shown to depend on the frequency dependent correlation-current-density effect. Furthermore, via adiabatic coupling constant (\ensuremath{\lambda}) perturbation theory, we prove: (iv) the exchange potential ${\ensuremath{\nu}}_{x}(\mathbf{r},t)$ is the work done in a conservative field representative of Pauli correlations and lowest-order $O(\ensuremath{\lambda})$ correlation-kinetic and correlation-current-density effects; (v) the correlation potential ${\ensuremath{\nu}}_{c}(\mathbf{r},t)$ commences in $O({\ensuremath{\lambda}}^{2}),$ and, at each order, it is the work done in a conservative field representative of Coulomb correlations and correlation-kinetic and correlation-current-density effects; (vi) we derive the integral virial theorem relating ${\ensuremath{\nu}}_{\mathrm{ee}}(\mathbf{r},t)$ to the electron-interaction and correlation-kinetic energy for arbitrary coupling constant strength \ensuremath{\lambda}, and show there are no explicit correlation-current-density contributions to the energy. From this integral virial theorem we (vii) obtain the fully interacting $(\ensuremath{\lambda}=1)$ and exchange-only $(\ensuremath{\lambda}=0)$ integral virial theorems as special cases, the latter showing there is no explicit correlation-kinetic contribution to the exchange energy; and (viii) write expressions for the electron-interaction and correlation-kinetic actions for arbitrary coupling constant \ensuremath{\lambda} in terms of the corresponding fields.
Given the time-evolution of an electron charge density, the local potential in Kohn-Sham time-dependent density functional theory (KS-TDDFT) can be modeled as a sum of instantaneous and dynamic contributions by assuming a certain form of the time-dependent Kohn-Sham (TD-KS) orbitals. The instantaneous part is obtained numerically using methods from ground-state density functional theory (DFT) and the dynamic part is expressed in terms of a velocity potential that depends on the electron current density. The suggested form of the TD-KS orbitals satisfies several known constraints (orthonormality, N-representability, J-representability), and the domain of validity is shown to depend on the evolution of the instantaneous quantities. Through this decomposition, we can relate time-dependent and ground -state V-representability. The resulting potentials are shown numerically to approximate the exact time-dependent Kohn-Sham potentials for a set of 3 non-singlet two-particle systems (a Kohn-mode, a Coulomb explosion, and a double quantum well) where the exact solutions and reference densities are known or obtained through configuration interaction (CI) approaches.
We consider an analytically solvable model of two interacting electrons that allows for the calculation of the exact exchange-correlation kernel of time-dependent density functional theory. This kernel, as well as the corresponding density response function, is studied in the limit of large repulsive interactions between the electrons and we give analytical results for these quantities as an asymptotic expansion in powers of the square root of the interaction strength. We find that in the strong interaction limit the three leading terms in the expansion of the kernel act instantaneously while memory terms only appear in the next orders. We further derive an alternative expansion for the kernel in the strong interaction limit on the basis of the theory developed previously [Phys. Chem. Chem. Phys. 18, 21092 (2016)] using the formalism of strictly correlated electrons in the adiabatic approximation. We find that the first two leading terms in this series, corresponding to the strictly correlated limit and its zero-point vibration correction, coincide with the two leading terms of the exact expansion. We finally analyze the spatial nonlocality of these terms and show when the adiabatic approximation breaks down. The ability to reproduce the exact kernel in the strong interaction limit indicates that the adiabatic strictly correlated electron formalism is useful for studying the density response and excitation properties of other systems with strong electronic interactions.
Most present applications of time-dependent density functional theory use adiabatic functionals, i.e., the effective potential at time t is determined solely by the density at the same time. This paper discusses a method that aims to go beyond this approximation, by incorporating "memory" effects: the derived exchange-correlation potential will depend not only on present densities but also on the past. In order to ensure the potentials are causal, we formulate the action on the Keldysh contour for electrons in electromagnetic fields, from which we derive suitable Kohn-Sham equations. The exchange-correlation action is now a functional of the electron density and velocity field. A specific action functional is constructed which is Galilean invariant and yields a causal exchange-correlation vector potential for the Kohn-Sham equations incorporating memory effects. We show explicitly that the net exchange-correlation Lorentz force is zero. The potential is consistent with known dynamical properties of the homogeneous electron gas (in the linear response limit).