comparison test-data/nvt.mdp @ 0:644ca8dfdf28 draft

planemo upload for repository https://github.com/galaxycomputationalchemistry/galaxy-tools-compchem/tools/gromacs commit 177cab098ac548bfdecba9d9f04614aec5f6c618
author chemteam
date Thu, 04 Oct 2018 17:39:10 -0400
parents
children
comparison
equal deleted inserted replaced
-1:000000000000 0:644ca8dfdf28
1 title = OPLS Lysozyme NVT equilibration
2 define = -DPOSRES ; position restrain the protein
3 ; Run parameters
4 integrator = md ; leap-frog integrator
5 nsteps = 500 ; 2 * 50000 = 100 ps
6 dt = 0.002 ; 2 fs
7 ; Output control
8 nstxout = 50 ; save coordinates every 1.0 ps
9 nstvout = 50 ; save velocities every 1.0 ps
10 nstenergy = 50 ; save energies every 1.0 ps
11 nstlog = 50 ; update log file every 1.0 ps
12 nstxout-compressed = 50 ; save compressed coordinates every 10.0 ps
13 ; Bond parameters
14 continuation = no ; first dynamics run
15 constraint_algorithm = lincs ; holonomic constraints
16 constraints = all-bonds ; all bonds (even heavy atom-H bonds) constrained
17 lincs_iter = 1 ; accuracy of LINCS
18 lincs_order = 4 ; also related to accuracy
19 ; Neighborsearching
20 cutoff-scheme = Verlet
21 ns_type = grid ; search neighboring grid cells
22 nstlist = 10 ; 20 fs, largely irrelevant with Verlet
23 rcoulomb = 1.0 ; short-range electrostatic cutoff (in nm)
24 rvdw = 1.0 ; short-range van der Waals cutoff (in nm)
25 ; Electrostatics
26 coulombtype = PME ; Particle Mesh Ewald for long-range electrostatics
27 pme_order = 4 ; cubic interpolation
28 fourierspacing = 0.16 ; grid spacing for FFT
29 ; Temperature coupling is on
30 tcoupl = V-rescale ; modified Berendsen thermostat
31 tc-grps = Protein Non-Protein ; two coupling groups - more accurate
32 tau_t = 0.1 0.1 ; time constant, in ps
33 ref_t = 300 300 ; reference temperature, one for each group, in K
34 ; Pressure coupling is off
35 pcoupl = no ; no pressure coupling in NVT
36 ; Periodic boundary conditions
37 pbc = xyz ; 3-D PBC
38 ; Dispersion correction
39 DispCorr = EnerPres ; account for cut-off vdW scheme
40 ; Velocity generation
41 gen_vel = yes ; assign velocities from Maxwell distribution
42 gen_temp = 300 ; temperature for Maxwell distribution
43 gen_seed = -1 ; generate a random seed