I’m trying to extract velocity around each atom in my crystal with ppp boundary condition. But when I print the velocity of the neighbors (v=atom->v, then get v[j], v[j], v[j] for each j in the neighbor list of i). I got zero for some atoms, which does not make sense to me. I’m guessing they are near-boundary atoms which have neighbors wrapped due to the periodic boundary condition. PS: The positions are all making sense.
Hope my description is clear. Appreciate any guidance in extracting such information correctly!