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