6.3 Møller-Plesset Perturbation Theory

Møller-Plesset Perturbation Theory⁶⁶⁵ 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 ${\mathrm{\Psi }}_0$ and energy $E_0$ 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 $H_0$ eigenvalue problem. If we assume that the Hartree-Fock wave function ${\mathrm{\Psi }}_0$ and energy $E_0$ lie near the exact wave function $\mathrm{\Psi }$ and energy $E$, we can now write the exact Hamiltonian operator as

$$H=H_0+\lambda V$$

(6.1)

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

$$E=E^{(0)}+\lambda E^{(1)}+{\lambda }^2{E}^{(2)}+{\lambda }^3{E}^{(3)}+\mathrm{\dots }$$

(6.2)

and

$$\mathrm{\Psi }={\mathrm{\Psi }}_0+\lambda {\mathrm{\Psi }}^{(1)}+{\lambda }^2{\mathrm{\Psi }}^{(2)}+{\lambda }^3{\mathrm{\Psi }}^{(3)}+\mathrm{\dots }$$

(6.3)

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

$$H_0{\mathrm{\Psi }}_0=E^{(0)}{\mathrm{\Psi }}_0$$

(6.4)

$$H_0{\mathrm{\Psi }}^{(1)}+V{\mathrm{\Psi }}_0=E^{(0)}{\mathrm{\Psi }}^{(1)}+E^{(1)}{\mathrm{\Psi }}_0$$

(6.5)

$$H_0{\mathrm{\Psi }}^{(2)}+V{\mathrm{\Psi }}^{(1)}=E^{(0)}{\mathrm{\Psi }}^{(2)}+E^{(1)}{\mathrm{\Psi }}^{(1)}+E^{(2)}{\mathrm{\Psi }}_0$$

(6.6)

and so forth. Multiplying each of the above equations by ${\mathrm{\Psi }}_0$ and integrating over all space yields the following expression for the $n$th-order (MP$n$) energy:

$$E^{(0)}=⟨{\mathrm{\Psi }}_0|H_0|{\mathrm{\Psi }}_0⟩$$

(6.7)

$$E^{(1)}=⟨{\mathrm{\Psi }}_0|V|{\mathrm{\Psi }}_0⟩$$

(6.8)

$$E^{(2)}=⟨{\mathrm{\Psi }}_0|V|{\mathrm{\Psi }}^{(1)}⟩$$

(6.9)

Thus, the Hartree-Fock energy

$$E_0=⟨{\mathrm{\Psi }}_0\left|H_0+V\right|{\mathrm{\Psi }}_0⟩$$

(6.10)

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

$$E_0=E^{(0)}+E^{(1)}$$

(6.11)

The correlation energy can then be written

$$E_{\mathrm{corr}}=E_0^{(2)}+E_0^{(3)}+E_0^{(4)}+\mathrm{\dots }$$

(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

$$E_0^{(2)}=-\frac{1}{4}\sum _{ab}^{\mathrm{virt}}\sum _{ij}^{\mathrm{occ}}\frac{{\left|⟨ab||ij⟩\right|}^2}{{\epsilon }_a+{\epsilon }_b-{\epsilon }_i-{\epsilon }_j}$$

(6.13)

where

$$⟨ab||ij⟩=⟨ab|ij⟩-⟨ab|ji⟩$$

(6.14)

and

$$⟨ab|cd⟩=\int \int {\psi }_a({𝐫}_1){\psi }_c({𝐫}_1)\frac{1}{r_{12}}{\psi }_b({𝐫}_2){\psi }_d({𝐫}_2)𝑑{𝐫}_1𝑑{𝐫}_2$$

(6.15)

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

$$⟨ab|cd⟩=\sum _{\mu }\sum _{\nu }\sum _{\lambda }\sum _{\sigma }{C}_{\mu a}{C}_{\nu c}{C}_{\lambda b}{C}_{\sigma d}\left(\mu \nu |\lambda \sigma \right)$$

(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.