24小时热门版块排行榜    

查看: 840  |  回复: 2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

明月照清渠

金虫 (小有名气)

[求助] 求助,gromacs跑nve能量为神马不恒定? 已有1人参与

求大神指教,为神马跑NVE的时候,能量不恒定,还降低呢?NVE设置如下,不知哪里有错?

constraints         = h-bonds
constraint_algorithm =  LINCS
lincs_iter        = 2                ; accuracy of LINCS
lincs_order        = 4                ; also related to accuracy
integrator          =  md
dt                  =  0.002    ; ps !
nsteps              =  250000000  ; 500 ns
nstxout             =  50000
nstvout             =  100000
nstfout             =  100000
nstlog              =  10000
nstenergy           =  50000
nstxtcout           =  50000
nstlist             =  10
ns_type             =  grid
coulombtype         =  PME-Switch
vdw_type            =  Shift
rvdw_switch         = 0.8
DispCorr            =  no
rlist               =  1.2
rcoulomb            =  1.0
rvdw                =  1.0
pbc                 = xyz
optimize_fft        = yes
nstcomm             = 1
comm-grps           = system
; Berendsen temperature coupling is on
tcoupl              =  no
tc-grps             =  system
tau_t               =  0.1 ;  0.1  0.1   
ref_t               =  300;  298  298   
; Energy monitoring
energygrps          =  system
; Semiisotropic pressure coupling is now on
;Pcoupl              =  Berendsen
;Pcoupltype          =  isotropic
;tau_p               =  1       ; 1      ; 5      5       5       5
;compressibility     =  4.5e-5  ;4.5e-5 ; 4.5e-5 0       0       0
;ref_p               =  1.0     ; 1.0    ; 1.0    0       0       0
; Generate velocites is on at 300K.
gen_vel             =  yes
gen_temp            =  300
gen_seed            =  173529
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

明月照清渠

金虫 (小有名气)

引用回帖:
2楼: Originally posted by fangsteel at 2014-12-31 15:36:36
mdp为啥把tcoupl设置成没有温度耦合呢?那不就是没有热浴了么???

nve是这样的,不能控温

[ 发自小木虫客户端 ]
3楼2014-12-31 19:40:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

fangsteel

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
mdp为啥把tcoupl设置成没有温度耦合呢?那不就是没有热浴了么???
gromacs
2楼2014-12-31 15:36:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 uKOSCocCb6Fm 2026-11-13 3/150 2026-11-15 04:05 by xSrQNSSjcgoX
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 uKOSCocCb6Fm 2026-11-12 7/350 2026-11-15 03:43 by xSrQNSSjcgoX
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 uKOSCocCb6Fm 2026-11-12 8/400 2026-11-15 02:26 by xSrQNSSjcgoX
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 FChKX8meu1ds 2026-11-09 9/450 2026-11-14 22:39 by Equinoxhua
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 FChKX8meu1ds 2026-11-09 8/400 2026-11-14 22:04 by Equinoxhua
[考博] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 TZqyNTBa07NP 2026-11-13 3/150 2026-11-14 15:58 by KXypaGgBBlw7
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 zSbvVX9JFl9J 2026-11-09 9/450 2026-11-13 11:13 by Equinoxhua
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 12:48 by ECLTXnJnRSXf
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 12:28 by ECLTXnJnRSXf
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 Z7qfCsdIiD1T 2026-11-10 7/350 2026-11-12 11:31 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 Z7qfCsdIiD1T 2026-11-10 7/350 2026-11-12 11:22 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 czAiolzR0IJC 2026-11-10 7/350 2026-11-12 11:09 by HajghKOMynmz
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 G2rhSvqvhNzl 2026-11-10 6/300 2026-11-12 10:49 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 G2rhSvqvhNzl 2026-11-10 6/300 2026-11-12 10:36 by HajghKOMynmz
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 YHiw1gou1Eda 2026-11-10 6/300 2026-11-12 10:08 by HajghKOMynmz
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 gRND4NVAZ8VG 2026-11-11 3/150 2026-11-12 07:42 by HajghKOMynmz
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 FChKX8meu1ds 2026-11-09 8/400 2026-11-12 05:30 by 9fOTqj074tIg
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 FChKX8meu1ds 2026-11-09 9/450 2026-11-12 05:19 by 9fOTqj074tIg
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 zSbvVX9JFl9J 2026-11-09 8/400 2026-11-12 04:57 by 9fOTqj074tIg
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 Z7qfCsdIiD1T 2026-11-10 4/200 2026-11-12 02:23 by 9fOTqj074tIg
信息提示
请填处理意见