Automatic k-mesh generation

The second method generates k-meshes automatically, and requires only the input of subdivisions of the Brillouin zone in each direction and the origin ('shift') for the k-mesh. There are three possible input formats. The simplest one is only supported by VASP.4.5 and newer versions:

Automatic mesh 0 ! number of k-points = 0 ->automatic generation scheme Auto ! fully automatic 10 ! length (l)As before, the first line is treated as a comment. On the second line a number smaller or equal 0 must be specified. In the previous section, this value supplied the number of k-points, a zero in this line activates the automatic generation scheme. The fully automatic scheme, selected by the first character in the third line ('a'), generates centered Monkhorst-Pack grids, where the numbers of subdivisions along each reciprocal lattice vector are given by

A slightly enhanced version, allows to supply the numbers for the subdivisions , and manually:

Automatic mesh 0 ! number of k-points = 0 ->automatic generation scheme Gamma ! generate a Gamma centered grid 4 4 4 ! subdivisions N_1, N_2 and N_3 along recipr. l. vectors 0. 0. 0. ! optional shift of the mesh (s_1, s_2, s_3)In this case, the third line (again, only the first character is significant) might start with 'G' or 'g' --for generating meshes with their origin at the point (as above)-- or 'M' or 'm', which selects the original Monkhorst-Pack scheme. In the latter case k-point grids, with

If the third line does not start with 'M', 'm', 'G' or 'g' an alternative input mode is selected. this mode is mainly for experts, and should not be used for casual VASP users. Here one can provide directly the generating basis vectors for the k-point mesh (in cartesian or reciprocal coordinates). The input file has the following format:

Automatic generation 0 Cartesian 0.25 0.00 0.00 0.00 0.25 0.00 0.00 0.00 0.25 0.00 0.00 0.00The entry in the third line switches between cartesian and reciprocal coordinates (as in the explicit input format described first - key characters 'C', 'c', 'K' or 'k' switch to cartesian coordinates). On the fifth, sixth and seventh line the generating basis vectors must be given and the eighth line supplies the shift (if one likes to shift the k-mesh off , default is to take the origin at , the shift is given in multiples of the generating basis vectors, only (0,0,0) and (1/2,1/2,1/2) and arbitrary combinations are usually usefull). This method can always be replaced by an appropriate Monkhorst-Pack setting. For instance for a fcc lattice the setting

cart 0.25 0 0 0 0.25 0 0 0 0.25 0.5 0.5 0.5is equivalent to

Monkhorst-pack 4 4 4 0 0 0This input scheme is especially interesting to build meshes, which are based on the conventional cell (for instance sc for fcc and bcc), or the primitive cell if a large super cell is used. In the example above the k-point mesh is based on the sc-cell. (for the second input file the tetrahedron method can not be used because the shift breaks the symmetry (see below), whereas the first input file can be used together with the tetrahedron method). For more hints please read section 8.6.

*Mind:*
The division scheme (or the generating basis of the k-mesh) must lead to a k-mesh
which belongs to the same class of Bravais lattice as the
reciprocal unit cell (Brillouin zone).
Any symmetry-breaking set-up of the mesh cannot be handled by VASP.
Hence such set-ups are not allowed --
if you break this rule an error message will be displayed.
Furthermore the symmetrisation of the k-mesh can lead to meshes which can not be
divided into tetrahedrons (at least not by the tetrahedron division scheme
implemented in VASP) --
if one uses meshes which do not have their origin at
(for certain lower symmetric types of Bravais lattices or certain
non-symmetry conserving shifts).
Therefore only very special shifts are allowed.
If a shift is selected which can not be handled you get an error message.
For reasons
of safety it might be a good choice to use only meshes
with their origin at
(switch 'G' or 'g' on third line or odd divisions)
if the tetrahedron method is used.

N.B. Requests for support are to be addressed to: vasp.materialphysik@univie.ac.at