1.画PARCHG
只有当预收敛的 WAVECAR 文件存在时,才能计算部分电荷密度。
具体内容在vasp wiki中band decomposed charge densities可以查到。
必须注意的是,VASP计算的密度(包括带分解电荷密度)总是使用VASP确定的空间群和点群对称性进行对称化。在某些情况下(从选定的k点计算电荷),这可能导致带分解电荷密度的不良结果。在这种情况下,需要关闭对称性以进行基态和连续带分解电荷密度计算。
其中包含7个参数:
1.LPARD:评估部分(带和/或k点)分解的电荷密度,设置为.TURE.时计算分解电荷密度。
如果设置LPARD为.TURE.而不设置其他参数,则计算总电荷密度。
2.IBAND:计算阵列 IBAND 中指定的所有波段的部分电荷密度。例如设置IBAND = 192 193,则计算第192 193条能带的电荷密度。
3.EINT:指定用于评估部分电荷密度的波段的能量范围。应给出两个实值,可以用来区分自旋向上和自旋向下,设置一个窄的能量范围,只算某一能带某一自旋的电荷密度。
4.NBMOD:
>0:数组 IBAND 中的值数。如果指定了 IBAND,则 NBMOD 将自动设置为正确的值(不需要再在 INCAR 文件中手动设置 NBMOD)
0:取所有波段来计算电荷密度,甚至不占用的波段也要考虑在内。
-1:像往常一样计算总电荷密度。如果未给出任何其他内容,则为默认值。
-2:计算具有EINT指定范围内特征值的电子的部分电荷密度。
-3:与以前相同,但给出了能量范围与费米能量的关系。
5.KPUSE:指定在评估部分 DOS 时使用哪些 k 点。例如KPUSE= 1 2 3 4
表示对前四个 k 点的电荷密度进行评估和求和。
6.LSEPB:指定是否单独计算每个波段的电荷密度并将其写入文件 PARCHG.nb.* 设为.TURE.默认为FALSE,表示写入一个文件。
7.LSEPK:指定是否将每个 k 点的电荷密度写入文件 PARCHG.*.nk,设为.TURE.默认为FALSE,表示写入一个文件。
计算某单独条件的电荷密度,后面两个参数可以不设置,直接取默认值,将结果写入一个大文件。
计算时复制WAVECAR,CONTCAR复制为POSCAR,将INCAR设置为静态计算的INCAR,并加入相关参数,使用一样的POTCAR和KPOINTS,即可计算。
LPARD = .TRUE.
IBAND = #需要画的轨道数
EINT = #预估能量值,为两个实数
NBMOD = #如果设置IBAND就不用设置了,如果设置EINT了之后设为-2
KPUSE = #计算的K点个数,1 2 3 4 为前四个K点
LSEPB = #放到一个文件里还是分开,分开是TURE,默认为FALSE
LSEPK = #K点的态密度写到一起还是分开,分开是TURE