| 查看: 787 | 回复: 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 |
回复此楼» 猜你喜欢
基金未中，这种答复是模板吗？
已经有3人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有5人回复
-
有多少人是今天查系统知道结果的？
已经有17人回复
-
怎么看青基中了没有啊
已经有5人回复
-
梦想
已经有3人回复
-
面上合作单位盖章
已经有3人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有5人回复
-
基金不中，共勉
已经有7人回复
-
为什么 国际(地区)合作与交流项目 没有放榜？
已经有11人回复
-
让我中一个面上吧！
已经有17人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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