saving work on Pade

This commit is contained in:
Pierre-Francois Loos 2020-11-25 08:48:07 +01:00
parent 8281a2e66b
commit 00f17592ea
3 changed files with 172 additions and 128 deletions

View File

@ -6,7 +6,7 @@
%Control: page (0) single
%Control: year (1) truncated
%Control: production of eprint (0) enabled
\begin{thebibliography}{131}%
\begin{thebibliography}{132}%
\makeatletter
\providecommand \@ifxundefined [1]{%
\@ifx{#1\undefined}
@ -1277,6 +1277,16 @@
{journal} {\bibinfo {journal} {J. Phys. A: Math. Theor.}\ }\textbf {\bibinfo
{volume} {40}},\ \bibinfo {pages} {581} (\bibinfo {year} {2007})}\BibitemShut
{NoStop}%
\bibitem [{\citenamefont {Surj{\'a}n}\ and\ \citenamefont
{Szabados}(2000)}]{Surjan_2000}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {P.~R.}\ \bibnamefont
{Surj{\'a}n}}\ and\ \bibinfo {author} {\bibfnamefont {{\'A}.}~\bibnamefont
{Szabados}},\ }\href {\doibase 10.1063/1.481006} {\bibfield {journal}
{\bibinfo {journal} {The Journal of Chemical Physics}\ }\textbf {\bibinfo
{volume} {112}},\ \bibinfo {pages} {4438} (\bibinfo {year} {2000})},\ \Eprint
{http://arxiv.org/abs/https://doi.org/10.1063/1.481006}
{https://doi.org/10.1063/1.481006} \BibitemShut {NoStop}%
\bibitem [{\citenamefont {Surj{\'a}n}\ \emph {et~al.}(2018)\citenamefont
{Surj{\'a}n}, \citenamefont {Mih{\'a}lka},\ and\ \citenamefont
{Szabados}}]{Surjan_2018}%

View File

@ -1,146 +1,178 @@
%% This BibTeX bibliography file was created using BibDesk.
%% http://bibdesk.sourceforge.net/
%% Created for Pierre-Francois Loos at 2020-11-24 09:46:03 +0100
%% Created for Pierre-Francois Loos at 2020-11-25 08:29:13 +0100
%% Saved with string encoding Unicode (UTF-8)
@article{Surjan_2000,
author = {Surj{\'a}n,P. R. and Szabados,{\'A}.},
date-added = {2020-11-25 08:29:04 +0100},
date-modified = {2020-11-25 08:29:13 +0100},
doi = {10.1063/1.481006},
eprint = {https://doi.org/10.1063/1.481006},
journal = {The Journal of Chemical Physics},
number = {10},
pages = {4438-4446},
title = {Optimized partitioning in perturbation theory: Comparison to related approaches},
url = {https://doi.org/10.1063/1.481006},
volume = {112},
year = {2000},
Bdsk-Url-1 = {https://doi.org/10.1063/1.481006}}
@article{Tsuchimochi_2019,
author = {Takashi Tsuchimochi and Seiichiro L. Ten-no},
doi = {10.1021/acs.jctc.9b00897},
journal = {J. Chem. Theory Comput.},
volume ={15},
pages = {6688},
title = {Second-order perturbation theory with spin-symmetry-projected Hartree--Fock},
volume = {15},
year = {2019},
doi ={10.1021/acs.jctc.9b00897},
}
Bdsk-Url-1 = {https://doi.org/10.1021/acs.jctc.9b00897}}
@article{Tsuchimochi_2014,
author = {Takashi Tsuchimochi and Troy {Van Voorhis}},
doi = {10.1063/1.4898804},
journal = {J. Chem. Phys.},
volume ={141},
pages = {164117},
title = {Extended {M\oller--Plesset} perturbation theory for dynamical and static correlations},
volume = {141},
year = {2014},
doi ={10.1063/1.4898804},
}
Bdsk-Url-1 = {https://doi.org/10.1063/1.4898804}}
@article{Knowles_1988a,
author = {Peter J. Knowles and Nicholas C. Handy},
doi = {10.1021/j100322a018},
journal = {J. Phys. Chem.},
volume ={92},
pages = {3097},
title = {Convergence of projected unrestricted {Hartree--Fock M\oller--Plesset} series},
volume = {92},
year = {1988},
doi ={10.1021/j100322a018},
}
Bdsk-Url-1 = {https://doi.org/10.1021/j100322a018}}
@article{Knowles_1988b,
author = {Peter J. Knowles and Nicholas C. Handy},
doi = {10.1063/1.454397},
journal = {J. Chem. Phys.},
volume ={88},
pages = {6991},
title = {Projected unrestricted {M\oller--Plesset} second-order energies},
volume = {88},
year = {1988},
doi ={10.1063/1.454397},
}
Bdsk-Url-1 = {https://doi.org/10.1063/1.454397}}
@article{Schlegel_1986,
author = {H. Bernhard Schlegel},
doi = {10.1063/1.450026},
journal = {J. Chem. Phys.},
volume ={84},
pages = {4530},
title = {Potential energy curves using unrestricted {M\oller--Plesset} perturbation theory with spin annihilation},
volume = {84},
year = {1986},
doi ={10.1063/1.450026},
}
Bdsk-Url-1 = {https://doi.org/10.1063/1.450026}}
@article{Schlegel_1988,
author = {H. Bernhard Schlegel},
doi = {10.1021/j100322a014},
journal = {J. Phys. Chem.},
volume ={91},
pages = {3075},
title = {{M\oller--Plesset} perturbation theory with spin projection},
volume = {91},
year = {1988},
doi ={10.1021/j100322a014},
}
Bdsk-Url-1 = {https://doi.org/10.1021/j100322a014}}
@article{Gill_1988a,
author = {P. M. W. Gill and M. W. Wong and R. H. Nobes and L. Radom},
doi = {10.1016/0009-2614(88)80328-2},
journal = {Chem. Phys. Lett.},
volume ={148},
pages = {541},
title = {How well can {RMP4} theory treat homolytic fragmentations?},
volume = {148},
year = {1988},
doi ={10.1016/0009-2614(88)80328-2},
}
Bdsk-Url-1 = {https://doi.org/10.1016/0009-2614(88)80328-2}}
@article{Nobes_1987,
author = {R. H. Nobes and J. A. Pople and L. Radom and N. C. Handy and P. J. Knowles},
doi = {10.1016/0009-2614(87)80545-6},
journal = {Chem. Phys. Lett.},
volume ={138},
pages = {481},
title = {Slow convergence of the {M\oller--Plesset} perturbation series: the dissociation energy of hydrogen cyanide and the electron affinity of the cyano radical},
volume = {138},
year = {1987},
doi ={10.1016/0009-2614(87)80545-6},
}
Bdsk-Url-1 = {https://doi.org/10.1016/0009-2614(87)80545-6}}
@article{Laidig_1987,
author = {William D. Laidig and Paul Saxe and Rodney J. Bartlett},
doi = {10.1063/1.452291},
journal = {J. Chem. Phys.},
volume ={86},
pages = {887},
title = {The description of \ce{N2} and \ce{F2} potential energy surfaces using multireference coupled cluster theory},
volume = {86},
year = {1987},
doi ={10.1063/1.452291},
}
Bdsk-Url-1 = {https://doi.org/10.1063/1.452291}}
@article{Bartlett_1975,
author = {R. J. Bartlett and D. M. Silver},
doi = {10.1063/1.430878},
journal = {J. Chem. Phys.},
volume ={62},
pages = {3258},
title = {Many-body perturbation theory applied to electron pair correlation energies. I. Closed-shell first-row diatomic hydrides},
volume = {62},
year = {1975},
doi ={10.1063/1.430878},
}
Bdsk-Url-1 = {https://doi.org/10.1063/1.430878}}
@article{Krishnan_1980,
author = {R. Krishnan and M. J. Frisch and J. A. Pople},
doi = {10.1063/1.439657},
journal = {J. Chem. Phys.},
volume ={72},
pages = {4244},
title = {Contribution of triple substitutions to the electron correlation energy in fourth order perturbation theory},
volume = {72},
year = {1980},
doi ={10.1063/1.439657},
}
Bdsk-Url-1 = {https://doi.org/10.1063/1.439657}}
@article{Pople_1978,
author = {J. A. Pople and R. Krishnan and H. B. Schlegel and J. S. Binkley},
doi = {10.1002/qua.560140503},
journal = {Int. J. Quantum Chem.},
volume ={14},
pages = {545},
title = {Electron correlation theories and their application to the study of simple reaction potential surfaces},
volume = {14},
year = {1978},
doi ={10.1002/qua.560140503},
}
Bdsk-Url-1 = {https://doi.org/10.1002/qua.560140503}}
@article{Pople_1976,
author = {John A. Pople and Stephen Binkley and Rolf Seeger},
doi = {10.1002/qua.560100802},
journal = {Int. J. Quantum Chem. Symp.},
volume ={10},
pages = {1},
title = {Theoretical models incorporating electron correlation},
volume = {10},
year = {1976},
doi ={10.1002/qua.560100802},
}
Bdsk-Url-1 = {https://doi.org/10.1002/qua.560100802}}
@article{Knowles_1985,
author = {P. J. Knowles and K. Somasundram and N. C. Handy and K. Hirao},
doi = {10.1016/0009-2614(85)85002-8},
journal = {Chem. Phys. Lett.},
volume ={113},
pages = {8},
title = {The Calculation of High-Order Energies in the Many-Body Perturbation Theory Series},
volume = {113},
year = {1985},
doi ={10.1016/0009-2614(85)85002-8},
}
Bdsk-Url-1 = {https://doi.org/10.1016/0009-2614(85)85002-8}}
@article{Laidig_1985,
author = {William D. Laidig and George Fitzgerald and Rodney J. Bartlett},
doi = {10.1016/0009-2614(85)80934-9},
journal = {Chem. Phys. Lett.},
volume ={113},
pages = {151},
title = {Is Fifth-Order MBPT Enough?},
volume = {113},
year = {1985},
doi ={10.1016/0009-2614(85)80934-9}
}
Bdsk-Url-1 = {https://doi.org/10.1016/0009-2614(85)80934-9}}
@article{Hall_1951,
abstract = { An analysis of the `linear combination of atomic orbitals' approximation using the accurate molecular orbital equations shows that it does not lead to equations of the form usually assumed in the semi-empirical molecular orbital method. A new semi-empirical method is proposed, therefore, in terms of equivalent orbitals. The equations obtained, which do have the usual form, are applicable to a large class of molecules and do not involve the approximations that were thought necessary. In this method the ionization potentials are calculated by treating certain integrals as semi-empirical parameters. The value of these parameters is discussed in terms of the localization of equivalent orbitals and some approximate rules are suggested. As an illustration the ionization potentials of the paraffin series are considered and good agreement between the observed and calculated values is found. },
author = {Hall, G. G. and Lennard-Jones, John Edward},

View File

@ -1164,17 +1164,19 @@ We believe that $\alpha$ singularities are connected to states with non-negligib
\section{Exploiting complex analysis}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
The inability of Taylor series to model properly the energy function $E(\lambda$) can be simply understood by the fact that one aims at modelling a complicated function with potentially poles and singularities by a simple polynomial of finite order.
%==========================================%
\subsection{Pad\'e approximant}
%==========================================%
According to Wikipedia, \textit{``a Pad\'e approximant is the best approximation of a function by a rational function of given order''}.
A $[d_A/d_B]$ Pad\'e approximant is defined as
A $[d_A/d_B]$ Pad\'e approximant is defined as \cite{Pade}
\begin{equation}
\label{eq:PadeApp}
E(\lambda) = \frac{A(\lambda)}{B(\lambda)} = \frac{\sum_{k=0}^{d_A} a_k \lambda^k}{\sum_{k=0}^{d_B} b_k \lambda^k}
E_{[d_A/d_B]}(\lambda) = \frac{A(\lambda)}{B(\lambda)} = \frac{\sum_{k=0}^{d_A} a_k \lambda^k}{\sum_{k=0}^{d_B} b_k \lambda^k}
\end{equation}
with $b_0 = 1$.
Pad\'e approximants are nice and they can model poles, but they cannot model functions with square-root branch points, and that sucks.
(with $b_0 = 1$), where the coefficients of the polynomials are determined by collecting terms according to power of $\lambda$.
Pad\'e approximants are nice and they can model poles (the poles of a $[d_A/d_B]$ Pad\'e approximant appears at the roots of the polynomial $B(\lambda)$), but they cannot model functions with square-root branch points, and that sucks because it is not consistent with the singularity structure of a typical perturbative treatment.
%==========================================%
\subsection{Quadratic approximant}
@ -1212,7 +1214,7 @@ However, by ramping up high enough the degree of the polynomials, one is able to
\subsection{Analytic continuation}
%==========================================%
Recently, Mih\'alka \textit{et al.} studied the partitioning effect on the convergence properties of Rayleigh-Schr\"odinger perturbation theory by considering the MP and the EN partitioning as well as an alternative partitioning. \cite{Mihalka_2017a}
Recently, Mih\'alka \textit{et al.} studied the partitioning effect on the convergence properties of Rayleigh-Schr\"odinger perturbation theory by considering the MP and the EN partitioning as well as an alternative partitioning \cite{Mihalka_2017a} (see also Ref.~\onlinecite{Surjan_2000}).
Taking as an example (in particular) the water molecule at equilibrium and at stretched geometries, they could estimate the radius of convergence via a quadratic Pad\'e approximant and convert divergent perturbation expansions to convergent ones in some cases thanks to a judicious choice of the level shift parameter.
In a subsequent study by the same group, \cite{Mihalka_2017b} they use analytic continuation techniques to resum divergent MP series \cite{Goodson_2011} taking again as an example the water molecule in a stretched geometry.
In a nutshell, their idea consists in calculating the energy of the system for several values of $\lambda$ for which the MP series is rapidly convergent (\ie, for $\lambda < r_c$), and to extrapolate the final energy to the physical system at $\lambda = 1$ via a polynomial- or Pad\'e-based fit.