哪位用 castep 做过DFT+U 或者 杂化泛函能带计算？

wmg166

                请教大师，我用 castep 做过DFT+U 或者 杂化泛函能带计算，我的64G 内存，用 HSE06 计算能带，截断能仅300  就提示内存不足，怎么办？  
                 我的体系是硅酸盐里面掺杂了一个 Ce，即在 1*2*2 的超胞里面掺杂了一个Ce，共 152 个原子。    DFT+U 是不是搞不定能带？ 我不晓得怎样确定 U值，求大师指点。

卡开发发

1、硅酸盐要看贡献价带顶和导带底的是什么原子，如果是主族元素大概率DFT+U解决不了什么问题。
2、152个原子算杂化泛函64G内存基本上不大可能，杂化泛函要昂贵很多。这样的体系目前来说也没啥特别好的办法。

wolfli369

卡开发发 发表于 2021-11-27 02:05
1、硅酸盐要看贡献价带顶和导带底的是什么原子，如果是主族元素大概率DFT+U解决不了什么问题。
2、152个原 ...

请问怎么确定是价带顶和导带底的贡献，有什么简单的方法？提前谢谢

卡开发发

wolfli369 发表于 2021-11-27 09:50
请问怎么确定是价带顶和导带底的贡献，有什么简单的方法？提前谢谢

看PDOS。

wolfli369

卡开发发 发表于 2021-11-27 11:27
看PDOS。

收到，谢谢

wmg166

卡开发发 发表于 2021-11-27 02:05
1、硅酸盐要看贡献价带顶和导带底的是什么原子，如果是主族元素大概率DFT+U解决不了什么问题。
2、152个原 ...

再请教一下： 局域密度近似（LDA）或广义梯度近似（GGA）往往会低估带隙值，那么如何使理论计算的带隙值与实验值一致呢？  杂化泛函需要极大的内存，很耗时。  一般怎样计算能带准确些？

卡开发发

wmg166 发表于 2021-11-29 11:31
再请教一下： 局域密度近似（LDA）或广义梯度近似（GGA）往往会低估带隙值，那么如何使理论计算的带隙值 ...

除了杂化泛函，可以考虑：
1、非全满d轨道的过渡金属且d电子作为带边的情况可以考虑DFT+U。
2、非表面体系可以考虑mBJ(TB09)，但是SCF极为难收敛，比较考验技巧。
3、DFT-1/2（或者叫LDA-1/2和GGA-1/2），但是需要做参数测试。

北农路2号

卡开发发 发表于 2021-11-29 15:26
除了杂化泛函，可以考虑：
1、非全满d轨道的过渡金属且d电子作为带边的情况可以考虑DFT+U。
2、非表面 ...

老师您好，请问我计算半导体TiO2可以使用DFT+U吗

北农路2号

卡开发发 发表于 2021-11-29 15:26
除了杂化泛函，可以考虑：
1、非全满d轨道的过渡金属且d电子作为带边的情况可以考虑DFT+U。
2、非表面 ...

还有就是U的值该如何确定呢？

卡开发发

北农路2号 发表于 2021-12-6 09:23
还有就是U的值该如何确定呢？

1、可以DFT+U。
2、U值确定方法可以是拟合实验数据 、高级方法、基于线性响应或者通过cRPA计算。

NewYaren

卡开发发 发表于 2021-12-6 10:48
1、可以DFT+U。
2、U值确定方法可以是拟合实验数据 、高级方法、基于线性响应或者通过cRPA计算。

请问含有稀土La时需要加U计算吗？
我有一个含La的配位聚合物，用GGA不加U计算能带结构和DOS时，计算任务总会无缘无故的停掉，把金属换成Cd时就没有事。
另外还有一个问题想请教老师，如何在计算晶体的能带结构时让有机配体多代一个负电荷使结构呈电中性？

NewYaren

请问楼主问题解决了吗？我近期遇到了类似的wen题，计算含有La的配位聚合物的带隙时，计算任务总会无缘无故停掉，很郁闷。

卡开发发

NewYaren 发表于 2021-12-7 23:07
请问含有稀土La时需要加U计算吗？
我有一个含La的配位聚合物，用GGA不加U计算能带结构和DOS时，计算任务 ...

1、需要，但不+U计算无故停止的原因我不清楚，你可能需要提供更多信息。
2、如果体系本身已经是电中性的情况，应该不必做其他调整，自洽场过程会按照能量更低的方式分布电子。反之，如果体系不是电中性，没有合适的修正可能会得到不太正确的结果。

北农路2号

卡开发发 发表于 2021-12-6 10:48
1、可以DFT+U。
2、U值确定方法可以是拟合实验数据 、高级方法、基于线性响应或者通过cRPA计算。

老师您好，还有一个问题想请教您，就是在castep里具体实施加U操作时，在什么时候+U呢？是在收敛性测试后，几何优化和性质计算时+U，还是只在性质计算时加上，还是从一开始的收敛性测试就要+U呢？希望老师能解答

NewYaren

卡开发发 发表于 2021-12-7 23:55
1、需要，但不+U计算无故停止的原因我不清楚，你可能需要提供更多信息。
2、如果体系本身已经是电中性的 ...

好的，谢谢老师。
从晶体结构看，化合物是带一个正电荷的。但从实验证据看有机配体呈还原态，应该多带一个负电荷从而使整体呈电中性。现在不知道怎么给有机分子添加一个额外的负电荷。这个负电荷应该是离域到整个有机分子上的。请问您有相关的经验吗？