Regularized and Opposite Spin-Scaled Functionals from Møller-Plesset Adiabatic Connection─Higher Accuracy at Lower Cost

Chemical Potentials and the One-Electron Hamiltonian of the Second-Order Perturbation Theory from the Functional Derivative Approach

We develop a functional derivative approach to calculate the chemical potentials of second-order perturbation theory (MP2). In the functional derivative approach, the correlation part of the MP2 chemical potential, which is the derivative of the MP2 correlation energy with respect to the occupation number of frontier orbitals, is obtained from the chain rule via the noninteracting Green's function. First, the MP2 correlation energy is expressed in terms of the noninteracting Green's function, and its functional derivative to the noninteracting Green's function is the second-order self-energy. Then, the derivative of the noninteracting Green's function to the occupation number is obtained by including the orbital relaxation effect. We show that the MP2 chemical potentials obtained from the functional derivative approach agree with that obtained from the finite difference approach. The one-electron Hamiltonian, defined as the derivative of the MP2 energy with respect to the one particle density matrix, is also derived using the functional derivative approach, which can be used in the self-consistent calculations of MP2 and double-hybrid density functionals. The developed functional derivative approach is promising for calculating the chemical potentials and the one-electron Hamiltonian of approximate functionals and many-body perturbation approaches dependent explicitly on the noninteracting Green's function.

Møller-Plesset Adiabatic Connection at Large Coupling Strengths for Open-Shell Systems

We study the adiabatic connection that has as weak-coupling expansion the Møller-Plesset perturbation series, generalizing to the open-shell case previous closed-shell results for the large-coupling limit. We first focus on the hydrogen atom with fractional spins, providing results along the adiabatic connection from small to large coupling strengths. We reveal an intriguing phase diagram and an equation for the large-coupling leading order that has closed-form solutions for specific choices of its relevant quantum numbers. We then show that the hydrogen atom results provide variational estimates for the large-coupling leading terms for the general many-electron open-shell case in terms of functionals of the Hartree-Fock α-spin and β-spin densities.

Reaching High Accuracy for Energetic Properties at Second-Order Perturbation Cost by Merging Self-Consistency and Spin-Opposite Scaling

Quantum chemical methods dealing with challenging systems while retaining low computational costs have attracted attention. In particular, many efforts have been devoted to developing new methods based on second-order perturbation that may be the simplest correlated method beyond Hartree-Fock. We have recently developed a self-consistent perturbation theory named one-body Møller-Plesset second-order perturbation theory (OBMP2) and shown that it can resolve issues caused by the noniterative nature of standard perturbation theory. In this work, we extend the method by introducing spin-opposite scaling to the double-excitation amplitudes, resulting in the O2BMP2 method. We assess the O2BMP2 performance on the triple-bond N2 dissociation, singlet-triplet gaps, and ionization potentials. O2BMP2 performs much better than standard MP2 and reaches the accuracy of coupled-cluster methods in all cases considered in this work.

Adiabatic connection interaction strength interpolation method made accurate for the uniform electron gas

The adiabatic connection interaction strength interpolation (ISI)-like method provides a high-level expression for the correlation energy, being, in principle, exact not only in the weak-interaction limit, where it recovers the second-order Görling-Levy perturbation term, but also in the strong-interaction limit that is described by the strictly correlated electron approach. In this work, we construct a genISI functional made accurate for the uniform electron gas, a solid-state physics paradigm that is a very difficult test for ISI-like correlation functionals. We assess the genISI functional for various jellium spheres with the number of electrons Z ≤ 912 and for the non-relativistic noble atoms with Z ≤ 290. For the jellium clusters, the genISI is remarkably accurate, while for the noble atoms, it shows a good performance, similar to other ISI-like methods. Then, the genISI functional can open the path using the ISI-like method in solid-state calculations.