<div dir="ltr">Hi Xavier,<br><br>I have a question about the accuracy of the correlation times obtained using gromacs tools. (I am getting these by integrating the output .xvg file of g_rotacf using g_analyze). Since the experimental values are not available to compare, I am not so sure about the correlation times I am getting. They are of the order of 0.4 to 0.8ns for various N-H vectors. I believe the values depends upon the range of time frames we use to fit with the legendre polynomials. I am using second order legendre polynomial to fit the frames. Can you give me any general suggestions to follow to obtain the reasonable correlation time values with MD simulations.<br>
<br>Ram. <br><br><div class="gmail_quote">On Thu, Sep 4, 2008 at 11:02 AM, Xavier Periole <span dir="ltr"><<a href="mailto:X.Periole@rug.nl">X.Periole@rug.nl</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
On Thu, 4 Sep 2008 09:45:34 -0400<div class="Ih2E3d"><br>
"rams rams" <<a href="mailto:rams.crux@gmail.com" target="_blank">rams.crux@gmail.com</a>> wrote:<br>
</div><div class="Ih2E3d"><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi Xavier,<br>
<br>
I am extremely sorry for incomplete information. But this is a follow up to<br>
my previous two emails. Unfortunately I havent received any suggestions for<br>
them and it appears like you too missed them. For more clarity I am pasting<br>
them here again:<br>
</blockquote></div>
Well if you do not get an answer that is probably that your post is not<br>
really asking a quick "answerable" question and you might want to<br>
re-formulate your post or dig a bit to make things easier.<div class="Ih2E3d"><br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
To calculate the rotational auto correlation functions, the command<br>
mentioned in the maual is:<br>
<br>
g_rotacf -P 1 -nparm 2 -fft -n index -o .xvg -fa -beginfit -endfit<br>
<br>
what are the -nparm and -fa options are meant for ?<br>
</blockquote></div>
This is an old exemple ... a bad one! Those options do not seem to exist<br>
anymore! Never knew they existed!<br>
you can grep "\-nparm" and "\-fa" in tools/*.c and you'll see they<br>
only appear in gmx_rotacf.c<br>
anyways These things are probably taken care of more invisibly.<br>
-nparm seems to give the number of parameters in the function to fit<br>
-fa the function to fit<br>
<br>
now: -nparm do not exist and -fa is -fitfn, most likely!<br>
<br>
This exemple should be removed and replaced! Just look at the options\listed when you type g_rotacf -h.<div class="Ih2E3d"><br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Also do we need to use -fitfn option to obtain the rotational auto<br>
correlation functions ??<br>
</blockquote></div>
No, -fitfn defines the type of function you'd use to fit the correlation<br>
function, which itself is written cause -o rotacf.xvg ...<br><font color="#888888">
<br>
XAvier.</font><div><div></div><div class="Wj3C7c"><br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
Ram.<br>
<br>
<br>
<br>
<br>
<br>
Dear users,<br>
<br>
I have given a command like the following to calculate the rotational auto<br>
correlation function:<br>
<br>
g_rotacf -f .trr -s .tpr -P 2 -fft -o .xvg -b 0000 -e 10000 -n .ndx -d<br>
<br>
I want to use the second order Legendre polynomial to fit. I integrated the<br>
resulting .xvg file, to obtain the correlation time using:<br>
<br>
g_analyze -f .xvg -integrate<br>
<br>
The output is the following. Could some one help me in understanding it<br>
betterly.<br>
<br>
Calculating the integral using the trapezium rule<br>
Integral 1 120.92840 +/- 0.00000<br>
std. dev. relative deviation of<br>
standard --------- cumulants from those of<br>
set average deviation sqrt(n-1) a Gaussian distribition<br>
cum. 3 cum. 4<br>
SS1 2.428352e-02 2.972849e-02 4.204243e-04 9.301 120.228<br>
<br>
<br>
Which value corresponds to correlation time and in which units (its in ps I<br>
suppose) ?<br>
<br>
Thanks in advance.<br>
<br>
Ram.<br>
<br>
<br>
On Thu, Sep 4, 2008 at 2:32 AM, Xavier Periole <<a href="mailto:X.Periole@rug.nl" target="_blank">X.Periole@rug.nl</a>> wrote:<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
On Wed, 3 Sep 2008 23:40:14 -0400<br>
"rams rams" <<a href="mailto:rams.crux@gmail.com" target="_blank">rams.crux@gmail.com</a>> wrote:<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi,<br>
<br>
I am so surprised for not finding any one who have better experience with<br>
g_rotacf. I have been playing around with it and the time correlation<br>
value<br>
I got by g_rotacf is so small in comparison to the time correlation value<br>
I<br>
calcualted using the hydrodynamic radius of the protein. The value is<br>
nearly<br>
10 times less. Can some one give me a better idea about g_rotacf.<br>
<br>
</blockquote>
Many people have certainly used g_rotacf to get ACFs of different<br>
observables.<br>
<br>
Anyways the way you describe your system, command line and your problem<br>
does<br>
does not help anyone to help you. Read your message above and think about<br>
what<br>
you would answer! You've played around with g_rotacf therefore you know it<br>
is not straightforward to give you the magic command.<br>
<br>
XAvier.<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
Ram.<br>
<br>
</blockquote>
<br>
-----------------------------------------------------<br>
XAvier Periole - PhD<br>
<br>
Molecular Dynamics Group / NMR and Computation<br>
University of Groningen<br>
The Netherlands<br>
-----------------------------------------------------<br>
_______________________________________________<br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.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 www interface<br>
or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">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><br>
<br>
</blockquote></blockquote>
<br>
-----------------------------------------------------<br>
XAvier Periole - PhD<br>
<br>
Molecular Dynamics Group / NMR and Computation<br>
University of Groningen<br>
The Netherlands<br>
-----------------------------------------------------<br>
_______________________________________________<br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.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 www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">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><br>
</div></div></blockquote></div><br></div>