# UnitCell¶

class UnitCell

An UnitCell represent the box containing the atoms, and its periodicity

A unit cell is represented by the cell matrix, containing the three cell vectors:

```| a_x   b_x   c_x |
| a_y   b_y   c_y |
| a_z   b_z   c_z |
```

Public Types

enum CellShape

Possible shapes for the unit cell.

Values:

enumerator ORTHORHOMBIC

Orthorhombic cell, with the three angles equals to 90°

enumerator TRICLINIC

Triclinic cell, with any values for the angles.

enumerator INFINITE

Infinite cell, to use when there is no cell.

Public Functions

UnitCell()

Construct an `INFINITE` unit cell, with all lengths set to 0

UnitCell(Vector3D lengths)

Construct an `ORTHORHOMBIC` unit cell with the given cell `lengths`

UnitCell(Vector3D lengths, Vector3D angles)

Construct a unit cell with the given cell `lengths` and `angles`

If all lengths are set to 0, then the cell is `INFINITE`. If at least one length is not zero and all angles are 90.0, then the cell is `ORTHORHOMBIC`. Else a `TRICLINIC` cell is created.

UnitCell(Matrix3D matrix)

Construct a unit cell from a cell matrix.

If `matrix` contains only zeros, then an infinite cell is created. If only the diagonal of the matrix is non-zero, then the cell is `ORTHORHOMBIC`. Else a `TRICLINIC` cell is created. The matrix entries should be in Angstroms.

inline Matrix3D matrix() const

Get the cell matrix

inline CellShape shape() const

Get the cell shape

void set_shape(CellShape shape)

Set the cell shape to `shape`

Throws:

Error – if `shape` is `ORTHORHOMBIC` and some angles are not 90°, or if `shape` is `INFINITE` and some lengths are not 0.0.

Vector3D lengths() const

Get the lengths of the cell’s vectors, in angstroms.

Vector3D angles() const

Get the angles between the cell’s vectors in degrees

void set_lengths(Vector3D lengths)

Set the lengths of the cell’s vectors. The values should be in angstroms.

This function reset cell orientation!

After the call, the cell is aligned such that the first cell vector is along the x axis, and the second cell vector is in the xy plane.

Throws:

Error – if the cell shape is `INFINITE`.

void set_angles(Vector3D angles)

Set the angles between the cell’s vectors. The values should be in degrees.

This function reset cell orientation!

After the call, the cell is aligned such that the first cell vector is along the x axis, and the second cell vector is in the xy plane.

Throws:

Error – if the cell shape is not `TRICLINIC`.

double volume() const

Get the unit cell volume

Vector3D wrap(const Vector3D &vector) const

Wrap the `vector` in the unit cell, using periodic boundary conditions.

For an orthorhombic unit cell, this make sure that all the vector components are between `-L/2` and `L/2` where `L` is the corresponding cell length.