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

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

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

Any argument can be changed after the initial construction of the<br />

calculator, simply by setting it with the method<br />

>>> calc.set( keyword=value )<br />

Volumetric Data Output<br />

The class<br />

class FHI-aims.AimsCube(origin=(0,0,0),edges=[(0.1,0.0,0.0),(0.0,0.1,0.0),(0.0,0.0,0.1)],points=(50,50,50),plots=None)<br />

describes an object that takes care of the volumetric output requests within FHI-aims. An object of this type can<br />

be attached to the main Aims() object as an option.<br />

The possible arguments for AimsCube are:<br />

keyword type<br />

origin list<br />

edges 3x3-array<br />

points list<br />

plots list<br />

The possible values for the entry of plots are discussed in detail in the FHI-aims manual, see below for an example.<br />

Example<br />

As an example, here is a possible setup to obtain the geometry of a water molecule:<br />

from ase import Atoms<br />

from ase.visualize import view<br />

from ase.calculators.aims import Aims, AimsCube<br />

from ase.optimize import QuasiNewton<br />

water = Atoms(’HOH’, [(1,0,0), (0,0,0), (0,1,0)])<br />

water_cube = AimsCube(points=(29,29,29),<br />

plots=(’total_density’,’delta_density’,<br />

’eigenstate 5’,’eigenstate 6’))<br />

calc=Aims(xc=’pbe’,<br />

sc_accuracy_etot=1e-6,<br />

sc_accuracy_eev=1e-3,<br />

sc_accuracy_rho=1e-6,<br />

sc_accuracy_forces=1e-4,<br />

species_dir=’/home/hanke/codes/fhi-aims/fhi-aims.workshop/species_defaults/light/’,<br />

run_command=’aims.workshop.serial.x’,<br />

cubes=water_cube)<br />

water.set_calculator(calc)<br />

dynamics = QuasiNewton(water,trajectory=’square_water.traj’)<br />

dynamics.run(fmax=0.01)<br />

view(water)<br />

LAMMPS<br />

Introduction<br />

LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator) is a classical molecular dynamics code.<br />

126 Chapter 7. Documentation for modules in <strong>ASE</strong>

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

Saved successfully!

Ooh no, something went wrong!