minor corrections up to the end

This commit is contained in:
Pierre-Francois Loos 2020-07-20 16:46:17 +02:00
parent ef4e33e878
commit cc3008b9f1

View File

@ -186,14 +186,14 @@ In the next section, we illustrate these concepts and the various tricks that ca
\label{sec:exact}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Let us consider a two-level quantum system made of two orbitals in its singlet ground state (\ie, the lowest orbital is doubly occupied). \cite{Romaniello_2009b}
Let us consider a two-level quantum system made of two electrons in its singlet ground state (\ie, the lowest orbital is doubly occupied). \cite{Romaniello_2009b}
We will label these two orbitals, $\MO{v}$ and $\MO{c}$, as valence ($v$) and conduction ($c$) orbitals with respective one-electron Hartree-Fock (HF) energies $\e{v}$ and $\e{c}$.
In a more quantum chemical language, these correspond to the HOMO and LUMO orbitals (respectively).
The ground state $\ket{0}$ has a one-electron configuration $\ket{v\bar{v}}$, while the doubly-excited state $\ket{D}$ has a configuration $\ket{c\bar{c}}$.
There is then only one single excitation possible which corresponds to the transition $v \to c$ with different spin-flip configurations.
As usual, this can produce a singlet singly-excited state $\ket{S} = (\ket{v\bar{c}} + \ket{c\bar{v}})/\sqrt{2}$, and a triplet singly-excited state $\ket{T} = (\ket{v\bar{c}} - \ket{c\bar{v}})/\sqrt{2}$. \cite{SzaboBook}
As usual, this produces a singlet singly-excited state $\ket{S} = (\ket{v\bar{c}} + \ket{c\bar{v}})/\sqrt{2}$, and a triplet singly-excited state $\ket{T} = (\ket{v\bar{c}} - \ket{c\bar{v}})/\sqrt{2}$. \cite{SzaboBook}
For the singlet manifold, the exact Hamiltonian in the basis of the (spin-adapted) configuration state functions reads
For the singlet manifold, the exact Hamiltonian in the basis of these (spin-adapted) configuration state functions reads
\begin{equation} \label{eq:H-exact}
\bH^{\updw} =
\begin{pmatrix}
@ -227,7 +227,7 @@ and $\Delta\e{} = \e{c} - \e{v}$.
The energy of the only triplet state is simply $\mel{T}{\hH}{T} = \EHF + \Delta\e{} - \ERI{vv}{cc}$.
For the sake of illustration, we will use the same numerical example throughout this study, and consider the singlet ground state of the \ce{He} atom in Pople's 6-31G basis set.
This system contains two orbitals and the numerical values of the various quantities defined above are
The numerical values of the various quantities defined above are
\begin{subequations}
\begin{align}
\e{v} & = -0.914\,127
@ -351,7 +351,7 @@ However, following Strinati's footsteps, \cite{Strinati_1982,Strinati_1984,Strin
Based on the very same two-level model that we employ here, Romaniello and coworkers \cite{Romaniello_2009b} clearly evidenced that one can genuinely access additional excitations by solving the non-linear, frequency-dependent BSE eigenvalue problem.
For this particular system, they showed that a BSE kernel based on the random-phase approximation (RPA) produces indeed double excitations but also unphysical excitations, \cite{Romaniello_2009b} attributed to the self-screening problem. \cite{Romaniello_2009a}
This issue was resolved in the subsequent work of Sangalli \textit{et al.} \cite{Sangalli_2011} via the design of a diagrammatic number-conserving approach based on the folding of the second-RPA Hamiltonian. \cite{Wambach_1988}
Thanks to a careful diagrammatic analysis of the dynamic kernel, they showed that their approach produces the correct number of optically active poles, and this was further illustrated by computing the polarizability of two unsaturated hydrocarbon chains (\ce{C8H2} and \ce{C4H6}).
Thanks to a careful diagrammatic analysis of the dynamical kernel, they showed that their approach produces the correct number of optically active poles, and this was further illustrated by computing the polarizability of two unsaturated hydrocarbon chains (\ce{C8H2} and \ce{C4H6}).
Within the so-called $GW$ approximation of MBPT, \cite{Golze_2019} one can easily compute the quasiparticle energies associated with the valence and conduction orbitals.
Assuming that $W$ has been calculated at the random-phase approximation (RPA) level and within the TDA, the expression of the $\GW$ quasiparticle energy is
@ -398,9 +398,9 @@ with
and
\begin{subequations}
\begin{gather}
W^{\co}_R(\omega) = \frac{4 \ERI{vv}{vc} \ERI{vc}{cc}}{\omega - \Omega - \Delta\eGW{}}
W^{\co}_R(\omega) = \frac{4 \ERI{vv}{vc} \ERI{vc}{cc}}{\omega - \Omega - \Delta\eGW{}}
\\
W^{\co}_C(\omega) = \ERI{vc}{cv} + \frac{4 \ERI{vc}{cv}^2}{\omega - \Omega}
W^{\co}_C(\omega) = \frac{4 \ERI{vc}{cv}^2}{\omega - \Omega}
\end{gather}
\end{subequations}
are the elements of the correlation part of the dynamically-screened Coulomb potential for the resonant and coupling blocks of the dBSE Hamiltonian.
@ -435,7 +435,7 @@ As mentioned in Ref.~\onlinecite{Romaniello_2009b}, spurious solutions appears d
Indeed, diagonalizing the exact Hamiltonian \eqref{eq:H-exact} produces only two singlet solutions corresponding to the singly- and doubly-excited states, and one triplet state (see Sec.~\ref{sec:exact}).
Therefore, there is one spurious solution for the singlet manifold ($\omega_{2}^{\dBSE,\updw}$) and two spurious solutions for the triplet manifold ($\omega_{2}^{\dBSE,\upup}$ and $\omega_{3}^{\dBSE,\upup}$).
It is worth mentioning that, around $\omega = \omega_1^{\dBSE,\sigma}$, the slope of the curves depicted in Fig.~\ref{fig:dBSE} is small, while the two other solutions, $\omega_2^{\dBSE,\sigma}$ and $\omega_3^{\dBSE,\sigma}$, stem from poles and consequently the slope is very large around these frequency values.
\titou{T2: add comment on how one can detect fake solutions?}
This makes these two latter solutions quite hard to locate with the Newton-Raphson method.
%%% TABLE I %%%
\begin{table*}
@ -473,6 +473,7 @@ It is worth mentioning that, around $\omega = \omega_1^{\dBSE,\sigma}$, the slop
Enforcing the TDA, which corresponds to neglecting the coupling term between the resonant and anti-resonant part of the dBSE Hamiltonian \eqref{eq:HBSE}, allows to remove some of these spurious excitations.
There is thus only one spurious excitation in the triplet manifold ($\omega_{3}^{\BSE,\upup}$), the two solutions of the singlet manifold corresponding now to the single and double excitations.
The spin blindness of the dBSE kernel is probably to blame for the existence of this spurious triplet excitation.
Figure \ref{fig:dBSE-TDA} shows the same curves as Fig.~\ref{fig:dBSE} but in the TDA.
%%% FIGURE 2 %%%
@ -487,7 +488,7 @@ Figure \ref{fig:dBSE-TDA} shows the same curves as Fig.~\ref{fig:dBSE} but in th
%%% %%% %%% %%%
In the static approximation, only one solution per spin manifold is obtained by diagonalizing $\bH_{\BSE}^{\sigma}$ (see Fig.~\ref{fig:dBSE} and Table \ref{tab:BSE}).
Therefore, the static BSE Hamiltonian does not produce spurious excitations but misses the (singlet) double excitation, and it shows that the physical single excitation stemming from the dBSE Hamiltonian is the lowest one for each spin manifold, \ie, $\omega_1^{\dBSE,\updw}$ and $\omega_1^{\dBSE,\upup}$.
Therefore, the static BSE Hamiltonian misses the (singlet) double excitation (as it should), and it shows that the physical single excitation stemming from the dBSE Hamiltonian is the lowest in energy for each spin manifold, \ie, $\omega_1^{\dBSE,\updw}$ and $\omega_1^{\dBSE,\upup}$.
Another way to access dynamical effects while staying in the static framework is to use perturbation theory, \cite{Rohlfing_2000,Ma_2009a,Ma_2009b,Baumeier_2012b} a scheme we label as perturbative BSE (pBSE).
To do so, one must decompose the dBSE Hamiltonian into a (zeroth-order) static part and a dynamical perturbation, such that
@ -617,17 +618,18 @@ and
f_{\dBSE2}^{\co,\upup} = - \frac{\ERI{cc}{cc} \ERI{vc}{cv} + \ERI{vv}{vv} \ERI{vc}{cv} }{2 \Delta\eGF{}}
\end{gather}
\end{subequations}
Note that, unlike the dBSE Hamiltonian [see Eq.~\eqref{eq:HBSE}], the BSE2 dynamical kernel is spin-aware with distinct expressions for singlets and triplets, and the coupling block $C_{\dBSE2}^{\sigma}$ is frequency independent.
Note that, unlike the dBSE Hamiltonian [see Eq.~\eqref{eq:HBSE}], the BSE2 dynamical kernel is spin aware with distinct expressions for singlets and triplets, and the coupling block $C_{\dBSE2}^{\sigma}$ is frequency independent.
This latter point has an important consequence as this lack of frequency dependence removes one of the spurious pole (see Fig.~\ref{fig:BSE2}).
The singlet manifold has then the right number of excitations.
However, one spurious triplet excitation remains.
It is mentioned in Ref.~\onlinecite{Rebolini_2016} that the BSE2 kernel has some similarities with the second-order polarization-propagator approximation \cite{Oddershede_1977,Nielsen_1980} (SOPPA) and second RPA kernels. \cite{Huix-Rotllant_2011,Huix-Rotllant_PhD,Sangalli_2011}
Numerical results for the two-level model are reported in Table \ref{tab:BSE2} with the usual approximations and perturbative treatments.
In the case of BSE2, the perturbative partitioning is simply
\begin{equation}
\bH_{\dBSE2}^{\sigma}(\omega)
= \underbrace{\bH_{\BSE2}^{\sigma}}_{\bH_{\pBSE2}^{(0)}}
+ \underbrace{\qty[ \bH_{\dBSE2}^{\sigma}(\omega) - \bH_{\BSE2}^{\sigma} ]}_{\bH_{\pBSE2}^{(1)}}
+ \underbrace{\qty[ \bH_{\dBSE2}^{\sigma}(\omega) - \bH_{\BSE2}^{\sigma} ]}_{\bH_{\pBSE2}^{(1)}(\omega)}
\end{equation}
%%% TABLE II %%%
@ -669,53 +671,53 @@ For the double excitation, dBSE2 yields a slightly better energy, yet still in q
%%% %%% %%% %%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\subsection{The forgotten kernel: Sangalli's kernel}
\label{sec:Sangalli}
%\subsection{The forgotten kernel: Sangalli's kernel}
%\label{sec:Sangalli}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\titou{This section is experimental...}
In Ref.~\onlinecite{Sangalli_2011}, Sangalli proposed a dynamical kernel (based on the second RPA) without (he claims) spurious excitations thanks to the design of a number-conserving approach which correctly describes particle indistinguishability and Pauli exclusion principle.
We will first start by writing down explicitly this kernel as it is given in obscure physicist notations in the original article.
The Hamiltonian with Sangalli's kernel is (I think)
\begin{equation}
\bH_\text{S}^{\sigma}(\omega) =
\begin{pmatrix}
\bR_\text{S}^{\sigma}(\omega) & \bC_\text{S}^{\sigma}(\omega)
\\
-\bC_\text{S}^{\sigma}(-\omega) & -\bR_\text{S}^{\sigma}(-\omega)
\end{pmatrix}
\end{equation}
with
\begin{subequations}
\begin{gather}
R_{ia,jb}^{\sigma}(\omega) = \delta_{ij} \delta_{ab} (\eGW{a} - \eGW{i}) + f_{ia,jb}^{\sigma} (\omega)
\\
C_{ia,jb}^{\sigma}(\omega) = f_{ia,bj}^{\sigma} (\omega)
\end{gather}
\end{subequations}
and
\begin{subequations}
\begin{gather}
f_{ia,jb}^{\sigma} (\omega) = \sum_{m \neq n} \frac{ c_{ia,mn} c_{jb,mn} }{\omega - ( \omega_{m} + \omega_{n})}
\\
c_{ia,mn}^{\sigma} = \sum_{jb,kc} \qty{ \qty[ \ERI{ij}{kc} \delta_{ab} + \ERI{kc}{ab} \delta_{ij} ] \qty[ R_{m,jc} R_{n,kb}
+ R_{m,kb} R_{n,jc} ] }
\end{gather}
\end{subequations}
where $R_{m,ia}$ are the elements of the RPA eigenvectors.
For the two-level model, Sangalli's kernel reads
\begin{align}
R(\omega) & = \Delta\eGW{} + f_R (\omega)
\\
C(\omega) & = f_C (\omega)
\end{align}
\begin{gather}
f_R (\omega) = 2 \frac{ [\ERI{vv}{vc} + \ERI{vc}{cc}]^2 }{\omega - 2\omega_1}
\\
f_C (\omega) = 0
\end{gather}
%\titou{This section is experimental...}
%In Ref.~\onlinecite{Sangalli_2011}, Sangalli proposed a dynamical kernel (based on the second RPA) without (he claims) spurious excitations thanks to the design of a number-conserving approach which correctly describes particle indistinguishability and Pauli exclusion principle.
%We will first start by writing down explicitly this kernel as it is given in obscure physicist notations in the original article.
%
%The Hamiltonian with Sangalli's kernel is (I think)
%\begin{equation}
% \bH_\text{S}^{\sigma}(\omega) =
% \begin{pmatrix}
% \bR_\text{S}^{\sigma}(\omega) & \bC_\text{S}^{\sigma}(\omega)
% \\
% -\bC_\text{S}^{\sigma}(-\omega) & -\bR_\text{S}^{\sigma}(-\omega)
% \end{pmatrix}
%\end{equation}
%with
%\begin{subequations}
%\begin{gather}
% R_{ia,jb}^{\sigma}(\omega) = \delta_{ij} \delta_{ab} (\eGW{a} - \eGW{i}) + f_{ia,jb}^{\sigma} (\omega)
% \\
% C_{ia,jb}^{\sigma}(\omega) = f_{ia,bj}^{\sigma} (\omega)
%\end{gather}
%\end{subequations}
%and
%\begin{subequations}
%\begin{gather}
% f_{ia,jb}^{\sigma} (\omega) = \sum_{m \neq n} \frac{ c_{ia,mn} c_{jb,mn} }{\omega - ( \omega_{m} + \omega_{n})}
% \\
% c_{ia,mn}^{\sigma} = \sum_{jb,kc} \qty{ \qty[ \ERI{ij}{kc} \delta_{ab} + \ERI{kc}{ab} \delta_{ij} ] \qty[ R_{m,jc} R_{n,kb}
% + R_{m,kb} R_{n,jc} ] }
%\end{gather}
%\end{subequations}
%where $R_{m,ia}$ are the elements of the RPA eigenvectors.
%
%For the two-level model, Sangalli's kernel reads
%\begin{align}
% R(\omega) & = \Delta\eGW{} + f_R (\omega)
% \\
% C(\omega) & = f_C (\omega)
%\end{align}
%
%\begin{gather}
% f_R (\omega) = 2 \frac{ [\ERI{vv}{vc} + \ERI{vc}{cc}]^2 }{\omega - 2\omega_1}
% \\
% f_C (\omega) = 0
%\end{gather}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Take-home messages}