| 查看: 2526 | 回复: 9 | ||||
[求助]
黄铜矿几何结构优化 按照文献参数不能收敛，且晶胞常数误差过大如何解决 已有1人参与
|
|
小弟按照文献所述参数对黄铜矿进行了几何结构的优化，但是却一直不能收敛是怎么回事，而且算得的晶胞常数与实验值相差超过了2%，甚至达到了4%或5%，是什么地方出现了问题吗？下面吧具体模拟过程及文献提供给各位大神分析，谢谢！ 文献中的设置如下：黄铜矿，理想黄铜矿的晶胞参数为：a=b=5.289 Å，c=10.423 Å，α=β=γ=90°。采用GG-PBE梯度修正函数，在进行几何优化时，选取合适的截断能和k点，收敛标准为：原子位移为0.005 Å，原子间作用力为0.01 eV/nm，原子间的内应力为0.2 GPa，体系总能量的变化不大于5.0×10-4 eV/atom。 截断能和k点选择过程及数据如图一和图二： 我按照245计算，计算到98步进行不下去了，受不了就停止了。 手头上有一个270ev的数据，收敛过程如图三图四，设置参数如图五图六。 另外我算了E点为310，350和400的，k点位5x5x2，很快收敛，但是优化后的晶胞常数误差超过2%，最大的达到4%，应该是不能用的。 请问1.为什么我按照文献设置不能收敛？ 2.为什么提高E点可以收敛，但是晶胞常数误差过大？是不是K点设置问题，还是 其他问题？ 谢谢指教！  360截图20140718223046593.jpg  360截图20140718223114573.jpg  360截图20140718223437784.jpg  360截图20140718223452372.jpg  360截图20140718223536488.jpg  360截图20140718223610451.jpg |
回复此楼» 猜你喜欢
0703化学336分求调剂
已经有4人回复
-
[复试调剂]西南科技大学国防/材料导师推荐
已经有6人回复
-
化学工程321分求调剂
已经有12人回复
-
211本，11408一志愿中科院277分，曾在中科院自动化所实习
已经有4人回复
-
材料专硕326求调剂
已经有5人回复
-
东南大学364求调剂
已经有5人回复
-
国自科面上基金字体
已经有7人回复
-
药学383 求调剂
已经有4人回复
-
286求调剂
已经有5人回复
-
085601求调剂
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- CASTEP几何结构优化的问题 虽然收敛了 感觉不是很对 请大家帮忙找一下原因
已经有11人回复
- 收敛测试和结构优化
已经有10人回复
- 【讨论】MS中sorption模块的MC计算
已经有24人回复
- 【求助】castep中的参数对计算结果有什么影响
已经有11人回复
huxiao.gl
木虫 (小有名气)
- 应助: 37 (小学生)
- 金币: 2939.4
- 红花: 4
- 帖子: 191
- 在线: 373.4小时
- 虫号: 1849937
- 注册: 2012-06-06
- 性别: GG
- 专业: 理论和计算化学
yeahgyz
2楼2014-07-19 02:48:21
3楼2014-07-19 09:01:42
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 1ST强帖: 2
- 应助: 47 (小学生)
- 金币: 3507.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
4楼2014-07-19 09:55:00
huxiao.gl
木虫 (小有名气)
- 应助: 37 (小学生)
- 金币: 2939.4
- 红花: 4
- 帖子: 191
- 在线: 373.4小时
- 虫号: 1849937
- 注册: 2012-06-06
- 性别: GG
- 专业: 理论和计算化学
5楼2014-07-19 16:42:34
6楼2014-07-19 16:55:59
7楼2014-07-19 17:01:39
huxiao.gl
木虫 (小有名气)
- 应助: 37 (小学生)
- 金币: 2939.4
- 红花: 4
- 帖子: 191
- 在线: 373.4小时
- 虫号: 1849937
- 注册: 2012-06-06
- 性别: GG
- 专业: 理论和计算化学
8楼2014-07-19 20:53:54
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 1ST强帖: 2
- 应助: 47 (小学生)
- 金币: 3507.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
|
只计算能量而不解析导数优化结构，即单点能（single point）。 能不能吻合实验在于XC functional和pseudopotential的选取，怎么选取不妨参考不少测试的文献。据我（经验）所知，一些离子型的结构PBEsol稍微好些，不过均匀结构（如纯金属）往往LDA就能给出好的结果，这些往往在一些泛函性能评估的文献能够查到。重点看后续工作要做什么，目前没有一种泛函能够保证任何性质都和实验吻合。另外我了解到的情况是含Cu的晶格常数似乎经常不准确，也许和关联效应有关（这点目前我不确定）。 测试过程并不是一味增加。事实上，楼上的结论正确，2%的偏差不算大。如我前面所说，如果泛函的性能不合适，提升精度反倒是会产生偏离，粗略的精度恰好更接近实验也只是巧合，这种准确对理论计算其实没啥意义。 |
9楼2014-07-19 23:05:18
10楼2016-06-01 22:06:46
5