计算化学公社

 找回密码 Forget password
 注册 Register
Views: 2070|回复 Reply: 1
打印 Print 上一主题 Last thread 下一主题 Next thread

[VASP] 请问大家做DOS图时sigma设置多少适合

[复制链接 Copy URL]

559

帖子

0

威望

6208

eV
积分
6767

Level 6 (一方通行)

做孤立分子的态密度图时,k点取gamma点,ismear=0,如果sigma设置的太小(如0.05)则会导致DOS峰不圆滑太过尖锐,甚至导致该有的一个分子轨道处出现两个DOS峰,如果sigma设置大些(如0.2)会好很多,请问这是为什么?sigma取多大合适?

0d05.jpg (56.23 KB, 下载次数 Times of downloads: 8)

SIGMA=0.05

SIGMA=0.05

0d2.jpg (34.91 KB, 下载次数 Times of downloads: 11)

SIGMA=0.2

SIGMA=0.2

3809

帖子

3

威望

1万

eV
积分
20344

Level 6 (一方通行)

围观吃瓜群众

2#
发表于 Post on 2024-8-28 22:17:59 | 只看该作者 Only view this author
由于又有人在其他帖子问了一遍相同的问题,我觉得还是应当加以重视,实际上我们可以去看OUTCAR的Bands方面的信息,因为对分子而言盒子要取足够大这时候k点仅需要Gamma点即可,展宽sigma需要足够小。

然后接下来是如下问题的讨论:
1、DOS峰太尖锐,这是展宽足够小导致的,这也是合理的。
2、某位置出现相近的两个DOS峰,这可能是GGA低估gap所导致的。
除了这些还有可能(目前能想到):
3、周期边界下破坏了平移对称性且周期镜像的分子之间会有不应出现的偶极-偶极作用。
4、赝势精度不够,例如要考虑硬赝势(尾缀_h)等。
5、个别色散修正导致结构变化反过来影响电子结构。

排查上面这几个问题最简单的办法就是用其他非PBC的程序一模一样算一遍,然后这里用DMol3(PBE/DNP)和Gaussian(PBE/cc-pVTZ,较大的基组因为常用基组根本就不是在PBE这个交换关联泛函下优化的)都测试了一下,对alpha spin的gap几乎都是0.01 a.u.左右,与VASP当前给出的结果是吻合的。

峰比较窄又比较近的情况需要非常小心,态密度需要足够高的分辨率(态密度x的点数目,亦即NEDOS,以及足够小的sigma),否则有可能有些峰显示不出或者靠近的两个峰混叠在一起,这甚至可能会导致成文的时候产生错误的结论。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-21 23:51 , Processed in 0.182631 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list