Linear interpolation method in ensemble Kohn-Sham and range-separated density-functional approximations for excited states

TL;DR

This paper introduces a linear interpolation method (LIM) within ensemble density-functional theory to improve the calculation of electronic excitation energies, especially for challenging cases like charge transfer and double excitations.

Contribution

The paper proposes a novel linear interpolation approach (LIM) for ensemble DFT that effectively incorporates weight dependence and improves excitation energy predictions.

Findings

LIM outperforms standard TDDFT for various excitations.

Promising results for systems like He, Be, H₂, and HeH⁺.

Effective for both single and double excitations.

Abstract

Gross-Oliveira-Kohn density functional theory (GOK-DFT) for ensembles is in principle very attractive, but has been hard to use in practice. A novel, practical model based on GOK-DFT for the calculation of electronic excitation energies is discussed. The new model relies on two modifications of GOK-DFT: use of range separation and use of the slope of the linearly-interpolated ensemble energy, rather than orbital energies. The range-separated approach is appealing as it enables the rigorous formulation of a multi-determinant state-averaged DFT method. In the exact theory, the short-range density functional, that complements the long-range wavefunction-based ensemble energy contribution, should vary with the ensemble weights even when the density is held fixed. This weight dependence ensures that the range-separated ensemble energy varies linearly with the ensemble weights. When the (weight-independent) ground-state short-range exchange-correlation functional is used in this context, curvature appears thus leading to an approximate weight-dependent excitation energy. In order to obtain unambiguous approximate excitation energies, we propose to interpolate linearly the ensemble energy between equiensembles. It is shown that such a linear interpolation method (LIM) can be rationalized and that it effectively introduces weight dependence effects. As proof of principle, LIM has been applied to He, Be, H$_2$ in both equilibrium and stretched geometries as well as the stretched HeH$^+$ molecule. Very promising results have been obtained for both single (including charge transfer) and double excitations with spin-independent short-range local and semi-local functionals. Even at the Kohn--Sham ensemble DFT level, that is recovered when the range-separation parameter is set to zero, LIM performs better than standard time-dependent DFT.