Yuchong: Can you send me an example of a small system that produces the
problem (similar to Andrew's)?
Andrew:
looking at your simulation, I see that you were using fix rigid single for
the immobile group, how about not using any time integrator for that
immobile group (e.g. by commenting out fix rigid); the potential energy will
be the same and you could save some redundant computing.
Good question. Glad to hear from you.
I noticed that "fix_rigid.cpp" modifies the virial (although I found
part of this code difficult to understand). So I am using fix "rigid"
because I also want to be able to run these kinds of simulations under
either NVT or NPT conditions with minimal changes to the input script.
(More on NPT later...)
I am hoping that this command:
"neigh_modify exclude group immobile immobile"
will have avoided the redundant computation. Please let me know if I'm wrong.
In this case, if you'd like to print out the temperature of the
mobile group, then define a temperature compute for it:
compute tempMobile mobile temp
thermo_style custom step c_tempMobile
Thanks. Although I am still using fix "rigid" and "neigh_modify
exclude", I tried printing out two different versions of the
temperature using:
"thermo_style custom step c_tempMobile temp"
There is a slight difference between the two temperatures.
"tempMobile" is allways about half a degree (0.1% or 0.2%) warmer than
"temp", however both temperatures appear to be fluctuating around the
target temperature (300 degrees +/- 7 degrees). I don't know if I
should worry about the discrepancy between these temperatures.
Overall I am relieved.
If, for some reasons, you must use fix rigid, then forces should also be
turned off, just like torques: force * off off off
Glad you brought this up...
For what it's worth, the rigid "immobile" group appears to be doing
what I asked it to do (translate but not rotate).
Please clarify: I've noticed you provided the user an option in
LAMMPS to turn off rotation and translation independently. Is there
any scenario where the user would turn off the torques, but leave the
forces on?
Thank you for the feedback by the way!
I appreciate it and I have more questions for you regarding NPT which
will be posted very soon.
Andrew