[gmx-users] Umbrella Sampling - good histogram but no result in profile
rose.rhmn93 at gmail.com
Thu Mar 15 20:57:10 CET 2018
Check the size of pullf-umbrellaX files(in a directory which gmx WHAM use them). MAYBE you refer for example a wrong pullf-umbrella0 file unintentionally. It happened to me.
Check the size, they should be almost close to each other.
Sent from my iPhone
> On Mar 14, 2018, at 18:14, Ben Tam <btam125 at hotmail.co.uk> wrote:
> Dear all,
> I have a rather strange problem. Currently I am trying to get a PMF of a water molecule through a hydrophobic materials. At the end of using "gmx wham", I get a rather good histogram and it cover the whole range of the system. However at the making profile.xvg. I get the answer "nan" inside the file. Therefore I really cannot think of why is it giving me that?
> Currently I set it as :
> ;Pull code
> pull = yes
> ;pull-cylinder-r = 0.5
> pull_ngroups = 2
> pull_ncoords = 1
> pull_group1_name = SOL
> pull_group2_name = MOL
> pull_coord1_type = umbrella ; harmonic biasing force
> pull_coord1_geometry = direction-periodic ; simple distance increase, direction-periodic
> pull_coord1_vec = 1 1 1
> pull_coord1_groups = 1 2
> pull_coord1_dim = Y Y Y
> pull_coord1_rate = 0.00 ; 0.01 nm per ps = 10 nm per ns
> pull_coord1_k = 10000 ; kJ mol^-1 nm^-2
> pull_coord1_start = yes ; define initial COM distance > 0
> Also this appears when I used the command "gmx wham":
> Switched to exact iteration in iteration 1
> 1) Maximum change -1.000000e+20
> Converged in 2 iterations. Final maximum change -1e+20
> Therefore I really don't understand what I have to do in order to get a profile.xvg. I understand the k value is high, but without that high value, the water won't settle within the assign window.
> Am I doing something wrong? because I got a good histogram but no result in the profile.xvg? Thank you very much for your help.
> Best regards,
> Gromacs Users mailing list
> * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-request at gromacs.org.
More information about the gromacs.org_gmx-users