[gmx-users] System drift with velocity verlet integrator
Kristoffer Enøe Johansson
kristoffer.johansson at sund.ku.dk
Wed Nov 19 10:43:11 CET 2014
Dear GMX users,
I am tryng to simulate a solid molecular crystal but the system builds up translational movement when I use the md-vv integrator with or without pressure coupling. With the leapfrog integrator this never happens. The md-vv integrator is interesting for me for comparison with another software.
Without symmetry constraints translation is indeed possible with close to zero friction for a crystallografic (periodic) cell but the center of mass removal should prevent this from happening.
Setting nstcomm=1 and comm-grps to a central part of the system had no effect at all. In using gromacs 4.6.7.
Any comments or ideas are welcome since I am stuck with this. I can provide any input/output files on request.
Computational Materials Science
Department of Pharmacy
University of Copenhagen
More information about the gromacs.org_gmx-users