24小时热门版块排行榜    

查看: 1900  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +7 Lanmanbaby 2026-08-09 11/550 2026-08-09 17:48 by zhanghaozhu
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 好奇怪的filecode +3 布布和一二 2026-08-08 4/200 2026-08-08 17:46 by zhanghaozhu
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:27 by oEVWOejN9taj
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 16:47 by oEVWOejN9taj
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:02 by oEVWOejN9taj
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:47 by oEVWOejN9taj
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 10:19 by 3OOjAIS77qg2
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见