sobereva 发表于 2023-2-11 20:28
有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，这点在置顶的新社员 ...

好的老师，我下次会注意社内发帖说明

各向异性控压，Z方向设很大压力

1416339770 发表于 2023-2-11 10:05
这是平衡用eq.mdp文件
define =
integrator = md

有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，这点在置顶的新社员必读贴里明确说了。

这是平衡用eq.mdp文件
define =
integrator = md
dt         = 0.002  ; ps
nsteps     = 500000 ; 1000ps
comm-mode  = Linear
comm-grps  = drop1 drop2 AU
energygrps =
;
nstxout = 1000
nstvout = 1000
nstfout = 1000
nstlog  = 1000
nstenergy = 500
nstxout-compressed = 500
compressed-x-grps  = system
;
annealing = single single single
annealing_npoints = 2 2 2
annealing_time = 0 600 0 600 0 600  ;ps
annealing_temp = 0 298.15 0 298.15 0 298.15
;
pbc = xyz
cutoff-scheme = Verlet
coulombtype   = PME
rcoulomb      = 1.0
vdwtype       = cut-off
rvdw          = 1.0
;
Tcoupl  = V-rescale
tau_t   = 0.2 0.2 0.2
tc_grps = SOL _ICZ AU
ref_t   = 298.15 298.15 298.15
;
freezegrps  =
freezedim   =
constraints = hbonds
;
acc-grps = top bottom
accelerate = 0 -0.00004236 0 0 0.00004236 0