Exact density-functional potentials for time-dependent quasiparticles

Inverse Kohn-Sham Density Functional Theory: Progress and Challenges

Inverse Kohn-Sham (iKS) methods are needed to fully understand the one-to-one mapping between densities and potentials on which density functional theory is based. They can contribute to the construction of empirical exchange-correlation functionals and to the development of techniques for density-based embedding. Unlike the forward Kohn-Sham problems, numerical iKS problems are ill-posed and can be unstable. We discuss some of the fundamental and practical difficulties of iKS problems with constrained-optimization methods on finite basis sets. Various factors that affect the performance are systematically compared and discussed, both analytically and numerically, with a focus on two of the most practical methods: the Wu-Yang method (WY) and the partial differential equation constrained optimization (PDE-CO). Our analysis of the WY and PDE-CO highlights the limitation of finite basis sets. We introduce new ideas to make iKS problems more tractable, provide an overall strategy for performing numerical density-to-potential inversions, and discuss challenges and future directions.

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.

Asymptotic behavior of the Hartree-exchange and correlation potentials in ensemble density functional theory

We show that the Hartree-exchange and correlation potentials of ensemble systems can have unexpected features, including non-vanishing asymptotic constants and non-trivial screening of electrons. These features are demonstrated here on Li, C, and F.

Time-Dependent Multicomponent Density Functional Theory for Coupled Electron-Positron Dynamics

Electron-positron interactions have been utilized in various fields of science. Here we develop time-dependent multicomponent density functional theory to study the coupled electron-positron dynamics from first principles. We prove that there are coupled time-dependent single-particle equations that can provide the electron and positron density dynamics, and derive the formally exact expression for their effective potentials. Introducing the adiabatic local density approximation to time-dependent electron-positron correlation, we apply the theory to the dynamics of a positronic lithium hydride molecule under a laser field. We demonstrate the significance of the coupling between electronic and positronic motion by revealing the complex positron detachment mechanism and the suppression of electronic resonant excitation by the screening effect of the positron.

Nonadiabatic Dynamics in Single-Electron Tunneling Devices with Time-Dependent Density-Functional Theory

We simulate the dynamics of a single-electron source, modeled as a quantum dot with on-site Coulomb interaction and tunnel coupling to an adjacent lead in time-dependent density-functional theory. Based on this system, we develop a time-nonlocal exchange-correlation potential by exploiting analogies with quantum-transport theory. The time nonlocality manifests itself in a dynamical potential step. We explicitly link the time evolution of the dynamical step to physical relaxation timescales of the electron dynamics. Finally, we discuss prospects for simulations of larger mesoscopic systems.

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.

Real-space grids and the Octopus code as tools for the development of new simulation approaches for electronic systems

Real-space grids are a powerful alternative for the simulation of electronic systems. One of the main advantages of the approach is the flexibility and simplicity of working directly in real space where the different fields are discretized on a grid, combined with competitive numerical performance and great potential for parallelization. These properties constitute a great advantage at the time of implementing and testing new physical models. Based on our experience with the Octopus code, in this article we discuss how the real-space approach has allowed for the recent development of new ideas for the simulation of electronic systems. Among these applications are approaches to calculate response properties, modeling of photoemission, optimal control of quantum systems, simulation of plasmonic systems, and the exact solution of the Schrödinger equation for low-dimensionality systems.

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

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

Time-dependent exchange and tunneling: detection at the same place of two electrons emitted simultaneously from different sources

Two-particle scattering probabilities in tunneling scenarios with exchange interaction are analyzed with quasi-particle wave packets. Two initial one-particle wave packets (with opposite central momentums) are spatially localized at each side of a barrier. After impinging upon a tunneling barrier, each wave packet splits into transmitted and reflected components. When the initial two-particle anti-symmetrical state is defined as a Slater determinant of any type of (normalizable) one-particle wave packet, it is shown that the probability of detecting two (identically injected) electrons at the same side of the barrier is different from zero in very common (single or double barrier) scenarios. In some particular scenarios, the transmitted and reflected components become orthogonal and the mentioned probabilities reproduce those values associated to distinguishable particles. These unexpected non-zero probabilities are still present when non-separable Coulomb interaction or non-symmetrical potentials are considered. On the other hand, for initial wave packets close to Hamiltonian eigenstates, the usual zero two-particle probability for electrons at the same side of the barrier found in the literature is recovered. The generalization to many-particle scattering probabilities with quasi-particle wave packets for low and high phase-space density are also analyzed. The far-reaching consequences of these non-zero probabilities in the accurate evaluation of quantum noise in mesoscopic systems are briefly indicated.

Existence, uniqueness, and construction of the density-potential mapping in time-dependent density-functional theory

In this work we review the mapping from densities to potentials in quantum mechanics, which is the basic building block of time-dependent density-functional theory and the Kohn-Sham construction. We first present detailed conditions such that a mapping from potentials to densities is defined by solving the time-dependent Schrödinger equation. We specifically discuss intricacies connected with the unboundedness of the Hamiltonian and derive the local-force equation. This equation is then used to set up an iterative sequence that determines a potential that generates a specified density via time propagation of an initial state. This fixed-point procedure needs the invertibility of a certain Sturm-Liouville problem, which we discuss for different situations. Based on these considerations we then present a discussion of the famous Runge-Gross theorem which provides a density-potential mapping for time-analytic potentials. Further we give conditions such that the general fixed-point approach is well-defined and converges under certain assumptions. Then the application of such a fixed-point procedure to lattice Hamiltonians is discussed and the numerical realization of the density-potential mapping is shown. We conclude by presenting an extension of the density-potential mapping to include vector-potentials and photons.

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.

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.

Strong correlation in Kohn-Sham density functional theory

We use the exact strong-interaction limit of the Hohenberg-Kohn energy density functional to approximate the exchange-correlation energy of the restricted Kohn-Sham scheme. Our approximation corresponds to a highly nonlocal density functional whose functional derivative can be easily constructed, thus transforming exactly, in a physically transparent way, an important part of the electron-electron interaction into an effective local one-body potential. We test our approach on quasi-one-dimensional systems, showing that it captures essential features of strong correlation that restricted Kohn-Sham calculations using the currently available approximations cannot describe.