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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 如何在vasp中TDOS费米能不归零
192/2返回列表上一页12
查看: 2447  |  回复: 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
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-29 18:59 by 4FFAWE8HcgUD
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-29 18:15 by 4FFAWE8HcgUD
[基金申请] 有没有仍没收到信息的 +5 德尚中行 2026-08-27 6/300 2026-08-29 18:10 by manplx
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-29 18:04 by 4FFAWE8HcgUD
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 17:31 by 4FFAWE8HcgUD
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +5 gdd2018 2026-08-28 10/500 2026-08-29 14:10 by Jacob678
[论文投稿] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 13:51 by jCd0dEvKHShX
[找工作] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 11:51 by jCd0dEvKHShX
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-29 11:18 by jCd0dEvKHShX
[基金申请] 中青基了要发朋友圈吗？ +3 349506619 2026-08-28 3/150 2026-08-29 10:49 by 木风9012
[基金申请] 我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？ +6 gouxfjh 2026-08-28 9/450 2026-08-29 09:59 by jklily
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 国合现在查不到了吗？ +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
信息提示
关闭
请填处理意见
关闭