On Mon, Oct 21, 2019, Hira Jabeen wrote:
>I am trying to run MD simulation in implicit solvent.
>In the initial nanoseconds, results are fine but after 90 nanoseconds I
>get an error of NaNs in coordinates and the simulation stops.
Welcome to MD! There are a variety of reasons why a simulation may
experience an instability, even after a long period where things seem
OK.
It's a good idea to visualize your trajectory, paying particular
attention to any non-standard resdidues or ligands you might have.
Also, study the temperature and energy profile that you can get from the
mdout files. Look for any differences between what you see in "good"
segments of the trajectory from "bad" ones.
It's possible that the difficulty goes back to the use of the CHARMM
force field, or to the conversion programs between CHARMM and Amber
formats. That seems somewhat unlikely, since the simulation ran OK for
a long time, and other users have not reported problems.
Aside: your input files are not using implicit solvent (not sure if that
is important or not.)
....dac
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Oct 23 2019 - 05:30:06 PDT
