DOSCAR
DOSCAR文件包含DOS和integarated DOS,单位分别是状态数/eV和状态数。通过除以单位晶胞的体积来获得The intensive DOS。
DOSCAR文件的前几行的组成为:
Number of Ions (including empty spheres), 原胞中的离子数, 0(no partial DOS) or 1(incl. partial DOS), NCDIJ (currently not used)//
原胞体积 基矢量长度(a,b,c [m]) POTIM[s]//
最初的温度TEBEG//
'CAR'//
INCAR中的SYSYTEM注释 //
E(max), E(min), NEDOS, E(fermi), 1.0000//
接下来lines的format:
-
不考虑自旋ISPIN=1
energy DOS integrated DOS
态密度(DOS)
由两个自旋的差分决定:
是两个再DOSCAR文件中的两个grid point的能量差;
是DOS的积分,就是小于某个能量的所有粒子数:
-
考虑自旋ISPIN=2
energy DOS(up) DOS(dwn) integrated DOS(up) integrated DOS(dwn)
如果在INCAR文件中设置了RWIGS或LORBIT(对Wigner Seitz半径很重要),则计算每个离子的lm和点投影的DOS并将其写入DOSCAR文件。此数据同样包含具有各种列的NEDOS行,具体取决于LORBIT,ISPIN和LNONCOLLINEAR(决定自旋线性/非线性)的选择。
-
Colinear计算中:
-
ISPIN=1,每个离子的实部部分的format如下:
energy s-DOS p-DOS d-DOS
而虚部Im部分的format为:
energy s p_y p_z p_x d_{xy} d_{yz} d_{z2-r2} d_{xz} d_{x2-y2},...
-
ISPIN = 2实数部分每个自旋打印额外的列:
energy s-DOS(up) s-DOS(down) p-DOS(up) p-DOS(dwn) d-DOS(up) d-DOS(dwn)
<u>相应地lm分解的密度由LORBIT选择决定</u>。
-
-
Non-collinear计算中:
-
仅能获取DOS点投影各个自旋分量部分的信息:
energy s-DOS(total) s-DOS(mx) s-DOS(my) s-DOS(mz) p-DOS(total) p-DOS(mx) ...
在这种情况下,不同方向上的磁化强度密度(
)也相应的获取了。、
-
注意
在所有情况下,实部和DOS点投影的单位是状态/原子/能量(states/atom/energy)。
在以下情况下,并行版本不估计DOS的投影:
-
vasp.4.5:NPAR
1时,没有DOS的点投影;
-
vasp.4.6:NPAR
1时,LORBIT = 0-5,则没有DOS点投影
在vasp.4.6中,<u>当LORBIT = 10-12时</u>,即使满足NPAR不等于1(与之前的版本相反),依旧可以估计DOS的点投影
-
vasp.5:不需要规定NPAR
-
Mind:<u>对于relaxation,DOSCAR通常是无用的。如果要获得准确的DOS,首先将CONTCAR复制到POSCAR并继续进行一次静态(ISTART = 1; NSW = 0)计算。</u>
网友评论