Stefan Hoorman wrote:
Well, I tried again, this time simulating 1ns each window. I thought
about simulating more, but thought of making a test before making longer
simulations for nothing. And again, my histogram looks like a
chromatographic peak, but the range this time is different. The rest of
You should have overlapping roughly Gaussian distributions. I can't tell if
that's what you have, or if you have independent peaks. If the latter is the
case, then your windows do not overlap sufficiently to run WHAM.
the histogram is a straight line parallel to the x axis set in zero. I
Parallel lines indicate discontinuities in the WHAM analysis. I really think
you need to post an image of this online so we can see. Without seeing what
you're dealing with, it's very hard to be helpful.
didn't understand what you meant by "getting what you think you are
setting up". I have checked the distance between my groups, both with
g_dist and by visual aid. The separate all right, there is no doubt
about that. Each window is what it should be, a window for the 1.5nm
distance starts at 1.5, decreases a little bit in the first picoseconds
and then stabilizes at 1.47 or 1.45. When I start g_wham analysis, the
verbose states that my initial distance is 0.99nm and then it says that
"Determined boundaries to 1.046511 and 2.140948". I tried this time
Right, because your windows appear to range between these two values, so the
closest window is restraining your structures at roughly 1 nm, and the furthest
window is restraining the strucutres at 2.14 nm COM separation. Is this what
you expect? I thought before you were dealing with distances out to 2.5 nm.
simulating longer periods, and got the same result as before. The
profile graphic stars as a descending curve from 0 to -5 kCal, then it
starts rising in a "up the hill" way until it reaches +45kCal. I really
don't know what I am doing wrong this time.
Can you re-post the .mdp options you're using for pulling? Also, images of both
your histogram and profile would be extremely helpful. I'm having a hard time
picturing what you're describing with respect to the histograms.
-Justin
------------------------------------------------------------------------
_______________________________________________
gmx-users mailing list gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
--
========================================
Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin
========================================
_______________________________________________
gmx-users mailing list gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php