计算化学公社

标题: gromacs npt出现问题 [打印本页]

作者
Author:
哆啦哆啦    时间: 2019-1-24 21:23
标题: gromacs npt出现问题
各位大神,我的体系为聚合物和小分子,在对体系进行npt 平衡时,盒子内部出现了空洞,请问如何解决

tcoupl          = V-rescale                 ; modified Berendsen thermostat
tc-grps         = System; Protein Non-Protein   ; two coupling groups - more accurate
tau_t           = 0.05               ; time constant, in ps
ref_t           = 310.0              ; reference temperature, one for each group, in K
; Pressure coupling is off
pcoupl          = Berendsen  ;Parrinello-Rahman                 ; no pressure coupling in NVT
pcoupltype     = isotropic            ; use the function in the pressure conpling
tau_p           = 2.0                 ; time constant,in ps
ref_p           = 1.0                 ; 1 bar pressure
compressibility = 4.5e-5
; Periodic boundary conditions
pbc             = xyz               ; 3-D PBC
; Dispersion correction
DispCorr        = EnerPres      ; account for cut-off vdW scheme
; Velocity generation
gen_vel         = no    ; assign velocities from Maxwell distribution
gen_temp        = 310.0 ; temperature for Maxwell distribution
gen_seed        = -1            ; generate a random seed


作者
Author:
少年爱吃地三鲜    时间: 2019-1-25 09:39
盒子内部出现空洞是什么意思?
作者
Author:
tomwong4253    时间: 2019-1-25 10:26
盒子内部出现空洞,我认为大可能是你加水的时候水分子加的不够密集,导致水分子与水分子之间空隙太大,NPT开始的时候水分子会以彼此的作用力首先聚集。导致体系出现空洞。解决方法:1. 等足够长的时间让系统自己把空洞挤紧。2. 重新建模,在盒子中加入足够的水分子。




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3