docs/old-html/online/mdp.html

   1 <title>mdp file format</title>
   2 <P> See the user guide
   3 for a detailed description of the options</a>.  </P>
   4
   5 <P> Below is a sample mdp file.
   6 The ordering of the items is not important, but if you enter the same
   7 thing twice, the <b>last</b> is used (grompp gives you a note when
   8 overriding values). Dashes and underscores on the
   9 left hand side are ignored.</P>
  10
  11 <P> The values of the options are reasonable values for a 1 nanosecond
  12 MD run of a protein in a box of water. </P>
  13
  14 <hr>
  15 <pre>
  16 title                    = Yo
  17 cpp                      = /lib/cpp
  18 include                  = -I../top
  19 define                   =
  20 integrator               = md
  21 dt                       = 0.002
  22 nsteps                   = 500000
  23 nstxout                  = 5000
  24 nstvout                  = 5000
  25 nstlog                   = 5000
  26 nstenergy                = 250
  27 nstxout-compressed       = 250
  28 compressed-x-grps        = Protein
  29 energygrps               = Protein  SOL
  30 nstlist                  = 10
  31 ns-type                  = grid
  32 rlist                    = 0.8
  33 coulombtype              = cut-off
  34 rcoulomb                 = 1.4
  35 rvdw                     = 0.8
  36 tcoupl                   = Berendsen
  37 tc-grps                  = Protein      SOL
  38 tau-t                    = 0.1  0.1
  39 ref-t                    = 300  300
  40 Pcoupl                   = Berendsen
  41 tau-p                    = 1.0
  42 compressibility          = 4.5e-5
  43 ref-p                    = 1.0
  44 gen-vel                  = yes
  45 gen-temp                 = 300
  46 gen-seed                 = 173529
  47 constraints              = all-bonds
  48 </pre>
  49 <hr>
  50
  51 <p>
  52 With this input <a href="../programs/gmx-grompp.html"><tt>grompp</tt></a> will produce
  53 an <tt>mdout.mdp</tt> with all the options and descriptions:
  54 </p>
  55
  56 <hr>
  57 <pre>
  58 ; VARIOUS PREPROCESSING OPTIONS =
  59 title                    = Yo
  60 cpp                      = /lib/cpp
  61 include                  = -I../top
  62 define                   =
  63
  64 ; RUN CONTROL PARAMETERS =
  65 integrator               = md
  66 ; start time and timestep in ps =
  67 tinit                    = 0
  68 dt                       = 0.002
  69 nsteps                   = 500000
  70 ; number of steps for center of mass motion removal =
  71 nstcomm                  = 1
  72 comm-grps                =
  73
  74 ; LANGEVIN DYNAMICS OPTIONS =
  75 ; Temperature, friction coefficient (amu/ps) and random seed =
  76 bd-temp                  = 300
  77 bd-fric                  = 0
  78 ld-seed                  = 1993
  79
  80 ; ENERGY MINIMIZATION OPTIONS =
  81 ; Force tolerance and initial step-size =
  82 emtol                    = 100
  83 emstep                   = 0.01
  84 ; Max number of iterations in relax-shells =
  85 niter                    = 20
  86 ; Frequency of steepest descents steps when doing CG =
  87 nstcgsteep               = 1000
  88
  89 ; OUTPUT CONTROL OPTIONS =
  90 ; Output frequency for coords (x), velocities (v) and forces (f) =
  91 nstxout                  = 5000
  92 nstvout                  = 5000
  93 nstfout                  = 0
  94 ; Output frequency for energies to log file and energy file =
  95 nstlog                   = 5000
  96 nstenergy                = 250
  97 ; Output frequency and precision for xtc file =
  98 nstxout-compressed       = 250
  99 compressed-x-precision   = 1000
 100 ; This selects the subset of atoms for the xtc file. You can =
 101 ; select multiple groups. By default all atoms will be written. =
 102 compressed-x-grps        = Protein
 103 ; Selection of energy groups =
 104 energygrps               = Protein  SOL
 105
 106 ; NEIGHBORSEARCHING PARAMETERS =
 107 ; nblist update frequency =
 108 nstlist                  = 10
 109 ; ns algorithm (simple or grid) =
 110 ns-type                  = grid
 111 ; Periodic boundary conditions: xyz or none =
 112 pbc                      = xyz
 113 ; nblist cut-off         =
 114 rlist                    = 0.8
 115
 116 ; OPTIONS FOR ELECTROSTATICS AND VDW =
 117 ; Method for doing electrostatics =
 118 coulombtype              = cut-off
 119 rcoulomb-switch          = 0
 120 rcoulomb                 = 1.4
 121 ; Dielectric constant (DC) for cut-off or DC of reaction field =
 122 epsilon-r                = 1
 123 ; Method for doing Van der Waals =
 124 vdw-type                 = Cut-off
 125 ; cut-off lengths        =
 126 rvdw-switch              = 0
 127 rvdw                     = 0.8
 128 ; Apply long range dispersion corrections for Energy and Pressure =
 129 DispCorr                 = No
 130 ; Spacing for the PME/PPPM FFT grid =
 131 fourierspacing           = 0.12
 132 ; FFT grid size, when a value is 0 fourierspacing will be used =
 133 fourier-nx               = 0
 134 fourier-ny               = 0
 135 fourier-nz               = 0
 136 ; EWALD/PME/PPPM parameters =
 137 pme-order                = 4
 138 ewald-rtol               = 1e-05
 139 epsilon-surface          = 0
 140
 141 ; OPTIONS FOR WEAK COUPLING ALGORITHMS =
 142 ; Temperature coupling   =
 143 tcoupl                   = Berendsen
 144 ; Groups to couple separately =
 145 tc-grps                  = Protein      SOL
 146 ; Time constant (ps) and reference temperature (K) =
 147 tau-t                    = 0.1  0.1
 148 ref-t                    = 300  300
 149 ; Pressure coupling      =
 150 Pcoupl                   = Berendsen
 151 Pcoupltype               = Isotropic
 152 ; Time constant (ps), compressibility (1/bar) and reference P (bar) =
 153 tau-p                    = 1.0
 154 compressibility          = 4.5e-5
 155 ref-p                    = 1.0
 156
 157 ; SIMULATED ANNEALING CONTROL =
 158 annealing                = no
 159 ; Time at which temperature should be zero (ps) =
 160 zero-temp-time           = 0
 161
 162 ; GENERATE VELOCITIES FOR STARTUP RUN =
 163 gen-vel                  = yes
 164 gen-temp                 = 300
 165 gen-seed                 = 173529
 166
 167 ; OPTIONS FOR BONDS     =
 168 constraints              = all-bonds
 169 ; Type of constraint algorithm =
 170 constraint-algorithm     = Lincs
 171 ; Do not constrain the start configuration =
 172 unconstrained-start      = no
 173 ; Relative tolerance of shake =
 174 shake-tol                = 0.0001
 175 ; Highest order in the expansion of the constraint coupling matrix =
 176 lincs-order              = 4
 177 ; Lincs will write a warning to the stderr if in one step a bond =
 178 ; rotates over more degrees than =
 179 lincs-warnangle          = 30
 180 ; Convert harmonic bonds to morse potentials =
 181 morse                    = no
 182
 183 ; NMR refinement stuff  =
 184 ; Distance restraints type: No, Simple or Ensemble =
 185 disre                    = No
 186 ; Force weighting of pairs in one distance restraint: Equal or Conservative =
 187 disre-weighting          = Equal
 188 ; Use sqrt of the time averaged times the instantaneous violation =
 189 disre-mixed              = no
 190 disre-fc                 = 1000
 191 disre-tau                = 0
 192 ; Output frequency for pair distances to energy file =
 193 nstdisreout              = 100
 194
 195 ; Free energy control stuff =
 196 free-energy              = no
 197 init-lambda              = 0
 198 delta-lambda             = 0
 199 sc-alpha                 = 0
 200 sc-sigma                 = 0.3
 201
 202 ; Non-equilibrium MD stuff =
 203 acc-grps                 =
 204 accelerate               =
 205 freezegrps               =
 206 freezedim                =
 207 cos-acceleration         = 0
 208 energygrp-excl           =
 209
 210 ; Electric fields       =
 211 ; Format is number of terms (int) and for all terms an amplitude (real) =
 212 ; and a phase angle (real) =
 213 E-x                      =
 214 E-xt                     =
 215 E-y                      =
 216 E-yt                     =
 217 E-z                      =
 218 E-zt                     =
 219
 220 ; User defined thingies =
 221 user1-grps               =
 222 user2-grps               =
 223 userint1                 = 0
 224 userint2                 = 0
 225 userint3                 = 0
 226 userint4                 = 0
 227 userreal1                = 0
 228 userreal2                = 0
 229 userreal3                = 0
 230 userreal4                = 0
 231 </pre>