On Fri, Jul 12, 2019, Gustaf Olsson wrote:
>
>I am running three different systems that contain pehnylphosdhonic acid,
>the neutral form (ppa) the anion (pp1) and the double negative anion
>(pp2).
I bet that one of the protons on the phosphorous is getting very close
to one of the other phosphate oxygens. At least, this is a common
problem with protonated phosphates, although it's not clear why ppa
should be different than ppa1.
Are you using gaff2? One of the important changes between gaff and
gaff2 is that atom type "ho" is given non-zero van der Waals radii,
which should prevent the hydrogen from getting too close to a
neighboring oxygen. You may also need to increase (say double) the
P-O-H angle constant.
Since the problem happens quickly, run a short trajectory with ntwx=1 or
2, and visualize the geometry of the phopshate group as it gets close to
the point where the system crashes.
...hope this helps...problem could be something different, of
course....dac
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Sun Jul 14 2019 - 06:00:02 PDT
