Graphite 发表于 2021-10-29 11:31
建议把初始构型中两个东西放近一点，因为长程截断设置的12 A，从图片看两者距离可能已经（部分）超过12 A了 ...

谢谢您的建议，我想请教一下PME不是除了ntb=0的情况下就默认打开的么？还需要写其他的控制条件么？我之前也尝试过摆的比较近，但是配体链像是质心固定了，只有构象的变化没有吸附过程。我拿同样的构型用gromacs模拟就可以进行吸附。不知道这之间有什么问题。
以下是gromacs的产出文件。


;VARIOUS PREPROCESSING OPYIONS
title         = nvt_prod
define        = -DPOSRES

;RUN CONTROL PARAMETERS
integrator    = md
dt            = 0.002
nsteps        = 50000000

;OUTPUT CONTROL OPTIONS
nstxout       = 5000
nstvout       = 5000
nstfout       = 5000
nstenergy     = 5000
nstlog        = 5000
energygrps    = SIW TEN SOL ION

;NEIGHBORSEARCHING PARAMETERS
nstlist       = 5
ns-type       = grid
pbc           = xyz
rlist         = 1.2

;options for electrostatics and vdw
coulombtype= PME
pme_order= 4
fourierspacing = 0.12
rcoulomb       = 1.2
vdw-type       = cut-off
rvdw           = 1.2
cutoff_scheme  = verlet
;dispersion coorection
DispCorr      = EnerPres

;pressure coupling
Pcoupl        = no

;generate velocities for startup run
gen_vel= yes
gen-temp = 300
gen-seed = -1

;fix
freezegrps = SIW
Freezedim = Y Y Y

;temperature coupling
Tcoupl        = berendsen
tc-grps       = SIW TEN SOL ION
tau_t         = 0.5 0.5 0.5 0.5
ref_t         = 300 300 300 300

建议把初始构型中两个东西放近一点，因为长程截断设置的12 A，从图片看两者距离可能已经（部分）超过12 A了，或者使用PME之类的长程作用处理方法。

sobereva 发表于 2021-10-29 06:21
直接发个结构图，一目了然，免得别人还得载进去看
上传多个文件的时候始终说清楚每个文件都对应什么，别让 ...

结构和文件说明已经添加了，请sob老师看一下，谢谢。

直接发个结构图，一目了然，免得别人还得载进去看
上传多个文件的时候始终说清楚每个文件都对应什么，别让别人打开后才能判断