[gmx-users] Exchange interval in REMD

Gianluca Interlandi gianluca at u.washington.edu
Thu Jan 26 07:48:27 CET 2012


Hi David,

I am actually performing Simulated Tempering which is the same as RE 
except that it is in series. But I guess that the same thinking concerning 
exchange interval still applies.

Gianluca

On Thu, 26 Jan 2012, David van der Spoel wrote:

> On 2012-01-26 06:37, Gianluca Interlandi wrote:
>> Hi Mark,
>> 
>> Thanks for the answer. Just to be clear: you think that a too short
>> interval might introduce artefacts? I agree with this. So, do you think
>> that an interval between 1 and 10 ps could be a good value?
> It also depends whether or not you use constant pressure, in which case it 
> makes sense to increase the time to long enough to let the volume relax.
> I still do not understand why people do NVT REMD, because it makes all but 
> one replica have a pressure that is not the ambient pressure.
>
>> 
>> Gianluca
>> 
>> 
>> On Thu, 26 Jan 2012, Mark Abraham wrote:
>> 
>>> On 26/01/2012 11:38 AM, Gianluca Interlandi wrote:
>>>> Hi,
>>>> 
>>>> I have been going through the How-to for REMD in gromacs:
>>>> 
>>>> http://www.gromacs.org/Documentation/How-tos/REMD
>>>> 
>>>> Is there a recommended value for the exchange interval and does it
>>>> depend on system size? For example, is an interval of 100 ps
>>>> appropriate for a sytem with 60,000 atoms? Is a too long interval
>>>> simply a waste of time or can it actually introduce artefacts?
>>> 
>>> There's literature by Sindhikara (sp?) and Roitberg that says you
>>> should have a short an interval as possible, and other literature that
>>> says you shouldn't go shorter than around the autocorrelation time of
>>> PE, ~1ps. I don't buy the former, and think some of their data support
>>> the latter, and have some of my own not-ready-for-publication data
>>> that I think suggests artefacts from the former. Work from Hansmann et
>>> al suggests that achieving steady-state replica flow is rather tricky
>>> even for seemingly trivial peptides, so increasing the number of
>>> exchange attempts might help a little. I'm not aware of any good case
>>> for a long interval between attempts.
>>> 
>>> Mark
>>> --
>>> gmx-users mailing list gmx-users at gromacs.org
>>> http://lists.gromacs.org/mailman/listinfo/gmx-users
>>> Please search the archive at
>>> http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
>>> Please don't post (un)subscribe requests to the list. Use the www
>>> interface or send it to gmx-users-request at gromacs.org.
>>> Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>>> 
>> 
>> -----------------------------------------------------
>> Gianluca Interlandi, PhD gianluca at u.washington.edu
>> +1 (206) 685 4435
>> http://artemide.bioeng.washington.edu/
>> 
>> Postdoc at the Department of Bioengineering
>> at the University of Washington, Seattle WA U.S.A.
>> -----------------------------------------------------
>
>
> -- 
> David van der Spoel, Ph.D., Professor of Biology
> Dept. of Cell & Molec. Biol., Uppsala University.
> Box 596, 75124 Uppsala, Sweden. Phone:	+46184714205.
> spoel at xray.bmc.uu.se    http://folding.bmc.uu.se
> -- 
> gmx-users mailing list    gmx-users at gromacs.org
> http://lists.gromacs.org/mailman/listinfo/gmx-users
> Please search the archive at 
> http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
> Please don't post (un)subscribe requests to the list. Use the www interface 
> or send it to gmx-users-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>

-----------------------------------------------------
Gianluca Interlandi, PhD gianluca at u.washington.edu
                     +1 (206) 685 4435
                     http://artemide.bioeng.washington.edu/

Postdoc at the Department of Bioengineering
at the University of Washington, Seattle WA U.S.A.
-----------------------------------------------------



More information about the gromacs.org_gmx-users mailing list