from ase import Atoms
from gpaw import GPAW

system = Atoms('O')
system.center(vacuum=4.0)

calc = GPAW(hund=True, xc='PBE')

system.set_calculator(calc)
system.get_potential_energy()