GROMACS version: 5.1.4
Recently I am learning how to choose md parameters and I encounter some problems.
1.According to gromacs manual of mdp options, “Run control - comm-mode” controls restriction on the center of mass motion. The default setting is “None”(so as the minim.mdp in “Lysozyme in water” tutorial), but I found some minim.mdp files in the paper which set “comm-mode = Linear” to remove center of mass translational velocity. I wonder when this approach should be adopted and why.
- “Pressure coupling - refcoord_scaling” controls whether the reference coordinates for position restraints are scaled with the scaling matrix of the pressure coupling.
According to the manual, "no: The reference coordinates for position restraints are not modified. Note that with this option the virial and pressure will depend on the absolute positions of the reference coordinates.
all: The reference coordinates are scaled with the scaling matrix of the pressure coupling.
com: Scale the center of mass of the reference coordinates with the scaling matrix of the pressure coupling. The vectors of each reference coordinate to the center of mass are not scaled. Only one COM is used, even when there are multiple molecules with position restraints. For calculating the COM of the reference coordinates in the starting configuration, periodic boundary conditions are not taken into account. "
I noticed that some npt.mdp files choose “com” while some do not. I wonder when should we scale the reference coordinates and why.
Thanks in advance!