temp commit

This commit is contained in:
Hugh Burton 2020-11-24 09:56:03 +00:00
parent 2625cabe0e
commit d14b5bc580
2 changed files with 39 additions and 1 deletions

View File

@ -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}).

View File

@ -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):
-----------------------------------------------------------