| 查看: 1139 | 回复: 5 | |||
paladinlee铁杆木虫 (正式写手)
|
[交流]
【求助】dlpoly运算error——392的解决方法 已有4人参与
|
|
DL_POLY terminated due to error 392error - too many link cells requested 我用NPT系综算出的问题，NVT算就没问题。 config和field应该没问题，都是一样的。 control文件如下： integrator leapfrog verlet temperature 300.00 pressure 0.0001 ensemble npt berendsen 1.0000 1.0000 steps 100000 equilibration 20000 multiple step 1 scale 10 print 10 stack 100 stats 10 trajectory 0 1000 0 rdf 10 timestep 0.0050 primary cutoff 5.000 cutoff 5.000 delr width 1.000 rvdw cutoff 5.000 no electrostatics shake tolerance 1.0E-5 quaternion tolerance 1.0E-5 ... 我看了手册，说要改参数，具体怎么改？ 谢谢！ |
回复此楼» 猜你喜欢
293求调剂
已经有12人回复
-
070300化学319求调剂
已经有4人回复
-
268求调剂
已经有9人回复
-
301求调剂
已经有4人回复
-
0703化学336分求调剂
已经有5人回复
-
301求调剂
已经有4人回复
-
材料专硕306英一数二
已经有9人回复
-
304求调剂
已经有8人回复
-
材料工程专硕调剂
已经有5人回复
-
344求调剂
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】dlpoly运算结果的数据信息
已经有3人回复
paladinlee
铁杆木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 17609.8
- 红花: 1
- 帖子: 837
- 在线: 189.3小时
- 虫号: 1050994
- 注册: 2010-07-02
- 性别: GG
- 专业: 工程热物理相关交叉领域
2楼2010-10-13 16:23:58
043114076
木虫 (正式写手)
- 应助: 49 (小学生)
- 金币: 3364.9
- 散金: 52
- 红花: 23
- 帖子: 788
- 在线: 369.5小时
- 虫号: 570592
- 注册: 2008-06-07
- 性别: GG
- 专业: 理论和计算化学
★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zh1987hs(金币+4):谢谢 2010-10-13 19:00:14
小木虫(金币+0.5):给个红包，谢谢回帖交流
zh1987hs(金币+4):谢谢 2010-10-13 19:00:14
|
错误的原因是link cell数超过了系统默认最大的link cell数，最简单的办法是加入no link 参数。 NVT下没问题，NPT有问题：说明NPT下，体积增大，导致link cell 过大。首先，你的压力设置为0.0001，相当于0.1个大气压，是否有误？ 你的初始盒子大小是否有误？应尽量接近实验值。 如果初始构型太差的话，分子之间排斥占主导，也会引起体积过度变动。模拟前可先进行opt 或0kmd 驰豫一下。 另外，控压驰豫时间可以改小一点试试。 link cell 是一种加速模拟的方法，具体可以看Frenkel 和 smit 分子模拟那本书的附录。 [ Last edited by 043114076 on 2010-10-13 at 18:50 ] |
3楼2010-10-13 18:47:06
paladinlee
铁杆木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 17609.8
- 红花: 1
- 帖子: 837
- 在线: 189.3小时
- 虫号: 1050994
- 注册: 2010-07-02
- 性别: GG
- 专业: 工程热物理相关交叉领域
4楼2010-10-14 09:38:52
5楼2015-03-10 16:57:11
6楼2015-07-14 20:23:05
