MUsmanJaved007 From your input file, it seems that you have four species in your system Ti, Cr, C, O
, however, there are 7 atoms that form your structure (look at atom.xyz
). In this case atom.element
should point to type of species for each atom in atom.xyz
. Let's assume the 7 atoms in your system in order of their coordinates in atom.xyz
are Tr, O, O, C, Cr, O, Tr
, then atom.element = [1, 4, 4, 3, 2, 4, 1]
. Please correct you input file accordingly.