长波红外超构透镜这个方向我盯了有一阵子。真正上手之后才发现最难的往往不是“偏振不敏感聚焦”这句口号本身而是每一步设计到底怎么落地以及你最后用FDTD仿真算出来的焦点到底是物理真实还是“仿真的幻觉”。这篇文章就把我从设计目标、材料选型、单元扫描到全镜仿真的完整过程摊开来讲重点放在偏振不敏感聚焦的实现逻辑和FDTD仿真设置里那些没人明说的坑适合正在做超构透镜、或者刚接触光学FDTD仿真又不想在软件里瞎耗时间的朋友参考。1. 项目做的到底是什么事1.1 长波红外的应用背景长波红外一般指8到14微米这个波段正好落在大气窗口里面。常温附近的物体热辐射峰值也集中在这个区间所以安防监控、热成像、气体检测、工业测温基本都绕不开这个波段。传统长波红外光学系统用的是锗透镜、硒化锌透镜这类厚重的曲面元件加工贵、体积大像差校正也麻烦。超构透镜的核心思路是用一层亚波长纳米结构阵列来替代传统透镜的曲面折射。每个纳米柱相当于一个微型相位调控单元不同尺寸的柱子让透射光获得不同相位延迟通过合理的空间排列就能把平面波汇聚成球面波。换句话说用“雕刻”代替“打磨”透镜可以做得非常薄并且能在同一片基底上实现多焦点、消色差这些传统元件很难做到的功能。做长波红外超构透镜还有个现实优势波长长意味着特征尺寸可以做得比较大。同样在10微米波长下单元周期可以放到2到3微米柱子的直径也有1到2微米普通紫外光刻甚至激光直写就能搞定不像可见光超构表面动不动就要电子束曝光加原子层沉积。这也是我把设计波长定在10微米附近的原因之一工艺容差和仿真资源都更加友好。1.2 偏振不敏感聚焦意味着什么很多超构透镜设计用的是方形或者矩形纳米柱这类结构对入射光的偏振方向是敏感的。圆偏振、线偏振方向一变有效折射率就跟着变焦点强度、效率都会波动。在真实系统里入射光的偏振状态往往不可控尤其是热成像场景目标辐射本身是宽谱非偏振的如果透镜对偏振敏感成像效果就会不稳定。所以这个项目的目标很明确设计一种对入射偏振方向不敏感的超构透镜。也就是说不管是横电波入射还是横磁波入射焦点位置、聚焦效率、光斑形态都要基本一致。这个需求听起来简单实际上直接决定了单元结构的选型后面我会专门解释为什么圆柱结构是最省事的选择。1.3 指标设定与设计路线我这次按照一个中等难度的工程指标来设定参数参数数值说明中心波长10 μm长波红外大气窗口焦距22 μm从柱顶面算起透镜口径25 μm半径 12.5 μm数值孔径 NA约 0.5R/f 决定单元周期2.5 μm小于 λ/(1NA)目标聚焦效率大于 60%焦斑内能量占比TE/TM 效率差小于 5%偏振不敏感性量化指标设计路线分四步走先做材料与单元结构选型然后通过FDTD扫描单元库得到相位-半径映射表接着按照聚焦相位分布离散化生成透镜版图最后做全镜FDTD仿真验证聚焦效果并分别跑TE和TM偏振来确认偏振不敏感。整条链路里单元库扫描的准确性决定了最终透镜的成败所以我会把大篇幅放在这块。2. 结构选型与材料选择2.1 为什么圆柱单元对偏振天然不敏感先看一个简单的物理图像。电磁波照到一组亚波长柱状结构上每个柱子相当于一段截断的介质波导透射光的相位主要由波导模式的有效折射率决定。对同样的柱子如果入射电场沿x方向还是沿y方向看到的“波导截面”不同那有效折射率就不同透射相位也就不同。方形柱就是典型例子边长一样的方柱电场沿棱边和沿对角线入射响应不一样所以它本质上是双折射的。而圆柱对任何绕中心轴的旋转都不改变截面形状TE波和TM波看到的是同一个几何模式响应自然一致。我把这个类比成水管方柱像带棱的方管水流方向一变阻力就变圆柱像圆水管怎么转都是一样的。所以偏振不敏感聚焦这个要求基本锁定了单元结构必须具有旋转对称性。圆柱是最直接的方案加工也最简单。后续如果非要提高自由度也可以考虑正六边形、圆环或者同心圆环结构但设计复杂度会指数级上升对长波红外来说性价比很低。2.2 长波红外的材料取舍材料选择要同时看折射率、吸收系数和加工兼容性。长波红外候选材料主要是这几种材料折射率10 μm吸收损耗工艺难度备注锗 Ge约 4.0很低中等长波红外老牌材料容易获得高阻硅 Si约 3.45较低但有声子吸收高需要高阻率制造工艺成熟硒化锌 ZnSe约 2.4极低中等折射率偏低相位覆盖不足硫化锌 ZnS约 2.2极低中等折射率同样偏低我最后选了锗。原因很直接折射率高同样的高度能实现更大的相位跨度。超构透镜单元的相位延迟大概正比于有效折射率差乘以高度折射率越高柱子可以做得越矮加工纵横比压力越小。锗在10微米附近的吸收系数很低透过率表现好而且半导体工艺里刻蚀锗也不算冷门。这里特意提醒一句别直接用普通单晶硅片。本征硅在可见光和近红外表现很好但到长波红外会有多声子吸收透过率明显下降。非要走硅工艺得用高阻硅或者做浮区熔融处理否则单元库扫描出来的透射效率会很难看。2.3 相位调控的物理逻辑单元相位调控的原理可以粗看成对固定高度的圆柱增大半径相当于提高了波导模式的等效折射率透射相位就连续变化。严格来说这不是简单的平板波导公式因为在亚波长周期内柱子之间还有耦合效应但趋势是正确的半径越大相位延迟越大。在一个周期单元里入射平面波经过柱子后除了透射的零级波可能还有高阶衍射波。为了保证只有零级波参与聚焦周期p必须足够小。判据通常取 p λ / (1 NA)对10微米波长、NA0.5的情况算出来是6.7微米左右我用2.5微米已经是相当保守的选择高阶衍射基本可以忽略。相位覆盖范围则由柱高和有效折射率的变化范围共同决定。我最早把柱高设为5微米结果扫描完半径后相位跨度只有大约4.5弧度不够覆盖0到2π后来加高到6.5微米才留出足够裕量。这个细节我后面还会在常见问题里展开因为相位覆盖不足是最容易被忽略却又最致命的问题。3. FDTD仿真建模全流程3.1 单元库扫描边界条件、网格与参数设计单元库扫描是整个设计的基石。简单说就是在一系列固定半径下做FDTD仿真记录透射幅度和透射相位最后得到一张“半径-相位”映射表。这个过程我建议用FDTD类的时域工具完成我用的是商用软件开源方案比如MEEP也能做核心设置逻辑是一样的。单元结构按周期2.5微米建模基底是硒化锌厚度取5微米圆柱材料是锗高度先按6.5微米来。x和y方向必须用周期性边界条件模拟无限周期阵列z方向上下两个面用PML吸收边界。光源用平面波垂直入射偏振方向设为x。需要注意的是扫描阶段必须用周期边界不能为了省事用PML包住一个单元那模拟的是孤立柱子和阵列里的耦合行为完全不同。光源带宽可以取宽谱比如6到14微米这样一次扫描能同时看到单波长相位和多波长响应如果内存紧张也可以直接在10微米波长用单频连续波光源仿真时间要拉长到稳定状态。网格是扫描精度最敏感的地方。圆柱内部电场变化剧烈不能只依赖全局网格。我习惯的做法是全局网格步长设为0.15微米再给圆柱和周围一圈空气区域加一个网格加密区步长压到0.06微米。这里有个逐步收敛的检查要做把网格加密一倍重新提取相位变化小于百分之一才算收敛否则参数扫描表本身就是错的。半径扫描范围我取0.3到1.1微米步长0.05微米。为什么上限不是1.2微米因为周期2.5微米半径1.1微米时相邻柱子之间的净间距只剩0.3微米再大就会带来很强的耦合效应仿真出来的相位曲线会出现奇怪的跳变而且工艺上也做不出来。扫描结束后先看两件事相位跨度有没有超过2π透射效率有没有过低。如果某个半径下透射效率低于60%这个点要么不采要么在后续映射时权重降低。3.2 相位提取与聚焦相位离散化大部分FDTD工具都提供频域监视器可以直接输出复电场。但提取透射相位有个细节不能简单取某个平面监视器上所有网格点的平均相位因为场在横截面上有空间分布平均相位会被抹平。我用的方法是在周期性单元的中心正上方固定一个点监视器记录电场时域波形做傅里叶变换后取出10微米频率对应的复振幅和相位角。因为所有半径扫描用同一个参考点所以相位相对值是一致的这就可以了。拿到相位曲线后要做两步处理。首先是相位展开因为复相位会从π跳到-π直接用拟合会出问题。其次是用插值函数把半径映射到目标相位做成一个“查相位-选半径”的函数。聚焦相位分布本身用标准球面波公式φ(x, y) 2π/λ × (√(x² y² f²) − f)其中f是焦距这里取22微米。透镜中心相位为0越往边缘相位延迟越大所以离散化时给每个柱位算一个目标相位再去半径映射表里找最接近的半径。口径半径12.5微米周期2.5微米所以横向是5个格子整个透镜大约几十个柱子这个规模做全镜FDTD仿真完全可行。有一点要说清楚生成的柱阵不是规则方形网格那么简单。圆形口径内每个格点对应一个柱格点坐标就是柱心坐标。口径外的格点不加柱子。边界上的柱子尺寸也要按各自坐标算出不同的相位不要因为看着像阵列就全部用同一个半径。3.3 全镜仿真边界条件与光源监视器设置全镜仿真和单元扫描最大的区别在于x和y方向不能再用周期边界必须全部换成PML吸收边界否则透镜会被无限复制成周期阵列焦点旁边会出现一排虚假的衍射峰。这是我见过最多人踩的坑。全镜模型的光源还是平面波从基底一侧朝柱阵方向入射。因为透镜口径只有25微米仿真区域横向取40微米×40微米就足够了四周留出PML吸收带。z方向从基底底部延伸到焦点后方长度约35微米。网格不需要像单元扫描那么细全局0.15微米就够柱阵区域加0.08微米的加密网格。监视器布置是结果是否可靠的关键。我在三个位置放了频域场剖面监视器焦平面监视器放在z27微米处对应柱顶面正向22微米的位置记录Ex、Ey、Ez的完整场分布。z轴纵向监视器沿光轴放一个二维场剖面用来观察焦点位置和焦深。透射功率监视器放在柱顶面上方2微米处记录通过透镜口径的总功率作为聚焦效率的分母。关于坐标系还要特别注意如果你把计算区域z0设在基底底面基底5微米、柱高6.5微米那么柱顶面在z11.5微米焦点在z33.5微米。很多人把焦距的22微米直接当成仿真z坐标结果焦点全跑到空气里判读错位了。建议在脚本里定义一个变量表示柱顶面位置之后所有监视器坐标都由它推导不要手写死数字。4. 聚焦性能分析与偏振验证4.1 焦点评价指标怎么定义仿真跑完原始结果是一堆场分布数据直接说“有个明亮焦点”是没用的得量化。我用三个指标聚焦效率的定义是焦点区域内的光功率除以通过透镜口径的入射功率。焦点区域半径取第一艾里斑半径公式是1.22λ/(2NA)对NA0.5、λ10微米算出来约12.2微米。实现时在焦平面场分布数据里把中心半径12.2微米以内的坡印廷矢量z分量做面积分除以透射监视器上的总功率。半高全宽FWHM用来评价光斑大小直接在焦平面上沿径向取强度曲线找峰值一半对应的宽度。理想衍射极限下FWHM大约等于0.514λ/NA也就是10.3微米左右。斯特列尔比是峰值强度与无像差理想聚焦系统峰值强度的比。小于0.8一般就认为有明显像差。这三个指标各有侧重聚焦效率看能量利用率FWHM看光斑锐度斯特列尔比看相位误差整体水平。我建议全镜仿真后一次性全算出来这样后面做误差分析也有对比基础。4.2 TE/TM双偏振对比仿真偏振不敏感不能靠嘴说得分别跑TE和TM两种偏振入射。做法很简单在全镜仿真里把平面波光源的偏振方向从x改到y其他设置完全不变重新跑一遍。然后对比两者的聚焦效率、焦点位置、FWHM。我在我的设计参数下TE偏振聚焦效率约65%TM偏振约63%相差约2%FWHM都在10.8微米左右和衍射极限10.3微米很接近。这个结果说明圆柱结构的旋转对称性确实有效。但这里有一个细节容易被忽略FDTD数值网格本身可能引入偏振不对称。因为笛卡尔网格在x和y方向上的离散化程度一样理论上不会造成明显差异但如果网格加密区边界不圆、或者柱体被阶梯近似TE和TM仍可能出现百分之几的差异。要区分是物理还是数值可以把圆柱结构换成同一个网格下的等效旋转对称结构再跑一次不过这属于进阶诊断一般项目里看TE和TM差在5%以内就可以接受。4.3 带宽与焦深表现单波长超构透镜本质上是有色差的因为相位分布公式里波长是固定的。我用宽谱光源跑全镜仿真然后提取8、10、12微米三个波长的焦点分布结果10微米处效率最高8和12微米处效率会掉到50%左右焦点位置也会往光轴方向移动大约2到3微米。这是单波长设计的正常表现设计任务里并没有要求消色差所以我不追求带宽但会把这个数据记录下来作为后续迭代的基线。焦深方面沿光轴方向看强度分布半高全宽对应的焦深大约在正负8微米左右。对于焦距22微米的透镜这个焦深算正常范围。实际系统如果有装配误差这个焦深还算友好。5. 制造误差敏感性分析5.1 半径偏差与侧壁倾斜的影响仿真里用的是理想圆柱但实际加工出来的柱体一定会有偏差。我做了两组规则误差分析一组是半径整体膨胀20纳米另一组是把垂直侧壁改成80度的锥形侧壁然后把误差结构重新建模跑全镜仿真。半径膨胀20纳米的结果是聚焦效率下降约4%焦点位置基本不变。这主要是因为半径变化改变了相位而边缘部位相位梯度大、对误差更敏感。侧壁倾斜的影响更明显效率下降了大概10%FWHM也粗了一圈。所以如果工艺能力只能做到80度侧壁设计阶段就应该直接按锥形柱来建模而不是等加工出来再对着理想结构找原因。5.2 工艺窗口如何划定误差分析的主要产出是“工艺窗口”。比如半径误差±30纳米以内效率退化小于7%那说明这套设计对半径不敏感加工控制在这个范围内就行。如果侧壁角度变化5度就导致效率掉15%那工艺上就要盯紧刻蚀参数。这个信息对实际项目很重要因为仿真做得再漂亮工艺做不出来都是零。我在做这个项目时会把误差分析表附在设计交付文档里标注出“哪些误差可以放开、哪些必须卡死”。纳米柱的刻蚀纵横比是另一个约束我的柱高6.5微米、直径最小约0.6微米纵横比超过10深反应离子刻蚀需要专门调配方这点要提前和工艺同事对齐。6. 常见问题与排查经验6.1 相位扫描结果乱跳现象是半径连续变化相位曲线却像毛刺一样来回蹦。我先排查监视器位置点监视器离结构太近时容易受表面倏逝波干扰太远又会因为周期边界内高阶模式叠加出问题。比较好的位置是柱顶面上方0.5到1微米。其次要检查网格收敛性网格粗糙时相位对半径变化过于敏感加密后曲线会变平滑。还有可能是个别半径下透射效率太低噪声主导了相位这时候把扫描步长加密看看趋势有没有恢复。6.2 全镜仿真出现多余衍射峰这是我最想强调的问题。全镜仿真里如果还在x、y方向用周期边界FDTD会默认把整个透镜作为单个周期无限复制于是焦点两侧会出现等间距的虚假亮点看起来像是严重像差。排查方法很简单把边界条件改成PML后重跑虚假峰消失就说明边界设错了。这个错误太隐蔽因为仿真能跑完、能出图、还能看到焦点新手很容易忽略边界条件的物理含义。6.3 焦距位置判读混乱FDTD里不能只看场剖面监视器上的一个z位置就断言焦点在那。因为焦平面监视器如果放偏了看到的只是一个较大光斑误判为聚焦性能差。正确做法是先放一个沿z轴的线监视器找出z峰值位置再在那个位置放二维焦平面监视器。我的经验是聚焦仿真应该分两步走先纵向找峰值再横向看光斑这样最稳。6.4 网格收敛性怎么快速检验网格是否收敛是仿真结果可信度的根本。快速检验法是把全场网格和加密网格都加细一半重跑全镜仿真然后看聚焦效率变化。如果变化小于1%基本可信如果超过3%说明网格还没有收敛结果不能用来下结论。全镜仿真比较费时我一般是在单元扫描阶段就把网格收敛性验证好全镜阶段只做一次抽查。6.5 内存不足怎么办长波红外波长长网格相对宽松通常内存压力不大。但如果还是爆内存我的建议是把口径缩小比如做一个半径7.5微米的小型验证镜焦距按比例缩短验证聚焦逻辑后再放大到最终口径。不要牺牲网格精度来换尺寸因为网格精度决定了相位提取的准确性相位错了透镜口径再大也没有意义。现象常见原因解决方法相位曲线跳变监视器离结构太近点监视器上移0.5~1μm相位覆盖不足柱高太低增加柱高重新扫描整镜出现伪衍射峰复用周期边界全镜仿真改用PML焦点位置不对焦距换算错误先沿z轴找峰值再放焦平面监视器效率随网格波动大网格未收敛网格减半重跑扫描做完整套流程我最深的体会是偏振不敏感聚焦本身并不是难点难点在于沿着整条设计链路保持清醒单元扫描的相位误差会被透镜口径放大全镜仿真和单元扫描的边界条件完全不同网格收敛性不能偷懒。这套流程跑通之后长波红外超构透镜从设计到验证基本就有了一套可以重复执行的范本之后再去加消色差、多焦点这些高阶功能也只是在现有框架上加砖添瓦。
长波红外超构透镜偏振不敏感聚焦：单元扫描与FDTD仿真全流程详解
长波红外超构透镜这个方向我盯了有一阵子。真正上手之后才发现最难的往往不是“偏振不敏感聚焦”这句口号本身而是每一步设计到底怎么落地以及你最后用FDTD仿真算出来的焦点到底是物理真实还是“仿真的幻觉”。这篇文章就把我从设计目标、材料选型、单元扫描到全镜仿真的完整过程摊开来讲重点放在偏振不敏感聚焦的实现逻辑和FDTD仿真设置里那些没人明说的坑适合正在做超构透镜、或者刚接触光学FDTD仿真又不想在软件里瞎耗时间的朋友参考。1. 项目做的到底是什么事1.1 长波红外的应用背景长波红外一般指8到14微米这个波段正好落在大气窗口里面。常温附近的物体热辐射峰值也集中在这个区间所以安防监控、热成像、气体检测、工业测温基本都绕不开这个波段。传统长波红外光学系统用的是锗透镜、硒化锌透镜这类厚重的曲面元件加工贵、体积大像差校正也麻烦。超构透镜的核心思路是用一层亚波长纳米结构阵列来替代传统透镜的曲面折射。每个纳米柱相当于一个微型相位调控单元不同尺寸的柱子让透射光获得不同相位延迟通过合理的空间排列就能把平面波汇聚成球面波。换句话说用“雕刻”代替“打磨”透镜可以做得非常薄并且能在同一片基底上实现多焦点、消色差这些传统元件很难做到的功能。做长波红外超构透镜还有个现实优势波长长意味着特征尺寸可以做得比较大。同样在10微米波长下单元周期可以放到2到3微米柱子的直径也有1到2微米普通紫外光刻甚至激光直写就能搞定不像可见光超构表面动不动就要电子束曝光加原子层沉积。这也是我把设计波长定在10微米附近的原因之一工艺容差和仿真资源都更加友好。1.2 偏振不敏感聚焦意味着什么很多超构透镜设计用的是方形或者矩形纳米柱这类结构对入射光的偏振方向是敏感的。圆偏振、线偏振方向一变有效折射率就跟着变焦点强度、效率都会波动。在真实系统里入射光的偏振状态往往不可控尤其是热成像场景目标辐射本身是宽谱非偏振的如果透镜对偏振敏感成像效果就会不稳定。所以这个项目的目标很明确设计一种对入射偏振方向不敏感的超构透镜。也就是说不管是横电波入射还是横磁波入射焦点位置、聚焦效率、光斑形态都要基本一致。这个需求听起来简单实际上直接决定了单元结构的选型后面我会专门解释为什么圆柱结构是最省事的选择。1.3 指标设定与设计路线我这次按照一个中等难度的工程指标来设定参数参数数值说明中心波长10 μm长波红外大气窗口焦距22 μm从柱顶面算起透镜口径25 μm半径 12.5 μm数值孔径 NA约 0.5R/f 决定单元周期2.5 μm小于 λ/(1NA)目标聚焦效率大于 60%焦斑内能量占比TE/TM 效率差小于 5%偏振不敏感性量化指标设计路线分四步走先做材料与单元结构选型然后通过FDTD扫描单元库得到相位-半径映射表接着按照聚焦相位分布离散化生成透镜版图最后做全镜FDTD仿真验证聚焦效果并分别跑TE和TM偏振来确认偏振不敏感。整条链路里单元库扫描的准确性决定了最终透镜的成败所以我会把大篇幅放在这块。2. 结构选型与材料选择2.1 为什么圆柱单元对偏振天然不敏感先看一个简单的物理图像。电磁波照到一组亚波长柱状结构上每个柱子相当于一段截断的介质波导透射光的相位主要由波导模式的有效折射率决定。对同样的柱子如果入射电场沿x方向还是沿y方向看到的“波导截面”不同那有效折射率就不同透射相位也就不同。方形柱就是典型例子边长一样的方柱电场沿棱边和沿对角线入射响应不一样所以它本质上是双折射的。而圆柱对任何绕中心轴的旋转都不改变截面形状TE波和TM波看到的是同一个几何模式响应自然一致。我把这个类比成水管方柱像带棱的方管水流方向一变阻力就变圆柱像圆水管怎么转都是一样的。所以偏振不敏感聚焦这个要求基本锁定了单元结构必须具有旋转对称性。圆柱是最直接的方案加工也最简单。后续如果非要提高自由度也可以考虑正六边形、圆环或者同心圆环结构但设计复杂度会指数级上升对长波红外来说性价比很低。2.2 长波红外的材料取舍材料选择要同时看折射率、吸收系数和加工兼容性。长波红外候选材料主要是这几种材料折射率10 μm吸收损耗工艺难度备注锗 Ge约 4.0很低中等长波红外老牌材料容易获得高阻硅 Si约 3.45较低但有声子吸收高需要高阻率制造工艺成熟硒化锌 ZnSe约 2.4极低中等折射率偏低相位覆盖不足硫化锌 ZnS约 2.2极低中等折射率同样偏低我最后选了锗。原因很直接折射率高同样的高度能实现更大的相位跨度。超构透镜单元的相位延迟大概正比于有效折射率差乘以高度折射率越高柱子可以做得越矮加工纵横比压力越小。锗在10微米附近的吸收系数很低透过率表现好而且半导体工艺里刻蚀锗也不算冷门。这里特意提醒一句别直接用普通单晶硅片。本征硅在可见光和近红外表现很好但到长波红外会有多声子吸收透过率明显下降。非要走硅工艺得用高阻硅或者做浮区熔融处理否则单元库扫描出来的透射效率会很难看。2.3 相位调控的物理逻辑单元相位调控的原理可以粗看成对固定高度的圆柱增大半径相当于提高了波导模式的等效折射率透射相位就连续变化。严格来说这不是简单的平板波导公式因为在亚波长周期内柱子之间还有耦合效应但趋势是正确的半径越大相位延迟越大。在一个周期单元里入射平面波经过柱子后除了透射的零级波可能还有高阶衍射波。为了保证只有零级波参与聚焦周期p必须足够小。判据通常取 p λ / (1 NA)对10微米波长、NA0.5的情况算出来是6.7微米左右我用2.5微米已经是相当保守的选择高阶衍射基本可以忽略。相位覆盖范围则由柱高和有效折射率的变化范围共同决定。我最早把柱高设为5微米结果扫描完半径后相位跨度只有大约4.5弧度不够覆盖0到2π后来加高到6.5微米才留出足够裕量。这个细节我后面还会在常见问题里展开因为相位覆盖不足是最容易被忽略却又最致命的问题。3. FDTD仿真建模全流程3.1 单元库扫描边界条件、网格与参数设计单元库扫描是整个设计的基石。简单说就是在一系列固定半径下做FDTD仿真记录透射幅度和透射相位最后得到一张“半径-相位”映射表。这个过程我建议用FDTD类的时域工具完成我用的是商用软件开源方案比如MEEP也能做核心设置逻辑是一样的。单元结构按周期2.5微米建模基底是硒化锌厚度取5微米圆柱材料是锗高度先按6.5微米来。x和y方向必须用周期性边界条件模拟无限周期阵列z方向上下两个面用PML吸收边界。光源用平面波垂直入射偏振方向设为x。需要注意的是扫描阶段必须用周期边界不能为了省事用PML包住一个单元那模拟的是孤立柱子和阵列里的耦合行为完全不同。光源带宽可以取宽谱比如6到14微米这样一次扫描能同时看到单波长相位和多波长响应如果内存紧张也可以直接在10微米波长用单频连续波光源仿真时间要拉长到稳定状态。网格是扫描精度最敏感的地方。圆柱内部电场变化剧烈不能只依赖全局网格。我习惯的做法是全局网格步长设为0.15微米再给圆柱和周围一圈空气区域加一个网格加密区步长压到0.06微米。这里有个逐步收敛的检查要做把网格加密一倍重新提取相位变化小于百分之一才算收敛否则参数扫描表本身就是错的。半径扫描范围我取0.3到1.1微米步长0.05微米。为什么上限不是1.2微米因为周期2.5微米半径1.1微米时相邻柱子之间的净间距只剩0.3微米再大就会带来很强的耦合效应仿真出来的相位曲线会出现奇怪的跳变而且工艺上也做不出来。扫描结束后先看两件事相位跨度有没有超过2π透射效率有没有过低。如果某个半径下透射效率低于60%这个点要么不采要么在后续映射时权重降低。3.2 相位提取与聚焦相位离散化大部分FDTD工具都提供频域监视器可以直接输出复电场。但提取透射相位有个细节不能简单取某个平面监视器上所有网格点的平均相位因为场在横截面上有空间分布平均相位会被抹平。我用的方法是在周期性单元的中心正上方固定一个点监视器记录电场时域波形做傅里叶变换后取出10微米频率对应的复振幅和相位角。因为所有半径扫描用同一个参考点所以相位相对值是一致的这就可以了。拿到相位曲线后要做两步处理。首先是相位展开因为复相位会从π跳到-π直接用拟合会出问题。其次是用插值函数把半径映射到目标相位做成一个“查相位-选半径”的函数。聚焦相位分布本身用标准球面波公式φ(x, y) 2π/λ × (√(x² y² f²) − f)其中f是焦距这里取22微米。透镜中心相位为0越往边缘相位延迟越大所以离散化时给每个柱位算一个目标相位再去半径映射表里找最接近的半径。口径半径12.5微米周期2.5微米所以横向是5个格子整个透镜大约几十个柱子这个规模做全镜FDTD仿真完全可行。有一点要说清楚生成的柱阵不是规则方形网格那么简单。圆形口径内每个格点对应一个柱格点坐标就是柱心坐标。口径外的格点不加柱子。边界上的柱子尺寸也要按各自坐标算出不同的相位不要因为看着像阵列就全部用同一个半径。3.3 全镜仿真边界条件与光源监视器设置全镜仿真和单元扫描最大的区别在于x和y方向不能再用周期边界必须全部换成PML吸收边界否则透镜会被无限复制成周期阵列焦点旁边会出现一排虚假的衍射峰。这是我见过最多人踩的坑。全镜模型的光源还是平面波从基底一侧朝柱阵方向入射。因为透镜口径只有25微米仿真区域横向取40微米×40微米就足够了四周留出PML吸收带。z方向从基底底部延伸到焦点后方长度约35微米。网格不需要像单元扫描那么细全局0.15微米就够柱阵区域加0.08微米的加密网格。监视器布置是结果是否可靠的关键。我在三个位置放了频域场剖面监视器焦平面监视器放在z27微米处对应柱顶面正向22微米的位置记录Ex、Ey、Ez的完整场分布。z轴纵向监视器沿光轴放一个二维场剖面用来观察焦点位置和焦深。透射功率监视器放在柱顶面上方2微米处记录通过透镜口径的总功率作为聚焦效率的分母。关于坐标系还要特别注意如果你把计算区域z0设在基底底面基底5微米、柱高6.5微米那么柱顶面在z11.5微米焦点在z33.5微米。很多人把焦距的22微米直接当成仿真z坐标结果焦点全跑到空气里判读错位了。建议在脚本里定义一个变量表示柱顶面位置之后所有监视器坐标都由它推导不要手写死数字。4. 聚焦性能分析与偏振验证4.1 焦点评价指标怎么定义仿真跑完原始结果是一堆场分布数据直接说“有个明亮焦点”是没用的得量化。我用三个指标聚焦效率的定义是焦点区域内的光功率除以通过透镜口径的入射功率。焦点区域半径取第一艾里斑半径公式是1.22λ/(2NA)对NA0.5、λ10微米算出来约12.2微米。实现时在焦平面场分布数据里把中心半径12.2微米以内的坡印廷矢量z分量做面积分除以透射监视器上的总功率。半高全宽FWHM用来评价光斑大小直接在焦平面上沿径向取强度曲线找峰值一半对应的宽度。理想衍射极限下FWHM大约等于0.514λ/NA也就是10.3微米左右。斯特列尔比是峰值强度与无像差理想聚焦系统峰值强度的比。小于0.8一般就认为有明显像差。这三个指标各有侧重聚焦效率看能量利用率FWHM看光斑锐度斯特列尔比看相位误差整体水平。我建议全镜仿真后一次性全算出来这样后面做误差分析也有对比基础。4.2 TE/TM双偏振对比仿真偏振不敏感不能靠嘴说得分别跑TE和TM两种偏振入射。做法很简单在全镜仿真里把平面波光源的偏振方向从x改到y其他设置完全不变重新跑一遍。然后对比两者的聚焦效率、焦点位置、FWHM。我在我的设计参数下TE偏振聚焦效率约65%TM偏振约63%相差约2%FWHM都在10.8微米左右和衍射极限10.3微米很接近。这个结果说明圆柱结构的旋转对称性确实有效。但这里有一个细节容易被忽略FDTD数值网格本身可能引入偏振不对称。因为笛卡尔网格在x和y方向上的离散化程度一样理论上不会造成明显差异但如果网格加密区边界不圆、或者柱体被阶梯近似TE和TM仍可能出现百分之几的差异。要区分是物理还是数值可以把圆柱结构换成同一个网格下的等效旋转对称结构再跑一次不过这属于进阶诊断一般项目里看TE和TM差在5%以内就可以接受。4.3 带宽与焦深表现单波长超构透镜本质上是有色差的因为相位分布公式里波长是固定的。我用宽谱光源跑全镜仿真然后提取8、10、12微米三个波长的焦点分布结果10微米处效率最高8和12微米处效率会掉到50%左右焦点位置也会往光轴方向移动大约2到3微米。这是单波长设计的正常表现设计任务里并没有要求消色差所以我不追求带宽但会把这个数据记录下来作为后续迭代的基线。焦深方面沿光轴方向看强度分布半高全宽对应的焦深大约在正负8微米左右。对于焦距22微米的透镜这个焦深算正常范围。实际系统如果有装配误差这个焦深还算友好。5. 制造误差敏感性分析5.1 半径偏差与侧壁倾斜的影响仿真里用的是理想圆柱但实际加工出来的柱体一定会有偏差。我做了两组规则误差分析一组是半径整体膨胀20纳米另一组是把垂直侧壁改成80度的锥形侧壁然后把误差结构重新建模跑全镜仿真。半径膨胀20纳米的结果是聚焦效率下降约4%焦点位置基本不变。这主要是因为半径变化改变了相位而边缘部位相位梯度大、对误差更敏感。侧壁倾斜的影响更明显效率下降了大概10%FWHM也粗了一圈。所以如果工艺能力只能做到80度侧壁设计阶段就应该直接按锥形柱来建模而不是等加工出来再对着理想结构找原因。5.2 工艺窗口如何划定误差分析的主要产出是“工艺窗口”。比如半径误差±30纳米以内效率退化小于7%那说明这套设计对半径不敏感加工控制在这个范围内就行。如果侧壁角度变化5度就导致效率掉15%那工艺上就要盯紧刻蚀参数。这个信息对实际项目很重要因为仿真做得再漂亮工艺做不出来都是零。我在做这个项目时会把误差分析表附在设计交付文档里标注出“哪些误差可以放开、哪些必须卡死”。纳米柱的刻蚀纵横比是另一个约束我的柱高6.5微米、直径最小约0.6微米纵横比超过10深反应离子刻蚀需要专门调配方这点要提前和工艺同事对齐。6. 常见问题与排查经验6.1 相位扫描结果乱跳现象是半径连续变化相位曲线却像毛刺一样来回蹦。我先排查监视器位置点监视器离结构太近时容易受表面倏逝波干扰太远又会因为周期边界内高阶模式叠加出问题。比较好的位置是柱顶面上方0.5到1微米。其次要检查网格收敛性网格粗糙时相位对半径变化过于敏感加密后曲线会变平滑。还有可能是个别半径下透射效率太低噪声主导了相位这时候把扫描步长加密看看趋势有没有恢复。6.2 全镜仿真出现多余衍射峰这是我最想强调的问题。全镜仿真里如果还在x、y方向用周期边界FDTD会默认把整个透镜作为单个周期无限复制于是焦点两侧会出现等间距的虚假亮点看起来像是严重像差。排查方法很简单把边界条件改成PML后重跑虚假峰消失就说明边界设错了。这个错误太隐蔽因为仿真能跑完、能出图、还能看到焦点新手很容易忽略边界条件的物理含义。6.3 焦距位置判读混乱FDTD里不能只看场剖面监视器上的一个z位置就断言焦点在那。因为焦平面监视器如果放偏了看到的只是一个较大光斑误判为聚焦性能差。正确做法是先放一个沿z轴的线监视器找出z峰值位置再在那个位置放二维焦平面监视器。我的经验是聚焦仿真应该分两步走先纵向找峰值再横向看光斑这样最稳。6.4 网格收敛性怎么快速检验网格是否收敛是仿真结果可信度的根本。快速检验法是把全场网格和加密网格都加细一半重跑全镜仿真然后看聚焦效率变化。如果变化小于1%基本可信如果超过3%说明网格还没有收敛结果不能用来下结论。全镜仿真比较费时我一般是在单元扫描阶段就把网格收敛性验证好全镜阶段只做一次抽查。6.5 内存不足怎么办长波红外波长长网格相对宽松通常内存压力不大。但如果还是爆内存我的建议是把口径缩小比如做一个半径7.5微米的小型验证镜焦距按比例缩短验证聚焦逻辑后再放大到最终口径。不要牺牲网格精度来换尺寸因为网格精度决定了相位提取的准确性相位错了透镜口径再大也没有意义。现象常见原因解决方法相位曲线跳变监视器离结构太近点监视器上移0.5~1μm相位覆盖不足柱高太低增加柱高重新扫描整镜出现伪衍射峰复用周期边界全镜仿真改用PML焦点位置不对焦距换算错误先沿z轴找峰值再放焦平面监视器效率随网格波动大网格未收敛网格减半重跑扫描做完整套流程我最深的体会是偏振不敏感聚焦本身并不是难点难点在于沿着整条设计链路保持清醒单元扫描的相位误差会被透镜口径放大全镜仿真和单元扫描的边界条件完全不同网格收敛性不能偷懒。这套流程跑通之后长波红外超构透镜从设计到验证基本就有了一套可以重复执行的范本之后再去加消色差、多焦点这些高阶功能也只是在现有框架上加砖添瓦。
