[lammps-users] about NVE

hi all users,

in my simulation, I am using NVE ensemble. Total
energy should be constant but, in my simulation, it is
not constant. What can create this problem? or should
I make change in neighbor and neighbor_modify command?

Thanks.

Dilek Cakiroglu
Mechanical and Aerospace Engineering
Oklahoma State University

Many things could cause this. A common one is using
a potential with a cutoff and not shifting it to 0 at the cutoff.

See pair_modify shift

Steve