电磁波基础：麦克斯韦方程组、波动方程、偏振态、相位与振幅的物理意义
超表面光学本质上是通过亚波长人工结构对电磁波的振幅、相位、偏振和波前进行精准调控。要理解超表面的设计逻辑与物理极限，首先必须牢固掌握电磁波的基本描述框架。本章从麦克斯韦方程组出发，推导波动方程，进而深入剖析偏振态、相位与振幅的物理意义，为后续超表面单元结构设计、相位离散化及偏振复用等核心内容奠定坚实的理论基础。
2.1 麦克斯韦方程组：电磁场的完备描述
麦克斯韦方程组是经典电磁理论的基石，它由四个方程构成，分别描述了电场与电荷、磁场与磁单极子、电磁感应以及电流与磁场之间的关系。在超表面光学中，我们通常处理的是无源、线性、各向同性的介质，此时方程组可以简化为如下形式：
- 高斯定律（电）： ∇·D = 0，表明在无自由电荷区域，电位移矢量的散度为零，电场线是连续的。
- 高斯定律（磁）： ∇·B = 0，说明磁场是无源场，磁力线总是闭合的，不存在孤立的磁单极子。
- 法拉第电磁感应定律： ∇×E = -∂B/∂t，揭示了时变磁场会产生旋涡状电场，这是电磁波传播的核心机制之一。
- 安培-麦克斯韦定律： ∇×H = ∂D/∂t，表明时变电场（位移电流）会产生旋涡状磁场，与法拉第定律共同支撑了电磁波的自我维持传播。
在超表面仿真中，我们通常使用频域求解器（如FDTD或RCWA）来求解这些方程。理解方程中每一项的物理含义，有助于判断仿真结果的合理性。例如，若仿真中出现了非物理的磁场散度不为零，则说明网格剖分或边界条件设置有误。
2.2 波动方程与平面波解
从麦克斯韦方程组出发，对法拉第定律取旋度并代入安培定律，可以推导出无源区域中的波动方程：
∇²E - με ∂²E/∂t² = 0
∇²H - με ∂²H/∂t² = 0其中 μ 和 ε 分别为介质的磁导率和介电常数。波动方程最简单的解是平面波：
E(r,t) = E₀ · e^(j(ωt - k·r))
H(r,t) = H₀ · e^(j(ωt - k·r))式中，ω 为角频率，k 为波矢，其大小 k = ω√(με) = 2π/λ。平面波解表明，电磁波在均匀介质中以速度 v = 1/√(με) 传播。在真空中，v = c ≈ 3×10⁸ m/s。
对于超表面设计而言，平面波解提供了“入射波”和“出射波”的理想化描述。实际超表面仿真中，我们通常用平面波作为激励源，观察透射或反射场的振幅、相位和偏振变化。值得注意的是，平面波是无限延伸的理想模型，而实际光束具有有限宽度和发散角，这在设计大数值孔径超表面透镜时需要特别考虑。
| 物理量 | 符号 | 单位 | 在超表面中的意义 |
|---|---|---|---|
| 介电常数 | ε | F/m | 决定材料对电场的响应，影响谐振频率 |
| 磁导率 | μ | H/m | 决定材料对磁场的响应，超表面中常利用等效μ |
| 波阻抗 | η = √(μ/ε) | Ω | 决定界面反射率，阻抗匹配是高效超表面的关键 |
| 波数 | k = 2π/λ | rad/m | 决定相位累积速率，与单元尺寸直接相关 |
2.3 偏振态：矢量场的几何描述
偏振态描述了电场矢量在垂直于传播方向的平面内的振动轨迹。根据轨迹形状，偏振态可分为线偏振、圆偏振和椭圆偏振。在数学上，偏振态可以用琼斯矢量（Jones Vector）或斯托克斯参数（Stokes Parameters）来完整描述。
对于沿 z 方向传播的平面波，电场矢量可以分解为 x 和 y 两个正交分量：
E_x = E₀x · cos(ωt - kz + φ_x)
E_y = E₀y · cos(ωt - kz + φ_y)当 φ_x = φ_y 时，合成电场为线偏振；当 E₀x = E₀y 且 φ_y - φ_x = ±π/2 时，为圆偏振；其他情况为椭圆偏振。超表面通过设计各向异性单元结构（如矩形纳米柱、椭圆介质柱），可以对 x 和 y 偏振分量施加不同的相位延迟，从而实现偏振转换或偏振复用。
2.4 相位与振幅：超表面调控的核心自由度
相位描述了电磁波在传播过程中某一时刻的振动状态，它决定了波前的形状。振幅则决定了电磁波的能量密度。在传统光学元件中，相位通过光程差（折射率×几何厚度）来累积，而振幅通过吸收或反射来调制。超表面的革命性之处在于，它可以在亚波长尺度内同时独立控制相位和振幅。
对于透射型超表面，单元结构的复透射系数可以表示为：
t = |t| · e^(jφ)其中 |t| 为透射振幅，φ 为透射相位。理想情况下，我们希望 |t| ≈ 1（高透射率）且 φ 可以覆盖 0~2π 的完整范围。实现这一目标的方法包括：
- 谐振相位调控：利用等离子体或Mie谐振，在谐振频率附近相位变化剧烈，但振幅损耗也较大。
- 几何相位调控：通过旋转各向异性单元，对圆偏振光施加 2α 的相位延迟（α 为旋转角），振幅几乎不变。
- 传播相位调控：通过改变单元结构的尺寸或占空比，调节等效折射率，从而累积不同的相位延迟。
2.5 案例拆解：利用相位与偏振调控实现涡旋光束生成
为了将上述概念融会贯通，我们以一个典型的超表面应用——涡旋光束生成器为例进行拆解。涡旋光束具有螺旋波前，其相位分布为 φ(l,θ) = l·θ，其中 l 为拓扑荷数，θ 为方位角。
设计目标：在 λ = 633 nm 波长下，将入射的线偏振高斯光束转换为拓扑荷 l = 1 的涡旋光束，转换效率 > 80%。
关键参数与设计过程：
- 基底材料：熔融石英（n = 1.457），单元结构：非晶硅纳米柱（n = 3.5），高度 H = 600 nm，周期 P = 400 nm。
- 通过扫描纳米柱的直径 D（从 80 nm 到 240 nm），在 FDTD 仿真中获得透射相位 φ 与透射率 |t| 的对应关系。
- 选择透射率 > 90% 且相位覆盖 0~2π 的直径集合，建立“直径-相位”查找表。
- 根据目标相位分布 φ(x,y) = arctan(y/x)，将超表面划分为 20×20 个单元，每个单元选择对应的纳米柱直径。
仿真验证与结果：在 FDTD 中仿真整个超表面，观察透射场的相位分布，确认在中心处存在相位奇点，且围绕中心一周相位变化为 2π。远场强度分布呈现环形，中心光强为零，符合涡旋光束的特征。实测转换效率为 82%，与仿真值 85% 接近。
风险与复盘：在加工过程中，纳米柱的侧壁倾斜角会导致相位偏移，实测效率低于仿真值约 3%。复盘发现，在仿真中未考虑侧壁倾角（约 2°），后续通过引入倾角参数重新优化，将效率提升至 84%。
2.6 总结与常见误区
本章从麦克斯韦方程组出发，系统梳理了电磁波的基本描述框架。核心要点总结如下：
| 概念 | 物理意义 | 超表面设计中的角色 |
|---|---|---|
| 麦克斯韦方程组 | 电磁场基本规律 | 仿真求解的理论基础 |
| 波动方程 | 电磁波传播的数学描述 | 确定传播常数、阻抗 |
| 偏振态 | 电场矢量振动轨迹 | 偏振复用、几何相位 |
| 相位 | 波前形状的决定因素 | 波前整形、聚焦、偏折 |
| 振幅 | 能量密度分布 | 效率优化、旁瓣抑制 |
常见误区包括：将相位与时间延迟混淆、忽略偏振态在非正入射时的变化、以及认为振幅调控不重要。实际上，在超表面设计中，振幅与相位是耦合的，单纯追求相位覆盖而忽略振幅损耗，往往会导致器件效率低下。正确的做法是在设计单元结构时，同时优化透射/反射振幅和相位，必要时采用多目标优化算法。