We point out an error in the argument [PRL 105, 013002 (2010)] that the time independence of the occupation numbers in the adiabatic approximation follows from the invariance of the ground-state interaction energy functional with respect to changes in the phases of the natural orbitals.
The adiabatic approximation is problematic in time-dependent density matrix functional theory. With pure density matrix functionals (invariant under phase change of the natural orbitals) it leads to lack of response in the occupation numbers, hence wrong frequency dependent responses, in particular α(ω→0)≠α(0) (the static polarizability). We propose to relinquish the requirement that the functional must be a pure one-body reduced density matrix (1RDM) functional, and to introduce additional variables which can be interpreted as phases of the one-particle states of the independent particle reference system formed with the natural orbitals, thus obtaining so-called phase-including natural orbital (PINO) functionals. We also stress the importance of the correct choice of the complex conjugation in the two-electron integrals in the commonly used functionals (they should not be of exchange type). We demonstrate with the Löwdin-Shull energy expression for two-electron systems, which is an example of a PINO functional, that for two-electron systems exact responses (polarizabilities, excitation energies) are obtained, while writing this energy expression in the usual way as a 1RDM functional yields erroneous responses.
Recently, we have demonstrated that the problems finding a suitable adiabatic approximation in time-dependent one-body reduced density matrix functional theory can be remedied by introducing an additional degree of freedom to describe the system: the phase of the natural orbitals [K. J. H. Giesbertz, O. V. Gritsenko, and E. J. Baerends, Phys. Rev. Lett. 105, 013002 (2010); K. J. H. Giesbertz, O. V. Gritsenko, and E. J. Baerends, J. Chem. Phys. 133, 174119 (2010)]. In this article we will show in detail how the frequency-dependent response equations give the proper static limit (ω → 0), including the perturbation in the chemical potential, which is required in static response theory to ensure the correct number of particles. Additionally we show results for the polarizability for H2 and compare the performance of two different two-electron functionals: the phase-including Löwdin-Shull functional and the density matrix form of the Löwdin-Shull functional.
We use the natural orbitals to define an independent particle system, from which the exact one-particle density matrix can be obtained with an ensemble of degenerate determinantal ground states. Also defining explicit phases for the orbitals, and admitting functionals that are dependent on those phases, time-dependent equations for the orbitals and occupation numbers are obtained from an action principle. The wrong polarizability and lack of double excitations of straightforward adiabatic time-dependent density matrix functional theory are then corrected, and the important symmetry χ(ω)=χ{*}(-ω), lost in previous ad hoc improvements, is restored. The extension of the response calculations beyond the occupied-virtual pairs, which are the only ones admitted in time-dependent density functional theory, leads to greatly improved response properties.