VASP处理split_dos

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

盼高手解答。

附：

1、晶体硅Si的计算模型图

2、VASP计算的INCAR文件、KPOINTS文件。

INCAR文件（非自洽计算） general:

System = Crystal Si cell

ISTART = 0 ; ICHARG=11 #自洽计算 ICHARG=2; 非自洽计算 ICHARG=11 PREC=Accurate ENCUT=320

ISMEAR = 0; SIGMA = 0.1 IALGO = 48 GGA=91 ISPIN=2

VOSKOWN = 1 AMIX=0.2 BMIX=0.0001 AMIX_MAG=0.8 BMIX_MAG=0.0001 MAGMOM= 8*0

LREAL = .FALSE. #option for system contain less than 20 atoms

NPAR = 1 LWAVE = .F.

LORBIT = 1 #非自洽计算选项，计算spd分波态密度 RWIGS = 1.312

非自洽计算gama点态密度DOS时的KPOINTS文件： K-Points 0 Gama 1 1 1 0 0 0

费米能级的移动是什么意思呢？

费米能级不是确定的吗？在自洽计算得到的OUTCAR文件中有的。在非自洽计算得到的OUTCAR文件中也有的。这两处得到的费米能级有些差别，有的似乎较大。

做split dos的时候，会从与DOSCAR放在一起的的OUTCAR文件中读费米能级E-fermi。然后就自动以这个费米能级为零点。

我个人觉得，费米能级算不准，不同软件，以及vasp中静态自洽和非自洽计算结果也有差别。所以有文献干脆不提费米能级，而只以价带顶为参考零点。

你的体系并不大，静态自洽的K点只取一个似乎不好，可以取7*7*7.

费米能级是否能算的准，比如费米能级是否能出现在浅受主能级与价带顶之间，或者本征绝然体带隙正中间？我没有算出来过。

1、费米能级是能够确定的。

费米能级由vdos脚本（或split_dos脚本）自动从DOSCAR中抽取出来的，然后得到上面的DOS图。

2、k点是故意选成单个gama点的。

由于我想看gama点的态密度，所以只设定了一个k点，因而k点个数是不必增加的。 3、费米能级的确会算不准。

之前查看过自洽计算和非自洽计算的OUTCAR，它们的费米能级有差别，因而在作图时应该考虑自洽计算和非自洽计算费米能级的差值。

我的困惑在于，现在图中标示的费米能级是错的，应该是在大约红线的位置，可是不知怎之前的自洽计算，k点取了多少？

自洽计算时的KPOINTS文件： K-Points 0 m

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说综合文库VASP处理split_dos在线全文阅读。