Identifying self-interstitials of bcc and fcc crystals in molecular dynamics

Identification of self-interstitials in molecular dynamics (MD) simulations is of critical importance. There exist several criteria for identifying the self-interstitial. Most of the existing methods use an assumed cut-off value for the displacement of an atom from its lattice position to identify the self-interstitial. The results obtained are affected by the chosen cut-off value. Moreover, these chosen cut-off values are independent of temperature. We have developed a novel unsupervised learning algorithm called Max-Space Clustering (MSC) to identify an appropriate cut-off value and its dependence on temperature. This method is compared with some widely used methods such as effective sphere (ES) method and nearest neighbor sphere (NNS) method. The cut-off radius obtained using our method shows a linear variation with temperature. The value of cut-off radius and its temperature dependence is derived for five bcc (Cr, Fe, Mo, Nb, W) and six fcc (Ag, Au, Cu, Ni, Pd, Pt) crystals. It is seen that the ratio of the cut-off values “r” to the lattice constant “a” lies between 0.23 and 0.3 at 300 K and this ratio is on an average smaller for the fcc crystals. Collision cascade simulations are carried out for Primary knock-on Atom (PKA) energies of 5 keV in Fe (at 300 K and 1000 K) and W (at 300 K and 2500 K) and the results are compared using the various methods.