999 lines
45 KiB
TeX
999 lines
45 KiB
TeX
\documentclass[aip,jcp,reprint,noshowkeys,superscriptaddress]{revtex4-1}
|
|
\usepackage{graphicx,dcolumn,bm,xcolor,microtype,multirow,amsmath,amssymb,amsfonts,physics,mhchem}
|
|
|
|
\usepackage[utf8]{inputenc}
|
|
\usepackage[T1]{fontenc}
|
|
\usepackage{txfonts}
|
|
|
|
\usepackage[
|
|
colorlinks=true,
|
|
citecolor=blue,
|
|
breaklinks=true
|
|
]{hyperref}
|
|
\urlstyle{same}
|
|
|
|
%useful stuff
|
|
\newcommand{\cdash}{\multicolumn{1}{c}{---}}
|
|
\newcommand{\mc}{\multicolumn}
|
|
\newcommand{\mr}{\multirow}
|
|
\newcommand{\fnm}{\footnotemark}
|
|
\newcommand{\fnt}{\footnotetext}
|
|
\newcommand{\tabc}[1]{\multicolumn{1}{c}{#1}}
|
|
\newcommand{\la}{\lambda}
|
|
\newcommand{\si}{\sigma}
|
|
|
|
% functionals, potentials, densities, etc
|
|
\newcommand{\eps}{\epsilon}
|
|
\newcommand{\e}[2]{\eps_\text{#1}^{#2}}
|
|
\renewcommand{\v}[2]{v_\text{#1}^{#2}}
|
|
\newcommand{\be}[2]{\bar{\eps}_\text{#1}^{#2}}
|
|
\newcommand{\bv}[2]{\bar{f}_\text{#1}^{#2}}
|
|
\newcommand{\n}[1]{n^{#1}}
|
|
\newcommand{\DD}[2]{\Delta_\text{#1}^{#2}}
|
|
\newcommand{\LZ}[2]{\Xi_\text{#1}^{#2}}
|
|
|
|
|
|
% energies
|
|
\newcommand{\EHF}{E_\text{HF}}
|
|
\newcommand{\Ec}{E_\text{c}}
|
|
\newcommand{\Ecat}{E_\text{cat}}
|
|
\newcommand{\Eneu}{E_\text{neu}}
|
|
\newcommand{\Eani}{E_\text{ani}}
|
|
\newcommand{\EPT}{E_\text{PT2}}
|
|
\newcommand{\EFCI}{E_\text{FCI}}
|
|
|
|
% matrices
|
|
\newcommand{\br}{\bm{r}}
|
|
\newcommand{\bw}{\bm{w}}
|
|
\newcommand{\bG}{\bm{G}}
|
|
\newcommand{\bS}{\bm{S}}
|
|
\newcommand{\bGamma}[1]{\bm{\Gamma}^{#1}}
|
|
\newcommand{\bH}{\bm{H}}
|
|
\newcommand{\bHc}{\bm{H}^\text{c}}
|
|
\newcommand{\bF}[1]{\bm{F}^{#1}}
|
|
\newcommand{\Ex}[1]{\Omega^{#1}}
|
|
\newcommand{\E}[1]{E^{#1}}
|
|
|
|
% elements
|
|
\newcommand{\ew}[1]{w_{#1}}
|
|
\newcommand{\eG}[1]{G_{#1}}
|
|
\newcommand{\eS}[1]{S_{#1}}
|
|
\newcommand{\eGamma}[2]{\Gamma_{#1}^{#2}}
|
|
\newcommand{\eHc}[1]{H_{#1}^\text{c}}
|
|
\newcommand{\eF}[2]{F_{#1}^{#2}}
|
|
|
|
% Numbers
|
|
\newcommand{\Nel}{N}
|
|
\newcommand{\Nbas}{K}
|
|
|
|
% Ao and MO basis
|
|
\newcommand{\MO}[2]{\phi_{#1}^{#2}}
|
|
\newcommand{\cMO}[2]{c_{#1}^{#2}}
|
|
\newcommand{\AO}[1]{\chi_{#1}}
|
|
|
|
|
|
% units
|
|
\newcommand{\IneV}[1]{#1~eV}
|
|
\newcommand{\InAU}[1]{#1~a.u.}
|
|
\newcommand{\InAA}[1]{#1~\AA}
|
|
|
|
\newcommand{\SI}{\textcolor{blue}{supplementary material}}
|
|
|
|
\newcommand{\LCPQ}{Laboratoire de Chimie et Physique Quantiques (UMR 5626), Universit\'e de Toulouse, CNRS, UPS, France}
|
|
\newcommand{\LCQ}{Laboratoire de Chimie Quantique, Institut de Chimie, CNRS, Universit\'e de Strasbourg, Strasbourg, France}
|
|
|
|
|
|
%%%% added by Manu %%%%%
|
|
\newcommand{\manu}[1]{{\textcolor{blue}{ Manu: #1 }} }
|
|
\newcommand{\beq}{\begin{eqnarray}}
|
|
\newcommand{\eeq}{\end{eqnarray}}
|
|
%
|
|
\newcommand{\bmk}{\bm{\kappa}} % orbital rotation vector
|
|
\newcommand{\bmg}{\bm{\gamma}} % orbital rotation vector
|
|
\newcommand{\bfx}{\bf{x}}
|
|
\newcommand{\bfr}{\bf{r}}
|
|
%%%%
|
|
|
|
\begin{document}
|
|
|
|
\title{Supplementary Material for ``Weight-dependent local density-functional approximations for ensembles''}
|
|
|
|
\author{Pierre-Fran\c{c}ois Loos}
|
|
\email{loos@irsamc.ups-tlse.fr}
|
|
\affiliation{\LCPQ}
|
|
\author{Emmanuel Fromager}
|
|
\email{fromagere@unistra.fr}
|
|
\affiliation{\LCQ}
|
|
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
\begin{abstract}
|
|
\end{abstract}
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
\maketitle
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
%\section{Ensemble Hartree--Fock method}
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
|
|
%The Hartree--Fock (HF) ensemble energy can be written
|
|
%as
|
|
%\beq\label{eq:eHF_ener}
|
|
%&&E^{\bw}_{\rm
|
|
%HF}({\bm\kappa})=
|
|
%\sum_{pq}\langle
|
|
%\varphi_p(\bmk)\vert\hat{h}\vert \varphi_q(\bmk)\rangle\sum^M_{K=0}w^{(K)}D^{(K)}_{pq}
|
|
%\nonumber\\
|
|
%&&+\dfrac{1}{2}\sum_{pqrs}\langle \varphi_p(\bmk)\varphi_q(\bmk)\vert\vert
|
|
%\varphi_r(\bmk)\varphi_s(\bmk)\rangle
|
|
%%\times
|
|
%\sum^M_{K=0}w^{(K)}D^{(K)}_{pr}D^{(K)}_{qs},
|
|
%\nonumber\\
|
|
%\eeq
|
|
%where the one- and antisymmetrized two-electron integrals read,
|
|
%\beq
|
|
%\langle
|
|
%\varphi_p({\bmk})\vert\hat{h}\vert
|
|
%\varphi_q({\bmk})\rangle=\int d{\bfx}\;
|
|
%\varphi_p({\bmk},{\bfx})\hat{h}\varphi_q({\bmk},{\bfx})
|
|
%\eeq
|
|
%with $\hat{h}\equiv-\frac{\nabla_{\bfr}^2}{2}+v_{\rm
|
|
%ext}(\bfr)$
|
|
%and
|
|
%\beq
|
|
%&&\langle \varphi_p(\bmk)\varphi_q(\bmk)\vert\vert
|
|
%\varphi_r(\bmk)\varphi_s(\bmk)\rangle=
|
|
%\nonumber\\
|
|
%&&\int d{\bfx}_1\int d{\bfx}_2\;
|
|
%\varphi_p({\bmk},{\bfx}_1)\varphi_q({\bmk},{\bfx}_2)\frac{1}{\vert
|
|
%{\bfr}_1-{\bfr}_2\vert}
|
|
%\nonumber
|
|
%\\
|
|
%&&\times\Big[\varphi_r({\bmk},{\bfx}_1)\varphi_s({\bmk},{\bfx}_2)
|
|
%-\varphi_s({\bmk},{\bfx}_1)\varphi_r({\bmk},{\bfx}_2)\Big]
|
|
%,
|
|
%\eeq
|
|
%respectively. Note that we use {\it real algebra} and the shorthand
|
|
%notation $\int d{\bfx}\equiv\int
|
|
%d{\bfr}\sum_{\sigma}$ for integration over space and
|
|
%summation over spin.
|
|
%% normalization condition
|
|
%%$\sum^M_{K=0}w^{(K)}=1$
|
|
%%.
|
|
%The antihermitian $\bmk\equiv\{\kappa_{pq}\}_{p>q}$ matrix which appears in the integrals controls
|
|
%the rotation of the spin-orbitals as follows,
|
|
%\beq\label{eq:orb_taylor_expansion}
|
|
%&&\varphi_p({\bmk},{\bfx})=\sum_q\left[e^{-{\bmk}}\right]_{qp}\varphi_q({\bfx})
|
|
%\nonumber\\
|
|
%&&=
|
|
%\varphi_p({\bfx})+\sum_{q<p}\kappa_{pq}\varphi_q({\bfx})-\sum_{q>p}\kappa_{qp}\varphi_q({\bfx})
|
|
%+\mathcal{O}\left({\bmk}^2\right).
|
|
%\eeq
|
|
%The (${\bmk}$-independent) one-electron reduced density matrices (1RDMs)
|
|
%in Eq.~(\ref{eq:eHF_ener}) are defined in the unrotated molecular
|
|
%spin-orbital basis for each (unrotated) determinant $\Phi^{(K)}$ belonging to
|
|
%the ensemble as follows: $D^{(K)}_{pr}=\delta_{pr}$ if $\varphi_p$ and
|
|
%$\varphi_r$ are both
|
|
%occupied in $\Phi^{(K)}$, otherwise $D^{(K)}_{pr}=0$. If the unrotated
|
|
%spin-orbitals are the minimizing ensemble HF ones, then the following
|
|
%stationarity condition is fulfilled,
|
|
%\beq\label{eq:station_cond}
|
|
%\left.\dfrac{\partial E^{\bw}_{\rm
|
|
%HF}({\bm\kappa})}{\partial \kappa_{lm}}
|
|
%\right|_{{\bmk}=0}=0,
|
|
%\eeq
|
|
%with $l>m$. Since, according to Eq.~(\ref{eq:orb_taylor_expansion}),
|
|
%\beq
|
|
%\left.\dfrac{\partial
|
|
%\varphi_p({\bmk},{\bfx})}{\partial \kappa_{lm}}
|
|
%\right|_{{\bmk}=0}=\delta_{lp}\varphi_m({\bfx})-\delta_{mp}\varphi_l({\bfx}),
|
|
%\eeq
|
|
%Eq.~(\ref{eq:station_cond}) can be written more explicitly as the
|
|
%following commutation relation,
|
|
%\iffalse%%%%%
|
|
%%%%%%% intermediate steps ... %%%%
|
|
%\beq
|
|
%\sum^M_{K=0}w^{(K)}\sum_qD^{(K)}_{mq}f^{(K)}_{lq}
|
|
%-
|
|
%\sum^M_{K=0}w^{(K)}\sum_qf^{(K)}_{mq}D^{(K)}_{lq}=0
|
|
%\eeq
|
|
%% original %%
|
|
%\iffalse%%%
|
|
%\beq
|
|
%&&
|
|
%2\sum_q\langle\varphi_m\vert\hat{h}\vert \varphi_q\rangle\sum^M_{K=0}w^{(K)}D^{(K)}_{lq}
|
|
%\nonumber\\
|
|
%&&-
|
|
%2\sum_q\langle\varphi_l\vert\hat{h}\vert \varphi_q\rangle\sum^M_{K=0}w^{(K)}D^{(K)}_{mq}
|
|
%\nonumber\\
|
|
%&&+2\sum_{qrs}\langle \varphi_m\varphi_q\vert\vert
|
|
%\varphi_r\varphi_s\rangle
|
|
%%\times
|
|
%\sum^M_{K=0}w^{(K)}D^{(K)}_{lr}D^{(K)}_{qs}
|
|
%\nonumber\\
|
|
%&&-2\sum_{qrs}\langle \varphi_l\varphi_q\vert\vert
|
|
%\varphi_r\varphi_s\rangle
|
|
%%\times
|
|
%\sum^M_{K=0}w^{(K)}D^{(K)}_{mr}D^{(K)}_{qs}
|
|
%\nonumber\\
|
|
%&&=0\eeq
|
|
%\fi%%%%
|
|
%%%%%%%%
|
|
%\beq
|
|
%&&
|
|
%2\sum_q\langle\varphi_m\vert\hat{h}\vert \varphi_q\rangle\sum^M_{K=0}w^{(K)}D^{(K)}_{lq}
|
|
%\nonumber\\
|
|
%&&-
|
|
%2\sum_q\langle\varphi_l\vert\hat{h}\vert \varphi_q\rangle\sum^M_{K=0}w^{(K)}D^{(K)}_{mq}
|
|
%\nonumber\\
|
|
%&&+2\sum_{qrs}\langle \varphi_m\varphi_r\vert\vert
|
|
%\varphi_q\varphi_s\rangle
|
|
%%\times
|
|
%\sum^M_{K=0}w^{(K)}D^{(K)}_{lq}D^{(K)}_{rs}
|
|
%\nonumber\\
|
|
%&&-2\sum_{qrs}\langle \varphi_l\varphi_r\vert\vert
|
|
%\varphi_q\varphi_s\rangle
|
|
%%\times
|
|
%\sum^M_{K=0}w^{(K)}D^{(K)}_{mq}D^{(K)}_{rs}
|
|
%\nonumber\\
|
|
%&&=0\eeq
|
|
%%%%%%%%%%%%%%%%%%%%%%
|
|
%%%%%%%%%%%%%%%%%%%%%%
|
|
%\fi%%%%
|
|
%%%%%%%%
|
|
%\beq\label{eq:stat_cond_commut_ind}
|
|
%\sum^M_{K=0}w^{(K)}\left[{\bm f}^{(K)},{\bm D}^{(K)}\right]=0,
|
|
%\eeq
|
|
%where the $K$th Fock matrix elements read
|
|
%\beq
|
|
%f^{(K)}_{mq}=\langle\varphi_m\vert\hat{h}\vert \varphi_q\rangle
|
|
%+
|
|
%\sum_{rs}\langle \varphi_m\varphi_r\vert\vert
|
|
%\varphi_q\varphi_s\rangle D^{(K)}_{rs}.
|
|
%\eeq
|
|
%In the minimizing ensemble HF spin-orbital basis, Eq.~(\ref{eq:stat_cond_commut_ind}) reads
|
|
%\beq
|
|
%\sum^M_{K=0}w^{(K)}\Big(\nu^{(K)}_m-\nu_l^{(K)}\Big)f^{(K)}_{lm}=0,
|
|
%\eeq
|
|
%where $\nu^{(K)}_m$ is the occupation of the spin-orbital $\varphi_m$ in the
|
|
%determinant $\Phi^{(K)}$.\\
|
|
%
|
|
%Note that, in more conventional ensemble calculations, the following HF
|
|
%energy expression is employed,
|
|
%\beq\label{eq:GI_ensHF_ener}
|
|
%&&\tilde{E}^{\bw}_{\rm
|
|
%HF}({\bm\kappa})=
|
|
%\sum_{pq}\langle
|
|
%\varphi_p(\bmk)\vert\hat{h}\vert \varphi_q(\bmk)\rangle D^{\bw}_{pq}
|
|
%\nonumber\\
|
|
%&&
|
|
%+\dfrac{1}{2}\sum_{pqrs}\langle \varphi_p(\bmk)\varphi_q(\bmk)\vert\vert
|
|
%\varphi_r(\bmk)\varphi_s(\bmk)\rangle
|
|
%%\times
|
|
%D^{\bw}_{pr}D^{\bw}_{qs},
|
|
%\eeq
|
|
%where ${\bm D}^{\bw}=\sum^M_{K=0}w^{(K)}{\bm D}^{(K)}$ is the ensemble
|
|
%1RDM. In this case, the stationarity condition simply reads
|
|
%\beq
|
|
%\left[{\bm f}^{\bw},{\bm D}^{\bw}\right]=0,
|
|
%\eeq
|
|
%where the ensemble Fock matrix elements are defined as follows,
|
|
%\beq
|
|
%f^{\bw}_{mq}=\langle\varphi_m\vert\hat{h}\vert \varphi_q\rangle
|
|
%+
|
|
%\sum_{rs}\langle \varphi_m\varphi_r\vert\vert
|
|
%\varphi_q\varphi_s\rangle D^{\bw}_{rs}.
|
|
%\eeq
|
|
%The major issue with the expression of the ensemble energy in
|
|
%Eq.~(\ref{eq:GI_ensHF_ener}) is the
|
|
%ghost-interaction error from which our expression (see
|
|
%Eq.~(\ref{eq:eHF_ener})) is free. Note also
|
|
%that, by construction, the ensemble energy in Eq.~(\ref{eq:GI_ensHF_ener}) is quadratic in the
|
|
%ensemble weights while ours, like the exact one, varies linearly with
|
|
%the weights.
|
|
%
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
%\section{Ensemble Hartree--Fock exchange and density-functional
|
|
%ghost-interaction correction}
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
%
|
|
%\beq
|
|
%F^{\bw}_{\rm HF}[n]&=&
|
|
%\underset{\hat{\gamma}^{{\bw}}\rightarrow n}{\rm min}\left\{{\rm
|
|
%Tr}\left[\hat{\gamma}^{{\bw}}\hat{T}\right]+W_{\rm
|
|
%HF}\left[{\bmg}^{\bw}\right]\right\}
|
|
%\nonumber\\
|
|
%&=&{\rm
|
|
%Tr}\left[\hat{\gamma}^{{\bw}}[n]\hat{T}\right]+W_{\rm
|
|
%HF}\left[{\bmg}^{\bw}[n]\right]
|
|
%\eeq
|
|
%where
|
|
%$\hat{\gamma}^{{\bw}}=\sum^M_{K=0}w^{(K)}\vert\Phi^{(K)}\rangle\langle\Phi^{(K)}\vert=\sum^M_{K=0}w^{(K)}\hat{\gamma}^{(K)}$ is an ensemble density matrix operator constructed
|
|
%from Slater determinants, the ensemble 1RDM elements are $\gamma_{pq}^{\bw}={\rm
|
|
%Tr}\left[\hat{\gamma}^{{\bw}}\hat{a}^\dagger_p\hat{a}_q\right]$,
|
|
%and $W_{\rm
|
|
%HF}\left[{\bmg}\right]=\frac{1}{2}\sum_{pqrs}\langle \varphi_p\varphi_q\vert\vert
|
|
%\varphi_r\varphi_s\rangle
|
|
%%\times
|
|
%\gamma_{pr}\gamma_{qs}$.\\
|
|
%
|
|
%In-principle-exact decomposition:
|
|
%
|
|
%\beq
|
|
%F^{\bw}[n]= F^{\bw}_{\rm HF}[n]+\overline{E}^{{\bw}}_{\rm
|
|
%Hx}[n]+\overline{E}^{{\bw}}_{\rm c}[n]
|
|
%\eeq
|
|
%
|
|
%The complementary ensemble Hx energy removes the ghost-interaction
|
|
%errors introduced in $W_{\rm
|
|
%HF}\left[{\bmg}^{\bw}[n]\right]$:
|
|
%\beq
|
|
%\overline{E}^{{\bw}}_{\rm
|
|
%Hx}[n]=\sum^M_{K=0}w^{(K)}W_{\rm
|
|
%HF}\left[{\bmg}^{(K)}[n]\right]
|
|
%-W_{\rm
|
|
%HF}\left[{\bmg}^{\bw}[n]\right],
|
|
%\eeq
|
|
%which gives in the canonical orbital basis
|
|
%\beq
|
|
%&&\overline{E}^{{\bw}}_{\rm
|
|
%Hx}[n]=
|
|
%\dfrac{1}{2}\sum_{pq}
|
|
%\langle \varphi^{{\bw}}_p[n]\varphi^{{\bw}}_q[n]\vert\vert
|
|
%\varphi^{{\bw}}_p[n]\varphi^{{\bw}}_q[n]\rangle
|
|
%\nonumber\\
|
|
%&&\times\left[\sum^M_{K=0}w^{(K)}\nu^{(K)}_p \left(\nu^{(K)}_q
|
|
%-\sum^M_{L=0}w^{(L)} \nu^{(L)}_q\right)\right]
|
|
%.\eeq
|
|
%\manu{I would guess that, in a uniform system, the GOK-DFT and our
|
|
%canonical orbitals are the same. This is nice since we can construct
|
|
%in a clean way density-functional approximations for both $\overline{E}^{{\bw}}_{\rm
|
|
%Hx}[n]$ and $E^{{\bw}}_{\rm c}[n]$ functionals. Am I right ?}
|
|
%
|
|
%Variational expression for the ensemble energy:
|
|
%\beq
|
|
%E^{{\bw}}=\underset{\hat{\gamma}^{{\bw}}}{\rm min}\Big\{
|
|
%&&{\rm
|
|
%Tr}\left[\hat{\gamma}^{{\bw}}\hat{T}\right]+W_{\rm
|
|
%HF}\left[{\bmg}^{\bw}\right]
|
|
%+
|
|
%\overline{E}^{{\bw}}_{\rm
|
|
%Hxc}\left[n_{\hat{\gamma}^{{\bw}}}\right]
|
|
%%+E^{{\bw}}_{\rm c}\left[n_{\hat{\gamma}^{{\bw}}}\right]
|
|
%\nonumber\\
|
|
%&&
|
|
%+\int d{\br}\;v_{\rm ext}({\bfr})n_{\hat{\gamma}^{{\bw}}}({\bfr})
|
|
%\Big\}
|
|
%\eeq
|
|
%
|
|
%Note that, if we use orbital rotations, the gradient of the DFT energy
|
|
%contributions can be expressed as follows,
|
|
%\beq
|
|
%\left.\dfrac{\partial
|
|
%\overline{E}^{{\bw}}_{\rm Hxc}\left[n^{{\bw}}({\bmk})\right]
|
|
%}{\partial \kappa_{lm}}
|
|
%\right|_{{\bmk}=0}=\int d{\br}\dfrac{\delta \overline{E}^{{\bw}}_{\rm
|
|
%Hxc}\left[n^{{\bw}}\right]}{\delta
|
|
%n({\br})}\left.\dfrac{\partial n^{{\bw}}({\bmk},{\br})}{\partial \kappa_{lm}}
|
|
%\right|_{{\bmk}=0},
|
|
%\eeq
|
|
%where
|
|
%\beq
|
|
%n^{{\bw}}({\bmk},{\br})=\sum_\sigma\sum_{pq}\varphi_p({\bmk},{\bfx})\varphi_q({\bmk},{\bfx})\gamma_{pq}^{\bw}
|
|
%\eeq
|
|
%thus leading to
|
|
%\beq
|
|
%&&\left.\dfrac{\partial
|
|
%\overline{E}^{{\bw}}_{\rm Hxc}\left[n^{{\bw}}({\bmk})\right]
|
|
%}{\partial \kappa_{lm}}
|
|
%\right|_{{\bmk}=0}=
|
|
%\sum_{pq}\gamma_{pq}^{\bw}
|
|
%\nonumber\\
|
|
%&&\times\left.\dfrac{\partial}
|
|
%{\partial \kappa_{lm}}
|
|
%\Big[\left\langle\varphi_p(\bmk)\middle\vert\hat{\overline{v}}^{{\bw}}_{\rm
|
|
%Hxc}
|
|
%\middle\vert \varphi_q(\bmk)\right\rangle
|
|
%\Big]
|
|
%\right|_{{\bmk}=0}.
|
|
%\eeq
|
|
%
|
|
%In conclusion, the minimizing canonical orbitals fulfill the following
|
|
%hybrid HF/GOK-DFT equation,
|
|
%\beq
|
|
%&&\left(-\frac{\nabla_{\bfr}^2}{2}+v_{\rm
|
|
%ext}({\bfr})+\hat{u}_{\rm HF}\left[\gamma^{\bw}\right]
|
|
%+\dfrac{\delta \overline{E}^{{\bw}}_{\rm Hxc}\left[n^{{\bw}}\right]}{\delta
|
|
%n({\br})}\right)\varphi^{{\bw}}_p({\bfx})
|
|
%\nonumber
|
|
%\\
|
|
%&&=\varepsilon^{{\bw}}_p\varphi^{{\bw}}_p({\bfx}).
|
|
%\eeq
|
|
%
|
|
%
|
|
%Since $\partial \gamma_{pq}^{\bw}/\partial
|
|
%w^{(I)}=\gamma_{pq}^{(I)}-\gamma_{pq}^{(0)}$, it comes
|
|
%
|
|
%\manu{just for me ...
|
|
%\beq
|
|
%&&+\dfrac{1}{2}
|
|
%\sum_{pqrs}\langle \varphi_p\varphi_q\vert\vert
|
|
%\varphi_r\varphi_s\rangle
|
|
%%\times
|
|
%\left(\gamma_{pr}^{(I)}-\gamma_{pr}^{(0)}\right)\gamma^{\bw}_{qs}
|
|
%\nonumber\\
|
|
%&&+\dfrac{1}{2}\sum_{pqrs}\langle \varphi_q\varphi_p\vert\vert
|
|
%\varphi_s\varphi_r\rangle
|
|
%%\times
|
|
% \gamma^{\bw}_{pr}\left(\gamma_{qs}^{(I)}-\gamma_{qs}^{(0)}\right)
|
|
%\nonumber\\
|
|
%&&=
|
|
%\sum_{pqrs}\langle \varphi_p\varphi_q\vert\vert
|
|
%\varphi_r\varphi_s\rangle
|
|
%%\times
|
|
%\left(\gamma_{pr}^{(I)}-\gamma_{pr}^{(0)}\right)\gamma^{\bw}_{qs}
|
|
%\nonumber\\
|
|
%&&=
|
|
%\sum_{pr}\left[\hat{u}_{\rm HF}\left[\gamma^{\bw}\right]\right]_{pr}\left(\gamma_{pr}^{(I)}-\gamma_{pr}^{(0)}\right)
|
|
%\nonumber\\
|
|
%&&=
|
|
%\sum_p\left[\hat{u}_{\rm
|
|
%HF}\left[\gamma^{\bw}\right]\right]_{pp}\left(\nu_p^{(I)}-\nu_p^{(0)}\right)
|
|
%\eeq
|
|
%}
|
|
%
|
|
%\beq
|
|
%\dfrac{dE^{\bw}}{dw^{(I)}}=\sum_p\varepsilon^{{\bw}}_p\left(\nu_p^{(I)}-\nu_p^{(0)}\right)+\left.\dfrac{\partial \overline{E}^{{\bw}}_{\rm
|
|
%Hxc}\left[n\right]}{\partial w^{(I)}}\right|_{n=n^{{\bw}}}.
|
|
%\eeq
|
|
%
|
|
%LZ shift in this context: $\varepsilon^{{\bw}}_p\rightarrow
|
|
%\overline{\varepsilon}^{{\bw}}_p=\varepsilon^{{\bw}}_p+\overline{\Delta}_{\rm
|
|
%LZ}^{{\bw}}$ where
|
|
%
|
|
%\beq
|
|
%N\overline{\Delta}_{\rm
|
|
%LZ}^{{\bw}}&=&\overline{E}^{{\bw}}_{\rm Hxc}\left[n^{{\bw}}\right]
|
|
%-\int d{\br}\dfrac{\delta \overline{E}^{{\bw}}_{\rm
|
|
%Hxc}\left[n^{{\bw}}\right]}{\delta
|
|
%n({\br})}n^{{\bw}}({\bfr})
|
|
%\nonumber\\
|
|
%&&
|
|
%-W_{\rm
|
|
%HF}\left[{\bmg}^{\bw}\right]
|
|
%\eeq
|
|
%
|
|
%such that
|
|
%\beq
|
|
%E^{{\bw}}=\sum^M_{K=0}w^{(K)}\sum_p\nu_p^{(K)}\overline{\varepsilon}^{{\bw}}_p.
|
|
%\eeq
|
|
%
|
|
%Thus we conclude that individual energies can be expressed in principle
|
|
%exactly as follows,
|
|
%
|
|
%\beq
|
|
%E^{(K)}=\sum_p\nu_p^{(K)}\overline{\varepsilon}^{{\bw}}_p+\sum^M_{I>0}\left(\delta_{IK}-w^{(I)}\right)\left.\dfrac{\partial \overline{E}^{{\bw}}_{\rm
|
|
%Hxc}\left[n\right]}{\partial w^{(I)}}\right|_{n=n^{{\bw}}}.
|
|
%\eeq
|
|
%%%%%%%%%%%%%%%
|
|
%\iffalse%%%%
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
%\section{Generalized GOK-DFT}
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
%
|
|
%The energy to be minimized in a generalized GOK-DFT approach can be
|
|
%written as
|
|
%\beq
|
|
%\eeq
|
|
%\fi%%%%%%%%%%
|
|
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
\section*{Construction of the density-functional approximations}
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
|
|
The density-functional approximations designed in this manuscript are based on highly-accurate energies for the ground state ($I=0$), the first singly-excited state ($I=1$), and the first doubly-excited state ($I=2$) of the (spin-polarized) two-electron ringium system.
|
|
We refer the interested reader to Refs.~\onlinecite{Loos_2012, Loos_2013a, Loos_2014b} for more details about this paradigm.
|
|
|
|
The reduced (i.e.~per electron) HF energy for these three states is:
|
|
\begin{subequations}
|
|
\begin{align}
|
|
\e{HF}{(0)}(n) & = \frac{\pi^2}{8} n^2 + n,
|
|
\\
|
|
\e{HF}{(1)}(n) & = \frac{\pi^2}{2}n^2 + \frac{4}{3} n,
|
|
\\
|
|
\e{HF}{(2)}(n) & = \frac{9\pi^2}{8}n^2 + \frac{23}{15} n.
|
|
\end{align}
|
|
\end{subequations}
|
|
All these states have the same (uniform) density $n = 2/(2\pi R)$ where $R$ is the radius of the ring on which the electrons are confined.
|
|
|
|
The total energy of the ground and doubly-excited states are given by the two lowest eigenvalues of the Hamiltonian $\bH$ with elements
|
|
\begin{equation}
|
|
\begin{split}
|
|
H_{ij}
|
|
& = \int_0^\pi \qty[ \frac{\psi_i(\omega)}{R} \frac{\psi_j(\omega)}{R} + \frac{\psi_i(\omega)\psi_j(\omega)}{2R\sin(\omega/2)} ] d\omega
|
|
\\
|
|
& = \frac{\sqrt{\pi}}{2 R} \qty[ \frac{\Gamma\qty(\frac{i+j}{2})}{\Gamma\qty(\frac{i+j+1}{2})} + \frac{ij}{4R} \frac{\Gamma\qty(\frac{i+j-1}{2})}{\Gamma\qty(\frac{i+j+2}{2})} ],
|
|
\end{split}
|
|
\end{equation}
|
|
where $\omega = \theta_1 - \theta_2$ is the interelectronic angle, $\Gamma(x)$ is the Gamma function, \cite{NISTbook} and
|
|
\begin{equation}
|
|
\psi_i(\omega) = \sin(\omega/2) \sin^{i-1}(\omega/2), \quad i=1,\ldots,M
|
|
\end{equation}
|
|
are (non-orthogonal) explicitly-correlated basis functions with overlap matrix elements
|
|
\begin{equation}
|
|
S_{ij}
|
|
= \int_0^\pi \psi_i(\omega)\psi_j(\omega) d\omega
|
|
= \sqrt{\pi} \frac{\Gamma\qty(\frac{i+j+1}{2})}{\Gamma\qty(\frac{i+j+2}{2})}.
|
|
\end{equation}
|
|
Thanks to this explicitly-correlated basis, the convergence rate of the energy is exponential with respect to $M$.
|
|
Therefore, high accuracy is reached with a very small number of basis functions.
|
|
Here, we typically use $M=10$.
|
|
For the singly-excited state, one has to modify the basis functions as
|
|
\begin{equation}
|
|
\psi_i(\omega) = \cos(\omega/2) \sin^{i-1} (\omega/2),
|
|
\end{equation}
|
|
and its energy is obtained by the lowest root of the Hamiltonian in this basis, and the matrix elements reads
|
|
\begin{align}
|
|
H_{ij} & = \frac{\sqrt{\pi}}{4 R} \qty[ \frac{\Gamma\qty(\frac{i+j}{2})}{\Gamma\qty(\frac{i+j+1}{3})} + \frac{3ij+i+j-1}{4R} \frac{\Gamma\qty(\frac{i+j-1}{2})}{\Gamma\qty(\frac{i+j+4}{2})} ],
|
|
\\
|
|
S_{ij} & = \frac{\sqrt{\pi}}{2} \frac{\Gamma\qty(\frac{i+j+1}{2})}{\Gamma\qty(\frac{i+j+4}{2})}.
|
|
\end{align}
|
|
The numerical values of the correlation energy for various $R$ are reported in Table \ref{tab:Ref} for the three states of interest.
|
|
|
|
%%% FIG 1 %%%
|
|
\begin{figure}
|
|
\includegraphics[width=\linewidth]{Ec}
|
|
\caption{
|
|
Reduced (i.e., per electron) correlation energy $\e{c}{(I)}$ [see Eq.~\eqref{eq:ec}] as a function of $R = 1/(\pi n)$ for the ground state ($I=0$), the first singly-excited state ($I=1$), and the first doubly-excited state ($I=2$) of the (spin-polarized) two-electron ringium system.
|
|
The data gathered in Table \ref{tab:Ref} are also reported.
|
|
}
|
|
\label{fig:Ec}
|
|
\end{figure}
|
|
%%% %%% %%%
|
|
|
|
%%% TABLE I %%%
|
|
\begin{turnpage}
|
|
\begin{squeezetable}
|
|
\begin{table*}
|
|
\caption{
|
|
\label{tab:Ref}
|
|
$-\e{c}{(I)}$ as a function of the radius of the ring $R$ for the ground state ($I=0$), the first singly-excited state ($I=1$), and the first doubly-excited state ($I=2$) of the (spin-polarized) two-electron ringium system.
|
|
}
|
|
\begin{ruledtabular}
|
|
\begin{tabular}{lcddddddddddd}
|
|
State & $I$ & \mc{11}{c}{Ring's radius $R = 1/(\pi n)$} \\
|
|
\cline{3-13}
|
|
& & \tabc{$0$} & \tabc{$1/10$} & \tabc{$1/5$} & \tabc{$1/2$} & \tabc{$1$} & \tabc{$2$} & \tabc{$5$} & \tabc{$10$} & \tabc{$20$} & \tabc{$50$} & \tabc{$100$} \\
|
|
\hline
|
|
Ground state & $0$ & 0.013708 & 0.012859 & 0.012525 & 0.011620 & 0.010374 & 0.008558 & 0.005673 & 0.003697 & 0.002226 & 0.001046 & 0.000567 \\
|
|
Singly-excited state & $1$ & 0.0238184 & 0.023392 & 0.022979 & 0.021817 & 0.020109 & 0.017371 & 0.012359 & 0.008436 & 0.005257 & 0.002546 & 0.001399 \\
|
|
Doubly-excited state & $2$ & 0.018715 & 0.018653 & 0.018576 & 0.018300 & 0.017743 & 0.016491 & 0.013145 & 0.009670 & 0.006365 & 0.003231 & 0.001816 \\
|
|
\end{tabular}
|
|
\end{ruledtabular}
|
|
\end{table*}
|
|
\end{squeezetable}
|
|
\end{turnpage}
|
|
|
|
Based on these highly-accurate calculations, one can write down, for each state, an accurate analytical expression of the reduced correlation energy \cite{Loos_2013a, Loos_2014a} via the following Pad\'e approximant
|
|
\begin{equation}
|
|
\label{eq:ec}
|
|
\e{c}{(I)}(n) = \frac{a_1^{(I)}\,n}{n + a_2^{(I)} \sqrt{n} + a_3^{(I)}},
|
|
\end{equation}
|
|
where $a_2^{(I)}$ and $a_3^{(I)}$ are state-specific fitting parameters, which are provided in Table I of the manuscript.
|
|
The value of $a_1^{(I)}$ is obtained via the exact high-density expansion of the correlation energy. \cite{Loos_2013a, Loos_2014a}
|
|
Equation \eqref{eq:ec} is depicted in Fig.~\ref{fig:Ec} for each state alongside the data gathered in Table \ref{tab:Ref}.
|
|
|
|
|
|
%%% FIG 2 %%%
|
|
\begin{figure*}
|
|
\includegraphics[height=0.325\linewidth]{EvsL_2}
|
|
\includegraphics[height=0.325\linewidth]{EvsL_3}
|
|
\includegraphics[height=0.325\linewidth]{EvsL_4}
|
|
\includegraphics[height=0.325\linewidth]{EvsL_5}
|
|
\includegraphics[height=0.325\linewidth]{EvsL_6}
|
|
\includegraphics[height=0.325\linewidth]{EvsL_7}
|
|
\caption{
|
|
Error with respect to FCI in single and double excitation energies of $N$-boxium as a function of the box length $L$ for various methods.
|
|
See main text for additional details.}
|
|
\label{fig:EvsL}
|
|
\end{figure*}
|
|
%%% %%% %%%
|
|
|
|
%%% FIG 3 %%%
|
|
%\begin{figure*}
|
|
% \includegraphics[height=0.325\linewidth]{EvsN_0125}
|
|
% \includegraphics[height=0.325\linewidth]{EvsN_025}
|
|
% \includegraphics[height=0.325\linewidth]{EvsN_05}
|
|
% \includegraphics[height=0.325\linewidth]{EvsN_1}
|
|
% \includegraphics[height=0.325\linewidth]{EvsN_2}
|
|
% \includegraphics[height=0.325\linewidth]{EvsN_4}
|
|
% \includegraphics[height=0.325\linewidth]{EvsN_8}
|
|
% \caption{
|
|
% Error with respect to FCI in single and double excitation energies of $N$-boxium as a function of the number of electrons $N$ for various methods and box length $L$.
|
|
% }
|
|
% \label{fig:EvsL}
|
|
%\end{figure*}
|
|
%%% %%% %%%
|
|
|
|
%%% TABLE II %%%
|
|
\begin{table*}
|
|
\caption{
|
|
\label{tab:OptGap_2}
|
|
Deviation from the FCI quantities (in hartree) of the individual energies, $\E{(I)}$, and the corresponding excitation energies, $\Ex{(I)}$, for the ground ($I=0$), singly-excited ($I=1$) and doubly-excited ($I=2$) states of 2-boxium (i.e.,~$\Nel = 2$ electrons in a box of length $L$).
|
|
The values of the derivative discontinuity $\DD{c}{(I)}$ are also reported.
|
|
(DNC = KS calculation does not converge.)
|
|
}
|
|
\begin{ruledtabular}
|
|
\begin{tabular}{lclddddddd}
|
|
& & & \mc{7}{c}{2-boxium with a box of length $L$} \\
|
|
\cline{4-10}
|
|
Method & $\bw$ & State & \pi/8 & \pi/4 & \pi/2 & \pi & 2\pi & 4\pi & 8\pi \\
|
|
\hline
|
|
FCI & & $\E{(0)}$ & 168.1946 & 44.0662 & 12.0035 & 3.4747 & 1.0896 & 0.3719 & 0.1367 \\
|
|
& & $\E{(1)}$ & 330.2471 & 85.0890 & 22.5112 & 6.2247 & 1.8355 & 0.5845 & 0.2006 \\
|
|
& & $\E{(2)}$ & 809.9972 & 204.9840 & 52.4777 & 13.7252 & 3.7248 & 1.0696 & 0.3300 \\
|
|
& & $\Ex{(1)}$ & 162.0525 & 41.0228 & 10.5078 & 2.7500 & 0.7458 & 0.2125 & 0.0639 \\
|
|
& & $\Ex{(2)}$ & 641.8026 & 160.9177 & 40.4743 & 10.2505 & 2.6352 & 0.6977 & 0.1933 \\
|
|
\hline
|
|
& & & \mc{7}{c}{Deviation from FCI} \\
|
|
\hline
|
|
CIS & & $\Ex{(1)}$ & 0.0104 & 0.0102 & 0.0099 & 0.0092 & 0.0077 & 0.0051 & 0.0022 \\
|
|
\\
|
|
TDHF & & $\Ex{(1)}$ & 0.0019 & 0.0021 & 0.0023 & 0.0027 & 0.0029 & 0.0023 & 0.0011 \\
|
|
\\
|
|
TDA-TDLDA& & $\Ex{(1)}$ & 0.0099 & 0.0088 & 0.0058 & -0.0041 & -0.0316 & -0.0467 & \tabc{DNC} \\
|
|
\\
|
|
TDLDA & & $\Ex{(1)}$ & 0.0015 & 0.0006 & -0.0018 & -0.0106 & -0.0370 & -0.0518 & \tabc{DNC} \\
|
|
\\
|
|
% eHF & $(0,0)$ & $\E{(0)}$ & 0.0137 & 0.0131 & 0.0119 & 0.0098 & 0.0067 & 0.0033 & 0.0010 \\
|
|
% & & $\E{(1)}$ & 0.0685 & 0.0678 & 0.0664 & 0.0637 & 0.0588 & 0.0500 & 0.0372 \\
|
|
% & & $\E{(2)}$ & -0.0052 & -0.0055 & -0.0058 & -0.0055 & -0.0022 & 0.0086 & 0.0271 \\
|
|
% & & $\Ex{(1)}$ & 0.1082 & 0.1069 & 0.1044 & 0.0998 & 0.0911 & 0.0736 & \tabc{DNC} \\
|
|
% & & $\Ex{(2)}$ & 0.0345 & 0.0336 & 0.0322 & 0.0306 & 0.0302 & 0.0321 & \tabc{DNC} \\
|
|
% \\
|
|
% eHF & $(1/3,1/3)$ & $\E{(0)}$ & 0.0565 & 0.0557 & 0.0540 & 0.0509 & 0.0455 & 0.0372 & 0.0269 \\
|
|
% & & $\E{(1)}$ & 0.0560 & 0.0553 & 0.0538 & 0.0510 & 0.0461 & 0.0384 & 0.0287 \\
|
|
% & & $\E{(2)}$ & 0.0371 & 0.0366 & 0.0357 & 0.0342 & 0.0316 & 0.0277 & 0.0223 \\
|
|
% & & $\Ex{(1)}$ & 0.0529 & 0.0517 & 0.0494 & 0.0456 & 0.0418 & 0.0409 & \tabc{DNC} \\
|
|
% & & $\Ex{(2)}$ & 0.0340 & 0.0330 & 0.0314 & 0.0288 & 0.0274 & 0.0303 & \tabc{DNC} \\
|
|
% \\
|
|
KS-eLDA & $(0,0)$ & $\E{(0)}$ & -0.0397 & -0.0391 & -0.0380 & -0.0361 & -0.0323 & -0.0236 & \tabc{DNC}\\
|
|
& & $\E{(1)}$ & 0.0215 & 0.0213 & 0.0210 & 0.0200 & 0.0159 & 0.0102 & \tabc{DNC}\\
|
|
& & $\E{(2)}$ & -0.0426 & -0.0425 & -0.0419 & -0.0387 & -0.0250 & -0.0045 & \tabc{DNC}\\
|
|
& & $\Ex{(1)}$ & 0.0612 & 0.0604 & 0.0590 & 0.0561 & 0.0483 & 0.0337 & \tabc{DNC}\\
|
|
& & $\Ex{(2)}$ & -0.0029 & -0.0034 & -0.0039 & -0.0025 & 0.0074 & 0.0191 & \tabc{DNC}\\
|
|
& & $\DD{c}{(0)}$ & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & \tabc{DNC} \\
|
|
& & $\DD{c}{(1)}$ & 0.0064 & 0.0056 & 0.0043 & 0.0022 & -0.0007 & -0.0037 & \tabc{DNC}\\
|
|
& & $\DD{c}{(2)}$ & 0.0159 & 0.0147 & 0.0126 & 0.0093 & 0.0046 & -0.0009 & \tabc{DNC}\\
|
|
\\
|
|
KS-eLDA & $(1/3,1/3)$ & $\E{(0)}$ & 0.0031 & 0.0036 & 0.0044 & 0.0054 & 0.0042 & -0.0025 & \tabc{DNC}\\
|
|
& & $\E{(1)}$ & 0.0090 & 0.0087 & 0.0083 & 0.0076 & 0.0070 & 0.0071 & \tabc{DNC}\\
|
|
& & $\E{(2)}$ & -0.0005 & -0.0009 & -0.0015 & -0.0023 & -0.0030 & -0.0026 & \tabc{DNC}\\
|
|
& & $\Ex{(1)}$ & 0.0058 & 0.0052 & 0.0039 & 0.0022 & 0.0028 & 0.0096 & \tabc{DNC}\\
|
|
& & $\Ex{(2)}$ & -0.0036 & -0.0045 & -0.0058 & -0.0077 & -0.0072 & 0.0000 & \tabc{DNC}\\
|
|
& & $\DD{c}{(0)}$ & -0.0074 & -0.0067 & -0.0055 & -0.0036 & -0.0010 & 0.0019 & \tabc{DNC}\\
|
|
& & $\DD{c}{(1)}$ & -0.0010 & -0.0011 & -0.0014 & -0.0017 & -0.0021 & -0.0022 & \tabc{DNC}\\
|
|
& & $\DD{c}{(2)}$ & 0.0084 & 0.0079 & 0.0069 & 0.0053 & 0.0031 & 0.0003 & \tabc{DNC}\\
|
|
\end{tabular}
|
|
\end{ruledtabular}
|
|
\end{table*}
|
|
|
|
%%% TABLE III %%%
|
|
\begin{table*}
|
|
\caption{
|
|
\label{tab:OptGap_3}
|
|
Deviation from the FCI quantities (in hartree) of the individual energies, $\E{(I)}$, and the corresponding excitation energies, $\Ex{(I)}$, for the ground ($I=0$), singly-excited ($I=1$) and doubly-excited ($I=2$) states of 3-boxium (i.e.,~$\Nel = 3$ electrons in a box of length $L$).
|
|
The values of the derivative discontinuity $\DD{c}{(I)}$ are also reported.
|
|
}
|
|
\begin{ruledtabular}
|
|
\begin{tabular}{lclddddddd}
|
|
& & & \mc{7}{c}{3-boxium with a box of length $L$} \\
|
|
\cline{4-10}
|
|
Method & $\bw$ & State & \pi/8 & \pi/4 & \pi/2 & \pi & 2\pi & 4\pi & 8\pi \\
|
|
\hline
|
|
FCI & & $\E{(0)}$ & 475.6891 & 125.7776 & 34.8248 & 10.3536 & 3.3766 & 1.2126 & 0.4721 \\
|
|
& & $\E{(1)}$ & 702.8330 & 183.3370 & 49.5922 & 14.2255 & 4.4269 & 1.5105 & 0.5606 \\
|
|
& & $\E{(2)}$ & 1379.3128 & 353.5967 & 92.7398 & 25.3135 & 7.3546 & 2.3203 & 0.7990 \\
|
|
& & $\Ex{(1)}$ & 227.1438 & 57.5594 & 14.7674 & 3.8720 & 1.0504 & 0.2979 & 0.0885 \\
|
|
& & $\Ex{(2)}$ & 903.6236 & 227.8191 & 57.9150 & 14.9599 & 3.9780 & 1.1077 & 0.3269 \\
|
|
\hline
|
|
& & & \mc{7}{c}{Deviation from FCI} \\
|
|
\hline
|
|
% CIS & & $\Ex{(1)}$ & 0.0163 & 0.0161 & 0.0157 & 0.0149 & 0.0133 & 0.0102 & 0.0057 \\
|
|
% \\
|
|
% TDHF & & $\Ex{(1)}$ & 0.0013 & 0.0013 & 0.0014 & 0.0014 & 0.0013 & 0.0010 & 0.0007 \\
|
|
\\
|
|
TDA-TDLDA& & $\Ex{(1)}$ & 0.0162 & 0.0157 & 0.0146 & 0.0110 & -0.0049 & -0.0344 & -0.0378 \\
|
|
\\
|
|
TDLDA & & $\Ex{(1)}$ & 0.0262 & 0.0264 & 0.0264 & 0.0269 & 0.0273 & 0.0206 & -0.0116 \\
|
|
\\
|
|
% eHF & $(0,0)$ & $\E{(0)}$ & 0.0327 & 0.0316 & 0.0297 & 0.0261 & 0.0203 & 0.0126 & 0.0055 \\
|
|
% & & $\E{(1)}$ & 0.1051 & 0.1038 & 0.1013 & 0.0965 & 0.0875 & 0.0708 & 0.0440 \\
|
|
% & & $\E{(2)}$ & 0.0841 & 0.0829 & 0.0806 & 0.0766 & 0.0704 & 0.0625 & 0.0517 \\
|
|
% & & $\Ex{(1)}$ & 0.1532 & 0.1516 & 0.1486 & 0.1428 & 0.1320 & 0.1095 & 0.0698 \\
|
|
% & & $\Ex{(2)}$ & 0.1322 & 0.1307 & 0.1279 & 0.1229 & 0.1149 & 0.1012 & 0.0774 \\
|
|
% \\
|
|
% eHF & $(1/3,1/3)$ & $\E{(0)}$ & 0.0886 & 0.0874 & 0.0850 & 0.0805 & 0.0725 & 0.0593 & 0.0404 \\
|
|
% & & $\E{(1)}$ & 0.0880 & 0.0865 & 0.0836 & 0.0783 & 0.0692 & 0.0552 & 0.0374 \\
|
|
% & & $\E{(2)}$ & 0.1209 & 0.1200 & 0.1182 & 0.1147 & 0.1078 & 0.0953 & 0.0747 \\
|
|
% & & $\Ex{(1)}$ & 0.0802 & 0.0785 & 0.0754 & 0.0699 & 0.0611 & 0.0529 & 0.0397 \\
|
|
% & & $\Ex{(2)}$ & 0.1131 & 0.1120 & 0.1100 & 0.1062 & 0.0997 & 0.0930 & 0.0770 \\
|
|
% \\
|
|
KS-eLDA & $(0,0)$ & $\E{(0)}$ & -0.0481 & -0.0478 & -0.0473 & -0.0463 & -0.0446 & -0.0387 & -0.0257 \\
|
|
& & $\E{(1)}$ & 0.0343 & 0.0336 & 0.0321 & 0.0292 & 0.0220 & 0.0084 & 0.0008 \\
|
|
& & $\E{(2)}$ & 0.0277 & 0.0267 & 0.0247 & 0.0216 & 0.0187 & 0.0208 & 0.0209 \\
|
|
& & $\Ex{(1)}$ & 0.0824 & 0.0814 & 0.0794 & 0.0755 & 0.0666 & 0.0471 & 0.0266 \\
|
|
& & $\Ex{(2)}$ & 0.0759 & 0.0745 & 0.0720 & 0.0679 & 0.0633 & 0.0595 & 0.0467 \\
|
|
& & $\DD{c}{(0)}$ & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\
|
|
& & $\DD{c}{(1)}$ & 0.0100 & 0.0092 & 0.0077 & 0.0051 & 0.0012 & -0.0034 & -0.0072 \\
|
|
& & $\DD{c}{(2)}$ & 0.0244 & 0.0231 & 0.0208 & 0.0168 & 0.0108 & 0.0029 & -0.0050 \\
|
|
\\
|
|
KS-eLDA & $(1/3,1/3)$ & $\E{(0)}$ & 0.0078 & 0.0080 & 0.0082 & 0.0085 & 0.0081 & 0.0024 & -0.0022 \\
|
|
& & $\E{(1)}$ & 0.0172 & 0.0162 & 0.0144 & 0.0112 & 0.0064 & 0.0019 & 0.0004 \\
|
|
& & $\E{(2)}$ & 0.0645 & 0.0636 & 0.0621 & 0.0590 & 0.0530 & 0.0420 & 0.0300 \\
|
|
& & $\Ex{(1)}$ & 0.0094 & 0.0083 & 0.0062 & 0.0027 & -0.0018 & -0.0004 & 0.0026 \\
|
|
& & $\Ex{(2)}$ & 0.0567 & 0.0557 & 0.0539 & 0.0506 & 0.0449 & 0.0397 & 0.0323 \\
|
|
& & $\DD{c}{(0)}$ & -0.0115 & -0.0107 & -0.0094 & -0.0072 & -0.0038 & 0.0005 & 0.0045 \\
|
|
& & $\DD{c}{(1)}$ & -0.0015 & -0.0016 & -0.0018 & -0.0022 & -0.0028 & -0.0033 & -0.0032 \\
|
|
& & $\DD{c}{(2)}$ & 0.0129 & 0.0123 & 0.0113 & 0.0094 & 0.0066 & 0.0028 & -0.0013 \\
|
|
\end{tabular}
|
|
\end{ruledtabular}
|
|
\end{table*}
|
|
|
|
|
|
%%% TABLE IV %%%
|
|
\begin{table*}
|
|
\caption{
|
|
\label{tab:OptGap_4}
|
|
Deviation from the FCI quantities (in hartree) of the individual energies, $\E{(I)}$, and the corresponding excitation energies, $\Ex{(I)}$, for the ground ($I=0$), singly-excited ($I=1$) and doubly-excited ($I=2$) states of 4-boxium (i.e.,~$\Nel = 4$ electrons in a box of length $L$).
|
|
The values of the derivative discontinuity $\DD{c}{(I)}$ are also reported.
|
|
}
|
|
\begin{ruledtabular}
|
|
\begin{tabular}{lclddddddd}
|
|
& & & \mc{7}{c}{4-boxium with a box of length $L$} \\
|
|
\cline{4-10}
|
|
Method & $\bw$ & State & \pi/8 & \pi/4 & \pi/2 & \pi & 2\pi & 4\pi & 8\pi \\
|
|
\hline
|
|
FCI & & $\E{(0)}$ & 1020.3778 & 270.0849 & 74.9426 & 22.3790 & 7.3595 & 2.6798 & 1.0633 \\
|
|
& & $\E{(1)}$ & 1312.2776 & 344.0184 & 93.8936 & 27.3398 & 8.7021 & 3.0600 & 1.1764 \\
|
|
& & $\E{(2)}$ & 2183.4399 & 563.5949 & 149.6753 & 41.7213 & 12.5052 & 4.1033 & 1.4749 \\
|
|
& & $\Ex{(1)}$ & 291.8998 & 73.9335 & 18.9510 & 4.9608 & 1.3426 & 0.3802 & 0.1131 \\
|
|
& & $\Ex{(2)}$ & 1163.0621 & 293.5099 & 74.7326 & 19.3423 & 5.1457 & 1.4235 & 0.4116 \\
|
|
\hline
|
|
& & & \mc{7}{c}{Deviation from FCI} \\
|
|
\hline
|
|
% CIS & & $\Ex{(1)}$ & 0.0203 & 0.0202 & 0.0200 & 0.0195 & 0.0187 & 0.0167 & 0.0116 \\
|
|
% \\
|
|
% TDHF & & $\Ex{(1)}$ & 0.0008 & 0.0008 & 0.0009 & 0.0009 & 0.0008 & 0.0008 & 0.0007 \\
|
|
% \\
|
|
TDA-TDLDA& & $\Ex{(1)}$ & 0.0203 & 0.0201 & 0.0195 & 0.0181 & 0.0106 & -0.0178 & -0.0369 \\
|
|
\\
|
|
TDLDA & & $\Ex{(1)}$ & 0.0008 & 0.0007 & 0.0004 & -0.0006 & -0.0074 & -0.0360 & -0.0653 \\
|
|
\\
|
|
% eHF & $(0,0)$ & $\E{(0)}$ & 0.0541 & 0.0529 & 0.0504 & 0.0459 & 0.0381 & 0.0266 & 0.0139 \\
|
|
% & & $\E{(1)}$ & 0.1359 & 0.1346 & 0.1320 & 0.1271 & 0.1186 & 0.1041 & 0.0804 \\
|
|
% & & $\E{(2)}$ & 0.1393 & 0.1374 & 0.1335 & 0.1260 & 0.1126 & 0.0910 & 0.0604 \\
|
|
% & & $\Ex{(1)}$ & 0.1900 & 0.1885 & 0.1857 & 0.1806 & 0.1715 & 0.1545 & 0.1189 \\
|
|
% & & $\Ex{(2)}$ & 0.1934 & 0.1913 & 0.1872 & 0.1795 & 0.1656 & 0.1414 & 0.0990 \\
|
|
% \\
|
|
% eHF & $(1/3,1/3)$ & $\E{(0)}$ & 0.1167 & 0.1153 & 0.1125 & 0.1073 & 0.0977 & 0.0818 & 0.0592 \\
|
|
% & & $\E{(1)}$ & 0.1110 & 0.1093 & 0.1059 & 0.0996 & 0.0888 & 0.0723 & 0.0517 \\
|
|
% & & $\E{(2)}$ & 0.1688 & 0.1669 & 0.1630 & 0.1553 & 0.1400 & 0.1118 & 0.0705 \\
|
|
% & & $\Ex{(1)}$ & 0.1025 & 0.1008 & 0.0975 & 0.0915 & 0.0816 & 0.0702 & 0.0532 \\
|
|
% & & $\Ex{(2)}$ & 0.1603 & 0.1584 & 0.1546 & 0.1471 & 0.1328 & 0.1097 & 0.0720 \\
|
|
% \\
|
|
KS-eLDA & $(0,0)$ & $\E{(0)}$ & -0.0541 & -0.0539 & -0.0537 & -0.0534 & -0.0529 & -0.0504 & -0.0386 \\
|
|
& & $\E{(1)}$ & 0.0413 & 0.0406 & 0.0390 & 0.0362 & 0.0304 & 0.0159 & 0.0008 \\
|
|
& & $\E{(2)}$ & 0.0642 & 0.0622 & 0.0586 & 0.0517 & 0.0399 & 0.0254 & 0.0149 \\
|
|
& & $\Ex{(1)}$ & 0.0954 & 0.0945 & 0.0927 & 0.0896 & 0.0833 & 0.0663 & 0.0394 \\
|
|
& & $\Ex{(2)}$ & 0.1182 & 0.1162 & 0.1123 & 0.1051 & 0.0928 & 0.0758 & 0.0534 \\
|
|
& & $\DD{c}{(0)}$ & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\
|
|
& & $\DD{c}{(1)}$ & 0.0136 & 0.0127 & 0.0111 & 0.0083 & 0.0038 & -0.0022 & -0.0080 \\
|
|
& & $\DD{c}{(2)}$ & 0.0330 & 0.0316 & 0.0291 & 0.0248 & 0.0178 & 0.0080 & -0.0028 \\
|
|
\\
|
|
KS-eLDA & $(1/3,1/3)$ & $\E{(0)}$ & 0.0085 & 0.0085 & 0.0084 & 0.0082 & 0.0072 & 0.0021 & -0.0015 \\
|
|
& & $\E{(1)}$ & 0.0164 & 0.0152 & 0.0129 & 0.0087 & 0.0020 & -0.0050 & -0.0044 \\
|
|
& & $\E{(2)}$ & 0.0936 & 0.0917 & 0.0880 & 0.0807 & 0.0664 & 0.0434 & 0.0300 \\
|
|
& & $\Ex{(1)}$ & 0.0079 & 0.0067 & 0.0045 & 0.0006 & -0.0051 & -0.0071 & -0.0029 \\
|
|
& & $\Ex{(2)}$ & 0.0851 & 0.0832 & 0.0796 & 0.0725 & 0.0593 & 0.0413 & 0.0315 \\
|
|
& & $\DD{c}{(0)}$ & -0.0155 & -0.0148 & -0.0134 & -0.0110 & -0.0071 & -0.0017 & 0.0040 \\
|
|
& & $\DD{c}{(1)}$ & -0.0020 & -0.0021 & -0.0023 & -0.0027 & -0.0034 & -0.0042 & -0.0044 \\
|
|
& & $\DD{c}{(2)}$ & 0.0175 & 0.0168 & 0.0157 & 0.0137 & 0.0105 & 0.0059 & 0.0004 \\
|
|
\end{tabular}
|
|
\end{ruledtabular}
|
|
\end{table*}
|
|
|
|
%%% TABLE V %%%
|
|
\begin{table*}
|
|
\caption{
|
|
\label{tab:OptGap_5}
|
|
Deviation from the FCI quantities (in hartree) of the individual energies, $\E{(I)}$, and the corresponding excitation energies, $\Ex{(I)}$, for the ground ($I=0$), singly-excited ($I=1$) and doubly-excited ($I=2$) states of 5-boxium (i.e.,~$\Nel = 5$ electrons in a box of length $L$).
|
|
The values of the derivative discontinuity $\DD{c}{(I)}$ are also reported.
|
|
}
|
|
\begin{ruledtabular}
|
|
\begin{tabular}{lclddddddd}
|
|
& & & \mc{7}{c}{5-boxium with a box of length $L$} \\
|
|
\cline{4-10}
|
|
Method & $\bw$ & State & \pi/8 & \pi/4 & \pi/2 & \pi & 2\pi & 4\pi & 8\pi \\
|
|
\hline
|
|
FCI & & $\E{(0)}$ & 1867.6344 & 493.6760 & 136.7020 & 40.7244 & 13.3763 & 4.8811 & 1.9492 \\
|
|
& & $\E{(1)}$ & 2224.11488 & 583.8981 & 159.7957 & 46.7553 & 15.0029 & 5.3399 & 2.0855 \\
|
|
& & $\E{(2)}$ & 3289.2022 & 852.4249 & 228.0415 & 64.3597 & 19.6613 & 6.6206 & 2.4547 \\
|
|
& & $\Ex{(1)}$ & 356.4804 & 90.2221 & 23.0937 & 6.0308 & 1.6266 & 0.4588 & 0.1363 \\
|
|
& & $\Ex{(2)}$ & 1421.56773 & 358.7489 & 91.3395 & 23.6352 & 6.2850 & 1.7395 & 0.5055 \\
|
|
\hline
|
|
& & & \mc{7}{c}{Deviation from FCI} \\
|
|
\hline
|
|
% CIS & & $\Ex{(1)}$ & 0.0230 & 0.0230 & 0.0229 & 0.0229 & 0.0230 & 0.0225 & 0.0182 \\
|
|
% \\
|
|
% TDHF & & $\Ex{(1)}$ & 0.0005 & 0.0005 & 0.0005 & 0.0005 & 0.0005 & 0.0005 & 0.0006 \\
|
|
% \\
|
|
TDA-TDLDA& & $\Ex{(1)}$ & 0.0230 & 0.0230 & 0.0228 & 0.0223 & 0.0192 & -0.0015 & -0.0309 \\
|
|
\\
|
|
TDLDA & & $\Ex{(1)}$ & 0.0005 & 0.0005 & 0.0004 & 0.0000 & -0.0033 & -0.0248 & -0.0650 \\
|
|
\\
|
|
% eHF & $(0,0)$ & $\E{(0)}$ & 0.0769 & 0.0755 & 0.0728 & 0.0676 & 0.0584 & 0.0439 & 0.0257 \\
|
|
% & & $\E{(1)}$ & 0.1642 & 0.1628 & 0.1603 & 0.1556 & 0.1475 & 0.1344 & 0.1120 \\
|
|
% & & $\E{(2)}$ & 0.1800 & 0.1777 & 0.1732 & 0.1645 & 0.1493 & 0.1260 & 0.0988 \\
|
|
% & & $\Ex{(1)}$ & 0.2228 & 0.2215 & 0.2189 & 0.2143 & 0.2066 & 0.1933 & 0.1626 \\
|
|
% & & $\Ex{(2)}$ & 0.2387 & 0.2364 & 0.2318 & 0.2233 & 0.2084 & 0.1850 & 0.1495 \\
|
|
% \\
|
|
% eHF & $(1/3,1/3)$ & $\E{(0)}$ & 0.1426 & 0.1411 & 0.1382 & 0.1326 & 0.1225 & 0.1059 & 0.0823 \\
|
|
% & & $\E{(1)}$ & 0.1347 & 0.1329 & 0.1295 & 0.1233 & 0.1124 & 0.0957 & 0.0737 \\
|
|
% & & $\E{(2)}$ & 0.2019 & 0.1996 & 0.1950 & 0.1860 & 0.1696 & 0.1429 & 0.1070 \\
|
|
% & & $\Ex{(1)}$ & 0.1277 & 0.1260 & 0.1228 & 0.1169 & 0.1071 & 0.0943 & 0.0785 \\
|
|
% & & $\Ex{(2)}$ & 0.1949 & 0.1926 & 0.1882 & 0.1797 & 0.1643 & 0.1414 & 0.1119 \\
|
|
% \\
|
|
KS-eLDA & $(0,0)$ & $\E{(0)}$ & -0.0587 & -0.0586 & -0.0587 & -0.0588 & -0.0591 & -0.0590 & -0.0506 \\
|
|
& & $\E{(1)}$ & 0.0457 & 0.0450 & 0.0435 & 0.0409 & 0.0362 & 0.0241 & 0.0033 \\
|
|
& & $\E{(2)}$ & 0.0861 & 0.0838 & 0.0793 & 0.0712 & 0.0571 & 0.0377 & 0.0196 \\
|
|
& & $\Ex{(1)}$ & 0.1044 & 0.1036 & 0.1022 & 0.0997 & 0.0953 & 0.0830 & 0.0540 \\
|
|
& & $\Ex{(2)}$ & 0.1447 & 0.1424 & 0.1380 & 0.1300 & 0.1162 & 0.0966 & 0.0703 \\
|
|
& & $\DD{c}{(0)}$ & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\
|
|
& & $\DD{c}{(1)}$ & 0.0172 & 0.0163 & 0.0147 & 0.0117 & 0.0067 & -0.0004 & -0.0080 \\
|
|
& & $\DD{c}{(2)}$ & 0.0416 & 0.0402 & 0.0376 & 0.0329 & 0.0253 & 0.0140 & 0.0005 \\
|
|
\\
|
|
KS-eLDA & $(1/3,1/3)$ & $\E{(0)}$ & 0.0070 & 0.0070 & 0.0068 & 0.0063 & 0.0053 & 0.0015 & -0.0049 \\
|
|
& & $\E{(1)}$ & 0.0162 & 0.0151 & 0.0128 & 0.0086 & 0.0018 & -0.0066 & -0.0095 \\
|
|
& & $\E{(2)}$ & 0.1080 & 0.1056 & 0.1011 & 0.0925 & 0.0772 & 0.0538 & 0.0325 \\
|
|
& & $\Ex{(1)}$ & 0.0092 & 0.0081 & 0.0060 & 0.0022 & -0.0035 & -0.0081 & -0.0047 \\
|
|
& & $\Ex{(2)}$ & 0.1010 & 0.0986 & 0.0943 & 0.0862 & 0.0719 & 0.0523 & 0.0373 \\
|
|
& & $\DD{c}{(0)}$ & -0.0196 & -0.0188 & -0.0174 & -0.0148 & -0.0106 & -0.0044 & 0.0029 \\
|
|
& & $\DD{c}{(1)}$ & -0.0024 & -0.0025 & -0.0027 & -0.0032 & -0.0040 & -0.0050 & -0.0056 \\
|
|
& & $\DD{c}{(2)}$ & 0.0220 & 0.0213 & 0.0201 & 0.0180 & 0.0146 & 0.0093 & 0.0027 \\
|
|
\end{tabular}
|
|
\end{ruledtabular}
|
|
\end{table*}
|
|
|
|
%%% TABLE VI %%%
|
|
\begin{table*}
|
|
\caption{
|
|
\label{tab:OptGap_6}
|
|
Deviation from the FCI quantities (in hartree) of the individual energies, $\E{(I)}$, and the corresponding excitation energies, $\Ex{(I)}$, for the ground ($I=0$), singly-excited ($I=1$) and doubly-excited ($I=2$) states of 6-boxium (i.e.,~$\Nel = 6$ electrons in a box of length $L$).
|
|
The values of the derivative discontinuity $\DD{c}{(I)}$ are also reported.
|
|
}
|
|
\begin{ruledtabular}
|
|
\begin{tabular}{lclddddddd}
|
|
& & & \mc{7}{c}{6-boxium with a box of length $L$} \\
|
|
\cline{4-10}
|
|
Method & $\bw$ & State & \pi/8 & \pi/4 & \pi/2 & \pi & 2\pi & 4\pi & 8\pi \\
|
|
\hline
|
|
FCI & & $\E{(0)}$ & 3082.5386 & 813.0910 & 224.3734 & 66.5257 & 21.7454 & 7.9136 & 3.1633 \\
|
|
& & $\E{(1)}$ & 3503.4911 & 919.5487 & 251.5842 & 73.6145 & 23.6504 & 8.4487 & 3.3217 \\
|
|
& & $\E{(2)}$ & 4762.0921 & 1236.8257 & 332.1993 & 94.3988 & 29.1455 & 9.9582 & 3.7572 \\
|
|
& & $\Ex{(1)}$ & 420.9525 & 106.4577 & 27.2108 & 7.0888 & 1.9050 & 0.5351 & 0.1583 \\
|
|
& & $\Ex{(2)}$ & 1679.5536 & 423.7347 & 107.8259 & 27.8731 & 7.4001 & 2.0446 & 0.5938 \\
|
|
\hline
|
|
& & & \mc{7}{c}{Deviation from FCI} \\
|
|
\hline
|
|
% CIS & & $\Ex{(1)}$ & 0.0249 & 0.0248 & 0.0250 & 0.0253 & 0.0261 & 0.0272 & 0.0248 \\
|
|
% \\
|
|
% TDHF & & $\Ex{(1)}$ & 0.0002 & 0.0000 & 0.0003 & 0.0003 & 0.0003 & 0.0003 & 0.0004 \\
|
|
% \\
|
|
TDA-TDLDA& & $\Ex{(1)}$ & 0.0249 & 0.0248 & 0.0250 & 0.0250 & 0.0242 & 0.0114 & -0.0223 \\
|
|
\\
|
|
TDLDA & & $\Ex{(1)}$ & 0.0002 & 0.0000 & 0.0002 & 0.0000 & -0.0016 & -0.0162 & -0.0612 \\
|
|
\\
|
|
% eHF & $(0,0)$ & $\E{(0)}$ & 0.1004 & 0.0989 & 0.0960 & 0.0903 & 0.0802 & 0.0634 & 0.0403 \\
|
|
% & & $\E{(1)}$ & 0.1908 & 0.1893 & 0.1870 & 0.1824 & 0.1747 & 0.1627 & 0.1422 \\
|
|
% & & $\E{(2)}$ & 0.2144 & 0.2120 & 0.2073 & 0.1985 & 0.1830 & 0.1595 & 0.1333 \\
|
|
% & & $\Ex{(1)}$ & 0.2534 & 0.2520 & 0.2498 & 0.2456 & 0.2387 & 0.2281 & 0.2035 \\
|
|
% & & $\Ex{(2)}$ & 0.2770 & 0.2747 & 0.2701 & 0.2617 & 0.2470 & 0.2249 & 0.1946 \\
|
|
% \\
|
|
% eHF & $(1/3,1/3)$ & $\E{(0)}$ & 0.1676 & 0.1660 & 0.1631 & 0.1574 & 0.1471 & 0.1300 & 0.1049 \\
|
|
% & & $\E{(1)}$ & 0.1580 & 0.1561 & 0.1528 & 0.1464 & 0.1352 & 0.1173 & 0.0924 \\
|
|
% & & $\E{(2)}$ & 0.2294 & 0.2270 & 0.2222 & 0.2130 & 0.1966 & 0.1702 & 0.1336 \\
|
|
% & & $\Ex{(1)}$ & 0.1534 & 0.1516 & 0.1485 & 0.1425 & 0.1321 & 0.1167 & 0.0991 \\
|
|
% & & $\Ex{(2)}$ & 0.2249 & 0.2225 & 0.2179 & 0.2091 & 0.1935 & 0.1696 & 0.1404 \\
|
|
% \\
|
|
KS-eLDA & $(0,0)$ & $\E{(0)}$ & -0.0626 & -0.0627 & -0.0628 & -0.0632 & -0.0641 & -0.0654 & -0.0612 \\
|
|
& & $\E{(1)}$ & 0.0486 & 0.0477 & 0.0465 & 0.0440 & 0.0400 & 0.0308 & 0.0078 \\
|
|
& & $\E{(2)}$ & 0.1017 & 0.0992 & 0.0946 & 0.0862 & 0.0718 & 0.0507 & 0.0271 \\
|
|
& & $\Ex{(1)}$ & 0.1112 & 0.1104 & 0.1093 & 0.1072 & 0.1041 & 0.0962 & 0.0690 \\
|
|
& & $\Ex{(2)}$ & 0.1643 & 0.1619 & 0.1575 & 0.1494 & 0.1358 & 0.1162 & 0.0884 \\
|
|
& & $\DD{c}{(0)}$ & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\
|
|
& & $\DD{c}{(1)}$ & 0.0208 & 0.0199 & 0.0182 & 0.0151 & 0.0098 & 0.0018 & -0.0075 \\
|
|
& & $\DD{c}{(2)}$ & 0.0503 & 0.0488 & 0.0460 & 0.0412 & 0.0330 & 0.0205 & 0.0046 \\
|
|
\\
|
|
KS-eLDA & $(1/3,1/3)$ & $\E{(0)}$ & 0.0046 & 0.0045 & 0.0043 & 0.0039 & 0.0031 & 0.0006 & -0.0067 \\
|
|
& & $\E{(1)}$ & 0.0157 & 0.0144 & 0.0123 & 0.0080 & 0.0009 & -0.0091 & -0.0160 \\
|
|
& & $\E{(2)}$ & 0.1167 & 0.1142 & 0.1095 & 0.1007 & 0.0853 & 0.0616 & 0.0355 \\
|
|
& & $\Ex{(1)}$ & 0.0112 & 0.0099 & 0.0080 & 0.0041 & -0.0022 & -0.0097 & -0.0093 \\
|
|
& & $\Ex{(2)}$ & 0.1121 & 0.1097 & 0.1051 & 0.0968 & 0.0822 & 0.0610 & 0.0423 \\
|
|
& & $\DD{c}{(0)}$ & -0.0237 & -0.0229 & -0.0214 & -0.0188 & -0.0142 & -0.0073 & 0.0013 \\
|
|
& & $\DD{c}{(1)}$ & -0.0029 & -0.0030 & -0.0032 & -0.0037 & -0.0045 & -0.0057 & -0.0066 \\
|
|
& & $\DD{c}{(2)}$ & 0.0266 & 0.0259 & 0.0246 & 0.0224 & 0.0187 & 0.0130 & 0.0053 \\
|
|
\end{tabular}
|
|
\end{ruledtabular}
|
|
\end{table*}
|
|
|
|
|
|
%%% TABLE VII %%%
|
|
\begin{table*}
|
|
\caption{
|
|
\label{tab:OptGap_7}
|
|
Deviation from the FCI quantities (in hartree) of the individual energies, $\E{(I)}$, and the corresponding excitation energies, $\Ex{(I)}$, for the ground ($I=0$), singly-excited ($I=1$) and doubly-excited ($I=2$) states of 7-boxium (i.e.,~$\Nel = 7$ electrons in a box of length $L$).
|
|
The values of the derivative discontinuity $\DD{c}{(I)}$ are also reported.
|
|
}
|
|
\begin{ruledtabular}
|
|
\begin{tabular}{lclddddddd}
|
|
& & & \mc{7}{c}{7-boxium with a box of length $L$} \\
|
|
\cline{4-10}
|
|
Method & $\bw$ & State & \pi/8 & \pi/4 & \pi/2 & \pi & 2\pi & 4\pi & 8\pi \\
|
|
\hline
|
|
FCI & & $\E{(0)}$ & 4729.98018 & 1244.7753 & 342.1796 & 100.8943 & 32.7728 & 11.8683 & 4.7359 \\
|
|
& & $\E{(1)}$ & 5215.3307 & 1367.4316 & 373.4897 & 109.0326 & 34.9524 & 12.4779 & 4.9156 \\
|
|
& & $\E{(2)}$ & 6667.18516 & 1733.3319 & 466.4133 & 132.9686 & 41.2715 & 14.2096 & 5.4146 \\
|
|
& & $\Ex{(1)}$ & 485.3505 & 122.6563 & 31.3101 & 8.1382 & 2.1796 & 0.6096 & 0.1797 \\
|
|
& & $\Ex{(2)}$ & 1937.2050 & 488.5566 & 124.2336 & 32.0743 & 8.4987 & 2.3413 & 0.6787 \\
|
|
\hline
|
|
& & & \mc{7}{c}{Deviation from FCI} \\
|
|
\hline
|
|
% CIS & & $\Ex{(1)}$ & 0.0262 & 0.0264 & 0.0265 & 0.0270 & 0.0283 & 0.0308 & 0.0309 \\
|
|
% \\
|
|
% TDHF & & $\Ex{(1)}$ & 0.0000 & 0.0001 & 0.0000 & 0.0000 & 0.0000 & 0.0001 & 0.0003 \\
|
|
% \\
|
|
TDA-TDLDA& & $\Ex{(1)}$ & 0.0262 & 0.0264 & 0.0264 & 0.0269 & 0.0273 & 0.0206 & -0.0116 \\
|
|
\\
|
|
TDLDA & & $\Ex{(1)}$ & 0.0000 & 0.0001 & 0.0000 & -0.0001 & -0.0009 & -0.0107 & -0.0539 \\
|
|
\\
|
|
% eHF & $(0,0)$ & $\E{(0)}$ & 0.1240 & 0.1224 & 0.1194 & 0.1135 & 0.1027 & 0.0842 & 0.0570 \\
|
|
% & & $\E{(1)}$ & 0.2162 & 0.2150 & 0.2125 & 0.2081 & 0.2006 & 0.1894 & 0.1711 \\
|
|
% & & $\E{(2)}$ & 0.2437 & 0.2419 & 0.2376 & 0.2291 & 0.2137 & 0.1916 & 0.1650 \\
|
|
% & & $\Ex{(1)}$ & 0.2473 & 0.2458 & 0.2425 & 0.2366 & 0.2263 & 0.2099 & 0.1825 \\
|
|
% & & $\Ex{(2)}$ & 0.1796 & 0.1780 & 0.1746 & 0.1685 & 0.1576 & 0.1406 & 0.1202 \\
|
|
% \\
|
|
% eHF & $(1/3,1/3)$ & $\E{(0)}$ & 0.1918 & 0.1902 & 0.1873 & 0.1817 & 0.1713 & 0.1538 & 0.1273 \\
|
|
% & & $\E{(1)}$ & 0.1810 & 0.1792 & 0.1758 & 0.1694 & 0.1579 & 0.1392 & 0.1122 \\
|
|
% & & $\E{(2)}$ & 0.2532 & 0.2514 & 0.2470 & 0.2381 & 0.2220 & 0.1969 & 0.1584 \\
|
|
% & & $\Ex{(1)}$ & 0.1796 & 0.1780 & 0.1746 & 0.1685 & 0.1576 & 0.1406 & 0.1202 \\
|
|
% & & $\Ex{(2)}$ & 0.2518 & 0.2501 & 0.2458 & 0.2372 & 0.2217 & 0.1982 & 0.1663 \\
|
|
% \\
|
|
KS-eLDA & $(0,0)$ & $\E{(0)}$ & -0.0664 & -0.0666 & -0.0667 & -0.0672 & -0.0684 & -0.0707 & -0.0702 \\
|
|
& & $\E{(1)}$ & 0.0502 & 0.0495 & 0.0482 & 0.0459 & 0.0423 & 0.0355 & 0.0131 \\
|
|
& & $\E{(2)}$ & 0.1122 & 0.1104 & 0.1061 & 0.0979 & 0.0836 & 0.0635 & 0.0360 \\
|
|
& & $\Ex{(1)}$ & 0.1165 & 0.1161 & 0.1149 & 0.1131 & 0.1108 & 0.1062 & 0.0834 \\
|
|
& & $\Ex{(2)}$ & 0.1785 & 0.1769 & 0.1728 & 0.1652 & 0.1520 & 0.1342 & 0.1063 \\
|
|
& & $\DD{c}{(0)}$ & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\
|
|
& & $\DD{c}{(1)}$ & 0.0244 & 0.0235 & 0.0218 & 0.0186 & 0.0130 & 0.0043 & -0.0065 \\
|
|
& & $\DD{c}{(2)}$ & 0.0589 & 0.0574 & 0.0546 & 0.0496 & 0.0410 & 0.0275 & 0.0095 \\
|
|
\\
|
|
KS-eLDA & $(1/3,1/3)$ & $\E{(0)}$ & 0.0014 & 0.0013 & 0.0012 & 0.0009 & 0.0003 & -0.0013 & -0.0079 \\
|
|
& & $\E{(1)}$ & 0.0149 & 0.0138 & 0.0115 & 0.0072 & -0.0001 & -0.0110 & -0.0209 \\
|
|
& & $\E{(2)}$ & 0.1217 & 0.1198 & 0.1154 & 0.1069 & 0.0917 & 0.0691 & 0.0389 \\
|
|
& & $\Ex{(1)}$ & 0.0135 & 0.0125 & 0.0103 & 0.0063 & -0.0005 & -0.0096 & -0.0130 \\
|
|
& & $\Ex{(2)}$ & 0.1203 & 0.1185 & 0.1142 & 0.1060 & 0.0914 & 0.0705 & 0.0469 \\
|
|
& & $\DD{c}{(0)}$ & -0.0278 & -0.0270 & -0.0255 & -0.0227 & -0.0180 & -0.0105 & -0.0007 \\
|
|
& & $\DD{c}{(1)}$ & -0.0034 & -0.0034 & -0.0037 & -0.0041 & -0.0050 & -0.0063 & -0.0076 \\
|
|
& & $\DD{c}{(2)}$ & 0.0311 & 0.0304 & 0.0291 & 0.0268 & 0.0230 & 0.0168 & 0.0083 \\
|
|
\end{tabular}
|
|
\end{ruledtabular}
|
|
\end{table*}
|
|
|
|
|
|
\bibliography{../eDFT}
|
|
|
|
\end{document}
|