[gmx-users] Zero Potential of Mean Force with g_wham

Mark Abraham Mark.Abraham at anu.edu.au
Wed Mar 9 16:00:21 CET 2011


On 10/03/2011 1:43 AM, Susana Tomasio wrote:
> Thank you Jianguo.
>
> So can I use g_wham 4.5.2 even though I ran the simulations on 4.5.1?

Quite likely. Check the change logs on the GROMCACS website for relevant 
issues.

Mark

>
> Best regards,
>
> Susana
>
>
>
> On Wed, Mar 9, 2011 at 2:31 PM, Jianguo Li <ljggmx at yahoo.com.sg 
> <mailto:ljggmx at yahoo.com.sg>> wrote:
>
>     I met the same problem when using 4.5.1 for some systems, the PMF
>     shows zero curve, while the hist file looks fine. The problem
>     disappears when using 4.5.2.
>     Jianguo
>
>     ------------------------------------------------------------------------
>     *From:* Susana Tomasio <susietomasio at gmail.com
>     <mailto:susietomasio at gmail.com>>
>     *To:* gmx-users at gromacs.org <mailto:gmx-users at gromacs.org>
>     *Sent:* Wednesday, 9 March 2011 20:07:25
>     *Subject:* [gmx-users] Zero Potential of Mean Force with g_wham
>
>     Dear all,
>
>     I am running umbrella sampling of a molecule through a lipid
>     bilayer with gromacs 4.5.1.
>     When I extracted the potential of mean force with g_wham I got
>     zero for all the windows.
>     Any ideas of why this is happening?
>
>     This is the command I used:
>     g_wham_4.5.1 -it tpr.dat -if pullf.dat -o -hist
>
>     And this is the pull code of my .mdp file:
>
>     ; Pull code
>     pull                =  umbrella
>     pull_geometry       =  position
>     pull_dim            =  N N Y
>     pull_start          =  yes
>     pull_ngroups        =  1
>     pull_group0         =  Bilayer
>     pull_group1         =  Protein
>     pull_vec1           =  0 0 0
>     pull_init1          =  0 0 0
>     pull_rate1          =  0
>     pull_pbcatom0       =  1453
>     pull_pbcatom1       =  13187
>     pull_k1             =  3000     ;  kJ mol^-1 nm^-2
>     pull_nstxout        =  1000     ; every 2 ps
>     pull_nstfout        =  1000     ; every 2 ps
>
>     Thank you.
>     Best regards,
>
>     Susana
>
>
>     --
>     gmx-users mailing list gmx-users at gromacs.org
>     <mailto:gmx-users at gromacs.org>
>     http://lists.gromacs.org/mailman/listinfo/gmx-users
>     Please search the archive at
>     http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
>     Please don't post (un)subscribe requests to the list. Use the
>     www interface or send it to gmx-users-request at gromacs.org
>     <mailto:gmx-users-request at gromacs.org>.
>     Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20110310/c5874237/attachment.html>


More information about the gromacs.org_gmx-users mailing list