The density-density response function in time-dependent density functional theory: mathematical foundations and pole shifting
Thiago Carvalho Corso, Mi-Song Dupuy, Gero Friesecke
TL;DR
This work provides a rigorous mathematical foundation for TDDFT in the RPA by proving existence and uniqueness of the Dyson equation solution for the density-density response function, and by establishing that the RPA poles are forward-shifted relative to the non-interacting poles. It introduces a meromorphic extension framework for the Fourier transform of the response, connects pole locations to an eigenvalue problem involving the Hartree operator, and delivers a rigorous Lehmann representation that remains valid in the presence of continuous spectra. The analysis shows that RPA poles can be computed via a reduced eigenproblem and quantifies pole ranks, offering a solid theoretical justification for standard Casida-type approaches used in quantum chemistry to determine excitation frequencies. Overall, the results explain why Kohn–Sham gaps typically underestimate true excitation energies and provide concrete, operator-theoretic criteria for locating and counting RPA poles within the physically relevant spectral window $(-\Omega,\Omega)$.
Abstract
We establish existence and uniqueness of the solution to the Dyson equation for the density-density response function in time-dependent density functional theory (TDDFT) in the random phase approximation (RPA). We show that the poles of the RPA density-density response function are forward-shifted with respect to those of the non-interacting response function, thereby explaining mathematically the well known empirical fact that the non-interacting poles (given by the spectral gaps of the time-independent Kohn-Sham equations) underestimate the true transition frequencies. Moreover we show that the RPA poles are solutions to an eigenvalue problem, justifying the approach commonly used in the physics community to compute these poles.
