diff --git a/BSEdyn.tex b/BSEdyn.tex index 57d2862..aa0c85d 100644 --- a/BSEdyn.tex +++ b/BSEdyn.tex @@ -200,9 +200,9 @@ \affiliation{\NEEL} \begin{abstract} -Similarly to the ubiquitous adiabatic approximation in time-dependent density-functional theory, the static approximation, which substitutes a dynamical (\ie, frequency-dependent) kernel by its static limit, is usually enforced in most implementations of the Bethe-Salpeter equation (BSE) formalism. -Here, going beyond the static approximation, we compute molecular excitation energies thanks to a renormalized first-order perturbative correction to the static BSE excitation energies. -The present correction goes beyond the plasmon-pole approximation as the dynamical screening of the Coulomb interaction is computed exactly. +Similar to the ubiquitous adiabatic approximation in time-dependent density-functional theory, the static approximation, which substitutes a dynamical (\ie, frequency-dependent) kernel by its static limit, is usually enforced in most implementations of the Bethe-Salpeter equation (BSE) formalism. +Here, going beyond the static approximation, we compute the dynamical correction in the electron-hole screening for molecular excitation energies thanks to a renormalized first-order perturbative correction to the static BSE excitation energies. +The present dynamical correction goes beyond the plasmon-pole approximation as the dynamical screening of the Coulomb interaction is computed exactly. Moreover, we investigate quantitatively the effect of the Tamm-Dancoff approximation by computing both the resonant and anti-resonant dynamical corrections to the BSE excitation energies. %\\ %\bigskip @@ -219,7 +219,7 @@ Moreover, we investigate quantitatively the effect of the Tamm-Dancoff approxima \label{sec:intro} %%%%%%%%%%%%%%%%%%%%%%%% The Bethe-Salpeter equation (BSE) formalism \cite{Salpeter_1951,Strinati_1988} is to the $GW$ approximation \cite{Hedin_1965,Golze_2019} of many-body perturbation theory (MBPT) \cite{Martin_2016} what time-dependent density-functional theory (TD-DFT) \cite{Runge_1984,Casida_1995} is to Kohn-Sham density-functional theory (KS-DFT), \cite{Hohenberg_1964,Kohn_1965} an affordable way of computing the neutral excitations of a given electronic system. -In recent years, it has been shown to be a valuable tool for computational theoretical chemists with a large number of systematic benchmark studies on large molecular systems appearing in the scientific literature \cite{Korbel_2014,Jacquemin_2015a,Bruneval_2015,Jacquemin_2015b,Hirose_2015,Jacquemin_2017a,Jacquemin_2017b,Krause_2017,Gui_2018} (see Ref.~\onlinecite{Blase_2018} for a recent review). +In recent years, it has been shown to be a valuable tool for computational chemists with a large number of systematic benchmark studies on large molecular systems appearing in the literature \cite{Korbel_2014,Jacquemin_2015a,Bruneval_2015,Jacquemin_2015b,Hirose_2015,Jacquemin_2017a,Jacquemin_2017b,Krause_2017,Gui_2018} (see Ref.~\onlinecite{Blase_2018} for a recent review). Taking the optical gap (\ie, the lowest optical excitation energy) as an example, BSE builds on top of a $GW$ calculation by adding up excitonic effects $\EB$ to the $GW$ HOMO-LUMO gap \begin{equation} @@ -249,14 +249,14 @@ They are, moreover, a real challenge for high-level computational methods. \cite Going beyond the static approximation is tricky and very few groups have dared to take the plunge. \cite{Strinati_1988,Rohlfing_2000,Sottile_2003,Ma_2009a,Ma_2009b,Romaniello_2009b,Sangalli_2011,Huix-Rotllant_2011,Zhang_2013,Rebolini_2016,Olevano_2019,Lettmann_2019} Nonetheless, it is worth mentioning the seminal work of Strinati, \cite{Strinati_1988} who \titou{bla bla bla.} Following Strinati's footsteps, Rohlfing and coworkers have developed an efficient way of taking into account, thanks to first-order perturbation theory, the dynamical effects via a plasmon-pole approximation combined with the Tamm-Dancoff approximation (TDA). \cite{Rohlfing_2000,Ma_2009a,Ma_2009b,Baumeier_2012b} -With such as scheme, they have been able to compute the excited states of biological chromophores, showing that taking into account the electron-hole dynamical screening is important for an accurate description of the lowest $n \ra \pi^*$ excitations. \cite{Ma_2009a,Ma_2009b,Baumeier_2012b} +With such a scheme, they have been able to compute the excited states of biological chromophores, showing that taking into account the electron-hole dynamical screening is important for an accurate description of the lowest $n \ra \pi^*$ excitations. \cite{Ma_2009a,Ma_2009b,Baumeier_2012b} Indeed, studying PYP, retinal and GFP chromophore models, Ma \textit{et al.}~found that \textit{``the influence of dynamical screening on the excitation energies is about $0.1$ eV for the lowest $\pi \ra \pis$ transitions, but for the lowest $n \ra \pis$ transitions the influence is larger, up to $0.25$ eV.''} \cite{Ma_2009b} A similar conclusion was reached in Ref.~\onlinecite{Ma_2009a}. Zhang \textit{et al.}~have studied the frequency-dependent second-order Bethe-Salpeter kernel and they have observed an appreciable improvement over configuration interaction with singles (CIS), time-dependent Hartree-Fock (TDHF), and adiabatic TD-DFT results. \cite{Zhang_2013} Rebolini and Toulouse have performed a similar investigation in a range-separated context, and they have reported a modest improvement over its static counterpart. \cite{Rebolini_2016,Rebolini_PhD} In these two latter studies, they also followed a (non-self-consistent) perturbative approach within the TDA with a renormalization of the first-order perturbative correction. -It is important to note that, although all the studies mentioned above are clearly going beyond the static approximation of BSE, they are not able to recover double excitations as the perturbative treatment makes ultimately the BSE kernel static. +It is important to note that, although all the studies mentioned above are clearly going beyond the static approximation of BSE, they are not able to recover additional excitations as the perturbative treatment makes ultimately the BSE kernel static. However, it does permit to recover, for transitions with a dominant single-excitation character, additional relaxation effects coming from higher excitations (and, in particular, non-interacting double excitations). These higher excitations would be explicitly present in the BSE Hamiltonian by ``unfolding'' the dynamical BSE kernel, and one would recover a linear eigenvalue problem with, nonetheless, a much larger dimension. Finally, let us also mentioned the work of Romaniello and coworkers, \cite{Romaniello_2009b,Sangalli_2011} in which the authors genuinely accessed additional excitations by solving the non-linear, frequency-dependent eigenvalue problem. @@ -282,7 +282,7 @@ We present, in a second step, the perturbative implementation of the dynamical c The two-particle correlation function $L(1,2; 1',2')$ --- a central quantity in the BSE formalism --- relates the variation of the one-body Green's function $G(1,1')$ with respect to an external non-local perturbation $U(2',2)$, \ie, \begin{equation} - iL(1,2; 1',2') = \pdv{G(1,1')}{U(2',2)} + iL(1,2; 1',2') = \pdv{G(1,1')}{U(2',2)}, \end{equation} where, \eg, $1 \equiv (\bx_1 t_1)$ is a space-spin plus time composite variable. The relation between $G$ and the one-body charge density $\rho(1) = -i G(1,1^+)$ provides a direct connection with the density-density susceptibility $\chi(1,2) = L(1,2;1^+,2^+)$ at the core of TD-DFT. @@ -349,7 +349,7 @@ where $\tau_{12} = t_1 - t_2$, $\theta$ is the Heaviside step function, and \end{subequations} The $\Om{s}{}$'s are the neutral excitation energies of interest. -Picking up the $e^{+i \Om{s}{} t_2 }$ component in $L(1,2; 1',2')$ and $L(6,2;5,2')$, simplifying further by $\tchi_s(\bx_2,\bx_{2'})$ on both side of the BSE [see Eq.~\eqref{eq:BSE}], we seek the $e^{-i \Om{s}{} t_1 }$ Fourier component associated with the right-hand side of a modified dynamical BSE, which reads +Picking up the $e^{+i \Om{s}{} t_2 }$ component in $L(1,2; 1',2')$ and $L(6,2;5,2')$, simplifying further by $\tchi_s(\bx_2,\bx_{2'})$ on both side of the BSE [see Eq.~\eqref{eq:BSE}], we are left with seeking the $e^{-i \Om{s}{} t_1 }$ Fourier component associated with the right-hand side of a modified dynamical BSE, which reads \begin{multline} \label{eq:BSE_2} \mel{N}{T [ \hpsi(\bx_1) \hpsi^{\dagger}(\bx_{1}') ] } {N,s} e^{ - i \Om{s}{} t_1 } \theta ( \tau_{12} ) @@ -413,7 +413,7 @@ leads to the following simplified BSE kernel \Xi(3,5;4,6) = v(3,6) \delta(3,4) \delta(5,6) - W(3^+,4) \delta(3,6) \delta(4,5), \end{equation} where $W$ is the dynamically-screened Coulomb operator. -The $GW$ quasiparticle energies $\eGW{p}$ are good approximations to the removal/addition energies $\e{p}$ introduced in Eq.~\eqref{eq:G-Lehman}. +The $GW$ quasiparticle energies $\eGW{p}$ are usually good approximations to the removal/addition energies $\e{p}$ introduced in Eq.~\eqref{eq:G-Lehman}. %Selecting $(p,q)=(j,b)$ yields the largest components %$X_{jb}^{s} = \mel{N}{\ha_j^{\dagger} \ha_b}{N,s}$, while $(p,q)=(b,j)$ yields much weaker %$Y_{jb}^{s} = \mel{N}{\ha_b^{\dagger} \ha_j}{N,s}$ contributions. @@ -455,7 +455,7 @@ is an effective dynamically-screened Coulomb potential, \cite{Romaniello_2009b} \subsection{Dynamical screening} %================================= -In the present study, we consider the exact spectral representation of $W(\omega)$ at the random-phase approximation (RPA) level, which reads +In the present study, we consider the exact spectral representation of $W$ at the random-phase approximation (RPA) level: \begin{multline} \label{eq:W} W_{ij,ab}(\omega) @@ -498,8 +498,8 @@ with \B{ia,jb}{\RPA} & = 2 \ERI{ia}{bj}, \end{align} \end{subequations} -where the $\e{p}$'s are taken as the HF orbital energies in the case of $G_0W_0$ or as the $GW$ quasiparticle energies in the case of self-consistent schemes such as ev$GW$. -The RPA matrices $\bA{\RPA}$ and $\bB{\RPA}$ in Eq.~\eqref{eq:LR-RPA} are of size $\Nocc \Nvir \times \Nocc \Nvir$, where $\Nocc$ and $\Nvir$ are the number of occupied and virtual orbitals (\ie, $\Norb = \Nocc + \Nvir$ is the total number of spatial orbitals), respectively, and $\bX{m}{}$, and $\bY{m}{}$ are (eigen)vectors of length $\Nocc \Nvir$. +where the $\e{p}$'s are taken as the HF orbital energies in the case of $G_0W_0$ \cite{Hybertsen_1985a, Hybertsen_1986} or as the $GW$ quasiparticle energies in the case of self-consistent schemes such as ev$GW$. \cite{Hybertsen_1986,Shishkin_2007,Blase_2011,Faber_2011} +The RPA matrices $\bA{\RPA}$ and $\bB{\RPA}$ in Eq.~\eqref{eq:LR-RPA} are of size $\Nocc \Nvir \times \Nocc \Nvir$, where $\Nocc$ and $\Nvir$ are the number of occupied and virtual orbitals (\ie, $\Norb = \Nocc + \Nvir$ is the total number of spatial orbitals), respectively, and $\bX{m}{\RPA}$, and $\bY{m}{\RPA}$ are (eigen)vectors of length $\Nocc \Nvir$. The analysis of the poles of the integrand in Eq.~\eqref{eq:wtilde} yields \begin{multline} @@ -511,7 +511,7 @@ The analysis of the poles of the integrand in Eq.~\eqref{eq:wtilde} yields Due to excitonic effects, the lowest BSE excitation energy, $\Om{1}{}$, stands lower than the lowest RPA excitation energy, $\Om{1}{\RPA}$, so that, $\Om{ib}{s} - \Om{m}{\RPA} < 0 $ and $\widetilde{W}_{ij,ab}(\Om{s}{})$ has no resonances. Furthermore, $\Om{ib}{s}$ and $\Om{ja}{s}$ are necessarily negative quantities for in-gap low-lying BSE excitations. Thus, we have $\abs*{\Omega_{ib}^{s} - \Om{m}{\RPA}} > \Omega_m^{\RPA}$. -As a consequence, we observe a reduction of the electron-hole screening, \ie, an enhancement of electron-hole stabilizing binding energy, as compared to the standard static BSE, and yields smaller (red-shifted) excitation energies. +As a consequence, we observe a reduction of the electron-hole screening, \ie, an enhancement of electron-hole binding energy, as compared to the standard static BSE, and yields smaller (red-shifted) excitation energies. %================================ \subsection{Perturbative dynamical correction}