|
| 1 | +; |
| 2 | +; File 'avatar.mdp' was generated |
| 3 | +; By user: simon (1000) |
| 4 | +; On host: laptop |
| 5 | +; At date: Thu Apr 17 11:46:50 2025 |
| 6 | +; |
| 7 | +; Created by: |
| 8 | +; :-) GROMACS - gmx grompp, 2025.1 (-: |
| 9 | +; |
| 10 | +; Executable: /home/simon/Softwares/gromacs-2025.1/build/bin/gmx |
| 11 | +; Data prefix: /home/simon/Softwares/gromacs-2025.1 (source tree) |
| 12 | +; Working dir: /home/simon/Git/GROMACS/inputs/tutorial1/avatar |
| 13 | +; Command line: |
| 14 | +; gmx grompp -f ../inputs/avatar.mdp -c conf.gro -o avatar -pp avatar -po avatar -maxwarn 1 |
| 15 | + |
| 16 | +; VARIOUS PREPROCESSING OPTIONS |
| 17 | +; Preprocessor information: use cpp syntax. |
| 18 | +; e.g.: -I/home/joe/doe -I/home/mary/roe |
| 19 | +include = |
| 20 | +; e.g.: -DPOSRES -DFLEXIBLE (note these variable names are case sensitive) |
| 21 | +define = |
| 22 | + |
| 23 | +; RUN CONTROL PARAMETERS |
| 24 | +integrator = md |
| 25 | +; Start time and timestep in ps |
| 26 | +tinit = 0 |
| 27 | +dt = 0.001 |
| 28 | +nsteps = 4000 |
| 29 | +; For exact run continuation or redoing part of a run |
| 30 | +init-step = 0 |
| 31 | +; Part index is updated automatically on checkpointing (keeps files separate) |
| 32 | +simulation-part = 1 |
| 33 | +; Multiple time-stepping |
| 34 | +mts = no |
| 35 | +; factor by which to increase the mass of the lightest atoms |
| 36 | +mass-repartition-factor = 1 |
| 37 | +; mode for center of mass motion removal |
| 38 | +comm_mode = linear |
| 39 | +; number of steps for center of mass motion removal |
| 40 | +nstcomm = 100 |
| 41 | +; group(s) for center of mass motion removal |
| 42 | +comm_grps = system |
| 43 | + |
| 44 | +; LANGEVIN DYNAMICS OPTIONS |
| 45 | +; Friction coefficient (amu/ps) and random seed |
| 46 | +bd-fric = 0 |
| 47 | +ld-seed = 48456 |
| 48 | + |
| 49 | +; ENERGY MINIMIZATION OPTIONS |
| 50 | +; Force tolerance and initial step-size |
| 51 | +emtol = 10 |
| 52 | +emstep = 0.01 |
| 53 | +; Max number of iterations in relax-shells |
| 54 | +niter = 20 |
| 55 | +; Step size (ps^2) for minimization of flexible constraints |
| 56 | +fcstep = 0 |
| 57 | +; Frequency of steepest descents steps when doing CG |
| 58 | +nstcgsteep = 1000 |
| 59 | +nbfgscorr = 10 |
| 60 | + |
| 61 | +; TEST PARTICLE INSERTION OPTIONS |
| 62 | +rtpi = 0.05 |
| 63 | + |
| 64 | +; OUTPUT CONTROL OPTIONS |
| 65 | +; Output frequency for coords (x), velocities (v) and forces (f) |
| 66 | +nstxout = 0 |
| 67 | +nstvout = 0 |
| 68 | +nstfout = 0 |
| 69 | +; Output frequency for energies to log file and energy file |
| 70 | +nstlog = 100 |
| 71 | +nstcalcenergy = 100 |
| 72 | +nstenergy = 100 |
| 73 | +; Output frequency and precision for .xtc file |
| 74 | +nstxout-compressed = 2 |
| 75 | +compressed-x-precision = 1000 |
| 76 | +; This selects the subset of atoms for the compressed |
| 77 | +; trajectory file. You can select multiple groups. By |
| 78 | +; default, all atoms will be written. |
| 79 | +compressed-x-grps = |
| 80 | +; Selection of energy groups |
| 81 | +energygrps = |
| 82 | + |
| 83 | +; NEIGHBORSEARCHING PARAMETERS |
| 84 | +; cut-off scheme (Verlet: particle based cut-offs) |
| 85 | +cutoff-scheme = Verlet |
| 86 | +; nblist update frequency |
| 87 | +nstlist = 10 |
| 88 | +; Periodic boundary conditions: xyz, no, xy |
| 89 | +pbc = xyz |
| 90 | +periodic-molecules = no |
| 91 | +; Allowed energy drift due to the Verlet buffer in kJ/mol/ps per atom, |
| 92 | +; a value of -1 means: use rlist |
| 93 | +verlet-buffer-tolerance = 0.005 |
| 94 | +; Allowed error in the average pressure due to the Verlet buffer for LJ interactions |
| 95 | +verlet-buffer-pressure-tolerance = 0.5 |
| 96 | +; nblist cut-off |
| 97 | +rlist = 1 |
| 98 | +; long-range cut-off for switched potentials |
| 99 | + |
| 100 | +; OPTIONS FOR ELECTROSTATICS AND VDW |
| 101 | +; Method for doing electrostatics |
| 102 | +coulombtype = pme |
| 103 | +coulomb-modifier = Potential-shift-Verlet |
| 104 | +rcoulomb-switch = 0 |
| 105 | +rcoulomb = 1.0 |
| 106 | +; Relative dielectric constant for the medium and the reaction field |
| 107 | +epsilon-r = 1 |
| 108 | +epsilon-rf = 0 |
| 109 | +; Method for doing Van der Waals |
| 110 | +vdw-type = Cut-off |
| 111 | +vdw-modifier = Potential-shift-Verlet |
| 112 | +; cut-off lengths |
| 113 | +rvdw-switch = 0 |
| 114 | +rvdw = 1.0 |
| 115 | +; Apply long range dispersion corrections for Energy and Pressure |
| 116 | +DispCorr = No |
| 117 | +; Extension of the potential lookup tables beyond the cut-off |
| 118 | +table-extension = 1 |
| 119 | +; Separate tables between energy group pairs |
| 120 | +energygrp-table = |
| 121 | +; Spacing for the PME/PPPM FFT grid |
| 122 | +fourierspacing = 0.1 |
| 123 | +; FFT grid size, when a value is 0 fourierspacing will be used |
| 124 | +fourier-nx = 0 |
| 125 | +fourier-ny = 0 |
| 126 | +fourier-nz = 0 |
| 127 | +; EWALD/PME/PPPM parameters |
| 128 | +pme-order = 4 |
| 129 | +ewald-rtol = 1e-05 |
| 130 | +ewald-rtol-lj = 0.001 |
| 131 | +lj-pme-comb-rule = Geometric |
| 132 | +ewald-geometry = 3d |
| 133 | +epsilon-surface = 0 |
| 134 | +implicit-solvent = no |
| 135 | + |
| 136 | +; OPTIONS FOR WEAK COUPLING ALGORITHMS |
| 137 | +ensemble-temperature-setting = auto |
| 138 | +ensemble-temperature = -1 |
| 139 | +; Temperature coupling |
| 140 | +tcoupl = v-rescale |
| 141 | +nsttcouple = -1 |
| 142 | +nh-chain-length = 10 |
| 143 | +print-nose-hoover-chain-variables = no |
| 144 | +; Groups to couple separately |
| 145 | +tc-grps = Water non-Water |
| 146 | +; Time constant (ps) and reference temperature (K) |
| 147 | +tau-t = 0.5 0.5 |
| 148 | +ref-t = 280 280 |
| 149 | +; pressure coupling |
| 150 | +pcoupl = No |
| 151 | +pcoupltype = Isotropic |
| 152 | +nstpcouple = -1 |
| 153 | +; Time constant (ps), compressibility (1/bar) and reference P (bar) |
| 154 | +tau-p = 5 |
| 155 | +compressibility = |
| 156 | +ref-p = |
| 157 | +; Scaling of reference coordinates, No, All or COM |
| 158 | +refcoord-scaling = No |
| 159 | + |
| 160 | +; OPTIONS FOR QMMM calculations |
| 161 | +QMMM = no |
| 162 | +; Groups treated with MiMiC |
| 163 | +QMMM-grps = |
| 164 | + |
| 165 | +; SIMULATED ANNEALING |
| 166 | +; Type of annealing for each temperature group (no/single/periodic) |
| 167 | +annealing = |
| 168 | +; Number of time points to use for specifying annealing in each group |
| 169 | +annealing-npoints = |
| 170 | +; List of times at the annealing points for each group |
| 171 | +annealing-time = |
| 172 | +; Temp. at each annealing point, for each group. |
| 173 | +annealing-temp = |
| 174 | + |
| 175 | +; GENERATE VELOCITIES FOR STARTUP RUN |
| 176 | +gen-vel = no |
| 177 | +gen-temp = 300 |
| 178 | +gen-seed = -675160514 |
| 179 | + |
| 180 | +; OPTIONS FOR BONDS |
| 181 | +constraints = hbonds |
| 182 | +; Type of constraint algorithm |
| 183 | +constraint-algorithm = lincs |
| 184 | +; Do not constrain the start configuration |
| 185 | +continuation = no |
| 186 | +; Use successive overrelaxation to reduce the number of shake iterations |
| 187 | +Shake-SOR = no |
| 188 | +; Relative tolerance of shake |
| 189 | +shake-tol = 0.0001 |
| 190 | +; Highest order in the expansion of the constraint coupling matrix |
| 191 | +lincs-order = 4 |
| 192 | +; Number of iterations in the final step of LINCS. 1 is fine for |
| 193 | +; normal simulations, but use 2 to conserve energy in NVE runs. |
| 194 | +; For energy minimization with constraints it should be 4 to 8. |
| 195 | +lincs-iter = 1 |
| 196 | +; Lincs will write a warning to the stderr if in one step a bond |
| 197 | +; rotates over more degrees than |
| 198 | +lincs-warnangle = 30 |
| 199 | +; Convert harmonic bonds to morse potentials |
| 200 | +morse = no |
| 201 | + |
| 202 | +; ENERGY GROUP EXCLUSIONS |
| 203 | +; Pairs of energy groups for which all non-bonded interactions are excluded |
| 204 | +energygrp-excl = |
| 205 | + |
| 206 | +; WALLS |
| 207 | +; Number of walls, type, atom types, densities and box-z scale factor for Ewald |
| 208 | +nwall = 0 |
| 209 | +wall-type = 9-3 |
| 210 | +wall-r-linpot = -1 |
| 211 | +wall-atomtype = |
| 212 | +wall-density = |
| 213 | +wall-ewald-zfac = 3 |
| 214 | + |
| 215 | +; COM PULLING |
| 216 | +pull = no |
| 217 | + |
| 218 | +; AWH biasing |
| 219 | +awh = no |
| 220 | + |
| 221 | +; ENFORCED ROTATION |
| 222 | +; Enforced rotation: No or Yes |
| 223 | +rotation = no |
| 224 | + |
| 225 | +; Group to display and/or manipulate in interactive MD session |
| 226 | +IMD-group = |
| 227 | + |
| 228 | +; NMR refinement stuff |
| 229 | +; Distance restraints type: No, Simple or Ensemble |
| 230 | +disre = No |
| 231 | +; Force weighting of pairs in one distance restraint: Conservative or Equal |
| 232 | +disre-weighting = Conservative |
| 233 | +; Use sqrt of the time averaged times the instantaneous violation |
| 234 | +disre-mixed = no |
| 235 | +disre-fc = 1000 |
| 236 | +disre-tau = 0 |
| 237 | +; Output frequency for pair distances to energy file |
| 238 | +nstdisreout = 100 |
| 239 | +; Orientation restraints: No or Yes |
| 240 | +orire = no |
| 241 | +; Orientation restraints force constant and tau for time averaging |
| 242 | +orire-fc = 0 |
| 243 | +orire-tau = 0 |
| 244 | +orire-fitgrp = |
| 245 | +; Output frequency for trace(SD) and S to energy file |
| 246 | +nstorireout = 100 |
| 247 | + |
| 248 | +; Free energy variables |
| 249 | +free-energy = no |
| 250 | +couple-moltype = |
| 251 | +couple-lambda0 = vdw-q |
| 252 | +couple-lambda1 = vdw-q |
| 253 | +couple-intramol = no |
| 254 | +init-lambda = -1 |
| 255 | +init-lambda-state = -1 |
| 256 | +delta-lambda = 0 |
| 257 | +nstdhdl = 50 |
| 258 | +fep-lambdas = |
| 259 | +mass-lambdas = |
| 260 | +coul-lambdas = |
| 261 | +vdw-lambdas = |
| 262 | +bonded-lambdas = |
| 263 | +restraint-lambdas = |
| 264 | +temperature-lambdas = |
| 265 | +calc-lambda-neighbors = 1 |
| 266 | +init-lambda-weights = |
| 267 | +init-lambda-counts = |
| 268 | +init-wl-histogram-counts = |
| 269 | +dhdl-print-energy = no |
| 270 | +sc-function = beutler |
| 271 | +sc-alpha = 0 |
| 272 | +sc-power = 1 |
| 273 | +sc-r-power = 6 |
| 274 | +sc-sigma = 0.3 |
| 275 | +sc-coul = no |
| 276 | +sc-gapsys-scale-linpoint-lj = 0.85 |
| 277 | +sc-gapsys-scale-linpoint-q = 0.3 |
| 278 | +sc-gapsys-sigma-lj = 0.3 |
| 279 | +separate-dhdl-file = yes |
| 280 | +dhdl-derivatives = yes |
| 281 | +dh_hist_size = 0 |
| 282 | +dh_hist_spacing = 0.1 |
| 283 | + |
| 284 | +; Non-equilibrium MD stuff |
| 285 | +acc-grps = |
| 286 | +accelerate = |
| 287 | +freezegrps = |
| 288 | +freezedim = |
| 289 | +cos-acceleration = 0 |
| 290 | +deform = |
| 291 | +deform-init-flow = no |
| 292 | + |
| 293 | +; simulated tempering variables |
| 294 | +simulated-tempering = no |
| 295 | +simulated-tempering-scaling = geometric |
| 296 | +sim-temp-low = 300 |
| 297 | +sim-temp-high = 300 |
| 298 | + |
| 299 | +; Ion/water position swapping for computational electrophysiology setups |
| 300 | +; Swap positions along direction: no, X, Y, Z |
| 301 | +swapcoords = no |
| 302 | +adress = no |
| 303 | + |
| 304 | +; User defined thingies |
| 305 | +user1-grps = |
| 306 | +user2-grps = |
| 307 | +userint1 = 0 |
| 308 | +userint2 = 0 |
| 309 | +userint3 = 0 |
| 310 | +userint4 = 0 |
| 311 | +userreal1 = 0 |
| 312 | +userreal2 = 0 |
| 313 | +userreal3 = 0 |
| 314 | +userreal4 = 0 |
| 315 | +; Electric fields |
| 316 | +; Format for electric-field-x, etc. is: four real variables: |
| 317 | +; amplitude (V/nm), frequency omega (1/ps), time for the pulse peak (ps), |
| 318 | +; and sigma (ps) width of the pulse. Omega = 0 means static field, |
| 319 | +; sigma = 0 means no pulse, leaving the field to be a cosine function. |
| 320 | +electric-field-x = 0 0 0 0 |
| 321 | +electric-field-y = 0 0 0 0 |
| 322 | +electric-field-z = 0 0 0 0 |
| 323 | + |
| 324 | +; Density guided simulation |
| 325 | +density-guided-simulation-active = false |
| 326 | + |
| 327 | +; QM/MM with CP2K |
| 328 | +qmmm-cp2k-active = false |
| 329 | + |
| 330 | +; Colvars bias |
| 331 | +colvars-active = false |
| 332 | + |
| 333 | +; Neural Network potential |
| 334 | +nnpot-active = false |
0 commit comments