1. 为什么长波红外超构透镜必须把偏振不敏感放在第一位热成像场景的现实约束1.1 你目标光源的偏振状态可能从一开始就把方案定死了长波红外LWIR8~12μm最典型的应用是热成像、气体检测、夜视辅助这类被动探测系统。这个波段里的目标物体比如人体、车辆发动机、飞行器蒙皮发出的都是热辐射本质上是大量随机取向原子的自发辐射叠加宏观上呈现的是天然的非偏振光。就算中间经过了金属反射面或者粗糙表面散射偏振度也会被完全打乱成为一个统计意义上各向同性的随机偏振场。这意味着一个红外物镜在设计时不能假设入射光是某个固定偏振。很多做可见光超表面的同行习惯用几何相位Pancharatnam–Berry相位因为单个柱子旋转角度就能连续覆盖0到2π相位设计自由度非常大。但几何相位有个绕不开的性质它对左右旋圆偏振的响应正好相反一个偏振获得φ相位另一个偏振获得-φ相位。入射光一旦变成线偏振可以分解成两个等权的圆偏振分量这两个分量经过透镜时一个被聚焦一个被发散最终焦面上只剩下一堆杂散背景焦点能量大打折扣。实际系统里如果要救回来就得在光路里加偏振片把一半能量砍掉还要忍受系统复杂度上升。对于讲究通过率、体积、成本的整机这基本不可接受。所以在项目开始敲方案之前第一个决策就是放弃几何相位改用传播相位方案并且把单元结构做成严格的旋转对称。只有这样的全介质超构透镜才能在没有任何偏振元器件辅助的情况下对任意偏振入射光给出一致的聚焦行为。这个结论在我做过的几个长波红外项目里反复验证过不是理论推导层面的区分而是实测焦斑和仿真吻合度上的本质区别。1.2 材料选型长波红外不是可见光的那套玩法长波红外波段可用的透明介质材料其实不多。可见光常用的二氧化钛、氮化镓、非晶硅在10μm附近要么吸收偏高要么折射率不够。高折射率低损耗的候选集中在硅和锗两个老牌红外材料上。硅在10μm的折射率约3.42色散很小带隙1.12eV对热辐射波段完全透明锗折射率约4.0相位调节能力更强但密度大、脆性高、热光系数也大环境温度一变化焦点漂移比硅明显。考虑到后续加工和封装我一般优先选硅除非设计上需要更大的相位延迟量才考虑锗。这里还要提醒一个特别容易踩的坑硅材料的损耗模型。FDTD软件材料库里自带的SiSilicon模型通常是从可见光到近红外波段的拟合数据在8~12μm区间未必可靠。更要命的是硅片的掺杂水平。高阻本征硅在长波红外是透明的但低阻硅比如重掺杂衬底片因为自由载流子吸收透过率会显著下降。建模前一定要确认实际流片用的硅片类型从Palik手册或者可靠文献里取该波段的n和k数据自定义材料导入。我见过有人直接用默认硅模型仿真效率算出来85%流片回来实测只有60%出头排查到最后发现就是材料吸收参数对不上。这类问题不会让FDTD报错但会实打实地毁掉你的实验计划。1.3 偏振不敏感单元的基本形态要求传播相位单元的原理是利用高折射率介质柱与周围环境之间的折射率差让透过柱体的光累积额外的相位延迟。这个相位延迟量取决于柱体的等效折射率而等效折射率又由柱径、柱高、周期共同决定。要让结构对x偏振和y偏振完全一致单元截面必须对称。最稳妥的是圆柱其次是正方形截面柱和正多边形柱。注意十字形或者带切角的形状虽然在某两个正交方向对称但在45度方向上响应可能不同FDTD里一验证就会露馅。工程经验是能选旋转对称的圆柱就选圆柱建模简单、网格剖分对称、加工时也容易保证侧壁均匀性。方形柱虽然制造上对光刻方向更友好但仿真和实测里经常因为轻微圆角或侧壁倾斜引入微小的偏振差异。我在长波红外项目里最终采用的单元就是硅圆柱阵列后面所有偏振验证都基于这个选择。2. 设计路径拆解从双曲相位分布到可查表单元库2.1 先把透镜的相位分布写清楚偏振不敏感聚焦透镜本质上还是一个大口径的衍射透镜。设计的第一步是确定工作波长、透镜直径、设计焦距这几个宏观参数。我这里以一个实际跑通的设计为例中心波长λ010μm透镜直径D400μm设计焦距f500μm数值孔径NAD/(2f)0.4。这个NA在长波红外超构透镜里属于中等偏上既能明显看到亚波长聚焦效果又不至于让边缘处相位梯度过陡导致单元库选不出合适的柱子。透镜上任意一个径向位置r处要让它发出的光与中心光线同相到达焦点需要补偿的相位是φ(r) - (2π/λ0) × [sqrt(f² r²) - f]这个双曲相位分布在r0处为0在边缘r200μm处相位绝对值约50rad折算过去接近8个2π周期。超构透镜不可能让每个像素都有连续相位它只能把连续的φ(r)对2π取模再映射到一组离散单元上。单元库的任务就是要提供一组透射率较高、相位覆盖0到2π的单元让离散化误差尽可能小。还有一个容易被忽略的点上面这个相位分布只和径向坐标r有关和方位角无关。也就是说透镜上的柱子在同一个半径圆环上相位设计应当完全一致。旋转对称的入射响应与旋转对称的相位分布结合在一起偏振不敏感才能成立。如果你设计的是偏折器、涡旋光产生器这类非旋转对称波前那偏振不敏感的实现条件完全是另一套逻辑不能套用这里的结论。2.2 传播相位为什么能在单个单元里覆盖0到2π介质柱可以理解成一段截断的波导。光穿过高度为H的柱子时相对周围空气获得一个额外相位Δφ ≈ (2π/λ) × (n_eff - n_air) × H其中n_eff是柱内传播模式的等效折射率它随柱径增大而增大但不会超过材料本体折射率。当柱径变化时n_eff在空气折射率约1和硅折射率3.42之间连续变化于是Δφ也连续变化。问题来了柱高H要多大才能把相位范围拉开到2π假设n_eff-n_air平均能到1.5~2在10μm波长下H如果只有2μm最大相位延迟只有约1.8~2.5rad远远不够。所以长波红外的介质柱必须做高典型设计在5~8μm量级。这个结论和可见光超透镜完全反过来可见光波段柱子高度只有几百纳米而长波红外一上来就是好几个微米的深柱。柱高越深深宽比越大加工侧壁的垂直度和光刻胶的形貌控制就越难所以这个参数不是想加就能加的要在相位覆盖范围和可制造性之间取平衡。2.3 单元库扫描的具体方案固定什么、扫描什么单元库扫描是做超构透镜最枯燥也是最重要的一步。我习惯固定周期p和柱高H扫描柱径d。以p5μm约λ0/2、H6μm为例d从1μm扫到4μm步长0.05~0.1μm能得到几十个单元数据点。每个点都用FDTD跑一次单元仿真记录透射率和透射相位。扫描范围要留两个安全边界d不能太接近p否则相邻柱子几乎连在一起等效折射率模型失效实际加工也会直接连成一片d也不能太小太细的柱子机械强度差在干法刻蚀过程中容易歪倒。一般把d/p控制在0.2~0.8区间内扫描既保证单元是个明确的柱体又留有余量。这一轮扫描跑完之后先别急着画相位曲线先检查x偏振和y偏振两组结果的一致性。如果单元结构真对称、网格剖分也对称相位差异应该在0.01rad量级。要是差异明显说明模型里存在非对称因素比如用多边形近似圆导致网格不对称或者CAD结构里混入了偏移。要在继续往下走之前修掉不然后面全透镜的偏振验证会一路偏离。3. FDTD单元仿真实操如何可靠地提取相位并筛选单元3.1 边界条件、光源和网格的关键设置我用Lumerical FDTD Solutions跑单元仿真MEEP或者其它FDTD工具的流程类似区别只在脚本接口。单元仿真模型是一个包含单根柱子的长方体盒子xy方向设置周期边界条件Periodic或Blochz方向设置PML吸收边界。光源用平面波从衬底一侧垂直入射。为了做偏振对比偏振方向分别设为x和y其他所有条件不变各跑一次。网格设置直接影响相位提取的准确性。10μm波长的场在直径只有几微米的介质柱内部变化非常剧烈默认的mesh accuracy 2往往不够。我一般手动在柱子区域添加mesh override最小步长压到0.05~0.1μm特别是在上下界面附近加密。这个成本在单元仿真阶段是可接受的因为一个单元模型很小算得再密也就几分钟。如果这一步偷懒网格太粗相位误差可能到0.1rad甚至更大每个单元误差一点点累积到几百上千个单元的透镜上焦点质量就会肉眼可见地恶化。监视器的放置也有讲究。透射监视器不能贴着柱子上表面那里倏逝场很强会把透射率算得虚高也不建议放太远白白增加计算域。我习惯放在柱顶上方1~2μm处记录频域电场再把透射功率归一化到入射平面波。这样得到的透射率和相位基本能反映远场行为。3.2 相位unwrap与偏振转换分量检查每次单元仿真结束后从频域监视器提取复电场E(f)透射相位是angle(E(f))透射率是|E(f)|²相对入射的比值。直接用angle函数拿到的相位在±π边界处会有锯齿跳变也就是相位wrap。如果直接把wrap之后的相位画出来看起来像单元库相位覆盖不足其实是后处理没做unwrap。Lumerical脚本里有unwrap函数但要注意它是按数组顺序做连续性判断的前提是d的扫描步长足够小让相邻单元的真实相位差不超过π。扫描步长如果太大unwrap也会失败这时要减小步长重新跑那一段。除了透射率和同偏振相位还要额外看一眼交叉偏振分量。对于理想对称单元x偏振入射不应该产生y偏振出射。如果仿真结果里交叉偏振的透射率明显非零说明结构里存在非对称性或者单元进入了某个强手性共振区。在偏振不敏感透镜的设计语境里交叉偏振分量就是损失和像差的来源必须记录下来。这一项的数值通常用偏振转化效率来表征即交叉偏振出射功率与入射功率之比。针对对称单元的偏振不敏感设计我们希望它越低越好如果你的目标恰好是做偏振转换器件那当然另当别论。3.3 单元筛选标准不止是覆盖0到2π单元库的质量直接决定全透镜性能上限。筛选时我按三条标准来透射率优先单个单元的透射率至少要大于50%最好稳定在70%以上。超构透镜的聚焦效率是多个单元透射率乘积级别的一个单元掉到30%整条光路的效率立刻被拖下去。相位覆盖连续所有可用单元的相位要能在0到2π区间内平滑连续分布不能出现大段空洞。如果某个柱径区间相位曲线斜率极陡说明那里存在强谐振加工尺寸稍微偏移相位就会剧烈漂移这样的单元宁可不用。相邻单元相位差平滑实际透镜里相邻柱子的直径是渐变过渡的。如果单元库里两个相邻柱径的相位跳变很大会在透镜局部引入明显的相位不连续形成杂散衍射。筛选时我会看一眼phase-d曲线局部斜率斜率过大的区间就标记为禁用区。筛选完成后把每个可用单元的直径、相位、透射率整理成一张映射表这就是后续全透镜建模的字典。整个过程里相位误差和透射率的取舍非常普遍。尤其在透镜边缘高NA区域有时候查表找不到透射率又高、相位又准的单元就得接受小范围的相位误差优先保证透射率因为透射率低造成的能量损失是全局性的而局部相位误差只会影响局部波前质量。4. 全透镜建模仿真从单元库到偏振验证和焦点评价4.1 用脚本批量生成全透镜几何模型全透镜包含几百到上千根柱子手动建模不现实必须写脚本。以Lumerical为例思路是遍历透镜孔径内的网格坐标计算每个坐标位置的目标相位φ(r)再根据单元库映射表查到对应的柱子直径然后用addcylinder这类命令逐根添加。由于透镜相位分布是旋转对称的只需要在半径方向按分辨率切出采样点然后按笛卡尔坐标填充到圆形孔径内即可。这里有一个建模细节如果把透镜划分成规则矩形网格圆孔边缘会有一圈被截断的柱子截面不再是完整的圆会引入额外的衍射。更稳妥的做法是只保留中心距离r小于D/2的网格点最外圈柱子按实际位置放置边缘处缺失的面积作为已知孔径误差接受下来。如果项目对边缘衍射特别敏感可以再对最外圈柱子的位置做径向微调让实际孔径尽量逼近设计圆。这一步对焦点主瓣影响不大但对旁瓣结构有影响夜间热成像这类对背景噪声敏感的应用值得多花这点功夫。4.2 偏振验证四件套x/y/左旋/右旋对比全透镜模型建好之后用覆盖整个孔径的平面波照射波长设为中心设计波长λ0。第一轮先用x偏振仿真记录焦平面上的横向强度分布和轴向光强曲线。然后改成y偏振再跑一次完全相同的监视器设置。理想的偏振不敏感透镜这两组结果的峰值强度、焦点位置、主瓣半高宽、第一旁瓣位置应当几乎重合。严谨起见再把圆偏振也补上左旋圆偏振一次、右旋圆偏振一次。对传播相位对称设计来说左右旋圆偏振的聚焦响应应当完全一致。实际仿真结果会因为数值噪声有微小差异但一般控制在1%以内。把x/y/L/R四组焦平面曲线叠在一起看如果四条曲线在焦点附近重合得很好偏振不敏感才算真正验证通过。这一步在单元仿真阶段没法替代因为全透镜的柱间耦合和边缘效应只有在整孔径仿真里才体现出来。我在项目里见过一个有意思的现象单元仿真x/y差异极小但全透镜仿真x/y焦点出现可分辨的偏移。排查下来是网格剖分的问题——全透镜模型网格比较大柱体边缘的网格剖分不再像单元仿真那样对称相当于数值上引入了各向异性。解决办法是对每根柱子单独添加mesh override或者把关键区域的网格步长压到单元仿真同等水平。内存压力确实大但偏振一致性校验不过关后面的效率数据都不可信。4.3 聚焦效率计算与波前质量判断聚焦效率是超构透镜最核心的量化指标定义是焦平面上以焦点为中心、一定半径范围内的能量与入射到透镜孔径总能量之比。积分半径不能随便取。我通常取衍射极限第一暗环半径r00.61λ/NA作为积分边界或者放宽到1.22λ/NA以包含主瓣的主要能量。取半径太小会把主瓣截掉一部分效率偏低半径太大则会把旁瓣和背景杂散光也算进来虚高。焦斑形状需要和理论艾里斑做对比。把焦平面横向强度归一化后观察第一暗环半径和旁瓣强度。如果旁瓣明显抬高说明相位离散化误差大、某个直径区间的透射率不均匀或者单元间耦合在起作用。更直观的做法是看焦斑的半高全宽FWHM衍射极限下FWHM约0.51λ/NA设计良好的超构透镜应当非常接近这一数值。如果FWHM偏宽超过10%大概率是相位分布离散化引入的波前误差回到单元库检查相位映射是否足够细。还有一种更系统化的评价方法从焦平面电场反推波前做Zernike分解或计算RMS波前误差。这个方法麻烦一些但能精确定位球差、彗差等低阶像差来源。超构透镜通常不引入大的离轴像差但高NA下边缘相位离散化会造成类似球差的波前偏差这时就能看见焦点沿轴向拉长、旁瓣不对称的现象。5. 长波红外仿真的坑与经验材料、网格、耦合、可制造性5.1 材料模型别直接用默认库长波红外的材料数据是个重灾区。软件自带材料库覆盖最全的是可见光和近红外到8~12μm往往只有少量拟合点甚至是外推值。我现在的习惯是任何长波红外项目启动前先花一天整理材料数据硅、锗、二氧化硅、氮化硅全部从文献里找该波段的实测n/k数据整理成自定义材料导入。这一天的投入能避免后面所有仿真结果的系统性偏差。另外注意温度。长波红外热成像系统往往在户外工作环境温度差异大。硅的热光系数约1.8×10⁻⁴/K左右锗约2.7×10⁻⁴/K。如果设计工作温度是20°C而目标是适应-20°C到60°C环境折射率变化带来的相位偏移会让焦点漂移几个波长量级。仿真阶段至少做一次高温和低温下的折射率敏感性检查评估焦点的轴向漂移是否在系统允许范围内。5.2 网格精度与内存的平衡策略单元仿真网格越细越好这个不用犹豫。但全透镜模型是另一个局面400μm直径、6μm高的柱子阵列如果把整个区域的网格都压到0.05μm内存轻松突破几十GB普通工作站根本跑不动。我的常见做法是区域细分柱体内部和上下界面用mesh override加密柱体之间的空气区域用默认较粗网格。如果模型还是太大利用结构对称性砍成1/4模型前提是光源和监视器也要满足对称条件。这里要提醒一句圆偏振光源不满足x/y对称条件用对称边界后没法做圆偏振验证。所以在我的流程里对称边界只用于前期快速估算正式偏振验证必须跑完整的3D模型。如果完整模型实在跑不动退一步的做法是先做2D切面仿真验证相位离散化方案在截面上的聚焦表现再用一个较小孔径比如直径100μm的3D模型验证偏振不敏感特性。小孔径透镜的主瓣更宽但偏振一致性和焦点位置仍然能说明问题。这个折中方案在长波红外波段特别实用因为大孔径全介质超构透镜的3D仿真资源需求确实惊人。5.3 单元间耦合超构透镜性能低于设计的头号元凶单元仿真里每根柱子都处于周期性无限阵列环境中相邻柱子直径相同。真实透镜里相邻柱子的直径是变化的近场耦合会让局部相位偏移偏离单元库标称值。这是所有超构透镜从设计到实测对不上的最常见原因在高NA区域尤其明显因为那里的相位梯度大、相邻柱径变化快。缓解耦合的方法一是增大周期p让柱子间距变大耦合减弱但p增大后要检查高阶衍射确保在工作数值孔径范围内不出现其他衍射级次的出射光。周期一般控制在λ/2附近这是一个兼顾耦合抑制和亚波长采样的经验值。方法二是做相邻单元组合验证把两个不同直径的柱子放进同一个周期边界仿真对比组合里的相位与单独仿真的相位差。这个检查不必全做只针对透镜中柱径变化最剧烈的位置采样验证。如果偏差持续超过0.1rad就需要调整单元库或透镜口径设计。5.4 从仿真到流片加工容差到底怎么预留长波红外全介质超构透镜的柱高动辄5~8μm深宽比接近甚至超过6:1这对工艺是实打实的挑战。干法刻蚀过程中侧壁倾斜、柱径缩窄、顶部圆角都是常见缺陷。设计阶段就该做容差扫描把柱径往上和往下各偏5%或0.1μm重新查相位和透射率的漂移。以柱径变化5%为例等效折射率可能变化0.05以上对应相位漂移可能超过0.5rad这种敏感单元就要标记出来避免使用。我在实践中发现硅圆柱在深宽比超过8:1后侧壁倾斜几乎无法避免。仿真时可以人为给柱体加一个1~2度的锥度看看偏振不敏感性能是否还能扛住。锥度会让x/y方向的等效截面产生微弱差异偏振不敏感性能出现退化。如果仿真显示锥度对性能影响不可接受就只能降低柱高、换更大周期的设计方案或者改用锗来以更小高度获取相同相位范围。最后说一个切身的体会长波红外的超构透镜设计仿真环节的成败往往不在算法而在于你有没有把整个流程里的每一个假设都当成可疑对象来验证。单元库阶段多花三天做网格收敛性检查、材料参数核对、x/y偏振一致性测试后面全透镜仿真就能省三周。我做这个项目时最值钱的一个习惯是把单元库的相位-直径表直接导出成CSV再用同一个脚本流程把相位分布反解成透镜版图数据一步到位。这样每次修改设计参数只需要改CSV和跑一遍脚本从单元库到加工文件的时间能压缩到半天以内。等你真正开始流片就会发现这个流程帮你省下的不仅是时间还有大把返工的成本。