<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hello,<div><br><div><div>On 16.02.2010, at 20:25, Gabriel Rocklin wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><span class="Apple-style-span" style="font-family: arial, sans-serif; font-size: 13px; border-collapse: collapse; "><div>Hi all,</div><div><br></div><div>I'm seeing some confusing behavior with Gromacs 3.3.4 and PME electrostatics. &nbsp;My potential energies seem to not to converge a single "correct" value as I make the PME parameters more and more exact.</div>

<div><br></div></span></blockquote><div><br></div><div>How can you be sure, that the PME parameters are more exact ? Of course if you increase your cutoff in real space, the real space part of the electrostatic interaction is getting more accurate. However changing rcut has an influence on the splitting parameter of the Ewald summation. So increasing the cutoff while keeping rtol constant shifts the splitting paramter to a smaller value, that is perhaps not suitable for a precise calculation of the electrostatics.&nbsp;</div><div><br></div><div>To answer you final questions. The energy will converge if you have an appropriate set of parameters.&nbsp;</div><div><br></div><div>Cheers,</div><div><br></div><div>Flo</div><div><br></div><div><br></div><br><blockquote type="cite"><span class="Apple-style-span" style="font-family: arial, sans-serif; font-size: 13px; border-collapse: collapse; "><div>I took 10 snapshots from a single precision, constant volume simulation I'm doing of a protein in water and saved them as .gro files with -ndec = 8. &nbsp;I then used mdrun -rerun to re-evaluate their energies using a number of different PME parameters. &nbsp;I found that the energies were not sensitive to ewald_rtol or fourierspacing (below certain values), but they were very sensitive to rlist and rcoulomb even out to very long distances.</div>

<div><br></div><div>For example, using pme_order=6, spacing=0.1, and tolerance=0.00000001, as I increase rlist and rcoulomb together from 0.7 to 2.5nm, I get these potential energies for my 10 snapshots:</div><div><br></div>

<div>0.7<span style="white-space: pre; ">        </span>-656927.375<span style="white-space: pre; ">        </span>-657225.4375<span style="white-space: pre; ">        </span>-656475.0<span style="white-space: pre; ">        </span>-656451.875<span style="white-space: pre; ">        </span>-657604.8125<span style="white-space: pre; ">        </span>-656090.6875<span style="white-space: pre; ">        </span>-657553.25<span style="white-space: pre; ">        </span>-657241.8125<span style="white-space: pre; ">        </span>-657747.9375<span style="white-space: pre; ">        </span>-656863.5625</div>

<div>0.9<span style="white-space: pre; ">        </span>-654605.375<span style="white-space: pre; ">        </span>-654916.1875<span style="white-space: pre; ">        </span>-654163.9375<span style="white-space: pre; ">        </span>-654132.125<span style="white-space: pre; ">        </span>-655300.875<span style="white-space: pre; ">        </span>-653771.5<span style="white-space: pre; ">        </span>-655239.875<span style="white-space: pre; ">        </span>-654936.875<span style="white-space: pre; ">        </span>-655420.75<span style="white-space: pre; ">        </span>-654539.25</div>

<div>1.0<span style="white-space: pre; ">        </span>-654553.5625<span style="white-space: pre; ">        </span>-654862.0625<span style="white-space: pre; ">        </span>-654108.9375<span style="white-space: pre; ">        </span>-654078.0<span style="white-space: pre; ">        </span>-655246.875<span style="white-space: pre; ">        </span>-653718.4375<span style="white-space: pre; ">        </span>-655182.6875<span style="white-space: pre; ">        </span>-654884.3125<span style="white-space: pre; ">        </span>-655365.8125<span style="white-space: pre; ">        </span>-654486.375</div>

<div>1.1<span style="white-space: pre; ">        </span>-654536.875<span style="white-space: pre; ">        </span>-654842.5<span style="white-space: pre; ">        </span>-654094.3125<span style="white-space: pre; ">        </span>-654061.25<span style="white-space: pre; ">        </span>-655227.9375<span style="white-space: pre; ">        </span>-653701.4375<span style="white-space: pre; ">        </span>-655166.5<span style="white-space: pre; ">        </span>-654864.3125<span style="white-space: pre; ">        </span>-655349.625<span style="white-space: pre; ">        </span>-654468.75</div>

<div>1.2<span style="white-space: pre; ">        </span>-654533.5625<span style="white-space: pre; ">        </span>-654842.625<span style="white-space: pre; ">        </span>-654087.1875<span style="white-space: pre; ">        </span>-654054.5<span style="white-space: pre; ">        </span>-655224.625<span style="white-space: pre; ">        </span>-653695.125<span style="white-space: pre; ">        </span>-655164.25<span style="white-space: pre; ">        </span>-654863.5625<span style="white-space: pre; ">        </span>-655340.5<span style="white-space: pre; ">        </span>-654463.125</div>

<div>1.3<span style="white-space: pre; ">        </span>-654531.9375<span style="white-space: pre; ">        </span>-654837.0625<span style="white-space: pre; ">        </span>-654087.5625<span style="white-space: pre; ">        </span>-654054.5625<span style="white-space: pre; ">        </span>-655223.6875<span style="white-space: pre; ">        </span>-653694.125<span style="white-space: pre; ">        </span>-655161.6875<span style="white-space: pre; ">        </span>-654857.5<span style="white-space: pre; ">        </span>-655345.0625<span style="white-space: pre; ">        </span>-654462.5</div>

<div>1.5<span style="white-space: pre; ">        </span>-654532.375<span style="white-space: pre; ">        </span>-654839.8125<span style="white-space: pre; ">        </span>-654086.1875<span style="white-space: pre; ">        </span>-654057.25<span style="white-space: pre; ">        </span>-655223.25<span style="white-space: pre; ">        </span>-653693.3125<span style="white-space: pre; ">        </span>-655164.4375<span style="white-space: pre; ">        </span>-654864.1875<span style="white-space: pre; ">        </span>-655340.1875<span style="white-space: pre; ">        </span>-654461.5625</div>

<div>1.8<span style="white-space: pre; ">        </span>-654527.125<span style="white-space: pre; ">        </span>-654840.3125<span style="white-space: pre; ">        </span>-654086.9375<span style="white-space: pre; ">        </span>-654054.0625<span style="white-space: pre; ">        </span>-655221.25<span style="white-space: pre; ">        </span>-653696.6875<span style="white-space: pre; ">        </span>-655155.75<span style="white-space: pre; ">        </span>-654860.0625<span style="white-space: pre; ">        </span>-655343.5<span style="white-space: pre; ">        </span>-654465.8125</div>

<div>2.0<span style="white-space: pre; ">        </span>-654532.5625<span style="white-space: pre; ">        </span>-654843.125<span style="white-space: pre; ">        </span>-654083.75<span style="white-space: pre; ">        </span>-654053.6875<span style="white-space: pre; ">        </span>-655223.375<span style="white-space: pre; ">        </span>-653695.375<span style="white-space: pre; ">        </span>-655158.8125<span style="white-space: pre; ">        </span>-654860.1875<span style="white-space: pre; ">        </span>-655342.3125<span style="white-space: pre; ">        </span>-654464.4375</div>

<div>2.5<span style="white-space: pre; ">        </span>-654535.875<span style="white-space: pre; ">        </span>-654838.75<span style="white-space: pre; ">        </span>-654089.875<span style="white-space: pre; ">        </span>-654055.0625<span style="white-space: pre; ">        </span>-655220.25<span style="white-space: pre; ">        </span>-653692.75<span style="white-space: pre; ">        </span>-655163.4375<span style="white-space: pre; ">        </span>-654857.375<span style="white-space: pre; ">        </span>-655337.125<span style="white-space: pre; ">        </span>-654461.375</div>

<div><br></div><div>So not only do the numbers keep changing, but the relative energies change as well, so that for example going from 2.0 to 2.5 nm, the standard deviation of the change in energy across the 10 snapshots is 4 kJ.</div>

<div><br></div><div>The system is a dodecahedral box with a volume of about 500 nm^3.</div><div><br></div><div>Is this normal behavior? &nbsp;Should these values converge to specific energies? &nbsp;And how far out would you normally have to set rcoulomb/rlist to see this?</div>

<div><br></div><div>Thanks,</div><div>Gabriel</div><div><br></div><div>Gabriel Rocklin<br>UCSF Shoichet/Dill Labs</div></span>
-- <br>gmx-developers mailing list<br><a href="mailto:gmx-developers@gromacs.org">gmx-developers@gromacs.org</a><br><a href="http://lists.gromacs.org/mailman/listinfo/gmx-developers">http://lists.gromacs.org/mailman/listinfo/gmx-developers</a><br>Please don't post (un)subscribe requests to the list. Use the <br>www interface or send it to gmx-developers-request@gromacs.org.</blockquote></div><br><div>
<span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div>--</div><div>Florian Dommert</div><div>Dipl.-Phys.<br><br>Institute for Computational Physics<br><br>University Stuttgart<br><br>Pfaffenwaldring 27<br>70569 Stuttgart<br><br>Phone: +49(0)711/685-6-3613<br>Fax: &nbsp;&nbsp;+49-(0)711/685-6-3658&nbsp;<br><br>EMail:&nbsp;<a href="mailto:dommert@fias.uni-frankfurt.de">dommert@icp.uni-stuttgart.de</a><br>Home:&nbsp;<a href="http://fias.uni-frankfurt.de/~simbio/Florian_Dommert">http://www.icp.uni-stuttgart.de/~icp/Florian_Dommert</a></div></div></div></div></div></div></div></span></div></span></div></span></div></span></span>
</div>
<br></div></body></html>