VASP计算介电函数 — 光学性质DFT计算的参数调优与精度验证

光学性质计算在VASP中似乎很简单——加个LOPTICS=.TRUE.，多设几个空带，等着收结果就完了。但如果你真的这么做过，大概率会遇到这种情况：算出来的吸收系数比实验值蓝移了0.5 eV以上；介电函数虚部的谱峰形状和实验完全对不上；或者计算出的静态介电常数ε₁(0)比实验值差20%。VASP计算介电函数这件事，从参数设置到物理近似，至少有五个会影响定量结果的环节。本文从工程角度逐一拆解。

一、介电函数计算的理论基础

频率依赖复介电函数ε(ω)=ε₁(ω)+iε₂(ω)是VASP光学性质计算的核心输出。虚部ε₂(ω)通过直接跃迁的Fermi黄金规则计算：

ε₂(ω) = (4π²e²/Ω) × lim(q→0) (1/q²) × Σ_{c,v,k} 2w_k δ(E_ck−E_vk−ħω) × |⟨u_{ck+q}|u_{vk}⟩|²

其中求和遍及所有k点、所有价带v和导带c。跃迁矩阵元⟨u_{ck+q}|u_{vk}⟩描述从价带到导带的光学跃迁概率，δ函数确保能量守恒。

实部ε₁(ω)通过Kramers-Kronig变换从虚部推导：

ε₁(ω) = 1 + (2/π) × P∫₀^∞ [ε₂(ω′)·ω′/(ω′²−ω²)] dω′

其中P表示Cauchy主值积分。一旦有了ε₁(ω)和ε₂(ω)，所有其他光学量都可以推出来：吸收系数α(ω)=ω×ε₂/(n×c)，折射率n(ω)=√[(|ε|+ε₁)/2]，消光系数κ(ω)=√[(|ε|−ε₁)/2]，反射率R(ω)=|(ñ−1)/(ñ+1)|²。

二、NBANDS：空带数是最重要也最容易被忽视的参数

VASP计算介电函数时，NBANDS参数的设置直接决定了光学计算的质量。原因是：介电函数需要计算从占据态到非占据态的跃迁，而非占据态的数量由NBANDS决定。如果NBANDS不够大（空带太少），高能量的光学跃迁就无法被计算，导致高能区的ε₂(ω)被人为截断。

实际经验：NBANDS至少是价带数（占据态）的3倍。对于包含32个电子的半导体体系，价带数约16条，NBANDS至少48条。如果想计算到20 eV能量的光学性质（覆盖真空紫外区），需要NBANDS达到价带数的4-5倍。

但NBANDS也不能无限增大——每增加一条空带，计算量和内存都线性增加。一个实用的判断方法：跑一次测试计算，检查OUTCAR末尾”optical properties”部分的最高跃迁能量是否覆盖了你的目标能量范围。如果最高跃迁能量低于你想要的ε₂(ω)最大能量，增加NBANDS再跑。

三、k点密度与展宽：光谱精细结构的控制

光学性质对k点密度的要求远高于总能量计算——因为跃迁矩阵元在k空间的变化比电子密度剧烈得多。对于Si这样的简单半导体，6×6×6 k点就能给出合理的带隙和吸收谱；但对钙钛矿这种带边有强色散的体系，8×8×8甚至12×12×12才够。

k点不够密的最明显症状是：ε₂(ω)的谱峰被人为平滑，精细结构（如E₁、E₂跃迁肩峰）被抹平。这在和实验椭偏光谱对比时尤其致命——实验谱上的特征峰在你的计算中消失了。

CSHIFT参数控制介电函数计算中δ函数的展宽（替代为Lorentz函数半宽）。默认CSHIFT=0.1 eV对大多数体系适用。减小CSHIFT可以让谱峰更尖锐（接近真实线形），但需要更密的k点来补偿——否则谱峰会因为k点不足出现锯齿状噪声。CSHIFT与k点密度的关系是：CSHIFT越小，需要的k点越密。

四、独立粒子近似与多体效应的补救

VASP的默认光学计算基于独立粒子近似（IPA）——电子和空穴被视为独立运动，不考虑它们之间的库仑相互作用（激子效应）。这对于宽禁带半导体（GaN、ZnO）或金属是合理的，但对窄带隙半导体、有机半导体和低维材料，激子效应可能大到0.3-1.0 eV，不能用IPA忽略。

Bethe-Salpeter方程（BSE）是VASP中处理激子效应的标准方法，但计算量比IPA大1-2个数量级。对于大多数工程应用（而不是基础物理研究），如果IPA+Tauc plot方法给出的光学带隙与实验一致，不需要上BSE。

另一个实际中有用的近似是scissor算符修正——手动将计算出的PBE带隙平移到实验带隙值，然后将导带整体上移同样数量。这不能替代更精确的泛函（如HSE06），但在定性讨论吸收谱形状和相对趋势时是一个实用的工程折中。

五、光学性质的衍生计算

有了ε(ω)，VASP可以输出很多实际可用的光学量。实际上VASP的OUTCAR中直接包含了大部分光学量的频率分布，只需在INCAR中设置即可输出完整数据集。

吸收系数的应用最为直接——与实验UV-Vis吸收光谱对比。但有一个容易混淆的点：实验吸收系数α(ω)通常用比尔-朗伯定律定义（从透射率推导），而VASP输出的是电动力学定义（从介电函数推导）。在带隙以下的弱吸收区，两者的差异可能很大。

折射率的频率分布n(ω)在光学薄膜和光伏器件的抗反射涂层设计中直接有用。静态折射率n(0)通常与实验值吻合较好（偏差<5%）——因为n(0)=√ε₁(0)，主要受带间跃迁的加权平均影响，对带隙的具体值不太敏感。

六、实际案例：TCO透明导电氧化物的光学分析

以ITO（Sn-doped In₂O₃）为例，这是一个典型的需要VASP计算介电函数的体系。ITO在可见光区透明（3.1-4.1 eV），但近红外区有自由载流子吸收（等离子体反射边）。

用VASP计算ITO的光学性质，需要在标准的带间介电函数基础上额外考虑Drude项（自由载流子贡献）。VASP本身不直接输出Drude模型参数，需要从能带结构中提取等离子体频率ω_p（从带内跃迁的联合态密度积分的低频极限估算），然后在后处理中加上Drude项：

ε(ω) = ε_interband(ω) − ω_p²/(ω² + iωγ)

加上Drude项后，近红外区的反射率提升和可见光区的透明度保持可以同时重现——这就是ITO作为透明电极的光学物理机制。

配图建议：

1. 介电函数频谱图（ALT：”VASP计算介电函数ε1ε2频谱半导体光学性质实部虚部”）：双层图，上为ε₁(ω)、下为ε₂(ω)，标注带隙起始位置和主要跃迁峰。
2. 计算vs实验吸收光谱对比图（ALT：”VASP光学吸收系数DFT计算与UV-Vis实验吸收光谱对比验证”）：计算吸收系数与实验UV-Vis谱叠加对比。
3. 光学常数合集图（ALT：”VASP计算折射率消光系数反射率频谱光学常数”）：展示n(ω)、κ(ω)、R(ω)的频谱曲线。