vasp计算光学

我做过一个半导体光学性质计算的项目，客户想预测他的材料在可见光到紫外波段的光吸收和折射率。我这个项目交付的核心就是复介电函数谱和吸收/折射率谱两张图。

光学性质的计算基础是介电函数 ε(ω) = ε1(ω) + iε2(ω)。ε2 描述材料对光的吸收能力，ε1 描述极化响应。两者通过 Kramers-Kronig 关系联系起来。左侧图给出的是 0-8 eV 光子能量范围内的 ε1 和 ε2。蓝色实线是 ε1，在 2 eV 以下从 1.3 上升到峰值约 4.0，然后逐渐下降，在 6 eV 附近穿过 1.0；红色虚线是 ε2，阈值约在 2 eV，主吸收峰在 3.5 eV，第二峰在 5 eV 附近。ε2 的阈值位置对应带边跃迁能量，对于间接带隙材料，阈值可能不完全等于带隙，因为需要声子参与。

右侧图是吸收谱 α 和折射率 n。吸收谱在 3-4 eV 达到主吸收峰，和 ε2 的峰位对应；折射率 n 在 2-3 eV 区间达到约 2.0 的平台，随后下降。可见光波段（1.6-3.1 eV）内 n 在 1.8-2.0 之间，说明材料对可见光有较强的折射能力，适合做光学涂层或光电器件。

计算上，我用 VASP 的线性响应理论（DFPT）或 RPA 近似算介电函数。具体设置：先做一个高精度的自洽计算，ENCUT=500 eV，K 点 12×12×12；然后非自洽算介电函数，NEDOS=2000，LOPTICS=.TRUE.，CSHIFT=0.1 eV。CSHIFT 控制 Lorentz 展宽，太小会导致峰很尖锐但可能有数值噪声，太大峰会变宽变矮。我对这个体系试了几个值，0.1 eV 在平滑度和分辨率之间比较平衡。

踩坑经验：第一个是带隙低估。PBE 算出来的带隙比实验小，导致 ε2 阈值偏低、吸收边红移。如果客户关心可见光吸收，我建议先用 HSE06 或 GW 修正带隙，再用 scissors operator 整体平移介电函数，而不是直接用 PBE。第二个是 K 点密度：光学跃迁涉及 k 空间积分，K 点太稀会丢失尖锐的 Van Hove 奇点。第三个是真空层：如果是 slab 或二维材料，垂直方向介电函数的处理要格外小心，需要加真空方向修正。

我的经验是，光学计算最容易被客户质疑的是峰位和实验对不上。这通常不是计算错了，而是泛函低估带隙。交付时我都会同时给出 PBE 结果和带隙修正后的估计，让客户心里有数。另外，吸收系数 α 的单位是 cm⁻¹，折射率 n 无量纲，这两个量对光伏器件设计最有用。

关于 VASP 光学性质和介电函数计算的更多细节，我整理在 [vasp计算光学](https://www.keyanxueshu.com/category/dft/) 栏目里。[vasp计算光学](https://www.keyanxueshu.com/)

关于实验对比再多说一点：左侧 ε2 的主峰位置通常和 PBE 带隙低估有关。比如图中主峰在 3.5 eV，如果实验测得吸收边在 2.8 eV，两者差 0.7 eV，很可能就是带隙误差。我的处理办法是先用 scissors operator 把介电函数整体 rigid shift 0.7 eV，再和实验对比峰位；峰形（相对强度、半高宽）则由跃迁矩阵元和展宽参数决定。如果 rigid shift 后峰位对上了但强度对不上，就要检查 K 点密度和 CSHIFT。还有一个实用技巧：把 ε1 的零点穿越能量和反射率最低能量对应起来，可以估算等离子体频率，这对金属/半导体光电响应设计很有帮助。