| 查看: 540 | 回复: 2 | ||||
ddwy木虫 (小有名气)
|
[求助]
分子动力学模拟出错 已有1人参与
|
|
大家好： 我用DL_POLY做分子动力学模拟，CONTROL输入文件是： temperature 298.15 # pressure 0.001 ensemble nvt hoover 0.5 integrat velocity restart steps 3000000 equilibration 1000000 mult 1 scale 5 stats 20 rdf 10 print rdf traj 1000001 100 0 densvar 20 timestep 0.002 cutoff 12.0 ewald precision 1.0E-6 minimise energy 50 5.0E-03 shake 1.0E-6 mxshak 10000 quaternion tolerance 1.0E-05 job time 99999999999999.0 close time 200000.0 finish 但是计算到第步的时候出现错误，OUTPUT文件为： ---------------------------------------------------------------------------------------------------------------------------------- 174600 -2.8822E+04 2.9815E+02 -3.1606E+04 3.9295E+03 -3.5876E+04 1.7169E+02 1.0398E+02 6.5269E+01 0.0000E+00 349.20000 -2.9405E+04 0.0000E+00 9.0503E+03 -9.0711E+04 3.5875E+04 8.3469E+01 0.0000E+00 6.3803E+04 0.0000E+00 119751.3 4.9579E+04 0.0000E+00 0.0000E+00 0.0000E+00 9.0000E+01 9.0000E+01 9.0000E+01 0.0000E+00 -1.2860E+00 rolling -2.6814E+04 2.9487E+02 -2.9563E+04 4.1468E+03 -3.4341E+04 3.6908E+02 1.8861E+02 7.4187E+01 0.0000E+00 averages -2.6525E+04 0.0000E+00 6.3658E+03 -9.3434E+04 3.4340E+04 1.4041E+03 0.0000E+00 6.4056E+04 0.0000E+00 4.9579E+04 0.0000E+00 0.0000E+00 0.0000E+00 9.0000E+01 9.0000E+01 9.0000E+01 0.0000E+00 -8.0166E-02 ---------------------------------------------------------------------------------------------------------------------------------- *** warning - global constraint number 6 , with particle numbers: 2885 & 2886 , converges to a length of 7.94 Angstroms with a factor 3.1039E+01 ,contributes towards next error !!! *** *** warning - global constraint number 7 , with particle numbers: 2885 & 2887 , converges to a length of 12.88 Angstroms with a factor 8.2412E+01 ,contributes towards next error !!! *** *** warning - global constraint number 8 , with particle numbers: 2886 & 2887 , converges to a length of 18.15 Angstroms with a factor 1.0006E+02 ,contributes towards next error !!! *** DL_POLY_4 terminated due to error 105 error - shake algorithm (constraints_shake) failed to converge "OUTPUT" 64939L, 6075639C 请各位指点一下，谢谢 |
回复此楼» 猜你喜欢
一志愿东华大学化学070300，求调剂
已经有8人回复
-
277材料科学与工程080500求调剂
已经有7人回复
-
317求调剂
已经有18人回复
-
293求调剂
已经有5人回复
-
280分求调剂 一志愿085802
已经有7人回复
-
0854电子信息求调剂
已经有3人回复
-
263求调剂
已经有4人回复
-
石河子大学（211、双一流）硕博研究生长期招生公告
已经有3人回复
-
一志愿西安交通大学材料工程专业 282分求调剂
已经有13人回复
-
北科281学硕材料求调剂
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 分子模拟的时间问题
已经有4人回复
- 分子动力学模拟的时间、体系大小和温度的争议问题！！
已经有23人回复
- 分子动力学模拟时间的意义？
已经有15人回复
- 分子动力学模拟轨迹文件格式的选择问题
已经有4人回复
- 分子动力学模拟“边界问题”含义
已经有10人回复
- 对接的分子用于分子动力学模拟的问题
已经有12人回复
- 拉伸分子动力学步骤
已经有8人回复
- Amber 分子动力学 求助！！！！
已经有10人回复
- 关于MS做分子动力学的问题
已经有8人回复
- 分子动力学结果分析
已经有7人回复
- 分子动力学模拟中的系综问题
已经有5人回复
- 分子动力学模拟完毕之后要怎么分析结果阿
已经有7人回复
- 分子动力学做冲击求助
已经有11人回复
- 分子动力学模拟长程库仑势的处理问题
已经有8人回复
- 请教分子力学与分子动力学的区别
已经有14人回复
- 请问各位熟悉分子动力学模拟的朋友，熟悉精通MS的分子动力学模拟需要多久？
已经有17人回复
- 【求助】求推荐分子动力学模拟的教程数据或ppt，谢谢
已经有6人回复
- 【求助】分子动力学模拟求扩散系数，最后一步不明白，请求给位大侠。
已经有21人回复
- 【求助】关于vasp进行分子动力学模拟的求助
已经有12人回复
- 【求助】分子动力学模拟尺度问题？
已经有11人回复
- 【求助】关于用分子动力学来做光谱的问题
已经有14人回复
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
【答案】应助回帖
感谢参与，应助指数 +1
|
http://www.uark.edu/ua/fengwang/DLPOLY2/node350.html Message 105: error - shake algorithm failed to converge The RD-SHAKE algorithm for bond constraints is iterative. If the maximum number of permitted iterations is exceeded, the program terminates. Possible causes include: a bad starting configuration; too large a time step used; incorrect force field specification; too high a temperature; inconsistent constraints involving shared atoms etc. Action: Corrective action depends on the cause. It is unlikely that simply increasing the iteration number will cure the problem, but you can try: follow the standard user response to increase the control parameter mxshak. But the trouble is much more likely to be cured by careful consideration of the physical system being simulated. For example, is the system stressed in some way? Too far from equilibrium? 试试从以上提出的可能性来查找错误。另外，怎么没看见shake tolerance关键字？ |
2楼2014-09-04 10:45:22
ddwy
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1713.9
- 散金: 15
- 红花: 3
- 帖子: 140
- 在线: 29.2小时
- 虫号: 2224794
- 注册: 2013-01-06
- 性别: MM
- 专业: 项目管理
3楼2014-09-05 16:13:18
