# Reading total energies energies
e1 = nlread('Si-alpha-100.hdf5', TotalEnergy)[0]
e2 = nlread('SiH4.hdf5', TotalEnergy)[0]

e3 = nlread('SiH4-Si.hdf5', TotalEnergy)[0]
e4 = nlread('SiH3-SiH.hdf5', TotalEnergy)[0]
e5 = nlread('SiH2-2SiH.hdf5', TotalEnergy)[0]

# computing adsorption energy
E3 = e3.evaluate().inUnitsOf(eV) - e1.evaluate().inUnitsOf(eV) - e2.evaluate().inUnitsOf(eV)
print("Adsorption energy SiH4 on Si(100) = %.2f eV" % E3)

E4 = e4.evaluate().inUnitsOf(eV) - e1.evaluate().inUnitsOf(eV) - e2.evaluate().inUnitsOf(eV)
print("Adsorption energy of SiH3 and H on Si(100) = %.2f eV" % E4)

E5 = e5.evaluate().inUnitsOf(eV) - e1.evaluate().inUnitsOf(eV) - e2.evaluate().inUnitsOf(eV)
print("Adsorption energy of SiH2 and 2 H on Si(100) = %.2f eV" % E5)