/doc/ase/calculators/dftb_ex1_relax.py

 1from ase.calculators.dftb import Dftb
 2from ase.optimize import QuasiNewton
 3from ase.io import write
 4
 5from ase.structure import molecule
 6test = molecule('H2O')
 7test.set_calculator(Dftb(label='h2o',
 8                         atoms=test,
 9                         Hamiltonian_MaxAngularMomentum_='',
10                         Hamiltonian_MaxAngularMomentum_O='"p"',
11                         Hamiltonian_MaxAngularMomentum_H='"s"',
12                         ))
13
14dyn = QuasiNewton(test, trajectory='test.traj')
15dyn.run(fmax=0.01)
16write('test.final.xyz', test)
17
18
19