Real-time electron dynamics with exact-exchange time-dependent density-functional theory

Size-dependent errors in real-time electron density propagation

Real-time (RT) electron density propagation with time-dependent density functional theory (TDDFT) or Hartree-Fock (TDHF) is one of the most popular methods to model the charge transfer in molecules and materials. However, both RT-TDHF and RT-TDDFT within the adiabatic approximation are known to produce inaccurate evolution of the electron density away from the ground state in model systems, leading to large errors in charge transfer and erroneous shifting of peaks in absorption spectra. Given the poor performance of these methods with small model systems and the widespread use of the methods with larger molecular and material systems, here we bridge the gap in our understanding of these methods and examine the size-dependence of errors in RT density propagation. We analyze the performance of RT density propagation for systems of increasing size during the application of a continuous resonant field to induce Rabi-like oscillations, during charge-transfer dynamics, and for peak shifting in simulated absorption spectra. We find that the errors in the electron dynamics are indeed size dependent for these phenomena, with the largest system producing the results most aligned with those expected from linear response theory. The results suggest that although the RT-TDHF and RT-TDDFT methods may produce severe errors for model systems, the errors in charge transfer and resonantly driven electron dynamics may be much less significant for more realistic, large-scale molecules and materials.

Making ab initio QED functional(s): Nonperturbative and photon-free effective frameworks for strong light-matter coupling

Strong light-matter coupling provides a promising path for the control of quantum matter where the latter is routinely described from first principles. However, combining the quantized nature of light with this ab initio tool set is challenging and merely developing as the coupled light-matter Hilbert space is conceptually different and computational cost quickly becomes overwhelming. In this work, we provide a nonperturbative photon-free formulation of quantum electrodynamics (QED) in the long-wavelength limit, which is formulated solely on the matter Hilbert space and can serve as an accurate starting point for such ab initio methods. The present formulation is an extension of quantum mechanics that recovers the exact results of QED for the zero- and infinite-coupling limit and the infinite-frequency as well as the homogeneous limit, and we can constructively increase its accuracy. We show how this formulation can be used to devise approximations for quantum-electrodynamical density-functional theory (QEDFT), which in turn also allows us to extend the ansatz to the full minimal-coupling problem and to nonadiabatic situations. Finally, we provide a simple local density-type functional that takes the strong coupling to the transverse photon degrees of freedom into account and includes the correct frequency and polarization dependence. This QEDFT functional accounts for the quantized nature of light while remaining computationally simple enough to allow its application to a large range of systems. All approximations allow the seamless application to periodic systems.

A reformulation of time-dependent Kohn-Sham theory in terms of the second time derivative of the density

The Kohn-Sham approach to time-dependent density-functional theory (TDDFT) can be formulated, in principle, exactly by invoking the force-balance equation for the density, which leads to an explicit expression for the exchange-correlation potential as an implicit density functional. It is shown that this suggests a reformulation of TDDFT in terms of the second time derivative of the density, rather than the density itself. The result is a time-local Kohn-Sham scheme of second-order in time whose causal structure is more transparent than that of the usual Kohn-Sham formalism. The scheme can be used to construct new approximations at the exchange-only level and beyond, and it offers a straightforward definition of the exact adiabatic approximation.

Time-dependent density functional theory applied to average atom opacity

We focus on studying the opacity of iron, chromium, and nickel plasmas at conditions relevant to experiments carried out at Sandia National Laboratories [J. E. Bailey et al., Nature (London) 517, 56 (2015)NATUAS0028-083610.1038/nature14048]. We calculate the photoabsorption cross sections and subsequent opacity for plasmas using linear-response time-dependent density functional theory (TD-DFT). Our results indicate that the physics of channel mixing accounted for in linear-response TD-DFT leads to an increase in the opacity in the bound-free quasicontinuum, where the Sandia experiments indicate that models underpredict iron opacity. However, the increase seen in our calculations is only in the range of 5%-10%. Further, we do not see any change in this trend for chromium and nickel. This behavior indicates that channel mixing effects do not explain the trends in opacity observed in the Sandia experiments.

Octopus, a computational framework for exploring light-driven phenomena and quantum dynamics in extended and finite systems

Over the last few years, extraordinary advances in experimental and theoretical tools have allowed us to monitor and control matter at short time and atomic scales with a high degree of precision. An appealing and challenging route toward engineering materials with tailored properties is to find ways to design or selectively manipulate materials, especially at the quantum level. To this end, having a state-of-the-art ab initio computer simulation tool that enables a reliable and accurate simulation of light-induced changes in the physical and chemical properties of complex systems is of utmost importance. The first principles real-space-based Octopus project was born with that idea in mind, i.e., to provide a unique framework that allows us to describe non-equilibrium phenomena in molecular complexes, low dimensional materials, and extended systems by accounting for electronic, ionic, and photon quantum mechanical effects within a generalized time-dependent density functional theory. This article aims to present the new features that have been implemented over the last few years, including technical developments related to performance and massive parallelism. We also describe the major theoretical developments to address ultrafast light-driven processes, such as the new theoretical framework of quantum electrodynamics density-functional formalism for the description of novel light-matter hybrid states. Those advances, and others being released soon as part of the Octopus package, will allow the scientific community to simulate and characterize spatial and time-resolved spectroscopies, ultrafast phenomena in molecules and materials, and new emergent states of matter (quantum electrodynamical-materials).

First-Principles Modeling of Electronic Stopping in Complex Matter under Ion Irradiation

Electronic stopping refers to the dynamical energy-transfer process to electrons in matter from highly energetic charged particles such as high-velocity protons. We discuss recent progress in theoretical studies of electronic stopping in condensed matter under ion irradiation, focusing on modern electronic structure theory's role in enabling the study of electronic excitation dynamics that result from the energy transfer. In the last few decades, first-principles simulation approaches based on real-time time-dependent density functional theory have greatly advanced the field. While linear response theory is widely used to study electronic stopping processes, especially for simple solids, novel first-principles dynamics approaches now allow us to study chemically complex systems and also yield detailed descriptions of electronic excitations at the molecular scale. Outstanding challenges for further advancement of electronic stopping modeling are also discussed from the viewpoint of electronic structure theory.

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.

Superadiabatic Forces via the Acceleration Gradient in Quantum Many-Body Dynamics

We apply the formally exact quantum power functional framework (J. Chem. Phys.2015, 143, 174108) to a one-dimensional Hooke’s helium model atom. The physical dynamics are described on the one-body level beyond the density-based adiabatic approximation. We show that gradients of both the microscopic velocity and acceleration field are required to correctly describe the effects due to interparticle interactions. We validate the proposed analytical forms of the superadiabatic force and transport contributions by comparison to one-body data from exact numerical solution of the Schrödinger equation. Superadiabatic contributions beyond the adiabatic approximation are important in the dynamics and they include effective dissipation.

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.

Time-Local Equation for the Exact Optimized Effective Potential in Time-Dependent Density Functional Theory

A long-standing challenge in the time-dependent density functional theory is to efficiently solve the exact time-dependent optimized effective potential (TDOEP) integral equation derived from orbital-dependent functionals, especially for the study of nonadiabatic dynamics in time-dependent external fields. In this Letter, we formulate a completely equivalent time-local TDOEP equation that admits a unique real-time solution in terms of time-dependent Kohn-Sham and effective memory orbitals. The time-local formulation is numerically implemented, with the incorporation of exponential memory loss to address the unaccounted for correlation component in the exact-exchange-only functional, to enable the study of the many-electron dynamics of a one-dimensional hydrogen chain. It is shown that the long time behavior of the electric dipole converges correctly and the zero-force theorem is fulfilled in the current implementation.

Functional theories of thermoelectric phenomena

We review the progress that has been recently made in the application of time-dependent density functional theory to thermoelectric phenomena. As the field is very young, we emphasize open problems and fundamental issues. We begin by introducing the formal structure of thermal density functional theory, a density functional theory with two basic variables-the density and the energy density-and two conjugate fields-the ordinary scalar potential and Luttinger's thermomechanical potential. The static version of this theory is contrasted with the familiar finite-temperature density functional theory, in which only the density is a variable. We then proceed to constructing the full time-dependent non equilibrium theory, including the practically important Kohn-Sham equations that go with it. The theory is shown to recover standard results of the Landauer theory for thermal transport in the steady state, while showing greater flexibility by allowing a description of fast thermal response, temperature oscillations and related phenomena. Several results are presented here for the first time, i.e. the proof of invertibility of the thermal response function in the linear regime, the full expression of the thermal currents in the presence of Luttinger's thermomechanical potential, an explicit prescription for the evaluation of the Kohn-Sham potentials in the adiabatic local density approximation, a detailed discussion of the leading dissipative corrections to the adiabatic local density approximation and the thermal corrections to the resistivity that follow from it.

Optimized Effective Potential for Quantum Electrodynamical Time-Dependent Density Functional Theory

We propose an orbital exchange-correlation functional for applying time-dependent density functional theory to many-electron systems coupled to cavity photons. The time nonlocal equation for the electron-photon optimized effective potential (OEP) is derived. In the static limit our OEP energy functional reduces to the Lamb shift of the ground state energy. We test the new approximation in the Rabi model. It is shown that the OEP (i) reproduces quantitatively the exact ground-state energy from the weak to the deep strong coupling regime and (ii) accurately captures the dynamics entering the ultrastrong coupling regime. The present formalism opens the path to a first-principles description of correlated electron-photon systems, bridging the gap between electronic structure methods and quantum optics for real material applications.

Nonperturbative time-dependent density functional theory (TDDFT) and time-dependent electron localization function (TDELF) study of the ionization of OCS and CS2 with ultrashort intense laser pulses — Orientational effects

The nonlinear nonperturbative response of OCS and CS2 to ultrashort (few cycles) intense laser pulses was studied numerically by time-dependent density functional theory (TDDFT) methods to understand molecular ionization as a function of laser–molecule orientation. A time-dependent electron localization function(TDELF) was used to visualize the nonlinear nonperturbative electron transfer occurring during the laser pulse. It was found that, for intensities I > 3.5 × 1014 W/cm2, the inner shell Kohn–Sham (KS) molecular orbitals contribute significantly to the ionization, whereas for the intensity I < 3.5 × 1014 W/cm2, the highest occupied molecular orbital (HOMO) shows the dominant response to the field. In general, the ionization rate maxima correspond to the alignment of maximum KS orbital densities with the laser pulse polarization instead of orbital ionization potentials (IP). These findings are corroborated through analysis of the TDELF images, where the ionization occurs from the lone pair or bond regions of the corresponding molecules.

Kohn-Sham self-interaction correction in real time

We present a solution scheme for the time-dependent Kohn-Sham self-interaction correction. Based on the generalized optimized effective potential approach, the multiplicative Kohn-Sham potential is constructed in real time and real space for the self-interaction corrected local density approximation. Excitations of different character, including charge-transfer excitations that had been regarded as prime examples for the failure of standard time-dependent density functionals, are described correctly by this approach. We analyze the time-dependent exchange-correlation potential and density, revealing features that are decisive for the correct description of the response.

Progress in time-dependent density-functional theory

The classic density-functional theory (DFT) formalism introduced by Hohenberg, Kohn, and Sham in the mid-1960s is based on the idea that the complicated N-electron wave function can be replaced with the mathematically simpler 1-electron charge density in electronic structure calculations of the ground stationary state. As such, ordinary DFT cannot treat time-dependent (TD) problems nor describe excited electronic states. In 1984, Runge and Gross proved a theorem making TD-DFT formally exact. Information about electronic excited states may be obtained from this theory through the linear response (LR) theory formalism. Beginning in the mid-1990s, LR-TD-DFT became increasingly popular for calculating absorption and other spectra of medium- and large-sized molecules. Its ease of use and relatively good accuracy has now brought LR-TD-DFT to the forefront for this type of application. As the number and the diversity of applications of TD-DFT have grown, so too has our understanding of the strengths and weaknesses of the approximate functionals commonly used for TD-DFT. The objective of this article is to continue where a previous review of TD-DFT in Volume 55 of the Annual Review of Physical Chemistry left off and highlight some of the problems and solutions from the point of view of applied physical chemistry. Because doubly-excited states have a particularly important role to play in bond dissociation and formation in both thermal and photochemistry, particular emphasis is placed on the problem of going beyond or around the TD-DFT adiabatic approximation, which limits TD-DFT calculations to nominally singly-excited states.

REAL-TIME APPROACH TO TIME-DEPENDENT DENSITY-FUNCTIONAL THEORY IN THE LINEAR AND NONLINEAR REGIME

The linear and nonlinear response of Si 4 and Na 4 to an external perturbation is investigated in the framework of time-dependent density-functional theory. The time-dependent Kohn–Sham equations, which are the central equations in this approach, are solved in real space and real time. A parallelized implementation to solve these nonlinear, one-particle Schrödinger equations is presented. In contrast to Na 4, Si 4 shows high-harmonic generation far beyond the cut-off predicted by the quasiclassical model and predictions for extended systems.

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.

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.

Time-dependent density-functional theory with a self-interaction correction

We discuss an implementation of the self-interaction correction for the local-density approximation to time-dependent density-functional theory. A variational formulation is given, taking care of the necessary constraints. A manageable and transparent propagation scheme using two sets of wave functions is proposed and applied to laser excitation with subsequent ionization of a dimer molecule.

Adiabatic approximation in nonperturbative time-dependent density-functional theory

We construct the exact exchange-correlation potential of time-dependent density-functional theory and the approximation to it that is adiabatic but exact otherwise. For the strong-field double ionization of the Helium atom these two potentials are virtually identical. Thus, memory effects play a negligible role in this paradigm process of nonlinear, nonperturbative electron dynamics. We identify the regime of high-frequency excitations where the adiabatic approximation breaks down and explicitly calculate the nonadiabatic contribution to the exchange-correlation potential.