Contents:
Get HF and orbital energies from a QCHEM output file
Return the HF energy of a QChem calculation
Raise exception if HF energy not found
Run Qchen for given molecular coordinates.
X_au: 3N array, molecular coordinates in bohr
if state == 1, EHF - e1 …