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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】VASP怎样精确计算出费米能级处的态密度值？
81/1返回列表
查看: 4095  |  回复: 22
查看全部回帖@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

xinqing98

木虫 (小有名气)

少将

[交流] 【求助】VASP怎样精确计算出费米能级处的态密度值？

大家好！
   我想精确计算出体系费米能级处的态密度值，
用VASP该怎么算呢？
   如果只是做一次态密度计算，从DOS图中得到的
费米能级处态密度值太粗略了，大家有好的办法吗？
谢谢了！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
青春就是勇于实践
1楼2010-09-08 14:56:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)
xinqing98(金币+1):谢谢你！ 2010-09-08 15:51:41
引用回帖:
Originally posted by xinqing98 at 2010-09-08 14:56:31:
大家好！
   我想精确计算出体系费米能级处的态密度值，
用VASP该怎么算呢？
   如果只是做一次态密度计算，从DOS图中得到的
费米能级处态密度值太粗略了，大家有好的办法吗？
谢谢了！

vasp不是给出费米能级么
赞一下 回复此楼
心之所向无惧无悔
2楼2010-09-08 15:43:59
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)
★
aylayl08(金币+1):感谢讨论 2010-09-08 19:03:59
引用回帖:
Originally posted by xinqing98 at 2010-09-08 15:53:09:


VASP虽然可以给出费米能级的值，但是做态密度计算给不出
费米能级处的精确态密度值啊！ 我需要费米能级处的精确
态密度值，该怎么办呢？

我不知道vasp还能怎样计算精确态密度，直接从DOSCAR里找到Fermi-level处的态密度不可以吗
赞一下(1人) 回复此楼
心之所向无惧无悔
4楼2010-09-08 15:57:38
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)
★
aylayl08(金币+1):谢谢回复 2010-09-08 19:03:47
引用回帖:
Originally posted by xinqing98 at 2010-09-08 16:18:37:


DOSCAR中能量的取点只有301个，在里边是找不到Fermi-level的准确位置的，那这样就读不出Fermi-level处的态密度。

INCAR里可以设置NEDOS，默认是300你可以加的多一点比如1500
赞一下(1人) 回复此楼
心之所向无惧无悔
6楼2010-09-08 18:11:33
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)
★
zzy870720z(金币+1):谢谢讨论 2010-09-09 08:25:15
引用回帖:
Originally posted by meinvren at 2010-09-08 20:20:21:


改变NEDOS有时侯是不科学的。比如BZ是５＊５＊５而NEDOS＝１５００，那么总共K点是多少？远远小于NEDOS！DOS计算时NEDOS不能大于K点，否则没有意义。

还没有这样想过，那你说默认的NEDOS就是合理的吗？而且算DOS是固定优化后的电荷密度算的，跟原来的K点设置没什么关系，就是一个固定的薛定谔方程求本征值的过程
赞一下(1人) 回复此楼
心之所向无惧无悔
10楼2010-09-08 22:15:17
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)
★ ★
zzy870720z(金币+2):谢谢讨论 2010-09-09 15:40:57
引用回帖:
Originally posted by meinvren at 2010-09-09 02:57:18:


根据DOS(E0)=sum_{i,k}{delta(E-E0)}。不管是你所谓的固定薛定谔方程求本征值还是不固定，DOS只与你进行DOS计算时的kmesh有关。如果能量点NEDOS远大于kmesh的k点数目，结果已经没有意义。

不敢苟同
据我所知，NEDOS是将能量下限到费米能级之间均匀的插入NEDOS个能量点。如果K点总数是100，那么它会在一个能带中标出100个不同的能量值（先不考虑简并）。而一个能带显然不能贯穿能量下限到费米能级这样大的能量范围，这之间通常包含很多能带，每个能带的能量范围其实是很窄的，也就是说这些K点经常可以将能量分的很细。默认的NEDOS我觉得远没有达到这个标准，再提高几倍没有问题
赞一下(2人) 回复此楼
心之所向无惧无悔
13楼2010-09-09 08:29:51
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)
★ ★
zzy870720z(金币+2):谢谢讨论交流 2010-09-09 15:41:38
引用回帖:
Originally posted by meinvren at 2010-09-09 14:48:29:


你不妨做个测试。比如３＊３＊３mesh，NEDOS＝３０００。DOS会怎样？许多许多的刺！如果提高一直提高kmesh直到足够密，最终不会有这样的情况。

我觉得你对NEDOS与K点关系的理解有问题。K点的作用在于将一个分子轨道扩展成一个能带，因此DOS图中K点只是起到分割一个能带的作用。通常能带图可以看到每一个能带特别是接近费米面时是很窄的，再进行分割其能量间隙会很小，所以适当的增加NEDOS减小统计能量间隔是合理的

有尖峰是很正常的，因为我们取的K点本来就是离散的，能量本征值自然不能覆盖所有位置，DOS出现尖峰正是对离散的K点一种反映，如果人为地减少NEDOS从而使DOS变得平滑你觉得合理吗？
赞一下(1人) 回复此楼
心之所向无惧无悔
16楼2010-09-09 15:17:12
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)
引用回帖:
Originally posted by meinvren at 2010-09-10 18:35:32:


应该是你理解有问题。我只是强调NEDOS的取值的合理范围，没有说人为平滑。你认为NEDOS取值过大得到许多带刺的DOS还有必要用于论文发表吗？这并不是简单地随机取值的问题，我已经在前面给出公式了。只是属于数 ...

呵呵，那我们就各持己见吧，我尊重你的看法
赞一下 回复此楼
心之所向无惧无悔
18楼2010-09-10 20:04:18
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xinqing98 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭