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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 如何在vasp中TDOS费米能不归零
191/2返回列表12下一页
查看: 2463  |  回复: 18
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

shala123456

铁杆木虫 (知名作家)

[求助] 如何在vasp中TDOS费米能不归零 已有3人参与

用split 分DOSCAR后，分出的Tdos费米能级就在0ev处，比如费米能是7.1ev时，那么其他的能量是7.1的相对值。
现在我的问题是：如何在DOS图中Oev处不代表费米能，改为它的实际能量？
谢谢大家啦！

如何在vasp中TDOS费米能不归零
5555555555.png
回复此楼

» 猜你喜欢
1楼 2016-01-20 10:40:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangpeng5226

新虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
shala123456: 金币+6, ★★★很有帮助 2016-01-21 15:29:36
把数据导入origin中，能量整体加上费米能，然后以新的能量作图

发自小木虫IOS客户端
赞一下 回复此楼
2楼2016-01-21 14:05:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shala123456

铁杆木虫 (知名作家)
引用回帖:
2楼: Originally posted by wangpeng5226 at 2016-01-21 14:05:42
把数据导入origin中，能量整体加上费米能，然后以新的能量作图

谢谢你的回复，你说的是set column value.
我还想问的是：在VASP中有没有设置哪个参数，使split_dos后的dos的值已经加上了费米能（我的结果有点多，想偷懒不想一个个的加）.拜托了。
赞一下 回复此楼
3楼2016-01-21 15:29:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shala123456

铁杆木虫 (知名作家)
沉没了
回复此楼
4楼2016-03-22 13:59:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小学英语本

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
shala123456: 金币+12, ★★★很有帮助, 感谢你的回复 2016-05-21 17:19:01
shala123456: 金币+6, ★★★★★最佳答案, 只能在这里加点金币了，谢谢你 2016-05-22 15:30:26
费米能级处在0eV，是由于split_dos脚本导致的，VASP本身没有参数使费米能未0eV。
不知道你用的是哪个split_dos脚本，一般在脚本开始的注释处有以下注释：
# Modified to also align the Fermi level from
# the ground state to E=0.0
具体代码也做了使费米能出于0eV.
你可以把代码改了。
赞一下 回复此楼
5楼2016-05-18 17:04:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shala123456

铁杆木虫 (知名作家)
引用回帖:
5楼: Originally posted by 小学英语本 at 2016-05-18 17:04:43
费米能级处在0eV，是由于split_dos脚本导致的，VASP本身没有参数使费米能未0eV。
不知道你用的是哪个split_dos脚本，一般在脚本开始的注释处有以下注释：
# Modified to also align the Fermi level from
# the ...

我有个问题想请教下：如何更改你所说的代码，新手求指教。
赞一下 回复此楼
6楼2016-05-21 18:12:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小学英语本

金虫 (小有名气)
引用回帖:
6楼: Originally posted by shala123456 at 2016-05-21 18:12:11
我有个问题想请教下：如何更改你所说的代码，新手求指教。
...

用写字板打开split_dos脚本，查找 $efermi 你应该能找到类似这样的几行代码
sed -n ''$start','$end' p' $dosfile | awk '{printf "%12.8f %12.8f %12.8f %12.8f %12.8f %12.8f %12.8f \n", $1+(-1)*'$efermi', $2, -1*$3, $4, -1*$5, $6, -1*$7}' >> DOS$i
该行代码是由sed和awk两个处理文件的命令组合而成，其中
$1+(-1)*'$efermi'
$1:表示指定行第一个变量,即DOSCAR中的第一列能量；
$efermi:表示从OUTCAR中读取的费米能；
你把代码中所有$1+(-1)*'$efermi'改为$1即可。
赞一下 回复此楼
7楼2016-05-22 15:15:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小学英语本

金虫 (小有名气)
你可以学简单学习一下Shell编程，看懂了就可以自己随意更改了。我也是因为要处理 f 电子，才去看的代码。不同的计算所得的DOSCAR文件列数不同，所以需要根据实际情况更改代码。不知道你是出于什么目的，但是我建议还是将费米能级设在0eV,这样其它位置都是以费米能级作为参照。比如：你题目中提到的费米能级7.1ev是参照那个位置得出来的呢？
赞一下 回复此楼
8楼2016-05-22 15:30:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shala123456

铁杆木虫 (知名作家)
引用回帖:
8楼: Originally posted by 小学英语本 at 2016-05-22 15:30:26
你可以学简单学习一下Shell编程，看懂了就可以自己随意更改了。我也是因为要处理 f 电子，才去看的代码。不同的计算所得的DOSCAR文件列数不同，所以需要根据实际情况更改代码。不知道你是出于什么目的，但是我建议还 ...

好的是得好好看看代码，DOSCAR的列数可以通过赋INCARE 中的LORBIT的值来控制，0/10时，分为10列，2/12时18列。
当时老师提了个想法：就是不同比例成分tdos的费米能不同，如果将0处改为费米能，那么我所有的比例成分就可以在同一范围内分析其低能级的变化趋势。
题目中的7.1是OUTCARE中读出来的，7.1就是费米能，0ev处
赞一下 回复此楼
9楼2016-05-22 15:46:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shala123456

铁杆木虫 (知名作家)
引用回帖:
7楼: Originally posted by 小学英语本 at 2016-05-22 15:15:18
用写字板打开split_dos脚本，查找 $efermi 你应该能找到类似这样的几行代码
sed -n ''$start','$end' p' $dosfile | awk '{printf "%12.8f %12.8f %12.8f %12.8f %12.8f %12.8f %12.8f \n", $1+(-1)*'$ ...

你所说的脚本我没找到，能麻烦你发我一份吗
这是我常用的split dos 软件
https://yunpan.cn/cSYEDe8jMKwZE （提取码：014c）
我一般用这个劈裂dos

还有就是VASP中有一个自带的split_dos命令
这两个好像都没有脚本
赞一下 回复此楼
10楼2016-05-22 16:03:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 shala123456 的主题更新
191/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:58 by ZPa0EcMwuECS
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:58 by ZPa0EcMwuECS
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 01:36 by ZPa0EcMwuECS
[论文投稿] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:25 by ZPa0EcMwuECS
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 01:03 by ZPa0EcMwuECS
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +16 爱看书的可乐 2026-08-26 18/900 2026-08-30 00:59 by winsaint
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 00:52 by ZPa0EcMwuECS
[论文投稿] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 00:41 by ZPa0EcMwuECS
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +7 ASdOkHsho7FD 2026-08-28 10/500 2026-08-30 00:41 by ZPa0EcMwuECS
[基金申请] 基金不中，共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
关闭
请填处理意见
关闭