Dear Lammps users,

I am tracking density/mass change in solid melting process. I did it using “compute chunk/atom” and “fix ave/chunk” but here is one issue I do not understand:

I got correct density when I used compute chunk/atom bin/3d and defined a box for calculation of density.

When I used “compute chunk/atom type”, the results got to 211054 which I do not know where it comes from (it is not atom number).

In the fix ave/chunk documentation, it says when a box is defined, use volume of box to calculate results, if not, used whole simulation box volume. This does not match my results either.

Anyone got any clue about how this happens?

Thanks for help.