企业官网建设流程全解析

长波红外超构透镜这个方向，我盯了有一阵子。真正上手之后才发现，最难的往往不是“偏振不敏感聚焦”这句口号本身，而是每一步设计到底怎么落地，以及你最后用FDTD仿真算出来的焦点，到底是物理真实还是“仿真的幻觉”。这篇文章就把我从设计目标、材料选型、单元扫描到全镜仿真的完整过程摊开来讲，重点放在偏振不敏感聚焦的实现逻辑和FDTD仿真设置里那些没人明说的坑，适合正在做超构透镜、或者刚接触光学FDTD仿真又不想在软件里瞎耗时间的朋友参考。

1. 项目做的到底是什么事

1.1 长波红外的应用背景

长波红外一般指8到14微米这个波段，正好落在大气窗口里面。常温附近的物体，热辐射峰值也集中在这个区间，所以安防监控、热成像、气体检测、工业测温基本都绕不开这个波段。传统长波红外光学系统用的是锗透镜、硒化锌透镜这类厚重的曲面元件，加工贵、体积大，像差校正也麻烦。

超构透镜的核心思路是用一层亚波长纳米结构阵列来替代传统透镜的曲面折射。每个纳米柱相当于一个微型相位调控单元，不同尺寸的柱子让透射光获得不同相位延迟，通过合理的空间排列就能把平面波汇聚成球面波。换句话说，用“雕刻”代替“打磨”，透镜可以做得非常薄，并且能在同一片基底上实现多焦点、消色差这些传统元件很难做到的功能。

做长波红外超构透镜还有个现实优势：波长长，意味着特征尺寸可以做得比较大。同样在10微米波长下，单元周期可以放到2到3微米，柱子的直径也有1到2微米，普通紫外光刻甚至激光直写就能搞定，不像可见光超构表面动不动就要电子束曝光加原子层沉积。这也是我把设计波长定在10微米附近的原因之一，工艺容差和仿真资源都更加友好。

1.2 偏振不敏感聚焦意味着什么

很多超构透镜设计用的是方形或者矩形纳米柱，这类结构对入射光的偏振方向是敏感的。圆偏振、线偏振方向一变，有效折射率就跟着变，焦点强度、效率都会波动。在真实系统里，入射光的偏振状态往往不可控，尤其是热成像场景，目标辐射本身是宽谱非偏振的，如果透镜对偏振敏感，成像效果就会不稳定。

所以这个项目的目标很明确：设计一种对入射偏振方向不敏感的超构透镜。也就是说，不管是横电波入射还是横磁波入射，焦点位置、聚焦效率、光斑形态都要基本一致。这个需求听起来简单，实际上直接决定了单元结构的选型，后面我会专门解释为什么圆柱结构是最省事的选择。

1.3 指标设定与设计路线

我这次按照一个中等难度的工程指标来设定参数：

设计路线分四步走：先做材料与单元结构选型，然后通过FDTD扫描单元库得到相位-半径映射表，接着按照聚焦相位分布离散化生成透镜版图，最后做全镜FDTD仿真验证聚焦效果，并分别跑TE和TM偏振来确认偏振不敏感。整条链路里，单元库扫描的准确性决定了最终透镜的成败，所以我会把大篇幅放在这块。

2. 结构选型与材料选择

2.1 为什么圆柱单元对偏振天然不敏感

先看一个简单的物理图像。电磁波照到一组亚波长柱状结构上，每个柱子相当于一段截断的介质波导，透射光的相位主要由波导模式的有效折射率决定。对同样的柱子，如果入射电场沿x方向还是沿y方向看到的“波导截面”不同，那有效折射率就不同，透射相位也就不同。

方形柱就是典型例子：边长一样的方柱，电场沿棱边和沿对角线入射，响应不一样，所以它本质上是双折射的。而圆柱对任何绕中心轴的旋转都不改变截面形状，TE波和TM波看到的是同一个几何，模式响应自然一致。我把这个类比成水管：方柱像带棱的方管，水流方向一变阻力就变；圆柱像圆水管，怎么转都是一样的。

所以偏振不敏感聚焦这个要求，基本锁定了单元结构必须具有旋转对称性。圆柱是最直接的方案，加工也最简单。后续如果非要提高自由度，也可以考虑正六边形、圆环或者同心圆环结构，但设计复杂度会指数级上升，对长波红外来说性价比很低。

2.2 长波红外的材料取舍

材料选择要同时看折射率、吸收系数和加工兼容性。长波红外候选材料主要是这几种：

我最后选了锗。原因很直接：折射率高，同样的高度能实现更大的相位跨度。超构透镜单元的相位延迟大概正比于有效折射率差乘以高度，折射率越高，柱子可以做得越矮，加工纵横比压力越小。锗在10微米附近的吸收系数很低，透过率表现好，而且半导体工艺里刻蚀锗也不算冷门。

这里特意提醒一句：别直接用普通单晶硅片。本征硅在可见光和近红外表现很好，但到长波红外会有多声子吸收，透过率明显下降。非要走硅工艺，得用高阻硅或者做浮区熔融处理，否则单元库扫描出来的透射效率会很难看。

2.3 相位调控的物理逻辑

单元相位调控的原理可以粗看成：对固定高度的圆柱，增大半径相当于提高了波导模式的等效折射率，透射相位就连续变化。严格来说这不是简单的平板波导公式，因为在亚波长周期内柱子之间还有耦合效应，但趋势是正确的：半径越大，相位延迟越大。

在一个周期单元里，入射平面波经过柱子后，除了透射的零级波，可能还有高阶衍射波。为了保证只有零级波参与聚焦，周期p必须足够小。判据通常取 p < λ / (1 + NA)，对10微米波长、NA=0.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微米的加密网格。

监视器布置是结果是否可靠的关键。我在三个位置放了频域场剖面监视器：

• 焦平面监视器：放在z=27微米处，对应柱顶面正向22微米的位置，记录Ex、Ey、Ez的完整场分布。
• z轴纵向监视器：沿光轴放一个二维场剖面，用来观察焦点位置和焦深。
• 透射功率监视器：放在柱顶面上方2微米处，记录通过透镜口径的总功率，作为聚焦效率的分母。

关于坐标系还要特别注意：如果你把计算区域z=0设在基底底面，基底5微米、柱高6.5微米，那么柱顶面在z=11.5微米，焦点在z=33.5微米。很多人把焦距的22微米直接当成仿真z坐标，结果焦点全跑到空气里判读错位了。建议在脚本里定义一个变量表示柱顶面位置，之后所有监视器坐标都由它推导，不要手写死数字。

4. 聚焦性能分析与偏振验证

4.1 焦点评价指标怎么定义

仿真跑完，原始结果是一堆场分布数据，直接说“有个明亮焦点”是没用的，得量化。我用三个指标：

聚焦效率的定义是焦点区域内的光功率除以通过透镜口径的入射功率。焦点区域半径取第一艾里斑半径，公式是1.22λ/(2NA)，对NA=0.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微米的小型验证镜，焦距按比例缩短，验证聚焦逻辑后再放大到最终口径。不要牺牲网格精度来换尺寸，因为网格精度决定了相位提取的准确性，相位错了，透镜口径再大也没有意义。

做完整套流程，我最深的体会是，偏振不敏感聚焦本身并不是难点，难点在于沿着整条设计链路保持清醒：单元扫描的相位误差会被透镜口径放大，全镜仿真和单元扫描的边界条件完全不同，网格收敛性不能偷懒。这套流程跑通之后，长波红外超构透镜从设计到验证基本就有了一套可以重复执行的范本，之后再去加消色差、多焦点这些高阶功能，也只是在现有框架上加砖添瓦。