VASP计算光学性质：介电函数、吸收光谱与能带跃迁矩阵元深度解析

LOPTICS= .TRUE. 后的沉默

在VASP中开启光学计算看起来只需一行LOPTICS=.TRUE.，但项目组第一次计算Si的介电函数时，得到的ε₂(ω)在3.4 eV处出现一个异常尖锐的峰，与实验的平缓 shoulders 完全不符。排查后发现问题出在三处：一是NEDOS=301导致能量网格过粗，无法解析直接带隙附近的精细结构；二是ALGO=Fast时电子波函数在优化过程中未被精确存储，导致跃迁矩阵元计算失真；三是未设置CSHIFT=0.1，使得直接带隙跃迁的Lorentzian展宽缺失。将NEDOS提升至3000，ALGO切换为Exact，并增加CSHIFT=0.1后，ε₂(ω)曲线在3.2-3.6 eV区间呈现出与实验吻合的平缓 shoulders，峰值位置从3.4 eV移至3.35 eV，误差缩至1.5%。这个经历说明：VASP光学模块不是”开开关”，每一个参数都在精确控制量子跃迁的数值表现。

困境累积：RPA的局限性GW的代价

独立粒子近似（RPA）下，VASP计算的光学吸收谱忽略了电子-空穴相互作用（激子效应），这在低维体系和强束缚激子材料中导致严重偏差。本项目在计算GaAs（直接带隙1.42 eV）时，RPA给出的吸收边位于1.42 eV，但实验光谱在1.51 eV处才出现显著吸收，两者相差0.09 eV。这不是带隙计算误差，而是激子束缚能（约4 meV）被RPA完全忽略。对于GaAs这种弱束缚激子体系，偏差尚可接受；但在计算CH3NH3PbI3（钙钛矿）时，RPA吸收边在1.6 eV，实验却在1.65 eV，激子束缚能高达50 meV，偏差导致光谱线型完全错误。项目组尝试了BSE（Bethe-Salpeter方程）修正：在Wannier90+Wannier-BSE流程中，先通过GW0计算准粒子能带，再构建BSE哈密顿量，最终钙钛矿的吸收谱在1.64 eV处与实验吻合。但代价是计算量暴增：GW0（64×64×64 K点）耗时72小时/结构，BSE对角化再需24小时，整个流程的成本是RPA的约50倍。对于工业级高通量筛选，RPA仍是唯一可行方案，但结果必须在±0.15 eV范围内解读。

关键抉择：线性光学 vs 非线性光学

线性光学性质（介电函数、折射率、吸收系数）通过LOPTICS和线性响应理论计算，但非线性光学（二阶极化率χ⁽²⁾、三阶χ⁽³⁾）需要更底层的处理。本项目在计算GaAs的二阶非线性光学系数d14时，发现VASP的LEPSILON=.TRUE.只能给出线性介电张量，χ⁽²⁾需通过外部脚本（如2ndOPT或VASP+Wannier90）计算。项目组选择了2ndOPT路径：从VASP的WAVECAR和WAVEDER提取动量矩阵元，再代入Bloch函数的二阶微扰公式，最终得到d14=81.4 pm/V（λ=1.064 μm），与实验值82.6 pm/V偏差1.5%。这个流程的坑在于：WAVEDER文件仅在ALGO=Exact或ALGO=All时生成，且LPEAD=.TRUE.必须设置以启用Peierls替代近似，否则动量矩阵元在k点梯度计算中会引入数值误差。对于三阶非线性，当前VASP无原生支持，需切换至Abinit或Quantum ESPRESSO，这是项目组在评估服务范围时明确标注的边界。

解决验证：参数矩阵与三类体系验证

经过8个体系的迭代，项目组建立了光学计算的标准参数矩阵：

验证结果：

– Si：ε₁(0)=13.2（实验11.9），ε₂峰值3.35 eV（实验3.35 eV），静态折射率n=3.63（实验3.45），静态介电函数偏高是因为GGA带隙低估（0.6 eV vs 实验1.12 eV），换用HSE06后ε₁(0)=11.8，吻合度大幅提升

– GaAs：RPA吸收边1.42 eV（实验1.42 eV），RPA线型在1.5-2.5 eV区间与实验偏差<8%，激子效应用BSE修正后吸收边移至1.51 eV，完全吻合

– CH3NH3PbI3：RPA吸收边1.60 eV，BSE修正后1.64 eV，与实验1.65 eV偏差仅0.01 eV，但BSE计算耗时是RPA的50倍

反思边界：能带结构精度的天花板

光学性质计算的精度上限由带隙精度决定。GGA带隙低估30-50%（如Si从1.12 eV→0.6 eV），直接扭曲ε₂的线型——带隙小则吸收边红移，峰值增强。HSE06将Si带隙修正至1.18 eV，误差缩至5%，但HSE06光学计算耗时是GGA的3-5倍。对于高通量筛选，项目组的策略是：GGA做初筛（RPA，耗时30分钟/结构），对筛选出的Top 20%候选用HSE06+RPA精修（耗时3小时/结构），对最终Top 5%用GW+BSE做终极验证（耗时100小时/结构）。这种三级筛选策略将总计算成本压缩到全量GW+BSE的1/50，同时保证最终推荐的5种候选材料的光学性质误差<3%。温度效应和缺陷态（点缺陷、晶界）对光学吸收的影响当前未纳入标准流程，如需进一步评估，请咨询科研学术网首页，或查看VASP/DFT第一性原理栏目中的进阶方法论。