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

24小时热门版块排行榜    

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

小学英语本

金虫 (小有名气)
引用回帖:
10楼: Originally posted by shala123456 at 2016-05-22 16:03:52
你所说的脚本我没找到，能麻烦你发我一份吗
这是我常用的split dos 软件
https://yunpan.cn/cSYEDe8jMKwZE （提取码：014c）
我一般用这个劈裂dos

还有就是VASP中有一个自带的split_dos命令
这两个好像都没 ...

软件的话，不是有一个P4v的软件吗？
split_dos是用于劈裂DOSCAR文件的，它将按照不同原子进行劈裂会得到DOS0、DOS1、DOS2、……DOSn文件。其中n为体系原子序数。DOS0为总态密度图。DOSn为所对应序号原子的分波态密度。
用法./split_dos (当前目录必须包含OUTCAR和DOSCAR文件)
sum_dos是用于将不同DOSn文件合在一起。
用法 ./sum_dos 0 1 n 表示合并DOS1直到DOSn文件。（第一个参数0表示连续合并）
用法./sum_dos 1 1 4 5 9表示合并指定的DOS1、DOS4、DOS5、DOS9文件。（第一个参数1表示指定合并）
sum_dos_np与sun_dos用法一样，用于合并非自旋极化（ISPIN=1）计算的DOS文件。
赞一下 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : split_dos
    • 2016-05-22 17:32:41, 2.87 K
  • 附件 2 : sum_dos
    • 2016-05-22 17:33:06, 2.04 K
  • 附件 3 : sum_dos_np
    • 2016-05-22 17:33:38, 1.72 K
11楼2016-05-22 17:53:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小学英语本

金虫 (小有名气)
忘记说，这个是只处理LORBIT设置大于等于10；费米能级与0eV对齐；而且也不处理f电子。
你看懂代码可以自己改。DOSCAR文件的组成你可以看VASP教程，我目前只知道
E（能量） s（上） s（下） p（上） p（下） d（上） d（下）7列的DOSCAR文件，（上、下表示自旋方向）
E（能量） s（上） s（下） p（上） p（下） d（上） d（下）f（上） f（下） 9列的DOSCAR文件，
对于十几列的DOSCAR文件，我只能确认第一列是能量。
赞一下 回复此楼
12楼2016-05-22 18:06:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shala123456

铁杆木虫 (知名作家)
引用回帖:
11楼: Originally posted by 小学英语本 at 2016-05-22 17:53:29
软件的话，不是有一个P4v的软件吗？
split_dos是用于劈裂DOSCAR文件的，它将按照不同原子进行劈裂会得到DOS0、DOS1、DOS2、……DOSn文件。其中n为体系原子序数。DOS0为总态密度图。DOSn为所对应序号原子的分波态密 ...

一直想找类似的命令，但没找到，想想当初用origin一个个的加就觉得心塞
我刚刚试了下，挺方便的，但这个命令不能适用于自旋极化（ISPIN=2)吗？还是说得出的结果不准，求指教
赞一下 回复此楼
13楼2016-05-22 20:07:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shala123456

铁杆木虫 (知名作家)
引用回帖:
12楼: Originally posted by 小学英语本 at 2016-05-22 18:06:44
忘记说，这个是只处理LORBIT设置大于等于10；费米能级与0eV对齐；而且也不处理f电子。
你看懂代码可以自己改。DOSCAR文件的组成你可以看VASP教程，我目前只知道
E（能量） s（上） s（下） p（上） p（下） d（上 ...

关键是当LORBIT>10时rwigs是不读的，会不会使结果不准呢？
赞一下 回复此楼
14楼2016-05-22 20:15:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小学英语本

金虫 (小有名气)
引用回帖:
13楼: Originally posted by shala123456 at 2016-05-22 20:07:01
一直想找类似的命令，但没找到，想想当初用origin一个个的加就觉得心塞
我刚刚试了下，挺方便的，但这个命令不能适用于自旋极化（ISPIN=2)吗？还是说得出的结果不准，求指教...

split_dos可以用于ISPIN=1，ISPIN=2两种情况。
sum_dos用于ISPIN=2的情况，sum_dos_np用于ISPIN=1的情况。
赞一下 回复此楼

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

shala123456
15楼2016-05-23 08:43:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小学英语本

金虫 (小有名气)
引用回帖:
14楼: Originally posted by shala123456 at 2016-05-22 20:15:28
关键是当LORBIT>10时rwigs是不读的，会不会使结果不准呢？...

这个和准不准没关系，关键看你计算什么性质了。有时候需要设置LORBIT>10。你可以根据自己情况把代码改了，脚本中读取LORBIT值的代码原本意义不大，只是为了说明在VASP4.*版本情况下LORBIT>10不可靠。
赞一下 回复此楼
16楼2016-05-23 08:52:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shala123456

铁杆木虫 (知名作家)
送红花一朵
引用回帖:
15楼: Originally posted by 小学英语本 at 2016-05-23 08:43:10
split_dos可以用于ISPIN=1，ISPIN=2两种情况。
sum_dos用于ISPIN=2的情况，sum_dos_np用于ISPIN=1的情况。...

感谢大神的解惑，这下放心啦
赞一下 回复此楼
17楼2016-05-23 13:54:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shala123456

铁杆木虫 (知名作家)
引用回帖:
16楼: Originally posted by 小学英语本 at 2016-05-23 08:52:51
这个和准不准没关系，关键看你计算什么性质了。有时候需要设置LORBIT>10。你可以根据自己情况把代码改了，脚本中读取LORBIT值的代码原本意义不大，只是为了说明在VASP4.*版本情况下LORBIT>10不可靠。...

好的，明白了，谢谢你
赞一下 回复此楼
18楼2016-05-23 13:56:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

前田庆次

木虫 (著名写手)

【答案】应助回帖

outcar得到费米能级 Origin第一列能量加上费米能级即可
赞一下 回复此楼
长枪一横花飘零，松风追月伴我行。无双人间世如梦，倾奇万世永留名。
19楼2017-04-08 10:30:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 shala123456 的主题更新
192/2返回列表上一页12
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +3 十三加油 2026-03-21 3/150 2026-03-21 13:18 by zhukairuo
[考研] 299求调剂 +3 某某某某位 2026-03-21 3/150 2026-03-21 13:11 by zhukairuo
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂！ +3 不想理你 2026-03-16 3/150 2026-03-21 03:59 by JourneyLucky
[考研] 初始318分求调剂（有工作经验） +3 1911236844 2026-03-17 3/150 2026-03-21 02:33 by JourneyLucky
[考研] 304求调剂 +6 曼殊2266 2026-03-18 6/300 2026-03-21 00:32 by JourneyLucky
[考研] 295求调剂 +4 一志愿京区211 2026-03-18 6/300 2026-03-20 23:41 by JourneyLucky
[考研] 294求调剂材料与化工专硕 +15 陌の森林 2026-03-18 15/750 2026-03-20 23:28 by JourneyLucky
[考研] 321求调剂 +9 何润采123 2026-03-18 11/550 2026-03-20 23:19 by JourneyLucky
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +4 葵梓卫队 2026-03-18 6/300 2026-03-20 23:02 by JourneyLucky
[考研] 323求调剂 +3 洼小桶 2026-03-18 3/150 2026-03-20 22:54 by JourneyLucky
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[考研] 北科281学硕材料求调剂 +5 tcxiaoxx 2026-03-20 5/250 2026-03-20 21:35 by laoshidan
[考研] 一志愿武理材料工程348求调剂 +3 ￣^￣゜汗 2026-03-19 4/200 2026-03-20 21:01 by zhukairuo
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +22 rare12345 2026-03-18 22/1100 2026-03-20 20:39 by zhukairuo
[考研] 260求调剂 +3 朱芷琳 2026-03-20 3/150 2026-03-20 20:35 by 学员8dgXkO
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
[考研] 材料与化工求调剂 +7 为学666 2026-03-16 7/350 2026-03-19 14:48 by 尽舜尧1
[考博] 26申博 +4 八6八68 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[论文投稿] 有没有大佬发小论文能带我个二作 +3 增锐漏人 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
信息提示
关闭
请填处理意见
关闭