| 查看: 789 | 回复: 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 |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有18人回复
-
基金系统什么内容也没有
已经有6人回复
-
国自然面上复盘~欢迎讨论
已经有12人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
面上合作单位盖章
已经有5人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有7人回复
-
申请删除本帖
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- gromacs全原子力场不能跑，求助！
已经有6人回复
fangsteel
木虫 (正式写手)
- 应助: 96 (初中生)
- 金币: 1215.1
- 散金: 32
- 红花: 22
- 帖子: 456
- 在线: 192.5小时
- 虫号: 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
投票:
10