Mpller-Plesset expansion

Arrays in which the calculated energy components will be stored are set to zero, etotal will contain the total energy whilst etwo will contain the two-body component. iz=l corresponds to the Hartree-Fock model zero-order Hamiltonian, that is the Mpller-Plesset expansion whereas iz=2 identifies the shifted denominator scheme which uses the Epstein-Nesbet zero-order Hamiltonian, ediag will be used to store the diagonal component. These energies are stored in the common block ptres together with the orbital energy (eorb (60)). 

The denominator factors are now required. They are stored in the arrays dl, d2, d3, d4, d5 corresponding to the different spin cases labelled by // = 1, 2, 3,4, 5. This program handles both the Mpller-Plesset and the Epstein-Nesbet perturbation series. For the Mpller-Plesset expansion, the denominators do not depend on the spin case and are given by... 

The elements of the arrays fl, f2, f3, f4, f5 can now be updated for both the Mpller-Plesset expansion and the Epstein-Nesbet series. 

In conclusion, the Mpller-Plesset expansion is a rather unpredictable one, which may diverge even for simple systems. Because of the special form of the zero-order operator, the M0ller-Plesset series will either alternate or otherwise oscillate with long periods. Nevertheless, the low-order corrections may be quite useful in many situations, as demonstrated in Chapter 15. 

Coupled cluster is closely connected with Mpller-Plesset perturbation theory, as mentioned at the start of this section. The infinite Taylor expansion of the exponential operator (eq. (4.46)) ensures that the contributions from a given excitation level are included to infinite order. Perturbation theory indicates that doubles are the most important, they are the only contributors to MP2 and MP3. At fourth order, there are contributions from singles, doubles, triples and quadruples. The MP4 quadruples... 

In SOPPA [5] a Mpller-Plesset perturbation theory expansion of the wave function [28,48] is employed ... 

The other area in which projection of an unrestricted result has received attention is projected Mpller-Plesset perturbation theory, the PUMPn methods [31], where n is the order of the perturbation theory. In cases in which the UHF approximation is a poor starting point (considerable spin contamination, for example), the convergence of the MP perturbation expansion can be slow and/or erratic. The PUMP methods apply projection operators to the perturbation expansion, although usually not full projection but simply annihilation of the leading contaminants. This approach has met with mixed success again, it represents a rather expensive modification to a technique that was originally chosen partly for its economy — seldom a recipe for success. 

MPn (Mpller-Plesset Perturbation Theory to Order n) 200, 206, 321 Mulliken polulation indices 182 Mulliken population analysis 229, 316 Multiple minima 52 Multipole expansion 270 Multipole moment 269 Mutual potential energy 27, 62... 

The Hartree—Fock model throws all its effort into obtaining the best possible one term expansion Do = 1, Dk = 0 for K > Q. The Configuration Interaction and Mpller—Plesset methods improve on this single-term model by... 

The evaluation of the derivatives of the energy of both the exact solutions of the Schrbdinger equation and variational approximations to it axe considerably simplified by the fact that, either the associated wavefunction solves the Schrbdinger equation, or the derivatives of the energy with respect to the variational peireime-ters are all zero. If we have an approximate solution to the Schrbdinger equation which is not variational then the situation is, in principle, much more complicated there are more terms to evaluate in the derivative expression. The problem is that many of the most useful and common approximations which are more accurate than the SCF HF wavefunction are exactly of this type perturbation expansions in terms of the SCF MOs, in particular, Mpller Plesset (many-body) perturbation expansions. 

The non-vanishing of the energy derivatives with respect to the linear expansion coefificients means that the full expression for the derivative of the energy, eqn ( 30.3) must be used, which demands a knowledge of many more derivatives in view of the large number of linear parameters involved in expansion of the orbitals. However, in studies of the evaluations of the derivatives of the second-order Mpller-Plesset energy an important technical result was obtained by Handy and Schaefer which drastically reduces the work involved in this (and other) models of electronic structure. 

As long as one stays with traditional methods for calculation of correlation energies, it is necessary to perform a transformation of the one- and two-electron integrals into the molecular spinor basis. While Mpller-Plesset perturbation expansions can be cast in a semi-direct form that does not require a complete integral transformation, it is... 

The use of low-order perturbation theory is probably the cheapest and conceptually simplest method for including correlation effects in a quantum-chemical calculation while maintaining a minimum of formal rigor. In particular, Mpller-Plesset perturbation expansions to various orders (commonly denoted MPn) have seen widespread use. For our purposes it is sufficient to discuss only the MP2 expansion, which is the lowest order that contributes beyond the mean-field approximation. 

As noted above, the theoretical challenges to be overcome include the validation of existing theories as well as the development of new theories. One of the surprises in electronic structure theory in the 1990s was the finding that Mpller-Plesset perturbation theory, the most widely used means to include electron correlation effects, does not lead to a convergent perturbation expansion series. This... 

A critical step in the development of any perturbation expansion is the division of the Hamiltonian 5f into a zero-order part Jfo and a perturbation JKi. In second quantized formalism, the Hamiltonian for the Mpller-Plesset theory is written as... 

---