| 查看: 786 | 回复: 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 |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有17人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有3人回复
-
基金不中，共勉
已经有5人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有3人回复
-
国自然面上复盘~欢迎讨论
已经有5人回复
-
我不理解！
已经有23人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有10人回复
-
系统查不到
已经有10人回复
-
怎么查啊
已经有5人回复
-
范进中举一文的中心思想
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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
投票:
