Dear all,
In my simulation, I am trying to use ‘compute temp/partial’ command under ‘npt’ to avoid the enforced
rigid velocity in z-direction.
the command i use for this purpose is like this:
compute new all temp/partial 1 1 0
fix 1 all npt 0.1 0.1 0.01 aniso 0.0 0.0 0.0 0.0 NULL NULL 10 drag 200
fix_modify 1 temp new
Then I found the temperature is keep increasing, rather than stay at my desired value and after a while,
‘lost atoms’ happens.
Is there something wrong with my setting temperature like this?
Thanks a lot for any suggestion!
Best
Ian