Dear Qian:

Debiasing equations for this can be found by analogy to eq. 16 and 17 in the the original Torrie and Valleau US paper: "Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling".

Alternatively, I used a simple method in doi:10.1016/j.cplett.2008.05.099 which was to run 2-dimensional WHAM with my quantity of interest as the second dimension, but assigning this second dimension a force constant of zero.

Be careful with either method, though, as apparent convergence of the free energy does not imply convergence of all other observables.

Chris.

-- original message --

Hi,

I am using Gromacs 4.0.5 for umbrella sampling. I know that I can use g_wham to calculate the free energy along the constrained distance. Is there a command in Gromacs 4.0.5 that can calculates the average of any other obserble (for example, potential energy) along the constrained distance? Thanks a lot.

Sincerely,
Qian


--
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

Reply via email to