[gmx-users] hi

Barnett, James W jbarnet4 at tulane.edu
Mon Sep 14 15:54:19 CEST 2015


On Sun, 2015-09-13 at 21:41 -0600, Idriz Seferi wrote:
> Continuation of email:
> 
> However, I still get a warning on pull-rate - which I don't
> understand.
> 
> Hope you can help me understand how to deal with the pull-rate so
> that
> grompp recognizes it. If you require the full list of parameters in
> the mdp
> file to solve the problem, then please let me know.
> Any help would be immensely appreciated.

'pull-rate1' should be 'pull-coord1-rate'

> On Sun, Sep 13, 2015 at 9:35 PM, Idriz Seferi <
> idrizseferi89 at gmail.com>
> wrote:
> 
> > Hi all,
> > 
> > I have a question regarding Umbrella Sampling in gromacs 5.1.
> > 
> > I am following the tutorial of Dr. Lemkul () and I am getting
> > errors about
> > the parameters when I am running gmx grompp.
> > Specifically, I am constantly getting an error on the pull-rate
> > option.
> > 
> > Dr. Lemkul uses the following options for the pull code:
> > 
> > ; Pull code
> > pull                    = yes
> > pull_ngroups            = 2
> > pull_ncoords            = 1
> > pull_group1_name        = Chain_A
> > pull_group2_name        = Chain_B
> > pull_coord1_type        = umbrella      ; harmonic biasing force
> > pull_coord1_geometry    = distance      ; simple distance increase
> > pull_coord1_groups      = 1 2
> > pull_coord1_dim         = N N Y
> > pull_start              = yes           ; define initial COM
> > distance > 0
> > pull_rate1              = 0.0
> > pull_k1                 = 1000          ; kJ mol^-1 nm^-2
> > 
> > However, when I run gmx grompp with these parameters I get a
> > warning on
> > pull_start, pull_rate1, and pull_k1. With the help of the gromacs
> > webpage
> > on mdp options, I modified the above list to the following:
> > 
> > pull                     = yes
> > pull_ngroups            = 2
> > pull_ncoords            = 1
> > pull_group1_name        = Chain_A
> > pull_group2_name        = Chain_B
> > pull_coord1_type        = umbrella      ; harmonic biasing force
> > pull_coord1_geometry    = distance      ; simple distance increase
> > pull_coord1_groups      = 1 2
> > pull_coord1_dim         = N N Y
> > pull_start              = yes           ; define initial COM
> > distance > 0
> > pull_coord1_rate        = 0.0          ; 0.01 nm per ps = 10 nm per
> > ns
> > pull_coord1_k           = 1000          ; kJ mol^-1 nm^-2
> > 
> > 
> > 
> > 

-- 
James “Wes” Barnett, Ph.D. Candidate
Louisiana Board of Regents Fellow

Chemical and Biomolecular Engineering
Tulane University
341-B Lindy Boggs Center for Energy and Biotechnology
6823 St. Charles Ave
New Orleans, Louisiana 70118-5674
jbarnet4 at tulane.edu


More information about the gromacs.org_gmx-users mailing list