<html><head><style type="text/css"><!-- DIV {margin:0px;} --></style></head><body><div style="font-family:times new roman,new york,times,serif;font-size:12pt"><div><pre>Thanks for your help, Mark.<br><br>I cannot demux the trajectories, because I'm running each box at with a different temperature<br>and with a different initial configuration (appropiate to that temperature)<br>(This, combined with the replica exchange, is usually called "parallel tempering").<br><br>Checking the g_msd -h I saw that I can save some time using g_msd -e and -b (end and begin <br>time for the MSD measurement), but that's the best I can see right now.<br><br>I do see the discontinuities you mention, thats what I would like to avoid. And, I agree: if I <br>realize I need statistics on blocks longer than 100ps, I'm gonna have to try something else,<br>at least reducing the replica exchange frequency<br><br>Thanks again, bye,<br><br>Julian<br><br>Mark Wrote:<br><br>&gt;Julián
 GC wrote:<br>&gt;&gt;<i> <br></i>&gt;&gt;<i> <br></i>&gt;&gt;<i> Hi all,<br></i>&gt;&gt;<i> <br></i>&gt;&gt;<i> I'm studying small water clusters with Gromacs 3.3 (I'm also performing <br></i>&gt;&gt;<i> MC with a self-writed MC program).<br></i>&gt;&gt;<i> I run 24 boxes in paralell, with -replex 100000 (it is important to use <br></i>&gt;&gt;<i> replica exchange).<br></i>&gt;&gt;<i> <br></i>&gt;&gt;<i> I need to get MSD results, but I don't think g_msd takes into account <br></i>&gt;&gt;<i> that every 100ps the boxes are being exchanged.<br></i><br>&gt;Constant-T REMD trajectories are intrinsically discontinuous, so g_msd <br>&gt;will also have discontinuities. If all of your replicas started from the <br>&gt;same structure, then the MSD calculation is measuring from the same <br>&gt;starting point. If you need statistics for longer continuous chunks, <br>&gt;then you need to reconsider your approach. If you can adjust suitably <br>&gt;for
 temperature, then generating demux-ed constant-replica REMD <br>&gt;trajectories will at least eliminate the discontinuities. (demux.pl is <br>&gt;either in the source distribution or on the contributions section of the <br>&gt;GROMACS webpage, I forget which)<br><br>&gt;&gt;<i> Does anyone knows how to fix that? Have anyone had to deal with <br></i>&gt;<i>&gt; something like this before?<br></i>&gt;&gt;<i> (The obvious "easy" solution is to chop the trajectory in blocks of <br></i>&gt;<i>&gt; 100ps, calculate MSD in each block, and then average all the results, <br></i>&gt;&gt;<i> but would be nice if this annoying "shortcut" could be avoided)<br></i><br>&gt;See g_msd -h... I think -trestart saves you some legwork here.<br><br>&gt;Mark</pre></div></div><br>
      <hr size=1><br><font face="Verdana" size="-2">Tarjeta de crédito Yahoo! de Banco Supervielle.
Solicitá tu nueva Tarjeta de crédito. De tu PC directo a tu casa. <br> Visitá www.tuprimeratarjeta.com.ar</font></body></html>