Two particular subjects tackled in the recent paper by Schirmer and Dreuw (SD) [J. Schirmer and A. Dreuw, Phys. Rev. A 75, 022513 (2007)] are commented on by us. By providing a convincing example, we answer positively the question posed by SD whether a local operator can be reconstructed from the knowledge of its particle-hole matrix elements (when the number of particle states exceeds one). In connection with the problem of the utmost reduction of the density-functional theory (DFT) to a solution of one equation for the density amplitude---the radical Kohn-Sham (rKS) approach---used by SD as a valuable tool to investigate the time-independent and time-dependent DFT, we demonstrate that the calculational scheme implementing the rKS formulation is rather uncompetitive with the traditional KS scheme.
Two particular subjects tackled in the recent paper by Schirmer and Dreuw (SD) [J. Schirmer and A. Dreuw, Phys. Rev. A 75, 022513 (2007)] are commented on by us. By providing a convincing example, we answer positively the question posed by SD whether a local operator can be reconstructed from the knowledge of its particle-hole matrix elements (when the number of particle states exceeds one). In connection with the problem of the utmost reduction of the density-functional theory (DFT) to a solution of one equation for the density amplitude---the radical Kohn-Sham (rKS) approach---used by SD as a valuable tool to investigate the time-independent and time-dependent DFT, we demonstrate that the calculational scheme implementing the rKS formulation is rather uncompetitive with the traditional KS scheme.
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.
In the preceding Comment, Maitra, van Leeuwen, and Burke (MLB) [Phys. Rev. A 78, 056501 (2008) attempt to refute our criticism of the foundations of time-dependent density functional theory [see J. Schirmer and A. Dreuw, Phys. Rev. A 75, 022513 (2007)]. However, their arguments miss the essence of our position. This is partly due to an ambiguity concerning the present status of the time-dependent (TD) Kohn-Sham (KS) equations. Based on the first Runge-Gross (mapping) theorem, the TD KS equations have to be established as a fixed-point iteration scheme for potential functionals (PF-FPI). Here, however, the convergence of the FPI is not assured. On the other hand, time propagation of the TD KS equations in their PF-FPI form fails because these equations do not represent physical equations of motion. MLB present a construction of the exact exchange-correlation potential, apparently contradicting one of our findings. We explain that their construction is explicitly based on the solution of the TD Schr\"odinger equation and, for this reason, does not address our objections, let alone refute them.
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.
Time-dependent [current]-density functional theory for open quantum systems (OQS) has emerged as a formalism that can incorporate dissipative effects in the dynamics of many-body quantum systems. Here, we review and clarify some formal aspects of these theories that have been recently questioned in the literature. In particular, we provide theoretical support for the following conclusions: (1) contrary to what we and others had stated before, within the master equation framework, there is in fact a one-to-one mapping between vector potentials and current densities for fixed initial state, particle-particle interaction, and memory kernel; (2) regardless of the first conclusion, all of our recently suggested Kohn-Sham (KS) schemes to reproduce the current and particle densities of the original OQS, and in particular, the use of a KS closed driven system, remains formally valid; (3) the Lindblad master equation maintains the positivity of the density matrix regardless of the time-dependence of the Hamiltonian or the dissipation operators; (4) within the stochastic Schrödinger equation picture, a one-to-one mapping from stochastic vector potential to stochastic current density for individual trajectories has not been proven so far, except in the case where the vector potential is the same for every member of the ensemble, in which case, it reduces to the Lindblad master equation picture; (5) master equations may violate certain desired properties of the density matrix, such as positivity, but they remain as one of the most useful constructs to study OQS when the environment is not easily incorporated explicitly in the calculation. The conclusions support our previous work as formally rigorous, offer new insights into it, and provide a common ground to discuss related theories.
Given the time-evolution of an electron charge density, the local potential in Kohn-Sham time-dependent density functional theory (KS-TDDFT) can be modeled as a sum of instantaneous and dynamic contributions by assuming a certain form of the time-dependent Kohn-Sham (TD-KS) orbitals. The instantaneous part is obtained numerically using methods from ground-state density functional theory (DFT) and the dynamic part is expressed in terms of a velocity potential that depends on the electron current density. The suggested form of the TD-KS orbitals satisfies several known constraints (orthonormality, N-representability, J-representability), and the domain of validity is shown to depend on the evolution of the instantaneous quantities. Through this decomposition, we can relate time-dependent and ground -state V-representability. The resulting potentials are shown numerically to approximate the exact time-dependent Kohn-Sham potentials for a set of 3 non-singlet two-particle systems (a Kohn-mode, a Coulomb explosion, and a double quantum well) where the exact solutions and reference densities are known or obtained through configuration interaction (CI) approaches.
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.