|
最近在计算单个原子的DOS,为了拉平两个原子的fermi level,我把一个O原子和一个金属原子放在了一个bulk里面,本来的bulk是18*19*20,两个原子之间的距离将近8埃,但是画出来的DOS图会发生劈裂,无论是O还是金属原子,对此我很是疑惑,但认为是可能二者之间的距离还没有很大。于是,我继续增大bulk.增大到24*25*26,发现还是DOS出现劈裂;我继续扩大bulk到40*41*42,上述情况依然出现,于是我想可能无法解决,我就把单个金属原子放在一个大bulk进行计算,用的赝势都是vasp官网推荐使用的赝势,通过看OUTCAR,发现电子占据是正确的,但是画出来的DOS图有劈裂,见下图。最奇怪的现象是我发现,当5个d轨道全部被自旋向下和自旋向上的两个电子占据时,画出来的DOS图是正确的,还有一种情况,就是像Mn,Cr这种5个d轨道被自旋向上电子单占据时画出来的DOS图也是正确的。其他的占据情况画出来的图都不正确,而且有点金属原子连占据都不是正确的,画出来的DOS图更不是正确的。注:上述算单个原子的bulk大小为24*25*26,另外我还计算了单个O原子的DOS,画出来的也是正确的,我就很疑惑,难道算过渡金属单个原子没法算准吗?想请问大家有什么解决办法嘛?是否也遇到过这类情况谢谢大家!
图片放不上来啊!!
|
|