[gmx-users] about COM pulling and g_dist
chris.neale at utoronto.ca
chris.neale at utoronto.ca
Mon Oct 4 15:54:47 CEST 2010
First, I see no reason why you would not expect your g_dist values to
be negative. If you want the positive distances, then reverse the
order of the groups as you pass then to g_dist. Also, what type of
values do you get directly out of mdrun and into the pull_x.xvg file?
Second, how did you go from values to a PMF? did you use wham? If so,
and if you gave a succession of increasing value positive numbers to
note restraint positions and then your values you passed were
increasingly negative, then you will certainly get the wrong PMF.
Third, your lack of symmetry may just be a lack of convergence. Have
you done and convergence tests?
If this doesn't solve it and you want more help, please post your
g_dist commands and post an image of your pmf online.
-- original message --
I have done a simulation about pulling a K+ throug a (9,9)CNT.
The results may be strange.
The length of the CNT is 1.477nm,the box is 5*5*5nm.I use the
Umbrella Sampling to calculate PMF ,following the tutorial by Justin
A. Lemkul. I have pulled (along Z direction) a K+ (pulling group)frome
one end (reference group ) of the CNT to another end to generate
Configurations for Umbrella Sampling,both ends of CNT were fixed. Then
I use g_dist to calculate the COM distance between K+ and the
reference group,the result is strange,the Z direction component -dz is
negative,it should be positive,because the Z coordinate of K+ is lager
than that of reference in the .xtc output file. I don't know why. the
pulling distance is not larger than half of the box size.the Z
direction component -dz is also negative in Umbrella Sampling windows.
Then do Umbrella Sampling,the pulling group and reference group is
as in COM pulling step.pull_geometry = distance,pull_rate = 0. The PMF
is negative and not symmetry about the middle of the CNT.I think it is
wrong.May be it is because the dz is negative,or the choice of
reference is not proper, it should be the middle the CNT as the
reference.Anybody could help? Thanks a lot!
; COM PULLING
; Pull type: no, umbrella, constraint or constant_force
pull = umbrella
; Pull geometry: distance, direction, cylinder or position
pull_geometry = distance
; Select components for the pull vector. default: Y Y Y
pull_dim = N N Y
; Cylinder radius for dynamic reaction force groups (nm)
pull_r1 = 1
; Switch from r1 to r0 in case of dynamic reaction force
pull_r0 = 1
pull_constr_tol = 1e-06
pull_start = yes
pull_nstxout = 500
pull_nstfout = 500
; Number of pull groups
pull_ngroups = 1
; Group name, weight (default all 1), vector, init, rate (nm/ps),
pull_group0 = refer
pull_pbcatom0 = 0
pull_group1 = K
pull_pbcatom1 = 0
pull_vec1 = 0 0 0
pull_init1 = 0
pull_rate1 = 0.0
pull_k1 = 1000
More information about the gromacs.org_gmx-users