It might be worse tan you think.
The amount of time you need to spend during what you call
'equilibration' depends very much on how fast you pull during that
pulling part of the simulation. For example, if you pull 1 A very
slowly, then you do not need to spend much time equilibrating. If you
pull very fast, then you better spend lots of time equilibrating
(technically, relaxing).
So, it is not only the relaxation time of your system you need to worry
about, but more importantly, you need to figure out how far out of
equilibrium you are after pulling.
Unfortunately, there are no formulas or set recipes for any of this, you
need to try different ways.
adrian
On 6/20/17 8:28 AM, David Case wrote:
> On Tue, Jun 20, 2017, Ramin Salimi wrote:
>> I am trying to pull my DNA molecule in increments of 1 A0, and then
>> run equilibration for 50ps, then do another pull for 1 A0, then do the
>> 50ps equilibration, and keep doing it until the molecule reaches my
>> target length. My primary purpose is to keep the system in quasi-static
>> condition as much as possible.
>>
>> So, based upon what I understood, to keep the system in quasi-static
>> during both pullings, and equilibration periods, the equilibration
>> periods between increments along with pulling periods should be larger
>> than the relaxation time of the system.
> Generally, there is no one number that qualifies as "the" relaxation time for
> a system. I'm no expert on pulling simulations, but my guess is that you will
> have to try different equilibration periods to examine convergence.
>> Does it have anything to do with correlation time, and the use of IRED
>> vectors, and IRED matrices in cpptraj?
> No.
>
> ....dac
>
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
--
Dr. Adrian E. Roitberg
University of Florida Research Foundation Professor
Department of Chemistry
University of Florida
roitberg.ufl.edu
352-392-6972
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Tue Jun 20 2017 - 06:00:03 PDT
