diff --git a/Manuscript/EPAWTFT.tex b/Manuscript/EPAWTFT.tex index 21f3cad..fd710cc 100644 --- a/Manuscript/EPAWTFT.tex +++ b/Manuscript/EPAWTFT.tex @@ -708,6 +708,33 @@ to the convergence properties of the MP expansion.} \subsection{M{\o}ller-Plesset Convergence in Molecular Systems} %=====================================================% +% GENERAL DESIRE FOR WELL-BEHAVED CONVERGENCE AND LOW-ORDER TERMS +\hugh{Among the most desirable properties of any electronic structure technique is the existence of +a systematic route to increasingly accurate energies. +In the context of MP theory, one would like a monotonic convergence of the perturbation +series towards the exact energy such that the accuracy increases as each term in the series is added. +If such well-behaved convergence can be established, then our ability to compute individual +terms in the series becomes the only barrier to computing the exact correlation in a finite basis set. +Unfortunately, the computational scaling of each term in the MP series increases with the perturbation +order, and practical calculations rely on fast convergence +to obtain high-accuracy results using only the lowest order terms. +} + +% INITIAL POSITIVITY AROUND THE CONVERGENCE PROPERTIES +\hugh{MP theory was first introduced to quantum chemistry through the pioneering +works of Bartlett \etal{},\cite{Bartlett_1975} and Pople \etal{}.\cite{Pople_1976,Pople_1978} +The first investigations into the convergence of the MP series followed soon after, and +it was quickly realised that this convergence was often very slow, or erratic. +} +Given the + + + +% WORRIES ABOUT ERRATIC OR SLOW CONVERGENCE + + +% CREMER AND HE + When one relies on MP perturbation theory (and more generally on any perturbative partitioning), it would be reasonable to ask for a systematic improvement of the energy with respect to the perturbative order, \ie, one would expect that the more terms of the perturbative series one can compute, the closer the result from the exact energy. In other words, one would like a monotonic convergence of the MP series. Assuming this, the only limiting process to get the exact correlation energy (in a finite basis set) would be our ability to compute the terms of this perturbation series. Unfortunately this is not as easy as one might think because i) the terms of the perturbative series become rapidly computationally cumbersome, and ii) erratic behaviour of the perturbative coefficients are not uncommon. For example, in the late 80's, Gill and Radom reported deceptive and slow convergences in stretched systems \cite{Gill_1986,Gill_1988} (see also Refs.~\onlinecite{Handy_1985,Lepetit_1988}). diff --git a/hugh_notes.txt b/hugh_notes.txt index 9bc8fed..b3c0503 100644 --- a/hugh_notes.txt +++ b/hugh_notes.txt @@ -56,7 +56,7 @@ low-order expansions give qualitatively wrong energetics (eg. unphysical barrier Schlegel, JCP, (1986): ---------------------- -Apply spin-projection to UHF and UMP to obtain improved potential enerrgy curves. Use a +Apply spin-projection to UHF and UMP to obtain improved potential energy curves. Use a post-perturbation projection to avoid mixing in higher energy states. Perturbation corrections do not significantly reduce spin contamination. PUHF has a gradient @@ -79,6 +79,14 @@ Worst barrier height estimate occurs at UMP4, after which there is very slow con They suspect that UMP problems can be attributed to spin-contamination. Conclude that incorrect "qualitative" nature of RMP is not as bad as spin-contamination in UMP. +Nobes, Pople, Radom, Handy and Knowles, CPL (1987): +--------------------------------------------------- +Look at the MP convergence in the cyanide anion (CN-) at this is a molecule with a symmetry-broken +solution at equilibrium. They observe very slow convergence in these UMP series too, with more than +MP4 being needed for accurate energies. This confirms that spin-contamination, rather than extended +bonds, provides the driving force for slow UMP convergence. + + Gill, Wong, Nobes, and Radom, CPL (1988): ----------------------------------------- Investigate performance of RMP expansions for homolytic bond breaking. @@ -113,3 +121,6 @@ The EN partitioning avoids this, but the PT terms then become undetermined (zero Single excitations can interact with the doubly-excited determinants. This matrix elements goes through a maximum at intermediate distances. This contribution enters at fourth-order. +Raghavachari, Pople, Replogle, and Head-Gordon, JPC, (1990): +----------------------------------------------------------- +