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

cbutadiene = {
{"\!\(\*SuperscriptBox[\(1\), \(\(\\ \\ \)\(3\)\)]\)(Subscript[A, 2  g])", 0.317,0.549,0.265,0.217,0.083,0.404},
{"\!\(\*SuperscriptBox[\(2\), \(\(\\ \\ \)\(1\)\)]\)(Subscript[A, 1  g])",2.650,1.783,1.658,1.123,1.621,1.685},
{"\!\(\*SuperscriptBox[\(1\), \(\(\\ \\ \)\(1\)\)]\)(Subscript[B, 2  g])",3.125,1.966,1.904,1.799,1.930,1.976}
};

ΔΕ=Table[cbutadiene[[;;,k]],{k,2,7}];
state=cbutadiene[[;;,1]];
method={"SF-CIS","SF-TDDFT","SF-ADC(2)-s","SF-ADC(2)-x","SF-ADC(3)","SF-BSE"};
color={Red,Blue,Green,Purple,Black,Red};
Show[{
ListPlot[Table[Callout[Table[{2m-3/2,ΔΕ[[m,k]]},{m,Length[ΔΕ]}],state[[k]],After],{k,3}],Joined->True,PlotStyle->{{LightGreen,Thickness[0.01]},{LightOrange,Thickness[0.01]},{LightRed,Thickness[0.01]},{LightPurple,Thickness[0.01]}}]
,
Show[
Table[Plot[If[2m-2≤x≤2m-1,ΔΕ[[m]]],{x,0,2Length[ΔΕ]-1},PlotStyle->Directive[Thick,color[[m]]]],{m,Length[ΔΕ]}]
]
},
PlotLabel->MaTeX["\\text{\\ce{C4H4}/cc-pVTZ (D4h)}",FontSize->30],
Frame->True,Axes->False,PlotRange->{All,{0.2,3}},GridLines->{{5.5,13.5,21.5},None},FrameTicks->{{Automatic,None},{Table[{2k-1,Rotate[MaTeX["\\text{"<>method[[k]]<>"}",FontSize->20],π/2]},{k,Length[method]}],None}},
BaseStyle->20,PlotRangePadding->0.5,FrameStyle->Directive[Thick,20,Black],AspectRatio->1/2,FrameLabel->{None,MaTeX["\\text{Excitation energy (eV)}",FontSize->20]},ImageSize->800]

cbutadieneD2h = {
{"\!\(\*SuperscriptBox[\(1\), \(\(\\ \\ \)\(3\)\)]\)(Subscript[B, 1  g])", 0.0,1.654,1.573,1.576,1.456},
{"\!\(\*SuperscriptBox[\(1\), \(\(\\ \\ \)\(1\)\)]\)(Subscript[B, 1  g])",0.0,3.416,3.208,3.141,3.285},
{"\!\(\*SuperscriptBox[\(2\), \(\(\\ \\ \)\(1\)\)]\)(Subscript[A, 1  g])",0.0,4.360,4.247,3.796,4.334}
};

ΔΕ=Table[cbutadieneD2h[[;;,k]],{k,2,6}];
state=cbutadieneD2h[[;;,1]];
method={"SF-CIS","SF-CCSD","SF-ADC(2)-s","SF-ADC(2)-x","SF-ADC(3)"};
color={Red,Blue,Green,Purple,Black};
Show[{
ListPlot[Table[Callout[Table[{2m-3/2,ΔΕ[[m,k]]},{m,Length[ΔΕ]}],state[[k]],After,LeaderSize->20],{k,3}],Joined->True,PlotStyle->{{LightGreen,Thickness[0.01]},{LightOrange,Thickness[0.01]},{LightRed,Thickness[0.01]},{LightPurple,Thickness[0.01]}}]
,
Show[
Table[Plot[If[2m-2≤x≤2m-1,ΔΕ[[m]]],{x,0,2Length[ΔΕ]-1},PlotStyle->Directive[Thick,color[[m]]]],{m,Length[ΔΕ]}]
]
},
PlotLabel->MaTeX["\\text{\\ce{C4H4}/cc-pVTZ (D2h)}",FontSize->30],
Frame->True,Axes->False,PlotRange->{All,{0,4.5}},GridLines->{{5.5,13.5,21.5},None},FrameTicks->{{Automatic,None},{Table[{2k-1,Rotate[MaTeX["\\text{"<>method[[k]]<>"}",FontSize->20],π/2]},{k,Length[method]}],None}},
BaseStyle->20,PlotRangePadding->0.5,FrameStyle->Directive[Thick,20,Black],AspectRatio->1/2,FrameLabel->{None,MaTeX["\\text{Excitation energy (eV)}",FontSize->20]},ImageSize->800]