From f9294732b4dce6c3264cb1f4bf7e5af971bc54ec Mon Sep 17 00:00:00 2001 From: Pierre-Francois Loos Date: Thu, 12 Dec 2019 18:13:35 +0100 Subject: [PATCH] Sec IID --- Manuscript/srDFT_SC.tex | 52 +++++++++++++++++++++-------------------- 1 file changed, 27 insertions(+), 25 deletions(-) diff --git a/Manuscript/srDFT_SC.tex b/Manuscript/srDFT_SC.tex index 3985711..03320cf 100644 --- a/Manuscript/srDFT_SC.tex +++ b/Manuscript/srDFT_SC.tex @@ -447,7 +447,7 @@ Because of the very definition of $\wbasis$, one has the following property in t which is again fundamental to guarantee the correct behavior of the theory in the CBS limit. \subsubsection{Frozen-core approximation} -As all WFT calculations in this work are performed within the frozen-core approximation, we use a valence-only version of the various quantities needed for the complementary basis functional introduced in Ref.~\onlinecite{LooPraSceTouGin-JCPL-19}. We partition the basis set as $\Bas = \Cor \bigcup \BasFC$, where $\Cor$ and $\BasFC$ are the sets of core and active orbitals, respectively, and define the valence-only local range-separation parameter as +As all WFT calculations in this work are performed within the frozen-core approximation, we use a valence-only version of the various quantities needed for the \titou{complementary basis functional} introduced in Ref.~\onlinecite{LooPraSceTouGin-JCPL-19}. We partition the basis set as $\Bas = \Cor \bigcup \BasFC$, where $\Cor$ and $\BasFC$ are the sets of core and active orbitals, respectively, and define the valence-only local range-separation parameter as \begin{equation} \label{eq:def_mur_val} \murpsival = \frac{\sqrt{\pi}}{2} \wbasiscoalval{}, @@ -475,32 +475,34 @@ is the valence-only effective interaction and One would note the restrictions of the sums to the set of active orbitals in Eqs.~\eqref{eq:fbasis_val} and \eqref{eq:twordm_val}. It is also noteworthy that, with the present definition, $\wbasisval$ still tends to the usual Coulomb interaction as $\Bas \to \CBS$. -\subsection{Generic form and properties of the approximations for $\efuncden{\den}$ } +\subsection{Short-range correlation density functional approximations} \label{sec:functional} -\subsubsection{Generic form of the approximate functionals} +\subsubsection{Generic form} \label{sec:functional_form} -As originally proposed and motivated in Ref.~\onlinecite{GinPraFerAssSavTou-JCP-18}, we approximate the \titou{complementary basis} 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 total density $\denr$, the spin polarization $\zeta(\br{})=[n_\uparrow(\br{})-n_\downarrow(\br{})]/\denr$, 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, a given approximation X of $\efuncden{\den}$ will have the following generic local form -\begin{equation} - \begin{aligned} +As originally proposed and motivated in Ref.~\onlinecite{GinPraFerAssSavTou-JCP-18}, we approximate the \titou{complementary basis 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 total density $\denr$, the spin polarization $\zeta(\br{})=[n_\uparrow(\br{})-n_\downarrow(\br{})]/\denr$, 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$. +\trashPFL{Therefore, a given approximation X of $\efuncden{\den}$ will have the following generic local form} +\titou{Therefore, $\efuncden{\den}$ has the following generic form} +\begin{multline} \label{eq:def_ecmdpbebasis} - &\efuncdenpbe{\argebasis} = \;\;\;\;\;\;\;\; \\ &\int \d\br{} \,\denr \ecmd(\argrebasis), - \end{aligned} -\end{equation} -where $\ecmd(\argecmd)$ is the correlation energy per particle taken as + \efuncdenpbe{\argebasis} = + \\ + \int \d\br{} \,\denr \ecmd(\argrebasis), +\end{multline} +where \begin{equation} \label{eq:def_ecmdpbe} - \ecmd(\argecmd) = \frac{\varepsilon_{\text{c}}^{\text{PBE}}(\argepbe)}{1+ \beta(\argepbe) \; \mu^3}, + \ecmd(\argecmd) = \frac{\varepsilon_{\text{c}}^{\text{PBE}}(\argepbe)}{1+ \beta(\argepbe,\titou{\ntwo}) \; \mu^3}, \end{equation} -with +is the correlation energy per particle, with \begin{equation} \label{eq:def_beta} - \beta(\argepbe) = \frac{3}{2\sqrt{\pi}(1 - \sqrt{2})}\frac{\varepsilon_{\text{c}}^{\text{PBE}}(\argepbe)}{\ntwo/\den}, + \beta(\argepbe,\titou{\ntwo}) = \frac{3}{2\sqrt{\pi}(1 - \sqrt{2})}\frac{\varepsilon_{\text{c}}^{\text{PBE}}(\argepbe)}{\ntwo/\den}, \end{equation} -where $\varepsilon_{\text{c}}^{\text{PBE}}(\argepbe)$ is the usual PBE correlation energy per particle~\cite{PerBurErn-PRL-96}. Before introducing the different flavors of approximate functionals that we will use here (see Section~\ref{sec:def_func}), we would like to give some motivations for this choice of functional form. +where $\varepsilon_{\text{c}}^{\text{PBE}}(\argepbe)$ is the usual PBE correlation energy per particle. \cite{PerBurErn-PRL-96} Before introducing the different flavors of approximate functionals that we will use here (see Sec.~\ref{sec:def_func}), we would like to give some motivations for this choice of functional form. -The functional form of $\ecmd(\argecmd)$ in Eq.~\ref{eq:def_ecmdpbe} has been originally proposed in Ref.~\onlinecite{FerGinTou-JCP-18} in the context of RSDFT. In the $\mu\to 0$ limit, it reduces to the usual PBE correlation functional +The functional form of $\ecmd(\argecmd)$ in Eq.~\ref{eq:def_ecmdpbe} has been originally proposed in Ref.~\onlinecite{FerGinTou-JCP-18} in the context of RSDFT. In the $\mu\to 0$ limit, it reduces to the usual PBE correlation functional, \ie, \begin{equation} \lim_{\mu \to 0} \ecmd(\argecmd) = \varepsilon_{\text{c}}^{\text{PBE}}(\argepbe), \end{equation} @@ -509,7 +511,7 @@ which is relevant in the weak-correlation (or high-density) limit. In the large- \label{eq:lim_mularge} \ecmd(\argecmd) \isEquivTo{\mu\to\infty} \frac{2\sqrt{\pi}(1 - \sqrt{2})}{3 \mu^3} \frac{\ntwo}{n}, \end{equation} -which is the exact large-$\mu$ behavior of the exact ECMD correlation energy~\cite{PazMorGorBac-PRB-06,FerGinTou-JCP-18}. Of course, for a specific system, the large-$\mu$ behavior will be exact only if one uses for $n_2$ the \textit{exact} on-top pair density of this system. This large-$\mu$ limit in Eq.~\eqref{eq:lim_mularge} is relevant in the strong-correlation (or low-density) limit. In the context of RSDFT, some of the present authors have illustrated in Ref.~\onlinecite{FerGinTou-JCP-18} that the on-top pair density involved in Eq.~\eqref{eq:def_ecmdpbe} plays indeed a crucial role when reaching the strong-correlation regime. The importance of the on-top pair density in the strong-correlation regime have been also recently acknowledged by Gagliardi and coworkers~\cite{CarTruGag-JPCA-17} and Pernal and coworkers\cite{GritMeePer-PRA-18}. +which is the exact large-$\mu$ behavior of the exact ECMD correlation energy. \cite{PazMorGorBac-PRB-06,FerGinTou-JCP-18} Of course, for a specific system, the large-$\mu$ behavior will be exact only if one uses for $n_2$ the \textit{exact} on-top pair density of this system. This large-$\mu$ limit in Eq.~\eqref{eq:lim_mularge} is relevant in the strong-correlation (or low-density) limit. In the context of RSDFT, some of the present authors have illustrated in Ref.~\onlinecite{FerGinTou-JCP-18} that the on-top pair density involved in Eq.~\eqref{eq:def_ecmdpbe} plays indeed a crucial role when reaching the strong-correlation regime. The importance of the on-top pair density in the strong-correlation regime have been also recently acknowledged by Gagliardi and coworkers \cite{CarTruGag-JPCA-17} and Pernal and coworkers.\cite{GritMeePer-PRA-18} Note also that $\ecmd(\argecmd)$ vanishes when $\ntwo$ vanishes \begin{equation} @@ -521,19 +523,19 @@ which is expected for systems with a vanishing on-top pair density, such as the \lim_{\mu \rightarrow \infty} \ecmd(\argecmd) = 0. \end{equation} -\subsubsection{Properties of approximate functionals} +\subsubsection{Properties} \label{sec:functional_prop} -Within the definitions of Eqs.~\eqref{eq:def_mur} and \eqref{eq:def_ecmdpbebasis}, any approximate complementary basis functional $\efuncdenpbe{\argebasis}$ satisfies two important properties. +Within the definitions of Eqs.~\eqref{eq:def_mur} and \eqref{eq:def_ecmdpbebasis}, any approximate \titou{complementary basis functional} $\efuncdenpbe{\argebasis}$ satisfies two important properties. -First, thanks to the properties in Eq.~\eqref{eq:cbs_mu} and~\eqref{eq:lim_muinf}, $\efuncdenpbe{\argebasis}$ vanishes in the CBS limit, whatever the wave function $\psibasis$ used to define the local range-separation parameter $\mu(\br{})$, +First, thanks to the properties in Eq.~\eqref{eq:cbs_mu} and~\eqref{eq:lim_muinf}, $\efuncdenpbe{\argebasis}$ vanishes in the CBS limit, independently of the wave function $\psibasis$ used to define the local range-separation parameter $\mu(\br{})$, \begin{equation} \label{eq:lim_ebasis} \lim_{\basis \rightarrow \text{CBS}} \efuncdenpbe{\argebasis} = 0, \quad \forall\, \psibasis, \end{equation} which guarantees an unaltered CBS limit. -Second, the fact that $\efuncdenpbe{\argebasis}$ vanishes for systems with vanishing on-top pair density guarantees the correct limit for one-electron systems and for the stretched H$_2$ molecule. This property is guaranteed independently by i) the definition of the effective interaction $\wbasis$ (see Eq.~\eqref{eq:wbasis}) together with the condition in Eq.~\eqref{eq:lim_muinf}, ii) the fact that $\ecmd(\argecmd)$ vanishes when the on-top pair density vanishes (see Eq.~\eqref{eq:lim_n2}). +Second, the fact that $\efuncdenpbe{\argebasis}$ vanishes for systems with vanishing on-top pair density guarantees the correct limit for one-electron systems and for the stretched H$_2$ molecule. This property is guaranteed independently by i) the definition of the effective interaction $\wbasis$ [see Eq.~\eqref{eq:wbasis}] together with the condition in Eq.~\eqref{eq:lim_muinf}, ii) the fact that $\ecmd(\argecmd)$ vanishes when the on-top pair density vanishes [see Eq.~\eqref{eq:lim_n2}]. \subsection{Requirements for the approximate functionals in the strong-correlation regime} \label{sec:requirements} @@ -601,21 +603,21 @@ We define the following functionals:\\ i) PBE-UEG-$\tilde{\zeta}$ which uses the UEG on-top pair density defined in Eq.~\eqref{eq:def_n2ueg} and the effective spin polarization of Eq.~\eqref{eq:def_effspin} \begin{equation} - \label{eq:def_pbeueg} + \label{eq:def_pbeueg_i} \begin{aligned} \bar{E}^\Bas_{\pbeuegXi} = \int \d\br{} \,\denr \ecmd(\argrpbeuegXi), \end{aligned} \end{equation} ii) PBE-ot-$\tilde{\zeta}$ which uses the on-top pair density of Eq.~\eqref{eq:def_n2extrap} \begin{equation} - \label{eq:def_pbeueg} + \label{eq:def_pbeueg_ii} \begin{aligned} \bar{E}^\Bas_{\pbeontXi} = \int \d\br{} \,\denr \ecmd(\argrpbeontXi), \end{aligned} \end{equation} iii) PBE-ot-$0{\zeta}$ where uses zero spin polarization and the on-top pair density of Eq.~\eqref{eq:def_n2extrap} \begin{equation} - \label{eq:def_pbeueg} + \label{eq:def_pbeueg_iii} \begin{aligned} \bar{E}^\Bas_{\pbeontns} = \int \d\br{} \,\denr \ecmd(\argrpbeontns). \end{aligned} @@ -646,7 +648,7 @@ The purpose of the present paper being the study of the basis-set correction in In the case of the C$_2$, N$_2$, O$_2$, and F$_2$ molecules, approximations to the FCI energies are obtained using converged frozen-core (1s orbitals are kept frozen) CIPSI calculations and the extrapolation scheme for the perturbative correction of Umrigar \textit{et. al.} (see Refs.~\onlinecite{HolUmrSha-JCP-17, SceGarCafLoo-JCTC-18, LooSceBloGarCafJac-JCTC-18, SceBenJacCafLoo-JCP-18, LooBogSceCafJac-JCTC-19, QP2} for more details) using the Quantum Package software~\cite{QP2}. The estimated exact potential energy curves are obtained from Ref.~\onlinecite{LieCle-JCP-74a}. For all geometries and basis sets, the error with respect to the exact FCI energies are estimated to be below 0.5 mH. In the case of the H$_{10}$ chain, the approximation to the FCI energies together with the estimated exact potential energy curves are obtained from the data of Ref.~\onlinecite{h10_prx} where the authors performed MRCI+Q calculations with a minimal valence active space as reference (see below for the description of the active space). -Regarding the complementary basis functional, we first perform full-valence complete-active-space self-consistent-field (CASSCF) calculations with the GAMESS-US software\cite{gamess} to obtain the wave function $\psibasis$. Then, all density-like quantities involved in the functional (density $n(\br{})$, spin polarization $\zeta(\br{})$, reduced density gradient $s(\br{})$, and on-top pair density $n_2(\br{})$) together with the local range-separation function $\mu(\br{})$ of Eq.~\eqref{eq:def_mur} are calculated with this full-valence CASSCF wave function. These CASSCF wave functions correspond to the following active spaces: 10 electrons in 10 orbitals for H$_{10}$, 8 electrons in 8 electrons for C$_2$, 10 electrons in 8 orbitals for N$_2$, 12 electrons in 8 orbitals for O$_2$, and 14 electrons in 8 orbitals for F$_2$. +Regarding the \titou{complementary basis functional}, we first perform full-valence complete-active-space self-consistent-field (CASSCF) calculations with the GAMESS-US software\cite{gamess} to obtain the wave function $\psibasis$. Then, all density-like quantities involved in the functional (density $n(\br{})$, spin polarization $\zeta(\br{})$, reduced density gradient $s(\br{})$, and on-top pair density $n_2(\br{})$) together with the local range-separation function $\mu(\br{})$ of Eq.~\eqref{eq:def_mur} are calculated with this full-valence CASSCF wave function. These CASSCF wave functions correspond to the following active spaces: 10 electrons in 10 orbitals for H$_{10}$, 8 electrons in 8 electrons for C$_2$, 10 electrons in 8 orbitals for N$_2$, 12 electrons in 8 orbitals for O$_2$, and 14 electrons in 8 orbitals for F$_2$. Also, as the frozen core approximation is used in all our CIPSI calculations, we use the corresponding valence-only complementary functionals. Therefore, all density-like quantities exclude any contribution from the core 1s orbitals, and the range-separation function is taken as the one defined in Eq.~\eqref{eq:def_mur_val}. @@ -765,7 +767,7 @@ We report in Figure \ref{fig:C2}, \ref{fig:N2}, \ref{fig:O2}, and \ref{fig:F2} t Just as in the case of \ce{H10}, the quality of $D_0$ are globally improved by adding the basis-set correction and it is remarkable that PBE-ot-$\tilde{\zeta}$ and PBE-ot-$0{\zeta}$ give very similar results. The latter observation confirms that the dependency on the on-top pair density allows one to remove the dependency of any kind of spin polarization for a quite wide range of electron density and also for a high-spin system like \ce{O2}. More quantitatively, an error below 1.0 mH on the estimated exact valence-only $D_0$ is found for \ce{N2}, \ce{O2}, and \ce{F2} with the aug-cc-pVTZ basis set using the PBE-ot-$0{\zeta}$ functional, whereas such a result is far from reach within the same basis set at the near-FCI level. In the case of \ce{C2} with the aug-cc-pVTZ basis set, an error of about 5.5 mH is found with respect to the estimated exact $D_0$. Such an error is remarkably large with respect to the other diatomic molecules studied here and might be associated to the level of strong correlation in the \ce{C2} molecule. -Regarding now the performance of the basis-set correction along the whole potential energy curve, it is interesting to notice that it fails to provide a noticeable improvement far from the equilibrium geometry. Acknowledging that the weak correlation effects in these regions are dominated by dispersion interactions which are long-range effects, the failure of the present approximations for the complementary basis functionals can be understood easily. Indeed, the whole scheme designed here is based on the physics of correlation near the electron-electron cusp: the local range-separation function $\mu(\br{})$ is designed by looking at the electron-electron coalescence point and the ECMD functionals are suited for short-range correlation effects. Therefore, the failure of the present basis-set correction to describe dispersion interactions is expected. +Regarding now the performance of the basis-set correction along the whole potential energy curve, it is interesting to notice that it fails to provide a noticeable improvement far from the equilibrium geometry. Acknowledging that the weak correlation effects in these regions are dominated by dispersion interactions which are long-range effects, the failure of the present approximations for the \titou{complementary basis functionals} can be understood easily. Indeed, the whole scheme designed here is based on the physics of correlation near the electron-electron cusp: the local range-separation function $\mu(\br{})$ is designed by looking at the electron-electron coalescence point and the ECMD functionals are suited for short-range correlation effects. Therefore, the failure of the present basis-set correction to describe dispersion interactions is expected. \section{Conclusion} \label{sec:conclusion}