Needs["MaTeX`"];
SetOptions[MaTeX,"Preamble"->{"\\usepackage{amssymb,amsmath,latexsym,amsfonts,amsthm,mathpazo,xcolor,bm,mhchem}"}];
SetDirectory["/Users/monino/Work/bse/sfBSE/output"];
coord=Import["seq","List"];
sfBSE={Import["sf_bse_first_excitation","List"],Import["sf_bse_second_excitation","List"],Import["sf_bse_third_excitation","List"],Import["sf_bse_fourth_excitation","List"],Import["sf_bse_sixth_excitation","List"],Import["sf_bse_eighth_excitation","List"]} ;

sfTDDFT={Import["sf_tddft_first_excitation","List"],Import["sf_tddft_second_excitation","List"],Import["sf_tddft_third_excitation","List"],Import["sf_tddft_fourth_excitation","List"],Import["sf_tddft_fifth_excitation","List"],Import["sf_tddft_sixth_excitation","List"],Import["sf_tddft_seventh_excitation","List"],Import["sf_tddft_eighth_excitation","List"]} ;

sfCIS={Import["sf_cis_first_excitation","List"],Import["sf_cis_second_excitation","List"],Import["sf_cis_third_excitation","List"],Import["sf_cis_fourth_excitation","List"],Import["sf_cis_fifth_excitation","List"],Import["sf_cis_sixth_excitation","List"],Import["sf_cis_seventh_excitation","List"],Import["sf_cis_eighth_excitation","List"]} ;

eomCCSD= {Import["eom_cssd_first_excitation","List"],Import["eom_cssd_second_excitation","List"],Import["eom_cssd_third_excitation","List"],Import["eom_cssd_sixth_excitation","List"]} ;

CIS= {Import["cis_first_excitation","List"],Import["cis_second_excitation","List"],Import["cis_third_excitation","List"],Import["cis_fourth_excitation","List"],Import["cis_fifth_excitation","List"],Import["cis_sixth_excitation","List"],Import["cis_seventh_excitation","List"]} ;

TDDFT= {Import["tddft_first_excitation","List"],Import["tddft_second_excitation","List"],Import["tddft_third_excitation","List"],Import["tddft_fourth_excitation","List"],Import["tddft_fifth_excitation","List"],Import["tddft_sixth_excitation","List"],Import["tddft_seventh_excitation","List"]} ;

BSE={Import["bse_first_excitation","List"],Import["bse_second_excitation","List"],Import["bse_third_excitation","List"],Import["bse_fourth_excitation","List"],Import["bse_fifth_excitation","List"],Import["bse_sixth_excitation","List"],Import["bse_seventh_excitation","List"]} ;

plot=Table[Transpose@{coord,sfBSE[[k+1]]-sfBSE[[1]]},{k,1,5}];
plot1=Table[Transpose@{coord,sfCIS[[k+1]]-sfCIS[[1]]},{k,1,7}];
plot2=Table[Transpose@{coord,eomCCSD[[k]]},{k,1,4}];
plot3=Table[Transpose@{coord,sfTDDFT[[k+1]]-sfTDDFT[[1]]},{k,1,7}];
plot4=Table[Transpose@{coord,CIS[[k]]},{k,1,7}];
plot5=Table[Transpose@{coord,BSE[[k]]},{k,1,7}];
plot6 = Table[Transpose@{coord,CIS[[k]]},{k,1,7}];

Show[{ListPlot[plot,Joined->True,PlotStyle-> Blue,PlotLegends-> {MaTeX["\\text{SF-BSE@G}_0\\text{W}_0",FontSize->24]}],ListPlot[plot2,Joined->True,PlotStyle-> Black,PlotLegends-> {MaTeX["\\text{EOM-CCSD}",FontSize->24]}]},Frame->True,PlotLabel-> MaTeX["\\ce{H2} ~/~ \\text{cc-pVTZ}",FontSize->24],LabelStyle->Directive[Black],FrameLabel-> {MaTeX["R \\text{(H-H)  (\\AA)}",FontSize->24],MaTeX["\\Delta E \\text{ (H-H)  (eV)}",FontSize->24]},FrameStyle->Directive[Gray,Medium],ImageSize->Large]

Show[{ListPlot[plot1,Joined->True,PlotMarkers->{Automatic,5},PlotStyle-> Red,PlotLegends-> {MaTeX["\\text{SF-CIS}"]}],ListPlot[plot2,Joined->True,PlotMarkers->{Automatic,5},PlotStyle-> Black,PlotLegends-> {MaTeX["\\text{EOM-CCSD}"]}]},Frame->True,PlotLabel-> MaTeX["\\text{H}_2 / \\text{CC-PVTZ}"],LabelStyle->Directive[Black],FrameLabel-> {MaTeX["R \\text{(H-H)  (angstrom)}"],MaTeX["E \\text{(H-H)  (eV)}"]},FrameStyle->Directive[Black,Medium]]

Show[{ListPlot[plot3,Joined->True,PlotMarkers->{Automatic,5},PlotStyle-> Purple,PlotLegends-> {MaTeX["\\text{SF-TDDFT}"]}],ListPlot[plot2,Joined->True,PlotMarkers->{Automatic,5},PlotStyle-> Black,PlotLegends-> {MaTeX["\\text{EOM-CCSD}"]}]},Frame->True,PlotLabel-> MaTeX["\\text{H}_2 / \\text{CC-PVTZ}"],LabelStyle->Directive[Black],FrameLabel-> {MaTeX["R \\text{(H-H)  (angstrom)}"],MaTeX["E \\text{(H-H)  (eV)}"]},FrameStyle->Directive[Black,Medium]]

Show[{ListPlot[plot5,Joined->True,PlotMarkers->{Automatic,5},PlotStyle-> Blue,PlotLegends-> {MaTeX["\\text{BSE@G0W0}"]}],ListPlot[plot2,Joined->True,PlotMarkers->{Automatic,5},PlotStyle-> Black,PlotLegends-> {MaTeX["\\text{EOM-CCSD}"]}]},Frame->True,PlotLabel-> MaTeX["\\text{H}_2 / \\text{CC-PVTZ}"],LabelStyle->Directive[Black],FrameLabel-> {MaTeX["R \\text{(H-H)  (angstrom)}"],MaTeX["E \\text{(H-H)  (eV)}"]},FrameStyle->Directive[Black,Medium]]

Show[{ListPlot[plot4,Joined->True,PlotMarkers->{Automatic,5},PlotStyle-> Red,PlotLegends-> {MaTeX["\\text{CIS}"]}],ListPlot[plot2,Joined->True,PlotMarkers->{Automatic,5},PlotStyle-> Black,PlotLegends-> {MaTeX["\\text{EOM-CCSD}"]}]},Frame->True,PlotLabel-> MaTeX["\\text{H}_2 / \\text{CC-PVTZ}"],LabelStyle->Directive[Black],FrameLabel-> {MaTeX["R \\text{(H-H)  (angstrom)}"],MaTeX["E \\text{(H-H)  (eV)}"]},FrameStyle->Directive[Black,Medium]]

test1 = ListPlot[Table[Transpose@{coord,sfBSE[[k+1]]-sfBSE[[1]]},{k,1,5}],Joined->False,PlotStyle->{{Blue,Thin,DotDashed}},PlotLegends-> {MaTeX["\\text{SF-BSE@G}_0\\text{W}_0",FontSize->24]}];
test2 = ListPlot[Table[Transpose@{coord,eomCCSD[[k]]},{k,1,4}],Joined->True,PlotStyle->Black,PlotLegends-> {MaTeX["\\text{EOM-CCSD}",FontSize->24]}];
test3 = ListPlot[Table[Transpose@{coord,BSE[[k]]},{k,1,7}],Joined->True,PlotStyle->Cyan,PlotLegends-> {MaTeX["\\text{BSE}",FontSize->24]}];
Show[test1,test2,test3,Frame->True,PlotLabel-> MaTeX["\\ce{H2} ~/~ \\text{cc-pVTZ}",FontSize->24],LabelStyle->Directive[Black],FrameLabel-> {MaTeX["R \\text{(H-H)  (\\AA)}",FontSize->24],MaTeX["\\Delta E \\text{ (H-H)  (eV)}",FontSize->24]},FrameStyle->Directive[Gray,Medium],ImageSize->Full]

test={"B\!\(\*SuperscriptBox[\(\\ \), \(1\)]\)" Superscript[Subscript[Σ,u],"+"],"EF\!\(\*SuperscriptBox[\(\\ \), \(1\)]\)" Superscript[Subscript[Σ,g],"+"],"EF\!\(\*SuperscriptBox[\(\\ \), \(1\)]\)" Superscript[Subscript[Σ,g],"+"],""};
test4 = ListPlot[Table[Callout[Table[{coord[[k]],eomCCSD[[l,k]]},{k,1,51}],test[[l]]],{l,4}],Joined->True,PlotStyle->Black,PlotLegends-> {MaTeX["\\text{EOM-CCSD}",FontSize->24]}];
Length[BSE[[1,;;]]];
Show[test4,test3,test1,Frame->True,PlotLabel-> MaTeX["\\ce{H2} ~/~ \\text{cc-pVTZ}",FontSize->24],LabelStyle->Directive[Black],FrameLabel-> {MaTeX["R \\text{(H-H)  (\\AA)}",FontSize->24],MaTeX["\\Delta E \\text{ (H-H)  (eV)}",FontSize->24]},FrameStyle->Directive[Gray,Medium],BaseStyle->20,ImageSize->Full,PlotRange->All]