<html><head><style type="text/css"><!-- DIV {margin:0px;} --></style></head><body><div style="font-family:times new roman,new york,times,serif;font-size:12pt"><div>&gt; The question is:<br>&gt; Do you trust you simulations when you/others do not find significant problems<br>&gt; or do you only trust your simulations when you are sure the integrators etc<br>&gt; are implemented exactly correctly or at least the magnitude of the (small)<br>&gt; errors are known.<br>&gt; <br>&gt; Personally I (and the majority of people) fall in the first category.<br>&gt; David is clearly in the second category.</div><div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"><br>I think you can recognise free energy calculation specialists like David and myself by a heightened level of paranoia, along with other signs like staring blankly at the wall wondering which part of your protocol could be causing that 5 kcal/mol deviation... so while I'm a big
 supporter of the first category I see where David's coming from, having spent a lot of time in the past running controls for almost certainly inconsequential variables like this.<br><span>If I can offer my two eurocents-worth, I would sleep easier at night with a barostat that's both well-behaved and produces a rigorous NPT ensemble... I'm in the process of migrating to GROMACS from NAMD and have never heard any complaints about their Langevin piston thermostat (<a target="_blank" href="http://dx.doi.org/10.1063/1.470648">http://dx.doi.org/10.1063/1.470648</a>), would there be any downsides to having this method in GROMACS?</span><br>best wishes,<br><br>Floris<br><br><br><div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"><font face="Tahoma" size="2"><hr size="1"><b><span style="font-weight: bold;">From:</span></b> Berk Hess &lt;gmx3@hotmail.com&gt;<br><b><span style="font-weight: bold;">To:</span></b> Discussion list for
 GROMACS users &lt;gmx-users@gromacs.org&gt;<br><b><span style="font-weight: bold;">Sent:</span></b> Friday, 7 August, 2009 16:16:19<br><b><span style="font-weight: bold;">Subject:</span></b> RE: [gmx-users] Re: Free Energy Calculation<br></font><br>

<style>
.hmmessage P
{
margin:0px;padding:0px;}
body.hmmessage
{
font-size:10pt;font-family:Verdana;}
</style>

All the discussion is mainly a matter about how paranoid you are.<br><br>A Berendsen barostat does not generate a proper ensemble.<br>But it does get the average pressure right. The questions is how large<br>the effects of the "ensemble error" is. I would say that in practice<br>this would be negligible compared to fluctuations, force field errors,<br>insufficient sampling etc. But strictly speaking it is there.<br><br>Parrinello-Rahman in principle gives you a proper NPT ensemble<br>(combined with a proper thermostat). But to have a real symplectic<br>or even only reversible integrator is difficult and might be computationally<br>expensive, especially when constraints are present, as will nearly always<br>be the case with most water models.<br>Again errors of a non-reversible integrator are, in most cases, extremely small<br>and smaller than the other errors you make.<br>In addition PR-coupling might produce nasty volume oscillations.<br><br>The
 question is:<br>Do you trust you simulations when you/others do not find significant problems<br>or do you only trust your simulations when you are sure the integrators etc<br>are implemented exactly correctly or at least the magnitude of the (small)<br>errors are known.<br><br>Personally I (and the majority of people) fall in the first category.<br>David is clearly in the second category.<br><br>In the past decade it has been shown that several algorithms that had been<br>used a lot caused serious artifacts. But I would say that we are currently<br>in a state where most things have been checked and although we know<br>there are some errors (BTW descretizing analytic differential equations also<br>produces errors), these errors have no effect on MOST of the things we are<br>interested in. Of course you should not derive compressibilities using<br>a Berendsen barostat, as a big name in the field has done even not to long ago. <br><br>Berk<br><br><hr
 id="stopSpelling">Date: Fri, 7 Aug 2009 07:32:33 -0700<br>From: floris_buelens@yahoo.com<br>Subject: Re: [gmx-users] Re: Free Energy Calculation<br>To: gmx-users@gromacs.org<br><br>

<style>
.ExternalClass DIV
{}
</style><div style="font-family: times new roman,new york,times,serif; font-size: 12pt;">Hi David and Berk,<br><br><div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"><div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"><div class="EC_gmail_quote"><blockquote class="EC_gmail_quote" style="padding-left: 1ex;"><div>I guess the best proper solution would be to use a Parrinello-Rahman barostat.<br>But in Gromacs this is currently not implemented in a reversible way.<br></div></blockquote><div><br>&gt; Yes, I would love to see a barostat that samples the correct distribution.<br><br> Can you clarify? The wording in the manual implies that Parrinello-Rahman coupling samples the correct distribution... When you say it's not reversible do I understand right that means you
 can't run your simulation backwards and get an identical result, if yes, is this ever relevant in practice for e.g. solvation free energy calculation? David, why don't you use Parrinello-Rahman coupling?<br>thanks,<br><br>Floris<br><br><br></div><blockquote class="EC_gmail_quote" style="padding-left: 1ex;">
<div><br><br>In practice I think the issues with the Berendsen thermostat are not relevant.<br>It does not generate the correct pressure and volume fluctuations,<br>but the average pressure and volume are correct. The effects of incorrect<br>
fluctuations is negligible and is surely much smaller than simulating at the wrong<br>volume, which can have a large effect on solvation free energies.<br></div></blockquote><div><br>It seems fairly like you are right (assuming your first line is meant to refer to the Berendsen barostat), but I think no one knows for sure at this point whether it matters and when it might matter. I agree that simulating at an incorrect volume can be a big problem, however. This is unlikely, but can happen, which is why we now do the affine transformation I just mentioned when looking at solvation free energies. <br>
&nbsp;<br>Anyway -- I'd just like to see a proper barostat implemented at some point, then the paranoid (like me) can run with a proper barostat instead of worrying about whether or not this matters, AND someone could test whether observables they are interested in come out different depending on the choice of barostat.<br>
<br>David<br><br><br></div><blockquote class="EC_gmail_quote" style="padding-left: 1ex;"><div><br>I also think that issues with Berendsen pressure coupling are much less severe<br>
than the issues with Berendsen temperature coupling. Berendsen tcoupl can have<br>significant effects on your ensemble. But in Gromacs 4.0 this can easily be avoided<br>by using the v-rescale thermostat.<br><br>Berk<br><br>
<hr>Subject: RE: [gmx-users] Re: Free Energy Calculation<br>Date: Thu, 6 Aug 2009 16:01:04 -0400<br>From: <a rel="nofollow" ymailto="mailto:PNAGY@UTNet.UToledo.Edu" target="_blank" href="mailto:PNAGY@UTNet.UToledo.Edu">PNAGY@UTNet.UToledo.Edu</a><br>To: <a rel="nofollow" ymailto="mailto:gmx-users@gromacs.org" target="_blank" href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><div>
<div class="h5"><br><br>










<br>

<font size="2">A question to the answer below:<br>
<br>
&nbsp;&nbsp; If the Berendsen barostat does not sample correctly the distribution of pressures<br>
then what guarantees that the the volume/density is correct<br>
at the end of the equilibration phase when the user is advised to switch to<br>
NVT?<br>
<br>
Peter Nagy<br>
The University of Toledo<br>
<br>
-----Original Message-----<br>
From: <a rel="nofollow" ymailto="mailto:gmx-users-bounces@gromacs.org" target="_blank" href="mailto:gmx-users-bounces@gromacs.org">gmx-users-bounces@gromacs.org</a> on behalf of David Mobley<br>
Sent: Thu 8/6/2009 2:38 PM<br>
To: Mauricio Carrillo Tripp<br>
Cc: Discussion list for GROMACS users<br>
Subject: [gmx-users] Re: Free Energy Calculation<br>
<br>
Hi,<br>
<br>
Sorry for the delay answering. These questions are better put on the&nbsp;<br>
GROMACS users list.<br>
<br>
1) Yes -- the Berendsen barostat does not sample the correct&nbsp;<br>
distribution of pressures.<br>
2) Regenerating velocities is fine for a couple of reasons: (a) we are&nbsp;<br>
after thermodynamics, not dynamics, and (b) regenerating velocities&nbsp;<br>
does not de-equilibrate the system (see the Andersen thermostat, for&nbsp;<br>
example).<br>
<br>
<br>
<br>
On Jun 8, 2009, at 3:02 PM, Mauricio Carrillo Tripp wrote:<br>
<br>
&gt; Hi,<br>
&gt;<br>
&gt; I've been following directions at <a rel="nofollow" target="_blank" href="http://www.dillgroup.ucsf.edu/group/wiki/index.php/Free_Energy:_Tutorial">http://www.dillgroup.ucsf.edu/group/wiki/index.php/Free_Energy:_Tutorial</a><br>

&gt; to calculate hydration free energies of a few small molecules, and I&nbsp;<br>
&gt; have a couple of questions:<br>
&gt; 1) Is there a reason for using NVT instead of NPT for the production&nbsp;<br>
&gt; step?<br>
&gt; 2) I noticed that between the equilibration and the production runs&nbsp;<br>
&gt; you generate velocities at each step,<br>
&gt;&nbsp;&nbsp;&nbsp;&nbsp; instead of using the ones from the previous step. Is there a&nbsp;<br>
&gt; reason for doing this? wouldn't this<br>
&gt;&nbsp;&nbsp;&nbsp;&nbsp; 'de-equilibrate' an already equilibrated system?<br>
&gt;<br>
&gt; Thank you.<br>
&gt;<br>
&gt; --<br>
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp; |&nbsp; Mauricio Carrillo Tripp, PhD<br>
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; /&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; |&nbsp;&nbsp;&nbsp;&nbsp; Department of Molecular Biology, TPC6<br>
&gt;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; |&nbsp;&nbsp;&nbsp;&nbsp; The Scripps Research Institute<br>
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; \&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; |&nbsp;&nbsp;&nbsp;&nbsp; 10550 North Torrey Pines Road<br>
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp; |&nbsp;&nbsp;&nbsp;&nbsp; La Jolla, California 92037<br>
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; /&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; |&nbsp;&nbsp;&nbsp;&nbsp; <a rel="nofollow" ymailto="mailto:trippm@scripps.edu" target="_blank" href="mailto:trippm@scripps.edu">trippm@scripps.edu</a><br>
&gt;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; |&nbsp; <a rel="nofollow" target="_blank" href="http://www.scripps.edu/%7Etrippm">http://www.scripps.edu/~trippm</a><br>
&gt;<br>
&gt; ** Aut tace aut loquere meliora silentio **<br>
<br>
<br>
</font>
<br><br></div></div><hr>What can you do with the new Windows Live? <a rel="nofollow" target="_blank" href="http://www.microsoft.com/windows/windowslive/default.aspx">Find out</a></div>
<br>_______________________________________________<br>
gmx-users mailing list &nbsp; &nbsp;<a rel="nofollow" ymailto="mailto:gmx-users@gromacs.org" target="_blank" href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
<a rel="nofollow" target="_blank" href="http://lists.gromacs.org/mailman/listinfo/gmx-users">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a rel="nofollow" target="_blank" href="http://www.gromacs.org/search">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 rel="nofollow" ymailto="mailto:gmx-users-request@gromacs.org" target="_blank" href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read <a rel="nofollow" target="_blank" href="http://www.gromacs.org/mailing_lists/users.php">http://www.gromacs.org/mailing_lists/users.php</a><br></blockquote><br></div><br>
</div></div></div><br><br><hr>See all the ways you can stay connected <a rel="nofollow" target="_blank" href="http://www.microsoft.com/windows/windowslive/default.aspx">to friends and family</a></div></div></div><br>



      </body></html>