Building Crystals and Slabs¶

The Crystal Structures Tool¶

If you are lucky your crystal structure is already included in AMSinput. Most of the common structures are there. NaCl is one of them.

Start AMSinput (SCM → New Input) and switch from to (or to DFTB, and switch the periodicity to ‘Bulk’)

Click on the snowflake-like pictogram on the toolbar.

Select a Cubic lattice and then NaCl

Next a builder menu will appear on the right-hand side, where you can change the parameters of the structure, such as lattice constants

In this case there is no need to change anything.

You can adjust the details like lattice constants by hand. For example LiF: it has the same crystal structure as NaCl, but other elements and a different lattice constant, namely 4.02:

→ Cubic → NaCl

Change the lattice constant and the elements as shown below

Click OK

The NaCl crystal structure will be replaced by the new LiF crystal structure.

Some crystal structures have more parameters.

→ Tetragonal → Rutile

Click OK

As you can see, it has two lattice parameters that need to be specified. But now also the positions of the oxygens have a parameter “x”

After clicking OK your screen will look like this

The Crystal Structures Database¶

All crystal structures known to the crystal structures tool are also available via the search tool.

Just as you can search for a molecule, you can search for a crystal. So lets make NaCl again:

You should get a couple of matches:

Select the NaCl crystal result (not the molecule!):

You now should have a NaCl crystal again.

Currently the structures of many simple crystals are included, as well as a few complex structure (optimized by MOPAC). The zeolite frameworks are also included as a starting point to make your zeolite structure.

Crystal builder (from space group information)¶

The structure database is very convenient but by necessity incomplete. A more powerful approach is to build 3D crystals from the space group information. Here we give an example on how to build TiO2 (Rutile) again, now from its space group.

Open the crystal builder:

Click on → New From Space Group …

And the following menu will appear on the right-hand side:

Rutile has the symmetry of space group 136

Enter 136 in the space group Number field, and press return/enter.

Set the lattice parameters a to 4.59 and b to 2.96

Click on next to Coordinates, select Ti atom and a site

Click on again,. Select O atom and f site, and set x=0.3

The final should look like this:

In the toolbar there is another tool that may be of use: the star button. Just like for molecules this button will try to detect the symmetry of your system and make small changes to adjust to perfect symmetry. It does this via spglib, see the help balloon for details.

You already had perfect symmetry, so nothing much changes. The detected symmetry should be reported:

You can view a larger part of the crystal and still see the unit cell clearly:

This will give you a view like this:

Slicer: building slabs, transform primitive to conventional cell¶

The slicer is a very easy, yet powerful tool to make slabs from any crystal structure. In this example we build a three layer slab of the Cu(111) surface.

→ Cubic → fcc | Select fcc from the “Cubic” crystals

The element and lattice constant are already correct for Cu.

Let us invoke the slicer tool to cut out the slab.

Click on the slice tool, or select → Generate Slab.

The right side of the AMSinput window will show the Generate Slab options.

If you hold your mouse over the Miller indices fields, you should get a balloon that explains how to set the proper Miller indices. In this particular example we have the primitive unit cell of Cu, so first we need to convert to conventional unit cell.

You will see (from the top), if the periodic view tool is activated,

You may want to rotate it or select View → View Direction → Along X-axis to convince yourself that it is a three layer slab.

Creating a supercell¶

Quite often you want to use a larger unit cell. You can do this for chains, slabs, and crystals, but we will demonstrate how it works for a slab.

We will continue to work with the 3-layer Cu slab:

Thus invoking the Super Cell Tool

Here you see how new lattice vectors are expressed in terms of old ones. Because we have a slab this is a 2x2 matrix.

and you get a bigger unit cell with three atoms per layer: