update

MidTermSummary.tex

@@ -20,7 +20,7 @@
 \section{Generalities}
 
 \subsection{Objectives}
-Our goal is to understand the physics of exceptional points (EPs) and perturbation series, i.e., where are EPs localize in the complex plane and how do they affect the convergence properties of the perturbation series. To do this we use the electronic ground-state of the spherium model (i.e., two opposite-spin electrons restricted to remain on a surface of a sphere of radius $R$) as a playground.
+Our goal is to understand the physics of exceptional points (EPs) and perturbation series, i.e., where are EPs localize in the complex plane and how do they affect the convergence properties of the perturbation series. To do this we use the electronic ground-state of the spherium model (i.e., two opposite-spin electrons restricted to remain on a surface of a sphere of radius $R$) as a playground. 
 
 \subsection{Variables}
 

RapportStage/Rapport.aux

@@ -84,8 +84,10 @@
 \@writefile{toc}{\contentsline {subsection}{\numberline {3.3}The singularity structure}{7}{subsection.3.3}\protected@file@percent }
 \citation{Goodson_2004}
 \citation{Olsen_2000}
+\citation{SzaboBook}
 \@writefile{toc}{\contentsline {subsection}{\numberline {3.4}The physics of quantum phase transition}{8}{subsection.3.4}\protected@file@percent }
 \@writefile{toc}{\contentsline {section}{\numberline {4}The spherium model}{8}{section.4}\protected@file@percent }
+\@writefile{toc}{\contentsline {subsection}{\numberline {4.1}Weak correlation regime}{8}{subsection.4.1}\protected@file@percent }
 \bibstyle{unsrt}
 \bibdata{Rapport}
 \bibcite{Bittner_2012}{1}

RapportStage/Rapport.log

@@ -1,4 +1,4 @@
-This is pdfTeX, Version 3.14159265-2.6-1.40.20 (TeX Live 2019/Arch Linux) (preloaded format=pdflatex 2020.5.9)  13 JUL 2020 10:12
+This is pdfTeX, Version 3.14159265-2.6-1.40.20 (TeX Live 2019/Arch Linux) (preloaded format=pdflatex 2020.5.9)  13 JUL 2020 15:18
 entering extended mode
  restricted \write18 enabled.
  %&-line parsing enabled.
@@ -754,7 +754,7 @@ Package epstopdf Info: Source file: <logoens.eps>
 df logoens.eps>
 (epstopdf)             \includegraphics on input line 72.
 Package epstopdf Info: Output file is already uptodate.
-<logoens-eps-converted-to.pdf, id=36, 263.98625pt x 97.36375pt>
+<logoens-eps-converted-to.pdf, id=40, 263.98625pt x 97.36375pt>
 File: logoens-eps-converted-to.pdf Graphic file (type pdf)
 <use logoens-eps-converted-to.pdf>
 Package pdftex.def Info: logoens-eps-converted-to.pdf  used on input line 72.
@@ -769,7 +769,7 @@ Package epstopdf Info: Source file: <logoucbl.eps>
 pdf logoucbl.eps>
 (epstopdf)             \includegraphics on input line 72.
 Package epstopdf Info: Output file is already uptodate.
-<logoucbl-eps-converted-to.pdf, id=37, 185.69376pt x 134.5025pt>
+<logoucbl-eps-converted-to.pdf, id=41, 185.69376pt x 134.5025pt>
 File: logoucbl-eps-converted-to.pdf Graphic file (type pdf)
 <use logoucbl-eps-converted-to.pdf>
 Package pdftex.def Info: logoucbl-eps-converted-to.pdf  used on input line 72.
@@ -784,7 +784,7 @@ Package epstopdf Info: Source file: <logounivlyon.eps>
 -to.pdf logounivlyon.eps>
 (epstopdf)             \includegraphics on input line 72.
 Package epstopdf Info: Output file is already uptodate.
-<logounivlyon-eps-converted-to.pdf, id=38, 301.125pt x 135.50626pt>
+<logounivlyon-eps-converted-to.pdf, id=42, 301.125pt x 135.50626pt>
 File: logounivlyon-eps-converted-to.pdf Graphic file (type pdf)
 <use logounivlyon-eps-converted-to.pdf>
 Package pdftex.def Info: logounivlyon-eps-converted-to.pdf  used on input line 
@@ -796,7 +796,7 @@ LaTeX Font Info:    Trying to load font information for T1+cmtt on input line 1
 (/usr/share/texmf-dist/tex/latex/base/t1cmtt.fd
 File: t1cmtt.fd 2014/09/29 v2.5h Standard LaTeX font definitions
 )
-<LCPQ_logo.pdf, id=41, 370.38374pt x 168.63pt>
+<LCPQ_logo.pdf, id=45, 370.38374pt x 168.63pt>
 File: LCPQ_logo.pdf Graphic file (type pdf)
 <use LCPQ_logo.pdf>
 Package pdftex.def Info: LCPQ_logo.pdf  used on input line 122.
@@ -811,13 +811,13 @@ Package epstopdf Info: Source file: <LogoCNRS.eps>
 pdf LogoCNRS.eps>
 (epstopdf)             \includegraphics on input line 122.
 Package epstopdf Info: Output file is already uptodate.
-<LogoCNRS-eps-converted-to.pdf, id=42, 57.21375pt x 57.21375pt>
+<LogoCNRS-eps-converted-to.pdf, id=46, 57.21375pt x 57.21375pt>
 File: LogoCNRS-eps-converted-to.pdf Graphic file (type pdf)
 <use LogoCNRS-eps-converted-to.pdf>
 Package pdftex.def Info: LogoCNRS-eps-converted-to.pdf  used on input line 122.
 
 (pdftex.def)             Requested size: 56.90804pt x 56.9055pt.
-<UPS_logo.jpg, id=43, 592.1322pt x 199.4652pt>
+<UPS_logo.jpg, id=47, 592.1322pt x 199.4652pt>
 File: UPS_logo.jpg Graphic file (type jpg)
 <use UPS_logo.jpg>
 Package pdftex.def Info: UPS_logo.jpg  used on input line 122.
@@ -832,7 +832,7 @@ Package pdftex.def Info: UPS_logo.jpg  used on input line 122.
 \openout4 = `Rapport.toc'.
 
  [2] [3] [4]
-<gill1986.png, id=188, 328.22626pt x 329.23pt>
+<gill1986.png, id=193, 328.22626pt x 329.23pt>
 File: gill1986.png Graphic file (type png)
 <use gill1986.png>
 Package pdftex.def Info: gill1986.png  used on input line 212.
@@ -842,44 +842,53 @@ Package pdftex.def Info: gill1986.png  used on input line 212.
 LaTeX Warning: `!h' float specifier changed to `!ht'.
 
 [5 <./gill1986.png>]
-<Nedivergence.png, id=209, 329.23pt x 234.8775pt>
+<Nedivergence.png, id=214, 329.23pt x 234.8775pt>
 File: Nedivergence.png Graphic file (type png)
 <use Nedivergence.png>
 Package pdftex.def Info: Nedivergence.png  used on input line 246.
 (pdftex.def)             Requested size: 219.35216pt x 156.49014pt.
-<HFdivergence.png, id=210, 322.20375pt x 231.86626pt>
+<HFdivergence.png, id=215, 322.20375pt x 231.86626pt>
 File: HFdivergence.png Graphic file (type png)
 <use HFdivergence.png>
 Package pdftex.def Info: HFdivergence.png  used on input line 248.
 (pdftex.def)             Requested size: 217.4634pt x 156.49014pt.
- [6 <./Nedivergence.png> <./HFdivergence.png>] [7] [8]
+ [6 <./Nedivergence.png> <./HFdivergence.png>] [7]
-(./Rapport.bbl [9] [10])
+LaTeX Font Info:    Trying to load font information for OMS+cmr on input line 3
-Package atveryend Info: Empty hook `BeforeClearDocument' on input line 311.
+06.
- [11]
+
-Package atveryend Info: Empty hook `AfterLastShipout' on input line 311.
+(/usr/share/texmf-dist/tex/latex/base/omscmr.fd
+File: omscmr.fd 2014/09/29 v2.5h Standard LaTeX font definitions
+)
+LaTeX Font Info:    Font shape `OMS/cmr/m/n' in size <10.95> not available
+(Font)              Font shape `OMS/cmsy/m/n' tried instead on input line 306.
+ [8] [9] (./Rapport.bbl
+[10] [11])
+Package atveryend Info: Empty hook `BeforeClearDocument' on input line 327.
+ [12]
+Package atveryend Info: Empty hook `AfterLastShipout' on input line 327.
  (./Rapport.aux)
-Package atveryend Info: Executing hook `AtVeryEndDocument' on input line 311.
+Package atveryend Info: Executing hook `AtVeryEndDocument' on input line 327.
-Package atveryend Info: Executing hook `AtEndAfterFileList' on input line 311.
+Package atveryend Info: Executing hook `AtEndAfterFileList' on input line 327.
 Package rerunfilecheck Info: File `Rapport.out' has not changed.
-(rerunfilecheck)             Checksum: 88987B6FDBFC32AB6DEBEB04F7925937;527.
+(rerunfilecheck)             Checksum: 0A559EFFB3496CDF4124AA3173D1DDF8;599.
 
 
 LaTeX Warning: There were multiply-defined labels.
 
-Package atveryend Info: Empty hook `AtVeryVeryEnd' on input line 311.
+Package atveryend Info: Empty hook `AtVeryVeryEnd' on input line 327.
  )
 (\end occurred inside a group at level 1)
 
 ### simple group (level 1) entered at line 207 ({)
 ### bottom level 
 Here is how much of TeX's memory you used:
- 21723 strings out of 492167
+ 21744 strings out of 492167
- 398015 string characters out of 6131559
+ 398339 string characters out of 6131559
  493579 words of memory out of 5000000
- 25816 multiletter control sequences out of 15000+600000
+ 25832 multiletter control sequences out of 15000+600000
  557560 words of font info for 90 fonts, out of 8000000 for 9000
  1141 hyphenation exceptions out of 8191
- 52i,12n,71p,1636b,994s stack positions out of 5000i,500n,10000p,200000b,80000s
+ 52i,12n,71p,1615b,994s stack positions out of 5000i,500n,10000p,200000b,80000s
 {/usr/share/texmf-dist/fonts/enc/dvips/cm-super/cm-super-t1.enc
 }</usr/share/texmf-dist/fonts/type1/public/amsfonts/cm/cmbsy10.pfb></usr/share/
 texmf-dist/fonts/type1/public/amsfonts/cm/cmbx10.pfb></usr/share/texmf-dist/fon
@@ -901,10 +910,10 @@ cm-super/sfrm1095.pfb></usr/share/texmf-dist/fonts/type1/public/cm-super/sfti10
 00.pfb></usr/share/texmf-dist/fonts/type1/public/cm-super/sfti1095.pfb></usr/sh
 are/texmf-dist/fonts/type1/public/cm-super/sftt1095.pfb></usr/share/texmf-dist/
 fonts/type1/public/cm-super/sfxc1095.pfb>
-Output written on Rapport.pdf (11 pages, 1005081 bytes).
+Output written on Rapport.pdf (12 pages, 1008118 bytes).
 PDF statistics:
- 353 PDF objects out of 1000 (max. 8388607)
+ 365 PDF objects out of 1000 (max. 8388607)
- 298 compressed objects within 3 object streams
+ 308 compressed objects within 4 object streams
- 78 named destinations out of 1000 (max. 500000)
+ 80 named destinations out of 1000 (max. 500000)
- 110 words of extra memory for PDF output out of 10000 (max. 10000000)
+ 118 words of extra memory for PDF output out of 10000 (max. 10000000)
 

RapportStage/Rapport.out

@@ -6,3 +6,4 @@
 \BOOKMARK [2][-]{subsection.3.3}{The singularity structure}{section.3}% 6
 \BOOKMARK [2][-]{subsection.3.4}{The physics of quantum phase transition}{section.3}% 7
 \BOOKMARK [1][-]{section.4}{The spherium model}{}% 8
+\BOOKMARK [2][-]{subsection.4.1}{Weak correlation regime}{section.4}% 9

RapportStage/Rapport.tex

@@ -91,7 +91,7 @@ hyperfigures=false]
 \vspace{1.5cm}
 
 \parbox{15cm}{\small
-\textbf{Abstract} : \it In this work, we explore the description of quantum chemistry in the complex plane. We see that the physic of the system can be connected to the position of the singularities of the energy in the complex plane. After briefly present the fundamental notions of quantum chemistry and perturbation theory in the complex plane, we perform an historical review of the researches that have been done on the physic of singularities. Then we make links between all those points of view on this problem using the spherium model (i.e., two opposite-spin electrons restricted to remain on a surface of a sphere of radius $R$) as a theoretical playground. In particular, we explore the effects of symmetry breaking of the wave functions on the singularity structure.
+\textbf{Abstract} : \it In this work, we explore the description of quantum chemistry in the complex plane. We see that the physics of the system can be 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 perform an historical review of the researches that have been done on the physic of singularities. Then we connect all those points of view on this problem using the spherium model (i.e., two opposite-spin electrons restricted to remain on the surface of a sphere of radius $R$) as a theoretical playground. In particular, we explore the effects of symmetry breaking of the wave functions on the singularity structure.
 }
 
 \vspace{0.5cm}
@@ -138,26 +138,26 @@ Laboratoire de Chimie et Physique Quantiques
 \section{Introduction}
 %============================================================%
 
-It has always been of great importance for theoretical chemists to better understand excited states and their properties because processes involving excited states are ubiquitous in nature (physics, chemistry and biology). One of the major challenges is to accurately compute energies of a chemical system (atoms, molecules, ..). Plenty of methods have been developed to this aim and each of them have its own qualities but also its own flaws. The fact that none of all those methods is successful for every molecule in every geometry encourage chemists to continue the development of new methodologies to get accurate energies and to try to understand deeply why each methods fails or not in each situation. All those methods rely on the notion of quantised energy levels of Hermitian quantum mechanics. In quantum chemistry, the ordering of the energy levels represents the different electronic states of a molecule, the lowest being the ground state while the higher ones are the so-called excited states. And we need methods to get accurately how those states are ordered.
+It has always been of great importance for theoretical chemists to better understand excited states and their properties because processes involving excited states are ubiquitous in nature (physics, chemistry and biology). One of the major challenges is to accurately compute energies of a chemical system (atoms, molecules, ..). Plenty of methods have been developed to this end and each of them have its own advantages but also its own flaws. The fact that none of all those methods is successful for every molecule in every geometry encourages chemists to continue the development of new methodologies to get accurate energies and to try to deeply understand why each method fails or not in each situation. All those methods rely on the notion of quantised energy levels of Hermitian quantum mechanics. In quantum chemistry, the ordering of the energy levels represents the different electronic states of a molecule, the lowest being the ground state while the higher ones are the so-called excited states. We need methods to accurately get how those states are ordered.
 
 Within this quantised paradigm, electronic states look completely disconnected from one another.
-However, one can gain a different perspective on quantisation if one extends quantum chemistry into the complex domain.
+However, one can gain a different perspective on quantisation extending quantum chemistry into the complex domain.
 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 quantised nature of conventional Hermitian quantum mechanics arises only from our limited perception of the more complex and profound structure of its non-Hermitian variant.
 
-Therefore, by analytically continuing the energy $E(\lambda)$ in the complex domain (where $\lambda$ is a coupling parameter), one can smoothly connect the ground and excited states of a molecule.
+Therefore, by analytically continuing the 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, one loses the ordering property of real numbers.
+This connection is possible because by extending real numbers to the complex domain, the ordering property of real numbers is lost.
-Hence, one can interchange electronic states away from the real axis, as the concept of ground and excited states has been 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 been recently realized experimentally 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) \cite{Heiss_1990, Heiss_1999, Heiss_2012, Heiss_2016} are non-Hermitian analogs of conical intersections (CIs) \cite{Yarkony_1996} where two states become exactly degenerate.
 CIs 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 CIs in the decay of bound excited states.
 Although Hermitian and non-Hermitian Hamiltonians are closely related, the behavior of their eigenvalues near degeneracies is starkly different.
-For example, by 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.
-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 (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 CIs introduces only a geometric phase while leaving the states unchanged.
+In contrast, encircling Hermitian degeneracies at CIs only introduces a geometric phase while leaving the states unchanged.
 More dramatically, whilst eigenvectors remain orthogonal at CIs, 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{Olsen_1996, Olsen_2000} 
 
@@ -190,13 +190,13 @@ H_0= \sum\limits_{i=1}^{n} f(i)
 f(i) = h(i) + \sum\limits_{j=1,j \neq i}^{n} \left[J_j(i) - K_j(i)\right]
 \end{equation}
 
-In Hartree-Fock theory the exact wave function is approximated as a Slater-determinant (which is an anti-symmetric combination of mono-electronic orbitals) and those wave functions are eigenvectors of the Fock operators. In the perturbation theory the energy is a power series of $\lambda$ and the physical energy is obtained by taking $\lambda$ equal to 1. We will refer to there energy up to the n-th order as the MPn energy. The MP0 energy overestimate the energy by double counting the electron-electron interaction, the MP1 correct this effect and the MP1 energy is equal to the Hartree-Fock energy. The MP2 energy start to recover a part of the correlation energy.
+In Hartree-Fock theory the exact wave function is approximated as a Slater-determinant (which is an anti-symmetric combination of mono-electronic orbitals) and those wave functions are eigenvectors of the Fock operators. In the perturbation theory the energy is a power series of $\lambda$ and the physical energy is obtained by taking $\lambda$ equal to 1. We will refer to the energy up to the n-th order as the MPn energy. The MP0 energy overestimates the energy by double counting the electron-electron interaction, the MP1 corrects this effect and the MP1 energy is equal to the Hartree-Fock energy. The MP2 energy starts to recover a part of the correlation energy.
 
 \begin{equation}
 E_{\text{MP\textsubscript{n}}}= \sum_{k=0}^n E^{(k)}
 \end{equation}
 
-But as mentioned before \textit{a priori} there are no reasons that this power series is always convergent for $\lambda$=1 when n goes to infinity. In fact, it is known that when the Hartree-Fock wave function is a bad approximation of the exact wave function, for example for multi-reference states, the M{\o}ller-Plesset will give bad results\cite{Gill_1986, Gill_1988, Handy_1985, Lepetit_1988}. A smart way to investigate the convergence properties of the MP series is to transform the coupling parameter $\lambda$ into a complex variable. By doing this the Hamiltonian and the energy become functions of this variable. By doing this the energy become a multivalued function on $n$ Riemann sheets. As mentioned above by searching the singularities of the function $E(\lambda)$ we can get information on the convergence properties of the MPPT. Those singularities of the energy are exactly the exceptional points connecting the electronic states mentioned in the introduction. The direct computation of the terms of the series is relatively easy up to the 4th order and the 5th and 6th order can be obtained at high cost. But to understand deeply the behavior of the MP series and how it is connected to the singularities, we need to have access to high order terms of the series. For small systems we can have access to the whole series using Full Configuration Interaction. If you diagonalize the Hamiltonian $H(\lambda)$ in the FCI basis set you get the exact energies (in this finite basis set) and expanding in $\lambda$ allow you to get the M{\o}ller-Plesset perturbation series at every order.
+But as mentioned before \textit{a priori} there are no reasons that this power series is always convergent for $\lambda$=1 when n goes to infinity. In fact, it is known that when the Hartree-Fock wave function is a bad approximation of the exact wave function, for example for multi-reference states, the M{\o}ller-Plesset will give bad results\cite{Gill_1986, Gill_1988, Handy_1985, Lepetit_1988}. A smart way to investigate the convergence properties of the MP series is to transform the coupling parameter $\lambda$ into a complex variable. By doing so the Hamiltonian and the energy become functions of this variable. The energy becomes a multivalued function on $n$ Riemann sheets. As mentioned above by searching the singularities of the function $E(\lambda)$ we can get information on the convergence properties of the MPPT. Those singularities of the energy are exactly the exceptional points connecting the electronic states mentioned in the introduction. The direct computation of the terms of the series is quite easy up to the 4th order and the 5th and 6th order can be obtained at high cost. But to deeply understand the behavior of the MP series and how it is connected to the singularities, we need to have access to high order terms of the series. For small systems we can have access to the whole series using Full Configuration Interaction. If the Hamiltonian $H(\lambda)$ is diagonalized in the FCI basis set we get the exact energies (in this finite basis set) and expanding in $\lambda$ allows to to get the M{\o}ller-Plesset perturbation series at every order.
 
 %============================================================%
 \section{Historical overview}
@@ -204,8 +204,8 @@ But as mentioned before \textit{a priori} there are no reasons that this power s
 
 \subsection{Behavior of the M{\o}ller-Plesset series}
 
-When you use M{\o}ller-Plesset perturbation theory it would be very convenient that each time you compute a term of higher order the result obtained is closer to exact energy. In other words that the M{\o]ller-Plesset series would be monotonically convergent. Assuming this the only limiting process to get the exact correlation energy in a finite basis is our ability to compute the terms of the perturbation series.
+When we use M{\o}ller-Plesset perturbation theory it would be very convenient that each time a higher order term is computed the result obtained is closer to exact energy. In other words, that the M{\o]ller-Plesset series would be monotonically convergent. Assuming this, the only limiting process to get the exact correlation energy in a finite basis set is our ability to compute the terms of the perturbation series.
-Unfortunately this not true in generic cases and rapidly some strange behaviors of the series have been exhibited. In the late 80's Gill et al. have reported deceptive and slow convergences in stretch systems\cite{Gill_1986, Gill_1988, Handy_1985, Lepetit_1988}. In the figure below we can see that the restricted M{\o}ller-Plesset series is convergent but oscillating which is not convenient if you are able to compute only few terms (for example here RMP5 is worse than RMP4). On the other hand, the unrestricted M{\o}ller-Plesset series is monotonically converging (except for the first few orders) but very slowly so we can't use it for systems where we can compute only the first terms.
+Unfortunately this is not true in generic cases and rapidly some strange behaviors of the series were exhibited. In the late 80's Gill et al. reported deceptive and slow convergences in stretch systems\cite{Gill_1986, Gill_1988, Handy_1985, Lepetit_1988}. In the figure below we can see that the restricted M{\o}ller-Plesset series is convergent but oscillating which is not convenient if you are only able to compute few terms (for example here RMP5 is worse than RMP4). On the other hand, the unrestricted M{\o}ller-Plesset series is monotonically converging (except for the first few orders) but very slowly so we can't use it for systems where we can only compute the first terms.
 
 \begin{figure}[h!]
     \centering
@@ -214,9 +214,9 @@ Unfortunately this not true in generic cases and rapidly some strange behaviors
     \label{fig:my_label}
 \end{figure}
 
-When a bond is stretched the exact function can undergo a symmetry breaking becoming multi-reference during this process (see for example the case of \ce{H_2} in \cite{SzaboBook}). A restricted HF Slater determinant is a poor approximation of a broken symmetry wave function but even in the unrestricted formalism, where the spatial orbitals of electrons $\alpha$ and $\beta$ are not restricted to be the same\cite{Fukutome_1981}, which allow a better description of symmetry broken system the series doesn't give accurate results at low orders. Even with this improvement of the zeroth order wave function the series doesn't have the smooth and rapidly converging behavior wanted. 
+When a bond is stretched the exact function can undergo a symmetry breaking becoming multi-reference during this process (see for example the case of \ce{H_2} in \cite{SzaboBook}). A restricted HF Slater determinant is a poor approximation of a broken symmetry wave function but even in the unrestricted formalism, where the spatial orbitals of electrons $\alpha$ and $\beta$ are not restricted to be the same\cite{Fukutome_1981}, which allows a better description of broken symmetry system, the series doesn't give accurate results at low orders. Even with this improvement of the zeroth order wave function the series doesn't have the smooth and rapidly converging behavior wanted. 
 
-In the unrestricted framework the ground state singlet wave function is allowed to mix with triplet states which lead to spin contamination. Gill et al. highlighted that there is a link between the slow convergence of the unrestricted MP series and the spin contamination of the wave function as it's shown in the Table 1 within the example of \ce{H_2} in a minimal basis.
+In the unrestricted framework the ground state singlet wave function is allowed to mix with triplet states which leads to spin contamination. Gill et al. highlighted the link between the slow convergence of the unrestricted MP series and the spin contamination of the wave function as shown in the Table 1 in the example of \ce{H_2} in a minimal basis.
 
 \begin{table}[h!]
     \centering
@@ -234,13 +234,13 @@ In the unrestricted framework the ground state singlet wave function is allowed
     \label{tab:my_label}
 \end{table}
 
-Handy and co-workers. exhibited the same behaviors of the series (oscillating and monotonically slowly) in stretched \ce{H_2O} and \ce{NH_2} systems \cite{Handy_1985}. Cremer and He performed the same analysis with 29 FCI systems \cite{Cremer_1996} and regrouped all the systems in two classes. The class A systems which have a monotonic convergence to the FCI value and the class B which converge erratically after initial oscillations. The sample of systems contains stretched molecules and also some at equilibrium geometry, there are also some systems in various basis sets. They highlighted that systems with class A convergence have well-separated electrons pairs whereas class B systems present electrons clustering.
+Handy and co-workers exhibited the same behaviors of the series (oscillating and monotonically slowly) in stretched \ce{H_2O} and \ce{NH_2} systems \cite{Handy_1985}. Cremer and He performed the same analysis with 29 FCI systems \cite{Cremer_1996} and regrouped all the systems in two classes. The class A systems which have a monotonic convergence to the FCI value and the class B which converge erratically after initial oscillations. The sample of systems contains stretched molecules and also some at equilibrium geometry, there are also some systems in various basis sets. They highlighted that systems with class A convergence have well-separated electrons pairs whereas class B systems present electrons clustering.
 
-This classification was encouraging in order to develop methods based on perturbation theory as it rationalize the two different convergence modes observed. If it is possible to predict if a system is class A or B, then one can use extrapolation method of the first terms adapted to the class of the systems \cite{Cremer_1996}.
+This classification was encouraging in order to develop methods based on perturbation theory as it rationalizes the two different observed convergence modes. If it is possible to predict if a system is class A or B, then one can use extrapolation method of the first terms adapted to the class of the systems \cite{Cremer_1996}.
 
 \subsection{Cases of divergence}
 
-But Olsen et al. have discovered even more preoccupying behavior of the MP series in the late 90's. They have shown 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 those two systems and classified them as class B systems. But Olsen and his co-workers have done the analysis in larger basis sets containing diffuse functions and in those basis sets the series become divergent at high order.
+However Olsen et al. have discovered an even more preoccupying behavior of the MP series in the late 90's. They have shown 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 those two systems and classified them as class B systems. But Olsen and his co-workers have done the analysis in larger basis sets containing diffuse functions and in those basis sets the series become divergent at high order.
 
 \begin{figure}[h!]
     \includegraphics[height=5.5cm]{Nedivergence.png}
@@ -251,25 +251,25 @@ But Olsen et al. have discovered even more preoccupying behavior of the MP serie
     \label{fig:my_label}
 \end{figure}
 
-The discovery of those divergent behavior was really worrying because in order to get more and more accurate results theoretical chemists need to work in large basis sets. So they investigated the causes of those divergences and in the same time the reasons of the different types of convergence. In order to do this they analyzed the relation between the dominant singularity (i.e. the closest singularity to the origin) and the convergence behavior of the series \cite{Olsen_2000}. Their analysis is based on Darboux's theorem: 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, one can explain the convergence patterns of the MP series by looking at the dominant singularity
+The discovery of this divergent behavior was really worrying because in order to get more and more accurate results theoretical chemists need to work in large basis sets. As a consequence they investigated the causes of those divergences and in the same time the reasons of the different types of convergence. To do this they analyzed the relation between the dominant singularity (i.e. the closest singularity to the origin) and the convergence behavior of the series \cite{Olsen_2000}. Their analysis is based on Darboux's theorem: 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.
 
-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). The method used is to do a scan of the real axis to identify the avoided crossing responsible for the pair of dominant singularity. Then by modeling this avoided crossing by a two-state Hamiltonian one can get an approximation of the dominant conjugate pair of singularity by finding the EPs of the 2x2 Hamiltonian. The diagonal matrix is the unperturbed Hamiltonian and the other matrix is the perturbation part of the 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). The method used is to do a scan of the real axis to identify the avoided crossing responsible for the pair of dominant singularity. Then by modeling this avoided crossing by a two-state Hamiltonian one can get an approximation of the dominant conjugate pair of singularity by finding the EPs of the $22\times$ Hamiltonian. The diagonal matrix is the unperturbed Hamiltonian and the other matrix is the perturbative part of the Hamiltonian.
 
 \begin{equation}
 \mqty(\alpha & \delta \\ \delta & \beta) = \mqty(\alpha + \alpha_s & 0 \\ 0 & \beta + \beta_s ) + \mqty(- \alpha_s & \delta \\ \delta & - \beta_s)
 \end{equation}
  
-They first studied molecules with low-lying doubly excited states of the same spatial and spin symmetry because in those systems the HF wave function is a bad approximation. The exact wave function have non-negligible contribution from the doubly excited states so those 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 have shown that the dominant singularity lies outside the unit circle but close to it causing the slow convergence.
+They first studied molecules with low-lying doubly excited states of the same spatial and spin symmetry 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 those 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.
 
-Then they have shown that the divergence for the \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 conclude that the divergence of the series was due to the interaction with a highly diffuse excited states. 
+Then they demonstrated that the divergence for the \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 formula of Cremer and He \cite{Cremer_1996} can't be used for all systems. Even more that those formula were not mathematically motivated when one look 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 result in alternated terms up to order 10 and then the series is monotonically convergent. This is due to the fact that two pairs of singularity are approximately at the same distance of the origin.
+Moreover they proved that the extrapolation formula of Cremer and He \cite{Cremer_1996} can't be used for all systems. Even more, that those formula 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 order ten and then the series is monotonically convergent. This is due to the fact that two pairs of singularity 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 his co-workers were trying to find the dominant singularity causing the divergence. They regroup singularities in two classes: the $\alpha$ singularities which have unit order imaginary parts and the $\beta$ singularities which have very small imaginary parts. The singularities $\alpha$ are related to large avoided crossing between the ground state and a low-lying excited states. Whereas the singularities $\beta$ come from a sharp avoided crossing between the ground state and a highly diffuse state. They succeed to explain the divergence of the series caused by $\beta$ singularities using a previous work of Stillinger \cite{Stillinger_2000}. 
+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 his co-workers were trying to find the dominant singularity causing the divergence. They regrouped singularities in two classes: the $\alpha$ singularities which have unit order imaginary parts and the $\beta$ singularities which have very small imaginary parts. The singularities $\alpha$ are related to large avoided crossing between the ground state and a low-lying excited states. Whereas the singularities $\beta$ 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 using a previous work of Stillinger \cite{Stillinger_2000}. 
 
-The M{\o}ller-Plesset Hamiltonian is defined as below and by reassembling the term we get the expression (10).
+The M{\o}ller-Plesset Hamiltonian is defined as below and by reassembling the term we get the expression (11).
 
 \begin{equation}
 H(\lambda)=H_0 + \lambda (H_\text{phys} - H_0)    
@@ -287,19 +287,35 @@ H(\lambda)=H_0 + \lambda (H_\text{phys} - H_0)
     H(\lambda)=\sum\limits_{j=1}^{n}\left[-\frac{1}{2}\grad_j^2 - \sum\limits_{k=1}^{N} \frac{Z_k}{|\vb{r}_j-\vb{R}_k|} + (1-\lambda)V_j^{(scf)}+\lambda\sum\limits_{j<l}^{n}\frac{1}{|\vb{r}_j-\vb{r}_l|} \right]
 \end{equation}
 
-The first two terms, the kinetic energy and the electron-nucleus attraction, form the mono-electronic core Hamiltonian which is independant of $\lambda$. The third term is the mean field repulsion of the Hartree-Fock calculation done to get $H_0$ and the last term is the Coulomb repulsion. If $\lambda$ is negative, the Coulomb interaction becomes attractive but the mean field stay repulsive as it is proportional to $(1-\lambda)$. If $\lambda$ becomes more and more negative the mean field becomes more and more repulsive so the nucleus can't bind anymore the electrons because the electron-nucleus attraction is not scaled with $\lambda$. But the repulsive mean field is localized around nucleus whereas the electrons interactions persist away from nucleus. So there is a real negative value $\lambda_c$ where the electron form a bound cluster and goes to infinity. According to Baker this value is a critical point of the system and by analogy with thermodynamics the energy $E(\lambda)$ exhibits a singularity at $\lambda_c$. Beyond $\lambda_c$ there is a continuum of eigenstates with electrons dissociated from the nucleus.
+The first two terms, the kinetic energy and the electron-nucleus attraction, form the mono-electronic core Hamiltonian which is independant of $\lambda$. The third term is the mean field repulsion of the Hartree-Fock calculation done to get $H_0$ and the last term is the Coulomb repulsion. If $\lambda$ is negative, the Coulomb interaction becomes attractive but the mean field stays repulsive as it is proportional to $(1-\lambda)$. If $\lambda$ becomes more and more negative the mean field becomes more and more repulsive so the nucleus can't bind the electrons anymore because the electron-nucleus attraction is not scaled with $\lambda$. The repulsive mean field is localized around nucleus whereas the electrons interactions persist away from nucleus. There is a real negative value $\lambda_c$ where the electrons form a bound cluster and goes to infinity. According to Baker this value is a critical point of the system and by analogy with thermodynamics the energy $E(\lambda)$ exhibits a singularity at $\lambda_c$. Beyond $\lambda_c$ there is a continuum of eigenstates with electrons dissociated from the nucleus.
 
-This reasoning is done on the exact Hamiltonian and energy, this is the exact energy which exhibits this singularity on the negative real axis. But in finite basis set, one can prove that for a Hermitian Hamiltonian the singularities of $E(\lambda)$ occurs in complex conjugate pair with non-zero imaginary parts. Sergeev and Goodson proved that, as predicted by Stillinger, in a finite basis set the critical point on the real axis modeled by a cluster of sharp avoided crossings with diffuse functions, equivalently by a cluster of $\beta$ singularities in the negative half plane. They explain that Olsen et al., because they used a 2x2 model, only observed the first singularity of this cluster of singularities causing the divergence.
+This reasoning is done on the exact Hamiltonian and energy, this is the exact energy which exhibits this singularity on the negative real axis. But in finite basis set, one can prove that for a Hermitian Hamiltonian the singularities of $E(\lambda)$ occurs in complex conjugate pair with non-zero imaginary parts. Sergeev and Goodson proved, as predicted by Stillinger, 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. They explain that Olsen et al., because they used a $2\times2$ model, only observed the first singularity of this cluster of singularities causing the divergence.
 
-Finally, it has been shown that $\beta$ singularities are very sensitive to the basis sets but not to 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 the singularity structure from molecules stretched from the equilibrium geometry is difficult \cite{Goodson_2004}, this is conform with the observation of Olsen and co-workers on the \ce{HF} molecule at equilibrium geometry and stretched geometry \cite{Olsen_2000}. To our knowledge the effect of bond stretching on singularities, its link with spin contamination and symmetry breaking of the wave function haven't been as well understood as the ionization effect and its link with diffuse function. In this work we 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 the basis sets 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 the singularity structure from molecules stretched from the equilibrium geometry is difficult \cite{Goodson_2004}, this is consistent with the observation of Olsen and co-workers on the \ce{HF} molecule at equilibrium geometry and stretched geometry \cite{Olsen_2000}. To our knowledge the effect of bond stretching on singularities, its link with spin contamination and symmetry breaking of the wave function hasn't been as well understood as the ionization effect and its link with diffuse function. In this work we 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.
 
 \subsection{The physics of quantum phase transition}
 
+At a sharp avoided crossings the two eigenstates undergo 
+
 %============================================================%
 \section{The spherium model}
 %============================================================%
 
+Simple systems that are analytically solvable (or at least quasi-exactly solvable) are of great importance in theoretical chemistry. Those systems are very useful benchmarks to test new methods as they are mathematically easy but retain much of the key physics. To investigate the physics of EPs we use one such system named spherium model. It consists of two electrons confined to the surface of a sphere interacting through the long-range Coulomb potential. The radius R of the sphere dictates the correlation regime, i.e., weak correlation regime at small $R$ where the kinetic energy dominates, or strong correlation regime where the electron repulsion term drives the physics. We will use this model to try to rationalize the effects of the variables that may influence the physics of EPs:
+\begin{itemize}
+	\item Partitioning of the Hamiltonian and the actual zeroth-order reference: weak correlation reference [restricted Hartree-Fock (RHF) or unrestricted Hartree-Fock (UHF) references, M{\o}ller-Plesset or Epstein-Nesbet (EN) partitioning], or strongly correlated reference.
+	\item Basis set: minimal basis or infinite (i.e., complete) basis made of localized or delocalized basis functions
+	\item Radius of the spherium that ultimately dictates the correlation regime.
+\end{itemize}
 
+\subsection{Weak correlation regime}
+
+In the restricted Hartree-Fock formalism, the wave function can't model properly the physics of the system at large R because the spatial orbitals are restricted to be the same. Then a fortiori it can't represent two electrons on opposite side of the sphere. In the unrestricted formalism
+At a critical value of R, called the Coulson-Fischer point, a second unrestricted Hartree-Fock solution appear. This solution is symmetry-broken as the two electrons tends to localize on opposite side of the sphere. By analogy with the case of \ce{H_2} \cite{SzaboBook}, the unrestricted Hartree-Fock wave function is defined as:
+
+
+
+Then the mono-electronic wave function are expand in the spatial basis set of the zonal spherical harmonics:
 
 \newpage
 

RapportStage/Rapport.toc

@@ -7,3 +7,4 @@
 \contentsline {subsection}{\numberline {3.3}The singularity structure}{7}{subsection.3.3}% 
 \contentsline {subsection}{\numberline {3.4}The physics of quantum phase transition}{8}{subsection.3.4}% 
 \contentsline {section}{\numberline {4}The spherium model}{8}{section.4}% 
+\contentsline {subsection}{\numberline {4.1}Weak correlation regime}{8}{subsection.4.1}%