from ASE import Atom, ListOfAtoms
x = ListOfAtoms([Atom('C'), Atom('O', (0, 0, 1.0))], cell = (5, 5, 5))
from Dacapo import Dacapo
y = Dacapo(nbands = 6)
x.SetCalculator(y)
y.Calculate()

a0 = 5.0
box0 = (a0, a0, a0)
r1 = (0.0, 0.0, 0.0) 
r2 = (0.0, 0.0, 1.0)
 
from ASE import Atom
atom1 = Atom('C', r1)
atom2 = Atom('O', r2)
molecule0 = [atom1, atom2]
 
from ASE import ListOfAtoms
model0 = ListOfAtoms(molecule0, cell = box0)
 
from Dacapo import Dacapo
method0 = Dacapo(nbands = 6)
 
model0.SetCalculator(method0)
method0.Calculate()

a0 = 5.0
b0 = 1.1
 
a1 = (a0 , 0.0, 0.0)
a2 = (0.0, a0 , 0.0)
a3 = (0.0, 0.0, a0 )
 
box0 = [a1, a2, a3]
 
r1 = (0.0, 0.0, 0.0)
r2 = (0.0, 0.0, b0/a0)
 
from ASE import Atom
atom1 = Atom('C', r1)
atom2 = Atom('O', r2)
 
molecule0 = []
molecule0.append(atom1)
molecule0.append(atom2)
 
from ASE import ListOfAtoms
model0 = ListOfAtoms(molecule0)
 
model1.SetUnitCell(box1)
 
from Dacapo import Dacapo 
method1 = Dacapo()
method1.SetNumberOfBands(6) 
 
model1.SetCalculator(method1)
method1.Calculate()

from ASE import Atom, ListOfAtoms
model0 = ListOfAtoms([])
from Dacapo import Dacapo 
method0 = Dacapo()
 
model0.append(Atom('C', (0, 0, 0)))
model0.SetUnitCell((5, 5, 5))
model0.append(Atom('O', (0, 0, 1.1)))
 
method1.SetNumberOfBands(6) 
 
model1.SetCalculator(method1)
method1.Calculate()