diff --git a/benzene.nb b/benzene.nb index 177ef26..daf204c 100644 --- a/benzene.nb +++ b/benzene.nb @@ -10,10 +10,10 @@ NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 158, 7] -NotebookDataLength[ 955487, 17961] -NotebookOptionsPosition[ 952722, 17906] -NotebookOutlinePosition[ 953115, 17922] -CellTagsIndexPosition[ 953072, 17919] +NotebookDataLength[ 967358, 18303] +NotebookOptionsPosition[ 962105, 18206] +NotebookOutlinePosition[ 962498, 18222] +CellTagsIndexPosition[ 962455, 18219] WindowFrame->Normal*) (* Beginning of Notebook Content *) @@ -81,8 +81,7 @@ Cell[BoxData[{ 3.807029969788039*^9}, {3.807243645653872*^9, 3.8072436465240803`*^9}, { 3.807261814034483*^9, 3.807261814866003*^9}, {3.807262188755821*^9, 3.807262188904146*^9}}, - CellLabel-> - "In[343]:=",ExpressionUUID->"fab9ee47-3e93-4201-bb7f-09ae5e949598"], + CellLabel->"In[73]:=",ExpressionUUID->"fab9ee47-3e93-4201-bb7f-09ae5e949598"], Cell[CellGroupData[{ @@ -280,8 +279,7 @@ Cell[BoxData[{ 3.807029983522843*^9, 3.807030011586466*^9}, {3.80715444508188*^9, 3.807154458718595*^9}, {3.807240495668721*^9, 3.807240496242681*^9}, { 3.8072632505601*^9, 3.807263267134137*^9}}, - CellLabel-> - "In[480]:=",ExpressionUUID->"070c4b28-8c43-4bc2-ad60-9265a5b4fc44"], + CellLabel->"In[75]:=",ExpressionUUID->"070c4b28-8c43-4bc2-ad60-9265a5b4fc44"], Cell[BoxData[ TagBox[ @@ -488,26 +486,7 @@ Cell[BoxData[ 16.8586758847529, -838.2441671000151}}, {{ 17.551823065312842`, -841.7734071000211}}, {{ 18.24497024587279, -844.5056871000247}}, {{ - 18.938117426432733`, -847.6478071000031}}}]}}, { - {GrayLevel[0], PointSize[0.011000000000000001`], Thickness[Large]}, - {RGBColor[0, - NCache[ - Rational[2, 3], 0.6666666666666666], 0], PointSize[ - 0.011000000000000001`], Thickness[Large]}, - {RGBColor[1, 0, 0], PointSize[0.011000000000000001`], Thickness[ - Large]}, - {RGBColor[0, 0, 1], PointSize[0.011000000000000001`], Thickness[ - Large]}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}}, { - {GrayLevel[0], PointSize[0.011000000000000001`], Thickness[Large]}, - {RGBColor[0, - NCache[ - Rational[2, 3], 0.6666666666666666], 0], PointSize[ - 0.011000000000000001`], Thickness[Large]}, - {RGBColor[1, 0, 0], PointSize[0.011000000000000001`], Thickness[ - Large]}, - {RGBColor[0, 0, 1], PointSize[0.011000000000000001`], Thickness[ - Large]}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}}}, {{}, \ -{}}}, InsetBox[ + 18.938117426432733`, -847.6478071000031}}}]}}}, {{}, {}}}, InsetBox[ TemplateBox[{ GraphicsBox[{ Thickness[0.00970591089973794], @@ -965,6 +944,7 @@ Pdgq4n+7vAN6fAMAzM/UBw== 109.74700000000001`, 11.8656}, {109.49200000000002`, 11.882799999999998`}, {109.07299999999998`, 11.9203}, { 109.07299999999998`, 16.6484}}}], + FilledCurveBox[{{{0, 2, 0}, {1, 3, 3}, {0, 1, 0}, {0, 1, 0}, {0, 1, 0}, {1, 3, 3}, {0, 1, 0}, {0, 1, 0}, {1, 3, 3}, {0, 1, 0}, {0, 1, 0}, {0, 1, 0}, {0, 1, 0}, {0, 1, 0}, {1, 3, 3}, {0, 1, @@ -1980,6 +1960,7 @@ etkE1OH856D4+msA50+fwF9lZm3ssGHu++XHDkc59Hi9YjGZaAJxbxCC/wJo /RoeBN/HvNMx4W6kgwY4fZk4MAODa1NppMPhtuXhpxYZwfng+GgxcYjb5cnD lB0J1w+0Xc1QARrfyibw9McAAg+MHWDp+fRhp7WZdQg+LH0DAOVUPzw= "]], + FilledCurveBox[{{{1, 4, 3}, {1, 3, 3}, {1, 3, 3}, {0, 1, 0}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {0, 1, 0}}, {{1, 4, 3}, {1, 3, 3}, {0, 1, 0}}}, {CompressedData[" @@ -2631,7 +2612,6 @@ aPotByJer+eg0y5289ziHIfrQp8czx8zcAi3BAZAX44Dr//6Kakahg6xwOzC Automatic, Automatic}, {{{6.907755278982137, FormBox["1000", TraditionalForm], {0.01, 0.}}, {9.210340371976184, FormBox[ - TemplateBox[{"10", "4"}, "Superscript", SyntaxForm -> SuperscriptBox], TraditionalForm], {0.01, 0.}}, { 11.512925464970229`, @@ -3310,52 +3290,8 @@ aPotByJer+eg0y5289ziHIfrQp8czx8zcAi3BAZAX44Dr//6Kakahg6xwOzC 16.8586758847529, -838.2441671000151}}, {{ 17.551823065312842`, -841.7734071000211}}, {{ 18.24497024587279, -844.5056871000247}}, {{ - 18.938117426432733`, -847.6478071000031}}}]}}, {{ - Directive[ - PointSize[0.011000000000000001`], - AbsoluteThickness[1.6], - Thickness[Large], - GrayLevel[0]]}, { - Directive[ - PointSize[0.011000000000000001`], - AbsoluteThickness[1.6], - Thickness[Large], - RGBColor[0, - Rational[2, 3], 0]]}, { - Directive[ - PointSize[0.011000000000000001`], - AbsoluteThickness[1.6], - Thickness[Large], - RGBColor[1, 0, 0]]}, { - Directive[ - PointSize[0.011000000000000001`], - AbsoluteThickness[1.6], - Thickness[Large], - RGBColor[ - 0, 0, 1]]}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}}, {{ - Directive[ - PointSize[0.011000000000000001`], - AbsoluteThickness[1.6], - Thickness[Large], - GrayLevel[0]]}, { - Directive[ - PointSize[0.011000000000000001`], - AbsoluteThickness[1.6], - Thickness[Large], - RGBColor[0, - Rational[2, 3], 0]]}, { - Directive[ - PointSize[0.011000000000000001`], - AbsoluteThickness[1.6], - Thickness[Large], - RGBColor[1, 0, 0]]}, { - Directive[ - PointSize[0.011000000000000001`], - AbsoluteThickness[1.6], - Thickness[Large], - RGBColor[ - 0, 0, 1]]}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}}}, \ -{{}, {}}}, {DisplayFunction -> Identity, GridLines -> { + 18.938117426432733`, -847.6478071000031}}}]}}}, {{}, {}}}, { + DisplayFunction -> Identity, GridLines -> { Charting`ScaledTickValues[{Log, Exp}], Automatic}, DisplayFunction -> Identity, DisplayFunction -> Identity, DisplayFunction -> Identity, DisplayFunction -> Identity, AspectRatio -> 1, Axes -> {False, False}, @@ -3562,8 +3498,8 @@ nyDsg/Fh7oHxS7aK/j5dZ+KAHp4AGg2gZQ== FilledCurve[{{{1, 4, 3}, {1, 3, 3}, {0, 1, 0}, {1, 3, 3}, {0, 1, 0}, {0, 1, 0}, {0, 1, 0}, {1, 3, 3}, {0, 1, 0}, {1, 3, 3}, {1, 3, 3}, {0, 1, 0}, {1, 3, 3}, {0, 1, 0}, {0, 1, 0}, {0, 1, 0}, { - 1, 3, 3}, {0, 1, 0}, {1, 3, 3}, {1, 3, 3}, {0, 1, 0}, {0, 1, 0}, - {1, 3, 3}, {0, 1, 0}}}, CompressedData[" + 1, 3, 3}, {0, 1, 0}, {1, 3, 3}, {1, 3, 3}, {0, 1, 0}, {0, 1, + 0}, {1, 3, 3}, {0, 1, 0}}}, CompressedData[" 1:eJxTTMoPSmVmYGBgBGJ9IGYC4v21shbpMyIcrgt9cjwvZuwwEwQWRjiov+Xd Z6Bp7HAGBDZEODxJXHjNRB8qfxDBXzTJitH3DYJ/+LJ2quSvCIfF9/fxzbls BOfz+q+fkrrCEM5/kaX9bXqtOpyfmgYEx1QcvrHFz/B5GeFw4WrYG31pFYdr @@ -3954,6 +3890,7 @@ g3d1s76PXZCDY9Oj4zN2GzvA0t/nDQHZs9wRfFh6BABafjoZ 16.156299999999998`}, {88.13749999999999, 18.421899999999994`}, {89.2578, 19.900000000000002`}, { 90.94999999999999, 19.900000000000002`}}}], + FilledCurve[{{{0, 2, 0}, {0, 1, 0}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {0, 1, 0}, {0, 1, 0}, {0, 1, 0}, {0, 1, 0}, {1, 3, 3}, {0, 1, 0}, {0, 1, 0}, { @@ -4133,6 +4070,7 @@ oPj6awDnT5/AX2VmbexQlDHxbY12hkOP1ysWk4kmDjdA7v2cDueDrH9/HsF3 MI3b5bkj3UEDnL5MHOYAg6uoJd3hcNvy8FOLjOB8cHy0mDgAU9d8qXqEfqDt upsiof5VNoGnPwYQeGDsAEvPpw87rc2sQ/Bh6RsATVY9FQ== "]], + FilledCurve[{{{1, 4, 3}, {1, 3, 3}, {0, 1, 0}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {0, 1, 0}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, { @@ -5083,9 +5021,8 @@ W1zuAg== 3.807240481625876*^9, 3.8072404967028103`*^9}, 3.807241599596191*^9, 3.807243658126053*^9, 3.8072615834249372`*^9, 3.807261817092359*^9, 3.807262092857091*^9, 3.807262190577567*^9, 3.807262358869186*^9, { - 3.807263251662565*^9, 3.807263268017001*^9}}, - CellLabel-> - "Out[489]=",ExpressionUUID->"d0d8ae32-ad10-4dcf-9f77-479e69e9ec0e"] + 3.807263251662565*^9, 3.807263268017001*^9}, 3.807356529851677*^9}, + CellLabel->"Out[84]=",ExpressionUUID->"3dc7aee2-dfb7-4266-8af9-123c83a80027"] }, Open ]], Cell[CellGroupData[{ @@ -8998,9 +8935,10 @@ Cell[BoxData[{ 3.807025314769651*^9}, {3.8070254940806017`*^9, 3.807025494241976*^9}, { 3.80702577364151*^9, 3.807025793957727*^9}, {3.807029968143608*^9, 3.807029969788039*^9}, {3.807243645653872*^9, 3.8072436465240803`*^9}, { - 3.807261814034483*^9, 3.807261814866003*^9}}, + 3.807261814034483*^9, 3.807261814866003*^9}, {3.807356542943404*^9, + 3.8073565589786997`*^9}}, CellLabel-> - "In[502]:=",ExpressionUUID->"88b83480-ff85-479c-9bd5-86e50cd48da6"], + "In[166]:=",ExpressionUUID->"88b83480-ff85-479c-9bd5-86e50cd48da6"], Cell[CellGroupData[{ @@ -9202,7 +9140,7 @@ Cell[BoxData[{ 3.807154458718595*^9}, {3.807240495668721*^9, 3.807240496242681*^9}, 3.807263260480081*^9}, CellLabel-> - "In[504]:=",ExpressionUUID->"f0f1eea9-20b9-48df-9572-c510784075b4"], + "In[168]:=",ExpressionUUID->"f0f1eea9-20b9-48df-9572-c510784075b4"], Cell[BoxData[ TagBox[ @@ -9411,26 +9349,7 @@ Cell[BoxData[ 16.8586758847529, -842.4038671000176}}, {{ 17.551823065312842`, -845.2500170999997}}, {{ 18.24497024587279, -848.0130471000109}}, {{ - 18.938117426432733`, -849.9845171000118}}}]}}, { - {GrayLevel[0], PointSize[0.011000000000000001`], Thickness[Large]}, - {RGBColor[0, - NCache[ - Rational[2, 3], 0.6666666666666666], 0], PointSize[ - 0.011000000000000001`], Thickness[Large]}, - {RGBColor[1, 0, 0], PointSize[0.011000000000000001`], Thickness[ - Large]}, - {RGBColor[0, 0, 1], PointSize[0.011000000000000001`], Thickness[ - Large]}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}}, { - {GrayLevel[0], PointSize[0.011000000000000001`], Thickness[Large]}, - {RGBColor[0, - NCache[ - Rational[2, 3], 0.6666666666666666], 0], PointSize[ - 0.011000000000000001`], Thickness[Large]}, - {RGBColor[1, 0, 0], PointSize[0.011000000000000001`], Thickness[ - Large]}, - {RGBColor[0, 0, 1], PointSize[0.011000000000000001`], Thickness[ - Large]}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}}}, {{}, \ -{}}}, InsetBox[ + 18.938117426432733`, -849.9845171000118}}}]}}}, {{}, {}}}, InsetBox[ TemplateBox[{ GraphicsBox[{ Thickness[0.00970591089973794], @@ -10833,6 +10752,7 @@ TJDW 12.5797}, {65.07969999999999, 12.674999999999999`}, { 65.07969999999999, 14.0828}, {66.0578, 14.940599999999998`}, { 68.9422, 15.9891}, {68.9422, 12.532799999999998`}}}], + FilledCurveBox[{{{0, 2, 0}, {0, 1, 0}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {0, 1, 0}, {0, 1, 0}, {1, 3, 3}, {1, 3, 3}, {0, 1, 0}, {1, 3, 3}, {1, 3, 3}, {0, 1, 0}, {0, 1, 0}, {0, 1, @@ -11324,6 +11244,7 @@ dDelIvjuD6tE1s2PdSg/vM115lpluPth8QXjl4Lk7yo5oKcXAJ2oB70= 19.900000000000002`}, {119.30599999999998`, 18.7781}, { 119.30599999999998`, 17.515600000000003`}, {119.30599999999998`, 12.0313}}}], + FilledCurveBox[{{{1, 4, 3}, {1, 3, 3}, {1, 3, 3}, {0, 1, 0}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {0, 1, 0}}, {{1, 4, 3}, {1, 3, 3}, {0, 1, 0}}}, {CompressedData[" @@ -11554,7 +11475,6 @@ aPotByJer+eg0y5289ziHIfrQp8czx8zcAi3BAZAX44Dr//6Kakahg6xwOzC Automatic, Automatic}, {{{6.907755278982137, FormBox["1000", TraditionalForm], {0.01, 0.}}, {9.210340371976184, FormBox[ - TemplateBox[{"10", "4"}, "Superscript", SyntaxForm -> SuperscriptBox], TraditionalForm], {0.01, 0.}}, { 11.512925464970229`, @@ -12235,52 +12155,8 @@ aPotByJer+eg0y5289ziHIfrQp8czx8zcAi3BAZAX44Dr//6Kakahg6xwOzC 16.8586758847529, -842.4038671000176}}, {{ 17.551823065312842`, -845.2500170999997}}, {{ 18.24497024587279, -848.0130471000109}}, {{ - 18.938117426432733`, -849.9845171000118}}}]}}, {{ - Directive[ - PointSize[0.011000000000000001`], - AbsoluteThickness[1.6], - Thickness[Large], - GrayLevel[0]]}, { - Directive[ - PointSize[0.011000000000000001`], - AbsoluteThickness[1.6], - Thickness[Large], - RGBColor[0, - Rational[2, 3], 0]]}, { - Directive[ - PointSize[0.011000000000000001`], - AbsoluteThickness[1.6], - Thickness[Large], - RGBColor[1, 0, 0]]}, { - Directive[ - PointSize[0.011000000000000001`], - AbsoluteThickness[1.6], - Thickness[Large], - RGBColor[ - 0, 0, 1]]}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}}, {{ - Directive[ - PointSize[0.011000000000000001`], - AbsoluteThickness[1.6], - Thickness[Large], - GrayLevel[0]]}, { - Directive[ - PointSize[0.011000000000000001`], - AbsoluteThickness[1.6], - Thickness[Large], - RGBColor[0, - Rational[2, 3], 0]]}, { - Directive[ - PointSize[0.011000000000000001`], - AbsoluteThickness[1.6], - Thickness[Large], - RGBColor[1, 0, 0]]}, { - Directive[ - PointSize[0.011000000000000001`], - AbsoluteThickness[1.6], - Thickness[Large], - RGBColor[ - 0, 0, 1]]}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}}}, \ -{{}, {}}}, {DisplayFunction -> Identity, GridLines -> { + 18.938117426432733`, -849.9845171000118}}}]}}}, {{}, {}}}, { + DisplayFunction -> Identity, GridLines -> { Charting`ScaledTickValues[{Log, Exp}], Automatic}, DisplayFunction -> Identity, DisplayFunction -> Identity, DisplayFunction -> Identity, DisplayFunction -> Identity, AspectRatio -> 1, Axes -> {False, False}, @@ -12384,7 +12260,6 @@ A3p8AgDLMM+X 19.707800000000002`}, {51.7141, 19.040599999999998`}, { 52.02340000000001, 16.967200000000002`}, {47.1609, 16.967200000000002`}}}], - FilledCurve[{{{0, 2, 0}, {1, 3, 3}, {1, 3, 3}, {0, 1, 0}, {1, 3, 3}, {0, 1, 0}, {0, 1, 0}, {0, 1, 0}, {1, 3, 3}, {0, 1, 0}, {0, 1, 0}, {1, 3, 3}, {0, 1, 0}}}, CompressedData[" @@ -12479,8 +12354,9 @@ nyDsg/Fh7oHxS7aK/j5dZ+KAHp4AGg2gZQ== 91.2688, 10.0297}, {90.38749999999999, 10.6016}, { 89.88589999999999, 11.5547}, {89.2188, 12.7953}, { 88.83749999999999, 14.464099999999997`}, {88.83749999999999, - 16.156299999999998`}, {88.83749999999999, 18.421899999999994`}, - {89.95780000000002, 19.900000000000002`}, {91.64999999999999, + 16.156299999999998`}, {88.83749999999999, + 18.421899999999994`}, {89.95780000000002, + 19.900000000000002`}, {91.64999999999999, 19.900000000000002`}}}], FilledCurve[{{{1, 4, 3}, {1, 3, 3}, {0, 1, 0}, {1, 3, 3}, {0, 1, @@ -14008,9 +13884,11 @@ W1zuAg== 3.807088359348612*^9}, {3.807154441293125*^9, 3.8071544593463383`*^9}, { 3.807240481625876*^9, 3.8072404967028103`*^9}, 3.807241599596191*^9, 3.807243658126053*^9, 3.8072615834249372`*^9, 3.807261817092359*^9, - 3.807262092857091*^9, 3.807263261194553*^9, 3.8072632916131144`*^9}, + 3.807262092857091*^9, 3.807263261194553*^9, 3.8072632916131144`*^9, + 3.807355433753173*^9, {3.807356544437235*^9, 3.807356571477729*^9}, + 3.807356695973937*^9}, CellLabel-> - "Out[513]=",ExpressionUUID->"0de2ef17-1213-4011-ac72-c6e9e0caf6db"] + "Out[177]=",ExpressionUUID->"f08a5516-ecb9-4b9d-b128-154a3ac2cb7b"] }, Open ]], Cell[CellGroupData[{ @@ -14179,7 +14057,7 @@ $E_\\\\text{rPT2}$}\>\"", ",", RowBox[{"{", RowBox[{ RowBox[{"-", "600"}], ",", - RowBox[{"-", "900"}]}], "}"}]}], "}"}]}], ",", + RowBox[{"-", "880"}]}], "}"}]}], "}"}]}], ",", "\[IndentingNewLine]", RowBox[{"Joined", "\[Rule]", "False"}]}], "\[IndentingNewLine]", "]"}], "\[IndentingNewLine]", ",", "\[IndentingNewLine]", @@ -14202,25 +14080,25 @@ $E_\\\\text{rPT2}$}\>\"", ",", RowBox[{"Export", "[", RowBox[{"\"\\"", ",", "%"}], "]"}], ";"}]}], "Input", CellChangeTimes->CompressedData[" -1:eJwdznss1AEAB/C7JFYWSSXsLqe7yKMzlN1woeRRlMeinVwPdsUIuUyNWB7N -GW263cijzIizuJmhhqmd6CyPW1wT53Bep4UdjWn9vr8/vvvs+893X9u7KeHx -+ygUig8R+Ld8Ztf2qs43ejnUmkFYr3c7BQWrK85wYvugHwwU0UkvcbWXobGV -NBjmhk3egMsVN+9BS9uwVCj7yS2FNK5hOZy27xFDB3bLAPxM2RyCcq7K7DSh -mXTpKBxtyz0JJStm1jCOSWHALcVrF7hw4YQb3G8x4gFbKwo4kO7dFQIn+tJi -IINGvQU7V9r5sKjjdjwc+becBnt01HTI3igqhYLoXTHUrugrocxO8Z78dayh -BfJUzaS5RTsOTOy3j52FmiqVGzSK43nBpRmb1H7CgkfrT6DcIuMZVPQKcqB7 -zfALWPtbnQ8bxW3F0CTrqwj6b4iHLcN0vlKh+SjMVkYc1vN0vt7yYCu4ZJrP -gCwTc0eYUK9zgY6u9CC4RqVfg92bBqGwXCKJhHlj6VGwQMiNgeHxHD6cHEp+ -AI2dx4WwNUSWBfMTPfLhYIK8BKbxJ2vhxZmBDrJzFd/g45cl32H120FSD7W7 -Eq42RU7BKtNEDcwM8/8DZXWCjC1Co4YgITzQtJcLz6uO5JE9O64Qahc6SfMa -WGVQeb3xFZw7Uy2BixHTpOI68374QzhLmhT8YXCbUMSZJe0SMUdgDd2d9OOd -WWufWGI35YoLfJN43w2GW815wl+xnjYRhKP6h3TYzNuZjyJ0ilEvwuLkTLaA -kN/NcoUsp4DpVEKlJod0bV6qgQGNQfNQaxfg8Imw7+kOqVYf2ArL+qgy6NNb -PAy/JE2RcvwMhHJC9XhSIbR3XCyGDP91Edxr85TAnPKBSvicxjYcIPTinjsE -me82zSA/k3Yc/gcxKftB +1:eJwdznss1AEAB/C7JFaWSyphjpOLPGJHmeEXSh5FeSzs5HqwK0bIZWrEOjRn +tF23G3mUGXksbmaoYWp3EctroYnzOq/TQkdjqt/398d3n33/+e5rcTspJHYf +jUYjyMA/4uldiytqr4jlIBMWabWGYw75qyv2cGz7oDf0EzEpLxKqS1DXuCEA +ZgePX4fLJTfuQCOL4GQo+04UQTNCWwynrDsl0MaxUQE/0jb7oJyYYJwiZTQs +HYVDzdknoXSFYQJjrGgsuNX/0gEunD/BgfsNB11gU0muG2R6tAfCse6USMgy +o0fBtpUWHsxvvRkLB/8up8BONT0VOm7kF0F+xK4EqlY0pVBm2f+W+nWsphFy +J+ops/N3bKyw3zJ8Bs6WTXCgTgzXHS5Nmyb3kOY+WH8E5YZpT2B/Fz8LOlcM +PIOVP5VCWCtpLoB6GZ9F0GdDMmAUrPZqEBgMwczR0MMartrLQx5gDJf0hSzI +1jOwhXHVagdo68T0h2t05lXYsakVBMVSaRjMGU4Nh7kCIhKGxLrx4Hhf4j2o +az8igE2BsgwojHcRwt44eSFM4Y1XwgvTilaqE/1f4MPnhV9h+eteShel8yhc +rQubhGX68bMwPdjnF5RV8dO2SHVq/AXwQN1eNjw3cSSH6pkxeVC10EaZU8Mu +hqPXal/AudPlUrgYOkUpqTLogd8EM5QJAe96t0lFbjOU7SKrQVjBdKZ8f2vG +xDOa3E267ABfxd/lwBDjOVf4I9rVNJR0SHOfCeu5O/PhpHaRykVYkJjuyCfl +dbCdINvOdyqZdHQ2i3JtvmEW+tb6z0OVpa/NB9LuxzuUKo1fEyzupsugZ1fB +APyUMEnp5q0lkJMqRxLyoLXtYgFk+ayL4F6zqxRmiRWl8KmZo7aC1J04ewha +vdlkQF662XEYVWoe9ps06Pk/yv8ulwPg "], CellLabel-> - "In[515]:=",ExpressionUUID->"0daea3d8-8f63-4ebd-ba0f-4b294d15916e"], + "In[179]:=",ExpressionUUID->"0daea3d8-8f63-4ebd-ba0f-4b294d15916e"], Cell[BoxData[ RowBox[{ @@ -14243,9 +14121,11 @@ Cell[BoxData[ 3.8072617557205467`*^9}, {3.80726179986965*^9, 3.807261818038783*^9}, { 3.807261896511039*^9, 3.807261918612938*^9}, 3.807262054107587*^9, 3.807262093637086*^9, 3.807262792191979*^9, {3.807262858065199*^9, - 3.807262897281643*^9}, 3.807263298877152*^9}, + 3.807262897281643*^9}, 3.807263298877152*^9, 3.807355434877316*^9, { + 3.807355555438017*^9, 3.807355560511257*^9}, {3.8073565469916763`*^9, + 3.807356571967432*^9}, 3.807356696395651*^9}, CellLabel-> - "Out[523]=",ExpressionUUID->"7a77a453-90c2-4388-b095-5d7536c2ebd6"], + "Out[187]=",ExpressionUUID->"70b3141b-5f73-46ce-a120-e28e26a497f7"], Cell[BoxData[ RowBox[{ @@ -14268,9 +14148,11 @@ Cell[BoxData[ 3.8072617557205467`*^9}, {3.80726179986965*^9, 3.807261818038783*^9}, { 3.807261896511039*^9, 3.807261918612938*^9}, 3.807262054107587*^9, 3.807262093637086*^9, 3.807262792191979*^9, {3.807262858065199*^9, - 3.807262897281643*^9}, 3.807263298880394*^9}, + 3.807262897281643*^9}, 3.807263298877152*^9, 3.807355434877316*^9, { + 3.807355555438017*^9, 3.807355560511257*^9}, {3.8073565469916763`*^9, + 3.807356571967432*^9}, 3.807356696396323*^9}, CellLabel-> - "Out[525]=",ExpressionUUID->"3d50c171-50af-40d7-83b3-f304118a6b8b"], + "Out[189]=",ExpressionUUID->"b622dec9-d93a-44a3-973b-53c6aa852652"], Cell[BoxData[ TagBox[ @@ -14388,17 +14270,7 @@ Cell[BoxData[ {{-0.06438604999999598, -775.6955871000173}}, {{-0.055509260000008, \ -786.8946071000096}}, {{-0.04794706999999221, -797.3029471000075}}, \ {{-0.04126446000000783, -806.7485871000031}}, {{-0.03514844000000039, \ --814.8360771000114}}}]}}, { - {GrayLevel[0], PointSize[0.012833333333333334`], Thickness[Large]}, - {RGBColor[0, 0, 1], PointSize[0.012833333333333334`], Thickness[ - Large]}, - {RGBColor[1, 0, 0], PointSize[0.012833333333333334`], Thickness[ - Large]}, {}, {}, {}, {}, {}, {}, {}}, { - {GrayLevel[0], PointSize[0.012833333333333334`], Thickness[Large]}, - {RGBColor[0, 0, 1], PointSize[0.012833333333333334`], Thickness[ - Large]}, - {RGBColor[1, 0, 0], PointSize[0.012833333333333334`], Thickness[ - Large]}, {}, {}, {}, {}, {}, {}, {}}}, {{}, {}}}, {{{}, {}, +-814.8360771000114}}}]}}}, {{}, {}}}, {{{}, {}, TagBox[ {RGBColor[0, 0, 1], Thickness[Large], Opacity[1.], Dashing[{Small, Small}], LineBox[CompressedData[" @@ -14431,7 +14303,7 @@ ZWhpoDNyovPLOwJNhrhe7IRtsIdeya5VE/hy2azu4Jwb5Hd+dYh7T2Bi6gT/ N60jOEFuuuYDgQVp8pR3c3bAdh72qJoi0JVc++0drTUoSh1VZz4SKB0PWPVM awTLKRevbNIQWDQmO/s0ybHhZO7j3ep5/weUAqmA "]]}, - Annotation[#, "Charting`Private`Tag$29186#1"]& ], + Annotation[#, "Charting`Private`Tag$16402#1"]& ], TagBox[ {RGBColor[1, 0, 0], Thickness[Large], Opacity[1.], Dashing[{Small, Small}], LineBox[CompressedData[" @@ -14464,7 +14336,7 @@ Oph0d2h3zhBcM8i2jYs+CEH2xvlJXwgOrVtYMiC3hPt1lJnzLMEY4SfqtmI3 6OdaOvXMEbzzuyRhUm4EHmtvDtz6SnC/Su1PRYrNkC6YCHH7RrD4o/PGlwoN 6F6wX9CdJ5jzoeZiG293w5rwPNK/2P8BZ3Cl8Q== "]]}, - Annotation[#, "Charting`Private`Tag$29186#2"]& ]}, {}}}, InsetBox[ + Annotation[#, "Charting`Private`Tag$16402#2"]& ]}, {}}}, InsetBox[ TemplateBox[{ GraphicsBox[{ Thickness[0.00970591089973794], @@ -16149,7 +16021,6 @@ uH2v23K7jX5LwPkzwUDcYT9Q+bbVURh8mHqddrGb5/yjHIzBQNIBFp8o9gH5 108.10199999999998`, 11.8656}, {106.01099999999997`, 11.8656}, { 105.64699999999999`, 11.8656}, {105.392, 11.882799999999998`}, { 104.97300000000001`, 11.9203}, {104.97300000000001`, 16.6484}}}], - FilledCurveBox[{{{0, 2, 0}, {1, 3, 3}, {0, 1, 0}, {0, 1, 0}, {0, 1, 0}, {1, 3, 3}, {0, 1, 0}, {0, 1, 0}, {1, 3, 3}, {0, 1, 0}, {0, 1, 0}, {0, 1, 0}, {0, 1, 0}, {0, 1, 0}, {1, 3, 3}, {0, 1, @@ -16202,7 +16073,7 @@ j38AhtbtwA== Part[#, 1]], Identity[ Part[#, 2]]}& )}}, - PlotRange->{{-0.2, 0.}, {-900, -600}}, + PlotRange->{{-0.2, 0.}, {-880, -600}}, PlotRangeClipping->True, PlotRangePadding->{{0, 0}, {0, 0}}, Ticks->{Automatic, Automatic}], @@ -16353,38 +16224,7 @@ j38AhtbtwA== {{-0.06438604999999598, -775.6955871000173}}, {{-0.055509260000008, \ -786.8946071000096}}, {{-0.04794706999999221, -797.3029471000075}}, \ {{-0.04126446000000783, -806.7485871000031}}, {{-0.03514844000000039, \ --814.8360771000114}}}]}}, {{ - Directive[ - PointSize[0.012833333333333334`], - AbsoluteThickness[1.6], - Thickness[Large], - GrayLevel[0]]}, { - Directive[ - PointSize[0.012833333333333334`], - AbsoluteThickness[1.6], - Thickness[Large], - RGBColor[0, 0, 1]]}, { - Directive[ - PointSize[0.012833333333333334`], - AbsoluteThickness[1.6], - Thickness[Large], - RGBColor[1, 0, 0]]}, {}, {}, {}, {}, {}, {}, {}}, {{ - Directive[ - PointSize[0.012833333333333334`], - AbsoluteThickness[1.6], - Thickness[Large], - GrayLevel[0]]}, { - Directive[ - PointSize[0.012833333333333334`], - AbsoluteThickness[1.6], - Thickness[Large], - RGBColor[0, 0, 1]]}, { - Directive[ - PointSize[0.012833333333333334`], - AbsoluteThickness[1.6], - Thickness[Large], - RGBColor[ - 1, 0, 0]]}, {}, {}, {}, {}, {}, {}, {}}}, {{}, {}}}, {{{{}, {}, +-814.8360771000114}}}]}}}, {{}, {}}}, {{{{}, {}, Annotation[{ Directive[ Opacity[1.], @@ -16422,7 +16262,7 @@ ZWhpoDNyovPLOwJNhrhe7IRtsIdeya5VE/hy2azu4Jwb5Hd+dYh7T2Bi6gT/ N60jOEFuuuYDgQVp8pR3c3bAdh72qJoi0JVc++0drTUoSh1VZz4SKB0PWPVM awTLKRevbNIQWDQmO/s0ybHhZO7j3ep5/weUAqmA "]]}, - "Charting`Private`Tag$29186#1"], + "Charting`Private`Tag$16402#1"], Annotation[{ Directive[ Opacity[1.], @@ -16460,7 +16300,7 @@ Oph0d2h3zhBcM8i2jYs+CEH2xvlJXwgOrVtYMiC3hPt1lJnzLMEY4SfqtmI3 6OdaOvXMEbzzuyRhUm4EHmtvDtz6SnC/Su1PRYrNkC6YCHH7RrD4o/PGlwoN 6F6wX9CdJ5jzoeZiG293w5rwPNK/2P8BZ3Cl8Q== "]]}, - "Charting`Private`Tag$29186#2"]}}, {}}}, { + "Charting`Private`Tag$16402#2"]}}, {}}}, { DisplayFunction -> Identity, DisplayFunction -> Identity, AspectRatio -> 1, Axes -> {False, False}, AxesLabel -> {None, None}, AxesOrigin -> {0, -600.}, BaseStyle -> 18, DisplayFunction :> Identity, @@ -16719,6 +16559,7 @@ Ef6D8WH+310CDFGnRAw+TD04fGUSHTxADlJXgYSXcCIkXnci+Cd27ehlE1CH 2UcUNtQlOGiA05eJAyMwuG56JTgcblsefmqREZwPjo8WEwdQ6mp0RejP3gM0 UQHqX2UTePpjAIEHxg6w9Hz6sNPazDoEH5a+ATd6Nk4= "]], + FilledCurve[{{{1, 4, 3}, {1, 3, 3}, {1, 3, 3}, {0, 1, 0}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {0, 1, 0}}, {{1, 4, 3}, {1, 3, 3}, {0, 1, 0}}}, {CompressedData[" @@ -17147,7 +16988,7 @@ j38AhtbtwA== Identity[ Part[#, 1]], Identity[ - Part[#, 2]]}& )}}, PlotRange -> {{-0.2, 0.}, {-900, -600}}, + Part[#, 2]]}& )}}, PlotRange -> {{-0.2, 0.}, {-880, -600}}, PlotRangeClipping -> True, PlotRangePadding -> {{0, 0}, {0, 0}}, Ticks -> {Automatic, Automatic}}], Placed[ @@ -17877,9 +17718,11 @@ WOz7Hw== 3.8072617557205467`*^9}, {3.80726179986965*^9, 3.807261818038783*^9}, { 3.807261896511039*^9, 3.807261918612938*^9}, 3.807262054107587*^9, 3.807262093637086*^9, 3.807262792191979*^9, {3.807262858065199*^9, - 3.807262897281643*^9}, 3.80726329907366*^9}, + 3.807262897281643*^9}, 3.807263298877152*^9, 3.807355434877316*^9, { + 3.807355555438017*^9, 3.807355560511257*^9}, {3.8073565469916763`*^9, + 3.807356571967432*^9}, 3.807356696503303*^9}, CellLabel-> - "Out[526]=",ExpressionUUID->"5b22d62c-c18f-4d0b-b5f0-2d98e0faf6a0"] + "Out[190]=",ExpressionUUID->"8d3b739a-d028-4230-b557-d5db482643b2"] }, Open ]], Cell[CellGroupData[{ @@ -17893,20 +17736,477 @@ Cell[BoxData[ CellChangeTimes->{{3.807265686484186*^9, 3.8072657251931047`*^9}, 3.8072658656281633`*^9, {3.807265980818079*^9, 3.807265980953951*^9}}, CellLabel-> - "In[531]:=",ExpressionUUID->"ba14da0a-87ae-46da-85c2-b124d0f9c105"], + "In[192]:=",ExpressionUUID->"ba14da0a-87ae-46da-85c2-b124d0f9c105"], Cell[BoxData["159744"], "Output", CellChangeTimes->{3.807265726567082*^9, 3.807265930898308*^9, - 3.807265981285715*^9}, + 3.807265981285715*^9, 3.8073554375827417`*^9, 3.8073565723491163`*^9, + 3.807356696769352*^9}, CellLabel-> - "Out[531]=",ExpressionUUID->"a9ea879d-0902-454b-a505-60f0edd2cffd"] + "Out[192]=",ExpressionUUID->"172c6c74-0ddd-4b88-801f-ed238e92e91d"] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"40", " ", "24", " ", "128"}]], "Input", + CellChangeTimes->{{3.807355755448935*^9, 3.807355761392128*^9}}, + CellLabel-> + "In[193]:=",ExpressionUUID->"3e45b97c-e088-4c82-99ef-15e7d70333ed"], + +Cell[BoxData["122880"], "Output", + CellChangeTimes->{{3.8073557418182497`*^9, 3.8073557622024517`*^9}, + 3.8073565723543777`*^9, 3.8073566967878428`*^9}, + CellLabel-> + "Out[193]=",ExpressionUUID->"482a00d7-a1a8-40be-82e4-7361552f8d24"] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData["SCIrPT2"], "Input", + CellLabel-> + "In[194]:=",ExpressionUUID->"12da2fba-f4a2-4f36-9b70-0977ea98aa96"], + +Cell[BoxData[ + RowBox[{"{", + RowBox[{ + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.3705474699999911`"}], ",", + RowBox[{"-", "379.85715710001955`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.3212134599999956`"}], ",", + RowBox[{"-", "439.44437710001694`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.27039005000000316`"}], ",", + RowBox[{"-", "502.81320710001864`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.2244133300000044`"}], ",", + RowBox[{"-", "561.4759271000196`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.18783665000000838`"}], ",", + RowBox[{"-", "608.3896670999991`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.16063102000001095`"}], ",", + RowBox[{"-", "644.1888871000003`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.1396871499999861`"}], ",", + RowBox[{"-", "671.0687671000244`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.12424601000000735`"}], ",", + RowBox[{"-", "692.313257100011`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.111580459999999`"}], ",", + RowBox[{"-", "710.1328171000034`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.10003670999998349`"}], ",", + RowBox[{"-", "725.1875671000221`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.08975225999998315`"}], ",", + RowBox[{"-", "739.1511171000218`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.08040832000000364`"}], ",", + RowBox[{"-", "751.9313671000134`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.07162208999997688`"}], ",", + RowBox[{"-", "764.0102371000239`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.06438604999999598`"}], ",", + RowBox[{"-", "775.6955871000173`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.055509260000008`"}], ",", + RowBox[{"-", "786.8946071000096`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.04794706999999221`"}], ",", + RowBox[{"-", "797.3029471000075`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.04126446000000783`"}], ",", + RowBox[{"-", "806.7485871000031`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.03514844000000039`"}], ",", + RowBox[{"-", "814.8360771000114`"}]}], "}"}]}], "}"}]], "Output", + CellChangeTimes->{3.80735596590476*^9, 3.807356572382283*^9, + 3.807356696792811*^9}, + CellLabel-> + "Out[194]=",ExpressionUUID->"a4fce446-6d7d-44bb-8eda-8363e14b342d"] +}, Open ]], + +Cell[BoxData[{ + RowBox[{ + RowBox[{"SCI", "=", + RowBox[{"Table", "[", + RowBox[{ + RowBox[{"{", + RowBox[{ + RowBox[{"dat", "\[LeftDoubleBracket]", + RowBox[{"k", ",", "1"}], "\[RightDoubleBracket]"}], ",", + RowBox[{ + SuperscriptBox["10", "3"], + RowBox[{"(", + RowBox[{ + RowBox[{"dat", "\[LeftDoubleBracket]", + RowBox[{"k", ",", "2"}], "\[RightDoubleBracket]"}], "-", "EHF"}], + ")"}]}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"k", ",", "start", ",", + RowBox[{"Length", "[", "dat", "]"}]}], "}"}]}], "]"}]}], + ";"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"SCIPT2", "=", + RowBox[{"Table", "[", + RowBox[{ + RowBox[{"{", + RowBox[{ + RowBox[{"dat", "\[LeftDoubleBracket]", + RowBox[{"k", ",", "1"}], "\[RightDoubleBracket]"}], ",", + RowBox[{ + SuperscriptBox["10", "3"], + RowBox[{"(", + RowBox[{ + RowBox[{"dat", "\[LeftDoubleBracket]", + RowBox[{"k", ",", "3"}], "\[RightDoubleBracket]"}], "-", "EHF"}], + ")"}]}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"k", ",", "start", ",", + RowBox[{"Length", "[", "dat", "]"}]}], "}"}]}], "]"}]}], + ";"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"SCIrPT2", "=", + RowBox[{"Table", "[", + RowBox[{ + RowBox[{"{", + RowBox[{ + RowBox[{"dat", "\[LeftDoubleBracket]", + RowBox[{"k", ",", "1"}], "\[RightDoubleBracket]"}], ",", + RowBox[{ + SuperscriptBox["10", "3"], + RowBox[{"(", + RowBox[{ + RowBox[{"dat", "\[LeftDoubleBracket]", + RowBox[{"k", ",", "5"}], "\[RightDoubleBracket]"}], "-", "EHF"}], + ")"}]}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"k", ",", "start", ",", + RowBox[{"Length", "[", "dat", "]"}]}], "}"}]}], "]"}]}], + ";"}]}], "Input", + CellLabel-> + "In[195]:=",ExpressionUUID->"3c63494d-0a28-4d52-b82f-5fa01ee388d7"], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"Table", "[", + RowBox[{ + RowBox[{"{", + RowBox[{ + RowBox[{"dat", "\[LeftDoubleBracket]", + RowBox[{"k", ",", "1"}], "\[RightDoubleBracket]"}], ",", + RowBox[{ + SuperscriptBox["10", "3"], + RowBox[{"(", + RowBox[{ + RowBox[{"dat", "\[LeftDoubleBracket]", + RowBox[{"k", ",", "2"}], "\[RightDoubleBracket]"}], "-", "EHF"}], + ")"}]}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"k", ",", "start", ",", + RowBox[{"Length", "[", "dat", "]"}]}], "}"}]}], "]"}]], "Input", + CellChangeTimes->{3.8073561412146273`*^9}, + CellLabel-> + "In[198]:=",ExpressionUUID->"bb7c7fa2-e87a-461f-bf8e-2d5b920e7e5f"], + +Cell[BoxData[ + RowBox[{"{", + RowBox[{ + RowBox[{"{", + RowBox[{"1280", ",", + RowBox[{"-", "379.85715710001955`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"2560", ",", + RowBox[{"-", "439.44437710001694`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"5120", ",", + RowBox[{"-", "502.81320710001864`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"10240", ",", + RowBox[{"-", "561.4759271000196`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"20480", ",", + RowBox[{"-", "608.3896670999991`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"40960", ",", + RowBox[{"-", "644.1888871000003`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"81920", ",", + RowBox[{"-", "671.0687671000244`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"163840", ",", + RowBox[{"-", "692.313257100011`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"327680", ",", + RowBox[{"-", "710.1328171000034`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"655360", ",", + RowBox[{"-", "725.1875671000221`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"1310720", ",", + RowBox[{"-", "739.1511171000218`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"2621440", ",", + RowBox[{"-", "751.9313671000134`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"5242880", ",", + RowBox[{"-", "764.0102371000239`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"10485760", ",", + RowBox[{"-", "775.6955871000173`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"20971520", ",", + RowBox[{"-", "786.8946071000096`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"41943040", ",", + RowBox[{"-", "797.3029471000075`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"83886080", ",", + RowBox[{"-", "806.7485871000031`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"167772160", ",", + RowBox[{"-", "814.8360771000114`"}]}], "}"}]}], "}"}]], "Output", + CellChangeTimes->{3.807356040568919*^9, 3.807356142565374*^9, + 3.807356572432342*^9, 3.807356696818439*^9}, + CellLabel-> + "Out[198]=",ExpressionUUID->"d086efb9-2f73-4004-98d7-80886efc5243"] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"863.3581802152587`", "-", "850.1618171000018`"}]], "Input", + CellChangeTimes->{{3.807356077062683*^9, 3.8073560772001667`*^9}}, + CellLabel-> + "In[199]:=",ExpressionUUID->"fb7d5558-fd2e-4127-9a0d-1eb25ec9ac64"], + +Cell[BoxData["13.196363115256872`"], "Output", + CellChangeTimes->{3.8073560782407312`*^9, 3.807356572438064*^9, + 3.80735669682316*^9}, + CellLabel-> + "Out[199]=",ExpressionUUID->"315dd49f-4ced-4ec2-93c8-e1f975c23fac"] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData["LinearModelFitError"], "Input", + CellChangeTimes->{{3.807356189470543*^9, 3.807356216543808*^9}}, + CellLabel-> + "In[200]:=",ExpressionUUID->"58dbc827-bfeb-4b41-bd51-761724b81a5c"], + +Cell[BoxData["LinearModelFitError"], "Output", + CellChangeTimes->{3.807356572466381*^9, 3.8073566968277807`*^9}, + CellLabel-> + "Out[200]=",ExpressionUUID->"6e064dd9-efc0-4032-85ea-c45452551709"] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData["LinearModelFitError"], "Input", + CellLabel-> + "In[201]:=",ExpressionUUID->"b21de86e-c1f5-48c5-8892-7a9be82adb1b"], + +Cell[BoxData["LinearModelFitError"], "Output", + CellChangeTimes->{3.807356572471368*^9, 3.8073566968328123`*^9}, + CellLabel-> + "Out[201]=",ExpressionUUID->"05abd436-38fc-4934-8a44-9b687c70ad4a"] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[{ + RowBox[{"fitPT2", "=", + RowBox[{"Fit", "[", + RowBox[{ + RowBox[{"Take", "[", + RowBox[{"SCIPT2", ",", + RowBox[{"-", "nfit"}]}], "]"}], ",", + RowBox[{"{", + RowBox[{"1", ",", "x"}], "}"}], ",", "x"}], + "]"}]}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"SCIrPT2", "=", + RowBox[{"Table", "[", + RowBox[{ + RowBox[{"{", + RowBox[{ + RowBox[{ + RowBox[{"dat", "\[LeftDoubleBracket]", + RowBox[{"k", ",", "5"}], "\[RightDoubleBracket]"}], "-", + RowBox[{"dat", "\[LeftDoubleBracket]", + RowBox[{"k", ",", "2"}], "\[RightDoubleBracket]"}]}], ",", + RowBox[{ + SuperscriptBox["10", "3"], + RowBox[{"(", + RowBox[{ + RowBox[{"dat", "\[LeftDoubleBracket]", + RowBox[{"k", ",", "2"}], "\[RightDoubleBracket]"}], "-", "EHF"}], + ")"}]}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"k", ",", "start", ",", + RowBox[{"Length", "[", "dat", "]"}]}], "}"}]}], "]"}]}], + ";"}], "\[IndentingNewLine]", + RowBox[{"fitrPT2", "=", + RowBox[{"Fit", "[", + RowBox[{ + RowBox[{"Take", "[", + RowBox[{"SCIrPT2", ",", + RowBox[{"-", "nfit"}]}], "]"}], ",", + RowBox[{"{", + RowBox[{"1", ",", "x"}], "}"}], ",", "x"}], "]"}]}]}], "Input", + CellLabel-> + "In[202]:=",ExpressionUUID->"7325d99b-25ff-4080-9c55-fe88f51bc7fa"], + +Cell[BoxData[ + RowBox[{ + RowBox[{"-", "843.1397249260968`"}], "-", + RowBox[{"4.563321652616199`*^-8", " ", "x"}]}]], "Output", + CellChangeTimes->{3.807356572507834*^9, 3.8073566968513002`*^9}, + CellLabel-> + "Out[202]=",ExpressionUUID->"4f4d5a30-cc54-48f6-895e-18717701761f"], + +Cell[BoxData[ + RowBox[{ + RowBox[{"-", "863.3581802152587`"}], "-", + RowBox[{"1376.8363964252949`", " ", "x"}]}]], "Output", + CellChangeTimes->{3.807356572507834*^9, 3.80735669685188*^9}, + CellLabel-> + "Out[204]=",ExpressionUUID->"7c31df32-8329-427d-986b-c782115d01e1"] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{ + RowBox[{"f", "=", + RowBox[{"LinearModelFit", "[", + RowBox[{ + RowBox[{"Take", "[", + RowBox[{"SCIrPT2", ",", + RowBox[{"-", "nfit"}]}], "]"}], ",", + RowBox[{"{", + RowBox[{"1", ",", "x"}], "}"}], ",", "x"}], "]"}]}], + "\[IndentingNewLine]"}]], "Input", + CellChangeTimes->{{3.807356266704557*^9, 3.807356273244464*^9}, { + 3.807356330208477*^9, 3.807356339740362*^9}}, + CellLabel-> + "In[205]:=",ExpressionUUID->"5d82dd32-c3b1-4401-93b3-6a427f759fd0"], + +Cell[BoxData[ + TagBox[ + RowBox[{"FittedModel", "[", + TagBox[ + PanelBox[ + TagBox[ + RowBox[{ + RowBox[{"-", "863.3581802152567`"}], "-", + RowBox[{"1376.8363964252455`", " ", "x"}]}], + Short[#, 2]& ], + FrameMargins->5], + Editable -> False], "]"}], + InterpretTemplate[ + FittedModel[{ + "Linear", {-863.3581802152567, -1376.8363964252455`}, {{$CellContext`x}, { + 1, $CellContext`x}}, {0, 0}}, {{1., 1., 1., + 1.}}, {{-0.055509260000008, -786.8946071000096}, {-0.04794706999999221, \ +-797.3029471000075}, {-0.04126446000000783, -806.7485871000031}, \ +{-0.03514844000000039, -814.8360771000114}}, {{1., -0.055509260000008}, { + 1., -0.04794706999999221}, {1., -0.04126446000000783}, { + 1., -0.03514844000000039}}, + Function[Null, + Internal`LocalizedBlock[{$CellContext`x}, #], {HoldAll}]]& ], + Editable->False, + SelectWithContents->True, + Selectable->True]], "Output", + CellChangeTimes->{3.807356273518469*^9, 3.8073563401447144`*^9, + 3.807356414213882*^9, 3.807356572516573*^9, 3.8073566968566933`*^9}, + CellLabel-> + "Out[205]=",ExpressionUUID->"24dbd7fd-4a7f-405f-bab3-00b41d675f8f"] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"f", "[", "\"\\"", "]"}]], "Input", + CellChangeTimes->{{3.807356449044463*^9, 3.807356450011277*^9}, { + 3.807356490396319*^9, 3.8073564913631067`*^9}}, + CellLabel-> + "In[206]:=",ExpressionUUID->"11276897-ac60-4098-bfe3-804e94f6b5a8"], + +Cell[BoxData[ + StyleBox[ + TagBox[GridBox[{ + {"\<\"\"\>", "\<\"Estimate\"\>", "\<\"Standard Error\"\>", "\<\"t\ +\[Hyphen]Statistic\"\>", "\<\"P\[Hyphen]Value\"\>"}, + {"1", + RowBox[{"-", "863.3581802152567`"}], "0.526587143072088`", + RowBox[{"-", "1639.5352442113572`"}], "3.720131089846035`*^-7"}, + {"x", + RowBox[{"-", "1376.8363964252455`"}], "11.54731952566032`", + RowBox[{"-", "119.23428578949907`"}], "0.00007033182256310777`"} + }, + AutoDelete->False, + GridBoxAlignment->{"Columns" -> {{Left}}, "Rows" -> {{Automatic}}}, + GridBoxDividers->{ + "ColumnsIndexed" -> {2 -> GrayLevel[0.7]}, + "RowsIndexed" -> {2 -> GrayLevel[0.7]}}, + GridBoxItemSize->{"Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, + GridBoxSpacings->{ + "ColumnsIndexed" -> {2 -> 1}, "RowsIndexed" -> {2 -> 0.75}}], + "Grid"], "DialogStyle", + StripOnInput->False]], "Output", + CellChangeTimes->{3.807356450508697*^9, 3.807356491593397*^9, + 3.807356572554665*^9, 3.8073566968910007`*^9}, + CellLabel-> + "Out[206]=",ExpressionUUID->"1d05d124-17d0-4b64-8237-0ac9d00311ee"] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"f", "[", "\"\\"", "]"}]], "Input", + CellChangeTimes->{{3.807356402451572*^9, 3.807356425645406*^9}}, + CellLabel-> + "In[207]:=",ExpressionUUID->"b22ca888-a8e8-4120-b6e9-9321f209da71"], + +Cell[BoxData[ + RowBox[{"{", + RowBox[{"0.526587143072088`", ",", "11.54731952566032`"}], "}"}]], "Output",\ + + CellChangeTimes->{{3.807356402623375*^9, 3.80735642591369*^9}, + 3.80735657256141*^9, 3.8073566968955183`*^9}, + CellLabel-> + "Out[207]=",ExpressionUUID->"f93e1c9b-3d6f-4087-bd6f-29699c9090c6"] }, Open ]] }, Open ]] }, Open ]] }, -WindowSize->{1280, 755}, +WindowSize->{1184, 1394}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, -FrontEndVersion->"12.1 for Mac OS X x86 (64-bit) (March 13, 2020)", +FrontEndVersion->"12.1 for Mac OS X x86 (64-bit) (June 19, 2020)", StyleDefinitions->"Default.nb", ExpressionUUID->"ad7ac4fb-d7dd-4a37-9960-87ba52ede3bd" ] @@ -17932,34 +18232,76 @@ Cell[CellGroupData[{ Cell[1738, 54, 150, 3, 72, "Title",ExpressionUUID->"b01a9543-ed00-4d3e-bf93-f93e73170b6d"], Cell[CellGroupData[{ Cell[1913, 61, 166, 3, 54, "Subsection",ExpressionUUID->"0d117b7e-62a1-424a-bf47-2d071fb26070"], -Cell[2082, 66, 825, 18, 52, "Input",ExpressionUUID->"fab9ee47-3e93-4201-bb7f-09ae5e949598"], +Cell[2082, 66, 821, 17, 52, "Input",ExpressionUUID->"fab9ee47-3e93-4201-bb7f-09ae5e949598"], Cell[CellGroupData[{ -Cell[2932, 88, 8487, 195, 780, "Input",ExpressionUUID->"070c4b28-8c43-4bc2-ad60-9265a5b4fc44"], -Cell[11422, 285, 255963, 4802, 584, "Output",ExpressionUUID->"d0d8ae32-ad10-4dcf-9f77-479e69e9ec0e"] +Cell[2928, 87, 8483, 194, 780, "Input",ExpressionUUID->"070c4b28-8c43-4bc2-ad60-9265a5b4fc44"], +Cell[11414, 283, 253634, 4741, 584, "Output",ExpressionUUID->"3dc7aee2-dfb7-4266-8af9-123c83a80027"] }, Open ]], Cell[CellGroupData[{ -Cell[267422, 5092, 10060, 226, 875, "Input",ExpressionUUID->"abf97065-2c12-498b-b03e-d7901e4cc68f"], -Cell[277485, 5320, 1503, 23, 34, "Output",ExpressionUUID->"9bdab6ac-e2d5-4ce8-a951-571c9f0a6425"], -Cell[278991, 5345, 1501, 23, 34, "Output",ExpressionUUID->"ac0955a2-ef61-42bb-abf3-f760a5c7a9b4"], -Cell[280495, 5370, 197256, 3604, 565, "Output",ExpressionUUID->"883fce0c-13d8-4e46-bb43-ba9dd7d6f602"] +Cell[265085, 5029, 10060, 226, 875, "Input",ExpressionUUID->"abf97065-2c12-498b-b03e-d7901e4cc68f"], +Cell[275148, 5257, 1503, 23, 34, "Output",ExpressionUUID->"9bdab6ac-e2d5-4ce8-a951-571c9f0a6425"], +Cell[276654, 5282, 1501, 23, 34, "Output",ExpressionUUID->"ac0955a2-ef61-42bb-abf3-f760a5c7a9b4"], +Cell[278158, 5307, 197256, 3604, 565, "Output",ExpressionUUID->"883fce0c-13d8-4e46-bb43-ba9dd7d6f602"] }, Open ]] }, Closed]], Cell[CellGroupData[{ -Cell[477800, 8980, 168, 3, 38, "Subsection",ExpressionUUID->"32a5cbc6-7ca0-469d-b531-3ef832f9f00b"], -Cell[477971, 8985, 782, 17, 52, "Input",ExpressionUUID->"88b83480-ff85-479c-9bd5-86e50cd48da6"], +Cell[475463, 8917, 168, 3, 38, "Subsection",ExpressionUUID->"32a5cbc6-7ca0-469d-b531-3ef832f9f00b"], +Cell[475634, 8922, 833, 18, 52, "Input",ExpressionUUID->"88b83480-ff85-479c-9bd5-86e50cd48da6"], Cell[CellGroupData[{ -Cell[478778, 9006, 8583, 198, 801, "Input",ExpressionUUID->"f0f1eea9-20b9-48df-9572-c510784075b4"], -Cell[487364, 9206, 256135, 4806, 584, "Output",ExpressionUUID->"0de2ef17-1213-4011-ac72-c6e9e0caf6db"] +Cell[476492, 8944, 8583, 198, 801, "Input",ExpressionUUID->"f0f1eea9-20b9-48df-9572-c510784075b4"], +Cell[485078, 9144, 253860, 4746, 584, "Output",ExpressionUUID->"f08a5516-ecb9-4b9d-b128-154a3ac2cb7b"] }, Open ]], Cell[CellGroupData[{ -Cell[743536, 14017, 8088, 205, 896, "Input",ExpressionUUID->"0daea3d8-8f63-4ebd-ba0f-4b294d15916e"], -Cell[751627, 14224, 1472, 23, 34, "Output",ExpressionUUID->"7a77a453-90c2-4388-b095-5d7536c2ebd6"], -Cell[753102, 14249, 1473, 23, 34, "Output",ExpressionUUID->"3d50c171-50af-40d7-83b3-f304118a6b8b"], -Cell[754578, 14274, 197439, 3607, 565, "Output",ExpressionUUID->"5b22d62c-c18f-4d0b-b5f0-2d98e0faf6a0"] +Cell[738975, 13895, 8100, 205, 896, "Input",ExpressionUUID->"0daea3d8-8f63-4ebd-ba0f-4b294d15916e"], +Cell[747078, 14102, 1618, 25, 70, "Output",ExpressionUUID->"70b3141b-5f73-46ce-a120-e28e26a497f7"], +Cell[748699, 14129, 1619, 25, 70, "Output",ExpressionUUID->"b622dec9-d93a-44a3-973b-53c6aa852652"], +Cell[750321, 14156, 195942, 3568, 70, "Output",ExpressionUUID->"8d3b739a-d028-4230-b557-d5db482643b2"] }, Open ]], Cell[CellGroupData[{ -Cell[952054, 17886, 420, 9, 30, "Input",ExpressionUUID->"ba14da0a-87ae-46da-85c2-b124d0f9c105"], -Cell[952477, 17897, 205, 4, 34, "Output",ExpressionUUID->"a9ea879d-0902-454b-a505-60f0edd2cffd"] +Cell[946300, 17729, 420, 9, 30, "Input",ExpressionUUID->"ba14da0a-87ae-46da-85c2-b124d0f9c105"], +Cell[946723, 17740, 278, 5, 70, "Output",ExpressionUUID->"172c6c74-0ddd-4b88-801f-ed238e92e91d"] +}, Open ]], +Cell[CellGroupData[{ +Cell[947038, 17750, 212, 4, 30, "Input",ExpressionUUID->"3e45b97c-e088-4c82-99ef-15e7d70333ed"], +Cell[947253, 17756, 238, 4, 70, "Output",ExpressionUUID->"482a00d7-a1a8-40be-82e4-7361552f8d24"] +}, Open ]], +Cell[CellGroupData[{ +Cell[947528, 17765, 116, 2, 30, "Input",ExpressionUUID->"12da2fba-f4a2-4f36-9b70-0977ea98aa96"], +Cell[947647, 17769, 2684, 78, 70, "Output",ExpressionUUID->"a4fce446-6d7d-44bb-8eda-8363e14b342d"] +}, Open ]], +Cell[950346, 17850, 1975, 59, 79, "Input",ExpressionUUID->"3c63494d-0a28-4d52-b82f-5fa01ee388d7"], +Cell[CellGroupData[{ +Cell[952346, 17913, 677, 19, 33, "Input",ExpressionUUID->"bb7c7fa2-e87a-461f-bf8e-2d5b920e7e5f"], +Cell[953026, 17934, 2089, 60, 70, "Output",ExpressionUUID->"d086efb9-2f73-4004-98d7-80886efc5243"] +}, Open ]], +Cell[CellGroupData[{ +Cell[955152, 17999, 234, 4, 30, "Input",ExpressionUUID->"fb7d5558-fd2e-4127-9a0d-1eb25ec9ac64"], +Cell[955389, 18005, 219, 4, 70, "Output",ExpressionUUID->"315dd49f-4ced-4ec2-93c8-e1f975c23fac"] +}, Open ]], +Cell[CellGroupData[{ +Cell[955645, 18014, 194, 3, 30, "Input",ExpressionUUID->"58dbc827-bfeb-4b41-bd51-761724b81a5c"], +Cell[955842, 18019, 195, 3, 70, "Output",ExpressionUUID->"6e064dd9-efc0-4032-85ea-c45452551709"] +}, Open ]], +Cell[CellGroupData[{ +Cell[956074, 18027, 128, 2, 30, "Input",ExpressionUUID->"b21de86e-c1f5-48c5-8892-7a9be82adb1b"], +Cell[956205, 18031, 195, 3, 70, "Output",ExpressionUUID->"05abd436-38fc-4934-8a44-9b687c70ad4a"] +}, Open ]], +Cell[CellGroupData[{ +Cell[956437, 18039, 1375, 41, 75, "Input",ExpressionUUID->"7325d99b-25ff-4080-9c55-fe88f51bc7fa"], +Cell[957815, 18082, 278, 6, 70, "Output",ExpressionUUID->"4f4d5a30-cc54-48f6-895e-18717701761f"], +Cell[958096, 18090, 272, 6, 70, "Output",ExpressionUUID->"7c31df32-8329-427d-986b-c782115d01e1"] +}, Open ]], +Cell[CellGroupData[{ +Cell[958405, 18101, 509, 14, 52, "Input",ExpressionUUID->"5d82dd32-c3b1-4401-93b3-6a427f759fd0"], +Cell[958917, 18117, 1150, 29, 70, "Output",ExpressionUUID->"24dbd7fd-4a7f-405f-bab3-00b41d675f8f"] +}, Open ]], +Cell[CellGroupData[{ +Cell[960104, 18151, 275, 5, 30, "Input",ExpressionUUID->"11276897-ac60-4098-bfe3-804e94f6b5a8"], +Cell[960382, 18158, 1112, 25, 70, "Output",ExpressionUUID->"1d05d124-17d0-4b64-8237-0ac9d00311ee"] +}, Open ]], +Cell[CellGroupData[{ +Cell[961531, 18188, 225, 4, 30, "Input",ExpressionUUID->"b22ca888-a8e8-4120-b6e9-9321f209da71"], +Cell[961759, 18194, 306, 7, 34, "Output",ExpressionUUID->"f93e1c9b-3d6f-4087-bd6f-29699c9090c6"] }, Open ]] }, Open ]] }, Open ]] diff --git a/benzene.tex b/benzene.tex index 1830785..34375f7 100644 --- a/benzene.tex +++ b/benzene.tex @@ -43,7 +43,7 @@ % Abstract \begin{abstract} -Following the recent work of Eriksen \textit{et al.}~[\href{https://arxiv.org/abs/2008.02678}{arXiv:2008.02678 [physics.chem-ph]}], we report the performance of the \textit{Configuration Interaction using a Perturbative Selection made Iteratively} (CIPSI) method on the non-relativistic frozen-core correlation energy of the benzene molecule in the cc-pVDZ basis. Following our usual protocol, we obtain a correlation energy of \titou{$-86x.x(x)$} m$E_h$ which agrees with the theoretical estimate of $-863$ m$E_h$ proposed by Eriksen \textit{et al.}~using an extensive array of highly-accurate new electronic structure methods. +Following the recent work of Eriksen \textit{et al.}~[\href{https://arxiv.org/abs/2008.02678}{arXiv:2008.02678 [physics.chem-ph]}], we report the performance of the \textit{Configuration Interaction using a Perturbative Selection made Iteratively} (CIPSI) method on the non-relativistic frozen-core correlation energy of the benzene molecule in the cc-pVDZ basis. Following our usual protocol, we obtain a correlation energy of $-863.4(5)$ m$E_h$ which agrees with the theoretical estimate of $-863$ m$E_h$ proposed by Eriksen \textit{et al.}~using an extensive array of highly-accurate new electronic structure methods. \end{abstract} % Title @@ -97,7 +97,7 @@ The outcome of this work is nicely summarized in the abstract of Ref.~\onlinecit \hline ph-AFQMC & $-864.3(4)$ & \onlinecite{Lee_2020} \\ \hline - CIPSI & \titou{$-86x.x(x)$} & This work \\ + CIPSI & $-863.4(5)$ & This work \\ \end{tabular} \end{ruledtabular} \end{table} @@ -142,25 +142,14 @@ higher-lying $\sigma$ [37,38,40,44,45,47,48,51,52,58--70,75--81,86,88--91,94--97 Like Pipek-Mezey, \cite{Pipek_1989} this choice of orbital windows allows to preserve a strict $\sigma$-$\pi$ separation in planar systems like benzene. As one can see from the energies of Table \ref{tab:NOvsLO}, for a given value of $N_\text{det}$, the variational energy as well as the PT2-corrected energies are much lower with localized orbitals than with natural orbitals. We, therefore, consider these energies more trustworthy, and we will base our best estimate of the correlation energy of benzene on these calculations. The convergence of the CIPSI correlation energy using localized orbitals is illustrated in Fig.~\ref{fig:CIPSI}, where one can see the behavior of the correlation energy, $\Delta E_\text{var.}$ and $\Delta E_\text{var.} + E_\text{(r)PT2}$, as a function of $N_\text{det}$ (left panel). -The right panel of Fig.~\ref{fig:CIPSI} is more instructive as it shows $\Delta E_\text{var.}$ as a function of $E_\text{(r)PT2}$, and their corresponding \titou{five}-point linear extrapolation curves that we have used to get our final estimate of the correlation energy. +The right panel of Fig.~\ref{fig:CIPSI} is more instructive as it shows $\Delta E_\text{var.}$ as a function of $E_\text{(r)PT2}$, and their corresponding four-point linear extrapolation curves that we have used to get our final estimate of the correlation energy. From this figure, one clearly sees that the rPT2-based correction behaves more linearly than its corresponding PT2 version, and is thus systematically employed in the following. % Results Our final number are gathered in Table \ref{tab:extrap_dist_table}, where, following the notations of Ref.~\onlinecite{Eriksen_2020}, we report, in addition to the final variational energies $\Delta E_{\text{var.}}$, the extrapolation distances, $\Delta E_{\text{dist}}$, defined as the difference between the final computed energy, $\Delta E_{\text{final}}$, and the extrapolated energy, $\Delta E_{\text{extrap.}}$ associated with ASCI, iCI, SHCI, DMRS, and CIPSI. -The three flavours of SCI fall into an interval ranging from $-860.0$ m$E_h$ (ASCI) to $-864.2$ m$E_h$ (SHCI), while the other non-SCI methods yield correlation energies ranging from $-863.7$ to $-862.8$ m$E_h$ (see Table \ref{tab:energy}). Our final CIPSI number (obtained with localized orbitals and rPT2 correction via a \titou{five}-point linear extrapolation) is \titou{$-86x.x(x)$} m$E_h$. +The three flavours of SCI fall into an interval ranging from $-860.0$ m$E_h$ (ASCI) to $-864.2$ m$E_h$ (SHCI), while the other non-SCI methods yield correlation energies ranging from $-863.7$ to $-862.8$ m$E_h$ (see Table \ref{tab:energy}). Our final CIPSI number (obtained with localized orbitals and rPT2 correction via a four-point linear extrapolation) is $-863.4(5)$ m$E_h$, where the error reported in parenthesis represents the fitting error (not the extrapolation error for which it is much harder to provide a theoretically sound estimate). For comparison, the best post blind test SHCI estimate is $-863.3$ m$E_h$, which agrees almost perfectly with our best CIPSI estimate. -\titou{Although it is not possible to provide a theoretically sound error bar, we estimate the extrapolation error by ... -We believe that it provides a very safe estimate of the extrapolation error.} - -% Timings -The present calculations have been performed on the AMD partition of GENCI's Irene supercomputer. -Each Irene's AMD node is a dual-socket AMD Rome (Epyc) CPU@2.60 GHz with 256GiB of RAM, with a total of 64 physical CPU cores per socket. -These nodes are connected via Infiniband HDR100. -The first step of the calculation, \ie, performing a CIPSI calculation up to $N_\text{det} \sim 10^7$ with Hartree-Fock orbitals in order to produce natural orbitals, takes roughly 24 hours on a single node, and reaching the same number of determinants with natural orbitals or localized orbitals takes roughly the same amount of time. -A second 24-hour run on 10 distributed nodes was performed to push the selection to $8 \times 10^7$ determinants, and a third distributed run using 40 nodes was used to reach $16 \times 10^7$ determinants. -A last 24-hour run on 60 distributed nodes was performed to obtain the largest wave function of $25 \times 10^7$ determinants and its PT2 correction. -In total, the present calculation has required \titou{340k} core hours, almost half of it being spent in the last stage of the computation. %%$ FIG. 1 %%% \begin{figure*} @@ -171,7 +160,7 @@ In total, the present calculation has required \titou{340k} core hours, almost h Convergence of the CIPSI correlation energy of benzene using localized orbitals. Left: $\Delta E_\text{var.}$, $\Delta E_\text{var.} + E_\text{PT2}$, and $\Delta E_\text{var.} + E_\text{rPT2}$ (in m$E_h$) as functions of the number of determinants in the variational space $N_\text{det}$. Right: $\Delta E_\text{var.}$ (in m$E_h$) as a function of $E_\text{PT2}$ or $E_\text{rPT2}$. - The \titou{five}-point linear extrapolation curves (dashed lines) are also reported. + The four-point linear extrapolation curves (dashed lines) are also reported. The theoretical estimate of $-863$ m$E_h$ from Ref.~\onlinecite{Eriksen_2020} is marked by a black line for comparison purposes. \label{fig:CIPSI} } @@ -216,7 +205,6 @@ The statistical error on $E_\text{(r)PT2}$, corresponding to one standard deviat 41\,943\,040 & $-231.487\,978$ & $-231.564\,529(382)$ & $-231.563\,593(377)$ & $-231.519\,122$ & $-231.567\,419(240)$ & $-231.567\,069(238)$ \\ 83\,886\,080 & $-231.501\,334$ & $-231.566\,994(317)$ & $-231.566\,325(314)$ & $-231.528\,568$ & $-231.570\,084(199)$ & $-231.569\,832(198)$ \\ 167\,772\,160 & $-231.514\,009$ & $-231.569\,939(273)$ & $-231.569\,467(271)$ & $-231.536\,655$ & $-231.571\,981(175)$ & $-231.571\,804(174)$ \\ - 251\,658\,240 & \\ \end{tabular} \end{ruledtabular} \end{table*} @@ -240,11 +228,18 @@ The statistical error on $E_\text{(r)PT2}$, corresponding to one standard deviat SHCI & $-827.2$ & $-852.8$ & $-864.2$ & $-11.4$ \\ DMRG & $-859.2$ & $-859.2$ & $-862.8$ & $-3.6$ \\ \hline - CIPSI & \titou{$-8xx.x$} & \titou{$-8xx.x$} & \titou{$-86x.x$} & \titou{$-xx.x$} \\ + CIPSI & $-814.8$ & $-850.2$ & $-863.4$ & $-13.2$ \\ \end{tabular} \end{ruledtabular} \end{table} +% Timings +The present calculations have been performed on the AMD partition of GENCI's Irene supercomputer. +Each Irene's AMD node is a dual-socket AMD Rome (Epyc) CPU@2.60 GHz with 256GiB of RAM, with a total of 64 physical CPU cores per socket. +These nodes are connected via Infiniband HDR100. +The first step of the calculation, \ie, performing a CIPSI calculation up to $N_\text{det} \sim 10^7$ with Hartree-Fock orbitals in order to produce natural orbitals, takes roughly 24 hours on a single node, and reaching the same number of determinants with natural orbitals or localized orbitals takes roughly the same amount of time. +A second 24-hour run on 10 distributed nodes was performed to push the selection to $8 \times 10^7$ determinants, and a third distributed run using 40 nodes was used to reach $16 \times 10^7$ determinants. +In total, the present calculation has required 150k core hours, most of it being spent in the last stage of the computation. % Acknowledgements diff --git a/fig1a.pdf b/fig1a.pdf index 0608819..4dd4454 100644 Binary files a/fig1a.pdf and b/fig1a.pdf differ diff --git a/fig1b.pdf b/fig1b.pdf index 506610b..769acea 100644 Binary files a/fig1b.pdf and b/fig1b.pdf differ