| 查看: 809 | 回复: 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 |
回复此楼» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有7人回复
-
浙江师范大学是怎么坑我的
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- gromacs全原子力场不能跑，求助！
已经有6人回复
fangsteel
木虫 (正式写手)
- 应助: 96 (初中生)
- 金币: 1206.1
- 散金: 32
- 红花: 22
- 帖子: 457
- 在线: 192.9小时
- 虫号: 2248690
- 注册: 2013-01-18
- 性别: GG
- 专业: 原子和分子物理
2楼2014-12-31 15:36:36
明月照清渠
金虫 (小有名气)
- 应助: 10 (幼儿园)
- 金币: 3550.4
- 散金: 1
- 红花: 3
- 帖子: 95
- 在线: 69.3小时
- 虫号: 1101813
- 注册: 2010-09-17
- 性别: MM
- 专业: 理论和计算化学
3楼2014-12-31 19:40:35
