starting results and adding stuff in theory
This commit is contained in:
parent
df184cd10e
commit
dd948a8a43
182
BSE-PES.tex
182
BSE-PES.tex
@ -41,9 +41,8 @@
|
|||||||
\newcommand{\dbr}[1]{d\br{#1}}
|
\newcommand{\dbr}[1]{d\br{#1}}
|
||||||
|
|
||||||
% methods
|
% methods
|
||||||
\newcommand{\GW}{$GW$}
|
\newcommand{\evGW}{ev$GW$}
|
||||||
\newcommand{\evGW}{ev\GW}
|
\newcommand{\qsGW}{qs$GW$}
|
||||||
\newcommand{\qsGW}{qs\GW}
|
|
||||||
\newcommand{\GOWO}{$G_0W_0$}
|
\newcommand{\GOWO}{$G_0W_0$}
|
||||||
\newcommand{\Hxc}{\text{Hxc}}
|
\newcommand{\Hxc}{\text{Hxc}}
|
||||||
\newcommand{\xc}{\text{xc}}
|
\newcommand{\xc}{\text{xc}}
|
||||||
@ -64,10 +63,9 @@
|
|||||||
\newcommand{\Ec}{E_\text{c}}
|
\newcommand{\Ec}{E_\text{c}}
|
||||||
\newcommand{\EHF}{E^\text{HF}}
|
\newcommand{\EHF}{E^\text{HF}}
|
||||||
\newcommand{\EBSE}{E^\text{BSE}}
|
\newcommand{\EBSE}{E^\text{BSE}}
|
||||||
\newcommand{\EcRPA}{E_\text{c}^\text{dRPA}}
|
\newcommand{\EcRPA}{E_\text{c}^\text{RPA}}
|
||||||
|
\newcommand{\EcRPAx}{E_\text{c}^\text{RPAx}}
|
||||||
\newcommand{\EcBSE}{E_\text{c}^\text{BSE}}
|
\newcommand{\EcBSE}{E_\text{c}^\text{BSE}}
|
||||||
\newcommand{\EcsBSE}{{}^1\EcBSE}
|
|
||||||
\newcommand{\EctBSE}{{}^3\EcBSE}
|
|
||||||
\newcommand{\IP}{\text{IP}}
|
\newcommand{\IP}{\text{IP}}
|
||||||
\newcommand{\EA}{\text{EA}}
|
\newcommand{\EA}{\text{EA}}
|
||||||
|
|
||||||
@ -77,7 +75,7 @@
|
|||||||
\newcommand{\eKS}[1]{\epsilon^\text{KS}_{#1}}
|
\newcommand{\eKS}[1]{\epsilon^\text{KS}_{#1}}
|
||||||
\newcommand{\eQP}[1]{\epsilon^\text{QP}_{#1}}
|
\newcommand{\eQP}[1]{\epsilon^\text{QP}_{#1}}
|
||||||
\newcommand{\eGOWO}[1]{\epsilon^\text{\GOWO}_{#1}}
|
\newcommand{\eGOWO}[1]{\epsilon^\text{\GOWO}_{#1}}
|
||||||
\newcommand{\eGW}[1]{\epsilon^\text{\GW}_{#1}}
|
\newcommand{\eGW}[1]{\epsilon^{GW}_{#1}}
|
||||||
\newcommand{\eevGW}[1]{\epsilon^\text{\evGW}_{#1}}
|
\newcommand{\eevGW}[1]{\epsilon^\text{\evGW}_{#1}}
|
||||||
\newcommand{\eGnWn}[2]{\epsilon^\text{\GnWn{#2}}_{#1}}
|
\newcommand{\eGnWn}[2]{\epsilon^\text{\GnWn{#2}}_{#1}}
|
||||||
\newcommand{\Om}[2]{\Omega_{#1}^{#2}}
|
\newcommand{\Om}[2]{\Omega_{#1}^{#2}}
|
||||||
@ -101,7 +99,7 @@
|
|||||||
\newcommand{\Wc}[1]{W^\text{c}_{#1}}
|
\newcommand{\Wc}[1]{W^\text{c}_{#1}}
|
||||||
\newcommand{\vc}[1]{v_{#1}}
|
\newcommand{\vc}[1]{v_{#1}}
|
||||||
\newcommand{\Sig}[1]{\Sigma_{#1}}
|
\newcommand{\Sig}[1]{\Sigma_{#1}}
|
||||||
\newcommand{\SigGW}[1]{\Sigma^\text{\GW}_{#1}}
|
\newcommand{\SigGW}[1]{\Sigma^{GW}_{#1}}
|
||||||
\newcommand{\Z}[1]{Z_{#1}}
|
\newcommand{\Z}[1]{Z_{#1}}
|
||||||
\newcommand{\MO}[1]{\phi_{#1}}
|
\newcommand{\MO}[1]{\phi_{#1}}
|
||||||
\newcommand{\ERI}[2]{(#1|#2)}
|
\newcommand{\ERI}[2]{(#1|#2)}
|
||||||
@ -124,9 +122,9 @@
|
|||||||
\newcommand{\bSig}{\mathbf{\Sigma}}
|
\newcommand{\bSig}{\mathbf{\Sigma}}
|
||||||
\newcommand{\bSigX}{\mathbf{\Sigma}^\text{x}}
|
\newcommand{\bSigX}{\mathbf{\Sigma}^\text{x}}
|
||||||
\newcommand{\bSigC}{\mathbf{\Sigma}^\text{c}}
|
\newcommand{\bSigC}{\mathbf{\Sigma}^\text{c}}
|
||||||
\newcommand{\bSigGW}{\mathbf{\Sigma}^\text{\GW}}
|
\newcommand{\bSigGW}{\mathbf{\Sigma}^{GW}}
|
||||||
\newcommand{\be}{\mathbf{\epsilon}}
|
\newcommand{\be}{\mathbf{\epsilon}}
|
||||||
\newcommand{\beGW}{\mathbf{\epsilon}^\text{\GW}}
|
\newcommand{\beGW}{\mathbf{\epsilon}^{GW}}
|
||||||
\newcommand{\beGnWn}[1]{\mathbf{\epsilon}^\text{\GnWn{#1}}}
|
\newcommand{\beGnWn}[1]{\mathbf{\epsilon}^\text{\GnWn{#1}}}
|
||||||
\newcommand{\bde}{\mathbf{\Delta\epsilon}}
|
\newcommand{\bde}{\mathbf{\Delta\epsilon}}
|
||||||
\newcommand{\bdeHF}{\mathbf{\Delta\epsilon}^\text{HF}}
|
\newcommand{\bdeHF}{\mathbf{\Delta\epsilon}^\text{HF}}
|
||||||
@ -192,22 +190,22 @@ Moreover, when performed on top of a (partially) self-consistently {\evGW} calcu
|
|||||||
However, similar to adiabatic TD-DFT, \cite{Levine_2006,Tozer_2000,Huix-Rotllant_2010,Elliott_2011} the static version of BSE cannot describe multiple excitations. \cite{Romaniello_2009a,Sangalli_2011}
|
However, similar to adiabatic TD-DFT, \cite{Levine_2006,Tozer_2000,Huix-Rotllant_2010,Elliott_2011} the static version of BSE cannot describe multiple excitations. \cite{Romaniello_2009a,Sangalli_2011}
|
||||||
|
|
||||||
A significant limitation of the BSE formalism, as compared to TD-DFT, lies in the lack of analytic nuclear gradients (\ie, the first derivatives of the energy with respect to the nuclear displacements) for both the ground and excited states, \cite{Furche_2002} preventing efficient applications to the study of chemoluminescence, fluorescence and other related processes \cite{Navizet_2011} associated with geometric relaxation of ground and excited states, and structural changes upon electronic excitation. \cite{Bernardi_1996,Olivucci_2010,Robb_2007}
|
A significant limitation of the BSE formalism, as compared to TD-DFT, lies in the lack of analytic nuclear gradients (\ie, the first derivatives of the energy with respect to the nuclear displacements) for both the ground and excited states, \cite{Furche_2002} preventing efficient applications to the study of chemoluminescence, fluorescence and other related processes \cite{Navizet_2011} associated with geometric relaxation of ground and excited states, and structural changes upon electronic excitation. \cite{Bernardi_1996,Olivucci_2010,Robb_2007}
|
||||||
While calculations of the {\GW} quasiparticle energies ionic gradients is becoming popular,
|
While calculations of the $GW$ quasiparticle energies ionic gradients is becoming popular,
|
||||||
\cite{Lazzeri_2008,Faber_2011b,Yin_2013,Faber_2015,Montserrat_2016,Zhenglu_2019} only one pioneering study of the excited-state BSE gradients has been published so far. \cite{Beigi_2003} In this study devoted to small molecules (\ce{CO} and \ce{NH3}), only the BSE excitation energy gradients were calculated, while computing the KS-DFT (LDA) forces as its ground-state analog.
|
\cite{Lazzeri_2008,Faber_2011b,Yin_2013,Faber_2015,Montserrat_2016,Zhenglu_2019} only one pioneering study of the excited-state BSE gradients has been published so far. \cite{Beigi_2003} In this study devoted to small molecules (\ce{CO} and \ce{NH3}), only the BSE excitation energy gradients were calculated, while computing the KS-DFT (LDA) forces as its ground-state analog.
|
||||||
|
|
||||||
Contrary to TD-DFT which relies on Kohn-Sham density-functional theory (KS-DFT) \cite{Hohenberg_1964,Kohn_1965,ParrBook} as its ground-state counterpart, the BSE ground-state energy is not a well-defined quantity, and no clear consensus has been found regarding its formal definition.
|
Contrary to TD-DFT which relies on Kohn-Sham density-functional theory (KS-DFT) \cite{Hohenberg_1964,Kohn_1965,ParrBook} as its ground-state counterpart, the BSE ground-state energy is not a well-defined quantity, and no clear consensus has been found regarding its formal definition.
|
||||||
It then remains in its infancy with very few available studies for atomic and molecular systems. \cite{Olsen_2014,Holzer_2018,Li_2019,Li_2020}
|
It then remains in its infancy with very few available studies for atomic and molecular systems. \cite{Olsen_2014,Holzer_2018,Li_2019,Li_2020}
|
||||||
As a matter of fact, in the largest recent available benchmark study \cite{Holzer_2018} of the total energies of the atoms H?Ne, the atomization energies of the 26 small molecules forming the HEAT test set, \cite{Harding_2008} and the bond lengths and harmonic vibrational frequencies of $3d$ transition-metal monoxides, the BSE correlation energy, as evaluated within the adiabatic-connection fluctuation-dissipation (ACFDT) framework, \cite{Furche_2005} was mostly discarded from the set of tested techniques due to instabilities (negative frequency modes in the BSE polarization propagator) and replaced by an approximate (RPAsX) approach where the screened-Coulomb potential matrix elements was removed from the resonant electron-hole contribution. \cite{Maggio_2016,Holzer_2018}
|
As a matter of fact, in the largest recent available benchmark study \cite{Holzer_2018} of the total energies of the atoms H?Ne, the atomization energies of the 26 small molecules forming the HEAT test set, \cite{Harding_2008} and the bond lengths and harmonic vibrational frequencies of $3d$ transition-metal monoxides, the BSE correlation energy, as evaluated within the adiabatic-connection fluctuation-dissipation (ACFDT) framework, \cite{Furche_2005} was mostly discarded from the set of tested techniques due to instabilities (negative frequency modes in the BSE polarization propagator) and replaced by an approximate (RPAsX) approach where the screened-Coulomb potential matrix elements was removed from the resonant electron-hole contribution. \cite{Maggio_2016,Holzer_2018}
|
||||||
Such a modified BSE polarization propagator was inspired by a previous study on the homogeneous electron gas. \cite{Maggio_2016}
|
Such a modified BSE polarization propagator was inspired by a previous study on the homogeneous electron gas. \cite{Maggio_2016}
|
||||||
With such an approximation, amounting to neglect excitonic effects in the electron-hole propagator, the question of using either KS-DFT or {\GW} eigenvalues in the construction of the propagator becomes further relevant, increasing accordingly the number of possible definitions for the ground-state correlation energy.
|
With such an approximation, amounting to neglect excitonic effects in the electron-hole propagator, the question of using either KS-DFT or $GW$ eigenvalues in the construction of the propagator becomes further relevant, increasing accordingly the number of possible definitions for the ground-state correlation energy.
|
||||||
Finally, renormalizing or not the Coulomb interaction by the interaction strength $\IS$ in the Dyson equation for the interacting polarizability leads to two different versions of the BSE correlation energy. \cite{Holzer_2018}
|
Finally, renormalizing or not the Coulomb interaction by the interaction strength $\IS$ in the Dyson equation for the interacting polarizability leads to two different versions of the BSE correlation energy. \cite{Holzer_2018}
|
||||||
|
|
||||||
Here, in analogy to random-phase approximation (RPA)-type formalisms \cite{Furche_2008,Toulouse_2009,Toulouse_2010,Angyan_2011,Ren_2012} and similarly to Refs.~\onlinecite{Olsen_2014,Maggio_2016,Holzer_2018}, the ground-state BSE energy is calculated in the adiabatic connection framework.
|
Here, in analogy to random-phase approximation (RPA)-type formalisms \cite{Furche_2008,Toulouse_2009,Toulouse_2010,Angyan_2011,Ren_2012} and similarly to Refs.~\onlinecite{Olsen_2014,Maggio_2016,Holzer_2018}, the ground-state BSE energy is calculated in the adiabatic connection framework.
|
||||||
Embracing this definition, the purpose of the present study is to investigate the quality of ground-state PES near equilibrium obtained within the BSE approach for several diatomic molecules.
|
Embracing this definition, the purpose of the present study is to investigate the quality of ground-state PES near equilibrium obtained within the BSE approach for several diatomic molecules.
|
||||||
The location of the minima on the ground-state PES is of particular interest.
|
The location of the minima on the ground-state PES is of particular interest.
|
||||||
This study is a first preliminary step towards the development of analytical nuclear gradients within the BSE@{\GW} formalism.
|
This study is a first preliminary step towards the development of analytical nuclear gradients within the BSE@$GW$ formalism.
|
||||||
Thanks to comparison with both similar and state-of-art computational approaches, we show that the ACFDT@BSE@{\GW} approach is surprisingly accurate, and can even compete with high-order coupled cluster (CC) methods in terms of absolute energies.
|
Thanks to comparison with both similar and state-of-art computational approaches, we show that the ACFDT@BSE@$GW$ approach is surprisingly accurate, and can even compete with high-order coupled cluster (CC) methods in terms of absolute energies.
|
||||||
However, we also observe, in some cases, unphysical irregularities on the ground-state PES, which are due to the appearance of a satellite resonance with a weight similar to that of the {\GW} quasiparticle peak. \cite{vanSetten_2015,Maggio_2017,Loos_2018,Veril_2018,Duchemin_2020}
|
However, we also observe, in some cases, unphysical irregularities on the ground-state PES, which are due to the appearance of a satellite resonance with a weight similar to that of the $GW$ quasiparticle peak. \cite{vanSetten_2015,Maggio_2017,Loos_2018,Veril_2018,Duchemin_2020}
|
||||||
|
|
||||||
%The paper is organized as follows.
|
%The paper is organized as follows.
|
||||||
%In Sec.~\ref{sec:theo}, we introduce the equations behind the BSE formalism.
|
%In Sec.~\ref{sec:theo}, we introduce the equations behind the BSE formalism.
|
||||||
@ -243,7 +241,7 @@ which takes into account the self-consistent variation of the Hartree potential
|
|||||||
\end{equation}
|
\end{equation}
|
||||||
(where $\vc{}$ is the bare Coulomb operator) and the exchange-correlation self-energy $\Sig{\xc}$.
|
(where $\vc{}$ is the bare Coulomb operator) and the exchange-correlation self-energy $\Sig{\xc}$.
|
||||||
In Eq.~\eqref{eq:BSE}, $\LBSE{0}(1,2,1',2') = - i \G{}(1,2')\G{}(2,1')$, and $(1)=(\br{1},\sigma_1,t_1)$ is a composite index gathering space, spin and time variables.
|
In Eq.~\eqref{eq:BSE}, $\LBSE{0}(1,2,1',2') = - i \G{}(1,2')\G{}(2,1')$, and $(1)=(\br{1},\sigma_1,t_1)$ is a composite index gathering space, spin and time variables.
|
||||||
In the {\GW} approximation, \cite{Hedin_1965,Aryasetiawan_1998,Onida_2002,Martin_2016,Reining_2017} we have
|
In the $GW$ approximation, \cite{Hedin_1965,Aryasetiawan_1998,Onida_2002,Martin_2016,Reining_2017} we have
|
||||||
\begin{equation}
|
\begin{equation}
|
||||||
\SigGW{\xc}(1,2) = i \G{}(1,2) \W{}{}(1^+,2),
|
\SigGW{\xc}(1,2) = i \G{}(1,2) \W{}{}(1^+,2),
|
||||||
\end{equation}
|
\end{equation}
|
||||||
@ -300,10 +298,10 @@ In the case of BSE, the specific expression of the matrix elements are
|
|||||||
\label{eq:LR_BSE}
|
\label{eq:LR_BSE}
|
||||||
\ABSE{ia,jb}{\IS} & = \delta_{ij} \delta_{ab} (\eGW{a} - \eGW{i}) + \IS \W{ia,bj}{\IS}(\omega = 0) - \lambda \ERI{ia}{jb},
|
\ABSE{ia,jb}{\IS} & = \delta_{ij} \delta_{ab} (\eGW{a} - \eGW{i}) + \IS \W{ia,bj}{\IS}(\omega = 0) - \lambda \ERI{ia}{jb},
|
||||||
\\
|
\\
|
||||||
\BBSE{ia,jb}{\IS} & = \W{ia,jb}{\IS}(\omega = 0) - \lambda \ERI{ib}{ja} ,
|
\BBSE{ia,jb}{\IS} & = \W{ia,jb}{\IS}(\omega = 0) - \lambda \ERI{ia}{bj} ,
|
||||||
\end{align}
|
\end{align}
|
||||||
\end{subequations}
|
\end{subequations}
|
||||||
where $\eGW{p}$ are the {\GW} quasiparticle energies,
|
where $\eGW{p}$ are the $GW$ quasiparticle energies,
|
||||||
\begin{multline}
|
\begin{multline}
|
||||||
\label{eq:W}
|
\label{eq:W}
|
||||||
\W{ia,jb}{\IS}(\omega) = 2 \ERI{ia}{jb}
|
\W{ia,jb}{\IS}(\omega) = 2 \ERI{ia}{jb}
|
||||||
@ -323,9 +321,9 @@ In Eq.~\eqref{eq:W}, $\eta$ is a positive infinitesimal, and $\OmRPA{m}{\IS}$ a
|
|||||||
\begin{subequations}
|
\begin{subequations}
|
||||||
\begin{align}
|
\begin{align}
|
||||||
\label{eq:LR_RPA}
|
\label{eq:LR_RPA}
|
||||||
\ARPA{ia,jb}{\IS} & = \delta_{ij} \delta_{ab} (\eHF{a} - \eHF{i}) + \IS \ERI{ia}{bj},
|
\ARPA{ia,jb}{\IS} & = \delta_{ij} \delta_{ab} (\eHF{a} - \eHF{i}) + 2 \IS \ERI{ia}{bj},
|
||||||
\\
|
\\
|
||||||
\BRPA{ia,jb}{\IS} & = \IS \ERI{ia}{jb},
|
\BRPA{ia,jb}{\IS} & = 2 \IS \ERI{ia}{jb},
|
||||||
\end{align}
|
\end{align}
|
||||||
\end{subequations}
|
\end{subequations}
|
||||||
where $\eHF{p}$ are the HF orbital energies.
|
where $\eHF{p}$ are the HF orbital energies.
|
||||||
@ -369,11 +367,24 @@ is the interaction kernel \cite{Angyan_2011, Holzer_2018} [with $\tA{ia,jb}{\IS}
|
|||||||
\end{pmatrix}
|
\end{pmatrix}
|
||||||
\end{equation}
|
\end{equation}
|
||||||
is the correlation part of the two-electron density matrix at interaction strength $\IS$, and $\Tr$ denotes the matrix trace.
|
is the correlation part of the two-electron density matrix at interaction strength $\IS$, and $\Tr$ denotes the matrix trace.
|
||||||
Note that the present definition of the BSE correlation energy [see Eq.~\eqref{eq:EcBSE}] has be labeled ``XBS'' for ``extended Bethe Salpeter'' by Holzer et al. \cite{Holzer_2018}
|
Note that the present definition of the BSE correlation energy [see Eq.~\eqref{eq:EcBSE}], which we refer to as BSE@$GW$@HF here, has be labeled ``XBS'' for ``extended Bethe Salpeter'' by Holzer et al. \cite{Holzer_2018}
|
||||||
Contrary to DFT where the electron density is fixed along the adiabatic path, in the present formalism, the density is not maintain with respect to $\IS$.
|
Contrary to DFT where the electron density is fixed along the adiabatic path, in the present formalism, the density is not maintain with respect to $\IS$.
|
||||||
Therefore, an additional contribution to Eq.~\eqref{eq:EcBSE} originating from the variation of the Green's function along the adiabatic connection should be added.
|
Therefore, an additional contribution to Eq.~\eqref{eq:EcBSE} originating from the variation of the Green's function along the adiabatic connection should be added.
|
||||||
However, as commonly done within RPA and RPAx (\ie, RPA with exchange), \cite{Toulouse_2009, Toulouse_2010, Holzer_2018} we shall neglect it in the present study.
|
However, as commonly done within RPA and RPAx (\ie, RPA with exchange), \cite{Toulouse_2009, Toulouse_2010, Holzer_2018} we shall neglect it in the present study.
|
||||||
|
|
||||||
|
Equation \eqref{eq:EcBSE} can also be straightforwardly applied to RPA and RPAx, the only difference being the expressions of $\bA{\IS}$ and $\bB{\IS}$ used to obtain the eigenvectors $\bX{\IS}$ and $\bY{\IS}$ entering the definition of $\bP{\IS}$ [see Eq.~\eqref{eq:2DM}].
|
||||||
|
For RPA, these expressions have been provided in Eq.~\eqref{eq:LR_RPA}, and their RPAx analogs read
|
||||||
|
\begin{subequations}
|
||||||
|
\begin{align}
|
||||||
|
\label{eq:LR_RPAx}
|
||||||
|
\ARPAx{ia,jb}{\IS} & = \delta_{ij} \delta_{ab} (\eHF{a} - \eHF{i}) + 2 \IS \ERI{ia}{bj} - \IS \ERI{ia}{jb},
|
||||||
|
\\
|
||||||
|
\BRPAx{ia,jb}{\IS} & = 2 \IS \ERI{ia}{jb} - \IS \ERI{ia}{bj}.
|
||||||
|
\end{align}
|
||||||
|
\end{subequations}
|
||||||
|
These two types of calculations will be refer to as RPA@HF and RPAx@HF respectively in the following.
|
||||||
|
Finally, we will also consider the RPA@$GW$@HF scheme which consists in replacing the HF orbital energies in Eq.~\eqref{eq:LR_RPA} by the $GW$ quasiparticles energies.
|
||||||
|
|
||||||
Several important comments are in order here.
|
Several important comments are in order here.
|
||||||
For spin-restricted closed-shell molecular systems around their equilibrium geometry (such as the ones studied here), it is rare to encounter singlet instabilities as these systems can be labeled as weakly correlated.
|
For spin-restricted closed-shell molecular systems around their equilibrium geometry (such as the ones studied here), it is rare to encounter singlet instabilities as these systems can be labeled as weakly correlated.
|
||||||
However, singlet instabilities may appear in the presence of strong correlation (\eg, when the bond is stretched).
|
However, singlet instabilities may appear in the presence of strong correlation (\eg, when the bond is stretched).
|
||||||
@ -388,15 +399,14 @@ However, to the best of our knowledge, such alternative plasmon expression does
|
|||||||
%\section{Computational details}
|
%\section{Computational details}
|
||||||
%\label{sec:comp_details}
|
%\label{sec:comp_details}
|
||||||
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
||||||
All the {\GW} calculations performed to obtain the screened Coulomb operator and the quasiparticle energies have been done using a (restricted) Hartree-Fock (HF) starting point, which is a very adequate choice in the case of the (small) systems that we have considered here.
|
All the $GW$ calculations performed to obtain the screened Coulomb operator and the quasiparticle energies have been done using a (restricted) Hartree-Fock (HF) starting point, which is a very adequate choice in the case of the (small) systems that we have considered here.
|
||||||
Perturbative {\GW} (or {\GOWO}) \cite{Hybertsen_1985a, Hybertsen_1986} calculations are employed as starting point to compute the BSE neutral excitations.
|
Perturbative $GW$ (or {\GOWO}) \cite{Hybertsen_1985a, Hybertsen_1986} calculations are employed as starting point to compute the BSE neutral excitations.
|
||||||
These will be labeled as BSE@{\GOWO}.
|
|
||||||
In the case of {\GOWO}, the quasiparticle energies have been obtained by linearizing the non-linear, frequency-dependent quasiparticle equation.
|
In the case of {\GOWO}, the quasiparticle energies have been obtained by linearizing the non-linear, frequency-dependent quasiparticle equation.
|
||||||
Further details about our implementation of {\GOWO} can be found in Refs.~\onlinecite{Loos_2018, Veril_2018}.
|
Further details about our implementation of {\GOWO} can be found in Refs.~\onlinecite{Loos_2018, Veril_2018}.
|
||||||
Finally, the infinitesimal $\eta$ has been set to zero for all calculations.
|
Finally, the infinitesimal $\eta$ has been set to zero for all calculations.
|
||||||
The numerical integration required to compute the correlation energy along the adiabatic path [see Eq.~\eqref{eq:EcBSE}] has been performed with a 21-point Gauss-Legendre quadrature.
|
The numerical integration required to compute the correlation energy along the adiabatic path [see Eq.~\eqref{eq:EcBSE}] has been performed with a 21-point Gauss-Legendre quadrature.
|
||||||
For comparison purposes, we have also computed the PES at the MP2, CC2 \cite{Christiansen_1995}, CCSD, \cite{Purvis_1982} and CC3 \cite{Christiansen_1995b} levels of theory using DALTON. \cite{dalton}
|
For comparison purposes, we have also computed the PES at the MP2, CC2 \cite{Christiansen_1995}, CCSD, \cite{Purvis_1982} and CC3 \cite{Christiansen_1995b} levels of theory using DALTON. \cite{dalton}
|
||||||
All the other calculations have been performed with our locally developed {\GW} software. \cite{Loos_2018,Veril_2018}
|
All the other calculations have been performed with our locally developed $GW$ software. \cite{Loos_2018,Veril_2018}
|
||||||
As one-electron basis sets, we employ the Dunning family (cc-pVXZ) defined with cartesian gaussian functions.
|
As one-electron basis sets, we employ the Dunning family (cc-pVXZ) defined with cartesian gaussian functions.
|
||||||
Unless, otherwise stated, the frozen-core approximation has not been enforced in our calculations in order to provide a fairer comparison between methods.
|
Unless, otherwise stated, the frozen-core approximation has not been enforced in our calculations in order to provide a fairer comparison between methods.
|
||||||
However, we have found that the conclusions drawn in the present study hold within the frozen-core approximation (see {\SI}).
|
However, we have found that the conclusions drawn in the present study hold within the frozen-core approximation (see {\SI}).
|
||||||
@ -410,27 +420,9 @@ However, we are currently pursuing different avenues to lower this cost by compu
|
|||||||
%\label{sec:PES}
|
%\label{sec:PES}
|
||||||
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
||||||
In order to illustrate the performance of the BSE-based adiabatic connection formulation, we have computed the ground-state PES of several closed-shell diatomic molecules around their equilibrium geometry: \ce{H2}, \ce{LiH}, \ce{LiF}, \ce{N2}, \ce{CO}, \ce{BF}, \ce{F2}, and \ce{HCl}.
|
In order to illustrate the performance of the BSE-based adiabatic connection formulation, we have computed the ground-state PES of several closed-shell diatomic molecules around their equilibrium geometry: \ce{H2}, \ce{LiH}, \ce{LiF}, \ce{N2}, \ce{CO}, \ce{BF}, \ce{F2}, and \ce{HCl}.
|
||||||
The PES of these molecules for various methods and Dunning's triple-$\zeta$ basis cc-pVTZ are represented in Fig.~\ref{fig:PES}, and the computed equilibrium distances are gathered in Table \ref{tab:Req}.
|
The PES of these molecules for various methods and Dunning's triple-$\zeta$ basis cc-pVTZ are represented in Figs.~\ref{fig:PES-H2-LiH}, \ref{fig:PES-LiF-HCl}, \ref{fig:PES-N2-CO-BF}, and \ref{fig:PES-F2}, while the computed equilibrium distances are gathered in Table \ref{tab:Req}.
|
||||||
Additional graphs for other basis sets can be found in the {\SI}.
|
Additional graphs for other basis sets can be found in the {\SI}.
|
||||||
|
|
||||||
%%% FIG 1 %%%
|
|
||||||
\begin{figure*}
|
|
||||||
\includegraphics[width=0.45\linewidth]{H2_GS_VTZ}
|
|
||||||
\includegraphics[width=0.45\linewidth]{LiH_GS_VTZ}
|
|
||||||
\includegraphics[width=0.45\linewidth]{LiF_GS_VTZ}
|
|
||||||
\includegraphics[width=0.45\linewidth]{N2_GS_VTZ}
|
|
||||||
\includegraphics[width=0.45\linewidth]{CO_GS_VTZ}
|
|
||||||
\includegraphics[width=0.45\linewidth]{BF_GS_VTZ}
|
|
||||||
\includegraphics[width=0.45\linewidth]{F2_GS_VTZ}
|
|
||||||
\includegraphics[width=0.45\linewidth]{HCl_GS_VTZ}
|
|
||||||
\caption{
|
|
||||||
Ground-state PES of diatomic molecules around their respective equilibrium geometry obtained at various levels of theory with the cc-pVTZ basis set.
|
|
||||||
Additional graphs for other basis sets and within the frozen-core approximation can be found in the {\SI}.
|
|
||||||
\label{fig:PES}
|
|
||||||
}
|
|
||||||
\end{figure*}
|
|
||||||
%%% %%% %%%
|
|
||||||
|
|
||||||
%%% TABLE I %%%
|
%%% TABLE I %%%
|
||||||
\begin{table*}
|
\begin{table*}
|
||||||
\caption{
|
\caption{
|
||||||
@ -455,17 +447,17 @@ Equilibrium distances (in bohr) of the ground state of diatomic molecules obtain
|
|||||||
MP2 & cc-pVDZ & 1.426 & 3.041 & 3.010 & 2.133 & 2.166 & 2.431 & 2.681 & 2.426 \\
|
MP2 & cc-pVDZ & 1.426 & 3.041 & 3.010 & 2.133 & 2.166 & 2.431 & 2.681 & 2.426 \\
|
||||||
& cc-pVTZ & 1.393 & 3.004 & 2.968 & 2.095 & 2.144 & 2.383 & 2.636 & 2.405 \\
|
& cc-pVTZ & 1.393 & 3.004 & 2.968 & 2.095 & 2.144 & 2.383 & 2.636 & 2.405 \\
|
||||||
& cc-pVQZ & 1.391 & 3.008 & 2.970 & 2.091 & 2.137 & 2.382 & 2.634 & 2.395 \\
|
& cc-pVQZ & 1.391 & 3.008 & 2.970 & 2.091 & 2.137 & 2.382 & 2.634 & 2.395 \\
|
||||||
BSE@{\GOWO}@HF & cc-pVDZ & 1.437 & 3.042 & 3.000 & 2.107 & 2.153 & 2.407 & 2.700 & >2.440 \\
|
BSE@{\GOWO}@HF & cc-pVDZ & 1.437 & 3.042 & 3.000 & 2.107 & 2.153 & 2.407 & 2.700 & >2.440 \\
|
||||||
& cc-pVTZ & 1.404 & 3.023 & glitch & & & <2.420 & & <2.410 \\
|
& cc-pVTZ & 1.404 & 3.023 & glitch & & & <2.420 & & <2.410 \\
|
||||||
& cc-pVQZ & 1.399 & & & & & & & \\
|
& cc-pVQZ & 1.399 & & & & & & & \\
|
||||||
RPA@{\GOWO}@HF & cc-pVDZ & 1.426 & 3.019 & 2.994 & 2.083 & 2.144 & 2.403 & 2.691 & 2.436 \\
|
RPA@{\GOWO}@HF & cc-pVDZ & 1.426 & 3.019 & 2.994 & 2.083 & 2.144 & 2.403 & 2.691 & 2.436 \\
|
||||||
& cc-pVTZ & 1.388 & 3.013 & glitch & & & <2.420 & & <2.410 \\
|
& cc-pVTZ & 1.388 & 3.013 & glitch & & & <2.420 & & <2.410 \\
|
||||||
& cc-pVQZ & 1.382 & & & & & & & \\
|
& cc-pVQZ & 1.382 & & & & & & & \\
|
||||||
RPAx@HF & cc-pVDZ & 1.428 & 3.040 & 2.998 & 2.077 & 2.130 & 2.417 & NaN & 2.424 \\
|
RPAx@HF & cc-pVDZ & 1.428 & 3.040 & 2.998 & 2.077 & 2.130 & 2.417 & NaN & 2.424 \\
|
||||||
& cc-pVTZ & 1.395 & 3.003 & <2.990 & & & <2.420 & & <2.410 \\
|
& cc-pVTZ & 1.395 & 3.003 & <2.990 & & & <2.420 & & <2.410 \\
|
||||||
& cc-pVQZ & 1.394 & & & & & & & \\
|
& cc-pVQZ & 1.394 & & & & & & & \\
|
||||||
RPA@HF & cc-pVDZ & 1.431 & 3.021 & 2.999 & 2.083 & 2.134 & & 2.623 & 2.424 \\
|
RPA@HF & cc-pVDZ & 1.431 & 3.021 & 2.999 & 2.083 & 2.134 & & 2.623 & 2.424 \\
|
||||||
& cc-pVTZ & 1.388 & 2.978 & <2.990 & & & 2.416 & & <2.410 \\
|
& cc-pVTZ & 1.388 & 2.978 & <2.990 & & & 2.416 & & <2.410 \\
|
||||||
& cc-pVQZ & 1.386 & & & & & <2.420 & & \\
|
& cc-pVQZ & 1.386 & & & & & <2.420 & & \\
|
||||||
% FROZEN CORE VERSION
|
% FROZEN CORE VERSION
|
||||||
% Method & Basis & \ce{H2} & \ce{LiH}& \ce{LiF}& \ce{N2} & \ce{CO} & \ce{BF} & \ce{F2} & \ce{HCl}\\
|
% Method & Basis & \ce{H2} & \ce{LiH}& \ce{LiF}& \ce{N2} & \ce{CO} & \ce{BF} & \ce{F2} & \ce{HCl}\\
|
||||||
@ -499,20 +491,76 @@ Equilibrium distances (in bohr) of the ground state of diatomic molecules obtain
|
|||||||
\end{ruledtabular}
|
\end{ruledtabular}
|
||||||
\end{table*}
|
\end{table*}
|
||||||
|
|
||||||
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
Let us first start with the two smallest molecules, \ce{H2} and \ce{LiH} which are both linked by covalent bonds (see Fig.~\ref{fig:PES-H2-LiH}).
|
||||||
%\subsection{Hydrogen molecule}
|
For \ce{H2}, we take as reference the full configuration interaction (FCI) energies and we also report the MP2 curve and its third-order variant (MP3), which improves upon MP2 towards FCI.
|
||||||
%\label{sec:H2}
|
RPA@HF and RPA@{\GOWO}@HF yield almost identical results, and significantly overestimate (in absolute value) the FCI energy, while RPAx@HF and BSE@{\GOWO}@HF slightly underestimate and overestimate the FCI energy, respectively, RPAx@HF being the best match in the case of \ce{H2}.
|
||||||
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
Interestingly though, the BSE@{\GOWO}@HF scheme yields a more accurate equilibrium bond length than any other method irrespectively of the basis set.
|
||||||
|
For example, with the cc-pVQZ basis set, BSE@{\GOWO}@HF is only off by $0.003$ bohr compared to FCI, while RPAx@HF, MP2, and CC2 underestimate the bond length by $0.008$, $0.011$, and $0.011$ bohr, respectively.
|
||||||
|
The RPA-based schemes are much less accurate, with even shorter equilibrium bond lengths.
|
||||||
|
|
||||||
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
The scenario is almost identical for \ce{LiH} for which we report the CC2, CCSD and CC3 energies in addition to MP2.
|
||||||
%\subsection{Lithium hydride and lithium fluoride}
|
In this case, RPAx@HF and BSE@{\GOWO}@HF nestle the CCSD and CC3 energy curves, and they are almost perfectly parallel.
|
||||||
%\label{sec:LiX}
|
|
||||||
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
%%% FIG 1 %%%
|
||||||
|
\begin{figure*}
|
||||||
|
\includegraphics[width=0.49\linewidth]{H2_GS_VTZ}
|
||||||
|
\includegraphics[width=0.49\linewidth]{LiH_GS_VTZ}
|
||||||
|
\caption{
|
||||||
|
Ground-state PES of \ce{H2} (left) and \ce{LiH} (right) around their respective equilibrium geometry obtained at various levels of theory with the cc-pVTZ basis set.
|
||||||
|
Additional graphs for other basis sets and within the frozen-core approximation can be found in the {\SI}.
|
||||||
|
\label{fig:PES-H2-LiH}
|
||||||
|
}
|
||||||
|
\end{figure*}
|
||||||
|
%%% %%% %%%
|
||||||
|
|
||||||
|
The cases of \ce{LiF} and \ce{HCl} (see Fig.~\ref{fig:PES-LiF-HCl}) are interesting as they corresponds to a strongly polarized bond towards the halogen atoms which are much more electronegative than the first row elements.
|
||||||
|
For these ionic bond, the performance of BSE@{\GOWO}@HF are terrific with an almost perfect match to the CC3 curve.
|
||||||
|
For \ce{LiF}, the two curves starting to deviate a few tenths of bohr after the equilibrium geometry, but they remain tightly ... for much longer in the case of \ce{HCl}.
|
||||||
|
Maybe surprisingly, BSE@{\GOWO}@HF outperforms both CC2 and CCSD, as well as RPAx@HF by a big margin for these two molecules exhibiting charge transfer.
|
||||||
|
However, in the case of \ce{LiF}, the attentive reader would have observed a small glitch in the $GW$-based curves very close to their minimum.
|
||||||
|
|
||||||
|
%%% FIG 2 %%%
|
||||||
|
\begin{figure*}
|
||||||
|
\includegraphics[width=0.49\linewidth]{LiF_GS_VTZ}
|
||||||
|
\includegraphics[width=0.49\linewidth]{HCl_GS_VTZ}
|
||||||
|
\caption{
|
||||||
|
Ground-state PES of \ce{LiF} (left) and \ce{HCl} (right) around their respective equilibrium geometry obtained at various levels of theory with the cc-pVTZ basis set.
|
||||||
|
Additional graphs for other basis sets and within the frozen-core approximation can be found in the {\SI}.
|
||||||
|
\label{fig:PES-LiF-HCl}
|
||||||
|
}
|
||||||
|
\end{figure*}
|
||||||
|
%%% %%% %%%
|
||||||
|
|
||||||
|
Let us now look at the isoelectronic series \ce{N2}, \ce{CO}, and \ce{BF}, which have a decreasing bond order (from triple bond to single bond).
|
||||||
|
In that case again, the performance of BSE@{\GOWO}@HF are outstanding as shown in Fig.~\ref{fig:PES-N2-CO-BF}.
|
||||||
|
|
||||||
|
|
||||||
|
%%% FIG 3 %%%
|
||||||
|
\begin{figure*}
|
||||||
|
\includegraphics[width=0.33\linewidth]{N2_GS_VTZ}
|
||||||
|
\includegraphics[width=0.33\linewidth]{CO_GS_VTZ}
|
||||||
|
\includegraphics[width=0.33\linewidth]{BF_GS_VTZ}
|
||||||
|
\caption{
|
||||||
|
Ground-state PES of the isoelectronic series \ce{N2} (left), \ce{CO} (center), and \ce{BF} (right) around their respective equilibrium geometry obtained at various levels of theory with the cc-pVTZ basis set.
|
||||||
|
Additional graphs for other basis sets and within the frozen-core approximation can be found in the {\SI}.
|
||||||
|
\label{fig:PES-N2-CO-BF}
|
||||||
|
}
|
||||||
|
\end{figure*}
|
||||||
|
%%% %%% %%%
|
||||||
|
|
||||||
|
The \ce{F2} molecule is a notoriously difficult case to treat due to the relative weakness of its covalent bond (see Fig.~\ref{fig:PES-F2}).
|
||||||
|
|
||||||
|
%%% FIG 4 %%%
|
||||||
|
\begin{figure}
|
||||||
|
\includegraphics[width=\linewidth]{F2_GS_VTZ}
|
||||||
|
\caption{
|
||||||
|
Ground-state PES of \ce{F2} around its equilibrium geometry obtained at various levels of theory with the cc-pVTZ basis set.
|
||||||
|
Additional graphs for other basis sets and within the frozen-core approximation can be found in the {\SI}.
|
||||||
|
\label{fig:PES-F2}
|
||||||
|
}
|
||||||
|
\end{figure}
|
||||||
|
%%% %%% %%%
|
||||||
|
|
||||||
%%%%%%%%%%%%%%%%%%%%%%%%
|
|
||||||
%\subsection{The isoelectronic sequence: \ce{N2}, \ce{CO}, and \ce{BF}}
|
|
||||||
%\label{sec:isoN2}
|
|
||||||
%%%%%%%%%%%%%%%%%%%%%%%%
|
|
||||||
|
|
||||||
%%%%%%%%%%%%%%%%%%%%%%%%
|
%%%%%%%%%%%%%%%%%%%%%%%%
|
||||||
%\section{Conclusion}
|
%\section{Conclusion}
|
||||||
@ -520,10 +568,10 @@ Equilibrium distances (in bohr) of the ground state of diatomic molecules obtain
|
|||||||
%%%%%%%%%%%%%%%%%%%%%%%%
|
%%%%%%%%%%%%%%%%%%%%%%%%
|
||||||
In this Letter, we have shown that calculating the BSE correlation energy in the ACFDT framework yield extremely accurate PES around equilibrium.
|
In this Letter, we have shown that calculating the BSE correlation energy in the ACFDT framework yield extremely accurate PES around equilibrium.
|
||||||
We have illustrated this for 8 diatomic molecules for which we have also computed reference ground-state energies using coupled cluster methods (CC2, CCSD, and CC3).
|
We have illustrated this for 8 diatomic molecules for which we have also computed reference ground-state energies using coupled cluster methods (CC2, CCSD, and CC3).
|
||||||
However, we have also observed that, in some cases, unphysical irregularities on the ground-state PES due to the appearance of a satellite resonance with a weight similar to that of the {\GW} quasiparticle peak.
|
However, we have also observed that, in some cases, unphysical irregularities on the ground-state PES due to the appearance of a satellite resonance with a weight similar to that of the $GW$ quasiparticle peak.
|
||||||
This shortcoming, which is entirely due to the quasiparticle nature of the underlying {\GW} calculation, has been thoroughly described in several previous studies.\cite{vanSetten_2015,Maggio_2017,Loos_2018,Veril_2018,Duchemin_2020}
|
This shortcoming, which is entirely due to the quasiparticle nature of the underlying $GW$ calculation, has been thoroughly described in several previous studies.\cite{vanSetten_2015,Maggio_2017,Loos_2018,Veril_2018,Duchemin_2020}
|
||||||
We believe that this central issue must be resolved if one wants to expand the applicability of the present methods.
|
We believe that this central issue must be resolved if one wants to expand the applicability of the present methods.
|
||||||
In the perspective of developing analytical nuclear gradients within the BSE@{\GW} formalism, we are currently investigating the accuracy of the ACFDT@BSE scheme for excited-state PES.
|
In the perspective of developing analytical nuclear gradients within the BSE@$GW$ formalism, we are currently investigating the accuracy of the ACFDT@BSE scheme for excited-state PES.
|
||||||
We hope to be able to report on this in the near future.
|
We hope to be able to report on this in the near future.
|
||||||
|
|
||||||
%%%%%%%%%%%%%%%%%%%%%%%%
|
%%%%%%%%%%%%%%%%%%%%%%%%
|
||||||
|
Loading…
Reference in New Issue
Block a user