求助用DFT+U算Y2O3带隙，加完U值带隙比不加U值更小了

bushiyueliang

我在用DFT+U算Y2O3时，U值加的越大，算出的带隙越小，用PBE直接算带隙4.12eV,加完变成3点几了，这是什么情况啊，求大佬看看
            

bushiyueliang

@卡开发发 求助大佬

卡开发发

bushiyueliang 发表于 2024-4-23 09:54
@卡开发发 求助大佬

信息不够充分：比如
1、这里开了ISPIN=2，但是最终总磁矩算出来是多少？是否收敛到合理的结果？
2、尽管开了DFT+U，但实际上+U后电子密度要发生驰豫，结构也会随着结构变形，其实没有原理保证+U绝对能够确保带隙增加。
3、其他参数的合理性，如k点的充分性等等，当然甚至有小概率必须要考虑半芯态结果才会正常。之前甚至遇到产生输入时赝势顺序弄反的。
所以暂时没办法下定论。

bushiyueliang

卡开发发 发表于 2024-4-23 14:40
信息不够充分：比如
1、这里开了ISPIN=2，但是最终总磁矩算出来是多少？是否收敛到合理的结果？
2、尽 ...

老师，全部文件太大了，我下了结构优化、单点能计算、还有能带计算的输入文件和OUTCAR，不知道够不够

卡开发发

bushiyueliang 发表于 2024-4-23 18:59
老师，全部文件太大了，我下了结构优化、单点能计算、还有能带计算的输入文件和OUTCAR，不知道够不够

我测试了一下，看来就是上面提到的第二点。

bushiyueliang

卡开发发 发表于 2024-4-23 20:41
我测试了一下，看来就是上面提到的第二点。

谢谢老师！！！那是说明DFT+U不适合这个体系吗？对于像这种稀土元素原胞原子数较多（40个原子）求能带除了杂化泛函和DFT+u,老师有别的方法推荐吗

卡开发发

bushiyueliang 发表于 2024-4-23 21:24
谢谢老师！！！那是说明DFT+U不适合这个体系吗？对于像这种稀土元素原胞原子数较多（40个原子）求能带除 ...

40原子不算多，杂化泛函其实还可以接受。其他可选的还有mBJ和DFT-1/2，表面上看mBJ更黑箱化，但SCF收敛会让人有点头疼。

bushiyueliang

卡开发发 发表于 2024-4-24 03:55
40原子不算多，杂化泛函其实还可以接受。其他可选的还有mBJ和DFT-1/2，表面上看mBJ更黑箱化，但SCF收敛会 ...

非常谢谢老师！！！我之前有看到mBJ方法，因为主要想求的是介电函数，不确定它对稀土体系是否适用

卡开发发

bushiyueliang 发表于 2024-4-24 08:55
非常谢谢老师！！！我之前有看到mBJ方法，因为主要想求的是介电函数，不确定它对稀土体系是否适用

可以用，效果怎样不确定。

bushiyueliang

卡开发发 发表于 2024-4-24 09:43
可以用，效果怎样不确定。

麻烦老师了！谢谢老师耐心回复！非常感谢老师！