24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 563  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 理学07化学 303求调剂 +12 睿08 2026-03-27 12/600 2026-04-01 22:36 by Dyhoer
[考研] 311求调剂 +13 zchqwer 2026-04-01 15/750 2026-04-01 22:30 by 2026材料调剂
[考研] 江苏科技大学招材料研究生 +4 Su032713. 2026-04-01 5/250 2026-04-01 22:03 by cccchenso
[考研] 266分求材料化工冶金矿业等专业的调剂 +7 哇呼哼呼哼 2026-03-26 7/350 2026-04-01 20:47 by cqupH
[考研] 307分求调剂 +14 (o~o) 2026-03-31 15/750 2026-04-01 20:43 by longlotian
[考研] 建环,能源,土木老师路过看一看!!! +4 嘿嘿uu 2026-04-01 4/200 2026-04-01 20:42 by 无懈可击111
[考研] 一志愿211,335分,0856,求调剂院校和导师 +11 倾____萧 2026-03-27 12/600 2026-04-01 15:27 by oooqiao
[考研] 一志愿 南京航空航天大学 ,080500材料科学与工程学硕 +7 @taotao 2026-03-30 7/350 2026-04-01 14:30 by chenqifeng666
[硕博家园] 博一被送出联培感觉不适应怎么办 +3 全村的狗 2026-03-31 3/150 2026-04-01 10:44 by 328838485
[考研] 一志愿 南京航空航天大学 ,080500材料科学与工程学硕 +10 @taotao 2026-03-31 11/550 2026-04-01 09:43 by xiayizhi
[考研] 求调剂,一志愿 南京航空航天大学 ,080500材料科学与工程学硕,总分289分 +10 @taotao 2026-03-29 10/500 2026-04-01 09:30 by oooqiao
[考研] 318一志愿吉林大学生物与医药 求调剂 +6 笃行致远. 2026-03-28 6/300 2026-04-01 09:28 by oooqiao
[考研] 289求调剂 +7 BrightLL 2026-03-29 7/350 2026-03-31 22:05 by 544594351
[考研] 352分-085602-一志愿985 +6 海纳百川Ly 2026-03-29 6/300 2026-03-31 21:06 by yuq
[考研] 合肥区域性重点一本招收调剂 +4 6266jl 2026-03-30 8/400 2026-03-31 18:43 by 6266jl
[考研] 总分322求生物学/生化与分子/生物信息学相关调剂 +6 星沉uu 2026-03-26 7/350 2026-03-31 10:19 by GdShizy
[考研] 11408总分309,一志愿东南大学求调剂,不挑专业 +5 天赋带到THU 2026-03-29 6/300 2026-03-30 20:49 by dick_runner
[考研] 085602 化工专硕 338分 求调剂 +12 路痴小琪 2026-03-27 12/600 2026-03-28 15:41 by L135790
[考研] 一志愿南京航空航天大学材料学硕求调剂 +3 @taotao 2026-03-28 3/150 2026-03-28 10:26 by JourneyLucky
[考研] 070300化学求调剂 +4 起个名咋这么难 2026-03-27 4/200 2026-03-27 21:39 by 83503孙老师
信息提示
请填处理意见