二、电磁学基础：麦克斯韦方程组、介质中的波动方程、边界条件及光学常数（介电函数）的物理意义
表面等离激元（Surface Plasmon Polaritons, SPPs）本质上是光与金属中自由电子气耦合形成的电磁模式。要理解其激发条件、传播特性与局域增强效应，必须回到电磁学的底层框架。本节从麦克斯韦方程组出发，推导介质中的波动方程，明确界面处的边界条件，并深入剖析介电函数（光学常数）如何决定金属的光学响应。这些内容不仅是后续理解SPP色散关系、穿透深度、共振波长的数学基础，更是设计等离激元器件时选择材料与结构的物理依据。
2.1 麦克斯韦方程组与介质中的波动方程
核心原理：从宏观电磁场到波动方程
表面等离激元的光学行为完全由经典电磁理论描述。麦克斯韦方程组是这一理论的基石，其微分形式为：
- 高斯定律：∇·D = ρ，表明电荷是电位移矢量的源。
- 高斯磁定律：∇·B = 0，自然界不存在孤立磁单极子。
- 法拉第定律：∇×E = -∂B/∂t，时变磁场产生旋涡电场。
- 安培-麦克斯韦定律：∇×H = J + ∂D/∂t，电流和时变电场共同产生旋涡磁场。
对于非磁性、无自由电荷、无外加电流的均匀介质，结合本构关系 D = ε₀εᵣE、B = μ₀μᵣH，对法拉第定律取旋度并代入安培定律，可消去磁场矢量，得到电场满足的波动方程：
∇²E - (εᵣμᵣ/c²) ∂²E/∂t² = 0对于时谐场（e-iωt），上式简化为亥姆霍兹方程：
∇²E + k₀²εᵣμᵣ E = 0其中 k₀ = ω/c 为真空波矢。在表面等离激元问题中，我们通常处理的是分层介质结构（如金属-介质界面），因此需要分别在各均匀层中求解波动方程，再通过边界条件将各层的解耦合起来。
判断方法：如何从波动方程判断模式类型
在分层介质中求解波动方程时，需要根据各层介电常数的符号和大小判断解的性质：
- 若 εᵣ > 0：波动方程的解为振荡型（传播波），对应体材料中的光传播。
- 若 εᵣ < 0：波动方程的解为指数衰减型（倏逝波），场无法在材料内部传播。
- 若界面一侧 ε₁ > 0，另一侧 ε₂ < 0 且 |ε₂| > ε₁：界面处可存在沿界面传播、垂直界面方向指数衰减的表面波——这正是SPP模式。
这一判断准则在后续设计等离激元波导、谐振器时至关重要。例如，在可见光波段，银的 ε ≈ -17 + 0.7i，而空气 ε = 1，满足 |ε_Ag| > ε_air，因此银-空气界面可以支持SPP传播。
2.2 边界条件与表面等离激元的场分布
核心原理：边界条件决定界面电磁场
在两种介质的分界面（设为 z = 0 平面，介质1为金属，介质2为电介质），电磁场必须满足以下边界条件：
| 边界条件 | 数学表达 | 物理含义 |
|---|---|---|
| 电场切向连续 | E₁ₜ = E₂ₜ | 切向电场无跳变 |
| 磁场切向连续（无表面电流） | H₁ₜ = H₂ₜ | 切向磁场无跳变 |
| 电位移法向连续（无表面电荷） | ε₁E₁ₙ = ε₂E₂ₙ | 法向电位移连续 |
| 磁感应法向连续 | B₁ₙ = B₂ₙ | 法向磁感应连续 |
对于横磁（TM）偏振的SPP模式，磁场只有 y 分量（H = H_y ŷ），电场有 x 和 z 分量。在金属（z < 0）和介质（z > 0）中分别写出场的表达式，并应用边界条件，可以得到SPP的色散关系：
k_spp = k₀ √(ε₁ε₂ / (ε₁ + ε₂))其中 ε₁ 为金属介电函数，ε₂ 为介质介电函数。该色散关系是表面等离激元光学中最核心的公式之一，它决定了SPP的传播常数、波长、传播长度和场束缚能力。
案例拆解：银-空气界面的SPP模式分析
考虑波长 λ₀ = 633 nm（He-Ne激光）下，银的介电函数 ε_Ag = -18.3 + 0.5i，空气 ε_air = 1。计算SPP的传播常数：
k₀ = 2π/λ₀ = 9.93 × 10⁶ m⁻¹
k_spp = k₀ √((-18.3+0.5i)×1 / (-18.3+0.5i+1))
= k₀ √((-18.3+0.5i)/(-17.3+0.5i))
≈ k₀ × 1.032 + i × 0.0014
≈ 1.025 × 10⁷ + i × 1.4 × 10⁴ m⁻¹由此可得：SPP波长 λ_spp = 2π/Re{k_spp} ≈ 613 nm，比自由空间波长略短；传播长度 L = 1/(2Im{k_spp}) ≈ 36 μm。这意味着SPP在银-空气界面上可以传播约36微米，之后能量因欧姆损耗衰减至1/e。这一传播长度对于设计等离激元波导、传感器等器件具有直接指导意义。
2.3 介电函数（光学常数）的物理意义
核心原理：介电函数如何描述材料光学响应
介电函数 ε(ω) = ε′(ω) + iε″(ω) 是描述材料对电磁场响应的核心光学常数。其物理意义可以从以下角度理解：
- 实部 ε′(ω)：决定材料对光的极化响应程度，影响光在材料中的传播速度（折射率 n² = (|ε|+ε′)/2）和色散特性。对于金属，ε′ < 0 意味着电子气的极化响应与外加电场反相，导致光无法在金属内部传播。
- 虚部 ε″(ω)：描述材料的吸收损耗，与消光系数 κ 相关（κ² = (|ε|-ε′)/2）。虚部越大，光在材料中传播时的能量损耗越严重。
- 介电函数与折射率的关系：ε = (n + iκ)²，其中 n 为折射率，κ 为消光系数。这一关系将介电函数与实验可直接测量的光学常数联系起来。
对于贵金属（金、银、铜），在可见光和近红外波段，介电函数可以用Drude模型近似描述：
ε(ω) = 1 - ωₚ² / (ω² + iγω)其中 ωₚ 为等离子体频率，γ 为电子散射速率（阻尼系数）。对于银，ωₚ ≈ 9.1 eV，γ ≈ 0.02 eV；对于金，ωₚ ≈ 9.0 eV，γ ≈ 0.07 eV。Drude模型在光子能量低于带间跃迁阈值时（银约3.9 eV，金约2.4 eV）具有较高精度。
判断方法：如何获取和验证介电函数数据
在实际工程中，获取准确的介电函数数据是仿真和设计的前提。常用方法包括：
- 查表法：使用Johnson和Christy（1972年）、Palik（1985年）等经典数据集，或使用在线数据库（如 refractiveindex.info）。
- 椭偏测量法：通过光谱椭偏仪测量薄膜样品的 Ψ 和 Δ 参数，反演得到 n 和 k，进而计算 ε。
- Kramers-Kronig变换：从反射率或吸收谱数据，通过K-K关系计算介电函数的实部和虚部。
- Drude-Lorentz拟合：将实验数据用多个Lorentz振子叠加拟合，获得解析表达式，便于在仿真软件中调用。
需要注意的是，薄膜的介电函数与块体材料存在差异，因为薄膜的晶粒尺寸、缺陷密度、表面氧化层等因素都会影响电子散射速率。因此，在精密设计中，建议对实际使用的薄膜样品进行椭偏测量，而非直接引用块体数据。
案例拆解：金薄膜介电函数的尺寸效应
假设我们需要设计一个工作在785 nm波长的等离激元生物传感器，计划使用50 nm厚的金膜。查阅块体金数据：ε_Au = -23.8 + 1.4i。然而，实际蒸镀的50 nm金膜由于晶界散射和表面粗糙度，其电子散射速率 γ 可能从块体的 0.07 eV 增加到 0.12 eV。根据Drude模型修正：
ε_film(ω) = 1 - ωₚ²/(ω² + iγ_film ω)
其中 γ_film = γ_bulk + A·v_F/d式中 A 为与散射机制相关的常数（约0.5-1），v_F 为费米速度（金约1.4×10⁶ m/s），d 为薄膜厚度。代入数值：
γ_film ≈ 0.07 + 0.75 × 1.4×10⁶ / 50×10⁻⁹ ≈ 0.07 + 0.021 eV ≈ 0.091 eV修正后的介电函数虚部约为 1.8i，比块体值增大约30%。这意味着SPP传播损耗增加约30%，传感器的品质因子（Q = Re{k_spp}/Im{k_spp}）相应下降。这一案例说明，在实际器件设计中，必须考虑薄膜厚度对介电函数的影响，否则仿真结果与实验测量会出现显著偏差。
常见误区与总结复盘
- 误区一：将介电函数视为常数。实际上 ε(ω) 强烈依赖于频率，尤其在金属中，从可见光到近红外波段，ε′ 的值随波长增大而增大。设计宽光谱器件时，必须考虑色散效应。
- 误区二：忽略虚部的影响。有些初学者在理论分析中只关注 ε′，将金属视为无损耗介质。这会导致SPP传播长度被高估数个数量级，无法解释实验中的实际损耗。
- 误区三：混淆介电函数与折射率。虽然两者可以互相换算，但在边界条件中必须使用介电函数 ε，而非折射率 n。因为边界条件中的法向分量关系是 ε₁E₁ₙ = ε₂E₂ₙ，而不是 n₁E₁ₙ = n₂E₂ₙ。
- 误区四：忽视表面粗糙度和氧化层。实际金属表面往往存在几纳米厚的自然氧化层（如铝表面的Al₂O₃），这会改变SPP的色散关系，导致共振峰偏移。在精密测量中，需考虑这一影响。