[gmx-users] Re: question about minimisation

reisingere at rostlab.informatik.tu-muenchen.de reisingere at rostlab.informatik.tu-muenchen.de
Thu Jul 26 10:47:25 CEST 2012


Ho,
first I minimize my structure. This is the corresponding mdp file:

define                  = -DPOSRES
integrator              = steep
emtol           = 10
nsteps          = 1500
nstenergy               = 1
energygrps              = System
coulombtype             = PME
rcoulomb                = 0.9
rvdw                    = 0.9
rlist           = 0.9
fourierspacing  = 0.12
pme_order               = 4
ewald_rtol              = 1e-5
pbc                     = xyz


and then I run a md run. This is the corresponding mdp file:

define          = -DPOSRES
integrator              = md
dt                      = 0.001
nsteps          = 5000
nstxout         = 100
nstvout         = 0
nstfout         = 0
nstlog          = 1000
nstxtcout               = 500
nstenergy               = 5
energygrps              = Protein Non-Protein
nstcalcenergy   = 5
nstlist         = 10
ns-type         = Grid
pbc                     = xyz
rlist           = 0.9
coulombtype             = PME
rcoulomb                = 0.9
rvdw            = 0.9
fourierspacing          = 0.12
pme_order               = 4
ewald_rtol              = 1e-5
gen_vel         = yes
gen_temp                = 200.0
gen_seed                = 9999
constraints             = all-bonds
tcoupl          = V-rescale
tc-grps         = Protein  Non-Protein
tau_t           = 0.1      0.1
ref_t           = 298      298
pcoupl          = no



In my topology file I include the restraint files like this:

; Include Position restraint file
#ifdef POSRES
#include "posre.itp"
#endif

#ifdef POSRES
#include "posre_memb.itp"
#endif


I just recognize that there is a "DPOSRES" in the mdp files and a "POSRES"
in my topology file. Is this the problem. Do I have to write it the same
way in the several files?

Thank you


> Hi,
>
> Did you call your genrestr output in your mdp file using define = * ? (for
> example define = -DPOSRES ; where posres was your file name.)
>
> Thanks
>
> Peterson J
>
>
>
> --
> View this message in context:
> http://gromacs.5086.n6.nabble.com/question-about-minimisation-tp4999681p4999685.html
> Sent from the GROMACS Users Forum mailing list archive at Nabble.com.
> --
> gmx-users mailing list    gmx-users at gromacs.org
> http://lists.gromacs.org/mailman/listinfo/gmx-users
> * Only plain text messages are allowed!
> * Please search the archive at
> http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
> * Please don't post (un)subscribe requests to the list. Use the
> www interface or send it to gmx-users-request at gromacs.org.
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>





More information about the gromacs.org_gmx-users mailing list