SetDirectory["/Users/monino/Work/bse/sfBSE/output"];

data=Import["h2_qchem_energy_sf_cis.txt","List"];
data1 = Import["h2_qchem_energy_sf_cis_double_exc.txt","List"];
data2=Import["h2_qchem_energy_triplet_cis.txt","List"] (*-2 *(-0.4998099076)*);
newdata = data (*-2 *(-0.4998099076)*);
newdata1 = data1  (*-2 *(-0.4998099076)*);
coord = Import["seq","List"]*1.89;
coord1 = Import["seq","List"]* 1.89;
plot = Transpose@{coord,newdata};(*or using Thread[{coord,newdata}]*)
plot1 = Transpose@{coord1,newdata1};
plot2 = Transpose@{coord1,data2};
ListPlot[{Labeled[plot,"GS"],Labeled[plot1,"Double excitation"],Labeled[plot2,"Simple excitation"]},Joined->True,AxesLabel->"E (au)"]