On 6/15/16 8:45 AM, Pedro Lacerda wrote:
Hi list,
RMSD analysis of trajectories differs very much from the intuition gained
with the visual inspection of trajectories. RMSD has almost 4nm while the
video seems correct.
I'm using Amber 14ffSB with an metallic nonbonded ion attached on it. It
was converted from Amber files using ACPYPE, the initial structure seems
correct and the problem doesn't seems originated in parametrization errors.
Please can you point me some directions to debug this defect? I guess that
isn't needed but do I need to attach some relevant file?
Look at the picture of the animation and RMSD of the first 2ns trajectory:
http://imgur.com/a/bFCr8
Please wait for full loading the animated gif, the remaining of the .xtc
look like this first 20s of the gif. The RMSD was taken with the -nopbc
option! If taken without -[no]pbc flag the plot breaks like if PBC was not
considered. It looks like quite the opposite from the expected, right?
I also measured the RMSD using VMD giving the same results. Maybe this is a
very simple error, what can I doing wrong?
Use trjconv to recenter and make your molecules whole. That's standard practice
before computing RMSD. I have had mixed results with -nopbc and do not rely on it.
-Justin
--
==================================================
Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow
Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 629
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201
jalem...@outerbanks.umaryland.edu | (410) 706-7441
http://mackerell.umaryland.edu/~jalemkul
==================================================
--
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.
