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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】VASP计算电子态密度
北京石油化工学院2026年研究生招生接收调剂公告
351/4返回列表1234下一页
查看: 8259  |  回复: 34
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

zzy870720z

荣誉版主 (文坛精英)

优秀版主优秀版主优秀版主优秀版主

[交流] 【求助】VASP计算电子态密度 已有14人参与

最近做点测试计算
可是DOS做出来效果很差
如下图


同时附上自己的INCAR
希望高手帮忙，如何做出漂亮的DOS图
并指出那些参数对DOS图有影响
先谢谢了


  ENCUT  = 300 eV
   ISTART = 0
   ISPIN  = 2
   ISYM   = 0
   GGA    = 91
   VOSKOWN = 1
   ISIF   = 0
   NELM   = 100
   ALGO   = FAST
   EDIFF  = 1E-4
   EDIFFG = 0.001
   LREAL  = .FALSE.
   NSW    = 0
   IBRION = -1
   ISMEAR = -5
   SIGMA  = 0.1
   NWRITE = 2
   NKPR   = 1
回复此楼

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

VASP 第一性原理 vasp 学习贴 模拟，理论
VASP 待办事宜 VASP报错建模问题集

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
博学、审问、慎思、明辨、笃学
1楼 2010-06-02 23:06:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

valenhou001

至尊木虫 (职业作家)
★ ★
zzy870720z(金币+2):你的意思是采用gaussian smearing？但是vasp手册上建议态密度用四面体方法取-5？ 2010-06-03 07:54:04
hometownlove(金币+2):谢谢交流 2010-06-09 13:40:08
1)增加k点网格的大小；
或者
2）采用ISMEAR=0, SIGMA = 0.1~0.2。
赞一下(3人) 回复此楼

» 本帖已获得的红花（最新10朵）

xiaobo1314
2楼2010-06-03 07:35:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)
★ ★ ★ ★
gavinliu7390(金币+3):非常感谢详细指导！ 2010-06-03 14:57:20
zzy870720z(金币+5):谢谢侯老师指点，我先试试，有问题再请教 2010-06-03 16:44:47
hometownlove(金币+1):谢谢交流 2010-06-09 13:44:57
下面的资料是从The FHIMD Toolkit User's Manual上拷贝下来的：
http://www.fhi-berlin.mpg.de/th/fhi98md/

见里面的公式3.2

1、先用vasp计算态密度，比如采用ISMEAR=-5，同时将NEDOS设置的大一些，默认值为301，将它调大，比如设置为1000左右；
2、得到了DOSCAR文件之后，取出total DOS. 可以简单的认为，在采用ISMEAR=-5时计算态密度，上面公式3.1中的F就是delta函数。
3、按上面公式3.1， 3.2，对计算出来的态密度做一个展宽，这个要自己编程：
公式3.1中的权重都为1.0了，不对k进行循环，只是对能量点的个数， 公式3.1里面|...|内含的量已经由vasp计算出来了的。

[ Last edited by valenhou001 on 2010-6-3 at 13:07 ]
赞一下(4人) 回复此楼
13楼2010-06-03 13:06:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主
★
zzy870720z(金币+2):感谢交流 2010-06-03 16:45:06
hometownlove(金币+1):谢谢交流 2010-06-09 13:45:52
引用回帖:
Originally posted by vasp001 at 2010-06-03 12:50:04:


展宽，是在你处理数据的脚本中来约定的，具体的细节你要看你的脚本，一般来说这个参数比约定SIGMA 对态密度的影响更大。如果展宽取得大一些的话，峰行不会那么尖锐，并且，态会相对弥散一些，我觉着这样才更符 ...

展宽不就是SIGMA吗?在INCAR中不是已经给出来了?为什么会有""在处理数据的脚本中来约定"这一说呢?

搜索到下面这一段话:"当体系的电子结构在Fermi 能级附近出现多个简并能级时，体系的自洽过程会相当慢，甚至不能收敛。这是由于能级简并，或能级之间相差不大，导致电子“不知道”应该在哪里分布，从而不断地在几个能级之间“跳来跳去”，但无论电子分布在哪个前线能级上整个体系都不能自洽。
所谓拖尾，是指在前线轨道上，电子可以以分数分布，而不再要求一定要以整数分布，从而使电子可以一分为多地分布在能量相近的各个前线轨道，而不需要在各轨道之间“跳来跳去”，使体系更容易自洽。使用该办法时要对拖尾效应设定一个能量值，使电子在能量相差在该值之内的轨道进行分数分布。设定值越小，计算结果越精确，但亦更难自洽。（也叫展宽）"
赞一下(3人) 回复此楼
14楼2010-06-03 15:08:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vasp001

金虫 (正式写手)
★
zzy870720z(金币+2):谢谢指导，能说更具体点吗？作图时有什么要注意的？ 2010-06-03 08:24:14
hometownlove(金币+1):谢谢交流 2010-06-09 13:41:21
引用回帖:
Originally posted by zzy870720z at 2010-06-02 23:06:33:
最近做点测试计算
可是DOS做出来效果很差
如下图


同时附上自己的INCAR
希望高手帮忙，如何做出漂亮的DOS图
并指出那些参数对DOS图有 ...

换用ISMEAR =-5 ，作图时，还有一个参数是可调的，就是DOS 里的展宽，很明显，你做的图不合理，因为展宽明显取得过小。
赞一下(2人) 回复此楼
一定要精通MaterialsStudios,VASP.
3楼2010-06-03 08:06:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)
★
zzy870720z(金币+2):谢谢专家指导，我试试 2010-06-03 09:48:31
hometownlove(金币+1):谢谢交流 2010-06-09 13:42:16
采用ISMEAR=-5时，需要将k-mesh增加的很大，出来的态密度才会光滑。

如果你想用适中的k-mesh，用ismear=0, SIGMA = 0.1~0.2也是可以接受的。这样弄出来的态密度，经过了Gussian展宽。

再要不然你用ISMEAR=-5计算出来的态密度，自己编个程序采用Guassian函数去展宽计算出来的态密度。
赞一下(2人) 回复此楼
6楼2010-06-03 08:40:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)
★
zzy870720z(金币+1):侯老师指导一下怎么做gaussian展宽？还有就是怎么gaussian处理后对DOS有没有影响？ 2010-06-03 11:07:45
hometownlove(金币+1):学习 2010-06-09 13:43:40
一般不从态密度上看能隙值，因为展宽的原因会使得价带顶和导带底的位置变得不准了。一般从计算的能带结构图上来看。
赞一下(2人) 回复此楼
10楼2010-06-03 10:32:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

zzy870720z

荣誉版主 (文坛精英)

优秀版主优秀版主优秀版主优秀版主
引用回帖:
Originally posted by valenhou001 at 2010-06-03 07:35:35:
1)增加k点网格的大小；
或者
2）采用ISMEAR=0, SIGMA = 0.1~0.2。

感谢侯老师指导，不过VASP手册上说计算态密度推荐采用四面体方法？
用Gaussian smearing方法好吗？
赞一下 回复此楼
博学、审问、慎思、明辨、笃学
4楼2010-06-03 08:25:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zzy870720z

荣誉版主 (文坛精英)

优秀版主优秀版主优秀版主优秀版主
引用回帖:
Originally posted by vasp001 at 2010-06-03 08:06:08:

换用ISMEAR =-5 ，作图时，还有一个参数是可调的，就是DOS 里的展宽，很明显，你做的图不合理，因为展宽明显取得过小。

首先谢谢你的解释
能说更详细些吗？作图时有什么要注意的？
麻烦了，很急
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

xiaobo1314
博学、审问、慎思、明辨、笃学
5楼2010-06-03 08:26:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dsagfds

银虫 (小有名气)
★ ★
zzy870720z(金币+1):谢谢提示，我做了静态计算 2010-06-03 09:48:47
hometownlove(金币+2):谢谢交流 2010-06-09 13:42:40
用VASP计算态密度需要执行两个步骤：
1、做静态计算
2.、将静态计算得出的CHGCAR文件，再加上（1）中的INCAR, KPOINTS, POSCAR和POTCAR一起作为输入文件。注意此时的INCAR中的ICHARG设为11。
增加DOS的平滑度可以增大NEDOS（约2000）来实现

你的问题是不是只做了第一步静态计算，然后利用静态计算的结果计算了态密度呢？
赞一下(2人) 回复此楼
让我真真实实的活吧
7楼2010-06-03 08:51:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zzy870720z

荣誉版主 (文坛精英)

优秀版主优秀版主优秀版主优秀版主
引用回帖:
Originally posted by valenhou001 at 2010-06-03 08:40:00:
采用ISMEAR=-5时，需要将k-mesh增加的很大，出来的态密度才会光滑。

如果你想用适中的k-mesh，用ismear=0, SIGMA = 0.1~0.2也是可以接受的。这样弄出来的态密度，经过了Gussian展宽。

再要不然你用ISMEAR=- ...

侯老师可以多提示一点吗？
没弄过gaussian展宽的东西？
有没有资料？谢谢
赞一下 回复此楼
博学、审问、慎思、明辨、笃学
8楼2010-06-03 10:02:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaojie7783

金虫 (小有名气)
引用回帖:
Originally posted by valenhou001 at 2010-06-03 08:40:00:
采用ISMEAR=-5时，需要将k-mesh增加的很大，出来的态密度才会光滑。

如果你想用适中的k-mesh，用ismear=0, SIGMA = 0.1~0.2也是可以接受的。这样弄出来的态密度，经过了Gussian展宽。

再要不然你用ISMEAR=- ...

请问如果采用了ISMEAR=0 SIGMA = 0.1~0.2的话，在态密度图上看到的能隙大小还准确吗？态密度图上看能隙是不是直接看费米能级附近态密度为0的两个点横坐标之间的差值呢？
赞一下 回复此楼
成功其实很简单，就是当你坚持不住的时候，再坚持一下！
9楼2010-06-03 10:19:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zzy870720z 的主题更新
351/4返回列表1234下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 362求调剂 +7 我要考大 2026-04-06 9/450 2026-04-07 09:04 by linbaixue
[考研] 一志愿211，0703化学305分求调剂 +10 严西西戏 2026-04-06 11/550 2026-04-07 08:38 by dick_runner
[考研] 327求调剂 +3 Xxjc1107. 2026-04-06 3/150 2026-04-06 22:34 by qlm5820
[考研] 调剂求助（生物与医药） +3 @6952 2026-04-06 3/150 2026-04-06 21:48 by Ecowxq666！
[考研] 266分，一志愿电气工程，本科材料，求材料专业调剂 +9 哇呼哼呼哼 2026-04-02 10/500 2026-04-06 19:47 by jean5056
[考研] 一志愿安徽某211 0703化学总分339求调剂 +7 晚风不晚 2026-04-04 7/350 2026-04-06 14:06 by houyaoxu
[考研] 一志愿 江南大学 085602 化工专硕 338分求调剂 +15 路痴小琪 2026-04-05 15/750 2026-04-06 09:27 by cql1109
[考研] 一志愿青科085500，初试295分，公共课213分 +3 遇到的人愿望都� 2026-04-05 3/150 2026-04-05 18:45 by 蓝云思雨
[考研] 0854求调剂 +4 assdll 2026-04-04 4/200 2026-04-05 09:44 by zhq0425
[考研] 材料调剂 +7 dxy调剂 2026-04-04 7/350 2026-04-05 09:15 by 陌秋26
[考研] 材料调剂 +12 一样YWY 2026-04-04 12/600 2026-04-05 08:24 by 544594351
[考研] 调剂 +11 JLLLLLLLLLL 2026-04-03 11/550 2026-04-04 22:21 by hemengdong
[考研] 292分，材料与化工，申请调剂 +22 程晴之 2026-04-01 26/1300 2026-04-04 22:03 by hemengdong
[考研] 277求调剂 +4 12A3 2026-04-02 5/250 2026-04-04 20:28 by 蓝云思雨
[考研] 调剂 +9 19945159693 2026-04-03 10/500 2026-04-04 20:16 by dongzh2009
[考研] 一志愿双非085502，267分，过四级求调剂 +3 再忙也要吃饭啊 2026-04-03 3/150 2026-04-04 05:03 by gswylq
[考研] 考研调剂 +3 Draa 2026-04-03 3/150 2026-04-03 17:37 by hgwz7468
[考研] 求调剂 +3 usbdndj 2026-04-03 3/150 2026-04-03 14:10 by dxiaoxin
[考研] 312 化工或制药调剂 +8 小小墨123 2026-04-02 9/450 2026-04-03 09:12 by zhouxiaoyu
[考研] 295求调剂 +7 愿旅途永远坦然 2026-04-02 7/350 2026-04-03 08:22 by fangshan711
信息提示
关闭
请填处理意见
关闭