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.

examples:<br />

m: int The length of the nanoribbon.<br />

type: str The orientation of the ribbon. Must be either ‘zigzag’ or ‘armchair’.<br />

saturated: bool If true, hydrogen atoms are placed along the edge.<br />

C_H: float Carbon-hydrogen bond length. Default: 1.09 Angstrom.<br />

C_C: float Carbon-carbon bond length. Default: 1.42 Angstrom.<br />

vacuum: float Amount of vacuum added to both sides. Default 2.5 Angstrom.<br />

magnetic: bool Make the edges magnetic.<br />

initial_mag: float Magnitude of magnetic moment if magnetic=True.<br />

sheet: bool If true, make an infinite sheet instead of a ribbon.<br />

>>> from ase.structure import graphene_nanoribbon<br />

>>> gnr1 = graphene_nanoribbon(3, 4, type=’armchair’, saturated=True)<br />

>>> gnr2 = graphene_nanoribbon(2, 6, type=’zigzag’, saturated=True,<br />

>>> C_H=1.1, C_C=1.4, vacuum=6.0,<br />

>>> magnetic=True, initial_mag=1.12)<br />

7.9.4 Special points in the Brillouin zone<br />

You can find the psecial points in the Brillouin zone:<br />

>>> from ase.structure import bulk<br />

>>> from ase.dft.kpoints import ibz_points, get_bandpath<br />

>>> si = bulk(’Si’, ’diamond’, a=5.459)<br />

>>> points = ibz_points[’fcc’]<br />

>>> G = points[’Gamma’]<br />

>>> X = points[’X’]<br />

>>> W = points[’W’]<br />

>>> K = points[’K’]<br />

>>> L = points[’L’]<br />

>>> kpts, x, X = get_bandpath([W, L, G, X, W, K], si.cell)<br />

>>> print len(kpts), len(x), len(X)<br />

50 50 6<br />

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

7.9. Creating atomic structures 85

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

Saved successfully!

Ooh no, something went wrong!