Density functional theory (DFT) has greatly expanded our ability to affordably compute and understand electronic ground states, by replacing intractable ab initio calculations by models based on paradigmatic physics from high- and low-density limits. But, a comparable treatment of excited states lags behind. Here, we solve this outstanding problem by employing a generalization of density functional theory to ensemble states (EDFT). We thus address important paradigmatic cases of all electronic systems in strongly (low-density) and weakly (high-density) correlated regimes. We show that the high-density limit connects to recent, exactly solvable EDFT results. The low-density limit reveals an unnoticed and most unexpected result-density functionals for strictly correlated ground states can be reused directly for excited states. Nontrivial dependence on excitation structure only shows up at third leading order. Overall, our results provide foundations for effective models of excited states that interpolate between exact low- and high-density limits, which we illustrate on the cases of singlet-singlet excitations in H_{2} and a ring of quantum wells.
Density functional theory (DFT) has greatly expanded our ability to affordably compute and understand electronic ground states, by replacing intractable ab initio calculations by models based on paradigmatic physics from high- and low-density limits. But, a comparable treatment of excited states lags behind. Here, we solve this outstanding problem by employing a generalization of density functional theory to ensemble states (EDFT). We thus address important paradigmatic cases of all electronic systems in strongly (low-density) and weakly (high-density) correlated regimes. We show that the high-density limit connects to recent, exactly solvable EDFT results. The low-density limit reveals an unnoticed and most unexpected result-density functionals for strictly correlated ground states can be reused directly for excited states. Nontrivial dependence on excitation structure only shows up at third leading order. Overall, our results provide foundations for effective models of excited states that interpolate between exact low- and high-density limits, which we illustrate on the cases of singlet-singlet excitations in H_{2} and a ring of quantum wells.
Ensemble density functional theory (eDFT) is an exact time-independent alternative to time-dependent DFT (TD-DFT) for the calculation of excitation energies. Despite its formal simplicity and advantages in contrast to TD-DFT (multiple excitations, for example, can be easily taken into account in an ensemble), eDFT is not standard, which is essentially due to the lack of reliable approximate exchange-correlation (xc) functionals for ensembles. Following Smith et al. [Phys. Rev. B 93, 245131 (2016)], we propose in this work to construct an exact eDFT for the nontrivial asymmetric Hubbard dimer, thus providing more insight into the weight dependence of the ensemble xc energy in various correlation regimes. For that purpose, an exact analytical expression for the weight-dependent ensemble exchange energy has been derived. The complementary exact ensemble correlation energy has been computed by means of Legendre-Fenchel transforms. Interesting features like discontinuities in the ensemble xc potential in the strongly correlated limit have been rationalized by means of a generalized adiabatic connection formalism. Finally, functional-driven errors induced by ground-state density-functional approximations have been studied. In the strictly symmetric case or in the weakly correlated regime, combining ensemble exact exchange with ground-state correlation functionals gives better ensemble energies than when calculated with the ground-state exchange-correlation functional. However, when approaching the asymmetric equiensemble in the strongly correlated regime, the former approximation leads to highly curved ensemble energies with negative slope which is unphysical. Using both ground-state exchange and correlation functionals gives much better results in that case. In fact, exact ensemble energies are almost recovered in some density domains. The analysis of density-driven errors is left for future work.
Ensemble density functional theory (EDFT) offers a promising path for computing excitation energies via a theory of the ensemble density --- a weighted sum of interacting excited-state densities. Recovering excitation energies is theoretically simple, but practical application requires approximations dependent on weight choice. Here, the authors derive exact conditions for EDFT. Using an exactly solvable model, they illustrate violations of exact conditions, weight-dependent derivative discontinuities in the strong-interaction limit, and ultimately the significance of weight dependence for future functional design.
Time-dependent density functional theory (TDDFT) is a standard approach for calculating optical excitations of molecules and solids, while ensemble DFT (EDFT) is a promising alternative under development. We introduce ensemble TDDFT (ETDDFT), a practical theory that combines the two, generalizing both; we ensemble generalize the Gross-Kohn equation and the exchange-correlation kernel of TDDFT, and generalize EDFT to time-dependent problems. We relate coordinate scaling to the adiabatic connection. The new theory provides multiple avenues for constructing and using approximations. We illustrate these on the two-site Hubbard model. We connect our results to the more general case of nonperturbative time dependence.
The recent use of a new ensemble in density functional theory (DFT) to produce direct corrections to the Kohn-Sham transitions yields the elusive double excitations that are missed by time-dependent DFT (TDDFT) with the standard adiabatic approximation. But accuracies are lower than for single excitations, and formal arguments about TDDFT suggest that a correction kernel is needed. In principle, ensemble DFT with direct corrections at the exchange level must yield accurate doubles in the weakly correlated limit. We illustrate with exact calculations and analytic results on the Hubbard dimer. We also explain the error in formal arguments in TDDFT.
Providing sufficient accuracy at low-computational cost, not surprisingly, density-functional theory ( DFT ) is a standard tool in theoretical material science. Similarly, quantum-electrodynamical density-functional theory ( QEDFT ) provides a promising route to describe experiments at the interface of quantum optics and quantum chemistry. Formally, in the Kohn-Sham ( KS ) formulation of DFT and QEDFT the exact many-particle Schrödinger equation is replaced by a set of coupled non-linear one-particle differential equations. The electron-electron, the electron-ion and the electron-photon interactions of the interacting system are mimicked by their corresponding effective KS potential with an, in principle unknown exchange-correlation ( xc ) term. In practice, DFT and QEDFT require approximations for this xc potential. The quality of which fully determines how much the DFT and QEDFT observables deviate from the exact ones. Understanding the shortcomings and improving the quality of approximate functionals makes the study of the exact functional of model systems indispensable. Here, we employ exact diagonalization for the Hamiltonian of different types of real-space lattice models in one spatial dimension to access the ground-state wave function, and thus to explicitly construct the exact functionals of DFT and QEDFT . To study the functional behavior of purely electronic systems, we consider two singlet-state electrons moving in an external potential with two connected atomic wells of variable depth and distance. To study coupled electron-photon systems, we consider the Rabi Hamiltonian with variable electron-photon coupling strength. By varying the model parameters, we systematically construct the functionals for systems from weak to strong coupling between the building blocks of the global system. Using the concept of the effective density-to- potential map, we discuss the intrasystem steepening which is an exact condition on the xc functional that captures the static correlation present in the system and carve out how it relates to the well-known intersystem derivative discontinuity. We further construct arbitrary observables as functionals of the internal variables which are directly affected by such features of the effective density-to-potential map. By constructing the reduced dipole map in the four-particle case of coupled electron-photon system, we show that the density-to-potential map helps to identify and understand recently discovered physical effects in the field of strongly coupled electron-photon systems. This work provides a contribution to the development of approximate function- als required in practical density-functional calculations of arbitrary experimentally accessible observables. Both the intrasystem steepening and the intersystem derivative discontinuity play a crucial role for systems with highly-localized electrons. Including both features in approximate functionals will improve the accu- racy of DFT evaluation of observables for such systems. Those results provide novel insight into the properties of the exact xc potentials in QEDFT . This work opens a new route to construct new functionals using coarse-grained or embed- ding approaches within the general and exact frameworks of DFT and QEDFT.