from qucochemistry.vqe import VQEexperiment from openfermion.hamiltonians import MolecularData
filename = 'molecules/H2_pyscf_equi.hdf5'
molecule = MolecularData(filename=filename)
vqe = VQEexperiment(molecule=molecule)
E = vqe.get_exact_gs()
vqe.start_vqe()
result = vqe.get_results()
print('Difference between VQE optimized and exact GS energy:')
print(str((result.fun-E)/0.0016) + 'kcal/mol')