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

玉米猫

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

granvia

PBE+U凑活一下

玉米猫

granvia 发表于 2022-2-1 14:11
PBE+U凑活一下

谢谢，纯泛函是不是会显著低估带隙

卡开发发

如果是纯有机体系DFT+U可能没用。bulk结构可以考虑使用DFT-1/2方法或者modified Becke-Johnson（mBJ或叫做TB09）的泛函算电子结构。

玉米猫

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

谢谢老师，我去找一找相关资料

KevinLee

大体系HSE计算好像要国内的PWmat才能搞定，U值有个线性响应算法可以确定

ycpan

可以考虑一下PWmat软件，算100个原子的HSE还是比较轻松的，并且速度会比vasp快很多

get-it

mBJ就挺好的

石不转

各位老师好！有一个挺基础的问题求教： vaspkit教程 用PBE泛函计算能带为何要先自洽再进行非自洽运算？与直接进行非自洽运算有何不同？
我针对同一体系进行了两次试验：
1）直接用vaspkit 303生成的KPATH.in作为KPOINTS进行非自洽运算；
2）按照vaspkit教程中的步骤，先用vaspkit 102生成的KPOINTS进行自洽运算，然后设置ICHARG=11，复制KPATH.in到KPOINTS进行非自洽运算。
试验结果：方法1计算时间约为方法2的两倍；二者画出的能带图几乎完全相同。
结论：对于未知构型，先用简单K点进行自洽运算，在此基础上更换精度更高的K点、读取已存在的波函数来计算能带，如此更方便快捷。
请问这个结论是否正确？如果正确，方法2更快捷的原因是什么？恳请老师们为我答疑解惑，谢谢！

乐平

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

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

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

876449830

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

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

heguangyuan

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

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