Δ=0.55 (*eV*)

γ = 0.370; (*eV*)
ωe = 2.15; (*eV*)

Ω=Range[0.02,0.4,0.01]; (*eV*)

Do[{
ω0 = 2.14, 
ω1m = 1.45+0.3,
ω1p = 2.76-0.3,

κ0 =2*0.0451,
κm = 2*0.0189-κ0,
κp = 2*0.0451-κ0,
 
δ = ω0-ωe,
Δm =ω1m-ω0,
Δp =ω1p-ω0,

M = (({{Δm -I κm/2, 0, 0, Ω[[i]]}, {0, 0, 0, Ω[[i]]}, {0, 0,  Δp -I κp/2, Ω[[i]]}, {Ω[[i]], Ω[[i]],  Ω[[i]], - δ +I(κ0-γ)/2}})),{vals,vecs}=Eigensystem[M],esorted=Sort[vals,Re[#1]<Re[#2]&],vecsorted[i]=SortBy[Transpose[{Re[vals],vecs}],First],vecsorted[i]=Flatten[vecsorted[i],1]
},{i,1,Length[Ω]}]

ListLinePlot[{
Table[{Ω[[i]],Abs[vecsorted[i][[4,1]]]^2},{i,1,Length[Ω]}],
Table[{Ω[[i]],Abs[vecsorted[i][[4,2]]]^2},{i,1,Length[Ω]}],
Table[{Ω[[i]],Abs[vecsorted[i][[4,3]]]^2},{i,1,Length[Ω]}],
Table[{Ω[[i]],Abs[vecsorted[i][[4,4]]]^2},{i,1,Length[Ω]}]},
Frame->True,PlotLabels->{"-1","0","+1","exc"},FrameLabel->{"Bare Rabi Coupling (eV)","Eigenstate Components"},PlotLabel->"FSR = "<>ToString[Δ]<>" eV"]