Response calculations based on an independent particle system with the exact one-particle density matrix: polarizabilities

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 [Phys. Rev. Lett. 105, 013002 (2010), 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 ($\omega\to0$), 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 H$_2$ and compare the performance of two different two-electron functionals: the phase-including L\"owdin-Shull functional and the density matrix form of the L\"owdin-Shull functional.