Hi Bharat and Matt,<br>
<br>
The tutorial was set up by a colleague of mine, who has left our group
a while ago. Given your remarks, I think it's better to replace the
current version with a new one. I hope I can find some time for the
revision and would welcome any suggestions for changes, additions
(topics to cover), which can be sent to me off the list. Consider the
tutorial in its current form deprecated (which I will add as a note in
the tutorial soon).<br>
<br>
Best regards,<br>
<br>
Tsjerk<br><br><div><span class="gmail_quote">On 6/2/06, <b class="gmail_sendername"><a href="mailto:mernst@tricity.wsu.edu">mernst@tricity.wsu.edu</a></b> &lt;<a href="mailto:mernst@tricity.wsu.edu">mernst@tricity.wsu.edu
</a>&gt; wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Hi Bharat,<br><br>This is the same tutorial that led me astray when I first began trying free energy
<br>calculations in Gromacs. I am not sure how close to the experimental value the author<br>intended to get with this hydration free energy tutorial. If you are using a recent<br>version of Gromacs, be aware that the soft core method has changed. The tutorial says
<br>sc-alpha should be 1.51, but that is more appropriate for the older method corresponding<br>to sc-power=2 rather than the sc-power=1 that is now recommended.<br><br>This tutorial also has you simultaneously mutate partial charges and atom types. I have
<br>recently been informed, and confirmed via observation, that performing the partial<br>charge mutation with soft core potentials on can lead to misleading and very noisy<br>results. Look in the recent mailing list archives for suggestions and papers about
<br>splitting the transformation into two series of calculations, one with partial charge<br>mutation and no soft core potentials and the other with soft core potentials and no<br>partial charge mutation. Even if you don't want to do this for the simple tutorial
<br>you'll probably want to do it for real systems in the future.<br><br>Finally, the tutorial gives the dubious advice to use the final structure from one<br>window as the starting structure for the next window. If you instead run each window's
<br>simulations completely independently, you can simultaneously run calculations on several<br>machines and don't have to worry about the output of previous windows somehow biasing or<br>contaminating the current window. Again, I don't know if it is necessary to operate this
<br>way to get good results on this simple tutorial, but it does become more important if<br>you want to run novel, non-trivial systems in the future.<br><br>Matt Ernst<br>Washington State University<br><br>&gt; Dear GMX users,
<br>&gt;<br>&gt; i am trying 'relative hydration free energy of p-cresol wrt toluene' tutorial<br>&gt; from <a href="http://md.chem.rug.nl/education/Free-Energy_Course/3.mutation-fe.html">http://md.chem.rug.nl/education/Free-Energy_Course/3.mutation-fe.html
</a>. i am using<br>&gt; all the mdp files given on<br>&gt; the site, the confout.gro from toluene/eq directory of the previous tutorial as starting<br>&gt; structure for this simulation.<br>&gt; and whatever relevant changes has to be made in top and gro have been made. then i
<br>&gt; performed the simulation for lamda<br>&gt; values from zero to one with an interval of 0.4, so total of 26 simulations. but i am<br>&gt; not getting the mentioned<br>&gt; experimental value of ddG-hyd.<br>&gt;<br>
&gt; i got the following values in rmsd.xvg<br>&gt;<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
IN
WATER&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;IN
VACUO<br>&gt; ==========================&nbsp;&nbsp; ============================<br>&gt;
0&nbsp;&nbsp;&nbsp;&nbsp;-61.874991&nbsp;&nbsp;
2.59103&nbsp;&nbsp;&nbsp;&nbsp;0&nbsp;&nbsp;&nbsp;&nbsp;-98.521845&nbsp;&nbsp;
8.80653<br>&gt; .04&nbsp;&nbsp;-57.180709&nbsp;&nbsp; 1.73975&nbsp;&nbsp;&nbsp;&nbsp;.04&nbsp;&nbsp;-67.562344&nbsp;&nbsp; 0.237903<br>&gt; .08&nbsp;&nbsp;-58.786951&nbsp;&nbsp; 0.266783&nbsp;&nbsp; .08&nbsp;&nbsp;-64.780817&nbsp;&nbsp; 1.23656<br>&gt; .12&nbsp;&nbsp;-60.764221&nbsp;&nbsp; 0.300319&nbsp;&nbsp; .12&nbsp;&nbsp;-63.280561&nbsp;&nbsp; 1.14712<br>&gt; .16&nbsp;&nbsp;-62.700774
&nbsp;&nbsp; 0.313339&nbsp;&nbsp; .16&nbsp;&nbsp;-60.123999&nbsp;&nbsp; 0.926409<br>&gt; .20&nbsp;&nbsp;-65.107642&nbsp;&nbsp; 1.26388&nbsp;&nbsp;&nbsp;&nbsp;.20&nbsp;&nbsp;-58.584394&nbsp;&nbsp; 0.103875<br>&gt; .24&nbsp;&nbsp;-63.749578&nbsp;&nbsp; 0.384702&nbsp;&nbsp; .24&nbsp;&nbsp;-56.668703&nbsp;&nbsp; 1.02234<br>&gt; .28&nbsp;&nbsp;-64.267572&nbsp;&nbsp; 1.22184&nbsp;&nbsp;&nbsp;&nbsp;.28&nbsp;&nbsp;-54.837984&nbsp;&nbsp; 
0.0807904<br>&gt; .32&nbsp;&nbsp;-61.524424&nbsp;&nbsp; 0.402878&nbsp;&nbsp; .32&nbsp;&nbsp;-53.550764&nbsp;&nbsp; 0.0760657<br>&gt; .36&nbsp;&nbsp;-55.897929&nbsp;&nbsp; 0.314078&nbsp;&nbsp; .36&nbsp;&nbsp;-51.510805&nbsp;&nbsp; 0.0622261<br>&gt; .40&nbsp;&nbsp;-49.084890&nbsp;&nbsp; 1.05137&nbsp;&nbsp;&nbsp;&nbsp;.40&nbsp;&nbsp;-49.736461&nbsp;&nbsp; 0.0627338<br>&gt; .44&nbsp;&nbsp;-38.339801
&nbsp;&nbsp; 1.39297&nbsp;&nbsp;&nbsp;&nbsp;.44&nbsp;&nbsp;-48.470839&nbsp;&nbsp; 0.0718284<br>&gt; .48&nbsp;&nbsp;-29.613058&nbsp;&nbsp; 1.38755&nbsp;&nbsp;&nbsp;&nbsp;.48&nbsp;&nbsp;-46.842750&nbsp;&nbsp; 0.0740531<br>&gt; .52&nbsp;&nbsp;-24.891702&nbsp;&nbsp; 1.04125&nbsp;&nbsp;&nbsp;&nbsp;.52&nbsp;&nbsp;-45.007666&nbsp;&nbsp; 0.0766915<br>&gt; .56&nbsp;&nbsp;-25.435188&nbsp;&nbsp; 1.70549&nbsp;&nbsp;&nbsp;&nbsp;.56&nbsp;&nbsp;-43.939682
&nbsp;&nbsp; 0.101024<br>&gt; .60&nbsp;&nbsp;-25.180946&nbsp;&nbsp; 1.82089&nbsp;&nbsp;&nbsp;&nbsp;.60&nbsp;&nbsp;-41.844138&nbsp;&nbsp; 0.111544<br>&gt; .64&nbsp;&nbsp;-28.395296&nbsp;&nbsp; 0.849711&nbsp;&nbsp; .64&nbsp;&nbsp;-41.250909&nbsp;&nbsp; 0.13644<br>&gt; .68&nbsp;&nbsp;-29.672103&nbsp;&nbsp; 2.26016&nbsp;&nbsp;&nbsp;&nbsp;.68&nbsp;&nbsp;-41.182423&nbsp;&nbsp; 0.139366<br>&gt; .72&nbsp;&nbsp;-34.441109
&nbsp;&nbsp; 2.00941&nbsp;&nbsp;&nbsp;&nbsp;.72&nbsp;&nbsp;-41.747919&nbsp;&nbsp; 0.127424<br>&gt; .76&nbsp;&nbsp;-39.148106&nbsp;&nbsp; 2.01051&nbsp;&nbsp;&nbsp;&nbsp;.76&nbsp;&nbsp;-44.426396&nbsp;&nbsp; 0.143008<br>&gt; .80&nbsp;&nbsp;-49.360365&nbsp;&nbsp; 2.34239&nbsp;&nbsp;&nbsp;&nbsp;.80&nbsp;&nbsp;-44.381975&nbsp;&nbsp; 0.135462<br>&gt; .84&nbsp;&nbsp;-57.031557&nbsp;&nbsp; 2.79298&nbsp;&nbsp;&nbsp;&nbsp;.84&nbsp;&nbsp;-45.564523
&nbsp;&nbsp; 0.123685<br>&gt; .88&nbsp;&nbsp;-64.059494&nbsp;&nbsp; 2.79089&nbsp;&nbsp;&nbsp;&nbsp;.88&nbsp;&nbsp;-46.074452&nbsp;&nbsp; 0.104784<br>&gt; .92&nbsp;&nbsp;-78.668957&nbsp;&nbsp; 1.97957&nbsp;&nbsp;&nbsp;&nbsp;.92&nbsp;&nbsp;-46.678656&nbsp;&nbsp; 0.0993529<br>&gt; .96&nbsp;&nbsp;-85.671014&nbsp;&nbsp; 1.92346&nbsp;&nbsp;&nbsp;&nbsp;.96&nbsp;&nbsp;-46.871521&nbsp;&nbsp; 0.0923201<br>&gt; 1.0&nbsp;&nbsp;-52.956186
&nbsp;&nbsp; 1.26136&nbsp;&nbsp;&nbsp;&nbsp;1.0&nbsp;&nbsp;-47.183714&nbsp;&nbsp; 0.0934921<br>&gt;<br>&gt;<br>&gt; upon integration in xmgr, following are the dG values:<br>&gt; dG-wat = -50.6556 kJ.mol-1<br>&gt; dG-vac = -51.1109 kJ.mol-1<br>&gt;<br>&gt; ddG-hyd = 0.455
 kJ.mol-1<br>&gt;<br>&gt; the experimental value mentioned on the site is -22.1 kJ.mol-1<br>&gt;<br>&gt; Please help...<br>&gt;<br>&gt; bharat<br><br><br>_______________________________________________<br>gmx-users mailing list&nbsp;&nbsp;&nbsp;&nbsp;
<a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br><a href="http://www.gromacs.org/mailman/listinfo/gmx-users">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>Please don't post (un)subscribe requests to the list. Use the
<br>www interface or send it to <a href="mailto: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">http://www.gromacs.org/mailing_lists/users.php
</a><br></blockquote></div><br><br clear="all"><br>-- <br><br>Tsjerk A. Wassenaar, M.Sc.<br>Groningen Biomolecular Sciences and Biotechnology Institute (GBB)<br>Dept. of Biophysical Chemistry<br>University of Groningen<br>
Nijenborgh 4<br>9747AG Groningen, The Netherlands<br>+31 50 363 4336<br>