咨询一下大体系能带计算的方法选择问题

玉米猫

最近在学习VASP，有一个100个原子的有机体系晶胞想进行能带计算，大家都说算能带用HSE06杂化泛函效果比较好，我用了vaspkit教程里的两种办法（1.pbe-scf后读取波函数进行HSE06能带计算；2.直接用HSE06进行能带计算），都算不动，所以想问一下各位老师这种大体系一般都怎么计算能带。是不是只能用纯泛函？该怎么选择方法呢？

卡开发发

ljh123 发表于 2025-5-22 18:15
老师我还想请教一下您，metaGGA（如mBJ）计算能带的流程是不是和杂化泛函是一样的呢？我看网上的教程，杂 ...

一样的，手动设置k点的做法肯定是最为通用的。

ljh123

卡开发发 发表于 2025-5-21 13:43
没意义，这在vaspwiki的METAGGA关键字处其实有说明：

这里其实讨论到你可以使用local mBJ（LMBJ）的方 ...

老师我还想请教一下您，metaGGA（如mBJ）计算能带的流程是不是和杂化泛函是一样的呢？我看网上的教程，杂化泛函相比纯泛函的计算方法就是多考虑了权重为0的K点。似乎metaGGA相比纯泛函也是多考虑了权重为0的K点，那是否可以按照HSE06的能带计算方法来算呢

ljh123

卡开发发 发表于 2025-5-21 13:43
没意义，这在vaspwiki的METAGGA关键字处其实有说明：

这里其实讨论到你可以使用local mBJ（LMBJ）的方 ...

十分感谢。我去看看LMBJ方法

卡开发发

ljh123 发表于 2025-5-21 13:06
是的。我目前的情况是，不主动设置CMBJ电子步可以正常收敛，收敛后的CMBJ值为2.0。我计算得到的Fe3O4@C的 ...

没意义，这在vaspwiki的METAGGA关键字处其实有说明：

In the presence of an extended vacuum region (e.g., surfaces) or an interface, the average of |▽ρ|/ρ has no meaning. Therefore, MBJ calculations should be done with a fixed value of c}, which can be done with the CMBJ tag., or alternatively with the LMBJ that was proposed for the purpose to be applicable to systems with vacuum or interfaces.

这里其实讨论到你可以使用local mBJ（LMBJ）的方法。不过我没有去读过LMBJ的原文，不太确定它的具体原理和适用范围。

ljh123

卡开发发 发表于 2025-5-21 09:15
且不说CMBJ调整，一般mBJ方法收敛都很困难。

是的。我目前的情况是，不主动设置CMBJ电子步可以正常收敛，收敛后的CMBJ值为2.0。我计算得到的Fe3O4@C的带隙在0.1eV左右。但是如果主动设置CMBJ=1.35的话一直不收敛（这个值是通过bulk计算收敛得到的）。所以我想请教下老师，“主动设置CMBJ电子步可以正常收敛，收敛后的CMBJ值为2.0“这个情况下，计算是否还有意义呢。因为2.0这个值似乎是CMBJ参数的上限值。

卡开发发

ljh123 发表于 2025-5-21 09:14
是的，不过我尝试强行调整后，收敛就变得极其困难。

且不说CMBJ调整，一般mBJ方法收敛都很困难。

ljh123

卡开发发 发表于 2025-5-20 21:02
我没记错对slab可以强行调整CMBJ，但这个意义我也很难说。

是的，不过我尝试强行调整后，收敛就变得极其困难。

卡开发发

ljh123 发表于 2025-5-20 16:50
老师您好，我之前尝试过mBJ方法用于slab。电子步能顺利收敛。但是确实发现CMBJ这个参数会从一开始就变为2 ...

我没记错对slab可以强行调整CMBJ，但这个意义我也很难说。

ljh123

卡开发发 发表于 2025-5-9 18:02
mBJ不能用于slab，数值稳定性存在问题。DFT-1/2原理上没有限制，但有可能需要调参。

老师您好，我之前尝试过mBJ方法用于slab。电子步能顺利收敛。但是确实发现CMBJ这个参数会从一开始就变为2.0并一直不变。请问这样的结果还有意义吗

卡开发发

ljh123 发表于 2025-5-9 14:09
老师您好，请问slab体系还能用DFT-1/2或者mBJ算电子结构吗

mBJ不能用于slab，数值稳定性存在问题。DFT-1/2原理上没有限制，但有可能需要调参。

ljh123

卡开发发 发表于 2022-2-1 17:28
如果是纯有机体系DFT+U可能没用。bulk结构可以考虑使用DFT-1/2方法或者modified Becke-Johnson（mBJ或叫做T ...

老师您好，请问slab体系还能用DFT-1/2或者mBJ算电子结构吗

heguangyuan

876449830 发表于 2022-9-4 21:36
我在读取WAVECAR和CHAGCAR和只读取CHAGCAR计算的结果不一样，后来听人建议只读取CHAGCAR会好些，至今不知 ...

因为计算能带的时候用的是新生成的K点，和之前生成wavecar的K点不同，读取了应该不太对

876449830

乐平 发表于 2022-9-4 17:15
个人觉得，“先用 PBE 泛函进行自洽计算，然后读取自洽计算的波函数和电荷密度再用 HSE06 杂化进行非自洽 ...

我在读取WAVECAR和CHAGCAR和只读取CHAGCAR计算的结果不一样，后来听人建议只读取CHAGCAR会好些，至今不知为何？

乐平

石不转 发表于 2022-5-30 09:04
各位老师好！有一个挺基础的问题求教： vaspkit教程 用PBE泛函计算能带为何要先自洽再进行非自洽运算？与直 ...

个人觉得，“先用 PBE 泛函进行自洽计算，然后读取自洽计算的波函数和电荷密度再用 HSE06 杂化进行非自洽计算”，这个方法类似于 Gaussian 里的先用 B3LYP泛函配合 6-31G(d) 基组进行几何优化，然后再换 meta-GGA 结合大基组，读取之前的 check 文件进行后续计算类似。

由于已经经过了自洽计算（往往这一步是比较耗时的），从自洽计算的波函数和电荷密度中读取信息一般会让后续更快。