Vkit计算光学性质详细步骤与参数设置指南.pdfVIP

Vaspkit计算光学性质具体步骤：

1.准备优化好的POSCAR如：cat./OPTIC/POSCAR(注：设当前计算目录为OPTIC)

systemSiC

4.35

0.50.50.0

0.00.50.5

0.50.00.5

SiC

11

cart

0.000.000.00

0.250.250.25

2在OPTIC目录运行：vaspkit-task101-inpST获取静态计算的INCAR

在OPTIC目录运行：vaspkit-task102-kpr0.04-kpsG获取静态计算KPOINTS和POTCAR

ls如下：

3输入vasp运行命令，如：mpirun-np20vasp_std.base(我的运行命令，你可能不同)

4运行完成后，继续输入：vaspkit-task101-inpOP或取光学计算INCAR，包含两个：1）

INCAR,2)INCAR.step2.lef

INCAR:INCAR.step2.lef

注：默认INCAR基于独立粒子近似计算

INCAR.step2.lef考虑局域场效应(ALGO=CHI)，其中LRPA=TRUE（默认）计算的是RPA层

级的局域场，只限于Hartree贡献，如果想计算DFT交换关联的贡献，请设置LRPA=.FALSE.

(去掉感叹号！)

5输入grepNBANDSOUTCAR获取静态计算的默认NBANDS数量（e.g.NBANDS=20）

6输入sed-i3cNBANDS=60INCAR（自动替换INCAR第二行NBANDS数量为60，取足够多空

带）

seds/NBANDS=/NBANDS=60/gINCAR.step2.lef

7输入vasp运行命令，如：mpirun-np20vasp_std.base(我的运行命令，你可能不同)

8运行optic.sh脚本（在vaspkitexamples/optic/optic.sh），建议在/home目录新建一个bin目录，并

把bin目录加入环境变量，然后把常用的小脚本放在/home/bin下，比如我把optic.sh放在

/home/bin/下，执行chomd+xoptic.sh赋予执行全，并在~/.bashrc文件中添加一行：

PATH=/home/bin/:$PATH，然后保存退出后输入：source~/.bashrc后，进入当前计算目录

（xx/OPTIC/）直接输入optic.sh,此时会产生IMAG.in（虚部）和REAL.in（实部）两个文件供

vaspkit后续使用。

其中IMAG.in每列数据物理意义参考如下图

REAL.in每列数据物理意义：

9输入vaspkit-task71获取ABSORB.dat（吸收），REFRACTIVE.dat（折射），

ENERGYLOSSSPECTRUM.dat(能量损失)，REFLECTIVITY.dat（反射），EXTINCTION.dat(灭绝)

以上数据都具如下图片的一个格式：

第一列能量区间（横轴数据），后六列，相应方向的物理量强度（或大小），各向同性或异性可通过xx

yyzz简单判断，如相同，各向同性，此时只需要用到第一列和第二（三，四）画图，如不同，

你随意，用第一列与第二、三、四画，或用第一列与第二、三、四列平均值画都有文献采用。

10用第二个INCAR.setp2.lef计算的流程与上一致，只是得到的最后结果和计算的时间上会有差别，此