VASP处理split_dos 下载本文

内容发布更新时间 : 2024/11/1 6:49:30星期一 下面是文章的全部内容请认真阅读。

实例1

split-dos的工作原理

我算了一个BaTiO3单胞,5个原子,加上总的态密度应该一共6组数据。但是计算结果的DOSCAR中却有1+5*7,36组数据,不知道每组数据对应什么?用split-dos脚本处理之后,确实生成了DOS0……DOS5,这是怎么生成的?脚本看不懂,望指教其原理。我对比过DOS4和DOSCAR,里面的数据没相同的。

PS:是不是一定要运行vp后才能运行split-dos?

split_dos是把费米能级算进去了的,所以能量会有不同 但是其他的应该只是从DOSCAR里面提取出来的。

可以给我发第一列是能量 后面的分别是s-up s-down p-up p-down d-up d-down

从 脚本上看 vp 不是必须的

实例2

电子态密度图是怎么做的?

我看好多文献在做电子结构计算时,都有态密度图,请问这是怎么作的?装个p4vasp,里面可以导入计算所得到的vasp.xml文件,选择electronic -> DOS+BAND就可以画出态密度图,然后Graph -> export 就可以到处数据到Origin中画图。p4vasp还可以处理各种PDOS,并自动把自旋向下的电子态密度数据设为负值,非常方便。

我不用P4VASP 行吗?直接从DOSCAR里找出数据画图或者用split-dos给分割一下之后画图 这样做也是正确的吧?

实例3

split_dos运行

我将split_dos 和vp放在vasp第一次运行后的同一个目录下,然后运行,出现结果./split_dos[94]: /vp: not found [No such file or directory] 请问是什么原因,在split_dos应该修改路径吗

赋予执行属性

是chmod 755 split_dos吧 我赋予了呀。split_dos中不用改路径吗,我的目录下都没bin文件

那就是你缺少一个vp

修正版的splitdos有三个文件:vp、 sumdos和split_dos.ksh

实例4

散尽金币求split_dos 运行

1 我将split_dos 与vp 放在 需要计算的文件夹下 然后在当前目录运行 sh split_dos 出现以下错误

zhangkun@zhangkun-desktop:~/band750$ sh split_dos : not found2: : not found3: : not found8: : not found11: : not found12: : not found18:

: No such file or directory set: 20: Illegal option -A

zhangkun@zhangkun-desktop:~/band750$ 请问是怎么回事 ? 怎样才能计算出呢

2 还有为什么我自己编译的vasp p4vasp 不能读入vasprun.xml文件呢 有没有别的方法 求的局域态密度的呢

我用了split-dos 通过读入procar文件 这几个需要输入的参数是什么含义啊 为什么会显示数组越界呢

我有自己些的分开DOS文件的shell,不知道和你用的split_dos功能一样不,如果需要站内PM我

运行shell时不用 sh split_dos直接 . /split_dos就行了

我有自己些的分开DOS文件的shell,不知道和你用的split_dos功能一样不,如果需要站内PM我

实例5

VASP计算出来的PDOS怎么画图啊?求指教

用origin 画图软件画,你把VASP计算出来的PDOS弄成txt格式(当然得从DOSCAR 中 split dos 得到,这个需要小程序,你自己搜一下),有两列,一列能量,一列态密度,直接用origin打开就行。这个软件很有用的,像画能带结构什么的都是用的到的。做分子模拟和材料计算 material studio 和 origin 是必不可少的,论坛上有这两个软件的教程,自己学习一下。还有就是计算

然后: ./split_dos这样他就会帮你执行啦

实例6

VASP 分波态密度 求助: 我算一个有Zn,O,H的体系,在用VASP计算态密度的时候,用的是PAW势,在INCAR里设置LORBIT=11计算, 用split-dos程序分割每个原子的分波态密度,竟然在O原子里看到了有d电子的成分,觉得很奇怪,后来以为是INCAR的设置有问题,于是舍弃LORBIT=11, 换成RWIGS ,将Zn,O,H三个元素的Weige-Size 半径写到INCAR里,可是最后得到的DOSCAR里连每个原子的dos 都没有,只有一个总的dos.请教大家一下,这

到底是为什么,要如何设置才能得到每个原子正确的分波态密度。

应该还是要设LORBIT=11这个参数的,至于O出现d的成份,很正常,因为PAW里面O默认的RWIGS偏大有可能把金属原子的d成份包括进去了,不过你可以把O原子所谓的p和d画在一起看看,你会发现d相对于p来说很小,几乎可以忽略不计,而且其所在的能量范围也比p高(可能就是金属原子的d所在能量范围)。

我查了一下O的dos,发现d电子的量还是非常大,在有的能级上能和O的p电子同量级,甚至有的比p电子的量还大。

这是一个O原子的PDOS,第二个是第一个在费米能级附近的放大图

曾经尝试过手动输入RWIGS,调节O的RWIGS,当然我的体系中只含一个氧原子,调节总的DOS影响不大,且尽可能保证OUTCAR中的volume of typ 接近100%。可以消除O的D态。

至于楼主说的没有分波态密度,猜测莫非不是静态计算,因为优化过程不产生DOS。

实例7

怎样移动态密度图DOS的费米能级Ef才是正确的呢?

问题:为了同其他模型对比,计算出了晶体硅Si在gama点的态密度图,如下图所示,但老师说费米能级是错的。

计算按标准计算DOS步骤进行,经过了自洽计算->非自洽计算->DOSCAR分割,使用vdos脚本分割DOSCAR时已经自动移动了费米能级。

可是经过这样做得到的费米能级似乎是错的,其位置应该是在红色点线附近才对。不知应该怎样移动费米能级才是正确的呢?从图中看应该移动大约2 eV左右,但这个值从哪里得到。

盼高手解答。