corrections toto

This commit is contained in:
Pierre-Francois Loos 2019-04-24 15:51:42 +02:00
parent 41b458639d
commit 844fc59cfa

View File

@ -169,12 +169,12 @@ By increasing the excitation degree of the CC expansion, one can systematically
One of the most fundamental drawbacks of conventional WFT methods is the slow convergence of energies and properties with respect to the size of the one-electron basis set.
This undesirable feature was put into light by Kutzelnigg more than thirty years ago. \cite{Kut-TCA-85}
To palliate this, following Hylleraas' footsteps, \cite{Hyl-ZP-29} Kutzelnigg proposed to introduce explicitly the interelectronic distance $r_{12} = \abs{\br{1} - \br{2}}$ to properly describe the electronic wave function around the coalescence of two electrons. \cite{Kut-TCA-85, KutKlo-JCP-91, NogKut-JCP-94}
The resulting F12 methods \trashAS{yields} \toto{yield} a prominent improvement of the energy convergence, and achieve chemical accuracy for small organic molecules with relatively small Gaussian basis sets. \cite{Ten-TCA-12, TenNog-WIREs-12, HatKloKohTew-CR-12, KonBisVal-CR-12}
The resulting F12 methods yield a prominent improvement of the energy convergence, and achieve chemical accuracy for small organic molecules with relatively small Gaussian basis sets. \cite{Ten-TCA-12, TenNog-WIREs-12, HatKloKohTew-CR-12, KonBisVal-CR-12}
For example, at the CCSD(T) level, one can obtain quintuple-$\zeta$ quality correlation energies with a triple-$\zeta$ basis, \cite{TewKloNeiHat-PCCP-07} although computational overheads are introduced by the large auxiliary basis used to resolve three- and four-electron integrals. \cite{BarLoo-JCP-17}
To reduce further the computational cost and/or ease the transferability of the F12 correction, approximated and/or universal schemes have recently emerged. \cite{TorVal-JCP-09, KonVal-JCP-10, KonVal-JCP-11, BooCleAlaTew-JCP-2012, IrmHumGru-arXiv-2019, IrmGru-arXiv-2019}
Present-day DFT calculations are almost exclusively done within the so-called Kohn-Sham (KS) formalism, which corresponds to an exact dressed one-electron theory. \cite{KohSha-PR-65}
\trashAS{DFT's attractiveness} \toto{The attractiveness of DFT} originates from its very favorable \trashAS{cost/accuracy} \toto{accuracy/cost} ratio as it often provides reasonably accurate energies and properties at a relatively low computational cost.
The attractiveness of DFT originates from its very favorable accuracy/cost ratio as it often provides reasonably accurate energies and properties at a relatively low computational cost.
Thanks to this, KS-DFT \cite{HohKoh-PR-64, KohSha-PR-65} has become the workhorse of electronic structure calculations for atoms, molecules and solids. \cite{ParYan-BOOK-89}
Although there is no clear way on how to systematically improve density-functional approximations, \cite{Bec-JCP-14} climbing Perdew's ladder of DFT is potentially the most satisfactory way forward. \cite{PerSch-AIPCP-01, PerRuzTaoStaScuCso-JCP-05}
In the context of the present work, one of the interesting feature of density-based methods is their much faster convergence with respect to the size of the basis set. \cite{FraMusLupTou-JCP-15}
@ -219,7 +219,7 @@ This implies that
\end{equation}
where $\E{\modY}{}$ is the energy associated with the method $\modY$ in the CBS limit.
In the case where $\modY = \FCI$ in Eq.~\eqref{eq:limitfunc}, we have a strict equality as $\E{\FCI}{} = \E{}{}$.
Provided that the functional $\bE{}{\Bas}[\n{}{}]$ is known exactly, the only \trashAS{source of error at this stage lies} \toto{sources of error at this stage lie} in the potential approximate nature of the methods $\modY$ and $\modZ$, \trashAS{and the lack of self-consistency} \toto{and to the absence of a self-consistent scheme}.
Provided that the functional $\bE{}{\Bas}[\n{}{}]$ is known exactly, the only sources of error at this stage lie in the potential approximate nature of the methods $\modY$ and $\modZ$, and the lack of self-consistency in the present scheme.
The functional $\bE{}{\Bas}[\n{}{}]$ is obviously \textit{not} universal as it depends on $\Bas$.
Moreover, as $\bE{}{\Bas}[\n{}{}]$ aims at fixing the incompleteness of $\Bas$, its main role is to correct
@ -229,7 +229,7 @@ Therefore, as we shall do later on, it feels natural to approximate $\bE{}{\Bas
Contrary to the conventional RS-DFT scheme which requires a range-separation \textit{parameter} $\rsmu{}{}$, here we use a range-separation \textit{function} $\rsmu{}{\Bas}(\br{})$ that automatically adapts to quantify the incompleteness of $\Bas$ in $\mathbb{R}^3$.
% https://english.stackexchange.com/questions/61600/consist-in-vs-consist-of
The first step of the present basis-set correction consists \trashAS{of} \toto{in} obtaining an effective two-electron interaction $\W{}{\Bas}(\br{1},\br{2})$ ``mimicking'' the Coulomb operator in an incomplete basis $\Bas$.
The first step of the present basis-set correction consists in obtaining an effective two-electron interaction $\W{}{\Bas}(\br{1},\br{2})$ ``mimicking'' the Coulomb operator in an incomplete basis $\Bas$.
In a second step, we shall link $\W{}{\Bas}(\br{1},\br{2})$ to $\rsmu{}{\Bas}(\br{})$.
As a final step, we employ short-range density functionals \cite{TouGorSav-TCA-05} with $\rsmu{}{\Bas}(\br{})$ as range-separation function.
@ -329,7 +329,7 @@ The local-density approximation (LDA) of the ECMD complementary functional is de
\label{eq:def_lda_tot}
\titou{\bE{\LDA}{\Bas}[\n{}{},\rsmu{}{\Bas}] = \int \n{}{}(\br{}) \be{\text{c,md}}{\sr,\LDA}\qty(\n{}{}(\br{}),\zeta(\br{}),\rsmu{}{\Bas}(\br{})) \dbr{},}
\end{equation}
where \titou{$\zeta = (\n{\uparrow}{} - \n{\downarrow}{})/\n{}{}$ is the spin-polarization and $\be{\text{c,md}}{\sr,\LDA}(\n{}{},\zeta,\rsmu{}{})$} is the ECMD correlation energy per electron of the uniform electron gas (UEG) \cite{LooGil-WIRES-16} \trashAS{parametrized} \toto{parameterized} in Ref.~\onlinecite{PazMorGorBac-PRB-06} as a function of the spin densities $\qty{\n{\sigma}{}}_{\sigma=\uparrow,\downarrow}$ and the range-separation parameter $\mu$.
where \titou{$\zeta = (\n{\uparrow}{} - \n{\downarrow}{})/\n{}{}$ is the spin-polarization and $\be{\text{c,md}}{\sr,\LDA}(\n{}{},\zeta,\rsmu{}{})$} is the ECMD correlation energy per electron of the uniform electron gas (UEG) \cite{LooGil-WIRES-16} parameterized in Ref.~\onlinecite{PazMorGorBac-PRB-06} as a function of the spin densities $\qty{\n{\sigma}{}}_{\sigma=\uparrow,\downarrow}$ and the range-separation parameter $\mu$.
The short-range LDA correlation functional relies on the transferability of the physics of the UEG which is certainly valid for large $\mu$ but is known to over correlate for small $\mu$.
In order to correct such a defect, we propose here a new Perdew-Burke-Ernzerhof (PBE)-based ECMD functional
\begin{equation}
@ -468,13 +468,13 @@ As a method $\modY$ we employ either CCSD(T) or exFCI.
Here, exFCI stands for extrapolated FCI energies computed with the CIPSI algorithm. \cite{HurMalRan-JCP-73, GinSceCaf-CJC-13, GinSceCaf-JCP-15}
We refer the interested reader to Refs.~\onlinecite{HolUmrSha-JCP-17, SceGarCafLoo-JCTC-18, LooSceBloGarCafJac-JCTC-18, SceBenJacCafLoo-JCP-18, LooBogSceCafJAc-JCTC-19} for more details.
In the case of the CCSD(T) calculations, we have $\modZ = \ROHF$ as we use the restricted open-shell HF (ROHF) one-electron density to compute the complementary basis-set correction energy.
In the case of exFCI, the one-electron density is computed from a very large CIPSI expansion containing several \trashAS{millions} \toto{million} determinants.
In the case of exFCI, the one-electron density is computed from a very large CIPSI expansion containing several million determinants.
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}
\trashAS{Except for} \toto{Apart from} 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 \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}.
Apart from 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 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.
The FC density-based correction is used consistently \titou{with the FC approximation in WFT methods.}
The FC density-based correction is used consistently 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.
As the exFCI atomization energies are converged with a precision of about 0.1 {\kcal}, we can label these as near FCI.