atomsk --nye <file1> <file2> [options] [<outputfile>] [formats]
This mode computes the per-atom Nye tensor, a quantity related to the strain gradient around an atom. The Nye tensor can be useful for visualizing dislocations. The calculation follows the method described in C.S. Hartley, Y. Mishin, Acta Mater. 53 (2005) 1313.
This mode requires two files: a <file1> containing the reference system, i.e. the system with no defect and all atoms in their perfect lattice sites; and a <file2> that can contain a defect (e.g. a dislocation). The two systems must have the exact same number of particles, and shall correspond to the same system (i.e. same base vectors, same atom types, etc.) for the calculation to make sense. The two files may be in different formats as long as these formats are supported by Atomsk.
The nine components of the Nye tensor (α11, α12, α13, α21, α22, α23, α31, α32, α33), are saved as auxiliary properties for each atom. Note that Atomsk will output these auxiliary properties only to the file formats that support it, like e.g. Atomeye CFG format (see this page for a list of formats supporting auxiliary properties).
If you use this mode with one or several options they will be applied to the <file1> and <file2> before the Nye tensor is computed.
Note that this mode assumes that the atom coordinates are wrapped, i.e. that all atoms are inside of the simulation box. If it is not the case then the calculation may be wrong. Coordinates can be wrapped thanks to the option -wrap
.
atomsk --nye reference.xsf dislocation.cfg nye.cfg
This will compute the Nye tensor due to the displacements in the file dislocation.cfg
, using reference.xsf
as a reference system. The positions and Nye tensors of atoms will be written to nye.cfg
.