We show that propagation with the best possible adiabatic approximation in time-dependent density-functional theory fails to properly transfer charge in an asymmetric two-site Hubbard model when beginning in the ground state. The approximation is adiabatic but exact otherwise, constructed from the exact ground-state exchange-correlation functional that we compute via constrained search. The model shares the essential features of charge-transfer dynamics in a real-space long-range molecule, so the results imply that the best possible adiabatic approximation, despite being able to capture nonlocal ground-state step features relevant to dissociation and charge-transfer excitations, cannot capture fully time-resolved charge-transfer dynamics out of the ground state.
We show that propagation with the best possible adiabatic approximation in time-dependent density-functional theory fails to properly transfer charge in an asymmetric two-site Hubbard model when beginning in the ground state. The approximation is adiabatic but exact otherwise, constructed from the exact ground-state exchange-correlation functional that we compute via constrained search. The model shares the essential features of charge-transfer dynamics in a real-space long-range molecule, so the results imply that the best possible adiabatic approximation, despite being able to capture nonlocal ground-state step features relevant to dissociation and charge-transfer excitations, cannot capture fully time-resolved charge-transfer dynamics out of the ground state.
We address and quantify the role of nonadiabaticity (``memory effects'') in the exchange-correlation (xc) functional of time-dependent density functional theory (TDDFT) for describing nonlinear dynamics of many-body systems. Time-dependent resonant processes are particularly challenging for available TDDFT approximations, due to their strong nonlinear and nonadiabatic character. None of the known approximate density functionals are able to cope with this class of problems in a satisfactory manner. In this work we look at the prototypical example of the resonant processes by considering Rabi oscillations within the exactly soluble two-site Hubbard model. We construct the exact adiabatic xc functional and show that (i) it does not reproduce correctly resonant Rabi dynamics, and (ii) there is a sizable nonadiabatic contribution to the exact xc potential, which turns out to be small only at the beginning and at the end of the Rabi cycle when the ground-state population is dominant. We then propose a ``two-level'' approximation for the time-dependent xc potential which can capture Rabi dynamics in the two-site problem. It works well both for resonant and for detuned Rabi oscillations and becomes essentially exact in the linear response regime.
This review explains the relationship between density functional theory and strongly correlated models using the simplest possible example, the two-site Hubbard model. The relationship to traditional quantum chemistry is included. Even in this elementary example, where the exact ground-state energy and site occupations can be found analytically, there is much to be explained in terms of the underlying logic and aims of density functional theory. Although the usual solution is analytic, the density functional is given only implicitly. We overcome this difficulty using the Levy-Lieb construction to create a parametrization of the exact function with negligible errors. The symmetric case is most commonly studied, but we find a rich variation in behavior by including asymmetry, as strong correlation physics vies with charge-transfer effects. We explore the behavior of the gap and the many-body Green's function, demonstrating the 'failure' of the Kohn-Sham (KS) method to reproduce the fundamental gap. We perform benchmark calculations of the occupation and components of the KS potentials, the correlation kinetic energies, and the adiabatic connection. We test several approximate functionals (restricted and unrestricted Hartree-Fock and Bethe ansatz local density approximation) to show their successes and limitations. We also discuss and illustrate the concept of the derivative discontinuity. Useful appendices include analytic expressions for density functional energy components, several limits of the exact functional (weak- and strong-coupling, symmetric and asymmetric), various adiabatic connection results, proofs of exact conditions for this model, and the origin of the Hubbard model from a minimal basis model for stretched H2.
Most applications of time-dependent density-functional theory (TDDFT) use the adiabatic local-density approximation (ALDA) for the dynamical exchange-correlation potential V(xc)(r,t). An exact (i.e., nonadiabatic) extension of the ground-state LDA into the dynamical regime leads to a V(xc)(r,t) with a memory, which causes the electron dynamics to become dissipative. To illustrate and explain this nonadiabatic behavior, this paper studies the dynamics of two interacting electrons on a two-dimensional quantum strip of finite size, comparing TDDFT within and beyond the ALDA with numerical solutions of the two-electron time-dependent Schrodinger equation. It is shown explicitly how dissipation arises through multiple particle-hole excitations, and how the nonadiabatic extension of the ALDA fails for finite systems but becomes correct in the thermodynamic limit.
We obtain the exact Kohn-Sham potentials VKS of time-dependent density-functional theory for one-dimensional Hubbard chains, driven by a dc external field, using the time-dependent electron density and current density obtained from exact many-body time evolution. The exact Vxc is compared to the adiabatically exact Vxcad and the ``instantaneous ground state'' Vxcigs. The latter is shown to work effectively in some cases when the former fails. Approximations for the exchange-correlation potential Vxc and its gradient, based on the local density and on the local current density, are also considered and both physical quantities are observed to be far outside the reach of any possible local approximation. Insight into the respective roles of ground-state and excited-state correlation in the time-dependent system, as reflected in the potentials, is provided by the pair correlation function.
This review explains the relationship between density functional theory and strongly correlated models using the simplest possible example, the two-site Hubbard model. The relationship to traditional quantum chemistry is included. Even in this elementary example, where the exact ground-state energy and site occupations can be found analytically, there is much to be explained in terms of the underlying logic and aims of Density Functional Theory. Although the usual solution is analytic, the density functional is given only implicitly. We overcome this difficulty using the Levy-Lieb construction to create a parametrization of the exact function with negligible errors. The symmetric case is most commonly studied, but we find a rich variation in behavior by including asymmetry, as strong correlation physics vies with charge-transfer effects. We explore the behavior of the gap and the many-body Green’s function, demonstrating the ‘failure ’ of the Kohn-Sham method to reproduce the fundamental gap. We perform benchmark calculations of the occupation and components of the KS potentials, the correlation kinetic energies, and the adiabatic connection. We test several approximate functionals (restricted and unrestricted Hartree-Fock and Bethe Ansatz Local Density Approximation) to show their successes and limitations. We also discuss and illustrate the concept of the derivative discontinuity. Useful appendices include analytic expressions for Density
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.