That is ok, but it was ok with fix NVT and without sorting id, which is the only change I made. Now, when I calculate RDF it is showing overlapping . Is it related with the sorting of the coordinates . I have added sorting just according to the particle index like 1,2,3,… instead of random .