Hi,

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

# Initial phase

timestep 0.00001

dump dump4 allatoms custom 100 10_100_2_init.voronoi type id c_compute3[1] c_compute3[2]

run 20000

undump dump4

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.

Thanks.

Tim