Your run is not converged. You could (i) run more slowly, (ii) run
lots of these non-equilibrium runs and treat them appropriately, or
(iii) use an equilibrium method.
In any event, I suggest that you begin by (iv) repeat the run for both
poses with different random seeds to see how much of the difference
has nothing to do with the starting structure. If the different
starting poses give you SMD results whose difference is no greater
than the difference between poses, then you're probably just looking
at statistical error between the two.
You asked if this error is ok... sorry, I can't judge that. But still
you need to know what kind of error it is (statistical, or
systematically based on starting pose). Note that you will always get
some difference between two runs, your goal is to have a small enough
difference that it doesn't change your conclusions.
Chris.
-- original message --
Hi,
I just want to know about the consistency we could obtain in plotting PMF.
I am using Gromacs 4.0.4 MPI version. Trying to find the PMF of the ligand
transport through the channel. I tried two SMDs (same pull rate, and force
constant) with different starting structures (only ligand poses and
position are slightly different). I performed the sampling of these two
pathways, the plots obtained are of similar overall shapes, but peaks (max
, min) have some variations?
Could you please give your comment on this attempt? Is that expectable?
Thank you,
--
Aswathy
--
gmx-users mailing list gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
