最近看了一些关于有类似膜体系的多层体系的模拟相关研究，总结了一下跑NPT系综时的.mdp的设置：
integrator               = md
dt                       = 0.002    ;ps
nsteps                   = 25000000  ;50 ns
comm-grps                = MEM SOL
nstlist                  = 10
pbc                      = xyz
periodic_molecules       = yes
rlist                    = 1.2
cutoff-scheme            = Verlet
verlet-buffer-tolerance  = 0.002
coulombtype              = PME
rcoulomb                 = 1.2
vdwtype                  = cut-off
rvdw-switch              = 1.0
rvdw                     = 1.2
DispCorr                 = EnerPres
tcoupl                   = V-rescale
tc-grps                  = system
tau-t                    = 0.1
ref-t                    = 293
pcoupl                   = C-rescale
pcoupltype               = semiisotropic
tau-p                    = 5
compressibility          = 0 4.5e-5
ref-p                    = 1.0 1.0
constraints =  hbonds