PTEROSOR/BSEdyn
10
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

done for 1st iteration

Browse Source
This commit is contained in:
Pierre-Francois Loos 2020-06-22 14:24:05 +02:00
parent 26cfdab924
commit 751aab3bb4
7 changed files with 212 additions and 178 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
BSEdyn.bib
View File

@ -1,13 +1,21 @@
%% This BibTeX bibliography file was created using BibDesk.
%% http://bibdesk.sourceforge.net/
%% Created for Pierre-Francois Loos at 2020-06-21 22:26:07 +0200
%% Created for Pierre-Francois Loos at 2020-06-22 10:36:09 +0200
%% Saved with string encoding Unicode (UTF-8)
@article{Loos_2020e,
Author = {X. Blase and Y. Duchemin and D. Jacquemin},
Date-Added = {2020-06-22 09:07:38 +0200},
Date-Modified = {2020-06-22 09:08:39 +0200},
Journal = {J. Phys. Chem. Lett.},
Title = {The Bethe-Salpeter Equation Formalism: From Physics to Chemistry},
Year = {submitted}}
@article{Loos_2020d,
Author = {P. F. Loos and E. Fromager},
Date-Added = {2020-06-21 21:41:42 +0200},
@ -1968,16 +1976,10 @@
@article{Romaniello_2009a,
Author = {Romaniello, P. and Guyot, S. and Reining, L.},
Date-Added = {2020-05-18 21:40:28 +0200},
Date-Modified = {2020-05-18 21:40:28 +0200},
Date-Modified = {2020-06-22 09:34:56 +0200},
Doi = {10.1063/1.3249965},
File = {/Users/loos/Zotero/storage/WCWVRT3M/Romaniello_2009.pdf},
Issn = {0021-9606, 1089-7690},
Journal = {J. Chem. Phys.},
Language = {en},
Month = oct,
Number = {15},
Pages = {154111},
Shorttitle = {The Self-Energy beyond {{GW}}},
Title = {The Self-Energy beyond {{GW}}: {{Local}} and Nonlocal Vertex Corrections},
Volume = {131},
Year = {2009},

2
Notes/BSEdyn-notes.rty
View File

@ -43,10 +43,12 @@
\newcommand{\HF}{\text{HF}}
\newcommand{\RPA}{\text{RPA}}
\newcommand{\RPAx}{\text{RPAx}}
\newcommand{\TDHF}{\text{TDHF}}
\newcommand{\dRPAx}{\text{dRPAx}}
\newcommand{\BSE}{\text{BSE}}
\newcommand{\TDABSE}{\text{BSE(TDA)}}
\newcommand{\dBSE}{\text{dBSE}}
\newcommand{\pBSE}{\text{pBSE}}
\newcommand{\TDAdBSE}{\text{dBSE(TDA)}}
\newcommand{\GW}{GW}
\newcommand{\GF}{\text{GF2}}

368
Notes/BSEdyn-notes.tex
View File

@ -67,9 +67,8 @@ where the explicit expressions of the resonant and coupling blocks, $\bR^{\sigma
Neglecting the coupling block [\ie, $\bC^{\sigma}(\omega) = 0$] between the resonant and anti-resonants parts, $\bR(\omega)$ and $-\bR(-\omega)$, is known as the Tamm-Dancoff approximation (TDA).
The non-linear eigenvalue problem defined in Eq.~\eqref{eq:LR} has particle-hole symmetry which means that it is invariant via the transformation $\omega \to -\omega$.
Therefore, without loss of generality, we will restrict our analysis to positive frequencies.
The central point here is that, thanks to their non-linear nature stemming from their frequency dependence, dynamical kernels potentially generate more than just single excitations.
In the one-electron basis of (real) spatial orbitals $\lbrace \MO{p} \rbrace$, the elements of the matrices defined in Eq.~\eqref{eq:LR} have the following generic forms: \cite{Dreuw_2005}
In the one-electron basis of (real) spatial orbitals $\lbrace \MO{p} \rbrace$, we will assume that the elements of the matrices defined in Eq.~\eqref{eq:LR} have the following generic forms: \cite{Dreuw_2005}
\begin{subequations}
\begin{gather}
R_{ia,jb}^{\sigma}(\omega) = (\e{a} - \e{i}) \delta_{ij} \delta_{ab} + f_{ia,jb}^{\Hxc,\sigma}(\omega)
@ -82,22 +81,23 @@ where $\delta_{pq}$ is the Kronecker delta, $\e{p}$ is the one-electron energy a
f_{ia,jb}^{\Hxc,\sigma}(\omega)
= \iint \MO{i}(\br) \MO{a}(\br) f^{\Hxc,\sigma}(\omega) \MO{j}(\br') \MO{b}(\br') d\br d\br'
\end{equation}
Here, $i$ and $j$ are occupied orbitals, $a$ and $b$ are unoccupied orbitals, and $p$ and $q$ indicate arbitrary orbitals.
In Eq.~\eqref{eq:kernel},
\begin{equation} \label{eq:kernel-Hxc}
f^{\Hxc,\sigma}(\omega) = f^{\Hx,\sigma} + f^{\co,\sigma}(\omega)
\end{equation}
is the (spin-resolved) Hartree-exchange-correlation (Hxc) dynamical kernel.
Here, $i$ and $j$ are occupied orbitals, $a$ and $b$ are unoccupied orbitals, and $p$ and $q$ indicate arbitrary orbitals.
In the case of a spin-independent kernel, we will drop the superscript $\sigma$.
As readily seen from Eq.~\eqref{eq:kernel-Hxc}, only the correlation (c) part of the kernel is frequency dependent and, in a wave function context, the static Hartree-exchange (Hx) matrix elements read
\begin{equation}
f_{ia,jb}^{\Hx,\sigma} = \sigma \ERI{ia}{jb} - \ERI{ib}{ja}
f_{ia,jb}^{\Hx,\sigma} = 2\sigma \ERI{ia}{jb} - \ERI{ib}{ja}
\end{equation}
where $\sigma = 1 $ or $0$ for singlet and triplet excited states (respectively), and
\begin{equation}
\ERI{ia}{jb} = \iint \MO{i}(\br) \MO{a}(\br) \frac{1}{\abs{\br - \br'}} \MO{j}(\br') \MO{b}(\br') d\br d\br'
\end{equation}
are the usual (bare) two-electron integrals.
In the case of a spin-independent kernel, we will drop the superscript $\sigma$.
are the usual two-electron integrals.
The central point here is that, thanks to its non-linear nature stemming from their frequency dependence, a dynamical kernel potentially generates more than just single excitations.
Unless otherwise stated, atomic units are used and we assume real quantities throughout this manuscript.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
@ -135,7 +135,7 @@ However, in practice, $N$ might be (very) large (\eg, equal to the total number
where the blocks $\bA_1$ and $\bA_2$, of sizes $N_1 \times N_1$ and $N_2 \times N_2$ (with $N_1 + N_2 = N$), can be associated with, for example, the single and double excitations of the system.
Note that this \textit{exact} decomposition does not alter, in any case, the values of the excitation energies, not their eigenvectors.
Solving separately each row of the system \eqref{eq:lin_sys_split} and assuming that $\omega \bI - \bA_2$ is invertible yields
Solving separately each row of the system \eqref{eq:lin_sys_split} and assuming that $\omega \bI - \bA_2$ is invertible, it follows that
\begin{subequations}
\begin{gather}
\label{eq:row1}
@ -160,7 +160,7 @@ For example, an operator $\Tilde{\bA}_1(\omega)$ built in the single-excitation
How have we been able to reduce the dimension of the problem while keeping the same number of solutions?
To do so, we have transformed a linear operator $\bA$ into a non-linear operator $\Tilde{\bA}_1(\omega)$ by making it frequency dependent.
In other words, we have sacrificed the linearity of the system in order to obtain a new, non-linear systems of equations of smaller dimension.
This procedure converting degrees of freedom into frequency or energy dependence is very general and can be applied in various contexts. \cite{Garniron_2018,QP2}
This procedure converting degrees of freedom into frequency or energy dependence is very general and can be applied in various contexts. \cite{Sottile_2003,Garniron_2018,QP2}
Thanks to its non-linearity, Eq.~\eqref{eq:non_lin_sys} can produce more solutions than its actual dimension.
However, because there is no free lunch, this non-linear system is obviously harder to solve than its corresponding linear analog given by Eq.~\eqref{eq:lin_sys}.
Nonetheless, approximations can be now applied to Eq.~\eqref{eq:non_lin_sys} in order to solve it efficiently.
@ -185,7 +185,7 @@ 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 \cite{Romaniello_2009b} in its singlet ground state (\ie, the lowest orbital is doubly occupied).
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}
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}}$.
@ -259,7 +259,9 @@ We are going to use these as reference for the remaining of this study.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\subsection{Maitra's dynamical kernel}
\label{sec:Maitra}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
The kernel proposed by Maitra and coworkers \cite{Maitra_2004,Cave_2004} in the context of dressed TDDFT (D-TDDFT) corresponds to an \textit{ad hoc} many-body theory correction to TDDFT.
More specifically, D-TDDFT adds to the static kernel a frequency-dependent part by reverse-engineering the exact Hamiltonian: a single and double excitations, assumed to be strongly coupled, are isolated from among the spectrum and added manually to the static kernel.
The very same idea was taking further by Huix-Rotllant, Casida and coworkers, \cite{Huix-Rotllant_2011} and tested on a large set of molecules.
@ -268,9 +270,9 @@ The static problem corresponds then to the TDHF Hamiltonian, while in the TDA, i
For the two-level model, the reverse-engineering process of the exact Hamiltonian \eqref{eq:H-exact} yields
\begin{equation} \label{eq:f-Maitra}
f_M^{\co,\updw}(\omega) = \frac{\abs*{\mel{S}{\hH}{D}}^2}{\omega - (\mel{D}{\hH}{D} - \mel{0}{\hH}{0}) }
f_\text{M}^{\co,\updw}(\omega) = \frac{\abs*{\mel{S}{\hH}{D}}^2}{\omega - (\mel{D}{\hH}{D} - \mel{0}{\hH}{0}) }
\end{equation}
while $f_M^{\co,\upup}(\omega) = 0$.
while $f_\text{M}^{\co,\upup}(\omega) = 0$.
The expression \eqref{eq:f-Maitra} can be easily obtained by folding the double excitation onto the single excitation, as explained in Sec.~\ref{sec:dyn}.
It is clear that one must know \textit{a priori} the structure of the Hamiltonian to construct such dynamical kernel, and this obviously hampers its applicability to realistic photochemical systems where it is sometimes hard to get a clear picture of the interplay between excited states. \cite{Boggio-Pasqua_2007}
@ -278,30 +280,30 @@ For the two-level model, the non-linear equations defined in Eq.~\eqref{eq:LR} p
\begin{equation} \label{eq:H-M}
\bH_\text{D-TDHF}^{\sigma}(\omega) =
\begin{pmatrix}
R_M^{\sigma}(\omega) & C_M^{\sigma}(\omega)
R_\text{M}^{\sigma}(\omega) & C_\text{M}^{\sigma}(\omega)
\\
-C_M^{\sigma}(-\omega) & -R_M^{\sigma}(-\omega)
-C_\text{M}^{\sigma}(-\omega) & -R_\text{M}^{\sigma}(-\omega)
\end{pmatrix}
\end{equation}
with
\begin{subequations}
\begin{gather}
\label{eq:R_M}
R_M^{\sigma}(\omega) = \Delta\e{} + 2 \sigma \ERI{vc}{vc} - \ERI{vc}{vc} + f_M^{\co,\sigma}(\omega)
R_\text{M}^{\sigma}(\omega) = \Delta\e{} + 2 \sigma \ERI{vc}{vc} - \ERI{vc}{vc} + f_\text{M}^{\co,\sigma}(\omega)
\\
\label{eq:C_M}
C_M^{\sigma}(\omega) = 2 \sigma \ERI{vc}{cv} - \ERI{vv}{cc} + f_M^{\co,\sigma}(\omega)
C_\text{M}^{\sigma}(\omega) = 2 \sigma \ERI{vc}{cv} - \ERI{vv}{cc} + f_\text{M}^{\co,\sigma}(\omega)
\end{gather}
\end{subequations}
yielding the excitation energies reported in Table \ref{tab:Maitra} when diagonalized.
The TDHF Hamiltonian is obtained from Eq.~\eqref{eq:H-M} by setting $f_M^{\co,\sigma}(\omega) = 0$ in Eqs.~\eqref{eq:R_M} and \eqref{eq:C_M}.
The TDHF Hamiltonian is obtained from Eq.~\eqref{eq:H-M} by setting $f_\text{M}^{\co,\sigma}(\omega) = 0$ in Eqs.~\eqref{eq:R_M} and \eqref{eq:C_M}.
In Fig.~\ref{fig:Maitra}, we plot $\det[\bH(\omega) - \omega \bI]$ as a function of $\omega$ for both the singlet (black and gray) and triplet (orange) manifolds.
The roots of $\det[\bH(\omega) - \omega \bI]$ indicate the excitation energies.
Because, there is nothing to dress for the triplet state, only the static TDHF excitation energy is reported.
%%% TABLE I %%%
\begin{table}
\caption{Singlet and triplet excitation energies at various levels of theory.
\caption{Singlet and triplet excitation energies (in hartree) at various levels of theory.
\label{tab:Maitra}
}
\begin{ruledtabular}
@ -309,7 +311,7 @@ Because, there is nothing to dress for the triplet state, only the static TDHF e
Singlets & CIS & TDHF & D-CIS & D-TDHF & Exact \\
\hline
$\omega_1^{\updw}$ & 1.91119 & 1.89758 & 1.90636 & 1.89314 & 1.92145 \\
$\omega_3^{\updw}$ & & & & 3.44865 & 3.47880 \\
$\omega_3^{\updw}$ & & & 3.44888 & 3.44865 & 3.47880 \\
\hline
Triplets & & & & & Exact \\
\hline
@ -323,7 +325,7 @@ Because, there is nothing to dress for the triplet state, only the static TDHF e
\begin{figure}
\includegraphics[width=\linewidth]{Maitra}
\caption{
$\det[\bH(\omega) - \omega \bI]$ as a function of $\omega$ for both the singlet (black and gray) and triplet (orange) manifolds.
$\det[\bH(\omega) - \omega \bI]$ as a function of $\omega$ for both the singlet (gray and black) and triplet (orange) manifolds.
The static TDHF Hamiltonian (dashed) and dynamic D-TDHF Hamiltonian (solid) are considered.
\label{fig:Maitra}
}
@ -331,26 +333,33 @@ Because, there is nothing to dress for the triplet state, only the static TDHF e
%%% %%% %%% %%%
Although not particularly accurate for the single excitations, Maitra's dynamical kernel allows to access the double excitation with good accuracy and provides exactly the right number of solutions (two singlets and one triplet).
Note that the correlation kernel is known to work best in the weak correlation regime (which is the case here) where the true excitations have a clear single and double excitation character, \cite{Loos_2019,Loos_2020d} but it is not intended to explore strongly correlated systems. \cite{Carrascal_2018}
Note that this correlation kernel is known to work best in the weak correlation regime (which is the case here) where the true excitations have a clear single and double excitation character, \cite{Loos_2019,Loos_2020d} but it is not intended to explore strongly correlated systems. \cite{Carrascal_2018}
Its accuracy for the single excitations could be certainly improved in a DFT context.
However, this is not the point of the present investigation.
Table \ref{tab:Maitra} also reports the slightly-improves (thanks to error compensation) CIS and D-CIS excitation energies.
Table \ref{tab:Maitra} also reports the slightly-improved (thanks to error compensation) CIS and D-CIS excitation energies.
In particular, single excitations are greatly improved without altering the accuracy of the double excitation.
Graphically, the curves obtained for CIS and D-CIS are extremely similar to the ones of TDHF and D-TDHF depicted in Fig.~\ref{fig:Maitra}.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\subsection{Dynamical BSE kernel}
\label{sec:BSE}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\titou{Titou will add a bit more background information about the BSE dynamical kernel.}
As mentioned in Sec.~\ref{sec:dyn}, most of BSE calculations performed nowadays are done within the static approximation. \cite{ReiningBook,Loos_2020e}
However, following Strinati's footsteps, \cite{Strinati_1982,Strinati_1984,Strinati_1988} several groups have explored this formalism beyond the static approximation by retaining (or reviving) the dynamical nature of the dynamically-screened Coulomb potential $W$ \cite{Sottile_2003,Romaniello_2009b,Sangalli_2011} or via a perturbative approach. \cite{Rohlfing_2000,Ma_2009a,Ma_2009b,Baumeier_2012b}
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}).
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 the dynamically-screened Coulomb potential $W$ has been calculated at the random-phase approximation (RPA) level and within the TDA, the expression of the $\GW$ quasiparticle energy is
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
\begin{equation}
\e{p}^{\GW} = \e{p} + Z_{p}^{\GW} \SigGW{p}(\e{p})
\end{equation}
where $p = v$ or $c$,
\begin{equation}
\label{eq:SigC}
\label{eq:SigGW}
\SigGW{p}(\omega) = \frac{2 \ERI{pv}{vc}^2}{\omega - \e{v} + \Omega} + \frac{2 \ERI{pc}{cv}^2}{\omega - \e{c} - \Omega}
\end{equation}
is the correlation part of the self-energy $\Sig{}$, and
@ -358,7 +367,7 @@ is the correlation part of the self-energy $\Sig{}$, and
Z_{p}^{\GW} = \qty( 1 - \left. \pdv{\SigGW{p}(\omega)}{\omega} \right|_{\omega = \e{p}} )^{-1}
\end{equation}
is the renormalization factor.
In Eq.~\eqref{eq:SigC}, $\Omega = \Delta\eGW{} + 2 \ERI{vc}{vc}$ is the sole (singlet) RPA excitation energy of the system, with $\Delta\eGW{} = \eGW{c} - \eGW{v}$.
In Eq.~\eqref{eq:SigGW}, $\Omega = \Delta\eGW{} + 2 \ERI{vc}{vc}$ is the sole (singlet) RPA excitation energy of the system, with $\Delta\eGW{} = \eGW{c} - \eGW{v}$.
Numerically, we get
\begin{align}
\Omega & = 2.769\,327
@ -368,7 +377,7 @@ Numerically, we get
\eGW{c} & = +1.373\,640
\end{align}
One can now build the dynamical BSE (dBSE) Hamiltonian \cite{Strinati_1988}
One can now build the dynamical BSE (dBSE) Hamiltonian \cite{Strinati_1988,Romaniello_2009b}
\begin{equation} \label{eq:HBSE}
\bH_{\dBSE}^{\sigma}(\omega) =
\begin{pmatrix}
@ -396,52 +405,7 @@ and
are the elements of the correlation part of the dynamically-screened Coulomb potential for the resonant and coupling blocks of the dBSE Hamiltonian.
Note that, in this case, the correlation kernel is spin blind.
It can be easily shown that solving the equation
\begin{equation}
\det[\bH_{\dBSE}^{\sigma}(\omega) - \omega \bI] = 0
\end{equation}
yields 3 solutions per spin manifold (see Fig.~\ref{fig:dBSE}).
Their numerical values are reported in Table \ref{tab:dBSE} alongside other variants discussed below.
As mentioned in Ref.~\onlinecite{Romaniello_2009b}, spurious solutions appears due to the approximate nature of the dBSE kernel.
Indeed, diagonalizing the exact Hamiltonian 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.
This evidences that dBSE reproduces qualitatively well the singlet and triplet single excitations, but quite badly the double excitation which is off by more than 1 hartree.
%%% TABLE I %%%
\begin{table*}
\caption{BSE singlet and triplet excitation energies at various levels of theory.
\label{tab:dBSE}
}
\begin{ruledtabular}
\begin{tabular}{|c|ccccccc|c|}
Singlets & BSE & BSE1 & BSE1(dTDA) & dBSE & BSE(TDA) & BSE1(TDA) & dBSE(TDA) & Exact \\
\hline
$\omega_1^{\updw}$ & 1.92778 & 1.90022 & 1.91554 & 1.90527 & 1.95137 & 1.94004 & 1.94005 & 1.92145 \\
$\omega_2^{\updw}$ & & & & 2.78377 & & & & \\
$\omega_3^{\updw}$ & & & & 4.90134 & & & 4.90117 & 3.47880 \\
\hline
Triplets & BSE & BSE1 & BSE1(dTDA) & dBSE & BSE(TDA) & BSE1(TDA) & dBSE(TDA) & Exact \\
\hline
$\omega_1^{\upup}$ & 1.48821 & 1.46860 & 1.46260 & 1.46636 & 1.49603 & 1.47070 & 1.47070 & 1.47085 \\
$\omega_2^{\upup}$ & & & & 2.76178 & & & & \\
$\omega_3^{\upup}$ & & & & 4.91545 & & & 4.91517 & \\
\end{tabular}
\end{ruledtabular}
\end{table*}
%%% %%% %%% %%%
%%% FIGURE 2 %%%
\begin{figure}
\includegraphics[width=\linewidth]{dBSE}
\caption{
$\det[\bH^{\dBSE}(\omega) - \omega \bI]$ as a function of $\omega$ for both the singlet (red) and triplet (blue) manifolds.
\label{fig:dBSE}
}
\end{figure}
%%% %%% %%% %%%
Within the usual static approximation, the BSE Hamiltonian is
Within the usual static approximation, the BSE Hamiltonian is simply
\begin{equation}
\bH_{\BSE}^{\sigma} =
\begin{pmatrix}
@ -458,41 +422,90 @@ with
C_{\BSE}^{\sigma} = 2 \sigma \ERI{vc}{vc} - W_C(\omega = 0)
\end{gather}
\end{subequations}
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:dBSE}).
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}$.
Enforcing the TDA, which corresponds to neglecting the coupling term between the resonant and anti-resonant part of the dBSE Hamiltonian, 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 to the single and double excitations.
It can be easily shown that solving the equation
\begin{equation}
\det[\bH_{\dBSE}^{\sigma}(\omega) - \omega \bI] = 0
\end{equation}
yields 3 solutions per spin manifold (see Fig.~\ref{fig:dBSE}).
Their numerical values are reported in Table \ref{tab:BSE} alongside other variants discussed below.
This evidences that dBSE reproduces qualitatively well the singlet and triplet single excitations, but quite badly the double excitation which is off by more than 1 hartree.
As mentioned in Ref.~\onlinecite{Romaniello_2009b}, spurious solutions appears due to the approximate nature of the dBSE kernel.
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.
%%% TABLE I %%%
\begin{table*}
\caption{BSE singlet and triplet excitation energies (in hartree) at various levels of theory.
\label{tab:BSE}
}
\begin{ruledtabular}
\begin{tabular}{|c|ccccccc|c|}
Singlets & BSE & pBSE & pBSE(dTDA) & dBSE & BSE(TDA) & pBSE(TDA) & dBSE(TDA) & Exact \\
\hline
$\omega_1^{\updw}$ & 1.92778 & 1.90022 & 1.91554 & 1.90527 & 1.95137 & 1.94004 & 1.94005 & 1.92145 \\
$\omega_2^{\updw}$ & & & & 2.78377 & & & & \\
$\omega_3^{\updw}$ & & & & 4.90134 & & & 4.90117 & 3.47880 \\
\hline
Triplets & BSE & pBSE & pBSE(dTDA) & dBSE & BSE(TDA) & pBSE(TDA) & dBSE(TDA) & Exact \\
\hline
$\omega_1^{\upup}$ & 1.48821 & 1.46860 & 1.46260 & 1.46636 & 1.49603 & 1.47070 & 1.47070 & 1.47085 \\
$\omega_2^{\upup}$ & & & & 2.76178 & & & & \\
$\omega_3^{\upup}$ & & & & 4.91545 & & & 4.91517 & \\
\end{tabular}
\end{ruledtabular}
\end{table*}
%%% %%% %%% %%%
%%% FIGURE 2 %%%
\begin{figure}
\includegraphics[width=\linewidth]{dBSE}
\caption{
$\det[\bH(\omega) - \omega \bI]$ as a function of $\omega$ for both the singlet (gray and black) and triplet (orange and red) manifolds.
The static BSE Hamiltonian (dashed) and dynamic dBSE Hamiltonian (solid) are considered.
\label{fig:dBSE}
}
\end{figure}
%%% %%% %%% %%%
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.
Figure \ref{fig:dBSE-TDA} shows the same curves as Fig.~\ref{fig:dBSE} but in the TDA.
%%% FIGURE 2 %%%
\begin{figure}
\includegraphics[width=\linewidth]{dBSE-TDA}
\caption{
$\det[\bH^{\TDAdBSE}(\omega) - \omega \bI]$ as a function of $\omega$ for both the singlet (red) and triplet (blue) manifolds within the TDA.
$\det[\bH(\omega) - \omega \bI]$ as a function of $\omega$ for both the singlet (gray and black) and triplet (orange and red) manifolds within the TDA.
The static BSE Hamiltonian (dashed) and dynamic dBSE Hamiltonian (solid) are considered.
\label{fig:dBSE-TDA}
}
\end{figure}
%%% %%% %%% %%%
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}$.
Another way to access dynamical effects while staying in the static framework is to use perturbation theory.
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
\begin{equation}
\bH_{\dBSE}^{\sigma}(\omega) = \underbrace{\bH_{\BSE}^{\sigma}}_{\bH^{(0)}} + \underbrace{\qty[ \bH_{\dBSE}^{\sigma}(\omega) - \bH_{\BSE}^{\sigma} ]}_{\bH^{(1)}}
\bH_{\dBSE}^{\sigma}(\omega)
= \underbrace{\bH_{\BSE}^{\sigma}}_{\bH_{\pBSE}^{(0)}}
+ \underbrace{\qty[ \bH_{\dBSE}^{\sigma}(\omega) - \bH_{\BSE}^{\sigma} ]}_{\bH_{\pBSE}^{(1)}(\omega)}
\end{equation}
Thanks to (renormalized) first-order perturbation theory, one gets
\begin{equation}
\begin{split}
\omega_{1}^{\BSE1,\sigma}
\omega_{1}^{\pBSE,\sigma}
& = \omega_{1}^{\BSE,\sigma}
\\
& + Z_{1}
& + Z_{1}^{\pBSE}
\T{\begin{pmatrix}
X_1 \\ Y_1
\end{pmatrix}
}
\cdot \qty[ \bH^{\dBSE}(\omega = \omega_{1}^{\BSE,\sigma}) - \bH^{\BSE} ] \cdot
\cdot \qty[ \bH_{\dBSE}^{\sigma}(\omega = \omega_{1}^{\BSE,\sigma}) - \bH_{\BSE}^{\sigma} ] \cdot
\begin{pmatrix}
X_1 \\ Y_1
\end{pmatrix}
@ -500,7 +513,7 @@ Thanks to (renormalized) first-order perturbation theory, one gets
\end{equation}
where
\begin{equation}
\bH^{\BSE}
\bH_{\BSE}^{\sigma}
\cdot
\begin{pmatrix}
X_1 \\ Y_1
@ -510,33 +523,40 @@ where
X_1 \\ Y_1
\end{pmatrix}
\end{equation}
are the eigenvectors of $\bH^{\BSE}$, and
and the renormalization factor is
\begin{equation}
Z_{1} = \qty{ 1 -
Z_{1}^{\pBSE} = \qty{ 1 -
\T{
\begin{pmatrix}
X_1 \\ Y_1
\end{pmatrix}
}
\cdot \left. \pdv{\bH^{\dBSE}(\omega)}{\omega} \right|_{\omega = \omega_{1}^{\BSE,\sigma}} \cdot
\cdot \left. \pdv{\bH_{\dBSE}^{\sigma}(\omega)}{\omega} \right|_{\omega = \omega_{1}^{\BSE,\sigma}} \cdot
\begin{pmatrix}
X_1 \\ Y_1
\end{pmatrix}
}^{-1}
\end{equation}
This corresponds to a dynamical correction to the static excitations, and the TDA can be applied to the dynamical correction, a scheme we label as dTDA in the following.
This corresponds to a dynamical perturbative correction to the static excitations.
Obviously, the TDA can be applied to the dynamical correction as well, a scheme we label as dTDA in the following.
The perturbatively-corrected values are also reported in Table \ref{tab:dBSE}, which shows that this scheme is very efficient at reproducing the dynamical value.
Note that, although the BSE1(dTDA) value is further from the dBSE value than BSE1, it is quite close to the exact excitation energy.
The perturbatively-corrected values are also reported in Table \ref{tab:BSE}, which shows that this scheme is very efficient at reproducing the dynamical value for the single excitations.
However, because the perturbative treatment is ultimately static, one cannot access double excitations with such a scheme.
Note that, although the pBSE(dTDA) value is further from the dBSE value than pBSE, it is quite close to the exact excitation energy.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\subsection{Second-order BSE kernel}
\label{sec:BSE2}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Here, we follow a different strategy and compute the dynamical second-order BSE kernel as illustrated by Yang and collaborators \cite{Zhang_2013}, and Rebolini and Toulouse \cite{Rebolini_2016}.
The third and final dynamical kernel that we consider here is the second-order BSE (BSE2) kernel as derived by Yang and collaborators in the TDA, \cite{Zhang_2013} and by Rebolini and Toulouse in a range-separated context \cite{Rebolini_2016,Rebolini_PhD} (see also Refs.~\onlinecite{Myohanen_2008,Sakkinen_2012}).
Note that a beyond-TDA BSE2 kernel was also derived in Ref.~\onlinecite{Rebolini_2016}, but was not tested.
In a nutshell, the BSE2 scheme applies second-order perturbation theory to optical excitations within the Green's function framework by taking the functional derivative of the second-order self-energy $\SigGF{}$ with respect to the one-body Green's function.
Because $\SigGF{}$ is a proper functional derivative, it was claimed in Ref.~\onlinecite{Zhang_2013} that BSE2 does not produce spurious excitations.
However, as we will show below, this is not always true.
First, let us compute the second-order quasiparticle energies, which reads
Like BSE requires $GW$ quasiparticle energies, BSE2 requires the second-order Green's function (GF2) quasiparticle energies, \cite{SzaboBook} which are defined as follows:
\begin{equation}
\eGF{p} = \e{p} + Z_{p}^{\GF} \SigGF{p}(\e{p})
\end{equation}
@ -549,75 +569,78 @@ and
\begin{equation}
Z_{p}^{\GF} = \qty( 1 - \left. \pdv{\SigGF{p}(\omega)}{\omega} \right|_{\omega = \e{p}} )^{-1}
\end{equation}
This expression can be easily obtained in the present case by the substitution $\Omega \to \e{c} - \e{v}$ which transforms the $GW$ self-energy into its GF2 analog.
The expression of the GF2 self-energy \eqref{eq:SigGF} can be easily obtained from its $GW$ counterpart \eqref{eq:SigGW} via the substitution $\Omega \to \e{c} - \e{v}$.
The static Hamiltonian for this theory is just the usual RPAx (or TDHF) Hamiltonian, \ie,
The static Hamiltonian of BSE2 is just the usual TDHF Hamiltonian where one substitutes the HF orbital energies by the GF2 quasiparticle energies, \ie,
\begin{equation}
\bH^{\RPAx} =
\bH_{\BSE2}^{\sigma} =
\begin{pmatrix}
A^{\stat} & B^{\stat}
R_{\BSE2}^{\sigma} & C_{\BSE2}^{\sigma}
\\
-B^{\stat} & -A^{\stat}
-C_{\BSE2}^{\sigma} & -R_{\BSE2}^{\sigma}
\end{pmatrix}
\end{equation}
with
\begin{subequations}
\begin{gather}
A^{\stat} = \Delta\eGF{} + 2 \sigma \ERI{vc}{vc} - \ERI{vv}{cc}
R_{\BSE2}^{\sigma} = \Delta\eGF{} + 2 \sigma \ERI{vc}{vc} - \ERI{vv}{cc}
\\
B^{\stat} = 2 \sigma \ERI{vc}{vc} - \ERI{vc}{cv}
C_{\BSE2}^{\sigma} = 2 \sigma \ERI{vc}{vc} - \ERI{vc}{cv}
\end{gather}
\end{subequations}
The dynamical part of the kernel for BSE2 (that we will call dRPAx for notational consistency) is a bit ugly but it simplifies greatly in the case of the present model to yield
To avoid any confusion with the results of Sec.~\ref{sec:Maitra} and for notational consistency with Sec.~\ref{sec:BSE}, we have labeled this static Hamiltonian as BSE2.
The correlation part of the dynamical kernel for BSE2 is a bit cumbersome \cite{Zhang_2013,Rebolini_2016,Rebolini_PhD} but it simplifies greatly in the case of the present model to yield
\begin{equation}
\bH^{\dRPAx} = \bH^{\RPAx} +
\bH_{\dBSE2}^{\sigma} = \bH_{\BSE2}^{\sigma} +
\begin{pmatrix}
A(\omega) & B
f_{\dBSE2}^{\co,\sigma}(\omega) & f_{\dBSE2}^{\co,\sigma}
\\
-B & -A(-\omega)
-f_{\dBSE2}^{\co,\sigma} & -f_{\BSE2}^{\co,\sigma}(-\omega)
\end{pmatrix}
\end{equation}
with
\begin{subequations}
\begin{gather}
A^{\updw}(\omega) = - \frac{4 \ERI{cv}{vv} \ERI{vc}{cc} - \ERI{vc}{cc}^2 - \ERI{cv}{vv}^2 }{\omega - 2 \Delta\eGF{}}
f_{\dBSE2}^{\co,\updw}(\omega) = - \frac{4 \ERI{cv}{vv} \ERI{vc}{cc} - \ERI{vc}{cc}^2 - \ERI{cv}{vv}^2 }{\omega - 2 \Delta\eGF{}}
\\
B^{\updw} = - \frac{4 \ERI{vc}{cv}^2 - \ERI{cc}{cc} \ERI{vc}{cv} - \ERI{vv}{vv} \ERI{vc}{cv} }{2 \Delta\eGF{}}
f_{\dBSE2}^{\co,\updw} = - \frac{4 \ERI{vc}{cv}^2 - \ERI{cc}{cc} \ERI{vc}{cv} - \ERI{vv}{vv} \ERI{vc}{cv} }{2 \Delta\eGF{}}
\end{gather}
\end{subequations}
and
\begin{subequations}
\begin{gather}
A^{\upup}(\omega) = - \frac{ \ERI{vc}{cc}^2 + \ERI{cv}{vv}^2 }{\omega - 2 \Delta\eGF{}}
f_{\dBSE2}^{\co,\upup}(\omega) = - \frac{ \ERI{vc}{cc}^2 + \ERI{cv}{vv}^2 }{\omega - 2 \Delta\eGF{}}
\\
B^{\upup} = - \frac{\ERI{cc}{cc} \ERI{vc}{cv} + \ERI{vv}{vv} \ERI{vc}{cv} }{2 \Delta\eGF{}}
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 the coupling blocks $B$ are frequency independent, as they should.
This has an important consequence as this lack of frequency dependence removes one of the spurious pole.
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 (see Fig.~\ref{fig:dBSE2}).
Numerical results for the two-level model are reported in Table \ref{tab:RPAx} with the usual approximations and perturbative treatments.
In the case of dRPAx, the perturbative partitioning is simply
However, one spurious triplet excitation remains.
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^{\dRPAx}(\omega) = \underbrace{\bH^{\RPAx}}_{\bH^{(0)}} + \underbrace{\qty[ \bH^{\dRPAx}(\omega) - \bH^{\RPAx} ]}_{\bH^{(1)}}
\bH_{\dBSE2}^{\sigma}(\omega)
= \underbrace{\bH_{\BSE2}^{\sigma}}_{\bH_{\pBSE2}^{(0)}}
+ \underbrace{\qty[ \bH_{\dBSE2}^{\sigma}(\omega) - \bH_{\BSE2}^{\sigma} ]}_{\bH_{\pBSE2}^{(1)}}
\end{equation}
This might not be the smartest way of decomposing the Hamiltonian though but it seems to work fine.
%%% TABLE II %%%
\begin{table*}
\caption{RPAx singlet and triplet excitation energies at various levels of theory.
\label{tab:RPAx}
\caption{BSE2 singlet and triplet excitation energies (in hartree) at various levels of theory.
\label{tab:BSE2}
}
\begin{ruledtabular}
\begin{tabular}{|c|ccccccc|c|}
Singlets & RPAx & RPAx1 & RPAx1(dTDA) & dRPAx & RPAx(TDA) & RPAx1(TDA) & dRPAx(TDA) & Exact \\
Singlets & BSE2 & pBSE2 & pBSE2(dTDA) & dBSE2 & BSE2(TDA) & pBSE2(TDA) & dBSE2(TDA) & Exact \\
\hline
$\omega_1$ & 1.84903 & 1.90940 & 1.90950 & 1.90362 & 1.86299 & 1.92356 & 1.92359 & 1.92145 \\
$\omega_2$ & & & & & & & & \\
$\omega_3$ & & & & 4.47124 & & & 4.47097 & 3.47880 \\
\hline
Triplets & RPAx & RPAx1 & RPAx1(dTDA) & dRPAx & RPAx(TDA) & RPAx1(TDA) & dRPAx(TDA) & Exact \\
Triplets & BSE2 & pBSE2 & pBSE2(dTDA) & dBSE2 & BSE2(TDA) & pBSE2(TDA) & dBSE2(TDA) & Exact \\
\hline
$\omega_1$ & 1.38912 & 1.44285 & 1.44304 & 1.42564 & 1.40765 & 1.46154 & 1.46155 & 1.47085 \\
$\omega_2$ & & & & & & & & \\
@ -627,74 +650,81 @@ This might not be the smartest way of decomposing the Hamiltonian though but it
\end{table*}
%%% %%% %%% %%%
As compared to dBSE, dBSE2 produces much larger corrections to the static excitation energies probably due to the poorer quality of its static reference (CIS or TDHF).
Overall, the accuracy of dBSE and dBSE2 are comparable (see Tables \ref{tab:BSE} and \ref{tab:BSE2}) for single excitations although their behavior is quite different.
For the double excitation, dBSE2 yields a slightly better energy, yet still in quite poor agreement with the exact value.
%%% FIGURE 3 %%%
\begin{figure}
\includegraphics[width=\linewidth]{dBSE2}
\caption{
$\det[\bH^{\dRPAx}(\omega) - \omega \bI]$ as a function of $\omega$ for both the singlet (red) and triplet (blue) manifolds.
\label{fig:dBSE2}
$\det[\bH(\omega) - \omega \bI]$ as a function of $\omega$ for both the singlet (gray and black) and triplet (orange and red) manifolds.
The static BSE2 Hamiltonian (dashed) and dynamic dBSE2 Hamiltonian (solid) are considered.
\label{fig:BSE2}
}
\end{figure}
%%% %%% %%% %%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\subsection{Sangalli's kernel}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\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.
The dynamical BSE Hamiltonian with Sangalli's kernel is
\begin{equation}
\bH^\text{NC}(\omega) =
\begin{pmatrix}
H(\omega) & K(\omega)
\\
-K(-\omega) & -H(-\omega)
\end{pmatrix}
\end{equation}
with
\begin{subequations}
\begin{gather}
H_{ia,jb}(\omega) = \delta_{ij} \delta_{ab} (\eGW{a} - \eGW{i}) + \Xi_{ia,jb} (\omega)
\\
K_{ia,jb}(\omega) = \Xi_{ia,bj} (\omega)
\end{gather}
\end{subequations}
and
\begin{subequations}
\begin{gather}
\Xi_{ia,jb} (\omega) = \sum_{m \neq n} \frac{ C_{ia,mn} C_{jb,mn} }{\omega - ( \omega_{m} + \omega_{n})}
\\
C_{ia,mn} = \frac{1}{2} \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}
H(\omega) & = \Delta\eGW{} + \Xi_H (\omega)
\\
K(\omega) & = \Xi_K (\omega)
\end{align}
\begin{gather}
\Xi_H (\omega) = 2 \frac{ [\ERI{vv}{vc} + \ERI{vc}{cc}]^2 }{\omega - 2\omega_1}
\\
\Xi_K (\omega) = 0
\end{gather}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%\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} = \frac{1}{2} \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}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
What have we learnt here?
What have we learned here?
%%%%%%%%%%%%%%%%%%%%%%%%
\acknowledgements{
%%%%%%%%%%%%%%%%%%%%%%%%
The author thanks the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (Grant agreement No.~863481) for financial support.
He also thanks Xavier Blase for numerous insightful discussions on dynamical kernels.}
%%%%%%%%%%%%%%%%%%%%%%%%
% BIBLIOGRAPHY
\bibliography{../BSEdyn}

BIN
Notes/Maitra.pdf
View File

Binary file not shown.

BIN
Notes/dBSE-TDA.pdf
View File

Binary file not shown.

BIN
Notes/dBSE.pdf
View File

Binary file not shown.

BIN
Notes/dBSE2.pdf
View File

Binary file not shown.