Skip to content

Common Neighbour Analysis

The common neighbour analysis is a signature attributed to all pairs of nearest neighbour atoms in the system used to classify the geometrical environment of the atoms. For each pair of neighbours we evaluate a triplet of integers: (r, s, t) as:

Number of shared neighbours of the pair, here the two green atoms share 4 neighbours (coloured in blue) so in this example the pair will have r = 4 Image title

Number of bonds between the r nieghbours (excluding pair), here you can see that there are two bonds between the 4 neighbours of the pair taken into consideration, so s = 1 Image title

Longest path along connected bonds between the r neighbours (excluding pair), here you can see that the two neighbours only form a path of length 1, so t = 1 Image title

Computing the pair CNA Signature

As said, the CNA signature is a pair property, in pySNOW we can compute the CNA for each pair of atom and obtain it as an \(3\cdot N_{pairs}\) array where for each pair the subarray contains the (r, s, t) values, optionally the function will also return a list containing the indeces of the atoms forming each pair for which the CNA signature has been computed, this can be used (and is used by other fnctions which we will present below) to identify patterns in the observed CNA signatures.

The method to compute the CNA involves, as always, as a first step, the reading of the coordinates of the atoms, then the function calculate_cna(...) can be used to obtain the desired output:

1
2
3
4
5
6
from snow.lodispp.pp_io import read_xyz
from snow.lodispp.cna import calculate_cna

el, coords = read_xyz("Au976To.xyz")
cut_off = 4.079 * 0.85 # (1)
n_pairs, cnas, pairs = calculate_cna(1, coords, cut_off, return_pair = True) # (2)
  1. Cutoff as 0.85*lattice parameter
  2. The first 1 is the index_frame placeholder. Setting return_pair to True will return the indeces of the atoms forming each pair
1
2
3
4
5
6
7
8
9
from ase.io import read
from snow.lodispp.cna import calculate_cna

nano = read("Au976To.xyz")

coords = nano.get_positions()

cut_off = 4.079 * 0.85 # (1)
n_pairs, cnas, pairs = calculate_cna(1, coords, cut_off, return_pair = True) # (2)
  1. Cutoff as 0.85*lattice parameter
  2. The first 1 is the index_frame placeholder. Setting return_pair to True will return the indeces of the atoms forming each pair
 1
 2
 3
 4
 5
 6
 7
 8
 9
10
from from ovito.io import import_file

from snow.lodispp.cna import calculate_cna

pipeline = import_file("Au976To.xyz")

coords = nano.get_positions()

cut_off = 4.079 * 0.85 # (1)
n_pairs, cnas, pairs = calculate_cna(1, coords, cut_off, return_pair = True) # (2)
  1. Cutoff as 0.85*lattice parameter
  2. The first 1 is the index_frame placeholder. Setting return_pair to True will return the indeces of the atoms forming each pair

CNAPs are described in the following table, also available in the readme of the library:

CNAp Index Descriptions

CNAp Description CNAp composition
1 Vertex between two (111) facets and a (100) facet [(1, (100)), (2, (211)), (1, (322)), (1, (422))]
2 Edge between (100) and a slightly distorted (111) [(1, (200)), (2, (211)), (2, (311)), (1, (421))]
3 Atoms lying on a (555) symmetry axis [(10, (422)), (2, (555))]
4 FCC bulk [(12, (421))]
5 Intersection of six five-fold axes [(12, (555))]
6 Edge between (100) facets [(2, (100)), (2, (211)), (2, (422))]
7 Vertex on twinning planes shared by (111) facets [(2, (200)), (1, (300)), (2, (311)), (1, (322)), (1, (422))]
8 Edge between (111) re-entrances and (111) facets [(2, (200)), (4, (311)), (1, (421))]
9 Re-entrance delimited by (111) facets [(2, (300)), (4, (311)), (2, (421)), (2, (422))]
10 Edge between (100) and (111) facets [(3, (211)), (2, (311)), (2, (421))]
11 Vertex shared by (100) and (111) facets [(4, (211)), (1, (421))]
12 (100) facet [(4, (211)), (4, (421))]
13 Five-fold symmetry axis (without center) [(4, (311)), (2, (322)), (2, (422))]
14 Five-fold vertex [(5, (322)), (1, (555))]
15 (111) face [(6, (311)), (3, (421))]
16 Twinning plane [(6, (421)), (6, (422))]

Ovito users can also make use of our Ovito modifier CNAPatterner which provides a similar kind of analysis. While Ovito pipelines are generally tricky to deal with in python, using it is advisable when treating with larger systems and/or systems with periodic boundary conditions.