fix mix up in basis names
This commit is contained in:
parent
7d59570047
commit
1488fc76fe
@ -39,8 +39,8 @@
|
||||
\newcommand{\EPT}{E_{\text{PT2}}}
|
||||
\newcommand{\EDMC}{E_{\text{FN-DMC}}}
|
||||
\newcommand{\Ndet}{N_{\text{det}}}
|
||||
\newcommand{\Nelec}{N_{\text{elec}}}
|
||||
\newcommand{\Nat}{N_{\text{atoms}}}
|
||||
\newcommand{\Nelec}{N}
|
||||
\newcommand{\Nat}{M}
|
||||
\newcommand{\hartree}{$E_h$}
|
||||
|
||||
\newcommand{\LCT}{Laboratoire de Chimie Th\'eorique (UMR 7616), Sorbonne Universit\'e, CNRS, Paris, France}
|
||||
@ -300,7 +300,7 @@ In RS-DFT, the Coulomb operator entering the electron-electron repulsion is spli
|
||||
\begin{equation}
|
||||
\frac{1}{r}
|
||||
= w_{\text{ee}}^{\text{sr}, \mu}(r)
|
||||
+ w_{\text{ee}}^{\text{lr}, \mu}(r)
|
||||
+ w_{\text{ee}}^{\text{lr}, \mu}(r),
|
||||
\end{equation}
|
||||
where
|
||||
\begin{align}
|
||||
@ -428,24 +428,24 @@ Geometries for each system are reported in the {\SI}.
|
||||
|
||||
All the calculations have been performed using Burkatzki-Filippi-Dolg (BFD)
|
||||
pseudopotentials \cite{Burkatzki_2007,Burkatzki_2008} with the associated double-,
|
||||
triple-, and quadruple-$\zeta$ basis sets (BFD-VXZ).
|
||||
triple-, and quadruple-$\zeta$ basis sets (VXZ-BFD).
|
||||
The small-core BFD pseudopotentials include scalar relativistic effects.
|
||||
Coupled cluster with singles, doubles, and perturbative triples [CCSD(T)] and KS-DFT energies have been computed with
|
||||
\emph{Gaussian09},\cite{g16} using the unrestricted formalism for open-shell systems.
|
||||
|
||||
All the CIPSI calculations have been performed with \emph{Quantum
|
||||
Package}.\cite{Garniron_2019,qp2_2020} We used the short-range version
|
||||
of the local-density approximation (LDA)\cite{Sav-INC-96a,TouSavFla-IJQC-04} and Perdew-Burke-Ernzerhof (PBE) \cite{PerBurErn-PRL-96} exchange
|
||||
and correlation functionals defined in
|
||||
Package}.\cite{Garniron_2019,qp2_2020} We consider the short-range version
|
||||
of the local-density approximation (LDA)\cite{Sav-INC-96a,TouSavFla-IJQC-04} and Perdew-Burke-Ernzerhof (PBE) \cite{PerBurErn-PRL-96}
|
||||
xc functionals defined in
|
||||
Ref.~\onlinecite{GolWerStoLeiGorSav-CP-06} (see also
|
||||
Refs.~\onlinecite{TouColSav-JCP-05,GolWerSto-PCCP-05}).
|
||||
Refs.~\onlinecite{TouColSav-JCP-05,GolWerSto-PCCP-05}) that we label srLDA and srPBE respectively in the following.
|
||||
In this work, we target chemical accuracy, so
|
||||
the convergence criterion for stopping the CIPSI calculations
|
||||
has been set to $\EPT < 10^{-3}$ \hartree{} or $ \Ndet > 10^7$.
|
||||
All the wave functions are eigenfunctions of the $\Hat{S}^2$ spin operator, as
|
||||
described in Ref.~\onlinecite{Applencourt_2018}.
|
||||
|
||||
QMC calculations have been performed with QMC=Chem,\cite{Scemama_2013}
|
||||
QMC calculations have been performed with \textit{QMC=Chem},\cite{Scemama_2013}
|
||||
in the determinant localization approximation (DLA),\cite{Zen_2019}
|
||||
where only the determinantal component of the trial wave
|
||||
function is present in the expression of the wave function on which
|
||||
@ -464,12 +464,12 @@ stochastic reconfiguration algorithm developed by Assaraf \textit{et al.},
|
||||
|
||||
%%% TABLE I %%%
|
||||
\begin{table}
|
||||
\caption{Fixed-node energy $\EDMC$ (in \hartree{}) and number of determinants $\Ndet$ in \ce{H2O} for various trial wave functions $\Psi^{\mu}$ obtained with the sr-PBE density functional.}
|
||||
\caption{FN-DMC energy $\EDMC$ (in \hartree{}) and number of determinants $\Ndet$ in \ce{H2O} for various trial wave functions $\Psi^{\mu}$ obtained with the srPBE density functional.}
|
||||
\label{tab:h2o-dmc}
|
||||
\centering
|
||||
\begin{ruledtabular}
|
||||
\begin{tabular}{crlrl}
|
||||
& \multicolumn{2}{c}{BFD-VDZ} & \multicolumn{2}{c}{BFD-VTZ} \\
|
||||
& \multicolumn{2}{c}{VDZ-BFD} & \multicolumn{2}{c}{VTZ-BFD} \\
|
||||
\cline{2-3} \cline{4-5}
|
||||
$\mu$ & $\Ndet$ & $\EDMC$ & $\Ndet$ & $\EDMC$ \\
|
||||
\hline
|
||||
@ -494,7 +494,7 @@ stochastic reconfiguration algorithm developed by Assaraf \textit{et al.},
|
||||
\begin{figure}
|
||||
\centering
|
||||
\includegraphics[width=\columnwidth]{h2o-dmc.pdf}
|
||||
\caption{Fixed-node energy of \ce{H2O} as a function
|
||||
\caption{FN-DMC energy of \ce{H2O} as a function
|
||||
of $\mu$ for various levels of theory to generate
|
||||
the trial wave function.}
|
||||
\label{fig:h2o-dmc}
|
||||
@ -502,14 +502,13 @@ stochastic reconfiguration algorithm developed by Assaraf \textit{et al.},
|
||||
%%% %%% %%% %%%
|
||||
|
||||
The first question we would like to address is the quality of the
|
||||
nodes of the wave function $\Psi^{\mu}$ obtained with an intermediate
|
||||
nodes of the wave function $\Psi^{\mu}$ obtained for intermediate values of the
|
||||
range separation parameter (\ie, $0 < \mu < +\infty$).
|
||||
For this purpose, we consider a weakly correlated molecular system, namely the water
|
||||
molecule near its equilibrium geometry. \cite{Caffarel_2016}
|
||||
molecule \titou{at its experimental geometry. \cite{Caffarel_2016}}
|
||||
We then generate trial wave functions $\Psi^\mu$ for multiple values of
|
||||
$\mu$, and compute the associated fixed-node energy keeping fixed all the
|
||||
parameters such as the CI coefficients and molecular orbitals impacting the
|
||||
nodal surface.
|
||||
$\mu$, and compute the associated FN-DMC energy keeping fixed all the
|
||||
parameters impacting the nodal surface, such as the CI coefficients and the molecular orbitals.
|
||||
|
||||
%======================================================
|
||||
\subsection{Fixed-node energy of $\Psi^\mu$}
|
||||
@ -520,7 +519,7 @@ one can clearly observe that relying on FCI trial
|
||||
wave functions ($\mu = \infty$) give FN-DMC energies lower
|
||||
than the energies obtained with a single KS determinant ($\mu=0$):
|
||||
a lowering of $3.2 \pm 0.6$~m\hartree{} at the double-$\zeta$ level and $7.2 \pm
|
||||
0.3$~m\hartree{} at the triple-$\zeta$ level are obtained.
|
||||
0.3$~m\hartree{} at the triple-$\zeta$ level are obtained with the srPBE functional.
|
||||
Coming now to the nodes of the trial wave functions $\Psi^{\mu}$ with
|
||||
intermediate values of $\mu$, Fig.~\ref{fig:h2o-dmc} shows that
|
||||
a smooth behaviour is obtained:
|
||||
@ -528,7 +527,7 @@ starting from $\mu=0$ (\ie, the KS determinant),
|
||||
the FN-DMC error is reduced continuously until it reaches a minimum
|
||||
for an optimal value of $\mu$ (which is obviously basis set and functional dependent),
|
||||
and then the FN-DMC error raises until it reaches the $\mu=\infty$ limit (\ie, the FCI wave function).
|
||||
For instance, with respect to the FN-DMC/VDZ-BFD energy at $\mu=\infty$,
|
||||
For instance, with respect to the fixed-node energy associated with the srPBE/VDZ-BFD trial wave function at $\mu=\infty$,
|
||||
one can obtain a lowering of the FN-DMC energy of $2.6 \pm 0.7$~m\hartree{}
|
||||
with an optimal value of $\mu=1.75$~bohr$^{-1}$.
|
||||
This lowering in FN-DMC energy is to be compared with the $3.2 \pm
|
||||
@ -537,12 +536,12 @@ and the FCI wave function ($\mu=\infty$). When the basis set is increased, the
|
||||
gain in FN-DMC energy with respect to the FCI trial wave function is reduced,
|
||||
and the optimal value of $\mu$ is slightly shifted towards large $\mu$.
|
||||
Last but not least, the nodes of the wave functions $\Psi^\mu$ obtained with the srLDA
|
||||
exchange-correlation functional give very similar FN-DMC energies with respect
|
||||
to those obtained with the srPBE functional, even if the
|
||||
functional give very similar FN-DMC energies with respect
|
||||
to those obtained with srPBE, even if the
|
||||
RS-DFT energies obtained with these two functionals differ by several
|
||||
tens of m\hartree{}.
|
||||
|
||||
An other important aspect here regards the compactness of the trial wave functions $\Psi^\mu$:
|
||||
Another important aspect here is the compactness of the trial wave functions $\Psi^\mu$:
|
||||
at $\mu=1.75$~bohr$^{-1}$, $\Psi^{\mu}$ has \textit{only} $54\,540$ determinants at the srPBE/VDZ-BFD level, while the FCI wave function contains $200\,521$ determinants (see Table \ref{tab:h2o-dmc}). Even at the srPBE/VTZ-BFD level, we observe a reduction by a factor two in the number of determinants between the optimal $\mu$ value and $\mu = \infty$.
|
||||
The take-home message of this numerical study is that RS-DFT trial wave functions can yield a lower fixed-node energy with more compact multi-determinant expansion as compared to FCI.
|
||||
|
||||
@ -552,22 +551,27 @@ The take-home message of this numerical study is that RS-DFT trial wave function
|
||||
%======================================================
|
||||
The data presented in Sec.~\ref{sec:fndmc_mu} evidence that, in a finite basis, RS-DFT can provide
|
||||
trial wave functions with better nodes than FCI wave function.
|
||||
Such behaviour can be directly compared to the common practice of
|
||||
As mentioned in Sec.~\ref{sec:SD}, such behavior can be directly compared to the common practice of
|
||||
re-optimizing the multi-determinant part of a trial wave function $\Psi$ (the so-called Slater part) in the presence of the exponentiated Jastrow factor $e^J$. \cite{Umrigar_2005,Scemama_2006,Umrigar_2007,Toulouse_2007,Toulouse_2008}
|
||||
Hence, in the present paragraph, we would like to elaborate further on the link between RS-DFT
|
||||
and wave function optimization in the presence of a Jastrow factor.
|
||||
For simplicity in the comparison, the molecular orbitals and the Jastrow
|
||||
factor are kept fixed: only the CI coefficients are modified.
|
||||
For the sake of simplicity, the molecular orbitals and the Jastrow
|
||||
factor are kept fixed; only the CI coefficients are varied.
|
||||
|
||||
Let us assume a fixed Jastrow factor $J(\br_1, \ldots , \br_N)$,
|
||||
Let us assume a fixed Jastrow factor $J(\br_1, \ldots , \br_N)$ (where $\br_i$ is the position of the $i$th electron),
|
||||
and a corresponding Slater-Jastrow wave function $\Phi = e^J \Psi$,
|
||||
where $\Psi = \sum_I c_I D_I$ is a general linear combination of Slater determinants $D_I$.
|
||||
where
|
||||
\begin{equation}
|
||||
\label{eq:Slater}
|
||||
\Psi = \sum_I c_I D_I
|
||||
\end{equation}
|
||||
is a general linear combination of Slater determinants $D_I$.
|
||||
The only remaining variational parameters in $\Phi$ are therefore the Slater part $\Psi$.
|
||||
Let us call $\Psi^J$ the linear combination of Slater determinant minimizing the variational energy
|
||||
Let us define $\Psi^J$ as the linear combination of Slater determinants minimizing the variational energy associated with $\Phi$, \ie,
|
||||
\begin{equation}
|
||||
\Psi^J = \argmin_{\Psi}\frac{ \mel{ \Psi }{ e^{J} \hat{H} e^{J} }{ \Psi } }{\mel{ \Psi }{ e^{2J} }{ \Psi } }.
|
||||
\end{equation}
|
||||
Such a wave function $\Psi^J$ satisfies the generalized Hermitian eigenvalue equation
|
||||
Such a wave function satisfies the generalized Hermitian eigenvalue equation
|
||||
\begin{equation}
|
||||
e^{J} \hat{H} \qty( e^{J} \Psi^J ) = E e^{2J} \Psi^J,
|
||||
\label{eq:ci-j}
|
||||
@ -578,27 +582,34 @@ but also the non-Hermitian transcorrelated eigenvalue problem\cite{BoyHan-PRSLA-
|
||||
e^{-J} \hat{H} \qty( e^{J} \Psi^J) = E \Psi^J,
|
||||
\end{equation}
|
||||
which is much easier to handle despite its non-Hermiticity.
|
||||
Of course, the FN-DMC energy of $\Phi$ depends therefore only on the nodes of $\Psi^J$.
|
||||
Of course, the FN-DMC energy of $\Phi$ depends only on the nodes of $\Psi^J$.
|
||||
In a finite basis set and with a quite accurate Jastrow factor, it is known that the nodes
|
||||
of $\Psi^J$ may be better than that of the FCI wave function, and therefore, we would like to compare $\Psi^J$ and $\Psi^\mu$.
|
||||
of $\Psi^J$ may be better than the nodes of the FCI wave function.
|
||||
Hence, we would like to compare $\Psi^J$ and $\Psi^\mu$.
|
||||
|
||||
To do so, we have made the following numerical experiment.
|
||||
First, we extract the 200 determinants with the largest weights in the FCI wave
|
||||
function out of a large CIPSI calculation. Within this set of determinants,
|
||||
function out of a large CIPSI calculation obtained with the VDZ-BFD basis. Within this set of determinants,
|
||||
we solve the self-consistent equations of RS-DFT [see Eq.~\eqref{rs-dft-eigen-equation}]
|
||||
with different values of $\mu$. This gives the CI expansions $\Psi^\mu$.
|
||||
Then, within the same set of determinants we optimize the CI coefficients in the presence of
|
||||
a simple one- and two-body Jastrow factor $e^J$ of the form $\exp(J_{eN} + J_{ee})$ with
|
||||
\begin{eqnarray}
|
||||
J_\text{eN} & = & - \sum_{A=1}^{\Nat} \sum_{i=1}^{\Nelec} \left( \frac{\alpha_A\, r_{iA}}{1 + \alpha_A\, r_{iA}} \right)^2
|
||||
Then, within the same set of determinants we optimize the CI coefficients $c_I$ [see Eq.~\eqref{eq:Slater}] in the presence of
|
||||
a simple one- and two-body Jastrow factor $e^J$ with $J = J_\text{eN} + J_\text{ee}$ and
|
||||
\begin{subequations}
|
||||
\begin{gather}
|
||||
J_\text{eN} = - \sum_{A=1}^{\Nat} \sum_{i=1}^{\Nelec} \qty( \frac{\alpha_A\, r_{iA}}{1 + \alpha_A\, r_{iA}} )^2,
|
||||
\label{eq:jast-eN} \\
|
||||
J_\text{ee} & = & \sum_{i=1}^{\Nelec} \sum_{j=1}^{i-1} \frac{a\, r_{ij}}{1 + b\, r_{ij}}. \label{eq:jast-ee}
|
||||
\end{eqnarray}
|
||||
$J_\text{eN}$ contains the electron-nucleus terms with a single parameter
|
||||
$\alpha_A$ per atom, and $J_\text{ee}$ contains the electron-electron terms
|
||||
where the indices $i$ and $j$ loop over all electron pairs. The parameters $a=1/2$
|
||||
J_\text{ee} = \sum_{i < j}^{\Nelec} \frac{a\, r_{ij}}{1 + b\, r_{ij}}.
|
||||
\label{eq:jast-ee}
|
||||
\end{gather}
|
||||
\end{subequations}
|
||||
The one-body Jastrow factor $J_\text{eN}$ contains the electron-nucleus terms (where $\Nat$ is the number of nuclei) with a single parameter
|
||||
$\alpha_A$ per nucleus.
|
||||
The two-body Jastrow factor $J_\text{ee}$ contains the electron-electron terms
|
||||
where the sum over $i < j$ loops over all electron pairs.
|
||||
In Eqs.~\eqref{eq:jast-eN} and \eqref{eq:jast-ee}, $r_{iA}$ is the distance between the $i$th electron and the $A$th nucleus while $r_{ij}$ is the interlectronic distance between electrons $i$ and $j$.
|
||||
The parameters $a=1/2$
|
||||
and $b=0.89$ were fixed, and the parameters $\gamma_{\text{O}}=1.15$ and $\gamma_{\text{H}}=0.35$
|
||||
were obtained by energy minimization with a single determinant.
|
||||
were obtained by energy minimization with a single \titou{HF?} determinant.
|
||||
The optimal CI expansion $\Psi^J$ is obtained by sampling the matrix elements
|
||||
of the Hamiltonian ($\mathbf{H}$) and the overlap ($\mathbf{S}$), in the
|
||||
basis of Jastrow-correlated determinants $e^J D_i$:
|
||||
@ -611,7 +622,7 @@ S_{ij} = \expval{ \frac{e^J D_i}{\Psi^J}\, \frac{e^J D_j}{\Psi^J} }
|
||||
and solving Eq.~\eqref{eq:ci-j}.\cite{Nightingale_2001}
|
||||
|
||||
We can easily compare $\Psi^\mu$ and $\Psi^J$ as they are developed
|
||||
on the same Slater determinant basis.
|
||||
on the same set of Slater determinants.
|
||||
In Fig.~\ref{fig:overlap}, we plot the overlaps
|
||||
$\braket*{\Psi^J}{\Psi^\mu}$ obtained for water,
|
||||
and in Fig.~\ref{dmc_small} the FN-DMC energy of the wave functions
|
||||
@ -621,10 +632,9 @@ $\Psi^\mu$ together with that of $\Psi^J$.
|
||||
\begin{figure}
|
||||
\centering
|
||||
\includegraphics[width=\columnwidth]{overlap.pdf}
|
||||
\caption{\ce{H2O}, double-zeta basis set, 200 most important
|
||||
determinants of the FCI expansion (see Sec.~\ref{sec:rsdft-j}).
|
||||
Overlap of the RS-DFT CI expansions $\Psi^\mu$ with the CI
|
||||
expansion optimized in the presence of a Jastrow factor $\Psi^J$.}
|
||||
\caption{Overlap between $\Psi^\mu$ and $\Psi^J$ as a function of $\mu$ for \ce{H2O}.
|
||||
For these two trial wave functions, the CI expansion consists of the 200 most important
|
||||
determinants of the FCI expansion obtained with the VDZ-BFD basis (see Sec.~\ref{sec:rsdft-j} for more details).}
|
||||
\label{fig:overlap}
|
||||
\end{figure}
|
||||
%%% %%% %%% %%%
|
||||
@ -916,7 +926,7 @@ impacted by this spurious effect, as opposed to FCI.
|
||||
|
||||
%%% TABLE IV %%%
|
||||
\begin{table}
|
||||
\caption{FN-DMC energies (in hartree) of the triplet carbon atom (BFD-VDZ) with
|
||||
\caption{FN-DMC energies (in hartree) of the triplet carbon atom (VDZ-BFD) with
|
||||
different values of $m_s$.}
|
||||
\label{tab:spin}
|
||||
\begin{ruledtabular}
|
||||
@ -1124,7 +1134,7 @@ while keeping the size-consistency.
|
||||
|
||||
Note that when $\mu=0$ the number of determinants is not equal to one because
|
||||
we have used the natural orbitals of a first CIPSI calculation, and
|
||||
not the sr-PBE orbitals.
|
||||
not the srPBE orbitals.
|
||||
So the Kohn-Sham determinant is expressed as a linear combination of
|
||||
determinants built with natural orbitals. It is possible to add
|
||||
an extra step to the algorithm to compute the natural orbitals from the
|
||||
|
Loading…
Reference in New Issue
Block a user