24小时热门版块排行榜    

查看: 11892  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

xiaomao_16

银虫 (小有名气)

[求助] CASTEP计算带隙比实验值小很多怎么解决

1. 如题,新手刚接触CASTEP,看着Tutorials里的教程自己算了个块体MgO的能带,感觉MgO应该没有太特殊的地方吧,但算出来只有4.1eV,实验值足足7.8eV,有文章说这是系统误差,但是不是有点差太多了,有没有解决的方法呢?
设置参数如下:GGA-PW91(其他像PBE也试过结果差不多),截断能500eV,USPP,k点取样之类都选择的是Fine,基本就是按照教程默认来的

2. 本来的目标是NiO,也试着找文献的参数算了,给Ni的d轨道施加Hubbard U10eV,GGA-PBEsol,USPP,其他还是基本设为Fine,能带结果导带穿过了费米能级,变成金属了,和没有加U的情况基本一样了,是不是哪里我设置的不对?求指教。如果需要什么图我再补,拜托了

另外,有文章用B3LYP算NiO能带,但用GGA优化结构之后换成杂化泛函算能带提示错误,把USPP之类的按提示修改之后可以运行,但刚运行一会儿就报错,提示ionic mismatch,赝势重启时发生改变之类,这又是哪里出问题了呢?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaomao_16

银虫 (小有名气)

引用回帖:
2楼: Originally posted by 卡开发发 at 2020-03-13 13:26:29
1、带隙偏小可能还是LDA/GGA本身的缺陷,这个没办法,要么DFT+U要么杂化泛函,对MgO而言DFT+U起不到效果,只能考虑杂化泛函。
2、NiO体系的话DFT+U应该是显著改善的,不过得确定你的磁矩是否按照文献正确设定。
3 ...

好的,感谢指导!还有个关于+U的问题,是只需要选中所有Ni原子Modify的时候加上Hubbard U,之后继续普通的GGA计算,就可以得到GGA+U的结果吗,还是要勾选LDA+U或者进行什么设置呢?
3楼2020-03-13 14:11:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★ ★ ★ ★
xiaomao_16(漫天飘雪代发): 金币+5, 谢谢交流 2020-03-13 23:17:12
1、带隙偏小可能还是LDA/GGA本身的缺陷,这个没办法,要么DFT+U要么杂化泛函,对MgO而言DFT+U起不到效果,只能考虑杂化泛函。
2、NiO体系的话DFT+U应该是显著改善的,不过得确定你的磁矩是否按照文献正确设定。
3、没记错的话CASTEP算杂化泛函应该使用模守恒赝势,这是必须的条件;另外不推荐B3LYP算固体,不符合自由电子极限,建议采用HSE06或者sX-LDA,相关资料你可以再看看文献。
不一定挂在论坛,计算问题问题欢迎留言。
2楼2020-03-13 13:26:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
3楼: Originally posted by xiaomao_16 at 2020-03-13 14:11:08
好的,感谢指导!还有个关于+U的问题,是只需要选中所有Ni原子Modify的时候加上Hubbard U,之后继续普通的GGA计算,就可以得到GGA+U的结果吗,还是要勾选LDA+U或者进行什么设置呢?...

需要勾选LDA+U和自旋极化,注意此处LDA+U其实就是DFT+U,并不区分结合使用LDA还是GGA甚至是MetaGGA。
不一定挂在论坛,计算问题问题欢迎留言。
4楼2020-03-13 16:09:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaomao_16

银虫 (小有名气)

引用回帖:
4楼: Originally posted by 卡开发发 at 2020-03-13 16:09:19
需要勾选LDA+U和自旋极化,注意此处LDA+U其实就是DFT+U,并不区分结合使用LDA还是GGA甚至是MetaGGA。...

好的 万分感激

发自小木虫IOS客户端
5楼2020-03-13 16:49:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 TZqyNTBa07NP 2026-11-13 3/150 2026-11-15 14:22 by 0aDbZAghdVq9
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 TZqyNTBa07NP 2026-11-13 3/150 2026-11-15 11:47 by xSrQNSSjcgoX
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 TZqyNTBa07NP 2026-11-13 3/150 2026-11-15 11:03 by xSrQNSSjcgoX
[考博] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 TZqyNTBa07NP 2026-11-13 4/200 2026-11-15 09:57 by xSrQNSSjcgoX
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 SX5mJmPWDfOO 2026-11-13 3/150 2026-11-15 07:56 by xSrQNSSjcgoX
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 SX5mJmPWDfOO 2026-11-13 4/200 2026-11-15 07:45 by xSrQNSSjcgoX
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 XJXMuRjD7cK7 2026-11-13 4/200 2026-11-15 07:12 by xSrQNSSjcgoX
[基金申请] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 XJXMuRjD7cK7 2026-11-13 3/150 2026-11-15 05:44 by xSrQNSSjcgoX
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 uKOSCocCb6Fm 2026-11-13 3/150 2026-11-15 04:05 by xSrQNSSjcgoX
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 C6CRn66RYYXS 2026-11-10 9/450 2026-11-13 11:02 by Equinoxhua
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 cdqGD6ftHmbl 2026-11-11 5/250 2026-11-12 12:28 by ECLTXnJnRSXf
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 czAiolzR0IJC 2026-11-10 8/400 2026-11-12 10:49 by HajghKOMynmz
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 G2rhSvqvhNzl 2026-11-10 6/300 2026-11-12 10:30 by HajghKOMynmz
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 C6CRn66RYYXS 2026-11-10 8/400 2026-11-12 10:27 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 C6CRn66RYYXS 2026-11-10 7/350 2026-11-12 10:16 by HajghKOMynmz
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 C6CRn66RYYXS 2026-11-10 6/300 2026-11-12 10:16 by HajghKOMynmz
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 FChKX8meu1ds 2026-11-09 8/400 2026-11-12 09:46 by HajghKOMynmz
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 gRND4NVAZ8VG 2026-11-11 3/150 2026-11-12 07:42 by HajghKOMynmz
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 Z7qfCsdIiD1T 2026-11-10 4/200 2026-11-12 02:23 by 9fOTqj074tIg
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 cqv1zKnzd7hr 2026-11-09 7/350 2026-11-12 00:44 by 9fOTqj074tIg
信息提示
请填处理意见