VASP处理split_dos(5)

来O有s p 两个轨道，但是得到的ldos中一共有5列，其中第一列是能量，请问我怎么区分其余四列各是什么啊？？

还有就是各位有什么更好的处理方式么？

用p4vasp作图，挺方便的。你可以到网上下载p4vasp。

vasp manual DOSCAR里边什么都有。。。

nanual中的情况跟我的结果对应不上啊

下了一个，但是安装完了之后提示错误说：AttributeRrror：'module? object has no attribute ?VisIsosurface Drawer_draw_as_points_set? 解决不了这个问题

对于画DOS我一般用的是split—dos来处理DOSCAR 的。至于DOSCAR里面有多少列，那么这个和你的设置有关系，建议仔细看看手册。

第2-5列分别是S、py、pz、px

实例14

态密度和能带问题 各位朋友：我现在遇到的问题是对同一体系计算的态密度和计算的能带对应不上，不知道为什么原因，我所想到的办法都试过了，还是那样子，请大家指教，谢谢你们了

能带是和态密度一一对应, 你认真查看你处理DOSCAR和EIGENVAL文件的方式, 有些软件的处理会出错的,比如split-dos或者p4vasp不可能对应不上的

你给出INCAR等更多信息，我们才好帮你解决。否则只能说，你一定算错什么了

我想这个计算可能出错的地方是，第一DOSCAR计算错了，第二，对于该体系的布里渊区的高对称点没有弄对，所以bands错了。仔细检查吧

检查一下高对称点，如果确认设置没错，可能像二楼说的后处理有问题，画能带的时候要考虑对称性,p4vasp等软件可能没有考虑，不知道楼主计算什么材料，有没有过渡金属？

在想来估计是两个fermi能级计算有点出入，所以就不能对应。

现在你可以用静态计算的Fermi能，统一一下，估计就没有问题了。

实例16

关于态密度问题

我现在做总态密度图和分拨态密度。但是我得出结果是总态密度文件夹DOS0中第二列出现有总的态密度值，但是分波态密度第二列什么也没有也就是都是0，我郁闷得很，不知道怎么回事，请大家帮助。

总态密度文件好分波态密度文件怎么能这样比较了。直接打开DOSCAR读~读数据不就很清

楚了

实例17

lev00处理PDOS

我用lev00处理DOS要得到PDOS图，首先是在做到输入： Chemical elment (character*2） of special 1> 我按我POSCAR顺序里面的La输入La (谰)元素，它提示：RROR! specials La has not been recognised! 但输入其他元素却可以? 然后是我的PROCAR 明明在目录文件下，它提示：FATAL! the file PROCAR is bad or absent! 我也碰到过这个情况，估计在Levoo中没有包含稀土元素。所以稀土的元素符号是不认的，其实你完全可以用其它levoo认可的元素符号代替。应该没有任何影响的。

我输入了其他元素替代稀土元素La 但老提示我的PROCAR有问题，我不知道是因为我输入的其他元素不行，还是其他什么的。但我计算应该没有问题，因为我算了好几次啊！

难道Lev00就不可以处理稀土元素了吗？有没有谁用 lev00 得到过正确的PDOS图?还是跟Lev00的版本有关系？

诶有关系的，你想这个只是一个符号的问题，计算出来的电子密度等文件已经是确定的了，那还有什么问题呢。我认为输入任何符号都对结果没有影响的。

那难道是我的PROCAR 错了，因为它老提示：FATAL! the file PROCAR is bad or absent! 可是我计算了好多次结果都这样啊！我INCAR里面设置计算PDOS的参数是 LORBIT = 1 ; NPAR = 1; RWIGS = 1.535 0.82 1.302 1.217

想问一下，用origin画d电子态密度图，自旋向上我把在DOSCAR中最后10列的前五列直接加起来，自旋向下我又把后五列加起来，横坐标能量减去费米能可以吗？

如果没记错的话，spin up & spin down应该是相间排列的

你可以再网上找个split-dos（忘记是个script还是f code），可以把DOSCAR按原子序号分开

如果没记错的话，spin up & spin down应该是相间排列的

你可以再网上找个split-dos（忘记是个script还是f code），可以把DOSCAR按原子序号分开

实例18

态密度图跟文献里面的怎么都不一样，求救 弛豫过程 System=L ENCUT=520 GGA=91

VOSKOWN=1 ISTART=0 ICHARG=2

ISMEAR=-5 SIGMA=0.2 ISPIN=2

MAGMOM=4*0 16*0 4*0 4 -4 -4 4 EDIFF=0.1E-05 EDIFFG=-0.03 IBRION=2 POTIM=0.5 NSW=200 ISIF=3 NPAR=1 NSIM=8

LPLANE=.TRUE. LSCALU=.FALSE. NBANDS=200 PREC=High LDAU=.TRUE. LDAUTYPE=2

LDAUL=-1 -1 -1 2 LDAUU=0.0 0.0 0.0 5.3 LDAUJ=0.0 0.0 0.0 1.0 LDAUPRINT=2 静态过程

System=L ENCUT=520 GGA=91

VOSKOWN=1 ISTART=0 ICHARG=2 ISMEAR=-5 SIGMA=0.2 ISPIN=2

MAGMOM=4*0 16*0 4*0 8 -8 -8 8 EDIFF=0.1E-05 EDIFFG=-0.03 IBRION=-1 POTIM=0.5 NSW=0 ISIF=3 PREC=High LDAU=.TRUE. LDAUTYPE=2

LDAUL=-1 -1 -1 2 LDAUU=0.0 0.0 0.0 5.3

LDAUJ=0.0 0.0 0.0 1.0 LDAUPRINT=2 态密度过程

System=L ENCUT=520 GGA=91

VOSKOWN=1 ISTART=1 ICHARG=11 ISMEAR=-5 SIGMA=0.2 ISPIN=2

MAGMOM=4*0 16*0 4*0 8 -8 -8 8 EDIFF=0.1E-05 EDIFFG=-0.03 #IBRION=-1 #POTIM=0.5 #NSW=0 ISIF=3 PREC=High LDAU=.TRUE. LDAUTYPE=2

LDAUL=-1 -1 -1 2 LDAUU=0.0 0.0 0.0 5.3 LDAUJ=0.0 0.0 0.0 1.0 LDAUPRINT=2

主要问题不是你说的那样，而是 bandgap 我一直都算不出文献的3.7eV,我DOS过程算出来都是1.5eV,但是静态过程可以得到3.6eV

文献DOS显示体系是磁性的， 你的DOS显示的是非磁性的。可能自洽计算就出了问题。检查一下磁量子数的初始值设置有没有问题。

哪位高手路过拔刀相助下呀。。。主要问题是这半导体材料的DOS应该不对称，结果对称了，

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