24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 564  |  回复: 1

tonylee33

新虫 (初入文坛)

[求助] 隐式水溶液不能用NPT吗 已有1人参与

做了个模拟,隐式水,加了个压强的coupling,用的是MTTK的方法。但是模拟跑出来的结果,盒子不停的缩小,最后体系就崩溃了。求高手给点意见?是不是说隐式水本身就不能用NPT模拟呢?

先谢谢诸位。

; OPTIONS FOR ELECTROSTATICS AND VDW
; Method for doing electrostatics
coulombtype              = cut-off
rcoulomb-switch          = 0.8
rcoulomb                 = 1.2
; Relative dielectric constant for the medium and the reaction field
epsilon_r                = 2
epsilon_rf               = 0
; Method for doing Van der Waals
vdw-type                 = cut-off
; cut-off lengths      
rvdw-switch              = 0.8
rvdw                     = 1.2
; SIMULATED ANNEALING  
; Type of annealing for each temperature group (no/single/periodic)
annealing                = single
; Number of time points to use for specifying annealing in each group
annealing-npoints        = 4
; List of times at the annealing points for each group
annealing-time           = 0 10 20 30
; Temp. at each annealing point, for each group.
annealing-temp           = 10 100 200 300
; Temperature coupling is on
tcoupl                   = nose-hoover  ; nose-hoover thermostat
tc-grps                  = System             ; coupling groups - more accurate
tau_t                     = 5                          ; time constant, in ps
ref_t                      = 300                     ; reference temperature, one for each group, in K
pcoupl                  = MTTK                 ; Pressure coupling on in NPT
pcoupltype         = isotropic           ; uniform scaling of box vectors
tau_p                    = 5                          ; time constant, in ps
ref_p                     = 1.0                      ; reference pressure, in bar
compressibility = 4.5e-5                ; isothermal compressibility of water, bar^-1
refcoord_scaling = com
; Periodic boundary conditions
pbc                         = xyz                      ; 3-D PBC
; Dispersion correction
DispCorr               = EnerPres          ; account for cut-off vdW scheme
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

宅得要死
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsloneil

专家顾问 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
implicit solvent一般不用NPT。假设压强为1 atm,这其中很大一部分是由水贡献的,implicit solvent在计算压力时没法包含水的这一部分贡献,软件所计算的压力仅仅是溶质的部分。

implicit solvent下一般用NVT。
2楼2015-03-13 11:21:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tonylee33 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 总分343,求生物学调剂 +6 深刻的凯撒 2026-03-26 6/300 2026-04-01 15:58 by 2007huanle
[考研] 321求调剂 +8 y-yh 2026-04-01 8/400 2026-04-01 14:44 by 星空星月
[考研] 一志愿中国科学院大学265求调剂 +8 恬淡ye 2026-03-31 9/450 2026-04-01 14:34 by 逆水乘风
[考研] 生物学296求调剂 +10 汤圆包 2026-03-29 14/700 2026-04-01 10:44 by 求调剂zz
[考研] 0856材料化工调剂 总分330 +18 zhubinhao 2026-03-27 18/900 2026-04-01 09:37 by oooqiao
[考研] 085900土木水利336分求调剂 +3 Zhangjiangj 2026-03-31 5/250 2026-04-01 09:14 by Zhangjiangj
[考研] 求调剂 生物学 377分 +6 zzll03 2026-03-31 6/300 2026-03-31 17:33 by 唐沐儿
[考研] 334求调剂 +7 Trying] 2026-03-31 7/350 2026-03-31 12:33 by 无际的草原
[考研] 266分,求材料相关专业调剂 +10 哇呼哼呼哼 2026-03-30 12/600 2026-03-31 11:00 by 熊一刀
[考研] 293分求调剂,外语为俄语 +5 加一一九 2026-03-31 5/250 2026-03-31 09:39 by zhshch
[考研] 吉大生物学326分求调剂 +3 sunnyupup 2026-03-31 3/150 2026-03-31 09:28 by longlotian
[考研] 285求调剂 +6 AZMK 2026-03-29 9/450 2026-03-30 21:02 by dophin1985
[考研] 328求调剂 +8 嗯滴的基本都 2026-03-27 8/400 2026-03-30 17:20 by Wang200018
[考研] 332求调剂 +6 @MZB382400 2026-03-28 6/300 2026-03-30 16:57 by 无际的草原
[考研] 303求调剂 +7 DLkz1314. 2026-03-30 7/350 2026-03-30 16:05 by shuang5186
[考研] 085602 化学工程专硕 340分求调剂 +4 qianbai11 2026-03-29 4/200 2026-03-30 11:34 by 唐沐儿
[考研] 085600,材料与化工321分求调剂 +10 大馋小子 2026-03-28 10/500 2026-03-29 23:35 by 飞行日记西
[考研] 304求调剂 +6 曼殊2266 2026-03-27 6/300 2026-03-28 14:10 by 唐沐儿
[考研] 285求调剂 +4 AZMK 2026-03-27 7/350 2026-03-27 20:59 by AZMK
[考研] 调剂 +4 柚柚yoyo 2026-03-26 4/200 2026-03-26 20:43 by fmesaito
信息提示
请填处理意见