We derive the basic formalism of density functional theory for time-dependent electron-nuclear systems. The basic variables of this theory are the electron density in body-fixed frame coordinates and the diagonal of the nuclear $N$-body density matrix. The body-fixed frame transformation is carried out in order to achieve an electron density that reflects the internal symmetry of the system. We discuss the implications of this body-fixed frame transformation and establish a Runge-Gross-type theorem and derive Kohn-Sham equations for the electrons and nuclei. We illustrate the formalism by performing calculations on a one-dimensional diatomic molecule for which the many-body Schr\"odinger equation can be solved numerically. These benchmark results are then compared to the solution of the time-dependent Kohn-Sham equations in the Hartree approximation. Furthermore, we analyze the excitation energies obtained from the linear response formalism in the single pole approximation. We find that there is a clear need for improved functionals that go beyond the simple Hartree approximation.
We derive the basic formalism of density functional theory for time-dependent electron-nuclear systems. The basic variables of this theory are the electron density in body-fixed frame coordinates and the diagonal of the nuclear N-body density matrix. The body-fixed frame transformation is carried out in order to achieve an electron density that reflects the internal symmetry of the system. We discuss the implications of this body-fixed frame transformation and establish a Runge-Gross-type theorem and derive Kohn-Sham equations for the electrons and nuclei. We illustrate the formalism by performing calculations on a one-dimensional diatomic molecule for which the many-body Schr"odinger equation can be solved numerically. These benchmark results are then compared to the solution of the time-dependent Kohn-Sham equations in the Hartree approximation. Furthermore, we analyze the excitation energies obtained from the linear response formalism in the single pole approximation. We find that there is a clear need for improved functionals that go beyond the simple Hartree approximation.
For the calculation of neutral excitations, time-dependent density functional theory (TDDFT) is an exact reformulation of the many-body time-dependent Schrödinger equation, based on knowledge of the density instead of the many-body wavefunction. The density can be determined in an efficient scheme by solving one-particle non-interacting Schrödinger equations—the Kohn–Sham equations. The complication of the problem is hidden in the—unknown—time-dependent exchange and correlation potential that appears in the Kohn–Sham equations and for which it is essential to find good approximations. Many approximations have been suggested and tested for finite systems, where even the very simple adiabatic local-density approximation (ALDA) has often proved to be successful. In the case of solids, ALDA fails to reproduce optical absorption spectra, which are instead well described by solving the Bethe–Salpeter equation of many-body perturbation theory (MBPT). On the other hand, ALDA can lead to excellent results for loss functions (at vanishing and finite momentum transfer). In view of this and thanks to recent successful developments of improved linear-response kernels derived from MBPT, TDDFT is today considered a promising alternative to MBPT for the calculation of electronic spectra, even for solids. After reviewing the fundamentals of TDDFT within linear response, we discuss different approaches and a variety of applications to extended systems.
Time-dependent density functional theory (TDDFT) provides a successful approach to calculate excitation energies of atomic and molecular systems. In part I of this work we present a double-pole approximation (DPA) to the response equations of TDDFT. The double-pole approximation provides an exact description of systems with two strongly coupled excitations which are isolated from the rest of the spectrum. In contrast to the traditional single- pole approximation of TDDFT the DPA also yields corrections to the Kohn-Sham oscillator strengths. Several critical pole separations can be identified, e.g. we find that the pole coupling can cause transitions to vanish entirely from the optical spectrum. We also demonstrate how to invert the double-pole solution which allows us to predict matrix elements of the exchange- correlation kernel fxc from experimental input. This can serve as benchmark for the construction of future approximations for the kernel fxc. Reduced density matrix functional theory (RDMFT) has emerged recently as promising candidate to treat strongly correlated electronic many-body systems beyond traditional density functional theory (DFT). The research within RDMFT was so far focussed on the static theory. In this work we attempt some first steps towards a time-dependent generalization of RDMFT. In part II we derive equations of motion for natural orbitals and occupation numbers. Using the equation of motion for the occupation numbers we show that an adiabatic extension of presently known ground-state functionals of static RDMFT always leads to occupation numbers which are constant in time. From the stationary conditions of the equations of motion for the N-body correlations (correlated parts of the N-body matrices) we derive a new class of ground-state functionals which can be used in static RDMFT. Applications are presented for a one-dimensional model system where the time-dependent many-body Schrödinger equation can be propagated numerically. We use optimal control theory to find optimized laser pulses for transitions in a model for atomic Helium. From the numerically exact correlated wavefunction we extract the exact time evolution of natural orbitals and occupation numbers for (i) laser-driven Helium and (ii) electron-ion scattering. Part III of this work considers time-dependent quantum transport within TDDFT. We present an algorithm for the calculation of extended eigenstates of single- particle Hamiltonians which is especially tailored to a finite-difference discretization of the Schrödinger equation. We consider the propagation of finite mesoscopic systems and demonstrate the limitations of such an approach. To overcome the shortcomings of a description of quantum transport in terms of finite systems we develop a time-propagation scheme for extended states which utilizes a mixed basis representation. Our discretization scheme allows to treat central device and lead regions on the same footing thus preventing artificial reflections at grid boundaries.
Since the early days of quantum mechanics, the wave function theory has proven to be very successful in describing many different quantum processes and phenomena. However, in many problems of quantum chemistry and solid-state physics, where the dimensionality of the systems studied is relatively high, ab initio calculations of the structure of atoms, molecules, clusters, and crystals, and their interactions are very often prohibitive. Hence, alternative formulations based on the direct use of the probability density, gathered under what is generally known as the density matrix theory [1], were also developed since the very beginning of the new mechanics. The independent electron approximation or Thomas-Fermi model, and the Hartree and Hartree-Fock approaches are former statistical models developed in that direction [2]. These models can be considered direct predecessors of the more recent density functional theory (DFT) [3], whose principles were established by Hohenberg, Kohn, and Sham [4,5] in the mid-1960s. According to this theory, the fundamental physical information about a many-body system is provided by single-particle densities in a three-dimensional space, which are obtained variationally within a time-independent framework. When compared with other previous formalisms, DFT presents two clear advantages: (1) it is able to treat many-body problems in a sufficiently accurate way and (2) it is computationally simple. This explains why it is one of the most widely used theories to deal with electronic structure-the electronic ground-state energy as a function of the position of the atomic nuclei determines molecular structures and solids, providing at the same time the forces acting on the atomic nuclei when they are not at their equilibrium positions. At present, DFT is used routinely to solve many problems in gas phase and condensed matter. Furthermore, it has made possible the development of accurate molecular dynamics schemes in which the forces are evaluated quantum mechanically ‘‘on the fly.’’ Nonetheless, DFT is a fundamental tool provided the systems studied are relatively large; for small systems, standard methods based on the use of the wave function render quite accurate results [6]. Moreover, it is also worth stressing that all practical applications of DFT rely on essentially uncontrolled approximations [7] (e.g., the local density approximation [4,5], the local spin-density approximation, or generalized gradient approximations [8]), and therefore the validity of DFT is conditioned to its capability to provide fairly good values of the experimental data.
In the present thesis, a rigorous mathematical formulation of time-dependent density-functional theory for lattice systems is derived, serving a formally exact approach towards solving time-dependent many-particle Schrödinger problems. After introducing the necessary mathematical foundations, a one-to-one mapping from external potentials onto electronic densities obtained by solving the Schrödinger equation is introduced. Its properties and the related Banach spaces of scalar potentials and electronic densities are discussed in detail. It allows to precisely formulate the Kohn-Sham Iteration scheme, mapping a given many-particle to an effective single-particle Schrödinger problem, both generating the exact same density. Employing the Banach fixed point theorem, the Kohn-Sham Iteration scheme is proven to be convergent depending on the chosen initial density.
The general expectation that, in principle, the time-dependent density-functional theory (TDDFT) is an exact formulation of the time evolution of an interacting N-electron system is critically reexamined. It is demonstrated that the previous TDDFT foundation, resting on four theorems by Runge and Gross (RG) [Phys. Rev. Lett. 52, 997 (1984)], is invalid because undefined phase factors corrupt the RG action integral functionals. Our finding confirms much of a previous analysis by van Leeuwen [Int. J. Mod. Phys. B 15, 1969 (2001)]. To analyze the RG theorems and other aspects of TDDFT, an utmost simplification of the Kohn-Sham (KS) concept has been introduced, in which the ground-state density is obtained from a single KS equation for one spatial (spinless) orbital. The time-dependent (TD) form of this radical Kohn-Sham (rKS) scheme, which has the same validity status as the ordinary KS version, has proved to be a valuable tool for analysis. The rKS concept is used to clarify also the alternative nonvariational formulation of TD KS theory. We argue that it is just a formal theory, allowing one to reproduce but not predict the time development of the exact density of the interacting N-electron system. Besides the issue of the formal exactness of TDDFT, it is shown that both the static and time-dependent KS linear response equations neglect the particle-particle (p-p) and hole-hole (h-h) matrix elements of the perturbing operator. For a local (multiplicative) operator this does not lead to a loss of information due to a remarkable general property of local operators. Accordingly, no logical inconsistency arises with respect to DFT, because DFT requires any external potential to be local. For a general nonlocal operator the error resulting from the neglected matrix elements is of second order in the electronic repulsion.