DOS图为啥特别密集和离散

qwe1832518773

按照网上说的拷贝静态计算的整个文件夹，然后做dos计算，同时保证dos计算和静态计算级别一致，但是为啥算完的dos图这么丑。。是展宽的原因吗？但是我用的是四面体方法（ISMEAR=-5），SIGMA的设置应该不影响的吧。dalao们能帮我看一下是出了什么问题吗？

静态计算INCAR：
SYSTEM = sp
ENCUT = 400
ISTART =0; ICHARG =2
ISMEAR = 0
LREAL = Auto
ALGO = Fast
IBRION = -1
NSW = 0
POTIM = 0.4
ISIF = 2
PREC = accurate
EDIFF = 1.E-6
EDIFFG = -0.002
NELMIN = 4
NPAR = 4
LVDW = .TRUE.
LWAVE = .TRUE.
LCHARG = .TRUE.


DOS的INCAR：
SYSTEM = DOS
ENCUT = 400
ISTART =1; ICHARG =11
ISMEAR = -5
LREAL = Auto
ALGO = Fast
ISIF = 2
PREC = accurate
EDIFF = 1.E-6
NELMIN = 4
LVDW = .TRUE.
SIGMA = 0.1
EMIN = -30
EMAX = 15
LORBIT = 11
NEDOS = 4500
LWAVE = .FALSE.

Penson

如果是含O的体系，ENCUT有点小，最好设到500以上。此外，设置EMIN和EMAX是为了什么？

卡开发发

非自洽做一下再分析，然后sigma设置大一些，sigma有影响，太小的话峰会不平滑。
或者你可以自己写程序进行展宽。

qwe1832518773

Penson 发表于 2021-1-29 09:25
如果是含O的体系，ENCUT有点小，最好设到500以上。此外，设置EMIN和EMAX是为了什么？

感谢回复呢，那我加大ENCUT试一下，EMIN和EMAX是因为我除了这个体系外还做了几个别的体系，想控制DOS的范围一致，这样好比较一些

qwe1832518773

卡开发发 发表于 2021-1-30 15:29
非自洽做一下再分析，然后sigma设置大一些，sigma有影响，太小的话峰会不平滑。
或者你可以自己写程序进行 ...

感谢卡老师，每次都会麻烦您给我回帖，真的十分感谢！
我看了网上说ISEAR=-5的话SIGMA就不会去读取，所以一直没管它，那我加大到0.5试一下！

卡开发发

qwe1832518773 发表于 2021-1-31 20:49
感谢卡老师，每次都会麻烦您给我回帖，真的十分感谢！
我看了网上说ISEAR=-5的话SIGMA就不会去读取，所 ...

非自洽不要ismear=-5，取个0，给个sigma。

qwe1832518773

卡开发发 发表于 2021-1-31 21:19
非自洽不要ismear=-5，取个0，给个sigma。

嗯呢好的！谢谢卡老师

qwe1832518773

卡开发发 发表于 2021-1-31 21:19
非自洽不要ismear=-5，取个0，给个sigma。

卡老师，忘了说了我算的是alpha-Fe2O3，后续加了展宽，并考虑铁磁性和使用半芯态的赝势，但是结果还是很乱，并且和文献里的DOS对不上，文献里的费米能级处是没有态密度的，而我的费米能级处态密度是最高的，这是为啥？

卡开发发

qwe1832518773 发表于 2021-2-3 13:18
卡老师，忘了说了我算的是alpha-Fe2O3，后续加了展宽，并考虑铁磁性和使用半芯态的赝势，但是结果还是很 ...

a-Fe2O3是反铁磁的，需要手工指定磁矩。
另外DFT+U以及表面态各方面都有影响，我不确定你具体的体系是什么。