sfBSE/sfBSE.tex

\documentclass[aip,jcp,reprint,noshowkeys,superscriptaddress]{revtex4-1}
\usepackage{graphicx,dcolumn,bm,xcolor,microtype,multirow,amscd,amsmath,amssymb,amsfonts,physics,longtable,wrapfig,txfonts}
\usepackage[version=4]{mhchem}

\usepackage[utf8]{inputenc}
\usepackage[T1]{fontenc}
\usepackage{txfonts}

\usepackage[
	colorlinks=true,
    citecolor=blue,
    breaklinks=true
	]{hyperref}
\urlstyle{same}

\begin{document}	

\title{Spin-Conserved and Spin-Flip Optical Excitations From the Bethe-Salpeter Equation Formalism}

\author{Enzo \surname{Monino}}
	\affiliation{\LCPQ}	 
\author{Pierre-Fran\c{c}ois \surname{Loos}}
	\email{loos@irsamc.ups-tlse.fr}
	\affiliation{\LCPQ}

\begin{abstract}
\alert{Here comes the abstract.}
%\bigskip
%\begin{center}
%	\boxed{\includegraphics[width=0.5\linewidth]{TOC}}
%\end{center}
%\bigskip
\end{abstract}

\maketitle

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Introduction}
\label{sec:intro}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\alert{Here comes the introduction.}
Unless otherwise stated, atomic units are used, and we assume real quantities throughout this manuscript.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Unrestricted $GW$ formalism}
\label{sec:UGW}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Let us consider an electronic system consisting of $n = n_\up + n_\dw$ electrons (where $n_\up$ and $n_\dw$ are the number of spin-up and spin-down electrons respectively) and $N$ one-electron basis functions.
The number of spin-up and spin-down occupied orbitals are $O_\up = n_\up$ and $O_\dw = n_\dw$, respectively, and there is $V_\up = N - O_\up$ and $V_\dw = N - O_\dw$ spin-up and spin-down virtual (\ie, unoccupied) orbitals.
The number of spin-conserved single excitations is then $S^\spc = S_{\up\up}^\spc + S_{\dw\dw}^\spc = O_\up V_\up + O_\dw V_\dw$, while the number of spin-flip excitations is $S^\spf = S_{\up\dw}^\spf + S_{\dw\up}^\spf = O_\up V_\dw + O_\dw V_\up$.
Let us denote as $\MO{p_\sig}$ the $p$th orbital of spin $\sig$ (where $\sig =$ $\up$ or $\dw$).
In the following, $i$ and $j$ are occupied orbitals, $a$ and $b$ are unoccupied orbitals, $p$, $q$, $r$, and $s$ indicate arbitrary orbitals, and $m$ labels single excitations.
It is important to understand that, in a spin-conserved excitation the hole orbital $\MO{i_\sig}$ and particle orbital $\MO{a_\sig}$ have the same spin $\sig$.
A bra and ket composed by these two orbitals will be denoted as $\rbra{ia\sig}$ and $\rket{ia\sig}$.
In a spin-flip excitation, the hole has a spin $\sig$ and the particle has the opposite spin $\bsig$.

%================================
\subsection{The dynamical screening}
%================================

Within the $GW$ formalism, the dynamical screening $W(\omega)$ is computed at the RPA level using the spin-conserved neutral excitations.
The matrix elements of $W(\omega)$ read
\begin{multline}
	W_{p_\sig q_\sig,r_\sigp s_\sigp}(\omega) = \ERI{p_\sig q_\sig}{r_\sigp s_\sigp} 
	+ \sum_m \ERI{p_\sig q_\sig}{m}\ERI{r_\sigp s_\sigp}{m} 
	\\
	\times \qty[ \frac{1}{\omega - \Om{m}{\spc,\RPA} + i \eta} - \frac{1}{\omega + \Om{m}{\spc,\RPA} - i \eta} ]
\end{multline}
where the two-electron integrals are
\begin{equation}
		\ERI{p_\sig q_\tau}{r_\sigp s_\taup} = \iint \MO{p_\sig}(\br) \MO{q_\tau}(\br) \frac{1}{\abs{\br - \br'}}  \MO{r_\sigp}(\br') \MO{s_\taup}(\br') d\br d\br'
\end{equation}

\begin{equation}
		\ERI{p_\sig q_\sig}{m} = \sum_{ia\sigp} \ERI{p_\sig q_\sig}{r_\sigp s_\sigp} (\bX{m}{\spc,\RPA}+\bY{m}{\spc,\RPA})_{i_\sigp a_\sigp}
\end{equation}

\begin{equation}
\label{eq:LR-RPA}
	\begin{pmatrix}
		\bA{}{\spc,\RPA}	&	\bB{}{\spc,\RPA}	\\
		-\bB{}{\spc,\RPA}	&	-\bA{}{\spc,\RPA}	\\
	\end{pmatrix}
	\cdot
	\begin{pmatrix}
		\bX{m}{\spc,\RPA}	\\
		\bY{m}{\spc,\RPA}	\\
	\end{pmatrix}
	=
	\Om{m}{\spc,\RPA}
	\begin{pmatrix}
		\bX{m}{\spc,\RPA}	\\
		\bY{m}{\spc,\RPA}	\\
	\end{pmatrix},
\end{equation}

The spin structure of these matrices are general and reads
\begin{align}
\label{eq:LR-RPA-AB}
	\bA{}{\spc} & = \begin{pmatrix}
		\bA{\upup,\upup}{}	&	\bA{\upup,\dwdw}{}	\\
		\bA{\dwdw,\upup}{}	&	\bA{\dwdw,\dwdw}{}	\\
	\end{pmatrix}
	&
	\bB{}{\spc} & = \begin{pmatrix}
		\bB{\upup,\upup}{}	&	\bB{\upup,\dwdw}{}	\\
		\bB{\dwdw,\upup}{}	&	\bB{\dwdw,\dwdw}{}	\\
	\end{pmatrix}
\\
\label{eq:LR-RPA-AB}
	\bA{}{\spf} & = \begin{pmatrix}
		\bA{\updw,\updw}{}	&	\bO	\\
		\bO					&	\bA{\dwup,\dwup}{}	\\
	\end{pmatrix}
	&
	\bB{}{\spf} & = \begin{pmatrix}
		\bO			&	\bB{\updw,\dwup}{}	\\
		\bB{\dwup,\updw}{}	&	\bO	\\
	\end{pmatrix}
\end{align}
with
\begin{subequations}
\begin{align}
	\label{eq:LR_RPA-A}
	\A{i_\sig a_\tau,j_\sigp b_\taup}{\RPA} & = \delta_{ij} \delta_{ab} \delta_{\sig \sigp}  \delta_{\tau \taup} (\e{a_\tau} - \e{i_\sig}) + \ERI{i_\sig a_\tau}{b_\sigp j_\taup}
	\\ 
	\label{eq:LR_RPA-B}
	\B{i_\sig a_\tau,j_\sigp b_\taup}{\RPA} & = \ERI{i_\sig a_\tau}{j_\sigp b_\taup}
\end{align}
\end{subequations}
from which we obtain, at the RPA level, the following expressions
\begin{subequations}
\begin{align}
	\label{eq:LR_RPA-Asc}
	\A{i_\sig a_\sig,j_\sigp b_\sigp}{\spc,\RPA} & = \delta_{ij} \delta_{ab} \delta_{\sig \sigp} (\e{a_\sig} - \e{i_\sig}) + \ERI{i_\sig a_\sig}{b_\sigp j_\sigp}
	\\ 
	\label{eq:LR_RPA-Bsc}
	\B{i_\sig a_\sig,j_\sigp b_\sigp}{\spc,\RPA} & = \ERI{i_\sig a_\sig}{j_\sigp b_\sigp}
\end{align}
\end{subequations}
for the spin-conserved excitations and 
\begin{subequations}
\begin{align}
	\label{eq:LR_RPA-Asf}
	\A{i_\sig a_\bsig,j_\sig b_\bsig}{\spf,\RPA} & = \delta_{ij} \delta_{ab} (\e{a_\bsig} - \e{i_\sig})
	\\ 
	\label{eq:LR_RPA-Bsf}
	\B{i_\sig a_\bsig,j_\bsig b_\sig}{\spf,\RPA} & = 0
\end{align}
\end{subequations}
for the spin-flip excitations.

%================================
\subsection{The $GW$ self-energy}
%================================

\begin{equation}
\begin{split}
	\SigC{p_\sig q_\sig}(\omega) 
	& = \sum_i \sum_m \frac{\ERI{p_\sig i_\sig}{m} \ERI{q_\sig i_\sig}{m}}{\omega - \e{i_\sig} + \Om{m}{\spc,\RPA} - i \eta}
	\\
	& + \sum_a \sum_m \frac{\ERI{p_\sig a_\sig}{m} \ERI{q_\sig a_\sig}{m}}{\omega - \e{a_\sig} - \Om{m}{\spc,\RPA} + i \eta}			 
\end{split}
\end{equation}

The quasiparticle energies $\eGW{p}$ are obtained by solving the frequency-dependent quasiparticle equation
\begin{equation}
	\omega = \eHF{p\sigma} + \SigC{p\sigma}(\omega)
\end{equation}

%================================
\subsection{The Bethe-Salpeter equation formalism}
%================================
Defining $W^{\stat}_{p_\sig q_\sig,r_\sigp s_\sigp} = W_{p_\sig q_\sig,r_\sigp s_\sigp}(\omega = 0)$, we have
\begin{subequations}
\begin{align}
	\label{eq:LR_BSE-A}
	\A{i_\sig a_\tau,j_\sigp b_\taup}{\BSE} & = \A{i_\sig a_\tau,j_\sigp b_\taup}{\RPA} - \delta_{\sig \sigp} W^{\stat}_{i_\sig j_\sigp,b_\taup a_\tau}
	\\ 
	\label{eq:LR_BSE-B}
	\B{i_\sig a_\tau,j_\sigp b_\taup}{\BSE} & = \B{i_\sig a_\tau,j_\sigp b_\taup}{\RPA} - \delta_{\sig \sigp} W^{\stat}_{i_\sig b_\taup,j_\sigp a_\tau}
\end{align}
\end{subequations}
from which we obtain, at the BSE level, the following expressions
\begin{subequations}
\begin{align}
	\label{eq:LR_BSE-Asc}
	\A{i_\sig a_\sig,j_\sigp b_\sigp}{\spc,\BSE} & = \A{i_\sig a_\sig,j_\sigp b_\sigp}{\spc,\RPA} - \delta_{\sig \sigp} W^{\stat}_{i_\sig j_\sigp,b_\sigp a_\sig}
	\\ 
	\label{eq:LR_BSE-Bsc}
	\B{i_\sig a_\sig,j_\sigp b_\sigp}{\spc,\BSE} & = \B{i_\sig a_\sig,j_\sigp b_\sigp}{\spc,\RPA} - \delta_{\sig \sigp} W^{\stat}_{i_\sig b_\sigp,j_\sigp a_\sig}
\end{align}
\end{subequations}
for the spin-conserved excitations and 
\begin{subequations}
\begin{align}
	\label{eq:LR_BSE-Asf}
	\A{i_\sig a_\bsig,j_\sig b_\bsig}{\spf,\BSE} & = \A{i_\sig a_\bsig,j_\sig b_\bsig}{\spf,\RPA} - W^{\stat}_{i_\sig j_\sig,b_\bsig a_\bsig}
	\\ 
	\label{eq:LR_BSE-Bsf}
	\B{i_\sig a_\bsig,j_\bsig b_\sig}{\spf,\BSE} & = - W^{\stat}_{i_\sig b_\sig,j_\bsig a_\bsig}
\end{align}
\end{subequations}
for the spin-flip excitations.

%================================
\subsection{The dynamical Bethe-Salpeter equation correction}

%================================
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Computational details}
\label{sec:compdet}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Conclusion}
\label{sec:ccl}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


%%%%%%%%%%%%%%%%%%%%%%%%
\acknowledgements{
We would like to thank Xavier Blase and Denis Jacquemin for insightful discussions. 
This project has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (Grant agreement No.~863481).}
%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section*{Data availability statement}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
The data that supports the findings of this study are available within the article and its supplementary material.

%%%%%%%%%%%%%%%%%%%%%%%%
\bibliography{sf-BSE}
%%%%%%%%%%%%%%%%%%%%%%%%

\end{document}