From ba6fc13e4112fcbb12cba94b9475f95d555081aa Mon Sep 17 00:00:00 2001 From: Pierre-Francois Loos Date: Wed, 13 Nov 2019 23:48:28 +0100 Subject: [PATCH] xc func --- FarDFT/FarDFT.nb | 140 +++++++++++++-------- FarDFT/Manuscript/FarDFT.bib | 2 +- FarDFT/Manuscript/FarDFT.tex | 237 ++++++++++++++++++++--------------- 3 files changed, 226 insertions(+), 153 deletions(-) diff --git a/FarDFT/FarDFT.nb b/FarDFT/FarDFT.nb index 6f1df33..86dc93e 100644 --- a/FarDFT/FarDFT.nb +++ b/FarDFT/FarDFT.nb @@ -10,10 +10,10 @@ NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 158, 7] -NotebookDataLength[ 61820, 1710] -NotebookOptionsPosition[ 52802, 1565] -NotebookOutlinePosition[ 53155, 1581] -CellTagsIndexPosition[ 53112, 1578] +NotebookDataLength[ 63150, 1752] +NotebookOptionsPosition[ 53923, 1604] +NotebookOutlinePosition[ 54259, 1619] +CellTagsIndexPosition[ 54216, 1616] WindowFrame->Normal*) (* Beginning of Notebook Content *) @@ -37,7 +37,7 @@ mathpazo,xcolor,bm,mhchem}\>\"", "}"}]}]}], "]"}], ";"}]}], "Input", InitializationCell->True, CellChangeTimes->{{3.7288240181604652`*^9, 3.728824027007351*^9}, { 3.733131339213026*^9, 3.733131352923026*^9}}, - CellLabel->"In[4]:=",ExpressionUUID->"8b7c0bbd-9f99-42dc-b8a5-a4dc3cb8cc89"], + CellLabel->"In[1]:=",ExpressionUUID->"8b7c0bbd-9f99-42dc-b8a5-a4dc3cb8cc89"], Cell[BoxData[ RowBox[{ @@ -45,7 +45,7 @@ Cell[BoxData[ InitializationCell->True, CellChangeTimes->{{3.7208031947801647`*^9, 3.7208032000677156`*^9}, { 3.7208034541742477`*^9, 3.720803455246439*^9}}, - CellLabel->"In[9]:=",ExpressionUUID->"b24e8243-3c38-4b03-af3a-708cb1afa004"] + CellLabel->"In[3]:=",ExpressionUUID->"b24e8243-3c38-4b03-af3a-708cb1afa004"] }, Closed]], Cell[CellGroupData[{ @@ -812,6 +812,46 @@ Cell[BoxData[ 3.782626306892585*^9}},ExpressionUUID->"e2512605-e53f-4895-bc29-\ b005c70f53b5"], +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{ + RowBox[{ + FractionBox["3", + RowBox[{"2", + SuperscriptBox["R", "2"]}]], "+", + FractionBox["176", + RowBox[{"105", " ", "\[Pi]", " ", "R"}]]}], "/.", + RowBox[{"R", "\[Rule]", + FractionBox["1", + RowBox[{ + SuperscriptBox["\[Pi]", + RowBox[{"2", "/", "3"}]], " ", + SuperscriptBox["\[Rho]", + RowBox[{"1", "/", "3"}]]}]]}]}]], "Input", + CellChangeTimes->{{3.782672032733528*^9, 3.782672036312338*^9}, + 3.782672105472003*^9}, + CellLabel->"In[6]:=",ExpressionUUID->"2caf6c5b-6ae4-41ea-a5bc-6e2553b311b0"], + +Cell[BoxData[ + RowBox[{ + FractionBox[ + RowBox[{"176", " ", + SuperscriptBox["\[Rho]", + RowBox[{"1", "/", "3"}]]}], + RowBox[{"105", " ", + SuperscriptBox["\[Pi]", + RowBox[{"1", "/", "3"}]]}]], "+", + RowBox[{ + FractionBox["3", "2"], " ", + SuperscriptBox["\[Pi]", + RowBox[{"4", "/", "3"}]], " ", + SuperscriptBox["\[Rho]", + RowBox[{"2", "/", "3"}]]}]}]], "Output", + CellChangeTimes->{3.782672036741086*^9, 3.7826721062008343`*^9}, + CellLabel->"Out[6]=",ExpressionUUID->"71a67885-8142-46b8-9292-410654424b59"] +}, Open ]], + Cell["This energy can be decomposed as", "Text", CellChangeTimes->{{3.782625326608221*^9, 3.7826253316050653`*^9}},ExpressionUUID->"32015b28-f5cc-4eff-8d8d-\ @@ -1066,7 +1106,6 @@ transfered to the ", Cell[BoxData[ FormBox[ SubscriptBox["C", "x"], TraditionalForm]], - FormatType->"TraditionalForm", FontColor->RGBColor[1, 0, 0],ExpressionUUID-> "fa14d797-4f81-4166-97c2-b14e12078e3f"], StyleBox[" coefficient.\nThis is obvious but kind of nice.", @@ -1565,8 +1604,7 @@ c0652866375e"] }, WindowSize->{1280, 755}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, -FrontEndVersion->"11.3 for Mac OS X x86 (32-bit, 64-bit Kernel) (March 5, \ -2018)", +FrontEndVersion->"12.0 for Mac OS X x86 (64-bit) (April 8, 2019)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) @@ -1655,60 +1693,64 @@ Cell[25417, 763, 310, 7, 35, "Text",ExpressionUUID->"137a2815-3e41-4c54-bb77-2ba Cell[25730, 772, 570, 16, 49, "Input",ExpressionUUID->"ee1a49a3-d78d-48fd-b18b-b45c556268d6"], Cell[26303, 790, 318, 6, 35, "Text",ExpressionUUID->"7f6139ea-87a7-485f-96a7-072785ed5a50"], Cell[26624, 798, 571, 14, 48, "Input",ExpressionUUID->"e2512605-e53f-4895-bc29-b005c70f53b5"], -Cell[27198, 814, 176, 3, 35, "Text",ExpressionUUID->"32015b28-f5cc-4eff-8d8d-b8691b888e8f"], -Cell[27377, 819, 847, 23, 48, "Input",ExpressionUUID->"1fcbd290-f3ea-43bd-b2ef-f31601137679"], -Cell[28227, 844, 201, 3, 35, "Text",ExpressionUUID->"9b1e04fb-55cb-4e81-8daf-ba707c78dfce"], -Cell[28431, 849, 827, 22, 45, "Input",ExpressionUUID->"c51ccbbc-5ee1-4cca-ade7-af230c7f6ba6"], -Cell[29261, 873, 153, 3, 35, "Text",ExpressionUUID->"8bfb35af-a4e4-4813-ae15-3703f4aeab96"], -Cell[29417, 878, 1240, 35, 50, "Input",ExpressionUUID->"60e340d1-6bd1-4928-9814-51c59353edd1"], -Cell[30660, 915, 250, 4, 35, "Text",ExpressionUUID->"f5f45693-a80c-4a46-ad73-07811620c0d9"], -Cell[30913, 921, 787, 23, 58, "Input",ExpressionUUID->"f99857fd-0430-402d-a733-929b6b359e23"] +Cell[CellGroupData[{ +Cell[27220, 816, 570, 17, 50, "Input",ExpressionUUID->"2caf6c5b-6ae4-41ea-a5bc-6e2553b311b0"], +Cell[27793, 835, 544, 16, 55, "Output",ExpressionUUID->"71a67885-8142-46b8-9292-410654424b59"] +}, Open ]], +Cell[28352, 854, 176, 3, 35, "Text",ExpressionUUID->"32015b28-f5cc-4eff-8d8d-b8691b888e8f"], +Cell[28531, 859, 847, 23, 48, "Input",ExpressionUUID->"1fcbd290-f3ea-43bd-b2ef-f31601137679"], +Cell[29381, 884, 201, 3, 35, "Text",ExpressionUUID->"9b1e04fb-55cb-4e81-8daf-ba707c78dfce"], +Cell[29585, 889, 827, 22, 45, "Input",ExpressionUUID->"c51ccbbc-5ee1-4cca-ade7-af230c7f6ba6"], +Cell[30415, 913, 153, 3, 35, "Text",ExpressionUUID->"8bfb35af-a4e4-4813-ae15-3703f4aeab96"], +Cell[30571, 918, 1240, 35, 50, "Input",ExpressionUUID->"60e340d1-6bd1-4928-9814-51c59353edd1"], +Cell[31814, 955, 250, 4, 35, "Text",ExpressionUUID->"f5f45693-a80c-4a46-ad73-07811620c0d9"], +Cell[32067, 961, 787, 23, 58, "Input",ExpressionUUID->"f99857fd-0430-402d-a733-929b6b359e23"] }, Open ]], Cell[CellGroupData[{ -Cell[31737, 949, 348, 5, 54, "Subsection",ExpressionUUID->"ef62b877-dbc0-41f4-bb50-f7bab5069a50"], -Cell[32088, 956, 252, 6, 35, "Text",ExpressionUUID->"5b10be0f-a732-422f-9446-681c78cbd103"], -Cell[32343, 964, 561, 17, 32, "Input",ExpressionUUID->"009be4dc-79a6-499d-b9af-ddbb9078f85a"], -Cell[32907, 983, 918, 28, 53, "Input",ExpressionUUID->"afd8bd15-999f-411c-a673-53d5b75e5201"], -Cell[33828, 1013, 297, 4, 35, "Text",ExpressionUUID->"1fe4835f-3b82-458a-9e2a-9272922b081a"], -Cell[34128, 1019, 1329, 40, 58, "Input",ExpressionUUID->"b5e6e435-424a-4f3e-8416-709ae4356a73"], -Cell[35460, 1061, 579, 15, 58, "Text",ExpressionUUID->"d05917ae-a767-4a78-834a-c4b18c1c0b13"] +Cell[32891, 989, 348, 5, 54, "Subsection",ExpressionUUID->"ef62b877-dbc0-41f4-bb50-f7bab5069a50"], +Cell[33242, 996, 252, 6, 35, "Text",ExpressionUUID->"5b10be0f-a732-422f-9446-681c78cbd103"], +Cell[33497, 1004, 561, 17, 32, "Input",ExpressionUUID->"009be4dc-79a6-499d-b9af-ddbb9078f85a"], +Cell[34061, 1023, 918, 28, 53, "Input",ExpressionUUID->"afd8bd15-999f-411c-a673-53d5b75e5201"], +Cell[34982, 1053, 297, 4, 35, "Text",ExpressionUUID->"1fe4835f-3b82-458a-9e2a-9272922b081a"], +Cell[35282, 1059, 1329, 40, 58, "Input",ExpressionUUID->"b5e6e435-424a-4f3e-8416-709ae4356a73"], +Cell[36614, 1101, 546, 14, 58, "Text",ExpressionUUID->"d05917ae-a767-4a78-834a-c4b18c1c0b13"] }, Open ]], Cell[CellGroupData[{ -Cell[36076, 1081, 419, 8, 54, "Subsection",ExpressionUUID->"bcd158fd-9aa7-4c7b-837e-afd37c54f55a"], -Cell[36498, 1091, 571, 11, 81, "Text",ExpressionUUID->"a0112f10-d580-4370-a6ed-bff98b24a4a6"], -Cell[37072, 1104, 2757, 86, 101, "Input",ExpressionUUID->"aabf1284-1a62-408f-8d1c-50c99b5acc35"], -Cell[39832, 1192, 278, 6, 35, "Text",ExpressionUUID->"88a4bd50-0ba7-4f7c-9fc0-e0feb01831fc"], -Cell[40113, 1200, 606, 19, 48, "Input",ExpressionUUID->"a889f096-9a50-426c-9cf3-81c363ec2eaf"], -Cell[40722, 1221, 355, 7, 35, "Text",ExpressionUUID->"d3c69495-05ef-49fe-8484-fa7fd62cc408"], -Cell[41080, 1230, 1887, 55, 58, "Input",ExpressionUUID->"0138b572-760b-4f82-a618-73197e776fa2"], -Cell[42970, 1287, 452, 7, 35, "Text",ExpressionUUID->"b83db319-04e5-42b8-ae4c-57f93abe09fd"] +Cell[37197, 1120, 419, 8, 54, "Subsection",ExpressionUUID->"bcd158fd-9aa7-4c7b-837e-afd37c54f55a"], +Cell[37619, 1130, 571, 11, 81, "Text",ExpressionUUID->"a0112f10-d580-4370-a6ed-bff98b24a4a6"], +Cell[38193, 1143, 2757, 86, 101, "Input",ExpressionUUID->"aabf1284-1a62-408f-8d1c-50c99b5acc35"], +Cell[40953, 1231, 278, 6, 35, "Text",ExpressionUUID->"88a4bd50-0ba7-4f7c-9fc0-e0feb01831fc"], +Cell[41234, 1239, 606, 19, 48, "Input",ExpressionUUID->"a889f096-9a50-426c-9cf3-81c363ec2eaf"], +Cell[41843, 1260, 355, 7, 35, "Text",ExpressionUUID->"d3c69495-05ef-49fe-8484-fa7fd62cc408"], +Cell[42201, 1269, 1887, 55, 58, "Input",ExpressionUUID->"0138b572-760b-4f82-a618-73197e776fa2"], +Cell[44091, 1326, 452, 7, 35, "Text",ExpressionUUID->"b83db319-04e5-42b8-ae4c-57f93abe09fd"] }, Open ]] }, Open ]], Cell[CellGroupData[{ -Cell[43471, 1300, 240, 4, 67, "Section",ExpressionUUID->"258ca12d-c8c5-4bc9-84ec-1711f7860786"], +Cell[44592, 1339, 240, 4, 67, "Section",ExpressionUUID->"258ca12d-c8c5-4bc9-84ec-1711f7860786"], Cell[CellGroupData[{ -Cell[43736, 1308, 347, 5, 54, "Subsection",ExpressionUUID->"ee1991f0-6627-4afa-871c-73710071f6fd"], -Cell[44086, 1315, 324, 7, 58, "Text",ExpressionUUID->"a6c10556-4f53-4da2-8cf0-79bbf3ec1242"], -Cell[44413, 1324, 1209, 36, 64, "Input",ExpressionUUID->"1885d9ca-6bb6-490d-a0bb-a3b96576a614"], -Cell[45625, 1362, 300, 6, 35, "Text",ExpressionUUID->"ead26072-79d4-4549-aab6-46fe8cf9b6db"] +Cell[44857, 1347, 347, 5, 54, "Subsection",ExpressionUUID->"ee1991f0-6627-4afa-871c-73710071f6fd"], +Cell[45207, 1354, 324, 7, 58, "Text",ExpressionUUID->"a6c10556-4f53-4da2-8cf0-79bbf3ec1242"], +Cell[45534, 1363, 1209, 36, 64, "Input",ExpressionUUID->"1885d9ca-6bb6-490d-a0bb-a3b96576a614"], +Cell[46746, 1401, 300, 6, 35, "Text",ExpressionUUID->"ead26072-79d4-4549-aab6-46fe8cf9b6db"] }, Open ]], Cell[CellGroupData[{ -Cell[45962, 1373, 353, 5, 54, "Subsection",ExpressionUUID->"ac951d60-fbce-44ba-8f12-9e5545df5beb"], -Cell[46318, 1380, 228, 5, 35, "Text",ExpressionUUID->"dca5790f-2358-43e0-b359-677842457a7d"], -Cell[46549, 1387, 1279, 38, 64, "Input",ExpressionUUID->"cce86f90-52f7-41f4-a63d-73a6c518b4c3"] +Cell[47083, 1412, 353, 5, 54, "Subsection",ExpressionUUID->"ac951d60-fbce-44ba-8f12-9e5545df5beb"], +Cell[47439, 1419, 228, 5, 35, "Text",ExpressionUUID->"dca5790f-2358-43e0-b359-677842457a7d"], +Cell[47670, 1426, 1279, 38, 64, "Input",ExpressionUUID->"cce86f90-52f7-41f4-a63d-73a6c518b4c3"] }, Open ]], Cell[CellGroupData[{ -Cell[47865, 1430, 348, 5, 54, "Subsection",ExpressionUUID->"9f2ce783-2740-403c-974c-46a7ef9321fe"], -Cell[48216, 1437, 310, 7, 35, "Text",ExpressionUUID->"6c7b6dc7-06c2-4f95-a398-fe8f7adea7db"], -Cell[48529, 1446, 610, 18, 32, "Input",ExpressionUUID->"c5451bea-41ca-4be3-9429-a6d8ddd06c4a"] +Cell[48986, 1469, 348, 5, 54, "Subsection",ExpressionUUID->"9f2ce783-2740-403c-974c-46a7ef9321fe"], +Cell[49337, 1476, 310, 7, 35, "Text",ExpressionUUID->"6c7b6dc7-06c2-4f95-a398-fe8f7adea7db"], +Cell[49650, 1485, 610, 18, 32, "Input",ExpressionUUID->"c5451bea-41ca-4be3-9429-a6d8ddd06c4a"] }, Open ]], Cell[CellGroupData[{ -Cell[49176, 1469, 468, 8, 54, "Subsection",ExpressionUUID->"ab8ea4d9-6951-44dc-96fc-10c74ccd1a4f"], -Cell[49647, 1479, 279, 6, 35, "Text",ExpressionUUID->"f0dac186-7114-4a5e-b53c-f938d8c4fd37"], -Cell[49929, 1487, 1342, 40, 33, "Input",ExpressionUUID->"2ca8802e-361f-4487-a3f0-67d18fc8fe5f"], -Cell[51274, 1529, 386, 8, 35, "Text",ExpressionUUID->"00fbe382-6c25-4893-a74d-546919a831b8"], -Cell[51663, 1539, 368, 8, 32, "Input",ExpressionUUID->"b7e98868-257f-43e4-b057-ff1cbb083046"], -Cell[52034, 1549, 728, 11, 104, "Text",ExpressionUUID->"336e8167-7e94-4334-9956-c0652866375e"] +Cell[50297, 1508, 468, 8, 54, "Subsection",ExpressionUUID->"ab8ea4d9-6951-44dc-96fc-10c74ccd1a4f"], +Cell[50768, 1518, 279, 6, 35, "Text",ExpressionUUID->"f0dac186-7114-4a5e-b53c-f938d8c4fd37"], +Cell[51050, 1526, 1342, 40, 33, "Input",ExpressionUUID->"2ca8802e-361f-4487-a3f0-67d18fc8fe5f"], +Cell[52395, 1568, 386, 8, 35, "Text",ExpressionUUID->"00fbe382-6c25-4893-a74d-546919a831b8"], +Cell[52784, 1578, 368, 8, 32, "Input",ExpressionUUID->"b7e98868-257f-43e4-b057-ff1cbb083046"], +Cell[53155, 1588, 728, 11, 104, "Text",ExpressionUUID->"336e8167-7e94-4334-9956-c0652866375e"] }, Open ]] }, Open ]] }, Open ]] diff --git a/FarDFT/Manuscript/FarDFT.bib b/FarDFT/Manuscript/FarDFT.bib index 758621a..012e98b 100644 --- a/FarDFT/Manuscript/FarDFT.bib +++ b/FarDFT/Manuscript/FarDFT.bib @@ -1,7 +1,7 @@ %% This BibTeX bibliography file was created using BibDesk. %% http://bibdesk.sourceforge.net/ -%% Created for Pierre-Francois Loos at 2019-11-12 22:50:55 +0100 +%% Created for Pierre-Francois Loos at 2019-11-13 20:40:48 +0100 %% Saved with string encoding Unicode (UTF-8) diff --git a/FarDFT/Manuscript/FarDFT.tex b/FarDFT/Manuscript/FarDFT.tex index 0e2a49c..be6aecb 100644 --- a/FarDFT/Manuscript/FarDFT.tex +++ b/FarDFT/Manuscript/FarDFT.tex @@ -1,15 +1,11 @@ \documentclass[aip,jcp,reprint,noshowkeys,superscriptaddress]{revtex4-1} -\usepackage{graphicx,dcolumn,bm,xcolor,microtype,multirow,amscd,amsmath,amssymb,amsfonts,physics,mhchem,longtable,pifont,wrapfig,txfonts,multirow} +\usepackage{graphicx,dcolumn,bm,xcolor,microtype,multirow,amsmath,amssymb,amsfonts,physics,mhchem} - -\usepackage{hyperref} -\hypersetup{ - colorlinks=true, - linkcolor=blue, - filecolor=blue, - urlcolor=blue, - citecolor=blue -} +\usepackage[ + colorlinks=true, + citecolor=blue, + breaklinks=true + ]{hyperref} \urlstyle{same} \newcommand{\alert}[1]{\textcolor{red}{#1}} @@ -22,15 +18,6 @@ \newcommand{\cmark}{\color{green}{\text{\ding{51}}}} \newcommand{\xmark}{\color{red}{\text{\ding{55}}}} -\usepackage{hyperref} -\hypersetup{ - colorlinks=true, - linkcolor=blue, - filecolor=blue, - urlcolor=blue, - citecolor=blue -} - %useful stuff \newcommand{\ie}{\textit{i.e.}} \newcommand{\eg}{\textit{e.g.}} @@ -55,6 +42,7 @@ % functionals, potentials, densities, etc \newcommand{\eps}{\epsilon} \newcommand{\e}[2]{\eps_\text{#1}^{#2}} +\newcommand{\kin}[2]{t_\text{#1}^{#2}} \newcommand{\E}[2]{E_\text{#1}^{#2}} \newcommand{\bE}[2]{\overline{E}_\text{#1}^{#2}} \newcommand{\be}[2]{\overline{\eps}_\text{#1}^{#2}} @@ -62,6 +50,7 @@ \newcommand{\n}[2]{n_{#1}^{#2}} \newcommand{\DD}[2]{\Delta_\text{#1}^{#2}} \newcommand{\LZ}[2]{\Xi_\text{#1}^{#2}} +\newcommand{\Cx}[1]{C_\text{x}^{#1}} % energies \newcommand{\EHF}{E_\text{HF}} @@ -71,6 +60,7 @@ \newcommand{\Eani}{E_\text{ani}} \newcommand{\EPT}{E_\text{PT2}} \newcommand{\EFCI}{E_\text{FCI}} +\newcommand{\LDA}{\text{LDA}} % matrices \newcommand{\br}{\bm{r}} @@ -132,29 +122,36 @@ Their accuracy is illustrated by computing on the prototypical H$_2$ molecule. %%% INTRODUCTION %%% %%%%%%%%%%%%%%%%%%%% \section{Introduction} -Here is a nice introduction. - -%Most of the standard local and semi-local DFAs rely on the infinite uniform electron gas (UEG) model (also known as jellium). \cite{ParrBook, Loos_2016} -%One major drawback of the jellium paradigm, when it comes to develop eDFAs, is that the ground and excited states cannot be easily identified like in a molecule. \cite{Gill_2012, Loos_2012, Loos_2014a, Loos_2014b, Agboola_2015, Loos_2017a} -%Moreover, because the infinite UEG model is a metal, it is gapless, which means that both the fundamental and optical gaps are zero. -%From this point of view, using finite UEGs \cite{Loos_2011b, Gill_2012} (which have, like an atom, discrete energy levels) to construct eDFAs can be seen as more relevant. \cite{Loos_2014a, Loos_2014b, Loos_2017a} -% -%As a finite uniform electron gas, we consider the glomium model in which electrons move on the surface of a four-dimensional sphere (also known as a glome). \cite{Loos_2011b, Agboola_2015} -%The most appealing feature of glomium (regarding the development of functionals in the context of eDFT) is the fact that both the ground- and excited-state densities have uniform densities. -%As a result, the ensemble density will remain constant (and uniform) as the ensemble weights vary. -%This is a necessary condition for being able to model derivative discontinuities. +Time-dependent density-functional theory (TD-DFT) has been the dominant force in the calculation of excitation energies of molecular systems in the last two decades. +At a relatively low computational cost (at least compared to the other excited-state methods), TD-DFT can provide accurate transition energies for low-lying excited states in organic molecules. +Importantly, setting up a TD-DFT calculation for a given system is an almost pain-free process from the user perspective as the only (yet essential) input variable is the choice of the so-called exchange-correlation (xc) functional. +Indeed, TD-DFT is a in-principle exact theory which recast the many-body problem by transferring its complexity to the xc functional. +However, TD-DFT is far from being perfect, and, in practice, approximations must be made for the xc functional. +One of its issues actually originates directly from the choice of the xc functional, and more specifically, the possible substantial variations in the quality of the excitation energy for two different choices of xc functionals. +Moreover, because it was so popular, it has been studied in excruciated details, and researchers have quickly unveiled various theoretical and practical deficiencies of approximate TD-DFT. +Practically, the TD-DFT xc kernel is usually considered as static instead of being frequency dependent. +One key consequence of this so-called adiabatic approximation is that double excitations are completely absent from the TD-DFT spectra. +Moreover, TD-DFT has problems with charge-transfer and Rydberg excited states (the excitation energies are usually drastically underestimated) due to the wrong asymptotic behaviour of the xc functional. +The paper is organised as follows. +In Sec.~\ref{sec:theo}, ... +Section \ref{sec:func} provides details about the construction of the weight-dependent exchange-correlation functional. +The results of our calculations for the prototypical \ce{H2} molecule are reported and discussed in Sec.~\ref{sec:res}. +Finally, we draw our conclusions in Sec.~\ref{sec:ccl}. +Unless otherwise stated, atomic units are used throughout. %%%%%%%%%%%%%%%%%%%% %%% THEORY %%% %%%%%%%%%%%%%%%%%%%% \section{Theory} +\label{sec:theo} Here is the theory. %%%%%%%%%%%%%%%%%% %%% FUNCTIONAL %%% %%%%%%%%%%%%%%%%%% \section{Functional} +\label{sec:func} We adopt the usual decomposition, and write down the weight-dependent exchange-correlation functional as \begin{equation} @@ -162,30 +159,85 @@ We adopt the usual decomposition, and write down the weight-dependent exchange-c \end{equation} where $\e{x}{\ew{}}(\n{}{})$ and $\e{c}{\ew{}}(\n{}{})$ are the weight-dependent exchange and correlation functionals, respectively. The construction of these two functionals is described below. -Here, we restrict our study to spin-unpolarized systems, i.e., closed-shell systems. +Here, we restrict our study to spin-unpolarized systems, \ie, $\n{\uparrow}{} = \n{\downarrow}{} = \n{}{}/2$. The present weight-dependent eDFA is specifically designed for the calculation of double excitations within eDFT. As mentioned previously, we consider a two-state ensemble including the ground state ($I=0$), and the first doubly-excited state ($I=1$) of the two-electron glomium system. -All these states have the same (uniform) density $\n{}{} = 2/(4\pi/3 R^3)$ where $R$ is the radius of the glome where the electrons are confined. +All these states have the same (uniform) density $\n{}{} = 2/(2\pi^2 R^3)$ where $R$ is the radius of the glome where the electrons are confined. We refer the interested reader to Refs.~\onlinecite{Loos_2011b} for more details about this paradigm. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \subsection{Weight-dependent exchange functional} \label{sec:Ex} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +The reduced (\ie, per electron) HF energy for these two states is +\begin{subequations} +\begin{align} + \e{HF}{(0)}(\n{}{}) & = \frac{4}{3\pi^{1/3}} \n{}{1/3}, + \\ + \e{HF}{(1)}(\n{}{}) & = \frac{3\pi^{4/3}}{2} \n{}{2/3} + \frac{176}{105\pi^{1/3}} \n{}{1/3}. +\end{align} +\end{subequations} +These two energies can be conveniently decomposed as +\begin{equation} + \e{HF}{(I)}(\n{}{}) = \kin{s}{(I)}(\n{}{}) + \e{H}{(0)}(\n{}{}) + \e{x}{(I)}(\n{}{}), +\end{equation} +with +\begin{subequations} +\begin{align} + \kin{s}{(0)}(\n{}{}) & = 0, + & + \kin{s}{(1)}(\n{}{}) & = \frac{3\pi^{4/3}}{2} \n{}{2/3}. + \\ + \e{H}{(0)}(\n{}{}) & = \frac{8}{3\pi^{1/3}} \n{}{1/3}, + & + \e{H}{(1)}(\n{}{}) & = \frac{352}{105\pi^{1/3}} \n{}{1/3}. + \\ + \e{x}{(0)}(\n{}{}) & = - \frac{4}{3\pi^{1/3}} \n{}{1/3}, + & + \e{x}{(1)}(\n{}{}) & = - \frac{176}{105\pi^{1/3}} \n{}{1/3}. +\end{align} +\end{subequations} + +Knowing that the exchange functional has the following form +\begin{equation} + \e{x}{(I)}(\n{}{}) = \Cx{(I)} \n{}{1/3} +\end{equation} +we obtain +\begin{align} + \Cx{(0)} & = - \frac{4}{3} \qty( \frac{2}{\pi} )^{1/3}, + & + \Cx{(1)} & = - \frac{176}{105} \qty( \frac{2}{\pi} )^{1/3} +\end{align} +We can now combine these two exchange functionals to create a weight-dependent exchange functional +\begin{equation} +\begin{split} + \e{x}{\ew{}}(\n{}{}) + & = (1-\ew{}) \e{x}{(0)}(\n{}{}) + \ew{} \e{x}{(1)}(\n{}{}) + \\ + & = \Cx{\ew{}} \n{}{1/3} +\end{split} +\end{equation} +with +\begin{equation} + \Cx{\ew{}} = (1-\ew{}) \Cx{(0)} + \ew{} \Cx{(1)} +\end{equation} +Amazingly, the weight dependence of the exchange functional can be transfered to the Subscript[C, x] coefficient. +This is obvious but kind of nice. + %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \subsection{Weight-dependent correlation functional} \label{sec:Ec} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -Based on highly-accurate calculations (see below), one can write down, for each state, an accurate analytical expression of the reduced (i.e., per electron) correlation energy \cite{Loos_2013a, Loos_2014a} via the following Pad\'e approximant +Based on 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{xc}{(I)}(\n{}{}) = \frac{c_1^{(I)}\,\n{}{}}{\n{}{} + c_2^{(I)} \sqrt{\n{}{}} + c_3^{(I)}}, + \e{c}{(I)}(\n{}{}) = \frac{a_1^{(I)}}{1 + a_2^{(I)} \n{}{-1/6} + a_3^{(I)} \n{}{-1/3}}, \end{equation} -where the $c_k^{(I)}$'s are state-specific fitting parameters, which are provided in Table \ref{tab:OG_func}. -The value of $c_1^{(I)}$ is obtained via the exact high-density expansion of the correlation energy. \cite{Loos_2013a, Loos_2014a} +where the $a_k^{(I)}$'s are state-specific fitting parameters, which are provided in Table \ref{tab:OG_func}. +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} provides two state-specific correlation DFAs based on a two-electron system. Combining these, one can build a two-state weight-dependent correlation eDFA: \begin{equation} @@ -193,46 +245,11 @@ Combining these, one can build a two-state weight-dependent correlation eDFA: \e{c}{\ew{}}(\n{}{}) = (1-\ew{}) \e{c}{(0)}(\n{}{}) + \ew{} \e{c}{(1)}(\n{}{}). \end{equation} -The reduced (i.e.~per electron) HF energy for these two states is: -\begin{subequations} -\begin{align} - \e{HF}{(0)}(n) & = \ldots, - \\ - \e{HF}{(1)}(n) & = \ldots. -\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)} ] \sin^2 \omega \, d\omega - \\ - & = \ldots, -\end{split} -\end{equation} -where $\omega$ is the interelectronic angle, $\Gamma(x)$ is the Gamma function, \cite{NISTbook} and -\begin{equation} - \psi_i(\omega) = \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) \sin^2 \omega \, d\omega - = \ldots. -\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 -The numerical values of the correlation energy for various $R$ are reported in Table \ref{tab:Ref} for the two 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. + Reduced (i.e., per electron) correlation energy $\e{c}{(I)}$ [see Eq.~\eqref{eq:ec}] as a function of $R = ...$ for the ground state ($I=0$), and the first doubly-excited state ($I=1$) of the (spin-unpolarized) two-electron glomium system. The data gathered in Table \ref{tab:Ref} are also reported. } \label{fig:Ec} @@ -240,32 +257,36 @@ The numerical values of the correlation energy for various $R$ are reported in T %%% %%% %%% %%% TABLE I %%% -\begin{turnpage} -\begin{squeezetable} -\begin{table*} +\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. + $-\e{c}{(I)}$ as a function of the radius of the glome $R$ for the ground state ($I=0$), and the first doubly-excited state ($I=1$) of the (spin-unpolarized) two-electron glomium 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$} \\ + \begin{tabular}{ldd} + $R$ & \mc{2}{c}{State } \\ + \cline{2-3} + & \tabc{Ground state} & \tabc{Doubly-excited state} \\ \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 \\ + $0$ & & \\ + $1/10$ & & \\ + $1/5$ & & \\ + $1/2$ & & \\ + $1$ & & \\ + $2$ & & \\ + $5$ & & \\ + $10$ & & \\ + $20$ & & \\ + $50$ & & \\ + $100$ & & \\ \end{tabular} \end{ruledtabular} -\end{table*} -\end{squeezetable} -\end{turnpage} +\end{table} 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{c_1^{(I)}\,n}{n + c_2^{(I)} \sqrt{n} + c_3^{(I)}}, + \e{c}{(I)}(n) = \frac{a^{(I)}\,n}{n + b^{(I)} \sqrt{n} + c^{(I)}}, \end{equation} where $c_2^{(I)}$ and $c_3^{(I)}$ are state-specific fitting parameters, which are provided in Table I of the manuscript. The value of $c_1^{(I)}$ is obtained via the exact high-density expansion of the correlation energy. \cite{Loos_2013a, Loos_2014a} @@ -274,19 +295,19 @@ Equation \eqref{eq:ec} is depicted in Fig.~\ref{fig:Ec} for each state alongside %%% TABLE 1 %%% -\begin{table} +\begin{table*} \caption{ \label{tab:OG_func} Parameters of the correlation DFAs defined in Eq.~\eqref{eq:ec}.} \begin{ruledtabular} \begin{tabular}{lcddd} - State & $I$ & \tabc{$c_1^{(I)}$} & \tabc{$c_2^{(I)}$} & \tabc{$c_3^{(I)}$} \\ + State & $I$ & \tabc{$a_1^{(I)}$} & \tabc{$a_2^{(I)}$} & \tabc{$a_3^{(I)}$} \\ \hline - Ground state & $0$ & & & \\ - Doubly-excited state & $1$ & & & \\ + Ground state & $0$ & -0.0238184 & +0.00575719 & +0.0830576 \\ + Doubly-excited state & $1$ & -0.0144633 & -0.0504501 & +0.0331287 \\ \end{tabular} \end{ruledtabular} -\end{table} +\end{table*} %%% %%% %%% %%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% @@ -299,21 +320,30 @@ In order to make the two-electron-based eDFA defined in Eq.~\eqref{eq:ecw} more \end{equation} where \begin{equation} - \be{xc}{(I)}(\n{}{}) = \e{xc}{(I)}(\n{}{}) + \e{xc}{\text{LDA}}(\n{}{}) - \e{xc}{(0)}(\n{}{}). + \be{xc}{(I)}(\n{}{}) = \e{xc}{(I)}(\n{}{}) + \e{xc}{\LDA}(\n{}{}) - \e{xc}{(0)}(\n{}{}). \end{equation} The local-density approximation (LDA) exchange-correlation functional is \begin{equation} - \e{xc}{\text{LDA}}(\n{}{}) = \e{x}{\text{LDA}}(\n{}{}) + \e{c}{\text{LDA}}(\n{}{}). + \e{xc}{\LDA}(\n{}{}) = \e{x}{\LDA}(\n{}{}) + \e{c}{\LDA}(\n{}{}). \end{equation} - +where we use here the Dirac exchange functional and the VWN5 correlation functional +\begin{subequations} +\begin{align} + \e{x}{\LDA}(\n{}{}) & = \Cx{\LDA} \n{}{1/3} + \\ + \e{c}{\LDA}(\n{}{}) & \equiv \e{c}{\text{VWN5}}(\n{}{}). +\end{align} +\end{subequations} +with $\Cx{\LDA} = -\frac{3}{2} \qty(\frac{3}{4\pi})^{1/3}$. + Equation \eqref{eq:becw} can be recast \begin{equation} \label{eq:eLDA} \be{xc}{\ew{}}(\n{}{}) - = \e{xc}{\text{LDA}}(\n{}{}) + \ew{} \qty[\e{xc}{(1)}(\n{}{})-\e{xc}{(0)}(\n{}{})], + = \e{xc}{\LDA}(\n{}{}) + \ew{} \qty[\e{xc}{(1)}(\n{}{})-\e{xc}{(0)}(\n{}{})], \end{equation} which nicely highlights the centrality of the LDA in the present eDFA. -In particular, $\be{xc}{(0)}(\n{}{}) = \e{xc}{\text{LDA}}(\n{}{})$. +In particular, $\be{xc}{(0)}(\n{}{}) = \e{xc}{\LDA}(\n{}{})$. Consequently, in the following, we name this correlation functional ``eLDA'' as it is a natural extension of the LDA for ensembles. This procedure can be theoretically justified by the generalized adiabatic connection formalism for ensembles (GACE) which was originally derived by Franck and Fromager. \cite{Franck_2014} @@ -324,18 +354,19 @@ Finally, we note that, by construction, \begin{equation} \left. \pdv{\be{xc}{\ew{}}[\n{}{}]}{\ew{J}}\right|_{\n{}{} = \n{}{\ew{}}(\br)} = \be{xc}{(J)}[\n{}{\ew{}}(\br)] - \be{xc}{(0)}[\n{}{\ew{}}(\br)]. \end{equation} -\alert{As shown by Gould and Pittalis, comment on density- and and state-driven errors. \cite{Gould_2019}} - - - -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 glomium model. -We refer the interested reader to Refs.~\onlinecite{} for more details about this paradigm. +%%%%%%%%%%%%%%% +%%% RESULTS %%% +%%%%%%%%%%%%%%% +\section{Results} +\label{sec:res} +Here, we do \ce{H2} because \ce{H2} is very interesting. %%%%%%%%%%%%%%%%%% %%% CONCLUSION %%% %%%%%%%%%%%%%%%%%% \section{Conclusion} +\label{sec:ccl} As concluding remarks, we would like to say that, what we have done is awesome. %%%%%%%%%%%%%%%%%%%%%%%%