3. 光与物质的相互作用：介电函数、极化率、吸收与色散，以及洛伦兹模型和德鲁德模型

光在介质中传播时，其速度、振幅、相位和偏振状态都会发生改变，这些宏观现象的背后是光与物质内部带电粒子（电子、离子）的微观相互作用。本章的核心任务是建立从微观振子模型到宏观光学常数的桥梁。我们将从介电函数和极化率的定义出发，推导出描述介质光学响应的核心参数，并深入剖析两个最基础的理论模型——洛伦兹模型（适用于绝缘体和半导体带间跃迁）与德鲁德模型（适用于金属和掺杂半导体中的自由载流子）。理解这两个模型，是后续分析微纳结构中光场调控、等离激元效应和超表面设计的前提。

上图展示了从微观振子到宏观光学常数的完整逻辑链。入射光场驱动介质内部的束缚电子或自由电子振荡，这些振荡的幅度和相位响应决定了介质的复介电函数，进而推导出折射率、消光系数等宏观参数。下面我们分模块展开。

3.1 核心原理：介电函数与极化率的物理内涵

当电磁波作用于介质时，电场分量会使介质中的正负电荷中心发生相对位移，形成电偶极矩。单位体积内的电偶极矩总和称为极化强度矢量 P。在线性、各向同性的介质中，极化强度与电场强度成正比：

P = ε₀ · χₑ · E

其中，ε₀ 是真空介电常数（8.854 × 10⁻¹² F/m），χₑ 是电极化率（无量纲）。电位移矢量 D 与电场 E 的关系为：

D = ε₀ · E + P = ε₀ · (1 + χₑ) · E = ε₀ · εᵣ · E

这里 εᵣ = 1 + χₑ 就是相对介电函数。在光频段，由于介质对电磁波的响应存在相位延迟，介电函数通常是复数：

ε(ω) = ε₁(ω) + i · ε₂(ω)

实部 ε₁ 描述介质对光的极化能力（决定光速和折射率），虚部 ε₂ 描述介质对光的吸收损耗。两者通过克拉默斯-克勒尼希关系相互联系——这意味着只要知道全频段的吸收谱，就可以计算出折射率的色散曲线。

从微观角度看，介电函数是介质内部各种极化机制（电子极化、离子极化、偶极子转向极化）在特定频率下的综合响应。在光频段（10¹⁴ ~ 10¹⁵ Hz），离子极化和偶极子极化因质量较大无法响应如此快的电场变化，只有质量极轻的电子能够跟随光频电场振荡。因此，微纳光学中我们主要关注电子极化机制。

3.2 洛伦兹模型：束缚电子的谐振响应

洛伦兹模型将介质中的电子视为被弹性力束缚在平衡位置附近的经典谐振子。电子受到三种力的作用：弹性恢复力（-mω₀²x）、阻尼力（-mγ·dx/dt）和电场驱动力（-eE）。根据牛顿第二定律，电子的运动方程为：

m · d²x/dt² + mγ · dx/dt + mω₀² · x = -e · E₀ · exp(-iωt)

求解该方程得到电子的稳态位移 x(ω)，进而得到单个电子的偶极矩 p = -e·x。若单位体积内有 N 个电子，则极化强度 P = N·p，代入 P = ε₀·χₑ·E 可得：

ε(ω) = 1 + (N·e²)/(ε₀·m) · 1/(ω₀² - ω² - iγω)

定义等离子体频率 ωₚ² = N·e²/(ε₀·m)，则洛伦兹模型的介电函数简化为：

ε(ω) = 1 + ωₚ²/(ω₀² - ω² - iγω)

当入射光频率 ω 接近共振频率 ω₀ 时，虚部 ε₂ 达到最大值，产生强烈的吸收峰；实部 ε₁ 在共振频率附近发生剧烈变化——低于共振频率时 ε₁ 较大（正常色散），高于共振频率时 ε₁ 迅速减小（反常色散）。远离共振频率时，ε₁ 趋近于 1，介质表现为透明。

3.2.1 洛伦兹模型的典型参数与行为

实际介质中往往存在多个共振频率（对应不同能级的电子跃迁），完整的介电函数是多个洛伦兹振子的叠加：

ε(ω) = 1 + Σⱼ ωₚⱼ²/(ω₀ⱼ² - ω² - iγⱼω)

例如，水在红外区有多个吸收峰（对应O-H键振动），在紫外区有电子跃迁吸收峰，因此水的折射率色散曲线需要用多个洛伦兹项来拟合。

3.3 德鲁德模型：自由电子的输运响应

对于金属和重掺杂半导体，价电子不再被束缚在特定原子周围，而是在晶格中自由运动。德鲁德模型将这些自由电子视为无恢复力的电子气，运动方程中弹性恢复力项为零（ω₀ = 0）：

m · d²x/dt² + mγ · dx/dt = -e · E₀ · exp(-iωt)

求解后得到德鲁德模型的介电函数：

ε(ω) = 1 - ωₚ²/(ω² + iγω) = 1 - ωₚ²/(ω² + γ²) + i · (ωₚ²·γ)/(ω·(ω² + γ²))

其中 γ = 1/τ，τ 是电子的平均散射时间（弛豫时间）。将实部和虚部分开：

ε₁(ω) = 1 - ωₚ²/(ω² + γ²)

ε₂(ω) = ωₚ²·γ/(ω·(ω² + γ²))

当 ω < ωₚ 时，ε₁ 为负值，这意味着金属对光具有强烈的反射特性（如铝在可见光区的反射率超过90%）；当 ω > ωₚ 时，ε₁ 变为正值，金属变得透明（如紫外区的铝）。这个临界频率 ωₚ 称为等离子体频率，是金属光学性质的分水岭。

3.4 案例拆解：金、银、铝在可见光区的光学响应对比

让我们通过三种典型金属的德鲁德模型参数来理解它们在微纳光学中的应用差异。下表列出了室温下三种金属的德鲁德参数：

案例分析过程：假设我们需要设计一个工作在 633 nm（红光）波长的等离激元纳米天线。首先计算该频率对应的光子能量 E = hc/λ ≈ 1.96 eV。对于金，ω = 1.96 eV，ωₚ = 9.0 eV，γ = 0.07 eV，代入德鲁德公式：

ε₁(Au) = 1 - 9.0²/(1.96² + 0.07²) ≈ 1 - 81/3.84 ≈ -20.1

ε₂(Au) = 9.0² × 0.07/(1.96 × (1.96² + 0.07²)) ≈ 81 × 0.07/(1.96 × 3.84) ≈ 0.75

对于银，γ 更小（0.02 eV），计算得 ε₁(Ag) ≈ -21.9，ε₂(Ag) ≈ 0.21。对比可见，银的虚部远小于金，意味着银纳米天线的欧姆损耗更低，共振品质因子 Q 更高。但银在空气中容易氧化，需要额外的保护层。铝的 ωₚ 更高，在可见光区 ε₁ 的值更大，但其 γ 也较大，损耗相对较高。

设计决策：在需要高灵敏度的传感应用中，优先选择银（损耗低，共振峰尖锐）；在需要化学稳定性或表面功能化的场景中，选择金（虽然损耗稍高，但生物相容性好）；在紫外波段（λ < 400 nm），铝是唯一可行的等离激元材料。

3.5 吸收与色散的联系：克拉默斯-克勒尼希关系

介电函数的实部和虚部并非独立，它们通过克拉默斯-克勒尼希关系（K-K关系）相互关联：

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

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

其中 P 表示柯西主值积分。这个关系的物理意义是：只要知道介质在全频段的吸收谱（ε₂），就可以唯一确定其色散曲线（ε₁）。反之亦然。这在实际工作中非常有用——例如，通过测量薄膜的透射光谱得到吸收系数，再利用K-K关系计算出折射率的色散。

在微纳光学器件设计中，K-K关系也意味着我们不能独立地优化吸收和色散——降低损耗（减小ε₂）必然会影响折射率的频率依赖性。例如，在超表面设计中，为了获得大的相位调控范围，需要选择低损耗材料（如TiO₂、GaN），因为高损耗材料的ε₂会导致K-K积分贡献增大，使ε₁的色散变得平缓，限制了相位调制效率。

3.6 常见误区与总结复盘

• 误区一：将洛伦兹模型和德鲁德模型视为完全独立的理论。实际上，德鲁德模型是洛伦兹模型在 ω₀ = 0 时的特例。更完整的描述是洛伦兹-德鲁德模型，它同时包含束缚电子项和自由电子项，适用于过渡金属和重掺杂半导体。
• 误区二：忽略阻尼项 γ 的频率依赖性。在低温或高纯度材料中，γ 可能随频率变化。对于金属纳米颗粒，当颗粒尺寸小于电子平均自由程时，γ 会因表面散射而增大，导致吸收峰展宽。
• 误区三：将介电函数的实部为负等同于“金属性”。某些介电材料在声子共振区（Reststrahlen带）也会出现 ε₁ < 0 的现象（如SiC在10-12 μm波段），此时材料表现出“光学金属”行为，支持表面声子极化激元。
• 误区四：混淆相速度和群速度的色散。洛伦兹模型在共振频率附近会出现反常色散（dn/dω < 0），此时群速度可能超过相速度甚至为负，但信息传播速度仍受光速限制。

总结复盘：本章的核心是建立“微观振子参数（ω₀, γ, ωₚ）→ 复介电函数 ε(ω) → 宏观光学常数（n, k, R, α）”的完整链路。洛伦兹模型适用于描述束缚电子的共振吸收，德鲁德模型适用于描述自由电子的输运行为。在实际微纳光学设计中，我们需要根据材料类型和波段选择合适的模型，并注意模型的适用范围和局限性。对于定量分析，建议结合实验数据或第一性原理计算，而非仅依赖理想化模型。