This modifier calculates the centrosymmetry parameter (CSP) [Kelchner, Plimpton, Hamilton, Phys. Rev. B, 58, 11085 (1998)] for each particle. In solid-state systems the centrosymmetry parameter is a useful measure of the local lattice disorder around an atom and can be used to characterize whether the atom is part of a perfect lattice, a local defect (e.g. a dislocation or stacking fault), or located at a surface.
The CSP value pCSP of an atom having N nearest neighbors (N = 12 for face-centered cubic, N = 8 for body-centered cubic lattices) is given by
where ri and ri+N/2 are vectors pointing from the central atom to a pair of opposite neighbors. For lattice sites in an ideal centrosymmetric crystal, the contributions of all neighbor pairs in this formula will cancel, and the resulting CSP value will hence be zero. Atomic sites within a defective crystal region, in contrast, typically have a disturbed, non-centrosymmetric neighborhood. In this case the CSP becomes positive. Using an appropriate threshold, to allow for small perturbations due to thermal displacements and elastic strains, the CSP can be used as an order parameter to filter out atoms that are part of crystal defects.
The calculated atomic CSP values are stored in the
Centrosymmetry output particle property by the modifier.
A histogram of the CSP values of the entire particle system is displayed in the modifier panel.
Number of neighbors
This parameter specifies the number of nearest neighbors that should be taken into account when computing the centrosymmetry value for an atom. This parameter value should match the ideal number of nearest neighbors in the crystal lattice at hand (12 in fcc crystals; 8 in bcc). More generally, it must be a positive, even integer.
The modifier supports two modes of operation:
- Conventional CSP
This mode implements the same algorithm as LAMMPS to calculate the centrosymmetry parameter. Weights are calculated between all N (N - 1) / 2 pairs of neighbor atoms, and the CSP is calculated as the summation of the N / 2 smallest weights.
- Minimum-weight matching CSP
The conventional CSP algorithm performs well on highly centrosymmetric structures. In acentrosymmetric structures, however, it often assigns similar CSP values to very different structures, which results from the "greedy" selection of neighbor pair weights. The Minimum-weight matching CSP [Larsen] ensures that neighbor relationships are reciprocal, which results in a better separation of CSP values between e.g. HCP atoms and surface defect atoms.
This algorithm has a significantly higher computational requirements.
Note that the modifier needs to see the complete set of particles to perform the computation. It should therefore be placed at the beginning of the data pipeline, preceding any modifiers that delete particles.
You can use the Color coding modifier to color atoms based on their CSP value or use the Expression selection modifier to select atoms having a CSP value below some threshold. These undisturbed atoms can then be hidden to reveal crystal defect atoms by using the Delete selected modifier.