全斯托克斯偏振成像 超构透镜模型 超表面 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。后来改用热氧化工艺平滑侧壁透射率图谱的毛刺明显减少。这种工程细节往往比理论模型更考验耐心但也正是超表面从论文走向应用的关键坎。