Surfaces

Common surfaces

A number of utility functions are provided to set up the most common surfaces, to add vacuum layers, and to add adsorbates to a surface. In general, all surfaces can be set up with the modules described in the section General crystal structures and surfaces, but these utility functions make common tasks easier.

Example

To setup an Al(111) surface with a hydrogen atom adsorbed in an on-top position:

from ase.lattice.surface import fcc111
slab = fcc111('Al', size=(2,2,3), vacuum=10.0)

This will produce a slab 2x2x3 times the minimal possible size, with a (111) surface in the z direction. A 10 Å vacuum layer is added on each side.

To set up the same surface with with a hydrogen atom adsorbed in an on-top position 1.5 Å above the top layer:

from ase.lattice.surface import fcc111, add_adsorbate
slab = fcc111('Al', size=(2,2,3))
add_adsorbate(slab, 'H', 1.5, 'ontop')
slab.center(vacuum=10.0, axis=2)

Note that in this case is is probably not meaningful to use the vacuum keyword to fcc111, as we want to leave 10 Å of vacuum after the adsorbate has been added. Instead, the center() method of the Atoms is used to add the vacuum and center the system.

The atoms in the slab will have tags set to the layer number: First layer atoms will have tag=1, second layer atoms will have tag=2, and so on. Adsorbates get tag=0:

>>> print atoms.get_tags()
[3 3 3 3 2 2 2 2 1 1 1 1 0]

This can be useful for setting up ase.constraints (see Diffusion of gold atom on Al(100) surface (NEB)).

Utility functions for setting up surfaces

All the functions setting up surfaces take the same arguments.

symbol:
The chemical symbol of the element to use.
size:
A tuple giving the system size in units of the minimal unit cell.
a:
(optional) The lattice constant. If specified, it overrides the expermental lattice constant of the element. Must be specified if setting up a crystal structure different from the one found in nature.
c:
(optional) Extra HCP lattice constant. If specified, it overrides the expermental lattice constant of the element. Can be specified if setting up a crystal structure different from the one found in nature and an ideal \(c/a\) ratio is not wanted (\(c/a=(8/3)^{1/2}\)).
vacuum:
The thickness of the vacuum layer. The specified amount of vacuum appears on both sides of the slab. Default value is None, meaning not to add any vacuum. In that case a “vacuum” layer equal to the interlayer spacing will be present on the upper surface of the slab. Specify vacuum=0.0 to remove it.
orthogonal:
(optional, not supported by all functions). If specified and true, forces the creation of a unit cell with orthogonal basis vectors. If the default is such a unit cell, this argument is not supported.

Each function defines a number of standard adsorption sites that can later be used when adding an adsorbate with ase.lattice.surface.add_adsorbate().

The following functions are provided

ase.lattice.surface.fcc100(symbol, size, a=None, vacuum=0.0)[source]
ase.lattice.surface.fcc110(symbol, size, a=None, vacuum=0.0)[source]
ase.lattice.surface.bcc100(symbol, size, a=None, vacuum=0.0)[source]
ase.lattice.surface.hcp10m10(symbol, size, a=None, c=None, vacuum=0.0)[source]
ase.lattice.surface.diamond100(symbol, size, a=None, vacuum=0.0)[source]

These allways give orthorhombic cells:

fcc100 fcc100
fcc110 fcc110
bcc100 bcc100
hcp10m10 hcp10m10
diamond100 diamond100
ase.lattice.surface.fcc111(symbol, size, a=None, vacuum=0.0, orthogonal=False)[source]
ase.lattice.surface.fcc111_root(symbol, root, size, a=None, vacuum=0.0, orthogonal=False, search_space=(20, 20))[source]
ase.lattice.surface.fcc211(symbol, size, a=None, vacuum=0.0, orthogonal=True)[source]
ase.lattice.surface.bcc110(symbol, size, a=None, vacuum=0.0, orthogonal=False)[source]
ase.lattice.surface.bcc111(symbol, size, a=None, vacuum=0.0, orthogonal=False)[source]
ase.lattice.surface.hcp0001(symbol, size, a=None, c=None, vacuum=0.0, orthogonal=False)[source]
ase.lattice.surface.diamond111(symbol, size, a=None, vacuum=0.0, orthogonal=False)[source]

These can give both non-orthorhombic and orthorhombic cells:

fcc111 fcc111 fcc111o
fcc111_root varies on root not_implemented
fcc211 not implemented fcc211o
bcc110 bcc110 bcc110o
bcc111 bcc111 bcc111o
hcp0001 hcp0001 hcp0001o
diamond111 diamond111 not implemented

The adsorption sites are marked with:

ontop hollow fcc hcp bridge shortbridge longbridge
ontop hollow fcc hcp bridge bridge bridge

This can be used for MX2 2D structures such as MoS2:

ase.lattice.surface.mx2(formula='MoS2', kind='2H', a=3.18, thickness=3.19, size=(1, 1, 1), vacuum=7.5)[source]

Create three-layer 2D materials with hexagonal structure.

For metal dichalcogenites, ect.

The kind argument accepts ‘2H’, which gives a mirror plane symmetry and ‘1T’, which gives an inversion symmetry.

../_images/mx2.png

Adding new utility functions

If you need other surfaces than the ones above, the easiest is to look in the source file surface.py, and adapt one of the existing functions. Please contribute any such function that you make either by checking it into SVN or by mailing it to the developers.

Adding adsorbates

After a slab has been created, a vacuum layer can be added. It is also possible to add one or more adsorbates.

ase.lattice.surface.add_adsorbate(slab, adsorbate, height, position=(0, 0), offset=None, mol_index=0)[source]

Add an adsorbate to a surface.

This function adds an adsorbate to a slab. If the slab is produced by one of the utility functions in ase.lattice.surface, it is possible to specify the position of the adsorbate by a keyword (the supported keywords depend on which function was used to create the slab).

If the adsorbate is a molecule, the atom indexed by the mol_index optional argument is positioned on top of the adsorption position on the surface, and it is the responsibility of the user to orient the adsorbate in a sensible way.

This function can be called multiple times to add more than one adsorbate.

Parameters:

slab: The surface onto which the adsorbate should be added.

adsorbate: The adsorbate. Must be one of the following three types:
A string containing the chemical symbol for a single atom. An atom object. An atoms object (for a molecular adsorbate).

height: Height above the surface.

position: The x-y position of the adsorbate, either as a tuple of
two numbers or as a keyword (if the surface is produced by one of the functions in ase.lattice.surfaces).
offset (default: None): Offsets the adsorbate by a number of unit
cells. Mostly useful when adding more than one adsorbate.
mol_index (default: 0): If the adsorbate is a molecule, index of
the atom to be positioned above the location specified by the position argument.

Note position is given in absolute xy coordinates (or as a keyword), whereas offset is specified in unit cells. This can be used to give the positions in units of the unit cell by using offset instead.

ase.lattice.surface.add_vacuum(atoms, vacuum)[source]

Add vacuum layer to the atoms.

Parameters:

atoms: An Atoms object most likely created by one of the ase.lattice modules.

vacuum: The thickness of the vacuum layer (in Angstrom).

Create specific non-common surfaces

New in version 3.5.2.

In addition to the most normal surfaces, a function has been constructed to create more uncommon surfaces that one could be interested in. It is constructed upon the Miller Indices defining the surface and can be used for both fcc, bcc and hcp structures. The theory behind the implementation can be found here: general_surface.pdf.

ase.lattice.surface.surface(lattice, indices, layers, vacuum=0.0, tol=1e-10)

Create surface from a given lattice and Miller indices.

lattice: Atoms object or str
Bulk lattice structure of alloy or pure metal. One can also give the chemical symbol as a string, in which case the correct bulk lattice will be generated automatically.
indices: sequence of three int
Surface normal in Miller indices (h,k,l).
layers: int
Number of equivalent layers of the slab.
vacuum: float
Amount of vacuum added on both sides of the slab.

Example

To setup a Au(211) surface with 9 layers and 10 Å of vacuum:

from ase.lattice.surface import surface
s1 = surface('Au', (2, 1, 1), 9)
s1.center(vacuum=10, axis=2)

This is the easy way, where you use the experimental lattice constant for gold bulk structure. You can write:

from ase.visualize import view
view(s1)

or simply s1.edit() if you want to see and rotate the structure.

../_images/s1.png

Next example is a molybdenum bcc(321) surface where we decide what lattice constant to use:

from ase.lattice import bulk
Mobulk = bulk('Mo', 'bcc', a=3.16, cubic=True)
s2 = surface(Mobulk, (3, 2, 1), 9)
s2.center(vacuum=10, axis=2)
../_images/s2.png

As the last example, creation of alloy surfaces is also very easily carried out with this module. In this example, two Pt3Rh fcc(211) surfaces will be created:

a = 4.0
from ase import Atoms
Pt3Rh = Atoms('Pt3Rh',
              scaled_positions=[(0, 0, 0),
                                (0.5, 0.5, 0),
                                (0.5, 0, 0.5),
                                (0, 0.5, 0.5)],
              cell=[a, a, a],
              pbc=True)
s3 = surface(Pt3Rh, (2, 1, 1), 9)
s3.center(vacuum=10, axis=2)

Pt3Rh.set_chemical_symbols('PtRhPt2')
s4 = surface(Pt3Rh , (2, 1, 1), 9)
s4.center(vacuum=10, axis=2)

s3 s4