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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】VASP计算电子态密度
351/4返回列表1234下一页
查看: 8167  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 293求调剂 +6 世界首富 2026-03-11 6/300 2026-03-17 17:04 by ruiyingmiao
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 26考研求调剂 +6 丶宏Sir 2026-03-13 6/300 2026-03-17 16:13 by 醉在风里
[考研] 求调剂，总分315，考的生物医药，一志愿湖南师范大学。调剂到任何专业都可以 +4 小丁想进步 2026-03-11 5/250 2026-03-17 16:05 by 外星文明
[考研] 289求调剂 +6 步川酷紫123 2026-03-11 6/300 2026-03-17 10:23 by Sammy2
[考研] 285化工学硕求调剂（081700） +9 柴郡猫_ 2026-03-12 9/450 2026-03-17 10:18 by Sammy2
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
[考研] 326求调剂 +4 诺贝尔化学奖觊� 2026-03-15 7/350 2026-03-16 17:11 by 诺贝尔化学奖觊�
[考研] 085600材料与化工 求调剂 +13 enenenhui 2026-03-13 14/700 2026-03-16 15:19 by 了了了了。。
[考研] 0703 物理化学调剂 +3 我可以上岸的对� 2026-03-13 5/250 2026-03-16 10:50 by 我可以上岸的对�
[考研] 0703化学调剂 290分有科研经历，论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 080500，材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 0703化学调剂 +4 快乐的香蕉 2026-03-11 4/200 2026-03-13 22:41 by JourneyLucky
[考研] 310求调剂 +3 【上上签】 2026-03-11 3/150 2026-03-13 16:16 by JourneyLucky
[考研] 求调剂 +3 程雨杭 2026-03-12 3/150 2026-03-13 15:06 by JourneyLucky
[考研] 270求调剂 085600材料与化工专硕 +3 YXCT 2026-03-11 3/150 2026-03-13 10:13 by houyaoxu
信息提示
关闭
请填处理意见
关闭