[lammps-users] fix setforce

Dear,

velocity lower set -1.0 0.0 0.0 units box
velocity upper set 1.0 0.0 0.0 units box

fix 1 all nvt 1.0 1.0 1.0

fix lf lower setforce 0.0 0.0 0.0
fix uf upper setforce 0.0 0.0 0.0

The fix setfoce command can’t fix the lower and upper atoms’ velocity in above codes, What’s wrong? If I change fix nvt to fix nve and fix temp/rescale, the fix setforce can work. SO, I think there is a bug in fix nvt.

Any suggestion? Thanks

Fix nvt will change the velocity of your wall atoms, which
probably isn't what you want. So don't do NVT on the wall
atoms.

Steve