全斯托克斯偏振成像技术：超构透镜模型与远场电场分布的FDTD仿真研究

全斯托克斯偏振成像 超构透镜模型 超表面 FDTD仿真 2018_Arbabi et al_Full-Stokes Imaging Polarimetry Using Dielectric Metasurfaces 主要包括硅纳米柱的单元结构仿真、相位和透射率的参数化扫描，偏振复用超构透镜的偏振解耦合相位计算代码，三组正交偏振成像的超构透镜模型，和对应的远场电场分布计算，以及六像素全斯托克斯超透镜的模型和远场计算结果

硅基超构透镜搞全斯托克斯偏振成像这事，本质上就是玩转纳米柱子的几何参数和相位操控。Arbabi团队2018年的骚操作至今让人眼前一亮——他们用一堆硅纳米柱子当乐高积木，硬是搭出了能同时处理六个偏振态的光学系统。

先说最基础的单元结构仿真。在FDTD Solutions里建模时，每个硅纳米柱的直径D和高度H直接决定相位响应。我习惯用Python脚本批量生成参数矩阵：

import numpy as np params = [(d, h) for d in np.linspace(100, 250, 20) for h in np.linspace(600, 800, 10)]

这种参数遍历要跑上万个仿真，建议搞个分布式任务队列。重点观察532nm波长下的相位延迟曲线是否出现突变——这往往意味着出现高阶模式影响线性度。

全斯托克斯偏振成像 超构透镜模型 超表面 FDTD仿真 2018_Arbabi et al_Full-Stokes Imaging Polarimetry Using Dielectric Metasurfaces 主要包括硅纳米柱的单元结构仿真、相位和透射率的参数化扫描，偏振复用超构透镜的偏振解耦合相位计算代码，三组正交偏振成像的超构透镜模型，和对应的远场电场分布计算，以及六像素全斯托克斯超透镜的模型和远场计算结果

偏振解耦阶段最烧脑。需要同时满足三组正交偏振态的相位需求，代码里其实就是解三个线性方程组：

A = [cos(2*theta_x), sin(2*theta_x), 1; cos(2*theta_y), sin(2*theta_y), 1; cos(2*theta_z), sin(2*theta_z), 1]; b = [phase_x; phase_y; phase_z]; x = A\b; % 关键在这步矩阵除法

这个操作相当于在琼斯矩阵空间里做投影，每个纳米柱的位置都要对应特定的偏振权重。实际调试时经常遇到矩阵奇异的情况，得加个正则化项λ*I来稳定解。

远场计算别直接用FDTD自带的傅里叶变换，内存会爆炸。改用角谱传播算法，核心代码段长这样：

def angular_spectrum(E0, wavelength, z): k = 2*np.pi/wavelength dx = 50e-9 # 网格尺寸 fx = np.fft.fftfreq(E0.shape[0], dx) fy = np.fft.fftfreq(E0.shape[1], dx) transfer = np.exp(1j*z*np.sqrt(k**2 - (2*np.pi*fx)**2 - (2*np.pi*fy)**2)) return np.fft.ifft2(np.fft.fft2(E0) * transfer)

注意处理倏逝波时要加个低通滤波，否则数值会飘。当六像素超透镜的远场出现六个独立聚焦点时，说明偏振解耦成功了，这时候的PSF半高宽通常在1.2λ左右。

实测中发现硅纳米柱的侧壁粗糙度影响比预期大，5nm的起伏就能让消光比掉10dB。后来改用热氧化工艺平滑侧壁，透射率图谱的毛刺明显减少。这种工程细节往往比理论模型更考验耐心，但也正是超表面从论文走向应用的关键坎。