Right, no, I don't think there is anything wrong with the script or so. It's just not clear to me exactly how this "0 eV" appears. Is it an eigenvalue of the MPSH spectrum? How does it look, exactly, in the output log?hi, Blom, I am glad that you give me reply rapidly.
if processIsMaster(): nlPrint(projected_hamiltonian_energy_spectrum)
if processIsMaster():
print "# MPSH Spectrum (eV)"
print "# -----------------------------------------------------------------------------"
for E in projected_hamiltonian_energy_spectrum.energies():
print E.inUnitsOf(eV)
I see nothing wrong here. You happen to have two molecular levels that coincide with the average Fermi level, but if there were more decimals I'm sure you would not have 0.000000000 eV but perhaps 0.0043 eV or so. They might not even be degenerate, but only so to within 2 decimals.Dear Anders Blom ,
To print the spectrum with more decimals (nlPrint doesn't do a supergood job there...), replace the lineCodeif processIsMaster(): nlPrint(projected_hamiltonian_energy_spectrum)
byCodeif processIsMaster():
print "# MPSH Spectrum (eV)"
print "# -----------------------------------------------------------------------------"
for E in projected_hamiltonian_energy_spectrum.energies():
print E.inUnitsOf(eV)
See http://quantumwise.com/forum/index.php?topic=103.0Thank you so much.