loos/srDFT_SC
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

theory

Browse Source
This commit is contained in:
Pierre-Francois Loos 2020-01-23 10:51:07 +01:00
parent ea74d0e2d5
commit 5427cab71e
1 changed files with 10 additions and 20 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

30
Manuscript/srDFT_SC.tex
View File

@ -314,7 +314,6 @@ In practice, WFT uses a finite one-electron basis set. The exact solution of the
Beside the difficulties of accurately describing the molecular electronic structure within a given basis set, a crucial limitation of WFT methods is the slow convergence of the energy (and related properties) with respect to the size of the basis set. As initially shown by the seminal work of Hylleraas \cite{Hyl-ZP-29} and further developed by Kutzelnigg and coworkers, \cite{Kut-TCA-85,KutKlo-JCP-91, NogKut-JCP-94} the main convergence problem originates from the divergence of the electron-electron Coulomb interaction at the coalescence point, which induces a discontinuity in the first derivative of the exact wave function (the so-called electron-electron cusp). Describing such a discontinuity with an incomplete one-electron basis set is impossible and, as a consequence, the convergence of the computed energies and properties are strongly affected. To alleviate this problem, extrapolation techniques have been developed, either based on a partial-wave expansion analysis, \cite{HelKloKocNog-JCP-97,HalHelJorKloKocOlsWil-CPL-98} or more recently based on perturbative arguments. \cite{IrmHulGru-PRL-19,IrmGru-JCP-2019} A more rigorous approach to tackle the basis-set convergence problem is provided by the so-called explicitly correlated F12 (or R12) methods \cite{Ten-TCA-12,TenNog-WIREs-12,HatKloKohTew-CR-12, KonBisVal-CR-12, GruHirOhnTen-JCP-17, MaWer-WIREs-18} which introduce a geminal function depending explicitly on the interelectronic distance. This ensures a correct representation of the Coulomb correlation hole around the electron-electron coalescence point, and leads to a much faster convergence of the energy than usual WFT methods. For instance, using the explicitly correlated version of coupled cluster with singles, doubles, and perturbative triples [CCSD(T)] in a triple-$\zeta$ basis set is equivalent to using a quintuple-$\zeta$ basis set with the usual CCSD(T) method, \cite{TewKloNeiHat-PCCP-07} although a computational overhead is introduced by the auxiliary basis set needed to compute the three- and four-electron integrals involved in F12 theory. \cite{BarLoo-JCP-17} In addition to the computational cost, a possible drawback of F12 theory is its rather complex formalism which requires non-trivial developments for adapting it to a new method. For strongly correlated systems, several multi-reference methods have been extended to explicit correlation (see, for example, Refs.~\onlinecite{Ten-CPL-07,ShiWer-JCP-10,TorKniWer-JCP-11,DemStanMatTenPitNog-PCCP-12,GuoSivValNee-JCP-17}), including approaches based on the so-called universal F12 theory which are potentially applicable to any electronic-structure computational methods. \cite{TorVal-JCP-09,KonVal-JCP-11,HauMaoMukKlo-CPL-12,BooCleAlaTew-JCP-12}
An alternative way to improve the convergence towards the complete basis set (CBS) limit is to treat the short-range correlation effects within DFT and to use WFT methods to deal only with the long-range and/or strong correlation effects. A rigorous approach achieving this mixing of DFT and WFT is range-separated DFT (RSDFT) (see Ref.~\onlinecite{TouColSav-PRA-04} and references therein) which relies on a decomposition of the electron-electron Coulomb interaction in terms of the interelectronic distance thanks to a range-separation parameter $\mu$. The advantage of this approach is at least two-fold: i) the DFT part deals primarily with the short-range part of the Coulomb interaction, and consequently the usual semilocal density-functional approximations are more accurate than for standard KS DFT; ii) the WFT part deals only with a smooth non-divergent interaction, and consequently the wave function has no electron-electron cusp \cite{GorSav-PRA-06} and the basis-set convergence is much faster. \cite{FraMusLupTou-JCP-15} A number of approximate RSDFT schemes have been developed involving single-reference \cite{AngGerSavTou-PRA-05, GolWerSto-PCCP-05, TouGerJanSavAng-PRL-09,JanHenScu-JCP-09, TouZhuSavJanAng-JCP-11, MusReiAngTou-JCP-15,KalTou-JCP-18,KalMusTou-JCP-19} and multi-reference \cite{LeiStoWerSav-CPL-97, FroTouJen-JCP-07, FroCimJen-PRA-10, HedKneKieJenRei-JCP-15, HedTouJen-JCP-18, FerGinTou-JCP-18} WFT methods. Nevertheless, there are still some open issues in RSDFT, such as remaining fractional-charge and fractional-spin errors in the short-range density functionals \cite{MusTou-MP-17} or the dependence of the quality of the results on the value of the range-separation parameter $\mu$.
% which can be seen as an empirical parameter.
Building on the development of RSDFT, a possible solution to the basis-set convergence problem has been recently proposed by some of the present authors~\cite{GinPraFerAssSavTou-JCP-18} in which RSDFT functionals are used to recover only the correlation effects outside a given basis set. The key point here is to realize that a wave function developed in an incomplete basis set is cuspless and could also originate from a Hamiltonian with a non-divergent long-range electron-electron interaction. Therefore, a mapping with RSDFT can be performed through the introduction of an effective non-divergent interaction representing the usual electron-electron Coulomb interaction projected in an incomplete basis set. First applications to weakly correlated molecular systems have been successfully carried out, \cite{LooPraSceTouGin-JCPL-19} together with extensions of this approach to the calculations of excitation energies \cite{GinSceTouLoo-JCP-19} and ionization potentials. \cite{LooPraSceGinTou-JCTC-20} The goal of the present work is to further develop this approach for the description of strongly correlated systems.
@ -387,11 +386,11 @@ As originally derived in Ref.~\onlinecite{GinPraFerAssSavTou-JCP-18} (see Sec.~I
\infty, & \text{otherwise,}
\end{cases}
\end{equation}
where $\twodmrdiagpsi$ is the opposite-spin pair density associated with $\wf{}{\Bas}$
where
\begin{equation}
\twodmrdiagpsi = \sum_{pqrs \in \Bas} \SO{p}{1} \SO{q}{2} \Gam{pq}{rs} \SO{r}{1} \SO{s}{2},
\end{equation}
and $\Gam{pq}{rs} = 2 \mel*{\wf{}{\Bas}}{ \aic{r_\downarrow}\aic{s_\uparrow}\ai{q_\uparrow}\ai{p_\downarrow}}{\wf{}{\Bas}}$ its associated tensor in a basis of spatial orthonormal orbitals $\{\SO{p}{}\}$,
is the opposite-spin pair density associated with $\wf{}{\Bas}$, and $\Gam{pq}{rs} = 2 \mel*{\wf{}{\Bas}}{ \aic{r_\downarrow}\aic{s_\uparrow}\ai{q_\uparrow}\ai{p_\downarrow}}{\wf{}{\Bas}}$ its associated tensor in a basis of spatial orthonormal orbitals $\{\SO{p}{}\}$,
\begin{equation}
\label{eq:fbasis}
\fbasis
@ -471,7 +470,7 @@ It is also noteworthy that, with the present definition, $\wbasisval$ still tend
\subsubsection{Generic approximate form}
\label{sec:functional_form}
As originally proposed and motivated in Ref.~\onlinecite{GinPraFerAssSavTou-JCP-18}, we approximate the complementary functional $\efuncden{\den}$ by using the so-called correlation energy functional with multideterminant reference (ECMD) introduced by Toulouse \textit{et al.}\cite{TouGorSav-TCA-05} Following the recent work in Ref.~\onlinecite{LooPraSceTouGin-JCPL-19}, we propose to use a Perdew-Burke-Ernzerhof (PBE)-like functional which uses the one-electron density $\denr$, the spin polarization $\zeta(\br{})=[n_\uparrow(\br{})-n_\downarrow(\br{})]/\denr$ (where $n_\uparrow(\br{})$ and $n_\downarrow(\br{})$ are the spin-up and spin-down densities), the reduced density gradient $s(\br{}) = \nabla \denr/\denr^{4/3}$, and the on-top pair density $\ntwo(\br{})\equiv \ntwo(\br{},\br{})$. In the present work, all these quantities are computed with the same wave function $\psibasis$ used to define $\mur \equiv\murpsi$.
As originally proposed and motivated in Ref.~\onlinecite{GinPraFerAssSavTou-JCP-18}, we approximate the complementary functional $\efuncden{\den}$ by using the so-called correlation energy functional with multideterminant reference (ECMD) introduced by Toulouse \textit{et al.}\cite{TouGorSav-TCA-05} Following the recent work in Ref.~\onlinecite{LooPraSceTouGin-JCPL-19}, we propose to consider a Perdew-Burke-Ernzerhof (PBE)-like functional which uses the one-electron density $\denr$, the spin polarization $\zeta(\br{})=[n_\uparrow(\br{})-n_\downarrow(\br{})]/\denr$ (where $n_\uparrow(\br{})$ and $n_\downarrow(\br{})$ are the spin-up and spin-down densities), the reduced density gradient $s(\br{}) = \nabla \denr/\denr^{4/3}$, and the on-top pair density $\ntwo(\br{})\equiv \ntwo(\br{},\br{})$. In the present work, all these quantities are computed with the same wave function $\psibasis$ used to define $\mur \equiv\murpsi$.
Therefore, $\efuncden{\den}$ has the following generic form
\begin{multline}
\label{eq:def_ecmdpbebasis}
@ -482,7 +481,7 @@ Therefore, $\efuncden{\den}$ has the following generic form
where
\begin{equation}
\label{eq:def_ecmdpbe}
\ecmd(\argecmd) = \frac{\varepsilon_{\text{c}}^{\text{PBE}}(\argepbe)}{1+ \beta(\argepbe,\ntwo) \; \mu^3},
\ecmd(\argecmd) = \frac{\varepsilon_{\text{c}}^{\text{PBE}}(\argepbe)}{1+ \beta(\argepbe,\ntwo) \; \mu^3}
\end{equation}
is the correlation energy per particle, with
\begin{equation}
@ -508,7 +507,6 @@ Note also that $\ecmd(\argecmd)$ vanishes when $\ntwo$ vanishes, \ie,
\lim_{\ntwo \to 0} \ecmd(\argecmd) = 0,
\end{equation}
which is expected for systems with a vanishing on-top pair density.
%, such as the totally dissociated H$_2$ molecule which is the archetype of strongly correlated systems.
Finally, the function $\ecmd(\argecmd)$ vanishes when $\mu \to \infty$ like all RSDFT short-range functionals, \ie,
\begin{equation}
\label{eq:lim_muinf}
@ -531,7 +529,7 @@ Second, $\efuncdenpbe{\argebasis}$ correctly vanishes for systems with uniformly
\label{eq:lim_ebasis}
\lim_{n_2 \to 0} \efuncdenpbe{\argebasis} = 0.
\end{equation}
This property is doubly guaranteed by i) the choice of setting the effective interaction $\wbasis$ at $\infty$ for a vanishing pair density [see Eq.~\eqref{eq:wbasis}] leading to $\mu(\br{}) \to \infty$ and thus a vanishing $\ecmd(\argecmd)$ according to Eq.~\eqref{eq:lim_muinf}, and ii) the fact that $\ecmd(\argecmd)$ vanishes anyway when the on-top pair density vanishes [see Eq.~\eqref{eq:lim_n2}].
This property is doubly guaranteed by i) the choice of setting $\wbasis = +\infty$ for a vanishing pair density [see Eq.~\eqref{eq:wbasis}], which leads to $\mu(\br{}) \to \infty$ and thus a vanishing $\ecmd(\argecmd)$ [see Eq.~\eqref{eq:lim_muinf}], and ii) the fact that $\ecmd(\argecmd)$ vanishes anyway when the on-top pair density vanishes [see Eq.~\eqref{eq:lim_n2}].
\subsection{Requirements on the complementary functional for strong correlation}
\label{sec:requirements}
@ -547,13 +545,9 @@ A sufficient condition to achieve spin-multiplet degeneracy is to eliminate all
To do so, it has been proposed to replace the dependence on the spin polarization by the dependence on the on-top pair density. Most often, it is done by introducing an effective spin polarization~\cite{MosSan-PRA-91,BecSavSto-TCA-95,Sav-INC-96a,Sav-INC-96,MieStoSav-MP-97,TakYamYam-CPL-02,TakYamYam-IJQC-04,GraCre-MP-05,TsuScuSav-JCP-10,LimCarLuoMaOlsTruGag-JCTC-14,GarBulHenScu-JCP-15,GarBulHenScu-PCCP-15,CarTruGag-JCTC-15,GagTruLiCarHoyBa-ACR-17} (see, also, Refs.~\onlinecite{PerSavBur-PRA-95,StaDav-CPL-01})
\begin{equation}
\label{eq:def_effspin}
\tilde{\zeta}(n,n_{2}) =
% \begin{cases}
\sqrt{ 1 - 2 \; n_{2}/n^2 },
% 0 & \text{otherwise.}
% \end{cases}
\tilde{\zeta}(n,n_{2}) = \sqrt{ 1 - 2 \; n_{2}/n^2 }
\end{equation}
expressed as a function of the density $n$ and the on-top pair density $n_2$, calculated from a given wave function. The advantage of this approach is that this effective spin polarization $\tilde{\zeta}$ is independent from $S_z$, since the on-top pair density is $S_z$-independent. Nevertheless, the use of $\tilde{\zeta}$ in Eq.~\eqref{eq:def_effspin} presents some disadvantages since this expression was derived for a single-determinant wave function. Hence, it does not appear justified to use it for a multideterminant wave function. More particularly, it may happen, in the multideterminant case, that $1 - 2 \; n_{2}/n^2 < 0 $ which results in a complex-valued effective spin polarization. \cite{BecSavSto-TCA-95}
expressed as a function of the density $n$ and the on-top pair density $n_2$ calculated from a given wave function. The advantage of this approach is that this effective spin polarization $\tilde{\zeta}$ is independent from $S_z$ since the on-top pair density is $S_z$-independent. Nevertheless, the use of $\tilde{\zeta}$ in Eq.~\eqref{eq:def_effspin} presents some disadvantages since this expression was derived for a single-determinant wave function. Hence, it does not appear justified to use it for a multideterminant wave function. More particularly, it may happen, in the multideterminant case, that $1 - 2 \; n_{2}/n^2 < 0 $ which results in a complex-valued effective spin polarization. \cite{BecSavSto-TCA-95}
Therefore, following other authors, \cite{MieStoSav-MP-97,LimCarLuoMaOlsTruGag-JCTC-14,GarBulHenScu-JCP-15} we use the following definition
\begin{equation}
\label{eq:def_effspin-0}
@ -570,13 +564,12 @@ An alternative way to eliminate the $S_z$ dependence is to simply set $\zeta=0$,
\subsubsection{Size consistency}
Since $\efuncdenpbe{\argebasis}$ is computed via a single integral over $\mathbb{R}^3$ [see Eq.~\eqref{eq:def_ecmdpbebasis}] which involves only local quantities [$n(\br{})$, $\zeta(\br{})$, $s(\br{})$, $n_2(\br{})$, and $\mu(\br{})$], in the case of non-overlapping fragments $\text{A}+\text{B}$, it can be written as the sum of two local contributions: one coming from the integration over the region of subsystem \ce{A} and the other one from the region of subsystem \ce{B}. Therefore, a sufficient condition for size consistency is that these local quantities coincide in the isolated systems and in the subsystems of the supersystem $\text{A}+\text{B}$. Since these local quantities are calculated from the wave function $\psibasis$, a sufficient condition is that the wave function is multiplicatively separable in the limit of non-interacting fragments, \ie, $\ket{\Psi_{\text{A}+\text{B}}^{\basis}} = \ket{\Psi_{\ce{A}}^{\basis}} \otimes \ket{\Psi_{\ce{B}}^{\basis}}$. We refer the interested reader to Appendix~\ref{app:sizeconsistency} for a detailed proof and discussion of the latter statement.
Since $\efuncdenpbe{\argebasis}$ is computed via a single integral over $\mathbb{R}^3$ [see Eq.~\eqref{eq:def_ecmdpbebasis}] which involves only local quantities [$n(\br{})$, $\zeta(\br{})$, $s(\br{})$, $n_2(\br{})$, and $\mu(\br{})$], in the case of non-overlapping fragments $\text{A}+\text{B}$, it can be written as the sum of two local contributions: one coming from the integration over the region of subsystem \ce{A} and the other one from the region of subsystem \ce{B}. Therefore, a sufficient condition for size consistency is that these local quantities coincide in the isolated systems and in the subsystems of the supersystem $\text{A}+\text{B}$. Since these local quantities are calculated from the wave function $\psibasis$, a sufficient condition is that the wave function is multiplicatively separable in the limit of non-interacting fragments, \ie, $\ket*{\Psi_{\text{A}+\text{B}}^{\basis}} = \ket*{\Psi_{\ce{A}}^{\basis}} \otimes \ket*{\Psi_{\ce{B}}^{\basis}}$. We refer the interested reader to Appendix~\ref{app:sizeconsistency} for a detailed proof and discussion of the latter statement.
In the case where the two subsystems \ce{A} and \ce{B} dissociate in closed-shell systems, a simple RHF wave function ensures this property, but when one or several covalent bonds are broken, a properly chosen CASSCF wave function is sufficient to recover this property. The underlying active space must however be chosen in such a way that it leads to size-consistent energies in the limit of dissociated fragments.
\subsection{Actual approximations used for the complementary functional}
\label{sec:def_func}
%\subsubsection{Definition of the protocol to design functionals}
As the present work focuses on the strong-correlation regime, we propose here to investigate only approximate functionals which are $S_z$ independent and size-consistent in the case of covalent bond breaking. Therefore, the wave functions $\psibasis$ used throughout this paper are CASSCF wave functions in order to ensure size consistency of all local quantities. The difference between two flavors of functionals are only due to the type of i) spin polarization, and ii) on-top pair density.
@ -622,11 +615,8 @@ We then define three complementary functionals:
\bar{E}^\Bas_{\pbeontns} = \int \d\br{} \,\denr \ecmd(\argrpbeontns).
\end{equation}
\end{itemize}
The performance of each of these functionals is tested in the following. Notice that we did not define a spin-unpolarized PBE-UEG functional which would be significantly inferior to the three other functionals. Indeed, without knowledge of the spin polarization or of the accurate on-top pair density, such a functional would be inaccurate.
The performance of each of these functionals is tested in the following. Note that we did not define a spin-unpolarized version of the PBE-UEG functional because it would have been significantly inferior (in terms of performance) compared to the three other functionals. Indeed, because to the lack of knowledge on the spin polarization or on the accurate on-top pair density, such a functional would be inaccurate. This assumption has been numerically confirmed by preliminary calculations.
%DFT: BLACK BOX and not CASSCF
%%%%%%%%%%%%%%%%%%%%%%%%
\section{Results}
\label{sec:results}
@ -732,7 +722,7 @@ Regarding the computational cost of the present approach, it should be stressed
The study of the \ce{H10} chain with equally distant atoms is a good prototype of strongly-correlated systems as it consists in the simultaneous breaking of 10 covalent $\sigma$ bonds which all interact with each other. Also, being a relatively small system, benchmark calculations at near-CBS values can be obtained (see Ref.~\onlinecite{h10_prx} for a detailed study of this system).
We report in Fig.~\ref{fig:H10} the potential energy curves computed using the cc-pVXZ (X $=$ D, T, and Q) basis sets for different levels of approximation. The computation of the atomization energies $D_0$ for each level of theory is reported in Table \ref{tab:d0}. A general trend that can be observed from these data is that, in a given basis set, the quality of the potential energy curves are globally improved by adding the basis-set correction, independently of the approximation level of $\efuncbasis$. Also, no erratic behavior is found when stretching the bonds, which shows that the present procedure (\textit{i.e.} the determination of the range-separation function and the definition of the functionals) is robust when reaching the strong-correlation regime.
We report in Fig.~\ref{fig:H10} the potential energy curves computed using the cc-pVXZ (X $=$ D, T, and Q) basis sets for different levels of approximation. The computation of the atomization energies $D_0$ for each level of theory is reported in Table \ref{tab:d0}. A general trend that can be observed from these data is that, in a given basis set, the quality of the potential energy curves are globally improved by adding the basis-set correction, independently of the approximation level of $\efuncbasis$. Also, no erratic behavior is found when stretching the bonds, which shows that the present procedure (\ie, the determination of the range-separation function and the definition of the functionals) is robust when reaching the strong-correlation regime.
In other words, smooth potential energy curves are obtained with the present basis-set correction.
More quantitatively, the values of $D_0$ are within chemical accuracy (\ie, an error below $1.4$ mHa) from the cc-pVTZ basis set when using the $\pbeontXi$ and $\pbeontns$ functionals, whereas such an accuracy is not even reached at the standard MRCI+Q/cc-pVQZ level of theory.