I am running simulations of bead-spring polymers in the glassy phase (LJ
interactions, FENE bonds) and I find that my simulation box acquires a
constant centre of mass velocity upon cooling. I don't know if this is a
bug exactly, but it is certainly very annoying. I've so far been
modifying the fix_msd code every time I download a new version, but I'm
hoping you can help me find a more permanent solution. I've attached the
data file I use, and the input file. Let me know if there is anything
else you need. Any suggestions would be much appreciated.
p.s. It might be useful to put a link to this mailing list in the
documentation on how to report bugs!
in.test (1.04 KB)
T03.data (76.6 KB)