saving work on Pade

This commit is contained in:
Pierre-Francois Loos 2020-11-25 09:58:06 +01:00
parent 00f17592ea
commit d0cf7f59e3
3 changed files with 141 additions and 17 deletions

View File

@ -6,7 +6,7 @@
%Control: page (0) single %Control: page (0) single
%Control: year (1) truncated %Control: year (1) truncated
%Control: production of eprint (0) enabled %Control: production of eprint (0) enabled
\begin{thebibliography}{132}% \begin{thebibliography}{138}%
\makeatletter \makeatletter
\providecommand \@ifxundefined [1]{% \providecommand \@ifxundefined [1]{%
\@ifx{#1\undefined} \@ifx{#1\undefined}
@ -1277,16 +1277,63 @@
{journal} {\bibinfo {journal} {J. Phys. A: Math. Theor.}\ }\textbf {\bibinfo {journal} {\bibinfo {journal} {J. Phys. A: Math. Theor.}\ }\textbf {\bibinfo
{volume} {40}},\ \bibinfo {pages} {581} (\bibinfo {year} {2007})}\BibitemShut {volume} {40}},\ \bibinfo {pages} {581} (\bibinfo {year} {2007})}\BibitemShut
{NoStop}% {NoStop}%
\bibitem [{\citenamefont {{Baker Jr.}}\ and\ \citenamefont
{Graves-Morris}(1996)}]{BakerBook}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {G.~A.}\ \bibnamefont
{{Baker Jr.}}}\ and\ \bibinfo {author} {\bibfnamefont {P.}~\bibnamefont
{Graves-Morris}},\ }\href@noop {} {\emph {\bibinfo {title} {Pad\'e
Approximants}}}\ (\bibinfo {publisher} {Cambridge University Press},\
\bibinfo {year} {1996})\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Pad\'e}(1892)}]{Pade_1892}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {H.}~\bibnamefont
{Pad\'e}},\ }in\ \href@noop {} {\emph {\bibinfo {booktitle} {Annales
scientifiques de l'{\'E}.N.S.}}},\ Vol.~\bibinfo {volume} {9},\ \bibinfo
{editor} {edited by\ \bibinfo {editor} {\bibnamefont {Gauthier-Villars}}}\
(\bibinfo {publisher} {{\'E}ditions scientifiques et m{\'e}dicales
Elsevier},\ \bibinfo {year} {1892})\ pp.\ \bibinfo {pages}
{3--93}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Loos}(2013)}]{Loos_2013}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {P.-F.}\ \bibnamefont
{Loos}},\ }\href {\doibase 10.1063/1.4790613} {\bibfield {journal} {\bibinfo
{journal} {J. Chem. Phys.}\ }\textbf {\bibinfo {volume} {138}},\ \bibinfo
{pages} {064108} (\bibinfo {year} {2013})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Gluzman}(2020)}]{Gluzman_2020}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {S.}~\bibnamefont
{Gluzman}},\ }\href {\doibase 10.3390/sym12101600} {\bibfield {journal}
{\bibinfo {journal} {Symmetry}\ }\textbf {\bibinfo {volume} {12}},\ \bibinfo
{pages} {1600} (\bibinfo {year} {2020})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Mayer}\ and\ \citenamefont
{Tong}(1985)}]{Mayer_1985}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {I.~L.}\ \bibnamefont
{Mayer}}\ and\ \bibinfo {author} {\bibfnamefont {B.~Y.}\ \bibnamefont
{Tong}},\ }\href {\doibase 10.1088/0022-3719/18/17/008} {\bibfield {journal}
{\bibinfo {journal} {J. Phys. C.: Solid State Phys.}\ }\textbf {\bibinfo
{volume} {18}},\ \bibinfo {pages} {3297} (\bibinfo {year}
{1985})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Goodson}(2019)}]{Goodson_2019}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {D.~Z.}\ \bibnamefont
{Goodson}},\ }in\ \href {\doibase
https://doi.org/10.1016/B978-0-12-813651-5.00009-7} {\emph {\bibinfo
{booktitle} {Mathematical Physics in Theoretical Chemistry}}},\ \bibinfo
{series and number} {Developments in Physical {\&} Theoretical Chemistry},\
\bibinfo {editor} {edited by\ \bibinfo {editor} {\bibfnamefont
{S.}~\bibnamefont {Blinder}}\ and\ \bibinfo {editor} {\bibfnamefont
{J.}~\bibnamefont {House}}}\ (\bibinfo {publisher} {Elsevier},\ \bibinfo
{year} {2019})\ pp.\ \bibinfo {pages} {295 -- 325}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Surj{\'a}n}\ and\ \citenamefont \bibitem [{\citenamefont {Surj{\'a}n}\ and\ \citenamefont
{Szabados}(2000)}]{Surjan_2000}% {Szabados}(2000)}]{Surjan_2000}%
\BibitemOpen \BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {P.~R.}\ \bibnamefont \bibfield {author} {\bibinfo {author} {\bibfnamefont {P.~R.}\ \bibnamefont
{Surj{\'a}n}}\ and\ \bibinfo {author} {\bibfnamefont {{\'A}.}~\bibnamefont {Surj{\'a}n}}\ and\ \bibinfo {author} {\bibfnamefont {{\'A}.}~\bibnamefont
{Szabados}},\ }\href {\doibase 10.1063/1.481006} {\bibfield {journal} {Szabados}},\ }\href {\doibase 10.1063/1.481006} {\bibfield {journal}
{\bibinfo {journal} {The Journal of Chemical Physics}\ }\textbf {\bibinfo {\bibinfo {journal} {J. Chem. Phys.}\ }\textbf {\bibinfo {volume} {112}},\
{volume} {112}},\ \bibinfo {pages} {4438} (\bibinfo {year} {2000})},\ \Eprint \bibinfo {pages} {4438} (\bibinfo {year} {2000})}\BibitemShut {NoStop}%
{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 \bibitem [{\citenamefont {Surj{\'a}n}\ \emph {et~al.}(2018)\citenamefont
{Surj{\'a}n}, \citenamefont {Mih{\'a}lka},\ and\ \citenamefont {Surj{\'a}n}, \citenamefont {Mih{\'a}lka},\ and\ \citenamefont
{Szabados}}]{Surjan_2018}% {Szabados}}]{Surjan_2018}%

View File

@ -1,24 +1,97 @@
%% This BibTeX bibliography file was created using BibDesk. %% This BibTeX bibliography file was created using BibDesk.
%% http://bibdesk.sourceforge.net/ %% http://bibdesk.sourceforge.net/
%% Created for Pierre-Francois Loos at 2020-11-25 08:29:13 +0100 %% Created for Pierre-Francois Loos at 2020-11-25 09:35:42 +0100
%% Saved with string encoding Unicode (UTF-8) %% Saved with string encoding Unicode (UTF-8)
@article{Loos_2013,
author = {Loos, Pierre-Fran{\c c}ois},
date-added = {2020-11-25 09:34:55 +0100},
date-modified = {2020-11-25 09:35:07 +0100},
doi = {10.1063/1.4790613},
journal = {J. Chem. Phys.},
pages = {064108},
title = {High-Density Correlation Energy Expansion of the One-Dimensional Uniform Electron Gas},
volume = {138},
year = {2013},
Bdsk-Url-1 = {https://dx.doi.org/10.1063/1.4790613}}
@article{Gluzman_2020,
author = {S. Gluzman},
date-added = {2020-11-25 09:32:54 +0100},
date-modified = {2020-11-25 09:33:54 +0100},
doi = {10.3390/sym12101600},
journal = {Symmetry},
pages = {1600},
title = {Pad\'e and Post-Pad\'e Approximations for Critical Phenomena},
volume = {12},
year = {2020}}
@incollection{Goodson_2019,
abstract = {The Schr{\"o}dinger equation for an atom or molecule includes parameters, such as bond lengths or nuclear charges, and the resulting energy eigenvalue can be treated as a function with the parameter values as continuous variables. Analysis of singular points of this function, at nonphysical parameter values, can explain and predict the success or failure of quantum chemical calculation methods. An introduction to the theory of singularities in functions of a complex variable is presented and examples of applications to quantum chemistry are described, including the calculation of molecular potential energy curves, the theoretical description of ionization, and the summation of perturbation theories.},
author = {David Z. Goodson},
booktitle = {Mathematical Physics in Theoretical Chemistry},
date-added = {2020-11-25 09:13:38 +0100},
date-modified = {2020-11-25 09:14:27 +0100},
doi = {https://doi.org/10.1016/B978-0-12-813651-5.00009-7},
editor = {S.M. Blinder and J.E. House},
isbn = {978-0-12-813651-5},
keywords = {Singularities, Avoided crossings, Quadratic approximants, Molecular potential energy curves, Ionization, Finite-size scaling, Perturbation theory, Series summation},
pages = {295 - 325},
publisher = {Elsevier},
series = {Developments in Physical {\&} Theoretical Chemistry},
title = {Chapter 9 - Singularity analysis in quantum chemistry},
year = {2019},
Bdsk-Url-1 = {http://www.sciencedirect.com/science/article/pii/B9780128136515000097},
Bdsk-Url-2 = {https://doi.org/10.1016/B978-0-12-813651-5.00009-7}}
@article{Mayer_1985,
abstract = {The quadratic Pade method-a new method for calculating the local density of states in various physical systems-is introduced and discussed. The method is based upon the use of Hermite-Pade polynomials and it makes the calculation of densities of states a straightforward and relatively simple matter. Its advantages over other methods with similar generality and complexity are outlined and numerical results for various systems, which illustrate the virtues of the new method, are presented and discussed.},
author = {I L Mayer and B Y Tong},
date-added = {2020-11-25 09:01:38 +0100},
date-modified = {2020-11-25 09:03:36 +0100},
doi = {10.1088/0022-3719/18/17/008},
journal = {J. Phys. C.: Solid State Phys.},
pages = {3297--3318},
title = {The quadratic Pade approximant method and its application for calculating densities of states},
volume = {18},
year = 1985,
Bdsk-Url-1 = {https://doi.org/10.1088%2F0022-3719%2F18%2F17%2F008},
Bdsk-Url-2 = {https://doi.org/10.1088/0022-3719/18/17/008}}
@book{BakerBook,
author = {G. A. {Baker Jr.} and P. Graves-Morris},
date-added = {2020-11-25 08:58:42 +0100},
date-modified = {2020-11-25 08:59:33 +0100},
publisher = {Cambridge University Press},
title = {Pad\'e Approximants},
year = {1996}}
@incollection{Pade_1892,
author = {H. Pad\'e},
booktitle = {Annales scientifiques de l'{\'E}.N.S.},
date-added = {2020-11-25 08:48:36 +0100},
date-modified = {2020-11-25 08:56:44 +0100},
editor = {Gauthier-Villars},
pages = {3--93},
publisher = {{\'E}ditions scientifiques et m{\'e}dicales Elsevier},
title = {Sur la repr{\'e}sentation approch{\'e}e d'une fonction par des fractions rationnelles},
volume = {9},
year = {1892}}
@article{Surjan_2000, @article{Surjan_2000,
author = {Surj{\'a}n,P. R. and Szabados,{\'A}.}, author = {Surj{\'a}n,P. R. and Szabados,{\'A}.},
date-added = {2020-11-25 08:29:04 +0100}, date-added = {2020-11-25 08:29:04 +0100},
date-modified = {2020-11-25 08:29:13 +0100}, date-modified = {2020-11-25 09:15:28 +0100},
doi = {10.1063/1.481006}, doi = {10.1063/1.481006},
eprint = {https://doi.org/10.1063/1.481006}, journal = {J. Chem. Phys.},
journal = {The Journal of Chemical Physics},
number = {10}, number = {10},
pages = {4438-4446}, pages = {4438-4446},
title = {Optimized partitioning in perturbation theory: Comparison to related approaches}, title = {Optimized partitioning in perturbation theory: Comparison to related approaches},
url = {https://doi.org/10.1063/1.481006},
volume = {112}, volume = {112},
year = {2000}, year = {2000},
Bdsk-Url-1 = {https://doi.org/10.1063/1.481006}} Bdsk-Url-1 = {https://doi.org/10.1063/1.481006}}

View File

@ -1161,27 +1161,31 @@ Those excited states have a non-negligible contribution to the exact FCI solutio
We believe that $\alpha$ singularities are connected to states with non-negligible contribution in the CI expansion thus to the dynamical part of the correlation energy, while $\beta$ singularities are linked to symmetry breaking and phase transitions of the wave function, \ie, to the multi-reference nature of the wave function thus to the static part of the correlation energy. We believe that $\alpha$ singularities are connected to states with non-negligible contribution in the CI expansion thus to the dynamical part of the correlation energy, while $\beta$ singularities are linked to symmetry breaking and phase transitions of the wave function, \ie, to the multi-reference nature of the wave function thus to the static part of the correlation energy.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Exploiting complex analysis} \section{Resummation techniques}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
As frequently claimed by Carl Bender, the most stupid thing that one can do with a series is to sum it.
Nonetheless, quantum chemists are basically doing exactly this on a daily basis.
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. 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} \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''}. According to Wikipedia, \textit{``a Pad\'e approximant is the best approximation of a function by a rational function of given order''}. \cite{BakerBook}
A $[d_A/d_B]$ Pad\'e approximant is defined as \cite{Pade} A $[d_A/d_B]$ Pad\'e approximant is defined as \cite{Pade_1892}
\begin{equation} \begin{equation}
\label{eq:PadeApp} \label{eq:PadeApp}
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} 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} \end{equation}
(with $b_0 = 1$), where the coefficients of the polynomials are determined by collecting terms according to power of $\lambda$. (with $b_0 = 1$), where the coefficients of the polynomials $A(\lambda)$ and $B(\lambda)$ 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. Pad\'e approximants are extremely useful in many areas of physics and chemistry \cite{Loos_2013,Gluzman_2020} as they can model poles [the poles of a Pad\'e approximant appears at the roots of the polynomial $B(\lambda)$].
However, they are unable to model functions with square-root branch points, which are ubiquitous in the singularity structure of a typical perturbative treatment.
%==========================================% %==========================================%
\subsection{Quadratic approximant} \subsection{Quadratic approximant}
%==========================================% %==========================================%
In a nutshell, the idea behind quadratic approximant is to model the singularity structure of the energy function $E(\lambda)$ via a generalized version of the square-root singularity expression In a nutshell, the idea behind quadratic approximant is to model the singularity structure of the energy function $E(\lambda)$ via a generalized version of the square-root singularity expression \cite{Mayer_1985,Goodson_2011,Goodson_2019}
\begin{equation} \begin{equation}
\label{eq:QuadApp} \label{eq:QuadApp}
E(\lambda) = \frac{1}{2 Q(\lambda)} \qty[ P(\lambda) \pm \sqrt{P^2(\lambda) - 4 Q(\lambda) R(\lambda)} ] E(\lambda) = \frac{1}{2 Q(\lambda)} \qty[ P(\lambda) \pm \sqrt{P^2(\lambda) - 4 Q(\lambda) R(\lambda)} ]
@ -1195,7 +1199,7 @@ where
& &
R(\lambda) & = \sum_{k=0}^{d_R} r_k \lambda^k R(\lambda) & = \sum_{k=0}^{d_R} r_k \lambda^k
\end{align} \end{align}
are polynomials, such that $d_P + d_Q + d_R = n - 1$ where $n$ is the highest-order series coefficient known from the Taylor expansion of $E(\lambda)$. are polynomials, such that $d_P + d_Q + d_R = n - 1$, where $n$ is the highest-order series coefficient known from the Taylor expansion of $E(\lambda)$.
Recasting Eq.~\eqref{eq:QuadApp} as a second-order expression in $E(\lambda)$, \ie, Recasting Eq.~\eqref{eq:QuadApp} as a second-order expression in $E(\lambda)$, \ie,
\begin{equation} \begin{equation}
Q(\lambda) E^2(\lambda) - P(\lambda) E(\lambda) + R(\lambda) \sim \order*{\lambda^{n+1}} Q(\lambda) E^2(\lambda) - P(\lambda) E(\lambda) + R(\lambda) \sim \order*{\lambda^{n+1}}