Perturbative analysis of the functional U[ n , ψ] that yields the correlation component U of the electron–electron repulsion energy in terms of the vectors ψ(1) and n of the natural spinorbitals and their occupation numbers (the 1-matrix functional) facilitates examination of the flaws inherent to the present implementations of the density matrix functional theory. Recognizing that the practical usefulness of any approximate 1-matrix functional hinges upon its capability of exactly reproducing the leading contribution to U at the limit of vanishing electron–electron interactions gives rise to asymptotic bilinear constraints for the (exact or model) 2-cumulant that enters the expression for U. The asymptotic behavior of certain blocks of is found to be equally important. These identities, which are obtained for both the single-determinantal and a model multideterminantal cases, take precedence over the linear constraints commonly enforced in the course of approximate construction of such functionals. This observation reveals the futility of designing sophisticated approximations tailored for the second-order contribution to while neglecting proper formulation of the respective first-order contribution that in the case of the so-called JKL-only functionals requires abandoning the JK-dependence altogether. It has its repercussions not only for the functionals of the PNOF family but also for the expressions involving only the L-type two-electron repulsion integrals (in the guise of their exchange counterparts) that account only for the correlation effects due to electrons with antiparallel spins and are well-defined only for spin-unpolarized and high-spin systems (yielding vanishing U for the latter).

If citing, it is advised that you check and use the publisher's definitive version for pagination, volume/issue, and date of publication details. And where the final published version is provided on the Research Portal, if citing you are again advised to check the publisher's website for any subsequent corrections.

Convergence features of the Rayleigh-Schrödinger perturbation theory (PT) strongly depend on the partitioning applied. We investigate the large order behavior of the Møller-Plesset and Epstein Nesbet partitionings in comparison with a less known partitioning obtained by level shift parameters minimizing the norm of operator Q^W^, with W^ being the perturbation operator while Q standing for the reduced resolvent of the zero order Hamiltonian H^. Numerical results, presented for molecular systems for the first time, indicate that it is possible to find level shift parameters in this way which convert divergent perturbation expansions to convergent ones in some cases. Besides numerical calculations of high-order PT terms, convergence radii of the corresponding perturbation expansions are also estimated using quadratic Padé approximants.

No abstract

Allowing triplet components of individual geminals, spin-contaminated strongly orthogonal geminal wave functions may emerge, which can be ameliorated by spin-projection techniques. Of the latter, halfprojection was previously shown to be useful, offering a compromise between the amount of remaining spin-contamination and the violation of size-consistency generated by projection. This paper investigates how a half-projected spin-contaminated geminal wave function can be improved by multi-configuration perturbation theory, to incorporate dynamical correlation effects.