[gmx-users] Variation in free energy between GROMACS versions?
Kenneth Huang
khuang8 at student.gsu.edu
Wed Jul 24 21:33:14 CEST 2019
Hi Mark,
Thanks for the reply- I guessed a number of things had changed in the last decade, but wasn't sure if it was even salvageable as a comparison point anymore.
Out of curiosity, could you elaborate on what sort of measurements I could make to scale the difference? Since I don’t have a run with the older version, wouldn’t anything I measure reflect the current version? I imagine it can't be as simple as adjusting the cutoffs in the electrostatics.
Best,
Kenneth
________________________________
From: gromacs.org_gmx-users-bounces at maillist.sys.kth.se<mailto:gromacs.org_gmx-users-bounces at maillist.sys.kth.se> <gromacs.org_gmx-users-bounces at maillist.sys.kth.se<mailto:gromacs.org_gmx-users-bounces at maillist.sys.kth.se>> on behalf of Mark Abraham <mark.j.abraham at gmail.com<mailto:mark.j.abraham at gmail.com>>
Sent: Wednesday, July 24, 2019 1:13:51 PM
To: Discussion list for GROMACS users <gmx-users at gromacs.org<mailto:gmx-users at gromacs.org>>
Subject: Re: [gmx-users] Variation in free energy between GROMACS versions?
Hi,
On Wed, 24 Jul 2019 at 17:25, Kenneth Huang <khuang8 at student.gsu.edu<mailto:khuang8 at student.gsu.edu>> wrote:
> Hi all,
>
> I'm trying to implement a variation of the simulated tempering procedure
> as described in a 2007 paper (Choosing weights for simulated tempering,
> Sanghyun Park and Vijay S. Pande, Phys. Rev. E 76, 016703) on GROMACS 2018,
> and am trying to verify that my implementation is correct compared to the
> system in the paper (Ala10). While I've managed to match a majority of the
> parameters that are described in the paper, the energy values I'm getting
> are all shifted about 3,000 kJ/mol relative to the ones in the paper.
>
> I'm not familiar enough with how mdrun itself has changed within GROMACS,
> and don't expect to fully be able to reproduce the results of the paper, so
> I was wondering if that large of a difference is expected? Considering that
> I'm using as many of the same settings as in the paper (ie, water model,
> forcefield, number of waters, thermostat, etc), is it unusual to see that
> magnitude of a change for a relatively simple system like Ala10 between an
> older and newer version of GROMACS? Or is this something that's expected,
> and I shouldn't worry to much about it?
>
There's a couple of philosophical questions here, but the important one for
the comparison is that such an old version of GROMACS probably didn't shift
the potential to zero at the cutoff (several MD packages still don't), and
even the values of the physical constants have been updated since then
(this was also a problem in comparing single-point energies between MD
packages in a SAMPL6? paper from Shirts & co). Given that the method used
reaction-field, you could probably make some measurements on things and
scale the potential shift by the number of short-range interactions...
Mark
Best,
>
> Kenneth
> --
> Gromacs Users mailing list
>
> * Please search the archive at
> https://nam05.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.gromacs.org%2FSupport%2FMailing_Lists%2FGMX-Users_List&data=02%7C01%7Ckhuang8%40student.gsu.edu%7Cd1f0add12cea4611b07708d7105a5928%7C704d822c358a47849a1649e20b75f941%7C0%7C0%7C636995852545341419&sdata=6beERxh%2B9boxZMCtnDGqLxS8wJsZ20nXnUTCGHC4oMU%3D&reserved=0 before
> posting!
>
> * Can't post? Read https://nam05.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.gromacs.org%2FSupport%2FMailing_Lists&data=02%7C01%7Ckhuang8%40student.gsu.edu%7Cd1f0add12cea4611b07708d7105a5928%7C704d822c358a47849a1649e20b75f941%7C0%7C0%7C636995852545341419&sdata=AJrhqpp4utnDv6RsfN8BZI5AmGttzGeCX0FnDgHikT0%3D&reserved=0
>
> * For (un)subscribe requests visit
> https://nam05.safelinks.protection.outlook.com/?url=https%3A%2F%2Fmaillist.sys.kth.se%2Fmailman%2Flistinfo%2Fgromacs.org_gmx-users&data=02%7C01%7Ckhuang8%40student.gsu.edu%7Cd1f0add12cea4611b07708d7105a5928%7C704d822c358a47849a1649e20b75f941%7C0%7C0%7C636995852545341419&sdata=O2I373t%2FCNd19gqOD6M9O%2FXGqDYtvc4yG66TM1Ou3Pc%3D&reserved=0 or
> send a mail to gmx-users-request at gromacs.org<mailto:gmx-users-request at gromacs.org>.
>
--
Gromacs Users mailing list
* Please search the archive at https://nam05.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.gromacs.org%2FSupport%2FMailing_Lists%2FGMX-Users_List&data=02%7C01%7Ckhuang8%40student.gsu.edu%7Cd1f0add12cea4611b07708d7105a5928%7C704d822c358a47849a1649e20b75f941%7C0%7C0%7C636995852545341419&sdata=6beERxh%2B9boxZMCtnDGqLxS8wJsZ20nXnUTCGHC4oMU%3D&reserved=0 before posting!
* Can't post? Read https://nam05.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.gromacs.org%2FSupport%2FMailing_Lists&data=02%7C01%7Ckhuang8%40student.gsu.edu%7Cd1f0add12cea4611b07708d7105a5928%7C704d822c358a47849a1649e20b75f941%7C0%7C0%7C636995852545341419&sdata=AJrhqpp4utnDv6RsfN8BZI5AmGttzGeCX0FnDgHikT0%3D&reserved=0
* For (un)subscribe requests visit
https://nam05.safelinks.protection.outlook.com/?url=https%3A%2F%2Fmaillist.sys.kth.se%2Fmailman%2Flistinfo%2Fgromacs.org_gmx-users&data=02%7C01%7Ckhuang8%40student.gsu.edu%7Cd1f0add12cea4611b07708d7105a5928%7C704d822c358a47849a1649e20b75f941%7C0%7C0%7C636995852545351412&sdata=1%2BIZceKQjOy%2BBIWqcUet8Qg0vRS15FICa8IwfIGO%2B%2Bo%3D&reserved=0 or send a mail to gmx-users-request at gromacs.org<mailto:gmx-users-request at gromacs.org>.
More information about the gromacs.org_gmx-users
mailing list