NPT_short_eq_STRONG.mdp

title                    = Martini
integrator               = md
define                   = -DSTRONG
dt                       = 0.005
nsteps                   = 1000 ;2000000 ;50 ns
nstcomm                  = 100
comm-grps		         = 

nstxtcout 		 = 1000 ;50 ps
nstxout                  = 0
nstvout                  = 0
nstfout                  = 0
nstlog                   = 1000 ;50 ps
nstenergy                = 0

cutoff-scheme            = Verlet
nstlist                  = 30
ns_type                  = grid
pbc                      = xyz
verlet-buffer-tolerance  = 0.005

coulombtype              = reaction-field 
rcoulomb                 = 1.1
epsilon_r                = 15	; 2.5 (with polarizable water)
epsilon_rf               = 0
vdw_type                 = cutoff  
vdw-modifier             = Potential-shift-verlet
rvdw                     = 1.1

tcoupl                   = v-rescale 
tc-grps                  = Protein DPPC POPC WN ION
tau_t                    = 1.0  1.0 1.0 1.0 1.0
ref_t                    = 312 312 312 312 312
Pcoupl                   = Berendsen ;parrinello-rahman 
Pcoupltype               = semiisotropic
tau_p                    = 12.0 ;parrinello-rahman is more stable with larger tau-p, DdJ, 20130422
compressibility          = 3e-4  3e-4
ref_p                    = 1.0  1.0

gen_vel                  = yes
gen_temp                 = 312
gen_seed                 = 473529

refcoord_scaling	 = com

constraints              = none 
constraint_algorithm     = Lincs
lincs-order=8
lincs-iter=2