From f3b5f49c150495a796544a7130b60efa9517b047 Mon Sep 17 00:00:00 2001 From: Pierre-Francois Loos Date: Fri, 14 Feb 2020 10:53:16 +0100 Subject: [PATCH] save work --- FarDFT.nb | 1551 +++++++++++++++++++++++++++++------------ Manuscript/FarDFT.tex | 4 +- 2 files changed, 1103 insertions(+), 452 deletions(-) diff --git a/FarDFT.nb b/FarDFT.nb index 7fc7bd1..42992d0 100644 --- a/FarDFT.nb +++ b/FarDFT.nb @@ -10,10 +10,10 @@ NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 158, 7] -NotebookDataLength[ 196078, 5262] -NotebookOptionsPosition[ 178259, 4982] -NotebookOutlinePosition[ 178596, 4997] -CellTagsIndexPosition[ 178553, 4994] +NotebookDataLength[ 219832, 5913] +NotebookOptionsPosition[ 200640, 5612] +NotebookOutlinePosition[ 200977, 5627] +CellTagsIndexPosition[ 200934, 5624] WindowFrame->Normal*) (* Beginning of Notebook Content *) @@ -98,7 +98,7 @@ Cell[BoxData[ InitializationCell->True, CellChangeTimes->{{3.790612553197805*^9, 3.790612558626453*^9}}, CellLabel->"In[7]:=",ExpressionUUID->"5dbbd40c-d911-445f-b476-85cd91f8d014"] -}, Open ]], +}, Closed]], Cell[CellGroupData[{ @@ -366,7 +366,7 @@ Cell[BoxData[ 3.790614210485475*^9}}, CellLabel->"In[11]:=",ExpressionUUID->"440fb372-2001-4d90-ade1-e5498676702f"] }, Open ]] -}, Open ]] +}, Closed]] }, Open ]], Cell[CellGroupData[{ @@ -2920,9 +2920,47 @@ Cell[BoxData[{ CellChangeTimes->{{3.79059838113128*^9, 3.7905985061910954`*^9}, 3.7906057050230303`*^9, {3.7906109479713078`*^9, 3.790610952959711*^9}, { 3.790617319564548*^9, 3.790617416879177*^9}, {3.790618439965136*^9, - 3.7906184689663363`*^9}}, - CellLabel-> - "In[154]:=",ExpressionUUID->"f51f78c8-fe8e-48f3-a562-1209bfc50f71"], + 3.7906184689663363`*^9}, {3.79066088543036*^9, + 3.790660891318091*^9}},ExpressionUUID->"f51f78c8-fe8e-48f3-a562-\ +1209bfc50f71"], + +Cell[BoxData[{ + RowBox[{ + RowBox[{ + RowBox[{"d\[Epsilon]xLDAdn", "[", "n_", "]"}], ":=", + RowBox[{"D", "[", + RowBox[{ + RowBox[{"\[Epsilon]xLDA", "[", "n", "]"}], ",", "n"}], "]"}]}], + ";"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{ + RowBox[{"d\[Epsilon]x0dn", "[", "n_", "]"}], ":=", + RowBox[{"D", "[", + RowBox[{ + RowBox[{"\[Epsilon]x0", "[", "n", "]"}], ",", "n"}], "]"}]}], + ";"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"d\[Epsilon]x1dn", "[", "n_", "]"}], ":=", + RowBox[{"D", "[", + RowBox[{ + RowBox[{"\[Epsilon]x1", "[", "n", "]"}], ",", "n"}], + "]"}]}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"d\[Epsilon]xeLDAdn", "[", + RowBox[{"n_", ",", "w_"}], "]"}], ":=", + RowBox[{ + RowBox[{"d\[Epsilon]xLDAdn", "[", "n", "]"}], "+", + RowBox[{"w", + RowBox[{"(", + RowBox[{ + RowBox[{"d\[Epsilon]x1dn", "[", "n", "]"}], "-", + RowBox[{"d\[Epsilon]x0dn", "[", "n", "]"}]}], ")"}]}]}]}]}], "Input", + InitializationCell->True, + CellChangeTimes->{{3.79059838113128*^9, 3.7905985061910954`*^9}, + 3.7906057050230303`*^9, {3.7906109479713078`*^9, 3.790610952959711*^9}, { + 3.790617319564548*^9, 3.790617416879177*^9}, {3.790618439965136*^9, + 3.7906184689663363`*^9}, {3.79066088543036*^9, 3.7906609462132797`*^9}}, + CellLabel->"In[85]:=",ExpressionUUID->"dd062623-5a91-460e-945c-d5f15e24abfb"], Cell["LDA correlation functionals", "Text", CellChangeTimes->{{3.790617745116544*^9, 3.790617768448262*^9}}, @@ -3096,8 +3134,7 @@ Cell[BoxData[ 3.783426658830944*^9}, {3.783427391840214*^9, 3.783427392389812*^9}, 3.7834287156166573`*^9, {3.7836773637524137`*^9, 3.7836775311274233`*^9}, { 3.7906176573902884`*^9, 3.790617676290105*^9}}, - CellLabel-> - "In[158]:=",ExpressionUUID->"c5e15fc2-7ce6-4ed3-bc7f-bcdc9085134d"], + CellLabel->"In[16]:=",ExpressionUUID->"c5e15fc2-7ce6-4ed3-bc7f-bcdc9085134d"], Cell[BoxData[{ RowBox[{ @@ -3168,8 +3205,45 @@ Cell[BoxData[{ 3.790617319564548*^9, 3.790617516913457*^9}, {3.790617702215055*^9, 3.790617719701961*^9}, {3.7906182605153427`*^9, 3.790618262031168*^9}, { 3.790618364188952*^9, 3.790618401731085*^9}}, - CellLabel-> - "In[164]:=",ExpressionUUID->"8db93606-ec73-4fe8-9228-51f79c4ef957"], + CellLabel->"In[22]:=",ExpressionUUID->"8db93606-ec73-4fe8-9228-51f79c4ef957"], + +Cell[BoxData[{ + RowBox[{ + RowBox[{ + RowBox[{"d\[Epsilon]cLDAdn", "[", "n_", "]"}], ":=", + RowBox[{"D", "[", + RowBox[{ + RowBox[{"\[Epsilon]cLDA", "[", "n", "]"}], ",", "n"}], "]"}]}], + ";"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{ + RowBox[{"d\[Epsilon]c0dn", "[", "n_", "]"}], ":=", + RowBox[{"D", "[", + RowBox[{ + RowBox[{"\[Epsilon]c0", "[", "n", "]"}], ",", "n"}], "]"}]}], + ";"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"d\[Epsilon]c1dn", "[", "n_", "]"}], ":=", + RowBox[{"D", "[", + RowBox[{ + RowBox[{"\[Epsilon]c1", "[", "n", "]"}], ",", "n"}], + "]"}]}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"d\[Epsilon]ceLDAdn", "[", + RowBox[{"n_", ",", "w_"}], "]"}], ":=", + RowBox[{ + RowBox[{"d\[Epsilon]cLDAdn", "[", "n", "]"}], "+", + RowBox[{"w", + RowBox[{"(", + RowBox[{ + RowBox[{"d\[Epsilon]c1dn", "[", "n", "]"}], "-", + RowBox[{"d\[Epsilon]c0dn", "[", "n", "]"}]}], ")"}]}]}]}]}], "Input", + InitializationCell->True, + CellChangeTimes->{{3.79059838113128*^9, 3.7905985061910954`*^9}, + 3.7906057050230303`*^9, {3.7906109479713078`*^9, 3.790610952959711*^9}, { + 3.790617319564548*^9, 3.790617416879177*^9}, {3.790618439965136*^9, + 3.7906184689663363`*^9}, {3.79066088543036*^9, 3.790660978615497*^9}}, + CellLabel->"In[90]:=",ExpressionUUID->"3c5be8aa-305f-4c0e-a51a-4bc92cdcf171"], Cell["LDA exchange-correlation functionals", "Text", CellChangeTimes->{{3.790617745116544*^9, 3.790617768448262*^9}, { @@ -3205,8 +3279,63 @@ Cell[BoxData[{ RowBox[{"n", ",", "w"}], "]"}]}]}]}], "Input", InitializationCell->True, CellChangeTimes->{{3.790617806453972*^9, 3.7906178574975023`*^9}}, + CellLabel->"In[26]:=",ExpressionUUID->"f6517707-c164-4ea9-b57c-8d4d27eafad3"], + +Cell[BoxData[{ + RowBox[{ + RowBox[{"d\[Epsilon]xcLDAdn", "[", "n_", "]"}], ":=", + RowBox[{ + RowBox[{"d\[Epsilon]xLDAdn", "[", "n", "]"}], "+", + RowBox[{"d\[Epsilon]cLDAdn", "[", "n", "]"}]}]}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{ + RowBox[{"d\[Epsilon]xc0dn", "[", "n_", "]"}], ":=", + RowBox[{ + RowBox[{"d\[Epsilon]x0dn", "[", "n", "]"}], "+", + RowBox[{"d\[Epsilon]c0dn", "[", "n", "]"}]}]}], + ";"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{ + RowBox[{"d\[Epsilon]xc1dn", "[", "n_", "]"}], ":=", + RowBox[{ + RowBox[{"d\[Epsilon]x1dn", "[", "n", "]"}], "+", + RowBox[{"d\[Epsilon]c1dn", "[", "n", "]"}]}]}], + ";"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"d\[Epsilon]xceLDAdn", "[", + RowBox[{"n_", ",", "w_"}], "]"}], ":=", + RowBox[{ + RowBox[{"d\[Epsilon]xeLDAdn", "[", + RowBox[{"n", ",", "w"}], "]"}], "+", + RowBox[{"d\[Epsilon]ceLDAdn", "[", + RowBox[{"n", ",", "w"}], "]"}]}]}]}], "Input", + InitializationCell->True, + CellChangeTimes->{{3.790617806453972*^9, 3.7906178574975023`*^9}, { + 3.790660991606427*^9, 3.790661018422545*^9}}, + CellLabel->"In[94]:=",ExpressionUUID->"33dfcd59-33cc-4f88-bc22-62c17b91dbeb"], + +Cell[BoxData[{ + RowBox[{ + RowBox[{"d\[Epsilon]xeLDAdw", "[", "n_", "]"}], ":=", + RowBox[{ + RowBox[{"\[Epsilon]x1", "[", "n", "]"}], "-", + RowBox[{"\[Epsilon]x0", "[", "n", "]"}]}]}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"d\[Epsilon]ceLDAdw", "[", "n_", "]"}], ":=", + RowBox[{ + RowBox[{"\[Epsilon]c1", "[", "n", "]"}], "-", + RowBox[{"\[Epsilon]c0", "[", "n", "]"}]}]}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"d\[Epsilon]xceLDAdw", "[", "n_", "]"}], ":=", + RowBox[{ + RowBox[{"\[Epsilon]xc1", "[", "n", "]"}], "-", + RowBox[{"\[Epsilon]xc0", "[", "n", "]"}]}]}]}], "Input", + InitializationCell->True, + CellChangeTimes->{{3.7906610362852507`*^9, 3.790661042510214*^9}, { + 3.790661095924047*^9, 3.790661108253295*^9}, {3.790661781802154*^9, + 3.790661784729022*^9}}, CellLabel-> - "In[168]:=",ExpressionUUID->"f6517707-c164-4ea9-b57c-8d4d27eafad3"], + "In[131]:=",ExpressionUUID->"9663902a-8667-481b-b613-f1c2b56b8e76"], Cell[TextData[StyleBox["Molecular coordinates", FontColor->RGBColor[1, 0, 0]]], "Text", @@ -3233,7 +3362,7 @@ Cell[BoxData[ InitializationCell->True, CellChangeTimes->{{3.790612571234482*^9, 3.790612583186882*^9}, { 3.790612652169979*^9, 3.790612655860166*^9}}, - CellLabel->"In[14]:=",ExpressionUUID->"36b985d7-30a2-430d-8dc3-30d6b4a45727"], + CellLabel->"In[30]:=",ExpressionUUID->"36b985d7-30a2-430d-8dc3-30d6b4a45727"], Cell[TextData[StyleBox["Basis functions", FontColor->RGBColor[1, 0, 0]]], "Text", @@ -3333,7 +3462,7 @@ Cell[BoxData[ 3.790612947951182*^9, 3.790612965035768*^9}, {3.790613311831658*^9, 3.790613358728168*^9}, {3.790613507527297*^9, 3.79061350977219*^9}, { 3.790614141437229*^9, 3.7906141534618187`*^9}}, - CellLabel->"In[15]:=",ExpressionUUID->"f0a2825f-92fc-42d2-aeee-57f53535b1d5"], + CellLabel->"In[31]:=",ExpressionUUID->"f0a2825f-92fc-42d2-aeee-57f53535b1d5"], Cell[TextData[StyleBox["Integrals over basis functions", FontColor->RGBColor[1, 0, 0]]], "Text", @@ -4032,7 +4161,7 @@ Cell[BoxData[ 3.790590259300644*^9, 3.79059026249852*^9}, {3.790590572948614*^9, 3.790590598010453*^9}, {3.79061268483996*^9, 3.790612743776194*^9}, { 3.790614708170142*^9, 3.790614712227317*^9}}, - CellLabel->"In[16]:=",ExpressionUUID->"341cb469-2a0b-4b3f-a5c7-4a3d8fb61885"], + CellLabel->"In[32]:=",ExpressionUUID->"341cb469-2a0b-4b3f-a5c7-4a3d8fb61885"], Cell[TextData[StyleBox["Spatial orbitals in minimal basis:", FontColor->RGBColor[1, 0, 0]]], "Text", @@ -4097,70 +4226,7 @@ Cell[BoxData[ 3.7905859157992983`*^9}, {3.790587433782695*^9, 3.790587438383739*^9}, { 3.790613976434304*^9, 3.79061399579806*^9}, {3.790614159389223*^9, 3.7906141606657143`*^9}}, - CellLabel->"In[17]:=",ExpressionUUID->"3bf77406-a77c-45d7-8c2f-817ad440cdff"], - -Cell[CellGroupData[{ - -Cell[BoxData[{ - RowBox[{"NIntegrate", "[", - RowBox[{ - SuperscriptBox[ - RowBox[{ - SubscriptBox["\[Phi]", "1"], "[", - RowBox[{ - RowBox[{"{", - RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "1"}], "]"}], "2"], ",", - RowBox[{"{", - RowBox[{"x", ",", - RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", - RowBox[{"{", - RowBox[{"y", ",", - RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", - RowBox[{"{", - RowBox[{"z", ",", - RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}]}], - "]"}], "\[IndentingNewLine]", - RowBox[{"NIntegrate", "[", - RowBox[{ - SuperscriptBox[ - RowBox[{ - SubscriptBox["\[Phi]", "2"], "[", - RowBox[{ - RowBox[{"{", - RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "1"}], "]"}], "2"], ",", - RowBox[{"{", - RowBox[{"x", ",", - RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", - RowBox[{"{", - RowBox[{"y", ",", - RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", - RowBox[{"{", - RowBox[{"z", ",", - RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}]}], - "]"}]}], "Input", - CellChangeTimes->{{3.790613009939126*^9, 3.790613015432294*^9}, { - 3.790613377187613*^9, 3.790613388996928*^9}, {3.790613521860344*^9, - 3.790613523881935*^9}, {3.79061416814648*^9, 3.7906141686251717`*^9}}, - CellLabel->"In[80]:=",ExpressionUUID->"92574105-d574-49ac-904d-a9f7b38ef1c3"], - -Cell[BoxData["0.9999999978278258`"], "Output", - CellChangeTimes->{ - 3.7906129813068457`*^9, 3.790613015763709*^9, 3.7906131807387753`*^9, - 3.790613237870916*^9, 3.7906133237388697`*^9, {3.790613363224881*^9, - 3.7906133898515778`*^9}, {3.79061350242299*^9, 3.7906135247380457`*^9}, - 3.790613611708521*^9, 3.790613943806418*^9, 3.7906139978828173`*^9, { - 3.790614169581675*^9, 3.790614190198777*^9}, 3.790614228992675*^9}, - CellLabel->"Out[80]=",ExpressionUUID->"aa0b75e7-d797-49d3-9a12-b14fac9e91cc"], - -Cell[BoxData["1.000000001455642`"], "Output", - CellChangeTimes->{ - 3.7906129813068457`*^9, 3.790613015763709*^9, 3.7906131807387753`*^9, - 3.790613237870916*^9, 3.7906133237388697`*^9, {3.790613363224881*^9, - 3.7906133898515778`*^9}, {3.79061350242299*^9, 3.7906135247380457`*^9}, - 3.790613611708521*^9, 3.790613943806418*^9, 3.7906139978828173`*^9, { - 3.790614169581675*^9, 3.790614190198777*^9}, 3.790614229974461*^9}, - CellLabel->"Out[81]=",ExpressionUUID->"9defd3f6-f5d8-462a-a8a4-c0bc80c8daea"] -}, Open ]], + CellLabel->"In[33]:=",ExpressionUUID->"3bf77406-a77c-45d7-8c2f-817ad440cdff"], Cell[TextData[StyleBox["One-electron densities", FontColor->RGBColor[1, 0, 0]]], "Text", @@ -4205,69 +4271,7 @@ Cell[BoxData[ 3.7905859157992983`*^9}, {3.790587433782695*^9, 3.790587438383739*^9}, { 3.790605575633658*^9, 3.790605631352661*^9}, {3.790611080408037*^9, 3.790611083617523*^9}}, - CellLabel->"In[18]:=",ExpressionUUID->"4bf40614-2147-47eb-9ce0-13fe949a9c00"], - -Cell[CellGroupData[{ - -Cell[BoxData[{ - RowBox[{"With", "[", - RowBox[{ - RowBox[{"{", - RowBox[{"R", "=", "1"}], "}"}], ",", - RowBox[{"NIntegrate", "[", - RowBox[{ - RowBox[{ - SubscriptBox["n", "1"], "[", - RowBox[{ - RowBox[{"{", - RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R"}], "]"}], ",", - RowBox[{"{", - RowBox[{"x", ",", - RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", - RowBox[{"{", - RowBox[{"y", ",", - RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", - RowBox[{"{", - RowBox[{"z", ",", - RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}]}], "]"}]}], - "]"}], "\[IndentingNewLine]", - RowBox[{"With", "[", - RowBox[{ - RowBox[{"{", - RowBox[{"R", "=", "1"}], "}"}], ",", - RowBox[{"NIntegrate", "[", - RowBox[{ - RowBox[{ - SubscriptBox["n", "2"], "[", - RowBox[{ - RowBox[{"{", - RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R"}], "]"}], ",", - RowBox[{"{", - RowBox[{"x", ",", - RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", - RowBox[{"{", - RowBox[{"y", ",", - RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", - RowBox[{"{", - RowBox[{"z", ",", - RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}]}], "]"}]}], - "]"}]}], "Input", - CellChangeTimes->{{3.790611087595409*^9, 3.790611117851706*^9}, { - 3.790611266223362*^9, 3.790611269262435*^9}}, - CellLabel->"In[85]:=",ExpressionUUID->"246fa925-65a3-41a1-b83a-54965c0376f6"], - -Cell[BoxData["1.9999999956556516`"], "Output", - CellChangeTimes->{{3.79061110341438*^9, 3.790611118269318*^9}, - 3.79061114928259*^9, 3.790611270030858*^9, 3.790612995128433*^9, - 3.7906132277625103`*^9, {3.790614234007246*^9, 3.7906142375844183`*^9}}, - CellLabel->"Out[85]=",ExpressionUUID->"e0356b94-048c-4f45-9ac7-db1db0236396"], - -Cell[BoxData["2.0000000029112863`"], "Output", - CellChangeTimes->{{3.79061110341438*^9, 3.790611118269318*^9}, - 3.79061114928259*^9, 3.790611270030858*^9, 3.790612995128433*^9, - 3.7906132277625103`*^9, {3.790614234007246*^9, 3.79061423853084*^9}}, - CellLabel->"Out[86]=",ExpressionUUID->"2b09c0f7-0efb-4358-85fb-939d40cdda47"] -}, Open ]], + CellLabel->"In[34]:=",ExpressionUUID->"4bf40614-2147-47eb-9ce0-13fe949a9c00"], Cell[BoxData[ RowBox[{ @@ -4299,8 +4303,7 @@ Cell[BoxData[ 3.7905859157992983`*^9}, {3.790587433782695*^9, 3.790587438383739*^9}, { 3.790605575633658*^9, 3.790605631352661*^9}, {3.790611080408037*^9, 3.790611083617523*^9}, {3.790619238747375*^9, 3.790619276696437*^9}}, - CellLabel-> - "In[192]:=",ExpressionUUID->"33c091a9-0792-4be2-bcd6-e02a256b26b5"], + CellLabel->"In[35]:=",ExpressionUUID->"33c091a9-0792-4be2-bcd6-e02a256b26b5"], Cell[TextData[StyleBox["Integrals over orbitals", FontColor->RGBColor[1, 0, 0]]], "Text", @@ -4358,7 +4361,7 @@ Cell[BoxData[ 3.790587455877968*^9, 3.790587477805868*^9}, {3.7905875469262457`*^9, 3.790587604542864*^9}, {3.790587670334567*^9, 3.790587690230748*^9}, { 3.790590272880838*^9, 3.790590291365608*^9}, 3.7905903766007347`*^9}, - CellLabel->"In[30]:=",ExpressionUUID->"41a0d75d-44b6-4653-b70e-858ab3e937b8"], + CellLabel->"In[36]:=",ExpressionUUID->"41a0d75d-44b6-4653-b70e-858ab3e937b8"], Cell[BoxData[{ RowBox[{ @@ -4450,7 +4453,7 @@ Cell[BoxData[{ 3.790588310918952*^9, 3.790588322533786*^9}, {3.790588364656127*^9, 3.7905883758140497`*^9}, {3.790588409551832*^9, 3.7905884242861843`*^9}, { 3.790588557126391*^9, 3.7905885830952806`*^9}}, - CellLabel->"In[20]:=",ExpressionUUID->"1430ab04-b05a-4d2d-bbf6-f9edc3cd5345"], + CellLabel->"In[37]:=",ExpressionUUID->"1430ab04-b05a-4d2d-bbf6-f9edc3cd5345"], Cell[BoxData[{ RowBox[{ @@ -4497,7 +4500,119 @@ Cell[BoxData[{ InitializationCell->True, CellChangeTimes->{{3.790588642589957*^9, 3.790588682679491*^9}, { 3.790588776310627*^9, 3.790588778758974*^9}}, - CellLabel->"In[23]:=",ExpressionUUID->"09f7d062-2fdd-45b1-8df2-0dd679b5cbc6"], + CellLabel->"In[40]:=",ExpressionUUID->"09f7d062-2fdd-45b1-8df2-0dd679b5cbc6"], + +Cell[TextData[StyleBox["Nuclear energy", + FontColor->RGBColor[1, 0, 0]]], "Text", + CellChangeTimes->{{3.790587253093514*^9, 3.790587257373096*^9}, { + 3.7906590777707577`*^9, + 3.790659080098378*^9}},ExpressionUUID->"94763911-d015-4b98-a645-\ +b27f631322bf"], + +Cell[BoxData[ + RowBox[{ + RowBox[{ + RowBox[{ + SubscriptBox["\[CapitalEpsilon]", "NN"], "[", "R_", "]"}], "=", + FractionBox[ + RowBox[{"ZA", " ", "ZB"}], + RowBox[{"Norm", "[", + RowBox[{ + RowBox[{"rA", "[", "R", "]"}], "-", + RowBox[{"rB", "[", "R", "]"}]}], "]"}]]}], ";"}]], "Input", + InitializationCell->True, + CellChangeTimes->{{3.790585364582788*^9, 3.790585422341731*^9}, { + 3.790588811486704*^9, 3.7905888669431877`*^9}, {3.7906190335224743`*^9, + 3.7906190337111177`*^9}, {3.790619093049782*^9, 3.790619093357505*^9}, { + 3.7906590832436953`*^9, 3.790659125362174*^9}}, + CellLabel->"In[66]:=",ExpressionUUID->"46e4593c-f9e1-4f4a-a999-d0f67e325be6"], + +Cell[TextData[StyleBox["CID energies", + FontColor->RGBColor[1, 0, 0]]], "Text", + CellChangeTimes->{{3.790587253093514*^9, 3.790587257373096*^9}, { + 3.79060545274035*^9, + 3.790605453095336*^9}},ExpressionUUID->"eeb33af8-ab09-4c5f-98a3-\ +684ae529c2b6"], + +Cell[BoxData[ + RowBox[{ + RowBox[{ + RowBox[{ + SubscriptBox["\[CapitalEpsilon]", "CID"], "[", "R_", "]"}], "=", + RowBox[{ + RowBox[{"{", + RowBox[{ + RowBox[{ + FractionBox[ + RowBox[{ + RowBox[{ + RowBox[{ + SubscriptBox["\[CapitalEpsilon]", "HF"], "[", "R", "]"}], + "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], "+", + RowBox[{ + RowBox[{ + SubscriptBox["\[CapitalEpsilon]", "HF"], "[", "R", "]"}], + "\[LeftDoubleBracket]", "2", "\[RightDoubleBracket]"}]}], "2"], "-", + RowBox[{ + FractionBox["1", "2"], + SqrtBox[ + RowBox[{ + SuperscriptBox[ + RowBox[{"(", + RowBox[{ + RowBox[{ + RowBox[{ + SubscriptBox["\[CapitalEpsilon]", "HF"], "[", "R", "]"}], + "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], "-", + RowBox[{ + RowBox[{ + SubscriptBox["\[CapitalEpsilon]", "HF"], "[", "R", "]"}], + "\[LeftDoubleBracket]", "2", "\[RightDoubleBracket]"}]}], ")"}], + "2"], "+", + RowBox[{"4", + SuperscriptBox[ + RowBox[{ + SubscriptBox["K", "12"], "[", "R", "]"}], "2"]}]}]]}]}], ",", + RowBox[{ + FractionBox[ + RowBox[{ + RowBox[{ + RowBox[{ + SubscriptBox["\[CapitalEpsilon]", "HF"], "[", "R", "]"}], + "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], "+", + RowBox[{ + RowBox[{ + SubscriptBox["\[CapitalEpsilon]", "HF"], "[", "R", "]"}], + "\[LeftDoubleBracket]", "2", "\[RightDoubleBracket]"}]}], "2"], "+", + RowBox[{ + FractionBox["1", "2"], + SqrtBox[ + RowBox[{ + SuperscriptBox[ + RowBox[{"(", + RowBox[{ + RowBox[{ + RowBox[{ + SubscriptBox["\[CapitalEpsilon]", "HF"], "[", "R", "]"}], + "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], "-", + RowBox[{ + RowBox[{ + SubscriptBox["\[CapitalEpsilon]", "HF"], "[", "R", "]"}], + "\[LeftDoubleBracket]", "2", "\[RightDoubleBracket]"}]}], ")"}], + "2"], "+", + RowBox[{"4", + SuperscriptBox[ + RowBox[{ + SubscriptBox["K", "12"], "[", "R", "]"}], "2"]}]}]]}]}]}], "}"}], + "+", + RowBox[{ + SubscriptBox["\[CapitalEpsilon]", "NN"], "[", "R", "]"}]}]}], + ";"}]], "Input", + InitializationCell->True, + CellChangeTimes->{{3.790605458040031*^9, 3.790605458543043*^9}, { + 3.790605491525733*^9, 3.79060554005932*^9}, {3.790614460047532*^9, + 3.790614470197546*^9}, 3.7906591751978273`*^9}, + CellLabel->"In[69]:=",ExpressionUUID->"ab8163b6-fc27-4314-ae7e-a7f0f7d95629"], Cell[TextData[StyleBox["Hartree-Fock energies", FontColor->RGBColor[1, 0, 0]]], "Text", @@ -4510,32 +4625,38 @@ Cell[BoxData[ RowBox[{ RowBox[{ SubscriptBox["\[CapitalEpsilon]", "HF"], "[", "R_", "]"}], "=", - RowBox[{"{", + RowBox[{ + RowBox[{"{", "\[IndentingNewLine]", + RowBox[{ + RowBox[{ + RowBox[{"2", + RowBox[{ + SubscriptBox["h", "1"], "[", "R", "]"}]}], "+", + RowBox[{"2", + RowBox[{ + SubscriptBox["J", "11"], "[", "R", "]"}]}], "-", + RowBox[{ + SubscriptBox["K", "11"], "[", "R", "]"}]}], ",", + "\[IndentingNewLine]", + RowBox[{ + RowBox[{"2", + RowBox[{ + SubscriptBox["h", "2"], "[", "R", "]"}]}], "+", + RowBox[{"2", + RowBox[{ + SubscriptBox["J", "22"], "[", "R", "]"}]}], "-", + RowBox[{ + SubscriptBox["K", "22"], "[", "R", "]"}]}]}], "\[IndentingNewLine]", + "}"}], "+", RowBox[{ - RowBox[{ - RowBox[{"2", - RowBox[{ - SubscriptBox["h", "1"], "[", "R", "]"}]}], "+", - RowBox[{"2", - RowBox[{ - SubscriptBox["J", "11"], "[", "R", "]"}]}], "-", - RowBox[{ - SubscriptBox["K", "11"], "[", "R", "]"}]}], ",", - RowBox[{ - RowBox[{"2", - RowBox[{ - SubscriptBox["h", "2"], "[", "R", "]"}]}], "+", - RowBox[{"2", - RowBox[{ - SubscriptBox["J", "22"], "[", "R", "]"}]}], "-", - RowBox[{ - SubscriptBox["K", "22"], "[", "R", "]"}]}]}], "}"}]}], ";"}]], "Input",\ - + SubscriptBox["\[CapitalEpsilon]", "NN"], "[", "R", "]"}]}]}], + ";"}]], "Input", InitializationCell->True, CellChangeTimes->{{3.790585364582788*^9, 3.790585422341731*^9}, { 3.790588811486704*^9, 3.7905888669431877`*^9}, {3.7906190335224743`*^9, - 3.7906190337111177`*^9}, {3.790619093049782*^9, - 3.790619093357505*^9}},ExpressionUUID->"78388547-17f1-4781-bf29-\ + 3.7906190337111177`*^9}, {3.790619093049782*^9, 3.790619093357505*^9}, { + 3.79065915742903*^9, 3.790659168377755*^9}, {3.790660617755475*^9, + 3.790660622626095*^9}},ExpressionUUID->"78388547-17f1-4781-bf29-\ 2487d17b5c53"], Cell[BoxData[ @@ -4582,96 +4703,96 @@ Cell[BoxData[ RowBox[{ SubscriptBox["J", "12"], "[", "R", "]"}]}], "-", RowBox[{ - SubscriptBox["K", "12"], "[", "R", "]"}]}], ")"}]}]}]}], + SubscriptBox["K", "12"], "[", "R", "]"}]}], ")"}]}], "+", + RowBox[{ + SubscriptBox["\[CapitalEpsilon]", "NN"], "[", "R", "]"}]}]}], ";"}]], "Input", InitializationCell->True, - CellLabel-> - "In[183]:=",ExpressionUUID->"1039bc21-357b-49d9-a514-fd3af5ae05a7"], - -Cell[TextData[StyleBox["CID energies", - FontColor->RGBColor[1, 0, 0]]], "Text", - CellChangeTimes->{{3.790587253093514*^9, 3.790587257373096*^9}, { - 3.79060545274035*^9, - 3.790605453095336*^9}},ExpressionUUID->"eeb33af8-ab09-4c5f-98a3-\ -684ae529c2b6"], + CellChangeTimes->{3.790659173613167*^9}, + CellLabel->"In[68]:=",ExpressionUUID->"1039bc21-357b-49d9-a514-fd3af5ae05a7"], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{ - SubscriptBox["\[CapitalEpsilon]", "CID"], "[", "R_", "]"}], "=", + SubscriptBox["\[CapitalDelta]\[CapitalEpsilon]w", "HF"], "[", + RowBox[{"R_", ",", "w_"}], "]"}], "=", + RowBox[{ + RowBox[{ + RowBox[{"(", + RowBox[{"1", "-", "w"}], ")"}], + RowBox[{ + RowBox[{ + SubscriptBox["\[CapitalEpsilon]", "HF"], "[", "R", "]"}], + "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}]}], "+", + RowBox[{"w", + RowBox[{ + RowBox[{ + SubscriptBox["\[CapitalEpsilon]", "HF"], "[", "R", "]"}], + "\[LeftDoubleBracket]", "2", "\[RightDoubleBracket]"}]}], "-", + RowBox[{ + SubscriptBox["\[CapitalEpsilon]w", "HF"], "[", + RowBox[{"R", ",", "w"}], "]"}]}]}], ";"}]], "Input", + InitializationCell->True, + CellChangeTimes->{{3.790660001969771*^9, 3.790660072444165*^9}}, + CellLabel->"In[79]:=",ExpressionUUID->"01749816-6d54-416a-ab4f-bab2e420acb4"], + +Cell[BoxData[ + RowBox[{ + RowBox[{ + RowBox[{ + SubscriptBox["\[Epsilon]w", "HF"], "[", + RowBox[{"R_", ",", "w_"}], "]"}], "=", RowBox[{"{", RowBox[{ RowBox[{ - FractionBox[ - RowBox[{ - RowBox[{ - RowBox[{ - SubscriptBox["\[CapitalEpsilon]", "HF"], "[", "R", "]"}], - "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], "+", - RowBox[{ - RowBox[{ - SubscriptBox["\[CapitalEpsilon]", "HF"], "[", "R", "]"}], - "\[LeftDoubleBracket]", "2", "\[RightDoubleBracket]"}]}], "2"], "-", - RowBox[{ - FractionBox["1", "2"], - SqrtBox[ + SubscriptBox["h", "1"], "[", "R", "]"}], "+", + RowBox[{ + RowBox[{"(", + RowBox[{"1", "-", "w"}], ")"}], + RowBox[{"(", RowBox[{ - SuperscriptBox[ - RowBox[{"(", - RowBox[{ - RowBox[{ - RowBox[{ - SubscriptBox["\[CapitalEpsilon]", "HF"], "[", "R", "]"}], - "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], "-", - RowBox[{ - RowBox[{ - SubscriptBox["\[CapitalEpsilon]", "HF"], "[", "R", "]"}], - "\[LeftDoubleBracket]", "2", "\[RightDoubleBracket]"}]}], ")"}], - "2"], "+", - RowBox[{"4", - SuperscriptBox[ - RowBox[{ - SubscriptBox["K", "12"], "[", "R", "]"}], "2"]}]}]]}]}], ",", + RowBox[{"2", + RowBox[{ + SubscriptBox["J", "11"], "[", "R", "]"}]}], "-", + RowBox[{ + SubscriptBox["K", "11"], "[", "R", "]"}]}], ")"}]}], "+", + RowBox[{"w", + RowBox[{"(", + RowBox[{ + RowBox[{"2", + RowBox[{ + SubscriptBox["J", "12"], "[", "R", "]"}]}], "-", + RowBox[{ + SubscriptBox["K", "12"], "[", "R", "]"}]}], ")"}]}]}], ",", RowBox[{ - FractionBox[ - RowBox[{ - RowBox[{ - RowBox[{ - SubscriptBox["\[CapitalEpsilon]", "HF"], "[", "R", "]"}], - "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], "+", - RowBox[{ - RowBox[{ - SubscriptBox["\[CapitalEpsilon]", "HF"], "[", "R", "]"}], - "\[LeftDoubleBracket]", "2", "\[RightDoubleBracket]"}]}], "2"], "+", - RowBox[{ - FractionBox["1", "2"], - SqrtBox[ + SubscriptBox["h", "2"], "[", "R", "]"}], "+", + RowBox[{ + RowBox[{"(", + RowBox[{"1", "-", "w"}], ")"}], + RowBox[{"(", RowBox[{ - SuperscriptBox[ - RowBox[{"(", - RowBox[{ - RowBox[{ - RowBox[{ - SubscriptBox["\[CapitalEpsilon]", "HF"], "[", "R", "]"}], - "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], "-", - RowBox[{ - RowBox[{ - SubscriptBox["\[CapitalEpsilon]", "HF"], "[", "R", "]"}], - "\[LeftDoubleBracket]", "2", "\[RightDoubleBracket]"}]}], ")"}], - "2"], "+", - RowBox[{"4", - SuperscriptBox[ - RowBox[{ - SubscriptBox["K", "12"], "[", "R", "]"}], "2"]}]}]]}]}]}], - "}"}]}], ";"}]], "Input", + RowBox[{"2", + RowBox[{ + SubscriptBox["J", "12"], "[", "R", "]"}]}], "-", + RowBox[{ + SubscriptBox["K", "12"], "[", "R", "]"}]}], ")"}]}], "+", + RowBox[{"w", + RowBox[{"(", + RowBox[{ + RowBox[{"2", + RowBox[{ + SubscriptBox["J", "22"], "[", "R", "]"}]}], "-", + RowBox[{ + SubscriptBox["K", "22"], "[", "R", "]"}]}], ")"}]}]}]}], "}"}]}], + ";"}]], "Input", InitializationCell->True, - CellChangeTimes->{{3.790605458040031*^9, 3.790605458543043*^9}, { - 3.790605491525733*^9, 3.79060554005932*^9}, {3.790614460047532*^9, - 3.790614470197546*^9}}, - CellLabel->"In[27]:=",ExpressionUUID->"ab8163b6-fc27-4314-ae7e-a7f0f7d95629"], + CellChangeTimes->{{3.790660423886979*^9, 3.790660530808771*^9}, { + 3.790661586617362*^9, 3.7906615882014847`*^9}}, + CellLabel-> + "In[114]:=",ExpressionUUID->"36f3a3dd-e13f-4bbf-a50c-b11a31e789f8"], Cell[TextData[StyleBox["LDA energies", FontColor->RGBColor[1, 0, 0]]], "Text", @@ -4685,82 +4806,87 @@ Cell[BoxData[ RowBox[{ RowBox[{ SubscriptBox["\[CapitalEpsilon]", "LDA"], "[", "R_", "]"}], ":=", - RowBox[{"{", - RowBox[{ + RowBox[{ + RowBox[{"{", "\[IndentingNewLine]", RowBox[{ - RowBox[{"2", - RowBox[{ - SubscriptBox["h", "1"], "[", "R", "]"}]}], "+", - RowBox[{"2", - RowBox[{ - SubscriptBox["J", "11"], "[", "R", "]"}]}], "+", - RowBox[{"NIntegrate", "[", - RowBox[{ + RowBox[{ + RowBox[{"2", RowBox[{ - RowBox[{"\[Epsilon]xcLDA", "[", + SubscriptBox["h", "1"], "[", "R", "]"}]}], "+", + RowBox[{"2", + RowBox[{ + SubscriptBox["J", "11"], "[", "R", "]"}]}], "+", + RowBox[{"NIntegrate", "[", + RowBox[{ + RowBox[{ + RowBox[{"\[Epsilon]xcLDA", "[", + RowBox[{ + SubscriptBox["n", "1"], "[", + RowBox[{ + RowBox[{"{", + RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R"}], "]"}], + "]"}], RowBox[{ SubscriptBox["n", "1"], "[", RowBox[{ RowBox[{"{", - RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R"}], "]"}], - "]"}], - RowBox[{ - SubscriptBox["n", "1"], "[", - RowBox[{ - RowBox[{"{", - RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R"}], "]"}]}], - ",", - RowBox[{"{", - RowBox[{"x", ",", - RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", - RowBox[{"{", - RowBox[{"y", ",", - RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", - RowBox[{"{", - RowBox[{"z", ",", - RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}]}], - "]"}]}], ",", - RowBox[{ - RowBox[{"2", - RowBox[{ - SubscriptBox["h", "2"], "[", "R", "]"}]}], "+", - RowBox[{"2", - RowBox[{ - SubscriptBox["J", "22"], "[", "R", "]"}]}], "+", - RowBox[{"NIntegrate", "[", - RowBox[{ + RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R"}], "]"}]}], + ",", + RowBox[{"{", + RowBox[{"x", ",", + RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", + RowBox[{"{", + RowBox[{"y", ",", + RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", + RowBox[{"{", + RowBox[{"z", ",", + RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}]}], + "]"}]}], ",", "\[IndentingNewLine]", + RowBox[{ + RowBox[{"2", RowBox[{ - RowBox[{"\[Epsilon]xcLDA", "[", + SubscriptBox["h", "2"], "[", "R", "]"}]}], "+", + RowBox[{"2", + RowBox[{ + SubscriptBox["J", "22"], "[", "R", "]"}]}], "+", + RowBox[{"NIntegrate", "[", + RowBox[{ + RowBox[{ + RowBox[{"\[Epsilon]xcLDA", "[", + RowBox[{ + SubscriptBox["n", "2"], "[", + RowBox[{ + RowBox[{"{", + RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R"}], "]"}], + "]"}], RowBox[{ SubscriptBox["n", "2"], "[", RowBox[{ RowBox[{"{", - RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R"}], "]"}], - "]"}], - RowBox[{ - SubscriptBox["n", "2"], "[", - RowBox[{ - RowBox[{"{", - RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R"}], "]"}]}], - ",", - RowBox[{"{", - RowBox[{"x", ",", - RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", - RowBox[{"{", - RowBox[{"y", ",", - RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", - RowBox[{"{", - RowBox[{"z", ",", - RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}]}], - "]"}]}]}], "}"}]}], ";"}]], "Input", + RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R"}], "]"}]}], + ",", + RowBox[{"{", + RowBox[{"x", ",", + RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", + RowBox[{"{", + RowBox[{"y", ",", + RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", + RowBox[{"{", + RowBox[{"z", ",", + RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}]}], + "]"}]}]}], "\[IndentingNewLine]", "}"}], "+", + RowBox[{ + SubscriptBox["\[CapitalEpsilon]", "NN"], "[", "R", "]"}]}]}], + ";"}]], "Input", InitializationCell->True, CellChangeTimes->{{3.790605458040031*^9, 3.790605458543043*^9}, { - 3.790605491525733*^9, 3.79060554005932*^9}, {3.7906056671315413`*^9, - 3.790605692291338*^9}, {3.790610838410942*^9, 3.790610918068091*^9}, { - 3.7906110542128553`*^9, 3.790611062835252*^9}, {3.790618152013009*^9, - 3.790618153069519*^9}, {3.790618673469589*^9, - 3.790618675001255*^9}},ExpressionUUID->"4419a1ba-532a-442b-95d2-\ -330c205bfd7e"], + 3.790605491525733*^9, 3.79060554005932*^9}, {3.7906056671315413`*^9, + 3.790605692291338*^9}, {3.790610838410942*^9, 3.790610918068091*^9}, { + 3.7906110542128553`*^9, 3.790611062835252*^9}, {3.790618152013009*^9, + 3.790618153069519*^9}, {3.790618673469589*^9, 3.790618675001255*^9}, + 3.790659180965103*^9, {3.790660627881838*^9, 3.7906606344245253`*^9}}, + CellLabel-> + "In[115]:=",ExpressionUUID->"4419a1ba-532a-442b-95d2-330c205bfd7e"], Cell[BoxData[ RowBox[{ @@ -4814,14 +4940,165 @@ Cell[BoxData[ RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", RowBox[{"{", RowBox[{"z", ",", - RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}]}], - "]"}]}]}], ";"}]], "Input", + RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}]}], "]"}], + "+", + RowBox[{ + SubscriptBox["\[CapitalEpsilon]", "NN"], "[", "R", "]"}]}]}], + ";"}]], "Input", InitializationCell->True, CellChangeTimes->{{3.7906191044624767`*^9, 3.7906191773964367`*^9}, { - 3.790619286782889*^9, 3.790619299369211*^9}, {3.790619337303144*^9, - 3.790619344809649*^9}, {3.790619402147251*^9, 3.790619402226974*^9}}, + 3.790619286782889*^9, 3.790619299369211*^9}, {3.790619337303144*^9, + 3.790619344809649*^9}, {3.790619402147251*^9, 3.790619402226974*^9}, + 3.790659183868813*^9}, CellLabel-> - "In[193]:=",ExpressionUUID->"c0b69368-df8e-46ab-a64b-3ea147722dd9"], + "In[116]:=",ExpressionUUID->"c0b69368-df8e-46ab-a64b-3ea147722dd9"], + +Cell[BoxData[ + RowBox[{ + RowBox[{ + RowBox[{ + SubscriptBox["\[CapitalDelta]\[CapitalEpsilon]w", "LDA"], "[", + RowBox[{"R_", ",", "w_"}], "]"}], ":=", + RowBox[{ + RowBox[{ + RowBox[{"(", + RowBox[{"1", "-", "w"}], ")"}], + RowBox[{ + RowBox[{ + SubscriptBox["\[CapitalEpsilon]", "LDA"], "[", "R", "]"}], + "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}]}], "+", + RowBox[{"w", + RowBox[{ + RowBox[{ + SubscriptBox["\[CapitalEpsilon]", "LDA"], "[", "R", "]"}], + "\[LeftDoubleBracket]", "2", "\[RightDoubleBracket]"}]}], "-", + RowBox[{ + SubscriptBox["\[CapitalEpsilon]w", "LDA"], "[", + RowBox[{"R", ",", "w"}], "]"}]}]}], ";"}]], "Input", + InitializationCell->True, + CellChangeTimes->{{3.790660001969771*^9, 3.790660091003895*^9}, + 3.790660542656966*^9}, + CellLabel-> + "In[117]:=",ExpressionUUID->"0dd5826b-1b72-4f2c-834b-e79b419ac43e"], + +Cell[BoxData[ + RowBox[{ + RowBox[{ + RowBox[{ + SubscriptBox["\[Epsilon]w", "LDA"], "[", + RowBox[{"R_", ",", "w_"}], "]"}], ":=", + RowBox[{"{", "\[IndentingNewLine]", + RowBox[{ + RowBox[{ + RowBox[{ + SubscriptBox["h", "1"], "[", "R", "]"}], "+", + RowBox[{"2", + RowBox[{"(", + RowBox[{"1", "-", "w"}], ")"}], + RowBox[{ + SubscriptBox["J", "11"], "[", "R", "]"}]}], "+", + RowBox[{"2", "w", " ", + RowBox[{ + SubscriptBox["J", "12"], "[", "R", "]"}]}], "+", + RowBox[{ + FractionBox["1", "2"], + RowBox[{"NIntegrate", "[", + RowBox[{ + RowBox[{ + RowBox[{"(", + RowBox[{ + RowBox[{ + RowBox[{"d\[Epsilon]xcLDAdn", "[", + RowBox[{"nw", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R", ",", + "w"}], "]"}], "]"}], + RowBox[{"nw", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R", ",", + "w"}], "]"}]}], "+", + RowBox[{"\[Epsilon]xcLDA", "[", + RowBox[{"nw", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R", ",", + "w"}], "]"}], "]"}]}], ")"}], + RowBox[{ + SubscriptBox["n", "1"], "[", + RowBox[{ + RowBox[{"{", + RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R"}], "]"}]}], + ",", + RowBox[{"{", + RowBox[{"x", ",", + RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", + RowBox[{"{", + RowBox[{"y", ",", + RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", + RowBox[{"{", + RowBox[{"z", ",", + RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}]}], + "]"}]}]}], ",", "\[IndentingNewLine]", + RowBox[{ + RowBox[{ + SubscriptBox["h", "2"], "[", "R", "]"}], "+", + RowBox[{"2", + RowBox[{"(", + RowBox[{"1", "-", "w"}], ")"}], + RowBox[{ + SubscriptBox["J", "12"], "[", "R", "]"}]}], "+", + RowBox[{"2", "w", " ", + RowBox[{ + SubscriptBox["J", "22"], "[", "R", "]"}]}], "+", + RowBox[{ + FractionBox["1", "2"], + RowBox[{"NIntegrate", "[", + RowBox[{ + RowBox[{ + RowBox[{"(", + RowBox[{ + RowBox[{ + RowBox[{"d\[Epsilon]xcLDAdn", "[", + RowBox[{"nw", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R", ",", + "w"}], "]"}], "]"}], + RowBox[{"nw", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R", ",", + "w"}], "]"}]}], "+", + RowBox[{"\[Epsilon]xcLDA", "[", + RowBox[{"nw", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R", ",", + "w"}], "]"}], "]"}]}], ")"}], + RowBox[{ + SubscriptBox["n", "2"], "[", + RowBox[{ + RowBox[{"{", + RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R"}], "]"}]}], + ",", + RowBox[{"{", + RowBox[{"x", ",", + RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", + RowBox[{"{", + RowBox[{"y", ",", + RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", + RowBox[{"{", + RowBox[{"z", ",", + RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}]}], + "]"}]}]}]}], "\[IndentingNewLine]", "}"}]}], ";"}]], "Input", + InitializationCell->True, + CellChangeTimes->{{3.790660423886979*^9, 3.790660614840477*^9}, { + 3.790660673368799*^9, 3.790660734031897*^9}, {3.790661577778819*^9, + 3.7906615810821047`*^9}, {3.79066238202234*^9, + 3.7906624156374083`*^9}},ExpressionUUID->"61cc9c57-6cf1-464b-bd07-\ +45e24fd7c9c3"], Cell[TextData[StyleBox["eLDA energies", FontColor->RGBColor[1, 0, 0]]], "Text", @@ -4836,85 +5113,90 @@ Cell[BoxData[ RowBox[{ RowBox[{ SubscriptBox["\[CapitalEpsilon]", "eLDA"], "[", "R_", "]"}], ":=", - RowBox[{"{", + RowBox[{ + RowBox[{"{", "\[IndentingNewLine]", + RowBox[{ + RowBox[{ + RowBox[{"2", + RowBox[{ + SubscriptBox["h", "1"], "[", "R", "]"}]}], "+", + RowBox[{"2", + RowBox[{ + SubscriptBox["J", "11"], "[", "R", "]"}]}], "+", + RowBox[{"NIntegrate", "[", + RowBox[{ + RowBox[{ + RowBox[{"\[Epsilon]xceLDA", "[", + RowBox[{ + RowBox[{ + SubscriptBox["n", "1"], "[", + RowBox[{ + RowBox[{"{", + RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R"}], "]"}], + ",", "0"}], "]"}], + RowBox[{ + SubscriptBox["n", "1"], "[", + RowBox[{ + RowBox[{"{", + RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R"}], "]"}]}], + ",", + RowBox[{"{", + RowBox[{"x", ",", + RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", + RowBox[{"{", + RowBox[{"y", ",", + RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", + RowBox[{"{", + RowBox[{"z", ",", + RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}]}], + "]"}]}], ",", "\[IndentingNewLine]", + RowBox[{ + RowBox[{"2", + RowBox[{ + SubscriptBox["h", "2"], "[", "R", "]"}]}], "+", + RowBox[{"2", + RowBox[{ + SubscriptBox["J", "22"], "[", "R", "]"}]}], "+", + RowBox[{"NIntegrate", "[", + RowBox[{ + RowBox[{ + RowBox[{"\[Epsilon]xceLDA", "[", + RowBox[{ + RowBox[{ + SubscriptBox["n", "2"], "[", + RowBox[{ + RowBox[{"{", + RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R"}], "]"}], + ",", "1"}], "]"}], + RowBox[{ + SubscriptBox["n", "2"], "[", + RowBox[{ + RowBox[{"{", + RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R"}], "]"}]}], + ",", + RowBox[{"{", + RowBox[{"x", ",", + RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", + RowBox[{"{", + RowBox[{"y", ",", + RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", + RowBox[{"{", + RowBox[{"z", ",", + RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}]}], + "]"}]}]}], "\[IndentingNewLine]", "}"}], "+", RowBox[{ - RowBox[{ - RowBox[{"2", - RowBox[{ - SubscriptBox["h", "1"], "[", "R", "]"}]}], "+", - RowBox[{"2", - RowBox[{ - SubscriptBox["J", "11"], "[", "R", "]"}]}], "+", - RowBox[{"NIntegrate", "[", - RowBox[{ - RowBox[{ - RowBox[{"\[Epsilon]xceLDA", "[", - RowBox[{ - RowBox[{ - SubscriptBox["n", "1"], "[", - RowBox[{ - RowBox[{"{", - RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R"}], "]"}], - ",", "0"}], "]"}], - RowBox[{ - SubscriptBox["n", "1"], "[", - RowBox[{ - RowBox[{"{", - RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R"}], "]"}]}], - ",", - RowBox[{"{", - RowBox[{"x", ",", - RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", - RowBox[{"{", - RowBox[{"y", ",", - RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", - RowBox[{"{", - RowBox[{"z", ",", - RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}]}], - "]"}]}], ",", - RowBox[{ - RowBox[{"2", - RowBox[{ - SubscriptBox["h", "2"], "[", "R", "]"}]}], "+", - RowBox[{"2", - RowBox[{ - SubscriptBox["J", "22"], "[", "R", "]"}]}], "+", - RowBox[{"NIntegrate", "[", - RowBox[{ - RowBox[{ - RowBox[{"\[Epsilon]xceLDA", "[", - RowBox[{ - RowBox[{ - SubscriptBox["n", "2"], "[", - RowBox[{ - RowBox[{"{", - RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R"}], "]"}], - ",", "1"}], "]"}], - RowBox[{ - SubscriptBox["n", "2"], "[", - RowBox[{ - RowBox[{"{", - RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R"}], "]"}]}], - ",", - RowBox[{"{", - RowBox[{"x", ",", - RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", - RowBox[{"{", - RowBox[{"y", ",", - RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", - RowBox[{"{", - RowBox[{"z", ",", - RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}]}], - "]"}]}]}], "}"}]}], ";"}]], "Input", + SubscriptBox["\[CapitalEpsilon]", "NN"], "[", "R", "]"}]}]}], + ";"}]], "Input", InitializationCell->True, CellChangeTimes->{{3.790605458040031*^9, 3.790605458543043*^9}, { - 3.790605491525733*^9, 3.79060554005932*^9}, {3.7906056671315413`*^9, - 3.790605692291338*^9}, {3.790610838410942*^9, 3.790610918068091*^9}, { - 3.7906110542128553`*^9, 3.790611062835252*^9}, {3.790615225694908*^9, - 3.7906152329132547`*^9}, {3.790618147868698*^9, 3.790618149719122*^9}, { - 3.790618670193342*^9, 3.79061867932382*^9}}, + 3.790605491525733*^9, 3.79060554005932*^9}, {3.7906056671315413`*^9, + 3.790605692291338*^9}, {3.790610838410942*^9, 3.790610918068091*^9}, { + 3.7906110542128553`*^9, 3.790611062835252*^9}, {3.790615225694908*^9, + 3.7906152329132547`*^9}, {3.790618147868698*^9, 3.790618149719122*^9}, { + 3.790618670193342*^9, 3.79061867932382*^9}, 3.790659187972827*^9, { + 3.790660643667115*^9, 3.790660649050318*^9}}, CellLabel-> - "In[174]:=",ExpressionUUID->"d4265afd-5fb1-4adc-ad27-da2413b01953"], + "In[119]:=",ExpressionUUID->"d4265afd-5fb1-4adc-ad27-da2413b01953"], Cell[BoxData[ RowBox[{ @@ -4972,15 +5254,363 @@ Cell[BoxData[ "]"}]}]}]], "Input", InitializationCell->True, CellChangeTimes->{{3.7906191044624767`*^9, 3.7906191943435097`*^9}, { - 3.790619315024371*^9, 3.79061931729541*^9}, {3.790619349519698*^9, - 3.790619356926655*^9}, {3.790619407198872*^9, 3.790619409338375*^9}, { - 3.7906194502646437`*^9, 3.790619451317535*^9}, {3.7906195478586807`*^9, - 3.790619547970214*^9}}, + 3.790619315024371*^9, 3.79061931729541*^9}, {3.790619349519698*^9, + 3.790619356926655*^9}, {3.790619407198872*^9, 3.790619409338375*^9}, { + 3.7906194502646437`*^9, 3.790619451317535*^9}, {3.7906195478586807`*^9, + 3.790619547970214*^9}, 3.790659191412819*^9}, CellLabel-> - "In[203]:=",ExpressionUUID->"9cd135b3-7599-482f-9e19-dfd18b561c56"] + "In[120]:=",ExpressionUUID->"9cd135b3-7599-482f-9e19-dfd18b561c56"], + +Cell[BoxData[ + RowBox[{ + RowBox[{ + RowBox[{ + SubscriptBox["\[CapitalDelta]\[CapitalEpsilon]w", "eLDA"], "[", + RowBox[{"R_", ",", "w_"}], "]"}], ":=", + RowBox[{ + RowBox[{ + RowBox[{"(", + RowBox[{"1", "-", "w"}], ")"}], + RowBox[{ + RowBox[{ + SubscriptBox["\[CapitalEpsilon]", "eLDA"], "[", "R", "]"}], + "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}]}], "+", + RowBox[{"w", + RowBox[{ + RowBox[{ + SubscriptBox["\[CapitalEpsilon]", "eLDA"], "[", "R", "]"}], + "\[LeftDoubleBracket]", "2", "\[RightDoubleBracket]"}]}], "-", + RowBox[{ + SubscriptBox["\[CapitalEpsilon]w", "eLDA"], "[", + RowBox[{"R", ",", "w"}], "]"}]}]}], ";"}]], "Input", + InitializationCell->True, + CellChangeTimes->{{3.790660001969771*^9, 3.790660105027638*^9}, + 3.790660544962287*^9}, + CellLabel-> + "In[121]:=",ExpressionUUID->"99ae0dcb-f21c-46d4-84aa-7f172765d9d8"], + +Cell[BoxData[ + RowBox[{ + RowBox[{ + RowBox[{ + SubscriptBox["\[Epsilon]w", "eLDA"], "[", + RowBox[{"R_", ",", "w_"}], "]"}], ":=", + RowBox[{"{", "\[IndentingNewLine]", + RowBox[{ + RowBox[{ + RowBox[{ + SubscriptBox["h", "1"], "[", "R", "]"}], "+", + RowBox[{"2", + RowBox[{"(", + RowBox[{"1", "-", "w"}], ")"}], + RowBox[{ + SubscriptBox["J", "11"], "[", "R", "]"}]}], "+", + RowBox[{"2", "w", " ", + RowBox[{ + SubscriptBox["J", "12"], "[", "R", "]"}]}], "+", + RowBox[{ + FractionBox["1", "2"], + RowBox[{"NIntegrate", "[", + RowBox[{ + RowBox[{ + RowBox[{"(", + RowBox[{ + RowBox[{ + RowBox[{"d\[Epsilon]xceLDAdn", "[", + RowBox[{ + RowBox[{"nw", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R", ",", + "w"}], "]"}], ",", "w"}], "]"}], + RowBox[{"nw", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R", ",", + "w"}], "]"}]}], "+", + RowBox[{"\[Epsilon]xceLDA", "[", + RowBox[{ + RowBox[{"nw", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R", ",", + "w"}], "]"}], ",", "w"}], "]"}]}], ")"}], + RowBox[{ + SubscriptBox["n", "1"], "[", + RowBox[{ + RowBox[{"{", + RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R"}], "]"}]}], + ",", + RowBox[{"{", + RowBox[{"x", ",", + RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", + RowBox[{"{", + RowBox[{"y", ",", + RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", + RowBox[{"{", + RowBox[{"z", ",", + RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}]}], + "]"}]}]}], ",", "\[IndentingNewLine]", + RowBox[{ + RowBox[{ + SubscriptBox["h", "2"], "[", "R", "]"}], "+", + RowBox[{"2", + RowBox[{"(", + RowBox[{"1", "-", "w"}], ")"}], + RowBox[{ + SubscriptBox["J", "12"], "[", "R", "]"}]}], "+", + RowBox[{"2", "w", " ", + RowBox[{ + SubscriptBox["J", "22"], "[", "R", "]"}]}], "+", + RowBox[{ + FractionBox["1", "2"], + RowBox[{"NIntegrate", "[", + RowBox[{ + RowBox[{ + RowBox[{"(", + RowBox[{ + RowBox[{ + RowBox[{"d\[Epsilon]xceLDAdn", "[", + RowBox[{ + RowBox[{"nw", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R", ",", + "w"}], "]"}], ",", "w"}], "]"}], + RowBox[{"nw", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R", ",", + "w"}], "]"}]}], "+", + RowBox[{"\[Epsilon]xceLDA", "[", + RowBox[{ + RowBox[{"nw", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R", ",", + "w"}], "]"}], ",", "w"}], "]"}]}], ")"}], + RowBox[{ + SubscriptBox["n", "2"], "[", + RowBox[{ + RowBox[{"{", + RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R"}], "]"}]}], + ",", + RowBox[{"{", + RowBox[{"x", ",", + RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", + RowBox[{"{", + RowBox[{"y", ",", + RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", + RowBox[{"{", + RowBox[{"z", ",", + RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}]}], + "]"}]}]}]}], "\[IndentingNewLine]", "}"}]}], ";"}]], "Input", + InitializationCell->True, + CellChangeTimes->{{3.790660423886979*^9, 3.790660614840477*^9}, { + 3.790660673368799*^9, 3.790660753527668*^9}, {3.7906615534098883`*^9, + 3.790661574889584*^9}}, + CellLabel-> + "In[139]:=",ExpressionUUID->"57a35adb-5f0e-4384-a38f-33a0653fb6b6"], + +Cell[BoxData[ + RowBox[{"With", "[", + RowBox[{ + RowBox[{"{", + RowBox[{ + RowBox[{"Rmin", "=", "0.1"}], ",", + RowBox[{"Rmax", "=", "5"}], ",", + RowBox[{"\[CapitalDelta]R", "=", "0.1"}]}], "}"}], ",", + "\[IndentingNewLine]", + RowBox[{ + RowBox[{"ScanEHF", "=", + RowBox[{"Table", "[", + RowBox[{ + RowBox[{"Flatten", "[", + RowBox[{"{", + RowBox[{"R", ",", + RowBox[{ + SubscriptBox["\[CapitalEpsilon]", "HF"], "[", "R", "]"}]}], "}"}], + "]"}], ",", + RowBox[{"{", + RowBox[{"R", ",", "Rmin", ",", "Rmax", ",", "\[CapitalDelta]R"}], + "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", + RowBox[{"ScanEwHF", "=", + RowBox[{"Table", "[", + RowBox[{ + RowBox[{"Flatten", "[", + RowBox[{"{", + RowBox[{"R", ",", + RowBox[{ + SubscriptBox["\[CapitalEpsilon]w", "HF"], "[", + RowBox[{"R", ",", + FractionBox["1", "2"]}], "]"}]}], "}"}], "]"}], ",", + RowBox[{"{", + RowBox[{"R", ",", "Rmin", ",", "Rmax", ",", "\[CapitalDelta]R"}], + "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", + RowBox[{"ScanECID", "=", + RowBox[{"Table", "[", + RowBox[{ + RowBox[{"Flatten", "[", + RowBox[{"{", + RowBox[{"R", ",", + RowBox[{ + SubscriptBox["\[CapitalEpsilon]", "CID"], "[", "R", "]"}]}], "}"}], + "]"}], ",", + RowBox[{"{", + RowBox[{"R", ",", "Rmin", ",", "Rmax", ",", "\[CapitalDelta]R"}], + "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", + RowBox[{"ScanELDA", "=", + RowBox[{"Table", "[", + RowBox[{ + RowBox[{"Flatten", "[", + RowBox[{"{", + RowBox[{"R", ",", + RowBox[{ + SubscriptBox["\[CapitalEpsilon]", "LDA"], "[", "R", "]"}]}], "}"}], + "]"}], ",", + RowBox[{"{", + RowBox[{"R", ",", "Rmin", ",", "Rmax", ",", "\[CapitalDelta]R"}], + "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", + RowBox[{"ScanEwLDA", "=", + RowBox[{"Table", "[", + RowBox[{ + RowBox[{"Flatten", "[", + RowBox[{"{", + RowBox[{"R", ",", + RowBox[{ + SubscriptBox["\[CapitalEpsilon]w", "LDA"], "[", + RowBox[{"R", ",", + FractionBox["1", "2"]}], "]"}]}], "}"}], "]"}], ",", + RowBox[{"{", + RowBox[{"R", ",", "Rmin", ",", "Rmax", ",", "\[CapitalDelta]R"}], + "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", + RowBox[{"ScanEeLDA", "=", + RowBox[{"Table", "[", + RowBox[{ + RowBox[{"Flatten", "[", + RowBox[{"{", + RowBox[{"R", ",", + RowBox[{ + SubscriptBox["\[CapitalEpsilon]", "eLDA"], "[", "R", "]"}]}], + "}"}], "]"}], ",", + RowBox[{"{", + RowBox[{"R", ",", "Rmin", ",", "Rmax", ",", "\[CapitalDelta]R"}], + "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", + RowBox[{"ScanEweLDA", "=", + RowBox[{"Table", "[", + RowBox[{ + RowBox[{"Flatten", "[", + RowBox[{"{", + RowBox[{"R", ",", + RowBox[{ + SubscriptBox["\[CapitalEpsilon]w", "eLDA"], "[", + RowBox[{"R", ",", + FractionBox["1", "2"]}], "]"}]}], "}"}], "]"}], ",", + RowBox[{"{", + RowBox[{"R", ",", "Rmin", ",", "Rmax", ",", "\[CapitalDelta]R"}], + "}"}]}], "]"}]}], ";"}]}], "\[IndentingNewLine]", "]"}]], "Input", + InitializationCell->True, + CellChangeTimes->{{3.790658215232338*^9, 3.790658236824192*^9}, { + 3.7906583897274837`*^9, 3.790658524654791*^9}, {3.790659704040538*^9, + 3.790659705870874*^9}, {3.790659735886373*^9, 3.790659813232091*^9}}, + CellLabel-> + "In[123]:=",ExpressionUUID->"93054b36-78ec-4bb1-83ea-8fbb31835fb9"], + +Cell[BoxData[ + RowBox[{ + RowBox[{"d\[CapitalEpsilon]weLDAdw", "[", + RowBox[{"R_", ",", "w_"}], "]"}], ":=", + RowBox[{"NIntegrate", "[", + RowBox[{ + RowBox[{ + RowBox[{"d\[Epsilon]xceLDAdw", "[", + RowBox[{"nw", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R", ",", "w"}], + "]"}], "]"}], + RowBox[{"nw", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", "R", ",", "w"}], + "]"}]}], ",", + RowBox[{"{", + RowBox[{"x", ",", + RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", + RowBox[{"{", + RowBox[{"y", ",", + RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", + RowBox[{"{", + RowBox[{"z", ",", + RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}]}], + "]"}]}]], "Input", + CellChangeTimes->{{3.7906608146381903`*^9, 3.790660820086027*^9}, { + 3.790661275357267*^9, 3.790661293506763*^9}, {3.790661435266225*^9, + 3.790661483849555*^9}, {3.790661526538616*^9, 3.790661526618237*^9}, { + 3.790661636466195*^9, 3.790661654523027*^9}, 3.790661792604311*^9}, + CellLabel-> + "In[134]:=",ExpressionUUID->"9eac41d8-7253-4bcb-9a2a-4a332b730ddb"], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{ + RowBox[{"2", + RowBox[{"(", + RowBox[{ + RowBox[{ + RowBox[{ + SubscriptBox["\[Epsilon]w", "LDA"], "[", + RowBox[{"R", ",", "w"}], "]"}], "\[LeftDoubleBracket]", "2", + "\[RightDoubleBracket]"}], "-", + RowBox[{ + RowBox[{ + SubscriptBox["\[Epsilon]w", "LDA"], "[", + RowBox[{"R", ",", "w"}], "]"}], "\[LeftDoubleBracket]", "1", + "\[RightDoubleBracket]"}]}], ")"}]}], "/.", + RowBox[{"{", + RowBox[{ + RowBox[{"R", "\[Rule]", "1.4"}], ",", + RowBox[{"w", "\[Rule]", "0"}]}], "}"}]}]], "Input", + CellChangeTimes->{{3.790662000166115*^9, 3.7906620185750113`*^9}, { + 3.790662138430471*^9, 3.790662147334878*^9}, 3.790662401812807*^9}, + CellLabel-> + "In[150]:=",ExpressionUUID->"cc5802c8-5164-4d35-a4aa-ac0f06412d48"], + +Cell[BoxData["1.4989461809762727`"], "Output", + CellChangeTimes->{{3.790662007156047*^9, 3.7906620189650087`*^9}, { + 3.790662139159606*^9, 3.7906621611887913`*^9}, {3.7906623969111633`*^9, + 3.7906624071288967`*^9}}, + CellLabel-> + "Out[150]=",ExpressionUUID->"6c0f8654-3af4-43e9-b5c1-7a47f72b02d4"] +}, Open ]], + +Cell[BoxData[ + RowBox[{ + RowBox[{ + RowBox[{"2", + RowBox[{"(", + RowBox[{ + RowBox[{ + RowBox[{ + SubscriptBox["\[Epsilon]w", "eLDA"], "[", + RowBox[{"R", ",", "w"}], "]"}], "\[LeftDoubleBracket]", "2", + "\[RightDoubleBracket]"}], "-", + RowBox[{ + RowBox[{ + SubscriptBox["\[Epsilon]w", "eLDA"], "[", + RowBox[{"R", ",", "w"}], "]"}], "\[LeftDoubleBracket]", "1", + "\[RightDoubleBracket]"}]}], ")"}]}], "+", + RowBox[{"d\[CapitalEpsilon]weLDAdw", "[", + RowBox[{"R", ",", "w"}], "]"}]}], "\[IndentingNewLine]"}]], "Input", + CellChangeTimes->{{3.790661475825005*^9, 3.790661489241352*^9}, { + 3.7906618427444572`*^9, 3.790661844966864*^9}, {3.790661881697806*^9, + 3.7906618992816668`*^9}, {3.7906619466147127`*^9, 3.790661955270863*^9}}, + CellLabel-> + "In[140]:=",ExpressionUUID->"8cbe4386-fbc6-4914-8e8c-e90d47f39535"] }, Open ]] }, -WindowSize->{1767, 1268}, +WindowSize->{1815, 1395}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, FrontEndVersion->"12.0 for Mac OS X x86 (64-bit) (April 8, 2019)", StyleDefinitions->"Default.nb" @@ -5010,9 +5640,9 @@ Cell[2842, 83, 343, 8, 46, "Input",ExpressionUUID->"058853fd-911e-44b3-8eac-1c3f InitializationCell->True], Cell[3188, 93, 279, 6, 46, "Input",ExpressionUUID->"5dbbd40c-d911-445f-b476-85cd91f8d014", InitializationCell->True] -}, Open ]], +}, Closed]], Cell[CellGroupData[{ -Cell[3504, 104, 157, 3, 98, "Title",ExpressionUUID->"98d632cb-7177-481c-9b9f-d6f48d10e284"], +Cell[3504, 104, 157, 3, 72, "Title",ExpressionUUID->"98d632cb-7177-481c-9b9f-d6f48d10e284"], Cell[CellGroupData[{ Cell[3686, 111, 207, 4, 67, "Section",ExpressionUUID->"9458067c-0b7a-4f70-9a58-1f771c4b3842"], Cell[CellGroupData[{ @@ -5023,10 +5653,10 @@ Cell[5774, 169, 1412, 35, 65, "Input",ExpressionUUID->"f4215905-d329-4e5c-925a-f InitializationCell->True], Cell[7189, 206, 2169, 59, 65, "Input",ExpressionUUID->"517a10dc-802a-4ead-9f44-0905089379ba", InitializationCell->True], -Cell[9361, 267, 3587, 99, 137, "Input",ExpressionUUID->"440fb372-2001-4d90-ade1-e5498676702f", +Cell[9361, 267, 3587, 99, 72, "Input",ExpressionUUID->"440fb372-2001-4d90-ade1-e5498676702f", InitializationCell->True] }, Open ]] -}, Open ]] +}, Closed]] }, Open ]], Cell[CellGroupData[{ Cell[13009, 373, 151, 3, 98, "Title",ExpressionUUID->"cdfcc3a6-7299-4a1b-a9d0-f29078225f1e"], @@ -5201,68 +5831,89 @@ Cell[100310, 2828, 448, 13, 55, "Output",ExpressionUUID->"27597e7b-3c41-49b8-990 Cell[CellGroupData[{ Cell[100819, 2848, 211, 4, 98, "Title",ExpressionUUID->"d8a7ce9c-a1d9-4ac0-b8bf-9ca66a94ac58"], Cell[101033, 2854, 197, 3, 35, "Text",ExpressionUUID->"d408aeeb-1654-4525-903e-88994da645ee"], -Cell[101233, 2859, 2135, 65, 171, "Input",ExpressionUUID->"f51f78c8-fe8e-48f3-a562-1209bfc50f71", +Cell[101233, 2859, 2158, 65, 171, "Input",ExpressionUUID->"f51f78c8-fe8e-48f3-a562-1209bfc50f71", InitializationCell->True], -Cell[103371, 2926, 198, 3, 35, "Text",ExpressionUUID->"683219a6-4bcc-4fb8-93b5-7c590bc511b5"], -Cell[103572, 2931, 5787, 168, 394, "Input",ExpressionUUID->"c5e15fc2-7ce6-4ed3-bc7f-bcdc9085134d", +Cell[103394, 2926, 1370, 36, 110, "Input",ExpressionUUID->"dd062623-5a91-460e-945c-d5f15e24abfb", InitializationCell->True], -Cell[109362, 3101, 2355, 70, 167, "Input",ExpressionUUID->"8db93606-ec73-4fe8-9228-51f79c4ef957", +Cell[104767, 2964, 198, 3, 35, "Text",ExpressionUUID->"683219a6-4bcc-4fb8-93b5-7c590bc511b5"], +Cell[104968, 2969, 5783, 167, 394, "Input",ExpressionUUID->"c5e15fc2-7ce6-4ed3-bc7f-bcdc9085134d", InitializationCell->True], -Cell[111720, 3173, 256, 4, 35, "Text",ExpressionUUID->"264d5181-32f7-492d-a01f-33c1ad90708e"], -Cell[111979, 3179, 1093, 29, 110, "Input",ExpressionUUID->"f6517707-c164-4ea9-b57c-8d4d27eafad3", +Cell[110754, 3138, 2351, 69, 167, "Input",ExpressionUUID->"8db93606-ec73-4fe8-9228-51f79c4ef957", InitializationCell->True], -Cell[113075, 3210, 214, 4, 35, "Text",ExpressionUUID->"d38e2fe6-260d-4638-b639-6984dd9d14e4"], -Cell[113292, 3216, 679, 19, 63, "Input",ExpressionUUID->"36b985d7-30a2-430d-8dc3-30d6b4a45727", +Cell[113108, 3209, 1368, 36, 110, "Input",ExpressionUUID->"3c5be8aa-305f-4c0e-a51a-4bc92cdcf171", InitializationCell->True], -Cell[113974, 3237, 210, 4, 35, "Text",ExpressionUUID->"e115f0cc-d9e8-4997-a5a4-88b0e9ed97e0"], -Cell[114187, 3243, 3586, 92, 187, "Input",ExpressionUUID->"f0a2825f-92fc-42d2-aeee-57f53535b1d5", +Cell[114479, 3247, 256, 4, 35, "Text",ExpressionUUID->"264d5181-32f7-492d-a01f-33c1ad90708e"], +Cell[114738, 3253, 1089, 28, 110, "Input",ExpressionUUID->"f6517707-c164-4ea9-b57c-8d4d27eafad3", InitializationCell->True], -Cell[117776, 3337, 276, 5, 35, "Text",ExpressionUUID->"5353bb94-3a0f-4a57-a4c7-bf24da4dc996"], -Cell[118055, 3344, 28451, 690, 843, "Input",ExpressionUUID->"341cb469-2a0b-4b3f-a5c7-4a3d8fb61885", +Cell[115830, 3283, 1180, 31, 110, "Input",ExpressionUUID->"33dfcd59-33cc-4f88-bc22-62c17b91dbeb", InitializationCell->True], -Cell[146509, 4036, 227, 4, 35, "Text",ExpressionUUID->"5b2593a5-509e-44f9-a6a7-5b7dcf896d99"], -Cell[146739, 4042, 1855, 57, 68, "Input",ExpressionUUID->"3bf77406-a77c-45d7-8c2f-817ad440cdff", +Cell[117013, 3316, 880, 21, 89, "Input",ExpressionUUID->"9663902a-8667-481b-b613-f1c2b56b8e76", InitializationCell->True], +Cell[117896, 3339, 214, 4, 35, "Text",ExpressionUUID->"d38e2fe6-260d-4638-b639-6984dd9d14e4"], +Cell[118113, 3345, 679, 19, 63, "Input",ExpressionUUID->"36b985d7-30a2-430d-8dc3-30d6b4a45727", + InitializationCell->True], +Cell[118795, 3366, 210, 4, 35, "Text",ExpressionUUID->"e115f0cc-d9e8-4997-a5a4-88b0e9ed97e0"], +Cell[119008, 3372, 3586, 92, 187, "Input",ExpressionUUID->"f0a2825f-92fc-42d2-aeee-57f53535b1d5", + InitializationCell->True], +Cell[122597, 3466, 276, 5, 35, "Text",ExpressionUUID->"5353bb94-3a0f-4a57-a4c7-bf24da4dc996"], +Cell[122876, 3473, 28451, 690, 843, "Input",ExpressionUUID->"341cb469-2a0b-4b3f-a5c7-4a3d8fb61885", + InitializationCell->True], +Cell[151330, 4165, 227, 4, 35, "Text",ExpressionUUID->"5b2593a5-509e-44f9-a6a7-5b7dcf896d99"], +Cell[151560, 4171, 1855, 57, 68, "Input",ExpressionUUID->"3bf77406-a77c-45d7-8c2f-817ad440cdff", + InitializationCell->True], +Cell[153418, 4230, 264, 5, 35, "Text",ExpressionUUID->"8580055f-686c-403d-8800-a4eab889fce4"], +Cell[153685, 4237, 1223, 36, 46, "Input",ExpressionUUID->"4bf40614-2147-47eb-9ce0-13fe949a9c00", + InitializationCell->True], +Cell[154911, 4275, 1097, 30, 46, "Input",ExpressionUUID->"33c091a9-0792-4be2-bcd6-e02a256b26b5", + InitializationCell->True], +Cell[156011, 4307, 267, 5, 35, "Text",ExpressionUUID->"6e302dfa-38b2-43a3-a507-71bf2633d8e8"], +Cell[156281, 4314, 1612, 49, 66, "Input",ExpressionUUID->"41a0d75d-44b6-4653-b70e-858ab3e937b8", + InitializationCell->True], +Cell[157896, 4365, 2945, 90, 153, "Input",ExpressionUUID->"1430ab04-b05a-4d2d-bbf6-f9edc3cd5345", + InitializationCell->True], +Cell[160844, 4457, 1366, 45, 109, "Input",ExpressionUUID->"09f7d062-2fdd-45b1-8df2-0dd679b5cbc6", + InitializationCell->True], +Cell[162213, 4504, 258, 5, 35, "Text",ExpressionUUID->"94763911-d015-4b98-a645-b27f631322bf"], +Cell[162474, 4511, 685, 16, 65, "Input",ExpressionUUID->"46e4593c-f9e1-4f4a-a999-d0f67e325be6", + InitializationCell->True], +Cell[163162, 4529, 253, 5, 35, "Text",ExpressionUUID->"eeb33af8-ab09-4c5f-98a3-684ae529c2b6"], +Cell[163418, 4536, 2845, 78, 63, "Input",ExpressionUUID->"ab8163b6-fc27-4314-ae7e-a7f0f7d95629", + InitializationCell->True], +Cell[166266, 4616, 214, 4, 35, "Text",ExpressionUUID->"04f91dc9-a818-4dc9-a431-8d20ab592198"], +Cell[166483, 4622, 1311, 37, 110, "Input",ExpressionUUID->"78388547-17f1-4781-bf29-2487d17b5c53", + InitializationCell->True], +Cell[167797, 4661, 1512, 50, 46, "Input",ExpressionUUID->"1039bc21-357b-49d9-a514-fd3af5ae05a7", + InitializationCell->True], +Cell[169312, 4713, 884, 24, 46, "Input",ExpressionUUID->"01749816-6d54-416a-ab4f-bab2e420acb4", + InitializationCell->True], +Cell[170199, 4739, 1693, 55, 46, "Input",ExpressionUUID->"36f3a3dd-e13f-4bbf-a50c-b11a31e789f8", + InitializationCell->True], +Cell[171895, 4796, 301, 5, 35, "Text",ExpressionUUID->"acf62eaa-9c57-4e21-98ab-4e5f44a5d9f7"], +Cell[172199, 4803, 3140, 85, 110, "Input",ExpressionUUID->"4419a1ba-532a-442b-95d2-330c205bfd7e", + InitializationCell->True], +Cell[175342, 4890, 2108, 63, 46, "Input",ExpressionUUID->"c0b69368-df8e-46ab-a64b-3ea147722dd9", + InitializationCell->True], +Cell[177453, 4955, 919, 26, 46, "Input",ExpressionUUID->"0dd5826b-1b72-4f2c-834b-e79b419ac43e", + InitializationCell->True], +Cell[178375, 4983, 4206, 117, 159, "Input",ExpressionUUID->"61cc9c57-6cf1-464b-bd07-45e24fd7c9c3", + InitializationCell->True], +Cell[182584, 5102, 352, 6, 35, "Text",ExpressionUUID->"1e4ce167-3a3b-4c53-9c2a-63a8f4ac6e4e"], +Cell[182939, 5110, 3268, 88, 110, "Input",ExpressionUUID->"d4265afd-5fb1-4adc-ad27-da2413b01953", + InitializationCell->True], +Cell[186210, 5200, 2089, 61, 46, "Input",ExpressionUUID->"9cd135b3-7599-482f-9e19-dfd18b561c56", + InitializationCell->True], +Cell[188302, 5263, 923, 26, 46, "Input",ExpressionUUID->"99ae0dcb-f21c-46d4-84aa-7f172765d9d8", + InitializationCell->True], +Cell[189228, 5291, 4345, 121, 159, "Input",ExpressionUUID->"57a35adb-5f0e-4384-a38f-33a0653fb6b6", + InitializationCell->True], +Cell[193576, 5414, 3751, 102, 279, "Input",ExpressionUUID->"93054b36-78ec-4bb1-83ea-8fbb31835fb9", + InitializationCell->True], +Cell[197330, 5518, 1251, 33, 30, "Input",ExpressionUUID->"9eac41d8-7253-4bcb-9a2a-4a332b730ddb"], Cell[CellGroupData[{ -Cell[148619, 4103, 1463, 40, 56, "Input",ExpressionUUID->"92574105-d574-49ac-904d-a9f7b38ef1c3"], -Cell[150085, 4145, 511, 7, 34, "Output",ExpressionUUID->"aa0b75e7-d797-49d3-9a12-b14fac9e91cc"], -Cell[150599, 4154, 510, 7, 34, "Output",ExpressionUUID->"9defd3f6-f5d8-462a-a8a4-c0bc80c8daea"] +Cell[198606, 5555, 803, 22, 30, "Input",ExpressionUUID->"cc5802c8-5164-4d35-a4aa-ac0f06412d48"], +Cell[199412, 5579, 301, 5, 34, "Output",ExpressionUUID->"6c0f8654-3af4-43e9-b5c1-7a47f72b02d4"] }, Open ]], -Cell[151124, 4164, 264, 5, 35, "Text",ExpressionUUID->"8580055f-686c-403d-8800-a4eab889fce4"], -Cell[151391, 4171, 1223, 36, 46, "Input",ExpressionUUID->"4bf40614-2147-47eb-9ce0-13fe949a9c00", - InitializationCell->True], -Cell[CellGroupData[{ -Cell[152639, 4211, 1565, 45, 52, "Input",ExpressionUUID->"246fa925-65a3-41a1-b83a-54965c0376f6"], -Cell[154207, 4258, 335, 4, 34, "Output",ExpressionUUID->"e0356b94-048c-4f45-9ac7-db1db0236396"], -Cell[154545, 4264, 332, 4, 34, "Output",ExpressionUUID->"2b09c0f7-0efb-4358-85fb-939d40cdda47"] -}, Open ]], -Cell[154892, 4271, 1101, 31, 46, "Input",ExpressionUUID->"33c091a9-0792-4be2-bcd6-e02a256b26b5", - InitializationCell->True], -Cell[155996, 4304, 267, 5, 35, "Text",ExpressionUUID->"6e302dfa-38b2-43a3-a507-71bf2633d8e8"], -Cell[156266, 4311, 1612, 49, 66, "Input",ExpressionUUID->"41a0d75d-44b6-4653-b70e-858ab3e937b8", - InitializationCell->True], -Cell[157881, 4362, 2945, 90, 153, "Input",ExpressionUUID->"1430ab04-b05a-4d2d-bbf6-f9edc3cd5345", - InitializationCell->True], -Cell[160829, 4454, 1366, 45, 109, "Input",ExpressionUUID->"09f7d062-2fdd-45b1-8df2-0dd679b5cbc6", - InitializationCell->True], -Cell[162198, 4501, 214, 4, 35, "Text",ExpressionUUID->"04f91dc9-a818-4dc9-a431-8d20ab592198"], -Cell[162415, 4507, 1018, 31, 46, "Input",ExpressionUUID->"78388547-17f1-4781-bf29-2487d17b5c53", - InitializationCell->True], -Cell[163436, 4540, 1392, 48, 46, "Input",ExpressionUUID->"1039bc21-357b-49d9-a514-fd3af5ae05a7", - InitializationCell->True], -Cell[164831, 4590, 253, 5, 35, "Text",ExpressionUUID->"eeb33af8-ab09-4c5f-98a3-684ae529c2b6"], -Cell[165087, 4597, 2670, 76, 63, "Input",ExpressionUUID->"ab8163b6-fc27-4314-ae7e-a7f0f7d95629", - InitializationCell->True], -Cell[167760, 4675, 301, 5, 35, "Text",ExpressionUUID->"acf62eaa-9c57-4e21-98ab-4e5f44a5d9f7"], -Cell[168064, 4682, 2803, 80, 68, "Input",ExpressionUUID->"4419a1ba-532a-442b-95d2-330c205bfd7e", - InitializationCell->True], -Cell[170870, 4764, 1996, 59, 68, "Input",ExpressionUUID->"c0b69368-df8e-46ab-a64b-3ea147722dd9", - InitializationCell->True], -Cell[172869, 4825, 352, 6, 35, "Text",ExpressionUUID->"1e4ce167-3a3b-4c53-9c2a-63a8f4ac6e4e"], -Cell[173224, 4833, 2953, 83, 68, "Input",ExpressionUUID->"d4265afd-5fb1-4adc-ad27-da2413b01953", - InitializationCell->True], -Cell[176180, 4918, 2063, 61, 68, "Input",ExpressionUUID->"9cd135b3-7599-482f-9e19-dfd18b561c56", - InitializationCell->True] +Cell[199728, 5587, 896, 22, 52, "Input",ExpressionUUID->"8cbe4386-fbc6-4914-8e8c-e90d47f39535"] }, Open ]] } ] diff --git a/Manuscript/FarDFT.tex b/Manuscript/FarDFT.tex index ae59599..afe8570 100644 --- a/Manuscript/FarDFT.tex +++ b/Manuscript/FarDFT.tex @@ -628,9 +628,9 @@ At the (ground-state) LDA level (\ie, we only consider ground-state functionals) \end{subequations} with \begin{align} - \n{}{(0)}(\br{}) & = 2 \MO{1}{2}(\br{}), + \n{}{(0)}(\br{}) & = 2 [\MO{1}{\ew{}}(\br{})]^2, & - \n{}{(1)}(\br{}) & = 2 \MO{2}{2}(\br{}), + \n{}{(1)}(\br{}) & = 2 [\MO{2}{\ew{}}(\br{})]^2, \end{align} Note that, contrary to the HF case, self-interaction is present in LDA. %The KS orbital energies are given by