This commit is contained in:
Pierre-Francois Loos 2019-04-23 23:37:58 +02:00
parent 3421d7ce07
commit 484fc3f0f8

View File

@ -146,9 +146,9 @@
\centering \centering
\includegraphics[width=\linewidth]{TOC} \includegraphics[width=\linewidth]{TOC}
\end{wrapfigure} \end{wrapfigure}
We report a universal density-based basis set incompleteness correction that can be applied to any wave function method. We report a universal density-based basis-set incompleteness correction that can be applied to any wave function method.
The present correction, which appropriately vanishes in the complete basis set (CBS) limit, relies on short-range correlation density functionals (with multi-determinant reference) from range-separated density-functional theory (RS-DFT) to estimate the basis set incompleteness error. The present correction, which appropriately vanishes in the complete basis set (CBS) limit, relies on short-range correlation density functionals (with multi-determinant reference) from range-separated density-functional theory (RS-DFT) to estimate the basis-set incompleteness error.
Contrary to conventional RS-DFT schemes which require an \textit{ad hoc} range-separation \textit{parameter} $\mu$, the key ingredient here is a range-separation \textit{function} $\mu(\bf{r})$ which automatically adapts to the spatial non-homogeneity of the basis set incompleteness error. Contrary to conventional RS-DFT schemes which require an \textit{ad hoc} range-separation \textit{parameter} $\mu$, the key ingredient here is a range-separation \textit{function} $\mu(\bf{r})$ that automatically adapts to the spatial non-homogeneity of the basis-set incompleteness error.
As illustrative examples, we show how this density-based correction allows us to obtain CCSD(T) atomization and correlation energies near the CBS limit for the G2 set of molecules with compact Gaussian basis sets. As illustrative examples, we show how this density-based correction allows us to obtain CCSD(T) atomization and correlation energies near the CBS limit for the G2 set of molecules with compact Gaussian basis sets.
\end{abstract} \end{abstract}
@ -180,7 +180,6 @@ Progress toward unifying WFT and DFT are on-going.
In particular, range-separated DFT (RS-DFT) (see Ref.~\onlinecite{TouColSav-PRA-04} and references therein) rigorously combines these two approaches via a decomposition of the electron-electron (e-e) interaction into a non-divergent long-range part and a (complementary) short-range part treated with WFT and DFT, respectively. In particular, range-separated DFT (RS-DFT) (see Ref.~\onlinecite{TouColSav-PRA-04} and references therein) rigorously combines these two approaches via a decomposition of the electron-electron (e-e) interaction into a non-divergent long-range part and a (complementary) short-range part treated with WFT and DFT, respectively.
As the WFT method is relieved from describing the short-range part of the correlation hole around the e-e coalescence points, the convergence with respect to the one-electron basis set is greatly improved. \cite{FraMusLupTou-JCP-15} As the WFT method is relieved from describing the short-range part of the correlation hole around the e-e coalescence points, the convergence with respect to the one-electron basis set is greatly improved. \cite{FraMusLupTou-JCP-15}
Therefore, a number of approximate RS-DFT schemes have been developed within single-reference \cite{AngGerSavTou-PRA-05, GolWerSto-PCCP-05, TouGerJanSavAng-PRL-09,JanHenScu-JCP-09, TouZhuSavJanAng-JCP-11, MusReiAngTou-JCP-15} or multi-reference \cite{LeiStoWerSav-CPL-97, FroTouJen-JCP-07, FroCimJen-PRA-10, HedKneKieJenRei-JCP-15, HedTouJen-JCP-18, FerGinTou-JCP-18} WFT approaches. Therefore, a number of approximate RS-DFT schemes have been developed within single-reference \cite{AngGerSavTou-PRA-05, GolWerSto-PCCP-05, TouGerJanSavAng-PRL-09,JanHenScu-JCP-09, TouZhuSavJanAng-JCP-11, MusReiAngTou-JCP-15} or multi-reference \cite{LeiStoWerSav-CPL-97, FroTouJen-JCP-07, FroCimJen-PRA-10, HedKneKieJenRei-JCP-15, HedTouJen-JCP-18, FerGinTou-JCP-18} WFT approaches.
Very recently, a major step forward has been taken by some of the present authors thanks to the development of a density-based basis-set correction for WFT methods. \cite{GinPraFerAssSavTou-JCP-18} Very recently, a major step forward has been taken by some of the present authors thanks to the development of a density-based basis-set correction for WFT methods. \cite{GinPraFerAssSavTou-JCP-18}
The present work proposes an extension of this new methodological development alongside the first numerical tests on molecular systems. The present work proposes an extension of this new methodological development alongside the first numerical tests on molecular systems.
@ -268,7 +267,7 @@ Because Eq.~\eqref{eq:int_eq_wee} can be recast as
\iint \W{}{\Bas}(\br{1},\br{2}) \n{2}{\Bas}(\br{1},\br{2}) \dbr{1} \dbr{2}, \iint \W{}{\Bas}(\br{1},\br{2}) \n{2}{\Bas}(\br{1},\br{2}) \dbr{1} \dbr{2},
\end{equation} \end{equation}
it intuitively motivates $\W{}{\Bas}(\br{1},\br{2})$ as a potential candidate for an effective interaction. it intuitively motivates $\W{}{\Bas}(\br{1},\br{2})$ as a potential candidate for an effective interaction.
Note that the divergence condition of $\W{}{\Bas}(\br{1},\br{2})$ in Eq.~\eqref{eq:def_weebasis} ensures that one-electron systems are free of correction as the present approach must only correct the basis set incompleteness error originating from the e-e cusp. Note that the divergence condition of $\W{}{\Bas}(\br{1},\br{2})$ in Eq.~\eqref{eq:def_weebasis} ensures that one-electron systems are free of correction as the present approach must only correct the basis-set incompleteness error originating from the e-e cusp.
As already discussed in Ref.~\onlinecite{GinPraFerAssSavTou-JCP-18}, $\W{}{\Bas}(\br{1},\br{2})$ is symmetric, \textit{a priori} non translational, nor rotational invariant if $\Bas$ does not have such symmetries. As already discussed in Ref.~\onlinecite{GinPraFerAssSavTou-JCP-18}, $\W{}{\Bas}(\br{1},\br{2})$ is symmetric, \textit{a priori} non translational, nor rotational invariant if $\Bas$ does not have such symmetries.
Thanks to its definition one can show that (see Appendix B of Ref.~\onlinecite{GinPraFerAssSavTou-JCP-18}) Thanks to its definition one can show that (see Appendix B of Ref.~\onlinecite{GinPraFerAssSavTou-JCP-18})
\begin{equation} \begin{equation}
@ -470,9 +469,9 @@ In the case of exFCI, the one-electron density is computed from a very large CIP
CCSD(T) energies are computed with Gaussian09 using standard threshold values, \cite{g09} while RS-DFT and exFCI calculations are performed with {\QP}. \cite{QP2} CCSD(T) energies are computed with Gaussian09 using standard threshold values, \cite{g09} while RS-DFT and exFCI calculations are performed with {\QP}. \cite{QP2}
For the numerical quadratures, we employ the SG-2 grid. \cite{DasHer-JCC-17} For the numerical quadratures, we employ the SG-2 grid. \cite{DasHer-JCC-17}
Except for the carbon dimer where we have taken the experimental equilibrium bond length (\InAA{1.2425}), all geometries have been extracted from Ref.~\onlinecite{HauJanScu-JCP-09} and have been obtained at the B3LYP/6-31G(2df,p) level of theory. Except for the carbon dimer where we have taken the experimental equilibrium bond length (\InAA{1.2425}), all geometries have been extracted from Ref.~\onlinecite{HauJanScu-JCP-09} and have been obtained at the B3LYP/6-31G(2df,p) level of theory.
Frozen-core calculations are defined as such: a \ce{He} core is frozen from \ce{Li} to \ce{Ne}, while a \ce{Ne} core is frozen from \ce{Na} to \ce{Ar}. Frozen-core calculations are \titou{systematically performed and} defined as such: a \ce{He} core is frozen from \ce{Li} to \ce{Ne}, while a \ce{Ne} core is frozen from \ce{Na} to \ce{Ar}.
In the context of the basis-set correction, the set of active MOs, $\BasFC$, involved in the definition of the effective interaction [see Eq.~\eqref{eq:WFC}] refers to the non-frozen MOs. In the context of the basis-set correction, the set of active MOs, $\BasFC$, involved in the definition of the effective interaction [see Eq.~\eqref{eq:WFC}] refers to the non-frozen MOs.
The FC density-based correction is used consistently when the FC approximation was applied in WFT methods. The FC density-based correction is used consistently \titou{with the FC approximation in WFT methods.}
To estimate the CBS limit of each method, following Ref.~\onlinecite{HalHelJorKloKocOlsWil-CPL-98}, we perform a two-point X$^{-3}$ extrapolation of the correlation energies using the quadruple- and quintuple-$\zeta$ data that we add up to the HF energies obtained in the largest (i.e.~quintuple-$\zeta$) basis. To estimate the CBS limit of each method, following Ref.~\onlinecite{HalHelJorKloKocOlsWil-CPL-98}, we perform a two-point X$^{-3}$ extrapolation of the correlation energies using the quadruple- and quintuple-$\zeta$ data that we add up to the HF energies obtained in the largest (i.e.~quintuple-$\zeta$) basis.
As the exFCI atomization energies are converged with a precision of about 0.1 {\kcal}, we can label these as near FCI. As the exFCI atomization energies are converged with a precision of about 0.1 {\kcal}, we can label these as near FCI.
@ -498,14 +497,14 @@ The raw data can be found in the {\SI}.
A statistical analysis of these data is also provided in Table \ref{tab:stats}, where we report the mean absolute deviation (MAD), root-mean-square deviation (RMSD), and maximum deviation (MAX) with respect to the CCSD(T)/CBS atomization energies. A statistical analysis of these data is also provided in Table \ref{tab:stats}, where we report the mean absolute deviation (MAD), root-mean-square deviation (RMSD), and maximum deviation (MAX) with respect to the CCSD(T)/CBS atomization energies.
Note that the MAD of our CCSD(T)/CBS atomization energies is only 0.37 {\kcal} compared to the values extracted from Ref.~\onlinecite{HauKlo-JCP-12} which corresponds to frozen-core non-relativistic atomization energies obtained at the CCSD(T)(F12)/cc-pVQZ-F12 level of theory corrected for higher-excitation contributions ($E_\text{CCSDT(Q)/cc-pV(D+d)Z} - E_\text{CCSD(T)/cc-pV(D+d)Z})$. Note that the MAD of our CCSD(T)/CBS atomization energies is only 0.37 {\kcal} compared to the values extracted from Ref.~\onlinecite{HauKlo-JCP-12} which corresponds to frozen-core non-relativistic atomization energies obtained at the CCSD(T)(F12)/cc-pVQZ-F12 level of theory corrected for higher-excitation contributions ($E_\text{CCSDT(Q)/cc-pV(D+d)Z} - E_\text{CCSD(T)/cc-pV(D+d)Z})$.
From double-$\zeta$ to quintuple-$\zeta$ basis, the MAD associated with the CCSD(T) atomization energies goes down slowly from 14.29 to 1.28 {\kcal}. From double-$\zeta$ to quintuple-$\zeta$ basis, the MAD associated with the CCSD(T) atomization energies goes down slowly from 14.29 to 1.28 {\kcal}.
For a commonly-used basis like cc-pVTZ, the MAD of CCSD(T) is still 6.06 {\kcal}. For a commonly used basis like cc-pVTZ, the MAD of CCSD(T) is still 6.06 {\kcal}.
Applying the basis-set correction drastically reduces the basis set incompleteness error. Applying the basis-set correction drastically reduces the basis-set incompleteness error.
Already at the CCSD(T)+LDA/cc-pVDZ and CCSD(T)+PBE/cc-pVDZ levels, the MAD is reduced to 3.24 and 1.96 {\kcal}. Already at the CCSD(T)+LDA/cc-pVDZ and CCSD(T)+PBE/cc-pVDZ levels, the MAD is reduced to 3.24 and 1.96 {\kcal}.
With the triple-$\zeta$ basis, the MAD of CCSD(T)+PBE/cc-pVTZ is already below 1 {\kcal} with 36 cases (out of 55) where we achieve chemical accuracy. With the triple-$\zeta$ basis, the MAD of CCSD(T)+PBE/cc-pVTZ is already below 1 {\kcal} with 36 cases (out of 55) where we achieve chemical accuracy.
CCSD(T)+LDA/cc-pVQZ and CCSD(T)+PBE/cc-pVQZ return MAD of 0.33 and 0.31 kcal/mol (respectively) while CCSD(T)/cc-pVQZ still yields a fairly large MAD of 2.50 {\kcal}. CCSD(T)+LDA/cc-pVQZ and CCSD(T)+PBE/cc-pVQZ return MAD of 0.33 and 0.31 kcal/mol (respectively) while CCSD(T)/cc-pVQZ still yields a fairly large MAD of 2.50 {\kcal}.
Therefore, similar to F12 methods, \cite{TewKloNeiHat-PCCP-07} we can safely claim that the present basis-set correction provides significant basis set reduction and recovers quintuple-$\zeta$ quality atomization and correlation energies with triple-$\zeta$ basis sets for a much cheaper computational cost. Therefore, similar to F12 methods, \cite{TewKloNeiHat-PCCP-07} we can safely claim that the present basis-set correction provides significant basis-set reduction and recovers quintuple-$\zeta$ quality atomization and correlation energies with triple-$\zeta$ basis sets for a much cheaper computational cost.
Encouraged by these promising results, we are currently pursuing various avenues toward basis set reduction for strongly correlated systems and electronically excited states. Encouraged by these promising results, we are currently pursuing various avenues toward basis-set reduction for strongly correlated systems and electronically excited states.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section*{Supporting information} \section*{Supporting information}