<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;">Thanks for the reply.&nbsp;<div>I used following command to extract first frame that has all the surfectants ( as suggested by first step for characterizing micelle clustering).</div><div><br></div><div>&nbsp;&nbsp; trjconv<span class="Apple-style-span" style="font-family: monospace; font-size: 14px; line-height: 17px; ">&nbsp;-f a.xtc -o a_cluster.gro -dump 0 -pbc cluster</span></div><div>&nbsp;Here is the output I got:</div><div><span class="Apple-style-span" style="font-family: monospace; font-size: 14px; line-height: 17px; "><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; padding-top: 0px; padding-right: 0px; padding-bottom: 0px; padding-left: 0px; line-height: 1.2em; outline-style: none; outline-width: initial; outline-color: initial; ">COM: &nbsp; &nbsp;0.000 &nbsp; &nbsp; 0.000 &nbsp; &nbsp;19.999 &nbsp;iter = 19526
 &nbsp;Isq = 12671033.000</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; padding-top: 0px; padding-right: 0px; padding-bottom: 0px; padding-left: 0px; line-height: 1.2em; outline-style: none; outline-width: initial; outline-color: initial; ">COM: &nbsp; &nbsp;0.000 &nbsp; &nbsp; 0.000 &nbsp; &nbsp; 0.000 &nbsp;iter = 19527 &nbsp;Isq = 12628369.000</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; padding-top: 0px; padding-right: 0px; padding-bottom: 0px; padding-left: 0px; line-height: 1.2em; outline-style: none; outline-width: initial; outline-color: initial; ">COM: &nbsp; &nbsp;0.000 &nbsp; &nbsp; 0.000 &nbsp; &nbsp;19.999 &nbsp;iter = 19528 &nbsp;Isq = 12671033.000</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; padding-top: 0px; padding-right: 0px; padding-bottom: 0px; padding-left: 0px; line-height: 1.2em; outline-style:
 none; outline-width: initial; outline-color: initial; ">COM: &nbsp; &nbsp;0.000 &nbsp; &nbsp; 0.000 &nbsp; &nbsp; 0.000 &nbsp;iter = 19529 &nbsp;Isq = 12628369.000</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; padding-top: 0px; padding-right: 0px; padding-bottom: 0px; padding-left: 0px; line-height: 1.2em; outline-style: none; outline-width: initial; outline-color: initial; ">COM: &nbsp; &nbsp;0.000 &nbsp; &nbsp; 0.000 &nbsp; &nbsp;19.999 &nbsp;iter = 19530 &nbsp;Isq = 12671033.000</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; padding-top: 0px; padding-right: 0px; padding-bottom: 0px; padding-left: 0px; line-height: 1.2em; outline-style: none; outline-width: initial; outline-color: initial; ">COM: &nbsp; &nbsp;0.000 &nbsp; &nbsp; 0.000 &nbsp; &nbsp; 0.000 &nbsp;iter = 19531 &nbsp;Isq = 12628369.000</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom:
 0px; margin-left: 0px; padding-top: 0px; padding-right: 0px; padding-bottom: 0px; padding-left: 0px; line-height: 1.2em; outline-style: none; outline-width: initial; outline-color: initial; ">COM: &nbsp; &nbsp;0.000 &nbsp; &nbsp; 0.000 &nbsp; &nbsp;19.999 &nbsp;iter = 19532 &nbsp;Isq = 12671033.000</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; padding-top: 0px; padding-right: 0px; padding-bottom: 0px; padding-left: 0px; line-height: 1.2em; outline-style: none; outline-width: initial; outline-color: initial; ">COM: &nbsp; &nbsp;0.000 &nbsp; &nbsp; 0.000 &nbsp; &nbsp; 0.000 &nbsp;iter = 19533 &nbsp;Isq = 12628369.000</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; padding-top: 0px; padding-right: 0px; padding-bottom: 0px; padding-left: 0px; line-height: 1.2em; outline-style: none; outline-width: initial; outline-color: initial; ">COM: &nbsp; &nbsp;0.000 &nbsp; &nbsp;
 0.000 &nbsp; &nbsp;19.999 &nbsp;iter = 19534 &nbsp;Isq = 12671033.000</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; padding-top: 0px; padding-right: 0px; padding-bottom: 0px; padding-left: 0px; line-height: 1.2em; outline-style: none; outline-width: initial; outline-color: initial; ">COM: &nbsp; &nbsp;0.000 &nbsp; &nbsp; 0.000 &nbsp; &nbsp; 0.000 &nbsp;iter = 19535 &nbsp;Isq = 12628369.000</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; padding-top: 0px; padding-right: 0px; padding-bottom: 0px; padding-left: 0px; line-height: 1.2em; outline-style: none; outline-width: initial; outline-color: initial; ">COM: &nbsp; &nbsp;0.000 &nbsp; &nbsp; 0.000 &nbsp; &nbsp;19.999 &nbsp;iter = 19536 &nbsp;Isq = 12671033.000</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; padding-top: 0px; padding-right: 0px; padding-bottom: 0px;
 padding-left: 0px; line-height: 1.2em; outline-style: none; outline-width: initial; outline-color: initial; ">COM: &nbsp; &nbsp;0.000 &nbsp; &nbsp; 0.000 &nbsp; &nbsp; 0.000 &nbsp;iter = 19537 &nbsp;Isq = 12628369.000</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; padding-top: 0px; padding-right: 0px; padding-bottom: 0px; padding-left: 0px; line-height: 1.2em; outline-style: none; outline-width: initial; outline-color: initial; ">COM: &nbsp; &nbsp;0.000 &nbsp; &nbsp; 0.000 &nbsp; &nbsp;19.999 &nbsp;iter = 19538 &nbsp;Isq = 12671033.000</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; padding-top: 0px; padding-right: 0px; padding-bottom: 0px; padding-left: 0px; line-height: 1.2em; outline-style: none; outline-width: initial; outline-color: initial; "><br></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; padding-top: 0px;
 padding-right: 0px; padding-bottom: 0px; padding-left: 0px; line-height: 1.2em; outline-style: none; outline-width: initial; outline-color: initial; "><span class="Apple-style-span" style="line-height: 17px; ">looks like lsq is only fluctuating between two numbers and the program is in an infinte loop.&nbsp;</span></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; padding-top: 0px; padding-right: 0px; padding-bottom: 0px; padding-left: 0px; line-height: 1.2em; outline-style: none; outline-width: initial; outline-color: initial; "><span class="Apple-style-span" style="line-height: 17px; "><br></span></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; padding-top: 0px; padding-right: 0px; padding-bottom: 0px; padding-left: 0px; line-height: 1.2em; outline-style: none; outline-width: initial; outline-color: initial; "><span class="Apple-style-span" style="line-height: 17px;
 ">Jagannath&nbsp;</span></div></span></div><div><br>--- On <b>Wed, 7/7/10, Mark Abraham <i>&lt;mark.abraham@anu.edu.au&gt;</i></b> wrote:<br><blockquote style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;"><br>From: Mark Abraham &lt;mark.abraham@anu.edu.au&gt;<br>Subject: Re: [gmx-users] problem with trjconv -pbc cluster<br>To: "Discussion list for GROMACS users" &lt;gmx-users@gromacs.org&gt;<br>Date: Wednesday, 7 July, 2010, 12:20 PM<br><br><div id="yiv1635627800"><br><br>----- Original Message -----<br>From: jagannath mondal &lt;jmondal2004@yahoo.co.in&gt;<br>Date: Wednesday, July 7, 2010 16:36<br>Subject: [gmx-users] problem with trjconv -pbc cluster<br>To: gmx-users@gromacs.org<br><br><table border="0" cellpadding="0" cellspacing="0"><tbody><tr><td
 style="font-family:inherit;font-style:inherit;font-variant:inherit;font-weight:inherit;font-size:inherit;line-height:inherit;font-size-adjust:inherit;font-stretch:inherit;" valign="top"><font style="font-style:normal;font-weight:normal;background-color:rgb(245, 248, 240);font-size:14px;">&gt; </font>Hi,<div>&nbsp;&nbsp;I had a system of surfectants which are started with an initial configuration where all of them are well dispersed. Visual study of trajectory shows they start aggregating and finally form two discrete micelles.&nbsp;</div><div><font style="font-style:normal;font-weight:normal;background-color:rgb(245, 248, 240);font-size:14px;">&gt; </font>To quantify this micelle clusterization, I tried to use the suggestions present in the gromacs documentations:</div><div><font style="font-style:normal;font-weight:normal;background-color:rgb(245, 248, 240);font-size:14px;">&gt;
 </font>&nbsp;http://www.gromacs.org/Documentation/How-tos/Micelle_Clustering</div><div><font style="font-style:normal;font-weight:normal;background-color:rgb(245, 248, 240);font-size:14px;">&gt; </font><br></div><div><font style="font-style:normal;font-weight:normal;background-color:rgb(245, 248, 240);font-size:14px;">&gt; </font>i.e&nbsp;</div><div><font style="font-style:normal;font-weight:normal;background-color:rgb(245, 248, 240);font-size:14px;">&gt; </font>&nbsp;&nbsp;&nbsp;</div><span class="Apple-style-span" style="font-family:Verdana, sans-serif;font-size:14px;"><ol><li style="line-height:1.25;"><font style="font-style:normal;font-weight:normal;background-color:rgb(245, 248, 240);font-size:14px;">&gt; </font>use&nbsp;<code><a rel="nofollow" target="_blank" href="http://www.gromacs.org/Documentation/Gromacs_Utilities/trjconv" style="color:rgb(0, 84, 166);">trjconv</a>&nbsp;-pbc cluster</code>&nbsp;to obtain a single frame that has all of the
 lipids in the  unit cell. This must be the first frame of your trajectory. A similar frame from some previous timepoint will not work.</li><li style="line-height:1.25;"><font style="font-style:normal;font-weight:normal;background-color:rgb(245, 248, 240);font-size:14px;">&gt; </font>use&nbsp;<a rel="nofollow" target="_blank" href="http://www.gromacs.org/Documentation/Gromacs_Utilities/grompp" style="color:rgb(0, 84, 166);">grompp</a>&nbsp;to make a new tpr file based on the frame that was output from the step above.</li><li style="line-height:1.25;"><font style="font-style:normal;font-weight:normal;background-color:rgb(245, 248, 240);font-size:14px;">&gt; </font>use&nbsp;<code><a rel="nofollow" target="_blank" href="http://www.gromacs.org/Documentation/Gromacs_Utilities/trjconv" style="color:rgb(0, 84, 166);">trjconv</a>&nbsp;-pbc nojump</code>&nbsp;to produce the desired trajectory using the newly produced tpr file.</li></ol></span><div><span
 class="Apple-style-span" style="font-family:Verdana, sans-serif;font-size:14px;"><ol><li style="line-height:1.25;"><code><a rel="nofollow" target="_blank" href="http://www.gromacs.org/Documentation/Gromacs_Utilities/trjconv" style="color:rgb(0, 84, 166);"><font style="font-style:normal;font-weight:normal;background-color:rgb(245, 248, 240);font-size:14px;">&gt; </font>trjconv</a>&nbsp;-f a.xtc -o  a_cluster.gro -e 0.001 -pbc cluster</code></li><li style="line-height:1.25;"><code><a rel="nofollow" target="_blank" href="http://www.gromacs.org/Documentation/Gromacs_Utilities/grompp" style="color:rgb(0, 84, 166);"><font style="font-style:normal;font-weight:normal;background-color:rgb(245, 248, 240);font-size:14px;">&gt; </font>grompp</a>&nbsp;-f a.mdp -c a_cluster.gro -o a_cluster.tpr</code></li><li style="line-height:1.25;"><code><a rel="nofollow" target="_blank" href="http://www.gromacs.org/Documentation/Gromacs_Utilities/trjconv" style="color:rgb(0, 84,
 166);"><font style="font-style:normal;font-weight:normal;background-color:rgb(245, 248, 240);font-size:14px;">&gt; </font>trjconv</a>&nbsp;-f a.xtc -o a_cluster.xtc -s a_cluster.tpr -pbc nojump</code></li></ol><div><font class="Apple-style-span" face="monospace"><span class="Apple-style-span" style="line-height:17px;"><br></span></font></div><div><font class="Apple-style-span" face="monospace"><span class="Apple-style-span" style="line-height:17px;"><font style="font-style:normal;font-weight:normal;background-color:rgb(245, 248, 240);font-size:14px;">&gt; </font>But, The first step i.e. trjconv -pbc cluster does not work. Looks like it is going through an infinite loop and is not stopping for convergence.<br><br>So what command did you give, what output did you get, and does the command make sense with respect to your trajectory contents?<br><br>Mark<br><br></span></font></div><div><font class="Apple-style-span" face="monospace"><span
 class="Apple-style-span" style="line-height:17px;"><font style="font-style:normal;font-weight:normal;background-color:rgb(245, 248, 240);font-size:14px;">&gt; </font><br></span></font></div><div><font class="Apple-style-span" face="monospace"><span class="Apple-style-span" style="line-height:17px;"><font style="font-style:normal;font-weight:normal;background-color:rgb(245, 248, 240);font-size:14px;">&gt; </font>I am using gromacs 4.0.7</span></font></div><div><font class="Apple-style-span" face="monospace"><span class="Apple-style-span" style="line-height:17px;"><font style="font-style:normal;font-weight:normal;background-color:rgb(245, 248, 240);font-size:14px;">&gt; </font><br></span></font></div><div><font class="Apple-style-span" face="monospace"><span class="Apple-style-span" style="line-height:17px;"><font style="font-style:normal;font-weight:normal;background-color:rgb(245, 248, 240);font-size:14px;">&gt; </font>Any help will be
 appreciated.</span></font></div><div><font class="Apple-style-span" face="monospace"><span class="Apple-style-span" style="line-height:17px;"><font style="font-style:normal;font-weight:normal;background-color:rgb(245, 248, 240);font-size:14px;">&gt; </font><br></span></font></div><div><font class="Apple-style-span" face="monospace"><span class="Apple-style-span" style="line-height:17px;"><font style="font-style:normal;font-weight:normal;background-color:rgb(245, 248, 240);font-size:14px;">&gt; </font>Jagannath</span></font></div></span></div></td></tr></tbody></table><font style="font-style:normal;font-weight:normal;background-color:rgb(245, 248, 240);font-size:14px;">&gt; </font><br>&gt; -- <br>&gt; gmx-users mailing list&nbsp;&nbsp;&nbsp; gmx-users@gromacs.org<br>&gt; http://lists.gromacs.org/mailman/listinfo/gmx-users<br>&gt; Please search the archive at http://www.gromacs.org/search <br>&gt; before posting!<br>&gt; Please don't post (un)subscribe
 requests to the list. Use the <br>&gt; www interface or send it to gmx-users-request@gromacs.org.<br>&gt; Can't post? Read http://www.gromacs.org/mailing_lists/users.php
</div><br>-----Inline Attachment Follows-----<br><br><div class="plainMail">-- <br>gmx-users mailing list&nbsp; &nbsp; <a ymailto="mailto:gmx-users@gromacs.org" href="/mc/compose?to=gmx-users@gromacs.org">gmx-users@gromacs.org</a><br><a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before posting!<br>Please don't post (un)subscribe requests to the list. Use the <br>www interface or send it to <a ymailto="mailto:gmx-users-request@gromacs.org" href="/mc/compose?to=gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>Can't post? Read <a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a></div></blockquote></div></td></tr></table><br>