24小时热门版块排行榜    

查看: 2024  |  回复: 18

小学英语本

金虫 (小有名气)

引用回帖:
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朵)

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 的主题更新
信息提示
请填处理意见