[lammps-users] fix rigid and timestep


There are some rigids in my system when I do DPD simulation. We know we can integrate the equation of motion with timestep 0.05 in reduced LJ units. BUT, I can’t confirm whether it is suitable to " fix rigid command" because it use constant-energy time integration.

Any suggestion?

Thanks in advance.


Try it and see what happens. If you conserve energy, good, if
not, see what happens if you reduce the timestep.