<br><br><div><span class="gmail_quote"></span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><br><br>---------- Пересланное письмо ----------
<br>From:&nbsp;David van der Spoel &lt;<a href="mailto:spoel@xray.bmc.uu.se">spoel@xray.bmc.uu.se</a>&gt;<br>To:&nbsp;Discussion list for GROMACS users &lt;<a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a>&gt;<br>Date:&nbsp;Fri, 26 Jan 2007 10:02:01 +0100
<br>Subject:&nbsp;Re: [gmx-users] g_gyrate<br>Dmitriy Golubobsky wrote:<br>&gt; Dear Gurus,<br>&gt; I&#39;ve got question on Radius of gyrtion calculation.<br>&gt; my system consists of 1 polymer and a solvent.<br>&gt; I did MD and decide to calculate Rg for this system
<br>&gt; I made index group consisting my residues<br>&gt; g_gyrate -f xtc_file -s tpr_file -n ndx_file -p -o xvg_output<br>&gt; ok. I&#39;ve got the output file with 4 colums: time Rg Rx Ry Rz<br>&gt; I use -p option to get components in principal axes, because I&#39; d like
<br>&gt; to calculate form-factors for my system.<br>&gt;<br>&gt; as I undersand Rg = sqrt ( sum_{i}^{N} r_i*m_i)/M)<br>&gt; ok. then fx=Rx/Rg^2, fy=Ry/Rg^2, fz=Rz/Rg^2<br>&gt; and I &nbsp;wish to get fx+fy+fz=1 but it is not so.
<br>&gt; where I&#39;m mistaken?<br>&gt; How does components Rxyz in principal axes are calculated?<br>&gt; I&#39;m using gromacs 3.3.1<br>&gt; Thanks a lot.<br>Rg = sqrt ( sum_{i}^{N} (r_i^2)*m_i)/M)<br><br>--<br>David.<br>
<br></blockquote></div>Dear David,<br>I understand this, thus for the <br>Rx=Rg = sqrt ( sum_{i}^{N} (x_i^2)*m_i)/M)<br><br>and we have to get Rx^2+Ry^2+Rz^2=Rg^2<br>but we didn&#39;t get this.<br>