应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】怎样移动态密度图DOS的费米能级Ef才是正确的呢？
北京石油化工学院2026年研究生招生接收调剂公告
71/1返回列表
查看: 5052  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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光学

» 猜你喜欢

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

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

查看全部散金贴
1楼 2011-04-14 10:25:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)
★ ★
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 ]
赞一下(1人) 回复此楼
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的回帖

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的回帖

菲顶顶

铜虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
3楼: Originally posted by qghrr at 2011-04-14 12:07:06
1、费米能级是能够确定的。
费米能级由vdos脚本（或split_dos脚本）自动从DOSCAR中抽取出来的，然后得到上面的DOS图。
2、k点是故意选成单个gama点的。
由于我想看gama点的态密度，所以只设定了一个k点，因而 ...

楼主的困惑解决了嘛？我也遇到了同样的问题
赞一下 回复此楼
6楼2014-09-14 17:19:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paopaotu326

金虫 (正式写手)
★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
franch: 金币+2, 谢谢回帖交流。。。 2014-09-30 17:00:42
不知道我的理解对不对，经过splidos的时候，程序已经自动把Ef放在0点了，没有不对啊？为什么要移动？你可以问下老师，不要怕，一切反动派都是纸老虎，嘿嘿：）
赞一下 回复此楼
7楼2014-09-26 09:46:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 qghrr 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学调剂 +13 艾志恒 2026-04-03 14/700 2026-04-05 10:59 by 风雨无晴
[考研] 材料调剂 +9 革微桂 2026-04-04 9/450 2026-04-05 08:27 by 544594351
[考研] 材料334求调剂 +15 Eecho# 2026-04-03 15/750 2026-04-04 23:05 by 无际的草原
[考研] 本科211 分数293请求调剂 +4 莲菜就是藕吧 2026-04-01 4/200 2026-04-04 22:32 by hemengdong
[考研] 277求调剂 +4 12A3 2026-04-02 5/250 2026-04-04 20:28 by 蓝云思雨
[考研] 340求调剂 +4 jhx777 2026-03-29 4/200 2026-04-04 20:08 by 无际的草原
[考研] 278求调剂 +6 Yy7400 2026-04-03 6/300 2026-04-04 09:53 by zhangdingwa
[考研] 考研调剂 +5 小sun要好运 2026-04-03 5/250 2026-04-03 21:43 by 啵啵啵0119
[考研] 266求调剂 +3 08电气工程 2026-04-03 3/150 2026-04-03 14:05 by 1753564080
[硕博家园] 求老师收留 +9 lllq123 2026-04-03 9/450 2026-04-03 13:48 by 呼吸都是减肥
[考研] 数一英一285求调剂 +7 AZMK 2026-04-03 9/450 2026-04-03 13:03 by ms629
[考研] 考研调剂 +3 李木子0120 2026-04-02 5/250 2026-04-02 21:45 by dongzh2009
[考研] 一志愿华南师范大学-22408计算机-292分-求华南师范大学调剂 +4 爱读书的小鳄鱼 2026-04-02 4/200 2026-04-02 18:35 by 求调剂zz
[考研] 一志愿上海海洋大学083200食品学硕，求调剂，接受其他专业 +6 what张 2026-04-01 7/350 2026-04-02 16:48 by zzsw+
[考研] 354求调剂 +4 lxb598 2026-03-31 5/250 2026-04-02 09:55 by Jaylen.
[考研] 考研调剂 +12 Amber00 2026-03-31 12/600 2026-04-02 09:04 by sanrepian
[考研] 一志愿北交材料工程总分358 +5 cs0106 2026-04-01 7/350 2026-04-01 11:45 by wangjy2002
[考研] 一志愿西交大080500材料学硕349 +6 jqx1258 2026-03-31 7/350 2026-03-31 21:08 by yuq
[考研] 085404 22408 315分 +5 zhuangyan123 2026-03-31 6/300 2026-03-31 13:48 by limeifeng
[考研] 085602 化学工程专硕 340分求调剂 +4 qianbai11 2026-03-29 4/200 2026-03-30 11:34 by 唐沐儿
信息提示
关闭
请填处理意见
关闭