you can do a few things:
- perform 2 runs, one starting from all replicas folded and 1 starting from
all replicas unfolded. All measures should match between the 2 runs.
anything that doesn't match is unreliable.
- do analysis on each replica (not each temperature) separately. Every
replica should have the same properties (rmsd histograms, fraction folded,
fraction of time at each temperature, and so on). If they don't, it's not
converged.
remember that "converged" is a bit vague since different properties
converge at different rates. Potential energy converges quickly, while
structure sampling does not. Every thing you measure needs error bars. Then
you know the precision, which is much more quantitative than a binary idea
of "converged" or not.
On Wed, Feb 6, 2013 at 7:21 AM, gargi borgohai <gargib2011.gmail.com> wrote:
> Dear AMBER users,
>
> I am using REMD method to study protein folding/unfolding behavior. I have
> run the system for 35ns. Can anyone help me regarding how can I confirm
> whether my system has converged or not?
>
>
> Thanking You..
> Sincerely
> Gargi Borgohain.
> PhD Scholar.
> Indian Institute of Technology, Guwahati.
> India.
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Feb 06 2013 - 05:00:02 PST
