| 查看: 570 | 回复: 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 请各位指点一下，谢谢 |
回复此楼» 猜你喜欢
献血感触
已经有4人回复
-
评审有感
已经有16人回复
-
反应很差，大量原料没有反应
已经有3人回复
-
Sci. Bull. 悲剧经验
已经有4人回复
-
找博士生导师
已经有7人回复
-
上海大学实验技术岗位非升即走
已经有11人回复
-
26/27申博自荐-锂/钠电池方向
已经有4人回复
-
同样的基金本子，换个专家直接从C变A！
已经有3人回复
-
别被青基扩招骗了！26年科研内卷才刚刚开始
已经有4人回复
-
26/27博士推荐
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 分子模拟的时间问题
已经有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
30