Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Correlation Energy Estimators based on M{\o}ller-Plesset Perturbation Theory

Herbert H. H. Homeier

Published Mar 5, 1996• 0 citations
OpenAlex

Abstract

Some methods for the convergence acceleration of the Moller-Plesset perturbation series for the correlation energy are discussed. The order-by-order summation is less effective than the Feenberg series. The latter is obtained by renormalizing the unperturbed Hamilton operator by a constant factor that is optimized for the third order energy. In the fifth order case, the Feenberg series can be improved by order-dependent optimization of the parameter. Alternatively, one may use Pad{e} approximants or a further method based on effective characteristic polynomials to accelerate the convergence of the perturbation series. Numerical evidence is presented that, besides the Feenberg-type approaches, suitable Pad{e} approximants, and also the effective second order characteristic polynomial, are excellent tools for correlation energy estimation.

Related papers

On the acceleration of convergence of many-body perturbation theory. I. General theory

Klaus Dietz, C Schmidt, Markus Warken +1 more

The theoretical background of schemes for accelerating convergence of Moller-Plesset-Goldstone perturbation theory is discussed. Aiming at practical applications the authors consider approximate formulations of the many-body theory of atoms and molecules in finite-dimensional Hilbert spaces. A key role is played by nearby singularities of energy levels and eigenfunctions considered as functions of a suitably chosen coupling parameter. Infinite subsummations in the perturbation series are performed such that the influence of these singularities is eliminated or, at least, sufficiently attenuated. The danger of double counting of diagrammatic contributions already included in these infinite subsummations is avoided by the very construction of the authors' scheme.

Journal of Physics B Atomic Molecular and Optical Physics•Jul 14, 1993•35 citations

Stochastic many-body perturbation theory for electron correlation energies

Zhendong Li

Treating electron correlation more accurately and efficiently is at the heart of the development of electronic structure methods. In the present work, we explore the use of stochastic approaches to evaluate high-order electron correlation energies, whose conventional computational scaling is unpleasantly steep, being O(Nn+3) with respect to the system size N and the perturbation order n for the Møller-Plesset (MP) series. To this end, starting from Goldstone’s time-dependent formulation of ab initio many-body perturbation theory (MBPT), we present a reformulation of MBPT, which naturally leads to a Monte Carlo scheme with O(nN2 + n2N + f(n)) scaling at each step, where f(n) is a function of n depending on the specific numerical scheme. Proof-of-concept calculations demonstrate that the proposed quantum Monte Carlo algorithm successfully extends the previous Monte Carlo approaches for MP2 and MP3 to higher orders by overcoming the factorial scaling problem. For the first time, Goldstone’s time-dependent formulation is made useful numerically for electron correlation energies, not only being purely as a theoretical tool.

arXiv:1909.00307•The Journal of Chemical Physics•Dec 27, 2019•12 citations

A GPU Implementation of the Second- and Third-Order Møller-Plesset Perturbation Theory

Satoki Tsuji, Yasuaki Ito, Nobuya Yokogawa +4 more

Post-Hartree-Fock methods are wavefunction-based quantum chemistry approaches that explicitly account for the electron correlation beyond the mean-field approximation, allowing accurate determination of electron orbitals and energies in a given system. However, their computational cost scales as the fourth power or higher with respect to the number of atoms, severely limiting their practical applications. This paper proposes an efficient GPU implementation of the Møller-Plesset method, a post-Hartree-Fock approach based on perturbation theory. We focus on the computation of second- and thirdorder perturbation energies (MP2 and MP3) for closed-shell systems, which incorporate electron correlation on top of the Hartree-Fock energy. To accelerate the extensive contraction operations arising in the Møller-Plesset method, our implementation employs a formulation that maximizes the use of highly optimized matrix multiplication routines. Furthermore, we effectively utilize computational capabilities of GPU architecture through a hierarchical reduction using warp-level primitives and shared memory, as well as contiguous coalesced access to matrix elements. Performance evaluation on an NVIDIA A100 GPU shows up to 255.5 and 71.1 times speedups for the MP2 correlation energy computation, compared with the widely used quantum chemistry software package PySCF and its GPU plugin GPU4PySCF, respectively. In the MP3 calculations, our DGEMM-based implementation achieves a maximum speedup of 96.2 times over PySCF.

International Symposium on Computing and Networking Workshops (CANDARW)•Nov 25, 2025•1 citation

Applications of Padé approximants to correlation functions

R. Paul

In this paper we have presented a general scheme whereby the method of Padé approximants may be applied to the power series expansion of time dependent correlation functions. It has been demonstrated by actual numerical calculations, using the computer simulated results of Berne and Harp [5], that even a few terms in the power series expansion are capable of reproducing the correlation curves to a fairly large extent. It is pointed out that the method can be extended to higher orders without the introduction of fitting parameters, and without violating the time parity requirements of the Liouville equation. In the final analysis the method is a truly microscopic method since it relies on the microscopic hamiltonian, and the underlying principles are very straightforward. The basic requirements are the spectral moments, but we do not need many of these as the Padé approximants are highly convergent.

Molecular Physics•Sep 1, 1980•11 citations

Regularized second-order correlation methods for extended systems

Elisabeth Keller, Theodoros Tsatsoulis, Karsten Reuter +1 more

Second-order Møller-Plesset perturbation theory (MP2) constitutes the simplest form of many-body wavefunction theory and often provides a good compromise between efficiency and accuracy. There are, however, well-known limitations to this approach. In particular, MP2 is known to fail or diverge for some prototypical condensed matter systems like the homogeneous electron gas (HEG) and to overestimate dispersion-driven interactions in strongly polarizable systems. In this paper, we explore how the issues of MP2 for metallic, polarizable, and strongly correlated periodic systems can be ameliorated through regularization. To this end, two regularized second-order methods (including a new, size-extensive Brillouin-Wigner approach) are applied to the HEG, the one-dimensional Hubbard model, and the graphene-water interaction. We find that regularization consistently leads to improvements over the MP2 baseline and that different regularizers are appropriate for the various systems.

The Journal of Chemical Physics•Dec 22, 2021•17 citations

A Regularized Second-Order Correlation Method from Green’s Function Theory

Christopher J. N. Coveney, David P. Tew

We present a scalable single-particle framework to treat electronic correlation in molecules and materials motivated by Green's function theory. We derive a size-extensive Brillouin-Wigner perturbation theory from the single-particle Green's function by introducing the Goldstone self-energy. This new ground state correlation energy, referred to as Quasi-Particle MP2 theory (QPMP2), avoids the characteristic divergences present in both second-order Møller-Plesset perturbation theory and Coupled Cluster Singles and Doubles within the strongly correlated regime. We show that the exact ground state energy and properties of the Hubbard dimer are reproduced by QPMP2 and demonstrate the advantages of the approach for larger Hubbard models where the metal-to-insulator transition is qualitatively reproduced, contrasting with the complete failure of traditional methods. We apply this formalism to characteristic strongly correlated molecular systems and show that QPMP2 provides an efficient, size-consistent regularization of MP2.

arXiv:2302.13296•Journal of Chemical Theory and Computation•Jun 27, 2023•22 citations
Some methods for the convergence acceleration of the Moller-Plesset perturbation series for the correlation energy are discussed. The order-by-order summation is less effective than the Feenberg series. The latter is obtained by renormalizing the unperturbed Hamilton operator by a constant factor that is optimized for the third order energy. In the fifth order case, the Feenberg series can be improved by order-dependent optimization of the parameter. Alternatively, one may use Pad{e} approximants or a further method based on effective characteristic polynomials to accelerate the convergence of the perturbation series. Numerical evidence is presented that, besides the Feenberg-type approaches, suitable Pad{e} approximants, and also the effective second order characteristic polynomial, are excellent tools for correlation energy estimation.
The theoretical background of schemes for accelerating convergence of Moller-Plesset-Goldstone perturbation theory is discussed. Aiming at practical applications the authors consider approximate formulations of the many-body theory of atoms and molecules in finite-dimensional Hilbert spaces. A key role is played by nearby singularities of energy levels and eigenfunctions considered as functions of a suitably chosen coupling parameter. Infinite subsummations in the perturbation series are performed such that the influence of these singularities is eliminated or, at least, sufficiently attenuated. The danger of double counting of diagrammatic contributions already included in these infinite subsummations is avoided by the very construction of the authors' scheme.
Treating electron correlation more accurately and efficiently is at the heart of the development of electronic structure methods. In the present work, we explore the use of stochastic approaches to evaluate high-order electron correlation energies, whose conventional computational scaling is unpleasantly steep, being O(Nn+3) with respect to the system size N and the perturbation order n for the Møller-Plesset (MP) series. To this end, starting from Goldstone’s time-dependent formulation of ab initio many-body perturbation theory (MBPT), we present a reformulation of MBPT, which naturally leads to a Monte Carlo scheme with O(nN2 + n2N + f(n)) scaling at each step, where f(n) is a function of n depending on the specific numerical scheme. Proof-of-concept calculations demonstrate that the proposed quantum Monte Carlo algorithm successfully extends the previous Monte Carlo approaches for MP2 and MP3 to higher orders by overcoming the factorial scaling problem. For the first time, Goldstone’s time-dependent formulation is made useful numerically for electron correlation energies, not only being purely as a theoretical tool.
Post-Hartree-Fock methods are wavefunction-based quantum chemistry approaches that explicitly account for the electron correlation beyond the mean-field approximation, allowing accurate determination of electron orbitals and energies in a given system. However, their computational cost scales as the fourth power or higher with respect to the number of atoms, severely limiting their practical applications. This paper proposes an efficient GPU implementation of the Møller-Plesset method, a post-Hartree-Fock approach based on perturbation theory. We focus on the computation of second- and thirdorder perturbation energies (MP2 and MP3) for closed-shell systems, which incorporate electron correlation on top of the Hartree-Fock energy. To accelerate the extensive contraction operations arising in the Møller-Plesset method, our implementation employs a formulation that maximizes the use of highly optimized matrix multiplication routines. Furthermore, we effectively utilize computational capabilities of GPU architecture through a hierarchical reduction using warp-level primitives and shared memory, as well as contiguous coalesced access to matrix elements. Performance evaluation on an NVIDIA A100 GPU shows up to 255.5 and 71.1 times speedups for the MP2 correlation energy computation, compared with the widely used quantum chemistry software package PySCF and its GPU plugin GPU4PySCF, respectively. In the MP3 calculations, our DGEMM-based implementation achieves a maximum speedup of 96.2 times over PySCF.
In this paper we have presented a general scheme whereby the method of Padé approximants may be applied to the power series expansion of time dependent correlation functions. It has been demonstrated by actual numerical calculations, using the computer simulated results of Berne and Harp [5], that even a few terms in the power series expansion are capable of reproducing the correlation curves to a fairly large extent. It is pointed out that the method can be extended to higher orders without the introduction of fitting parameters, and without violating the time parity requirements of the Liouville equation. In the final analysis the method is a truly microscopic method since it relies on the microscopic hamiltonian, and the underlying principles are very straightforward. The basic requirements are the spectral moments, but we do not need many of these as the Padé approximants are highly convergent.
Second-order Møller-Plesset perturbation theory (MP2) constitutes the simplest form of many-body wavefunction theory and often provides a good compromise between efficiency and accuracy. There are, however, well-known limitations to this approach. In particular, MP2 is known to fail or diverge for some prototypical condensed matter systems like the homogeneous electron gas (HEG) and to overestimate dispersion-driven interactions in strongly polarizable systems. In this paper, we explore how the issues of MP2 for metallic, polarizable, and strongly correlated periodic systems can be ameliorated through regularization. To this end, two regularized second-order methods (including a new, size-extensive Brillouin-Wigner approach) are applied to the HEG, the one-dimensional Hubbard model, and the graphene-water interaction. We find that regularization consistently leads to improvements over the MP2 baseline and that different regularizers are appropriate for the various systems.
We present a scalable single-particle framework to treat electronic correlation in molecules and materials motivated by Green's function theory. We derive a size-extensive Brillouin-Wigner perturbation theory from the single-particle Green's function by introducing the Goldstone self-energy. This new ground state correlation energy, referred to as Quasi-Particle MP2 theory (QPMP2), avoids the characteristic divergences present in both second-order Møller-Plesset perturbation theory and Coupled Cluster Singles and Doubles within the strongly correlated regime. We show that the exact ground state energy and properties of the Hubbard dimer are reproduced by QPMP2 and demonstrate the advantages of the approach for larger Hubbard models where the metal-to-insulator transition is qualitatively reproduced, contrasting with the complete failure of traditional methods. We apply this formalism to characteristic strongly correlated molecular systems and show that QPMP2 provides an efficient, size-consistent regularization of MP2.