6 Wave Function-Based Correlation Methods

6.3 Møller-Plesset Perturbation Theory

Møller-Plesset Perturbation Theory665 is a widely used method for approximating the correlation energy of molecules. In particular, second-order Møller-Plesset perturbation theory (MP2) is one of the simplest and most useful levels of theory beyond the Hartree-Fock approximation. Conventional and local MP2 methods available in Q-Chem are discussed in detail in Sections 6.4 and 6.5 respectively. The MP3 method is still occasionally used, while MP4 calculations are quite commonly employed as part of the G2 and G3 thermochemical methods.198, 197 In the remainder of this section, the theoretical basis of Møller-Plesset theory is reviewed.

The Hartree-Fock wave function Ψ0 and energy E0 are approximate solutions (eigenfunction and eigenvalue) to the exact Hamiltonian eigenvalue problem or Schrödinger’s electronic wave equation, Eq. (4.5). The HF wave function and energy are, however, exact solutions for the Hartree-Fock Hamiltonian H0 eigenvalue problem. If we assume that the Hartree-Fock wave function Ψ0 and energy E0 lie near the exact wave function Ψ and energy E, we can now write the exact Hamiltonian operator as

H=H0+λV (6.1)

where V is the small perturbation and λ is a dimensionless parameter. Expanding the exact wave function and energy in terms of the HF wave function and energy yields

E=E(0)+λE(1)+λ2E(2)+λ3E(3)+ (6.2)

and

Ψ=Ψ0+λΨ(1)+λ2Ψ(2)+λ3Ψ(3)+ (6.3)

Substituting these expansions into the Schrödinger equation and collecting terms according to powers of λ yields

H0Ψ0=E(0)Ψ0 (6.4)
H0Ψ(1)+VΨ0=E(0)Ψ(1)+E(1)Ψ0 (6.5)
H0Ψ(2)+VΨ(1)=E(0)Ψ(2)+E(1)Ψ(1)+E(2)Ψ0 (6.6)

and so forth. Multiplying each of the above equations by Ψ0 and integrating over all space yields the following expression for the nth-order (MPn) energy:

E(0)=Ψ0|H0|Ψ0 (6.7)
E(1)=Ψ0|V|Ψ0 (6.8)
E(2)=Ψ0|V|Ψ(1) (6.9)

Thus, the Hartree-Fock energy

E0=Ψ0|H0+V|Ψ0 (6.10)

is simply the sum of the zeroth- and first- order energies

E0=E(0)+E(1) (6.11)

The correlation energy can then be written

Ecorr=E0(2)+E0(3)+E0(4)+ (6.12)

of which the first term is the MP2 energy.

It can be shown that the MP2 energy can be written (in terms of spin-orbitals) as

E0(2)=-14abvirtijocc|ab||ij|2εa+εb-εi-εj (6.13)

where

ab||ij=ab|ij-ab|ji (6.14)

and

ab|cd=ψa(𝐫1)ψc(𝐫1)1r12ψb(𝐫2)ψd(𝐫2)𝑑𝐫1𝑑𝐫2 (6.15)

which can be written in terms of the two-electron repulsion integrals

ab|cd=μνλσCμaCνcCλbCσd(μν|λσ) (6.16)

Expressions for higher order terms follow similarly, although with much greater algebraic and computational complexity. MP3 and particularly MP4 (the third and fourth order contributions to the correlation energy) are both occasionally used, although they are increasingly supplanted by the coupled-cluster methods described in the following sections. The disk and memory requirements for MP3 are similar to the self-consistent pair correlation methods discussed in Section 6.8 while the computational cost of MP4 is similar to the (T) corrections discussed in Section 6.9.