from ase import Atom, Atoms
from ase.calculators.jacapo import Jacapo
molecule1 = Atoms([Atom('C'), Atom('O', [1.2, 0, 0])], cell=(6, 6, 6), pbc = True)
solver1 = Jacapo(nbands = 6, atoms = molecule1)   
print solver1.get_potential_energy()