Hello all,<br><br>I'm running the same process from the free energy tutorial by Justin Lemkul...<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/gmx-tutorials/free_energy/index.html">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/gmx-tutorials/free_energy/index.html</a>.<br>
<br>How did the number of solvent particles get chosen (in the tutorial, 210 molecules were chosen)? I seem to be getting slightly different results (ranging from as small as 200 J/mol to about 5 kJ/mol depending on how many molecules I choose (ranging for example from 210 ethanol molecules to about 610 ethanol molecules for the largest energy difference change of about 5 kJ). I keep running tests to see if there is some sort of minimum atom number to get steady consistent numbers but I can't seem to find it. When I plot for example the bar.xvg & barint.xvg for both sets to see where the lines don't match up, its usually one or two points that differ slightly which cause the free energies in the end to be slightly different. I seem to be noticing too that the more atoms I use, the free energy gets a little bit lower.<br>
<br>Does anybody have any experience with this?<br><br>Thanks.<br clear="all"><br>-- <br><div><b><font face="arial,helvetica,sans-serif">Best regards,</font></b></div>
<div><b></b><font face="arial,helvetica,sans-serif"> </font></div>
<div><b><font face="arial,helvetica,sans-serif" size="4">Fabian F. Casteblanco</font></b></div>
<div><b><font face="arial,helvetica,sans-serif">Rutgers University -- </font></b></div>
<div><b><font face="arial,helvetica,sans-serif">C: +908 917 0723</font></b></div>
<div><b><font face="arial,helvetica,sans-serif">E: </font></b><a href="mailto:fabian.casteblanco@gmail.com" target="_blank"><b><font face="arial,helvetica,sans-serif">fabian.casteblanco@gmail.com</font></b></a></div><br>