[交流] 求助：vasp计算rutile TiO2的电子DOS时，为何单胞和supercell的结果不一样？ 已有2人参与

我一般都是用单胞来计算电子dos。这次因为想看一个特定的晶格distortion对rutile的影响，而这个distortion需要在一个2*2*2的supercell实现，所以我就尝试用一个supercell来算dos了。 优化的单胞展开成2*2*2的supercell（这个supercell没有任何distortion，我是想看一看单胞和8倍的supercell得到的东西是不是一样），开始计算，想着结果应该和单胞出来的没有什么区别啊（我觉得dos的值应该大约是8倍，除此之外分布应该没啥区别），但是奇怪的是supercell产生的dos很明显的少了一个高能级的conduction band peak，我不知道是咋回事。比较单胞和supercell的区别，只有kpoint的密度低了一倍，这应该是没有问题的，以下是我对这个超胞的计算文件 （1）先以4*4*4的 K point 优化结构，进而得到chgcar INCAR: SYSTEM = rutile Relax ISTART = 0 ICHARG = 2 ENCUT  = 500 PREC   = Accurate NELMIN = 4           #  do a minimum of four electronic steps NELM   = 40          #  maximum of 40 electronic steps EDIFF  = 1E-8 EDIFFG = -1E-6        #  high accuracy ADDGRID = .TRUE. IBRION =  -1 ISMEAR  = 0 SIGMA = 0.05 LREAL  = .False. LWAVE  = .False. LPLANE = .TRUE.    NPAR   = 8       # twice the number of nodes ------------------------------------------------------------------------------ (2) 拷贝chgcar到新directory，增加KPOINT 到 11*11*11 INCAR SYSTEM = rutile Relax ISTART = 0 ICHARG = 11 ENCUT  = 500 PREC   = Accurate NELMIN = 4           #  do a minimum of four electronic steps NELM   = 40          #  maximum of 40 electronic steps EDIFF  = 1E-8 EDIFFG = -1E-6        #  high accuracy ADDGRID = .TRUE. IBRION =  -1 ISMEAR  = 0 SIGMA = 0.05 LORBIT = 10 #EMIN = -15 #EMAX = 15 #NEDOS = 1000 LREAL  = .False. LWAVE  = .False. LPLANE = .TRUE.       NPAR   = 8       # twice the number of nodes ----------------------------------------------------------------------- 单胞的计算和上面一样，除了KPOINT 分别是8*8*8 和 21*21*21，我也改过ISMEAR 到-1，木有啥变化。两个系统内力都收敛的很小，是优化的很好的结构。二者的能量是7.999倍，非常好的8倍关系。附件中给出了元胞和supercell的比较，很明显，最上面的一个peak在超胞计算中不见了 本人菜鸟，实在是想不出来什么地方不对了，望高人指点，谢谢啦! dos_rutile.png

回复此楼