| 查看: 1124 | 回复: 6 | |||
[交流]
【求助】DL-POLY的error185出现应该怎么处理已有3人参与
|
|
各位大家好， 我DL－POLY运行时出现了error185，我不知道是哪里错了，哪里须改进，知道的高手们请赐教，谢谢各位了！ |
回复此楼» 猜你喜欢
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有6人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
三甲基碘化亚砜的氧化反应
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
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