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.

seq: A sequence of three int‘s.<br />

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

For parameters in the list without default value given, VASP will set the default value. Most of the parameters<br />

used in the VASP INCAR file are allowed keywords. See the official VASP manual for more details.<br />

Note: Parameters can be changed after the calculator has been constructed by using the set() method:<br />

>>> calc.set(prec=’Accurate’, ediff=1E-5)<br />

This would set the precision to Accurate and the break condition for the electronic SC-loop to 1E-5 eV.<br />

Spin-polarized calculation<br />

If the atoms object has non-zero magnetic moments, a spin-polarized calculation will be performed by default.<br />

Here follows an example how to calculate the total magnetic moment of a sodium chloride molecule.<br />

from ase import Atoms, Atom<br />

from ase.calculators.vasp import Vasp<br />

a = [6.5, 6.5, 7.7]<br />

d = 2.3608<br />

NaCl = Atoms([Atom(’Na’, [0, 0, 0], magmom=1.928),<br />

Atom(’Cl’, [0, 0, d], magmom=0.75)],<br />

cell=a)<br />

calc = Vasp(prec = ’Accurate’,<br />

xc = ’PBE’,<br />

lreal = False)<br />

NaCl.set_calculator(calc)<br />

print NaCl.get_magnetic_moment()<br />

In this example the initial magnetic moments are assigned to the atoms when defining the Atoms object. The calculator<br />

will detect that at least one of the atoms has a non-zero magnetic moment and a spin-polarized calculation<br />

will automatically be performed. The <strong>ASE</strong> generated INCAR file will look like:<br />

INCAR created by Atomic Simulation Environment<br />

PREC = Accurate<br />

LREAL = .FALSE.<br />

ISPIN = 2<br />

MAGMOM = 1*1.9280 1*0.7500<br />

Note: It is also possible to manually tell the calculator to perform a spin-polarized calculation:<br />

>>> calc.set(ispin=2)<br />

This can be useful for continuation jobs, where the initial magnetic moment is read from the WAVECAR file.<br />

Restart old calculation<br />

To continue an old calculation which has been performed without the interface use the restart parameter when<br />

constructing the calculator<br />

>>> calc = Vasp(restart=True)<br />

Then the calculator will read atomic positions from the CONTCAR file, physical quantities from the OUTCAR file,<br />

k-points from the KPOINTS file and parameters from the INCAR file.<br />

7.14. Calculators 123

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

Saved successfully!

Ooh no, something went wrong!