Time-dependent optimized effective potential

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.

Molecular excitations from meta-generalized gradient approximations in the Kohn-Sham scheme

Meta-Generalized Gradient Approximations (meta-GGAs) can, in principle, include spatial and temporal nonlocality in time-dependent density functional theory at a much lower computational cost than functionals that use exact exchange. We here test whether a meta-GGA that has recently been developed with a focus on capturing nonlocal response properties and the particle number discontinuity can realize such features in practice. To this end, we extended the frequency-dependent Sternheimer formalism to the meta-GGA case. Using the Krieger-Li-Iafrate (KLI) approximation, we calculate the optical response for the selected paradigm molecular systems and compare the meta-GGA Kohn-Sham response to the one found with exact exchange and conventional (semi-)local functionals. We find that the new meta-GGA captures important properties of the nonlocal exchange response. The KLI approximation, however, emerges as a limiting factor in the evaluation of charge-transfer excitations.

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.

Many-Body Quantum Dynamics by the Reduced Density Matrix Based on Time-Dependent Density-Functional Theory

We evaluate the density matrix of an arbitrary quantum mechanical system in terms of the quantities pertinent to the solution of the time-dependent density functional theory (TDDFT) problem. Our theory utilizes the adiabatic connection perturbation method of Görling and Levy, from which the expansion of the many-body density matrix in powers of the coupling constant λ naturally arises. We then find the reduced density matrix ρ_{λ}(r,r^{'},t), which, by construction, has the λ independent diagonal elements ρ_{λ}(r,r,t)=n(r,t), n(r,t) being the particle density. The off-diagonal elements of ρ_{λ}(r,r^{'},t) contribute importantly to the processes unaccessible via the density, directly or by the use of the known TDDFT functionals. Of those, we consider the momentum-resolved photoemission, doing this to the first order in λ, i.e., on the level of the exact exchange theory. In illustrative calculations of photoemission from the quasi-2D electron gas and isolated atoms, we find quantitatively strong and conceptually far-reaching differences with the independent-particle Fermi's golden rule formula.

A theoretical study of Ar8+-acetylene collisions at 1.2 MeV: Ionization and dissociation dynamics

We theoretically study Ar⁸⁺-induced dissociation of C₂H₂ molecule at 1.2 MeV using the time-dependent density-functional theory non-adiabatically coupled to nuclear dynamics. We find that molecular dissociation depends strongly on the ionization at the initial stage and the collision configuration. A detailed analysis shows a correspondence between the charge state of [C₂H₂]^q+ and the final fragments. A remarkable impact parameter effect provides deep insights of bond breakup and electronic transport. We analyze two typical sequential dissociation channels reported in experiments by tracking structural and electronic dynamics in real time. Our results provide better understanding of experiments. Moreover, the comparison between various exchange-correlation functionals reveals that electrons' correlation and self-interaction do not significantly impact the initial ionization and fragment distribution in the present study.

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.

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.

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.

Time-dependent density functional theory for many-electron systems interacting with cavity photons

Time-dependent (current) density functional theory for many-electron systems strongly coupled to quantized electromagnetic modes of a microcavity is proposed. It is shown that the electron-photon wave function is a unique functional of the electronic (current) density and the expectation values of photonic coordinates. The Kohn-Sham system is constructed, which allows us to calculate the above basic variables by solving self-consistent equations for noninteracting particles. We suggest possible approximations for the exchange-correlation potentials and discuss implications of this approach for the theory of open quantum systems. In particular we show that it naturally leads to time-dependent density functional theory for systems coupled to the Caldeira-Leggett bath.

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.

FIRST PRINCIPLES MOLECULAR DYNAMICS INVOLVING EXCITED STATES AND NONADIABATIC TRANSITIONS

Extensions of traditional molecular dynamics to excited electronic states and non-Born–Oppenheimer dynamics are reviewed focusing on applicability to chemical reactions of large molecules, possibly in condensed phases. The latter imposes restrictions on both the level of accuracy of the underlying electronic structure theory and the treatment of nonadiabaticity. This review, therefore, exclusively deals with ab initio "on the fly" molecular dynamics methods. For the same reason, mainly mixed quantum-classical approaches to nonadiabatic dynamics are considered.

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.

An introduction to “Computational Crystallography”

The currently available methods for the computation of structures and their properties are reviewed. After a brief introduction into some common technical aspects, the capabilities and limitations of the most commonly used approaches are discussed. Examples are given to show the state of the art in Computational “Crystallography”, and possible future developments are outlined

Self-interaction Free Relativistic Spin-density Functional Theory

We review the progress in the recent development of the relativistic optimized effective potential (ROEP) method within spin-density functional theory. The ROEP-equations for spin-polarized systems are derived and their application to open-shell atoms using the exact exchange approximation for the exchange correlation-functional is presented. Further, we expand the ROEP framework to treat extended systems within the KKR-multiple scattering formalism. We illustrate the application of the theory to open-shell free atoms and solids.

Time-dependent density-functional theory and strongly correlated systems: insight from numerical studies

We illustrate the scope of time-dependent density-functional theory for strongly correlated (lattice) models out of equilibrium. Using the exact many-body time evolution, we reverse engineer the exact exchange correlation (xc) potential v_(xc) for small Hubbard chains exposed to time-dependent fields. We introduce an adiabatic local density approximation to v_(xc) for the 1D Hubbard model and compare it to exact results, to gain insight about approximate xc potentials. Finally, we provide some remarks on the v-representability for the 1D Hubbard model.

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.

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

The exact exchange potential in time-dependent density-functional theory is defined as an orbital functional through the time-dependent optimized effective potential (TDOEP) method. We numerically solve the TDOEP integral equation for the real-time nonlinear intersubband electron dynamics in a semiconductor quantum well with two occupied subbands. It is found that memory effects become significant in the vicinity of intersubband resonances.

Time-dependent density functional theory: past, present, and future

Time-dependent density functional theory (TDDFT) is presently enjoying enormous popularity in quantum chemistry, as a useful tool for extracting electronic excited state energies. This article discusses how TDDFT is much broader in scope, and yields predictions for many more properties. We discuss some of the challenges involved in making accurate predictions for these properties.

Time-dependent density-functional theory beyond the adiabatic approximation: Insights from a two-electron model system

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.

Linear response time-dependent density functional theory for van der Waals coefficients

A linear response time-dependent density functional theory is described and used to calculate the dynamic polarizabilities and van der Waals C(6) coefficients of complex atom pairs. We present values of C(6) for dimers of main group atoms and the first row of transition metal atoms.

Protein-bound chromophores astaxanthin and phytochromobilin: excited state quantum chemical studies

We present an overview of excited state quantum chemical calculations aimed at elucidating controversial issues regarding the photochemistry of the protein-bound chromophores astaxanthin and phytochromobilin. In particular, we show how the application of time-dependent density functional theory and other single-reference quantum chemical excited state methods have contributed to shed new light on the origin of the >0.5 eV bathochromic shift of the electronic absorption by the carotenoid astaxanthin in the protein macromolecular complex crustacyanin, and the mechanism for C15-Z,syn --> C15-E,anti isomerization of the tetrapyrrole phytochromobilin that underlies the photoactivation of the plant photoreceptor phytochrome. Within the approximation that exciton coupling is neglected, the calculations on astaxanthin provide support for the notion that the bathochromic shift, which is responsible for the slate-blue coloration of lobster shell, is due to polarization rather than a conformational change of the chromophore in the protein-bound state. Furthermore, the polarization is attributed to a hydrogen-bonded protonated histidine residue. The calculations on phytochromobilin, in turn, suggest that a stepwise C15-Z,syn --> C15-E,syn (photochemical), C15-E,syn --> C15-E,anti (thermal) mechanism is much more favorable than a concerted, fully photochemical mechanism, and that neutral forms of the chromophore are much less likely to photoisomerize than the parent, protonated form. Accordingly, the calculations indirectly support the view that the photoactivation of phytochrome does not involve a proton transfer from the chromophore to the surrounding protein.

Derivative discontinuities in time-dependent density-functional theory

Based on the Runge-Gross theorem for ensembles we investigate the influence of particle-number-changes on the exchange-correlation potential in time-dependent density-functional theory. We show that the potential changes discontinuously when the particle number crosses an integer value. Real-time simulations of an atomic ionization process demonstrate that this discontinuity appears naturally in the theory of the time-dependent optimized effective potential. The importance of such a discontinuity for physical processes, even such ones in which the particle number is a constant, is discussed.