18.11.2012 Views

ASE Manual Release 3.6.1.2825 CAMd - CampOS Wiki

ASE Manual Release 3.6.1.2825 CAMd - CampOS Wiki

ASE Manual Release 3.6.1.2825 CAMd - CampOS Wiki

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

# Optimize the Path:<br />

relax = BFGS(band, trajectory=’neb.traj’)<br />

relax.run(fmax=0.05)<br />

Example3: Diffusion run using NEB, Restart old calculation<br />

#!/usr/bin/env python<br />

import os<br />

from ase.io import read<br />

from ase.neb import NEB<br />

from ase.calculators.turbomole import Turbomole<br />

from ase.optimize import BFGS<br />

initial = read(’initial.coord’)<br />

final = read(’final.coord’)<br />

os.system(’rm -f coord; cp initial.coord coord’)<br />

#restart<br />

configs = read(’neb.traj@-5:’)<br />

band = NEB(configs, climb=True)<br />

#Set calculators<br />

for config in configs:<br />

config.set_calculator(Turbomole())<br />

# Optimize:<br />

relax = BFGS(band, trajectory=’neb.traj’)<br />

relax.run(fmax=0.05)<br />

For developers: python files affected by the turbomole interface<br />

ase:<br />

__init__.py<br />

new lines 17from<br />

ase.calculators import LennardJones, \<br />

EMT, ASAP, Siesta, Dacapo, Vasp, Turbomole<br />

ase/io:<br />

__init__.py<br />

new line 44<br />

TURBOMOLE coord file (filename===’coord’)<br />

new lines 145if<br />

format == ’turbomole’:<br />

from ase.io.turbomole import read_turbomole<br />

return read_turbomole(filename)<br />

new line 184<br />

TURBOMOLE coord file turbomole<br />

new lines 267elif<br />

format == ’tmol’:<br />

from ase.io.turbomole import write_turbomole<br />

write_turbomole(filename, images)<br />

return<br />

new lines 349if<br />

lines[0].startswith(’$coord’):<br />

return ’turbomole’<br />

<strong>ASE</strong> <strong>Manual</strong>, <strong>Release</strong> 3.6.1.2828<br />

7.14. Calculators 121

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!