diff --git a/Manuscript/FarDFT.tex b/Manuscript/FarDFT.tex index 8935ee7..d15e516 100644 --- a/Manuscript/FarDFT.tex +++ b/Manuscript/FarDFT.tex @@ -40,7 +40,7 @@ \newcommand{\hH}{\Hat{H}} \newcommand{\hHc}{\Hat{h}} \newcommand{\hT}{\Hat{T}} -\newcommand{\bH}{\Hat{T}} +\newcommand{\bH}{\bm{H}} \newcommand{\hVext}{\Hat{V}_\text{ext}} \newcommand{\vext}{v_\text{ext}} \newcommand{\hWee}{\Hat{W}_\text{ee}} @@ -68,6 +68,8 @@ \newcommand{\EFCI}{E_\text{FCI}} \newcommand{\HF}{\text{HF}} \newcommand{\LDA}{\text{LDA}} +\newcommand{\eLDA}{\text{eLDA}} +\newcommand{\CID}{\text{CID}} \newcommand{\Hxc}{\text{Hxc}} \newcommand{\Ha}{\text{H}} \newcommand{\ex}{\text{x}} @@ -82,7 +84,7 @@ \newcommand{\bGamma}[1]{\bm{\Gamma}^{#1}} \newcommand{\bHc}{\bm{h}} \newcommand{\bF}[1]{\bm{F}^{#1}} -\newcommand{\Ex}[1]{\Omega^{#1}} +\newcommand{\Ex}[2]{\Omega_{#1}^{#2}} % elements \newcommand{\ew}[1]{w_{#1}} @@ -91,6 +93,8 @@ \newcommand{\eGamma}[2]{\Gamma_{#1}^{#2}} \newcommand{\hGamma}[2]{\Hat{\Gamma}_{#1}^{#2}} \newcommand{\eHc}[1]{h_{#1}} +\newcommand{\eJ}[1]{J_{#1}} +\newcommand{\eK}[1]{K_{#1}} \newcommand{\eF}[2]{F_{#1}^{#2}} \newcommand{\ON}[2]{f_{#1}^{#2}} @@ -104,6 +108,8 @@ \newcommand{\cMO}[2]{c_{#1}^{#2}} \newcommand{\AO}[1]{\chi_{#1}} +\newcommand{\RHH}{R_{\ce{H-H}}} + % units \newcommand{\IneV}[1]{#1 eV} \newcommand{\InAU}[1]{#1 a.u.} @@ -200,7 +206,7 @@ Degeneracies can be easily handled. One of the key feature of eDFT in the present context is that one can easily extract individual excitation energies from the ensemble energy via differentiation with respect to individual weights: \begin{equation} - \pdv{\E{}{\bw}}{\ew{I}} = \E{}{(I)} - \E{}{(0)}, + \pdv{\E{}{\bw}}{\ew{I}} = \E{}{(I)} - \E{}{(0)} = \Ex{}{(I)}, \end{equation} where the weights are normalised by setting $\ew{0} = 1 - \sum_{I \ne 0} \ew{I}$. @@ -469,7 +475,159 @@ In the case of a homogeneous system (or equivalently within the LDA), substituti %%%%%%%%%%%%%%% \section{Results} \label{sec:res} -Here, we do \ce{H2} because \ce{H2} is very interesting. +Here, we consider as testing ground the minimal-basis \ce{H2} molecule. +We select STO-3G as minimal basis, and study the behaviour of the total energy of \ce{H2} as a function of the internuclear distance $\RHH$ (in bohr). +The bonding and antibonding orbitals of the \ce{H2} molecule are given by +\begin{subequations} +\begin{align} + \MO{1}{}(\br{}) & = \qty[ \AO{A}(\br{}) + \AO{B}(\br{}) ]/\sqrt{2 + S_{AB}}, + \\ + \MO{2}{}(\br{}) & = \qty[ \AO{A}(\br{}) - \AO{B}(\br{}) ]/\sqrt{2 - S_{AB}}, +\end{align} +\end{subequations} +where $\AO{A}$ and $\AO{B}$ are the two contracted Gaussian basis functions centred on each of the nucleus, and $S_{AB} = \braket{\AO{A}}{\AO{B}}$. + +As reference results, we consider CID (configuration interaction with doubles) computed in the same (minimal) basis set. +The CID energies of the ground state and doubly-excited states are provided by the eigenvalues of the following CID matrix: +\begin{equation} + \bH_\CID = + \begin{pmatrix} + \E{\HF}{(0)} & \eK{12} + \\ + \eK{12} & \E{\HF}{(1)} + \end{pmatrix}, +\end{equation} +with +\begin{subequations} +\begin{align} + \label{eq:HF0} + \E{\HF}{(0)} & = \eHc{1} + 2 \eJ{11} - \eK{11}, + \\ + \label{eq:HF1} + \E{\HF}{(1)} & = \eHc{2} + 2 \eJ{22} - \eK{22}, +\end{align} +\end{subequations} +and +\begin{subequations} +\begin{align} + \eHc{p} & = \int \MO{p}{}(\br{}) \qty[-\frac{\nabla^2}{2} + \vext(\br{})] \MO{p}{}(\br{})d\br{}, + \\ + \eJ{pq} & = \iint \frac{\MO{p}{}(\br{})\MO{p}{}(\br{}) \MO{q}{}(\br{}')\MO{q}{}(\br{}')}{\abs{\br{} - \br{}'}} d\br{} d\br{}', + \\ + \eK{pq} & = \iint \frac{\MO{p}{}(\br{})\MO{q}{}(\br{}) \MO{q}{}(\br{}')\MO{p}{}(\br{}')}{\abs{\br{} - \br{}'}} d\br{} d\br{}'. +\end{align} +\end{subequations} +Note that, in the HF case, there is no self-interaction error as $\eJ{pp} = \eK{pp}$. + +The CID energies are explicitly given by +\begin{subequations} +\begin{align} + \E{\CID}{(0)} & = \frac{\E{\HF}{(0)} + \E{\HF}{(1)}}{2} - \frac{1}{2} \sqrt{\qty(\E{\HF}{(1)} - \E{\HF}{(0)})^2 + 4 \eK{12}^2}, + \\ + \E{\CID}{(1)} & = \frac{\E{\HF}{(0)} + \E{\HF}{(1)}}{2} + \frac{1}{2} \sqrt{\qty(\E{\HF}{(1)} - \E{\HF}{(0)})^2 + 4 \eK{12}^2}, +\end{align} +\end{subequations} +and the CID excitation energy reads +\begin{equation} + \Ex{\CID}{(1)} = \sqrt{\qty(\E{\HF}{(1)} - \E{\HF}{(0)})^2 + 4 \eK{12}^2} \ge \Ex{\HF}{(1)}. +\end{equation} + +At the (ground-state) LDA level (\ie, we only consider ground-state functionals), these energies reads +\begin{subequations} +\begin{align} + \label{eq:LDA0} + \E{\LDA}{(0)} & = \eHc{1} + 2 \eJ{11} + \int \e{\xc}{\LDA}[\n{}{(0)}(\br{})] \n{}{(0)}(\br{}) d\br{}, + \\ + \label{eq:LDA1} + \E{\LDA}{(1)} & = \eHc{2} + 2 \eJ{22} + \int \e{\xc}{\LDA}[\n{}{(1)}(\br{})] \n{}{(1)}(\br{}) d\br{}, +\end{align} +\end{subequations} +with +\begin{align} + \n{}{(0)}(\br{}) & = 2 \MO{1}{2}(\br{}), + & + \n{}{(1)}(\br{}) & = 2 \MO{2}{2}(\br{}), +\end{align} +Note that, contrary to the HF case, self-interaction is present in LDA. + +At the eLDA, we have +\begin{subequations} +\begin{align} + \label{eq:eLDA0} + \E{\eLDA}{(0)} & = \eHc{1} + 2 \eJ{11} + \int \be{\xc}{(0)}[\n{}{(0)}(\br{})] \n{}{(0)}(\br{}) d\br{}, + \\ + \label{eq:eLDA1} + \E{\eLDA}{(1)} & = \eHc{2} + 2 \eJ{22} + \int \be{\xc}{(1)}[\n{}{(1)}(\br{})] \n{}{(1)}(\br{}) d\br{}, +\end{align} +\end{subequations} +with $\be{\xc}{(0)} \equiv \e{\xc}{\LDA}$ and $\be{\xc}{(1)}(\n{}{}) = \e{\xc}{\LDA}(\n{}{}) + \e{\xc}{(1)}(\n{}{}) - \e{\xc}{(0)}(\n{}{})$. + +%\titou{Note that we do not consider symmetry-broken solutions.} + +Interestingly here, there is a strong connection between the LDA and eLDA excitation energies: +\begin{equation} + \Ex{\eLDA}{(1)} = \Ex{\LDA}{(1)} + \int \qty( \e{\xc}{(1)} - \e{\xc}{(0)} )[\n{}{(1)}(\br{})] \n{}{(1)}(\br{}) d\br{}. +\end{equation} + +These equations can be combined to define three ensemble energies +\begin{subequations} +\begin{align} + \label{eq:EwHF} + \E{\HF}{\ew{}} & = (1-\ew{}) \E{\HF}{(0)} + \ew{} \E{\HF}{(1)}, + \\ + \label{eq:EwLDA} + \E{\LDA}{\ew{}} & = (1-\ew{}) \E{\LDA}{(0)} + \ew{} \E{\LDA}{(1)}, + \\ + \label{eq:EweLDA} + \E{\eLDA}{\ew{}} & = (1-\ew{}) \E{\eLDA}{(0)} + \ew{} \E{\eLDA}{(1)}, +\end{align} +\end{subequations} +which are all, by construction, linear with respect to $\ew{}$. + +These energies given in Eqs.~\eqref{eq:EwHF}, \eqref{eq:EwLDA} and \eqref{eq:EweLDA} can also be obtained directly from the ensemble density $\n{}{\ew{}} = (1-\ew{}) \n{}{(0)} + \ew{} \n{}{(1)}$. +(This is what one would do in practice, \ie, by performing a KS ensemble calculation.) +We will label these energies as $\bE{}{\ew{}}$. + +For HF, we have +\begin{equation} +\begin{split} + \bE{\HF}{\ew{}} + & = (1-\ew{}) \eHc{1} + \ew{} \eHc{2} + + \frac{1}{2} \iint \frac{\n{}{\ew{}}(\br{})\n{}{\ew{}}(\br{}')}{\abs{\br{} - \br{}'}} d\br{} d\br{}' + \\ + & = \ldots +\end{split} +\end{equation} +which is clearly quadratic with respect to $\ew{}$. + +For LDA, we have +\begin{equation} +\begin{split} + \bE{\LDA}{\ew{}} + & = (1-\ew{}) \eHc{1} + \ew{} \eHc{2} + + \iint \frac{\n{}{\ew{}}(\br{})\n{}{\ew{}}(\br{}')}{\abs{\br{} - \br{}'}} d\br{} d\br{}' + \\ + & + \int \e{\xc}{\LDA}[\n{}{\ew{}}(\br{})] \n{}{\ew{}}(\br{}) d\br{} + \\ + & = \ldots , +\end{split} +\end{equation} +which is also clearly quadratic with respect to $\ew{}$ because the (weight-independent) LDA functional cannot compensate the ``quadraticity'' of the Hartree term. + +For eLDA, we have +\begin{equation} +\begin{split} + \bE{\eLDA}{\ew{}} + & = (1-\ew{}) \eHc{1} + \ew{} \eHc{2} + + \iint \frac{\n{}{\ew{}}(\br{})\n{}{\ew{}}(\br{}')}{\abs{\br{} - \br{}'}} d\br{} d\br{}' + \\ + & + \int \be{\xc}{\ew{}}[\n{}{\ew{}}(\br{})] \n{}{\ew{}}(\br{}) d\br{} + \\ + & = \ldots , +\end{split} +\end{equation} +which *could* be linear with respect to the weight if the weight-dependent xc functional is very well constructed. +This would be, for example, the case with the exact xc functional. %%%%%%%%%%%%%%%%%% %%% CONCLUSION %%%