2. 介电函数与光学常数的物理基础

在DFT框架下讨论光学性质，核心出发点是复介电函数 ε(ω)=ε₁(ω)+iε₂(ω)。它把电子结构（能带、态密度、跃迁矩阵元）与宏观可测量的光学量连接起来。理解这一节的关键，不是记住一堆公式，而是弄清一条“数据流”：从能带与波函数出发，先得到虚部 ε₂，再由Kramers-Kronig关系得到实部 ε₁，最后由 ε₁、ε₂派生出折射率 n、消光系数 k、吸收系数 α、反射率 R 和能量损失函数 L。只要这条链路清楚，后续VASP计算中读什么文件、改什么参数、结果如何判读都会变得明确。

一、复介电函数的两个分量

介电函数的虚部 ε₂(ω) 描述电子从占据态到未占据态的跃迁吸收，其表达式为：

ε₂(ω) = (4π²e²/Ω) Σ_{c,v,k} w_k |<c,k|u·r|v,k>|² δ(E_c - E_v - ħω)

其中 Ω 是原胞体积，w_k 是k点权重，<c,k|u·r|v,k> 是跃迁矩阵元，δ 函数保证能量守恒。这个式子说明：ε₂ 的峰位对应可能的带间跃迁能量，峰强取决于该能量处有多少对占据态与空态、以及跃迁是否被偶极选择定则允许。因此 ε₂ 谱本质上是“联合态密度 × 跃迁强度”的加权结果。

实部 ε₁(ω) 不能独立给定，它由 ε₂ 通过Kramers-Kronig变换得到：

ε₁(ω) = 1 + (2/π) P ∫₀^∞ [ω' ε₂(ω') / (ω'² - ω²)] dω'

P 表示柯西主值积分。这一关系意味着：ε₁ 在某一频率的值，取决于 ε₂ 在全频段的分布，而不是只由该频率附近的跃迁决定。实际计算中，ε₂ 通常由DFT直接算出，ε₁ 则由程序内部做KK变换得到，因此 ε₁ 的静态值 ε₁(0) 可以作为检验计算是否合理的一个指标——它应与材料的静态介电常数实验值量级一致。

二、从介电函数到光学常数

有了 ε₁ 和 ε₂，其余光学量都是代数派生。复折射率定义为 N = n + ik，满足 N² = ε，因此：

n(ω) = sqrt( (sqrt(ε₁²+ε₂²) + ε₁) / 2 )
k(ω) = sqrt( (sqrt(ε₁²+ε₂²) - ε₁) / 2 )

其中 n 是折射率，k 是消光系数。k 不为零表示介质对光有衰减，k 越大衰减越强。吸收系数 α 与 k 的关系为 α = 2ωk/c = 4πk/λ。反射率 R 在正入射下由菲涅尔公式给出：

R(ω) = [(n-1)² + k²] / [(n+1)² + k²]

能量损失函数 L(ω) = Im(-1/ε) = ε₂/(ε₁²+ε₂²)，它描述电子穿过材料时损失能量的谱分布。当 ε₁ 和 ε₂ 同时趋于零时，L 出现峰，对应体等离激元共振。因此 L 谱的峰位常被用来判断材料的等离激元能量，也是检验介电函数是否自洽的一个敏感量。

三、与电子结构的内在联系

上述公式中真正决定谱形的是电子结构。具体来说，ε₂ 的峰来自特定k点上的带间跃迁。例如，对半导体，吸收边位置由最小直接带隙决定；对金属，低能区由带内Drude贡献主导，ε₁ 在低频为负，反射率接近1。若把DFT计算得到的能带与 ε₂ 谱对照，可以定位每个峰对应的跃迁路径，例如 Γ 点价带顶到导带底的跃迁、或 X 点附近的跃迁。这种“谱-带对应”分析是判断光学性质计算是否合理的重要手段。

需要特别注意的是，标准DFT（如LDA、GGA）通常会低估带隙，导致 ε₂ 的吸收边偏小、ε₁(0) 偏大。若要与实验对比，常需引入剪刀修正（scissor shift）或使用杂化泛函。此外，计算中必须包含足够的空带数目，否则高能区 ε₂ 不完整，KK变换得到的 ε₁ 也会失真。一般建议空带数至少为占据带数的2～3倍，并在收敛测试中检查 ε₁(0) 是否稳定。

四、实际计算中的操作要点

在VASP中，介电函数计算通常分两步：先做自洽计算得到收敛的电荷密度，再用独立粒子近似（LOPTICS=.TRUE.）或含时DFT（ALGO=TDHF）计算 ε₂。关键INCAR参数包括：

LOPTICS = .TRUE.    ! 开启光学性质计算
CSHIFT  = 0.1       ! 复数位移，避免分母为零
NBANDS  = 96        ! 空带数，需覆盖目标能量范围
NEDOS   = 2000      ! 态密度和介电函数能量网格点数

计算完成后，OUTCAR中会输出 ε₁、ε₂ 以及派生量。检查时先看 ε₂ 是否在低能区非负、高能区是否衰减到零；再看 ε₁(0) 是否与材料静态介电常数合理接近；最后看 L 谱是否出现合理的等离激元峰。若 ε₁(0) 异常大或 L 谱出现非物理尖峰，通常意味着空带不足或k点采样不够密。

总结来说，复介电函数是连接电子结构与宏观光学响应的枢纽。ε₂ 由跃迁直接决定，ε₁ 由KK关系约束，n、k、α、R、L 都是 ε₁ 和 ε₂ 的代数组合。掌握这条数据流，就能在计算中快速定位问题：谱形不对，先查能带和跃迁矩阵元；ε₁(0) 不对，先查空带收敛和带隙修正；派生量异常，先检查 ε₁、ε₂ 是否自洽。这样，后续更复杂的介电与光学计算就有了稳固的物理基础。