<div dir="ltr">Hi,<div><br></div><div>I think Igor has a point that we should make our implementation of Amber as accurate as possible, but let&#39;s make sure we only change things once.</div><div><br></div><div>There has also been at least once instance where Amber changed they way they defined improper torsions internally without documentation, so we went from perfect matches in Amber ~8 to mismatching, and then back to matching.</div><div><br></div><div>1. I agree the Amber implementation should be our reference for Amber, but do the parameter files in the Amber package match the published force field exactly? Let&#39;s double-check that paper.</div><div><br></div><div>2. This If we should change, I would suggest making a quick audit of _all_ our Amber parameters and see if there are more rounding errors/deviations. If so, let&#39;s change all of them at the same time.</div><div><br></div><div>3. I also agree with Berk that this is only worth doing if we really get everything to match pretty close to perfectly. It&#39;s better to have a known insignificant discrepancy than a set of parameters that keeps changing!</div><div><br></div><div><br></div><div>Cheers,</div><div><br></div><div>Erik</div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, May 29, 2017 at 9:36 AM, Berk Hess <span dir="ltr">&lt;<a href="mailto:hess@kth.se" target="_blank">hess@kth.se</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="auto">Hi,<div dir="auto"><br></div><div dir="auto">My opinion is that we should minimize the number of time we change force field parameters, as this is confusing for users. The difference here is very small and the Coulomb constant also differs from Amber. I would only want to change this when we can get the total energy to match the Amber code up to all decimals.</div><div dir="auto"><br></div><div dir="auto">Cheers,</div><div dir="auto"><br></div><div dir="auto">Berk</div></div><br>--<br>
Gromacs Developers mailing list<br>
<br>
* Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/GMX-developers_List" rel="noreferrer" target="_blank">http://www.gromacs.org/<wbr>Support/Mailing_Lists/GMX-<wbr>developers_List</a> before posting!<br>
<br>
* Can&#39;t post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" rel="noreferrer" target="_blank">http://www.gromacs.org/<wbr>Support/Mailing_Lists</a><br>
<br>
* For (un)subscribe requests visit<br>
<a href="https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-developers" rel="noreferrer" target="_blank">https://maillist.sys.kth.se/<wbr>mailman/listinfo/gromacs.org_<wbr>gmx-developers</a> or send a mail to <a href="mailto:gmx-developers-request@gromacs.org">gmx-developers-request@<wbr>gromacs.org</a>.<br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr">--<div>Erik Lindahl &lt;<a href="mailto:erik.lindahl@gmail.com" target="_blank">erik.lindahl@gmail.com</a>&gt;</div><div>Professor of Biophysics, Dept. Biochemistry &amp; Biophysics, Stockholm University</div><div>Science for Life Laboratory, Box 1031, 17121 Solna, Sweden</div></div></div></div></div>
</div>