From 77fa6c4dcc53b8478d9e82218a38d1bbe6848704 Mon Sep 17 00:00:00 2001
From: Pierre-Francois Loos <pierrefrancois.loos@gmail.com>
Date: Tue, 24 Nov 2020 09:25:43 +0100
Subject: [PATCH 1/3] minor corrections in QPT

---
 Manuscript/EPAWTFT.bbl | 12 ++++++------
 Manuscript/EPAWTFT.tex |  7 ++++---
 2 files changed, 10 insertions(+), 9 deletions(-)

diff --git a/Manuscript/EPAWTFT.bbl b/Manuscript/EPAWTFT.bbl
index 18ecaee..0398b28 100644
--- a/Manuscript/EPAWTFT.bbl
+++ b/Manuscript/EPAWTFT.bbl
@@ -530,12 +530,6 @@
   {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 {Goodson}(2011)}]{Goodson_2011}%
-  \BibitemOpen
-  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {D.~Z.}\ \bibnamefont
-  {Goodson}},\ }\href {\doibase 10.1002/wcms.92} {\bibfield  {journal}
-  {\bibinfo  {journal} {{WIREs} Comput. Mol. Sci.}\ }\textbf {\bibinfo {volume}
-  {2}},\ \bibinfo {pages} {743} (\bibinfo {year} {2011})}\BibitemShut {NoStop}%
 \bibitem [{\citenamefont {Cejnar}\ \emph {et~al.}(2007)\citenamefont {Cejnar},
   \citenamefont {Heinze},\ and\ \citenamefont {Macek}}]{Cejnar_2007}%
   \BibitemOpen
@@ -545,6 +539,12 @@
   {\doibase 10.1103/PhysRevLett.99.100601} {\bibfield  {journal} {\bibinfo
   {journal} {Phys. Rev. Lett.}\ }\textbf {\bibinfo {volume} {99}},\ \bibinfo
   {pages} {100601} (\bibinfo {year} {2007})}\BibitemShut {NoStop}%
+\bibitem [{\citenamefont {Goodson}(2011)}]{Goodson_2011}%
+  \BibitemOpen
+  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {D.~Z.}\ \bibnamefont
+  {Goodson}},\ }\href {\doibase 10.1002/wcms.92} {\bibfield  {journal}
+  {\bibinfo  {journal} {{WIREs} Comput. Mol. Sci.}\ }\textbf {\bibinfo {volume}
+  {2}},\ \bibinfo {pages} {743} (\bibinfo {year} {2011})}\BibitemShut {NoStop}%
 \bibitem [{\citenamefont {Szabo}\ and\ \citenamefont
   {Ostlund}(1989)}]{SzaboBook}%
   \BibitemOpen
diff --git a/Manuscript/EPAWTFT.tex b/Manuscript/EPAWTFT.tex
index 21f3cad..228acf1 100644
--- a/Manuscript/EPAWTFT.tex
+++ b/Manuscript/EPAWTFT.tex
@@ -1015,7 +1015,7 @@ To the best of our knowledge, the effect of bond stretching on singularities, it
 %====================================================
 
 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. 
+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,Borisov_2015,Sindelka_2017,CarrBook,Vojta_2003,SachdevBook,GilmoreBook} 
 In quantum mechanics, the Hamiltonian is almost always dependent of, at least, one parameter. 
@@ -1023,7 +1023,7 @@ In some cases the variation of a parameter can lead to abrupt changes at a criti
 These QPTs exist both for ground and excited states as shown by Cejnar and coworkers. \cite{Cejnar_2005,Cejnar_2007,Caprio_2008,Cejnar_2009,Sachdev_2011,Cejnar_2015,Cejnar_2016, Macek_2019,Cejnar_2020} 
 A ground-state QPT is characterised 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. 
+Otherwise, it is called continuous and of $m$th order (with $m \ge 2$) 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. 
@@ -1033,7 +1033,7 @@ One of the major obstacles that one faces in order to achieve this resides in th
 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} 
+More recently Stransky and coworkers proved that the distribution of EPs is characteristic of the QPT order. \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 interacting 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.
@@ -1046,6 +1046,7 @@ For example, without interaction the ground state is the spherical phase (a cond
 In particular, we see that the transition from the spherical phase to the axially symmetric one is analog to the symmetry breaking of the wave function of the hydrogen molecule when the bond is stretched. \cite{SzaboBook}
 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 modelling 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. 

From a96a52c464998acf4dd94297dccdc1abd0343ca1 Mon Sep 17 00:00:00 2001
From: Pierre-Francois Loos <pierrefrancois.loos@gmail.com>
Date: Tue, 24 Nov 2020 09:40:33 +0100
Subject: [PATCH 2/3] add Antoine example and degreenified Hugh comments

---
 Manuscript/EPAWTFT.tex | 58 ++++++++++++++++++++----------------------
 1 file changed, 28 insertions(+), 30 deletions(-)

diff --git a/Manuscript/EPAWTFT.tex b/Manuscript/EPAWTFT.tex
index 228acf1..3532575 100644
--- a/Manuscript/EPAWTFT.tex
+++ b/Manuscript/EPAWTFT.tex
@@ -195,36 +195,36 @@ We refer the interested reader to the excellent book of Moiseyev for a general o
 \label{sec:TDSE}
 %%%%%%%%%%%%%%%%%%%%%%%
 Within the Born-Oppenheimer approximation, the exact molecular Hamiltonian with $\Ne$ electrons and 
-$\Nn$ (clamped) nuclei is defined \hugh{for a given nuclear framework} as
+$\Nn$ (clamped) nuclei is defined for a given nuclear framework as
 \begin{equation}\label{eq:ExactHamiltonian}
-    \hugh{\hH(\vb{R})} = 
+    \hH(\vb{R}) = 
     - \frac{1}{2} \sum_{i}^{\Ne} \grad_i^2 
     - \sum_{i}^{\Ne} \sum_{A}^{\Nn} \frac{Z_A}{\abs{\vb{r}_i-\vb{R}_A}} 
     + \sum_{i<j}^{\Ne}\frac{1}{\abs{\vb{r}_i-\vb{r}_j}},
 \end{equation}
 where $\vb{r}_i$ defines the position of the $i$-th electron, $\vb{R}_{A}$ and $Z_{A}$ are the position
-and charge of the $A$-th nucleus respectively, \hugh{and $\vb{R} = (\vb{R}_{1}, \dots, \vb{R}_{\Nn})$ is a
-collective vector for the nuclear positions.}
+and charge of the $A$-th nucleus respectively, and $\vb{R} = (\vb{R}_{1}, \dots, \vb{R}_{\Nn})$ is a
+collective vector for the nuclear positions.
 The first term represents the kinetic energy of the electrons, while 
 the two following terms account for the electron-nucleus attraction and the electron-electron repulsion.
 
 % EXACT SCHRODINGER EQUATION
-The exact many-electron wave function \hugh{at a given nuclear geometry} $\Psi(\vb{R})$ corresponds 
+The exact many-electron wave function at a given nuclear geometry $\Psi(\vb{R})$ corresponds 
 to the solution of the (time-independent) Schr\"{o}dinger equation
 \begin{equation} 
-    \hugh{\hH(\vb{R})\, \Psi(\vb{R}) = E(\vb{R})\, \Psi(\vb{R}),}
+    \hH(\vb{R})\, \Psi(\vb{R}) = E(\vb{R})\, \Psi(\vb{R}),
     \label{eq:SchrEq}
 \end{equation} 
 with the eigenvalues $E(\vb{R})$ providing the exact energies.
-\hugh{The energy $E(\vb{R})$ can be considered as a ``one-to-many'' function since each input nuclear geometry
-yields several eigenvalues corresponding to the ground and excited states of the exact spectrum.}
+The energy $E(\vb{R})$ can be considered as a ``one-to-many'' function since each input nuclear geometry
+yields several eigenvalues corresponding to the ground and excited states of the exact spectrum.
 However, exact solutions to Eq.~\eqref{eq:SchrEq} are only possible in the simplest of systems, such as 
 the one-electron hydrogen atom and some specific two-electron systems with well-defined mathematical 
 properties.\cite{Taut_1993,Loos_2009b,Loos_2010e,Loos_2012}
-\hugh{In practice, approximations to the exact Schr\"{o}dinger equation must be introduced, including
+In practice, approximations to the exact Schr\"{o}dinger equation must be introduced, including
 the perturbation theories and Hartree--Fock approximation considered in this review
 In what follows, we will drop the parametric dependence on the nuclear geometry and, 
-unless otherwise stated, atomic units will be used throughout.}
+unless otherwise stated, atomic units will be used throughout.
 
 %===================================%
 \subsection{Exceptional Points in the Hubbard Dimer}
@@ -249,10 +249,8 @@ unless otherwise stated, atomic units will be used throughout.}
 \end{figure*}
 
 To illustrate the concepts discussed throughout this article, we consider the symmetric Hubbard dimer at half filling, \ie\ with two opposite-spin fermions.
-Analytic\trashHB{ally solvable} model systems are essential in theoretical chemistry and physics as their \hugh{mathematical} simplicity \trashHB{of the
-mathematics} compared to realistic systems (e.g., atoms and molecules) allows new concepts and methods to be 
-easily \trashHB{illustrated and} tested while retaining the key physical phenomena.
-\hugh{(HGAB: This sentence felt too long to me. Feel free to re-instate words if you think they are neccessary)}
+Analytically solvable models are essential in theoretical chemistry and physics as their mathematical simplicity compared to realistic systems (e.g., atoms and molecules) allows new concepts and methods to be 
+easily tested while retaining the key physical phenomena.
 
 Using the (localised) site basis, the Hilbert space of the Hubbard dimer comprises the four configurations
 \begin{align}
@@ -274,12 +272,11 @@ where $t$ is the hopping parameter and $U$ is the on-site Coulomb repulsion.
 We refer the interested reader to Refs.~\onlinecite{Carrascal_2015,Carrascal_2018} for more details about this system.
 The parameter $U$ controls the strength of the electron correlation.
 In the weak correlation regime (small $U$), the kinetic energy dominates and the electrons are delocalised over both sites.
-In the large-$U$ (or strong correlation) regime, the electron repulsion term \hugh{becomes dominant} \trashHB{drives the physics}
+In the large-$U$ (or strong correlation) regime, the electron repulsion term becomes dominant 
 and the electrons localise on opposite sites to minimise their Coulomb repulsion. 
 This phenomenon is often referred to as Wigner crystallisation. \cite{Wigner_1934}
 
-To illustrate the formation of an EP, we scale the off-diagonal coupling strength by introducing the complex parameter $\lambda$ through the transformation $t\rightarrow \lambda t$
-\hugh{to give the parameterised Hamiltonian $\hH(\lambda)$.}
+To illustrate the formation of an EP, we scale the off-diagonal coupling strength by introducing the complex parameter $\lambda$ through the transformation $t\rightarrow \lambda t$ to give the parameterised Hamiltonian $\hH(\lambda)$.
 When $\lambda$ is real, the Hamiltonian~\eqref{eq:H_FCI} is Hermitian with the distinct (real-valued) (eigen)energies
 \begin{subequations}
 \begin{align}
@@ -321,7 +318,7 @@ As a result, completely encircling an EP leads to the interconversion of the two
 Additionally, the wave functions pick up a geometric phase in the process, and four complete loops are required to recover their starting forms.\cite{MoiseyevBook}
 
 % LOCATING EPS
-\hugh{To locate EPs in practice, one must simultaneously solve
+To locate EPs in practice, one must simultaneously solve
 \begin{subequations}
 \begin{align}
 	\label{eq:PolChar}
@@ -336,15 +333,15 @@ Equation \eqref{eq:PolChar} is the well-known secular equation providing the (ei
 If the energy is also solution of Eq.~\eqref{eq:DPolChar}, then this energy value is at least two-fold 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$.}
+same symmetry for complex values of $\lambda$.
 
 
 %============================================================%
 \subsection{Rayleigh-Schr\"odinger Perturbation Theory}
 %============================================================%
 
-\hugh{One of the most common routes to approximately solving the Schr\"odinger equation
-is to introduce a perturbative expansion of the exact energy.}
+One of the most common routes to approximately solving the Schr\"odinger equation
+is to introduce a perturbative expansion of the exact energy.
 % SUMMARY OF RS-PT
 Within Rayleigh-Schr\"odinger perturbation theory, the time-independent Schr\"odinger equation 
 is recast as 
@@ -401,10 +398,12 @@ However, this series diverges for $x \ge 1$.
 This divergence occurs because $f(x)$ has four singularities in the complex 
 ($\e^{\i\pi/4}$, $\e^{-\i\pi/4}$, $\e^{\i3\pi/4}$, and $\e^{-\i3\pi/4}$) with a modulus equal to $1$, demonstrating
 that complex singularities are essential to fully understand the series convergence on the real axis.\cite{BenderBook}
-\titou{Include Antoine's example $\sum_{n=1}^\infty \lambda^n/n$ which is divergent at $\lambda = 1$ but convergent at $\lambda = -1$.}
 
-The radius of convergence for the perturbation series Eq.~\eqref{eq:E_expansion} is therefore dictated by the magnitude \titou{$\abs{\lambda_0}$} of the
+The radius of convergence for the perturbation series Eq.~\eqref{eq:E_expansion} is therefore dictated by the magnitude $\abs{\lambda_c}$ of the
 singularity in $E(\lambda)$ that is closest to the origin.
+Note that when $\lambda = \lambda_c$, one cannot \textit{a priori} predict if the series is convergent or not.
+For example, the series $\sum_{k=1}^\infty \lambda^k/k$ diverges at $\lambda = 1$ but converges at $\lambda = -1$.
+
 Like the exact system in Sec.~\ref{sec:example}, the perturbation energy $E(\lambda)$ represents
 a ``one-to-many'' function with the output elements representing an approximation to both the ground and excited states.
 The most common singularities on $E(\lambda)$ therefore correspond to non-analytic EPs in the complex 
@@ -423,7 +422,7 @@ This Slater determinant is defined as an antisymmetric combination of $\Ne$ (rea
 \begin{equation}\label{eq:FockOp}
     \Hat{f}(\vb{x}) \phi_p(\vb{x}) = \qty[ \Hat{h}(\vb{x}) + \Hat{v}_\text{HF}(\vb{x}) ] \phi_p(\vb{x}) = \epsilon_p \phi_p(\vb{x}).
 \end{equation}
-Here the \hugh{(one-electron)} core Hamiltonian is
+Here the (one-electron) core Hamiltonian is
 \begin{equation}
 \label{eq:Hcore}
 	\Hat{h}(\vb{x}) = -\frac{\grad^2}{2} + \sum_{A}^{M} \frac{Z_A}{\abs{\vb{r}-\vb{R}_A}}
@@ -566,8 +565,8 @@ Time-reversal symmetry dictates that this UHF wave function must be degenerate w
 by swapping $\ta^{\text{UHF}}$ and $\tb^{\text{UHF}}$ in Eqs.~\eqref{eq:ta_uhf} and \eqref{eq:tb_uhf}.
 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}
-\hugh{Symmetry breaking can also occur in RHF theory when a charge-density wave is formed from an oscillation 
-between the two closed-shell configurations with both electrons localised on one site or the other.\cite{StuberPaldus,Fukutome_1981}}
+Symmetry breaking can also occur in RHF theory when a charge-density wave is formed from an oscillation 
+between the two closed-shell configurations with both electrons localised on one site or the other.\cite{StuberPaldus,Fukutome_1981}
 
 %============================================================%
 \subsection{Self-Consistency as a Perturbation} %OR {Complex adiabatic connection}
@@ -580,7 +579,7 @@ Alternatively, the non-linear terms arising from the Coulomb and exchange operat
 be considered as a perturbation from the core Hamiltonian \eqref{eq:Hcore} by introducing the
 transformation $U \rightarrow \lambda\, U$, giving the parametrised Fock operator 
 \begin{equation}
-    \Hat{f}(\vb{x} \hugh{; \lambda}) = \Hat{h}(\vb{x}) + \lambda\, \Hat{v}_\text{HF}(\vb{x}).
+    \Hat{f}(\vb{x} ; \lambda) = \Hat{h}(\vb{x}) + \lambda\, \Hat{v}_\text{HF}(\vb{x}).
 \end{equation}
 The orbitals in the reference problem $\lambda=0$ correspond to the symmetry-pure eigenfunctions of the one-electron core
 Hamiltonian, while self-consistent solutions at $\lambda = 1$ represent the orbitals of the true HF solution.
@@ -700,9 +699,9 @@ i) the Epstein-Nesbet (EN) partitioning which consists in taking the diagonal el
 %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}
-\hugh{While an in-depth comparison of these different approaches can offer insight into 
+While an in-depth comparison of these different approaches can offer insight into 
 their relative strengths and weaknesses for various situations, we will restrict our current discussion
-to the convergence properties of the MP expansion.}
+to the convergence properties of the MP expansion.
 
 %=====================================================%
 \subsection{M{\o}ller-Plesset Convergence in Molecular Systems}
@@ -928,7 +927,6 @@ To do so, they analysed the relation between the dominant singularity (\ie, the
 \begin{quote}
 	\textit{``In the limit of large order, the series coefficients become equivalent to the Taylor series coefficients of the singularity closest to the origin. ''}
 \end{quote}
-\titou{T2: should we move this theorem earlier?}
 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). 
 

From abc5f2c96e32bf5099d3e0915cdcabcc3c8072b6 Mon Sep 17 00:00:00 2001
From: Pierre-Francois Loos <pierrefrancois.loos@gmail.com>
Date: Tue, 24 Nov 2020 09:47:29 +0100
Subject: [PATCH 3/3] add Hugh references

---
 Manuscript/EPAWTFT.bbl | 24 +++++++++++++++++++++++-
 Manuscript/EPAWTFT.bib | 42 +++++++++++++++++++++++++++++++++++++++++-
 Manuscript/EPAWTFT.tex |  4 ++--
 3 files changed, 66 insertions(+), 4 deletions(-)

diff --git a/Manuscript/EPAWTFT.bbl b/Manuscript/EPAWTFT.bbl
index 0398b28..56d9b73 100644
--- a/Manuscript/EPAWTFT.bbl
+++ b/Manuscript/EPAWTFT.bbl
@@ -6,7 +6,7 @@
 %Control: page (0) single
 %Control: year (1) truncated
 %Control: production of eprint (0) enabled
-\begin{thebibliography}{113}%
+\begin{thebibliography}{116}%
 \makeatletter
 \providecommand \@ifxundefined [1]{%
  \@ifx{#1\undefined}
@@ -587,6 +587,12 @@
   {Kryachko}}}\ (\bibinfo  {publisher} {Kluwer Academic},\ \bibinfo {address}
   {Dordrecht},\ \bibinfo {year} {2003})\ p.~\bibinfo {pages} {67}\BibitemShut
   {NoStop}%
+\bibitem [{\citenamefont {Slater}(1951)}]{Slater_1951}%
+  \BibitemOpen
+  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {J.~C.}\ \bibnamefont
+  {Slater}},\ }\href {\doibase 10.1103/PhysRev.82.538} {\bibfield  {journal}
+  {\bibinfo  {journal} {Phys. Rev.}\ }\textbf {\bibinfo {volume} {82}},\
+  \bibinfo {pages} {538} (\bibinfo {year} {1951})}\BibitemShut {NoStop}%
 \bibitem [{\citenamefont {Coulson}\ and\ \citenamefont
   {Fischer}(1949)}]{Coulson_1949}%
   \BibitemOpen
@@ -610,6 +616,22 @@
   \bibfield  {author} {\bibinfo {author} {\bibfnamefont {H.}~\bibnamefont
   {Fukutome}},\ }\href {\doibase 10.1002/qua.560200502} {\ \textbf {\bibinfo
   {volume} {20}},\ \bibinfo {pages} {955}}\BibitemShut {NoStop}%
+\bibitem [{\citenamefont {Roothaan}(1951)}]{Roothaan_1951}%
+  \BibitemOpen
+  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {C.~C.~J.}\
+  \bibnamefont {Roothaan}},\ }\href {\doibase 10.1103/RevModPhys.23.69}
+  {\bibfield  {journal} {\bibinfo  {journal} {Rev. Mod. Phys.}\ }\textbf
+  {\bibinfo {volume} {23}},\ \bibinfo {pages} {69} (\bibinfo {year}
+  {1951})}\BibitemShut {NoStop}%
+\bibitem [{\citenamefont {Hall}\ and\ \citenamefont
+  {Lennard-Jones}(1951)}]{Hall_1951}%
+  \BibitemOpen
+  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {G.~G.}\ \bibnamefont
+  {Hall}}\ and\ \bibinfo {author} {\bibfnamefont {J.~E.}\ \bibnamefont
+  {Lennard-Jones}},\ }\href {\doibase 10.1098/rspa.1951.0048} {\bibfield
+  {journal} {\bibinfo  {journal} {Proc. R. Soc. Lond. A}\ }\textbf {\bibinfo
+  {volume} {205}},\ \bibinfo {pages} {541} (\bibinfo {year}
+  {1951})}\BibitemShut {NoStop}%
 \bibitem [{\citenamefont {Hiscock}\ and\ \citenamefont
   {Thom}(2014)}]{Hiscock_2014}%
   \BibitemOpen
diff --git a/Manuscript/EPAWTFT.bib b/Manuscript/EPAWTFT.bib
index 9b5ae24..6b6ae9b 100644
--- a/Manuscript/EPAWTFT.bib
+++ b/Manuscript/EPAWTFT.bib
@@ -1,13 +1,53 @@
 %% This BibTeX bibliography file was created using BibDesk.
 %% http://bibdesk.sourceforge.net/
 
-%% Created for Pierre-Francois Loos at 2020-11-23 11:07:33 +0100 
+%% Created for Pierre-Francois Loos at 2020-11-24 09:46:03 +0100 
 
 
 %% Saved with string encoding Unicode (UTF-8) 
 
 
 
+@article{Hall_1951,
+	abstract = { An analysis of the `linear combination of atomic orbitals' approximation using the accurate molecular orbital equations shows that it does not lead to equations of the form usually assumed in the semi-empirical molecular orbital method. A new semi-empirical method is proposed, therefore, in terms of equivalent orbitals. The equations obtained, which do have the usual form, are applicable to a large class of molecules and do not involve the approximations that were thought necessary. In this method the ionization potentials are calculated by treating certain integrals as semi-empirical parameters. The value of these parameters is discussed in terms of the localization of equivalent orbitals and some approximate rules are suggested. As an illustration the ionization potentials of the paraffin series are considered and good agreement between the observed and calculated values is found. },
+	author = {Hall, G. G. and Lennard-Jones, John Edward},
+	date-added = {2020-11-24 09:45:15 +0100},
+	date-modified = {2020-11-24 09:45:50 +0100},
+	doi = {10.1098/rspa.1951.0048},
+	journal = {Proc. R. Soc. Lond. A},
+	pages = {541-552},
+	title = {The molecular orbital theory of chemical valency VIII. A method of calculating ionization potentials},
+	volume = {205},
+	year = {1951},
+	Bdsk-Url-1 = {https://royalsocietypublishing.org/doi/abs/10.1098/rspa.1951.0048},
+	Bdsk-Url-2 = {https://doi.org/10.1098/rspa.1951.0048}}
+
+@article{Roothaan_1951,
+	author = {Roothaan, C. C. J.},
+	date-added = {2020-11-24 09:43:57 +0100},
+	date-modified = {2020-11-24 09:44:09 +0100},
+	doi = {10.1103/RevModPhys.23.69},
+	journal = {Rev. Mod. Phys.},
+	pages = {69--89},
+	title = {New Developments in Molecular Orbital Theory},
+	volume = {23},
+	year = {1951},
+	Bdsk-Url-1 = {https://link.aps.org/doi/10.1103/RevModPhys.23.69},
+	Bdsk-Url-2 = {https://doi.org/10.1103/RevModPhys.23.69}}
+
+@article{Slater_1951,
+	author = {Slater, J. C.},
+	date-added = {2020-11-24 09:42:40 +0100},
+	date-modified = {2020-11-24 09:42:58 +0100},
+	doi = {10.1103/PhysRev.82.538},
+	journal = {Phys. Rev.},
+	pages = {538--541},
+	title = {Magnetic Effects and the Hartree-Fock Equation},
+	volume = {82},
+	year = {1951},
+	Bdsk-Url-1 = {https://link.aps.org/doi/10.1103/PhysRev.82.538},
+	Bdsk-Url-2 = {https://doi.org/10.1103/PhysRev.82.538}}
+
 @article{Loos_2019d,
 	author = {P. F. Loos and B. Pradines and A. Scemama and J. Toulouse and E. Giner},
 	date-added = {2020-11-23 11:07:32 +0100},
diff --git a/Manuscript/EPAWTFT.tex b/Manuscript/EPAWTFT.tex
index 3532575..9f4477d 100644
--- a/Manuscript/EPAWTFT.tex
+++ b/Manuscript/EPAWTFT.tex
@@ -124,7 +124,7 @@
 \newcommand{\UOX}{Physical and Theoretical Chemical Laboratory, Department of Chemistry, University of Oxford, Oxford, OX1 3QZ, U.K.}
 \begin{document}	
 
-\title{Perturbation theory in the complex plane: 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}
@@ -574,7 +574,7 @@ between the two closed-shell configurations with both electrons localised on one
 
 % INTRODUCE PARAMETRISED FOCK HAMILTONIAN
 The inherent non-linearity in the Fock eigenvalue problem arises from self-consistency 
-in the HF approximation, and is usually solved through an iterative approach.\cite{Roothaan1951,Hall1951}
+in the HF approximation, and is usually solved through an iterative approach.\cite{Roothaan_1951,Hall_1951}
 Alternatively, the non-linear terms arising from the Coulomb and exchange operators can 
 be considered as a perturbation from the core Hamiltonian \eqref{eq:Hcore} by introducing the
 transformation $U \rightarrow \lambda\, U$, giving the parametrised Fock operator