熟悉DFT计算的小伙伴们都已熟知GGA-PBE泛函所计算出的能带带隙往往小于材料实验值,而使用杂化泛函如HSE则可能计算得出相对接近的带隙值,但是这种方法在计算资源上的占用非常巨大,往往需要计算非常久。这里我们介绍一下mBj交换势这种在计算资源上相对较少,速度快且也能得到相对准确带隙的计算方法,当然这种方法的缺点是...
2, MBJ能带:SYSTEM=NbX2-hex ENCUT=600 ISMEAR=0 SIGMA=0.01 NSW=0 IBRION=-1 ISIF=2 NELMD...
可如果用这个泛函得到的结构去做mBJ计算,最终得到的能隙却真是一言难尽,黑磷块体材料说好的0.3 eV的能隙被活生生算没有了(如下图)。测试了一堆泛函之后,孔博士和乔博士发现:如果用optB88-vdW弛豫得到的结构,再用mBJ去算个能隙,这种组合得到的结果(0.31 eV)似乎跟实验符合的挺不错的。这么看来optB88-vdW这...
可以看到单独的GGA-PBE肯定是低估带隙了,并且+U确实有效果,带隙提高到了2.22 eV左右,mBJ对带隙的计算依旧发挥稳定,3.04 eV。还尝试了一下在mBJ的基础上再+U(有点邪门,但是确实可以这么操作),带隙继续提高到了3.44 eV,令人意外的是常用来矫正能量的SCAN和R2SCAN泛函也有一定改善带隙计算的能力,效果和+U接近,...
能带计算中涉及的k点的权重均设置为0 Step 5:最后一步计算1. INCAR、KPOINTS、POSCAR、POTCAR & WAVECAR(Step 4) --> Step 52. 修改INCAR的参数:ISTART = 1 # 读取WAVECARICHGCAR = 0NELM = 400 # MBJ修正的收敛性比较糟糕,所以电子步上限建议设高一些ALGO = All # 该算法计算得偏慢,但是相对来说比...
ALGO = All # 该算法计算得偏慢,但是相对来说比较好收敛 METAGGA = MBJ #在vaspwiki中有其他的选择,比如SCAN等 LASPH = .T. GGA = CAiii) 如果研究体系是体材料,可以不用在INCAR中设置CMBJ的数值;如果研究的的体系是二维材料,建议在INCAR中,将CMBJ的值设置为在对应体系体材料计算时的CMBJ对应的值,该数...
本帖内容被屏蔽
计算COF的能带,原子数70多,实验带隙较高,通过vasp计算带隙较低。我想使用MBJ的方法提高带隙,在网上...
调整INCAR参数:ISTART = 1,读取WAVECAR;ICHGCAR = 0,NELM = 400(电子步上限建议设高一些),ALGO = All(算法计算较慢,但收敛性较好),METAGGA = MBJ(选择其他选项如SCAN等),LASPH = .T.,GGA = CA。对于体材料,无需在INCAR中设置CMBJ值;对于二维材料,建议在INCAR中将CMBJ值设为...
大部分是算能带的,说是先要算普通的PBE,因为mBJ计算需要读取WAVECAR,然后KPOINTS根据杂化算能带的方式...