| 查看: 2066 | 回复: 8 | |||
[求助]
smearing到底选多大合适？
|
|
看过一些帖子，前辈们认为smearing选0.005Ha 最好。 现在我在计算FeOOH的体系，用dmol软件包，稳定构型选0.005Ha 有时根本不收敛，而且选0.05Ha算出的能量比选0.005Ha 还要低，也就是说选0.05Ha构型更稳定。更要命的是搜索过渡态时，smearing选0.005Ha 一个点都跑不了，0.05Ha还勉强可以找个过渡态。到底我该怎么选smearing,请有经验的虫友们不吝赐教！！！ |
回复此楼» 猜你喜欢
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
考博求助
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- origin图中的字号选多大合适？
已经有16人回复
2楼2013-01-16 09:23:15
3楼2013-01-16 11:53:12
hisap-nano
新虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 924.1
- 红花: 2
- 帖子: 412
- 在线: 226.3小时
- 虫号: 297318
- 注册: 2006-11-18
- 专业: 金属材料的力学行为
4楼2013-01-17 06:03:17
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
hyliu0618: 金币+8, ★★★很有帮助 2013-01-17 10:24:14
liliangfang: 金币+3, 谢谢交流 2013-01-20 16:07:45
感谢参与，应助指数 +1
hyliu0618: 金币+8, ★★★很有帮助 2013-01-17 10:24:14
liliangfang: 金币+3, 谢谢交流 2013-01-20 16:07:45
|
“选0.05Ha算出的能量比选0.005Ha 还要低，也就是说选0.05Ha构型更稳定”，这个不一定。0.05的结果中，smearing contribution 的energy会很大，很多时候0.05Ha比0.005Ha的能量更低很正常，但0.005Ha才是对的，因为其smearing contribution要小很多。 过渡态开始选0.005Ha跑不动很正常，可能某一个或几个image因为初始degauss过低，无法收敛。可以先用0.05跑，在收敛之后，restart，使用更低的degauss，会对收敛有帮助。 对Fe这类过渡金属体系，不知道你是否考虑了d电子的强相关作用，如果没有考虑，不收敛很正常。 |
5楼2013-01-17 08:30:02
6楼2013-01-17 10:52:12
【答案】应助回帖
|
DMOL没用过，记得DMOL用的不是平面波吧，对DMOL不熟。 VASP里面DFT+U，Hybrid应该是现在比较成熟的解决办法，个人愚见，仅供参考。 另外，只是因为你的体系很难收敛，我才想到可能是强关联作用造成的，可能fe的d电子巡游地很厉害，没有localize，这是不能收敛的一种可能性，但并不是唯一的可能性。 如果你用dmol的话，可以试试加大smearing，或者加大ecut，先令体系收敛，然后降低smearing，再restart。不过smearing也和你要计算的band数量以及smearing方法有很大关系，解决不收敛的问题很confusing，需要多次尝试，特别是对这种强关联体系。 建议你先看文献，文献没有借鉴看类似体系，这样可以给你很多启发。 |
7楼2013-01-17 11:00:09
zmrright
金虫 (著名写手)
- 应助: 49 (小学生)
- 金币: 3460.8
- 散金: 706
- 红花: 9
- 帖子: 1975
- 在线: 299.1小时
- 虫号: 828781
- 注册: 2009-08-16
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
8楼2013-09-28 21:45:17
zmrright
金虫 (著名写手)
- 应助: 49 (小学生)
- 金币: 3460.8
- 散金: 706
- 红花: 9
- 帖子: 1975
- 在线: 299.1小时
- 虫号: 828781
- 注册: 2009-08-16
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
|
“选0.05Ha算出的能量比选0.005Ha 还要低，也就是说选0.05Ha构型更稳定”，这个不一定。0.05的结果中，smearing contribution 的energy会很大，很多时候0.05Ha比0.005Ha的能量更低很正常，但0.005Ha才是对的，因为其smearing contribution要小很多。 过渡态开始选0.005Ha跑不动很正常，可能某一个或几个image因为初始degauss过低，无法收敛。可以先用0.05跑，在收敛之后，restart，使用更低的degauss，会对收敛有帮助。 对Fe这类过渡金属体系，不知道你是否考虑了d电子的强相关作用，如果没有考虑，不收敛很正常。 你好，你说的restart是不是把smearing变小再优化一次啊？ |
9楼2013-09-28 21:46:37
10