24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1673  |  回复: 8

hyliu0618

铁虫 (正式写手)

[求助] smearing到底选多大合适?

看过一些帖子,前辈们认为smearing选0.005Ha 最好。
现在我在计算FeOOH的体系,用dmol软件包,稳定构型选0.005Ha 有时根本不收敛,而且选0.05Ha算出的能量比选0.005Ha 还要低,也就是说选0.05Ha构型更稳定。更要命的是搜索过渡态时,smearing选0.005Ha 一个点都跑不了,0.05Ha还勉强可以找个过渡态。到底我该怎么选smearing,请有经验的虫友们不吝赐教!!!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

计算虫子
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hookah

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
hyliu0618: 金币+1 2013-01-16 11:52:43
hyliu0618: 金币+5, ★★★很有帮助 2013-01-17 10:50:33
liliangfang: 金币+1, 谢谢交流 2013-01-20 16:07:37
这个选择,没有严格的标准,一般都是在做一个体系前对于smearing参数做测试的。也分你要计算什么东西。
2楼2013-01-16 09:23:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hyliu0618

铁虫 (正式写手)

引用回帖:
2楼: Originally posted by hookah at 2013-01-16 09:23:15
这个选择,没有严格的标准,一般都是在做一个体系前对于smearing参数做测试的。也分你要计算什么东西。

我在计算FeOOH的体系
计算虫子
3楼2013-01-16 11:53:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hisap-nano

新虫 (正式写手)

引用回帖:
3楼: Originally posted by hyliu0618 at 2013-01-16 11:53:12
我在计算FeOOH的体系...

你可以参考下面的文献
Metals at finite temperature a modified smearing scheme
你用google搜一下,希望对你有帮助!
4楼2013-01-17 06:03:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wxqdelphi

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hyliu0618

铁虫 (正式写手)

引用回帖:
5楼: Originally posted by wxqdelphi at 2013-01-17 08:30:02
“选0.05Ha算出的能量比选0.005Ha 还要低,也就是说选0.05Ha构型更稳定”,这个不一定。0.05的结果中,smearing contribution 的energy会很大,很多时候0.05Ha比0.005Ha的能量更低很正常,但0.005Ha才是对的,因为其 ...

请问这位前辈在DMOL中如何考虑d电子的强相关,如果您没有用过这个软件包,能否简单介绍一下这个强相关该如何处理,谢谢!
计算虫子
6楼2013-01-17 10:52:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wxqdelphi

金虫 (正式写手)

【答案】应助回帖

引用回帖:
6楼: Originally posted by hyliu0618 at 2013-01-17 10:52:12
请问这位前辈在DMOL中如何考虑d电子的强相关,如果您没有用过这个软件包,能否简单介绍一下这个强相关该如何处理,谢谢!...

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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zmrright

金虫 (著名写手)

引用回帖:
5楼: Originally posted by wxqdelphi at 2013-01-17 08:30:02
“选0.05Ha算出的能量比选0.005Ha 还要低,也就是说选0.05Ha构型更稳定”,这个不一定。0.05的结果中,smearing contribution 的energy会很大,很多时候0.05Ha比0.005Ha的能量更低很正常,但0.005Ha才是对的,因为其 ...

你说的restart是什么意思?收敛后,再把smearing变成0.005再优化吗?
做好自我,不枉此生
8楼2013-09-28 21:45:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zmrright

金虫 (著名写手)

引用回帖:
5楼: Originally posted by wxqdelphi at 2013-01-17 08:30:02
“选0.05Ha算出的能量比选0.005Ha 还要低,也就是说选0.05Ha构型更稳定”,这个不一定。0.05的结果中,smearing contribution 的energy会很大,很多时候0.05Ha比0.005Ha的能量更低很正常,但0.005Ha才是对的,因为其 ...

“选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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hyliu0618 的主题更新
信息提示
请填处理意见