from ase.structure import molecule
from gpaw import GPAW

system = molecule('H2O')
system.center(vacuum=4.0)

calc = GPAW(mode='pw', parallel=dict(band=4), nbands=8, basis='dzp',
            #dtype=complex
            )
system.set_calculator(calc)

system.get_potential_energy()