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

minor changes in abstract and intro

Browse Source
This commit is contained in:
Pierre-Francois Loos 2021-01-20 22:25:22 +01:00
parent 59ad3ae420
commit 2c14770632
1 changed files with 6 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
Manuscript/sfBSE.tex
View File

@ -201,11 +201,11 @@
\affiliation{\LCPQ} \affiliation{\LCPQ}
\begin{abstract} \begin{abstract}
Like adiabatic time-dependent density-functional theory (TD-DFT), the Bethe-Salpeter equation (BSE) formalism of many-body perturbation theory, in its static approximation, is ``blind'' to double (and higher) excitations, which are, for example, ubiquitous in conjugated molecules like polyenes. Like adiabatic time-dependent density-functional theory (TD-DFT), the Bethe-Salpeter equation (BSE) formalism of many-body perturbation theory, in its static approximation, is ``blind'' to double (and higher) excitations, which are ubiquitous, for example, in conjugated molecules like polyenes.
Here, we apply the spin-flip \textit{ansatz} (which considers the lowest triplet state as the reference configuration instead of the singlet ground state) to the BSE formalism in order to access double excitations. Here, we apply the spin-flip \textit{ansatz} (which considers the lowest triplet state as the reference configuration instead of the singlet ground state) to the BSE formalism in order to access, in particular, double excitations.
The present scheme is based on a spin-unrestricted version of the $GW$ approximation employed to compute the charged excitations and screened Coulomb potential required for the BSE calculations. The present scheme is based on a spin-unrestricted version of the $GW$ approximation employed to compute the charged excitations and screened Coulomb potential required for the BSE calculations.
Dynamical corrections to the static BSE optical excitations are taken into account via an unrestricted generalization of our recently developed (renormalized) perturbative treatment. Dynamical corrections to the static BSE optical excitations are taken into account via an unrestricted generalization of our recently developed (renormalized) perturbative treatment.
The performance of the present spin-flip BSE formalism is illustrated by computing the vertical excitation energies of the beryllium atom, the hydrogen molecule at various bond lengths, and cyclobutadiene in its rectangular and square-planar geometries. The performance of the present spin-flip BSE formalism is illustrated by computing excited-state energies of the beryllium atom, the hydrogen molecule at various bond lengths, and cyclobutadiene in its rectangular and square-planar geometries.
%\bigskip %\bigskip
%\begin{center} %\begin{center}
% \boxed{\includegraphics[width=0.5\linewidth]{TOC}} % \boxed{\includegraphics[width=0.5\linewidth]{TOC}}
@ -225,17 +225,17 @@ Accurately predicting ground- and excited-state energies (hence excitation energ
An armada of theoretical and computational methods have been developed to this end, each of them being plagued by its own flaws. \cite{Roos_1996,Piecuch_2002,Dreuw_2005,Krylov_2006,Sneskov_2012,Gonzales_2012,Laurent_2013,Adamo_2013,Ghosh_2018,Blase_2020,Loos_2020d,Casanova_2020} An armada of theoretical and computational methods have been developed to this end, each of them being plagued by its own flaws. \cite{Roos_1996,Piecuch_2002,Dreuw_2005,Krylov_2006,Sneskov_2012,Gonzales_2012,Laurent_2013,Adamo_2013,Ghosh_2018,Blase_2020,Loos_2020d,Casanova_2020}
The fact that none of these methods is successful in every chemical scenario has encouraged chemists to carry on the development of new excited-state methodologies, their main goal being to get the most accurate excitation energies (and properties) at the lowest possible computational cost in the most general context. \cite{Loos_2020d} The fact that none of these methods is successful in every chemical scenario has encouraged chemists to carry on the development of new excited-state methodologies, their main goal being to get the most accurate excitation energies (and properties) at the lowest possible computational cost in the most general context. \cite{Loos_2020d}
Originally developed in the framework of nuclear physics, \cite{Salpeter_1951} and popularized in condensed-matter physics, \cite{Sham_1966,Strinati_1984,Delerue_2000} one of the new emerging method in the computational chemistry landscape is the Bethe-Salpeter equation (BSE) formalism \cite{Salpeter_1951,Strinati_1988,Albrecht_1998,Rohlfing_1998,Benedict_1998,vanderHorst_1999,Blase_2018,Blase_2020} from many-body perturbation theory \cite{Onida_2002,Martin_2016} which, based on an underlying $GW$ calculation to compute accurate charged excitations and the dynamically-screened Coulomb potential, \cite{Hedin_1965,Golze_2019} is able to provide accurate optical (\ie, neutral) excitations for molecular systems at a rather modest computational cost.\cite{Rohlfing_1999a,Horst_1999,Puschnig_2002,Tiago_2003,Boulanger_2014,Jacquemin_2015a,Bruneval_2015,Jacquemin_2015b,Hirose_2015,Jacquemin_2017a,Jacquemin_2017b,Rangel_2017,Krause_2017,Gui_2018,Blase_2018,Liu_2020,Blase_2020,Holzer_2018a,Holzer_2018b,Loos_2020e} Originally developed in the framework of nuclear physics, \cite{Salpeter_1951} and popularized in condensed-matter physics, \cite{Sham_1966,Strinati_1984,Delerue_2000} one of the new emerging method in the computational chemistry landscape is the Bethe-Salpeter equation (BSE) formalism \cite{Salpeter_1951,Strinati_1988,Albrecht_1998,Rohlfing_1998,Benedict_1998,vanderHorst_1999,Blase_2018,Blase_2020} from many-body perturbation theory \cite{Onida_2002,Martin_2016} which, based on an underlying $GW$ calculation to compute accurate charged excitations (\ie, ionization potentials and electron affinities) and the dynamically-screened Coulomb potential, \cite{Hedin_1965,Golze_2019} is able to provide accurate optical (\ie, neutral) excitations for molecular systems at a rather modest computational cost.\cite{Rohlfing_1999a,Horst_1999,Puschnig_2002,Tiago_2003,Boulanger_2014,Jacquemin_2015a,Bruneval_2015,Jacquemin_2015b,Hirose_2015,Jacquemin_2017a,Jacquemin_2017b,Rangel_2017,Krause_2017,Gui_2018,Blase_2018,Liu_2020,Blase_2020,Holzer_2018a,Holzer_2018b,Loos_2020e}
Most of BSE implementations rely on the so-called static approximation, \cite{Blase_2018,Bruneval_2016,Krause_2017,Liu_2020} which approximates the dynamical (\ie, frequency-dependent) BSE kernel by its static limit. Most of BSE implementations rely on the so-called static approximation, \cite{Blase_2018,Bruneval_2016,Krause_2017,Liu_2020} which approximates the dynamical (\ie, frequency-dependent) BSE kernel by its static limit.
Like adiabatic time-dependent density-functional theory (TD-DFT), \cite{Runge_1984,Casida_1995,Petersilka_1996,UlrichBook} the static BSE formalism is plagued by the lack of double (and higher) excitations, which are, for example, ubiquitous in conjugated molecules like polyenes \cite{Maitra_2004,Cave_2004,Saha_2006,Watson_2012,Shu_2017,Barca_2018a,Barca_2018b,Loos_2019} or the ground state of open-shell molecules. \cite{Casida_2005,Huix-Rotllant_2011,Loos_2020f} Like adiabatic time-dependent density-functional theory (TD-DFT), \cite{Runge_1984,Casida_1995,Petersilka_1996,UlrichBook} the static BSE formalism is plagued by the lack of double (and higher) excitations, which are, for example, ubiquitous in conjugated molecules like polyenes \cite{Maitra_2004,Cave_2004,Saha_2006,Watson_2012,Shu_2017,Barca_2018a,Barca_2018b,Loos_2019} or the ground state of open-shell molecules. \cite{Casida_2005,Huix-Rotllant_2011,Loos_2020f}
Indeed, both adiabatic TD-DFT \cite{Levine_2006,Tozer_2000,Elliott_2011,Maitra_2012,Maitra_2016} and static BSE \cite{ReiningBook,Romaniello_2009b,Sangalli_2011,Loos_2020h,Authier_2020} can only access (singlet and triplet) single excitations with respect to the reference determinant usually taken as the closed-shell singlet ground state. Indeed, both adiabatic TD-DFT \cite{Levine_2006,Tozer_2000,Elliott_2011,Maitra_2012,Maitra_2016} and static BSE \cite{ReiningBook,Romaniello_2009b,Sangalli_2011,Loos_2020h,Authier_2020} can only access (singlet and triplet) single excitations with respect to the reference determinant usually taken as the closed-shell singlet ground state.
Double excitations are even challenging for state-of-the-art methods, \cite{Loos_2018a,Loos_2019,Loos_2020c,Loos_2020d,Veril_2020} like the approximate third-order coupled-cluster (CC3) method \cite{Christiansen_1995b,Koch_1997} or equation-of-motion coupled-cluster with singles, doubles and triples (EOM-CCSDT). \cite{Kucharski_1991,Kallay_2004,Hirata_2000,Hirata_2004} Double excitations are even challenging for state-of-the-art methods, \cite{Loos_2018a,Loos_2019,Loos_2020c,Loos_2020d,Veril_2020} like the approximate third-order coupled-cluster (CC3) method \cite{Christiansen_1995b,Koch_1997} or equation-of-motion coupled-cluster with singles, doubles and triples (EOM-CCSDT). \cite{Kucharski_1991,Kallay_2004,Hirata_2000,Hirata_2004}
One way to access double excitations is via the spin-flip formalism established by Krylov in 2001, \cite{Krylov_2001a,Krylov_2001b,Krylov_2002} with earlier attempts by Bethe, \cite{Bethe_1931} as well as Shibuya and McKoy. \cite{Shibuya_1970} One way to access double excitations is via the spin-flip formalism established by Krylov in 2001, \cite{Krylov_2001a,Krylov_2001b,Krylov_2002} with earlier attempts by Bethe, \cite{Bethe_1931} as well as Shibuya and McKoy. \cite{Shibuya_1970}
The idea behind the spin-flip \textit{ansatz} is rather simple: instead of considering the singlet ground state as reference, the reference is taken as the lowest triplet state. The idea behind the spin-flip \textit{ansatz} is rather simple: instead of considering the singlet ground state as reference, the reference configuration is taken as the lowest triplet state.
In such a way, one can access the singlet ground state and the singlet doubly-excited state via a spin-flip deexcitation and excitation (respectively), the difference of these two excitation energies providing an estimate of the double excitation. In such a way, one can access the singlet ground state and the singlet doubly-excited state via a spin-flip deexcitation and excitation (respectively), the difference of these two excitation energies providing an estimate of the double excitation.
We refer the interested reader to Refs.~\onlinecite{Krylov_2006,Krylov_2008,Casanova_2020} for detailed reviews on spin-flip methods. We refer the interested reader to Refs.~\onlinecite{Krylov_2006,Krylov_2008,Casanova_2020} for detailed reviews on spin-flip methods.
Note that a similar idea has been exploited by the group of Weito Yang to access double excitations in the context of the particle-particle random-phase approximation. \cite{Peng_2013,Yang_2013b,Yang_2014a,Peng_2014,Zhang_2016,Sutton_2018} Note that a similar idea has been exploited by the group of Yang to access double excitations in the context of the particle-particle random-phase approximation. \cite{Peng_2013,Yang_2013b,Yang_2014a,Peng_2014,Zhang_2016,Sutton_2018}
One obvious issue of spin-flip methods is that not all double excitations are accessible in such a way. One obvious issue of spin-flip methods is that not all double excitations are accessible in such a way.
Moreover, spin-flip methods are usually hampered by spin contamination \cite{Casanova_2020} (\ie, artificial mixing with configurations of different spin multiplicities) due to spin incompleteness of the configuration interaction expansion as well as the possible spin contamination of the reference configuration. \cite{Krylov_2000b} Moreover, spin-flip methods are usually hampered by spin contamination \cite{Casanova_2020} (\ie, artificial mixing with configurations of different spin multiplicities) due to spin incompleteness of the configuration interaction expansion as well as the possible spin contamination of the reference configuration. \cite{Krylov_2000b}