[gmx-users] pullx (gromacs 4.6 vs 2016)
gromacs query
gromacsquery at gmail.com
Fri Aug 25 12:00:14 CEST 2017
Hi All,
I am using 'equivalent' mdp settings for pull run in Gromacs 4.6 and 2016.
With gromacs 4.6 I get time, 0Z and 1dZ; which I can understand as
explained here: https://www.mail-archive.com/gmx-users@gromacs.org/
msg22453.html
With Gromacs 2016 I get only time and some distance and xvg file says is
Position. Is it same as 1dZ? Just wondering why 0Z is not printed.
#####> pull settings for 4.6<####
pull = umbrella
pull_geometry = direction ;
pull-vec1 = 0 0 -1
pull_dim = N N Y ;
pull_start = yes ; define initial COM distance > 0
pull_ngroups = 1
pull_group0 = Mem ; restrained group
pull_group1 = Peptide; group to pull
pull_init1 = 0
pull_rate1 = 0.0; nm per ps
pull_k1 = 1000 ; kJ mol^-1 nm^-2
pull-nstxout = 1000;
pull-nstfout = 1000;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
#####> pull settings for 2016<####
pull = yes
pull_coord1_type = umbrella ; equivalent to pull = umbrella
pull_coord1_geometry = direction ; equivalent to pull_geometry =
direction ;
pull-coord1-vec = 0 0 -1 ; equivalent to pull-vec1 = 0 0 -1
pull_coord1_dim = N N Y ; equivalent to pull_dim = N N Y ;
pull_coord1_start = yes; equivalent to pull_start = yes ; define
initial COM distance > 0
pull_ngroups = 2 ; equivalent to pull_ngroups = 1
pull_coord1_groups = 1 2
pull_group1_name = Mem ; equivalent to pull_group0
pull_group2_name = Peptide ; equivalent to
pull_group1
pull-coord1-init = 0 ; pull_init1 = 0
pull_coord1_rate = 0.00000 ; equivalent to pull_rate1
pull_coord1_k = 1000 ; equivalent to pull_k1 =
1000 ; kJ mol^-1 nm^-2
Thanks
JIom
More information about the gromacs.org_gmx-users
mailing list