| 查看: 692 | 回复: 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 |
回复此楼» 猜你喜欢
0856材料与化工调剂，339
已经有9人回复
-
一志愿西安交大材料学硕（英一数二）347，求调剂到高分子/材料相关专业
已经有7人回复
-
324分 085600材料与化工
已经有9人回复
-
考研调剂
已经有6人回复
-
材料专硕306英一数二
已经有3人回复
-
材料调剂
已经有5人回复
-
309分085801求调剂
已经有6人回复
-
一志愿武理材料工程302调剂环化或化工
已经有10人回复
-
311（085601）求调剂
已经有13人回复
-
335求调剂
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- gromacs全原子力场不能跑，求助！
已经有6人回复
fangsteel
超级版主
- 应助: 96 (初中生)
- 金币: 1213.6
- 散金: 32
- 红花: 22
- 帖子: 456
- 在线: 192.1小时
- 虫号: 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
130