计算化学公社

标题: DOS图为啥特别密集和离散 [打印本页]

作者
Author: qwe1832518773 时间: 2021-1-28 23:19
标题: DOS图为啥特别密集和离散
按照网上说的拷贝静态计算的整个文件夹，然后做dos计算，同时保证dos计算和静态计算级别一致，但是为啥算完的dos图这么丑。。是展宽的原因吗？但是我用的是四面体方法（ISMEAR=-5），SIGMA的设置应该不影响的吧。dalao们能帮我看一下是出了什么问题吗？
(, 下载次数 Times of downloads: 33)

静态计算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.

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

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

作者
Author: qwe1832518773 时间: 2021-1-31 20:48

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

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

作者
Author: qwe1832518773 时间: 2021-1-31 20:49

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

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

作者
Author: 卡开发发 时间: 2021-1-31 21:19

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

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

作者
Author: qwe1832518773 时间: 2021-1-31 21:47

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

嗯呢好的！谢谢卡老师

作者
Author: qwe1832518773 时间: 2021-2-3 13:18

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

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

作者
Author: 卡开发发 时间: 2021-2-3 13:34

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)