Initialization

Needs["MaTeX`"]
SetOptions[MaTeX,"Preamble"->{"\\usepackage{amssymb,amsmath,latexsym,amsfonts,amsthm,mathpazo,xcolor,bm,mhchem,inputenc,fontenc,txfonts}"}];

HaToeV=27.21138602;

SetDirectory["~/Dropbox/Manuscripts/Ec/Data"];

Hilbert spaces

Five-membered ring

(* cyclopentadiene, furan, imidazole, pyrrole, and thiophene *)

{((Binomial[100-5,18-5])^2),((Binomial[90-5,18-5])^2),((Binomial[90-5,18-5])^2),((Binomial[95-5,18-5])^2),((Binomial[94-14,22-14])^2)}//N

(* benzene, pyrazine, pyridazine, pyridine, pyrimidine, tetrazine, and triazine *)

{((Binomial[114-6,21-6])^2),((Binomial[104-6,21-6])^2),((Binomial[104-6,21-6])^2),((Binomial[109-6,21-6])^2),((Binomial[104-6,21-6])^2),((Binomial[94-6,21-6])^2),((Binomial[99-6,21-6])^2)}//N

TableForm[{
{MoleculePlot3D[Molecule["cyclopentadiene"]],MoleculePlot3D[Molecule["furan"]],MoleculePlot3D[Molecule["imidazole"]],MoleculePlot3D[Molecule["pyrrole"]],MoleculePlot3D[Molecule["thiophene"]]}
,{"cyclopentadiene","furan","imidazole","pyrrole","thiophene"}
},TableAlignments->Center]
TableForm[{
{MoleculePlot3D[Molecule["benzene"]],MoleculePlot3D[Molecule["pyrazine"]],MoleculePlot3D[Molecule["pyridazine"]],MoleculePlot3D[Molecule["pyridine"]],MoleculePlot3D[Molecule["pyrimidine"]],MoleculePlot3D[Molecule["1,2,4,5-tetrazine"]],MoleculePlot3D[Molecule["1,3,5-triazine"]]}
,{"benzene","pyrazine","pyridazine","pyridine","pyrimidine","tetrazine","triazine"}
},TableAlignments->Center]

Cyclopentadiene

MP & CC calculations

CyclopentadieneHF=-192.808319784581;
CyclopentadieneMP={{5,-663.374798421},{6,-701.038770308868},{7,-734.485164342225},{8,-733.4333154}};
CyclopentadieneCCn={{5,-669.867722683},{7,-736.574222273275},{9,-738.366684546823}};
CyclopentadieneCC={{6,-707.237579224683},{8,-736.246708418799},{10,-738.202355812839}};

SizeTitle=30;SizeLabel=20;SizeLegend=20;
PlotOptions={
Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1/2,GridLines->Automatic,PlotMarkers->{"OpenMarkers",Large},
PlotStyle->{
{Thickness[0.01],Black},{Thickness[0.01],Darker[Green]},{Thickness[0.01],Red},{Thickness[0.01],Blue}
}
};
ListPlot[{CyclopentadieneMP,CyclopentadieneCCn,CyclopentadieneCC},Joined->True,PlotOptions,
PlotLabel->MaTeX["\\text{Cyclopentadiene/cc-pVDZ}",FontSize->SizeLabel],
FrameTicks->{{Automatic,None},{{5,6,7,8,9,10},None}},
FrameLabel->{
MaTeX["\\text{Computational cost $\\mathcal{O}(N^\\alpha)$}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{MP$n$}",FontSize->SizeLegend],MaTeX["\\text{CC$n$}",FontSize->SizeLegend],
MaTeX["\\text{CC}",FontSize->SizeLegend]
}],{Right,Top}]]

Shanks

s=CyclopentadieneMP[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

s=CyclopentadieneCCn[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

s=CyclopentadieneCC[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

CIPSI  calculations

Cyclopentadiene=Import["cyclopentadiene.dat"];

SizeTitle=30;SizeLabel=24;SizeLegend=24;
PlotOptions={
Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1,GridLines->Automatic,PlotMarkers->"OpenMarkers",
PlotStyle->{
{Thick,Black},{Thick,Darker[Green]},{Thick,Red},{Thick,Blue}
}
};
start=9;
CCSDT={{0.01,CyclopentadieneCC[[2,2]]},{(10^9),CyclopentadieneCC[[2,2]]}};
CCSDTQ={{0.01,CyclopentadieneCC[[3,2]]},{(10^9),CyclopentadieneCC[[3,2]]}};
SCI=Table[{Cyclopentadiene[[k,1]],(10^3)(Cyclopentadiene[[k,2]]-CyclopentadieneHF)},{k,start,Length[Cyclopentadiene]}];
SCIPT2=Table[{Cyclopentadiene[[k,1]],(10^3)(Cyclopentadiene[[k,3]]-CyclopentadieneHF)},{k,start,Length[Cyclopentadiene]}];
SCIrPT2=Table[{Cyclopentadiene[[k,1]],(10^3)(Cyclopentadiene[[k,5]]-CyclopentadieneHF)},{k,start,Length[Cyclopentadiene]}];
ListLogLinearPlot[{
CCSDTQ,SCI,SCIPT2,SCIrPT2
},
FrameLabel->{
MaTeX["\\text{Number of determinants}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{CCSDTQ}",FontSize->SizeLegend],MaTeX["\\text{$\\Delta E_\\text{var.}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.} + E_\\text{PT2}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.} + E_\\text{rPT2}$}",FontSize->SizeLegend]
}],{Right,Top}],
PlotRange->{{(10^3),3 (10^8)},{-780,-550}},
PlotOptions,
Joined->False,
PlotRange->All
]
(*Export["furanEvsNdet.pdf",%];*)

SizeTitle=30;SizeLabel=24;SizeLegend=24;
PlotOptions={Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1,GridLines->Automatic,PlotMarkers->"OpenMarkers",
PlotStyle->{
{Thick,Gray},{Thick,Black},{Thick,Blue},{Thick,Red},{Thick,Dashed,Green}
}
};
start=9;
nfit=4;
CCSDT={{-1000,CyclopentadieneCC[[2,2]]},{0,CyclopentadieneCC[[2,2]]}};
CCSDTQ={{-1000,CyclopentadieneCC[[3,2]]},{0,CyclopentadieneCC[[3,2]]}};
SCIPT2=Table[{(10^3)(Cyclopentadiene[[k,3]]-Cyclopentadiene[[k,2]]),(10^3)(Cyclopentadiene[[k,2]]-CyclopentadieneHF)},{k,start,Length[Cyclopentadiene]}];
fitPT2=Fit[Take[SCIPT2,-nfit],{1,x},x]
SCIrPT2=Table[{(10^3)(Cyclopentadiene[[k,5]]-Cyclopentadiene[[k,2]]),(10^3)(Cyclopentadiene[[k,2]]-CyclopentadieneHF)},{k,start,Length[Cyclopentadiene]}];
fitrPT2=Fit[Take[SCIrPT2,-nfit],{1,x},x]
Show[{
ListPlot[{
CCSDT,CCSDTQ,SCIPT2,SCIrPT2
},
PlotLabel->MaTeX["\\text{Cyclopentadiene/cc-pVDZ}",FontSize->SizeLabel],
FrameLabel->{
MaTeX["\\text{$E_\\text{PT2}$ or $E_\\text{rPT2}$ ($mE_h$)}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{CCSDT}",FontSize->SizeLegend],
MaTeX["\\text{CCSDTQ}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.}$ vs $E_\\text{PT2}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.}$ vs $E_\\text{rPT2}$}",FontSize->SizeLegend]
}],{Right,Top}],
PlotOptions,
PlotRange->{{-200,0},{-500,-780}},
Joined->False
]
,
Plot[{fitPT2,fitrPT2},{x,-200,0},PlotStyle->{{Thick,Dashed,Blue},{Thick,Dashed,Red}}]
}]
(*Export["furanEcvsPT2.pdf",%];*)

Furan

Shanks for the 6-31G basis

FuranHF=-228.523990058384;
FuranCCn={{5,-229.001730644282},{7,-229.040591463170},{9,-229.041751058791},{11,-229.041751139916}};
FuranCC={{6,-229.022381696587},{8,-229.039971317883},{10,-229.041597770327}};

s=(FuranCCn[[;;,2]]-FuranHF)(10^3)
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

s=(FuranCC[[;;,2]]-FuranHF)(10^3)
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

CIPSI  calculations

Furan=Import["furan_6-31g.dat"];

SizeTitle=30;SizeLabel=24;SizeLegend=24;
PlotOptions={Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1,GridLines->Automatic,PlotMarkers->"OpenMarkers",
PlotStyle->{
{Thick,Gray},{Thick,Black},{Thick,Blue},{Thick,Red},{Thick,Dashed,Green}
}
};
start=9;
nfit=3;
CCSDT={{-1000,(FuranCC[[2,2]]-FuranHF)(10^3)},{0,(FuranCC[[2,2]]-FuranHF)(10^3)}};
CCSDTQ={{-1000,(FuranCC[[3,2]]-FuranHF)(10^3)},{0,(FuranCC[[3,2]]-FuranHF)(10^3)}};
SCIPT2=Table[{(10^3)(Furan[[k,3]]-Furan[[k,2]]),(10^3)(Furan[[k,2]]-FuranHF)},{k,start,Length[Furan]}];
fitPT2=Fit[Take[SCIPT2,-nfit],{1,x},x]
SCIrPT2=Table[{(10^3)(Furan[[k,5]]-Furan[[k,2]]),(10^3)(Furan[[k,2]]-FuranHF)},{k,start,Length[Furan]}];
fitrPT2=Fit[Take[SCIrPT2,-nfit],{1,x},x]
Show[{
ListPlot[{
CCSDT,CCSDTQ,SCIPT2,SCIrPT2
},
PlotLabel->MaTeX["\\text{Furan/cc-pVDZ}",FontSize->SizeLabel],
FrameLabel->{
MaTeX["\\text{$E_\\text{PT2}$ or $E_\\text{rPT2}$ ($mE_h$)}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{CCSDT}",FontSize->SizeLegend],
MaTeX["\\text{CCSDTQ}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.}$ vs $E_\\text{PT2}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.}$ vs $E_\\text{rPT2}$}",FontSize->SizeLegend]
}],{Right,Top}],
PlotOptions,
PlotRange->{{-200,0},Automatic},
Joined->False
]
,
Plot[{fitPT2,fitrPT2},{x,-200,0},PlotStyle->{{Thick,Dashed,Blue},{Thick,Dashed,Red}}]
}]
(*Export["furanEcvsPT2.pdf",%];*)

MP & CC calculations

FuranHF=-228.643312572484575;
FuranMP={{5,-707.450373824},{6,-727.820193088592},{7,-766.627890777360},{8,-759.8932275}};
FuranCCn={{5,-717.183336927},{7,-765.730048755297},{9,-766.866256888079}};
FuranCC={{6,-734.972685409664},{8,-764.332120385917},{10,-766.65420877679}};

SizeTitle=30;SizeLabel=20;SizeLegend=20;
PlotOptions={
Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1/2,GridLines->Automatic,PlotMarkers->{"OpenMarkers",Large},
PlotStyle->{
{Thickness[0.01],Black},{Thickness[0.01],Darker[Green]},{Thickness[0.01],Red},{Thickness[0.01],Blue}
}
};
ListPlot[{FuranMP,FuranCCn,FuranCC},Joined->True,PlotOptions,
FrameTicks->{{Automatic,None},{{5,6,7,8,9,10},None}},
PlotLabel->MaTeX["\\text{Furan/cc-pVDZ}",FontSize->SizeLabel],
FrameLabel->{
MaTeX["\\text{Computational cost $\\mathcal{O}(N^\\alpha)$}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{MP$n$}",FontSize->SizeLegend],MaTeX["\\text{CC$n$}",FontSize->SizeLegend],
MaTeX["\\text{CC}",FontSize->SizeLegend]
}],{Right,Top}]]

Shanks

s=FuranMP[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

s=FuranCCn[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

s=FuranCC[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

CIPSI  calculations

Furan=Import["furan.dat"];

SizeTitle=30;SizeLabel=24;SizeLegend=24;
PlotOptions={
Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1,GridLines->Automatic,PlotMarkers->"OpenMarkers",
PlotStyle->{
{Thick,Black},{Thick,Darker[Green]},{Thick,Red},{Thick,Blue}
}
};
start=9;
CCSDT={{0.01,FuranCC[[2,2]]},{(10^9),FuranCC[[2,2]]}};
CCSDTQ={{0.01,FuranCC[[3,2]]},{(10^9),FuranCC[[3,2]]}};
SCI=Table[{Furan[[k,1]],(10^3)(Furan[[k,2]]-FuranHF)},{k,start,Length[Furan]}];
SCIPT2=Table[{Furan[[k,1]],(10^3)(Furan[[k,3]]-FuranHF)},{k,start,Length[Furan]}];
SCIrPT2=Table[{Furan[[k,1]],(10^3)(Furan[[k,5]]-FuranHF)},{k,start,Length[Furan]}];
ListLogLinearPlot[{
CCSDTQ,SCI,SCIPT2,SCIrPT2
},
FrameLabel->{
MaTeX["\\text{Number of determinants}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{CCSDTQ}",FontSize->SizeLegend],MaTeX["\\text{$\\Delta E_\\text{var.}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.} + E_\\text{PT2}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.} + E_\\text{rPT2}$}",FontSize->SizeLegend]
}],{Right,Top}],
PlotRange->{{(10^3),3 (10^8)},{-780,-550}},
PlotOptions,
Joined->False,
PlotRange->All
]
(*Export["furanEvsNdet.pdf",%];*)

SizeTitle=30;SizeLabel=24;SizeLegend=24;
PlotOptions={Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1,GridLines->Automatic,PlotMarkers->"OpenMarkers",
PlotStyle->{
{Thick,Gray},{Thick,Black},{Thick,Blue},{Thick,Red},{Thick,Dashed,Green}
}
};
start=9;
nfit=3;
CCSDT={{-1000,FuranCC[[2,2]]},{0,FuranCC[[2,2]]}};
CCSDTQ={{-1000,FuranCC[[3,2]]},{0,FuranCC[[3,2]]}};
SCIPT2=Table[{(10^3)(Furan[[k,3]]-Furan[[k,2]]),(10^3)(Furan[[k,2]]-FuranHF)},{k,start,Length[Furan]}];
fitPT2=Fit[Take[SCIPT2,-nfit],{1,x},x]
SCIrPT2=Table[{(10^3)(Furan[[k,5]]-Furan[[k,2]]),(10^3)(Furan[[k,2]]-FuranHF)},{k,start,Length[Furan]}];
fitrPT2=Fit[Take[SCIrPT2,-nfit],{1,x},x]
Show[{
ListPlot[{
CCSDT,CCSDTQ,SCIPT2,SCIrPT2
},
PlotLabel->MaTeX["\\text{Furan/cc-pVDZ}",FontSize->SizeLabel],
FrameLabel->{
MaTeX["\\text{$E_\\text{PT2}$ or $E_\\text{rPT2}$ ($mE_h$)}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{CCSDT}",FontSize->SizeLegend],
MaTeX["\\text{CCSDTQ}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.}$ vs $E_\\text{PT2}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.}$ vs $E_\\text{rPT2}$}",FontSize->SizeLegend]
}],{Right,Top}],
PlotOptions,
PlotRange->{{-200,0},{-500,-780}},
Joined->False
]
,
Plot[{fitPT2,fitrPT2},{x,-200,0},PlotStyle->{{Thick,Dashed,Blue},{Thick,Dashed,Red}}]
}]
(*Export["furanEcvsPT2.pdf",%];*)

Imidazole

MP & CC calculations

ImidazoleHF=-224.835356714467;
ImidazoleMP={{5,-720.429915262},{6,-737.828542914060},{7,-777.223689935675},{8,-770.7756855}};
ImidazoleCCn={{5,-729.015795304},{7,-776.111512572204},{9,-777.204281879068}};
ImidazoleCC={{6,-744.198886306645},{8,-774.551173949665},{10,-776.949476371663}};

SizeTitle=30;SizeLabel=20;SizeLegend=20;
PlotOptions={
Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1/2,GridLines->Automatic,PlotMarkers->{"OpenMarkers",Large},
PlotStyle->{
{Thickness[0.01],Black},{Thickness[0.01],Darker[Green]},{Thickness[0.01],Red},{Thickness[0.01],Blue}
}
};
ListPlot[{ImidazoleMP,ImidazoleCCn,ImidazoleCC},Joined->True,PlotOptions,
FrameTicks->{{Automatic,None},{{5,6,7,8,9,10},None}},
PlotLabel->MaTeX["\\text{Imidazole/cc-pVDZ}",FontSize->SizeLabel],
FrameLabel->{
MaTeX["\\text{Computational cost $\\mathcal{O}(N^\\alpha)$}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{MP$n$}",FontSize->SizeLegend],MaTeX["\\text{CC$n$}",FontSize->SizeLegend],
MaTeX["\\text{CC}",FontSize->SizeLegend]
}],{Right,Top}]]

Shanks

s=ImidazoleMP[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

s=ImidazoleCCn[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

s=ImidazoleCC[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

CIPSI  calculations

Imidazole=Import["imidazole.dat"];

SizeTitle=30;SizeLabel=24;SizeLegend=24;
PlotOptions={
Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1,GridLines->Automatic,PlotMarkers->"OpenMarkers",
PlotStyle->{
{Thick,Black},{Thick,Darker[Green]},{Thick,Red},{Thick,Blue}
}
};
start=9;
CCSDT={{0.01,ImidazoleCC[[2,2]]},{(10^9),ImidazoleCC[[2,2]]}};
CCSDTQ={{0.01,ImidazoleCC[[3,2]]},{(10^9),ImidazoleCC[[3,2]]}};
SCI=Table[{Imidazole[[k,1]],(10^3)(Imidazole[[k,2]]-ImidazoleHF)},{k,start,Length[Imidazole]}];
SCIPT2=Table[{Imidazole[[k,1]],(10^3)(Imidazole[[k,3]]-ImidazoleHF)},{k,start,Length[Imidazole]}];
SCIrPT2=Table[{Imidazole[[k,1]],(10^3)(Imidazole[[k,5]]-ImidazoleHF)},{k,start,Length[Imidazole]}];
ListLogLinearPlot[{
CCSDTQ,SCI,SCIPT2,SCIrPT2
},
FrameLabel->{
MaTeX["\\text{Number of determinants}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{CCSDTQ}",FontSize->SizeLegend],MaTeX["\\text{$\\Delta E_\\text{var.}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.} + E_\\text{PT2}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.} + E_\\text{rPT2}$}",FontSize->SizeLegend]
}],{Right,Top}],
PlotRange->{{(10^3),3 (10^8)},{-780,-550}},
PlotOptions,
Joined->False,
PlotRange->All
]
(*Export["furanEvsNdet.pdf",%];*)

SizeTitle=30;SizeLabel=24;SizeLegend=24;
PlotOptions={Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1,GridLines->Automatic,PlotMarkers->"OpenMarkers",
PlotStyle->{
{Thick,Gray},{Thick,Black},{Thick,Blue},{Thick,Red},{Thick,Dashed,Green}
}
};
start=9;
nfit=3;
CCSDT={{-1000,ImidazoleCC[[2,2]]},{0,ImidazoleCC[[2,2]]}};
CCSDTQ={{-1000,ImidazoleCC[[3,2]]},{0,ImidazoleCC[[3,2]]}};
SCIPT2=Table[{(10^3)(Imidazole[[k,3]]-Imidazole[[k,2]]),(10^3)(Imidazole[[k,2]]-ImidazoleHF)},{k,start,Length[Imidazole]}];
fitPT2=Fit[Take[SCIPT2,-nfit],{1,x},x]
SCIrPT2=Table[{(10^3)(Imidazole[[k,5]]-Imidazole[[k,2]]),(10^3)(Imidazole[[k,2]]-ImidazoleHF)},{k,start,Length[Imidazole]}];
fitrPT2=Fit[Take[SCIrPT2,-nfit],{1,x},x]
Show[{
ListPlot[{
CCSDT,CCSDTQ,SCIPT2,SCIrPT2
},
PlotLabel->MaTeX["\\text{Imidazole/cc-pVDZ}",FontSize->SizeLabel],
FrameLabel->{
MaTeX["\\text{$E_\\text{PT2}$ or $E_\\text{rPT2}$ ($mE_h$)}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{CCSDT}",FontSize->SizeLegend],
MaTeX["\\text{CCSDTQ}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.}$ vs $E_\\text{PT2}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.}$ vs $E_\\text{rPT2}$}",FontSize->SizeLegend]
}],{Right,Top}],
PlotOptions,
PlotRange->{{-200,0},{-500,-780}},
Joined->False
]
,
Plot[{fitPT2,fitrPT2},{x,-200,0},PlotStyle->{{Thick,Dashed,Blue},{Thick,Dashed,Red}}]
}]
(*Export["furanEcvsPT2.pdf",%];*)

Pyrrole

MP & CC calculations

PyrroleHF=-208.828614693633;
PyrroleMP={{5,-695.721950709},{6,-720.603108743368},{7,-756.459033548742},{8,-752.2589064}};
PyrroleCCn={{5,-702.531904930},{7,-756.270764152895},{9,-757.576382167585}};
PyrroleCC={{6,-725.712852643528},{8,-755.183721462796},{10,-757.360899120397}};

SizeTitle=30;SizeLabel=20;SizeLegend=20;
PlotOptions={
Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1/2,GridLines->Automatic,PlotMarkers->{"OpenMarkers",Large},
PlotStyle->{
{Thickness[0.01],Black},{Thickness[0.01],Darker[Green]},{Thickness[0.01],Red},{Thickness[0.01],Blue}
}
};
ListPlot[{PyrroleMP,PyrroleCCn,PyrroleCC},Joined->True,PlotOptions,
FrameTicks->{{Automatic,None},{{5,6,7,8,9,10},None}},
PlotLabel->MaTeX["\\text{Pyrrole/cc-pVDZ}",FontSize->SizeLabel],
FrameLabel->{
MaTeX["\\text{Computational cost $\\mathcal{O}(N^\\alpha)$}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{MP$n$}",FontSize->SizeLegend],MaTeX["\\text{CC$n$}",FontSize->SizeLegend],
MaTeX["\\text{CC}",FontSize->SizeLegend]
}],{Right,Top}]]

Shanks

s=PyrroleMP[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

s=PyrroleCCn[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

s=PyrroleCC[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

CIPSI  calculations

Pyrrole=Import["pyrrole.dat"];

SizeTitle=30;SizeLabel=24;SizeLegend=24;
PlotOptions={
Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1,GridLines->Automatic,PlotMarkers->"OpenMarkers",
PlotStyle->{
{Thick,Black},{Thick,Darker[Green]},{Thick,Red},{Thick,Blue}
}
};
start=9;
CCSDT={{0.01,PyrroleCC[[2,2]]},{(10^9),PyrroleCC[[2,2]]}};
CCSDTQ={{0.01,PyrroleCC[[3,2]]},{(10^9),PyrroleCC[[3,2]]}};
SCI=Table[{Pyrrole[[k,1]],(10^3)(Pyrrole[[k,2]]-PyrroleHF)},{k,start,Length[Pyrrole]}];
SCIPT2=Table[{Pyrrole[[k,1]],(10^3)(Pyrrole[[k,3]]-PyrroleHF)},{k,start,Length[Pyrrole]}];
SCIrPT2=Table[{Pyrrole[[k,1]],(10^3)(Pyrrole[[k,5]]-PyrroleHF)},{k,start,Length[Pyrrole]}];
ListLogLinearPlot[{
CCSDTQ,SCI,SCIPT2,SCIrPT2
},
FrameLabel->{
MaTeX["\\text{Number of determinants}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{CCSDTQ}",FontSize->SizeLegend],MaTeX["\\text{$\\Delta E_\\text{var.}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.} + E_\\text{PT2}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.} + E_\\text{rPT2}$}",FontSize->SizeLegend]
}],{Right,Top}],
PlotRange->{{(10^3),3 (10^8)},{-780,-550}},
PlotOptions,
Joined->False,
PlotRange->All
]
(*Export["furanEvsNdet.pdf",%];*)

SizeTitle=30;SizeLabel=24;SizeLegend=24;
PlotOptions={Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1,GridLines->Automatic,PlotMarkers->"OpenMarkers",
PlotStyle->{
{Thick,Gray},{Thick,Black},{Thick,Blue},{Thick,Red},{Thick,Dashed,Green}
}
};
start=9;
nfit=3;
CCSDT={{-1000,PyrroleCC[[2,2]]},{0,PyrroleCC[[2,2]]}};
CCSDTQ={{-1000,PyrroleCC[[3,2]]},{0,PyrroleCC[[3,2]]}};
SCIPT2=Table[{(10^3)(Pyrrole[[k,3]]-Pyrrole[[k,2]]),(10^3)(Pyrrole[[k,2]]-PyrroleHF)},{k,start,Length[Pyrrole]}];
fitPT2=Fit[Take[SCIPT2,-nfit],{1,x},x]
SCIrPT2=Table[{(10^3)(Pyrrole[[k,5]]-Pyrrole[[k,2]]),(10^3)(Pyrrole[[k,2]]-PyrroleHF)},{k,start,Length[Pyrrole]}];
fitrPT2=Fit[Take[SCIrPT2,-nfit],{1,x},x]
Show[{
ListPlot[{
CCSDT,CCSDTQ,SCIPT2,SCIrPT2
},
PlotLabel->MaTeX["\\text{Pyrrole/cc-pVDZ}",FontSize->SizeLabel],
FrameLabel->{
MaTeX["\\text{$E_\\text{PT2}$ or $E_\\text{rPT2}$ ($mE_h$)}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{CCSDT}",FontSize->SizeLegend],
MaTeX["\\text{CCSDTQ}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.}$ vs $E_\\text{PT2}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.}$ vs $E_\\text{rPT2}$}",FontSize->SizeLegend]
}],{Right,Top}],
PlotOptions,
PlotRange->{{-200,0},{-500,-780}},
Joined->False
]
,
Plot[{fitPT2,fitrPT2},{x,-200,0},PlotStyle->{{Thick,Dashed,Blue},{Thick,Dashed,Red}}]
}]
(*Export["furanEcvsPT2.pdf",%];*)

Thiophene

MP & CC calculations

ThiopheneHF=-551.321002444800;
ThiopheneMP={{5,-661.487676839},{6,-689.433481554693},{7,-726.576354093504},{8,-721.5988552}};
ThiopheneCCn={{5,-669.535846587},{7,-726.268305928940},{9,-727.662544525096}};
ThiopheneCC={{6,-694.455334004781},{8,-725.135024376793},{10,-727.450489687365}};

SizeTitle=30;SizeLabel=20;SizeLegend=20;
PlotOptions={
Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1/2,GridLines->Automatic,PlotMarkers->{"OpenMarkers",Large},
PlotStyle->{
{Thickness[0.01],Black},{Thickness[0.01],Darker[Green]},{Thickness[0.01],Red},{Thickness[0.01],Blue}
}
};
ListPlot[{ThiopheneMP,ThiopheneCCn,ThiopheneCC},Joined->True,PlotOptions,
FrameTicks->{{Automatic,None},{{5,6,7,8,9,10},None}},
PlotLabel->MaTeX["\\text{Thiophene/cc-pVDZ}",FontSize->SizeLabel],
FrameLabel->{
MaTeX["\\text{Computational cost $\\mathcal{O}(N^\\alpha)$}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{MP$n$}",FontSize->SizeLegend],MaTeX["\\text{CC$n$}",FontSize->SizeLegend],
MaTeX["\\text{CC}",FontSize->SizeLegend]
}],{Right,Top}]]

Shanks

s=ThiopheneMP[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

s=ThiopheneCCn[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

s=ThiopheneCC[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

CIPSI  calculations

Thiophene=Import["thiophene.dat"];

SizeTitle=30;SizeLabel=24;SizeLegend=24;
PlotOptions={
Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1,GridLines->Automatic,PlotMarkers->"OpenMarkers",
PlotStyle->{
{Thick,Black},{Thick,Darker[Green]},{Thick,Red},{Thick,Blue}
}
};
start=9;
CCSDT={{0.01,ThiopheneCC[[2,2]]},{(10^9),ThiopheneCC[[2,2]]}};
CCSDTQ={{0.01,ThiopheneCC[[3,2]]},{(10^9),ThiopheneCC[[3,2]]}};
SCI=Table[{Thiophene[[k,1]],(10^3)(Thiophene[[k,2]]-ThiopheneHF)},{k,start,Length[Thiophene]}];
SCIPT2=Table[{Thiophene[[k,1]],(10^3)(Thiophene[[k,3]]-ThiopheneHF)},{k,start,Length[Thiophene]}];
SCIrPT2=Table[{Thiophene[[k,1]],(10^3)(Thiophene[[k,5]]-ThiopheneHF)},{k,start,Length[Thiophene]}];
ListLogLinearPlot[{
CCSDTQ,SCI,SCIPT2,SCIrPT2
},
FrameLabel->{
MaTeX["\\text{Number of determinants}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{CCSDTQ}",FontSize->SizeLegend],MaTeX["\\text{$\\Delta E_\\text{var.}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.} + E_\\text{PT2}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.} + E_\\text{rPT2}$}",FontSize->SizeLegend]
}],{Right,Top}],
PlotRange->{{(10^3),3 (10^8)},{-780,-550}},
PlotOptions,
Joined->False,
PlotRange->All
]
(*Export["furanEvsNdet.pdf",%];*)

SizeTitle=30;SizeLabel=24;SizeLegend=24;
PlotOptions={Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1,GridLines->Automatic,PlotMarkers->"OpenMarkers",
PlotStyle->{
{Thick,Gray},{Thick,Black},{Thick,Blue},{Thick,Red},{Thick,Dashed,Green}
}
};
start=9;
nfit=4;
CCSDT={{-1000,ThiopheneCC[[2,2]]},{0,ThiopheneCC[[2,2]]}};
CCSDTQ={{-1000,ThiopheneCC[[3,2]]},{0,ThiopheneCC[[3,2]]}};
SCIPT2=Table[{(10^3)(Thiophene[[k,3]]-Thiophene[[k,2]]),(10^3)(Thiophene[[k,2]]-ThiopheneHF)},{k,start,Length[Thiophene]}];
fitPT2=Fit[Take[SCIPT2,-nfit],{1,x},x]
SCIrPT2=Table[{(10^3)(Thiophene[[k,5]]-Thiophene[[k,2]]),(10^3)(Thiophene[[k,2]]-ThiopheneHF)},{k,start,Length[Thiophene]}];
fitrPT2=Fit[Take[SCIrPT2,-nfit],{1,x},x]
Show[{
ListPlot[{
CCSDT,CCSDTQ,SCIPT2,SCIrPT2
},
PlotLabel->MaTeX["\\text{Thiophene/cc-pVDZ}",FontSize->SizeLabel],
FrameLabel->{
MaTeX["\\text{$E_\\text{PT2}$ or $E_\\text{rPT2}$ ($mE_h$)}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{CCSDT}",FontSize->SizeLegend],
MaTeX["\\text{CCSDTQ}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.}$ vs $E_\\text{PT2}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.}$ vs $E_\\text{rPT2}$}",FontSize->SizeLegend]
}],{Right,Top}],
PlotOptions,
PlotRange->{{-200,0},{-500,-780}},
Joined->False
]
,
Plot[{fitPT2,fitrPT2},{x,-200,0},PlotStyle->{{Thick,Dashed,Blue},{Thick,Dashed,Red}}]
}]
(*Export["furanEcvsPT2.pdf",%];*)

Benzene

MP & CC calculations

BenzeneHF=-230.722244985640;
BenzeneMP={{5,-782.332077228},{6,-816.387142630248},{7,-858.539911321773},{8,-853.7830144}};
BenzeneCCn={{5,-789.412298096},{7,-859.142698599615},{9,-860.586074871286}};
BenzeneCC={{6,-821.788776613854},{8,-857.926313401076},{10,-860.365191225800}};

SizeTitle=30;SizeLabel=20;SizeLegend=20;
PlotOptions={
Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1/2,GridLines->Automatic,PlotMarkers->{"OpenMarkers",Large},
PlotStyle->{
{Thickness[0.01],Black},{Thickness[0.01],Darker[Green]},{Thickness[0.01],Red},{Thickness[0.01],Blue}
}
};
ListPlot[{BenzeneMP,BenzeneCCn,BenzeneCC},Joined->True,PlotOptions,
FrameTicks->{{Automatic,None},{{5,6,7,8,9,10},None}},
FrameLabel->{
MaTeX["\\text{Computational cost $\\mathcal{O}(N^\\alpha)$}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{MP$n$}",FontSize->SizeLegend],MaTeX["\\text{CC$n$}",FontSize->SizeLegend],
MaTeX["\\text{CC}",FontSize->SizeLegend]
}],{Right,Top}]]

Shanks

s=BenzeneMP[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

s=BenzeneCCn[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

s=BenzeneCC[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

CIPSI  calculations

Benzene=Import["benzene.dat"];

SizeTitle=30;SizeLabel=24;SizeLegend=24;
PlotOptions={
Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1,GridLines->Automatic,PlotMarkers->"OpenMarkers",
PlotStyle->{
{Thick,Black},{Thick,Darker[Green]},{Thick,Red},{Thick,Blue}
}
};
start=9;
CCSDT={{0.01,BenzeneCC[[2,2]]},{(10^9),BenzeneCC[[2,2]]}};
CCSDTQ={{0.01,BenzeneCC[[3,2]]},{(10^9),BenzeneCC[[3,2]]}};
SCI=Table[{Benzene[[k,1]],(10^3)(Benzene[[k,2]]-BenzeneHF)},{k,start,Length[Benzene]}];
SCIPT2=Table[{Benzene[[k,1]],(10^3)(Benzene[[k,3]]-BenzeneHF)},{k,start,Length[Benzene]}];
SCIrPT2=Table[{Benzene[[k,1]],(10^3)(Benzene[[k,5]]-BenzeneHF)},{k,start,Length[Benzene]}];
ListLogLinearPlot[{
CCSDTQ,SCI,SCIPT2,SCIrPT2
},
FrameLabel->{
MaTeX["\\text{Number of determinants}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{CCSDTQ}",FontSize->SizeLegend],MaTeX["\\text{$\\Delta E_\\text{var.}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.} + E_\\text{PT2}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.} + E_\\text{rPT2}$}",FontSize->SizeLegend]
}],{Right,Top}],
PlotRange->{{(10^3),3 (10^8)},{-880,-550}},
PlotOptions,
Joined->False,
PlotRange->All
]
(*Export["furanEvsNdet.pdf",%];*)

SizeTitle=30;SizeLabel=24;SizeLegend=24;
PlotOptions={Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1,GridLines->Automatic,PlotMarkers->"OpenMarkers",
PlotStyle->{
{Thick,Gray},{Thick,Black},{Thick,Blue},{Thick,Red},{Thick,Dashed,Green}
}
};
start=9;
nfit=4;
CCSDT={{-1000,BenzeneCC[[2,2]]},{0,BenzeneCC[[2,2]]}};
CCSDTQ={{-1000,BenzeneCC[[3,2]]},{0,BenzeneCC[[3,2]]}};
SCIPT2=Table[{(10^3)(Benzene[[k,3]]-Benzene[[k,2]]),(10^3)(Benzene[[k,2]]-BenzeneHF)},{k,start,Length[Benzene]}];
fitPT2=Fit[Take[SCIPT2,-nfit],{1,x},x]
SCIrPT2=Table[{(10^3)(Benzene[[k,5]]-Benzene[[k,2]]),(10^3)(Benzene[[k,2]]-BenzeneHF)},{k,start,Length[Benzene]}];
fitrPT2=Fit[Take[SCIrPT2,-nfit],{1,x},x]
Show[{
ListPlot[{
CCSDT,CCSDTQ,SCIPT2,SCIrPT2
},
PlotLabel->MaTeX["\\text{Benzene/cc-pVDZ}",FontSize->SizeLabel],
FrameLabel->{
MaTeX["\\text{$E_\\text{PT2}$ or $E_\\text{rPT2}$ ($mE_h$)}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{CCSDT}",FontSize->SizeLegend],
MaTeX["\\text{CCSDTQ}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.}$ vs $E_\\text{PT2}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.}$ vs $E_\\text{rPT2}$}",FontSize->SizeLegend]
}],{Right,Top}],
PlotOptions,
PlotRange->{{-200,0},{-500,-880}},
Joined->False
]
,
Plot[{fitPT2,fitrPT2},{x,-200,0},PlotStyle->{{Thick,Dashed,Blue},{Thick,Dashed,Red}}]
}]
(*Export["furanEcvsPT2.pdf",%];*)

Pyrazine

MP & CC calculations

PyrazineHF=-262.703010944756;
PyrazineMP={{5,-834.617733280},{6,-853.655354014131},{7,-902.907089183363},{8,-893.8350552}};
PyrazineCCn={{5,-844.463447870},{7,-901.497912273981},{9,-902.636842327822}};
PyrazineCC={{6,-861.028517352127},{8,-899.350733684886},{10,-902.330775312939}};

SizeTitle=30;SizeLabel=20;SizeLegend=20;
PlotOptions={
Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1/2,GridLines->Automatic,PlotMarkers->{"OpenMarkers",Large},
PlotStyle->{
{Thickness[0.01],Black},{Thickness[0.01],Darker[Green]},{Thickness[0.01],Red},{Thickness[0.01],Blue}
}
};
ListPlot[{PyrazineMP,PyrazineCCn,PyrazineCC},Joined->True,PlotOptions,
FrameTicks->{{Automatic,None},{{5,6,7,8,9,10},None}},
PlotLabel->MaTeX["\\text{Pyrazine/cc-pVDZ}",FontSize->SizeLabel],
FrameLabel->{
MaTeX["\\text{Computational cost $\\mathcal{O}(N^\\alpha)$}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{MP$n$}",FontSize->SizeLegend],MaTeX["\\text{CC$n$}",FontSize->SizeLegend],
MaTeX["\\text{CC}",FontSize->SizeLegend]
}],{Right,Top}]]

Shanks

s=PyrazineMP[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

s=PyrazineCCn[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

s=PyrazineCC[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

CIPSI  calculations

Pyrazine=Import["pyrazine.dat"];

SizeTitle=30;SizeLabel=24;SizeLegend=24;
PlotOptions={
Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1,GridLines->Automatic,PlotMarkers->"OpenMarkers",
PlotStyle->{
{Thick,Black},{Thick,Darker[Green]},{Thick,Red},{Thick,Blue}
}
};
start=9;
CCSDT={{0.01,PyrazineCC[[2,2]]},{(10^9),PyrazineCC[[2,2]]}};
CCSDTQ={{0.01,PyrazineCC[[3,2]]},{(10^9),PyrazineCC[[3,2]]}};
SCI=Table[{Pyrazine[[k,1]],(10^3)(Pyrazine[[k,2]]-PyrazineHF)},{k,start,Length[Pyrazine]}];
SCIPT2=Table[{Pyrazine[[k,1]],(10^3)(Pyrazine[[k,3]]-PyrazineHF)},{k,start,Length[Pyrazine]}];
SCIrPT2=Table[{Pyrazine[[k,1]],(10^3)(Pyrazine[[k,5]]-PyrazineHF)},{k,start,Length[Pyrazine]}];
ListLogLinearPlot[{
CCSDTQ,SCI,SCIPT2,SCIrPT2
},
FrameLabel->{
MaTeX["\\text{Number of determinants}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{CCSDTQ}",FontSize->SizeLegend],MaTeX["\\text{$\\Delta E_\\text{var.}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.} + E_\\text{PT2}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.} + E_\\text{rPT2}$}",FontSize->SizeLegend]
}],{Right,Top}],
PlotRange->{{(10^3),3 (10^8)},{-880,-550}},
PlotOptions,
Joined->False,
PlotRange->All
]
(*Export["furanEvsNdet.pdf",%];*)

SizeTitle=30;SizeLabel=24;SizeLegend=24;
PlotOptions={Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1,GridLines->Automatic,PlotMarkers->"OpenMarkers",
PlotStyle->{
{Thick,Gray},{Thick,Black},{Thick,Blue},{Thick,Red},{Thick,Dashed,Green}
}
};
start=9;
nfit=4;
CCSDT={{-1000,PyrazineCC[[2,2]]},{0,PyrazineCC[[2,2]]}};
CCSDTQ={{-1000,PyrazineCC[[3,2]]},{0,PyrazineCC[[3,2]]}};
SCIPT2=Table[{(10^3)(Pyrazine[[k,3]]-Pyrazine[[k,2]]),(10^3)(Pyrazine[[k,2]]-PyrazineHF)},{k,start,Length[Pyrazine]}];
fitPT2=Fit[Take[SCIPT2,-nfit],{1,x},x]
SCIrPT2=Table[{(10^3)(Pyrazine[[k,5]]-Pyrazine[[k,2]]),(10^3)(Pyrazine[[k,2]]-PyrazineHF)},{k,start,Length[Pyrazine]}];
fitrPT2=Fit[Take[SCIrPT2,-nfit],{1,x},x]
Show[{
ListPlot[{
CCSDT,CCSDTQ,SCIPT2,SCIrPT2
},
PlotLabel->MaTeX["\\text{Pyrazine/cc-pVDZ}",FontSize->SizeLabel],
FrameLabel->{
MaTeX["\\text{$E_\\text{PT2}$ or $E_\\text{rPT2}$ ($mE_h$)}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{CCSDT}",FontSize->SizeLegend],
MaTeX["\\text{CCSDTQ}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.}$ vs $E_\\text{PT2}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.}$ vs $E_\\text{rPT2}$}",FontSize->SizeLegend]
}],{Right,Top}],
PlotOptions,
PlotRange->{{-200,0},{-500,-880}},
Joined->False
]
,
Plot[{fitPT2,fitrPT2},{x,-200,0},PlotStyle->{{Thick,Dashed,Blue},{Thick,Dashed,Red}}]
}]
(*Export["furanEcvsPT2.pdf",%];*)

Pyridazine

MP & CC calculations

PyridazineHF=-262.669873445804;
PyridazineMP={{5,-838.740583969},{6,-857.275882097803},{7,-907.925790862876},{8,-898.2587542}};
PyridazineCCn={{5,-848.908853611},{7,-906.214393900430},{9,-907.453759605642}};
PyridazineCC={{6,-864.874528282081},{8,-903.988565095408},{10,-907.116636262810}};

SizeTitle=30;SizeLabel=20;SizeLegend=20;
PlotOptions={
Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1/2,GridLines->Automatic,PlotMarkers->{"OpenMarkers",Large},
PlotStyle->{
{Thickness[0.01],Black},{Thickness[0.01],Darker[Green]},{Thickness[0.01],Red},{Thickness[0.01],Blue}
}
};
ListPlot[{PyridazineMP,PyridazineCCn,PyridazineCC},Joined->True,PlotOptions,
FrameTicks->{{Automatic,None},{{5,6,7,8,9,10},None}},
PlotLabel->MaTeX["\\text{Pyridazine/cc-pVDZ}",FontSize->SizeLabel],
FrameLabel->{
MaTeX["\\text{Computational cost $\\mathcal{O}(N^\\alpha)$}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{MP$n$}",FontSize->SizeLegend],MaTeX["\\text{CC$n$}",FontSize->SizeLegend],
MaTeX["\\text{CC}",FontSize->SizeLegend]
}],{Right,Top}]]

Shanks

s=PyridazineMP[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

s=PyridazineCCn[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

s=PyridazineCC[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

CIPSI  calculations

Pyridazine=Import["pyridazine.dat"];

SizeTitle=30;SizeLabel=24;SizeLegend=24;
PlotOptions={
Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1,GridLines->Automatic,PlotMarkers->"OpenMarkers",
PlotStyle->{
{Thick,Black},{Thick,Darker[Green]},{Thick,Red},{Thick,Blue}
}
};
start=9;
CCSDT={{0.01,PyridazineCC[[2,2]]},{(10^9),PyridazineCC[[2,2]]}};
CCSDTQ={{0.01,PyridazineCC[[3,2]]},{(10^9),PyridazineCC[[3,2]]}};
SCI=Table[{Pyridazine[[k,1]],(10^3)(Pyridazine[[k,2]]-PyridazineHF)},{k,start,Length[Pyridazine]}];
SCIPT2=Table[{Pyridazine[[k,1]],(10^3)(Pyridazine[[k,3]]-PyridazineHF)},{k,start,Length[Pyridazine]}];
SCIrPT2=Table[{Pyridazine[[k,1]],(10^3)(Pyridazine[[k,5]]-PyridazineHF)},{k,start,Length[Pyridazine]}];
ListLogLinearPlot[{
CCSDTQ,SCI,SCIPT2,SCIrPT2
},
FrameLabel->{
MaTeX["\\text{Number of determinants}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{CCSDTQ}",FontSize->SizeLegend],MaTeX["\\text{$\\Delta E_\\text{var.}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.} + E_\\text{PT2}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.} + E_\\text{rPT2}$}",FontSize->SizeLegend]
}],{Right,Top}],
PlotRange->{{(10^3),3 (10^8)},{-880,-550}},
PlotOptions,
Joined->False,
PlotRange->All
]
(*Export["furanEvsNdet.pdf",%];*)

SizeTitle=30;SizeLabel=24;SizeLegend=24;
PlotOptions={Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1,GridLines->Automatic,PlotMarkers->"OpenMarkers",
PlotStyle->{
{Thick,Gray},{Thick,Black},{Thick,Blue},{Thick,Red},{Thick,Dashed,Green}
}
};
start=9;
nfit=4;
CCSDT={{-1000,PyridazineCC[[2,2]]},{0,PyridazineCC[[2,2]]}};
CCSDTQ={{-1000,PyridazineCC[[3,2]]},{0,PyridazineCC[[3,2]]}};
SCIPT2=Table[{(10^3)(Pyridazine[[k,3]]-Pyridazine[[k,2]]),(10^3)(Pyridazine[[k,2]]-PyridazineHF)},{k,start,Length[Pyridazine]}];
fitPT2=Fit[Take[SCIPT2,-nfit],{1,x},x]
SCIrPT2=Table[{(10^3)(Pyridazine[[k,5]]-Pyridazine[[k,2]]),(10^3)(Pyridazine[[k,2]]-PyridazineHF)},{k,start,Length[Pyridazine]}];
fitrPT2=Fit[Take[SCIrPT2,-nfit],{1,x},x]
Show[{
ListPlot[{
CCSDT,CCSDTQ,SCIPT2,SCIrPT2
},
PlotLabel->MaTeX["\\text{Pyridazine/cc-pVDZ}",FontSize->SizeLabel],
FrameLabel->{
MaTeX["\\text{$E_\\text{PT2}$ or $E_\\text{rPT2}$ ($mE_h$)}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{CCSDT}",FontSize->SizeLegend],
MaTeX["\\text{CCSDTQ}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.}$ vs $E_\\text{PT2}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.}$ vs $E_\\text{rPT2}$}",FontSize->SizeLegend]
}],{Right,Top}],
PlotOptions,
PlotRange->{{-200,0},{-500,-880}},
Joined->False
]
,
Plot[{fitPT2,fitrPT2},{x,-200,0},PlotStyle->{{Thick,Dashed,Blue},{Thick,Dashed,Red}}]
}]
(*Export["furanEcvsPT2.pdf",%];*)

Pyridine

MP & CC calculations

PyridineHF=-246.715184707013;
PyridineMP={{5,-807.471827567},{6,-834.035340573813},{7,-879.923597854822},{8,-872.875093}};
PyridineCCn={{5,-816.328319360},{7,-879.599675234389},{9,-880.821972634781}};
PyridineCC={{6,-840.701382579433},{8,-877.874903251700},{10,-880.821972634781}};

SizeTitle=30;SizeLabel=20;SizeLegend=20;
PlotOptions={
Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1/2,GridLines->Automatic,PlotMarkers->{"OpenMarkers",Large},
PlotStyle->{
{Thickness[0.01],Black},{Thickness[0.01],Darker[Green]},{Thickness[0.01],Red},{Thickness[0.01],Blue}
}
};
ListPlot[{PyridineMP,PyridineCCn,PyridineCC},Joined->True,PlotOptions,
FrameTicks->{{Automatic,None},{{5,6,7,8,9,10},None}},
PlotLabel->MaTeX["\\text{Pyridine/cc-pVDZ}",FontSize->SizeLabel],,
FrameLabel->{
MaTeX["\\text{Computational cost $\\mathcal{O}(N^\\alpha)$}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{MP$n$}",FontSize->SizeLegend],MaTeX["\\text{CC$n$}",FontSize->SizeLegend],
MaTeX["\\text{CC}",FontSize->SizeLegend]
}],{Right,Top}]]

Shanks

s=PyridineMP[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

s=PyridineCCn[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

s=PyridineCC[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

CIPSI  calculations

Pyridine=Import["pyridine.dat"];

SizeTitle=30;SizeLabel=24;SizeLegend=24;
PlotOptions={
Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1,GridLines->Automatic,PlotMarkers->"OpenMarkers",
PlotStyle->{
{Thick,Black},{Thick,Darker[Green]},{Thick,Red},{Thick,Blue}
}
};
start=9;
CCSDT={{0.01,PyridineCC[[2,2]]},{(10^9),PyridineCC[[2,2]]}};
CCSDTQ={{0.01,PyridineCC[[3,2]]},{(10^9),PyridineCC[[3,2]]}};
SCI=Table[{Pyridine[[k,1]],(10^3)(Pyridine[[k,2]]-PyridineHF)},{k,start,Length[Pyridine]}];
SCIPT2=Table[{Pyridine[[k,1]],(10^3)(Pyridine[[k,3]]-PyridineHF)},{k,start,Length[Pyridine]}];
SCIrPT2=Table[{Pyridine[[k,1]],(10^3)(Pyridine[[k,5]]-PyridineHF)},{k,start,Length[Pyridine]}];
ListLogLinearPlot[{
CCSDTQ,SCI,SCIPT2,SCIrPT2
},
FrameLabel->{
MaTeX["\\text{Number of determinants}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{CCSDTQ}",FontSize->SizeLegend],MaTeX["\\text{$\\Delta E_\\text{var.}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.} + E_\\text{PT2}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.} + E_\\text{rPT2}$}",FontSize->SizeLegend]
}],{Right,Top}],
PlotRange->{{(10^3),3 (10^8)},{-880,-550}},
PlotOptions,
Joined->False,
PlotRange->All
]
(*Export["furanEvsNdet.pdf",%];*)

SizeTitle=30;SizeLabel=24;SizeLegend=24;
PlotOptions={Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1,GridLines->Automatic,PlotMarkers->"OpenMarkers",
PlotStyle->{
{Thick,Gray},{Thick,Black},{Thick,Blue},{Thick,Red},{Thick,Dashed,Green}
}
};
start=9;
nfit=4;
CCSDT={{-1000,PyridineCC[[2,2]]},{0,PyridineCC[[2,2]]}};
CCSDTQ={{-1000,PyridineCC[[3,2]]},{0,PyridineCC[[3,2]]}};
SCIPT2=Table[{(10^3)(Pyridine[[k,3]]-Pyridine[[k,2]]),(10^3)(Pyridine[[k,2]]-PyridineHF)},{k,start,Length[Pyridine]}];
fitPT2=Fit[Take[SCIPT2,-nfit],{1,x},x]
SCIrPT2=Table[{(10^3)(Pyridine[[k,5]]-Pyridine[[k,2]]),(10^3)(Pyridine[[k,2]]-PyridineHF)},{k,start,Length[Pyridine]}];
fitrPT2=Fit[Take[SCIrPT2,-nfit],{1,x},x]
Show[{
ListPlot[{
CCSDT,CCSDTQ,SCIPT2,SCIrPT2
},
PlotLabel->MaTeX["\\text{Pyridine/cc-pVDZ}",FontSize->SizeLabel],
FrameLabel->{
MaTeX["\\text{$E_\\text{PT2}$ or $E_\\text{rPT2}$ ($mE_h$)}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{CCSDT}",FontSize->SizeLegend],
MaTeX["\\text{CCSDTQ}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.}$ vs $E_\\text{PT2}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.}$ vs $E_\\text{rPT2}$}",FontSize->SizeLegend]
}],{Right,Top}],
PlotOptions,
PlotRange->{{-200,0},{-500,-880}},
Joined->False
]
,
Plot[{fitPT2,fitrPT2},{x,-200,0},PlotStyle->{{Thick,Dashed,Blue},{Thick,Dashed,Red}}]
}]
(*Export["furanEcvsPT2.pdf",%];*)

Pyrimidine

MP & CC calculations

PyrimidineHF=-262.713652796899;
PyrimidineMP={{5,-830.057493278},{6,-849.612212994749},{7,-899.266633591386},{8,-889.9912031}};
PyrimidineCCn={{5,-841.308837134},{7,-898.365293520950},{9,-899.263868742531}};
PyrimidineCC={{6,-857.976308572695},{8,-896.080430191809},{10,-898.976078520330}};

SizeTitle=30;SizeLabel=20;SizeLegend=20;
PlotOptions={
Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1/2,GridLines->Automatic,PlotMarkers->{"OpenMarkers",Large},
PlotStyle->{
{Thickness[0.01],Black},{Thickness[0.01],Darker[Green]},{Thickness[0.01],Red},{Thickness[0.01],Blue}
}
};
ListPlot[{PyrimidineMP,PyrimidineCCn,PyrimidineCC},Joined->True,PlotOptions,
FrameTicks->{{Automatic,None},{{5,6,7,8,9,10},None}},
PlotLabel->MaTeX["\\text{Pyrimidine/cc-pVDZ}",FontSize->SizeLabel],
FrameLabel->{
MaTeX["\\text{Computational cost $\\mathcal{O}(N^\\alpha)$}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{MP$n$}",FontSize->SizeLegend],MaTeX["\\text{CC$n$}",FontSize->SizeLegend],
MaTeX["\\text{CC}",FontSize->SizeLegend]
}],{Right,Top}]]

Shanks

s=PyrimidineMP[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

s=PyrimidineCCn[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

s=PyrimidineCC[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

CIPSI  calculations

Pyrimidine=Import["pyrimidine.dat"];

SizeTitle=30;SizeLabel=24;SizeLegend=24;
PlotOptions={
Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1,GridLines->Automatic,PlotMarkers->"OpenMarkers",
PlotStyle->{
{Thick,Black},{Thick,Darker[Green]},{Thick,Red},{Thick,Blue}
}
};
start=9;
CCSDT={{0.01,PyrimidineCC[[2,2]]},{(10^9),PyrimidineCC[[2,2]]}};
CCSDTQ={{0.01,PyrimidineCC[[3,2]]},{(10^9),PyrimidineCC[[3,2]]}};
SCI=Table[{Pyrimidine[[k,1]],(10^3)(Pyrimidine[[k,2]]-PyrimidineHF)},{k,start,Length[Pyrimidine]}];
SCIPT2=Table[{Pyrimidine[[k,1]],(10^3)(Pyrimidine[[k,3]]-PyrimidineHF)},{k,start,Length[Pyrimidine]}];
SCIrPT2=Table[{Pyrimidine[[k,1]],(10^3)(Pyrimidine[[k,5]]-PyrimidineHF)},{k,start,Length[Pyrimidine]}];
ListLogLinearPlot[{
CCSDTQ,SCI,SCIPT2,SCIrPT2
},
FrameLabel->{
MaTeX["\\text{Number of determinants}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{CCSDTQ}",FontSize->SizeLegend],MaTeX["\\text{$\\Delta E_\\text{var.}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.} + E_\\text{PT2}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.} + E_\\text{rPT2}$}",FontSize->SizeLegend]
}],{Right,Top}],
PlotRange->{{(10^3),3 (10^8)},{-880,-550}},
PlotOptions,
Joined->False,
PlotRange->All
]
(*Export["furanEvsNdet.pdf",%];*)

SizeTitle=30;SizeLabel=24;SizeLegend=24;
PlotOptions={Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1,GridLines->Automatic,PlotMarkers->"OpenMarkers",
PlotStyle->{
{Thick,Gray},{Thick,Black},{Thick,Blue},{Thick,Red},{Thick,Dashed,Green}
}
};
start=9;
nfit=4;
CCSDT={{-1000,PyrimidineCC[[2,2]]},{0,PyrimidineCC[[2,2]]}};
CCSDTQ={{-1000,PyrimidineCC[[3,2]]},{0,PyrimidineCC[[3,2]]}};
SCIPT2=Table[{(10^3)(Pyrimidine[[k,3]]-Pyrimidine[[k,2]]),(10^3)(Pyrimidine[[k,2]]-PyrimidineHF)},{k,start,Length[Pyrimidine]}];
fitPT2=Fit[Take[SCIPT2,-nfit],{1,x},x]
SCIrPT2=Table[{(10^3)(Pyrimidine[[k,5]]-Pyrimidine[[k,2]]),(10^3)(Pyrimidine[[k,2]]-PyrimidineHF)},{k,start,Length[Pyrimidine]}];
fitrPT2=Fit[Take[SCIrPT2,-nfit],{1,x},x]
Show[{
ListPlot[{
CCSDT,CCSDTQ,SCIPT2,SCIrPT2
},
PlotLabel->MaTeX["\\text{Pyrimidine/cc-pVDZ}",FontSize->SizeLabel],
FrameLabel->{
MaTeX["\\text{$E_\\text{PT2}$ or $E_\\text{rPT2}$ ($mE_h$)}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{CCSDT}",FontSize->SizeLegend],
MaTeX["\\text{CCSDTQ}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.}$ vs $E_\\text{PT2}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.}$ vs $E_\\text{rPT2}$}",FontSize->SizeLegend]
}],{Right,Top}],
PlotOptions,
PlotRange->{{-200,0},{-500,-880}},
Joined->False
]
,
Plot[{fitPT2,fitrPT2},{x,-200,0},PlotStyle->{{Thick,Dashed,Blue},{Thick,Dashed,Red}}]
}]
(*Export["furanEcvsPT2.pdf",%];*)

Tetrazine

MP & CC calculations

TetrazineHF=-294.615743073901;
TetrazineMP={{5,-895.922989130},{6,-899.456957539510},{7,-958.570330486010},{8,-944.2694261}};
TetrazineCCn={{5,-908.969490669},{7,-954.853218744419},{9,-955.880388848436}};
TetrazineCC={{6,-909.100948888639},{8,-951.585113405389},{10,-955.428165912650}};

SizeTitle=30;SizeLabel=20;SizeLegend=20;
PlotOptions={
Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1/2,GridLines->Automatic,PlotMarkers->{"OpenMarkers",Large},
PlotStyle->{
{Thickness[0.01],Black},{Thickness[0.01],Darker[Green]},{Thickness[0.01],Red},{Thickness[0.01],Blue}
}
};
ListPlot[{TetrazineMP,TetrazineCCn,TetrazineCC},Joined->True,PlotOptions,
FrameTicks->{{Automatic,None},{{5,6,7,8,9,10},None}},
PlotLabel->MaTeX["\\text{Tetrazine/cc-pVDZ}",FontSize->SizeLabel],
FrameLabel->{
MaTeX["\\text{Computational cost $\\mathcal{O}(N^\\alpha)$}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{MP$n$}",FontSize->SizeLegend],MaTeX["\\text{CC$n$}",FontSize->SizeLegend],
MaTeX["\\text{CC}",FontSize->SizeLegend]
}],{Right,Top}]]

Shanks

s=TetrazineMP[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

s=TetrazineCCn[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

s=TetrazineCC[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

CIPSI  calculations

Tetrazine=Import["tetrazine.dat"];

SizeTitle=30;SizeLabel=24;SizeLegend=24;
PlotOptions={
Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1,GridLines->Automatic,PlotMarkers->"OpenMarkers",
PlotStyle->{
{Thick,Black},{Thick,Darker[Green]},{Thick,Red},{Thick,Blue}
}
};
start=9;
CCSDT={{0.01,TetrazineCC[[2,2]]},{(10^9),TetrazineCC[[2,2]]}};
CCSDTQ={{0.01,TetrazineCC[[3,2]]},{(10^9),TetrazineCC[[3,2]]}};
SCI=Table[{Tetrazine[[k,1]],(10^3)(Tetrazine[[k,2]]-TetrazineHF)},{k,start,Length[Tetrazine]}];
SCIPT2=Table[{Tetrazine[[k,1]],(10^3)(Tetrazine[[k,3]]-TetrazineHF)},{k,start,Length[Tetrazine]}];
SCIrPT2=Table[{Tetrazine[[k,1]],(10^3)(Tetrazine[[k,5]]-TetrazineHF)},{k,start,Length[Tetrazine]}];
ListLogLinearPlot[{
CCSDTQ,SCI,SCIPT2,SCIrPT2
},
FrameLabel->{
MaTeX["\\text{Number of determinants}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{CCSDTQ}",FontSize->SizeLegend],MaTeX["\\text{$\\Delta E_\\text{var.}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.} + E_\\text{PT2}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.} + E_\\text{rPT2}$}",FontSize->SizeLegend]
}],{Right,Top}],
PlotRange->{{(10^3),3 (10^8)},{-880,-550}},
PlotOptions,
Joined->False,
PlotRange->All
]
(*Export["furanEvsNdet.pdf",%];*)

SizeTitle=30;SizeLabel=24;SizeLegend=24;
PlotOptions={Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1,GridLines->Automatic,PlotMarkers->"OpenMarkers",
PlotStyle->{
{Thick,Gray},{Thick,Black},{Thick,Blue},{Thick,Red},{Thick,Dashed,Green}
}
};
start=9;
nfit=4;
CCSDT={{-1000,TetrazineCC[[2,2]]},{0,TetrazineCC[[2,2]]}};
CCSDTQ={{-1000,TetrazineCC[[3,2]]},{0,TetrazineCC[[3,2]]}};
SCIPT2=Table[{(10^3)(Tetrazine[[k,3]]-Tetrazine[[k,2]]),(10^3)(Tetrazine[[k,2]]-TetrazineHF)},{k,start,Length[Tetrazine]}];
fitPT2=Fit[Take[SCIPT2,-nfit],{1,x},x]
SCIrPT2=Table[{(10^3)(Tetrazine[[k,5]]-Tetrazine[[k,2]]),(10^3)(Tetrazine[[k,2]]-TetrazineHF)},{k,start,Length[Tetrazine]}];
fitrPT2=Fit[Take[SCIrPT2,-nfit],{1,x},x]
Show[{
ListPlot[{
CCSDT,CCSDTQ,SCIPT2,SCIrPT2
},
PlotLabel->MaTeX["\\text{Tetrazine/cc-pVDZ}",FontSize->SizeLabel],
FrameLabel->{
MaTeX["\\text{$E_\\text{PT2}$ or $E_\\text{rPT2}$ ($mE_h$)}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{CCSDT}",FontSize->SizeLegend],
MaTeX["\\text{CCSDTQ}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.}$ vs $E_\\text{PT2}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.}$ vs $E_\\text{rPT2}$}",FontSize->SizeLegend]
}],{Right,Top}],
PlotOptions,
PlotRange->{{-200,0},{-500,-880}},
Joined->False
]
,
Plot[{fitPT2,fitrPT2},{x,-200,0},PlotStyle->{{Thick,Dashed,Blue},{Thick,Dashed,Red}}]
}]
(*Export["furanEcvsPT2.pdf",%];*)

Triazine

MP & CC calculations

TriazineHF=-278.717319389695;
TriazineMP={{5,-850.519282476},{6,-863.593139716035},{7,-916.717801063422},{8,-905.4360103}};
TriazineCCn={{5,-864.406418639},{7,-915.583810723581},{9,-916.075626870853}};
TriazineCC={{6,-873.819918012385},{8,-912.716493155011},{10,-915.770573661020}};

SizeTitle=30;SizeLabel=20;SizeLegend=20;
PlotOptions={
Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1/2,GridLines->Automatic,PlotMarkers->{"OpenMarkers",Large},
PlotStyle->{
{Thickness[0.01],Black},{Thickness[0.01],Darker[Green]},{Thickness[0.01],Red},{Thickness[0.01],Blue}
}
};
ListPlot[{TriazineMP,TriazineCCn,TriazineCC},Joined->True,PlotOptions,
FrameTicks->{{Automatic,None},{{5,6,7,8,9,10},None}},
PlotLabel->MaTeX["\\text{Triazine/cc-pVDZ}",FontSize->SizeLabel],
FrameLabel->{
MaTeX["\\text{Computational cost $\\mathcal{O}(N^\\alpha)$}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{MP$n$}",FontSize->SizeLegend],MaTeX["\\text{CC$n$}",FontSize->SizeLegend],
MaTeX["\\text{CC}",FontSize->SizeLegend]
}],{Right,Top}]]

Shanks

s=TriazineMP[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

s=TriazineCCn[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

s=TriazineCC[[;;,2]]
s=Table[((s[[n+1]]s[[n-1]]-((s[[n]])^2))/(s[[n+1]]-2s[[n]]+s[[n-1]])),{n,2,Length[s]-1}]

CIPSI  calculations

Triazine=Import["triazine.dat"];

SizeTitle=30;SizeLabel=24;SizeLegend=24;
PlotOptions={
Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1,GridLines->Automatic,PlotMarkers->"OpenMarkers",
PlotStyle->{
{Thick,Black},{Thick,Darker[Green]},{Thick,Red},{Thick,Blue}
}
};
start=9;
CCSDT={{0.01,TriazineCC[[2,2]]},{(10^9),TriazineCC[[2,2]]}};
CCSDTQ={{0.01,TriazineCC[[3,2]]},{(10^9),TriazineCC[[3,2]]}};
SCI=Table[{Triazine[[k,1]],(10^3)(Triazine[[k,2]]-TriazineHF)},{k,start,Length[Triazine]}];
SCIPT2=Table[{Triazine[[k,1]],(10^3)(Triazine[[k,3]]-TriazineHF)},{k,start,Length[Triazine]}];
SCIrPT2=Table[{Triazine[[k,1]],(10^3)(Triazine[[k,5]]-TriazineHF)},{k,start,Length[Triazine]}];
ListLogLinearPlot[{
CCSDTQ,SCI,SCIPT2,SCIrPT2
},
FrameLabel->{
MaTeX["\\text{Number of determinants}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{CCSDTQ}",FontSize->SizeLegend],MaTeX["\\text{$\\Delta E_\\text{var.}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.} + E_\\text{PT2}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.} + E_\\text{rPT2}$}",FontSize->SizeLegend]
}],{Right,Top}],
PlotRange->{{(10^3),3 (10^8)},{-880,-550}},
PlotOptions,
Joined->False,
PlotRange->All
]
(*Export["furanEvsNdet.pdf",%];*)

SizeTitle=30;SizeLabel=24;SizeLegend=24;
PlotOptions={Joined->True,PlotTheme->"Detailed",BaseStyle->18,ImageSize->Large,AspectRatio->1,GridLines->Automatic,PlotMarkers->"OpenMarkers",
PlotStyle->{
{Thick,Gray},{Thick,Black},{Thick,Blue},{Thick,Red},{Thick,Dashed,Green}
}
};
start=9;
nfit=4;
CCSDT={{-1000,TriazineCC[[2,2]]},{0,TriazineCC[[2,2]]}};
CCSDTQ={{-1000,TriazineCC[[3,2]]},{0,TriazineCC[[3,2]]}};
SCIPT2=Table[{(10^3)(Triazine[[k,3]]-Triazine[[k,2]]),(10^3)(Triazine[[k,2]]-TriazineHF)},{k,start,Length[Triazine]}];
fitPT2=Fit[Take[SCIPT2,-nfit],{1,x},x]
SCIrPT2=Table[{(10^3)(Triazine[[k,5]]-Triazine[[k,2]]),(10^3)(Triazine[[k,2]]-TriazineHF)},{k,start,Length[Triazine]}];
fitrPT2=Fit[Take[SCIrPT2,-nfit],{1,x},x]
Show[{
ListPlot[{
CCSDT,CCSDTQ,SCIPT2,SCIrPT2
},
PlotLabel->MaTeX["\\text{Triazine/cc-pVDZ}",FontSize->SizeLabel],
FrameLabel->{
MaTeX["\\text{$E_\\text{PT2}$ or $E_\\text{rPT2}$ ($mE_h$)}",FontSize->SizeLabel],
MaTeX["\\text{Correlation energy (m$E_h$)}",FontSize->SizeLabel]},
PlotLegends->Placed[LineLegend[{
MaTeX["\\text{CCSDT}",FontSize->SizeLegend],
MaTeX["\\text{CCSDTQ}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.}$ vs $E_\\text{PT2}$}",FontSize->SizeLegend],
MaTeX["\\text{$\\Delta E_\\text{var.}$ vs $E_\\text{rPT2}$}",FontSize->SizeLegend]
}],{Right,Top}],
PlotOptions,
PlotRange->{{-200,0},{-500,-880}},
Joined->False
]
,
Plot[{fitPT2,fitrPT2},{x,-200,0},PlotStyle->{{Thick,Dashed,Blue},{Thick,Dashed,Red}}]
}]
(*Export["furanEcvsPT2.pdf",%];*)