I'm having a consistent problem in trying to run a Brownian Dynamics simulation of rigid particles on the order of 10 atoms. Basically after fixing the rigid bodies, i then apply:
fix integrator all nve
fix thermostat all langevin 1.000000 1.000000 1.000000 48279
as the simulation runs, the temperature drops from 1.0 to close to zero. I have used neigh_modify to turn off the force calculations between atoms within a rigid body. Basically, the simulation is not maintaining temperature. I have varied the damping coefficient from 1 to 2000 to no avail. Any suggestions on how to solve this problem?
Thanks in advance!