|
关于带电缺陷的缺陷形成能,我查了蛮多文献依然很难完全理解计算过程,主要问题出在q(EVBM+EF)这一项中
问题一是带电量q是如何取值? 例如我计算CaCO3中一个二价钙离子缺陷,是否计算缺陷的时候设置NELECT=NELECT+2然后计算时设q=2?因为体系多了两个电子?
问题二是公式中的EVBM具体如何得到? 是否需要计算功函数,然后参考真空能级计算?,或者计算dos用pymatgen读取这个vbm?
问题三是公式中的Ef到底怎么算呢?文献中说EFis the Fermi level, referenced to the valence-band maximum in the bulk,费米能级是参考价带的?vasp中使用的grep "E-fermi" OUTCAR 是参考价带的能量吗?因为我用pymatgen得到的vbm和Efermi相差无几,所以感觉不太对,想知道OUTCAR里面的fermi能级是相对于谁的?如果不对应该怎样确定这个Efermi呢?
最后一个问题是化学势如何确定,我知道化学势是不同条件下不同的,算缺陷形成能时怎么确定呢?如计算碳酸钙钙缺陷时可以使用Materials Project找到Ca的最稳定单胞的每原子能量作为化学势吗?(根据我查找的结果是-2eV)? |
|