24小时热门版块排行榜    

查看: 10063  |  回复: 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的回帖

卡开发发

专家顾问 (著名写手)

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的回帖

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的回帖

卡开发发

专家顾问 (著名写手)

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的回帖

gengxiy

新虫 (初入文坛)

你好,我想问一下,我设置了正负2的磁矩,U值设置的是8,为什么算出来带隙只有零点几

发自小木虫IOS客户端
6楼2026-03-23 00:02:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

awmc2008

至尊木虫 (文坛精英)

学习者

MgO加剪刀差即可。NiO中,O的P轨道也要加U值。
xuexizhe
7楼2026-04-05 11:57:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaomao_16 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026国自然函评费到账 +9 羊腰板 2026-08-21 9/450 2026-08-22 00:03 by anji999
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 14/700 2026-08-21 16:51 by applepetty
[基金申请] 人气不行了 +7 fansofjerry 2026-08-21 7/350 2026-08-21 14:13 by Vivilian
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 什么时候开奖? +7 CrisMessi 2026-08-18 7/350 2026-08-21 12:15 by weiyin
[基金申请] 时间戳又变了8-15 +15 archvillain 2026-08-15 27/1350 2026-08-21 11:32 by tim76
[基金申请] 今天基金会出结果吗?20260819 +15 kkkl_v 2026-08-19 16/800 2026-08-21 11:22 by 365372687
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 时间戳今天,20号变了 +4 archvillain 2026-08-20 4/200 2026-08-21 08:10 by gatelove
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
信息提示
请填处理意见