Hi Omakar, LAMMPS mailing list is the best place to ask such questions. I think I saw your post. Just checked the issue you have raised here. It is clearly mentioned in lammps manual. Use following commands and you will get exactly same result you r expecting.

compute peratom all pe/atom

compute pe all reduce sum c_peratom

thermo_style custom step temp etotal press pe c_pe

Hope it helps. 

Thanks you Dr Tiwari for your kind suggestions. But as you have already understand my problem so I am asking you again because if I go for LAMMPS mailing list I have to again explain everything. If it is possible from your side then please answer. From next time I will opt for LAMMPS mailing list.

Here I am writing the out put from the commands

fix 1 all nve

fix 2 all temp/rescale 20 2000.0 2000.0 10.0 1.0

timestep 1.0

thermo 1

run 1000

unfix 2

compute 3 all pe/atom

compute 4 all reduce sum c_3

thermo_style custom step temp pe etotal c_4

dump 20 all custom 10 Trajectory_nve_2000 id xu yu zu vx vy vz c_3

run 1000

OUTPUT of thermo_style

Step Temp PotEng TotEng 4

1000 1996.2558 -1459674.8 -1445209.2 -748980.04

1001 1997.0884 -1459680.7 -1445209.1 -748987.39

Now just see the 3rd and the last column of the output. The 3rd column is the PE of entire system and the last column is the sum of the PE of individual atoms. These two quantity should be equal but they differs by a multiplicative factor of approximately two. I don't understand why ?

Omkar

The whole input file.

units real

atom_style full

read_data inp_nve.xyz

kspace_style ewald/n 0.00001

kspace_modify gewald 0.24287

pair_style buck/coul/long 11.0 11.0

pair_coeff * * 0.0 1.0 0.0

pair_coeff  1  4 

pair_coeff 2   4

pair_coeff 3   4

pair_coeff 4   4

fix 1 all nve

fix 2 all temp/rescale 20 2000.0 2000.0 10.0 1.0

timestep 1.0

thermo 1

run 1000

unfix 2

compute 3 all pe/atom

compute 4 all reduce sum c_3

thermo_style custom step temp pe etotal c_4

dump 20 all custom 10 Trajectory_nve_2000 id xu yu zu vx vy vz c_3

run 1000

Hello Sir

Can you give me the input file that you have used ?