从黄昆固体物理第九章第一节可知，击中厚度为x的材料的透射强度与x的关系如下
这里，r是反射系数，是归一化吸收系数，j(x )/J0是透射率。 光照射到材料表面时，一些光已经被反射，因此一般在计算透射光谱时，可以用阻挡函数t(x )表示，可以表示厚度x不同的材料对一定强度J0的光的阻挡能力。
t(x )=1-j(x )/j0(() ) ) () ) ) ) ) ) ) ) ) ) (但不限于此) ) ) ) )
本章以石墨为例。 笔者为了用笔记本电脑跑vasp，减少计算量，用MS对石墨单元进行了以下处理。 之后，简化的单元只包含两个原子，计算非常快。 当然，交换相关泛函是GGA-PBE，没有采用PBE HSE。 采用HSE的过程也基本相同，只是向INCAR追加HSE06的参数
首先在scf计算中，然后在grep NBANDS OUTCAR中检查NBANDS，在光学计算时将INCAR设定为该值2。
计算完成后，将vasp工具包安装单元的optical.sh文件复制到工作文件夹下，并在chmod x optical.sh中授予可执行权限。 然后，运行./optial.sh 该文件主要从vasprun.xml中提取介电函数的实部和虚部，执行vaspkit选择71，711，可以得到吸收系数、反射系数、折射率和光能的关系。
以下是为了吸收光谱，在这里取第一个峰值，将=1标准化的例子。
反射光谱如下，在与吸收光谱的峰值对应的能量下，反射系数为0.81。
有了以上数据，笔者通过直接写两段简单的python代码，可以得到阻止函数t(x )和x的关系。
最后，将Attenuation.dat拖动到orginlab进行绘制，即可获得石墨对13 eV光的透射光谱。 如下所示。
参考：
[1]黄昆.固体物理
-------------请参见
PS :因为很着急，所以可能有错误的地方。 欢迎同行们交流和学习。 另外，如有疑问、计算技术的讨论、协助，请联系笔者的邮箱地址。 1924311399@qq.com谢谢。