[gmx-users] g_msd, REMD, trestart...

[Julián GC]

Thanks for your help, Mark.

I cannot demux the trajectories, because I'm running each box at with a different temperature
and with a different initial configuration (appropiate to that temperature)
(This, combined with the replica exchange, is usually called "parallel tempering").

Checking the g_msd -h I saw that I can save some time using g_msd -e and -b (end and begin 
time for the MSD measurement), but that's the best I can see right now.

I do see the discontinuities you mention, thats what I would like to avoid. And, I agree: if I 
realize I need statistics on blocks longer than 100ps, I'm gonna have to try something else,
at least reducing the replica exchange frequency

Thanks again, bye,

Julian

Mark Wrote:

>Julián GC wrote:
>> 
>> 
>> Hi all,
>> 
>> I'm studying small water clusters with Gromacs 3.3 (I'm also performing 
>> MC with a self-writed MC program).
>> I run 24 boxes in paralell, with -replex 100000 (it is important to use 
>> replica exchange).
>> 
>> I need to get MSD results, but I don't think g_msd takes into account 
>> that every 100ps the boxes are being exchanged.

>Constant-T REMD trajectories are intrinsically discontinuous, so g_msd 
>will also have discontinuities. If all of your replicas started from the 
>same structure, then the MSD calculation is measuring from the same 
>starting point. If you need statistics for longer continuous chunks, 
>then you need to reconsider your approach. If you can adjust suitably 
>for temperature, then generating demux-ed constant-replica REMD 
>trajectories will at least eliminate the discontinuities. (demux.pl is 
>either in the source distribution or on the contributions section of the 
>GROMACS webpage, I forget which)

>> Does anyone knows how to fix that? Have anyone had to deal with 
>> something like this before?
>> (The obvious "easy" solution is to chop the trajectory in blocks of 
>> 100ps, calculate MSD in each block, and then average all the results, 
>> but would be nice if this annoying "shortcut" could be avoided)

>See g_msd -h... I think -trestart saves you some legwork here.

>Mark



      Yahoo! Encuentros.

Ahora encontrar pareja es mucho más fácil, probá el nuevo Yahoo! Encuentros http://yahoo.cupidovirtual.com/servlet/NewRegistration
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20080315/fa67f59e/attachment.html>