| 查看: 1475 | 回复: 6 | |||
[交流]
【求助】DL-POLY的error185出现应该怎么处理 已有3人参与
|
|
各位大家好， 我DL－POLY运行时出现了error185，我不知道是哪里错了，哪里须改进，知道的高手们请赐教，谢谢各位了！ |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有7人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有6人回复
-
各位大神，目前国内有哪些比较好用的逆合成软件？
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
2楼2010-05-31 08:17:17
3楼2010-05-31 15:23:07
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
4楼2010-06-01 08:27:49
|
各位下面是我的CONTROL文件： integration leapfrog verlet ensemble npt hoover 1.0 5.0 temperature 373.0 pressure 0.001 steps 50000 equilibration 40000 mult 1 scale 10 print 100 stack 100 stats 100 trajectory 1 100 1 rdf 10 timestep 0.0005 cutoff 18.0 delr width 1.2 ewald precision 1.0e-9 cap 1.0000E+04 quaternion tolerance 1.0000E-05 zdensity 0.1 80.0 print rdf job time 99999900.0 close time 200000.0 finish 我改了很多次有关静电加和的参数ewald precision和cutoff 的数值，还是出现错误185，说明中让增大cutoff，在Ewald sum中用更少的k矢量和重新设置参数kmaxb, kmaxc，我增大了cutoff，不知道该怎么减少k矢量，我的CONTROL中好像没有用到kmaxb, kmaxc参数，我就改变了ewald precision的数值但是没用，请各位好心人帮忙看一下提些建议好吗？我将不胜感激！！！ |
5楼2010-06-01 11:16:53
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
★ ★ ★ ★ ★
ghcacj(金币+5):谢谢 2010-06-01 16:02:48
tang327(金币+5): 2010-06-02 09:44:06
ghcacj(金币+5):谢谢 2010-06-01 16:02:48
tang327(金币+5): 2010-06-02 09:44:06
|
A。 ewald precision，这个值你可以放大一点，看看效果如何。你搞到1e-9，精度太高了一点吧。1e-6，甚至更低也是可以的。 B。如果这一招不行，你可以试试手动设置ewald sum alpha k1 k2 k3。对了，你对Ewald了解多少？如果不了解，要仔细阅读Manual的2.4节，尤其是2.4.6节，如何设置这些值，要读第3.2.5节。 我对有机分子没有经验，如果是熔盐，一般alpha＝0.3~0.6就差不多了吧。K1，K2，K3之类，一般在5，6或者7，你也可以搞小一点，比如4或者5。你如果是方格，三个数一般是相等的。 先多试试方法A。反复试验。cutoff一般不要乱动，因为它和你的势函数还有一定的关联，这个值你最不要大过元胞的半边长。 |
6楼2010-06-01 15:46:17
7楼2010-06-02 20:25:58
