Hello GMX users,
I added another ligand as a new residue in the topology database, and
then build the topology of the liand with pdb2gmx,
everything worked well. But I got this warning from grompp for a nvt
equilibration after successful minimizations:
There are atoms at both ends of an angle, connected by constraints and
with masses that differ by more than a factor of 13. This means that
there are likely dynamic modes that are only very weakly coupled. To
ensure good equipartitioning, you need to either not use constraints on
all bonds (but, if possible, only on bonds involving hydrogens) or use
integrator = sd or decrease one or more tolerances:
verlet-buffer-tolerance <= 0.0001, LINCS iterations >= 2, LINCS order >=
4 or SHAKE tolerance <= 1e-05
It seems that there might be some problem of the topology file. But I
can run grompp without any warning if I change the
constraints on all bonds to only bonds involving hydrogens. Then I ran a
500 ps NVT equilibration, I checked the trajectory,
I did not found any wired angle. Could some one give me some tips? What
is the possible reason that I got the warning?
Thanks a lot!
PS: to be consistent, I also need to apply constraints on all bonds.
All the best,
Gromacs Users mailing list
* Please search the archive at
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a
mail to gmx-users-requ...@gromacs.org.