In time-dependent density-functional theory, an important quantity is the exchange-correlation kernel whose knowledge allows for the calculations of the excitation spectrum of electronic systems. In most applications, the exchange-correlation kernel is approximated by the adiabatic local-density approximation, ignoring the frequency dependence and nonlocality. To gain insight into the nature of the exact exchange-correlation kernel, we have considered a two-dimensional one-band Hubbard model. The calculated exact exchange-correlation kernel reveals a number of striking features. It has a weak energy dependence up to the main excitation energy. On the other hand, the exact kernel shows a strong energy dependence in energy regions, where there are many-body excitations that are not contained in the local-density approximation, even though these excitations have small weight. It is found that a static approximation can well reproduce the main excitation peaks, but the satellite structures arising from many-body interaction are unlikely to be accounted for with a static approximation.
Time-dependent density-functional theory is a rather accurate and efficient way to compute electronic excitations for finite systems. However, in the macroscopic limit (systems of increasing size), for the usual adiabatic random-phase, local-density, or generalized-gradient approximations, one recovers the Kohn-Sham independent-particle picture, and thus the incorrect band gap. To clarify this trend, we investigate the macroscopic limit of the exchange-correlation kernel in such approximations by means of an algebraical analysis complemented with numerical studies of a one-dimensional tight-binding model. We link the failure to shift the Kohn-Sham spectrum of these approximate kernels to the fact that the corresponding operators in the transition space act only on a finite subspace.
The exact time-dependent correlation energy in time-dependent density-functional theory is shown to sometimes become positive, which is impossible with most present TDDFT approximations. Both the correlation potential and energy can be strongly nonlocal in time. A new inequality is derived for the time-dependent exchange-correlation energy. The correlation energy appears to scale to a constant function of scaled time in the high-density limit. In the linear response regime, the correlation energy is shown to become purely adiabatic, but the correlation potential is generally nonadiabatic. The usefulness of the virial theorem as a test of numerical accuracy is demonstrated. All results are found or inspired by exact numerical solution of a simple model system (Hooke’s atom), and inversion of the corresponding Kohn–Sham equations.
The exact exchange-correlation (xc) kernel ${f}_{\text{xc}}(x,{x}^{\ensuremath{'}},\ensuremath{\omega})$ of linear response time-dependent density functional theory is computed over a wide range of frequencies for three canonical one-dimensional finite systems. Methods used to ensure the numerical robustness of
fxc are set out. The frequency dependence of
fxc is found to be largely due to its analytic structure, i.e., its singularities at certain frequencies, which are required in order to capture particular transitions, including those of double excitation character. However, within the frequency range of the first few interacting excitations,
fxc is approximately $\ensuremath{\omega}$ independent, meaning the exact adiabatic approximation ${f}_{\text{xc}}(\ensuremath{\omega}=0)$ remedies the failings of the local density approximation and random phase approximation for these lowest transitions. The key differences between the exact
fxc and its common approximations are analyzed, and cannot be eliminated by exploiting the limited gauge freedom in
fxc. The optical spectrum benefits from using as accurate as possible an
fxc and ground-state xc potential, while maintaining exact compatibility between the two is of less importance.
In the context of inhomogeneous one-dimensional finite systems, recent numerical advances [Phys. Rev. B 103, 125155 (2021)] allow us to compute the exact coupling-constant dependent exchange-correlation kernel ${f}_{\text{xc}}^{\ensuremath{\lambda}}(x,{x}^{\ensuremath{'}},\ensuremath{\omega})$ within linear response time-dependent density-functional theory. This permits an improved understanding of ground-state total energies derived from the adiabatic-connection fluctuation-dissipation theorem (ACFDT). We consider both one-shot and self-consistent ACFDT calculations, and demonstrate that chemical accuracy is reliably preserved when the frequency dependence in the exact functional ${f}_{\text{xc}}[n](\ensuremath{\omega}=0)$ is neglected. This performance is understood on the grounds that the exact
fxc[n] varies slowly over the most relevant $\ensuremath{\omega}$ range (but not in general), and hence the spatial structure in ${f}_{\text{xc}}[n](\ensuremath{\omega}=0)$ is able to largely remedy the principal issue in the present context: self-interaction (examined from the perspective of the exchange-correlation hole). Moreover, we find that the implicit orbitals contained within a self-consistent ACFDT calculation utilizing the adiabatic exact kernel ${f}_{\text{xc}}[n](\ensuremath{\omega}=0)$ are remarkably similar to the exact Kohn-Sham orbitals, thus further establishing that the majority of the physics required to capture the ground-state total energy resides in the spatial dependence of
fxc[n] at $\ensuremath{\omega}=0$.
We investigate a number of formal properties of the adiabatic strictly-correlated electrons (SCE) functional, relevant for time-dependent potentials and for kernels in linear response time-dependent density functional theory. Among the former, we focus on the compliance to constraints of exact many-body theories, such as the generalised translational invariance and the zero-force theorem. Within the latter, we derive an analytical expression for the adiabatic SCE Hartree exchange-correlation kernel in one dimensional systems, and we compute it numerically for a variety of model densities. We analyse the non-local features of this kernel, particularly the ones that are relevant in tackling problems where kernels derived from local or semi-local functionals are known to fail.