4. 材料建模与色散处理

在FDTD仿真中，材料建模的准确程度直接决定结果的物理可信度。时域有限差分法在时域推进电场与磁场，因此材料响应必须以时域形式参与迭代。对于无色散介质，介电常数是一个实数常量，实现简单；但对于金属、等离子体、极性液体、半导体以及大多数光学材料，介电常数随频率变化，即存在色散。若忽略色散，仿真在宽带激励下会出现明显的波形失真、相位误差甚至非物理结果。本节围绕介电常数模型、Drude模型、Lorentz模型、Debye模型及拟合方法展开，重点说明各模型的物理来源、时域实现方式以及工程中如何选择与拟合参数。

一、介电常数模型的基本框架

在频域中，材料的本构关系可写为 D(ω) = ε₀εᵣ(ω)E(ω)，其中 εᵣ(ω) 为复相对介电常数。色散模型的核心任务，是用有理函数或极点形式逼近 εᵣ(ω)，使其在FDTD中能够转化为辅助微分方程或辅助变量递推。常见的三类模型为：

• Drude模型：描述自由载流子响应，适用于金属、重掺杂半导体、等离子体。其介电函数为 εᵣ(ω) = ε∞ − ωp²/(ω² + iγω)，其中 ωp 为等离子体频率，γ 为碰撞频率。
• Lorentz模型：描述束缚电子谐振响应，适用于介电材料、激子吸收、光学谐振。其形式为 εᵣ(ω) = ε∞ + Δε·ω₀²/(ω₀² − ω² − iγω)，ω₀ 为谐振频率，Δε 为振子强度。
• Debye模型：描述极性分子的弛豫响应，适用于水、生物组织、聚合物等。其形式为 εᵣ(ω) = ε∞ + Δε/(1 + iωτ)，τ 为弛豫时间。

三者的物理机制不同：Drude对应自由电子集体振荡，Lorentz对应束缚电子受迫振动，Debye对应偶极子取向弛豫。在FDTD中，它们都可以通过引入辅助极化变量 P 或电流密度 J，将频域关系转换为时域一阶微分方程，再以显式或半隐式格式递推。

二、Drude模型：金属与等离子体的时域实现

Drude模型将金属中的自由电子视为在电场作用下加速、并与晶格碰撞损失能量的经典气体。其频域介电函数为 εᵣ(ω) = ε∞ − ωp²/(ω² + iγω)。在光频段，金属的 ε∞ 通常取1，ωp 和 γ 由材料决定。例如银在可见光附近常用 ωp ≈ 1.37×10¹⁶ rad/s、γ ≈ 3.2×10¹³ rad/s 进行近似，但实际工程中需根据目标波段重新拟合。

在FDTD中，Drude模型通过引入自由电流密度 J 实现。由频域关系可得时域方程：

dJ/dt + γJ = ε₀ωp²E

采用中心差分或半隐式格式离散，得到递推式：

J^{n+1} = α J^n + β (E^{n+1} + E^n)/2
α = (1 − γΔt/2)/(1 + γΔt/2)
β = ε₀ωp²Δt/(1 + γΔt/2)

随后电场更新中，D 与 E 的关系变为 D = ε₀ε∞E + Δt·J 的累积形式。实际实现时，通常将 J 与 E 交错半个时间步，以保证二阶精度。判断Drude实现是否正确，可检查金属在低于 ωp 的频率下是否呈现负实部介电常数，以及反射率是否接近1。

三、Lorentz模型：谐振吸收与色散

Lorentz模型将束缚电子视为受恢复力和阻尼力作用的谐振子。其频域形式为 εᵣ(ω) = ε∞ + Δε·ω₀²/(ω₀² − ω² − iγω)。当 ω 接近 ω₀ 时，介电常数出现共振峰，对应材料吸收带。在FDTD中，引入极化强度 P，时域方程为：

d²P/dt² + γ dP/dt + ω₀²P = Δε·ε₀ω₀²E

将其拆为一阶系统，令 J = dP/dt，则：

dP/dt = J
dJ/dt = Δε·ε₀ω₀²E − γJ − ω₀²P

时间离散后，P 和 J 均可用显式格式更新。需要注意的是，当 ω₀Δt 较大时，显式格式可能不稳定，此时应减小时间步或采用半隐式处理。Lorentz模型常用于描述二氧化硅、氮化硅等介质在紫外波段的吸收，也可用于激子共振建模。

四、Debye模型：弛豫与宽频响应

Debye模型描述极性分子在电场作用下的取向弛豫。其频域形式为 εᵣ(ω) = ε∞ + Δε/(1 + iωτ)。与Lorentz不同，Debye没有谐振峰，而是呈现从静态介电常数到高频介电常数的单调过渡。时域方程为：

τ dP/dt + P = Δε·ε₀E

离散后得到：

P^{n+1} = (τ/(τ+Δt)) P^n + (Δε·ε₀·Δt/(τ+Δt)) E^{n+1}

Debye模型广泛用于水、生物组织、聚合物等。以水为例，室温下 τ ≈ 8.2 ps，Δε ≈ 78.4，ε∞ ≈ 5.2。在太赫兹和微波频段，水的Debye弛豫是影响传播损耗的关键因素。若仿真频段跨越多个弛豫过程，可采用多极Debye模型叠加。

五、拟合方法与工程实践

实际材料参数往往来自实验测量的折射率 n 和消光系数 k，或椭圆偏振光谱数据。拟合的目标是找到一组模型参数，使模型计算的 εᵣ(ω) 与实验数据在目标频段内误差最小。常用流程如下：

1. 将实验 n、k 转换为复介电常数：ε₁ = n² − k²，ε₂ = 2nk。
2. 选择模型类型与极点数量。金属优先Drude，介质吸收峰优先Lorentz，宽频弛豫优先Debye。
3. 设定目标函数，例如在频段内对 ε₁ 和 ε₂ 的加权最小二乘误差。
4. 采用Levenberg-Marquardt、遗传算法或粒子群算法优化参数。
5. 验证拟合结果是否满足Kramers-Kronig关系，确保因果性。
6. 在FDTD中设置时间步 Δt，检查 ωpΔt、ω₀Δt、Δt/τ 是否满足稳定性条件。

总结来说，材料建模与色散处理是FDTD仿真从“能跑”到“可信”的关键一步。Drude、Lorentz、Debye分别对应自由载流子、束缚谐振和偶极弛豫三种物理机制，它们的时域实现都归结为辅助微分方程的递推。工程中应根据材料类型和目标频段选择模型，并通过实验数据拟合参数，最后用稳定性条件和物理一致性检验确保仿真结果可靠。