求助vasp加入SOC计算态密度问题

wei97zheng

主要计算的是铈离子。
最开始是直接调用vasp_std计算态密度，有f层分量。
再加入SOC后，调用vasp_ncl进行态密度计算，最后使用vaspkit得到态密度分布，但是结果的DOSCAR里只有spd的轨道分量（只有37行，分别是energy，s（1*4），p（3*4），d（5*4）），没有f轨道。
并且vaspkit给出的态密度分布有问题，Energy只有两个点，不知道该怎么处理，有没有人可以指导一下，如下图所示：

主要问题有：（1）SOC计算后就没有f轨道了吗，我想看的就是f轨道怎么办？
（2）态密度分布如果vaspkit不能处理的话，我自己在DOSCAR里找相应的列进行绘图可以吗？

wei97zheng

SOC计算的代码如下：

Spin-Orbit Coupling Calculation
LSORBIT    = .TRUE.    (Activate SOC)
GGA_COMPAT = .FALSE.   (Apply spherical cutoff on gradient field)
VOSKOWN    =  1        (Enhances the magnetic moments and the magnetic energies)
LMAXMIX    =  6        (For d elements increase LMAXMIX to 4, f: LMAXMIX = 6)
ISYM       =  0       (Switch symmetry off)
# SAXIS    =  0 0 1    (Direction of the magnetic field)
# MAGMOM   =  0 0 3    (Set this parameters manually, Local magnetic moment parallel to SAXIS, 3*NIONS*1.0 for non-collinear magnetic systems)
NBANDS     =  720        (Set this parameters manually, 2 * number of bands of collinear-run)
LNONCOLLINEAR = .TRUE.

直接使用的vaspkit处理态密度功能，得到数个文件

sobereva

wei97zheng 发表于 2022-5-30 09:28
SOC计算的代码如下：

Spin-Orbit Coupling Calculation

有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，这点在置顶的新社员必读贴里明确说了。

wei97zheng

sobereva 发表于 2022-5-30 10:00
有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，这点在置顶的新社员 ...

好的，下次一定注意，我是怕图顺序乱掉，下次一定放一起。

柠檬鬼lemon

您好，请问这个问题您现在解决了吗，我最近在做的体系中也有f轨道，加上soc之后算dos就没有f轨道了，请问这是vasp程序的问题+soc后确实计算不出f轨道，还是可以通过修改什么标签获得f轨道呀

lxylalalala

请问解决了吗，如何解决？