Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2507.19464 · Physical Review Letters

Ensemble Time-Dependent Density Functional Theory

Kimberly J. Daas, Steven Crisostomo, Kieron Burke

Published May 12, 2026• 3 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/q7dh-ccch

Related papers

Accurate double excitations from ensemble density functional calculations

Francisca Sagredo, Kieron Burke

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.

arXiv:1806.03392•The Journal of Chemical Physics•Oct 2, 2018•46 citations

Exact ensemble density functional theory for excited states in a model system: Investigating the weight dependence of the correlation energy

Killian Deur, Laurent Mazouin, Emmanuel Fromager

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.

arXiv:1610.03308•Physical review. B./Physical review. B•Jan 13, 2017•53 citations

Time-dependent density-functional theory for extended systems

Silvana Botti, Arno Schindlmayr, Rodolfo Del Sole +1 more

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.

Reports on Progress in Physics•Feb 13, 2007•237 citations

Thermal Density Functional Theory: Time-Dependent Linear Response and Approximate Functionals from the Fluctuation-Dissipation Theorem

Aurora Pribram−Jones, Paul Grabowski, Kieron Burke

The van Leeuwen proof of linear-response time-dependent density functional theory (TDDFT) is generalized to thermal ensembles. This allows generalization to finite temperatures of the Gross-Kohn relation, the exchange-correlation kernel of TDDFT, and fluctuation dissipation theorem for DFT. This produces a natural method for generating new thermal exchange-correlation approximations.

arXiv:1509.03071•Physical Review Letters•Jun 8, 2016•80 citations

Diagrammatic techniques for time-dependent density-functional theory

Ralf Stubner

The density-functional theory (DFT) is one of the most important ab initio methods for calculating properties of many-body systems. The success of DFT is based on the Kohn-Sham (KS) method and the approximations for the exchange-correlation (xc) effects. These approximations are well understood. However, DFT is first and foremost a theory for the ground state. It is very difficult to obtain properties of excited states. The time-dependent density-functional theory (TDDFT) can overcome this restriction and in the linear-response regime can give access to the properties of excited states. However, the xc effects, which are contained in the xc kernel f_xc , have to be approximated. The approximations for f_xc are not well understood up to now. Green-functions-based conventional many-body theory (MBT) is an alternative approach to the many-body problem. Especially when combined with diagrammatic techniques, MBT allows to derive physically motivated approximations. This work investigates the use of diagrammatic techniques from MBT for TDDFT . The aim of the investigations is to contribute to the development, testing, and understanding of the approximations involved in the time-dependent case. It is shown that it is possible to expand the xc kernel in a diagrammatic series. Explicit diagrammatic rules for constructing f_xc in any order of the perturbation theory are given. The derived diagrammatic rules enable us to study the analytic properties of f_xc. It is shown that f_xc is finite at KS excitation energies in every order of the perturbation theory. If this were not the case, finite-order approximations to f_xc would not give sensible results for the kernel which plays the role of a mass operator connecting the KS density--density response function with the exact density--density response function. However, there is one important difference between f_xc and the self-energy, the latter playing the role of the mass operator connecting free and exact Green function. This becomes apparent when the perturbative expansion for f_xc is used to calculate the correction to the KS excitation energies. This correction is obtained up to the second-order with respect to the interaction and with respect to one- and two-particle irreducible elements. It is shown that if a perturbative approximation for f_xc is used, the correction to the KS excitation energies has to be calculated using a consistent approach based on the same order of perturbation theory for both quantities. One of the most promising applications of TDDFT in the linear response regime is the calculation of electronic excitation spectra of solids, where TDDFT could replace the computationally expensive solution of the BSE , which is normally used for describing the excitonic correlations. It is shown that the xc kernel can be separated into a quasiparticle and a excitonic part in an exact and unambiguous way. The diagrammatic rules derived for f_xc apply to these parts separately. This allows us to express the excitonic part in terms of the three-point function, for which an integral equation is derived. This integral equation is similar to the BSE and gives the exact translation of the BSE into the TDDFT language. There exists the possibility of cancellation effects in the kernel of the integral equation. If these effects are effective, a first-order approximation to the excitonic f_xc is sufficient. The cancellation effects in the equation for Lambda are explicitly studied for a two-band model semiconductor. Using a general short-range interaction between the quasiparticles (QPs) all relevant equations are solved analytically. It is shown that the cancellation is most effective directly at the band-gap energy. Consequently, the first-order and the exact excitonic f_xc are very similar close to the band-gap energy for all interaction strengths. For a weak interaction the first-order approximation to the excitonic f_xc describes the different excitonic correlations in the density©density response function very well. Qualitatively similar results are found when the QPs interact via the long-ranged Coulomb interaction.

•Jan 1, 2005•0 citations

Time-dependent density-functional theory beyond the adiabatic approximation: Insights from a two-electron model system

Carsten A. Ullrich

Most applications of time-dependent density-functional theory (TDDFT) use the adiabatic local-density approximation (ALDA) for the dynamical exchange-correlation potential V(xc)(r,t). An exact (i.e., nonadiabatic) extension of the ground-state LDA into the dynamical regime leads to a V(xc)(r,t) with a memory, which causes the electron dynamics to become dissipative. To illustrate and explain this nonadiabatic behavior, this paper studies the dynamics of two interacting electrons on a two-dimensional quantum strip of finite size, comparing TDDFT within and beyond the ALDA with numerical solutions of the two-electron time-dependent Schrodinger equation. It is shown explicitly how dissipation arises through multiple particle-hole excitations, and how the nonadiabatic extension of the ALDA fails for finite systems but becomes correct in the thermodynamic limit.

arXiv:cond-mat/0610341•The Journal of Chemical Physics•Dec 21, 2006•121 citations
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.
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 (xcxcxc) 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 xcxcxc 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 xcxcxc 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.
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.
The density-functional theory (DFT) is one of the most important ab initio methods for calculating properties of many-body systems. The success of DFT is based on the Kohn-Sham (KS) method and the approximations for the exchange-correlation (xc) effects. These approximations are well understood. However, DFT is first and foremost a theory for the ground state. It is very difficult to obtain properties of excited states. The time-dependent density-functional theory (TDDFT) can overcome this restriction and in the linear-response regime can give access to the properties of excited states. However, the xc effects, which are contained in the xc kernel f_xc , have to be approximated. The approximations for f_xc are not well understood up to now. Green-functions-based conventional many-body theory (MBT) is an alternative approach to the many-body problem. Especially when combined with diagrammatic techniques, MBT allows to derive physically motivated approximations. This work investigates the use of diagrammatic techniques from MBT for TDDFT . The aim of the investigations is to contribute to the development, testing, and understanding of the approximations involved in the time-dependent case. It is shown that it is possible to expand the xc kernel in a diagrammatic series. Explicit diagrammatic rules for constructing f_xc in any order of the perturbation theory are given. The derived diagrammatic rules enable us to study the analytic properties of f_xc. It is shown that f_xc is finite at KS excitation energies in every order of the perturbation theory. If this were not the case, finite-order approximations to f_xc would not give sensible results for the kernel which plays the role of a mass operator connecting the KS density--density response function with the exact density--density response function. However, there is one important difference between f_xc and the self-energy, the latter playing the role of the mass operator connecting free and exact Green function. This becomes apparent when the perturbative expansion for f_xc is used to calculate the correction to the KS excitation energies. This correction is obtained up to the second-order with respect to the interaction and with respect to one- and two-particle irreducible elements. It is shown that if a perturbative approximation for f_xc is used, the correction to the KS excitation energies has to be calculated using a consistent approach based on the same order of perturbation theory for both quantities. One of the most promising applications of TDDFT in the linear response regime is the calculation of electronic excitation spectra of solids, where TDDFT could replace the computationally expensive solution of the BSE , which is normally used for describing the excitonic correlations. It is shown that the xc kernel can be separated into a quasiparticle and a excitonic part in an exact and unambiguous way. The diagrammatic rules derived for f_xc apply to these parts separately. This allows us to express the excitonic part in terms of the three-point function, for which an integral equation is derived. This integral equation is similar to the BSE and gives the exact translation of the BSE into the TDDFT language. There exists the possibility of cancellation effects in the kernel of the integral equation. If these effects are effective, a first-order approximation to the excitonic f_xc is sufficient. The cancellation effects in the equation for Lambda are explicitly studied for a two-band model semiconductor. Using a general short-range interaction between the quasiparticles (QPs) all relevant equations are solved analytically. It is shown that the cancellation is most effective directly at the band-gap energy. Consequently, the first-order and the exact excitonic f_xc are very similar close to the band-gap energy for all interaction strengths. For a weak interaction the first-order approximation to the excitonic f_xc describes the different excitonic correlations in the density©density response function very well. Qualitatively similar results are found when the QPs interact via the long-ranged Coulomb interaction.
Most applications of time-dependent density-functional theory (TDDFT) use the adiabatic local-density approximation (ALDA) for the dynamical exchange-correlation potential V(xc)(r,t). An exact (i.e., nonadiabatic) extension of the ground-state LDA into the dynamical regime leads to a V(xc)(r,t) with a memory, which causes the electron dynamics to become dissipative. To illustrate and explain this nonadiabatic behavior, this paper studies the dynamics of two interacting electrons on a two-dimensional quantum strip of finite size, comparing TDDFT within and beyond the ALDA with numerical solutions of the two-electron time-dependent Schrodinger equation. It is shown explicitly how dissipation arises through multiple particle-hole excitations, and how the nonadiabatic extension of the ALDA fails for finite systems but becomes correct in the thermodynamic limit.