atomsk --central-symmetry <file> [options]
This mode computes a central symmetry parameter for each atom.
The central symmetry parameter was first introduced in C.L. Kelchner, Phys. Rev. B 58 (1998) 11085, and later modified in its implementation in the visualization software Atomeye. This mode follows the latter, i.e. the central symmetry parameter is unitless and defined by:
ci = ( ∑ |dj+dk|2 ) / ( 2 ∑ |dj|2 )
where the atoms j and k are opposite neighbors of the central atom i, and the sum runs over all pairs of opposite neighbors.
This central symmetry parameters was designed for mono-atomic materials. In fcc metals each atom has 12 nearest neighbors; in bcc metals each atom has 8 nearest neighbors. In perfect lattices each neighbor has a counterpart, so that the central symmetry is perfect and for each atom ci=0. When the symmetry is broken, e.g. in surfaces where some neighbors are missing, in stacking faults or in dislocations where the environment is distorted, the atoms will have a parameter ci>0.
Atomsk can also compute this parameter in other types of materials, when they possess some kind of central symmetry. It is the case, for instance, of some binary materials with the rock-salt lattice (NaCl, MgO, LiF...), or ternary materials with the perovskite lattice (SrTiO3, BaTiO3...). In such lattices Atomsk computes the central symmetry parameter using only the neighbors that belong to the same species as the first neighbor. For instance in NaCl, for a central atom i of natrium (Na) only chloride (Cl) neighbors are used in the computation of ci, and vice-versa. In SrTiO3, the first neighbors are oxygen ions for Sr, oxygen ions for Ti, and oxygen ions for oxygen.
Note that this parameter cannot be expected to be meaningful in systems that do not have some sort of central symmetry, for instance in diamond, zincblende, or hexagonal lattices. In non-cubic perovskites, the lattice may not be centrosymmetric, making this parameter less relevant.
If this mode is used with one or several options they will be applied to the system before computing the central symmetry parameter.
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 --central-symmetry Al_dislocation.cfg Al_centrosymm.cfg -wrap
This will read the file Al_dislocation.cfg, and compute the central symmetry parameter for each atom. The final result will be output to Al_centrosymm.cfg.
atomsk --central-symmetry SrTiO3.xsf SrTiO3_centrosymm.cfg -wrap
This will read the file SrTiO3.xsf, and compute the central symmetry parameter for each atom. The final result will be output to SrTiO3_centrosymm.cfg.