While using g_wham to obtain a PMF from umbrella sampling simulations I came up with the following problem. I have simulated my windows with a 0.1nm difference between structures starting from the initial position where my dimer is equilibrated and ending 0.5nm farther than what the lennard jones and coulomb interactions between these structures exist. When I include all these windows including the last 5 ones where my structures do not interact, my profile.xvg comes out as a straight line at zero energy. The histogram file looks fine, but the profile comes out wrong. When I do not include these last 5 windows my profile comes out as expected but there seems to be no final plateau. Is this normal in g_wham? Any advice on the matter would be great. Thank you
_______________________________________________ gmx-users mailing list [email protected] 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 [email protected]. Can't post? Read http://www.gromacs.org/mailing_lists/users.php
