Ok Justin,<br><br>So which gmx tool i can use to obtain the angle between two vectors defined by two atom pairs ?<br><br>SA<br><br><br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">

sa wrote:<br>
&gt; Hi all<br>
&gt;<br>
&gt; This message is related with my previous message (see below) about the<br>
&gt; calculation of the order value for the DPC alkyl chain in a micelle. So<br>
&gt; if I understand well the previous angel&#39;s message, I need to compute the<br>
&gt; theta angle between the 2 vectors defined by the (i-1,i+1) and (i, i+2).<br>
&gt; So for carbon atom C14, I have created an index file with two groups<br>
&gt; defined as following<br>
&gt;<br>
&gt; aC14 | aC16  and aC14 | aC16 with make_ndx_mpi<br>
&gt;<br>
&gt; and used the command g_sgangle_mpi -s run_1.tpr -f *.xtc -b 100000 -e<br>
&gt; 101000 -n index.ndx and choose the C13_C15 and C14_C16 groups.<br>
&gt; Unfortunately i obtain the following error<br>
&gt;<br>
&gt; Group C13_C15 contains the following atoms:<br>
&gt; Atomname 0: C13<br>
&gt; Atomname 1: C15<br>
&gt; Atomname 2: C13<br>
&gt; Atomname 3: C15<br>
&gt; Atomname 4: C13<br>
&gt; Atomname 5: C15<br>
&gt; Atomname 6: C13<br>
&gt; Atomname 7: C15<br>
&gt; Atomname 8: C13<br>
&gt; Atomname 9: C1<br>
&gt; ...<br>
&gt;<br>
&gt; Group C14_C16 contains the following atoms:<br>
&gt; Atomname 0: C14<br>
&gt; Atomname 1: C16<br>
&gt; Atomname 2: C14<br>
&gt; Atomname 3: C16<br>
&gt; Atomname 4: C14<br>
&gt;<br>
&gt; ...<br>
&gt; Atomname 105: C16<br>
&gt; Atomname 106: C14<br>
&gt; Atomname 107: C16<br>
&gt;<br>
&gt; Careful: distance only makes sense in some situations.<br>
&gt;<br>
&gt; Reading frame       0 time 100000.000<br>
&gt; Back Off! I just backed up sg_angle.xvg to ./#sg_angle.xvg.2#<br>
&gt;<br>
&gt; -------------------------------------------------------<br>
&gt; Program g_sgangle_mpi, VERSION 4.5.3<br>
&gt; Source code file: gmx_sgangle.c, line: 127<br>
&gt;<br>
&gt; Fatal error:<br>
&gt; Something wrong with contents of index file.<br>
&gt;<br>
&gt; For more information and tips for troubleshooting, please check the GROMACS<br>
&gt; website at <a href="http://www.gromacs.org/Documentation/Errors" target="_blank">http://www.gromacs.org/Documentation/Errors</a><br>
&gt; -------------------------------------------------------<br>
&gt;<br>
&gt; Did I make a error, what is the correct approach. to obtain the angle<br>
&gt; between 2 vectors ?<br>
&gt;<br>
<br>
Per the description in g_sgangle -h, groups can only be two or three atoms in<br>
size.  It seems as if this tool would be very inefficient for any more than a<br>
few calculations.<br>
<br>
-Justin<br>
<br>
&gt; Thank you in advance<br>
&gt;<br>
&gt; SA<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;     ------------------------------<br>
&gt;<br>
&gt;     Message: 3<br>
&gt;     Date: Wed, 16 Mar 2011 10:06:44 +0100<br>
&gt;     From: ?ngel Pi?eiro &lt;<a href="mailto:angel.pineiro@usc.es">angel.pineiro@usc.es</a> &lt;mailto:<a href="mailto:angel.pineiro@usc.es">angel.pineiro@usc.es</a>&gt;&gt;<br>
&gt;     Subject: Re: [gmx-users] Re: gmx-users Digest, Vol 83, Issue 106<br>
&gt;     To: Discussion list for GROMACS users &lt;<a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
&gt;     &lt;mailto:<a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a>&gt;&gt;<br>
&gt;     Message-ID: &lt;1300266404.6843.29.camel@arginine&gt;<br>
&gt;     Content-Type: text/plain; charset=&quot;utf-8&quot;<br>
&gt;<br>
&gt;     Then, please, let us know how it works for your systems. The results for<br>
&gt;     my systems were exactly as expected. This allows to evaluate the order<br>
&gt;     of the C-chains regardless their orientation... but I do not know if<br>
&gt;     there is a better method to do this. I would be happy to know the<br>
&gt;     opinion of anyone else who want to try this or to propose an alternative<br>
&gt;     method.<br>
&gt;<br>
&gt;     Cheers,<br>
&gt;<br>
&gt;     Ã ngel.<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;     On Wed, 2011-03-16 at 09:42 +0100, sa wrote:<br>
&gt;<br>
&gt;      &gt; Thank you Angel, I will try your suggestion.<br>
&gt;      &gt;<br>
&gt;      &gt; Cheers<br>
&gt;      &gt;<br>
&gt;      &gt; SA<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt;         Hi<br>
&gt;      &gt;         very recently I faced the same problem with a system that<br>
&gt;      &gt;         gives micelles<br>
&gt;      &gt;         of different geometries and, as far as I saw, g_order<br>
&gt;     don&#39;t do<br>
&gt;      &gt;         that.<br>
&gt;      &gt;         Then I decided to compute a kind of local order parameters<br>
&gt;      &gt;         defined as:<br>
&gt;      &gt;<br>
&gt;      &gt;         S_i=(3 cos(\theta)-1)/2<br>
&gt;      &gt;<br>
&gt;      &gt;         where theta is the angle between the segments joining the<br>
&gt;      &gt;         carbon atoms<br>
&gt;      &gt;         (i-1,i+1) and (i, i+2) in a linear C-chain. Perhaps you find<br>
&gt;      &gt;         this<br>
&gt;      &gt;         reasonable for your analysis...<br>
&gt;      &gt;<br>
&gt;      &gt;         Cheers,<br>
&gt;      &gt;<br>
&gt;      &gt;         Ã ngel.<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;<br>
&gt;<br>
<br>
--<br>
========================================<br>
<br>
Justin A. Lemkul<br>
Ph.D. Candidate<br>
ICTAS Doctoral Scholar<br>
MILES-IGERT Trainee<br>
Department of Biochemistry<br>
Virginia Tech<br>
Blacksburg, VA<br>
jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> | (540) 231-9080<br>
<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br>
<br>
========================================<br>
<font color="#888888"><br>
</font></blockquote></div><br>