Merge branch 'master' of git.irsamc.ups-tlse.fr:scemama/RSDFT-CIPSI-QMC
This commit is contained in:
commit
56095fff4b
@ -1,13 +1,48 @@
|
||||
%% This BibTeX bibliography file was created using BibDesk.
|
||||
%% http://bibdesk.sourceforge.net/
|
||||
|
||||
%% Created for Pierre-Francois Loos at 2020-08-07 13:35:24 +0200
|
||||
%% Created for Pierre-Francois Loos at 2020-08-07 20:17:07 +0200
|
||||
|
||||
|
||||
%% Saved with string encoding Unicode (UTF-8)
|
||||
|
||||
|
||||
|
||||
@article{Assaraf_2000,
|
||||
Author = {Assaraf, Roland and Caffarel, Michel and Khelif, Anatole},
|
||||
Date-Added = {2020-08-07 20:12:45 +0200},
|
||||
Date-Modified = {2020-08-07 20:12:45 +0200},
|
||||
Doi = {10.1103/physreve.61.4566},
|
||||
Issn = {1095-3787},
|
||||
Journal = {Phys. Rev. E},
|
||||
Month = {Apr},
|
||||
Number = {4},
|
||||
Pages = {4566--4575},
|
||||
Publisher = {American Physical Society (APS)},
|
||||
Title = {Diffusion Monte Carlo methods with a fixed number of walkers},
|
||||
Url = {http://dx.doi.org/10.1103/PhysRevE.61.4566},
|
||||
Volume = {61},
|
||||
Year = {2000},
|
||||
Bdsk-Url-1 = {http://dx.doi.org/10.1103/PhysRevE.61.4566},
|
||||
Bdsk-Url-2 = {http://dx.doi.org/10.1103/physreve.61.4566}}
|
||||
|
||||
@article{Assaraf_2007,
|
||||
Author = {Assaraf, Roland and Caffarel, Michel and Scemama, Anthony},
|
||||
Date-Added = {2020-08-07 20:12:45 +0200},
|
||||
Date-Modified = {2020-08-07 20:12:45 +0200},
|
||||
Doi = {10.1103/physreve.75.035701},
|
||||
Issn = {1550-2376},
|
||||
Journal = {Phys. Rev. E},
|
||||
Month = {Mar},
|
||||
Number = {3},
|
||||
Publisher = {American Physical Society (APS)},
|
||||
Title = {Improved Monte Carlo estimators for the one-body density},
|
||||
Url = {http://dx.doi.org/10.1103/PhysRevE.75.035701},
|
||||
Volume = {75},
|
||||
Year = {2007},
|
||||
Bdsk-Url-1 = {http://dx.doi.org/10.1103/PhysRevE.75.035701},
|
||||
Bdsk-Url-2 = {http://dx.doi.org/10.1103/physreve.75.035701}}
|
||||
|
||||
@article{Burkatzki_2007,
|
||||
Author = {Burkatzki, M. and Filippi, C. and Dolg, M.},
|
||||
Date-Added = {2020-08-07 13:35:15 +0200},
|
||||
|
@ -65,7 +65,7 @@
|
||||
By combining density-functional theory (DFT) and wave function theory (WFT) via the range separation (RS) of the interelectronic Coulomb operator, we obtain accurate fixed-node diffusion Monte Carlo (FN-DMC) energies with compact multideterminant trial wave functions.
|
||||
These compact trial wave functions are generated via the diagonalization of the RS-DFT Hamiltonian.
|
||||
In particular, we combine here short-range correlation functionals with selected configuration interaction (SCI).
|
||||
As the WFT method is relieved from describing the short-range part of the correlation hole around the electron?electron coalescence points, the number of determinants in the trial wave function required to reach a given accuracy is significantly reduced as compared to a conventional SCI calculation.
|
||||
As the WFT method is relieved from describing the short-range part of the correlation hole around the electron-electron coalescence points, the number of determinants in the trial wave function required to reach a given accuracy is significantly reduced as compared to a conventional SCI calculation.
|
||||
\titou{T2: work in progress.}
|
||||
\end{abstract}
|
||||
|
||||
@ -400,9 +400,10 @@ the pseudopotential is localized. Hence, in the DLA the fixed-node
|
||||
energy is independent of the Jastrow factor, as in all-electron
|
||||
calculations. Simple Jastrow factors were used to reduce the
|
||||
fluctuations of the local energy.
|
||||
\titou{time step $5 \times 10^{-4}$.
|
||||
All-electron move DMC.}
|
||||
|
||||
The FN-DMC simulations are performed with the stochastic reconfiguration
|
||||
algorithm developed by Assaraf \textit{et al.}, \cite{Assaraf_2000}
|
||||
with a time step of $5 \times 10^{-4}$ a.u.
|
||||
\titou{All-electron move DMC?}
|
||||
|
||||
|
||||
\section{Influence of the range-separation parameter on the fixed-node
|
||||
@ -497,7 +498,7 @@ where $\Psi = \sum_I c_I D_I$ is a general linear combination of Slater determin
|
||||
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
|
||||
\begin{equation}
|
||||
\Psi^J = \text{argmin}_{\Psi}\frac{ \langle \Psi | e^{J} H e^{J} |\Psi \rangle}{\langle \Psi | e^{2J} |\Psi \rangle}.
|
||||
\Psi^J = \text{argmin}_{\Psi}\frac{ \mel{ \Psi }{ e^{J} H e^{J} }{ \Psi } }{\mel{ \Psi }{ e^{2J} }{ \Psi } }.
|
||||
\end{equation}
|
||||
Such a wave function $\Psi^J$ satisfies the generalized hermitian eigenvalue equation
|
||||
\begin{equation}
|
||||
@ -558,16 +559,16 @@ introducing short-range correlation with DFT has
|
||||
an impact on the CI coefficients similar to the Jastrow factor.
|
||||
In the case of F$_2$, the Jastrow factor has
|
||||
very little effect on the CI coefficients, as the overlap
|
||||
$\braket{\Psi^J}{\Psi^{\mu=\infty}}$ is very close to
|
||||
$\braket*{\Psi^J}{\Psi^{\mu=\infty}}$ is very close to
|
||||
$1$.
|
||||
Nevertheless, a slight maximum is obtained for
|
||||
$\mu=5$~bohr$^{-1}$.
|
||||
In order to refine the comparison between $\Psi^\mu$ and $\Psi^J$,
|
||||
we report, in the case of the water molecule in the double-zeta basis set,
|
||||
several quantities related to the one- and two-body density of $\Psi^J$ and $\Psi^\mu$ with different values of $\mu$.
|
||||
First, we report in table~\ref{table_on_top} the integrated on-top pair density $\langle n_2({\bf r},{\bf r}) \rangle$
|
||||
First, we report in table~\ref{table_on_top} the integrated on-top pair density $\expval{ n_2({\bf r},{\bf r}) }$
|
||||
\begin{equation}
|
||||
\langle n_2({\bf r},{\bf r}) \rangle = \int \text{d}{\bf r} \,\,n_2({\bf r},{\bf r})
|
||||
\expval{ n_2({\bf r},{\bf r}) } = \int \text{d}{\bf r} \,\,n_2({\bf r},{\bf r})
|
||||
\end{equation}
|
||||
where $n_2({\bf r}_1,{\bf r}_2)$ is the two-body density (normalized to $N(N-1)$ where $N$ is the number of electrons)
|
||||
obtained for both $\Psi^\mu$ and $\Psi^J$.
|
||||
@ -601,12 +602,12 @@ Therefore, one can understand the similarity between the eigenfunctions of $H^\m
|
||||
they both deal with an effective non-divergent interaction but still produce reasonable one-body density.
|
||||
|
||||
\begin{table}
|
||||
\caption{H$_2$O, double-zeta basis set. Integrated on-top pair density $\langle n_2({\bf r},{\bf r}) \rangle$
|
||||
\caption{H$_2$O, double-zeta basis set. Integrated on-top pair density $\expval{ n_2({\bf r},{\bf r}) }$
|
||||
for $\Psi^J$ and $\Psi^\mu$ with different values of $\mu$. }
|
||||
\label{table_on_top}
|
||||
\begin{ruledtabular}
|
||||
\begin{tabular}{cc}
|
||||
$\mu$ & $\langle n_2({\bf r},{\bf r}) \rangle$ \\
|
||||
$\mu$ & $\expval{ n_2({\bf r},{\bf r}) }$ \\
|
||||
\hline
|
||||
0.00 & 1.443 \\
|
||||
0.25 & 1.438 \\
|
||||
@ -624,7 +625,7 @@ they both deal with an effective non-divergent interaction but still produce rea
|
||||
\begin{figure}
|
||||
\centering
|
||||
\includegraphics[width=\columnwidth]{on-top-mu.pdf}
|
||||
\caption{H$_2$O, double-zeta basis set. Integrated on-top pair density $\langle n_2({\bf r},{\bf r}) \rangle$ along the OH axis, for $\Psi^J$ and $\Psi^\mu$ with different values of $\mu$. }
|
||||
\caption{H$_2$O, double-zeta basis set. Integrated on-top pair density $\expval{ n_2({\bf r},{\bf r}) }$ along the OH axis, for $\Psi^J$ and $\Psi^\mu$ with different values of $\mu$. }
|
||||
\label{fig:n2}
|
||||
\end{figure}
|
||||
\begin{figure}
|
||||
|
Loading…
Reference in New Issue
Block a user