I’m studying damage cascades in FCC Ni-Fe-Cr. I’m trying to calculate vacancy and interstitial distributions using the “compute voronoi/atom occupation” command. Below is part of my code (showing only commands related to voronoi calculation)
After creating perfect FCC lattice
compute compute3 all voronoi/atom occupation
Initialize velocity at 800 C
Equilibrate for 2 ps
Set velocity PKA in the center = 10 keV
dump dump4 allatoms custom 100 10_100_2_init.voronoi type id c_compute3 c_compute3
I was able to successfully visualize vacancy and interstitial cascades using OVITO’s Wigner-Seitz defect analysis tool. But my output file for the voronoi data contains two columns of “1” only, which means that all the atoms are in the original voronoi cells as defined with the perfect lattice.
Why am I seeing only columns of 1’s? I was expecting 0’s for vacancies and >1 for interstitials.