saving work before pull
This commit is contained in:
parent
9078bbea75
commit
e12c8214c0
@ -502,6 +502,16 @@ It is worth noting the close similarity of the first-order elements with the one
|
||||
\subsection{Second-order matrix elements}
|
||||
% ///////////////////////////%
|
||||
|
||||
%%% FIG 1 %%%
|
||||
\begin{figure*}
|
||||
\centering
|
||||
\includegraphics[width=\linewidth]{fig1.pdf}
|
||||
\caption{
|
||||
Functional form of the qs$GW$ self-energy (left) for $\eta = 1$ and the SRG-qs$GW$ self-energy (right) for $s = 1/(2\eta^2) = 1/2$.
|
||||
\label{fig:plot}}
|
||||
\end{figure*}
|
||||
%%% %%% %%% %%%
|
||||
|
||||
The second-order renormalized quasiparticle equation is given by
|
||||
\begin{equation}
|
||||
\label{eq:GW_renorm}
|
||||
@ -546,16 +556,6 @@ while the dynamic part of the self-energy [see Eq.~\eqref{eq:srg_sigma}] tends t
|
||||
Therefore, the SRG flow continuously transforms the dynamical self-energy $\widetilde{\bSig}(\omega; s)$ into a static correction $\widetilde{\bF}^{(2)}(s)$.
|
||||
This transformation is done gradually starting from the states that have the largest denominators in Eq.~\eqref{eq:static_F2}.
|
||||
|
||||
%%% FIG 1 %%%
|
||||
\begin{figure*}
|
||||
\centering
|
||||
\includegraphics[width=\linewidth]{fig1.pdf}
|
||||
\caption{
|
||||
Functional form of the qs$GW$ self-energy (left) for $\eta = 1$ and the SRG-qs$GW$ self-energy (right) for $s = 1/(2\eta^2) = 1/2$.
|
||||
\label{fig:plot}}
|
||||
\end{figure*}
|
||||
%%% %%% %%% %%%
|
||||
|
||||
%///////////////////////////%
|
||||
\subsection{Alternative form of the static self-energy}
|
||||
% ///////////////////////////%
|
||||
@ -699,17 +699,6 @@ However, as we will see in the next subsection these are just particular molecul
|
||||
Also, the SRG-qs$GW_\TDA$ is better than qs$GW_\TDA$ in the three cases of Fig.~\ref{fig:fig3} but this is the other way around.
|
||||
Therefore, it seems that the effect of the TDA can not be systematically predicted.
|
||||
|
||||
|
||||
|
||||
%%%%%%%%%%%%%%%%%%%%%%
|
||||
\subsection{Statistical analysis}
|
||||
\label{sec:SRG_vs_Sym}
|
||||
%%%%%%%%%%%%%%%%%%%%%%
|
||||
|
||||
The test set considered in this study is composed of the GW20 set of molecules introduced by Lewis and Berkelbach. \cite{Lewis_2019}
|
||||
This set is made of the 20 smallest atoms and molecules of the GW100 benchmark set. \cite{vanSetten_2015}
|
||||
In addition, the MgO and O3 molecules (which are part of GW100 as well) has been added to the test set because they are known to possess intruder states. \cite{vanSetten_2015,Forster_2021}
|
||||
|
||||
\begin{table}
|
||||
\centering
|
||||
\caption{First ionization potential in eV calculated using $\Delta$CCSD(T) (reference), HF, $G_0W_0$@HF, qs$GW$ and SRG-qs$GW$. The statistical descriptors are computed for the errors with respect to the reference.}
|
||||
@ -754,6 +743,11 @@ In addition, the MgO and O3 molecules (which are part of GW100 as well) has been
|
||||
\end{tabular}
|
||||
\end{table}
|
||||
|
||||
%%%%%%%%%%%%%%%%%%%%%%
|
||||
\subsection{Statistical analysis}
|
||||
\label{sec:SRG_vs_Sym}
|
||||
%%%%%%%%%%%%%%%%%%%%%%
|
||||
|
||||
%%% FIG 4 %%%
|
||||
\begin{figure*}
|
||||
\centering
|
||||
@ -764,6 +758,10 @@ In addition, the MgO and O3 molecules (which are part of GW100 as well) has been
|
||||
\end{figure*}
|
||||
%%% %%% %%% %%%
|
||||
|
||||
The test set considered in this study is composed of the GW20 set of molecules introduced by Lewis and Berkelbach. \cite{Lewis_2019}
|
||||
This set is made of the 20 smallest atoms and molecules of the GW100 benchmark set. \cite{vanSetten_2015}
|
||||
In addition, the MgO and O3 molecules (which are part of GW100 as well) has been added to the test set because they are known to possess intruder states. \cite{vanSetten_2015,Forster_2021}
|
||||
|
||||
Table~\ref{tab:tab1} shows the principal IP of the 22 molecules considered in this work computed at various level of theories.
|
||||
As mentioned previously the HF IPs are overestimated with a mean signed error (MSE) and mean absolute error (MAE) of \SI{0.64}{\electronvolt} and \SI{0.74}{\electronvolt}, respectively.
|
||||
Performing a one-shot $G_0W_0$ calculation on top of it allows to divided by more than 2 the MSE and MAE, \SI{0.26}{\electronvolt} and \SI{0.32}{\electronvolt}, respectively.
|
||||
@ -773,23 +771,30 @@ In addition, the MSE and MAE (\SI{0.24}{\electronvolt}/\SI{0.25}{\electronvolt})
|
||||
|
||||
Now turning to the new results of this manuscript, \ie the alternative self-consistent scheme SRG-qs$GW$.
|
||||
Table~\ref{tab:tab1} shows the SRG-qs$GW$ values for $s=100$.
|
||||
For this value of the flow parameter, the MAE is converged to \SI{10d-3}{\electronvolt} (see Supplementary Material).
|
||||
For this value of the flow parameter, the MAE is converged to \SI{d-3}{\electronvolt} (see Supplementary Material).
|
||||
The statistical descriptors corresponding to the alternative static self-energy are all improved with respect to qs$GW$.
|
||||
Of course these are slight improvements but this is done with no additional computational cost and can be implemented really quickly just by changing the form of the static approximation.
|
||||
The evolution of statistical descriptors with respect to the various methods considered in Table~\ref{tab:tab1} is graphically illustrated by Fig.~\ref{fig:fig4}.
|
||||
The decrease of the MSE and SDE correspond to a shift of the maximum of the distribution toward 0 and a shrink of the width of the distribution, respectively.
|
||||
|
||||
In addition to this improvement of the accuracy, the SRG-qs$GW$ scheme is also much easier to converge than its qs$GW$ counterpart.
|
||||
\ANT{TO CONTINUE, waiting for s=10000}
|
||||
Indeed, up to $s=10^3$ self-consistency of the SRG-qs$GW$ scheme can be converged without any problems.
|
||||
For $s=10^4$, convergence could not be attained for the following molecules.
|
||||
On the other hand, the qs$GW$ convergence is much more erratic, the 22 molecules could be converged for $\eta=0.1$.
|
||||
However, if we decrease $\eta$ then convergence could not be attained for the whole set of molecules using the black-box convergence parameters (see Sec.~\ref{sec:comp_det}).
|
||||
Unfortunately, the convergence of the IP is not as tight as in the SRG case because for $\eta=0.01$ the IP that could be converged can vary by $10^{-2}$ to $10^{-1}$ with respect to $\eta=0.1$.
|
||||
|
||||
We will now gauge the effect of the TDA for the screening on the accuracy of the various methods considered previously.
|
||||
|
||||
Part on approximation and correction for W:
|
||||
TDHF G0W0 not that bad in GW100 but bad in GW22, qsGW TDHF even worse even with SRG,
|
||||
Maybe that would be nice to add SRG G0W0 to see if it mitigates the outliers of GW20 cf Bruneval 2021,
|
||||
That would be nice to understand clearly why qsGWTDHF is worse (screening, gap, etc)
|
||||
|
||||
\begin{table}
|
||||
\centering
|
||||
\caption{First ionization potential in eV calculated using \ANT{TO COMPLETE}. The statistical descriptors are computed for the errors with respect to the reference.}
|
||||
\caption{First ionization potential in eV calculated using $G_0W_{\text{TDA},0}$@HF, qs$GW_{\text{TDA}}$ and SRG-qs$GW_{\text{TDA}}$. The statistical descriptors are computed for the errors with respect to the reference.}
|
||||
\label{tab:tab1}
|
||||
-24.454162, -20.845919, -16.530702, -5.453359, -8.144555, -15.641055, -15.602545, -12.417807, -10.751003, -12.696334, -9.32538, -14.603221, -17.364737, -14.667361, -13.662675, -10.910541, -11.379826, -11.850654, -10.467744, -15.55351, -8.098392, -13.68235
|
||||
\begin{tabular}{lccc}
|
||||
\hline
|
||||
\hline
|
||||
@ -829,16 +834,9 @@ We will now gauge the effect of the TDA for the screening on the accuracy of the
|
||||
\end{tabular}
|
||||
\end{table}
|
||||
|
||||
Part on approximation and correction for W:
|
||||
TDHF G0W0 not that bad in GW100 but bad in GW22, qsGW TDHF even worse even with SRG,
|
||||
Maybe that would be nice to add SRG G0W0 to see if it mitigates the outliers of GW20 cf Bruneval 2021,
|
||||
That would be nice to understand clearly why qsGWTDHF is worse (screening, gap, etc)
|
||||
|
||||
Part on EA:
|
||||
MgO- does not converge yet but when we have it same analysis as Table 1 and Fig 4 but for the EA
|
||||
|
||||
Not sure that we need anything more, let's see once this is written
|
||||
|
||||
%=================================================================%
|
||||
\section{Conclusion}
|
||||
\label{sec:conclusion}
|
||||
|
Loading…
Reference in New Issue
Block a user