Hi, I am trying to calculate the PMF of two peptides using constraint in gromacs 4.0.5. But, I am finding a very large oscillation( from -1000 to +1000) in the force that is being printed in pullf.xvg file . Is this oscillation normal ? In that case at any distance, the mean force will be zero ( I assume).
Here is the excerpt of the mdp file: ; pull code pull = constraint pull_geometry = distance pull_dim = Y Y Y pull_init1 = 1.5 pull_k1 = pull_group0 = CA1 pull_ngroups = 1 pull_group1 = CA2 pull_nstxout = 50 ;step pull_nstfout = 20 Thanks Jagannath
-- 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