内容发布更新时间 : 2025/7/2 21:58:47星期一 下面是文章的全部内容请认真阅读。

实例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左右，但这个值从哪里得到。

盼高手解答。