A density-matrix formalism is developed based on the one-particle density-matrix of a single-determinantal reference state. Unlike traditional density-functional-theory approaches, the
v-representable problem does not appear in the proposed method, nor the need to introduce functionals defined by a constrained search. The correlation-energy functionals are nonuniversal, in the sense that they depend on the external potential. Nevertheless, model systems can still be used to derive universal energy functionals. Variational and nonvariational energy functionals are introduced that yield the target-state energy when the reference state---or its corresponding one-particle density matrix---is constructed from Brueckner orbitals. Nonvariational energy functionals yield generalized Hartree-Fock equations involving a nonlocal correlation potential and the Hartree-Fock exchange operator; these equations are obtained by imposing the Brillouin-Brueckner condition. The same equations---for the most part---are obtained from variational energy functionals using functional minimization, yielding the (kernel of the) correlation potential as the functional derivative of correlation-energy functionals. Approximations for the correlation-energy functions are introduced, including a one-particle-density-matrix variant of the local-density approximation, a variant of the Colle-Salvetti functional, and a linear combination of the two that is a variant of the correlation-energy functional within the hybrid, three-parameter, Becke-Lee-Yang-Parr density functional (B3LYP).
Using an approach based on many body perturbation theory, the correlation energy $\cEco$ is expressed as an explicit functional of
ρ1,
v, and
vs, where
ρ1 is the one-particle density matrix from the noninteracting, or reference, determinantal-state;
v is the external potential from the interacting, or target, state;
vs is the (kernel of the) external potential from the noninteracting determinantal-state. In other words we have $\cEco[ρ_1,v,v_s]$. Anther possibility is the following explicit functional: $\cEco[ρ_1,v_{\text{co}},v_s]$, where
vco is the (kernel of the) correlation potential from the noninteracting Hamiltonian. The proposed method can, in principle, be used to compute $\cEco$ in a very accurate and efficient manner, since, like the Kohn--Sham approach, there are no virtual orbitals to consider. However, in contrast to the Kohn--Sham approach, $\cEco$ is a known, explicit functional that can be approximated in a systematic manner. For simplicity, we only consider noninteracting closed-shell states and target states that are nondegenerate, singlet ground-states; so, in that case,
ρ1 denotes the spin-less one-particle density matrix from the determinantal reference state.
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.