这个标题读起来像实验记录但其实是个非常典型的全链路项目从物理设计到纳米结构再到FDTD仿真验证。长波红外、超构透镜、偏振复用聚焦三个关键词叠在一起做的人知道这里每一步都有坑不做的人可能觉得不就是“跑个仿真”。我这次把这个“仿真之旅”完整拆开从方案选型、单单元扫描、圆偏振光源设置、全透镜建模到结果后处理和踩坑记录全部倒出来。无论是刚接触FDTD的学生还是想快速上手超构透镜仿真的人照着这条路线走能省掉很多试错。1. 项目整体设计与思路拆解1.1 为什么偏偏是长波红外超构透镜长波红外LWIR8~14 μm是一个很特殊的波段。一方面它是大气透射窗口热成像、气体检测比如CO2激光的10.6 μm谱线、自由空间光通信都在这个区间另一方面传统红外光学元件又笨又贵——锗透镜、硒化锌透镜需要研磨抛光非球面加工更贵要做成多功能阵列难度很大。超构透镜的出现给了另一个思路用微纳柱状结构在一个平面基底上实现聚焦、偏折、偏振控制元件厚度只有几微米重量和体积直接缩小几个量级。但到了长波红外超构透镜的设计尺度跟可见光完全不是一回事。可见光波段单元尺寸是几百纳米而长波红外波长是10 μm上下单元周期可以做到3~5 μm这反而带来一个好处加工容差更宽松普通紫外光刻甚至激光直写可能就能做出来。仿真上也相对友好因为最小特征尺寸没到深亚波长网格压力没可见光那么恐怖。偏振复用聚焦说的是用一个超构透镜同时响应不同偏振态的光并让它们聚焦到不同位置。这意味着光学系统里可以省掉偏振分束器、波片这些额外的分立元件用一片平面结构就完成偏振分光与聚焦。这个能力在偏振成像、红外探测前端、激光加工多通道聚焦里都很有价值。1.2 偏振复用路线选型为什么要走几何相位偏振复用超构透镜的主流设计路线有两条一条是依赖共振相位与各向异性单元的“线偏振复用”另一条是基于几何相位也叫Pancharatnam–Berry相位PB相位的“圆偏振复用”。我这次选的是圆偏振几何相位路线原因有几个。几何相位的核心逻辑很简单当一束圆偏振光透过一个旋转了角度的各向异性单元时透射光中与原偏振相同的手性分量会获得一个大小等于两倍旋转角度的相位偏移。也就是说单元旋转角度θ → 相位变化2θ旋转0到180度相位就能完整覆盖0到2π。这个映射关系是纯几何的对波长不敏感只要单元的半波片条件满足一个旋转角度的扫描就能完成全相位调控设计效率很高。更关键的是左旋圆偏振和右旋圆偏振入射时几何相位会取相反符号2θ和-2θ。这意味着只要把透镜所有单元的旋转角度按某个相位分布设计LCP光得到一个波前RCP光自动获得互补的另一个波前。双焦点的偏振复用聚焦本质就是利用这种“一个角度两套相位”的性质让两种圆偏振分别聚焦到两个不同空间位置。相比线偏振复用需要在每个单元里独立设计两个正交偏振的响应几何相位方案的结构简单、参数少、扫描工作量小特别适合FDTD仿真验证。1.3 双焦点相位函数怎么定偏振复用的相位分布不是直接拿一个聚焦公式算的而是要把两个目标焦点对应的相位分布做融合。具体来说对于设计波长为λ、透镜中心坐标为(0, 0)、焦距为f目标焦点1位于(x1, 0, f)焦点2位于(x2, 0, f)理想相位分布为Φ₁(x, y) (2π/λ) · [√((x - x1)² y² f²) - f]Φ₂(x, y) (2π/λ) · [√((x - x2)² y² f²) - f]这里面减掉一个常数相f是为了让中心相位归零不影响波前形状。接下来关键一步我们实际能控制的是单元的旋转角度LCP光获得2θ的相位RCP光获得-2θ的相位。要让LCP光看到Φ₁、RCP光看到Φ₂理想情况是θ₁ Φ₁/2θ₂ -Φ₂/2但一个单元只有一个旋转角度不能同时满足两个条件。工程上常用的折中方案是把两个期望相位叠加成复数场再取相位也就是计算Φ_combined(x, y) arg[exp(iΦ₁) exp(iΦ₂)]然后令旋转角 θ Φ_combined / 2。这个方法在标量全息里叫双相位编码它会把能量近似均分到两个通道每个通道的理论衍射效率大约40%左右代价是总效率低于单焦点透镜。但如果两个焦点分得足够开串扰可以压到很低的水平。我在仿真实例里用的参数是λ10.6 μm透镜直径100 μm焦距120 μm两个焦点中心分别位于x30 μm和x-30 μm处。这套参数下双焦点横向分离60 μm远大于焦斑尺寸通道隔离度足够评估偏振复用效果。2. 材料选型与单单元结构设计2.1 长波红外材料怎么选超构透镜仿真很容易忽略的一个前置问题材料的折射率数据和损耗到底靠不靠谱。在长波红外材料选项跟可见光完全不同不能随手拿二氧化钛或氮化镓来用。我筛选材料时重点关注三个指标目标波段透明性消光系数要小、折射率对比度决定相位覆盖能力、工艺成熟度。候选材料里硅在1~7 μm表现很好但到了10 μm以上本征吸收明显上升做高透射超构透镜不合适锗在8~12 μm窗口非常透明折射率高达4.0左右是长波红外元件的传统材料硒化锌ZnSe折射率约2.4透明范围覆盖整个长波红外色散也小只是折射率偏低要达到同样的相位延迟需要更高的结构柱还有硫系玻璃折射率约2.6~3.2红外吸收低可以用纳米压印批量加工在长波红外超表面研究里也很流行。材料选择会直接影响单元的纵横比。折射率越高实现相同相位延迟所需的柱高越矮加工越容易。我优先选了锗原因很直接10.6 μm处折射率可参考约4.0柱高可以做到一个波长左右FDTD建模仿真时内存占用也相对可控。基底同样用锗这样柱子和基底之间没有异质材料界面的额外反射。这里要提醒一点Lumerical材料库里有Ge的Palik数据覆盖这个波段可以直接用比手动输入的近似值靠谱得多。材料10.6 μm处折射率长波红外吸收优点需要注意的问题Ge约4.0低折射率高工艺成熟高温下吸收增加ZnSe约2.4极低色散低稳定折射率低柱高较长Si约3.4偏高工艺最成熟长波端吸收损耗明显硫系玻璃约2.6–3.2低可压印成本低折射率可调但批次差异大2.2 单元结构建模圆柱不行必须打破旋转对称这是整个项目里最容易踩的第一个坑。很多人第一次做几何相位超构透镜想当然地在仿真里建立一个圆柱、旋转圆柱角度结果发现透过率和相位完全不变——因为圆柱是旋转对称的旋转它等于没旋转。几何相位要求单元必须有明确的本征光学轴方向最常见的就是椭圆柱、矩形柱、十字形或者带缺口的环形结构。我用的单元是锗椭圆柱。椭圆的长轴方向定义为结构旋转角度的参考方向。设计变量有四个椭圆长半轴a、短半轴b、柱高H、单元周期P。单元周期决定了工作波长下的衍射级次是否会被打开一般取亚波长量级我取P4 μm约为0.38λ保证零级透射为主。长短轴的比值直接决定单元的双折射大小。在FDTD里单独扫这个单元时要分别查看x偏振和y偏振入射下的透射幅度与相位。半波片条件要求两个正交偏振方向的透射相位差接近π同时两个偏振方向的透射率都尽量高。只有满足这个条件圆偏振入射时才能得到高强度且相位精确的相反圆偏振透射分量。换句话说这个单元本质上是一个“微小的波片”先扫出满足半波片条件的几何参数再谈旋转和相位调控。2.3 单单元FDTD扫描怎么把参数空间跑清楚单单元扫描建议用严格的周期边界Periodic跑一个单元光源用线偏振平面波垂直入射。要特别注意的是两次仿真中光源偏振方向需要分别与椭圆长短轴对齐。可以分别记录透射光的复振幅用透过率监视器加相位监视器得到振幅T_x、T_y和相位φ_x、φ_y。扫描策略要分层先固定一个初始高度比如H10 μm接近λ扫描长短轴组合观察相位差Δφ |φ_x - φ_y|然后在这个基础上扫描柱高因为高度其实是最快改变相位差的参数之一。以λ10.6 μm、P4 μm为例Ge椭圆柱的长轴扫描范围大致2.0~2.8 μm短轴1.0~1.8 μm高度8~12 μm。这种量级下FDTD网格最大步长设为0.2 μm左右即可一个单频仿真跑几秒钟就完成了完全不构成性能瓶颈。扫描结果里要同时看三个条件Δφ接近π、T_x和T_y都大于70%并且两者差距不要太大。因为几何相位转换效率跟这两个本征通道的透射振幅一致性密切相关若T_x和T_y相差过大即使相位差是π最后圆偏振通道也会出现不需要的零级透射分量导致偏振串扰。实际数据里我找到一组很稳的参数a2.6 μmb1.2 μmH9 μm。这组条件下两个正交偏振方向的透射率都在75%以上相位差在10.6 μm处约0.97π对于单波长设计已经足够用了。2.4 半波片条件没满足时会发生什么这里值得多说一句如果单元的双折射不够或者相位差不满足π圆偏振入射产生的几何相位分量幅度就会变小而且会混入一个“共极化”的背景透射分量。这个背景分量相当于一个不聚焦的简单透射基底打在人眼看到的模拟结果里就是在焦点旁边出现一个很大的中心亮斑偏振复用的对比度和效率都严重下降。所以单单元扫描阶段宁可多花时间也不能带着残次单元去做全透镜阵列仿真。我在初版仿真就是吃了这个亏短轴取太大相位差只有0.7π最后全透镜仿真得到的LCP焦点能量里混了一堆背景光一开始还以为是边界条件的问题排查半天才发现是单元本身没达到半波片条件。这个教训值得记下单单元是地基地基不牢后面所有优化都是白费。3. FDTD仿真建模与关键参数配置3.1 仿真工具与模型搭建流程我用的仿真工具是Ansys Lumerical FDTD老版本叫FDTD Solutions装的是2020 R1版本。长波红外超构透镜仿真的工作量不算小但好在结构简单不需要导入复杂CAD模型用内置的脚本函数就能完成建模这也是我推荐初学者走这条路的原因——你不需要去啃建模软件几行脚本就能生成几百上千个椭圆柱阵列。建模流程大致分四步。第一步是定义材料和基底基底用一块厚度3 μm的锗平板材料直接调用库里的Ge (Palik)。第二步是定义单个椭圆柱用addellipse命令设置中心位置、半径、旋转角度、高度和材料。第三步是写循环遍历透镜口径内所有网格点根据每个位置的相位值计算旋转角度并添加单元。第四步是设置光源、边界条件、监视器并运行仿真。整套流程看起来常规实际操作中一不留神就会在小地方出问题比如单位换算错误Lumerical默认单位是米但很多脚本习惯用微米这就够让人调一晚上的。3.2 圆偏振平面波光源很多人第一步就卡住FDTD仿真里没有直接叫“圆偏振平面波”的按钮这是新手最常卡住的地方。实现圆偏振入射的核心思路是用两个线偏振分量叠加振幅相等、相位差90度。在Lumerical中我采用的方法是添加一个平面波光源在光源属性里设置偏振相关参数将x偏振分量的振幅设为1y偏振分量的振幅设为1i或-1i这样在频域上就严格等效于左旋或右旋圆偏振入射。要提一下的是不同版本的Lumerical界面里这个设置路径不完全一样。2020 R1版本里可以在平面波光源的属性面板找到polarization angle或者类似参数同时在模型视图里看它的电场矢量方向更稳妥的做法是在脚本里直接设置光源的幅度与相位参数逐项指定Ex、Ey的实部和虚部。跑LCP和RCP两种情况时只需要把Ey的相位分别设为90度和-90度其他设置保持不变这样两个仿真之间的差异是纯粹可控的方便后续对比。还有一个细节光源应当放在结构的入射侧且距离结构表面至少半个波长保证结构近场散射不会反向耦合到光源平面造成反射仿真不真实。对于10.6 μm设计光源距结构表面取6~8 μm比较保险。3.3 边界条件与监视器的正确姿势单单元扫描时边界条件用Periodic周期边界因为我们要模拟的是无限周期阵列中的单个单元周期边界会自动在横向复制这个单元。注意斜入射才需要Bloch边界正入射情况下Periodic就够了。全透镜阵列仿真时没有周期性可言四周必须用PML完美匹配层吸收边界。PML层数我习惯设在8层以上并且PML到结构边缘留出至少一个工作波长的距离也就是10 μm以上。曾经我贪图节省计算区域把PML压到离结构只有3 μm结果反射杂波混进焦平面焦点附近多了好几个假次级峰。这个距离宁可多留后面的数据干净程度完全不同。监视器设置上我用了三类一是透射总功率监视器放在透镜基底之后用于归一化聚焦效率二是焦平面监视器放在zf处记录频域电场分布这是分析焦点位置和光斑大小的数据来源三是可以用一个x-z截面的场分布监视器直观看到两个偏振的聚焦波前形态。焦平面监视器的尺寸要覆盖两个焦点如果两个焦点间距60 μm监视器横向范围至少要120 μm宽否则另一个焦点被截掉一半都不知道。3.4 网格设置与收敛性验证FDTD的网格策略直接影响仿真可信度。网格太粗相位和透过率都不准网格太细内存和时间成本急剧上升。长波红外的优势在于波长长网格压力相对可见光小很多。我的经验是单单元扫描阶段最大网格步长设为0.2 μm约λ/53全透镜仿真阶段最大网格步长0.25 μm约λ/42。在锗内部折射率4.0等效波长2.65 μm这个最大网格仍然是等效波长的1/10以上精度可以接受。Lumerical的自动非均匀网格Auto non-uniform mesh建议开启同时最好打开共形网格conformal mesh对椭圆柱这类曲面结构共形网格能让边界的介电常量过渡更平滑透射谱的噪声会明显下降。网格收敛性验证必须做方法很朴素取同一结构把最大网格步长分别设为0.3、0.2、0.1 μm跑透射率和相位看结果是否稳定。如果0.2和0.1的结果相差很小说明0.2已经收敛。我实测下来从0.3降到0.2相位数据变化了约0.05 rad从0.2降到0.1变化不到0.02 rad所以选定0.2/0.25 μm这一档是安全的。如果完全不做收敛性验证审稿人或同事问一句“你的网格收敛了吗”就会很被动。4. 全透镜建模与偏振复用聚焦结果4.1 相位分布计算与柱体旋转角度映射全透镜建模前先把相位分布算出来。用前面提到的双相位编码方法算出每个单元位置(x, y)需要的旋转角度θ Φ_combined / 2。透镜口径设为半径50 μm圆心在(0,0)f120 μmx130x2-30。把这些数据算出来后注意一个细节由于几何相位对旋转角是模2π的所以θ需要通过取模操作映射到0到180度范围内这样每个椭圆柱的旋转角度不会超过半圈加工和建模都更友好。脚本建模时我用for循环遍历口径内的网格坐标点每个网格中心放置一个椭圆柱长轴a、短轴b、高度H取前面扫描得到的最优值旋转角θ按相位分布赋值。这里有一个提升效率的小技巧可以提前把旋转角度按位置存成一个二维数组然后用结构化遍历一次生成全部椭圆柱避免每加一个柱子都重新计算一次相位函数几百个单元时感觉不明显上到2500个单元时能省不少时间。需要注意单元位置与透镜口径的边界裁剪落在口径边缘内的单元旋转角正常赋值半落在口径外的单元直接跳过不要用“缩到口径边界内”的方式强行放置单元那会破坏局部相位分布。4.2 全尺寸仿真的性能压力全透镜阵列有大约2500个椭圆柱仿真区域约100 μm × 100 μm × 40 μm含上下PML。这个体量对FDTD来说不算小尤其要跑两个偏振态的入射综合时间和内存开销需要提前规划。我是在24核工作站上跑的开启并行加速后单次圆偏振入射仿真大概跑了4~6小时。如果你只有笔记本建议先做缩尺验证把透镜口径缩到60 μm直径仿一遍确认双焦点现象和光斑质量没问题再跑全口径。盲目追求“一步到位”容易跑了一晚上结果发现设错了监视器位置整个仿真白费。3D FDTD还有个技巧如果结构在x-y平面具有镜像对称性可以劈半仿真用对称/反对称边界条件把计算量减半。但这个透镜的相位分布没有对称性所以没法用。异质结构或参数异常导致的不收敛问题另说单就这个设计老老实实全尺寸仿真就是最稳妥的。4.3 焦平面光场分析双焦点、效率与光斑质量全透镜仿真结束后从焦平面监视器提取频域电场数据计算光强分布|E(x,y)|²。LCP入射时我在z120 μm处看到清晰的焦点峰值位置在x≈30 μmRCP入射时峰值位置跑到x≈-30 μm。两个焦点横向分离约60 μm与设计值吻合得非常好说明相位映射和旋转角度计算没有出错。光斑质量上我统计了焦点的半高全宽。x方向FWHM大约为15~16 μm与衍射极限基本一致。这个水平说明设计光路没有引入明显像差相位离散化和单元尺寸误差的影响在可接受范围内。聚焦效率是另一个关键指标。我用的是“焦点积分功率/透射总功率”的定义在焦平面以峰值位置为中心取一个半径10 μm的圆形积分区域计算落入这个区域的功率占总透射功率的比例。实测下来单个圆偏振通道的聚焦效率约35%左右。这个数值看起来不高但别忘了这是双相位编码带来的代价两个通道均分能量理论上限本来就只有单通道接近50%。跟文献里双焦点偏振复用超构透镜的典型效率相比35%是一个正常的工程水平。指标实测值说明LCP焦点位置(30.1, 0, 120) μm与设计值偏差0.1 μm量级RCP焦点位置(-29.8, 0, 120) μm与设计值偏差0.2 μm量级x方向FWHM15~16 μm接近衍射极限单通道聚焦效率约35%双相位编码的固有能量分配4.4 偏振串扰两个通道会不会互相“串门”偏振复用器最怕的就是LCP入射时在RCP焦点位置也看到亮点。我用两组仿真对比来评估串扰第一组LCP入射分别看两个目标焦点位置的能量第二组RCP入射同样看两个位置。结果发现两个通道的交叉位置能量大约只有主焦点峰值强度的3%~5%换算成功率串扰比低于-13 dB。这意味着两个偏振通道之间不会产生明显的“指纹重叠”。做到这个指标除了单元半波片条件达标之外还有一个因素是双焦点间距足够大——60 μm的间距大约是焦斑尺寸的四倍即使两个点光源稍微模糊也不会重叠。如果发现串扰偏高排查顺序建议是先回到单单元检查半波片条件再检查相位映射时是否做了正确的模运算最后才考虑缩小焦点间距导致的通道重叠。我在这个项目里没做进一步的串扰优化因为-13 dB对概念验证来说已经够用但如果目标是要做偏振成像这个指标还得通过调整相位组合算法或改用迭代优化来进一步压低。5. 常见问题与排查技巧实录5.1 问题速查表走到这一步把项目里遇到过的典型问题整理成一张表方便以后照单排查。现象可能原因处理方法旋转圆柱角度相位不变单元旋转对称没有光学轴换椭圆柱/矩形柱等各向异性结构透过率很低50%材料吸收数据不对或网格太粗检查材料库数据细化网格确认PML位置焦平面有中心背景亮斑单元半波片条件不满足返回单单元扫描优化高度与长短轴比双焦点间距与设计不符相位计算中λ单位或坐标单位错误代入数值复核用中心点相位做基准验证焦点旁边出现假次级峰PML离结构太近增大PML距结构的间隔至一个波长以上仿真内存爆掉网格过度细化或仿真区域过大合理放宽网格至收敛阈值先缩尺验证串扰过高10%半波片条件差或焦点间距太小优化单元参数增大焦点间距5.2 网格与内存的取舍心得FDTD仿真最常出现的争执就是“网格到底要多细”。我的原则是结构小特征尺寸决定网格下限但不需要无脑收紧。椭圆柱的最短轴1.2 μm这个特征尺寸相对网格0.2 μm有6个网格覆盖已经足够解析形状。再往上加密比如0.1 μm计算量翻好几倍结果差别却小到可以忽略。所以合理的流程是做一次定量收敛性测试而不是凭感觉设一个特别小的网格。内存方面3D FDTD的内存瓶颈主要来自仿真区域内的网格总数。我在仿真区域高度上做了一些压缩基底和柱子的上方空间留12 μm下方基底空间留8 μm加上PML总高度约40 μm。这个高度设置平衡了计算量和物理空间需求。如果后续需要把仿真波段展宽到8~14 μm整个窗口网格数会大幅增加那时建议先做单波长验证再考虑宽谱仿真。5.3 焦点位置偏移相位处理中的隐形雷我遇到过一次比较隐蔽的焦点偏移问题焦点位置与设计值偏差了近3 μm不是系统误差而是相位映射时没有做模运算导致的。具体来说旋转角θ取的是Φ/2但几何相位只有2θΦmod 2π如果直接把超过2π的Φ值除以2赋给θ就会造成柱体的旋转角度超过180度进而引入额外的相位歧义。解决办法是在算θ之前先把相位Φ映射到[0, 2π)区间然后再除以2。这个操作顺序不能反。另外由于几何相位对旋向的符号差异LCP和RCP通道看到的相位谱一个是Φ一个是-Φ在双相位编码的复数叠加中要特别注意符号否则两个焦点会跑到同一侧整个偏振复用没有实现。我建议在写脚本时就把两个通道的分量分开存先画出来看相位分布是否正确再生成模型省得全透镜跑一遍才发现问题。5.4 一个实用的校验技巧小口径先验证再放大口径这是最想分享的一条经验全透镜仿真前先跑一个口径减半比如直径50 μm的版本验证双焦点现象、焦点间距和大致效率。小口径版本的单元数只有625个仿真时间可以压缩到一个小时以内改参数和排查逻辑错误都非常快。等小口径结果确认无误再修整套参数跑全口径因为全口径版本运行时间长一旦发现错误返工成本太高。我在这个项目里之所以最后能一次跑出干净的双焦点结果靠的就是先做了口径60 μm的快速验证。当时小口径版本已经能看到两个清晰的焦斑间距也是预期的约40 μm小口径对应的缩放间距这时候我心里就有底了放大口径后大概率能成。做仿真最忌讳“直接上全尺寸”这种赌徒心态分阶段验证反而总耗时更短。写在最后这个项目从最开始的一无所知到最终跑出双焦点差不多花了两周时间大部分时间其实都耗在单元参数扫描和那些不起眼的小细节上。回头想最关键的几步认知转变就是几何相位必须用各向异性单元、单单元的半波片条件决定一切、圆偏振光源的本质是两个正交线偏振的叠加、全透镜建模前先做四分之一口径验证。把这四件事想清楚剩下的就是机械执行FDTD流程。下一步如果继续扩展我可能会把波长范围从单波长展宽到8~12 μm整段针对宽谱做色散优化也可以尝试在同一个透镜上实现更多偏振通道比如四个偏振态对应四个焦点相应的相位编码会复杂一些但思路是相通的。长波红外超构透镜的仿真路线其实没有想象中那么高不可攀只要把物理逻辑理清楚FDTD只是一个验证工具。工具永远是工具真正值钱的是你对机理的理解和对细节的把控。