Hello all,
I am trying to run a simulation with two of the molecules being restrained
to their original positions. However after I visualize the trajectory with
VMD, the restrained molecules show a drift from their current position over
time.
The prmtop, inpcrd, and in files can be found here:
https://www.dropbox.com/sh/ecng1mjobglgy6i/AAAIydLyGWe2O-e-2PnS5M-sa?dl=0
My procedure is the following:
1) minimize, 2) heat at NPT conditions
Here are the commands I used :
sander -O -i min.in -o min-file.out -p file.prmtop -c file.inpcrd \
-r min-file.rst -inf min-file.mdinfo -ref file.inpcrd
sander -O -i heat.in -o heat-file.out -p file.prmtop -c min-file.rst -r
heat-file.rst \
-x heat-file.mdcrd -inf heat-file.mdinfo -ref file.inpcrd
I have found that someone was having similar issues before:
http://archive.ambermd.org/201302/0115.html
I attempted to use the suggestion by setting nscm to be larger than nstlim
but this problem still persists.
Does anybody have an explanation for what is the cause of the drift and/or
how to solve this issue? Thank-you.
Best wishes,
Sanmeet
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Mar 19 2015 - 16:00:02 PDT