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.

data.cpk_colors<br />

data.reference_states<br />

data.vdw_radii<br />

All of these are lists that should be indexed with an atomic number:<br />

>>> from ase.data import *<br />

>>> atomic_names[92]<br />

’Uranium’<br />

>>> atomic_masses[2]<br />

4.0026000000000002<br />

data.atomic_numbers<br />

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

If you don’t know the atomic number of some element, then you can look it up in the atomic_numbers<br />

dictionary:<br />

>>> atomic_numbers[’Cu’]<br />

29<br />

>>> covalent_radii[29]<br />

1.1699999999999999<br />

The covalent radii are taken from [Cordeo08]. The source of the van der Waals radii is given in vdw.py.<br />

7.24.2 Molecular data<br />

The G1, G2, and G3-databases are available in the molecules module.<br />

Example:<br />

>>> from ase.data.molecules import molecule<br />

>>> atoms = molecule(’H2O’)<br />

All molecular members of each database is conveniently contained in a list of strings (g1, g2, g3), and one can<br />

look up the experimental atomization energy for each molecule. This is extrapolated from experimental heats of<br />

formation at room temperature, using calculated zero-point energies and thermal corrections.<br />

Example:<br />

>>> from ase.data.molecules import molecule, g2, get_atomization_energy<br />

>>> g2[11]<br />

’H2O’<br />

>>> atoms = molecule(g2[11])<br />

>>> get_atomization_energy(g2[11])<br />

232.57990000000001<br />

>>> from ase.units import kcal,mol<br />

>>> get_atomization_energy(g2[11])*kcal/mol<br />

10.08562144637833<br />

where the last line converts the experimental atomization energy of H2O from units of kcal/mol to eV.<br />

7.24.3 S22, s26, and s22x5 data<br />

The s22, s26, and s22x5 databases are available in the s22 module.<br />

Each weakly bonded complex is identified as an entry in a list of strings (s22, s26, s22x5), and is fully created by<br />

a ‘create’-function:<br />

>>> from ase.data.s22 import s22, create_s22_system<br />

>>> sys = s22[0]<br />

>>> sys<br />

’Ammonia_dimer’<br />

7.24. The data module 167

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

Saved successfully!

Ooh no, something went wrong!