I have a question about granular package in LAMMPS.
I wrote a code to determine the number of contacts of
each grain (in granular mode), by finding the distance
between grains from their xyz coordination at the
LAST time step. But the total number of contacts and
the mean values are very different from the values
appeared at LAMMPS output files. I wounder wether any
criterion other than the distance has been used for
finding neighbors in LAMMPS.
What values in the output are you referring to? LAMMPS
would always compute contact from the current timestep,
and it would mean 2 particles are overlapping, i.e. there
is a spring force and/or tangential force between them.