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 upon the idea that the complicated N-electron wavefunction can be replaced with the mathematically simpler 1-electron charge density in electronic struc- ture calculations of the ground stationary state. As such, ordinary DFT is neither able to 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. Begin- ning 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 has grown, so too has grown 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 this series [Annu. Rev. Phys. Chem. 55: 427 (2004)] left off and highlight some of the problems and solutions from the point of view of applied physical chemistry. Since doubly-excited states have a particularly important role to play in bond dissociation and formation in both thermal and photochemistry, particular emphasis will be placed upon the problem of going beyond or around the TD-DFT adiabatic approximation which limits TD-DFT calculations to nominally singly-excited states. Posted with permission from the Annual Review of Physical Chemistry, Volume 63 \c{opyright} 2012 by Annual Reviews, http://www.annualreviews.org

Fast computation of the Kohn-Sham susceptibility of large systems

For hybrid systems, such as molecules grafted onto solid surfaces, the calculation of linear response in time dependent density functional theory is slowed down by the need to calculate, in N^4 operations, the susceptibility of N non interacting Kohn-Sham reference electrons. We show how this susceptibility can be calculated N times faster within finite precision. By itself or in combination with previous methods, this should facilitate the calculation of TDDFT response and optical spectra of hybrid systems.Comment: submitted 25/1/200

Rydberg transition frequencies from the Local Density Approximation

A method is given that extracts accurate Rydberg excitations from LDA density functional calculations, despite the short-ranged potential. For the case of He and Ne, the asymptotic quantum defects predicted by LDA are in less than 5% error, yielding transition frequency errors of less than 0.1eV.Comment: 4 pages, 6 figures, submitted to Phys. Rev. Let

A joint time-dependent density-functional theory for excited states of electronic systems in solution

We present a novel joint time-dependent density-functional theory for the description of solute-solvent systems in time-dependent external potentials. Starting with the exact quantum-mechanical action functional for both electrons and nuclei, we systematically eliminate solvent degrees of freedom and thus arrive at coarse-grained action functionals which retain the highly accurate \emph{ab initio} description for the solute and are, in principle, exact. This procedure allows us to examine approximations underlying popular embedding theories for excited states. Finally, we introduce a novel approximate action functional for the solute-water system and compute the solvato-chromic shift of the lowest singlet excited state of formaldehyde in aqueous solution, which is in good agreement with experimental findings.Comment: 11 page

A new and efficient approach to time-dependent density-functional perturbation theory for optical spectroscopy

Using a super-operator formulation of linearized time-dependent density-functional theory, the dynamical polarizability of a system of interacting electrons is given a matrix continued-fraction representation whose coefficients can be obtained from the non-symmetric block-Lanczos method. The resulting algorithm allows for the calculation of the {\em full spectrum} of a system with a computational workload which is only a few times larger than that needed for {\em static} polarizabilities within time-independent density-functional perturbation theory. The method is demonstrated with the calculation of the spectrum of benzene, and prospects for its application to the large-scale calculation of optical spectra are discussed.Comment: 4 pages, 2 figure

Spin gaps and spin-flip energies in density-functional theory

Energy gaps are crucial aspects of the electronic structure of finite and extended systems. Whereas much is known about how to define and calculate charge gaps in density-functional theory (DFT), and about the relation between these gaps and derivative discontinuities of the exchange-correlation functional, much less is know about spin gaps. In this paper we give density-functional definitions of spin-conserving gaps, spin-flip gaps and the spin stiffness in terms of many-body energies and in terms of single-particle (Kohn-Sham) energies. Our definitions are as analogous as possible to those commonly made in the charge case, but important differences between spin and charge gaps emerge already on the single-particle level because unlike the fundamental charge gap spin gaps involve excited-state energies. Kohn-Sham and many-body spin gaps are predicted to differ, and the difference is related to derivative discontinuities that are similar to, but distinct from, those usually considered in the case of charge gaps. Both ensemble DFT and time-dependent DFT (TDDFT) can be used to calculate these spin discontinuities from a suitable functional. We illustrate our findings by evaluating our definitions for the Lithium atom, for which we calculate spin gaps and spin discontinuities by making use of near-exact Kohn-Sham eigenvalues and, independently, from the single-pole approximation to TDDFT. The many-body corrections to the Kohn-Sham spin gaps are found to be negative, i.e., single particle calculations tend to overestimate spin gaps while they underestimate charge gaps.Comment: 11 pages, 1 figure, 3 table

Time-dependent Density Functional calculation of e-H scattering

Phase shifts for single-channel elastic electron-atom scattering are derived from time-dependent density functional theory. The H$^-$ ion is placed in a spherical box, its discrete spectrum found, and phase shifts deduced. Exact-exchange yields an excellent approximation to the ground-state Kohn-Sham potential, while the adiabatic local density approximation yields good singlet and triplet phase shifts.Comment: 5 pages, 4 figures, 1 tabl