Znn3bq.jpeg
查看: 5130  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

qghrr

铁虫 (初入文坛)


[交流] 【求助】怎样移动态密度图DOS的费米能级Ef才是正确的呢?

问题:为了同其他模型对比,计算出了晶体硅Si在gama点的态密度图,如下图所示,但老师说费米能级是错的。



计算按标准计算DOS步骤进行,经过了自洽计算->非自洽计算->DOSCAR分割,使用vdos脚本分割DOSCAR时已经自动移动了费米能级。

可是经过这样做得到的费米能级似乎是错的,其位置应该是在红色点线附近才对。不知应该怎样移动费米能级才是正确的呢?从图中看应该移动大约2 eV左右,但这个值从哪里得到。

盼高手解答。

附:
1、晶体硅Si的计算模型图



2、VASP计算的INCAR文件、KPOINTS文件。

INCAR文件(非自洽计算)
general:
System = Crystal Si cell
ISTART = 0 ; ICHARG=11 #自洽计算 ICHARG=2; 非自洽计算 ICHARG=11
PREC=Accurate
ENCUT=320
ISMEAR = 0; SIGMA = 0.1
IALGO  =  48
GGA=91
ISPIN=2
VOSKOWN = 1
AMIX=0.2
BMIX=0.0001
AMIX_MAG=0.8
BMIX_MAG=0.0001
MAGMOM= 8*0
LREAL  =  .FALSE. #option for system contain less than 20 atoms
NPAR  =  1
LWAVE  = .F.
LORBIT = 1        #非自洽计算选项,计算spd分波态密度
RWIGS  = 1.312

非自洽计算gama点态密度DOS时的KPOINTS文件:
K-Points
0
Gama
1 1 1
0 0 0

[ Last edited by qghrr on 2011-4-14 at 12:11 ]
回复此楼

» 收录本帖的淘帖专辑推荐

vasp光学

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qghrr

铁虫 (初入文坛)


引用回帖:
Originally posted by ioribest at 2011-04-14 13:39:47:
之前的自洽计算,k点取了多少?

自洽计算时的KPOINTS文件:
K-Points
0
m
11 11 11
0  0  0
5楼2011-04-14 15:01:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答
★ ★
qghrr(金币+10): 2011-04-14 11:48:41
zzy870720z(金币+2): 谢谢指教 2011-04-14 20:20:19
费米能级的移动是什么意思呢?

费米能级不是确定的吗?在自洽计算得到的OUTCAR文件中有的。在非自洽计算得到的OUTCAR文件中也有的。这两处得到的费米能级有些差别,有的似乎较大。

做split dos的时候,会从与DOSCAR放在一起的的OUTCAR文件中读费米能级E-fermi。然后就自动以这个费米能级为零点。

我个人觉得,费米能级算不准,不同软件,以及vasp中静态自洽和非自洽计算结果也有差别。所以有文献干脆不提费米能级,而只以价带顶为参考零点。

你的体系并不大,静态自洽的K点只取一个似乎不好,可以取7*7*7.

费米能级是否能算的准,比如费米能级是否能出现在浅受主能级与价带顶之间,或者本征绝然体带隙正中间?我没有算出来过。希望高手提供帮助。

[ Last edited by cenwanglai on 2011-4-14 at 10:50 ]
2楼2011-04-14 10:48:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qghrr

铁虫 (初入文坛)


引用回帖:
Originally posted by cenwanglai at 2011-04-14 10:48:33:
费米能级的移动是什么意思呢?

费米能级不是确定的吗?在自洽计算得到的OUTCAR文件中有的。在非自洽计算得到的OUTCAR文件中也有的。这两处得到的费米能级有些差别,有的似乎较大。

做split dos的时候,会从 ...

1、费米能级是能够确定的。
费米能级由vdos脚本(或split_dos脚本)自动从DOSCAR中抽取出来的,然后得到上面的DOS图。
2、k点是故意选成单个gama点的。
由于我想看gama点的态密度,所以只设定了一个k点,因而k点个数是不必增加的。
3、费米能级的确会算不准。
之前查看过自洽计算和非自洽计算的OUTCAR,它们的费米能级有差别,因而在作图时应该考虑自洽计算和非自洽计算费米能级的差值。

我的困惑在于,现在图中标示的费米能级是错的,应该是在大约红线的位置,可是不知怎么如何移动到那里,以及移动量的依据。

同盼高手...

[ Last edited by qghrr on 2011-4-14 at 12:09 ]
3楼2011-04-14 12:07:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ioribest

银虫 (小有名气)


qghrr(金币+5): 2011-04-14 15:02:30
之前的自洽计算,k点取了多少?
4楼2011-04-14 13:39:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 qs8kylt1wr 2026-05-01 4/200 2026-05-03 00:00 by 火炎焱燚炙烤
[论文投稿] 有人投过CCC中国控制会议吗? 19+3 初九凡 2026-04-26 4/200 2026-05-02 21:11 by Equinoxhua
[考博] 考博 +4 Eternitf 2026-04-29 4/200 2026-05-02 16:10 by yejenny
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 qs8kylt1wr 2026-05-01 3/150 2026-05-02 15:51 by zahslv263y
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 qs8kylt1wr 2026-05-01 3/150 2026-05-02 15:51 by zahslv263y
[教师之家] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 qs8kylt1wr 2026-05-01 3/150 2026-05-02 15:36 by zahslv263y
[博后之家] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 qs8kylt1wr 2026-05-01 3/150 2026-05-02 15:06 by zahslv263y
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 d8tbypao1q 2026-05-01 3/150 2026-05-02 14:21 by zahslv263y
[博后之家] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 d8tbypao1q 2026-05-01 3/150 2026-05-02 14:06 by zahslv263y
[教师之家] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 d8tbypao1q 2026-05-01 3/150 2026-05-02 14:06 by zahslv263y
[博后之家] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k17zww5nna 2026-05-01 3/150 2026-05-02 10:21 by ovidyohpyz
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k17zww5nna 2026-05-01 4/200 2026-05-02 10:21 by ovidyohpyz
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k17zww5nna 2026-05-01 4/200 2026-05-02 10:06 by ovidyohpyz
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k17zww5nna 2026-05-01 5/250 2026-05-02 09:21 by dz7y1p1i98
[考博] 药化及相关博士的申请 5+3 中药学渣 2026-04-29 3/150 2026-05-01 18:14 by babero
[论文投稿] 一篇MDPI论文改变了学习工作和生活 +3 helsinglee 2026-04-29 4/200 2026-05-01 11:45 by bobvan
[文学芳草园] 太白金星有点烦 +3 yexuqing 2026-04-27 3/150 2026-04-30 21:56 by myrtle
[考博] 中国地质大学(北京)博士招生补录,数理学院材料科学与工程专业和材料与化工专业 +3 liuxh329 2026-04-29 3/150 2026-04-30 15:42 by 顺利毕业zs
[基金申请] 收到国自然专家邀请后几年才会有本子送过来评 +3 barley5 2026-04-30 3/150 2026-04-30 14:58 by aspect3000
[论文投稿] 有没有快的中文核心比较快录用的,纳米材料光催化 5+3 2914439618 2026-04-27 4/200 2026-04-29 08:59 by 北京莱茵润色
信息提示
请填处理意见