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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 如何在vasp中TDOS费米能不归零
192/2返回列表上一页12
查看: 2558  |  回复: 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一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-01 07:01 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 07:01 by YXM1FDCSE8Bn
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 06:53 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 X5bHxOULxFqb 2026-10-30 17/850 2026-11-01 06:48 by YXM1FDCSE8Bn
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 X5bHxOULxFqb 2026-10-30 21/1050 2026-11-01 06:46 by YXM1FDCSE8Bn
[基金申请] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:48 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:43 by YXM1FDCSE8Bn
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 05:41 by YXM1FDCSE8Bn
[论文投稿] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 03:59 by YXM1FDCSE8Bn
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 03:56 by YXM1FDCSE8Bn
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 10/500 2026-11-01 03:51 by YXM1FDCSE8Bn
[论文投稿] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 X5bHxOULxFqb 2026-10-30 25/1250 2026-11-01 03:30 by YXM1FDCSE8Bn
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 02:25 by YXM1FDCSE8Bn
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 12/600 2026-11-01 02:19 by YXM1FDCSE8Bn
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 4/200 2026-11-01 00:41 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 X5bHxOULxFqb 2026-10-30 12/600 2026-11-01 00:15 by YXM1FDCSE8Bn
[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 X5bHxOULxFqb 2026-10-30 13/650 2026-10-31 17:12 by hKe7U1kad5Sg
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +8 zuesli 2026-10-27 14/700 2026-10-30 22:44 by 好运连连公主
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
[考博] 申博求助 20+3 Katrinae 2026-10-26 4/200 2026-10-27 09:56 by 北京莱茵编辑
信息提示
关闭
请填处理意见
关闭