This commit is contained in:
Antoine Marie 2020-07-16 09:58:54 +02:00
parent b811e7e123
commit 8009e2d444
8 changed files with 115 additions and 91 deletions

View File

@ -47,7 +47,7 @@
\citation{Olsen_2000}
\@writefile{toc}{\contentsline {section}{\numberline {1}Introduction}{3}{section.1}\protected@file@percent }
\@writefile{toc}{\contentsline {section}{\numberline {2}Perturbation theory}{3}{section.2}\protected@file@percent }
\citation{Goodson_2011}
\citation{Goodson_2012}
\citation{Moller_1934}
\@writefile{lof}{\contentsline {figure}{\numberline {1}{\ignorespaces \centering A generic EP with the square root branch point topology. A loop around the EP interconvert the states.}}{4}{figure.1}\protected@file@percent }
\newlabel{fig:TopologyEP}{{1}{4}{\centering A generic EP with the square root branch point topology. A loop around the EP interconvert the states}{figure.1}{}}
@ -117,6 +117,8 @@
\newlabel{eq: RHFWF}{{16}{10}{Restricted and unrestricted equation for the spherium model}{equation.4.16}{}}
\newlabel{eq:ERHF}{{17}{10}{Restricted and unrestricted equation for the spherium model}{equation.4.17}{}}
\@writefile{toc}{\contentsline {subsubsection}{\numberline {4.1.2}The minimal basis example}{10}{subsubsection.4.1.2}\protected@file@percent }
\newlabel{eq:uhfstatio}{{20}{11}{The minimal basis example}{equation.4.20}{}}
\newlabel{eq:rhfstatio}{{21}{11}{The minimal basis example}{equation.4.21}{}}
\@writefile{toc}{\contentsline {subsubsection}{\numberline {4.1.3}Symmetry-broken solutions}{11}{subsubsection.4.1.3}\protected@file@percent }
\@writefile{toc}{\contentsline {subsection}{\numberline {4.2}Strongly correlated regime}{11}{subsection.4.2}\protected@file@percent }
\@writefile{toc}{\contentsline {section}{\numberline {5}Radius of convergence and exceptional points}{11}{section.5}\protected@file@percent }
@ -153,33 +155,32 @@
\bibcite{MoiseyevBook}{25}
\bibcite{Olsen_1996}{26}
\bibcite{Olsen_2000}{27}
\bibcite{Goodson_2011}{28}
\bibcite{Moller_1934}{29}
\bibcite{Gill_1986}{30}
\bibcite{Gill_1988}{31}
\bibcite{Handy_1985}{32}
\bibcite{Lepetit_1988}{33}
\bibcite{SzaboBook}{34}
\bibcite{Fukutome_1981}{35}
\bibcite{Cremer_1996}{36}
\bibcite{Christiansen_1996}{37}
\bibcite{Sergeev_2005}{38}
\bibcite{Sergeev_2006}{39}
\bibcite{Stillinger_2000}{40}
\bibcite{Baker_1971}{41}
\bibcite{Goodson_2004}{42}
\bibcite{Heiss_1988}{43}
\bibcite{Heiss_2002}{44}
\bibcite{Cejnar_2005}{45}
\bibcite{Cejnar_2007}{46}
\bibcite{Cejnar_2009}{47}
\bibcite{Borisov_2015}{48}
\bibcite{Sindelka_2017}{49}
\bibcite{Sachdev_2011}{50}
\bibcite{Cejnar_2015}{51}
\bibcite{Cejnar_2016}{52}
\bibcite{Caprio_2008}{53}
\bibcite{Macek_2019}{54}
\bibcite{Stransky_2018}{55}
\bibcite{Coulson_1949}{56}
\bibcite{Moller_1934}{28}
\bibcite{Gill_1986}{29}
\bibcite{Gill_1988}{30}
\bibcite{Handy_1985}{31}
\bibcite{Lepetit_1988}{32}
\bibcite{SzaboBook}{33}
\bibcite{Fukutome_1981}{34}
\bibcite{Cremer_1996}{35}
\bibcite{Christiansen_1996}{36}
\bibcite{Sergeev_2005}{37}
\bibcite{Sergeev_2006}{38}
\bibcite{Stillinger_2000}{39}
\bibcite{Baker_1971}{40}
\bibcite{Goodson_2004}{41}
\bibcite{Heiss_1988}{42}
\bibcite{Heiss_2002}{43}
\bibcite{Cejnar_2005}{44}
\bibcite{Cejnar_2007}{45}
\bibcite{Cejnar_2009}{46}
\bibcite{Borisov_2015}{47}
\bibcite{Sindelka_2017}{48}
\bibcite{Sachdev_2011}{49}
\bibcite{Cejnar_2015}{50}
\bibcite{Cejnar_2016}{51}
\bibcite{Caprio_2008}{52}
\bibcite{Macek_2019}{53}
\bibcite{Stransky_2018}{54}
\bibcite{Coulson_1949}{55}
\@writefile{toc}{\contentsline {section}{\numberline {A}ERHF and EUHF}{15}{appendix.A}\protected@file@percent }

View File

@ -154,11 +154,6 @@ Jeppe Olsen, Poul Jørgensen, Trygve Helgaker, and Ove Christiansen.
two-state model.
\newblock {\em J. Chem. Phys.}, 112(22):9736--9748.
\bibitem{Goodson_2011}
David~Z. Goodson.
\newblock Resummation methods.
\newblock 2(5):743--761.
\bibitem{Moller_1934}
Chr. Møller and M.~S. Plesset.
\newblock Note on an approximation treatment for many-electron systems.

View File

@ -143,15 +143,15 @@
}
@article{Goodson_2011,
Title = {Resummation methods},
Volume = {2},
@article{goodson_resummation_2012,
title = {Resummation methods},
volume = {2},
doi = {10.1002/wcms.92},
pages = {743--761},
number = {5},
journalTitle = {{WIREs} Computational Molecular Science},
journaltitle = {{WIREs} Computational Molecular Science},
author = {Goodson, David Z.},
Date = {2011},
date = {2012},
}
@article{Katz_1962,

View File

@ -3,10 +3,9 @@ Capacity: max_strings=100000, hash_size=100000, hash_prime=85009
The top-level auxiliary file: Rapport.aux
The style file: unsrt.bst
Database file #1: Rapport.bib
Warning--I didn't find a database entry for "Goodson_2012"
Warning--empty year in Olsen_1996
Warning--empty year in Olsen_2000
Warning--empty journal in Goodson_2011
Warning--empty year in Goodson_2011
Warning--empty year in Moller_1934
Warning--empty year in Gill_1986
Warning--empty year in Gill_1988
@ -50,45 +49,45 @@ Warning--empty journal in Stransky_2018
Warning--empty year in Stransky_2018
Warning--empty journal in Coulson_1949
Warning--empty year in Coulson_1949
You've used 56 entries,
You've used 55 entries,
1791 wiz_defined-function locations,
809 strings with 12646 characters,
and the built_in function-call counts, 13722 in all, are:
= -- 1130
> -- 651
806 strings with 12602 characters,
and the built_in function-call counts, 13500 in all, are:
= -- 1109
> -- 647
< -- 2
+ -- 240
- -- 184
* -- 1239
:= -- 2181
add.period$ -- 169
call.type$ -- 56
change.case$ -- 54
+ -- 238
- -- 183
* -- 1220
:= -- 2149
add.period$ -- 166
call.type$ -- 55
change.case$ -- 53
chr.to.int$ -- 0
cite$ -- 103
duplicate$ -- 520
empty$ -- 1321
format.name$ -- 184
if$ -- 2902
cite$ -- 100
duplicate$ -- 511
empty$ -- 1297
format.name$ -- 183
if$ -- 2856
int.to.chr$ -- 0
int.to.str$ -- 56
missing$ -- 59
newline$ -- 283
num.names$ -- 56
pop$ -- 140
int.to.str$ -- 55
missing$ -- 58
newline$ -- 278
num.names$ -- 55
pop$ -- 136
preamble$ -- 1
purify$ -- 0
quote$ -- 0
skip$ -- 187
skip$ -- 184
stack$ -- 0
substring$ -- 1146
substring$ -- 1121
swap$ -- 49
text.length$ -- 2
text.prefix$ -- 0
top$ -- 0
type$ -- 0
warning$ -- 47
while$ -- 130
width$ -- 58
write$ -- 572
(There were 47 warnings)
warning$ -- 45
while$ -- 127
width$ -- 57
write$ -- 563
(There were 46 warnings)

View File

@ -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) 15 JUL 2020 17:49
This is pdfTeX, Version 3.14159265-2.6-1.40.20 (TeX Live 2019/Arch Linux) (preloaded format=pdflatex 2020.5.9) 16 JUL 2020 09:51
entering extended mode
restricted \write18 enabled.
%&-line parsing enabled.
@ -838,8 +838,12 @@ Package pdftex.def Info: TopologyEP.pdf used on input line 158.
LaTeX Warning: `!h' float specifier changed to `!ht'.
[3] [4 <./TopologyEP.pdf>]
<gill1986.png, id=248, 328.22626pt x 329.23pt>
[3]
LaTeX Warning: Citation `Goodson_2012' on page 4 undefined on input line 180.
[4 <./TopologyEP.pdf>]
<gill1986.png, id=246, 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 211.
@ -852,12 +856,12 @@ LaTeX Warning: `!h' float specifier changed to `!ht'.
LaTeX Warning: `!h' float specifier changed to `!ht'.
[5]
<Nedivergence.png, id=277, 329.23pt x 234.8775pt>
<Nedivergence.png, id=275, 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 243.
(pdftex.def) Requested size: 219.35216pt x 156.49014pt.
<HFdivergence.png, id=278, 322.20375pt x 231.86626pt>
<HFdivergence.png, id=276, 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 245.
@ -877,25 +881,29 @@ 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 316.
[9] [10] [11] (./Rapport.bbl
[12] [13]) [14]
Package atveryend Info: Empty hook `BeforeClearDocument' on input line 432.
Package atveryend Info: Empty hook `BeforeClearDocument' on input line 453.
[15]
Package atveryend Info: Empty hook `AfterLastShipout' on input line 432.
Package atveryend Info: Empty hook `AfterLastShipout' on input line 453.
(./Rapport.aux)
Package atveryend Info: Executing hook `AtVeryEndDocument' on input line 432.
Package atveryend Info: Executing hook `AtEndAfterFileList' on input line 432.
Package atveryend Info: Executing hook `AtVeryEndDocument' on input line 453.
Package atveryend Info: Executing hook `AtEndAfterFileList' on input line 453.
Package rerunfilecheck Info: File `Rapport.out' has not changed.
(rerunfilecheck) Checksum: 747F8F3DB1A902E47B11902AB4777D3F;1342.
Package atveryend Info: Empty hook `AtVeryVeryEnd' on input line 432.
LaTeX Warning: There were undefined references.
Package atveryend Info: Empty hook `AtVeryVeryEnd' on input line 453.
)
(\end occurred inside a group at level 1)
### simple group (level 1) entered at line 206 ({)
### bottom level
Here is how much of TeX's memory you used:
21816 strings out of 492167
399412 string characters out of 6131559
494570 words of memory out of 5000000
25866 multiletter control sequences out of 15000+600000
21818 strings out of 492167
399435 string characters out of 6131559
494643 words of memory out of 5000000
25867 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,1615b,998s stack positions out of 5000i,500n,10000p,200000b,80000s
@ -920,10 +928,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 (15 pages, 1197396 bytes).
Output written on Rapport.pdf (15 pages, 1199296 bytes).
PDF statistics:
496 PDF objects out of 1000 (max. 8388607)
432 compressed objects within 5 object streams
117 named destinations out of 1000 (max. 500000)
118 named destinations out of 1000 (max. 500000)
203 words of extra memory for PDF output out of 10000 (max. 10000000)

Binary file not shown.

Binary file not shown.

View File

@ -177,7 +177,7 @@ The energy can then be written as a power series of $\lambda$
\begin{equation}
E(\lambda) = \sum_{k=0}^\infty \lambda^k E^{(k)}
\end{equation}
where $\lambda$ is a coupling parameter set equal to 1 at the end of the calculation. However it is not guaranteed that the series $E(\lambda)$ has a radius of convergence $\abs{\lambda_0} < 1$. It means that the series is divergent for the physical system at $\lambda=1$. One can prove that $\abs{\lambda_0}$ can be obtained by extending $\lambda$ in the complex plane and looking for the singularities of $E(\lambda)$. This is due to the following theorem \cite{Goodson_2011}: The Taylor series about a point $z_0$ of a function over the complex $z$ plane will converge at a value $z_1$ if the function is non-singular at all values of $z$ in the circular region centered at $z_0$ with radius $\abs{z_1 z_0}$. If the function has a singular point $z_s$ such that $\abs{z_s z_0} < \abs{z_1 z_0}$, then the series will diverge when evaluated at $z_1$. This theorem means that the radius of convergence of the perturbation series is equal to distance to the origin of the closest singularity of $E(\lambda)$.
where $\lambda$ is a coupling parameter set equal to 1 at the end of the calculation. However it is not guaranteed that the series $E(\lambda)$ has a radius of convergence $\abs{\lambda_0} < 1$. It means that the series is divergent for the physical system at $\lambda=1$. One can prove that $\abs{\lambda_0}$ can be obtained by extending $\lambda$ in the complex plane and looking for the singularities of $E(\lambda)$. This is due to the following theorem \cite{Goodson_2012}: The Taylor series about a point $z_0$ of a function over the complex $z$ plane will converge at a value $z_1$ if the function is non-singular at all values of $z$ in the circular region centered at $z_0$ with radius $\abs{z_1 z_0}$. If the function has a singular point $z_s$ such that $\abs{z_s z_0} < \abs{z_1 z_0}$, then the series will diverge when evaluated at $z_1$. This theorem means that the radius of convergence of the perturbation series is equal to distance to the origin of the closest singularity of $E(\lambda)$.
The discovery of a partitioning of the Hamiltonian that allowed chemists to recover a part of the correlation energy (i.e. the difference between the exact energy and the Hartree-Fock energy) using perturbation theory has been a major step in the development of post-Hartree-Fock methods. This case of the Rayleigh-Schrödinger perturbation theory is called the M{\o}ller-Plesset perturbation theory \cite{Moller_1934}. In the MPPT the unperturbed Hamiltonian is the sum of the $n$ mono-electronic Fock operators which are the sum of the one-electron core Hamiltonian $h(i)$, the Coulomb $J_j(i)$ and Exchange $K_j(i)$ operators.
@ -295,9 +295,9 @@ Finally, it was shown that $\beta$ singularities are very sensitive to the basis
In the previous section, we saw that a reasoning on 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 singularity $\beta$. It is now well-known that this phenomenon is a specific case of a more general phenomenon. Indeed, theoretical physicists proved that EPs are connected to quantum phase transitions \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 a parameter, in some cases the variation of a parameter can lead to abrupt changes at a critical point. Those quantum phase transitions exist both for ground and excited states \cite{Cejnar_2009, Sachdev_2011, Cejnar_2015, Cejnar_2016, Caprio_2008, Macek_2019}. A ground-state quantum phase transition is characterized by the successive derivative 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 n-th order if the n-th derivative is discontinuous. A quantum phase transition can also be identify by the discontinuity of an appropriate order parameter (or one of its derivative).
The presence of an EP close to the real axis is characteristic of a sharp avoided crossings. Yet at such an avoided crossings eigenstates change abruptly. Although it is now well understood that EPs are closely related to quantum phase transitions, the link between the type of QPT (ground state or excited state, first or superior order) and EPs still need to be clarify. One of the major challenge in order to do this reside in our 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. Cejnar et al. developped 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 co-workers proved that the distribution of EPs is not the same around a QPT of first or second order \cite{Stransky_2018}. Moreover, that when the dimension of the system increases they tends towards the real axis in a different manner, meaning respectively exponentially and algebraically.
The presence of an EP close to the real axis is characteristic of a sharp avoided crossing. Yet at such an avoided crossing eigen states change abruptly. Although it is now well understood that EPs are closely related to quantum phase transitions, the link between the type of QPT (ground state or excited state, first or superior order) and EPs still need to be clarify. One of the major challenge in order to do this reside in our 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. Cejnar 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 co-workers proved that the distribution of EPs is not the same around a QPT of first or second order \cite{Stransky_2018}. Moreover, that when the dimension of the system increases they tends towards the real axis in a different manner, meaning respectively exponentially and algebraically.
It seems like our understanding of the physics of spatial and/or spin symmetry breaking in the Hartree-Fock theory can be enlightened by quantum phase transition theory. Indeed, the second derivative of the energy is discontinuous at the Coulson-Fischer point which mean that the system undergo a second order quantum phase transition. The $\beta$ singularities introduced by Sergeev to describe the EPs modeling the formation of a bound cluster of electrons are actually a more general class of singularities.
It seems like our understanding of the physics of spatial and/or spin symmetry breaking in the Hartree-Fock theory can be enlightened by quantum phase transition theory. Indeed, the second derivative of the energy is discontinuous at the Coulson-Fischer point which means that the system undergo a second order quantum phase transition. The $\beta$ singularities introduced by Sergeev to describe the EPs modeling the formation of a bound cluster of electrons are actually a more general class of singularities.
Singularity $\alpha$ and quantum phase transition ?
@ -383,8 +383,29 @@ Because the transformation between two basis sets needs to be unitary, we get th
\phi(\theta_1)= -\sin(\chi)\frac{Y_{00}(\Omega_1)}{R} + \cos(\chi)\frac{Y_{10}(\Omega_1)}{R}
\end{equation}
The minimization of the UHF and RHF ground state energies with respect to $\chi$ respectively yiels the following equations:
\begin{equation}
\frac{2(75+6R-56R\cos(2\chi))\sin(2\chi))}{75R^2}=0
\label{eq:uhfstatio}
\end{equation}
\begin{equation}
\frac{2(75+6R+44R\cos(2\chi))\sin(2\chi))}{75R^2}=0
\label{eq:rhfstatio}
\end{equation}
It is evident that those equations are both verified for $\chi=\frac{n\pi}{2}$ with $n$ an integer. Those value of $\chi$ lead to the
\subsubsection{Symmetry-broken solutions}
Apparition of the sb-UHF solution
Negative R possible, sb-UHF now a maximum
Less know sb-RHF solution, maximum and minimum inversed
Charge density wave vs Spin density wave
NRJ graphics
\subsection{Strongly correlated regime}
\section{Radius of convergence and exceptional points}