计算化学公社

 找回密码 Forget password
 注册 Register
Views: 19725|回复 Reply: 22
打印 Print 上一主题 Last thread 下一主题 Next thread

[CASTEP/Dmol3/MS] the Hubbard-U correction如何操作?

[复制链接 Copy URL]

81

帖子

0

威望

593

eV
积分
674

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
文献中写道We have recently reported that inclusion of the Hubbard-U correction (i.e., GGA+U) leads to improvement of the description。。。
这个“+U”如何操作,恳请指点一二。我是小白。
万分感谢!!

65

帖子

0

威望

662

eV
积分
727

Level 4 (黑子)

23#
发表于 Post on 2018-6-11 09:27:10 | 只看该作者 Only view this author
卡开发发 发表于 2018-6-10 21:10
前面几楼的帖子提到的文献有数据以及其他的测试你也可以看D G Truhlar和S Grimme的文献,他们一般发展新 ...

嗯嗯,谢谢你的回答

3754

帖子

3

威望

1万

eV
积分
19676

Level 6 (一方通行)

围观吃瓜群众

22#
发表于 Post on 2018-6-10 21:10:30 | 只看该作者 Only view this author
wudazhuang509 发表于 2018-6-10 18:24
嗯嗯 好的,谢谢你,有没这种误差的分析的文献可以推荐下?

前面几楼的帖子提到的文献有数据以及其他的测试你也可以看D G Truhlar和S Grimme的文献,他们一般发展新的泛函总要大量测试。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

65

帖子

0

威望

662

eV
积分
727

Level 4 (黑子)

21#
发表于 Post on 2018-6-10 18:24:14 | 只看该作者 Only view this author
卡开发发 发表于 2018-5-29 12:00
实验值都±0.3的误差你也没啥好纠结的。这些近似的泛函本身也总会有误差,而且这种误差并不是说所有的体 ...

嗯嗯 好的,谢谢你,有没这种误差的分析的文献可以推荐下?

3754

帖子

3

威望

1万

eV
积分
19676

Level 6 (一方通行)

围观吃瓜群众

20#
发表于 Post on 2018-5-29 12:00:14 | 只看该作者 Only view this author
wudazhuang509 发表于 2018-5-29 11:29
完全搞不懂这种计算多少误差能够接受,希望你能理解,以前做理论的。纠结体质啊

实验值都±0.3的误差你也没啥好纠结的。这些近似的泛函本身也总会有误差,而且这种误差并不是说所有的体系系统性低估或高估某个值,到了不同的体系表现是不一样的。比如上面的TiO2,HSE06能给出和实验相当接近的带隙,而文献给出CdSe、ZnSe也都能低估0.4eV左右。

日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

65

帖子

0

威望

662

eV
积分
727

Level 4 (黑子)

19#
发表于 Post on 2018-5-29 11:29:17 | 只看该作者 Only view this author
kyuu 发表于 2018-5-28 22:34
Vasp Paw算出来的TiO2带隙2.67, GGA没加U 还纠结不?

完全搞不懂这种计算多少误差能够接受,希望你能理解,以前做理论的。纠结体质啊

357

帖子

0

威望

2069

eV
积分
2426

Level 5 (御坂)

真 掘墓者

18#
发表于 Post on 2018-5-28 22:34:23 | 只看该作者 Only view this author
本帖最后由 kyuu 于 2018-5-28 22:37 编辑
wudazhuang509 发表于 2018-5-24 15:23
嗯嗯 谢谢你了。
实验的带隙是2.8+-0.3,PBE0计算出来是2.36。不知道是否可以达到可以作为文章的一个使 ...

Vasp Paw算出来的TiO2带隙2.67, GGA没加U 还纠结不?
圣诞刨坟忙

65

帖子

0

威望

662

eV
积分
727

Level 4 (黑子)

17#
发表于 Post on 2018-5-24 19:30:12 | 只看该作者 Only view this author
卡开发发 发表于 2018-5-24 17:00
应该可以被接受,http://bbs.keinsci.com/thread-3000-1-1.html这里给的文章HSE06也有低估这么多的体系。

太感谢你了

3754

帖子

3

威望

1万

eV
积分
19676

Level 6 (一方通行)

围观吃瓜群众

16#
发表于 Post on 2018-5-24 17:00:44 | 只看该作者 Only view this author
wudazhuang509 发表于 2018-5-24 15:23
嗯嗯 谢谢你了。
实验的带隙是2.8+-0.3,PBE0计算出来是2.36。不知道是否可以达到可以作为文章的一个使 ...

应该可以被接受,http://bbs.keinsci.com/thread-3000-1-1.html这里给的文章HSE06也有低估这么多的体系。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

65

帖子

0

威望

662

eV
积分
727

Level 4 (黑子)

15#
发表于 Post on 2018-5-24 15:23:44 | 只看该作者 Only view this author
卡开发发 发表于 2018-5-24 13:37
2、
(1)看你材料本身的带隙吧,如果本身带隙就很窄的话相对偏差就很大,但一般即便是比较好的方法低估 ...

嗯嗯 谢谢你了。
实验的带隙是2.8+-0.3,PBE0计算出来是2.36。不知道是否可以达到可以作为文章的一个使用数据。所以很纠结

3754

帖子

3

威望

1万

eV
积分
19676

Level 6 (一方通行)

围观吃瓜群众

14#
发表于 Post on 2018-5-24 13:37:17 | 只看该作者 Only view this author
wudazhuang509 发表于 2018-5-24 12:13
谢谢你的回答,
1,第二种是指你之前回答的 2. CASTEP的Calculation开启自旋极化和LDA+U(注:这里的LDA ...

2、
(1)看你材料本身的带隙吧,如果本身带隙就很窄的话相对偏差就很大,但一般即便是比较好的方法低估0.几也是很可能的,论坛内应该发过相关的文献有测试数据。
(2)赝势的质量会影响计算结果,但有多大影响我也很难估计,CASTEP的老的模守恒赝势质量确实糟糕点。
(3)GW指的是多体Green函数方法,如Yambo或者BerkeleyGW可以做,VASP也能做。
如果实验数据无误且计算参数选择可靠,也只好这样,如果为了去凑一个带隙意义就不是很大了,

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
sobereva + 2

查看全部评分 View all ratings

日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

65

帖子

0

威望

662

eV
积分
727

Level 4 (黑子)

13#
发表于 Post on 2018-5-24 12:13:13 | 只看该作者 Only view this author
卡开发发 发表于 2018-5-24 05:19
1、Ag的组态是4d10,+U可能用处不大。你指的第二种是什么?
2、得看低多少,如果只是低0.几那结果就算不 ...

谢谢你的回答,
1,第二种是指你之前回答的 2. CASTEP的Calculation开启自旋极化和LDA+U(注:这里的LDA+U指的就是DFT+U,XC Functional实际上可以选择LDA或者GGA,并非只针对LDA)。
2,能带小0.34eV,是不是就可以认为是合理的呢?另外sX-LDA怎么好像得到的比PBE0还要小,还有GW是指的什么呢?

3754

帖子

3

威望

1万

eV
积分
19676

Level 6 (一方通行)

围观吃瓜群众

12#
发表于 Post on 2018-5-24 05:19:46 | 只看该作者 Only view this author
wudazhuang509 发表于 2018-5-23 22:13
请问,
1. 在计算AgF能带的时候,我使用第一种方法+U,d=5, 10, 15之后,计算得到的结果和PBE计算出来一 ...

1、Ag的组态是4d10,+U可能用处不大。你指的第二种是什么?
2、得看低多少,如果只是低0.几那结果就算不错了。PBE0如果算出来带隙还比较小的话你可能得先确定一下计算参数是否经过充分测试,保证你的计算参数选择合理,能带计算对这些还是比较敏感的,如果还不行也只好试试sX-LDA,再不行也就只好去跑GW了。
3、三种泛函大致上是这样:
(1)CASTEP的LDA本来是Slater交换泛函+PZ(Perdew-Zunger)的关联泛函;
(2)HF-LDA指的是HF的交换部分额外加上PZ的关联泛函构造的杂化泛函;
(3)sX-LDA则是介于这两者之间,类似于HSE,通过Thomas-Fermi屏蔽来使得低密度情形渐进为HF-LDA,高密度则被屏蔽,渐进为LDA。

评分 Rate

参与人数
Participants 1
eV +3 收起 理由
Reason
sobereva + 3

查看全部评分 View all ratings

日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

65

帖子

0

威望

662

eV
积分
727

Level 4 (黑子)

11#
发表于 Post on 2018-5-23 22:13:23 | 只看该作者 Only view this author
卡开发发 发表于 2018-2-7 19:16
1、Modify-Electronic Configuration-Hubbard U,然后针对某个角动量设置U值。
2、CASTEP的Calculation开 ...

请问,
1. 在计算AgF能带的时候,我使用第一种方法+U,d=5, 10, 15之后,计算得到的结果和PBE计算出来一样,是不是说明这个体系不适合+U呢?还是说要同时使用第二种?
2. 我用PBE0计算出来的带隙偏小,有没其他方法提高呢?
3. sX-LDA和HF-LDA与LDA有什么区别呢?或者说在计算能带中能不能提高带隙?
谢谢

183

帖子

2

威望

4931

eV
积分
5154

Level 6 (一方通行)

能科top50

10#
发表于 Post on 2018-3-6 19:41:25 | 只看该作者 Only view this author
卡开发发 发表于 2018-3-6 18:21
是否DFT+U,U值取多少一直是固体和表面计算中有争议的问题。
这里提供一个反面的观点的文章,Phys. Chem ...

谢谢卡开发发老师的解答,您提供的文献给我提供了很大帮助,我打算参考其中的观点回复审稿意见
In defeat, malice. In victory, revenge.

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-16 02:57 , Processed in 0.235724 second(s), 22 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list