From 566de9c3ca2ac73b5913e71ba81b91576a43d5ce Mon Sep 17 00:00:00 2001
From: Pierre-Francois Loos <pierrefrancois.loos@gmail.com>
Date: Fri, 3 Jun 2022 10:10:31 +0200
Subject: [PATCH] CASSCF fig

---
 Data/CASPT3.nb                      | 912 +++++++++++++++++++++++++++-
 Manuscript/CASPT3.tex               |  12 +-
 Manuscript/fig2a.pdf                | Bin 0 -> 20847 bytes
 Manuscript/{fig2.pdf => fig2b.pdf}  | Bin
 Response_Letter/Response_Letter.tex |  21 +-
 5 files changed, 914 insertions(+), 31 deletions(-)
 create mode 100644 Manuscript/fig2a.pdf
 rename Manuscript/{fig2.pdf => fig2b.pdf} (100%)

diff --git a/Data/CASPT3.nb b/Data/CASPT3.nb
index 237a76c..07224b5 100644
--- a/Data/CASPT3.nb
+++ b/Data/CASPT3.nb
@@ -10,10 +10,10 @@
 NotebookFileLineBreakTest
 NotebookFileLineBreakTest
 NotebookDataPosition[       158,          7]
-NotebookDataLength[   1289234,      25065]
-NotebookOptionsPosition[   1285992,      24999]
-NotebookOutlinePosition[   1286390,      25015]
-CellTagsIndexPosition[   1286347,      25012]
+NotebookDataLength[   1328597,      25927]
+NotebookOptionsPosition[   1325125,      25857]
+NotebookOutlinePosition[   1325522,      25873]
+CellTagsIndexPosition[   1325479,      25870]
 WindowFrame->Normal*)
 
 (* Beginning of Notebook Content *)
@@ -79,6 +79,864 @@ Cell[BoxData[
 
 Cell[CellGroupData[{
 
+Cell[BoxData[{
+ RowBox[{
+  RowBox[{"WFT", "=", 
+   RowBox[{"{", 
+    RowBox[{"\"\<CASSCF\>\"", ","}], "}"}]}], ";"}], "\[IndentingNewLine]", 
+ RowBox[{
+  RowBox[{"MAE", "=", 
+   RowBox[{
+    RowBox[{"Import", "[", "\"\<FCI-all.xlsx\>\"", "]"}], 
+    "\[LeftDoubleBracket]", 
+    RowBox[{"Sheet", ",", 
+     RowBox[{"3", ";;", "286"}], ",", 
+     RowBox[{"41", "+", 
+      RowBox[{"{", "18", "}"}]}]}], "\[RightDoubleBracket]"}]}], 
+  ";"}], "\[IndentingNewLine]", 
+ RowBox[{"graph", "=", "\[IndentingNewLine]", 
+  RowBox[{"Show", "[", 
+   RowBox[{
+    RowBox[{"{", "\[IndentingNewLine]", 
+     RowBox[{
+      RowBox[{"Histogram", "[", 
+       RowBox[{
+        RowBox[{"MAE", "\[LeftDoubleBracket]", 
+         RowBox[{";;", ",", "1"}], "\[RightDoubleBracket]"}], ",", "20", ",", 
+        "\"\<PDF\>\"", ",", 
+        RowBox[{"PlotTheme", "\[Rule]", "\"\<Scientific\>\""}], ",", 
+        RowBox[{"BaseStyle", "\[Rule]", "14"}], ",", 
+        RowBox[{"PlotRange", "\[Rule]", 
+         RowBox[{"{", 
+          RowBox[{
+           RowBox[{"{", 
+            RowBox[{
+             RowBox[{"-", "2.2"}], ",", 
+             RowBox[{"+", "2.2"}]}], "}"}], ",", "Automatic"}], "}"}]}], ",", 
+        
+        RowBox[{"FrameLabel", "\[Rule]", 
+         RowBox[{"{", 
+          RowBox[{"\"\<Error (eV)\>\"", ",", "None"}], "}"}]}], ",", 
+        RowBox[{"ChartStyle", "\[Rule]", 
+         RowBox[{"{", 
+          RowBox[{"Directive", "[", 
+           RowBox[{
+            RowBox[{"Opacity", "[", "0.5", "]"}], ",", "Black"}], "]"}], 
+          "}"}]}], ",", 
+        RowBox[{"FrameTicks", "\[Rule]", 
+         RowBox[{"{", 
+          RowBox[{"Automatic", ",", "None"}], "}"}]}], ",", 
+        RowBox[{"FrameStyle", "\[Rule]", 
+         RowBox[{"Directive", "[", 
+          RowBox[{"Thick", ",", "20", ",", "Black"}], "]"}]}], ",", 
+        RowBox[{"PlotLabel", "\[Rule]", 
+         RowBox[{"Style", "[", 
+          RowBox[{
+           RowBox[{
+            RowBox[{
+            "WFT", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], 
+            "<>", "\"\< MAE: \>\"", "<>", 
+            RowBox[{"ToString", "[", 
+             RowBox[{"SetAccuracy", "[", 
+              RowBox[{
+               RowBox[{"Mean", "[", 
+                RowBox[{"Abs", "[", 
+                 RowBox[{"DeleteCases", "[", 
+                  RowBox[{
+                   RowBox[{"MAE", "\[LeftDoubleBracket]", 
+                    RowBox[{";;", ",", "1"}], "\[RightDoubleBracket]"}], ",", 
+                   "\"\<\>\""}], "]"}], "]"}], "]"}], ",", "3"}], "]"}], 
+             "]"}], "<>", "\"\< eV\>\""}], ",", "20"}], "]"}]}]}], "]"}], 
+      "\[IndentingNewLine]", ",", "\[IndentingNewLine]", 
+      RowBox[{"SmoothHistogram", "[", 
+       RowBox[{
+        RowBox[{"DeleteCases", "[", 
+         RowBox[{
+          RowBox[{"MAE", "\[LeftDoubleBracket]", 
+           RowBox[{";;", ",", "1"}], "\[RightDoubleBracket]"}], ",", 
+          "\"\<\>\""}], "]"}], ",", "Automatic", ",", "\"\<PDF\>\"", ",", 
+        RowBox[{"PlotRange", "\[Rule]", 
+         RowBox[{"{", 
+          RowBox[{
+           RowBox[{"{", 
+            RowBox[{
+             RowBox[{"-", "2.2"}], ",", "2.2"}], "}"}], ",", "Automatic"}], 
+          "}"}]}], ",", 
+        RowBox[{"PlotStyle", "\[Rule]", 
+         RowBox[{"Directive", "[", 
+          RowBox[{"Black", ",", 
+           RowBox[{"Thickness", "[", "0.01", "]"}]}], "]"}]}], ",", 
+        RowBox[{"Filling", "\[Rule]", "Bottom"}], ",", 
+        RowBox[{"PlotTheme", "\[Rule]", "\"\<Web\>\""}]}], "]"}]}], 
+     "\[IndentingNewLine]", "}"}], ",", 
+    RowBox[{"ImageSize", "\[Rule]", "300"}]}], "]"}]}], "\[IndentingNewLine]", 
+ RowBox[{
+  RowBox[{"Export", "[", 
+   RowBox[{"\"\<CASSCF.pdf\>\"", ",", "%"}], "]"}], ";"}]}], "Input",
+ CellLabel->"In[64]:=",ExpressionUUID->"31fd4aa9-bc03-4d09-8e5a-c1a6e87694ae"],
+
+Cell[BoxData[
+ GraphicsBox[{{
+    {RGBColor[0.92, 0.488, 0.2432], EdgeForm[{Opacity[0.567], Thickness[
+     Small]}], {}, 
+     {RGBColor[0.92, 0.488, 0.2432], EdgeForm[{Opacity[0.567], Thickness[
+      Small]}], 
+      TagBox[
+       TooltipBox[
+        TagBox[
+         {GrayLevel[0], Opacity[0.5], 
+          TagBox[
+           DynamicBox[{
+             FEPrivate`If[
+              CurrentValue["MouseOver"], 
+              EdgeForm[{
+                GrayLevel[0.5], 
+                AbsoluteThickness[1.5], 
+                Opacity[0.66]}], {}, {}], 
+             
+             RectangleBox[{-1.2, 0}, {-1., 0.07547169811320754}, 
+              "RoundingRadius" -> 0]},
+            
+            ImageSizeCache->{{71.49000000000002, 84.95000000000002}, {
+             57.293578568398516`, 74.80609465871879}}],
+           "DelayedMouseEffectStyle"]},
+         StatusArea[#, 0.07547169811320754]& ,
+         TagBoxNote->"0.07547169811320754"],
+        StyleBox[
+        "0.07547169811320754`", {FontFamily -> "Times"}, StripOnInput -> 
+         False]],
+       Annotation[#, 
+        Style[0.07547169811320754, {FontFamily -> "Times"}], "Tooltip"]& ], 
+      TagBox[
+       TooltipBox[
+        TagBox[
+         {GrayLevel[0], Opacity[0.5], 
+          TagBox[
+           DynamicBox[{
+             FEPrivate`If[
+              CurrentValue["MouseOver"], 
+              EdgeForm[{
+                GrayLevel[0.5], 
+                AbsoluteThickness[1.5], 
+                Opacity[0.66]}], {}, {}], 
+             
+             RectangleBox[{-1., 0}, {-0.8, 0.1509433962264151}, 
+              "RoundingRadius" -> 0]},
+            
+            ImageSizeCache->{{84.45000000000002, 97.91000000000003}, {
+             40.28106247807825, 74.80609465871879}}],
+           "DelayedMouseEffectStyle"]},
+         StatusArea[#, 0.1509433962264151]& ,
+         TagBoxNote->"0.1509433962264151"],
+        StyleBox[
+        "0.1509433962264151`", {FontFamily -> "Times"}, StripOnInput -> 
+         False]],
+       Annotation[#, 
+        Style[0.1509433962264151, {FontFamily -> "Times"}], "Tooltip"]& ], 
+      TagBox[
+       TooltipBox[
+        TagBox[
+         {GrayLevel[0], Opacity[0.5], 
+          TagBox[
+           DynamicBox[{
+             FEPrivate`If[
+              CurrentValue["MouseOver"], 
+              EdgeForm[{
+                GrayLevel[0.5], 
+                AbsoluteThickness[1.5], 
+                Opacity[0.66]}], {}, {}], 
+             
+             RectangleBox[{-0.8, 0}, {-0.6, 0.4339622641509434}, 
+              "RoundingRadius" -> 0]},
+            
+            ImageSizeCache->{{97.41000000000003, 
+             110.87000000000003`}, {-23.51587286062278, 74.80609465871879}}],
+           "DelayedMouseEffectStyle"]},
+         StatusArea[#, 0.4339622641509434]& ,
+         TagBoxNote->"0.4339622641509434"],
+        StyleBox[
+        "0.4339622641509434`", {FontFamily -> "Times"}, StripOnInput -> 
+         False]],
+       Annotation[#, 
+        Style[0.4339622641509434, {FontFamily -> "Times"}], "Tooltip"]& ], 
+      TagBox[
+       TooltipBox[
+        TagBox[
+         {GrayLevel[0], Opacity[0.5], 
+          TagBox[
+           DynamicBox[{
+             FEPrivate`If[
+              CurrentValue["MouseOver"], 
+              EdgeForm[{
+                GrayLevel[0.5], 
+                AbsoluteThickness[1.5], 
+                Opacity[0.66]}], {}, {}], 
+             
+             RectangleBox[{-0.6, 0}, {-0.4, 0.2830188679245283}, 
+              "RoundingRadius" -> 0]},
+            
+            ImageSizeCache->{{110.37000000000003`, 123.83000000000003`}, {
+             10.509159320017773`, 74.80609465871879}}],
+           "DelayedMouseEffectStyle"]},
+         StatusArea[#, 0.2830188679245283]& ,
+         TagBoxNote->"0.2830188679245283"],
+        StyleBox[
+        "0.2830188679245283`", {FontFamily -> "Times"}, StripOnInput -> 
+         False]],
+       Annotation[#, 
+        Style[0.2830188679245283, {FontFamily -> "Times"}], "Tooltip"]& ], 
+      TagBox[
+       TooltipBox[
+        TagBox[
+         {GrayLevel[0], Opacity[0.5], 
+          TagBox[
+           DynamicBox[{
+             FEPrivate`If[
+              CurrentValue["MouseOver"], 
+              EdgeForm[{
+                GrayLevel[0.5], 
+                AbsoluteThickness[1.5], 
+                Opacity[0.66]}], {}, {}], 
+             
+             RectangleBox[{-0.4, 0}, {-0.2, 0.4528301886792453}, 
+              "RoundingRadius" -> 0]},
+            
+            ImageSizeCache->{{123.33000000000003`, 
+             136.79000000000002`}, {-27.769001883202847`, 74.80609465871879}}],
+           "DelayedMouseEffectStyle"]},
+         StatusArea[#, 0.4528301886792453]& ,
+         TagBoxNote->"0.4528301886792453"],
+        StyleBox[
+        "0.4528301886792453`", {FontFamily -> "Times"}, StripOnInput -> 
+         False]],
+       Annotation[#, 
+        Style[0.4528301886792453, {FontFamily -> "Times"}], "Tooltip"]& ], 
+      TagBox[
+       TooltipBox[
+        TagBox[
+         {GrayLevel[0], Opacity[0.5], 
+          TagBox[
+           DynamicBox[{
+             FEPrivate`If[
+              CurrentValue["MouseOver"], 
+              EdgeForm[{
+                GrayLevel[0.5], 
+                AbsoluteThickness[1.5], 
+                Opacity[0.66]}], {}, {}], 
+             
+             RectangleBox[{-0.2, 0}, {0., 0.7547169811320755}, 
+              "RoundingRadius" -> 0]},
+            
+            ImageSizeCache->{{136.29000000000002`, 
+             149.75000000000003`}, {-95.81906624448393, 74.80609465871879}}],
+           "DelayedMouseEffectStyle"]},
+         StatusArea[#, 0.7547169811320755]& ,
+         TagBoxNote->"0.7547169811320755"],
+        StyleBox[
+        "0.7547169811320755`", {FontFamily -> "Times"}, StripOnInput -> 
+         False]],
+       Annotation[#, 
+        Style[0.7547169811320755, {FontFamily -> "Times"}], "Tooltip"]& ], 
+      TagBox[
+       TooltipBox[
+        TagBox[
+         {GrayLevel[0], Opacity[0.5], 
+          TagBox[
+           DynamicBox[{
+             FEPrivate`If[
+              CurrentValue["MouseOver"], 
+              EdgeForm[{
+                GrayLevel[0.5], 
+                AbsoluteThickness[1.5], 
+                Opacity[0.66]}], {}, {}], 
+             
+             RectangleBox[{0., 0}, {0.2, 0.7735849056603774}, 
+              "RoundingRadius" -> 0]},
+            
+            ImageSizeCache->{{149.25000000000003`, 
+             162.71000000000004`}, {-100.07219526706402`, 74.80609465871879}}],
+           "DelayedMouseEffectStyle"]},
+         StatusArea[#, 0.7735849056603774]& ,
+         TagBoxNote->"0.7735849056603774"],
+        StyleBox[
+        "0.7735849056603774`", {FontFamily -> "Times"}, StripOnInput -> 
+         False]],
+       Annotation[#, 
+        Style[0.7735849056603774, {FontFamily -> "Times"}], "Tooltip"]& ], 
+      TagBox[
+       TooltipBox[
+        TagBox[
+         {GrayLevel[0], Opacity[0.5], 
+          TagBox[
+           DynamicBox[{
+             FEPrivate`If[
+              CurrentValue["MouseOver"], 
+              EdgeForm[{
+                GrayLevel[0.5], 
+                AbsoluteThickness[1.5], 
+                Opacity[0.66]}], {}, {}], 
+             
+             RectangleBox[{0.2, 0}, {0.4, 0.5283018867924528}, 
+              "RoundingRadius" -> 0]},
+            
+            ImageSizeCache->{{162.21000000000004`, 
+             175.67000000000004`}, {-44.78151797352311, 74.80609465871879}}],
+           "DelayedMouseEffectStyle"]},
+         StatusArea[#, 0.5283018867924528]& ,
+         TagBoxNote->"0.5283018867924528"],
+        StyleBox[
+        "0.5283018867924528`", {FontFamily -> "Times"}, StripOnInput -> 
+         False]],
+       Annotation[#, 
+        Style[0.5283018867924528, {FontFamily -> "Times"}], "Tooltip"]& ], 
+      TagBox[
+       TooltipBox[
+        TagBox[
+         {GrayLevel[0], Opacity[0.5], 
+          TagBox[
+           DynamicBox[{
+             FEPrivate`If[
+              CurrentValue["MouseOver"], 
+              EdgeForm[{
+                GrayLevel[0.5], 
+                AbsoluteThickness[1.5], 
+                Opacity[0.66]}], {}, {}], 
+             
+             RectangleBox[{0.4, 0}, {0.6, 0.5283018867924528}, 
+              "RoundingRadius" -> 0]},
+            
+            ImageSizeCache->{{175.17000000000004`, 
+             188.63000000000002`}, {-44.78151797352311, 74.80609465871879}}],
+           "DelayedMouseEffectStyle"]},
+         StatusArea[#, 0.5283018867924528]& ,
+         TagBoxNote->"0.5283018867924528"],
+        StyleBox[
+        "0.5283018867924528`", {FontFamily -> "Times"}, StripOnInput -> 
+         False]],
+       Annotation[#, 
+        Style[0.5283018867924528, {FontFamily -> "Times"}], "Tooltip"]& ], 
+      TagBox[
+       TooltipBox[
+        TagBox[
+         {GrayLevel[0], Opacity[0.5], 
+          TagBox[
+           DynamicBox[{
+             FEPrivate`If[
+              CurrentValue["MouseOver"], 
+              EdgeForm[{
+                GrayLevel[0.5], 
+                AbsoluteThickness[1.5], 
+                Opacity[0.66]}], {}, {}], 
+             
+             RectangleBox[{0.6, 0}, {0.8, 0.4528301886792453}, 
+              "RoundingRadius" -> 0]},
+            
+            ImageSizeCache->{{188.13000000000002`, 
+             201.59000000000003`}, {-27.769001883202847`, 74.80609465871879}}],
+           "DelayedMouseEffectStyle"]},
+         StatusArea[#, 0.4528301886792453]& ,
+         TagBoxNote->"0.4528301886792453"],
+        StyleBox[
+        "0.4528301886792453`", {FontFamily -> "Times"}, StripOnInput -> 
+         False]],
+       Annotation[#, 
+        Style[0.4528301886792453, {FontFamily -> "Times"}], "Tooltip"]& ], 
+      TagBox[
+       TooltipBox[
+        TagBox[
+         {GrayLevel[0], Opacity[0.5], 
+          TagBox[
+           DynamicBox[{
+             FEPrivate`If[
+              CurrentValue["MouseOver"], 
+              EdgeForm[{
+                GrayLevel[0.5], 
+                AbsoluteThickness[1.5], 
+                Opacity[0.66]}], {}, {}], 
+             
+             RectangleBox[{0.8, 0}, {1., 0.24528301886792453`}, 
+              "RoundingRadius" -> 0]},
+            
+            ImageSizeCache->{{201.09000000000003`, 214.55000000000004`}, {
+             19.015417365177903`, 74.80609465871879}}],
+           "DelayedMouseEffectStyle"]},
+         StatusArea[#, 0.24528301886792453`]& ,
+         TagBoxNote->"0.24528301886792453"],
+        StyleBox[
+        "0.24528301886792453`", {FontFamily -> "Times"}, StripOnInput -> 
+         False]],
+       Annotation[#, 
+        Style[0.24528301886792453`, {FontFamily -> "Times"}], "Tooltip"]& ], 
+      TagBox[
+       TooltipBox[
+        TagBox[
+         {GrayLevel[0], Opacity[0.5], 
+          TagBox[
+           DynamicBox[{
+             FEPrivate`If[
+              CurrentValue["MouseOver"], 
+              EdgeForm[{
+                GrayLevel[0.5], 
+                AbsoluteThickness[1.5], 
+                Opacity[0.66]}], {}, {}], 
+             
+             RectangleBox[{1., 0}, {1.2, 0.18867924528301888`}, 
+              "RoundingRadius" -> 0]},
+            
+            ImageSizeCache->{{214.05000000000004`, 227.51000000000005`}, {
+             31.77480443291811, 74.80609465871879}}],
+           "DelayedMouseEffectStyle"]},
+         StatusArea[#, 0.18867924528301888`]& ,
+         TagBoxNote->"0.18867924528301888"],
+        StyleBox[
+        "0.18867924528301888`", {FontFamily -> "Times"}, StripOnInput -> 
+         False]],
+       Annotation[#, 
+        Style[0.18867924528301888`, {FontFamily -> "Times"}], "Tooltip"]& ], 
+      TagBox[
+       TooltipBox[
+        TagBox[
+         {GrayLevel[0], Opacity[0.5], 
+          TagBox[
+           DynamicBox[{
+             FEPrivate`If[
+              CurrentValue["MouseOver"], 
+              EdgeForm[{
+                GrayLevel[0.5], 
+                AbsoluteThickness[1.5], 
+                Opacity[0.66]}], {}, {}], 
+             
+             RectangleBox[{1.2, 0}, {1.4, 0.018867924528301886`}, 
+              "RoundingRadius" -> 0]},
+            
+            ImageSizeCache->{{227.01000000000005`, 240.47000000000003`}, {
+             70.05296563613872, 74.80609465871879}}],
+           "DelayedMouseEffectStyle"]},
+         StatusArea[#, 0.018867924528301886`]& ,
+         TagBoxNote->"0.018867924528301886"],
+        StyleBox[
+        "0.018867924528301886`", {FontFamily -> "Times"}, StripOnInput -> 
+         False]],
+       Annotation[#, 
+        Style[0.018867924528301886`, {FontFamily -> "Times"}], "Tooltip"]& ], 
+      
+      TagBox[
+       TooltipBox[
+        TagBox[
+         {GrayLevel[0], Opacity[0.5], 
+          TagBox[
+           DynamicBox[{
+             FEPrivate`If[
+              CurrentValue["MouseOver"], 
+              EdgeForm[{
+                GrayLevel[0.5], 
+                AbsoluteThickness[1.5], 
+                Opacity[0.66]}], {}, {}], 
+             
+             RectangleBox[{1.4, 0}, {1.6, 0.03773584905660377}, 
+              "RoundingRadius" -> 0]},
+            
+            ImageSizeCache->{{239.97000000000003`, 253.43000000000006`}, {
+             65.79983661355865, 74.80609465871879}}],
+           "DelayedMouseEffectStyle"]},
+         StatusArea[#, 0.03773584905660377]& ,
+         TagBoxNote->"0.03773584905660377"],
+        StyleBox[
+        "0.03773584905660377`", {FontFamily -> "Times"}, StripOnInput -> 
+         False]],
+       Annotation[#, 
+        Style[0.03773584905660377, {FontFamily -> "Times"}], "Tooltip"]& ], 
+      TagBox[
+       TooltipBox[
+        TagBox[
+         {GrayLevel[0], Opacity[0.5], 
+          TagBox[
+           DynamicBox[{
+             FEPrivate`If[
+              CurrentValue["MouseOver"], 
+              EdgeForm[{
+                GrayLevel[0.5], 
+                AbsoluteThickness[1.5], 
+                Opacity[0.66]}], {}, {}], 
+             
+             RectangleBox[{2., 0}, {2.2, 0.07547169811320754}, 
+              "RoundingRadius" -> 0]},
+            
+            ImageSizeCache->{{278.85, 292.31000000000006`}, {
+             57.293578568398516`, 74.80609465871879}}],
+           "DelayedMouseEffectStyle"]},
+         StatusArea[#, 0.07547169811320754]& ,
+         TagBoxNote->"0.07547169811320754"],
+        StyleBox[
+        "0.07547169811320754`", {FontFamily -> "Times"}, StripOnInput -> 
+         False]],
+       Annotation[#, 
+        Style[0.07547169811320754, {FontFamily -> "Times"}], 
+        "Tooltip"]& ]}, {}, {}}, {{}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, \
+{}, {}, {}, {}}}, {{}, GraphicsComplexBox[CompressedData["
+1:eJztm+VXVV3UtykRUUIpRRBQUEAkBBHUm3kAQaQVpEQBRaRDESkFKaVTkO4G
+KQlBWHR3d3d3x3ueMZ7nf3g/sL+ccY3fGmvvc8Y+a8754WJ6Y/L8HR4ODg4j
+IQ7O/3yGBP/PhVuK879XTGefk6L2Hvo/Lo+JkrOW30MJ7214fYPFYcz0PW3U
+oz2UdSNUV+eLBOAC51QV6x7i4NUTsXWSBEaSrYwFyj2kE38pgTpLGjADRdYX
+cffQjJWHr+oZOdBIdhC7v7SLPgp4qi+UyYPd56fkr/p2EQ6Vt7IgzXOIECMf
+cKjaRfj/NOlvOClACUVPfFLWLiJIjK8KOvcChsbCTZvDd5HIrzPrD5OV4ChD
+++Gm6y6qw9kcfqCtAnRfbxPSWuwioxEMMfMjNXgkvd4Kb3aR+me2sM931UGd
+9m/oO9ldpCCwkzgh8RpsZu103B/sooZysTsvMZoQkifOk3VzF20cPWtcMdKC
+QieSw+5Lu6gER5NlNP0N9D/vrD483kGwWqHHhasNe4yhvtcXdhDG/G57uO47
+uLyipS7Rs4MuexZHPpnSAYFi1lvGFTtIRXeEbvizLqi4r6z5Z+ygCX53ZMOk
+D59V8/79Dd1B50wp7uFPGEDgrS/fR77voB6ld8+Di40gb0v0+RnzHfTbkiuY
+KtsEuiuI6W9r7iDOaQHj+lFT2PJtm5GX3kFniv776LFtBpSav7ItBHaQtGOc
+8AjtR+Dj1PgSxryDqu/nFz+RNQfFQxaJcvIdpESyL/rS9xN8rF+8NHu4jVzj
+FqwjJi3A71fOEMncNoJc2R/hEpaQrWOdxNu1jfJzOVQXi62gnU/4o2rZNhKb
+HyGnFrOBNTwiIbv0bXTFRvcwc9gWyNuaieKDt1H4Vy1mBfevwB35s6PeeRuR
+n4Qb17DYg7yResSq2TbyyYkz3CX7BiYPb+hRv95GMsYzea4EDuB1bp73keQ2
+algaIUT4jvC7J/NYi38b8fjbxk2SOEFz/Oe679e30XO/t5RiN5xh6aNQQDrp
+NvqXz36J87ELXBA5o9Gxv4Vmme587zX9DhzkjWx701tIGSN5yT/5B0gP+21e
+69hCw6/0r8Ytu4JBmip6jLbQTfc99xcYd3CzZnTTT91CX/1y345HeECKxIyi
+T9AWIqPRpfYj9oJ66t8MeY5bSI1yWz/B0RvmJs3nB0y2UEti5jiXsw8Q5TzM
+xVXfQsXu3M5f5XyB9Rue/S2JLTT4/jrxKKMfSMjVScrwbSHSqAZ94z0/eE/v
+Q/WRcQvdO+iof9LnDy4LSqO/Lmwh/68qxj9KAyDhL31qye4m0m35JCKV8ROq
+v09+mpzcROmLfLdS4gNh6kUqhrhtE22OeZqVxQYBAfOH89zFmwgzaEIYl/wL
+mNcFul8kb6KnN+S2DfOD4XHpSZTNz02kkV365k5TCLz1qjaI/raJzmVEhO/M
+h4KDuid/jdEmSlp5zjlGHg4x7Iq4S6qbyDxxwutAKALKd2kbL4lvok7KJ9NK
+nyJhrHosUODuJlLOf+mFqxYFOD+TtF5f20TagelBLkdRwPjWhMOJeBPxX/Oj
+xsRHA4aHfyd5ewPZEqRR8D2LAY2Tw7KW8Q1UlRbzVgc3Fr42VXhsNW+gj1tm
+BwO5sRAR6qZ8tWgDcfksHQQax0GJ3rPrwokbSPX1xyeht+Nh6P7lJR3/DcTJ
+VEoztxgPh2dG8j3sNtAXxxcuDtkJQNcZ75BtsIE+3cJ5pWubCI9iDGV6lTeQ
+rHasYZhkErw05b18LLqBljxvNVyhTwZrof3xG9wbaPwozn9iIxlCLpSlP6Xb
+QCGfHSk2mlOgsP+7pQnRBor7UjH/5Hcq9CXJiv7cXEcXzzMvDvumwa4FFWnR
+6DpSHe+tzbNKh8tig72jjetIzGCNuObdbxCgiI0l/LuOhiAsm0QpA5TH9Iw5
+4teR1bekVWGWTLDI4BZ87ruONnnvCvtFZ0Lglx18yy/rSNakq2eTIQvypEqa
+w/XW0QezLUudqCzouuIcXPFiHRmH2RdNXM+GzRkp7TnhdSR94mWrl5gNlHmX
+uMg411FSf57Zxp0c4HPq2+OjXUcco9+vfM3LAYXnUZVqhOvoSlTQozPCf+AD
+43tv+/U1RHNFStG56Q/4Ld9RSxheQ45lMY17L3Mh+98mc2P9GtpO4LzzejEX
+2tyKVtby1lDYML3j7695sKriUEgTu4b+M/+rOn0pH8hvPXX+z3sN1fsJuuEm
+5wP3Fpn8W5s1VPuzgeMAUwByFd20ru/XkNM8LmdLfwEY+4ZP/VZYQ9mjdNJW
+Fn/BS0M7sxPW0HG7p8AWRSH8vnPbZv/2Guo+3GcXyimEpoM1McbLa+jXuNeo
+vGIRLNUVkIsTrKGPSllfb+4UwYVfdgMGq6tIcGIt5l/oP+DQEU/wHVxFNx/E
+LV0QKQZpPhKz/NpVxPlnsoJyvhgM8DofDv1ZRZ6pYspN/iXg1hpCiB+9ikJv
+JvziBQQpEVptrJ6rKESN2U54EUG9IWuYrNUqulOYHmslUwpzD1Z0zN+tIvTy
+mVJ0fCkQncvjCXm2in7NGhf/PC4F1h7bQ/TfKroyrByooFIGEvGiNVNsq0i6
+8dWZvuwyeP+R2O889SrCj2hMZyItBxfhNnUevFWkLHGF6LZ+OUgzcMv9HlxB
+EnWvPq5Vl4NbUfHk3p8VxDe5E/yBuQLqlaWsxTxX0CxNM0myQwUQbfaS+b5b
+QVoSrRd8xypAwkcnfvC/FfR48VDzpnAluHBsPmClXkEpZ7m/a0dXQnXtt9aP
+y8toXL5XSxyvCgjekemg6mWURUluW/22Ch7jhh8QRy6jFuYr90aqqsAhnN1X
+6fMy0uZu63FjrYZywYKbMXLLyIzlomO5ezXgdIv9W7q1jC6+nmv8ulINmA8d
+zwRxlpEkb3ZAgUINfCXVmnHqXUJvbbZ13hfUQEnKsm1r5hKa0zOicaSvhUNx
+20t0rksogShUH8+xFh5NnEt6r7WEHD0HaAdna8HaLui/HMElFKgU+eKMXB0U
+XmXpOL64hKhpNi6Y59bBbn62ruT8Ihq5/puSjq4eBBQxxz/LF5HIAgv7ukM9
+WKw2+Y+FLKK3c+uqc3P1kOfxku3Ox0VEGumitCXfAJuscyWWUovI2LtL62xB
+A/BVWShW3lhEaTVDFpcYGuGDFsE82eECyhvqjSJ0aYTsI1+7l50LSF7vKlfX
+YiOsBjNQJaYtIEFNQ01jhSbg5k9PWXdaQM1ayRcq/zaBcfsDjNArbB4VRFjO
+0Ay/jWu7XO8toNxwJjcV52ZYJFYy6CJZQIGi8SYW883AkTiBwzQ9j5pTb7gQ
+y7WAvqhZoGHJPIqSCaPCy2mBlJHj2wWB86gl3R3/5vlWmLXxKMM3mUc3Bg4+
+6j9qBdbLtMpyT+ZRhMWMtJhRK+j8SVwMYZhH4VfZdNPDWyFB/p7D9M4cAk2S
+9OCmVphcLKe52zqHpARqNnCOWoHZVT79S9IcMpa4J9J/uw3esAyL1NnPIQU8
+vDw6tTaIKTPopVSdQ5tGeN/R9zYYfbVnpMkzh/78d1Gy8E8bMO674Kedm0MP
+F2jHCcba4HUgZfDO2CxaLd5+EXShHSLuxnCKFs4iidaWWb377TDYzFXp5TeL
+hrqoAj9ptQOdQbFqv/4salyWuZnj1g5qZ6VWWERnkUDmYSxdTjuExPY6mV2d
+RcWNZV6p/e3QDzq0xRszqL6zvl0NtwMuD25kEDXOoBO2hwkstzpAxfKbmGLc
+DMqkcaMhku6AQEqygUjbGWREJ1VxZNIB3ZlhpguKM2gthVqWwL8DKGXYCe/f
+mUEmrkmtl3M7QHEuP9ThzAzqvNon9KC7A/ycxXiah6aRpk3J+ffbHdDO1FF9
+JW8a1RTUS0VSdQJ5iab6O69pZN6e6TrB2wnyastrmTrTSPB84P17zzrBa9vm
++6HQNPKzY2MNNOqEZr9z9BI00+hH652sc66dcIErKNt/ZQqR2xe1+sV2gnQD
+s8RIzRSSc7loxV/cCW7vs4fYo6aQr7njq72uTqjHx3y0sJxCDLq0sWNLnUAU
+1URULj+FpEi21rYIukDi0csIErYp1Klrsvbf1S5w6Z3lVcWdQhZOBsFl3F1Q
+bW5RF9c3iVYNTOncxLqA4CKBxmrWJDq7XVsVpdoFj9N9Nx+6TaIzMre6yIy6
+wOEpg9v3N5Mo8QMV/7hdF5RPpTF0PJhEWiaU3Df8ugDH4UHuNYpJpP9LYn84
+tgsw12ol9RcmkINcwLvruV3wtfDFaG7FBHp3n99jv6oLSpQmPuGGTaChPPsb
+5t1dcLhuel7GfAKp/RRIDZ/ugkfex1G/pCdQIMmesvN2F1jf9uCfZJ5Ahk94
+jwUJu6Gw5koj19E48rTEvVdP1Q27bxO1bLrGEU+9eJAYSzcI4NzbqU4fR3fJ
+rdsK+brBIqzc45LLOHrh6bos+Lgb8gTkr79+PY4EMFE7bQrdsNk5lJ/MP47y
+mX2r/N52A5+ZgcwW6ThiOqFZdv/YDR9I9sYxM2PIJ1FapsmxG7KTXSw90Bj6
+VkGc/y6gG1bFKEl7g8ZQUFb/X+34buAej469YTqG8GgxRn153WD8lUvQRGIM
+5YhRbbbXdsNv2uLmQsYxxBwupq0x0A2LeZLahHujyHlZhcVzuRs4FHr3nrWN
+orsvva7r4/WA/so77/DkURQtGLqJT9MDKe4bzHPfRhGrqP2n9xw9MHvrWyGf
+2ijiK/P9EC3SA6yVpPL2d0cR/rmUnGLVHtDRDJtqIB5FVtKXr9Wa9YD0nsI2
+Z/UIonzL8qjBrQdqeUyYqi1GkBxloFJbXA9I6LtJv7o1gq6kb4tOoh6ojIn/
+vNkzjAaPqnBJBnvg8UBpjPuPYVTDoounuNcDpRSDTdcFh5EF/SPaMppewEjv
+7P6dG0I9nZfuvbzfC/+cLjE/CxlCqQXeKqwqvfCo+I7crOQQqplW4Oe07oWC
+LQlru4NB9L0xMMYsvBcEOLXjqdMGkc+P5NTtsl74o2PXmq4+iKpsTD+UzvQC
+X2TIwWOSQWQxRhrdRtoHmT25NweLB1BK0UD17ft9wE3e9uyj8QCaUt4ma9Xs
+gzSJRVtihgFUJtruUOLeBxzfziZFt/SjT/dSLuzk90HS3+sdAvb9qMJHKPDz
+VB+wrv933MLdj6qzeS2EKPshjl2V7f1YH5p58WBC4nE/ML81Vzz27UMcQaI9
+/p/6ISrU2+6nSB/6aaTIT5vUD4ydKSkcG71o/uTJSs9AP4RdqO6qiO1F1MSa
+kdXkA0AnNobzUrEXtceJho6KD8CvL4e31wl60f7eofu1rwNwOY9G2TW3B42Q
+iTy3yhuAgOW7Dow6PSjveo307MoAUN6STc+n7kHm8Txsb9kHwUdDr1e2phvV
+a391G303COS/nPCnP3cjvLtW+0oxg+DRGsn5hbUbTbdElxSNDMKFc0WqlH1d
+KJZfmZng2hD8EO52SnXtQqlih0fMr4aAyHotQ+RBF9LScpq6FD4ETtkXBvrm
+O1EKXWB82dAQECzcIjQL7UQiI+RlzAzDYH9DlIdIuhNxfexN59IaBhz11+qR
+hx1oxy1evyt2GGwDrL7zp3egvKKUWqKZYThsDMhuetWB2uGB8D/2EbA8kzmk
+TdqBfEg8uMaNR2D3vwaiw5J29PasVq1hzgiYW0zz+pu0IwsbXHHZ3RHY/I2r
+wc7YjuzU6u/9uzYKpjN0bmWtbUiHrSvRWHgUVhkEclW+tSEfIUe1+9qjYKii
+MLrC04Zuo3vmZN9HYdHH+Pz38VZUq7/uu5I8Crp1rvzX/FuR9JTBcGvjKMzi
+xmvlirYimQ0rvbSVUag1t78SPNuCHMaUdGwujUHtVJ5JzdMWJPHO9ID/3hiU
+Ki1VbaU0o4jtvoheZSzX3KBjvtCMAq6SmKhaj0GBgNqH50ZNiIDj862cMCwn
++9TaNzeiXnb85yMlY5BJW3Mtg6sRPV+mr+8ZxbL7kfmQTwNaCL7zIxBvHJIO
+eRvOr9ejg8zdNCpmLBvpMz1QqEehDg3OcmLjEDUc9Vn3Tx3acBxME9LBslxP
+UyBVHcqhtvUdchmHX6UkzFUWtSgpFO/BrUQs8zy23uipQbYCJpnUNePgE2Pd
+yiRYgwKOVTQSp7FMkXVTPqQaEZ+n/tV5ZgJ+OM3Yfj2oQj9unjsIZsbyFn1H
+mnoVkuqudtkQmQB7HUW2geJKNGlOoDqkieUeN7tzDJVI09GgSuHrBFhKlHXd
+t69ApeXurxRCsfx357bOWDnSmaFs682fAFN2TocAkXL0iz2PaawDy6HaveWx
+ZajCbeGz3gqWL4RyrhGUocv4uhQmxJOg+6XNiUGnFFHml9qsME9C6yBZXqgX
+Qp1OgZqzQpNQWmL7Zle+BA36s68qq2A5ao70BUUxUnMm7nhoNgmZDkpFWV1F
+aNwgR9nPFcvaFe9JfxWiv2LvOlWjJyFKnJvSQO0v+u/fdplPAZZZw0tr6AoQ
+9aDr33stk+BDTGzEPJKHXrwj/fZkCsuLFle+Reci79sPUeX+JNg3T1QNvf2D
+yLTfZaSRTYF9pvyHBzdzELvZyvbRjSkw9Su+FjSbhVZmvXrz72PZnL1hIyUT
+mV/YetEnOQWlaZwrUf2/0duSN4Var7BsSiShppGGancyb8ibTEEU33gUxWQy
+Gp6SXoi1x/Ju0V6jbiJiKwkTUffF3u/fz+cuS3EozUzRyjIay/YmqfAhBmkH
+EI3tZGL3++FHtL4UiTDpayt9CMvk4oXaOyHo/N5kCX0zdv2vPYMenEB0X9Xo
+T/kAlhnT6SWJfZE8MG/WzGJZvizC6tYPVEb+Zp5jawpwyuQJSsxtETuRXtUq
+zjTglP4XePGeFmq1z4m7dAHLODcHVex0gIJWstWXGsujjxPHrtjBWZptEyNG
+LGs9ND1v7QozHhTJcWzTYM+qT+if4wuelHpR/HexHJWnyJgSCH1tSVbXH2D5
+Mn5selQoUCQEarwXxnJtomW4VBRsz45+IJDAMkEzw2R1DIie6y3alcEyZrOa
+XSQeaDLvGooqYNmW1tisOBE6ylYSJ5SnobQAQ1UgkAJPDPTjel9ieVPn33FO
+GogRR/gyaWCZ2/OtGFcGmJZdLyzWmgbNPNUz4l2ZwHbb5m3aW+x+nYsqc5ey
+wXDEnWJJG8vrdmke8jlQ0/35mt07LJNT4HB7/QGqBW7RN//DnAkKHQ25cCJ5
+6Uwgdn2UtGCixbl8OMi0f0mL3S9Kv3H/ypMC6CaBo2VNLP/QkC12+guLI0fH
+FK+xnLAerVleCJflxR1d1bDPV+m8hY/zD5jiiFjklbA8fvlp4n/FwHbymFT3
+GZZx0sIkbUqA6iSKtl0Ky9dgdakAQQXjOSJvsWkwvTYXLWFTCu8L4q9FCWE5
+PVuKlrYMXthdkiC8j+VHtlsLBWWAFtXZKjix3CAWWaxcDsPMdhK9LNjvo0b2
+1Hu7HG7LyCY8ocPyXO+65s8KEJh6RE1xCcuWMWF3+Sqh6onX0qOzWD5rKE7Q
+UQl2X92Daw6mwCfw3mqXWRVcvBFRl7uCZZaT4ETyavjU9egp8QSW/9SKWmVU
+w19ZHY7aLiyL+i1JytZAH1d+yHoN9v1vfxlEt1QDUvvpFJ5/sazFIrzsXos9
+l/z5QlOwvLo8j9jr4KZzdShjKJbtCgJ86+rg4vMoOSb3KcgkdRB6q1sPIdM6
+3NHWWA6XmuU72wAMGoEmcXpY5qDyI0xowM4dTEpcKlguGn7Y+7gRyIsi5cXF
+sf8fyaSp5IlGaPjFXLPEi+U+M28bhyYwYk9UYWfCsu5DQRkm7Jxn4u+FQ4rl
+HYKJa6XNcI917cwn7HnR6tLssfq6Ba6bzSl7TGOZ6hd/+VELMNTLfZVqx553
+7dyrT2Vb4SfH5Ym84kmYDaaqysbOdQZ4g/c7kibBVGs/+OpSKyxVh5HG+0/C
+KuuIsdOjNjh/ZVyI5ys2X60QXXJvg0shN8/Z6mLz/KTLSgNtgGu1Nun6fBIs
+7TyXStjboSR0hV/j0STsin8ov2XdDh9kxIRxbmJzUuUgn7p2sPnaSmtGhs27
+HhruXe6A8LEHC7m72PM+nFH4jW4H5D+M8WofmwCcd2eoG/I7gAwnrra+Hptz
+zM/znsXOSQbWj2JysPlmMwpT6oSE+lBd9TBsPSnKCTiT0Ak/F0Q8D50mgMjx
+l57xVifsZjcQuBhhc8kvQj2Pu4Ceziv/4AU2v/SGAhPQBZtkx0hTaAJ8+sRn
+kya6oLCefiX/5gSQR98uvsjbDXhptm74ZNhcl9zP2qEblpylHojvjAM595bO
+RFs3kEj1XnYYwdbDnb6H0kzYPpu/8d1fbD28jErIc017oH7Zn2s5A5u7xE7R
+l2L73uI2neu/sLnsj0IXsl4441NTr2yPrbdURt4rr3vhr00BlY/uODAOPdNW
++d0L5Evy9Y3y2DyOX7DsqBdkrj40JxPE5oZXSdll+oApJ0pLlQlbz3lxJvzC
++sCzoWI/9dw4sB5M5h8s9MGNkQF9ovUxSCqv89B+2A89jNl9Jv1jwOr2W6vJ
+rR8eNT2VnCjH9gvP/Pn5+/th7qCEVDt1DLivWJ6PZBuA8KKiyXV/bD6qPnrW
+agBWDx9T+Npi8yThXNPaAbjJ4/0O3mH7EZObbn00g5D1t7juWGYMBO6f1xB5
+PwiYMxJnW/mx+fEKb2reIPTIhor/YcDm1Z1ElIRDUEZRI5hGhO13PP8O2b4Y
+gsXi9t6CtVHAvIjInoobApMbxev9/aNQSuf4XXZzCApf3dKhqMTmk+/V80WH
+QcH6Hc7bdGx/lSrNw+g/DOeN3VzqA0dB4iMPoev4MLgfE1VI2mPzB9QDazwj
+IPyl3GdSbxQk8Q4y1L6NgCdXafwvhVForBtxqmgdAV+clor3/41CosQ4BzX5
+KNjoGbk/v4Xt79IUg+rvj0LrmH6W6sVRYCOvwbXTGIUO8abkg+0RmP8oaMiL
+7QvF6H7Pfm0YgdSe1O6Z36Pg99yXVihqBAwfXhMO6x6FKv+FE8FPI3An0idV
+/ngU0nOpX36UHIFlPHzqMzfHgOB70I95hhHI0Plk/xf7O3JtXCqP3BoG0/qZ
+eaNPY3D7vBmZf8Mw8HCqvbgePgatRXKqldHDsO7biLorx4D82dbyXcthyNkS
+YndfHAP24aS8YdlhMFfNCgDKcaAuieKqZhmGe8U3TjYejsM/9WKf2cMh2GYM
+1Et6Ow5U55bbRTuHIN+JqFPdfRxOivGTu1OHwHLWWuhizji0a7l9iXAcAkHp
+paSq/nEwHvhwK/jlEOxnaFBY402AxH/M/hW8Q1BE0f6Fk30CQi56kF69MAS2
+nx/Pjj+bgBkSgerQyUH4byDveZDVBNhpSINI8SAcC7EVS0VPwP3fjYTEgYOA
+YkJv4dRNwPWjN5YrxoNgT0jq92d1AmSTL4QuPRkEYX37Q93LkzCldun3CeMg
+4DWv69BjJkHl51Ncxv0BqOB519b2fhJE7mBIJTsGwOlnz0MX70ko51wnskwb
+ALG9pwkP8idBOGGUMs55AAhf/SNfGcb2XdHfuGtfD0BNKadNLOEUFPpGvxu+
+PwA/mKOnlDmnYP1OS/Aodv56+oNC/oLSFMzcD5OsnesH4kXnwtIv2L5mX4nC
+u7wfGuR2mT/FT8FToTetnKH94JGj783WNAXUxxIa0R/7QYZmaG9ocwrkQy2f
+jkr1A6mNnLYftu5pNzwImrrRDy3DZc3ij6dBZzycKPWgD3xE+AQPDKbBS4a+
+kb2jD54lJMRm+E/DcDERnkZKH1wivkKqXTQN/YRTNiLf+qDDyN3y8sQ0cPzF
+jqnKfRDQdjTeSDwDe9Nu07icfaB4z1Tm290ZyPAclOrE7wPq4PH8e2ozMJHt
+80+krxe6DxWvz3+bAUbV2HkJ7Dn0S7PGIyJ5BpSY5LsmHHpBpVJw53nbDByW
+xVWQKfcCLWua1tm9GWh9Ux9ez94L/e7XGosYZ0Ej0a7l4nEPhK348JtKzILD
+GxufidYeUFfAj2Y2nYXicMj4L7YHruV/Ot8XNAsrHtE3mD71wDDt7CdPNAvj
+S2mdHuI9EPVVbVR4Zha+Hng02dP0gOZ4o+Q26Rw4bv5zOprphuvikJvCPwc/
+eQU8CQq6YTw5i0Hj9RxQcA9fCvreDXEkzG4ULnOAG9LZna/UDdpmgZs16XMw
+H0kop8fSDTe7iDRsu+ZAwt7G6fdGF0wL2NRxH81BgG/GqltZFySFYasy8zyM
+aIaE4Hl3gS6OZkSw9DzcP+s5eEW9C9i124lkzedhNPy1XRdrF8zXPP6IFzYP
+25+s8h9g61Da7fyhvIp5oHltPytb1gmG3mwSBgvz8NL9TQmNZydwboRmM1As
+gMLHuBehKp2wrERK3/lgASo6btCN3uiEzEL77z/eLEBc7arm4nIHtt/aWHvk
+tgB9Htb3a/52wF2Hd+prWQuw6bxw/rNTB6xP9VTH9y3AeguNHr5cB/x5Ksmj
+hrsIw/ZMxx+udIB5+r9QUrZFWDHh6mucaAf+i1yEFfKLYEGa003zux22zaNN
+P1sugpIPSdNLy3Yo6KUYuB21CANPsxkiRdrB8pGL2GjNIjTellSZu4CdwqN2
+MwJWFmHIpqoI09MG+/gGtE9pliDlFk9MWnQb/Hs/5HQktARM3MxPeA3bwLZB
+biVLZwnuC309N8DfBkJc5ao6XkvwuS9pNRm3DY79+Cpp85Yg/3nZzYTGVijd
+TuBsGVqCkRjBkc6gVmy/eCXY8cwyREWsXnv0thVEStzxBe4sA+mvA+NRzlbA
+u35stKi4DIoTete9tlqgwtm0N8p2Ge5bltjOeraA09y4yIu4ZRC4VYrpvdkC
+YjIv0s81LgOBunaUNravIsyqoSnZWAZHPqm1ALVmqKF84PDh6gocy493Gm81
+wQ/LtMWboivg5p2sv+vTBE8HrykP6K+AmPX4Z/47TUCM8S3z9luB94NhhHfq
+sX1dLD7H48IVGLwXNDfyvhE8zloE7o6tAP0bOr9nhI2QcM/cvptgFe4cyjY5
+xjfA9/fpuZmXVyEk+slTS7EG0A2ennfjWIXPfgRBPNP1INHAwPgOswqvBiMX
+fn+vB9YjlReguAoSQ3wfVtnq4RyXn9sV3VXoyyR4utVYB/OaDWjDZhU49F/T
+F5vWQb0fwVaT9yqo98lniFHVQUrlf+xJsaswN03l7VVYC+7bFhoO+asAnv52
+AZq1YMiaGaDesAratY94lc7WgrTaXB3/yCq0b5Q6tv6uAQ6P6yfkG6tQyDab
+iaNcAyQlL/kWCNeg8HumyMRJNSyvBOhV0a7Bn+RYA4vkamhmao6I5FwDm0rt
+T3kK1fBb4WynlcgaDHeoiEUdV4G3M+acotIaXPAifsqdUgWm+VZCnPpr4FTC
+/UlbqQrk57I/En3Frn8w+UkQvwq4ry4mjfti13cl6KRlVsJFGZbhf/Fr8ITj
+jnHh60pY//qaIujvGngz4lprkVRCe2aQhFnTGnwa6uqN+FcB2eOtX6TG1kAP
+M0BrZlgB/pTEOSxba0CHwW/soqsAc3HR2ROidaC3jZdvbCoHRUtb+n66dcjy
+ynRVsCsHvpTc53+412H78+ALfZ5yoBpc/u71GMvflRXIJ8tgm4S1WFdlHS5U
+dvOJBZVBN2itixiuY4ezhEUiqTLIMwu5RW+/DgyF9JzqJ6UQFNuhvuO/Duc/
+GIlh/pSCZdcFv7bEddCU2zzO1SsFlbPiNalF64DDYlJTwlAKAoJ2h84t62AU
+a9VRkIbgikEBj+bEOlThNvpdv4pgP2xN58HOOtxPe7lR4lYC/c3sYZTnN2Cn
+cmPP/6AYCnG025avbcA7so6EBKNiCL0bTlh3dwM2aBY2j0b/ga1298NY8Q0I
+CC+rjnjxD9QDycy+qG2AbIMnj39DETyqlUhQNt6A91QySaMiRUC//22Ax2ED
+Orm+xnsVFcLx7SLyC4Eb4HLWgjH2XiEMv9oUm07eAJPFND22rL9Q4n3HprR4
+A6Q6FCdZOP9CZJlOZkjbBtzwj/gRm1YA9huRU+ZTG3DgthWWyFEAmix9tHJ7
+G3DLyN324e98wChfkmcj2QSHZM0VbZ58YHKVcsZn2oR2Fq84jrw8wCtyKhzi
+24T8n1a4wY/yYHyxeCVfYhNKzRRkMytzofzaDrOf+iZ4RJj+/CybC7Hy3GqG
+ppvA8JM2cbvvDzg56HmLO23CWMpHzbvv/4D2n5hKxl+bQL9oPce5nQOPpwf2
+9lM34R11IeOqSw6wXKbi6kKbEPFysdnqSg4QSspqZ3RsYvvRwPKW9GyYtvke
+7DqzCT9XyjI2RLOhOr20+e3BJjhrtujPDGRB4sgevhDZFsjUDrakfsqCHxd5
+BS/f2IL456POwhezQFfU0Hidfwvmuw2WE39ngsSn+NhGyS2YSm+9PCKTCWyJ
+w70Jr7cgvzLZsGYoA4j7aEi/fdgCzFK5m/6NDJgnfib60mULqPlIjsINf0P9
+IzfLeyFbkNmy3aL+Nx1SjSvSyX5vAQH+l8HQc+ngEXU4Ple2Ba8FFiWVX6WB
+Yfu9y5VdW+B825r+259UkCYwkYmY2wJHE4tMErJUuMOf5GB5tAWvxFV0Vg1T
+gFR3LP/5xW0QVkguvNqcDMvBtEscLNuwoBh+ze5uMjQ3KFw/K7gNr4cNmclD
+kiDjyEN5THob8KgY4ooIksCHq9qjSHMbPlj7yn38kIidV0/KfppvQ3+Krz7L
+RALI+wvsmPzYhkE+14RSpQTgqTLjkAzbhld3LN14m+Lh4k6KFnPmNqyYF703
+FI+HddbJwOOKbfDNyW5WLo+DdjX6xt6ebWCZCCiexsRBtocSbs7CNlixMO0Q
+lseCf4k3v+fJNhz+C7uUIBYL5qu1Bu8pdqB590lMZkMMKF7Hixa+tQPnAhR0
+LivGAJ/iw+6rD3dAGpZeNw1HA5WL+flt2R2YmN73LzaIhu38dEzrmx0Q5q8V
+7NuPgu656U8pFjvA9Kqd8qJHFORdZUx1ctsBQqJuCiWGKAiSUR19HbED7suf
+KO7pRWLnbT8qwewdMHtgNRITFAEqWQ2SFNU74KMhxTRQFw6CEwT2S3078EGR
+i7T+JAyuUAnl1iztgOjysdcLwTDYF/88H427C4e1n7WNLUKh3zKTwZZqF55a
+qP08zAuBopQ5RSW2XfDgU+fb3g+G0MHrbtz/7ULOhR5SCZFgsCVVR8TPdrHv
+yZ2lSY9foI75uTmpvQuGjmNyZf1B8N+HZjZkid2PbPdNx+0goI87qxHssQuW
+V4s7ye0D4bgLE/AxahcwV5iOPvf8hOGz1nUyf3ZhV8zj4JDnJyDBnONbtbtA
+86E7Lt47ACINFnnxBneBQWjfVW/VH+zDWfQGV3ahQ13ntrSiP2i2vI7Iw98D
+0V2PAckiPxDG/dXhQ7MHnzW55A1Y/ICJt43I4PYepNKwkfz28wW8d8RCYrAH
+esXCP2nxfWE8UPQjg8IevBq6vppr4QMVtbZJezp7kEkxWHdhwBti93OHOqz3
+IL5+cPhaixc4caxc+u21B8fR0QpldZ6g/ZpV4kfMHpDHteDhNHiAmI/Wlzd5
+e3DrG2v8Zrs7sJSHZD+q3wODdebkrDE3INzsmKEe3gPPOnFfuR1XmGYhoV9b
+2wPW8uqFeQpXqFEWf95wZh9s5706A+//gERXu+/xV/bh89FikPGb7/CjqOCf
+3Z19IAohaHcLcAHdpbU1VeF9oC7620Tc7AxPGW7f4nuxDzLtJTuEZM7A9kxb
+nVRvH076XIRClJ2A2DHcd9Z2HwZvCl8bTXSE+T/d1eU++xD9bPr2wYkDNEyT
+HYbF7UNDRPwdAk0HSL38lOdzwT6wLuJu0NR+Aw9JB51njfsgSMD6VEHgGxja
+FoXeHt0Hxzc58+2Z9iDze7P1zCb2+SclTXN47OHO6B3C0bMHUMna4Hjdwg5I
+L71/WHj1ABi/E128GPIVlkWjTAO4DiAqTQGj3/IFWj71xRuLHsDBrT8U76i/
+QEbipQEJ5QOIb6kvvG9iCz59UuQ3DA7gfvTyW7J+GzA97yx29PUAKHJYX3Mr
+28Cz/0qse/wO4FF4YE31tDXwmOxkZCUcgC5/PwG+qzVcjOaeci88gPNA+uOs
+kDWst+vR6jQfQCLF5tBVQmtoJ4iVw4wfwPPdzjKdMSvI4R90ot0+AL9XxhUk
+rVbgr0tVuHnuEPAnVJk5263APER2pZn+EKTq9bpPZqxAsfE7czLPIUzKCTIl
+XbSGe8elqo5ih5A6d7MX5KyBinvf65XqIcQlffYeiLKGbS3eyvtGh5BDeTIQ
+SWgD3f6Gexe/HQLT+uUfafY2kF8Vz7kYcAjcTXIq58hsIWhn+G110iH0Klnq
+VWbZgiXb5eCof4eQOFxLMaP7BVRePmu2bj2EtZI59098X0HQ0w3/xeQheB/Y
+2BtT2sEVVCHAtXsI0zUrD7QJ7GF/9dDo3IUjYEk27v37wx76r/PHTjAcgVwI
+nSgvzTcoUjTpLeY9AgZKwYJ/Gd8g1CWJ5NeTI9Df026Reu4AtgVjIh9eHgHd
+Tb+46mMHUJ+ntZQ2OYLWJzG2FLmO8B+dYvpNxyNYHHPmovzkBPSynuM4QUdw
+LPJ+M17IGY7tqmkGUo5gT8yKK/qiCwxnnUjnlhxBF8eIcN+SC6AJAQfv9iPQ
+e8KiTtfxHSKpPuTrTR8Bc3qFAm/ZD7B/krooun8ECnHU34byXUHTapLpGukx
+sErtfhrNcwPhVHrlXaZjCIs917ta7A5MQ0oe7feOYboSjlGjB+CR+ZSlPT2G
+9s8fH26Me8I4pm7b5dUx4LU+bHqE4w0VH/A4tMyOYczy+nn6Sz4QG/dQ66Hz
+MVx0qr4jnOgDTt3mgVTBx3A7srNTEXxBm+h3w0raMTCZcRAz9PuC2IMZnPrS
+YxD5rDD19rMfsBgy8sd1HkOTnk7wJLU/EEaoGnydPYZUZYLvGgX+MN3iF6Vy
+eAwmruLPfV4GQA1uY9dd8hOIH2w75ML5CYm8Z86TMJ9AbPhZgq34n/DjnRBm
+5v4JKO4pOvpLB4Ju0OdPZVInkIEXrFa8EQhP6zJTQjVOAE9QR5YxNAjYDuZG
+Pn08gbWeq3V6Ir+A+M4NKvnvJ+AxO4YrM/cL5l+rS7KHnkAg3WSPkXcwNPj8
+tCPIOIHlyC/XbflCILW8+c9w+Qncn/8rQdcbAh6bZ+cLuk9A9CEv4bZVKBje
+FGbwnz+BCltN/8wrYSCjYq1odHwCpVeK59fzwoC6qq/yxjkcjGf4LXLtZ+Fw
+9rXFCDkFDuZZ+A211Nlw2Nm+tH9Eh4PZkk74Z2obAbPeGZQLN3EwI67kVK8u
+REIvqzRXLzcOJv7algXtr0ioLZt9WvUAB7OpEPsohSwKCtSctbMf42CMxsuF
+B+WjIGmDyS5SFgfznF6UxMY7Cn55lAR7qOBglO9yB2o0RIEry8s/Vm9wMIfD
+vne0CaLBqmSnWccQB3OtfCVG+0E06CsHzClY4GA08hreCxlFg9oqN4GwPQ6G
+m3w4qzUsGiRdm65xuuFgRAdHQw9ro+HBdX3BqwE4mD99BY/j16KBvYhQkSgC
++3yLjYsZ1DFAqxhrvJWIg1Fl1LhIIBAD55fAdTwLBwP7/5zdX8TAgfNgbEsR
+DuZhauQxj0kMLF6zKvlXhYNJEMtqnHaKgcF8qr7kFhzMzVG5ldDAGGiUz94I
+7MPBvH8mnPI4Lgb+zcmSOk3gYDCIIK4nPQbSHBZYzZZwMBEVQg9lc2Ig7OoP
+0dc7OJjHxM884/7EgOcf5tdSuLiY1IA9q76sGPgiU2YpcB4XM026EbucEgPG
+06/8WahwMZlSAcazkTHw2m4//RIDLkYhxWy3yScGZC8H1Z6w4mJKO89SJH+J
+AaEs3onFu7gYTJ/NfRedGOCUbD3qe4SLecj+zdJcKgauTRherhHHxcTpW5s5
+c8QAme053j/yuBgK1uW7Tedi4IQyQSZaDRdDmmTI8WYiGlbTRXS9tHExBYVn
+KV/9jYZR8REHG2NcTD/pt4Ze92hoHbEJ17XExTg2KImuq0VDqeXlghcO2Pux
+9L2svRkNmRdz20U8cDHbshaVNitREJXybIkrEBdzxcywjzM3Cv7Xd8D8n99w
+6j+c+g+n/sOp/3DqP5z6D6f+w6n/cOo/nPoPp/7Dqf9w6j+c+g+n/sOp/3Dq
+P5z6D6f+w6n/cOo/nPoPp/7Dqf9w6j+c+g+n/sOp/3DqP5z6D6f+w6n/cOo/
+nPoPp/7Dqf9w6j+c+g+n/sOp/3DqP5z6D6f+w6n/cOo/nPoPp/7Dqf/w/7v/
+8P8AUD6mog==
+     "], {{{}, 
+       {GrayLevel[0], Opacity[0.3], EdgeForm[None], 
+        GraphicsGroupBox[PolygonBox[CompressedData["
+1:eJwl1kPUGIYSBtA/tu2ksW2bjW3btts4jdnYtm3btm0n7+Z1cc/MZmY1Z86X
+qGG7im0DBwQEBCJCkICA333EIP/14QlHWMIQmlCEJATBCUZQgvyeJRAB/LLo
+Jz/4zje+8oXPfOIjH3jPO97yhte84iUveM4znvKExzziIQ+4zz3ucofb3OIm
+N7jONa5yhctc4iIXOM85znKG05ziJCc4zjGOcoTDHOIgB9jPPvayh93sYic7
+2M42trKFzWxiIxtYzzrWsobVrGIlK1jOMpayhMUsYiELmM885jKH2cxiJjP4
+l+lMYypTmMwkJjKB8YxjLGP4h9GMYiQjGM4whjKEv/mLwQxiIAPoTz/60ofe
+9KInPehON7rShc50oiMdaE87fh9mG1rTipa0oDnNaEoTGtOIhjSgPvWoSx1q
+U4ua1KA61ahKFSpTiYpUoDzlKMuflKE0pShJCYpTjKIUoTCFKEgB8pOPvOQh
+N7nISQ6yk42sZCEzmchIBtKTjrSkITWpSEkKkpOMpCQhMX+QiIQkID7xiEsc
+YhOLmMQgOtGIShQiE4mIRCA84QhLGEITipCEIDjBCEoQAhOI30/kl4fykx98
+5xtf+cJnPvGRD7znHW95w2te8ZIXPOcZT3nCYx7xkAfc5x53ucNtbnGTG1zn
+Gle5wmUucZELnOccZznDaU5xkhMc5xhHOcJhDnGQA+xnH3vZw252sZMdbGcb
+W9nCZjaxkQ2sZx1rWcNqVrGSFSxnGUtZwmIWsZAFzGcec5nDbGYxkxn8y3Sm
+MZUpTGYSE5nAeMYxljH8w2hGMZIRDGcYQxnC3/zFYAYxkAH0px996UNvetGT
+HnSnG13pQmc60ZEOtKcdbWlDa1rRkhY0pxlNaUJjGtGQBtSnHnWpQ21qUZMa
+VKcaValCZSpRkQqUpxxl+ZMylKYUJSlBcYpRlCIUphAFKUB+8pGXPOQmFznJ
+QXaykZUsZCYTGclAetKRljSkJhUpSUFykpGUJCTmDxKRkATEJx5xiUNsYhGT
+GEQnGlGJQmQiEZEIhCccYQlDaEIRkhAEJxhBCUJgAgX6L0z8Un7yg+984ytf
++MwnPvKB97zjLW94zSte8oLnPOMpT3jMIx7ygPvc4y53uM0tbnKD61zjKle4
+zCUucoHznOMsZzjNKU5yguMc4yhHOMwhDnKA/exjL3vYzS52soPtbGMrW9jM
+JjaygfWsYy1rWM0qVrKC5SxjKUtYzCIWsoD5zGMuc5jNLGYyg3+ZzjSmMoXJ
+TGIiExjPOMYyhn8YzShGMoLhDGMoQ/ibvxjMIAYygP70oy996E0vetKD7nSj
+K13oTCc60oH2tKMtbWhNK1rSguY0oylNaEwjGtKA+tSjLnWoTS1qUoPqVKMq
+VahMJSpSgfKUoyx/UobSlKIkJShOMYpShMIUoiAFyE8+8pKH3OQiJznITjay
+koXMZCIjGUhPOtKShtSkIiUpSE4ykpKExPxBIhKSgPjEIy5xiE0sYhKD6EQj
+KlGITCQiEoHwhCMsYQhNKEISguAEI2jA/yN/wP8AXVqEEQ==
+          
+          "]]]}, {}, {}}, {{}, {}, 
+       {GrayLevel[0], PointSize[
+         NCache[
+          Rational[1, 120], 0.008333333333333333]], Thickness[0.01], CapForm[
+        "Butt"], LineBox[CompressedData["
+1:eJwl10O0GIYSANAX21ZjNrbd2LZt20Zj23Ya27Zt28a/OX9xz8xmlqOEDdpW
+aBMoICAgYpCAgD8xMNKAoAQjOCEISShCE4awhCM8Ef7UE4nIRCEq0YhODGIS
+i9jEIS7x+Iv4JCAhiUhMEpKSjOSkICWp+JvUpCEt6UhPBjKSicxkISvZyE4O
+cpKL3OQhL/nITwEKUojC/EMRilKM4pSgJKUoTRnKUo7yVKAilahMFapSjerU
+oCa1qE0d6lKP+jSgIY1oTBOa0ozmtKAlrWhNG9rSjvZ0oCOd6EwXutKN7vSg
+J73oTR/60o/+DGAggxjMEIYyjOGMYCT/MorRjGEs4xjPBCYyiclMYSrTmM4M
+ZjKL2cxhLvOYzwIWsojFLGEpy1jOClayitWsYS3rWM9/bGAjm9jMFrayje3s
+YCe72M0e9rKP/RzgIIc4zBGOcozjnOAkpzjNGc5yjvNc4CKXuMwVrnKN69zg
+Jre4zR3uco/7POAhj3jME57yjOe84CWveM0b3vKO93zgI5/4zBe+8o3v/OAn
+v/jNn+YPRGCCEJRgBCcEIQlFaMIQlnCEJwIRiURkohCVaEQnBjGJRWziEJd4
+/EV8EpCQRCQmCUlJRnJSkJJU/E1q0pCWdKQnAxnJRGaykJVsZCcHOclFbvKQ
+l3zkpwAFKURh/qEIRSlGcUpQklKUpgxlKUd5KlCRSlSmClWpRnVqUJNa1KYO
+dalHfRrQkEY0pglNaUZzWtCSVrTmz/BuSzva04GOdKIzXehKN7rTg570ojd9
+6Es/+jOAgQxiMEMYyjCGM4KR/MsoRjOGsYxjPBOYyCQmM4WpTGM6M5jJLGYz
+h7nMYz4LWMgiFrOEpSxjOStYySpWs4a1rGM9/7GBjWxiM1vYyja2s4Od7GI3
+e9jLPvZzgIMc4jBHOMoxjnOCk5ziNGc4yznOc4GLXOIyV7jKNa5zg5vc4jZ3
+uMs97vOAhzziMU94yjOe84KXvOI1b3jLO97zgY984jNf+Mo3vvODn/ziN38W
+fyACE4SgBCM4IQhJKEIThrCEIzwRiEgkIhOFqEQjOjGISSxiE4e4xOMv4pOA
+hCQiMUlISjKSk4KUpOJvUpOGtKQjPRnISCYyk4WsZCM7OchJLnKTh7zkIz8F
+KEghCvMPRShKMYpTgpKUojRlKEs5ylOBilSiMlWoSjWqU4Oa1KI2dahLPerT
+gIY0ojFNaEozmtOClrSiNW1oSzva04GOdKIzXehKN7rTg570ojd96Es/+jOA
+gQxiMEMYyjCGM4KR/MsoRjOGsYxjPBOYyCQmM4WpTGM6M5jJLGYzh7nMYz4L
+WMgiFrOEpSxjOStYySpWs4a1rGM9/7GBjWxiM1vYyja2s4Od7GI3e9jLPvZz
+gIMc4jBHOMoxjnOCk5ziNGc4yznOc4GLXOIyV7jKNa5zg5vc4jZ3uMs97vOA
+hzziMU94yjOe84KXvOI1b3jLO97zgY984jNf+Mo3vvODn/ziN3+O/kAEJghB
+CUZwQhCSUIQmDGEJR3giBPn///A/0nOEEQ==
+         "]]}}}], {{}, {}}}},
+  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
+  Axes->{False, False},
+  AxesLabel->{None, None},
+  AxesOrigin->{-2.2880000000000003`, 0},
+  BaseStyle->14,
+  Frame->True,
+  FrameLabel->{{None, None}, {
+     FormBox["\"Error (eV)\"", TraditionalForm], None}},
+  FrameStyle->Directive[
+    Thickness[Large], 20, 
+    GrayLevel[0]],
+  FrameTicks->{{None, None}, {Automatic, Automatic}},
+  GridLines->{{0}, {0}},
+  GridLinesStyle->Directive[
+    GrayLevel[0.5, 0.4]],
+  ImageSize->300,
+  LabelStyle->{FontFamily -> "Times"},
+  PlotLabel->FormBox[
+    StyleBox["\"CASSCF MAE: 0.47 eV\"", 20, StripOnInput -> False], 
+    TraditionalForm],
+  PlotRange->{{-2.2, 2.2}, {All, All}},
+  PlotRangePadding->{{
+     Scaled[0.02], 
+     Scaled[0.02]}, {None, 
+     Scaled[0.05]}},
+  Ticks->{Automatic, Automatic}]], "Output",
+ CellLabel->"Out[66]=",ExpressionUUID->"6388ff56-ba75-4aa8-9483-c26182754d88"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
 Cell[BoxData[{
  RowBox[{
   RowBox[{"WFT", "=", 
@@ -24998,8 +25856,8 @@ b5JAxM18hw2qT5rnrVWF88HuEdCA89HzC4wPAHciYt4=
 }, Open  ]]
 },
 WindowSize->{1440, 847},
-WindowMargins->{{1441, Automatic}, {Automatic, 0}},
-FrontEndVersion->"13.0 for Mac OS X ARM (64-bit) (December 2, 2021)",
+WindowMargins->{{Automatic, 126}, {Automatic, 0}},
+FrontEndVersion->"13.0 for Mac OS X x86 (64-bit) (December 2, 2021)",
 StyleDefinitions->"Default.nb",
 ExpressionUUID->"85b7a10d-2875-445d-bac3-ec6fc5ce6624"
 ]
@@ -25028,43 +25886,47 @@ Cell[2603, 70, 89, 0, 98, "Title",ExpressionUUID->"7bca69e7-ea88-42f1-8cb1-a3017
 Cell[2695, 72, 245, 5, 46, "Input",ExpressionUUID->"717b96b6-16dc-49c5-8f1c-5691e5e5da0f",
  InitializationCell->True],
 Cell[CellGroupData[{
-Cell[2965, 81, 5142, 122, 304, "Input",ExpressionUUID->"919b4508-e93a-4c58-acbd-bccf8883719c"],
-Cell[8110, 205, 126303, 2334, 531, "Output",ExpressionUUID->"015b06bf-1e91-4e1a-9cff-a6bd49cc491b"]
+Cell[2965, 81, 3811, 94, 262, "Input",ExpressionUUID->"31fd4aa9-bc03-4d09-8e5a-c1a6e87694ae"],
+Cell[6779, 177, 35282, 757, 272, "Output",ExpressionUUID->"6388ff56-ba75-4aa8-9483-c26182754d88"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[134450, 2544, 1392, 34, 94, "Input",ExpressionUUID->"f9587fd2-1e23-4030-9f68-19ee3273f3fd"],
-Cell[135845, 2580, 910, 23, 34, "Output",ExpressionUUID->"ddd7a8f9-59bf-40aa-b7d9-bc4f1657909d"]
+Cell[42098, 939, 5142, 122, 304, "Input",ExpressionUUID->"919b4508-e93a-4c58-acbd-bccf8883719c"],
+Cell[47243, 1063, 126303, 2334, 531, "Output",ExpressionUUID->"015b06bf-1e91-4e1a-9cff-a6bd49cc491b"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[136792, 2608, 33857, 784, 1858, "Input",ExpressionUUID->"d3cfcc37-7870-45d8-9dde-1f2d1fb25587"],
-Cell[170652, 3394, 1668, 42, 34, "Output",ExpressionUUID->"20eac2ff-49eb-454b-8032-5f9dde62980f"],
-Cell[172323, 3438, 965099, 18691, 2059, "Output",ExpressionUUID->"6574969c-81b8-467d-9bfb-ccfcce64c9d4"]
+Cell[173583, 3402, 1392, 34, 94, "Input",ExpressionUUID->"f9587fd2-1e23-4030-9f68-19ee3273f3fd"],
+Cell[174978, 3438, 910, 23, 34, "Output",ExpressionUUID->"ddd7a8f9-59bf-40aa-b7d9-bc4f1657909d"]
+}, Open  ]],
+Cell[CellGroupData[{
+Cell[175925, 3466, 33857, 784, 1858, "Input",ExpressionUUID->"d3cfcc37-7870-45d8-9dde-1f2d1fb25587"],
+Cell[209785, 4252, 1668, 42, 34, "Output",ExpressionUUID->"20eac2ff-49eb-454b-8032-5f9dde62980f"],
+Cell[211456, 4296, 965099, 18691, 2059, "Output",ExpressionUUID->"6574969c-81b8-467d-9bfb-ccfcce64c9d4"]
 }, Open  ]]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[1137471, 22135, 86, 0, 98, "Title",ExpressionUUID->"a87d4be2-7864-4d6c-a2e3-36b4dcb9831f"],
-Cell[1137560, 22137, 244, 5, 46, "Input",ExpressionUUID->"106aeec9-b1fc-431a-b4b9-d3a7cdcac07d",
+Cell[1176604, 22993, 86, 0, 98, "Title",ExpressionUUID->"a87d4be2-7864-4d6c-a2e3-36b4dcb9831f"],
+Cell[1176693, 22995, 244, 5, 46, "Input",ExpressionUUID->"106aeec9-b1fc-431a-b4b9-d3a7cdcac07d",
  InitializationCell->True],
 Cell[CellGroupData[{
-Cell[1137829, 22146, 4452, 128, 283, "Input",ExpressionUUID->"0f8cc5fb-1a0b-4172-b346-6b0a3699cd11"],
-Cell[1142284, 22276, 31728, 630, 318, "Output",ExpressionUUID->"4a3f128c-d621-4a5d-9883-46aa65f2e128"]
+Cell[1176962, 23004, 4452, 128, 283, "Input",ExpressionUUID->"0f8cc5fb-1a0b-4172-b346-6b0a3699cd11"],
+Cell[1181417, 23134, 31728, 630, 318, "Output",ExpressionUUID->"4a3f128c-d621-4a5d-9883-46aa65f2e128"]
 }, Open  ]]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[1174061, 22912, 152, 3, 98, "Title",ExpressionUUID->"adbad97c-5cc1-459d-b662-b85d0910cb88"],
+Cell[1213194, 23770, 152, 3, 98, "Title",ExpressionUUID->"adbad97c-5cc1-459d-b662-b85d0910cb88"],
 Cell[CellGroupData[{
-Cell[1174238, 22919, 3328, 70, 94, "Input",ExpressionUUID->"5d41447f-af63-414c-a72f-279d73505b87"],
-Cell[1177569, 22991, 53160, 880, 5156, "Output",ExpressionUUID->"b7b5c954-23f1-476d-83c3-e63def1b3658"]
+Cell[1213371, 23777, 3328, 70, 94, "Input",ExpressionUUID->"5d41447f-af63-414c-a72f-279d73505b87"],
+Cell[1216702, 23849, 53160, 880, 5156, "Output",ExpressionUUID->"b7b5c954-23f1-476d-83c3-e63def1b3658"]
 }, Open  ]]
 }, Closed]],
 Cell[CellGroupData[{
-Cell[1230778, 23877, 150, 3, 72, "Title",ExpressionUUID->"61cd51bf-e012-4d32-bf6d-b9b9c71fc72c"],
-Cell[1230931, 23882, 2518, 68, 229, "Input",ExpressionUUID->"ad61fbc1-6c66-408c-8f64-9aeb20e11f88",
+Cell[1269911, 24735, 150, 3, 72, "Title",ExpressionUUID->"61cd51bf-e012-4d32-bf6d-b9b9c71fc72c"],
+Cell[1270064, 24740, 2518, 68, 229, "Input",ExpressionUUID->"ad61fbc1-6c66-408c-8f64-9aeb20e11f88",
  InitializationCell->True],
-Cell[1233452, 23952, 2761, 73, 283, "Input",ExpressionUUID->"d48b5f18-7800-4eb1-91de-a74892343606"],
+Cell[1272585, 24810, 2761, 73, 283, "Input",ExpressionUUID->"d48b5f18-7800-4eb1-91de-a74892343606"],
 Cell[CellGroupData[{
-Cell[1236238, 24029, 2169, 44, 115, "Input",ExpressionUUID->"88c5defa-ce5a-47a9-9104-0124f20e76a8"],
-Cell[1238410, 24075, 47554, 920, 345, "Output",ExpressionUUID->"b7f2b372-fe4a-440e-bb84-c204b0cb3a05"]
+Cell[1275371, 24887, 2169, 44, 115, "Input",ExpressionUUID->"88c5defa-ce5a-47a9-9104-0124f20e76a8"],
+Cell[1277543, 24933, 47554, 920, 345, "Output",ExpressionUUID->"b7f2b372-fe4a-440e-bb84-c204b0cb3a05"]
 }, Open  ]]
 }, Open  ]]
 }
diff --git a/Manuscript/CASPT3.tex b/Manuscript/CASPT3.tex
index e7d8b0e..8043cf3 100644
--- a/Manuscript/CASPT3.tex
+++ b/Manuscript/CASPT3.tex
@@ -159,8 +159,8 @@ We underline that, although a third-order version of NEVPT has been developed \c
 \end{figure}
 %%% %%% %%% %%%
 
-For each compound represented in Fig.~\ref{fig:mol}, we have computed the CASPT2 and CASPT3 vertical excitation energies with Dunning's aug-cc-pVTZ 
-basis set. \cite{Kendall_1992} 
+For each compound represented in Fig.~\ref{fig:mol}, we have computed the \alert{single-state} CASPT2 and CASPT3 vertical excitation energies.
+\alert{All calculations reported in the present manuscript have been performed with Dunning's aug-cc-pVTZ basis set.} \cite{Kendall_1992} 
 Geometries and reference theoretical best estimates (TBEs) for the vertical excitation energies have been extracted from the QUEST database \cite{Veril_2021} and can be downloaded at \url{https://lcpq.github.io/QUESTDB_website}.
 
 All the CASPT2 and CASPT3 calculations have been carried out in the frozen-core approximation and within the RS2 and RS3 contraction schemes as implemented in MOLPRO and described in Refs.~\onlinecite{Werner_1996} and \onlinecite{Werner_2020}. 
@@ -501,8 +501,11 @@ TBEs listed as ``safe'' are assumed to be chemically accurate (\ie, absolute err
 
 %%% FIGURE 2 %%%
 \begin{figure}
-	\includegraphics[width=\linewidth]{fig2}
-	\caption{Histograms of the errors (in \si{\eV}) obtained for CASPT2 and CASPT3 with and without IPEA shift.
+	\includegraphics[width=0.5\linewidth]{fig2a}
+	\\
+	\vspace{0.01\textwidth}
+	\includegraphics[width=\linewidth]{fig2b}
+	\caption{Histograms of the errors (in \si{\eV}) obtained for \alert{CASSCF as well as} CASPT2 and CASPT3 with and without IPEA shift.
 	Raw data are given in Table \ref{tab:BigTab}.}
 	\label{fig:PT2_vs_PT3}
 \end{figure}
@@ -658,6 +661,7 @@ The two principal take-home messages of this study are that:
 (ii) CASPT2(IPEA) and CASPT3 have a very similar accuracy.
 These global trends are also true for specific sets of excitations and various system sizes.
 Therefore, if one can afford the extra computation cost associated with the third-order energy (which is only several times more than its second-order counterpart), one can eschew the delicate choice of the IPEA value in CASPT2, and rely solely on the CASPT3(NOIPEA) excitation energies.
+\alert{Of course, it is worth stressing that the present conclusions are only valid for small- and medium-sized molecules and that the present study does not cover transition metal complexes.}
 
 %%%%%%%%%%%%%%%%%%%%%%
 \section*{Supplementary Material}
diff --git a/Manuscript/fig2a.pdf b/Manuscript/fig2a.pdf
new file mode 100644
index 0000000000000000000000000000000000000000..d28ab14d47fe829182e4625c9ef041af02e1be1d
GIT binary patch
literal 20847
zcma&N19W8D)&?5eb}F`Q+qP}nwr$%^Cmna(vC(nI>agRNKIfkM-yd(hS5>R_o?4iD
z?=ilw<|;BpQE@sZdNvrcwu|}~7*+ryz`@uWhL;z>AZuoC;c5wB`&22z00006aVuL_
zGw08vt&yvlsF{g_sTmAEKa7j3vzd_{jA!nJo=gLql+)&8!!a4{;r#t+ga?+9j3n(j
zwpWN20Jm{gEPWLAy_bc!xvQq;+1Eq!1!_QcxQ<^>FTq#ToXkl7TwojgrJ(WXN#LFP
zZKl;X^~?LC!H3uTFk$-B%!iHTnAiJ8>-WB<k9V`@A6Iol%pXUo{Ki=UJP~FEB3u0J
zol^m7HNK%(hcizvCzt_H&vA5kat)Sxdj#@edYJLYqXHW~B>7akCJ9~$RlWzQNBbmo
zRtUdH2&~?JCEOOg0^?e?Ci%_bl1Sl5X<wN1o^K_X!xI{HD8!7Hm;4503P|O1FL^3H
z;!0+CxffZr)sxQbxG<q=>oYoAPT1MNp8@t=vh)1-<ph?~Lgrg*H3a%I^Y*z3@9m{;
z?U(ll=8xyct5Jzf+sKps)w5(RlgYQofbRXL_lB<fD*^ZTg>)#L3j3_zPPA`8y~keh
zJCpuzNRC~`zgH*LzWCy8Ecd+Q_DpL8w9H(uH#Vi%t=(5DoVP1e-8H<s&z#*G+~bBL
zX-G$(0dLUGi2Qb+v4Q%dN<2?^c}8X5!Q^Wo@?CPk5@z>GlS`&2=aBP1jWUHCCza#~
zXXw%o`Q5zF7;TKpE=ucu(c?(njqF>pasS}pj<fakFrl&T!Rg6_u8IFwlkJy)rw$+e
ztB>k!MlhX%<ihVDVr>MPR72fFzTHP`p$*UXB9V!i<vsz?*0059T#*et>?LEK#B2O|
z?K#Q@wzyYu$`l}UgaXdj?)11pEkn%*QBb`rg+B|mL7GdWNhxh|2w!^>e>#3ESQmJ3
z%K@XS5wSrKeW2q^4k6f>Dn}M3gCj**<&S6kj+#Wq*-FOOb`e2hwuWEjXF7w`Y!@$B
zCr{V9E^}%!oU~-hmsIi~m4nDFD#6Htqcoq9Rg1uHtwTn1wwr)Iob(lAO5y?~g@C)A
zm8*5<+^G7#{(BPgDjwy}*_&^BlG5G~;rmCTl(gQ|=J1Sbsb@_G21&L|?nFM$q!MpP
z4`2afvPB{Ss=g;AVMK}doOoL5M5ENAw<ykg2ERqjb<l;Xjw5U-#EE=ylWb23$GW%!
zcR^PpWdl2Tp>rk^9r4-~y~oJslA5Lc1<j#$-J7)5dhms_BSq~q45I_KIf$*c@y6d(
zx$qjF5_P%n5-++k8<Clp!U{cYQ|~9+d77c1Y68=I@Fl;UBQaG;lzxkfyr8g;8@!-2
zg%ujR6jV>BQw3nWP2_W)1*TVK5WcfKu;Fso`KJf<a{!^JgJizv32x$zamk7*yMDui
z0PclLVgVwBgk+=#A>&HM;4^J0pJ|O#E-v(QBgOZ&l^aiLtU_D_P7GrXAL~*;7%$K|
zJ73l3ey!>uc5*wQj5kCG#k)04)GB;&&P_zJ)g&57LOSd2JbsxzAu|@q1FKUt#=1Yv
zGnN$wIkM$oYCqYx(b2Z}DPY$PDm_49xoeTh*NW*jQAM0dO+TWPRLHfch1gveAf-Lf
zS&x8jaCd&ud7?4s%-L)j=Ro!BtO-%mzEBirjkJWqp64D+L}WE1r5nNbC7#MZ87|Vh
z9L!EtMD}b<^rj&iZC;R+F;EiTs5FN}S;W1)ksX!jq(~S#Qa*@eKPnbM@@RGCcmdNo
z)7lZL=_EaHtd}$9EX6n8qKLl;iOBFWr8R0wa0`FnS6QHDXK7%mysSQ1um6UgO6HuC
zs%SLP92f>j4<m4NHsBVehBy%cLb<G$bl2R+%;lsm&qxye7rTq@K%Jtk%Fw}q<mQj(
zU!%>icZLB}@q#v?wk{g&<3;^LUloZ1<~NfgISAA5K#`(;51YjRb8<1-W7>t$N&DiK
z4l@EdbK=wcwgZ&%+`Vgre}A)oP;dsAbmJ!b*$0Pq5WVk(ZJ`65I}>HrwvJy5;?Cm>
za*<0=Zo>cKGaMuCM}rBxPjuuE7@P?q;=B^Po{>@<jIKzFkbq(n)IMw1^RNPHc{jxd
zF7mQ~68%j|S);&RP>djn@+w>jl8lYQ3O?BN%MfI|Mm|nx5r#@-xmsmUVbUH+n9oRZ
z4y>*QO)Ofu2bG9A9PNpd8=Wc~?PL_13R(bd^BPEa$`DE(Er3#nCG$#&$IL3+b5Tv7
zO&idQXOoN**)4w~8X4XHc~fsiMFZ?}jfhmikQV@?Hzd1B&JPsAXfMbzrB{f~3&mC+
zYF8a{ZndJpsbNH`W=vlga|0D8a8qqCCp<whtPDkVLW4aW)szsV{6rZ0gm0SHJB#Cl
z_nDxK1&onE+_T>Vbtf|Pl=&BjRNeERk=!$@Hh-F>q;gQ6wZO;&0UHQ|Qtl6l*C0fu
z8bK*DrU38?6swU~0)_bxtm9INlL+MZr;t{gAyUa)df|#h5l?)5(s7YD$sq)z2!z>4
zq|-5t#)y&CQ9{2?WpHuUAjAl{lE68oSwMw-^<qJfN<KMEr$w0#8+0x9xP%-4Ot^Yt
z&DWES5i)8F%oGu|zGG%aBS@AqBqs-AM7g$9{uBc$vg;Gnzlv8JV^|V`ekx2Y0Qw?i
zg=FW`lFn8K(I3w-q_IEA6ZTG|A`bN3yTwTkQZ62C10tOEy@Y28BvxN5UMnCtXtFIh
zCu+!h%m^LHa4tLUq<B3R+V?!ErVfN>WW*A~Y&{yZoy(9`+J8{GQD}&9vsI{QLvB6$
zNj-o>h2$KPLNAcd_@b5sf-V2h&jT`y+J_0;C~_;S67HeL+UnSN<Lk?zJjQ*!9EWY4
z919i|_IF8UGY+ocg{}e!*wk7nV3tT<zulrWE>^(+L{{L31%;r6FPY42He+$(l)5QL
z-R`>-=oDbk-|@y_Kmk@$l~F<4*$C2RniNXWTvh~zv87YMuIz$aL~2!e`Yp`v@IQbg
zO<O28)x)zlZC4z0+Nk8z24FZ$h@^zDfOfloS|?xuU3?zINxAULZ{_+eEyFA#)%y>D
z%x7HRtK+Pwp;-{(2ueU3O9o=jBC~~=i+F&c`8A(n%2j_Zfd&e)>1mP*%i;><<LsO+
z17XdTAIQea#jl!cYcLNKsLGR47)d)I2&UsSO2P^Ht+DHn<?~m4{avJ}dP5F3U}NVS
zCV<mqW}!lt1qZ#3R8p>mz$FLK4-{6-#cCEVI^Swvr-eI;Xe~CxC<i&Y2DIG`Avl3y
z8+$I;tOD&<WF#iVaSUrEx)91}psNSfFCoNXIS3c{ovsm;wLWdWHPw1W{E4<~>^lJH
z#JYLyUPZU+3O=z66PDZuWQQGFwWcl@T62X(X}SYi;mM0hr8eE=1sAAEd*|XQlx$-5
ztIIgJU<htvX{lU|aBz0(O^7v+yAD5D4wV)vehQ=rm>{?MObAe^3I7x=C;)$#;diXn
zXz-QfjW+zi8n9e33%Hn_ijWV!oyr0h@<DiMJn*#)5#=(BsiQ=ACl$79wNYlny~K5!
zOq-o-+z|hajGiATz-Xap4b#{L>Pt~t3p=Sj_o)-pFk%#JjLty_Fp?Hyr<*K%AauTb
zj};K@HCTlPs2-Ioyd%pHkPi6jy!M7*v0*q>kIc8hE=DRi741pG{z3QTI6PB2ky$+S
zK*ezMIz9-O14f?yM72@S0VN*HLaqCzES}YMrNb~_#;-;yGKV;Z#2f<wxnJR_E{dgr
z4zwm-0~K%>kTy<CGYofZeZQ%-1Cx(Z0E$$>AUc3!LYvok3Q>rze9cJ2_-tLFH9-L<
zz@~)S2LpkksclFgVtj$?BPiu6qyv%C-o>69?AN0G5jbOZz^iQWVJ$H$(&~dF(Rfun
zh(cXBB9JMhU?V)e@+z&Ay+K!q=z;tLXkkOZ$EdQzS{gv?KDjMlm+5C2HK9L}@<I{5
z7EEm*r7ZKeu0sg!l8TO*=?4a}|450>83Pnv012r{Nni)M$q!qMtzUyGSCIlbEf{s9
zFs>=41S5GhB0(Gji<sFeVmJvmUg?)HlQ>RN8)Ar6tV*mP3l$5?>g*vEr^|;byF)-{
zX`owuoB)OiRX*j(QpBWSf<xFl+u%Il#>6WE8wNxYq&b+9R*?jA6oQqjR59@|Xwozt
z`nB*>!-hkLKzPv17|BuQC+jVY9i%gL@NrP`niGD=Jb4P3lR!CU&CP9~Bv?C9uN<Pk
zL#WfdS%3(I!d+Z|sP|MKWD0?jgle4V79(eopB#a3kj<MJZiM496}AU}IElV+*P7cc
zfCnkf5bV?!F=KY<4CsoOKpcJ4%7tO;5o_l>Gb?9iCV(%mL+A?5?rm80UM&`&1uT!O
zr&p)e#*K+l**oAhUNIJ@Eeilwt581b+0DvbqaZq{bYBvfbE$jNHnPF=@PiiG!moaj
zLG+qDdy)#Pc8}x6UjAm^Ue0=z#Fc18H{BlG&`#G@Z*QiJtG#s2`cW)uZe8aMdC)d>
zQSF%ebm2-1bzG{2d$SZfT+4B~WNt7K&4aem$ljL9wwoUqlKkjIS=y|Y=f|jQBDc0I
zQRw`m5ZS%0NJpr!?%hM@mlX|tQ@?#q>AE3wQ;*<U5gqRQ`fQE4yA4CgioB)scJM+t
zVar9b+I8SLq7IvzJAs*E{DI;_hz~k1nFi_;hNCrx1yW*QQ)7?I%@kpA0_r9vXcasn
z1!aFAJAUqQZDE^TUH|69?J5Ts=eQaUl*6n4UT!1&3K^H5u^cx~no%3O9u{aPMq%&n
z3)=RN7is)1%<Z=qYK_X@hhXsbFs>frK7$5@f{#laF>NQ*L>GW6_y_E(vf^zrp@*re
z5<R1QXwe>ifQr{;&Lq62dpXapgm=KIe<t}Qy0<q^kNmHfhZDtI2uqj2UopOdu90m;
z?2jTrG_s{9vZ17uMRw&y%pMS5$$@n!d(Ul!w+%H&*^2O*tnSW&x{`FPHE>Y{h%yJK
zV2$y5J+Ei9Dvu7LYTy;lp)S>k{EiJhhQ$&UTF2cia9<c_+m6n)QDH8SqCT`d=;ZN!
zfL`yFPdfrMss%v-?eWI~5t2;PERXg0(3UKs<(#W>L=NBTPq$-AK|wJ&6<^QfrM^mC
z33oN#N{g9KW%k4+nN5JO7mWaFLF%wU3r7T9W;wtr2mqU|mK}rGzzoGluxePs-e|^v
z&6Fu0vY1-<zWr;*jdUw!9Ti$`JnIGh<-_LSXbM;n7x$9$0?7r%LOYl=Z5?Td;t5gM
z;;hY~!UH)wqFX6}#LBkXxM9|$Ua?aTFeNJRcF8)o-{hbMeUUcGz5N(@qBZHkz%s!1
z?X+(GylW0sI#w!`z=59gbev@nVoBgTM^=co+JnPf?FZv=g<yg;49#`Ws?%}<@<_h{
zOPEZIhF2oFrZ7r)JZ6HEQG_i7^5f`o?3FK+#tbQQAT9||5VE{>b#5LvcDG{!nh;=f
z89gA(YnHRleM)&^YyGO<0S5zk{jv{_wJ-zy{uRI=5BQ=0TmNi>ns&@r>^I5}DNHIQ
zw>qt&8^!$?2*yPe_Qwtl>}^wTaMEg$u1Iu7Ki2V@(lVG}z1U_;q2n95QkD`QAS?-B
zrJkBNMhIWqtoE@-7~*lWsCs6=c=p1)gn7B88?^EDxuLy@d&4jgr-C6;XI0%4r)2G~
zR5v9ouu?DquvqS50K}LzM<5dt3T~^|WDyqB7*af5D6}x3X`Wb&<`V*HECYp$?(OQ!
zuKm0KO?m;62M_X96hP}Xw`Ed}XQh$GoxLCcK(kT{HhcP(O17lOWrIaVE;W09C4EK<
zN=%6UHb%u3IG3Ed*$%2Q)X0zj1pxLOoJ16D=y)YXB^eYiqr^^IqTI?DoyrEP(QQsg
zBX1nJO!%G0#J45~9)gh_L7^Zb_TfxO@C(S!?hq4l5;b!U=3XRksXktn&Y0?aC*_1K
zQ3f`?AsSiI;ofX;G3AM)5F8G4kxM79c8nl8>P_&7a2iOo+Q4J;ShG=N3)ENyNN#Z}
z%zk#MUE6BRAPod%JX|wBWpec*IB2wSUGsSJkyp^|E<&(C5<B8ranA}I{S*E&x!HK5
zgP65^20LV?PuUoyGx%Qq+`~b*4hGebZf;uhVB+9T0GsVDjMpAKqN+Y$ZGN8l_ZkWh
zGmq*_m}I_4c^#gH{Ef7<UxE>};(IMXWlm-sfzDIX5mOaU)ZfJoF@UUsYIXI<VH7n@
ze@(H)BnJv#hBfB~vXPT5B%+7`0X*5^!J*!!cD&%F0*!BKCi!F^>(QdISS{vAZfIQY
zpcJ)ot(Y}c70{GSgR>m-QHO?t{e+}OY7_M$l^fMPNb=T51Z$WW@{`GG!U92l)AU-S
zG<)8}*xDuf>CQ~OVTHnH@-}k)_~z0|9X^!FXbF}}0%&e>#Ja~N#+yIP(qQA4e-O6w
zcHsOKQ;%;-C6nYM+$qtYMT?+F9}>Znx<olH=R69m&1t;E%4nt-C_2xHb}A#?4(J!T
zH!x{7O(}VjFZHER?nR~Aj!L}}o$w(3$&+-GC;30rGyh%Bef|F{>c0p5=NKOeqSC63
zw6%s&Q%&AQr`ZV_K5hwouXop%6?1NcjJ|t}&>dFU;Z?r*XZB3bOWxWu_Tv$9a}GH#
zBiZ3okH$XkxSp$f3KfMsltW71oVHU=x#1))o|T00eiT^ZYrg>c5WtI3l{rW`5JHc_
zXzqi}8RZH$?nvW#vzKMha_0&m)}E>aCGOaZGnGY`+_fuCh|<^d$tn#J1!xx`<OHhs
zh2MGH0$3mjWW$y-GrdFh!=#w8=<>8!wD}Q7<_!7WAeN`VTuN%077i&y>>2L~Vh6VL
z@nk0;OIlcD9KvgIn&X^!1akWFf_O7$71SRLBR$HW4Zm-H_wxY4$uRVH#DVD88~=68
zmon{k;+F$?RN<6x%-BY-_<C&QMz48Q!sq3bSHjU5j5uw)FYl~dH_9d8SsF6Q&DmD6
z{j!7OF@3wEzm4aJaDG@iGh^?2b+jE7m_&m$9vWnebsQPq2ncf`VLP-~Hvx!qpPHI#
z5I|)H6`xS%#_aBP+T+cVNG!TYdc<1+ioD~iNg6L-)g=Og@N*8+@z2&;^7Xot96app
zyIMSMQNMT`y>nYP5+USC)_&c49pm288t1WCL&}zS)Q>%S?T>tjJrS(5-wDF)y&oBU
z@Vx*n*qHB0xgEg7-@F~>E8v}-qQJg5xEofVL9~4kbbr+S9k1IQLwmJ)2ie|3vaeI`
zUv2v?cKB%ffdzkD0-ngRfG*!TikrWO?Dw@e&O(;px^oZTgtrtn<8CgMW_=LaVUX$8
zq_sTZaI=ID4nj_YwK4{SPc)^=aDfhkjY`vPq&i$uo$#n@bYPuxjdev(8&~3N-}aGH
z7qeUi&B1!eQzO^G9|m-*BvcR0$s3BL-D!+Ter=R#gY$^BAZeBcX<b6PThN-&!VJ&7
z3*0zI@<N2vBleCxX{Tx=I!DooP=>%{%8yoiT+gEM8Dq*n*8Oqrnp~TI$@3Xqo)F|R
zf#`|Iaz-awPPY&KZD!^Ru*gk+n!aKRR^dQ4Tqu)-k)$!!JRVY(0}9c)YSfb$vsi(h
zQ<;gSuI!p63;RfJJpIp-8K_;pVsZ<xOJO-mFG&t;6y^2Qv=geDvT>g*1JV(=DweNt
zi>51#ALUYgUXrsXES~jf1xKWX>~fa=Ns_KJtlVfGR?bP8^{zAOI4IfSsTDL-r(6!F
z>Dp|4lD((%Ey2lNIdJskUO!yEKGQZ5$(?^E@kz2;CNG&{m|yA&W$%d1$>k^?{?W1*
z{_P8WTfJPbPga2j>hD8LXn#EzfD2d!xm2Q6=p*eQ(Ojd9n5{FqfjNslbOOR1c9P!2
zK{;3=c@2tJ#$)Rw5mSrI58D(TgR%g9=wbv{s#3prsuaZv9}H@=b2<e#NwXe`?)J4t
zhFi0SRBaHE`3n(kdcN;Ehx8+!We*7BeInkDIfAdmqn&ah?uL@8mn>OQ%N$>(#n$vR
z*j)=4-Vzo4%5?k2ChTw8d-zD*J*u|5bH9-D<dSdF`T0sRUMbIRklt%poof_#gNjhU
z(rM&>xl>7zZE$^I-|&|l3y@svOMWg`c;AQtMsU@D44(v-@?X7p2fhqr<3J0Sfg56r
zcp&~3(~WZB1_K!mi#e1Q3wrM6ViSVNv+jeBs3FibNpx|al))ieFGC3wZbQ8<hN174
zwcbr%#Lc;4mw7oT>rqtR<9K@`Kka50(eK!Fw0k-OjK8vpuy<c{_Au+&UD@ZisH-rz
z)cFfS#b-Q9FaY*+Hycf$`HR%>r22(k%`W${qft%I{Qb@HiW}4SpvzOufR6L)9XX$u
z*nyn0;C0WKi(loqyy;av2o*do7inK2Y~0?ax;HMqjuSOmU!eEC=WlyrUb`DOI!9`Z
z=E#KYPYqRUL`{K<Z>8ZeaJgd=6<kWYKT=g3@nl@|e4iQZ&Ouz^@n~vXL7&~=Q<4a;
zDSfkBco054fD)nZkg54GwfyVh>@EEKVAe6_!F}%Ggx=o6;PFY>!Q=Mvm%F1!$nPb#
z{SCc?(U*xId51?aGe0hlF*zz7I%*3Vj6ELP2(Mf{9(xV+5c+-q0=}AU{SYn4cc|&$
zSQj94=y~MWuQ)opy}*B&sqpzdq~`bIP^0pqYdS#h$)m6VdR7D85*AiA-V|8lBB+R0
z$Y7wb!3imCenVBL6I0rgzHC-Qr&G<eQVb`%3qJm?4S(bru>%MBxk!f$TH@InRYWM-
zX9`A=zrl76qusOGe@<hpQ)4cZe?}YkQ5jc{FTTiB!eGgyx=rBF=Om$*e@s!>l#&_6
zjI|?7vdfx`fjaw6O+mSC)GbZ?ONN3MbtCh#;g=<=sxqByYkJX<Xa)r}qwFpwl{z|>
zjxw|G8x{q0>q|9T0+}Co_l_gcTqbh(%RiJ3$J9J4>ZWg>GzX5VCJ;$@<){kl^=FS_
zvV~246wxgbR!ge%MUH9iA6Frupbafz$aJZxI-AyZEux!uE!k(|O=yu7PCv`*L01+@
zsv1aGdsDM7sg-)uYqO=cabrz5NE@{MY{EOHv4j_Ew|DC#kD%PlG1fB?YltW|u+w9N
z!cjxa+1N8O@oQ=VS?iHr(;T<0b$`%As%yv^>9Stf#85x0Zs1e5#I5OUn$JmgwXW)@
zC)kY1cJkBvLQ7CZ$G@Uh>7z7UPvfd#Jzmm+u~nzP^;%rIrRmmC^}T}Fr+RxO{mQ5M
z7soQNtu^>*1iJ=>UEDy^>5gVrW0k&UfF#Z#wG(HmkZ*PBs=7+<basr3j<Zwk<ZspP
zrn)Q6JT<^kIn;ql(O7TtLC^BRBkK9ZlY>fb087@`udJ%9)HdUdrMLK8zG#ChP97*h
zFPV?28VN`?zx7An1GnXe)7S=h{)MI9u3SMUjEbi|=8sWA^4m$-yTZ4;cl2y`1|;3X
z=+!>5jTV7dx98!4UY%Z@L&lHmiW7_Nx7Uci9l15aZ+$^hTk+t}g#nDYh~m7z&j_<W
zLLNlj>#5$gg3eO<Rq^c#P&=PZe;(@k$?u%1@IwB~?t#(uOI&X4xZ(~f;P#7*%ra4V
z-+I!?qDMBo_?q80J3}b@G4KBft-Bo>8a?~x#1DyVTY2FKte$_YM9S{+7c}dJgB68B
z9VW|)T|wB#%kx6U`+^Yo7$>iR2k}&(Alaz<BIf*^G|Vwuk9;Ygj&z(IqJ~(uKA8$0
zN00nLg&g?WL#v?t8@`TE{PP8_d?|L`lMb`Oj6TKuo$&Ic^Gi6Yl{XQgQO`YhQVk|6
z&jL*Yi$Ae|9lzk&Fxz=Rh2@E;kB(+1GE8!z+O@ZfoA~eKiZ&6um9YlD0-0JBcb}K%
zIO+!Lg!Nu;RVy{0#fPPmiE5X^txr+vkH|Jpjr+&yUm_je64B8>phmw)(B^kvsa*6>
zd0{S}VV-HLYr;Eo{qzXKI<tc6`%opgP?5jgq5!_$cE8|5R}U&z;O#!ck8U}lJF~V{
z0Wx*9sN(*&SJS6gHvfHSo*%X2J=e_27c!dJ^L1|w`jnxdGml;;@&R{RSUJ%|<fmVf
z@6V%zA1@NsUHUT2=lLLcz3)#la_=uQ`R8zOvv8ijnHvn=Mgx5O8C-fXlijvGuVx6h
zVhU3SfeC#+y0+iWHTquPUku(BX8K~xp4J=!p3VYh`v#Eb`_?@DKZYCt1)t?vfl7br
z78o;o(|=WmKhOTumH!lz{~SNd$sFulTz?$?h%EmuJO5Q^7V~tKP;vb%G6Q&d|0q6R
z3`UHf0)Rn82*60s3SclX`s?^ve*VYK^#5UJ{x3V*zv9^d{}IpfA9i+*e>^_V{=<j$
zzkFE!<-ziwcD8@)On+IK{$=~S2L^Emd)H5ge`ah>$4JKrVE?DQ{%4WqbWEQrj(=32
zai5=m<%s#;IRY@KdO4Z_7!-{x%m8EnroVEgXyk16**()gu?%u%rdCG64xRv=PZcA8
znTZ|1#K{WyL@A#GiF{7w)7=H|r{Mp0bzj2S!OihMruX+4e`)`jvC8Mj&PMhwj(^57
z@%l$C1#oqCGy6vp`ScPsbGI@vQ<f0^qx{$EmCaln+?-9!TmXNrw(@_hjPq0eADa%i
z{ynY#YtzA>75u$BfeFCM!OZ#JcOhsP*s82x3miQw<bF88BMjj39F7DHuaMhFK^a<}
z+EA2HEkYZ0W3J!h@Pq?n2PtQ}mrNaVD=XxwE=E*qS4f_fTK(MJG?k2%$3avZE|ssQ
zzArn|r!KN%NC{_G)6uZ#yVBo)Ph_-PpDcZv;P1Qg?Q@xaJj%V&WtqmZF>d)jBx*!0
z@U|PblE5-@gsp;rF`bXxT`NFf{QVZOQMv4dM8zWJXqx#D-?Wy^_u$(?nZioIZ&KHc
z>$l0dabKsmIF=)9*h2s!p4e!Y>Mc#;5x|2mkDFzoJ>mn7`4^9TF73e@?m5<hK{C3^
zQ<hb4j~&CoS>o{B$@ocq?mX}lVxoaD%(-d4Rsi%edt<*@Uj|DVgTDbIPOMnSoEGAX
zt-2BVA;a56o0dV8lY3CWPM!t<*mhqf;}+)8+<T?R`@*(Q_isoALuE(lG>%m0?Gjoj
zc@_$Jk;n4<on<>vPwt+9_~}OGIkhiddx}<fzVs_vU%VY1ey%ie2=tD;P&*jp%3r7o
zJe_-_t!($5Pj|rH4%ghad0O?xL-aYIn!~z(YY*f3Mfpo)lqTTv`4@l7jXGzBkH4#O
zl=@cp!>^ns)Qs&TaIrVhA6ePGzu8?Z?ieSOVQ;BlX6+5}kF88*n#g(!j-KEG(hXF%
zl=xBjiWSb0GKz6C(kcCI^AlWPpIy}apq?+>-cAYQ+4%Pwei7af?q)UcU+sLEZfNM&
zFjTNR_z?KmC@9!Dv3To^z`o79``UF^5sjSP&9%$Z?mFZ#;Zj8|Qw}|jfTg7pIY(1#
zXRG@BNXUAy9Jy($O{T<eC`p$7Ar=U(SyV_6FItkEm3(4t4%{u;F%Lb?MpdgUOr1>$
z<3&f&!v~4UeoM}$GbSCd94W2v9Z$eibjh;R7AkA9$Yv}Sic8XRmNo36`Bv-;ww{{Q
z82|hu@~!qXQe$~8W;2zA<}q@q4Kb^rWq-70xs^C+LQXlg)^sN8ROp>0$gQV|m*(Mu
zOZZ)t)rrQ!^-$`jv}8m1AiYwnU&KN%tEtnswm^TO35?PT4wy3fiQ-uBbAAIj`FLz~
z`eZ^5_@&rvVT22g!k9#xrF+!*ig5_O-oTWfq7#Tpay5~i+0!yV9i?tSFcN#WLTzWk
zves)c>pj<ugjmLiCP|~yoMb;CC;XJ`Ou0|V{VIob$6@h!Uf8HnFhU<)5?zhgL<}ke
zMp33Ta!?W*O%gs^fm4R1zFFt_2qnhTX4}sYZPwAEutRdhWPGI(a|c!;9X1`M=H|>9
zQqvrWVx(rVfXHyrN-LH2WLAVKbY<aolWYVEQV?WHL%p?fonl+)weUhR(+GBMqtz|r
zk@NJrlv7z@h_*Q<6W6{C?5JXagS+jLwOTy^!_|umNke|Rjr~@Orj!MD*tga_i5o`d
z-Kd{%Il+()3uXry=@)U;Sn5N_8VVwgXqBU(6eNuSU_Bk>)@fktv0<vBG8wjyAV-_y
zYLX`GD$S_jwi*k>G&xOA*^P~97#Qh&cyP|Th$pUV%~!uhf5!Y86A!h5_k|rsb*?ja
zu^x(KY0JyLD00s;_o=#=w<T;lW7|-$6K%n*d<!;*tb8jr=TLqzldBund!`fBM3)*T
zECEeU?3qm6?Yc>=iDgiESWK(X`w0;@zWh62XFFqO8)B7ghp<iQx5S0`Z;3jXkIIYi
zfW?td72kEjy4+`*x=zrRL#49FA!nt<Z4H$tKbEI<x?^3hk&+|dV-8a31cx*ftxXnV
z6R1=c`$|J1y!E)GGCAZd({QnG(*gm5&GkF&XVMoHzXtMWe4i`O1dA(OunL#q;i`hL
zfJ@bn@Ahu9y%`N;NA1Kwjh<#9PSiIc&}$83JSxy!*A11c1K_hAL=?1ms8USA6ogHi
zD)b4JB$Vh7^-!wnsUg7AEA>K>kV?ga&7iCd)#l%e#>M&gwdAXXOD8fNUh2@2T#u6V
z?*gex#l%c3ziX<7N{C{l5zKhA$c3>6X4i7tx?xFZ_Pz0HtA5|~5t*wirsnhRg=#B4
zBI(CKZ{8R^Fv6JV=8iS3)-=L;FLF{J@TH+G?jl==p$;Pb-Pq#m@jV)Qyd={E%-_Yn
z>eZL-z+Zw60i>rl-x_ZbO}MqCJ{ed-pf-NCIJHNW&-Xx<ZP0d5G8(O3x4CznYOG!x
z=9Nj(Cw=?U{ziR$c}cdawZnt1MYN%8qF5tvL$AKfg9An^xtNod$CU(gAdO1dd|v*B
zjkEA^yBI&=z<v&)j`P8ugkbi-@YNo+ob3yrez;w0-_BySDk&>f3RV#h{lZr5&(QX=
z=%M6qTrL(<g_k=pRmdN4a=bU&>m3K+X_cwb)UCD=qEppzT8qmAo^%L9iS45$)@&$h
z3a9c3{U>dgz@um2)6SGs9ecb2ozfo!vZE`ZrUKU@*_+|!se9>Sfy0cK^V+2Z%a*W{
zO8q_SAkpak{6O+M$aX<Id|W^k7O|XKflhl!s|t7^sNaYj1$9K1e1MK?1UH2Y7BIeS
zSO9(TBmc!&6!OhwmS~UbSn0DdX{$i~=B|_4aneo+6th7YL$$!#VRV8*ixer8!~;o*
z*Rv{2Ho+_+iPsB_=XOVa4*wKP{s=to$lSZDVwNyWTS^K$9Qaa~5WFXYFTfR4$MLLS
zCusF^)JLtbQMaK(yK2+ea9dgh_fith;<;`_ah?-lq6GuIYJ&^(q=PJ<t)pLV>FA}M
zgzL_Qp7gopv$<etk$vw<gxYo@LY@xkuHoOljX#*}I}<MM+Yv4fKqRRNh5Zc7g^ju{
zZeKlCvu3C{>4cM&R`9SdcSoeas(DwjxO!*3JKevi>T9kKbZL4&OzPMwd+6xP_HJMf
zxK*u0;f9<G;qZM7aq6w0#$LvUQ9b;qkD&H=&3N#}&bhZU#!Su7g5K6U)*1C)<Mb)-
zs?T^?^Ny(5A$-Mn@GDfQGpc|Ln6FRRk4yUi>ZsM}tgu^7)IZ*<&FD9l353)|d*Go=
z>^;RE{1(?AcT!Uu%|~nEar=d7F@bC70gc5SQAz`}CDte65XN%_0{o-u`*1su!%5Fz
zhJNB!hU2XThM$<gq^aSw0I<SBO2q$0=YLSwztA}|3oFMz$dl<0e*7N*jp@ID=KqDz
zJ~7|Fw}t8{cH4ggAe}!&(=iM~Pg`A*no@PK(%&fN`T~Uy3q=FhIqmi#CM<|a4=_(o
z$A<r0+4Z`n2M(wxl#n&d;|R4`l_Zhn<e8Dj$EsMUt5IX1X3ihp?<LoTZ~#p~`7J?E
zALY);o5i>b3JUZt4l(|OjbO0S#22R%E6`b}0MOhH5z4vHSM^nc1!>sZH@w*_d&@Kn
zh@Ighc5M@AljO$IEt!{`m>DZY;*+%T3Ez6n;bR{*z3(6<=Sj&DKA_om89EAPfH$Zf
zjnmtl=(u9K7f(f9{;_CDmU`=HU;x)tdHMwIa_(~YzZ3dbmj0blwtpw|e`EPS+(6XK
z#l+dl(bd8E4{ZL6F320%eFA+IF*!{!Wg0mnS4%THBUdXEqklVy8`)XedI2c^!$$SD
z(bmY~6KXU6+f(>&E<ne`!33aVV&`H7Ftc)TeSWd9|BDnzxf<D8nF!fi*qQ+t|B4fG
zG5JF(*g0AL68}u_k4*PTCjL^07&%IsSy@=R{?*N2Div2VJ9Pl(-|>G0$$x}DOZ>c7
z%$Nax0RMjyfIl4QpZn`?+VRKt@9x};UH>9Gs-G0&kL@q5u#t<|Uu^3CH8obwF0LY$
zM$Z3aK-TDAGAr9(J^c541DO6?k-vG6nEmI*W-EIO0E32=y^y_&)jx{A;xw#GT`gVy
z#i{;07yrCjm>ECE`bT1C{TI{vWK0ff_EvxPN(26Sj{bAqO#j7u{~<I0rvGAO|JQ@{
zpVWtoi;exiQy(uYPgS+{<8_Nl7kRb3Y-On_F10GvGU~Ka6l%~UI0`8;_)x)bV5m%)
zKth5^LO4QMsw0({IO|18>b4Q~6%u4{9gYXkk?pJNkJ0N0KM*Q;$u4?+1YW^Zd7tIN
z6h1n=JUsS3{xEnv&+S>(x5?{S(f6tAS*ZsC2BND4$NAcVaYy&u_&ehTxLF#w(#D)l
z=hV+p-;Yi3;RlfHWv<u?!|Ini0{FmvY#{fAdye`{7T*2!0wcuL_S&3{ZH|{59*1mD
zKpIJ%Gd;?8lJ%z5UgV7@7ff0`{w&6`#vb&IeFt668is5gPV|i{6Q0gOhbQ9RoENPA
z7`bbD?fJgrnT7^`L?u#71YmG}yhs0b^tPXcRt#d~joGAG<lRBL<g<llS<tL<T;T8=
zL4hHdy}#OI+a!on6joZ}7G>G<>a1Ae76tvLB;!t4)y>M7T4r<u8xr)=c^~;Up5J;&
z-*%sX;^hl50{uhs;OQkNQa~`ANj60f&E{IFUScicEb8EF7BySqfN4;7>nbsucv`uO
z7v@<xt-6QK8_!xpmo(WpuvTr)OyseB>G^ubUxW9(q9GA+GI8%ZOGf20r}GKIy)7=7
zx>4P#3E6dQm%iCnwX4!(p13WX;n_G7P|?_F6Ln*w67GoFfkSEc#Iwz#i_Y3mLzxW1
z5F?{7EoE90CvA-H5gcpj4yQ!l6bd9PID)&8vXfJ6NF#4&PEg}s)R1u?Pt}Lz@RTHt
zWjb&OJiW{)<Bs9>q)g_*&IDMK(2So`1ibH@SCb89_{3j%(%!qDwkO;BodyJ1=giEp
z1G_^U^b4;weeW8O2Cqv~mlj}7`5;Yr@5tJu@IRw4PTpiqWn0K5tubq<te^^?G53;p
z=vd1|%Ofp&x$~YgYUvY94@uya_cnB`b4()VcpUBfmR2o++M0*WuP`l2Y*rT7GBM7s
z<6^tk8_lh3q)uM3OAW$gO`g?dd9WHC^^ka5?c1z4<K8m(ManKLmPbRy3;s><1S?_u
z6Aj|io)~#2UVwuh46<`wVc*2P!_xSHZG|xe(ucsbYnM{|%GUc*York`Y1d3j<)5I@
zJ)ZRq{Z@L;q#{Dj!f81CEMkcb8<xOtmY@;hDUin!hEZrP{0wwW6?~UBvPt}#D}o4k
zB5f@gCOx{{2sMfl#`%1TR;zE!m#=|0M|0V>HFOl_O}oWTsSFJOOLjH2kxDg=r_+4z
z>a7e43VwqHgU5wM3W;UgY=0eYK^&<yg8*^3s)m5~j2VuyRi_Vc%L7Z|Hiwx>Q{JL{
zT%|+%Ax{42U#P*1uj;);eTG^C&mim9+0$l=C>PQpN5B=6$AEZa#0yXj%G8<Py1ru$
ze*K3aJWZH7eCnVVi8tUrg@Q6ikhht%pbl`C(=^c7GM+Z-NfaC$Xopouga``=FUmz~
z-DtrVLe!b1tuMqFRNrT%q6uQrim}rec{g01?@n-J+|W+X&<kh7<%Whz-VNfC6zc@Z
z(E44=Bj5_6MVTF2wbZ4(tC@Qh6?Z3ZGLZTh6jJQeH+jv~bB-CA%EVqG?VMEhQN#}?
z@?0Lvd+Ddi_b>4ZRJ5&e10+r*HwHQk@|1tLTxy-J2ljeX^_U5v63Y5DR{cPgR+H-`
zL;qSvH@ti(yPf9^B3C&Q6ex~g>QTbRF`j~Olw8eA90vUZ?6Q$Vpl=-i=GjfZT>os@
zj^(G+7gtgNfrO&0#-)^>dS^uX3RMv|n`cbutDTOL-OEzBNwvYRlsHkPkC)<Cf;{(y
zMe>;lkCbz>?t81(djzZ!ZssY9ami5#^C7?HT;?#Wbr+MWC#spa2$*qRaRTxO6^SdU
zHz;&<P$?IlX_qGLOsF}vd&T0%&4LVd;sfM;6sRXhGul!_TDY}FEa=dh&$Sm!eze%g
z44UX*tdv?;uqp8mL$Ib+2$UHV6@cKc4g?{q-Bn+Q3I4En7MrGDSm3f~wMEdx&Z$$`
zFfwk~xnP1E)-wJ&orLdkyh<u|;eov`zR8t>hV9BV$KcJ_OS?3rME1>^nMoatovDyw
z1h9I14q$?fG;@l??k6DSf>S52BZE?hSI^ixb>C!^OhBbT?7V%tEd=3&qk;7n5###V
zKuyT?Ol>sX=g#Ps?rOmFl0+Aye!hxPcuApEb$dl9=O>suIxg~a;e$1m`=Q%FKxSl-
z@)o{KVV=dYS8Hd)?iJB_hS$H>zv$=4`9>%DYu&TjAzQlIvfA1ayhWTbg+Rm=8=-l;
z_LfmE8b8ZT+t<RjH7g+*Xl?C~4Ia7Uf?1pnpLB|tjq<0o@Ojoq?gTzDHg8!6BiZ7h
z2g$n#@?maed-Gc%NY1+K3doC)fhO8@6T)$@7#!B}6}1v{9Nm6%$XubC;dEFoX{aNV
zh&;Agp)sazWC~3taXBPFG8a#(0PDs@G|8^~*(90t&MMru0tvNL9b+6`Dc^zM!oKM2
zr|B{83mK2^8`@`x0nf*<4)0GtEv&yn+=Nc6xw|mEEf(=nY=5aT=sWb=<G<z}FPJsp
za&y;O0gJ^pcqZ+J2bBy){$kZ@(woxEW$$F`h&RPY-IMcDf9JRSs4g^yQX5Yj!Ip~N
zZ0#wu)5@I(gethx9DR{qpb3_@B<rj<FP)Jq>yS7JRxxDb!5K{(T<F4;XB3({d9s5g
zWif-I>z>!&Bu_O7lDa^I<@YV(WiB7GMN&N(Z}qlNFqAZlHfGD+s!ntq1vp9)k`nmq
z+yWmTpZuugjAX)_3>%n@wY4~+B^d-qdzvU*jjY6K9GN?M6z4@b3uCads2(;CQy%s(
zdo~=~L<~=Q^IVE;++HjfJMd6T@mImMUNjVm*V0flRmC{FOXdKoq$q(&oJl+K1a`@}
z<V8-Z;c+|EwPgUZtXaIxVS?(MY%Q6EZhJn=Gq}Y-vx`j@myWikz~y7p9ibK~5A16q
z;>_5zKxbd8pM6F5#qD{g*wLzHhI!lfZeP3%Qg~+q_DYFhu3V;3^SO+J7b}@4nX#FS
z&fW++hK;S}xpvaBHB?r2NvUXjI0#h-+523@4hJ&zDQrwIiwmiWw>w5AE}Mu8d07W>
zzWIx3)T#VH21|F`8o}>F6R9>;AR!tCV&WkW!M9aB?&$_j?YxXR0CtE+X597Bt0}y&
zUf!+`n~(}!-0|peTn)kihI_gses52ajkW>L50-|ChrZ^c;gZsP6AL!)TjifQJOt8w
zDc&6SM~gvHNi`u0IQGv+BBTYPTb%Jsc`;r;_bn-4R1^rSk_w9BTE~TWuwv$nou%j<
zLgN{Y<k<A4RBfzFaU3ez7u-4gI9j5%0h?4?R-Fr;6({y5_E!krkS%(f_N(@r_Ul9!
z*N3CSok4q^`k?K(c*1KJzq#>3o#y(MGj+x1%%H6NI0NT$_6yHh4QuHrzV7^F%O-|}
zfzh^G4K}dGbB$Nu@v9BubZ*5Nxg_IN6l3>EcAfARdVuD=V8gN8+6E|UGbWM<j(O%H
z#M$R!^2PD<8f!3D;|z!8NyoD(o7A%#a7xVPao|W3({?cmcSnn=QZ&@J(=fsekw(|Y
zvGu|<&wZ=sjwR3=%Ce^CVxL1U#^WUwcOJLQS0{)?DF~v3j-$(&-3afKpUP$~%izT)
zuF|~tope=C)uGKTA>S?&HvOwKkXi93$*A}{3euskIATG!nN=D)(vIEmh*X$_{*`;W
z&V4ozJx0B{#WSpTW*%&sqFAx>?72Zvj~quK?rcEX=kKyMKfBxVfTb&bg5<R2495i_
zwq1dBPLS@JHW{~V>ajP(Qu~xvek)TXJ+k)Fp3cRiK1woFrvzAIi&Ixk6lW!)2+@0}
zk&g4OLDfsomBJr{>FTBp{S;m2@1R*(3j2kc5~oJ5Vs)?~iIZ2>)l9=6rT$KmDJ_>;
z<!yVn%EMO^uCOY1byvu=p$@;Hu#q>^8B-I5kbgvO%PCO;l>lyWQAI!&N>3b<QHiX2
zNoP*?6#i8nX;7EKD5mK~yYYkrf}2Rd5a~(SjsiLYtgbxr1gt>oRFxj$v`w&%9zsh!
zo=$=f3f9q$@jLYd4U-fWJdbS7%_5j8F;}}x7U)sMc-Kw3&g0Z3Gr#x8CDmL_z@zFU
zvdHWZCPqfw)9;aZwX(7|p90(XkNe@XE#^a>%suTJ@HH+*x5<gH77B3-S}fGFmy?h_
zo9V|ga{M)}8k>pL)2=X=hD+`Dv6uZAJW~cg+7t~cFKfbl<frRW-zl?&g{Ym0kG=(!
z@uQoF*DTVsbmv<pzMjmCXD{mQ8J_cP*t)WH_lpf}?Hbauff%W(t|GJKd3IUprHxTb
zacw&SIdLXs;twI5=K_i%EckZxeZ~*lz3P?>{yc8%_N3Y2*W0(lGLi<cZ_oF4bXnhS
z^VL0m8}Csu)MJ8R+LmxdwL|w$-h5^73mqa+*mpm+y=<)4UfIh(?A?0P5voAH9;0pU
z)U~ezgwlQalCPVAaWzGJG=Bz`57HOc-sQZR{|b~dHECZ!!BdKu<sv-&n)odOgh*D|
zm4G8Zo8Ds8R9^XeZ?J9^r3tlbmXZopccil%pWB2h49lZRg>}jFO0Iil$5QO-;sP^C
zzjHZMhz__-bf%|SML+6w@-W!m<jC7DBnJ;K*wm@}Zc82Zc60Kk3_n0(3{(aOHzW3I
zbF=$FtB+~v_ijcq@Nbs}b2o$Ja@U(4u}B9w7`X;tD+xH@>~8VzK#$THS*H8$b21eO
z6(3?o7q*k9?piWy)(eH*Rl;~J=vz<lJW1Lxa9qP^Z5&teBYclA^*f^xa>T{2itSFr
z5K|?`LZ#W=fCsJtXW22*X9vl}hY>gjmaSz?(HLj@Vr205bGK0=mQE@nQc5vh{)A-;
z$r-!wi$s(bO<~-5ktMfbsUi#nH*_D6ZaeKVA?az%0(j&dO;k8`wv=r&g-keh!E-+N
zerQ+wXpxQj>FEoh=xniqbHw=|96F*&<iy8T;;^C1kuxb}{XzMgeF`!y+CA)2&St;_
zsM1*qSkY3+69ISDrc?CrbzTGsC?v;JMzfZR0oBJ`KSDZ?zt(2SujE!S8OU<gaHNZ3
z@iGHBl{9yC7HcuN&~7+9<oKbhT2SZsE_nX93i!Al6Egq&6QK8xr?fK_g_P<0YB^y~
z3Jw32xgnG_zZ`dx4@!xd{q!5qP41vwaD7A&hDNx!$716(^CHT@Su{GsX9cVWrdbmi
z#EC0v2qZ<VkV?v$cmKlIo0KQ6;qv3kO3rM`8ms1Hn?zezb~ufn9Ldzr?~4a+2V0A!
zz0(i@AqR#!e%pk`@e?kx-P*WBEj)9Duh!*|;n}f_n6i>g{Wy{-HIjq`+co3tr)sAp
z4zoI^J%cZ?3=Dbd{xiXPrqm~5@6}5F7&kQxHy`^EJ6MEC_@)$IP?Pm@H`M7XV>SV|
zGW|;`m$a$XO;Tl<OY7GpRCK+b8Q>5W_d^n|!1XFN?+)=Jdg0|tHV*=Gw|#_$86aUy
zI|D27+X;^{k^%8ojJ*aGvyW){r^9#MtmZY18eUbZe(+$w2)bYh`qlKo1HR$b!wX+?
zRSiApd*Z%>rR?Bv!lCY3-JoIxjw`~$O41#q1Ps<zVG{&rT*sg{BAMyzSnLcn<jeF0
zw@Q2BBs*ff4T84KhV^X1_em0l>p7xKml8?^z%ywkpo}gQ2%#+#XMj-TOHVNpVk8+`
zTX~Z#?s_8jyIe~or__hoQbyoJ4&5_iV)V;3B6u+jI(a5^I<XHcM(?-wFF9h`M!XOK
zH=smqiu;d4oqRcYVuWYQOxU&?TKfubFwlp>;+2Ib*S+Aw<bTZ@(llFn<1xsvi~Ku{
z8I?2oCm4CHMpBKHfL=bvg30`m;#Q=3N)dYfFat)vbnZB)BbdiejvWO~>qymP`R2%4
zGSCCcB}oe=JRUATNXmCYQ92tVn9~3j>n`+5G2C~E{9$LPcWgQaB!MNtx-$6%b8Yi_
z(~Po}QjZD)>Ao<ULE}xe{mQydY8fLg+u*6nEA~~yhW?jM^$!q|P!e)u5&VbbgmDfu
zMRV%dp0{*;NyOQ1Y(sZu!WOslkww%h)RVJ?$MS8-4~Ja7uqMp@w1z=ujM2&nnqBJS
z@HNXbT-A@?;OF5&oPq5;2!#Ekj3A$kj)qCafpjHNXT0QqkxRU{gmOYuiqkxEU6L>M
zC%`S1Tb{=X5x;oy82ireZs^V-^Wxmro$1BJUGWyagNluYji-&ajrm@<i_I_KIqObn
zaU^1$(uSL$;UYyb2SIm$%dy7#40t$0v4F1;angVB)XNh%zR12{(CY5XttF&TTL?a#
zl618cynZsS6lzoobZMF_5q>AXMhg9=!_mHfS9a~t4=0KzGB2>BCd^6T8(y!FuX0;M
z&cV5ZH~7QKcxDx^s)@-}>a-drCne@eE>&KOAg2Vh3X0uTJ*vqNUGersFO8(fA$LyC
ze9sDZeN6J^R4%x;0$C>CA-_W(GJMzis*};2cg`B8r&6a??j_q2(H-QI{sI&}JMsZ~
zbT?#yQCoOf(i^d~sK7YwxLxiq6E~WsDfW;!I;)&evYl@71Qioo1{Z@diElY1ObfsS
z1AThuEzp>_@6?lz(3AS!LTjEf2RLQg-sgfRNJ0xLHV=?#;6Of?V-))M{cBiyV>=?W
zmk`MtRiCu&Rg--8LT&h~Kzgc@^K@TIDogT>n3VC&vY?z(Qpehr+*j;JpOk}anr@^%
zip#49@^H1Ua>E_fospe05rm#O7cTo-U?=zN6x!)i_?DM(CLWNctl90IP<fzrMJn9#
zf-2k9U);f>MO_9@5;!RE-a98sK&Sj#P%kvx*+%vF4LqEbf?20-W@H2Vn2cg${p)>?
z#OA8>INhD*c*O!bWdt?^Ug}HXSbWfHuh#3qeyrpbdZ#gcPYRtPR9Llt8}zU80=DUW
zf&31JZh3M>XeqR?$;Q8MuBPqhVN!EX&EPImcXcV9c>e>NO!Ws#%kd)zP3qRuR_Kny
z-3)%_b_{z&m7a+F+S>V5j0Z(wkIfhn&!V`~9PrDj6F=)j&|df+1bF=Inuf)WhKKp_
zhpS@61pYgi(b-~*JJE#XJ9MS_Add|Gw}h^Ewd(gKMn)Uj=-mZ5$MjniVU?@E6PQ=(
z0SH<DtA6bMup7jJXQncT)M9ALKDb=v!-%SB@C3nk%KXsa;$>2(3%!ugR}qio(iIA#
zORm5QIKLlorV(f!nV6#uO~~}Z!bONpgTxN7)<24Cb3FQa7>ynRJBE2MF>|3Sk+9!D
ze<<ZntLfb~ULN0&9oVmkwf1$tz?@SQ-wjp1!MU0s-b(mnHarC4#}+CoDbZw8u5P)R
zX&BpWrX}wGutd2l2)WRU3b7o%v`?k(t)v=MyFD-A&M2H74vQ3J>Cp}IGsOG>3YwL^
zox#b;xPzm(2Uw6t`TYpym6{0F2I`tds<lzyy!`A(AGY)O<X5Tgffqq7_Sn_<Rf_Uu
z1B2jJ1a?$$gQ+@f@<TwYbMLCj@#nVLZuO*&C+(?jnP%s6-}IoyYyY4%ecTNH1-wi>
z%W^K!p{(-%Q^~c!Lz%Vl+ET5GTBFEq_$cH0zH=WGV&op<(#j=7VqM3`H7(VYVT33|
z)Jjo_TDit&ql8+DVw5E8mJOv)F1zQQF|PI7@0)kdJm);;@;vW3&zbrC=Kp`s<iU3g
z(MXiW7q|Se$wL7rZm&tM46}BM+F~dxCvEg8H>^%-Gi9^%+Hh3FIvv%jtk*xe>{u=F
z3AN_r*(YN<hOcftK{KJVW#Ss|ivy&S^O*NBl3V(`$C-sUjLD`S>ealo`i0mZJ(^Z{
zCo-R9^bdC5nKYgDolb5}ZTViMP)xK!m~d%0-Z?rTwLYs=^0&@qx_&I@CIh!$hHY2F
zhZ4S&-0giOBu6^gD;epK|J13Nm`6wmylF1$(^qf1{_|NT1~t0kC9UE58}xbel_jxw
z`3>RX>fZRj47^FMC@oXWjp*w&W!EmC>hEr<EkReEaozYK)lNmfVP_8a1M-);%?UU*
z&?<e|*+e{u#ZjuW-G1A+%@(yQ)>S$GwmeHlu5wcHbM&u?;9uXz%CYMouT6FSwZ=5T
zu%nrYd2G`AePu<ATkPJKvJ#a&GIx4bI|an<ad<7dPcB$j6YFdJlbi&nZmhE|_fS!l
zW93DuA<k=mg8t6NrwQY(hNjHgeJ5&DWCQFBRfJwWJDQpv+xyIGDzHP-a?<zWI7k0|
zR7}s8n4aq*brccmOZ|Xc`=;PnwUSi%t78%+;w4*jbr?go0iCuvnH^S}Qnn=Z-k5Nu
z<y@f|Cq2}YqOKW~d>gXYtA1;VUaIBv)uK(mC4aYPr^9h&!#u&Wr)eoi8ZGc~SueD#
z$2LUQ=>=_TmL%!KG-PO_G^1LJ^rUb(tD<TRCL~t`I7NZ#&O}GD^~c;+=dA)3JHIo0
zbRxB3RL@lOcZuWa7vj;ZTj7>146;#uLNX^gr#|wH10$btTmAvM&(tVs*RUtqX!~~g
z!;(~VXO8eMxE$?Pvt*(D!Y4G$4b^w~C7)NOX@_8^qaSNduDW$k>ZW;^R=H!CO1a$^
zo6LqTQLBlH8+rm`%*4`toVXXy9g^gqFmn`O_m!E96{F<pD>muJd!=q#Dq-@I3RpYT
z9l;}{W)I@<vIWuaTmEQ$GikF)LW`!?`k%I^cQEYO*QU%Iuf7r4JT2UK>JLF-?m?$>
z)|j8_P!~OvY!phCPmI0vXrOGM)VzRowcEH1y}6rC5%%a)t&AZJhdDKlR5Cq6KH7Js
zn;wzdmE2wIWh&SV89%&I8bk9^66r*}!Z1|(P=*QN^fPUp1FAv&Z_kXo=XbSe(97<h
zE)@>`Z>`n&!CEl^qwsQj5m94W;Lj1%s)O+|=zU7r5}i9f*o;`G48{frciM}b=2Rs0
zHui@axi#I*uyl!xS=q$4$R@yM?^XkXg3QlC0?yjj3dd3+o86UdiXMs=DUL?F1_&s+
zZbRp@sCgp!?dMqS-MRjO5s;)Uwf%?o9G17M32Gxg#9)LriQgf{&Zlu!Vr)LCwQX&U
z*Yes`B*tp3LnYS0FiI?rZQ0Rc^4G<6;a<$`HtxHtYt<eU2D`UcBpkMr_*+|RihU&O
z`S~}(T~$^Ou;Eh=r^ZLO4(?^feB`7#4Rya0j?c}i=}NoRuT*$L`+at`_W9BmGGjvK
zFGIyg&5v(?AtQ4v^P!aEQO&;T@n%M7d3wfRhyU+A{W9$0&pkiJ4wadDj45WS+TP<d
zM~DLepesAB?Hckf*6qrUztnU+w9AG~+ShOQ<w@3~($D?k!$+2-ifF5R1)Aj_bp#MK
z0ER-d*esu-I`0T1K%L~HgI~ks7J%J}6~1?&b=2U)0{0nfCQnMKGo#RsS|6jzr<PGK
z(v}yi)Hifjt1GO+G}rWR3{uyn)rCe#XR<9SPf*RoYJR-DqM%V#@z`3Y>moMmIr$;l
zebwogHVdR(_F<DxJyN+ibYVQZk-%JkXxkl!VWVX_?>FIB4HOhkXzb2N_wD+-J5FsU
zhPfQSMQKcK>T1SDmUq{w^g*BJqiz#c5x&wF9{s_xr^ViDGdZqXuFR@ZS^Xzle5J<J
zP{G@FI{emtRf~*Q(isJT9#)yL1NLmK*}1gnKQ*5Q8y&%?VK&_APsP83aRMKt+rtf;
zKAp3!?u$4#a?S6i6K<gOXs@(<06Pt1QS~a~bm{)V4~1iSkK5lImY%RK*>EuB^TrXO
zC!^Qy-q)O{+qP#qa=Fga6rZoG3LgpvtAc}Yvx|aIVwWPs-vrF}+UZ&vSt(n4`nU(G
z()@f}d}l#lK48wR`xmGj3*peJcq|zWlOPCS!~i+x(PzJjl#>WN;Cl`!pHa*q<;XBx
zZ~hmgoQn|8hm5~I5xnvK1uFm6q<j-9hhhBew9QR5n2g!oz=y{FPpq7ek|PpJKt2gz
ze))-)<7>HXNQlX~i-bksdLj`C79v2^2?QiJcag~mlZAirYk9r++NBbYN8tFy7vdFi
z;cz(K;x9n-d55FTk+cz3%)&TDNZNtiX*BV4Lx=-hrok)?I|$)y&*MVF^U2T_bU(MC
z{ou$o#r-axbU!q#1{2jFwBkVq!#_|%9gyIG)`RZif6#M(pqd}u1H%hwJ#Zi!V3DA;
z45$VM5JUjcO9Yar2NDs2NLV0Ykn|x4WFb6|NSTBOk_-~KXQXacM*(>puZ*<EfwnkM
z#!`U<<s}m65IcxX#MXQp$T-|G;@3=`0Bj;U5|B6~fxbu?JCnx&zp#iuKuH9?Abw(j
z<hQ{iej`4E@n(HT+5#QlKA>hmSaAPHA0&_X&6VI@uoUi71nS&r^Yru8N1i-T#6nms
zL<9{fcoG4-mp8L!es15WofiXq8Wr*j!lSe3xd3X7p9eP#xB)T4)1T`ULU%tP2zCP`
z2ppSvp@{?nj(|RZ=ILMxi2~vad7*u0bPySaK<qEnk;n*%WT6g%U=XVdb)==^!7z>j
z_S%AWFa<|iVuOssF17(fB;sNlKt~1<y>L7jhQKCRtivwp3ls56bYvoV@f^Tp611cp
z1tKrK4`6t)T?!fb?pWXp_|}K!^<~gq0BMQNJrdE%Gt3=XBZP%R^Ya5+ckV$2Bi{pH
aht5))0vRrJ26t7#ctl80L4ms4K=2;}4>Sz`

literal 0
HcmV?d00001

diff --git a/Manuscript/fig2.pdf b/Manuscript/fig2b.pdf
similarity index 100%
rename from Manuscript/fig2.pdf
rename to Manuscript/fig2b.pdf
diff --git a/Response_Letter/Response_Letter.tex b/Response_Letter/Response_Letter.tex
index 128db85..73ec0a3 100644
--- a/Response_Letter/Response_Letter.tex
+++ b/Response_Letter/Response_Letter.tex
@@ -35,14 +35,19 @@ A few things need to be modified -- minor -- for the paper to be accepted.
 }
 \\
 \alert{
+We thank the reviewer for these positive comments and supporting publication of the present manuscript.
 }
 
 \begin{enumerate}
 
 \item 
-{Table I should be moved to the SI. Key entries could be kept if explicitly discussed in the manuscript.}
+{Table I should be moved to the SI. 
+Key entries could be kept if explicitly discussed in the manuscript.}
 \\
 \alert{
+We would prefer to keep Table I in the main text as it is the core table of the present manuscript.
+We believe that it would be unfortunate to move (parts of) the table in the SI as it would not be ideal for the reader to have to switch between the main article and the SI depending on the system he/she is interested in.
+We hope that the reviewer will understand our motivations.
 }
 
 \item 
@@ -50,18 +55,21 @@ A few things need to be modified -- minor -- for the paper to be accepted.
 The very important classes of transition metal complexes are not covered.}
 \\
 \alert{
+We have modified the conclusion to stress these two points.
 }
 
 \item 
 {Make it VERY clear that the excitation energies are Single-State CASPTX excitation energies.}
 \\
 \alert{
+The reviewer is right. We have added this additional information in the "Computational Details" section of the manuscript.
 }
 
 \item 
 {The authors should point out the difference between dynamic sigma polarization and so-called left-right polarization (https://doi.org/10.1021/jp9528020).}
 \\
 \alert{
+This has been clarified in the discussion section of our revised manuscript.
 }
 
 \end{enumerate}
@@ -75,6 +83,8 @@ I have however some comments that I would like the authors consider.
 }
 \\
 \alert{
+Thank you for recommending publication of the present manuscript.
+Our response to Reviewer \#2's comments are given below.
 }
 
 \begin{enumerate}
@@ -84,6 +94,7 @@ I have however some comments that I would like the authors consider.
 CASSCF has largest errors, sure, but it would be interesting to see how broad is the distribution in comparison to caspt3, represented in the same way.}
 \\
 \alert{
+The histogram for CASSCF has been added (and discussed) in the revised version of the manuscript.
 }
 
 \item 
@@ -103,12 +114,18 @@ I would say the latter, so while the TBEs are useful, a comparison with the expe
 }
 
 \item 
-{Can the authors also discuss if there is any basis set effect? Are the TBEs obtained with the same basis set as their calculations? 
+{Can the authors also discuss if there is any basis set effect? 
+Are the TBEs obtained with the same basis set as their calculations? 
 Probably not. 
 The authors could be more a bit more specific about the method/basis set of each of the TBEs in the SI, so that everyone can compare for a specific molecule (if desired) the setups in which the values where obtained. 
 Again, this should be discussed against the experimental value in order to convey a clear recommendation to the readers.}
 \\
 \alert{
+All calculations reported in the present manuscript have been performed with the very same basis as the TBEs (aug-cc-pVTZ), as mentioned in the "Computational Details" section. 
+We have modified the corresponding sentence to stress this further.
+Having an augmented triple-$\zeta$ basis allows us to describe faithfully all excited states (including the Rydberg states).
+Therefore, there is no basis set effect between the CASPT2/CASPT3 and TBE excitation energies.
+Moreover, because the basis set effects are extremely transferable within wave function methods, we can safely assume that the present conclusions would be identical in a small or larger basis set.
 }
 
 \end{enumerate}