PTEROSOR
/
EPAWTFT
10
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

update review

This commit is contained in:
Pierre-Francois Loos 2020-11-12 12:07:10 +01:00
parent e6b06be8e5
commit 8d9d97231d
158 changed files with 49148 additions and 1612 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

245
Manuscript/EPAWTFT.aux
View File

@ -17,28 +17,29 @@
\providecommand\HyField@AuxAddToFields[1]{}
\providecommand\HyField@AuxAddToCoFields[2]{}
\citation{MoiseyevBook,BenderPTBook}
\citation{Burton_2019}
\citation{Burton_2019,Burton_2019a}
\citation{Bittner_2012,Chong_2011,Chtchelkatchev_2012,Doppler_2016,Guo_2009,Hang_2013,Liertzer_2012,Longhi_2010,Peng_2014,Peng_2014a,Regensburger_2012,Ruter_2010,Schindler_2011,Szameit_2011,Zhao_2010,Zheng_2013,Choi_2018,El-Ganainy_2018}
\citation{Heiss_1990,Heiss_1999,Heiss_2012,Heiss_2016}
\citation{MoiseyevBook,Heiss_1988,Heiss_1990,Heiss_1999,Berry_2011,Heiss_2012,Heiss_2016,Benda_2018}
\citation{Yarkony_1996}
\citation{Benda_2018}
\citation{MoiseyevBook,Heiss_2016,Benda_2018}
\citation{Berry_1984}
\citation{MoiseyevBook}
\citation{BenderBook,Olsen_1996,Olsen_2000,Olsen_2019}
\citation{BenderBook,Olsen_1996,Olsen_2000,Olsen_2019,Mihalka_2017a,Mihalka_2017b,Mihalka_2019}
\citation{SzaboBook}
\newlabel{FirstPage}{{}{1}{}{section*.1}{}}
\@writefile{toc}{\contentsline {title}{Perturbation theory in the complex plane}{1}{section*.2}\protected@file@percent }
\@writefile{toc}{\contentsline {title}{Perturbation theory in the complex plane: Exceptional points and where to find them}{1}{section*.2}\protected@file@percent }
\@writefile{toc}{\contentsline {abstract}{Abstract}{1}{section*.1}\protected@file@percent }
\@writefile{toc}{\contentsline {section}{\numberline {I}Introduction}{1}{section*.3}\protected@file@percent }
\newlabel{sec:intro}{{I}{1}{}{section*.3}{}}
\@writefile{toc}{\contentsline {subsection}{\numberline {A}Background}{1}{section*.4}\protected@file@percent }
\citation{SzaboBook}
\@writefile{toc}{\contentsline {subsection}{\numberline {B}An illustrative example}{1}{section*.5}\protected@file@percent }
\newlabel{eq:H_2x2}{{1}{1}{}{equation.1.1}{}}
\citation{MoiseyevBook}
\citation{Berry_1984}
\citation{MoiseyevBook}
\citation{Mihalka_2017b}
\citation{Goodson_2012}
\citation{BenderBook}
\@writefile{toc}{\contentsline {subsection}{\numberline {B}An illustrative example}{2}{section*.5}\protected@file@percent }
\newlabel{eq:H_2x2}{{1}{2}{}{equation.1.1}{}}
\newlabel{eq:E_2x2}{{2}{2}{}{equation.1.2}{}}
\newlabel{eq:lambda_EP}{{3}{2}{}{equation.1.3}{}}
\newlabel{eq:E_EP}{{4}{2}{}{equation.1.4}{}}
@ -49,12 +50,18 @@
\newlabel{eq:SchrEq-PT}{{9}{2}{}{equation.2.9}{}}
\newlabel{eq:Elambda}{{10}{2}{}{equation.2.10}{}}
\newlabel{eq:DivExample}{{11}{2}{}{equation.2.11}{}}
\@writefile{toc}{\contentsline {subsection}{\numberline {B}The Hartree-Fock Hamiltonian}{2}{section*.8}\protected@file@percent }
\citation{SzaboBook}
\citation{Moller_1934}
\citation{SzaboBook}
\citation{Gill_1994}
\citation{Gill_1986,Gill_1988,Handy_1985,Lepetit_1988}
\citation{JensenBook}
\@writefile{toc}{\contentsline {subsection}{\numberline {B}The Hartree-Fock Hamiltonian}{3}{section*.8}\protected@file@percent }
\citation{Nesbet_1955,Epstein_1926}
\citation{Seidl_2018}
\citation{Gill_1986,Gill_1988}
\citation{Handy_1985,Lepetit_1988}
\citation{Gill_1986}
\newlabel{eq:FockOp}{{12}{3}{}{equation.2.12}{}}
\newlabel{eq:CoulOp}{{15a}{3}{}{equation.2.15alph1}{}}
\newlabel{eq:ExcOp}{{15b}{3}{}{equation.2.15alph2}{}}
@ -62,10 +69,6 @@
\@writefile{toc}{\contentsline {subsection}{\numberline {C}M{\o }ller-Plesset perturbation theory}{3}{section*.9}\protected@file@percent }
\newlabel{eq:MPHamiltonian}{{17}{3}{}{equation.2.17}{}}
\newlabel{eq:EMP2}{{19}{3}{}{equation.2.19}{}}
\citation{Seidl_2018}
\citation{Gill_1986,Gill_1988}
\citation{Handy_1985,Lepetit_1988}
\citation{Gill_1986}
\citation{SzaboBook}
\citation{Gill_1988}
\citation{Handy_1985}
@ -74,12 +77,10 @@
\citation{Cremer_1996}
\citation{Olsen_1996}
\citation{Olsen_1996,Christiansen_1996}
\@writefile{toc}{\contentsline {subsection}{\numberline {D}Alternative partitioning}{4}{section*.10}\protected@file@percent }
\newlabel{sec:AlterPart}{{II\tmspace +\thinmuskip {.1667em}D}{4}{}{section*.10}{}}
\@writefile{toc}{\contentsline {section}{\numberline {III}Historical overview}{4}{section*.11}\protected@file@percent }
\@writefile{toc}{\contentsline {subsection}{\numberline {A}Behavior of the M{\o }ller-Plesset series}{4}{section*.12}\protected@file@percent }
\@writefile{toc}{\contentsline {subsection}{\numberline {B}Cases of divergence}{4}{section*.13}\protected@file@percent }
\citation{Olsen_2000}
\@writefile{toc}{\contentsline {section}{\numberline {III}Historical overview}{4}{section*.10}\protected@file@percent }
\@writefile{toc}{\contentsline {subsection}{\numberline {A}Behavior of the M{\o }ller-Plesset series}{4}{section*.11}\protected@file@percent }
\@writefile{toc}{\contentsline {subsection}{\numberline {B}Cases of divergence}{4}{section*.12}\protected@file@percent }
\citation{Olsen_2019}
\citation{Cremer_1996}
\citation{Sergeev_2005,Sergeev_2006}
@ -88,151 +89,93 @@
\citation{Sergeev_2005}
\citation{Stillinger_2000}
\citation{Olsen_2000}
\@writefile{toc}{\contentsline {subsection}{\numberline {C}The singularity structure}{5}{section*.14}\protected@file@percent }
\newlabel{eq:HamiltonianStillinger}{{22}{5}{}{equation.3.22}{}}
\citation{Goodson_2004}
\citation{Olsen_2000}
\citation{Heiss_1988,Heiss_2002,Cejnar_2005,Cejnar_2007,Cejnar_2009,Borisov_2015,Sindelka_2017}
\citation{Cejnar_2009,Sachdev_2011,Cejnar_2015,Cejnar_2016,Caprio_2008,Macek_2019}
\citation{Cejnar_2009,Sachdev_2011}
\@writefile{toc}{\contentsline {subsection}{\numberline {C}The singularity structure}{5}{section*.13}\protected@file@percent }
\newlabel{eq:HamiltonianStillinger}{{22}{5}{}{equation.3.22}{}}
\citation{Cejnar_2005,Cejnar_2007}
\citation{Stransky_2018}
\citation{Lipkin_1965}
\citation{Ushveridze_1994}
\citation{Thompson_2005,Seidl_2007,Loos_2009b}
\@writefile{toc}{\contentsline {subsection}{\numberline {D}The physics of quantum phase transition}{6}{section*.15}\protected@file@percent }
\@writefile{toc}{\contentsline {section}{\numberline {IV}The spherium model}{6}{section*.16}\protected@file@percent }
\newlabel{sec:spherium}{{IV}{6}{}{section*.16}{}}
\citation{Loos_2009}
\citation{Wigner_1934}
\citation{Coulson_1949}
\citation{Loos_2009}
\citation{AngularBook}
\citation{SzaboBook}
\newlabel{eq:H-sph-omega}{{24}{7}{}{equation.4.24}{}}
\newlabel{eq:RHF_WF}{{25}{7}{}{equation.4.25}{}}
\newlabel{eq:UHF_WF}{{26}{7}{}{equation.4.26}{}}
\newlabel{eq:EHF}{{28}{7}{}{equation.4.28}{}}
\citation{GiulianiBook}
\citation{GiulianiBook}
\citation{Hiscock_2014,Burton_2019,Burton_2019a}
\citation{Cejnar_2007}
\citation{Yarkony_1996}
\newlabel{eq:EsbUHF}{{33}{8}{}{equation.4.33}{}}
\@writefile{lof}{\contentsline {figure}{\numberline {1}{\ignorespaces Energies of the five solutions of the HF equations (multiplied by $R^2$). The dotted curves correspond to the analytic continuation of the symmetry-broken solutions.}}{8}{figure.1}\protected@file@percent }
\newlabel{fig:SpheriumNrj}{{1}{8}{Energies of the five solutions of the HF equations (multiplied by $R^2$). The dotted curves correspond to the analytic continuation of the symmetry-broken solutions}{figure.1}{}}
\@writefile{toc}{\contentsline {section}{\numberline {V}Radius of convergence and exceptional points}{8}{section*.17}\protected@file@percent }
\@writefile{toc}{\contentsline {subsection}{\numberline {A}Evolution of the radius of convergence}{8}{section*.18}\protected@file@percent }
\newlabel{eq:PolChar}{{35a}{8}{}{equation.5.35alph1}{}}
\newlabel{eq:DPolChar}{{35b}{8}{}{equation.5.35alph2}{}}
\citation{AngularBook,SlaterBook,Loos_2009}
\citation{Goodson_2004}
\@writefile{lot}{\contentsline {table}{\numberline {I}{\ignorespaces RHF and sb-UHF energies in the minimal basis and exact energies (in the complete basis) for various $R$. }}{9}{table.1}\protected@file@percent }
\newlabel{tab:ERHFvsEUHF}{{I}{9}{RHF and sb-UHF energies in the minimal basis and exact energies (in the complete basis) for various $R$}{table.1}{}}
\newlabel{eq:rhfbasis}{{36}{9}{}{equation.5.36}{}}
\@writefile{lof}{\contentsline {figure}{\numberline {2}{\ignorespaces Energies $E(\lambda )$ in the restricted basis set \textup {\hbox {\mathsurround \z@ \normalfont (\ignorespaces \ref {eq:rhfbasis}\unskip \@@italiccorr )}} with $R=10$. One can clearly see the avoided crossing between the s\textsuperscript {2} and p\textsubscript {z}\textsuperscript {2} states around $\lambda = 1$.}}{9}{figure.2}\protected@file@percent }
\newlabel{fig:RHFMiniBas}{{2}{9}{Energies $E(\lambda )$ in the restricted basis set \eqref {eq:rhfbasis} with $R=10$. One can clearly see the avoided crossing between the s\textsuperscript {2} and p\textsubscript {z}\textsuperscript {2} states around $\lambda = 1$}{figure.2}{}}
\@writefile{lof}{\contentsline {figure}{\numberline {3}{\ignorespaces Radius of convergence $R_{\text {CV}}$ for two (left) and three (right) basis functions for various partitionings.}}{10}{figure.3}\protected@file@percent }
\newlabel{fig:RadiusPartitioning}{{3}{10}{Radius of convergence $R_{\text {CV}}$ for two (left) and three (right) basis functions for various partitionings}{figure.3}{}}
\@writefile{lot}{\contentsline {table}{\numberline {II}{\ignorespaces Dominant singularity in the CSF basis set ($K=8$) for various value of $R$ in the MP and WC partitioning.}}{10}{table.2}\protected@file@percent }
\newlabel{tab:SingAlpha}{{II}{10}{Dominant singularity in the CSF basis set ($K=8$) for various value of $R$ in the MP and WC partitioning}{table.2}{}}
\@writefile{toc}{\contentsline {subsection}{\numberline {B}Exceptional points in the UHF formalism}{10}{section*.19}\protected@file@percent }
\newlabel{sec:uhfSing}{{V\tmspace +\thinmuskip {.1667em}B}{10}{}{section*.19}{}}
\newlabel{eq:uhfbasis}{{39}{10}{}{equation.5.39}{}}
\@writefile{lof}{\contentsline {figure}{\numberline {4}{\ignorespaces Radius of convergence $R_{\text {CV}}$ in the CSF basis with $K$ basis functions for the MP (left) and WC (right) partitioning.}}{10}{figure.4}\protected@file@percent }
\newlabel{fig:RadiusBasis}{{4}{10}{Radius of convergence $R_{\text {CV}}$ in the CSF basis with $K$ basis functions for the MP (left) and WC (right) partitioning}{figure.4}{}}
\newlabel{eq:uhforbitals}{{41a}{10}{}{equation.5.41alph1}{}}
\newlabel{eq:MatrixElem}{{42}{10}{}{equation.5.42}{}}
\citation{Sergeev_2005}
\citation{Burton_2019a}
\citation{BenderPTBook}
\@writefile{lof}{\contentsline {figure}{\numberline {5}{\ignorespaces Energies $E(\lambda )$ in the unrestricted basis set \textup {\hbox {\mathsurround \z@ \normalfont (\ignorespaces \ref {eq:uhfbasis}\unskip \@@italiccorr )}} with $R=10$.}}{11}{figure.5}\protected@file@percent }
\newlabel{fig:UHFMiniBas}{{5}{11}{Energies $E(\lambda )$ in the unrestricted basis set \eqref {eq:uhfbasis} with $R=10$}{figure.5}{}}
\@writefile{lof}{\contentsline {figure}{\numberline {6}{\ignorespaces Energies $E(\lambda )$ in the unrestricted basis set \textup {\hbox {\mathsurround \z@ \normalfont (\ignorespaces \ref {eq:uhfbasis}\unskip \@@italiccorr )}} for $R=1.5$ (left) and $R=1.51$ (right).}}{11}{figure.6}\protected@file@percent }
\newlabel{fig:UHFEP}{{6}{11}{Energies $E(\lambda )$ in the unrestricted basis set \eqref {eq:uhfbasis} for $R=1.5$ (left) and $R=1.51$ (right)}{figure.6}{}}
\citation{Seidl_2018}
\@writefile{toc}{\contentsline {section}{\numberline {VI}Conclusion}{12}{section*.20}\protected@file@percent }
\@writefile{lof}{\contentsline {figure}{\numberline {7}{\ignorespaces Real part (left) and imaginary part (right) of $E(\lambda )$ in the unrestricted basis set \textup {\hbox {\mathsurround \z@ \normalfont (\ignorespaces \ref {eq:uhfbasis}\unskip \@@italiccorr )}} for $R=1$.}}{12}{figure.7}\protected@file@percent }
\newlabel{fig:UHFPT}{{7}{12}{Real part (left) and imaginary part (right) of $E(\lambda )$ in the unrestricted basis set \eqref {eq:uhfbasis} for $R=1$}{figure.7}{}}
\@writefile{toc}{\contentsline {subsection}{\numberline {D}The physics of quantum phase transitions}{6}{section*.14}\protected@file@percent }
\@writefile{toc}{\contentsline {section}{\numberline {IV}Conclusion}{6}{section*.15}\protected@file@percent }
\bibdata{EPAWTFTNotes,EPAWTFT}
\bibcite{MoiseyevBook}{{1}{2011}{{Moiseyev}}{{}}}
\bibcite{BenderPTBook}{{2}{2019}{{Bender}}{{}}}
\bibcite{Burton_2019}{{3}{2019{}}{{Burton\ \emph {et~al.}}}{{Burton, Thom,\ and\ Loos}}}
\bibcite{Bittner_2012}{{4}{2012}{{Bittner\ \emph {et~al.}}}{{Bittner, Dietz, G\"unther, Harney, Miski-Oglu, Richter,\ and\ Sch\"afer}}}
\bibcite{Chong_2011}{{5}{2011}{{Chong\ \emph {et~al.}}}{{Chong, Ge,\ and\ Stone}}}
\bibcite{Chtchelkatchev_2012}{{6}{2012}{{Chtchelkatchev\ \emph {et~al.}}}{{Chtchelkatchev, Golubov, Baturina,\ and\ Vinokur}}}
\bibcite{Doppler_2016}{{7}{2016}{{Doppler\ \emph {et~al.}}}{{Doppler, Mailybaev, B\"ohm, Kuhl, Girschik, Libisch, Milburn, Rabl, Moiseyev,\ and\ Rotter}}}
\bibcite{Guo_2009}{{8}{2009}{{Guo\ \emph {et~al.}}}{{Guo, Salamo, Duchesne, Morandotti, {Volatier-Ravat}, Aimez, Siviloglou,\ and\ Christodoulides}}}
\bibcite{Hang_2013}{{9}{2013}{{Hang\ \emph {et~al.}}}{{Hang, Huang,\ and\ Konotop}}}
\bibcite{Liertzer_2012}{{10}{2012}{{Liertzer\ \emph {et~al.}}}{{Liertzer, Ge, Cerjan, Stone, T\"ureci,\ and\ Rotter}}}
\bibcite{Longhi_2010}{{11}{2010}{{Longhi}}{{}}}
\bibcite{Peng_2014}{{12}{2014{}}{{Peng\ \emph {et~al.}}}{{Peng, Ozdemir, Rotter, Yilmaz, Liertzer, Monifi, Bender, Nori,\ and\ Yang}}}
\bibcite{Peng_2014a}{{13}{2014{}}{{Peng\ \emph {et~al.}}}{{Peng, \"Ozdemir, Lei, Monifi, Gianfreda, Long, Fan, Nori, Bender,\ and\ Yang}}}
\bibcite{Regensburger_2012}{{14}{2012}{{Regensburger\ \emph {et~al.}}}{{Regensburger, Bersch, Miri, Onishchukov, Christodoulides,\ and\ Peschel}}}
\bibcite{Ruter_2010}{{15}{2010}{{R\"uter\ \emph {et~al.}}}{{R\"uter, Makris, {El-Ganainy}, Christodoulides, Segev,\ and\ Kip}}}
\bibcite{Schindler_2011}{{16}{2011}{{Schindler\ \emph {et~al.}}}{{Schindler, Li, Zheng, Ellis,\ and\ Kottos}}}
\bibcite{Szameit_2011}{{17}{2011}{{Szameit\ \emph {et~al.}}}{{Szameit, Rechtsman, {Bahat-Treidel},\ and\ Segev}}}
\bibcite{Zhao_2010}{{18}{2010}{{Zhao\ \emph {et~al.}}}{{Zhao, Schaden,\ and\ Wu}}}
\bibcite{Zheng_2013}{{19}{2013}{{Zheng\ \emph {et~al.}}}{{Zheng, Hao,\ and\ Long}}}
\bibcite{Choi_2018}{{20}{2018}{{Choi\ \emph {et~al.}}}{{Choi, Hahn, Yoon,\ and\ Song}}}
\bibcite{El-Ganainy_2018}{{21}{2018}{{El-Ganainy\ \emph {et~al.}}}{{El-Ganainy, Makris, Khajavikhan, Musslimani, Rotter,\ and\ Christodoulides}}}
\bibcite{Heiss_1990}{{22}{1990}{{Heiss\ and\ Sannino}}{{}}}
\bibcite{Heiss_1999}{{23}{1999}{{D.}}{{}}}
\bibcite{Heiss_2012}{{24}{2012}{{Heiss}}{{}}}
\bibcite{Heiss_2016}{{25}{2016}{{Heiss}}{{}}}
\bibcite{Yarkony_1996}{{26}{1996}{{Yarkony}}{{}}}
\bibcite{Benda_2018}{{27}{2018}{{Benda\ and\ Jagau}}{{}}}
\bibcite{Berry_1984}{{28}{1984}{{Berry}}{{}}}
\bibcite{BenderBook}{{29}{1978}{{Bender\ and\ Orszag}}{{}}}
\bibcite{Olsen_1996}{{30}{}{{Olsen\ \emph {et~al.}}}{{Olsen, Christiansen, Koch,\ and\ J{\o }rgensen}}}
\bibcite{Olsen_2000}{{31}{}{{Olsen\ \emph {et~al.}}}{{Olsen, J{\o }rgensen, Helgaker,\ and\ Christiansen}}}
\bibcite{Olsen_2019}{{32}{2019}{{Olsen\ and\ J{\o }rgensen}}{{}}}
\bibcite{SzaboBook}{{33}{1989}{{Szabo\ and\ Ostlund}}{{}}}
\bibcite{Goodson_2012}{{34}{2012}{{Goodson}}{{}}}
\bibcite{Moller_1934}{{35}{1934}{{M{\o }ller\ and\ Plesset}}{{}}}
\bibcite{Gill_1986}{{36}{1986}{{Gill\ and\ Radom}}{{}}}
\bibcite{Gill_1988}{{37}{1988}{{Gill\ \emph {et~al.}}}{{Gill, Pople, Radom,\ and\ Nobes}}}
\bibcite{Handy_1985}{{38}{1985}{{Handy\ \emph {et~al.}}}{{Handy, Knowles,\ and\ Somasundram}}}
\bibcite{Lepetit_1988}{{39}{1988}{{Lepetit\ \emph {et~al.}}}{{Lepetit, P{\'e}lissier,\ and\ Malrieu}}}
\bibcite{JensenBook}{{40}{2017}{{Jensen}}{{}}}
\bibcite{Seidl_2018}{{41}{2018}{{Seidl\ \emph {et~al.}}}{{Seidl, Giarrusso, Vuckovic, Fabiano,\ and\ Gori-Giorgi}}}
\bibcite{Cremer_1996}{{42}{1996}{{Cremer\ and\ He}}{{}}}
\bibcite{Christiansen_1996}{{43}{1996}{{Christiansen\ \emph {et~al.}}}{{Christiansen, Olsen, J{\o }rgensen, Koch,\ and\ Malmqvist}}}
\bibcite{Sergeev_2005}{{44}{2005}{{Sergeev\ \emph {et~al.}}}{{Sergeev, Goodson, Wheeler,\ and\ Allen}}}
\bibcite{Sergeev_2006}{{45}{2006}{{Sergeev\ and\ Goodson}}{{}}}
\bibcite{Stillinger_2000}{{46}{2000}{{Stillinger}}{{}}}
\bibcite{Baker_1971}{{47}{1971}{{{Baker}}}{{}}}
\bibcite{Goodson_2004}{{48}{2004}{{Goodson\ and\ Sergeev}}{{}}}
\bibcite{Heiss_1988}{{49}{1988}{{Heiss}}{{}}}
\bibcite{Heiss_2002}{{50}{2002}{{Heiss\ and\ M{\"u}ller}}{{}}}
\bibcite{Cejnar_2005}{{51}{2005}{{Cejnar\ \emph {et~al.}}}{{Cejnar, Heinze,\ and\ Dobe{\v s}}}}
\bibcite{Cejnar_2007}{{52}{2007}{{Cejnar\ \emph {et~al.}}}{{Cejnar, Heinze,\ and\ Macek}}}
\bibcite{Cejnar_2009}{{53}{2009}{{Cejnar\ and\ Jolie}}{{}}}
\@writefile{toc}{\contentsline {section}{\numberline {}Acknowledgments}{13}{section*.21}\protected@file@percent }
\@writefile{toc}{\contentsline {section}{\numberline {}References}{13}{section*.22}\protected@file@percent }
\bibcite{Borisov_2015}{{54}{2015}{{Borisov\ \emph {et~al.}}}{{Borisov, Ru{\v z}i{\v c}ka,\ and\ Znojil}}}
\bibcite{Sindelka_2017}{{55}{2017}{{{\v S}indelka\ \emph {et~al.}}}{{{\v S}indelka, Santos,\ and\ Moiseyev}}}
\bibcite{Sachdev_2011}{{56}{2011}{{Sachdev}}{{}}}
\bibcite{Cejnar_2015}{{57}{2015}{{Cejnar\ \emph {et~al.}}}{{Cejnar, Str{\'a}nsk{\'y},\ and\ Kloc}}}
\bibcite{Cejnar_2016}{{58}{2016}{{Cejnar\ and\ Str{\'a}nsk{\'y}}}{{}}}
\bibcite{Caprio_2008}{{59}{2008}{{Caprio\ \emph {et~al.}}}{{Caprio, Cejnar,\ and\ Iachello}}}
\bibcite{Macek_2019}{{60}{2019}{{Macek\ \emph {et~al.}}}{{Macek, Str{\'a}nsk{\'y}, Leviatan,\ and\ Cejnar}}}
\bibcite{Stransky_2018}{{61}{2018}{{Str{\'a}nsk{\'y}\ \emph {et~al.}}}{{Str{\'a}nsk{\'y}, Dvo{\v r}{\'a}k,\ and\ Cejnar}}}
\bibcite{Lipkin_1965}{{62}{1965}{{Lipkin\ \emph {et~al.}}}{{Lipkin, Meshkov,\ and\ Glick}}}
\bibcite{Ushveridze_1994}{{63}{1994}{{Ushveridze}}{{}}}
\bibcite{Thompson_2005}{{64}{2005}{{Thompson\ and\ Alavi}}{{}}}
\bibcite{Seidl_2007}{{65}{2007}{{Seidl}}{{}}}
\bibcite{Loos_2009b}{{66}{2009{}}{{Loos\ and\ Gill}}{{}}}
\bibcite{Loos_2009}{{67}{2009{}}{{Loos\ and\ Gill}}{{}}}
\bibcite{Wigner_1934}{{68}{1934}{{Wigner}}{{}}}
\bibcite{Coulson_1949}{{69}{}{{Coulson\ and\ Fischer}}{{}}}
\bibcite{AngularBook}{{70}{1996}{{Edmonds}}{{}}}
\bibcite{GiulianiBook}{{71}{2005}{{Giuliani\ and\ Vignale}}{{}}}
\bibcite{Hiscock_2014}{{72}{2014}{{Hiscock\ and\ Thom}}{{}}}
\bibcite{Burton_2019a}{{73}{2019{}}{{Burton\ \emph {et~al.}}}{{Burton, Thom,\ and\ Loos}}}
\bibcite{SlaterBook}{{74}{1960}{{Slater}}{{}}}
\bibcite{Burton_2019a}{{4}{2019{}}{{Burton\ \emph {et~al.}}}{{Burton, Thom,\ and\ Loos}}}
\bibcite{Bittner_2012}{{5}{2012}{{Bittner\ \emph {et~al.}}}{{Bittner, Dietz, G\"unther, Harney, Miski-Oglu, Richter,\ and\ Sch\"afer}}}
\bibcite{Chong_2011}{{6}{2011}{{Chong\ \emph {et~al.}}}{{Chong, Ge,\ and\ Stone}}}
\bibcite{Chtchelkatchev_2012}{{7}{2012}{{Chtchelkatchev\ \emph {et~al.}}}{{Chtchelkatchev, Golubov, Baturina,\ and\ Vinokur}}}
\bibcite{Doppler_2016}{{8}{2016}{{Doppler\ \emph {et~al.}}}{{Doppler, Mailybaev, B\"ohm, Kuhl, Girschik, Libisch, Milburn, Rabl, Moiseyev,\ and\ Rotter}}}
\bibcite{Guo_2009}{{9}{2009}{{Guo\ \emph {et~al.}}}{{Guo, Salamo, Duchesne, Morandotti, {Volatier-Ravat}, Aimez, Siviloglou,\ and\ Christodoulides}}}
\bibcite{Hang_2013}{{10}{2013}{{Hang\ \emph {et~al.}}}{{Hang, Huang,\ and\ Konotop}}}
\bibcite{Liertzer_2012}{{11}{2012}{{Liertzer\ \emph {et~al.}}}{{Liertzer, Ge, Cerjan, Stone, T\"ureci,\ and\ Rotter}}}
\bibcite{Longhi_2010}{{12}{2010}{{Longhi}}{{}}}
\bibcite{Peng_2014}{{13}{2014{}}{{Peng\ \emph {et~al.}}}{{Peng, Ozdemir, Rotter, Yilmaz, Liertzer, Monifi, Bender, Nori,\ and\ Yang}}}
\bibcite{Peng_2014a}{{14}{2014{}}{{Peng\ \emph {et~al.}}}{{Peng, \"Ozdemir, Lei, Monifi, Gianfreda, Long, Fan, Nori, Bender,\ and\ Yang}}}
\bibcite{Regensburger_2012}{{15}{2012}{{Regensburger\ \emph {et~al.}}}{{Regensburger, Bersch, Miri, Onishchukov, Christodoulides,\ and\ Peschel}}}
\bibcite{Ruter_2010}{{16}{2010}{{R\"uter\ \emph {et~al.}}}{{R\"uter, Makris, {El-Ganainy}, Christodoulides, Segev,\ and\ Kip}}}
\bibcite{Schindler_2011}{{17}{2011}{{Schindler\ \emph {et~al.}}}{{Schindler, Li, Zheng, Ellis,\ and\ Kottos}}}
\bibcite{Szameit_2011}{{18}{2011}{{Szameit\ \emph {et~al.}}}{{Szameit, Rechtsman, {Bahat-Treidel},\ and\ Segev}}}
\bibcite{Zhao_2010}{{19}{2010}{{Zhao\ \emph {et~al.}}}{{Zhao, Schaden,\ and\ Wu}}}
\bibcite{Zheng_2013}{{20}{2013}{{Zheng\ \emph {et~al.}}}{{Zheng, Hao,\ and\ Long}}}
\bibcite{Choi_2018}{{21}{2018}{{Choi\ \emph {et~al.}}}{{Choi, Hahn, Yoon,\ and\ Song}}}
\bibcite{El-Ganainy_2018}{{22}{2018}{{El-Ganainy\ \emph {et~al.}}}{{El-Ganainy, Makris, Khajavikhan, Musslimani, Rotter,\ and\ Christodoulides}}}
\bibcite{Heiss_1988}{{23}{1988}{{Heiss}}{{}}}
\bibcite{Heiss_1990}{{24}{1990}{{Heiss\ and\ Sannino}}{{}}}
\bibcite{Heiss_1999}{{25}{1999}{{Heiss}}{{}}}
\bibcite{Berry_2011}{{26}{2011}{{Berry\ and\ Uzdin}}{{}}}
\bibcite{Heiss_2012}{{27}{2012}{{Heiss}}{{}}}
\bibcite{Heiss_2016}{{28}{2016}{{Heiss}}{{}}}
\bibcite{Benda_2018}{{29}{2018}{{Benda\ and\ Jagau}}{{}}}
\bibcite{Yarkony_1996}{{30}{1996}{{Yarkony}}{{}}}
\bibcite{Berry_1984}{{31}{1984}{{Berry}}{{}}}
\bibcite{BenderBook}{{32}{1978}{{Bender\ and\ Orszag}}{{}}}
\bibcite{Olsen_1996}{{33}{1996}{{Olsen\ \emph {et~al.}}}{{Olsen, Christiansen, Koch,\ and\ J{\o }rgensen}}}
\bibcite{Olsen_2000}{{34}{2000}{{Olsen\ \emph {et~al.}}}{{Olsen, J{\o }rgensen, Helgaker,\ and\ Christiansen}}}
\bibcite{Olsen_2019}{{35}{2019}{{Olsen\ and\ J{\o }rgensen}}{{}}}
\bibcite{Mihalka_2017a}{{36}{2017}{{Mih\'alka\ and\ Surj\'an}}{{}}}
\bibcite{Mihalka_2017b}{{37}{2017}{{Mih{\'a}lka\ \emph {et~al.}}}{{Mih{\'a}lka, Szabados,\ and\ Surj{\'a}n}}}
\bibcite{Mihalka_2019}{{38}{2019}{{Mih{\'a}lka\ \emph {et~al.}}}{{Mih{\'a}lka, Szabados,\ and\ Surj{\'a}n}}}
\bibcite{SzaboBook}{{39}{1989}{{Szabo\ and\ Ostlund}}{{}}}
\bibcite{Goodson_2012}{{40}{2012}{{Goodson}}{{}}}
\bibcite{Moller_1934}{{41}{1934}{{M{\o }ller\ and\ Plesset}}{{}}}
\bibcite{Gill_1994}{{42}{1994}{{Gill}}{{}}}
\bibcite{Gill_1986}{{43}{1986}{{Gill\ and\ Radom}}{{}}}
\bibcite{Gill_1988}{{44}{1988}{{Gill\ \emph {et~al.}}}{{Gill, Pople, Radom,\ and\ Nobes}}}
\bibcite{Handy_1985}{{45}{1985}{{Handy\ \emph {et~al.}}}{{Handy, Knowles,\ and\ Somasundram}}}
\bibcite{Lepetit_1988}{{46}{1988}{{Lepetit\ \emph {et~al.}}}{{Lepetit, P{\'e}lissier,\ and\ Malrieu}}}
\bibcite{JensenBook}{{47}{2017}{{Jensen}}{{}}}
\bibcite{Nesbet_1955}{{48}{1955}{{Nesbet\ and\ Hartree}}{{}}}
\bibcite{Epstein_1926}{{49}{1926}{{Epstein}}{{}}}
\bibcite{Seidl_2018}{{50}{2018}{{Seidl\ \emph {et~al.}}}{{Seidl, Giarrusso, Vuckovic, Fabiano,\ and\ Gori-Giorgi}}}
\bibcite{Cremer_1996}{{51}{1996}{{Cremer\ and\ He}}{{}}}
\bibcite{Christiansen_1996}{{52}{1996}{{Christiansen\ \emph {et~al.}}}{{Christiansen, Olsen, J{\o }rgensen, Koch,\ and\ Malmqvist}}}
\bibcite{Sergeev_2005}{{53}{2005}{{Sergeev\ \emph {et~al.}}}{{Sergeev, Goodson, Wheeler,\ and\ Allen}}}
\bibcite{Sergeev_2006}{{54}{2006}{{Sergeev\ and\ Goodson}}{{}}}
\bibcite{Stillinger_2000}{{55}{2000}{{Stillinger}}{{}}}
\bibcite{Baker_1971}{{56}{1971}{{{Baker}}}{{}}}
\bibcite{Goodson_2004}{{57}{2004}{{Goodson\ and\ Sergeev}}{{}}}
\bibcite{Heiss_2002}{{58}{2002}{{Heiss\ and\ M{\"u}ller}}{{}}}
\bibcite{Cejnar_2005}{{59}{2005}{{Cejnar\ \emph {et~al.}}}{{Cejnar, Heinze,\ and\ Dobe{\v s}}}}
\bibcite{Cejnar_2007}{{60}{2007}{{Cejnar\ \emph {et~al.}}}{{Cejnar, Heinze,\ and\ Macek}}}
\bibcite{Cejnar_2009}{{61}{2009}{{Cejnar\ and\ Jolie}}{{}}}
\bibcite{Borisov_2015}{{62}{2015}{{Borisov\ \emph {et~al.}}}{{Borisov, Ru{\v z}i{\v c}ka,\ and\ Znojil}}}
\bibcite{Sindelka_2017}{{63}{2017}{{{\v S}indelka\ \emph {et~al.}}}{{{\v S}indelka, Santos,\ and\ Moiseyev}}}
\bibcite{Sachdev_2011}{{64}{2011}{{Sachdev}}{{}}}
\@writefile{toc}{\contentsline {section}{\numberline {}Acknowledgments}{7}{section*.16}\protected@file@percent }
\@writefile{toc}{\contentsline {section}{\numberline {}References}{7}{section*.17}\protected@file@percent }
\bibcite{Cejnar_2015}{{65}{2015}{{Cejnar\ \emph {et~al.}}}{{Cejnar, Str{\'a}nsk{\'y},\ and\ Kloc}}}
\bibcite{Cejnar_2016}{{66}{2016}{{Cejnar\ and\ Str{\'a}nsk{\'y}}}{{}}}
\bibcite{Caprio_2008}{{67}{2008}{{Caprio\ \emph {et~al.}}}{{Caprio, Cejnar,\ and\ Iachello}}}
\bibcite{Macek_2019}{{68}{2019}{{Macek\ \emph {et~al.}}}{{Macek, Str{\'a}nsk{\'y}, Leviatan,\ and\ Cejnar}}}
\bibcite{Stransky_2018}{{69}{2018}{{Str{\'a}nsk{\'y}\ \emph {et~al.}}}{{Str{\'a}nsk{\'y}, Dvo{\v r}{\'a}k,\ and\ Cejnar}}}
\bibcite{Lipkin_1965}{{70}{1965}{{Lipkin\ \emph {et~al.}}}{{Lipkin, Meshkov,\ and\ Glick}}}
\bibstyle{apsrev4-1}
\citation{REVTEX41Control}
\citation{apsrev41Control}
\newlabel{LastBibItem}{{74}{14}{}{section*.22}{}}
\newlabel{LastPage}{{}{14}{}{}{}}
\newlabel{LastBibItem}{{70}{8}{}{section*.17}{}}
\newlabel{LastPage}{{}{8}{}{}{}}

217
Manuscript/EPAWTFT.bbl
View File

@ -6,7 +6,7 @@
%Control: page (0) single
%Control: year (1) truncated
%Control: production of eprint (0) enabled
\begin{thebibliography}{74}%
\begin{thebibliography}{70}%
\makeatletter
\providecommand \@ifxundefined [1]{%
\@ifx{#1\undefined}
@ -73,6 +73,17 @@
{journal} {\bibinfo {journal} {J. Chem. Phys.}\ }\textbf {\bibinfo {volume}
{150}},\ \bibinfo {pages} {041103} (\bibinfo {year}
{2019}{\natexlab{a}})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Burton}\ \emph
{et~al.}(2019{\natexlab{b}})\citenamefont {Burton}, \citenamefont {Thom},\
and\ \citenamefont {Loos}}]{Burton_2019a}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {H.~G.~A.}\
\bibnamefont {Burton}}, \bibinfo {author} {\bibfnamefont {A.~J.~W.}\
\bibnamefont {Thom}}, \ and\ \bibinfo {author} {\bibfnamefont {P.-F.}\
\bibnamefont {Loos}},\ }\href {\doibase 10.1021/acs.jctc.9b00289} {\bibfield
{journal} {\bibinfo {journal} {J. Chem. Theory Comput.}\ }\textbf {\bibinfo
{volume} {15}},\ \bibinfo {pages} {4374} (\bibinfo {year}
{2019}{\natexlab{b}})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Bittner}\ \emph {et~al.}(2012)\citenamefont
{Bittner}, \citenamefont {Dietz}, \citenamefont {G\"unther}, \citenamefont
{Harney}, \citenamefont {Miski-Oglu}, \citenamefont {Richter},\ and\
@ -305,6 +316,13 @@
{\doibase 10.1038/nphys4323} {\bibfield {journal} {\bibinfo {journal} {Nat.
Phys.}\ }\textbf {\bibinfo {volume} {14}},\ \bibinfo {pages} {11} (\bibinfo
{year} {2018})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Heiss}(1988)}]{Heiss_1988}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {W.~D.}\ \bibnamefont
{Heiss}},\ }\href {\doibase 10.1007/BF01283767} {\bibfield {journal}
{\bibinfo {journal} {Z. Physik A - Atomic Nuclei}\ }\textbf {\bibinfo
{volume} {329}},\ \bibinfo {pages} {133} (\bibinfo {year}
{1988})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Heiss}\ and\ \citenamefont
{Sannino}(1990)}]{Heiss_1990}%
\BibitemOpen
@ -314,12 +332,21 @@
{journal} {\bibinfo {journal} {J. Phys. Math. Gen.}\ }\textbf {\bibinfo
{volume} {23}},\ \bibinfo {pages} {1167} (\bibinfo {year}
{1990})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {D.}(1999)}]{Heiss_1999}%
\bibitem [{\citenamefont {Heiss}(1999)}]{Heiss_1999}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {W.}~\bibnamefont
{D.}},\ }\href {\doibase 10.1007/s100530050339} {\bibfield {journal}
\bibfield {author} {\bibinfo {author} {\bibfnamefont {W.~D.}\ \bibnamefont
{Heiss}},\ }\href {\doibase 10.1007/s100530050339} {\bibfield {journal}
{\bibinfo {journal} {Eur. Phys. J. D}\ }\textbf {\bibinfo {volume} {7}},\
\bibinfo {pages} {1} (\bibinfo {year} {1999})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Berry}\ and\ \citenamefont
{Uzdin}(2011)}]{Berry_2011}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {M.~V.}\ \bibnamefont
{Berry}}\ and\ \bibinfo {author} {\bibfnamefont {R.}~\bibnamefont {Uzdin}},\
}\href {\doibase 10.1088/1751-8113/44/43/435303} {\bibfield {journal}
{\bibinfo {journal} {J. Phys. A Math. Theor.}\ }\textbf {\bibinfo {volume}
{44}},\ \bibinfo {pages} {435303} (\bibinfo {year} {2011})}\BibitemShut
{NoStop}%
\bibitem [{\citenamefont {Heiss}(2012)}]{Heiss_2012}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {W.~D.}\ \bibnamefont
@ -333,13 +360,6 @@
{Heiss}},\ }\href {\doibase 10.1038/nphys3864} {\bibfield {journal}
{\bibinfo {journal} {Nat. Phys.}\ }\textbf {\bibinfo {volume} {12}},\
\bibinfo {pages} {823} (\bibinfo {year} {2016})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Yarkony}(1996)}]{Yarkony_1996}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {D.~R.}\ \bibnamefont
{Yarkony}},\ }\href {\doibase 10.1103/RevModPhys.68.985} {\bibfield
{journal} {\bibinfo {journal} {Rev. Mod. Phys.}\ }\textbf {\bibinfo {volume}
{68}},\ \bibinfo {pages} {985} (\bibinfo {year} {1996})}\BibitemShut
{NoStop}%
\bibitem [{\citenamefont {Benda}\ and\ \citenamefont
{Jagau}(2018)}]{Benda_2018}%
\BibitemOpen
@ -349,6 +369,13 @@
{journal} {\bibinfo {journal} {J. Phys. Chem. Lett.}\ }\textbf {\bibinfo
{volume} {9}},\ \bibinfo {pages} {6978} (\bibinfo {year} {2018})}\BibitemShut
{NoStop}%
\bibitem [{\citenamefont {Yarkony}(1996)}]{Yarkony_1996}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {D.~R.}\ \bibnamefont
{Yarkony}},\ }\href {\doibase 10.1103/RevModPhys.68.985} {\bibfield
{journal} {\bibinfo {journal} {Rev. Mod. Phys.}\ }\textbf {\bibinfo {volume}
{68}},\ \bibinfo {pages} {985} (\bibinfo {year} {1996})}\BibitemShut
{NoStop}%
\bibitem [{\citenamefont {Berry}(1984)}]{Berry_1984}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {M.~V.}\ \bibnamefont
@ -365,9 +392,9 @@
Methods for Scientists and Engineers: Asymptotics Methods and Perturbation
Theory}}}\ (\bibinfo {publisher} {Springer},\ \bibinfo {year}
{1978})\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Olsen}\ \emph {et~al.}({\natexlab{a}})\citenamefont
{Olsen}, \citenamefont {Christiansen}, \citenamefont {Koch},\ and\
\citenamefont {J{\o}rgensen}}]{Olsen_1996}%
\bibitem [{\citenamefont {Olsen}\ \emph {et~al.}(1996)\citenamefont {Olsen},
\citenamefont {Christiansen}, \citenamefont {Koch},\ and\ \citenamefont
{J{\o}rgensen}}]{Olsen_1996}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {J.}~\bibnamefont
{Olsen}}, \bibinfo {author} {\bibfnamefont {O.}~\bibnamefont {Christiansen}},
@ -375,10 +402,10 @@
{author} {\bibfnamefont {P.}~\bibnamefont {J{\o}rgensen}},\ }\href {\doibase
10.1063/1.472352} {\bibfield {journal} {\bibinfo {journal} {J. Chem.
Phys.}\ }\textbf {\bibinfo {volume} {105}},\ \bibinfo {pages} {5082}
({\natexlab{a}})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Olsen}\ \emph {et~al.}({\natexlab{b}})\citenamefont
{Olsen}, \citenamefont {J{\o}rgensen}, \citenamefont {Helgaker},\ and\
\citenamefont {Christiansen}}]{Olsen_2000}%
(\bibinfo {year} {1996})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Olsen}\ \emph {et~al.}(2000)\citenamefont {Olsen},
\citenamefont {J{\o}rgensen}, \citenamefont {Helgaker},\ and\ \citenamefont
{Christiansen}}]{Olsen_2000}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {J.}~\bibnamefont
{Olsen}}, \bibinfo {author} {\bibfnamefont {P.}~\bibnamefont {J{\o}rgensen}},
@ -386,7 +413,7 @@
\bibinfo {author} {\bibfnamefont {O.}~\bibnamefont {Christiansen}},\ }\href
{\doibase 10.1063/1.481611} {\bibfield {journal} {\bibinfo {journal} {J.
Chem. Phys.}\ }\textbf {\bibinfo {volume} {112}},\ \bibinfo {pages} {9736}
({\natexlab{b}})}\BibitemShut {NoStop}%
(\bibinfo {year} {2000})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Olsen}\ and\ \citenamefont
{J{\o}rgensen}(2019)}]{Olsen_2019}%
\BibitemOpen
@ -395,6 +422,37 @@
{J{\o}rgensen}},\ }\href {\doibase 10.1063/1.5110554} {\bibfield {journal}
{\bibinfo {journal} {J. Chem. Phys.}\ }\textbf {\bibinfo {volume} {151}},\
\bibinfo {pages} {084108} (\bibinfo {year} {2019})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Mih\'alka}\ and\ \citenamefont
{Surj\'an}(2017)}]{Mihalka_2017a}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {Z.~E.}\ \bibnamefont
{Mih\'alka}}\ and\ \bibinfo {author} {\bibfnamefont {P.~R.}\ \bibnamefont
{Surj\'an}},\ }\href {\doibase 10.1103/PhysRevA.96.062106} {\bibfield
{journal} {\bibinfo {journal} {Phys. Rev. A}\ }\textbf {\bibinfo {volume}
{96}},\ \bibinfo {pages} {062106} (\bibinfo {year} {2017})}\BibitemShut
{NoStop}%
\bibitem [{\citenamefont {Mih{\'a}lka}\ \emph {et~al.}(2017)\citenamefont
{Mih{\'a}lka}, \citenamefont {Szabados},\ and\ \citenamefont
{Surj{\'a}n}}]{Mihalka_2017b}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {Z.~{\'E}.}\
\bibnamefont {Mih{\'a}lka}}, \bibinfo {author} {\bibfnamefont
{{\'A}.}~\bibnamefont {Szabados}}, \ and\ \bibinfo {author} {\bibfnamefont
{P.~R.}\ \bibnamefont {Surj{\'a}n}},\ }\href {\doibase 10.1063/1.4978898}
{\bibfield {journal} {\bibinfo {journal} {J. Chem. Phys.}\ }\textbf
{\bibinfo {volume} {146}},\ \bibinfo {pages} {124121} (\bibinfo {year}
{2017})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Mih{\'a}lka}\ \emph {et~al.}(2019)\citenamefont
{Mih{\'a}lka}, \citenamefont {Szabados},\ and\ \citenamefont
{Surj{\'a}n}}]{Mihalka_2019}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {Z.~{\'E}.}\
\bibnamefont {Mih{\'a}lka}}, \bibinfo {author} {\bibfnamefont
{{\'A}.}~\bibnamefont {Szabados}}, \ and\ \bibinfo {author} {\bibfnamefont
{P.~R.}\ \bibnamefont {Surj{\'a}n}},\ }\href {\doibase 10.1063/1.5083191}
{\bibfield {journal} {\bibinfo {journal} {J. Chem. Phys.}\ }\textbf
{\bibinfo {volume} {150}},\ \bibinfo {pages} {031101} (\bibinfo {year}
{2019})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Szabo}\ and\ \citenamefont
{Ostlund}(1989)}]{SzaboBook}%
\BibitemOpen
@ -417,6 +475,13 @@
{Plesset}},\ }\href {\doibase 10.1103/PhysRev.46.618} {\bibfield {journal}
{\bibinfo {journal} {Phys. Rev.}\ }\textbf {\bibinfo {volume} {46}},\
\bibinfo {pages} {618} (\bibinfo {year} {1934})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Gill}(1994)}]{Gill_1994}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {P.~M.~W.}\
\bibnamefont {Gill}},\ }\href {\doibase 10.1016/S0065-3276(08)60019-2}
{\bibfield {journal} {\bibinfo {journal} {Adv. Quantum Chem.}\ }\textbf
{\bibinfo {volume} {25}},\ \bibinfo {pages} {141} (\bibinfo {year}
{1994})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Gill}\ and\ \citenamefont {Radom}(1986)}]{Gill_1986}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {P.~M.~W.}\
@ -462,6 +527,21 @@
{Jensen}},\ }\href@noop {} {\emph {\bibinfo {title} {Introduction to
computational chemistry}}}\ (\bibinfo {publisher} {Wiley},\ \bibinfo {year}
{2017})\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Nesbet}\ and\ \citenamefont
{Hartree}(1955)}]{Nesbet_1955}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {R.~K.}\ \bibnamefont
{Nesbet}}\ and\ \bibinfo {author} {\bibfnamefont {D.~R.}\ \bibnamefont
{Hartree}},\ }\href {\doibase 10.1098/rspa.1955.0134} {\bibfield {journal}
{\bibinfo {journal} {Proc. R. Soc. Lond. A}\ }\textbf {\bibinfo {volume}
{230}},\ \bibinfo {pages} {312} (\bibinfo {year} {1955})}\BibitemShut
{NoStop}%
\bibitem [{\citenamefont {Epstein}(1926)}]{Epstein_1926}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {P.~S.}\ \bibnamefont
{Epstein}},\ }\href {\doibase 10.1103/PhysRev.28.695} {\bibfield {journal}
{\bibinfo {journal} {Phys. Rev.}\ }\textbf {\bibinfo {volume} {28}},\
\bibinfo {pages} {695} (\bibinfo {year} {1926})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Seidl}\ \emph {et~al.}(2018)\citenamefont {Seidl},
\citenamefont {Giarrusso}, \citenamefont {Vuckovic}, \citenamefont
{Fabiano},\ and\ \citenamefont {Gori-Giorgi}}]{Seidl_2018}%
@ -535,13 +615,6 @@
{\bibinfo {booktitle} {Adv. Quantum Chem.}}},\ Vol.~\bibinfo {volume} {47}\
(\bibinfo {publisher} {Academic Press},\ \bibinfo {year} {2004})\ pp.\
\bibinfo {pages} {193--208}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Heiss}(1988)}]{Heiss_1988}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {W.~D.}\ \bibnamefont
{Heiss}},\ }\href {\doibase 10.1007/BF01283767} {\bibfield {journal}
{\bibinfo {journal} {Z. Physik A - Atomic Nuclei}\ }\textbf {\bibinfo
{volume} {329}},\ \bibinfo {pages} {133} (\bibinfo {year}
{1988})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Heiss}\ and\ \citenamefont
{M{\"u}ller}(2002)}]{Heiss_2002}%
\BibitemOpen
@ -664,96 +737,4 @@
{\doibase 10.1016/0029-5582(65)90862-X} {\bibfield {journal} {\bibinfo
{journal} {Nucl. Phys.}\ }\textbf {\bibinfo {volume} {62}},\ \bibinfo {pages}
{188} (\bibinfo {year} {1965})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Ushveridze}(1994)}]{Ushveridze_1994}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {A.~G.}\ \bibnamefont
{Ushveridze}},\ }\href@noop {} {\emph {\bibinfo {title} {Quasi-exactly
solvable models in quantum mechanics}}}\ (\bibinfo {publisher} {Institute of
Physics Publishing},\ \bibinfo {year} {1994})\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Thompson}\ and\ \citenamefont
{Alavi}(2005)}]{Thompson_2005}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {D.~C.}\ \bibnamefont
{Thompson}}\ and\ \bibinfo {author} {\bibfnamefont {A.}~\bibnamefont
{Alavi}},\ }\href {\doibase 10.1063/1.1869978} {\bibfield {journal}
{\bibinfo {journal} {J. Chem. Phys.}\ }\textbf {\bibinfo {volume} {122}},\
\bibinfo {pages} {124107} (\bibinfo {year} {2005})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Seidl}(2007)}]{Seidl_2007}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {M.}~\bibnamefont
{Seidl}},\ }\href {\doibase 10.1103/PhysRevA.75.062506} {\bibfield {journal}
{\bibinfo {journal} {Phys. Rev. A}\ }\textbf {\bibinfo {volume} {75}},\
\bibinfo {pages} {062506} (\bibinfo {year} {2007})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Loos}\ and\ \citenamefont
{Gill}(2009{\natexlab{a}})}]{Loos_2009b}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {P.-F.}\ \bibnamefont
{Loos}}\ and\ \bibinfo {author} {\bibfnamefont {P.~M.~W.}\ \bibnamefont
{Gill}},\ }\href {\doibase 10.1103/PhysRevLett.103.123008} {\bibfield
{journal} {\bibinfo {journal} {Phys. Rev. Let.}\ }\textbf {\bibinfo {volume}
{103}},\ \bibinfo {pages} {123008} (\bibinfo {year}
{2009}{\natexlab{a}})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Loos}\ and\ \citenamefont
{Gill}(2009{\natexlab{b}})}]{Loos_2009}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {P.-F.}\ \bibnamefont
{Loos}}\ and\ \bibinfo {author} {\bibfnamefont {P.~M.~W.}\ \bibnamefont
{Gill}},\ }\href {\doibase 10.1103/PhysRevA.79.062517} {\bibfield {journal}
{\bibinfo {journal} {Phys. Rev. A}\ }\textbf {\bibinfo {volume} {79}},\
\bibinfo {pages} {062517} (\bibinfo {year} {2009}{\natexlab{b}})}\BibitemShut
{NoStop}%
\bibitem [{\citenamefont {Wigner}(1934)}]{Wigner_1934}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {E.}~\bibnamefont
{Wigner}},\ }\href {\doibase 10.1103/PhysRev.46.1002} {\bibfield {journal}
{\bibinfo {journal} {Phys. Rev.}\ }\textbf {\bibinfo {volume} {46}},\
\bibinfo {pages} {1002} (\bibinfo {year} {1934})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Coulson}\ and\ \citenamefont
{Fischer}()}]{Coulson_1949}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {C.~A.}\ \bibnamefont
{Coulson}}\ and\ \bibinfo {author} {\bibfnamefont {I.}~\bibnamefont
{Fischer}},\ }\href {\doibase 10.1080/14786444908521726} {\bibfield
{journal} {\bibinfo {journal} {1949}\ }\textbf {\bibinfo {volume} {40}},\
\bibinfo {pages} {386}}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Edmonds}(1996)}]{AngularBook}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {A.~R.}\ \bibnamefont
{Edmonds}},\ }\href@noop {} {\emph {\bibinfo {title} {Angular {Momentum} in
{Quantum} {Mechanics}}}}\ (\bibinfo {publisher} {Princeton University
Press},\ \bibinfo {year} {1996})\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Giuliani}\ and\ \citenamefont
{Vignale}(2005)}]{GiulianiBook}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {G.}~\bibnamefont
{Giuliani}}\ and\ \bibinfo {author} {\bibfnamefont {G.}~\bibnamefont
{Vignale}},\ }\href {\doibase 10.1017/CBO9780511619915} {\emph {\bibinfo
{title} {Quantum {Theory} of the {Electron} {Liquid}}}}\ (\bibinfo
{publisher} {Cambridge University Press},\ \bibinfo {year}
{2005})\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Hiscock}\ and\ \citenamefont
{Thom}(2014)}]{Hiscock_2014}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {H.~G.}\ \bibnamefont
{Hiscock}}\ and\ \bibinfo {author} {\bibfnamefont {A.~J.~W.}\ \bibnamefont
{Thom}},\ }\href {\doibase 10.1021/ct5007696} {\bibfield {journal} {\bibinfo
{journal} {J. Chem. Theory Comput.}\ }\textbf {\bibinfo {volume} {10}},\
\bibinfo {pages} {4795} (\bibinfo {year} {2014})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Burton}\ \emph
{et~al.}(2019{\natexlab{b}})\citenamefont {Burton}, \citenamefont {Thom},\
and\ \citenamefont {Loos}}]{Burton_2019a}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {H.~G.~A.}\
\bibnamefont {Burton}}, \bibinfo {author} {\bibfnamefont {A.~J.~W.}\
\bibnamefont {Thom}}, \ and\ \bibinfo {author} {\bibfnamefont {P.-F.}\
\bibnamefont {Loos}},\ }\href {\doibase 10.1021/acs.jctc.9b00289} {\bibfield
{journal} {\bibinfo {journal} {J. Chem. Theory Comput.}\ }\textbf {\bibinfo
{volume} {15}},\ \bibinfo {pages} {4374} (\bibinfo {year}
{2019}{\natexlab{b}})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Slater}(1960)}]{SlaterBook}%
\BibitemOpen
\bibfield {author} {\bibinfo {author} {\bibfnamefont {J.~C.}\ \bibnamefont
{Slater}},\ }\href@noop {} {\emph {\bibinfo {title} {{Quantum Theory of
Atomic Structure}}}}\ (\bibinfo {publisher} {McGraw-Hill},\ \bibinfo {year}
{1960})\BibitemShut {NoStop}%
\end{thebibliography}%

1536
Manuscript/EPAWTFT.bib
View File

File diff suppressed because it is too large Load Diff

71
Manuscript/EPAWTFT.blg
View File

@ -23,46 +23,45 @@ Control: production of article title (-1) disabled
Control: page (0) single
Control: year (1) truncated
Control: production of eprint (0) enabled
Warning--empty year in Coulson_1949
You've used 76 entries,
You've used 72 entries,
5847 wiz_defined-function locations,
2147 strings with 28164 characters,
and the built_in function-call counts, 72588 in all, are:
= -- 4650
> -- 1882
< -- 494
+ -- 601
- -- 459
* -- 11148
:= -- 7202
add.period$ -- 74
call.type$ -- 76
change.case$ -- 298
chr.to.int$ -- 80
cite$ -- 77
duplicate$ -- 6786
empty$ -- 5460
format.name$ -- 1102
if$ -- 14470
int.to.chr$ -- 5
int.to.str$ -- 83
missing$ -- 887
newline$ -- 274
num.names$ -- 222
pop$ -- 2827
2133 strings with 28024 characters,
and the built_in function-call counts, 70033 in all, are:
= -- 4473
> -- 1850
< -- 481
+ -- 585
- -- 453
* -- 10828
:= -- 6899
add.period$ -- 70
call.type$ -- 72
change.case$ -- 282
chr.to.int$ -- 78
cite$ -- 72
duplicate$ -- 6517
empty$ -- 5216
format.name$ -- 1078
if$ -- 13975
int.to.chr$ -- 3
int.to.str$ -- 79
missing$ -- 860
newline$ -- 262
num.names$ -- 210
pop$ -- 2686
preamble$ -- 1
purify$ -- 367
purify$ -- 350
quote$ -- 0
skip$ -- 2666
skip$ -- 2570
stack$ -- 0
substring$ -- 1918
swap$ -- 6353
text.length$ -- 222
substring$ -- 1863
swap$ -- 6176
text.length$ -- 219
text.prefix$ -- 0
top$ -- 10
type$ -- 1036
warning$ -- 2
while$ -- 223
type$ -- 993
warning$ -- 1
while$ -- 215
width$ -- 0
write$ -- 633
(There were 2 warnings)
write$ -- 606
(There was 1 warning)

263
Manuscript/EPAWTFT.log
View File

@ -1,4 +1,4 @@
This is pdfTeX, Version 3.14159265-2.6-1.40.21 (TeX Live 2020) (preloaded format=pdflatex 2020.6.14) 24 AUG 2020 11:00
This is pdfTeX, Version 3.14159265-2.6-1.40.21 (TeX Live 2020) (preloaded format=pdflatex 2020.6.14) 12 NOV 2020 12:04
entering extended mode
restricted \write18 enabled.
file:line:error style messages enabled.
@ -587,55 +587,55 @@ Package uniquecounter Info: New unique counter `rerunfilecheck' on input line 2
\Hy@SectionHShift=\skip78
)
Class revtex4-1 Info: Unrecognized array package. Please update this document c
lass! (Proceeding with fingers crossed.) on input line 73.
lass! (Proceeding with fingers crossed.) on input line 74.
Class revtex4-1 Info: Patching unrecognized longtable package. (Proceeding with
fingers crossed) on input line 73.
Class revtex4-1 Info: cite was not loaded (OK!) on input line 73.
Class revtex4-1 Info: mcite was not loaded (OK!) on input line 73.
Class revtex4-1 Info: multicol was not loaded (OK!) on input line 73.
fingers crossed) on input line 74.
Class revtex4-1 Info: cite was not loaded (OK!) on input line 74.
Class revtex4-1 Info: mcite was not loaded (OK!) on input line 74.
Class revtex4-1 Info: multicol was not loaded (OK!) on input line 74.
(./EPAWTFT.aux)
\openout1 = `EPAWTFT.aux'.
LaTeX Font Info: Checking defaults for OML/txmi/m/it on input line 73.
LaTeX Font Info: Checking defaults for OML/txmi/m/it on input line 74.
LaTeX Font Info: Trying to load font information for OML+txmi on input line
73.
74.
(/usr/local/texlive/2020/texmf-dist/tex/latex/txfonts/omltxmi.fd
File: omltxmi.fd 2000/12/15 v3.1
)
LaTeX Font Info: ... okay on input line 73.
LaTeX Font Info: Checking defaults for OMS/txsy/m/n on input line 73.
LaTeX Font Info: ... okay on input line 74.
LaTeX Font Info: Checking defaults for OMS/txsy/m/n on input line 74.
LaTeX Font Info: Trying to load font information for OMS+txsy on input line
73.
74.
(/usr/local/texlive/2020/texmf-dist/tex/latex/txfonts/omstxsy.fd
File: omstxsy.fd 2000/12/15 v3.1
)
LaTeX Font Info: ... okay on input line 73.
LaTeX Font Info: Checking defaults for OT1/cmr/m/n on input line 73.
LaTeX Font Info: ... okay on input line 73.
LaTeX Font Info: Checking defaults for T1/cmr/m/n on input line 73.
LaTeX Font Info: ... okay on input line 73.
LaTeX Font Info: Checking defaults for TS1/cmr/m/n on input line 73.
LaTeX Font Info: ... okay on input line 73.
LaTeX Font Info: Checking defaults for OMX/txex/m/n on input line 73.
LaTeX Font Info: ... okay on input line 74.
LaTeX Font Info: Checking defaults for OT1/cmr/m/n on input line 74.
LaTeX Font Info: ... okay on input line 74.
LaTeX Font Info: Checking defaults for T1/cmr/m/n on input line 74.
LaTeX Font Info: ... okay on input line 74.
LaTeX Font Info: Checking defaults for TS1/cmr/m/n on input line 74.
LaTeX Font Info: ... okay on input line 74.
LaTeX Font Info: Checking defaults for OMX/txex/m/n on input line 74.
LaTeX Font Info: Trying to load font information for OMX+txex on input line
73.
74.
(/usr/local/texlive/2020/texmf-dist/tex/latex/txfonts/omxtxex.fd
File: omxtxex.fd 2000/12/15 v3.1
)
LaTeX Font Info: ... okay on input line 73.
LaTeX Font Info: Checking defaults for U/txexa/m/n on input line 73.
LaTeX Font Info: ... okay on input line 74.
LaTeX Font Info: Checking defaults for U/txexa/m/n on input line 74.
LaTeX Font Info: Trying to load font information for U+txexa on input line 7
3.
4.
(/usr/local/texlive/2020/texmf-dist/tex/latex/txfonts/utxexa.fd
File: utxexa.fd 2000/12/15 v3.1
)
LaTeX Font Info: ... okay on input line 73.
LaTeX Font Info: Checking defaults for PD1/pdf/m/n on input line 73.
LaTeX Font Info: ... okay on input line 73.
LaTeX Font Info: ... okay on input line 74.
LaTeX Font Info: Checking defaults for PD1/pdf/m/n on input line 74.
LaTeX Font Info: ... okay on input line 74.
(/usr/local/texlive/2020/texmf-dist/tex/context/base/mkii/supp-pdf.mkii
[Loading MPS to PDF converter (version 2006.09.02).]
@ -659,14 +659,14 @@ Package epstopdf-base Info: Redefining graphics rule for `.eps' on input line 4
File: epstopdf-sys.cfg 2010/07/13 v1.3 Configuration of (r)epstopdf for TeX Liv
e
))
LaTeX Info: Redefining \microtypecontext on input line 73.
LaTeX Info: Redefining \microtypecontext on input line 74.
Package microtype Info: Generating PDF output.
Package microtype Info: Character protrusion enabled (level 2).
Package microtype Info: Using default protrusion set `alltext'.
Package microtype Info: Automatic font expansion enabled (level 2),
(microtype) stretch: 20, shrink: 20, step: 1, non-selected.
Package microtype Info: Using default expansion set `basictext'.
LaTeX Info: Redefining \showhyphens on input line 73.
LaTeX Info: Redefining \showhyphens on input line 74.
Package microtype Info: No adjustment of tracking.
Package microtype Info: No adjustment of interword spacing.
Package microtype Info: No adjustment of character kerning.
@ -675,10 +675,10 @@ Package microtype Info: No adjustment of character kerning.
File: mt-ptm.cfg 2006/04/20 v1.7 microtype config. file: Times (RS)
)
Package chemgreek Info: Activating mapping `default' on line 73.
Package chemgreek Info: Activating mapping `default' on line 74.
\AtBeginShipoutBox=\box73
Package hyperref Info: Link coloring ON on input line 73.
Package hyperref Info: Link coloring ON on input line 74.
(/usr/local/texlive/2020/texmf-dist/tex/latex/hyperref/nameref.sty
Package: nameref 2019/09/16 v2.46 Cross-referencing by name of section
@ -692,20 +692,20 @@ Package: gettitlestring 2019/12/15 v1.6 Cleanup title references (HO)
)
\c@section@level=\count298
)
LaTeX Info: Redefining \ref on input line 73.
LaTeX Info: Redefining \pageref on input line 73.
LaTeX Info: Redefining \nameref on input line 73.
LaTeX Info: Redefining \ref on input line 74.
LaTeX Info: Redefining \pageref on input line 74.
LaTeX Info: Redefining \nameref on input line 74.
(./EPAWTFT.out) (./EPAWTFT.out)
\@outlinefile=\write4
\openout4 = `EPAWTFT.out'.
Class revtex4-1 Info: Taking \maketitle back from hyperref on input line 73.
Class revtex4-1 Info: Taking \maketitle back from hyperref on input line 74.
\openout3 = `EPAWTFTNotes.bib'.
LaTeX Font Info: Font shape `T1/txr/b/n' in size <12> not available
(Font) Font shape `T1/txr/bx/n' tried instead on input line 93.
(Font) Font shape `T1/txr/bx/n' tried instead on input line 95.
LaTeX Font Info: Trying to load font information for OT1+txr on input line 9
3.
5.
(/usr/local/texlive/2020/texmf-dist/tex/latex/txfonts/ot1txr.fd
File: ot1txr.fd 2000/12/15 v3.1
@ -715,194 +715,69 @@ File: mt-cmr.cfg 2013/05/19 v2.2 microtype config. file: Computer Modern Roman
(RS)
)
LaTeX Font Info: Trying to load font information for U+txsya on input line 9
3.
5.
(/usr/local/texlive/2020/texmf-dist/tex/latex/txfonts/utxsya.fd
File: utxsya.fd 2000/12/15 v3.1
)
LaTeX Font Info: Trying to load font information for U+txsyb on input line 9
3.
5.
(/usr/local/texlive/2020/texmf-dist/tex/latex/txfonts/utxsyb.fd
File: utxsyb.fd 2000/12/15 v3.1
)
LaTeX Font Info: Trying to load font information for U+txmia on input line 9
3.
5.
(/usr/local/texlive/2020/texmf-dist/tex/latex/txfonts/utxmia.fd
File: utxmia.fd 2000/12/15 v3.1
)
LaTeX Font Info: Trying to load font information for U+txsyc on input line 9
3.
5.
(/usr/local/texlive/2020/texmf-dist/tex/latex/txfonts/utxsyc.fd
File: utxsyc.fd 2000/12/15 v3.1
)
LaTeX Font Info: Font shape `T1/txr/b/n' in size <9> not available
(Font) Font shape `T1/txr/bx/n' tried instead on input line 96.
(Font) Font shape `T1/txr/bx/n' tried instead on input line 98.
[1{/usr/local/texlive/2020/texmf-var/fonts/map/pdftex/updmap/pdftex.map}] [2]
[3]
LaTeX Font Info: Trying to load font information for TS1+txr on input line 3
52.
(/usr/local/texlive/2020/texmf-dist/tex/latex/txfonts/ts1txr.fd
File: ts1txr.fd 2000/12/15 v3.1
) [4] [5]
[6] [7]
<EsbHF.pdf, id=300, 578.16pt x 373.395pt>
File: EsbHF.pdf Graphic file (type pdf)
<use EsbHF.pdf>
Package pdftex.def Info: EsbHF.pdf used on input line 566.
(pdftex.def) Requested size: 246.0pt x 158.87024pt.
[8 <./EsbHF.pdf>]
<EMP_RHF_R10.pdf, id=341, 520.94624pt x 512.91624pt>
File: EMP_RHF_R10.pdf Graphic file (type pdf)
<use EMP_RHF_R10.pdf>
Package pdftex.def Info: EMP_RHF_R10.pdf used on input line 607.
(pdftex.def) Requested size: 246.0pt x 242.2055pt.
<PartitioningRCV2.pdf, id=349, 578.16pt x 380.42125pt>
File: PartitioningRCV2.pdf Graphic file (type pdf)
<use PartitioningRCV2.pdf>
Package pdftex.def Info: PartitioningRCV2.pdf used on input line 624.
(pdftex.def) Requested size: 249.9028pt x 164.43123pt.
<PartitioningRCV3.pdf, id=350, 578.16pt x 380.42125pt>
File: PartitioningRCV3.pdf Graphic file (type pdf)
<use PartitioningRCV3.pdf>
Package pdftex.def Info: PartitioningRCV3.pdf used on input line 625.
(pdftex.def) Requested size: 249.9028pt x 164.43123pt.
Overfull \hbox (3.9028pt too wide) in paragraph at lines 624--626
[][]
[]
Overfull \hbox (3.9028pt too wide) in paragraph at lines 624--626
[]
[]
[9 <./EMP_RHF_R10.pdf>]
LaTeX Warning: A float is stuck (cannot be placed); try class option [floatfix]
on input line 636.
<MPlargebasis.pdf, id=383, 578.16pt x 388.45125pt>
File: MPlargebasis.pdf Graphic file (type pdf)
<use MPlargebasis.pdf>
Package pdftex.def Info: MPlargebasis.pdf used on input line 638.
(pdftex.def) Requested size: 249.9028pt x 167.90208pt.
<WCElargebasis.pdf, id=384, 578.16pt x 374.39874pt>
File: WCElargebasis.pdf Graphic file (type pdf)
<use WCElargebasis.pdf>
Package pdftex.def Info: WCElargebasis.pdf used on input line 639.
(pdftex.def) Requested size: 249.9028pt x 161.82811pt.
Overfull \hbox (3.9028pt too wide) in paragraph at lines 638--640
[][]
[]
Overfull \hbox (3.9028pt too wide) in paragraph at lines 638--640
[]
[]
Overfull \hbox (115.91017pt too wide) in paragraph at lines 647--658
[]$[]$
[]
[10 <./PartitioningRCV2.pdf> <./PartitioningRCV3.pdf
pdfTeX warning: /Library/TeX/texbin/pdflatex (file ./PartitioningRCV3.pdf): PDF
inclusion: multiple pdfs with page group included in a single page
> <./MPlargebasis.pdf
pdfTeX warning: /Library/TeX/texbin/pdflatex (file ./MPlargebasis.pdf): PDF inc
lusion: multiple pdfs with page group included in a single page
> <./WCElargebasis.pdf
pdfTeX warning: /Library/TeX/texbin/pdflatex (file ./WCElargebasis.pdf): PDF in
clusion: multiple pdfs with page group included in a single page
>]
<EMP_UHF_R10.pdf, id=435, 520.94624pt x 512.91624pt>
File: EMP_UHF_R10.pdf Graphic file (type pdf)
<use EMP_UHF_R10.pdf>
Package pdftex.def Info: EMP_UHF_R10.pdf used on input line 703.
(pdftex.def) Requested size: 246.0pt x 242.2055pt.
<UHFCI.pdf, id=441, 520.94624pt x 501.875pt>
File: UHFCI.pdf Graphic file (type pdf)
<use UHFCI.pdf>
Package pdftex.def Info: UHFCI.pdf used on input line 711.
(pdftex.def) Requested size: 229.49844pt x 221.0937pt.
<UHFEP.pdf, id=442, 520.94624pt x 501.875pt>
File: UHFEP.pdf Graphic file (type pdf)
<use UHFEP.pdf>
Package pdftex.def Info: UHFEP.pdf used on input line 712.
(pdftex.def) Requested size: 229.49844pt x 221.0937pt.
[11 <./EMP_UHF_R10.pdf> <./UHFCI.pdf
pdfTeX warning: /Library/TeX/texbin/pdflatex (file ./UHFCI.pdf): PDF inclusion:
multiple pdfs with page group included in a single page
> <./UHFEP.pdf
pdfTeX warning: /Library/TeX/texbin/pdflatex (file ./UHFEP.pdf): PDF inclusion:
multiple pdfs with page group included in a single page
>]
<ReNRJPT.pdf, id=496, 520.94624pt x 527.9725pt>
File: ReNRJPT.pdf Graphic file (type pdf)
<use ReNRJPT.pdf>
Package pdftex.def Info: ReNRJPT.pdf used on input line 724.
(pdftex.def) Requested size: 229.49844pt x 232.59058pt.
<ImNRJPT.pdf, id=497, 520.94624pt x 527.9725pt>
File: ImNRJPT.pdf Graphic file (type pdf)
<use ImNRJPT.pdf>
Package pdftex.def Info: ImNRJPT.pdf used on input line 725.
(pdftex.def) Requested size: 229.49844pt x 232.59058pt.
LaTeX Warning: `!h' float specifier changed to `!ht'.
[12 <./ReNRJPT.pdf> <./ImNRJPT.pdf
pdfTeX warning: /Library/TeX/texbin/pdflatex (file ./ImNRJPT.pdf): PDF inclusio
n: multiple pdfs with page group included in a single page
>] (./EPAWTFT.bbl [13])
[3] [4]
[5] [6] (./EPAWTFT.bbl [7])
Class revtex4-1 Warning: Environment acknowledgements is obsolete;
Use acknowledgments instead.
[14]
Package atveryend Info: Empty hook `AfterLastShipout' on input line 768.
[8]
Package atveryend Info: Empty hook `AfterLastShipout' on input line 744.
(./EPAWTFT.aux)
Package atveryend Info: Executing hook `AtVeryEndDocument' on input line 768.
Package atveryend Info: Executing hook `AtEndAfterFileList' on input line 768.
Package atveryend Info: Executing hook `AtVeryEndDocument' on input line 744.
Package atveryend Info: Executing hook `AtEndAfterFileList' on input line 744.
Package rerunfilecheck Info: File `EPAWTFT.out' has not changed.
(rerunfilecheck) Checksum: 92DE4D91FCC62EBAF1E5ED2DCAC4EAE9;1594.
(rerunfilecheck) Checksum: 81A90960DB67E307FACD10337C344D5D;1233.
)
Here is how much of TeX's memory you used:
14631 strings out of 480478
221331 string characters out of 5899595
570775 words of memory out of 5000000
29659 multiletter control sequences out of 15000+600000
631037 words of font info for 549 fonts, out of 8000000 for 9000
14313 strings out of 480478
216976 string characters out of 5899595
574099 words of memory out of 5000000
29499 multiletter control sequences out of 15000+600000
618445 words of font info for 469 fonts, out of 8000000 for 9000
1141 hyphenation exceptions out of 8191
44i,18n,58p,1927b,602s stack positions out of 5000i,500n,10000p,200000b,80000s
44i,18n,58p,1757b,626s stack positions out of 5000i,500n,10000p,200000b,80000s
{/usr/local/texlive/2020/texmf-dist/fonts/enc/dvips/base/8r.enc}</usr/local/t
exlive/2020/texmf-dist/fonts/type1/public/txfonts/rtcxr.pfb></usr/local/texlive
/2020/texmf-dist/fonts/type1/public/txfonts/rtxi.pfb></usr/local/texlive/2020/t
exmf-dist/fonts/type1/public/txfonts/rtxmi.pfb></usr/local/texlive/2020/texmf-d
ist/fonts/type1/public/txfonts/rtxr.pfb></usr/local/texlive/2020/texmf-dist/fon
ts/type1/public/txfonts/txbsy.pfb></usr/local/texlive/2020/texmf-dist/fonts/typ
e1/public/txfonts/txex.pfb></usr/local/texlive/2020/texmf-dist/fonts/type1/publ
ic/txfonts/txexa.pfb></usr/local/texlive/2020/texmf-dist/fonts/type1/public/txf
onts/txsy.pfb></usr/local/texlive/2020/texmf-dist/fonts/type1/public/txfonts/tx
sya.pfb></usr/local/texlive/2020/texmf-dist/fonts/type1/public/txfonts/txsyb.pf
b></usr/local/texlive/2020/texmf-dist/fonts/type1/urw/times/utmb8a.pfb></usr/lo
cal/texlive/2020/texmf-dist/fonts/type1/urw/times/utmr8a.pfb></usr/local/texliv
e/2020/texmf-dist/fonts/type1/urw/times/utmri8a.pfb>
Output written on EPAWTFT.pdf (14 pages, 849408 bytes).
exlive/2020/texmf-dist/fonts/type1/public/txfonts/rtxi.pfb></usr/local/texlive/
2020/texmf-dist/fonts/type1/public/txfonts/rtxmi.pfb></usr/local/texlive/2020/t
exmf-dist/fonts/type1/public/txfonts/rtxr.pfb></usr/local/texlive/2020/texmf-di
st/fonts/type1/public/txfonts/txbsy.pfb></usr/local/texlive/2020/texmf-dist/fon
ts/type1/public/txfonts/txex.pfb></usr/local/texlive/2020/texmf-dist/fonts/type
1/public/txfonts/txexa.pfb></usr/local/texlive/2020/texmf-dist/fonts/type1/publ
ic/txfonts/txsy.pfb></usr/local/texlive/2020/texmf-dist/fonts/type1/public/txfo
nts/txsyb.pfb></usr/local/texlive/2020/texmf-dist/fonts/type1/urw/times/utmb8a.
pfb></usr/local/texlive/2020/texmf-dist/fonts/type1/urw/times/utmr8a.pfb></usr/
local/texlive/2020/texmf-dist/fonts/type1/urw/times/utmri8a.pfb>
Output written on EPAWTFT.pdf (8 pages, 134017 bytes).
PDF statistics:
743 PDF objects out of 1000 (max. 8388607)
656 compressed objects within 7 object streams
170 named destinations out of 1000 (max. 500000)
59112 words of extra memory for PDF output out of 61914 (max. 10000000)
457 PDF objects out of 1000 (max. 8388607)
430 compressed objects within 5 object streams
120 named destinations out of 1000 (max. 500000)
47241 words of extra memory for PDF output out of 51595 (max. 10000000)

23
Manuscript/EPAWTFT.out
View File

@ -1,4 +1,4 @@
\BOOKMARK [0][-]{section*.2}{Perturbation theory in the complex plane}{}% 2
\BOOKMARK [0][-]{section*.2}{Perturbation theory in the complex plane: Exceptional points and where to find them}{}% 2
\BOOKMARK [1][-]{section*.1}{Abstract}{section*.2}% 1
\BOOKMARK [1][-]{section*.3}{Introduction}{section*.2}% 3
\BOOKMARK [2][-]{section*.4}{Background}{section*.3}% 4
@ -7,16 +7,11 @@
\BOOKMARK [2][-]{section*.7}{Rayleigh-Schr\366dinger perturbation theory}{section*.6}% 7
\BOOKMARK [2][-]{section*.8}{The Hartree-Fock Hamiltonian}{section*.6}% 8
\BOOKMARK [2][-]{section*.9}{M\370ller-Plesset perturbation theory}{section*.6}% 9
\BOOKMARK [2][-]{section*.10}{Alternative partitioning}{section*.6}% 10
\BOOKMARK [1][-]{section*.11}{Historical overview}{section*.2}% 11
\BOOKMARK [2][-]{section*.12}{Behavior of the M\370ller-Plesset series}{section*.11}% 12
\BOOKMARK [2][-]{section*.13}{Cases of divergence}{section*.11}% 13
\BOOKMARK [2][-]{section*.14}{The singularity structure}{section*.11}% 14
\BOOKMARK [2][-]{section*.15}{The physics of quantum phase transition}{section*.11}% 15
\BOOKMARK [1][-]{section*.16}{The spherium model}{section*.2}% 16
\BOOKMARK [1][-]{section*.17}{Radius of convergence and exceptional points}{section*.2}% 17
\BOOKMARK [2][-]{section*.18}{Evolution of the radius of convergence}{section*.17}% 18
\BOOKMARK [2][-]{section*.19}{Exceptional points in the UHF formalism}{section*.17}% 19
\BOOKMARK [1][-]{section*.20}{Conclusion}{section*.2}% 20
\BOOKMARK [1][-]{section*.21}{Acknowledgments}{section*.2}% 21
\BOOKMARK [1][-]{section*.22}{References}{section*.2}% 22
\BOOKMARK [1][-]{section*.10}{Historical overview}{section*.2}% 10
\BOOKMARK [2][-]{section*.11}{Behavior of the M\370ller-Plesset series}{section*.10}% 11
\BOOKMARK [2][-]{section*.12}{Cases of divergence}{section*.10}% 12
\BOOKMARK [2][-]{section*.13}{The singularity structure}{section*.10}% 13
\BOOKMARK [2][-]{section*.14}{The physics of quantum phase transitions}{section*.10}% 14
\BOOKMARK [1][-]{section*.15}{Conclusion}{section*.2}% 15
\BOOKMARK [1][-]{section*.16}{Acknowledgments}{section*.2}% 16
\BOOKMARK [1][-]{section*.17}{References}{section*.2}% 17

BIN
Manuscript/EPAWTFT.pdf
View File

Binary file not shown.

BIN
Manuscript/EPAWTFT.synctex.gz
View File

Binary file not shown.

746
Manuscript/EPAWTFT.tex
View File

@ -31,6 +31,7 @@
% operators
\newcommand{\bH}{\mathbf{H}}
\newcommand{\bV}{\mathbf{V}}
\newcommand{\bh}{\mathbf{h}}
\newcommand{\bQ}{\mathbf{Q}}
\newcommand{\bSig}{\mathbf{\Sigma}}
@ -72,32 +73,35 @@
\newcommand{\VU}{Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, FEW, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands.}
\begin{document}
\title{Perturbation theory in the complex plane}
%\title{Exceptional points and where to find them}
\title{Perturbation theory in the complex plane: Exceptional points and where to find them}
\author{Antoine \surname{Marie}}
\affiliation{\LCPQ}
\author{Hugh G.~A.~\surname{Burton}}
\affiliation{\UCAM}
\author{Paola \surname{Gori-Giorgi}}
\affiliation{\VU}
%\author{Paola \surname{Gori-Giorgi}}
%\affiliation{\VU}
\author{Pierre-Fran\c{c}ois \surname{Loos}}
\email{loos@irsamc.ups-tlse.fr}
\affiliation{\LCPQ}
\begin{abstract}
In this work, we explore the extension of quantum chemistry in the complex plane. We observe that the physics of a quantum system is intimately connected to the position of the energy singularities in the complex plane. After a brief presentation of the fundamental notions of quantum chemistry and perturbation theory in the complex plane, we provide a historical overview of the various research activities that have been performed on the physic of singularities. Then we connect and further discuss these different aspects using the spherium model (\ie, two opposite-spin electrons restricted to remain on the surface of a sphere) as a theoretical playground. In particular, we explore various perturbative partitioning strategies and the effects of symmetry breaking on the singularity structure of the electronic energy. This provides fundamental insights on the location of these singularities in the complex plane (that one calls exceptional points) and, specifically, on the magnitude of the radius of convergence associated with the perturbative treatment.
In this review, we explore the extension of quantum chemistry in the complex plane.
We observe that the physics of a quantum system is intimately connected to the position of the energy singularities in the complex plane.
After a presentation of the fundamental notions of quantum chemistry and perturbation theory in the complex plane, we provide a historical overview of the various research activities that have been performed on the physic of singularities.
\end{abstract}
\maketitle
%============================================================%
%%%%%%%%%%%%%%%%%%%%%%%
\section{Introduction}
\label{sec:intro}
%============================================================%
%%%%%%%%%%%%%%%%%%%%%%%
%======================%
\subsection{Background}
%=======================%
Due to the ubiquitous influence of processes involving electronic excited states in physics, chemistry, and biology, their faithful description from first principles has been one of the grand challenges faced by theoretical chemists since the dawn of computational chemistry.
Accurately predicting ground- and excited-state energies (hence excitation energies) is particularly valuable in this context, and it has concentrated most of the efforts within the community.
@ -111,24 +115,26 @@ However, one can gain a different perspective on quantization extending quantum
In a non-Hermitian complex picture, the energy levels are \textit{sheets} of a more complicated topological manifold called \textit{Riemann surface}, and they are smooth and continuous \textit{analytic continuation} of one another.
In other words, our view of the quantized nature of conventional Hermitian quantum mechanics arises only from our limited perception of the more complex and profound structure of its non-Hermitian variant. \cite{MoiseyevBook,BenderPTBook}
The realization that ground and excited states both emerge from one single mathematical structure with equal importance suggests that excited-state energies can be computed from first principles in their own right.
One could then exploit the structure of these Riemann surfaces to develop methods that directly target excited-state energies without needing ground-state information. \cite{Burton_2019}
One could then exploit the structure of these Riemann surfaces to develop methods that directly target excited-state energies without needing ground-state information. \cite{Burton_2019,Burton_2019a}
By analytically continuing the electronic energy $E(\lambda)$ in the complex domain (where $\lambda$ is a coupling parameter), the ground and excited states of a molecule can be smoothly connected.
This connection is possible because by extending real numbers to the complex domain, the ordering property of real numbers is lost.
Hence, electronic states can be interchanged away from the real axis since the concept of ground and excited states has been lost.
Amazingly, this smooth and continuous transition from one state to another has recently been experimentally realized in physical settings such as electronics, microwaves, mechanics, acoustics, atomic systems and optics. \cite{Bittner_2012, Chong_2011, Chtchelkatchev_2012, Doppler_2016, Guo_2009, Hang_2013, Liertzer_2012, Longhi_2010, Peng_2014, Peng_2014a, Regensburger_2012, Ruter_2010, Schindler_2011, Szameit_2011, Zhao_2010, Zheng_2013, Choi_2018, El-Ganainy_2018}
Amazingly, this smooth and continuous transition from one state to another has recently been experimentally realized in physical settings such as electronics, microwaves, mechanics, acoustics, atomic systems and optics. \cite{Bittner_2012,Chong_2011,Chtchelkatchev_2012,Doppler_2016,Guo_2009,Hang_2013,Liertzer_2012,Longhi_2010,Peng_2014, Peng_2014a,Regensburger_2012,Ruter_2010,Schindler_2011,Szameit_2011,Zhao_2010,Zheng_2013,Choi_2018,El-Ganainy_2018}
Exceptional points (EPs) are branch point singularities where two (or more) states become exactly degenerate. \cite{Heiss_1990,Heiss_1999,Heiss_2012,Heiss_2016}
Exceptional points (EPs) are branch point singularities where two (or more) states become exactly degenerate. \cite{MoiseyevBook,Heiss_1988,Heiss_1990,Heiss_1999,Berry_2011,Heiss_2012,Heiss_2016,Benda_2018}
They are the non-Hermitian analogs of conical intersections, \cite{Yarkony_1996} which are ubiquitous in non-adiabatic processes and play a key role in photo-chemical mechanisms.
In the case of auto-ionizing resonances, EPs have a role in deactivation processes similar to conical intersections in the decay of bound excited states. \cite{Benda_2018}
Although Hermitian and non-Hermitian Hamiltonians are closely related, the behavior of their eigenvalues near degeneracies is starkly different.
For example, encircling non-Hermitian degeneracies at EPs leads to an interconversion of states, and two loops around the EP are necessary to recover the initial energy.
For example, encircling non-Hermitian degeneracies at EPs leads to an interconversion of states, and two loops around the EP are necessary to recover the initial energy. \cite{MoiseyevBook,Heiss_2016,Benda_2018}
Additionally, the wave function picks up a geometric phase (also known as Berry phase \cite{Berry_1984}) and four loops are required to recover the initial wave function.
In contrast, encircling Hermitian degeneracies at conical intersections only introduces a geometric phase while leaving the states unchanged.
More dramatically, whilst eigenvectors remain orthogonal at conical intersections, at non-Hermitian EPs the eigenvectors themselves become equivalent, resulting in a \textit{self-orthogonal} state. \cite{MoiseyevBook}
More importantly here, although EPs usually lie off the real axis, these singular points are intimately related to the convergence properties of perturbative methods and avoided crossing on the real axis are indicative of singularities in the complex plane. \cite{BenderBook,Olsen_1996,Olsen_2000,Olsen_2019}
More importantly here, although EPs usually lie off the real axis, these singular points are intimately related to the convergence properties of perturbative methods and avoided crossing on the real axis are indicative of singularities in the complex plane. \cite{BenderBook,Olsen_1996,Olsen_2000,Olsen_2019,Mihalka_2017a,Mihalka_2017b,Mihalka_2019}
%===================================%
\subsection{An illustrative example}
%===================================%
In order to highlight the general properties of EPs mentioned above, we propose to consider the following $2 \times 2$ Hamiltonian commonly used in quantum chemistry
\begin{equation}
@ -180,53 +186,16 @@ and we have
E_{\pm}(4\pi) & = E_{\pm}(0). \notag
\end{align}
This evidences that encircling non-Hermitian degeneracies at EPs leads to an interconversion of states, and two loops around the EP are necessary to recover the initial energy.
Additionally, the wave function picks up a geometric phase (also known as Berry phase \cite{Berry_1984}) and four loops are required to recover the starting wave function.
Additionally, the wave function picks up a geometric phase and four loops are required to recover the starting wave function.
\cite{MoiseyevBook}
%The eigenvectors associated to the eigenenergies \eqref{eq:E_2x2} are
%\begin{equation}\label{eq:phi_2x2}
%\begin{split}
% \phi_{\pm}
% & =
% \begin{pmatrix}
% (\epsilon_1 - \epsilon_2 \pm \sqrt{(\epsilon_1 - \epsilon_2)^2 + 4\lambda^2})/2\lambda
% \\
% 1
% \end{pmatrix}
% \\
% =
% \begin{pmatrix}
% (\epsilon_1-\epsilon_2)/2\lambda \pm \sqrt{2\lambda_\text{EP}} \sqrt{\lambda - \lambda_\text{EP}}/\lambda
% \\
% 1
% \end{pmatrix},
%\end{split}
%\end{equation}
%and, for $\lambda=\lambda_\text{EP}$, they become
%\begin{equation}
% \phi_{\pm} = \begin{pmatrix} \mp i \\ 1\end{pmatrix},
%\end{equation}
%which are clearly self-orthogonal, \ie, their norm is equal to zero.
%
%As branch point (square-root) singularities, EPs inherit their topology from the multi-valued function $\sqrt{\lambda - \lambda_\text{EP}}$.
%Then, if the eigenvectors are properly normalized, they behave as $(\lambda - \lambda_\text{EP})^{-1/4}$, which gives the following pattern when looping around one EP:
%\begin{align}
% \phi_{\pm}(2\pi) & = \phi_{\mp}(0),
% &
% \phi_{\pm}(4\pi) & = -\phi_{\pm}(0), \notag
% \\
% \phi_{\pm}(6\pi) & = -\phi_{\mp}(0),
% &
% \phi_{\pm}(8\pi) & = \phi_{\pm}(0).\notag
%\end{align}
%In plain words, four loops around the EP are necessary to recover the initial state.
%We can also see that looping the other way around leads to a different pattern.
%============================================================%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Perturbation theory}
%============================================================%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%============================================================%
\subsection{Rayleigh-Schr\"odinger perturbation theory}
%============================================================%
Within the Born-Oppenheimer approximation,
\begin{equation}\label{eq:ExactHamiltonian}
@ -246,7 +215,7 @@ is recast as
\end{equation}
where $\hH^{(0)}$ is the zeroth-order Hamiltonian and $\hV = \hH - \hH^{(0)}$ is the so-called perturbation.
The ``physical'' system of interest is recovered by setting the coupling parameter $\lambda$ to unity.
This decomposition is obviously non-unique and motivated by several factors as discussed below.
This decomposition is obviously non-unique and motivated by several factors as discussed below. \cite{Mihalka_2017b}
Accordingly to Eq.~\eqref{eq:SchrEq-PT}, the energy can then be written as a power series of $\lambda$
\begin{equation} \label{eq:Elambda}
@ -265,8 +234,9 @@ This theorem means that the radius of convergence of the perturbation series is
\end{equation}
This function is smooth for $x \in \mathbb{R}$ and infinitely differentiable in $\mathbb{R}$. One would expect that the Taylor series of such a function would be convergent $\forall x \in \mathbb{R}$. However this series is divergent for $x \ge 1$. This is because the function has four singularities in the complex plane ($x = e^{i\pi/4}$, $e^{-i\pi/4}$, $e^{i3\pi/4}$, and $e^{-i3\pi/4}$) with a modulus equal to $1$. This simple yet powerful example emphasizes the importance of the singularities in the complex plane to understand the convergence properties on the real axis.
%============================================================%
\subsection{The Hartree-Fock Hamiltonian}
%============================================================%
In the Hartree-Fock (HF) approximation, the many-electron wave function is approximated as a single Slater determinant $\Psi^{\text{HF}}(\vb{x}_1,\ldots,\vb{x}_n)$ [where $\vb{x} = (\sigma,\vb{r})$ is a composite vector gathering spin and spatial coordinates] defined as an antisymmetric combination of $n$ (real-valued) one-electron spin-orbitals $\phi_p(\vb{x})$, which are, by definition, eigenfunctions of the one-electron Fock operator
\begin{equation}\label{eq:FockOp}
@ -291,14 +261,13 @@ is the HF mean-field potential with
\end{gather}
\end{subequations}
being the Coulomb and exchange operators (respectively) in the spin-orbital basis. \cite{SzaboBook}
If the spatial part of the spin-orbitals are restricted to be the same for spin-up and spin-down electrons, one talks about restricted HF (RHF) theory, whereas if one does not enforce this constrain it leads to the so-called unrestricted HF (UHF) theory.
If the spatial part of the spin-orbitals are restricted to be the same for spin-up and spin-down electrons, one talks about restricted HF (RHF) theory, whereas if one does not enforce this constrain it leads to the so-called unrestricted HF (UHF) theory.
From hereon, $i$ and $j$ are occupied orbitals, $a$ and $b$ are unoccupied (or virtual) orbitals, while $p$, $q$, $r$, and $s$ indicate arbitrary orbitals.
Rather than solving Eq.~\eqref{eq:SchrEq}, HF theory uses the variational principle to find an approximation of $\Psi$ as a single Slater determinant. Hence a Slater determinant is not an eigenfunction of the exact Hamiltonian $\hH$. However, it is, by definition, an eigenfunction of the so-called (approximated) HF many-electron Hamiltonian defined as the sum of the one-electron Fock operators
\begin{equation}\label{eq:HFHamiltonian}
\hH^{\text{HF}} = \sum_{i} f(\vb{x}_i).
\end{equation}
%Note that the lowest eigenvalue of the Hamiltonian \eqref{eq:HFHamiltonian} is not the HF energy but the sum of the eigenvalues associated with the occupied eigenvalues (see below).
\subsection{M{\o}ller-Plesset perturbation theory}
@ -329,9 +298,9 @@ The MP2 energy reads
\end{equation}
with $\mel{pq}{}{rs} = \braket{pq}{rs} - \braket{pq}{sr}$, and where
\begin{equation}
\braket{pq}{rs} = \iint \dd\vb{x}_1\dd\vb{x}_2\phi_p(\vb{x}_1)\phi_q(\vb{x}_2)\frac{1}{\abs{\vb{r}_1 - \vb{r}_2}}\phi_r(\vb{x}_1)\phi_s(\vb{x}_2)
\braket{pq}{rs} = \iint \dd\vb{x}_1\dd\vb{x}_2\frac{\phi_p(\vb{x}_1)\phi_q(\vb{x}_2)\phi_r(\vb{x}_1)\phi_s(\vb{x}_2)}{\abs{\vb{r}_1 - \vb{r}_2}}
\end{equation}
is a two-electron integral in the spin-orbital basis.
are two-electron integral in the spin-orbital basis. \cite{Gill_1994}
As mentioned earlier, there is, \textit{a priori}, no guarantee that the MP$m$ series converges to the exact energy when $m \to \infty$.
In fact, it is known that when the HF wave function is a poor approximation to the exact wave function, for example in multi-reference systems, the MP method yields deceptive results. \cite{Gill_1986, Gill_1988, Handy_1985, Lepetit_1988}
@ -344,22 +313,19 @@ In order to better understand the behavior of the MP series and how it is connec
For small systems, one can access the whole terms of the series using full configuration interaction (FCI).
If the Hamiltonian $H(\lambda)$ is diagonalized in the FCI space, one gets the exact energies (in this finite Hilbert space) and the Taylor expansion with respect to $\lambda$ allows to access the MP perturbation series at any order.
\subsection{Alternative partitioning}\label{sec:AlterPart}
Obviously, although practically convenient for electronic structure calculations, the MP partitioning is not the only possibility, and alternative partitioning have been proposed in the literature:
i) the Epstein-Nesbet (EN) partitioning which consists in taking the diagonal elements of $\hH$ as the zeroth-order Hamiltonian. \cite{Nesbet_1955,Epstein_1926}
Hence, the off-diagonal elements of $\hH$ are the perturbation operator,
ii) the weak correlation partitioning in which the one-electron part is consider as the unperturbed Hamiltonian $\hH^{(0)}$ and the two-electron part is the perturbation operator $\hV$, and
iii) the strong coupling partitioning where the two operators are inverted compared to the weak correlation partitioning. \cite{Seidl_2018}
Obviously, although practically convenient for electronic structure calculations, the MP partitioning is not the only possibility.
Here, we will consider three alternative partitioning schemes:
\begin{itemize}
\item The Epstein-Nesbet (EN) partitioning which consists in taking the diagonal elements of $\hH$ as the zeroth-order Hamiltonian.
Hence, the off-diagonal elements of $\hH$ are the perturbation operator.
\item The weak correlation (WC) partitioning in which the one-electron part is consider as the unperturbed Hamiltonian $\hH^{(0)}$ and the two-electron part is the perturbation operator $\hV$.
\item The strong coupling (SC) partitioning where the two operators are inverted compared to the WC partitioning. \cite{Seidl_2018}
\end{itemize}
%============================================================%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Historical overview}
%============================================================%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%=====================================================%
\subsection{Behavior of the M{\o}ller-Plesset series}
%=====================================================%
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, each time a higher-order term is computed, one would like to obtained an overall result closer to the exact energy.
@ -369,23 +335,23 @@ In Ref.~\onlinecite{Gill_1986}, the authors showed that the RMP series is conver
On the other hand, the UMP series is monotonically convergent (except for the first few orders) but very slowly.
Thus, one cannot practically use it for systems where only the first terms can be computed.
When a bond is stretched, in most cases the exact wave function becomes more and more multi-reference. Yet the HF wave function is restricted to be a single Slater determinant.
Thus, it is inappropriate to model (even qualitatively) stretched systems. Nevertheless, the HF wave function can undergo a symmetry breaking to lower its energy by sacrificing one of the symmetry of the exact wave function during the process (see for example the case of \ce{H2} in Ref.~\onlinecite{SzaboBook}).
One could then potentially claim that the RMP series exhibits deceptive convergence properties as the RHF Slater determinant is a poor approximation of the exact wave function for stretched system. However, even in the unrestricted formalism which clearly represents a better description of a stretched system, the UMP series does not have the smooth and rapidly converging behavior that one would wish for.
When a bond is stretched, in most cases the exact wave function becomes more and more of multi-reference nature.
Yet the HF wave function is restricted to be a single Slater determinant.
It is then inappropriate to model (even qualitatively) stretched systems. Nevertheless, the HF wave function can undergo symmetry breaking to lower its energy by sacrificing one of the symmetry of the exact wave function in the process (see for example the case of \ce{H2} in Ref.~\onlinecite{SzaboBook}).
One could then potentially claim that the RMP series exhibits deceptive convergence properties as the RHF Slater determinant is a poor approximation of the exact wave function for stretched system. However, even in the unrestricted formalism which clearly represents a better description of a stretched system, the UMP series does not have the smooth and rapidly convergent behavior that one would wish for.
In the unrestricted framework the singlet ground state wave function is allowed to mix with triplet wave functions, leading to the so-called spin contamination issue. Gill \textit{et al.}~highlighted the link between slow convergence of the UMP series and spin contamination for \ce{H2} in a minimal basis. \cite{Gill_1988}
Handy and coworkers reported the same behavior of the series (oscillatory and slowly monotonically convergent) in stretched \ce{H2O} and \ce{NH2} systems. \cite{Handy_1985} Lepetit \textit`{et al.}~analyzed the difference between the MP and EN partitioning for the UHF reference. \cite{Lepetit_1988}
Handy and coworkers reported the same behavior of the series (oscillatory and slowly monotonically convergent) in stretched \ce{H2O} and \ce{NH2}. \cite{Handy_1985} Lepetit \textit{et al.}~analyzed the difference between the MP and EN partitioning for the UHF reference. \cite{Lepetit_1988}
They concluded that the slow convergence is due to the coupling of the singly- and doubly-excited configurations.
Cremer and He analyzed 29 atomic and molecular systems at the FCI level \cite{Cremer_1996} and grouped them in two classes: i) the \textit{class A} systems where one observes a monotonic convergence to the FCI energy, and ii) the \textit{class B} systems for which convergence is erratic after initial oscillations. Their system set contains stretched molecules as well as molecules at their equilibrium geometry for various basis sets. They highlighted that \cite{Cremer_1996}
\begin{quote}
\textit{``Class A systems are characterized by electronic structures with well-separated electron pairs while class B systems are characterized by electronic structures with electron clustering in one or more regions.''}
\end{quote}
\textit{``Class A systems are characterized by electronic structures with well-separated electron pairs while class B systems are characterized by electronic structures with electron clustering in one or more regions.''}
As one can only compute the first terms of the MP series, a smart way of getting more accurate results is to use extrapolation formula, \ie, estimating the limit of the series with only few terms.
Cremer and He proved that using specific extrapolation formulas of the MP series for class A and class B systems improves the precision of the results compared to the formula used without classes. The mean deviation from FCI correlation energies is $0.3$ millihartree with the adapted formula whereas with the formula that do not distinguish the systems it is 12 millihartree.
Even if there were still shaded areas and that this classification was incomplete, this work showed that understanding the origin of the different modes of convergence would lead to a more rationalized use of the MP perturbation theory and to more accurate correlation energies.
Cremer and He proved that using specific extrapolation formulas of the MP series for class A and class B systems improves the precision of the results compared to the formula used without resorting to classes. The mean absolute deviation taking the FCI correlation energies as reference is $0.3$ millihartree with the class-specific formula whereas the deviation increases to 12 millihartree using the general formula.
Even if there were still shaded areas in their analysis and that their classification was incomplete, the work of Cremer and He clearly evidenced that understanding the origin of the different modes of convergence could potentially lead to a more rationalized use of MP perturbation theory and, hence, to more accurate correlation energy estimates.
%===============================%
\subsection{Cases of divergence}
%===============================%
In the late 90's, Olsen \textit{et al.}~discovered an even more preoccupying behavior of the MP series. \cite{Olsen_1996} They showed that the series could be divergent even in systems that they considered as well understood like \ce{Ne} and \ce{HF}. \cite{Olsen_1996, Christiansen_1996} Cremer and He had already studied these two systems and classified them as \textit{class B} systems. However, the analysis of Olsen and coworkers was performed in larger basis sets containing diffuse functions. In these basis sets, they found that the series become divergent at (very) high order.
@ -394,24 +360,24 @@ The discovery of this divergent behavior is worrying as in order to get meaningf
\textit{``In the limit of large order, the series coefficients become equivalent to the Taylor series coefficients of the singularity closest to the origin. Following the result of this theorem, the convergence patterns of the MP series can be explained by looking at the dominant singularity.''}
\end{quote}
A singularity in the unit circle is designated as an intruder state, more precisely as a front-door (respectively back-door) intruder state if the real part of the singularity is positive (respectively negative). Their method consists in performing a scan of the real axis to detect the avoided crossing responsible for the pair of dominant singularities. Then, by modeling this avoided crossing via a two-state Hamiltonian one can get an approximation of the dominant conjugate pair of singularities by finding the EPs of the following $2\times2$ Hamiltonian
A singularity in the unit circle is designated as an intruder state, more precisely as a front-door (respectively back-door) intruder state if the real part of the singularity is positive (respectively negative). Their method consists in performing a scan of the real axis to detect the avoided crossing responsible for the pair of dominant singularities. Then, by modeling this avoided crossing via a two-state Hamiltonian one can get an approximation of the dominant conjugate pair of singularities by finding the EPs of the following $2\times2$ matrix
\begin{equation}
\underbrace{\mqty(\alpha & \delta \\ \delta & \beta)}_{\hH} = \underbrace{\mqty(\alpha + \alpha_s & 0 \\ 0 & \beta + \beta_s )}_{\hH^{(0)}} + \underbrace{\mqty(- \alpha_s & \delta \\ \delta & - \beta_s)}_{\hV},
\underbrace{\mqty(\alpha & \delta \\ \delta & \beta)}_{\bH} = \underbrace{\mqty(\alpha + \alpha_s & 0 \\ 0 & \beta + \beta_s )}_{\bH^{(0)}} + \underbrace{\mqty(- \alpha_s & \delta \\ \delta & - \beta_s)}_{\bV},
\end{equation}
where the diagonal matrix is the unperturbed Hamiltonian and the other matrix is the perturbation.
where the diagonal matrix is the unperturbed Hamiltonian matrix $\bH^{(0)}$ and the second matrix in the right-hand-side $\bV$ is the perturbation.
See also Ref.~\onlinecite{Olsen_2019} where the present model is generalized to a non-symmetric (i.e, non-Hermitian) Hamiltonian.
They first studied molecules with low-lying doubly-excited states of the same spatial and spin symmetry.
% \titou{because in those systems the HF wave function is a bad approximation.}
The exact wave function has a non-negligible contribution from the doubly-excited states, so these low-lying excited states were good candidates for being intruder states. For \ce{CH_2} in a large basis set, the series is convergent up to the 50th order. They showed that the dominant singularity lies outside the unit circle but close to it causing the slow convergence.
The exact wave function has a non-negligible contribution from the doubly-excited states, so these low-lying excited states were good candidates for being intruder states. \titou{For \ce{CH_2} in a large basis set, the series is convergent up to the 50th order. They showed that the dominant singularity lies outside the unit circle but close to it causing the slow convergence.}
Then they demonstrated that the divergence for \ce{Ne} is due to a back-door intruder state. When the basis set is augmented with diffuse functions, the ground state undergo sharp avoided crossings with highly diffuse excited states leading to a back-door intruder state. They used their two-state model on this avoided crossings and the model was actually predicting the divergence of the series.
%They concluded that the divergence of the series was due to the interaction with a highly diffuse excited state.
Moreover they proved that the extrapolation formulas of Cremer and He \cite{Cremer_1996} cannot be used for all systems, and that these formulas were not mathematically motivated when looking at the singularity causing the divergence. For example, the hydrogen fluoride molecule contains both back-door intruder states and low-lying doubly-excited states which results in alternated terms up to 10th order. For higher orders, the series is monotonically convergent. This surprising behavior is due to the fact that two pairs of singularities are approximately at the same distance from the origin.
%=======================================
\subsection{The singularity structure}
%=======================================
In the 2000's, Sergeev and Goodson \cite{Sergeev_2005, Sergeev_2006} analyzed this problem from a more mathematical point of view by looking at the whole singularity structure where Olsen and collaborators were trying to find the dominant singularity causing the divergence. They regrouped singularities in two classes: i) $\alpha$ singularities which have ``large'' imaginary parts, and ii) $\beta$ singularities which have very small imaginary parts. Singularities of type $\alpha$ are related to large avoided crossing between the ground and low-lying excited states, whereas $\beta$ singularities come from a sharp avoided crossing between the ground state and a highly diffuse state. They succeeded to explain the divergence of the series caused by $\beta$ singularities following previous work of Stillinger. \cite{Stillinger_2000}
To understand the convergence properties of the perturbation series at $\lambda=1$, one must look at the whole complex plane, in particular, for negative (\ie, real) values of $\lambda$. If $\lambda$ is negative, the Coulomb interaction becomes attractive but the mean field (which has been computed at $\lambda = 1$) remains repulsive as it is proportional to $(1-\lambda)$:
@ -421,311 +387,321 @@ To understand the convergence properties of the perturbation series at $\lambda=
\hH(\lambda) =
\sum_{i}^{n} \Bigg[
\overbrace{-\frac{1}{2}\grad_i^2
- \sum_{A}^{N} \frac{Z_A}{\abs{\vb{r}_i-\vb{R}_A}}}^{\text{Independent of $\lambda$}}
- \sum_{A}^{N} \frac{Z_A}{\abs{\vb{r}_i-\vb{R}_A}}}^{\text{independent of $\lambda$}}
\\
+ \underbrace{(1-\lambda)v^{\text{HF}}(\vb{x}_i)}_{\text{Repulsive for $\lambda < 1$}}
+ \underbrace{\lambda\sum_{i<j}^{n}\frac{1}{|\vb{r}_i-\vb{r}_j|}}_{\text{Attractive for $\lambda < 0$}}
+ \underbrace{(1-\lambda)v^{\text{HF}}(\vb{x}_i)}_{\text{repulsive for $\lambda < 1$}}
+ \underbrace{\lambda\sum_{i<j}^{n}\frac{1}{|\vb{r}_i-\vb{r}_j|}}_{\text{attractive for $\lambda < 0$}}
\Bigg].
\end{multline}
The major difference between those two terms is that the repulsive mean field is localized around the nuclei whereas the interelectronic interaction persist away from the nuclei. If $\lambda$ becomes more and more negative the mean field becomes more and more repulsive so there exists a critical (negative) value of $\lambda$, $\lambda_\text{c}$, for which the Coulombic field created by the nuclei cannot bind the electrons anymore because of the $\lambda$-independent nature of the the electron-nucleus attraction. For $\lambda = \lambda_c$, the electrons dissociate from the nuclei and form a bound cluster which is infinitely separated from the nuclei. According to Baker, \cite{Baker_1971} this value is a critical point of the system and, by analogy with thermodynamics, the energy $E(\lambda)$ exhibits a singularity at $\lambda_\text{c}$. At this point the system undergo a phase transition and a symmetry breaking. Beyond $\lambda_c$ there is a continuum of eigenstates with electrons dissociated from the nucleus.
The major difference between these two terms is that the repulsive mean field is localized around the nuclei whereas the interelectronic interaction persist away from the nuclei. If $\lambda$ becomes more and more negative the mean field becomes more and more repulsive so there exists a critical (negative) value of $\lambda$, $\lambda_\text{c}$, for which the Coulombic field created by the nuclei cannot bind the electrons anymore because of the $\lambda$-independent nature of the the electron-nucleus attraction. For $\lambda = \lambda_c$, the electrons dissociate from the nuclei and form a bound cluster which is infinitely separated from the nuclei. According to Baker, \cite{Baker_1971} this value is a critical point of the system and, by analogy with thermodynamics, the energy $E(\lambda)$ exhibits a singularity at $\lambda_\text{c}$. At this point the system undergo a phase transition and a symmetry breaking.
Beyond $\lambda_c$ there is a continuum of eigenstates thanks to which the electrons dissociated from the nuclei.
This reasoning is done on the exact Hamiltonian and energy, \ie, the Hamiltonian in the complete Hilbert space, this is the exact energy which exhibits this singularity on the negative real axis. However, in a finite basis set which does not span the complete Hilbert space, one can prove that, for a Hermitian Hamiltonian, the singularities of $E(\lambda)$ occurs in complex conjugate pairs with non-zero imaginary parts. Sergeev and Goodson proved, \cite{Sergeev_2005} as predicted by Stillinger, \cite{Stillinger_2000} that in a finite basis set the critical point on the real axis is modeled by a cluster of sharp avoided crossings with diffuse functions, equivalently by a cluster of $\beta$ singularities in the negative half plane. This explains that Olsen \textit{et al.}, because they used a simple $2\times2$ model, only observed the first singularity of this cluster of singularities causing the divergence. \cite{Olsen_2000}
Finally, it was shown that $\beta$ singularities are very sensitive to changes in the basis set but not to the stretching of the system. On the contrary, $\alpha$ singularities are relatively insensitive to the basis sets but very sensitive to bond stretching. According to Goodson, \cite{Goodson_2004} the singularity structure of stretched molecules is difficult because there is more than one significant singularity. This is consistent with the observation of Olsen and coworkers \cite{Olsen_2000} on the \ce{HF} molecule at equilibrium geometry and stretched geometry. To the best our knowledge, the effect of bond stretching on singularities, its link with spin contamination and symmetry breaking of the wave function has not been as well understood as the ionization phenomenon and its link with diffuse functions. In this work, we shall try to improve our understanding of the effect of symmetry breaking on the singularities of $E(\lambda)$ and we hope that it will lead to a deeper understanding of perturbation theory.
Finally, it was shown that $\beta$ singularities are very sensitive to changes in the basis set but not to the stretching of the system. On the contrary, $\alpha$ singularities are relatively insensitive to the basis sets but very sensitive to bond stretching.
According to Goodson, \cite{Goodson_2004} the singularity structure of stretched molecules is difficult because there is more than one significant singularity.
This is consistent with the observation of Olsen and coworkers \cite{Olsen_2000} on the \ce{HF} molecule at equilibrium geometry and stretched geometry.
To the best of our knowledge, the effect of bond stretching on singularities, its link with spin contamination and symmetry breaking of the wave function has not been as well understood as the ionization phenomenon and its link with diffuse functions.
\subsection{The physics of quantum phase transition}
%====================================================
\subsection{The physics of quantum phase transitions}
%====================================================
In the previous section, we saw that a careful analysis of the structure of the Hamiltonian allows us to predict the existence of a critical point. In a finite basis set this critical point is model by a cluster of $\beta$ singularities. It is now well known that this phenomenon is a special case of a more general phenomenon. Indeed, theoretical physicists proved that EPs close to the real axis are connected to \textit{quantum phase transitions} (QPTs). \cite{Heiss_1988, Heiss_2002, Cejnar_2005, Cejnar_2007, Cejnar_2009, Borisov_2015, Sindelka_2017} In quantum mechanics, the Hamiltonian is almost always dependent of, at least, one parameter. In some cases the variation of a parameter can lead to abrupt changes at a critical point. These QPTs exist both for ground and excited states as shown by Cejnar and coworkers. \cite{Cejnar_2009, Sachdev_2011, Cejnar_2015, Cejnar_2016, Caprio_2008, Macek_2019} A ground-state QPT is characterized by the derivatives of the ground-state energy with respect to a non-thermal control parameter. \cite{Cejnar_2009, Sachdev_2011} The transition is called discontinuous and of first order if the first derivative is discontinuous at the critical parameter value. Otherwise, it is called continuous and of $m$th order if the $m$th derivative is discontinuous. A QPT can also be identify by the discontinuity of an appropriate order parameter (or one of its derivatives).
In the previous section, we saw that a careful analysis of the structure of the Hamiltonian allows us to predict the existence of a critical point. In a finite basis set this critical point is model by a cluster of $\beta$ singularities. It is now well known that this phenomenon is a special case of a more general phenomenon. Indeed, theoretical physicists proved that EPs close to the real axis are connected to \textit{quantum phase transitions} (QPTs). \cite{Heiss_1988,Heiss_2002,Cejnar_2005, Cejnar_2007, Cejnar_2009, Borisov_2015, Sindelka_2017} In quantum mechanics, the Hamiltonian is almost always dependent of, at least, one parameter. In some cases the variation of a parameter can lead to abrupt changes at a critical point. These QPTs exist both for ground and excited states as shown by Cejnar and coworkers. \cite{Cejnar_2009, Sachdev_2011, Cejnar_2015, Cejnar_2016, Caprio_2008, Macek_2019} A ground-state QPT is characterized by the derivatives of the ground-state energy with respect to a non-thermal control parameter. \cite{Cejnar_2009, Sachdev_2011} The transition is called discontinuous and of first order if the first derivative is discontinuous at the critical parameter value. Otherwise, it is called continuous and of $m$th order if the $m$th derivative is discontinuous. A QPT can also be identify by the discontinuity of an appropriate order parameter (or one of its derivatives).
The presence of an EP close to the real axis is characteristic of a sharp avoided crossing. Yet, at such an avoided crossing, eigenstates change abruptly. Although it is now well understood that EPs are closely related to QPTs, the link between the type of QPT (ground state or excited state, first or higher order) and EPs still need to be clarified. One of the major obstacles that one faces in order to achieve this resides in the ability to compute the distribution of EPs. The numerical assignment of an EP to two energies on the real axis is very difficult in large dimensions. Hence, the design of specific methods are required to get information on the location of EPs. Following this idea, Cejnar \textit{et al.}~developed a method based on a Coulomb analogy giving access to the density of EP close to the real axis. \cite{Cejnar_2005, Cejnar_2007} More recently Stransky and coworkers proved that the distribution of EPs is characteristic on the order of the QPT. \cite{Stransky_2018} In the thermodynamic limit, some of the EPs converge towards a critical point $\lambda_\text{c}$ on the real axis. They showed that, within the interaction boson model, \cite{Lipkin_1965} EPs associated to first- and second-order QPT behave differently when the number of particles increases. The position of these singularities converge towards the critical point on the real axis at different rates (exponentially and algebraically for the first and second orders, respectively) with respect to the number of particles.
It seems like our understanding of the physics of spatial and/or spin symmetry breaking in HF theory can be enlightened by QPT theory. Indeed, the second derivative of the HF ground-state energy is discontinuous at the point of spin symmetry-breaking which means that the system undergo a second-order QPT. Moreover, the $\beta$ singularities introduced by Sergeev and coworkers to describe the EPs modeling the formation of a bound cluster of electrons are actually a more general class of singularities. The EPs close to the real axis (the so-called $\beta$ singularities) are connected to QPT because they result from a sharp avoided crossings at which the eigenstates change quickly. However, the $\alpha$ singularities arise from large avoided crossings. Thus, they cannot be connected to QPT. The avoided crossings generating $\alpha$ singularities generally involve the ground state and low-lying doubly-excited states. Those excited states have a non-negligible contribution to the exact FCI solution because they have (usually) the same spatial and spin symmetry as the ground state. 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{The spherium model}\label{sec:spherium}
%============================================================%
Simple systems that are analytically solvable (or at least quasi-exactly solvable, \ie, models for which it is possible to obtain a finite portion of the exact solutions of the Schr{\"o}dinger equation \cite{Ushveridze_1994}) are of great importance in theoretical chemistry.
These systems are very useful to perform benchmark studies in order to test new methods as the mathematics are easier than in realistic systems (such as molecules or solids) but retain much of the key physics.
To investigate the physics of EPs we consider one such system named \textit{spherium}.
It consists of two electrons confined to the surface of a sphere interacting through the long-range Coulomb potential. \cite{Thompson_2005, Seidl_2007, Loos_2009b}
Thus, the Hamiltonian is
\begin{equation}
\hH = -\frac{\grad_1^2 + \grad_2^2}{2} + \frac{1}{r_{12}},
\end{equation}
or
\begin{equation} \label{eq:H-sph-omega}
\hH = -\frac{1}{R^2} \qty( \pdv[2]{}{\omega} + \cot \omega \pdv{}{\omega}) + \frac{1}{R \sqrt{2 - 2 \cos \omega}},
\end{equation}
where $\omega$ is the interelectronic angle.
The Laplace operators are the kinetic operators for each electrons and $r_{12}^{-1} = \abs{\vb{r}_1 - \vb{r}_2}^{-1}$ is the Coulomb operator.
Note that, as readily seen by the definition of the interelectronic distance $r_{12}$, the electrons interact through the sphere.
The radius of the sphere $R$ dictates the correlation regime. \cite{Loos_2009}
In the weak correlation regime (\ie, small $R$), the kinetic energy (which scales as $R^{-2}$) dominates and the electrons are delocalized over the sphere.
For large $R$ (or strong correlation regime), the electron repulsion term (which scales as $R^{-1}$) drives the physics and the electrons localize on opposite side of the sphere to minimize their Coulomb repulsion.
This phenomenon is sometimes referred to as a Wigner crystallization. \cite{Wigner_1934}
We will use this model in order to rationalize the effects of the parameters that may influence the physics of EPs:
\begin{itemize}
\item Partitioning of the Hamiltonian and the actual zeroth-order reference: weak correlation reference (RHF or UHF references, MP or EN partitioning), or strongly correlated reference.
\item Basis set: minimal basis or infinite (\ie, complete) basis.
\item Radius of the sphere that ultimately dictates the correlation regime.
\end{itemize}
In the RHF formalism, the two electrons are restricted to ``live'' in the same spatial orbital.
The spatial part of the RHF wave function is then
\begin{equation}\label{eq:RHF_WF}
\Psi_{\text{RHF}}(\theta_1,\theta_2) = Y_0(\theta_1) Y_0(\theta_2),
\end{equation}
where $\theta_i$ is the polar angle of the $i$th electron and $Y_{\ell}(\theta)$ is a zonal spherical harmonic.
Because $Y_0(\theta) = 1/\sqrt{4\pi}$, it is clear that the RHF wave function yields a uniform one-electron density.
The RHF wave function cannot model properly the physics of the system at large $R$ because the spatial orbitals are restricted to be the same, and, \textit{a fortiori}, it cannot represent two electrons on opposite side of the sphere.
Within the UHF formalism, there is a critical value of $R$, called Coulson-Fischer point, \cite{Coulson_1949} at which a UHF solution appears and is lower in energy than the RHF one.
The UHF solution has broken symmetry because the two electrons tends to localize on opposite sides of the sphere.
The spatial part of the UHF wave function is defined as
\begin{equation}\label{eq:UHF_WF}
\Psi_{\text{UHF}}(\theta_1,\theta_2)=\phi_\alpha(\theta_1)\phi_\beta(\theta_2),
\end{equation}
where $\phi_\sigma(\theta)$ is the spatial orbital associated with the spin-$\sigma$ electrons ($\sigma = \alpha$ for spin-up electrons and $\sigma = \beta$ for spin-down electrons).
These one-electron orbitals are expanded in the basis of zonal spherical harmonics
\begin{equation}
\phi_\sigma(\theta)=\sum_{\ell=0}^{\infty}C_{\sigma,\ell}Y_{\ell}(\theta)
\end{equation}
It is possible to obtain the formula for the HF energy in this basis set: \cite{Loos_2009}
\begin{equation}
E_{\text{HF}} = T_{\text{HF}} + V_{\text{HF}},
\label{eq:EHF}
\end{equation}
where the kinetic and potential energies are, respectively,
\begin{align}
T_{\text{HF}} & = \frac{1}{R^2} \sum_{\sigma=\alpha,\beta} \sum_{\ell=0}^{\infty} C_{\sigma,\ell}^2 \, \ell(\ell+1),
&
V_{\text{HF}} & = \frac{1}{R} \sum_{L=0}^{\infty}
v^\alpha_{L} v^\beta_{L},
\end{align}
and
\begin{equation}
v^\sigma_{L}
= \sum_{\ell_1,\ell_2} \sqrt{(2\ell_1+1)(2\ell_2+1)} C_{\sigma,\ell_1}C_{\sigma,\ell_2}
\begin{pmatrix}
\ell_1 & \ell_2 & L
\\
0 & 0 & 0
\end{pmatrix}^2
\end{equation}
is expressed in terms of the Wigner 3j-symbols. \cite{AngularBook}
The general method is to use a self-consistent field procedure as described in Ref.~\onlinecite{SzaboBook} to get the coefficients of the HF wave function corresponding to stationary solutions with respect to the coefficients $C_{\sigma,\ell}$, \ie,
\begin{equation}
\pdv{E_{\text{HF}}}{C_{\sigma,\ell}} = 0.
\end{equation}
Here, we work in a minimal basis, composed of $Y_{0}$ and $Y_{1}$, or equivalently, a s and p\textsubscript{z} orbital, to illustrate the difference between the RHF and UHF solutions. In this basis there is a shortcut to find the stationary solutions and ensure normalization of the orbitals. One can define the one-electron orbitals as
\begin{equation}
\phi_\sigma(\theta)= \cos(\chi_\sigma)Y_{0}(\theta) + \sin(\chi_\sigma)Y_{1}(\theta),
\end{equation}
using a mixing angle between the two basis functions for each spin manifold.
Hence, one has just to minimize/maximize the energy with respect to the two mixing angles $\chi_\alpha$ and $\chi_\beta$.
This process provides the three following solutions valid for all value of $R$, which are respectively a minimum, a maximum and a saddle point of the HF equations:
\begin{itemize}
\item The two electrons are in the s orbital which is a RHF solution. This solution is associated with the energy $1/R$.
\item The two electrons are in the p\textsubscript{z} orbital which is a RHF solution. This solution is associated with the energy $2/R^{2}+ 29/(25R)$.
\item One electron is in the s orbital and the other is in the p\textsubscript{z} orbital which is a UHF solution. This solution is associated with the energy $1/R^{2} + 1/R$.
\end{itemize}
In addition, the minimization process gives also the well-known symmetry-broken UHF (sb-UHF) solution. In this case the Coulson-Fischer point associated to this solution is $R=3/2$. For $R>3/2$ the sb-UHF solution is the global minimum of the HF equations and the RHF solution presented before is a local minimum. This solution corresponds to the configuration with the spin-up electron in an orbital on one side of the sphere and the spin-down electron in a miror-image orbital on the opposite side and the configuration the other way round. The electrons can be on opposite sides of the sphere because the choice of p\textsubscript{z} as a basis function induced a privileged axis on the sphere for the electrons. For $R>3/2$, this solution has the energy
\begin{equation}\label{eq:EsbUHF}
E_{\text{sb-UHF}}=-\frac{75}{112R^3}+\frac{25}{28R^2}+\frac{59}{84R}.
\end{equation}
The exact solution for the ground state is a singlet. The spherical harmonics are eigenvectors of $\hS^2$ (the spin operator) and they are associated to different eigenvalues. Yet, the symmetry-broken orbitals are linear combinations of $Y_0$ and $Y_1$. Hence, the symmetry-broken orbitals are not eigenvectors of $\hS^2$. However, this solution gives lower energies than the RHF one at large $R$ as shown in Table \ref{tab:ERHFvsEUHF} even if it does not have the exact spin symmetry. In fact, at the Coulson-Fischer point, it becomes more effective to minimize the Coulomb repulsion than the kinetic energy in order to minimize the total energy. Thus, within the HF approximation, the variational principle is allowed to break the spin symmetry because it yields a more effective minimization of the Coulomb repulsion. This type of symmetry breaking is called a spin density wave because the system oscillates between the two symmetry-broken configurations. \cite{GiulianiBook}
\begin{table*}
\caption{RHF and sb-UHF energies in the minimal basis and exact energies (in the complete basis) for various $R$.
\label{tab:ERHFvsEUHF}}
\begin{ruledtabular}
\begin{tabular}{ccccccccc}
$R$ & 0.1 & 1 & 2 & 3 & 5 & 10 & 100 & 1000 \\
\hline
RHF & 10.00000 & 1.000000 & 0.500000 & 0.333333 & 0.200000 & 0.100000 & 0.010000 & 0.001000 \\
UHF & 10.00000 & 1.000000 & 0.490699 & 0.308532 & 0.170833 & 0.078497 & 0.007112 & 0.000703 \\
Exact & 9.783874 & 0.852781 & 0.391959 & 0.247898 & 0.139471 & 0.064525 & 0.005487 & 0.000515 \\
\end{tabular}
\end{ruledtabular}
\end{table*}
There is also another symmetry-broken solution for $R>75/38$ but this one corresponds to a maximum of the HF equations. This solution is associated with another type of symmetry breaking somewhat less known. It corresponds to a configuration where both electrons are on the same side of the sphere, in the same spatial orbital. This solution is called symmetry-broken RHF (sb-RHF). The reasoning is counter-intuitive because the electrons tends to maximize their energy. The p\textsubscript{z} orbital is symmetric with respect to the center of the sphere. If the orbitals are symmetric, the maximum is when the two electrons are in the p\textsubscript{z} orbital because it maximizes the kinetic energy. At the critical value of $R$, placing the two electrons in the same symmetry-broken orbital \ie, on the same side of the sphere gives a superior energy than the p\textsubscript{z}\textsuperscript{2} state. Adding a s orbital on one side of the p\textsubscript{z} orbital to form a symmetry-broken orbital reduce the kinetic energy but increase the repulsion energy as the two electrons are more localized on one side of the sphere. It becomes more efficient to maximize the repulsion energy than the kinetic energy for $R>75/38$. This configuration breaks the spatial symmetry of charge. Hence this symmetry breaking is associated with a charge density wave, the system oscillates between the situations with the two electrons on each side. \cite{GiulianiBook}
The energy associated with this sb-RHF solution reads
\begin{equation}
E_{\text{sb-RHF}}=\frac{75}{88R^3}+\frac{25}{22R^2}+\frac{91}{66R}.
\end{equation}
\begin{figure}
\includegraphics[width=\linewidth]{EsbHF.pdf}
\caption{Energies of the five solutions of the HF equations (multiplied by $R^2$). The dotted curves correspond to the analytic continuation of the symmetry-broken solutions.}
\label{fig:SpheriumNrj}
\end{figure}
We can also consider negative values of $R$, which corresponds to the situation where one of the electrons is replaced by a positron (the anti-particle of the electron) as readily seen in Eq.~\eqref{eq:H-sph-omega}. For negative $R$ values, there are also a sb-RHF ($R<-3/2$) and a sb-UHF ($R<-75/38$) solution for negative values of $R$ (see Fig.~\ref{fig:SpheriumNrj}) but in this case the sb-RHF solution is a minimum and the sb-UHF is a maximum of the HF equations. Indeed, the sb-RHF state minimizes the attraction energy by placing the electron and the positron on the same side of the sphere. And the sb-UHF state maximizes the energy because the two attracting particles are on opposite sides of the sphere.
In addition, we can also consider the symmetry-broken solutions beyond their respective Coulson-Fischer points by analytically continuing their respective energies leading to the so-called holomorphic solutions. \cite{Hiscock_2014, Burton_2019, Burton_2019a} All those energies are plotted in Fig.~\ref{fig:SpheriumNrj}. The dotted curves corresponds to the holomorphic domain of the energies.
\section{Radius of convergence and exceptional points}
\subsection{Evolution of the radius of convergence}
In this subsection, we investigate how the partitioning of $\hH(\lambda)$ influence the radius of convergence of the perturbation series. Let us remind the reader that the radius of convergence is equal to the distance of the closest singularity to the origin of $E(\lambda)$. Hence, we have to determine the locations of the EPs to obtain information on the convergence properties of the perturbative series. To find them we solve simultaneously the following equations: \cite{Cejnar_2007}
\begin{subequations}
\begin{align}
\label{eq:PolChar}
\det[E\hI-\hH(\lambda)] & = 0,
\\
\label{eq:DPolChar}
\pdv{E}\det[E\hI-\hH(\lambda)] & = 0,
\end{align}
\end{subequations}
where $\hI$ is the identity operator.
Equation \eqref{eq:PolChar} is the well-known secular equation providing us with the (eigen)energies of the system. If an energy is also solution of Eq.~\eqref{eq:DPolChar}, then this energy is, at least two-fold degenerate. In this case the energies obtained are dependent of $\lambda$. Thus, solving these equations with respect to $E$ and $\lambda$ gives the value of $\lambda$ where two energies are degenerate. These degeneracies can be conical intersections between two states with different symmetry for real values of $\lambda$ \cite{Yarkony_1996} or EPs between two states with the same symmetry for complex values of $\lambda$.
Let us assume that electron 1 is spin-up and electron 2 is spin-down.
Hence, we can forget about the spin part of the spin-orbitals and from now on we will work with spatial orbitals. In the restricted formalism the spatial orbitals are the same so the two-electron basis set can be defined as
\begin{align}\label{eq:rhfbasis}
\psi_1 & =Y_{0}(\theta_1)Y_{0}(\theta_2),
&
\psi_2 & =Y_{0}(\theta_1)Y_{1}(\theta_2),\\
\psi_3 & =Y_{1}(\theta_1)Y_{0}(\theta_2),
&
\psi_4 & =Y_{1}(\theta_1)Y_{1}(\theta_2).
\end{align}
The Hamiltonian $\hH(\lambda)$ is block diagonal in this basis because of its symmetry, \ie, $\psi_1$ only interacts with $\psi_4$, and $\psi_2$ with $\psi_3$. The two singly-excited states yield, after diagonalization, a spatially anti-symmetric singlet $sp_z$ and a spatially symmetric triplet sp\textsubscript{z} state. Hence those states do not have the same symmetry as the spatially symmetric singlet ground state. Thus, these states cannot be involved in an avoided crossing with the ground state as can be seen in Fig.~\ref{fig:RHFMiniBas} and, \textit{a fortiori} cannot be involved in an EP with the ground state. However there is an avoided crossing between the s\textsuperscript{2} and p\textsubscript{z}\textsuperscript{2} states which gives two EPs in the complex plane.
\begin{figure}
\includegraphics[width=\linewidth]{EMP_RHF_R10.pdf}
\caption{Energies $E(\lambda)$ in the restricted basis set \eqref{eq:rhfbasis} with $R=10$.
One can clearly see the avoided crossing between the s\textsuperscript{2} and p\textsubscript{z}\textsuperscript{2} states around $\lambda = 1$.}
\label{fig:RHFMiniBas}
\end{figure}
To simplify the problem, it is convenient to only consider basis functions of a given symmetry. Such basis functions are called configuration state functions (CSFs). It simplifies greatly the problem because, with such a basis set, one only gets the degeneracies of interest associated with the convergence properties, \ie, the EPs between states with the same symmetry as the ground state. In the present context, the ground state is a totally symmetric singlet. According to angular momentum theory, \cite{AngularBook, SlaterBook, Loos_2009} we expand the exact wave function in the following two-electron basis:
\begin{equation}
\Phi_\ell(\omega)=\frac{\sqrt{2\ell+1}}{4\pi R^2}P_\ell(\cos\omega),
\end{equation}
where $P_\ell$ are Legendre polynomials.
Then, using this basis set we can compare the different partitioning of Sec.~\ref{sec:AlterPart}. Figure \ref{fig:RadiusPartitioning} shows the evolution of the radius of convergence $R_{\text{CV}}$ as a function of $R$ for the MP, the EN, the WC and the SC partitioning in a minimal basis (\ie, consisting of $P_0$ and $P_1$) of size $K = 2$, and in the same basis augmented with $P_2$ ($K = 3$). We see that, for the SC partitioning, $R_{\text{CV}}$ increases with $R$ whereas it is decreasing for the three others partitioning. This result is expected because the MP, EN, and WC partitioning use a weakly correlated reference so $\hH^{(0)}$ is a good approximation for small $R$. On the contrary, the SC partitioning consider naturally a strongly correlated reference so the SC series converges far better when the electron are strongly correlated, \ie, when $R$ is large in the spherium model.
Interestingly, the radius of convergence associated with the SC partitioning is greater than one for a greater range of radii for $K = 2$ than $K = 3$.
\begin{figure}
\includegraphics[width=0.49\textwidth]{PartitioningRCV2.pdf}
\includegraphics[width=0.49\textwidth]{PartitioningRCV3.pdf}
\caption{Radius of convergence $R_{\text{CV}}$ for two (left) and three (right) basis functions for various partitionings.}
\label{fig:RadiusPartitioning}
\end{figure}
The MP partitioning is always better than WC in Fig.~\ref{fig:RadiusPartitioning}. In the WC partitioning the powers of $R$ (the zeroth-order scales as $R^{-2}$ while the perturbation scales as $R^{-1}$) are well-separated so each term of the series has a well-defined power of $R$. This is not the case for the MP series.
Interestingly, it can be proved that the $m$th order energy of the WC series can be obtained as a Taylor series of MP$m$ with respect to $R$.
It seems that the EN is better than MP for very small $R$ in the minimal basis. In fact, it is just an artefact of the minimal basis because, for $K = 3$, the MP series has a greater radius of convergence for all values of $R$. It holds true for $K>3$.
Figure \ref{fig:RadiusBasis} shows that the radius of convergence is not very sensitive to the size of the basis set. The CSFs have all the same spin and spatial symmetries so we expect that the singularities obtained within this basis set will be $\alpha$ singularities. Table \ref{tab:SingAlpha} shows that the singularities considered in this case are indeed $\alpha$ singularities. This is consistent with the observation of Goodson and Sergeev \cite{Goodson_2004} who stated that $\alpha$ singularities are relatively insensitive to the basis set size. The discontinuities observed in Fig.~\ref{fig:RadiusBasis} for the MP partitioning are due to changes in dominant singularity. We can observe this change in Table \ref{tab:SingAlpha}, the value for $R=1$ and $R=2$ are respectively in the positive and negative planes.
\begin{figure}
\includegraphics[width=0.49\textwidth]{MPlargebasis.pdf}
\includegraphics[width=0.49\textwidth]{WCElargebasis.pdf}
\caption{Radius of convergence $R_{\text{CV}}$ in the CSF basis with $K$ basis functions for the MP (left) and WC (right) partitioning.}
\label{fig:RadiusBasis}
\end{figure}
\begin{table}
\footnotesize
\caption{Dominant singularity in the CSF basis set ($K=8$) for various value of $R$ in the MP and WC partitioning.}
\begin{tabular}{cccccccc}
\hline
\hline
$R$ & 0.1 & 1 & 2 & 3 & 5 & 10 & 100 \\
\hline
MP & $+14.1-10.9\,i$ & $+2.38-1.47\,i$ & $-0.67-1.30\,i$ & $-0.49-0.89\,i$ & $-0.33-0.55\,i$ & $-0.22-0.31\,i$ & $+0.03-0.05\,i$ \\
WC & $-9.6-10.7\,i$ & $-0.96-1.07\,i$ & $-0.48-0.53\,i$ & $-0.32-0.36\,i$ & $-0.19-0.21\,i$ & $-0.10-0.11\,i$ & $-0.01-0.01\,i$ \\
\hline
\hline
\end{tabular}
\label{tab:SingAlpha}
\end{table}
\subsection{Exceptional points in the UHF formalism}\label{sec:uhfSing}
Now, we investigate the differences in the singularity structure between the RHF and UHF formalism. To do so, we use the symmetry-broken orbitals discussed in Sec.~\ref{sec:spherium}. Thus, the UHF two-electron basis is
\begin{align}\label{eq:uhfbasis}
\psi_1 & =\phi_{\alpha,1}(\theta_1)\phi_{\beta,1}(\theta_2),
&
\psi_2 & =\phi_{\alpha,1}(\theta_1)\phi_{\beta,2}(\theta_2),\\
\psi_3 & =\phi_{\alpha,2}(\theta_1)\phi_{\beta,1}(\theta_2),
&
\psi_4 & =\phi_{\alpha,2}(\theta_1)\phi_{\beta,2}(\theta_2).
\end{align}
with the symmetry-broken orbitals
\begin{subequations}
\begin{align}\label{eq:uhforbitals}
\phi_{\alpha,1}(\theta)
& =\frac{\sqrt{75+62R}}{4\sqrt{7R}} Y_{0}(\theta)
+ \frac{5\sqrt{-3+2R}}{4\sqrt{7R}} Y_{1}(\theta),
\\
\phi_{\beta,1}(\theta)
& =\frac{\sqrt{75+62R}}{4\sqrt{7R}} Y_{0}(\theta)
- \frac{5\sqrt{-3+2R}}{4\sqrt{7R}} Y_{1}(\theta),
\\
\phi_{\alpha,2}(\theta)
& = - \frac{5\sqrt{-3+2R}}{4\sqrt{7R}} Y_{0}(\theta)
+ \frac{\sqrt{75+62R}}{4\sqrt{7R}} Y_{1}(\theta),
\\
\phi_{\beta,2}(\theta)
& =\frac{5\sqrt{-3+2R}}{4\sqrt{7R}} Y_{0}(\theta)
+\frac{\sqrt{75+62R}}{4\sqrt{7R}} Y_{1}(\theta).
\end{align}
\end{subequations}
In the UHF formalism the Hamiltonian $\hH(\lambda)$ is no more block diagonal, $\psi_4$ can interact with $\psi_2$ and $\psi_3$. The matrix elements of the Hamiltonian corresponding to this interaction are
\begin{equation}\label{eq:MatrixElem}
H_{24}=H_{34}=H_{42}=H_{43}=\sqrt{R-\frac{3}{2}}\sqrt{R+\frac{75}{62}}\qty(R+\frac{25}{2})\frac{\sqrt{31}}{70R^3}
\end{equation}
For $R=3/2$ the Hamiltonian is block diagonal because the matrix elements \eqref{eq:MatrixElem} are equal to zero so this is equivalent to the RHF case but for $R>3/2$ the matrix elements become real. This interaction corresponds to the spin contamination of the wave function. For $R<3/2$ the matrix elements are complex, this corresponds to the holomorphic solution of Fig.~\ref{fig:SpheriumNrj}, the singularities in this case will be treated later. The matrix elements become real again for $R<-75/62$, this corresponds to the sb-UHF solution for negative value of $R$ observed in Sec.~\ref{sec:spherium}. We will refer to the domain where the matrix elements are complex as the holomorphic domain.
The singularity structure in this case is more complex because of the spin contamination of the wave function. We can not use CSFs in this case. So when one compute all the degeneracies using Eqs.~\eqref{eq:PolChar} and \eqref{eq:DPolChar} some correspond to EPs and some correspond to conical intersections. The numerical distinction of those singularities is very difficult. We will first look at the energies $E(\lambda)$ obtained with this basis set to attribute a physical signification to the singularities obtained numerically.
Figure \ref{fig:UHFMiniBas} is the analog of Fig.~\ref{fig:RHFMiniBas} in the UHF formalism. We see that in this case the sp\textsubscript{z} triplet interacts with the s\textsuperscript{2} and the p\textsubscript{z}\textsuperscript{2} singlets. Those avoided crossings are due to the spin contamination of the wave function.
\begin{figure}
\includegraphics[width=\linewidth]{EMP_UHF_R10.pdf}
\caption{Energies $E(\lambda)$ in the unrestricted basis set \eqref{eq:uhfbasis} with $R=10$.}
\label{fig:UHFMiniBas}
\end{figure}
Within the RHF formalism, we have observed only $\alpha$ singularities and large avoided crossings but one can see in Fig.~\ref{fig:UHFMiniBas} that in the UHF case there are sharp avoided crossings which are known to be connected to $\beta$ singularities. For example at $R=10$ the pair of singularities connected to the avoided crossing between s\textsuperscript{2} and sp\textsubscript{z} $^{3}P$ is $0.999\pm0.014\,i$. And the one between sp\textsubscript{z} $^{3}P$ and p\textsubscript{z}\textsuperscript{2} is connected with the singularities $2.207\pm0.023\,i$. However, in spherium, the electrons cannot be ionized so those singularities cannot be the same $\beta$ singularities as the ones highlighted by Sergeev and Goodson. \cite{Sergeev_2005} We can see in Fig.~\ref{fig:UHFEP} that the s\textsuperscript{2} singlet and the sp\textsubscript{z} triplet states are degenerated for $R=3/2$. For $R>3/2$, it becomes an avoided crossing on the real axis and the degeneracies are $moved$ in the complex plane. The wave function is spin contaminated by $Y_1$ for $R>3/2$ this is why the s\textsuperscript{2} singlet energy cannot cross the sp\textsubscript{z} triplet curves anymore. When $R$ increases this avoided crossing becomes sharper. As presented before $\beta$ singularities are linked to quantum phase transition so it seems that this singularity is linked to the spin symmetry breaking of the UHF wave function. The fact that a similar pair of $\beta$ singularities appears for $R<-75/62$ confirms this assumption. The sharp avoided crossing between sp\textsubscript{z} $^{3}P$ and p\textsubscript{z}\textsuperscript{2} is not present on Fig. \ref{fig:UHFEP}. The second pair of $\beta$ singularities resulting from this avoided crossing appears for $R\gtrsim 2.5$, this is probably due to an excited-state quantum phase transition but this still need to be investigated.
\begin{figure}
\includegraphics[width=0.45\textwidth]{UHFCI.pdf}
\includegraphics[width=0.45\textwidth]{UHFEP.pdf}
\caption{Energies $E(\lambda)$ in the unrestricted basis set \eqref{eq:uhfbasis} for $R=1.5$ (left) and $R=1.51$ (right).}
\label{fig:UHFEP}
\end{figure}
As shown before, some matrix elements of the Hamiltonian become complex in the holomorphic domain. Therefore the Hamiltonian becomes non-Hermitian for these values of $R$. In Ref.~\onlinecite{Burton_2019a}, Burton \textit{et al.}~proved that for the \ce{H_2} molecule the unrestricted Hamiltonian is not \pt -symmetric in the holomorphic domain. An analog reasoning can be done with the spherium model to prove the same result. The \pt -symmetry (invariance with respect to combined space reflection $\mathcal{P}$ and time reversal $\mathcal{T}$) is a property which ensures that a non-Hermitian Hamiltonian has a real energy spectrum. \cite{BenderPTBook} Thus \pt -symmetric Hamiltonians can be seen as an intermediate class between Hermitian and non-Hermitian Hamiltonians.
Figure \ref{fig:UHFPT} shows that for the spherium model a part of the energy spectrum becomes complex when $R$ is in the holomorphic domain. The domain of values where the energy becomes complex is called the broken \pt-symmetry region. This is consistent with the fact that in the holomorphic domain the Hamiltonian is no more \pt -symmetric.
For a non-Hermitian Hamiltonian the EPs can lie on the real axis. In particular, at the point of {\pt} transition (the point where the energies become complex) the two energies are degenerate resulting in such an EP on the real axis. This degeneracy can be seen in Fig.~\ref{fig:UHFPT}.
\begin{figure}[h!]
\includegraphics[width=0.45\textwidth]{ReNRJPT.pdf}
\includegraphics[width=0.45\textwidth]{ImNRJPT.pdf}
\caption{Real part (left) and imaginary part (right) of $E(\lambda)$ in the unrestricted basis set \eqref{eq:uhfbasis} for $R=1$.}
\label{fig:UHFPT}
\end{figure}
%%============================================================%
%\section{The spherium model}
%\label{sec:spherium}
%%============================================================%
%
%Simple systems that are analytically solvable (or at least quasi-exactly solvable, \ie, models for which it is possible to obtain a finite portion of the exact solutions of the Schr{\"o}dinger equation \cite{Ushveridze_1994}) are of great importance in theoretical chemistry.
%These systems are very useful to perform benchmark studies in order to test new methods as the mathematics are easier than in realistic systems (such as molecules or solids) but retain much of the key physics.
%To investigate the physics of EPs we consider one such system named \textit{spherium}.
%It consists of two electrons confined to the surface of a sphere interacting through the long-range Coulomb potential. \cite{Thompson_2005, Seidl_2007, Loos_2009b}
%Thus, the Hamiltonian is
%\begin{equation}
% \hH = -\frac{\grad_1^2 + \grad_2^2}{2} + \frac{1}{r_{12}},
%\end{equation}
%or
%\begin{equation} \label{eq:H-sph-omega}
% \hH = -\frac{1}{R^2} \qty( \pdv[2]{}{\omega} + \cot \omega \pdv{}{\omega}) + \frac{1}{R \sqrt{2 - 2 \cos \omega}},
%\end{equation}
%in term of the interelectronic angle $\omega$.
%The Laplace operators are the kinetic operators for each electron and $r_{12}^{-1} = \abs{\vb{r}_1 - \vb{r}_2}^{-1}$ is the Coulomb operator.
%Note that, as readily seen by the definition of the interelectronic distance $r_{12}$, the electrons interact through the sphere.
%The radius of the sphere $R$ dictates the correlation regime. \cite{Loos_2009}
%In the weak correlation regime (\ie, small $R$), the kinetic energy (which scales as $R^{-2}$) dominates and the electrons are delocalized over the sphere.
%For large $R$ (or strong correlation regime), the electron repulsion term (which scales as $R^{-1}$) drives the physics and the electrons localize on opposite side of the sphere to minimize their Coulomb repulsion.
%This phenomenon is sometimes referred to as a Wigner crystallization. \cite{Wigner_1934}
%
%We will use this model in order to rationalize the effects of the parameters that may influence the physics of EPs:
%i) Partitioning of the Hamiltonian and the actual zeroth-order reference: weak correlation reference (RHF or UHF references, MP or EN partitioning), or strongly correlated reference.
%ii) Basis set: minimal basis or infinite (\ie, complete) basis.
%iii) Radius of the sphere that ultimately dictates the correlation regime.
%
%In the RHF formalism, the two electrons are restricted to ``live'' in the same spatial orbital.
%The spatial part of the RHF wave function is then
%\begin{equation}\label{eq:RHF_WF}
% \Psi_{\text{RHF}}(\theta_1,\theta_2) = Y_0(\theta_1) Y_0(\theta_2),
%\end{equation}
%where $\theta_i$ is the polar angle of the $i$th electron and $Y_{\ell}(\theta)$ is a zonal spherical harmonic.
%Because $Y_0(\theta) = 1/\sqrt{4\pi}$, it is clear that the RHF wave function yields a uniform one-electron density.
%
%The RHF wave function cannot model properly the physics of the system at large $R$ because the spatial orbitals are restricted to be the same, and, \textit{a fortiori}, it cannot represent two electrons on opposite side of the sphere.
%Within the UHF formalism, there is a critical value of $R$, called Coulson-Fischer point, \cite{Coulson_1949} at which a UHF solution appears and is lower in energy than the RHF one.
%The UHF solution has broken symmetry because the two electrons tends to localize on opposite sides of the sphere.
%The spatial part of the UHF wave function is defined as
%\begin{equation}\label{eq:UHF_WF}
% \Psi_{\text{UHF}}(\theta_1,\theta_2)=\phi_\alpha(\theta_1)\phi_\beta(\theta_2),
%\end{equation}
%where $\phi_\sigma(\theta)$ is the spatial orbital associated with the spin-$\sigma$ electrons ($\sigma = \alpha$ for spin-up electrons and $\sigma = \beta$ for spin-down electrons).
%These one-electron orbitals are expanded in the basis of zonal spherical harmonics
%\begin{equation}
% \phi_\sigma(\theta)=\sum_{\ell=0}^{\infty}C_{\sigma,\ell}Y_{\ell}(\theta).
%\end{equation}
%It is possible to obtain the formula for the HF energy in this basis set: \cite{Loos_2009}
%\begin{equation}
% E_{\text{HF}} = T_{\text{HF}} + V_{\text{HF}},
%\label{eq:EHF}
%\end{equation}
%where the kinetic and potential energies are, respectively,
%\begin{align}
% T_{\text{HF}} & = \frac{1}{R^2} \sum_{\sigma=\alpha,\beta} \sum_{\ell=0}^{\infty} C_{\sigma,\ell}^2 \, \ell(\ell+1),
% &
% V_{\text{HF}} & = \frac{1}{R} \sum_{L=0}^{\infty}
% v^\alpha_{L} v^\beta_{L},
%\end{align}
%and
%\begin{equation}
% v^\sigma_{L}
% = \sum_{\ell_1,\ell_2} \sqrt{(2\ell_1+1)(2\ell_2+1)} C_{\sigma,\ell_1}C_{\sigma,\ell_2}
% \begin{pmatrix}
% \ell_1 & L & \ell_2
% \\
% 0 & 0 & 0
% \end{pmatrix}^2
%\end{equation}
%is expressed in terms of the Wigner 3j-symbols. \cite{AngularBook}
%
%The general method is to use a self-consistent field procedure as described in Ref.~\onlinecite{SzaboBook} to get the coefficients of the HF wave function corresponding to stationary solutions with respect to the coefficients $C_{\sigma,\ell}$, \ie,
%\begin{equation}
% \pdv{E_{\text{HF}}}{C_{\sigma,\ell}} = 0.
%\end{equation}
%Here, we work in a minimal basis, composed of $Y_{0}$ and $Y_{1}$, or equivalently, a $s$ and $p_{z}$ orbital, to illustrate the difference between the RHF and UHF solutions. In this basis there is a shortcut to find the stationary solutions and ensure normalization of the orbitals. One can define the one-electron orbitals as
%\begin{equation}
% \phi_\sigma(\theta)= \cos(\chi_\sigma)Y_{0}(\theta) + \sin(\chi_\sigma)Y_{1}(\theta),
%\end{equation}
%using a mixing angle $\chi_\sigma$ between the two basis functions for each spin manifold.
%Hence, one has just to minimize/maximize the energy with respect to the two mixing angles $\chi_\alpha$ and $\chi_\beta$.
%
%This process provides the three following solutions valid for all value of $R$, which are respectively a minimum, a maximum and a saddle point of the HF equations:
%i) the two electrons are in the $s$ orbital which is a RHF solution. This solution is associated with the energy $1/R$;
%ii) the two electrons are in the $p_{z}$ orbital which is a RHF solution. This solution is associated with the energy $2/R^{2}+ 29/(25R)$;
%iii) one electron is in the $s$ orbital and the other is in the $p_{z}$ orbital which is a UHF solution. This solution is associated with the energy $1/R^{2} + 1/R$.
%
%In addition, the minimization process gives also the well-known symmetry-broken UHF (sb-UHF) solution. In this case the Coulson-Fischer point associated to this solution is $R=3/2$.
%For $R>3/2$, the sb-UHF solution is the global minimum of the HF equations and the RHF solution presented before is a local minimum. This solution corresponds to the configuration with the spin-up electron in an orbital on one side of the sphere and the spin-down electron in a miror-image orbital on the opposite side and the configuration the other way round. The electrons can be on opposite sides of the sphere because the choice of $p_{z}$ as a basis function induced a privileged axis on the sphere for the electrons. For $R>3/2$, this solution has the energy
%\begin{equation}\label{eq:EsbUHF}
% E_{\text{sb-UHF}}=-\frac{75}{112R^3}+\frac{25}{28R^2}+\frac{59}{84R}.
%\end{equation}
%
%The exact solution for the ground state is a singlet. The spherical harmonics are eigenvectors of $\hS^2$ (the spin operator) and they are associated to different eigenvalues.
%Yet, the symmetry-broken orbitals are linear combinations of $Y_0$ and $Y_1$.
%Hence, the symmetry-broken orbitals are not eigenvectors of $\hS^2$.
%However, this solution gives lower energies than the RHF one at large $R$, even if it does not have the exact spin symmetry.
%In fact, at the Coulson-Fischer point, it becomes more effective to minimize the Coulomb repulsion than the kinetic energy in order to minimize the total energy.
%Thus, within the HF approximation, the variational principle is allowed to break the spin symmetry because it yields a more effective minimization of the Coulomb repulsion.
%This type of symmetry breaking is also called a spin-density wave in the physics community as the system ``oscillates'' between the two symmetry-broken configurations. \cite{GiulianiBook}
%
%There is also another symmetry-broken solution for $R>75/38$ but this one corresponds to a maximum of the HF equations.
%This solution is associated with another type of symmetry breaking somewhat less known.
%It corresponds to a configuration where both electrons are on the same side of the sphere, in the same spatial orbital.
%This solution is called symmetry-broken RHF (sb-RHF). The reasoning is counter-intuitive because the electrons tends to maximize their energy.
%The $sp_{z}$ orbital is symmetric with respect to the center of the sphere.
%If the orbitals are symmetric, the maximum is when the two electrons are in the $p_{z}$ orbital because it maximizes the kinetic energy.
%At the critical value of $R$, placing the two electrons in the same symmetry-broken orbital \ie, on the same side of the sphere gives a superior energy than the $p_{z}^2$ state. Adding a s orbital on one side of the $p_{z}$ orbital to form a symmetry-broken orbital reduce the kinetic energy but increase the repulsion energy as the two electrons are more localized on one side of the sphere.
%It becomes more efficient to maximize the repulsion energy than the kinetic energy for $R>75/38$.
%This configuration breaks the spatial symmetry of charge.
%Hence this symmetry breaking is associated with a charge-density wave as the system oscillates between the situations where the two electrons are one side or the other. \cite{GiulianiBook}
%The energy associated with this sb-RHF solution reads
%\begin{equation}
%E_{\text{sb-RHF}}=\frac{75}{88R^3}+\frac{25}{22R^2}+\frac{91}{66R}.
%\end{equation}
%\begin{figure}
% \includegraphics[width=\linewidth]{EsbHF.pdf}
% \caption{Energies of the five solutions of the HF equations (multiplied by $R^2$). The dotted curves correspond to the analytic continuation of the symmetry-broken solutions.}
% \label{fig:SpheriumNrj}
%\end{figure}
%
%We can also consider negative values of $R$, which corresponds to the situation where one of the electrons is replaced by a positron as readily seen in Eq.~\eqref{eq:H-sph-omega}.
%For negative $R$ values, there are also a sb-RHF ($R<-3/2$) and a sb-UHF ($R<-75/38$) solution for negative values of $R$ (see Fig.~\ref{fig:SpheriumNrj}) but in this case the sb-RHF solution is a minimum and the sb-UHF is a maximum of the HF equations.
%Indeed, the sb-RHF state minimizes the attraction energy by placing the electron and the positron on the same side of the sphere.
%And the sb-UHF state maximizes the energy because the two attracting particles are on opposite sides of the sphere.
%
%In addition, we can also consider the symmetry-broken solutions beyond their respective Coulson-Fischer points by analytically continuing their respective energies leading to the so-called holomorphic solutions. \cite{Hiscock_2014, Burton_2019, Burton_2019a} All those energies are plotted in Fig.~\ref{fig:SpheriumNrj}. The dotted curves corresponds to the holomorphic domain of the energies.
%
%
%\section{Radius of convergence and exceptional points}
%
%\subsection{Evolution of the radius of convergence}
%
%In this subsection, we investigate how the partitioning of $\hH(\lambda)$ influence the radius of convergence of the perturbation series. Let us remind the reader that the radius of convergence is equal to the distance of the closest singularity to the origin of $E(\lambda)$. Hence, we have to determine the locations of the EPs to obtain information on the convergence properties of the perturbative series. To find them we solve simultaneously the following equations: \cite{Cejnar_2007}
%\begin{subequations}
%\begin{align}
% \label{eq:PolChar}
% \det[E\hI-\hH(\lambda)] & = 0,
% \\
% \label{eq:DPolChar}
% \pdv{E}\det[E\hI-\hH(\lambda)] & = 0,
%\end{align}
%\end{subequations}
%where $\hI$ is the identity operator.
%Equation \eqref{eq:PolChar} is the well-known secular equation providing us with the (eigen)energies of the system. If an energy is also solution of Eq.~\eqref{eq:DPolChar}, then this energy is, at least, two-fold degenerate. In this case the energies obtained are $\lambda$-dependent.
%Thus, solving these equations with respect to $E$ and $\lambda$ gives the value of $\lambda$ where two energies are degenerate.
%These degeneracies can be conical intersections between two states with different symmetries for real values of $\lambda$ \cite{Yarkony_1996} or EPs between two states with the same symmetry for complex values of $\lambda$.
%
%Let us assume that electron 1 is spin-up and electron 2 is spin-down.
%Hence, we can forget about the spin part of the spin-orbitals and from now on we will work with spatial orbitals. In the restricted formalism the spatial orbitals are the same so the two-electron basis set can be defined as
%
%\begin{align}\label{eq:rhfbasis}
% \psi_1 & =Y_{0}(\theta_1)Y_{0}(\theta_2),
% &
% \psi_2 & =Y_{0}(\theta_1)Y_{1}(\theta_2),\\
% \psi_3 & =Y_{1}(\theta_1)Y_{0}(\theta_2),
% &
% \psi_4 & =Y_{1}(\theta_1)Y_{1}(\theta_2).
%\end{align}
%The Hamiltonian $\hH(\lambda)$ is block diagonal in this basis because of its symmetry, \ie, $\psi_1$ only interacts with $\psi_4$, and $\psi_2$ with $\psi_3$. The two singly-excited states yield, after diagonalization, a spatially anti-symmetric singlet $sp_z$ and a spatially symmetric triplet $sp_{z}$ state.
%Hence those states do not have the same symmetry as the spatially symmetric singlet ground state.
%Thus, these states cannot be involved in an avoided crossing with the ground state as can be seen in Fig.~\ref{fig:RHFMiniBas} and, \textit{a fortiori} cannot be involved in an EP with the ground state.
%However there is an avoided crossing between the $s^{2}$ and $p_{z}^{2}$ states which gives two EPs in the complex plane.
%
%\begin{figure}
% \includegraphics[width=\linewidth]{EMP_RHF_R10.pdf}
% \caption{Energies $E(\lambda)$ in the restricted basis set \eqref{eq:rhfbasis} with $R=10$.
% One can clearly see the avoided crossing between the $s^{2}$ and $p_{z}^{2}$ states around $\lambda = 1$.}
% \label{fig:RHFMiniBas}
%\end{figure}
%
%To simplify the problem, it is convenient to only consider basis functions of a given symmetry. Such basis functions are called configuration state functions (CSFs). It simplifies greatly the problem because, with such a basis set, one only gets the degeneracies of interest associated with the convergence properties, \ie, the EPs between states with the same symmetry as the ground state. In the present context, the ground state is a totally symmetric singlet. According to angular momentum theory, \cite{AngularBook, SlaterBook, Loos_2009} we expand the exact wave function in the following two-electron basis:
%\begin{equation}
%\Phi_\ell(\omega)=\frac{\sqrt{2\ell+1}}{4\pi R^2}P_\ell(\cos\omega),
%\end{equation}
%where $P_\ell$ are Legendre polynomials.
%
%Then, using this basis set we can compare the different partitioning of Sec.~\ref{sec:AlterPart}. Figure \ref{fig:RadiusPartitioning} shows the evolution of the radius of convergence $R_{\text{CV}}$ as a function of $R$ for the MP, the EN, the WC and the SC partitioning in a minimal basis (\ie, consisting of $P_0$ and $P_1$) of size $K = 2$, and in the same basis augmented with $P_2$ ($K = 3$). We see that, for the SC partitioning, $R_{\text{CV}}$ increases with $R$ whereas it is decreasing for the three others partitioning. This result is expected because the MP, EN, and WC partitioning use a weakly correlated reference so $\hH^{(0)}$ is a good approximation for small $R$. On the contrary, the SC partitioning consider naturally a strongly correlated reference so the SC series converges far better when the electron are strongly correlated, \ie, when $R$ is large in the spherium model.
%
%Interestingly, the radius of convergence associated with the SC partitioning is greater than one for a greater range of radii for $K = 2$ than $K = 3$.
%
%\begin{figure}
% \includegraphics[width=0.49\textwidth]{PartitioningRCV2.pdf}
% \includegraphics[width=0.49\textwidth]{PartitioningRCV3.pdf}
% \caption{Radius of convergence $R_{\text{CV}}$ for two (left) and three (right) basis functions for various partitionings.}
% \label{fig:RadiusPartitioning}
%\end{figure}
%
%The MP partitioning is always better than WC in Fig.~\ref{fig:RadiusPartitioning}. In the WC partitioning the powers of $R$ (the zeroth-order scales as $R^{-2}$ while the perturbation scales as $R^{-1}$) are well-separated so each term of the series has a well-defined power of $R$. This is not the case for the MP series.
%Interestingly, it can be proved that the $m$th order energy of the WC series can be obtained as a Taylor series of MP$m$ with respect to $R$.
%It seems that the EN is better than MP for very small $R$ in the minimal basis. In fact, it is just an artefact of the minimal basis because, for $K = 3$, the MP series has a greater radius of convergence for all values of $R$. It holds true for $K>3$.
%
%
%Figure \ref{fig:RadiusBasis} shows that the radius of convergence is not very sensitive to the size of the basis set. The CSFs have all the same spin and spatial symmetries so we expect that the singularities obtained within this basis set will be $\alpha$ singularities. Table \ref{tab:SingAlpha} shows that the singularities considered in this case are indeed $\alpha$ singularities. This is consistent with the observation of Goodson and Sergeev \cite{Goodson_2004} who stated that $\alpha$ singularities are relatively insensitive to the basis set size. The discontinuities observed in Fig.~\ref{fig:RadiusBasis} for the MP partitioning are due to changes in dominant singularity. We can observe this change in Table \ref{tab:SingAlpha}, the value for $R=1$ and $R=2$ are respectively in the positive and negative planes.
%
%\begin{figure}
% \includegraphics[width=0.49\textwidth]{MPlargebasis.pdf}
% \includegraphics[width=0.49\textwidth]{WCElargebasis.pdf}
% \caption{Radius of convergence $R_{\text{CV}}$ in the CSF basis with $K$ basis functions for the MP (left) and WC (right) partitioning.}
% \label{fig:RadiusBasis}
%\end{figure}
%
%\begin{table*}
%\caption{Dominant singularity in the CSF basis set ($K=8$) for various value of $R$ in the MP and WC partitioning.}
%\begin{ruledtabular}
%\begin{tabular}{cccccccc}
%$R$ & 0.1 & 1 & 2 & 3 & 5 & 10 & 100 \\
%\hline
%MP & $+14.1-10.9\,i$ & $+2.38-1.47\,i$ & $-0.67-1.30\,i$ & $-0.49-0.89\,i$ & $-0.33-0.55\,i$ & $-0.22-0.31\,i$ & $+0.03-0.05\,i$ \\
%WC & $-9.6-10.7\,i$ & $-0.96-1.07\,i$ & $-0.48-0.53\,i$ & $-0.32-0.36\,i$ & $-0.19-0.21\,i$ & $-0.10-0.11\,i$ & $-0.01-0.01\,i$ \\
%\end{tabular}
%\end{ruledtabular}
%\label{tab:SingAlpha}
%\end{table*}
%
%\subsection{Exceptional points in the UHF formalism}\label{sec:uhfSing}
%
%Now, we investigate the differences in the singularity structure between the RHF and UHF formalism. To do so, we use the symmetry-broken orbitals discussed in Sec.~\ref{sec:spherium}. Thus, the UHF two-electron basis is
%\begin{align}\label{eq:uhfbasis}
% \psi_1 & =\phi_{\alpha,1}(\theta_1)\phi_{\beta,1}(\theta_2),
% &
% \psi_2 & =\phi_{\alpha,1}(\theta_1)\phi_{\beta,2}(\theta_2),\\
% \psi_3 & =\phi_{\alpha,2}(\theta_1)\phi_{\beta,1}(\theta_2),
% &
% \psi_4 & =\phi_{\alpha,2}(\theta_1)\phi_{\beta,2}(\theta_2).
%\end{align}
%with the symmetry-broken orbitals
%\begin{subequations}
%\begin{align}\label{eq:uhforbitals}
% \phi_{\alpha,1}(\theta)
% & =\frac{\sqrt{75+62R}}{4\sqrt{7R}} Y_{0}(\theta)
% + \frac{5\sqrt{-3+2R}}{4\sqrt{7R}} Y_{1}(\theta),
% \\
% \phi_{\beta,1}(\theta)
% & =\frac{\sqrt{75+62R}}{4\sqrt{7R}} Y_{0}(\theta)
% - \frac{5\sqrt{-3+2R}}{4\sqrt{7R}} Y_{1}(\theta),
% \\
% \phi_{\alpha,2}(\theta)
% & = - \frac{5\sqrt{-3+2R}}{4\sqrt{7R}} Y_{0}(\theta)
% + \frac{\sqrt{75+62R}}{4\sqrt{7R}} Y_{1}(\theta),
% \\
% \phi_{\beta,2}(\theta)
% & =\frac{5\sqrt{-3+2R}}{4\sqrt{7R}} Y_{0}(\theta)
% +\frac{\sqrt{75+62R}}{4\sqrt{7R}} Y_{1}(\theta).
%\end{align}
%\end{subequations}
%
%In the UHF formalism the Hamiltonian $\hH(\lambda)$ is no more block diagonal, $\psi_4$ can interact with $\psi_2$ and $\psi_3$. The matrix elements of the Hamiltonian corresponding to this interaction are
%\begin{equation}\label{eq:MatrixElem}
% H_{24}=H_{34}=H_{42}=H_{43}=\sqrt{R-\frac{3}{2}}\sqrt{R+\frac{75}{62}}\qty(R+\frac{25}{2})\frac{\sqrt{31}}{70R^3}
%\end{equation}
%
%For $R=3/2$ the Hamiltonian is block diagonal because the matrix elements \eqref{eq:MatrixElem} are equal to zero so this is equivalent to the RHF case but for $R>3/2$ the matrix elements become real. This interaction corresponds to the spin contamination of the wave function. For $R<3/2$ the matrix elements are complex, this corresponds to the holomorphic solution of Fig.~\ref{fig:SpheriumNrj}, the singularities in this case will be treated later. The matrix elements become real again for $R<-75/62$, this corresponds to the sb-UHF solution for negative value of $R$ observed in Sec.~\ref{sec:spherium}. We will refer to the domain where the matrix elements are complex as the holomorphic domain.
%
%The singularity structure in this case is more complex because of the spin contamination of the wave function. We can not use CSFs in this case. So when one compute all the degeneracies using Eqs.~\eqref{eq:PolChar} and \eqref{eq:DPolChar} some correspond to EPs and some correspond to conical intersections. The numerical distinction of those singularities is very difficult. We will first look at the energies $E(\lambda)$ obtained with this basis set to attribute a physical signification to the singularities obtained numerically.
%Figure \ref{fig:UHFMiniBas} is the analog of Fig.~\ref{fig:RHFMiniBas} in the UHF formalism. We see that in this case the $sp_{z}$ triplet interacts with the $s^{2}$ and the $p_{z}^{2}$ singlets. Those avoided crossings are due to the spin contamination of the wave function.
%
%\begin{figure}
% \includegraphics[width=\linewidth]{EMP_UHF_R10.pdf}
% \caption{Energies $E(\lambda)$ in the unrestricted basis set \eqref{eq:uhfbasis} with $R=10$.}
% \label{fig:UHFMiniBas}
%\end{figure}
%
%Within the RHF formalism, we have observed only $\alpha$ singularities and large avoided crossings but one can see in Fig.~\ref{fig:UHFMiniBas} that in the UHF case there are sharp avoided crossings which are known to be connected to $\beta$ singularities. For example at $R=10$ the pair of singularities connected to the avoided crossing between $s^{2}$ and $sp_{z}$ $^{3}P$ is $0.999\pm0.014\,i$. And the one between $sp_{z}$ $^{3}P$ and $p_{z}^{2}$ is connected with the singularities $2.207\pm0.023\,i$. However, in spherium, the electrons cannot be ionized so those singularities cannot be the same $\beta$ singularities as the ones highlighted by Sergeev and Goodson. \cite{Sergeev_2005} We can see in Fig.~\ref{fig:UHFEP} that the $s^{2}$ singlet and the $sp_{z}$ triplet states are degenerated for $R=3/2$. For $R>3/2$, it becomes an avoided crossing on the real axis and the degeneracies are $moved$ in the complex plane. The wave function is spin contaminated by $Y_1$ for $R>3/2$ this is why the $s^{2}$ singlet energy cannot cross the $sp_{z}$ triplet curves anymore. When $R$ increases this avoided crossing becomes sharper. As presented before $\beta$ singularities are linked to quantum phase transition so it seems that this singularity is linked to the spin symmetry breaking of the UHF wave function. The fact that a similar pair of $\beta$ singularities appears for $R<-75/62$ confirms this assumption. The sharp avoided crossing between $sp_{z}$ $^{3}P$ and $p_{z}^{2}$ is not present on Fig. \ref{fig:UHFEP}. The second pair of $\beta$ singularities resulting from this avoided crossing appears for $R\gtrsim 2.5$, this is probably due to an excited-state quantum phase transition but this still need to be investigated.
%
%\begin{figure}
% \includegraphics[width=0.45\textwidth]{UHFCI.pdf}
% \includegraphics[width=0.45\textwidth]{UHFEP.pdf}
% \caption{Energies $E(\lambda)$ in the unrestricted basis set \eqref{eq:uhfbasis} for $R=1.5$ (left) and $R=1.51$ (right).}
% \label{fig:UHFEP}
%\end{figure}
%
%As shown before, some matrix elements of the Hamiltonian become complex in the holomorphic domain. Therefore the Hamiltonian becomes non-Hermitian for these values of $R$. In Ref.~\onlinecite{Burton_2019a}, Burton \textit{et al.}~proved that for the \ce{H_2} molecule the unrestricted Hamiltonian is not \pt -symmetric in the holomorphic domain. An analog reasoning can be done with the spherium model to prove the same result. The \pt -symmetry (invariance with respect to combined space reflection $\mathcal{P}$ and time reversal $\mathcal{T}$) is a property which ensures that a non-Hermitian Hamiltonian has a real energy spectrum. \cite{BenderPTBook} Thus \pt -symmetric Hamiltonians can be seen as an intermediate class between Hermitian and non-Hermitian Hamiltonians.
%
%Figure \ref{fig:UHFPT} shows that for the spherium model a part of the energy spectrum becomes complex when $R$ is in the holomorphic domain. The domain of values where the energy becomes complex is called the broken \pt-symmetry region. This is consistent with the fact that in the holomorphic domain the Hamiltonian is no more \pt -symmetric.
%
%For a non-Hermitian Hamiltonian the EPs can lie on the real axis. In particular, at the point of {\pt} transition (the point where the energies become complex) the two energies are degenerate resulting in such an EP on the real axis. This degeneracy can be seen in Fig.~\ref{fig:UHFPT}.
%
%\begin{figure}[h!]
% \includegraphics[width=0.45\textwidth]{ReNRJPT.pdf}
% \includegraphics[width=0.45\textwidth]{ImNRJPT.pdf}
% \caption{Real part (left) and imaginary part (right) of $E(\lambda)$ in the unrestricted basis set \eqref{eq:uhfbasis} for $R=1$.}
% \label{fig:UHFPT}
%\end{figure}
\section{Conclusion}
@ -743,13 +719,13 @@ The $\beta$ singularities, which constitutes the second family, are artifacts ge
These singularities are close to the real axis and connected with sharp avoided crossing between the ground state and a highly diffuse state.
We have found that the $\beta$ singularities modeling the ionization phenomenon described by Sergeev and Goodson are actually part of a more general class of singularities. Indeed, those singularities close to the real axis are connected to quantum phase transition and symmetry breaking, and theoretical physics have demonstrated that the behavior of the EPs depends of the type of transitions from which the EPs result (first or higher orders, ground state or excited state transitions).
In this work, we have shown that $\beta$ singularities are involved in the spin symmetry breaking of the UHF wave function.
This confirms that $\beta$ singularities can occur for other types of transition and symmetry breaking than just the formation of a bound cluster of electrons.
It would be interesting to investigate the difference between the different type of symmetry breaking and how it affects the singularity structure.
Moreover the singularity structure in the non-Hermitian case still need to be investigated.
In the holomorphic domain, some singularities lie on the real axis and it would also be interesting to look at the differences between the different symmetry breaking and their respective holomorphic domain.
Furthermore, in this study we have used spherical harmonics (or combination of spherical harmonics) as basis functions which have a delocalized nature. It would also be interesting to investigate the use of localized basis functions \cite{Seidl_2018} (for example gaussians) because these functions would be more adapted to describe the strongly correlated regime. %More generally, to investigate the effect of the type of basis on the physics of EPs.
To conclude, this work shows that our understanding of the singularity structure of the energy is still incomplete but we hope that it opens new perspectives for the understanding of the physics of EPs in electronic structure theory.
%In this work, we have shown that $\beta$ singularities are involved in the spin symmetry breaking of the UHF wave function.
%This confirms that $\beta$ singularities can occur for other types of transition and symmetry breaking than just the formation of a bound cluster of electrons.
%It would be interesting to investigate the difference between the different type of symmetry breaking and how it affects the singularity structure.
%Moreover the singularity structure in the non-Hermitian case still need to be investigated.
%In the holomorphic domain, some singularities lie on the real axis and it would also be interesting to look at the differences between the different symmetry breaking and their respective holomorphic domain.
%Furthermore, in this study we have used spherical harmonics (or combination of spherical harmonics) as basis functions which have a delocalized nature. It would also be interesting to investigate the use of localized basis functions \cite{Seidl_2018} (for example gaussians) because these functions would be more adapted to describe the strongly correlated regime. %More generally, to investigate the effect of the type of basis on the physics of EPs.
%To conclude, this work shows that our understanding of the singularity structure of the energy is still incomplete but we hope that it opens new perspectives for the understanding of the physics of EPs in electronic structure theory.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

BIN
References/Altland_1997.pdf Normal file
View File

Binary file not shown.

BIN
References/Anandan_1990.pdf Normal file
View File

Binary file not shown.

BIN
References/Baker_1971.pdf Normal file
View File

Binary file not shown.

BIN
References/Benda_2017.pdf Normal file
View File

Binary file not shown.

BIN
References/Benda_2018.pdf Normal file
View File

Binary file not shown.

BIN
References/BenderBook.pdf Normal file
View File

Binary file not shown.

BIN
References/Bender_1998.pdf Normal file
View File

Binary file not shown.

BIN
References/Bender_1999.pdf Normal file
View File

Binary file not shown.

BIN
References/Bender_2002a.pdf Normal file
View File

Binary file not shown.

BIN
References/Bender_2002b.pdf Normal file
View File

Binary file not shown.

BIN
References/Bender_2003.pdf Normal file
View File

Binary file not shown.

BIN
References/Bender_2004.pdf Normal file
View File

Binary file not shown.

BIN
References/Bender_2005.pdf Normal file
View File

Binary file not shown.

BIN
References/Bender_2006.pdf Normal file
View File

Binary file not shown.

BIN
References/Bender_2007a.pdf Normal file
View File

Binary file not shown.

BIN
References/Bender_2007b.pdf Normal file
View File

Binary file not shown.

BIN
References/Bender_2008.pdf Normal file
View File

Binary file not shown.

BIN
References/Bender_2008a.pdf Normal file
View File

Binary file not shown.

BIN
References/Bender_2008b.pdf Normal file
View File

Binary file not shown.

BIN
References/Bender_2011.pdf Normal file
View File

Binary file not shown.

BIN
References/Bender_2014.pdf Normal file
View File

Binary file not shown.

BIN
References/Bender_2015.pdf Normal file
View File

Binary file not shown.

BIN
References/Bender_2016.pdf Normal file
View File

Binary file not shown.

BIN
References/Bender_2017.pdf Normal file
View File

Binary file not shown.

BIN
References/Bender_2017a.pdf Normal file
View File

Binary file not shown.

BIN
References/Bender_2017b.pdf Normal file
View File

Binary file not shown.

BIN
References/Bender_2018.pdf Normal file
View File

Binary file not shown.

BIN
References/Bender_slides.pdf Normal file
View File

Binary file not shown.

BIN
References/Berry_1984.pdf Normal file
View File

Binary file not shown.

BIN
References/Berry_2004.pdf Normal file
View File

Binary file not shown.

BIN
References/Berry_2011.pdf Normal file
View File

Binary file not shown.

BIN
References/Beygi_2018.pdf Normal file
View File

Binary file not shown.

BIN
References/Beygi_2018a.pdf Normal file
View File

Binary file not shown.

BIN
References/Beygi_2018b.pdf Normal file
View File

Binary file not shown.

BIN
References/Brody_2003.pdf Normal file
View File

Binary file not shown.

BIN
References/Burton_2015.pdf Normal file
View File

Binary file not shown.

BIN
References/Burton_2018.pdf Normal file
View File

Binary file not shown.

BIN
References/Burton_2019.pdf Normal file
View File

Binary file not shown.

BIN
References/Carlini_2006.pdf Normal file
View File

Binary file not shown.

BIN
References/Cederbaum_2004.pdf Normal file
View File

Binary file not shown.

BIN
References/Cejnar_2000.pdf Normal file
View File

Binary file not shown.

BIN
References/Cejnar_2005.pdf Normal file
View File

Binary file not shown.

BIN
References/Cejnar_2007.pdf Normal file
View File

Binary file not shown.

BIN
References/Cejnar_2008.pdf Normal file
View File

Binary file not shown.

BIN
References/Cejnar_2009.pdf Normal file
View File

Binary file not shown.

BIN
References/Cejnar_2015.pdf Normal file
View File

Binary file not shown.

BIN
References/Chaudhuri_1997.pdf Normal file
View File

Binary file not shown.

BIN
References/Cherbal_2012.pdf Normal file
View File

Binary file not shown.

BIN
References/Dasarathy_2013.pdf Normal file
View File

Binary file not shown.

BIN
References/Doppler_2016.pdf Normal file
View File

Binary file not shown.

BIN
References/Ernzerhof_2020.pdf Normal file
View File

Binary file not shown.

7574
References/Estrada_1986.pdf Normal file
View File

File diff suppressed because one or more lines are too long

BIN
References/Feuerbacher_2004.pdf Normal file
View File

Binary file not shown.

BIN
References/Finley_2000.pdf Normal file
View File

Binary file not shown.

BIN
References/Frydryszak_2008.pdf Normal file
View File

Binary file not shown.

BIN
References/Fukutome_1981.pdf Normal file
View File

Binary file not shown.

BIN
References/Gill_1986.pdf Normal file
View File

Binary file not shown.

11854
References/Gill_1988.pdf Normal file
View File

File diff suppressed because it is too large Load Diff

BIN
References/Goldzak_2018.pdf Normal file
View File

Binary file not shown.

BIN
References/Goodson_2002.pdf Normal file
View File

Binary file not shown.

BIN
References/Goodson_2004.pdf Normal file
View File

Binary file not shown.

BIN
References/Goodson_2011.pdf Normal file
View File

Binary file not shown.

BIN
References/Handy_1985.pdf Normal file
View File

Binary file not shown.

BIN
References/Hassan_2017.pdf Normal file
View File

Binary file not shown.

BIN
References/Hatano_1996.pdf Normal file
View File

Binary file not shown.

BIN
References/Hatano_1997.pdf Normal file
View File

Binary file not shown.

BIN
References/Hatano_1998.pdf Normal file
View File

Binary file not shown.

BIN
References/Heiss_1971.pdf Normal file
View File

Binary file not shown.

BIN
References/Heiss_1990.pdf Normal file
View File

Binary file not shown.

BIN
References/Heiss_1991.pdf Normal file
View File

Binary file not shown.

BIN
References/Heiss_1999.pdf Normal file
View File

Binary file not shown.

BIN
References/Heiss_2001.pdf Normal file
View File

Binary file not shown.

BIN
References/Heiss_2005.pdf Normal file
View File

Binary file not shown.

BIN
References/Heiss_2012.pdf Normal file
View File

Binary file not shown.

BIN
References/Heiss_2016.pdf Normal file
View File

Binary file not shown.

BIN
References/Henderson_2018.pdf Normal file
View File

Binary file not shown.

BIN
References/Hernandez_2006.pdf Normal file
View File

Binary file not shown.

BIN
References/Hiscock_2014.pdf Normal file
View File

Binary file not shown.

BIN
References/Jake_2018.pdf Normal file
View File

Binary file not shown.

BIN
References/Jimenez-Hoyos_2011.pdf Normal file
View File

Binary file not shown.

BIN
References/Jimenez-Hoyos_2012.pdf Normal file
View File

Binary file not shown.

BIN
References/Jones-Smith_2010.pdf Normal file
View File

Binary file not shown.

BIN
References/Jones-Smith_2013.pdf Normal file
View File

Binary file not shown.

BIN
References/Jones-Smith_2014.pdf Normal file
View File

Binary file not shown.

BIN
References/Jones-Smith_PhD.pdf Normal file
View File

Binary file not shown.

BIN
References/Kats_2018.pdf Normal file
View File

Binary file not shown.

BIN
References/Katz_1962.pdf Normal file
View File

Binary file not shown.

BIN
References/Kawabata_2019(SI).pdf Normal file
View File

Binary file not shown.

8977
References/Kawabata_2019.pdf Normal file
View File

File diff suppressed because one or more lines are too long

BIN
References/Kjonstad_2017a.pdf Normal file
View File

Binary file not shown.

BIN
References/Kjonstad_2017b.pdf Normal file
View File

Binary file not shown.

Some files were not shown because too many files have changed in this diff Show More