| 查看: 4905 | 回复: 6 | ||||
[交流]
【求助】怎样移动态密度图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光学 |
» 猜你喜欢
求个博导看看
已经有17人回复
-
青基代表作，AAAI之类的A会的special track在国内认可度高吗？还是归为workshop之流？
已经有3人回复
-
上海工程技术大学【激光智能制造】课题组招收硕士
已经有6人回复
-
带资进组求博导收留
已经有11人回复
-
自荐读博
已经有5人回复
-
上海工程技术大学张培磊教授团队招收博士生
已经有4人回复
-
求助院士们，这个如何合成呀
已经有4人回复
-
临港实验室与上科大联培博士招生1名
已经有9人回复
-
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助关于荧光共振能量转移的问题？
已经有8人回复
- 怎样从甲醇中提取出乙醛肟！
已经有4人回复
- 核磁碳谱中的氯的化学位移
已经有6人回复
- 紧急求助：H谱化学位移一样、峰面积不同，是否说明是不同混合物？
已经有3人回复
- 新手问一个算态密度的问题，谢谢！
已经有3人回复
- 怎么根据态密度图积分算区间的电子数
已经有8人回复
- 恒电流充放电时，为什么电流密度越大测得的容量越小
已经有19人回复
- 费米能级的状态密度低说明什么？（具体见内容）
已经有4人回复
- CASTEP 中DOS 的读法
已经有15人回复
- 电子密度图 标原子
已经有17人回复
- 3.22文献解读：三元材料在高能量密度锂离子电池中的应用
已经有236人回复
- 求低温等离子体电子密度测量方法
已经有4人回复
- DOS图中怎么判断非金属性，半金属性，金属性
已经有7人回复
- CASTEP几何优化，能带结构和态密度分析
已经有11人回复
- 等效位置的态密度分析
已经有9人回复
- 【求助】VASP怎样精确计算出费米能级处的态密度值？
已经有22人回复
» 抢金币啦！回帖就可以得到:
-
南京都市圈高校大龄离异博士征友
+2/520
-
双一流南京医科大学招计算机、AI、统计、生物信息等方向26年9月入学博士
+1/175
-
诚征另一半
+1/173
-
江苏科技大学能源材料化学课题组张俊豪教授招收博士研究生1-2名
+1/81
-
广州，真诚男征女
+1/77
-
广州
+1/69
-
最新看到一个观点：说高校教师的斩杀线是青基和面上
+1/65
-
坐标北京不异地
+1/63
-
南科大薛亚辉课题组诚聘离子输运、低维器件、原子力显微镜等方向“快响行动”博士生
+1/61
-
2026博士申请——有机化学\计算化学\药物化学方向
+1/52
-
上海大学昝鹏教授、军事医学研究院伯晓晨研究员/倪铭副研究员 课题组招聘博士生
+2/46
-
暨南大学理工学院 光子技术研究院段宣明团队申请制读博招生
+1/29
-
盐湖所镁基储氢材料课题组招聘
+1/26
-
天津大学化学系吴立朋课题组申请考核制博士招生/博后招聘-有机化学，金属有机
+1/10
-
天津大学化学系吴立朋课题组申请考核制博士招生/博后招聘
+1/7
-
求博导收留
+1/6
-
中国科学院大连化学物理研究所DNL0902研究组招聘博士后和职工
+1/5
-
酰胺水解求助
+1/3
-
探究TGF-β在癌症免疫调控中的作用机制|肿瘤
+1/1
-
博士后网站系统登录
+1/1
★ ★
qghrr(金币+10): 2011-04-14 11:48:41
zzy870720z(金币+2): 谢谢指教 2011-04-14 20:20:19
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
|
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
4楼2011-04-14 13:39:47
5楼2011-04-14 15:01:53
6楼2014-09-14 17:19:49
7楼2014-09-26 09:46:37