====== Lesson 2 : 「変数」を使う ====== 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()