minor corrections at proof stage

This commit is contained in:
Pierre-Francois Loos 2020-10-20 21:30:07 +02:00
parent ad8d048848
commit 6d76db7fe4
3 changed files with 174 additions and 78 deletions

View File

@ -10,10 +10,10 @@
NotebookFileLineBreakTest NotebookFileLineBreakTest
NotebookFileLineBreakTest NotebookFileLineBreakTest
NotebookDataPosition[ 158, 7] NotebookDataPosition[ 158, 7]
NotebookDataLength[ 961285, 17989] NotebookDataLength[ 965485, 18085]
NotebookOptionsPosition[ 958752, 17938] NotebookOptionsPosition[ 962949, 18034]
NotebookOutlinePosition[ 959145, 17954] NotebookOutlinePosition[ 963343, 18050]
CellTagsIndexPosition[ 959102, 17951] CellTagsIndexPosition[ 963300, 18047]
WindowFrame->Normal*) WindowFrame->Normal*)
(* Beginning of Notebook Content *) (* Beginning of Notebook Content *)
@ -8905,8 +8905,7 @@ Cell[BoxData[{
3.807261814034483*^9, 3.807261814866003*^9}, {3.807356542943404*^9, 3.807261814034483*^9, 3.807261814866003*^9}, {3.807356542943404*^9,
3.8073565589786997`*^9}, {3.8112324241990013`*^9, 3.811232425322535*^9}, { 3.8073565589786997`*^9}, {3.8112324241990013`*^9, 3.811232425322535*^9}, {
3.811232472089983*^9, 3.811232472146283*^9}}, 3.811232472089983*^9, 3.811232472146283*^9}},
CellLabel-> CellLabel->"In[4]:=",ExpressionUUID->"88b83480-ff85-479c-9bd5-86e50cd48da6"],
"In[419]:=",ExpressionUUID->"88b83480-ff85-479c-9bd5-86e50cd48da6"],
Cell[CellGroupData[{ Cell[CellGroupData[{
@ -9107,8 +9106,7 @@ Cell[BoxData[{
3.807029983522843*^9, 3.807030011586466*^9}, {3.80715444508188*^9, 3.807029983522843*^9, 3.807030011586466*^9}, {3.80715444508188*^9,
3.807154458718595*^9}, {3.807240495668721*^9, 3.807240496242681*^9}, 3.807154458718595*^9}, {3.807240495668721*^9, 3.807240496242681*^9},
3.807263260480081*^9, {3.811232475114682*^9, 3.8112324751735573`*^9}}, 3.807263260480081*^9, {3.811232475114682*^9, 3.8112324751735573`*^9}},
CellLabel-> CellLabel->"In[6]:=",ExpressionUUID->"f0f1eea9-20b9-48df-9572-c510784075b4"],
"In[421]:=",ExpressionUUID->"f0f1eea9-20b9-48df-9572-c510784075b4"],
Cell[BoxData[ Cell[BoxData[
TagBox[ TagBox[
@ -9317,7 +9315,26 @@ Cell[BoxData[
16.8586758847529, -842.4038671000176}}, {{ 16.8586758847529, -842.4038671000176}}, {{
17.551823065312842`, -845.2500170999997}}, {{ 17.551823065312842`, -845.2500170999997}}, {{
18.24497024587279, -848.0130471000109}}, {{ 18.24497024587279, -848.0130471000109}}, {{
18.938117426432733`, -849.9845171000118}}}]}}}, {{}, {}}}, InsetBox[ 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[
TemplateBox[{ TemplateBox[{
GraphicsBox[{ GraphicsBox[{
Thickness[0.00970591089973794], Thickness[0.00970591089973794],
@ -10720,7 +10737,6 @@ TJDW
12.5797}, {65.07969999999999, 12.674999999999999`}, { 12.5797}, {65.07969999999999, 12.674999999999999`}, {
65.07969999999999, 14.0828}, {66.0578, 14.940599999999998`}, { 65.07969999999999, 14.0828}, {66.0578, 14.940599999999998`}, {
68.9422, 15.9891}, {68.9422, 12.532799999999998`}}}], 68.9422, 15.9891}, {68.9422, 12.532799999999998`}}}],
FilledCurveBox[{{{0, 2, 0}, {0, 1, 0}, {1, 3, 3}, {1, 3, 3}, {1, 3, 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}, {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, 3}, {0, 1, 0}, {1, 3, 3}, {1, 3, 3}, {0, 1, 0}, {0, 1, 0}, {0, 1,
@ -11212,7 +11228,6 @@ dDelIvjuD6tE1s2PdSg/vM115lpluPth8QXjl4Lk7yo5oKcXAJ2oB70=
19.900000000000002`}, {119.30599999999998`, 18.7781}, { 19.900000000000002`}, {119.30599999999998`, 18.7781}, {
119.30599999999998`, 17.515600000000003`}, {119.30599999999998`, 119.30599999999998`, 17.515600000000003`}, {119.30599999999998`,
12.0313}}}], 12.0313}}}],
FilledCurveBox[{{{1, 4, 3}, {1, 3, 3}, {1, 3, 3}, {0, 1, 0}, {1, 3, 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, 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[" 1, 0}}, {{1, 4, 3}, {1, 3, 3}, {0, 1, 0}}}, {CompressedData["
@ -11443,6 +11458,7 @@ aPotByJer+eg0y5289ziHIfrQp8czx8zcAi3BAZAX44Dr//6Kakahg6xwOzC
Automatic, Automatic}, {{{6.907755278982137, Automatic, Automatic}, {{{6.907755278982137,
FormBox["1000", TraditionalForm], {0.01, 0.}}, {9.210340371976184, FormBox["1000", TraditionalForm], {0.01, 0.}}, {9.210340371976184,
FormBox[ FormBox[
TemplateBox[{"10", "4"}, "Superscript", SyntaxForm -> TemplateBox[{"10", "4"}, "Superscript", SyntaxForm ->
SuperscriptBox], TraditionalForm], {0.01, 0.}}, { SuperscriptBox], TraditionalForm], {0.01, 0.}}, {
11.512925464970229`, 11.512925464970229`,
@ -12123,8 +12139,52 @@ aPotByJer+eg0y5289ziHIfrQp8czx8zcAi3BAZAX44Dr//6Kakahg6xwOzC
16.8586758847529, -842.4038671000176}}, {{ 16.8586758847529, -842.4038671000176}}, {{
17.551823065312842`, -845.2500170999997}}, {{ 17.551823065312842`, -845.2500170999997}}, {{
18.24497024587279, -848.0130471000109}}, {{ 18.24497024587279, -848.0130471000109}}, {{
18.938117426432733`, -849.9845171000118}}}]}}}, {{}, {}}}, { 18.938117426432733`, -849.9845171000118}}}]}}, {{
DisplayFunction -> Identity, GridLines -> { 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 -> {
Charting`ScaledTickValues[{Log, Exp}], Automatic}, DisplayFunction -> Charting`ScaledTickValues[{Log, Exp}], Automatic}, DisplayFunction ->
Identity, DisplayFunction -> Identity, DisplayFunction -> Identity, Identity, DisplayFunction -> Identity, DisplayFunction -> Identity,
DisplayFunction -> Identity, AspectRatio -> 1, Axes -> {False, False}, DisplayFunction -> Identity, AspectRatio -> 1, Axes -> {False, False},
@ -12228,6 +12288,7 @@ A3p8AgDLMM+X
19.707800000000002`}, {51.7141, 19.040599999999998`}, { 19.707800000000002`}, {51.7141, 19.040599999999998`}, {
52.02340000000001, 16.967200000000002`}, {47.1609, 52.02340000000001, 16.967200000000002`}, {47.1609,
16.967200000000002`}}}], 16.967200000000002`}}}],
FilledCurve[{{{0, 2, 0}, {1, 3, 3}, {1, 3, 3}, {0, 1, 0}, {1, 3, 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, 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[" 1, 0}, {1, 3, 3}, {0, 1, 0}}}, CompressedData["
@ -12322,9 +12383,8 @@ nyDsg/Fh7oHxS7aK/j5dZ+KAHp4AGg2gZQ==
91.2688, 10.0297}, {90.38749999999999, 10.6016}, { 91.2688, 10.0297}, {90.38749999999999, 10.6016}, {
89.88589999999999, 11.5547}, {89.2188, 12.7953}, { 89.88589999999999, 11.5547}, {89.2188, 12.7953}, {
88.83749999999999, 14.464099999999997`}, {88.83749999999999, 88.83749999999999, 14.464099999999997`}, {88.83749999999999,
16.156299999999998`}, {88.83749999999999, 16.156299999999998`}, {88.83749999999999, 18.421899999999994`},
18.421899999999994`}, {89.95780000000002, {89.95780000000002, 19.900000000000002`}, {91.64999999999999,
19.900000000000002`}, {91.64999999999999,
19.900000000000002`}}}], 19.900000000000002`}}}],
FilledCurve[{{{1, 4, 3}, {1, 3, 3}, {0, 1, 0}, {1, 3, 3}, {0, 1, FilledCurve[{{{1, 4, 3}, {1, 3, 3}, {0, 1, 0}, {1, 3, 3}, {0, 1,
@ -13854,9 +13914,9 @@ W1zuAg==
3.807243658126053*^9, 3.8072615834249372`*^9, 3.807261817092359*^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.807355433753173*^9, {3.807356544437235*^9, 3.807356571477729*^9},
3.807356695973937*^9, 3.8112231486510353`*^9, 3.8112325556465607`*^9}, 3.807356695973937*^9, 3.8112231486510353`*^9, 3.8112325556465607`*^9,
CellLabel-> 3.812206375101117*^9},
"Out[430]=",ExpressionUUID->"94181d5f-8832-4765-b54e-766c93d91980"] CellLabel->"Out[15]=",ExpressionUUID->"7d8b90dc-d960-4d55-affb-657fffbd0185"]
}, Open ]], }, Open ]],
Cell[CellGroupData[{ Cell[CellGroupData[{
@ -13904,7 +13964,7 @@ Cell[BoxData[{
RowBox[{ RowBox[{
RowBox[{"{", RowBox[{"{",
RowBox[{ RowBox[{
RowBox[{"-", "1"}], ",", RowBox[{"-", "1000"}], ",",
RowBox[{"-", "863"}]}], "}"}], ",", RowBox[{"-", "863"}]}], "}"}], ",",
RowBox[{"{", RowBox[{"{",
RowBox[{"0", ",", RowBox[{"0", ",",
@ -14101,9 +14161,8 @@ $E_\\\\text{rPT2}$}\>\"", ",",
3.811223364466147*^9}, {3.811223406473289*^9, 3.811223407520194*^9}, { 3.811223364466147*^9}, {3.811223406473289*^9, 3.811223407520194*^9}, {
3.8112234398007393`*^9, 3.811223522415801*^9}, {3.811223557239375*^9, 3.8112234398007393`*^9, 3.811223522415801*^9}, {3.811223557239375*^9,
3.811223581794196*^9}, 3.8112236218642673`*^9, {3.811223696961484*^9, 3.811223581794196*^9}, 3.8112236218642673`*^9, {3.811223696961484*^9,
3.811223753295334*^9}}, 3.811223753295334*^9}, {3.8122064089901867`*^9, 3.812206417622491*^9}},
CellLabel-> CellLabel->"In[31]:=",ExpressionUUID->"0daea3d8-8f63-4ebd-ba0f-4b294d15916e"],
"In[432]:=",ExpressionUUID->"0daea3d8-8f63-4ebd-ba0f-4b294d15916e"],
Cell[BoxData[ Cell[BoxData[
RowBox[{ RowBox[{
@ -14132,9 +14191,8 @@ Cell[BoxData[
3.811223185687756*^9, {3.811223309105949*^9, 3.811223364760091*^9}, { 3.811223185687756*^9, {3.811223309105949*^9, 3.811223364760091*^9}, {
3.811223479081833*^9, 3.811223582933836*^9}, 3.811223700996222*^9, { 3.811223479081833*^9, 3.811223582933836*^9}, 3.811223700996222*^9, {
3.8112237386115932`*^9, 3.811223753676*^9}, 3.811224082055791*^9, 3.8112237386115932`*^9, 3.811223753676*^9}, 3.811224082055791*^9,
3.811232556096879*^9}, 3.811232556096879*^9, 3.812206376976618*^9, 3.812206418637349*^9},
CellLabel-> CellLabel->"Out[39]=",ExpressionUUID->"6aabe757-b82c-493a-bfe6-83e278b98fd5"],
"Out[440]=",ExpressionUUID->"874fb105-0c41-4067-8df9-4692f067686e"],
Cell[BoxData[ Cell[BoxData[
RowBox[{ RowBox[{
@ -14163,15 +14221,14 @@ Cell[BoxData[
3.811223185687756*^9, {3.811223309105949*^9, 3.811223364760091*^9}, { 3.811223185687756*^9, {3.811223309105949*^9, 3.811223364760091*^9}, {
3.811223479081833*^9, 3.811223582933836*^9}, 3.811223700996222*^9, { 3.811223479081833*^9, 3.811223582933836*^9}, 3.811223700996222*^9, {
3.8112237386115932`*^9, 3.811223753676*^9}, 3.811224082055791*^9, 3.8112237386115932`*^9, 3.811223753676*^9}, 3.811224082055791*^9,
3.811232556097681*^9}, 3.811232556096879*^9, 3.812206376976618*^9, 3.812206418639378*^9},
CellLabel-> CellLabel->"Out[41]=",ExpressionUUID->"ac7444a1-9cea-44e7-bee8-9d5f30ed62ec"],
"Out[442]=",ExpressionUUID->"a7a8dc3c-2a42-42ae-98e4-bf4b227fb1b0"],
Cell[BoxData[ Cell[BoxData[
TagBox[ TagBox[
GraphicsBox[{{{{}, {{{}, {}, GraphicsBox[{{{{}, {{{}, {},
{GrayLevel[0], PointSize[0.012833333333333334`], Thickness[Large], {GrayLevel[0], PointSize[0.012833333333333334`], Thickness[Large],
LineBox[{{-1., -863.}, {0., -863.}}]}, LineBox[{{-1000., -863.}, {0., -863.}}]},
{RGBColor[0, 0, 1], PointSize[0.012833333333333334`], Thickness[ {RGBColor[0, 0, 1], PointSize[0.012833333333333334`], Thickness[
Large], LineBox[{{-203.62591770977707`, -600.}, \ Large], LineBox[{{-203.62591770977707`, -600.}, \
{-196.22449999999958`, -608.3896670999991}, {-166.2804400000084, \ {-196.22449999999958`, -608.3896670999991}, {-166.2804400000084, \
@ -14210,7 +14267,7 @@ Cell[BoxData[
Thickness[Large], Thickness[Large],
GrayLevel[0]]}], GrayLevel[0]]}],
TraditionalForm], {0., 0.}, Automatic, Scaled[ TraditionalForm], {0., 0.}, Automatic, Scaled[
9.75]], {{{-1., -863.}}, {{0., -863.}}}]}, 9.75]], {{{-1000., -863.}}, {{0., -863.}}}]},
{RGBColor[0, 0, 1], PointSize[0.012833333333333334`], Thickness[ {RGBColor[0, 0, 1], PointSize[0.012833333333333334`], Thickness[
Large], GeometricTransformationBox[InsetBox[ Large], GeometricTransformationBox[InsetBox[
FormBox[ FormBox[
@ -14282,7 +14339,17 @@ Cell[BoxData[
{{-64.38604999999598, -775.6955871000173}}, {{-55.509260000008, \ {{-64.38604999999598, -775.6955871000173}}, {{-55.509260000008, \
-786.8946071000096}}, {{-47.94706999999221, -797.3029471000075}}, \ -786.8946071000096}}, {{-47.94706999999221, -797.3029471000075}}, \
{{-41.26446000000783, -806.7485871000031}}, {{-35.14844000000039, \ {{-41.26446000000783, -806.7485871000031}}, {{-35.14844000000039, \
-814.8360771000114}}}]}}}, {{}, {}}}, {{{}, {}, -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]}, {}, {}, {}, {}, {}, {}, {}}}, {{}, {}}}, {{{}, {},
TagBox[ TagBox[
{RGBColor[0, 0, 1], Thickness[Large], Opacity[1.], {RGBColor[0, 0, 1], Thickness[Large], Opacity[1.],
Dashing[{Small, Small}], LineBox[CompressedData[" Dashing[{Small, Small}], LineBox[CompressedData["
@ -14315,7 +14382,7 @@ fu0b5+NQkdtZNxPLRwiI7FT1g8xOMrt5xl44SkBH58++C8qnpCOZlTo+RuC3
lti3n5j1ZCKjd+Pt9wRm8iKot5TVZONlB/WhDwSmhvLc2pXF5AKr+DMrxwlY lti3n5j1ZCKjd+Pt9wRm8iKot5TVZONlB/WhDwSmhvLc2pXF5AKr+DMrxwlY
pf578ezip1URWfVeI1r+D/vPmz8= pf578ezip1URWfVeI1r+D/vPmz8=
"]]}, "]]},
Annotation[#, "Charting`Private`Tag$28428#1"]& ], Annotation[#, "Charting`Private`Tag$8081#1"]& ],
TagBox[ TagBox[
{RGBColor[1, 0, 0], Thickness[Large], Opacity[1.], {RGBColor[1, 0, 0], Thickness[Large], Opacity[1.],
Dashing[{Small, Small}], LineBox[CompressedData[" Dashing[{Small, Small}], LineBox[CompressedData["
@ -14348,7 +14415,7 @@ pqKVjvBck08NkyAYrTW9vFd0QbmSyf5EwtDwj+4cSkcvzXX2bRohkVQf3zPB
q6GDPj/9JmOUxFRemNEtqorO1PZHcf8lMfYhj9VIldCNE54TNmMklp2azP3N q6GDPj/9JmOUxFRemNEtqorO1PZHcf8lMfYhj9VIldCNE54TNmMklp2azP3N
Tle5UJpHtkzzfwUymE8= Tle5UJpHtkzzfwUymE8=
"]]}, "]]},
Annotation[#, "Charting`Private`Tag$28428#2"]& ]}, {}}}, InsetBox[ Annotation[#, "Charting`Private`Tag$8081#2"]& ]}, {}}}, InsetBox[
TemplateBox[{ TemplateBox[{
GraphicsBox[{ GraphicsBox[{
Thickness[0.00970591089973794], Thickness[0.00970591089973794],
@ -15886,7 +15953,6 @@ mxg5fFi0XuFshpLDwbbl4acOmTt023juStsEDT9PC4i8hTLEf4kWkHzzWBXO
fwJKdwm6cP6L4q2iv7n1HcSnXuHMEELwwfoWm8H5YP+FmEHy5389iH0/TSHp fwJKdwm6cP6L4q2iv7n1HcSnXuHMEELwwfoWm8H5YP+FmEHy5389iH0/TSHp
/rOew1+QeRNNIenii77DTWBwGrWaOqDnfwB6Osfc /rOew1+QeRNNIenii77DTWBwGrWaOqDnfwB6Osfc
"]], "]],
FilledCurveBox[{{{1, 4, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3, FilledCurveBox[{{{1, 4, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3,
3}, {0, 1, 0}, {0, 1, 0}, {1, 3, 3}, {0, 1, 0}, {1, 3, 3}, {0, 3}, {0, 1, 0}, {0, 1, 0}, {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, 4, 3}, { 1, 0}, {0, 1, 0}, {0, 1, 0}, {1, 3, 3}, {0, 1, 0}}, {{1, 4, 3}, {
@ -15920,7 +15986,6 @@ D+fD3Ot9gt12tq4jnL+6+3YGg7kTnH9/H98c4ywnB423vPsMbho4HNe0mnR6
vpNDwC3pmsRNug6fNwRkz7qO4IOdedLJoXAN0KAPhlD/WsP5PiD7TK3h6m8A vpNDwC3pmsRNug6fNwRkz7qO4IOdedLJoXAN0KAPhlD/WsP5PiD7TK3h6m8A
KaOlaHxRG4f9tbIW6VcMIO6tt4G7ZwbI/E47jPCG8QGen8Le KaOlaHxRG4f9tbIW6VcMIO6tt4G7ZwbI/E47jPCG8QGen8Le
"]], "]],
FilledCurveBox[{{{0, 2, 0}, {1, 3, 3}, {0, 1, 0}, {0, 1, 0}, {1, 3, FilledCurveBox[{{{0, 2, 0}, {1, 3, 3}, {0, 1, 0}, {0, 1, 0}, {1, 3,
3}, {1, 3, 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}, {0, 1, 0}, {1, 3, 3}, {1, 3,
3}, {1, 3, 3}, {0, 1, 0}, {0, 1, 0}, {0, 1, 0}, {0, 1, 3}, {1, 3, 3}, {0, 1, 0}, {0, 1, 0}, {0, 1, 0}, {0, 1,
@ -16033,6 +16098,7 @@ uH2v23K7jX5LwPkzwUDcYT9Q+bbVURh8mHqddrGb5/yjHIzBQNIBFp8o9gH5
108.10199999999998`, 11.8656}, {106.01099999999997`, 11.8656}, { 108.10199999999998`, 11.8656}, {106.01099999999997`, 11.8656}, {
105.64699999999999`, 11.8656}, {105.392, 11.882799999999998`}, { 105.64699999999999`, 11.8656}, {105.392, 11.882799999999998`}, {
104.97300000000001`, 11.9203}, {104.97300000000001`, 16.6484}}}], 104.97300000000001`, 11.9203}, {104.97300000000001`, 16.6484}}}],
FilledCurveBox[{{{0, 2, 0}, {1, 3, 3}, {0, 1, 0}, {0, 1, 0}, {0, 1, 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}, {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, 0}, {0, 1, 0}, {0, 1, 0}, {0, 1, 0}, {1, 3, 3}, {0, 1,
@ -16194,7 +16260,7 @@ n3VZe0oep8PDF+b+KlD4cos7oMcPAMzCv+I=
AbsoluteThickness[1.6], AbsoluteThickness[1.6],
Thickness[Large], Thickness[Large],
GrayLevel[0]], GrayLevel[0]],
Line[{{-1., -863.}, {0., -863.}}]}, { Line[{{-1000., -863.}, {0., -863.}}]}, {
Hue[0.9060679774997897, 0.6, 0.6], Hue[0.9060679774997897, 0.6, 0.6],
Directive[ Directive[
PointSize[0.012833333333333334`], PointSize[0.012833333333333334`],
@ -16250,7 +16316,7 @@ n3VZe0oep8PDF+b+KlD4cos7oMcPAMzCv+I=
AbsoluteThickness[1.6], AbsoluteThickness[1.6],
Thickness[Large], Thickness[Large],
GrayLevel[0]]}], {0., 0.}, Automatic, GrayLevel[0]]}], {0., 0.}, Automatic,
Scaled[9.75]], {{{-1., -863.}}, {{0., -863.}}}]}, { Scaled[9.75]], {{{-1000., -863.}}, {{0., -863.}}}]}, {
Directive[ Directive[
PointSize[0.012833333333333334`], PointSize[0.012833333333333334`],
AbsoluteThickness[1.6], AbsoluteThickness[1.6],
@ -16333,7 +16399,38 @@ n3VZe0oep8PDF+b+KlD4cos7oMcPAMzCv+I=
{{-64.38604999999598, -775.6955871000173}}, {{-55.509260000008, \ {{-64.38604999999598, -775.6955871000173}}, {{-55.509260000008, \
-786.8946071000096}}, {{-47.94706999999221, -797.3029471000075}}, \ -786.8946071000096}}, {{-47.94706999999221, -797.3029471000075}}, \
{{-41.26446000000783, -806.7485871000031}}, {{-35.14844000000039, \ {{-41.26446000000783, -806.7485871000031}}, {{-35.14844000000039, \
-814.8360771000114}}}]}}}, {{}, {}}}, {{{{}, {}, -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]]}, {}, {}, {}, {}, {}, {}, {}}}, {{}, {}}}, {{{{}, {},
Annotation[{ Annotation[{
Directive[ Directive[
Opacity[1.], Opacity[1.],
@ -16370,8 +16467,7 @@ Wq45fEr7sppkWOSUiIe0+5cpTn8zB8gmQ4a5jpqA7jWh/2llN6nTF948PUzg
fu0b5+NQkdtZNxPLRwiI7FT1g8xOMrt5xl44SkBH58++C8qnpCOZlTo+RuC3 fu0b5+NQkdtZNxPLRwiI7FT1g8xOMrt5xl44SkBH58++C8qnpCOZlTo+RuC3
lti3n5j1ZCKjd+Pt9wRm8iKot5TVZONlB/WhDwSmhvLc2pXF5AKr+DMrxwlY lti3n5j1ZCKjd+Pt9wRm8iKot5TVZONlB/WhDwSmhvLc2pXF5AKr+DMrxwlY
pf578ezip1URWfVeI1r+D/vPmz8= pf578ezip1URWfVeI1r+D/vPmz8=
"]]}, "]]}, "Charting`Private`Tag$8081#1"],
"Charting`Private`Tag$28428#1"],
Annotation[{ Annotation[{
Directive[ Directive[
Opacity[1.], Opacity[1.],
@ -16409,7 +16505,7 @@ pqKVjvBck08NkyAYrTW9vFd0QbmSyf5EwtDwj+4cSkcvzXX2bRohkVQf3zPB
q6GDPj/9JmOUxFRemNEtqorO1PZHcf8lMfYhj9VIldCNE54TNmMklp2azP3N q6GDPj/9JmOUxFRemNEtqorO1PZHcf8lMfYhj9VIldCNE54TNmMklp2azP3N
Tle5UJpHtkzzfwUymE8= Tle5UJpHtkzzfwUymE8=
"]]}, "]]},
"Charting`Private`Tag$28428#2"]}}, {}}}, { "Charting`Private`Tag$8081#2"]}}, {}}}, {
DisplayFunction -> Identity, DisplayFunction -> Identity, AspectRatio -> DisplayFunction -> Identity, DisplayFunction -> Identity, AspectRatio ->
1, Axes -> {False, False}, AxesLabel -> {None, None}, 1, Axes -> {False, False}, AxesLabel -> {None, None},
AxesOrigin -> {0, -600.}, BaseStyle -> 18, DisplayFunction :> Identity, AxesOrigin -> {0, -600.}, BaseStyle -> 18, DisplayFunction :> Identity,
@ -16851,6 +16947,7 @@ mtnsDZrmqOHw4ZJvkkBBrsMZEFijDuEH5Drsza95O1NVwaFnep5Qs0euw5u2
3G6j3TKQ8HXJhaTrSik4H2z/TEk4v02BXfXMFnGHPys/XvL1y3WovP/jlnG3 3G6j3TKQ8HXJhaTrSik4H2z/TEk4v02BXfXMFnGHPys/XvL1y3WovP/jlnG3
mMMeUIDEIfiw/ALjg+2fkOtwEJRfmyUg6W8zND9ek3NAz98ArL7TnQ== mMMeUIDEIfiw/ALjg+2fkOtwEJRfmyUg6W8zND9ek3NAz98ArL7TnQ==
"]], "]],
FilledCurve[{{{1, 4, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {0, 1, FilledCurve[{{{1, 4, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {0, 1,
0}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {0, 1, 0}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {1, 3, 3}, {0, 1,
0}}}, CompressedData[" 0}}}, CompressedData["
@ -17929,16 +18026,15 @@ WOz7Hw==
3.811223185687756*^9, {3.811223309105949*^9, 3.811223364760091*^9}, { 3.811223185687756*^9, {3.811223309105949*^9, 3.811223364760091*^9}, {
3.811223479081833*^9, 3.811223582933836*^9}, 3.811223700996222*^9, { 3.811223479081833*^9, 3.811223582933836*^9}, 3.811223700996222*^9, {
3.8112237386115932`*^9, 3.811223753676*^9}, 3.811224082055791*^9, 3.8112237386115932`*^9, 3.811223753676*^9}, 3.811224082055791*^9,
3.811232556439719*^9}, 3.811232556096879*^9, 3.812206376976618*^9, 3.812206418796273*^9},
CellLabel-> CellLabel->"Out[42]=",ExpressionUUID->"1153d9f3-f21c-4cda-b1a2-659ea76a05a4"]
"Out[443]=",ExpressionUUID->"7c2de34b-d43d-497b-813a-04f2c9fe547d"] }, Open ]]
}, Open ]]
}, Open ]] }, Open ]]
}, Closed]]
}, Closed]]
}, },
WindowSize->{1184, 1394}, WindowSize->{1184, 1394},
WindowMargins->{{0, Automatic}, {Automatic, 0}}, WindowMargins->{{0, Automatic}, {Automatic, 0}},
FrontEndVersion->"12.1 for Mac OS X x86 (64-bit) (June 19, 2020)", FrontEndVersion->"12.1 for Mac OS X x86 (64-bit) (March 13, 2020)",
StyleDefinitions->"Default.nb", StyleDefinitions->"Default.nb",
ExpressionUUID->"ad7ac4fb-d7dd-4a37-9960-87ba52ede3bd" ExpressionUUID->"ad7ac4fb-d7dd-4a37-9960-87ba52ede3bd"
] ]
@ -17978,19 +18074,19 @@ Cell[278399, 5314, 195660, 3564, 70, "Output",ExpressionUUID->"1ef10c7e-62ec-42e
}, Closed]], }, Closed]],
Cell[CellGroupData[{ Cell[CellGroupData[{
Cell[474108, 8884, 168, 3, 38, "Subsection",ExpressionUUID->"32a5cbc6-7ca0-469d-b531-3ef832f9f00b"], Cell[474108, 8884, 168, 3, 38, "Subsection",ExpressionUUID->"32a5cbc6-7ca0-469d-b531-3ef832f9f00b"],
Cell[474279, 8889, 930, 19, 52, "Input",ExpressionUUID->"88b83480-ff85-479c-9bd5-86e50cd48da6"], Cell[474279, 8889, 925, 18, 52, "Input",ExpressionUUID->"88b83480-ff85-479c-9bd5-86e50cd48da6"],
Cell[CellGroupData[{ Cell[CellGroupData[{
Cell[475234, 8912, 8631, 198, 801, "Input",ExpressionUUID->"f0f1eea9-20b9-48df-9572-c510784075b4"], Cell[475229, 8911, 8626, 197, 801, "Input",ExpressionUUID->"f0f1eea9-20b9-48df-9572-c510784075b4"],
Cell[483868, 9112, 253908, 4746, 70, "Output",ExpressionUUID->"94181d5f-8832-4765-b54e-766c93d91980"] Cell[483858, 9110, 256303, 4808, 584, "Output",ExpressionUUID->"7d8b90dc-d960-4d55-affb-657fffbd0185"]
}, Open ]], }, Open ]],
Cell[CellGroupData[{ Cell[CellGroupData[{
Cell[737813, 13863, 10814, 242, 896, "Input",ExpressionUUID->"0daea3d8-8f63-4ebd-ba0f-4b294d15916e"], Cell[740198, 13923, 10861, 241, 896, "Input",ExpressionUUID->"0daea3d8-8f63-4ebd-ba0f-4b294d15916e"],
Cell[748630, 14107, 1881, 29, 70, "Output",ExpressionUUID->"874fb105-0c41-4067-8df9-4692f067686e"], Cell[751062, 14166, 1921, 28, 34, "Output",ExpressionUUID->"6aabe757-b82c-493a-bfe6-83e278b98fd5"],
Cell[750514, 14138, 1881, 29, 70, "Output",ExpressionUUID->"a7a8dc3c-2a42-42ae-98e4-bf4b227fb1b0"], Cell[752986, 14196, 1921, 28, 34, "Output",ExpressionUUID->"ac7444a1-9cea-44e7-bee8-9d5f30ed62ec"],
Cell[752398, 14169, 206314, 3764, 70, "Output",ExpressionUUID->"7c2de34b-d43d-497b-813a-04f2c9fe547d"] Cell[754910, 14226, 207999, 3803, 578, "Output",ExpressionUUID->"1153d9f3-f21c-4cda-b1a2-659ea76a05a4"]
}, Open ]]
}, Open ]]
}, Open ]] }, Open ]]
}, Closed]]
}, Closed]]
} }
] ]
*) *)

View File

@ -57,20 +57,20 @@ The same comment applies to the excited-state benchmark set of Thiel and coworke
Following a similar goal, we have recently proposed a large set of highly-accurate vertical transition energies for various types of excited states thanks to the renaissance of selected configuration interaction (SCI) methods \cite{Bender_1969,Huron_1973,Buenker_1974} which can now routinely produce near full configuration interaction (FCI) quality excitation energies for small- and medium-sized organic molecules. \cite{Loos_2018a,Loos_2019,Loos_2020a,Loos_2020b,Loos_2020c} Following a similar goal, we have recently proposed a large set of highly-accurate vertical transition energies for various types of excited states thanks to the renaissance of selected configuration interaction (SCI) methods \cite{Bender_1969,Huron_1973,Buenker_1974} which can now routinely produce near full configuration interaction (FCI) quality excitation energies for small- and medium-sized organic molecules. \cite{Loos_2018a,Loos_2019,Loos_2020a,Loos_2020b,Loos_2020c}
% The context % The context
In a recent preprint, \cite{Eriksen_2020} Eriksen \textit{et al.}~have proposed a blind test for a particular electronic structure problem inviting several groups around the world to contribute to this endeavour. In a recent article, \cite{Eriksen_2020} Eriksen \textit{et al.}~have proposed a blind test for a particular electronic structure problem inviting several groups around the world to contribute to this endeavour.
In addition to coupled cluster theory with singles, doubles, triples, and quadruples (CCSDTQ), \cite{Oliphant_1991,Kucharski_1992} a large panel of highly-accurate, emerging electronic structure methods were considered: In addition to coupled cluster theory with singles, doubles, triples, and quadruples (CCSDTQ), \cite{Oliphant_1991,Kucharski_1992} a large panel of highly-accurate, emerging electronic structure methods were considered:
(i) the many-body expansion FCI (MBE-FCI), \cite{Eriksen_2017,Eriksen_2018,Eriksen_2019a,Eriksen_2019b} (i) the many-body expansion FCI (MBE-FCI), \cite{Eriksen_2017,Eriksen_2018,Eriksen_2019a,Eriksen_2019b}
(ii) three SCI methods including a second-order perturbative correction \alert{[adaptive sampling CI (ASCI), \cite{Tubman_2016,Tubman_2018,Tubman_2020} iterative CI (iCI), \cite{Liu_2014,Liu_2016,Lei_2017,Zhang_2020} and semistochastic heat-bath CI (SHCI) \cite{Holmes_2016,Holmes_2017,Sharma_2017}]}, (ii) three SCI methods including a second-order perturbative correction [adaptive sampling CI (ASCI), \cite{Tubman_2016,Tubman_2018,Tubman_2020} iterative CI (iCI), \cite{Liu_2014,Liu_2016,Lei_2017,Zhang_2020} and semistochastic heat-bath CI (SHCI) \cite{Holmes_2016,Holmes_2017,Sharma_2017}],
(iii) \alert{the full coupled-cluster reduction (FCCR) \cite{Xu_2018,Xu_2020} which also includes a second-order perturbative correction}, (iii) the full coupled-cluster reduction (FCCR) \cite{Xu_2018,Xu_2020} which also includes a second-order perturbative correction,
(iv) the density-matrix renornalization group approach (DMRG), \cite{White_1992,White_1993,Chan_2011} and (iv) the density-matrix renornalization group (DMRG) approach, \cite{White_1992,White_1993,Chan_2011} and
(v) two flavors of FCI quantum Monte Carlo (FCIQMC), \cite{Booth_2009,Cleland_2010} namely AS-FCIQMC \cite{Ghanem_2019} and CAD-FCIQMC. \cite{Deustua_2017,Deustua_2018} (v) two flavors of FCI quantum Monte Carlo (FCIQMC), \cite{Booth_2009,Cleland_2010} namely AS-FCIQMC \cite{Ghanem_2019} and CAD-FCIQMC. \cite{Deustua_2017,Deustua_2018}
We refer the interested reader to Ref.~\onlinecite{Eriksen_2020} and its supporting information for additional details on each method and the complete list of references. We refer the interested reader to Ref.~\onlinecite{Eriksen_2020} and its supporting information for additional details on each method and the complete list of references.
Soon after, Lee \textit{et al.}~reported phaseless auxiliary-field quantum Monte Carlo \cite{Motta_2018} (ph-AFQMC) correlation energies for the very same problem. \cite{Lee_2020} Soon after, Lee \textit{et al.}~reported phaseless auxiliary-field quantum Monte Carlo \cite{Motta_2018} (ph-AFQMC) correlation energies for the very same problem. \cite{Lee_2020}
% The system % The system
The target application is the non-relativistic frozen-core correlation energy of the ground state of the benzene molecule in the cc-pVDZ basis. The target application is the non-relativistic frozen-core correlation energy of the ground state of the benzene molecule in the cc-pVDZ basis.
\alert{The geometry of benzene has been optimized at the MP2/6-31G* level \cite{Schreiber_2008} and its coordinates} can be found in the supporting information of Ref.~\onlinecite{Eriksen_2020} alongside its nuclear repulsion and Hartree-Fock energies. The geometry of benzene has been optimized at the MP2/6-31G* level \cite{Schreiber_2008} and its coordinates can be found in the supporting information of Ref.~\onlinecite{Eriksen_2020} alongside its nuclear repulsion and Hartree-Fock energies.
This corresponds to an active space of 30 electrons and 108 orbitals, \ie, the Hilbert space of benzene is of the order of $10^{35}$ Slater determinants. This corresponds to an active space of 30 electrons and 108 orbitals, \ie, the Hilbert space is of the order of $10^{35}$ Slater determinants.
Needless to say that this size of Hilbert space cannot be tackled by exact diagonalization with current architectures. Needless to say that this size of Hilbert space cannot be tackled by exact diagonalization with current architectures.
The correlation energies reported in Ref.~\onlinecite{Eriksen_2020} are gathered in Table \ref{tab:energy} alongside the best ph-AFQMC estimate from Ref.~\onlinecite{Lee_2020} based on a CAS(6,6) trial wave function. The correlation energies reported in Ref.~\onlinecite{Eriksen_2020} are gathered in Table \ref{tab:energy} alongside the best ph-AFQMC estimate from Ref.~\onlinecite{Lee_2020} based on a CAS(6,6) trial wave function.
The outcome of this work is nicely summarized in the abstract of Ref.~\onlinecite{Eriksen_2020}: The outcome of this work is nicely summarized in the abstract of Ref.~\onlinecite{Eriksen_2020}:
@ -79,12 +79,12 @@ The outcome of this work is nicely summarized in the abstract of Ref.~\onlinecit
%%% TABLE 1 %%% %%% TABLE 1 %%%
\begin{table} \begin{table}
\caption{ \caption{
The frozen-core correlation energy (in m$E_h$) of benzene in the cc-pVDZ basis set using various methods. The frozen-core correlation energy $\Delta E$ (in m$E_h$) of benzene in the cc-pVDZ basis set using various methods.
\label{tab:energy} \label{tab:energy}
} }
\begin{ruledtabular} \begin{ruledtabular}
\begin{tabular}{llc} \begin{tabular}{llc}
Method & \tabc{$E_c$} & Ref. \\ Method & \tabc{$\Delta E$} & Ref. \\
\hline \hline
ASCI & $-860.0$ & \onlinecite{Eriksen_2020} \\ ASCI & $-860.0$ & \onlinecite{Eriksen_2020} \\
iCI & $-861.1$ & \onlinecite{Eriksen_2020} \\ iCI & $-861.1$ & \onlinecite{Eriksen_2020} \\
@ -114,7 +114,7 @@ The outcome of this work is nicely summarized in the abstract of Ref.~\onlinecit
Right: $\Delta E_\text{var.}$ (in m$E_h$) as a function of $E_\text{PT2}$ or $E_\text{rPT2}$. Right: $\Delta E_\text{var.}$ (in m$E_h$) as a function of $E_\text{PT2}$ or $E_\text{rPT2}$.
The four-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. The theoretical estimate of $-863$ m$E_h$ from Ref.~\onlinecite{Eriksen_2020} is marked by a black line for comparison purposes.
\alert{The statistical error bars associated with $E_\text{PT2}$ or $E_\text{rPT2}$ (not shown) are of the order of the size of the markers.} The statistical error bars associated with $E_\text{PT2}$ or $E_\text{rPT2}$ (not shown) are of the order of the size of the markers.
\label{fig:CIPSI} \label{fig:CIPSI}
} }
\end{figure*} \end{figure*}
@ -193,9 +193,9 @@ The statistical error on $E_\text{(r)PT2}$, corresponding to one standard deviat
\end{table} \end{table}
% CIPSI % CIPSI
For the sake of completeness and our very own curiosity, we report in this Note the frozen-core correlation energy obtained with a fourth flavor of SCI known as \textit{Configuration Interaction using a Perturbative Selection made Iteratively} (CIPSI), \cite{Huron_1973} which also includes a second-order perturbative (PT2) correction. For the sake of completeness and our very own curiosity, we report in this Note the frozen-core correlation energy obtained with a fourth flavor of SCI known as \textit{Configuration Interaction using a Perturbative Selection made Iteratively} (CIPSI), \cite{Huron_1973} which also includes a PT2 correction.
In short, the CIPSI algorithm belongs to the family of SCI+PT2 methods. In short, the CIPSI algorithm belongs to the family of SCI+PT2 methods.
The idea behind such methods is to avoid the exponential increase of the size of the CI expansion by retaining the most energetically relevant determinants only, thanks to the use of a second-order energetic criterion to select perturbatively determinants in the FCI space. The idea behind such methods is to slow down the exponential increase of the size of the CI expansion by retaining the most energetically relevant determinants only, thanks to the use of a second-order energetic criterion to select perturbatively determinants in the FCI space.
However, performing SCI calculations rapidly becomes extremely tedious when one increases the system size as one hits the exponential wall inherently linked to these methods. However, performing SCI calculations rapidly becomes extremely tedious when one increases the system size as one hits the exponential wall inherently linked to these methods.
From a historical point of view, CIPSI is probably one of the oldest SCI algorithm. From a historical point of view, CIPSI is probably one of the oldest SCI algorithm.
@ -204,14 +204,14 @@ Recently, the determinant-driven CIPSI algorithm has been efficiently implemente
In particular, we were able to compute highly-accurate ground- and excited-state energies for small- and medium-sized molecules (including benzene). \cite{Loos_2018a,Loos_2019,Loos_2020a,Loos_2020b,Loos_2020c} In particular, we were able to compute highly-accurate ground- and excited-state energies for small- and medium-sized molecules (including benzene). \cite{Loos_2018a,Loos_2019,Loos_2020a,Loos_2020b,Loos_2020c}
CIPSI is also frequently used to provide accurate trial wave function for QMC calculations. \cite{Caffarel_2014,Caffarel_2016a,Caffarel_2016b,Giner_2013,Giner_2015,Scemama_2015,Scemama_2016,Scemama_2018,Scemama_2018b,Scemama_2019,Dash_2018,Dash_2019} CIPSI is also frequently used to provide accurate trial wave function for QMC calculations. \cite{Caffarel_2014,Caffarel_2016a,Caffarel_2016b,Giner_2013,Giner_2015,Scemama_2015,Scemama_2016,Scemama_2018,Scemama_2018b,Scemama_2019,Dash_2018,Dash_2019}
The particularity of the current implementation is that the selection step and the PT2 correction are computed \textit{simultaneously} via a hybrid semistochastic algorithm \cite{Garniron_2017,Garniron_2019} (which explains the statistical error associated with the PT2 correction in the following). The particularity of the current implementation is that the selection step and the PT2 correction are computed \textit{simultaneously} via a hybrid semistochastic algorithm \cite{Garniron_2017,Garniron_2019} (which explains the statistical error associated with the PT2 correction in the following).
\alert{Moreover, a renormalized version of the PT2 correction (dubbed rPT2 below) has been recently implemented and tested for a more efficient extrapolation to the FCI limit thanks to a partial resummation of the higher-order of perturbation. \cite{Garniron_2019} Moreover, a renormalized version of the PT2 correction (dubbed rPT2 below) has been recently implemented and tested for a more efficient extrapolation to the FCI limit thanks to a partial resummation of the higher-order of perturbation. \cite{Garniron_2019}
We refer the interested reader to Ref.~\onlinecite{Garniron_2019} where one can find all the details regarding the implementation of the rPT2 correction and the CIPSI algorithm.} We refer the interested reader to Ref.~\onlinecite{Garniron_2019} where one can find all the details regarding the implementation of the rPT2 correction and the CIPSI algorithm.
% Computational details % Computational details
Being late to the party, we obviously cannot report blindly our CIPSI results. Being late to the party, we obviously cannot report blindly our CIPSI results.
However, following the philosophy of Eriksen \textit{et al.} \cite{Eriksen_2020} and Lee \textit{et al.}, \cite{Lee_2020} we will report our results with the most neutral tone, leaving the freedom to the reader to make up his/her mind. However, following the philosophy of Eriksen \textit{et al.} \cite{Eriksen_2020} and Lee \textit{et al.}, \cite{Lee_2020} we will report our results with the most neutral tone, leaving the freedom to the reader to make up his/her mind.
We then follow our usual ``protocol'' \cite{Scemama_2018,Scemama_2018b,Scemama_2019,Loos_2018a,Loos_2019,Loos_2020a,Loos_2020b,Loos_2020c} by performing a preliminary SCI calculation using Hartree-Fock orbitals in order to generate a SCI wave function with at least $10^7$ determinants. We then follow our usual ``protocol'' \cite{Scemama_2018,Scemama_2018b,Scemama_2019,Loos_2018a,Loos_2019,Loos_2020a,Loos_2020b,Loos_2020c} by performing a preliminary SCI calculation using Hartree-Fock orbitals in order to generate a SCI wave function with at least $10^7$ determinants.
Natural orbitals are then computed based on this wave function, and a new, larger SCI calculation is performed with this new natural set of orbitals. Natural orbitals are then computed based on this wave function, and a new SCI calculation is performed with this new natural set of orbitals.
This has the advantage to produce a smoother and faster convergence of the SCI energy toward the FCI limit. This has the advantage to produce a smoother and faster convergence of the SCI energy toward the FCI limit.
The total SCI energy is defined as the sum of the variational energy $E_\text{var.}$ (computed via diagonalization of the CI matrix in the reference space) and a second-order perturbative correction $E_\text{(r)PT2}$ which takes into account the external determinants, \ie, the determinants which do not belong to the variational space but are linked to the reference space via a nonzero matrix element. The magnitude of $E_\text{(r)PT2}$ provides a qualitative idea of the ``distance'' to the FCI limit. The total SCI energy is defined as the sum of the variational energy $E_\text{var.}$ (computed via diagonalization of the CI matrix in the reference space) and a second-order perturbative correction $E_\text{(r)PT2}$ which takes into account the external determinants, \ie, the determinants which do not belong to the variational space but are linked to the reference space via a nonzero matrix element. The magnitude of $E_\text{(r)PT2}$ provides a qualitative idea of the ``distance'' to the FCI limit.
As mentioned above, SCI+PT2 methods rely heavily on extrapolation, especially when one deals with medium-sized systems. As mentioned above, SCI+PT2 methods rely heavily on extrapolation, especially when one deals with medium-sized systems.
@ -232,24 +232,24 @@ higher-lying $\pi$ [39,41--43,46,49,50,53--57,71--74,82--85,87,92,93,98];
higher-lying $\sigma$ [37,38,40,44,45,47,48,51,52,58--70,75--81,86,88--91,94--97,99--114].} higher-lying $\sigma$ [37,38,40,44,45,47,48,51,52,58--70,75--81,86,88--91,94--97,99--114].}
Like Pipek-Mezey, \cite{Pipek_1989} this choice of orbital windows allows to preserve a strict $\sigma$-$\pi$ separation in planar systems like benzene. 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. 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.
\alert{Indeed, localized orbitals significantly speed up the convergence of SCI calculations by taking benefit of the local character of electron correlation.\cite{Angeli_2003,Angeli_2009,BenAmor_2011,Suaud_2017,Chien_2018,Eriksen_2020}} Indeed, localized orbitals significantly speed up the convergence of SCI calculations by taking benefit of the local character of electron correlation.\cite{Angeli_2003,Angeli_2009,BenAmor_2011,Suaud_2017,Chien_2018,Eriksen_2020}
We, therefore, consider these energies more trustworthy, and we will base our best estimate of the correlation energy of benzene on these calculations. 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 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 four-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.
\alert{(In other words, the four largest variational wave functions are considered to perform the linear extrapolation.)} (In other words, the four largest variational wave functions are considered to perform the linear extrapolation.)
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. 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 % 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 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. 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, DMRG, 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 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). 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).
\footnote{\alert{Using the last 3, 4, 5, and 6 largest wave functions to perform the linear extrapolation yield the following correlation energy estimates: $-863.1(11)$, $-863.4(5)$, $-862.1(8)$, and $-863.5(11)$ mE$_h$, respectively. \footnote{Using the last 3, 4, 5, and 6 largest wave functions to perform the linear extrapolation yield the following correlation energy estimates: $-863.1(11)$, $-863.4(5)$, $-862.1(8)$, and $-863.5(11)$ mE$_h$, respectively.
These numbers vary by $1.4$ mE$_h$. These numbers vary by $1.4$ mE$_h$.
The four-point extrapolated value of $-863.4(5)$ mE$_h$ that we have chosen to report as our best estimate corresponds to the smallest fitting error. The four-point extrapolated value of $-863.4(5)$ mE$_h$ that we have chosen to report as our best estimate corresponds to the smallest fitting error.
Quadratic fits yield much larger variations and are discarded in practice. Quadratic fits yield much larger variations and are discarded in practice.
Due to the stochastic nature of $E_\text{rPT2}$, the fifth point is slightly off as compared to the others. Due to the stochastic nature of $E_\text{rPT2}$, the fifth point is slightly off as compared to the others.
Taking into account this fifth point yield a slightly smaller estimate of the correlation energy [$-862.1(8)$ mE$_h$], while adding a sixth point settles down the correlation energy estimate at $-863.5(11)$ mE$_h$ Taking into account this fifth point yield a slightly smaller estimate of the correlation energy [$-862.1(8)$ mE$_h$], while adding a sixth point settles down the correlation energy estimate at $-863.5(11)$ mE$_h$
}} }
For comparison, the best post blind test SHCI estimate is $-863.3$ m$E_h$, which agrees almost perfectly with our best CIPSI estimate, while the best post blind test ASCI and iCI correlation energies are $-861.3$ and $-864.15$ m$E_h$, respectively (see Table \ref{tab:extrap_dist_table}). For comparison, the best post blind test SHCI estimate is $-863.3$ m$E_h$, which agrees almost perfectly with our best CIPSI estimate, while the best post blind test ASCI and iCI correlation energies are $-861.3$ and $-864.15$ m$E_h$, respectively (see Table \ref{tab:extrap_dist_table}).
% Timings % Timings
@ -266,7 +266,7 @@ This work was performed using HPC resources from GENCI-TGCC (2020-gen1738) and f
PFL and AS have received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (Grant agreement No.~863481). PFL and AS have received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (Grant agreement No.~863481).
% Data availability statement % Data availability statement
\alert{The data that support the findings of this study are openly available in Zenodo at http://doi.org/10.5281/zenodo.4075286.} The data that support the findings of this study are openly available in Zenodo at \href{http://doi.org/10.5281/zenodo.4075286}{http://doi.org/10.5281/zenodo.4075286}.
%merlin.mbs apsrev4-1.bst 2010-07-25 4.21a (PWD, AO, DPC) hacked %merlin.mbs apsrev4-1.bst 2010-07-25 4.21a (PWD, AO, DPC) hacked
%Control: key (0) %Control: key (0)
@ -1114,7 +1114,7 @@ PFL and AS have received funding from the European Research Council (ERC) under
(\bibinfo {year} {2018})}\BibitemShut {NoStop}% (\bibinfo {year} {2018})}\BibitemShut {NoStop}%
\bibitem [{Note2()}]{Note2}% \bibitem [{Note2()}]{Note2}%
\BibitemOpen \BibitemOpen
\bibinfo {note} {\protect \leavevmode {\protect \color {black}Using the last \bibinfo {note} {\protect \leavevmode {\protect Using the
3, 4, 5, and 6 largest wave functions to perform the linear extrapolation 3, 4, 5, and 6 largest wave functions to perform the linear extrapolation
yield the following correlation energy estimates: $-863.1(11)$, $-863.4(5)$, yield the following correlation energy estimates: $-863.1(11)$, $-863.4(5)$,
$-862.1(8)$, and $-863.5(11)$ mE$_h$, respectively. These numbers vary by $-862.1(8)$, and $-863.5(11)$ mE$_h$, respectively. These numbers vary by

BIN
fig1b.pdf

Binary file not shown.