Solving the ultranonlocality problem in time-dependent spin-density-functional theory

It has been known for some time that the exchange-correlation potential in time-dependent density-functional theory is an intrinsically nonlocal functional of the density as soon as one goes beyond the adiabatic approximation. In this paper we show that a much more severe nonlocality problem, with a completely different physical origin, plagues the exchange-correlation potentials in time-dependent spin-density functional theory. We show how the use of the spin current density as a basic variable solves this problem, and we provide an explicit local expression for the exchange-correlation fields as functionals of the spin currents.     

Time-dependent density-functional theory beyond the local-density approximation

Approximations for the ground-state exchange-correlation potential of density-functional theory have reached a high level of sophistication. By contrast, time- or frequency-dependent exchange-correlation potentials are still being treated in a local approximation. Here we propose a novel approximation scheme, which effectively brings the power of the generalized gradient approximation (GGA) and meta-GGA to time-dependent density-functional theory. The theory should allow a more accurate treatment of strongly inhomogeneous electronic systems (e.g. molecular junctions) while remaining essentially exact for slowly varying densities and slowly varying external potentials.     

Time-dependent generalized Kohn–Sham theory

Generalized Kohn–Sham (GKS) theory extends the realm of density functional theory (DFT) by providing a rigorous basis for non-multiplicative potentials, the use of which is outside original Kohn–Sham theory. GKS theory is of increasing importance as it underlies commonly used approximations, notably (conventional or range-separated) hybrid functionals and meta-generalized-gradient-approximation (meta-GGA) functionals. While this approach is often extended in practice to time-dependent DFT (TDDFT), the theoretical foundation for this extension has been lacking, because the Runge–Gross theorem and the van Leeuwen theorem that serve as the basis of TDDFT have not been generalized to non-multiplicative potentials. Here, we provide the necessary generalization. Specifically, we show that with one simple but non-trivial additional caveat – upholding the continuity equation in the GKS electron gas – the Runge–Gross and van Leeuwen theorems apply to time-dependent GKS theory. We also discuss how this is manifested in common GKS-based approximations.     

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.     

Climbing the density functional ladder: nonempirical meta-generalized gradient approximation designed for molecules and solids

The electron density, its gradient, and the Kohn-Sham orbital kinetic energy density are the local ingredients of a meta-generalized gradient approximation (meta-GGA). We construct a meta-GGA density functional for the exchange-correlation energy that satisfies exact constraints without empirical parameters. The exchange and correlation terms respect two paradigms: one- or two-electron densities and slowly varying densities, and so describe both molecules and solids with high accuracy, as shown by extensive numerical tests. This functional completes the third rung of "Jacob's ladder" of approximations, above the local spin density and GGA rungs.     

Meta-generalized gradient approximation: non-empirical construction and performance of a density functional

The local ingredients of a meta-generalized gradient approximation (meta-GGA) include the electron density, its gradient, and the Kohn–Sham orbital kinetic energy density. We discuss the strategy of constructing a successful meta-GGA density functional for the exchange-correlation energy that satisfies exact constraints without empirical parameters. The new feature of this functional is that it simultaneously respects the two paradigms of electronic structure theory: one- or two-electron densities and slowly-varying densities, and so is uniformly accurate for atoms, molecules and solids. Results of extensive numerical tests of the new functional are summarized and evaluated.     

Time-dependent density-functional theory meets dynamical mean-field theory: real-time dynamics for the 3D Hubbard model

We introduce a new class of exchange-correlation potentials for a static and time-dependent density-functional theory of strongly correlated systems in 3D. The potentials are obtained via dynamical mean-field theory and, for strong enough interactions, exhibit a discontinuity at half-filling density, a signature of the Mott transition. For time-dependent perturbations, the dynamics is described in the adiabatic local density approximation. Results from the new scheme compare very favorably to exact ones in clusters. As an application, we study Bloch oscillations in the 3D Hubbard model.     

Dynamical corrections to the DFT-LDA electron conductance in nanoscale systems

Using time-dependent current-density functional theory, we derive analytically the dynamical exchange-correlation correction to the dc conductance of nanoscale junctions. The correction pertains to the conductance calculated in the zero-frequency limit of time-dependent density functional theory within the adiabatic local-density approximation. In particular, we show that in linear response, the correction depends nonlinearly on the gradient of the electron density; thus, it is more pronounced for molecular junctions than for quantum point contacts. We provide specific numerical examples to illustrate these findings.     

Charge transfer, double and bond-breaking excitations with time-dependent density matrix functional theory

Time-dependent density functional theory (TDDFT) in its current adiabatic implementations exhibits three remarkable failures: (a) completely wrong behavior of the excited state surface along a bond-breaking coordinate; (b) lack of doubly excited configurations; (c) much too low charge transfer excitation energies. These TDDFT failure cases are all strikingly exhibited by prototype two-electron systems such as dissociating H2 and HeH+. We find for these systems with time-dependent density matrix functional theory that: (a) Within previously formulated simple adiabatic approximations, the bonding-to-antibonding excited state surface as well as charge transfer excitations are described without problems, but not the double excitations; (b) An adiabatic approximation is formulated in which also the double excitations are fully accounted for.     

Optical saturation driven by exciton confinement in molecular chains: a time-dependent density-functional theory approach

We identify excitonic confinement in one-dimensional molecular chains (i.e., polyacetylene and H2) as the main driving force for the saturation of the chain polarizability as a function of the number of molecular units. This conclusion is based on first principles time-dependent density-functional theory calculations using a recently developed exchange-correlation kernel that accounts for excitonic effects. The failure of simple local and semilocal functionals is shown to be linked to the lack of memory effects, spatial ultranonlocality, and self-interaction corrections. These effects get smaller as the gap reduces, in which case such simple approximations do perform better.     

Formally exact quantum variational principles for collective motion based on the invariance principle of the Schrödinger equation

Utilizing the concept of invariant collective subspace of a many-body system (or invariance principle of the time-dependent Schrödinger equation), we derive a number of formally exact variational principles to characterize the subspace. Previous studies based on time-dependent or adiabatic time-dependent Hartree-Fock theory are, in principle, contained as approximations.     

Signatures of short-range many-body effects in the dielectric function of silicon for finite momentum transfer

We present an investigation of the dynamic structure factor and of the dielectric function epsilonM(Q,omega) of the prototypical semiconductor silicon for finite momentum transfer, combining inelastic x-ray scattering experiments and ab initio calculations. In contrast with optical spectra, for finite momentum transfer time-dependent density-functional theory in the adiabatic local-density approximation together with lifetime broadening describes the physics of valence excitations correctly. Major structures in the spectra, governed by short-range crystal and exchange-correlation local-field effects, are strongly influenced by a mixing of transitions of positive and negative energies, in striking difference to spectra for vanishing momentum transfer. This mixing gives rise to a pronounced Fano asymmetry.     

Local vertex corrections from exchange-correlation kernels with a discontinuity

The fundamental gap of an interacting many-electron system is given by the sum of the single-particle Kohn-Sham gap and the derivative discontinuity. The latter can be generated by advanced approximations to the exchange-correlation (XC) energy and is the key quantity to capture strong correlation with density functional theory (DFT). In this work we derive an expression for the derivative discontinuity in terms of the XC kernel of time-dependent density functional theory and demonstrate the crucial role of a discontinuity in the XC kernel itself. By relating approximate XC kernels to approximate local vertex corrections we then generate beyond-GW self-energies that include a discontinuity in the local vertex function. The quantitative importance of this result is illustrated with a numerical study of the local exchange vertex on model systems.     

Excitonic effects in solids described by time-dependent density-functional theory

Starting from the many-body Bethe-Salpeter equation we derive an exchange-correlation kernel f(xc) that reproduces excitonic effects in bulk materials within time-dependent density functional theory. The resulting f(xc) accounts for both self-energy corrections and the electron-hole interaction. It is static, nonlocal, and has a long-range Coulomb tail. Taking the example of bulk silicon, we show that the -alpha/q(2) divergency is crucial and can, in the case of continuum excitons, even be sufficient for reproducing the excitonic effects and yielding excellent agreement between the calculated and the experimental absorption spectrum.     

Exact density-functional potentials for time-dependent quasiparticles

We calculate the exact Kohn-Sham potential that describes, within time-dependent density-functional theory, the propagation of an electron quasiparticle wave packet of nonzero crystal momentum added to a ground-state model semiconductor. The potential is observed to have a highly nonlocal functional dependence on the charge density, in both space and time, giving rise to features entirely lacking in local or adiabatic approximations. The dependence of the nonequilibrium part of the Kohn-Sham electric field on the local current and charge density is identified as a key element of the correct Kohn-Sham functional.     

Semiclassical neutral atom as a reference system in density functional theory

We use the asymptotic expansions of the semiclassical neutral atom as a reference system in density functional theory to construct accurate generalized gradient approximations (GGAs) for the exchange-correlation and kinetic energies without any empiricism. These asymptotic functionals are among the most accurate GGAs for molecular systems, perform well for solid state, and overcome current GGA state of the art in frozen density embedding calculations. Our results also provide evidence for the conjointness conjecture between exchange and kinetic energies of atomic systems.     

Unifying description of the optical properties of InN from first principles

The optical properties of hexagonal InN have been studied using the all-electron approach based on density functional theory (DFT). The full-potential augmented plane wave method is employed with two different exchange-correlation potentials, the Perdew–Wang (PW) and the Engel–Vosko (EV) approximations. In addition, both non-relativistic and relativistic approximations are considered. We found that the PW and relativistic approximations give a metallic ground state; whereas using the EV and non-relativistic approximations a semiconductor phase is obtained, opening the gap up to 0.83 eV. Besides, the calculated interband transitions of the complex dielectric function up to 13 eV show favourable agreement with the recent spectroscopic ellipsometry results.     

Nonempirical construction of current-density functionals from conventional density-functional approximations

Density-functional approximations for the exchange-correlation energy Exc[n] of a many-electron ground state are highly developed and widely useful. When a paramagnetic current jp(r) is present, Vignale and Rasolt have extended the Kohn-Sham theorems and presented an additive correction valid to second order in the gauge-invariant vorticity nu=Delta x (jp/n):Exc[n, jp]=Exc[n, jp=0] + DeltaE(VR)(xc)[n, nu]. Apart from spin-polarization effects, their correction is unambiguous for a generalized gradient approximation (GGA). But for a meta-GGA (MGGA), one needs to know how to go back from the orbital kinetic energy density tau([n, jp];r) to tau([n,0];r); we show how to do this here. Numerical tests on the degeneracies for open-shell atoms show that current-density corrections reduce the error of GGA from 2 to 1 kcal/mol, and of MGGA from 5 to 2 kcal/mol.     

Exact exchange-correlation treatment of dissociated H(2) in density functional theory

Whereas the wave function of the dissociated H(2) molecule is easily obtained, an exact treatment of this system has not yet been given in density functional theory. We propose an orbital dependent exchange-correlation functional that yields this exact Kohn-Sham solution. This type of functional offers an alternative to local-density and generalized gradient approximations and the prospects are indicated for treating with it some outstanding problems in density functional theory.