constraints = all-angles does not convert angles to constraints - Redmine #3067
Archive from user: Marvin Bernhardt
I try to simulate united-atom propane gas (no LJ) with fixed C-C-C angle. For this I simply try to set constraints = all-angles in the mdp file, as this should:
Convert all bonds to constraints and all angles to bond-constraints.
(from http://manual.gromacs.org/2019.3/user-guide/mdp-options.html#bonds)
However It does not have an effect. This i conclude from measuring the RDF between the CH3 groups, and the above option does not change the width of the peak. Also i can see the peaks in the vibrational density of states.
In the md.log i find:
The number of constraints is 4000
but it should be 6000 because there are 2000 molecules.
grompp outputs
turning all bonds and angles into constraints...
but apparently it doesn’t happen. My current guess is that this is a grompp bug.
(from redmine: issue id 3067, created on 2019-08-26 by gmxdefault, closed on 2019-08-30)
- Changesets:
- Revision 9fee6b57 by Berk Hess on 2019-08-26T15:58:25Z:
Fix grompp constraints=all-angles with constraints
When using the mdp option constraints=all-angles, angles involving
bonds supplied as constraints in the topology would be removed,
but not replaced by angle constraints.
Fixes #3067
Change-Id: I7db390785a23c67abb30ae98e12272f88f5d5af7
- Uploads: