From 9e70ba237c641374b2ff704c09f97f5858f244d6 Mon Sep 17 00:00:00 2001 From: Pierre-Francois Loos Date: Thu, 30 Jul 2020 11:13:10 +0200 Subject: [PATCH 1/3] saving work in sec 4 --- RapportStage/Rapport.tex | 38 ++++++++++++++++++++------------------ 1 file changed, 20 insertions(+), 18 deletions(-) diff --git a/RapportStage/Rapport.tex b/RapportStage/Rapport.tex index b741062..4c8f388 100644 --- a/RapportStage/Rapport.tex +++ b/RapportStage/Rapport.tex @@ -641,28 +641,28 @@ and \titou{T2: please define Wigner 3j symbols.} \titou{STOPPED HERE.} -We obtained Eq.~\eqref{eq:EUHF} for the general form of the wave function \eqref{eq:UHF_WF}, but to be associated with a physical wave function the energy needs to be stationary with respect to the coefficients. The general method is to use the Hartree-Fock self-consistent field method \cite{SzaboBook} to get the coefficients of the wave functions corresponding to physical solutions. We will work in a minimal basis, composed of a $Y_{00}$ and a $Y_{10}$ spherical harmonic, or equivalently a s and a p\textsubscript{z} orbital, to illustrate the difference between the RHF and UHF solutions. In this basis there is a shortcut to find the stationary solutions. One can define the one-electron wave functions $\phi(\theta)$ using a mixing angle between the two basis functions. Hence we just need to minimize the energy with respect to the two mixing angles $\chi_\alpha$ and $\chi_\beta$. +Equation \eqref{eq:EUHF} is obtained from the general form of the wave function \eqref{eq:UHF_WF}, but to be associated with a physical wave function the energy has to be stationary with respect to the coefficients. The general method is to use a self-consistent field procedure \cite{SzaboBook} to get the coefficients of the wave functions corresponding to physical solutions. We will work in a minimal basis, composed of $Y_{0}$ and $Y_{1}$, or equivalently, a s and p\textsubscript{z} orbitals, to illustrate the difference between the RHF and UHF solutions. In this basis there is a shortcut to find the stationary solutions. One can define the one-electron wave functions $\phi(\theta)$ using a mixing angle between the two basis functions. Hence we just need to minimize the energy with respect to the two mixing angles $\chi_\alpha$ and $\chi_\beta$. \begin{equation} -\phi_\alpha(\theta_1)= \cos(\chi_\alpha)\frac{Y_{00}(\theta_1)}{R} + \sin(\chi_\alpha)\frac{Y_{10}(\theta_1)}{R} + \phi_\sigma(\theta)= \cos(\chi_\sigma)Y_{0}(\theta) + \sin(\chi_\sigma)Y_{1}(\theta) \end{equation} -The minimization gives the three following solutions valid for all value of R: +The minimization gives the three following solutions valid for all value of $R$: \begin{itemize} -\item The two electrons are in the s orbital which is a RHF solution. This solution is associated with the energy $R^{-2}$. -\item The two electrons are in the p\textsubscript{z} orbital which is a RHF solution. This solution is associated with the energy $R^{-2} + R^{-1}$ -\item One electron is in the s orbital and the other is in the p\textsubscript{z} orbital which is a UHF solution. This solution is associated with the energy $2R^{-2}+\frac{29}{25}R^{-1}$ +\item The two electrons are in the s orbital which is a RHF solution. This solution is associated with the energy $1/R^{2}$. +\item The two electrons are in the p\textsubscript{z} orbital which is a RHF solution. This solution is associated with the energy $1/R^{2} + 1/R$ +\item One electron is in the s orbital and the other is in the p\textsubscript{z} orbital which is a UHF solution. This solution is associated with the energy $2/R^{2}+ 29/(25R)$ \end{itemize} Those solutions are respectively a minimum, a maximum and a saddle point of the HF equations.\\ -In addition, there is also the well-known symmetry-broken UHF (sb-UHF) solution. For $R>3/2$ an other stationary UHF solution appears, this solution is a minimum of the HF equations. This solution corresponds to the configuration with the electron $\alpha$ on one side of the sphere and the electron $\beta$ on the opposite side and the configuration the other way round. The electrons can be on opposite sides of the sphere because the choice of p\textsubscript{z} as a basis function induced a privileged axis on the sphere for the electrons. This solution have the energy \eqref{eq:EsbUHF} for $R>3/2$. +In addition, there is also the well-known symmetry-broken UHF (sb-UHF) solution. For $R>3/2$ an other stationary UHF solution appears, this solution is a minimum of the HF equations. This solution corresponds to the configuration with the electron $\alpha$ on one side of the sphere and the electron $\beta$ on the opposite side and the configuration the other way round. The electrons can be on opposite sides of the sphere because the choice of p\textsubscript{z} as a basis function induced a privileged axis on the sphere for the electrons. This solution has the energy \eqref{eq:EsbUHF} for $R>3/2$. \begin{equation}\label{eq:EsbUHF} E_{\text{sb-UHF}}=-\frac{75}{112R^3}+\frac{25}{28R^2}+\frac{59}{84R} \end{equation} -The exact solution for the ground state is a singlet so this wave function does not have the true symmetry. Indeed, the spherical harmonics are eigenvectors of $S^2$ but the symmetry-broken solution is a linear combination of the two eigenvectors and is not an eigenvector of $S^2$. However this solution gives more accurate results for the energy at large R as shown in Table \ref{tab:ERHFvsEUHF}. In fact at the Coulson-Fischer point, it becomes more efficient to minimize the Coulomb repulsion than the kinetic energy in order to minimize the total energy. Thus the wave function break the spin symmetry because it allows a more efficient minimization of the Coulomb repulsion. This type of symmetry breaking is called a spin density wave because the system oscillates between the two symmetry-broken configurations \cite{GiulianiBook}. +The exact solution for the ground state is a singlet so this wave function does not have the true symmetry. Indeed, the spherical harmonics are eigenvectors of $S^2$ but the symmetry-broken solution is a linear combination of the two eigenvectors and is not an eigenvector of $S^2$. However this solution gives more accurate results for the energy at large $R$ as shown in Table \ref{tab:ERHFvsEUHF}. In fact at the Coulson-Fischer point, it becomes more efficient to minimize the Coulomb repulsion than the kinetic energy in order to minimize the total energy. Thus the wave function break the spin symmetry because it allows a more efficient minimization of the Coulomb repulsion. This type of symmetry breaking is called a spin density wave because the system oscillates between the two symmetry-broken configurations \cite{GiulianiBook}. \begin{table}[h!] \centering @@ -681,15 +681,6 @@ Exact & 9.783874 & 0.852781 & 0.391959 & 0.247898 & 0.139471 & 0.064525 & 0.005 \label{tab:ERHFvsEUHF} \end{table} -There is also another symmetry-broken solution for $R>75/38$ but this one corresponds to a maximum of the HF equations. This solution is associated with another type of symmetry breaking somewhat less known. It corresponds to a configuration where both electrons are on the same side of the sphere, in the same spatial orbital. This solution is called symmetry-broken RHF (sb-RHF). At the critical value of R, the repulsion of the two electrons on the same side of the sphere maximizes more the energy than the kinetic energy of the p\textsubscript{z} orbitals. This symmetry breaking is associated with a charge density wave, the system oscillates between the situations with the electrons on each side \cite{GiulianiBook}. - -\begin{equation} -E_{\text{sb-RHF}}=\frac{75}{88R^3}+\frac{25}{22R^2}+\frac{91}{66R} -\end{equation} - -We can also consider negative values of R. This corresponds to the situation where one of the electrons is replaced by a positron. There are also a sb-RHF ($R<-3/2$) and a sb-UHF ($R<-75/38$) solution for negative values of R (see Fig.\ref{fig:SpheriumNrj} but in this case the sb-RHF solution is a minimum and the sb-UHF is a maximum. Indeed, the sb-RHF states minimize the energy by placing the electron and the positron on the same side of the sphere. And the sb-UHF states maximize the energy because the two attracting particles are on opposite sides of the sphere. - -In addition, we can also consider the symmetry-broken solutions beyond their respective Coulson-Fischer points by analytically continuing their respective energies leading to the so-called holomorphic solutions \cite{Hiscock_2014, Burton_2019, Burton_2019a}. All those energies are plotted in Fig.~\ref{fig:SpheriumNrj}. The dotted curves corresponds to the holomorphic domain of the energies. \begin{wrapfigure}{r}{0.5\textwidth} \centering @@ -698,11 +689,22 @@ In addition, we can also consider the symmetry-broken solutions beyond their res \label{fig:SpheriumNrj} \end{wrapfigure} +There is also another symmetry-broken solution for $R>75/38$ but this one corresponds to a maximum of the HF equations. This solution is associated with another type of symmetry breaking somewhat less known. It corresponds to a configuration where both electrons are on the same side of the sphere, in the same spatial orbital. This solution is called symmetry-broken RHF (sb-RHF). At the critical value of $R$, the repulsion of the two electrons on the same side of the sphere maximizes more the energy than the kinetic energy of the p\textsubscript{z} orbitals. This symmetry breaking is associated with a charge density wave, the system oscillates between the situations with the electrons on each side \cite{GiulianiBook}. + +\begin{equation} +E_{\text{sb-RHF}}=\frac{75}{88R^3}+\frac{25}{22R^2}+\frac{91}{66R} +\end{equation} + +We can also consider negative values of $R$. This corresponds to the situation where one of the electrons is replaced by a positron. There are also a sb-RHF ($R<-3/2$) and a sb-UHF ($R<-75/38$) solution for negative values of $R$ (see Fig.~\ref{fig:SpheriumNrj} but in this case the sb-RHF solution is a minimum and the sb-UHF is a maximum. Indeed, the sb-RHF states minimize the energy by placing the electron and the positron on the same side of the sphere. And the sb-UHF states maximize the energy because the two attracting particles are on opposite sides of the sphere. + +In addition, we can also consider the symmetry-broken solutions beyond their respective Coulson-Fischer points by analytically continuing their respective energies leading to the so-called holomorphic solutions \cite{Hiscock_2014, Burton_2019, Burton_2019a}. All those energies are plotted in Fig.~\ref{fig:SpheriumNrj}. The dotted curves corresponds to the holomorphic domain of the energies. + + \section{Radius of convergence and exceptional points} \subsection{Evolution of the radius of convergence} -In this part, we will try to investigate how some parameters of $\hH(\lambda)$ influence the radius of convergence of the perturbation series. The radius of convergence is equal to the distance of the closest singularity to the origin of $E(\lambda)$. Hence we need to determine the locations of the exceptional points to obtain information on the convergence properties. To find them we solve simultaneously Eqs.~\eqref{eq:PolChar} and \eqref{eq:DPolChar}. Equation \eqref{eq:PolChar} is the well-known secular equation giving the energies of the system where $\hI$ is the identity operator. If an energy is also solution of Eq.~\eqref{eq:DPolChar} then this energy is degenerate. In this case the energies obtained are dependent of $\lambda$ so solving those equations with respect to $E$ and $\lambda$ gives the value of $\lambda$ where two energies are degenerate. These degeneracies can be conical intersections between two states with different symmetry for real value of $\lambda$ or exceptional points between two states with the same symmetry for complex value of $\lambda$. +In this part, we will try to investigate how some parameters of $\hH(\lambda)$ influence the radius of convergence of the perturbation series. The radius of convergence is equal to the distance of the closest singularity to the origin of $E(\lambda)$. Hence, we have to determine the locations of the exceptional points to obtain information on the convergence properties. To find them we solve simultaneously Eqs.~\eqref{eq:PolChar} and \eqref{eq:DPolChar}. Equation \eqref{eq:PolChar} is the well-known secular equation giving the energies of the system where $\hI$ is the identity operator. If an energy is also solution of Eq.~\eqref{eq:DPolChar} then this energy is degenerate. In this case the energies obtained are dependent of $\lambda$ so solving those equations with respect to $E$ and $\lambda$ gives the value of $\lambda$ where two energies are degenerate. These degeneracies can be conical intersections between two states with different symmetry for real value of $\lambda$ or exceptional points between two states with the same symmetry for complex value of $\lambda$. \begin{subequations} \begin{align} From 418e648d2c6c9b74bc4e9a7f32fc307ce89d4ea0 Mon Sep 17 00:00:00 2001 From: Pierre-Francois Loos Date: Thu, 30 Jul 2020 16:05:16 +0200 Subject: [PATCH 2/3] HF part for spherium --- RapportStage/Rapport.tex | 36 ++++++++++++++++++++---------------- 1 file changed, 20 insertions(+), 16 deletions(-) diff --git a/RapportStage/Rapport.tex b/RapportStage/Rapport.tex index 41ebe34..bde598d 100644 --- a/RapportStage/Rapport.tex +++ b/RapportStage/Rapport.tex @@ -616,23 +616,22 @@ These one-electron orbitals are expanded in the basis of zonal spherical harmoni \begin{equation} \phi_\sigma(\theta)=\sum_{\ell=0}^{\infty}C_{\sigma,\ell}Y_{\ell}(\theta) \end{equation} -It is possible to obtain the formula for the ground state UHF energy in this basis set \cite{Loos_2009}: - +It is possible to obtain the formula for the HF energy in this basis set \cite{Loos_2009}: \begin{equation} -E_{\text{UHF}} = T_{\alpha} + T_{\beta} + V -\label{eq:EUHF} + E_{\text{HF}} = T_{\text{HF}} + V_{\text{HF}} +\label{eq:EHF} \end{equation} -with -\begin{gather} - T_{\sigma} = \frac{1}{R^2} \sum_{\ell=0}^{\infty} C_{\sigma,\ell}^2 \, \ell(\ell+1) - \\ - V = \frac{1}{R} \sum_{\ell_1,\ell_2,\ell_3,\ell_4=0}^{\infty} \sum_{L=0}^{\infty} - (-1)^{\ell_3+\ell_4} v^\alpha_{\ell_1,\ell_2,L} v^\beta_{\ell_3,\ell_4,L} -\end{gather} +where the kinetic and potential energies are +\begin{align} + T_{\text{HF}} & = \sum_{\sigma=\alpha,\beta} \frac{1}{R^2} \sum_{\ell=0}^{\infty} C_{\sigma,\ell}^2 \, \ell(\ell+1) + & + V_{\text{HF}} & = \frac{1}{R} \sum_{L=0}^{\infty} + v^\alpha_{L} v^\beta_{L} +\end{align} and \begin{equation} - v^\sigma_{\ell_1,\ell_2,L} - = \sqrt{(2\ell_1+1)(2\ell_2+1)} C_{\sigma,\ell_1}C_{\sigma,\ell_2} + v^\sigma_{L} + = \sum_{\ell_1,\ell_2} \sqrt{(2\ell_1+1)(2\ell_2+1)} C_{\sigma,\ell_1}C_{\sigma,\ell_2} \begin{pmatrix} \ell_1 & \ell_2 & L \\ @@ -640,13 +639,18 @@ and \end{pmatrix}^2 \end{equation} is expressed in terms of the Wigner 3j-symbols \cite{AngularBook}. -\titou{STOPPED HERE.} - -Equation \eqref{eq:EUHF} is obtained from the general form of the wave function \eqref{eq:UHF_WF}, but to be associated with a physical wave function the energy has to be stationary with respect to the coefficients. The general method is to use a self-consistent field procedure \cite{SzaboBook} to get the coefficients of the wave functions corresponding to physical solutions. We will work in a minimal basis, composed of $Y_{0}$ and $Y_{1}$, or equivalently, a s and p\textsubscript{z} orbitals, to illustrate the difference between the RHF and UHF solutions. In this basis there is a shortcut to find the stationary solutions. One can define the one-electron wave functions $\phi(\theta)$ using a mixing angle between the two basis functions. Hence we just need to minimize the energy with respect to the two mixing angles $\chi_\alpha$ and $\chi_\beta$. +The general method is to use a self-consistent field procedure as described in Ref.~\cite{SzaboBook} to get the coefficients of the wave functions corresponding to stationary solutions with respect to the coefficients $C_{\sigma,\ell}$, i.e., +\begin{equation} + \pdv{E_{\text{HF}}}{C_{\sigma,\ell}} = 0. +\end{equation}. +Here, we work in a minimal basis, composed of $Y_{0}$ and $Y_{1}$, or equivalently, a s and p\textsubscript{z} orbital, to illustrate the difference between the RHF and UHF solutions. In this basis there is a shortcut to find the stationary solutions. One can define the one-electron wave functions \begin{equation} \phi_\sigma(\theta)= \cos(\chi_\sigma)Y_{0}(\theta) + \sin(\chi_\sigma)Y_{1}(\theta) \end{equation} +using a mixing angle between the two basis functions for each spin manifold. +Hence we just minimize the energy with respect to the two mixing angles $\chi_\alpha$ and $\chi_\beta$. +\titou{STOPPED HERE.} The minimization gives the three following solutions valid for all value of $R$: \begin{itemize} From 22afac608bef92c62d2ed407eae61b3275118f3a Mon Sep 17 00:00:00 2001 From: Pierre-Francois Loos Date: Thu, 30 Jul 2020 16:08:43 +0200 Subject: [PATCH 3/3] HF part for spherium --- RapportStage/Rapport.tex | 10 ++++------ 1 file changed, 4 insertions(+), 6 deletions(-) diff --git a/RapportStage/Rapport.tex b/RapportStage/Rapport.tex index bde598d..325a2ec 100644 --- a/RapportStage/Rapport.tex +++ b/RapportStage/Rapport.tex @@ -650,16 +650,14 @@ Here, we work in a minimal basis, composed of $Y_{0}$ and $Y_{1}$, or equivalent \end{equation} using a mixing angle between the two basis functions for each spin manifold. Hence we just minimize the energy with respect to the two mixing angles $\chi_\alpha$ and $\chi_\beta$. -\titou{STOPPED HERE.} -The minimization gives the three following solutions valid for all value of $R$: +This process provides the three following solutions valid for all value of $R$, which are respectively a minimum, a maximum and a saddle point of the HF equations: \begin{itemize} \item The two electrons are in the s orbital which is a RHF solution. This solution is associated with the energy $1/R^{2}$. -\item The two electrons are in the p\textsubscript{z} orbital which is a RHF solution. This solution is associated with the energy $1/R^{2} + 1/R$ -\item One electron is in the s orbital and the other is in the p\textsubscript{z} orbital which is a UHF solution. This solution is associated with the energy $2/R^{2}+ 29/(25R)$ +\item The two electrons are in the p\textsubscript{z} orbital which is a RHF solution. This solution is associated with the energy $1/R^{2} + 1/R$. +\item One electron is in the s orbital and the other is in the p\textsubscript{z} orbital which is a UHF solution. This solution is associated with the energy $2/R^{2}+ 29/(25R)$. \end{itemize} - -Those solutions are respectively a minimum, a maximum and a saddle point of the HF equations.\\ +\titou{STOPPED HERE.} In addition, there is also the well-known symmetry-broken UHF (sb-UHF) solution. For $R>3/2$ an other stationary UHF solution appears, this solution is a minimum of the HF equations. This solution corresponds to the configuration with the electron $\alpha$ on one side of the sphere and the electron $\beta$ on the opposite side and the configuration the other way round. The electrons can be on opposite sides of the sphere because the choice of p\textsubscript{z} as a basis function induced a privileged axis on the sphere for the electrons. This solution has the energy \eqref{eq:EsbUHF} for $R>3/2$.