[gmx-users] Sampling windows
gmelaugh01 at qub.ac.uk
Wed Jan 5 18:03:01 CET 2011
Many thanks. Two more quick questions (both refering to barrier region):
Does the value of the force constant affect the height of the barrier
providing that there is reasonable overlap? Obviously there is a
compromise between the size of the constant and the number of windows.
If I increase the force constant and number of windows to obtain a high
density of sampling, does this give a more accurate curve?
XAvier Periole wrote:
> On Jan 5, 2011, at 5:01 PM, Gavin Melaugh wrote:
>> Hi Xavier
>> Many thanks for the reply. When I compare the two curves
>> Curve 1) lack of overlap in barrier region; force constant = 1000
>> Curve 2) with increased number of sampling windows with increased force
>> constant (note I have only increased the force constant of the windows
>> in the barrier region).
>> Force constant =1000, Force constant in barrier region = 5000
>> The barrier increases due to better sampling in this region (expected),
>> but the rest of the curve in 2) is shifted upward compared to curve 1).
>> Why is this so when the only difference is the increased number of
>> windows and force constant in one region?
> the relative position of the regions on the left and right of your
> region is completely depending on the barrier region!
> you may have overestimated the difference between bound/unbound states
> because the energy barrier region was not sampled.
> Note also that cure one is simply wrong so comparing to it is dangerous.
> You have to be careful that you might still not have a converged
> you might have overlap of your histograms on the reaction coordinate but
> other degrees of freedom are not fully sampled ... that might result in a
> funny and wrong profile. You might want to look at your simulations for
> funny things ...
>> XAvier Periole wrote:
>>> On Jan 5, 2011, at 4:08 PM, Gavin Melaugh wrote:
>>>> Dear all
>>>> I have generated a PMF curve for introducing a guest molecule into a
>>>> host molecule using the umbrella sampling method. The curve seems
>>>> but there are no overlaps between histograms in the barrier region. To
>>>> circumvent this I have added more windows in this region and to ensure
>>>> sampling have increased the force constant only in this region. The
>>>> windows in this problematic region are now not equally spaced. Does
>>>> matter for g_WHAM?
>>> Nope it does not matter.
>>>> Could someone please reply to me on this matter as it
>>>> has been a sticking point for me for quite sometime?
>>>> Many thanks
>>>> gmx-users mailing list gmx-users at gromacs.org
>>>> 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.
>>>> Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>> gmx-users mailing list gmx-users at gromacs.org
>> 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.
>> Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
More information about the gromacs.org_gmx-users