简介针对光学超表面与超透镜领域的研究需求这份Matlab仿真资源围绕基于圆偏振复用实现双焦点超透镜的设计与验证展开重点解决在单一超透镜结构中编码并解码两个独立焦点的问题。相比传统方案圆偏振复用设计的集成度更高。资源包共1个文件为m脚本大小仅2KB代码紧凑涵盖理论模型构建、透镜几何参数设置、圆偏振态切换仿真及焦点分布可视化输出适合光学工程、微纳光学及信息光学方向的研究生和工程师快速复现并做参数化拓展。利用Matlab的光学仿真能力可在不加工透镜的情况下评估结构参数对双焦点位置和强度的影响透镜材料、几何尺寸与圆偏振态均可自由调整有助于优化设计通过运行脚本可直观对比不同圆偏振态下的焦点变化理解空间复用技术提升超透镜信息容量的机理。目前已有129人学习浏览可作为双焦点超透镜课题的起步模板为后续多焦点设计、实验验证或光学信息处理应用提供参考。1. 双焦点超透镜要解决什么问题为什么非要把两个焦点压进一层镜片先说个实际场景。做共聚焦显微或者双光镊实验的时候经常需要在两个不同深度同时聚焦传统做法要么加一条分光光路要么用多个透镜组累在一起要么上液晶空间光调制器动态切换。每一套方案都不便宜而且光学平台上一旦要同时架两套光路调试效率直线下降。双焦点超透镜的出现本质上就是想用一层平面元件替代这一堆东西入射圆偏振光的手性不同焦点位置就不同切换只需要动一下偏振片不涉及任何机械移动。这个标题里的圆偏振复用四个字才是整个方案的精髓。所谓复用就是在同一个器件上编码多套光学功能用光场的某个维度作为开关。常见的复用维度有波长、角度、偏振其中偏振复用最大的好处是不需要改变器件工作环境只需改变入射线偏振的方向或者左右旋圆偏振的手性就能瞬间切换功能响应速度只受偏振调控器件的限制。相比波长复用受制于光源波长范围偏振复用天然适合超表面这种对偏振敏感的结构。超透镜的本质是亚波长尺度的纳米天线阵列。每一个纳米柱也叫meta-atom都能独立调控入射光的振幅、相位和偏振所以我们可以把双焦点当作两套独立相位分布编码在同一片超表面上。在可见光波段想要做到这一点目前最成熟的手段就是几何相位也叫Pancharatnam-Berry相位。这个方案最吸引人的地方在于一套结构参数、一次刻蚀工艺完成两个功能整片透镜厚度只有几百纳米光路紧凑程度是传统系统无法比的。这篇文章我主要围绕圆偏振复用双焦点超透镜的仿真研究展开从PB相位原理讲到FDTD建模流程再到结果分析和排错经验。内容偏工程实操适合正在做超表面器件仿真、或者刚接手金属透镜课题的研究生也适合光学工程背景的同行想快速了解这一类器件的设计套路。2. 圆偏振复用的底层机制PB几何相位如何同时写进两套相位分布2.1 旋转纳米柱带来的几何相位先解释PB相位的基础。假设一个矩形纳米柱被放置在x-y平面内绕z轴旋转了角度σ。当一束圆偏振光垂直入射时透过光中与入射手性相反的交叉偏振分量会携带一个额外的相位大小恰好等于2σ或-2σ。具体符号取决于入射光的手性左旋圆偏振LCP入射透射的右旋圆偏振RCP分量相位增加 2σ右旋圆偏振RCP入射透射的左旋圆偏振LCP分量相位增加 -2σ。这个相位的来源不是光程差而是光在介质中传播时偏振态沿庞加莱球表面演化所积累的几何相位。它的好处非常直接只要把纳米柱旋转一个角度就能对交叉偏振分量精确赋予一个相位而且这个相位在很宽的波段范围内几乎无色散。正因为如此PB相位非常适用于需要同时实现多套相位分布的场景。2.2 双焦点编码的核心公式要让同一个超透镜对LCP聚焦在f1、对RCP聚焦在f2我们需要两套独立的相位分布[ \phi_{LCP}(r)-\frac{2\pi}{\lambda}\left(\sqrt{r^2f_1^2}-f_1\right), \quad \phi_{RCP}(r)-\frac{2\pi}{\lambda}\left(\sqrt{r^2f_2^2}-f_2\right) ]其中r是纳米柱中心到透镜中心的距离λ是设计波长。对于每个坐标位置纳米柱提供的总相位由两部分构成传播相位取决于纳米柱本身的尺寸形状和几何相位取决于旋转角。对LCP来说总相位是φ_prop 2σ对RCP来说总相位是φ_prop - 2σ。于是可以反推出每个位置需要的传播相位和旋转角[ \phi_{prop}(r)\frac{\phi_{LCP}(r)\phi_{RCP}(r)}{2}, \quad \sigma(r)\frac{\phi_{LCP}(r)-\phi_{RCP}(r)}{4} ]这个公式是整个设计流程的起点。实际操作中我习惯先写一个脚本对透镜面上每个像素坐标算出这两个数值存成两份矩阵后面查表反推结构参数用。下面是一段Python风格的伪代码逻辑和我在MATLAB里用的完全一致import numpy as np lam 633e-9 f1 30e-6 f2 45e-6 D 30e-6 # 透镜直径 N 512 # 网格数 x np.linspace(-D/2, D/2, N) X, Y np.meshgrid(x, x) R np.sqrt(X**2 Y**2) phi_lcp -2*np.pi/lam * (np.sqrt(R**2 f1**2) - f1) phi_rcp -2*np.pi/lam * (np.sqrt(R**2 f2**2) - f2) phi_prop (phi_lcp phi_rcp) / 2 sigma (phi_lcp - phi_rcp) / 4这里有个容易被忽略的问题σ的取值范围是无限的而几何相位2σ本身是周期性的纳米柱旋转角在0到π范围内编码的相位就已经覆盖0到2π。所以查表之前需要把σ归一化到0到π区间否则后面会出现同一个旋转角对应不同相位、查表不唯一的错误。我最初在这个问题上吃过亏后面排错章节会详细说。2.3 为什么选择圆偏振而不是线偏振可能有同行会问用线偏振复用不行吗理论上可以但实际限制很多。线偏振入射到各向异性纳米柱上透射的交叉偏振分量相位关系不仅跟旋转角有关还跟纳米柱在两个主轴方向的透射系数复振幅有关设计自由度虽然多了但对结构参数误差更敏感而且线偏振的方向角在工艺上不如圆偏振手性那么好精确切换。圆偏振复用把两个功能编码在同一个结构参数σ但符号相反的关系上结构简单相位稳定性更好这也是为什么目前文献里偏振复用超透镜大多数采用圆偏振方案。3. 从单元到整镜超透镜仿真建模的完整流程3.1 单元结构设计与参数扫描双焦点超透镜能不能成关键在于单元结构库的质量。我用的结构是标准的纳米鳍nanofin也就是矩形截面的TiO₂柱高度600nm周期350nm宽度和长度在80nm到240nm之间扫描。工作波长定在633nm衬底用熔融石英TiO₂折射率在632.8nm处约2.42。选TiO₂的理由很简单可见光波段损耗低、折射率高纳米柱高度不需要做得很高就能积累足够传播相位。如果是近红外波段换成非晶硅a-Si更常见折射率更高结构可以更矮。单元仿真的工具我推荐用RCWA或FDTD先做Bulk单元扫描RCWA速度快适合一次扫几千组参数。每组参数记录两个核心数据交叉偏振透射系数的相位和振幅。要注意的是扫描时必须明确传播相位的定义。很多文献用φ_prop arg(t_cross)来表示也就是交叉偏振透射系数的相位角。由于PB相位被分离出来后这部分不随旋转角变化所以φ_prop只由纳米柱的宽度、长度决定。扫描结果应满足两个筛选条件相位覆盖范围达到2π以上交叉偏振透射幅度尽量高我一般要求不低于0.8。如果相位覆盖不足2π说明柱子高度不够或者折射率太低需要增加高度如果透射率太低则要调整宽度-长度扫描范围。这一步没有捷径就是典型的参数网格扫描计算量也不大RCWA通常几十分钟能跑完。3.2 从相位库到透镜布局拿到单元相位库之后接下来就是反推设计根据第2节算出的φ_prop矩阵在相位库里查找匹配的宽度w, 长度l组合再根据σ矩阵给每个柱乘以旋转。查表时有两个细节一是要建立相位差最小的匹配策略。由于离散扫描步长的限制实际结构实现的传播相位和理想值之间总是有偏差一般控制在全库相位均方根误差小于0.1 rad就是比较好的水平。步长5nm的扫描密度通常能到这个精度如果不够就局部加密。二是要处理相位跳变。φ_prop的取值范围映射到0到2π后在边界处会出现从2π跳变到0的问题查表时如果不做unwrap操作会在透镜面上产生一圈明显的相位错位焦点质量直接崩掉。务必对φ_prop做一维unwrap之后再匹配。3.3 整镜仿真设置工具、边界与光源整镜仿真我推荐用有限时域差分法FDTD虽然慢但能直接把近场和远场传播算清楚结果更直观。Lumerical FDTD、Tidy3d、CST都可以我自己的习惯是单元结构用RCWA批量筛选整镜验证用Lumerical或者Tidy3dTidy3d的GPU加速在整镜大尺寸仿真上优势明显。整镜仿真中有几项参数直接决定结果的可靠性。光源用平面波波长633nm入射模式需要设置成圆偏振。早年间我直接在Lumerical里用两个正交线偏振平面波叠加、相位差90°来实现后来发现更稳的做法是用脚本定义注入偏振椭圆度并且必须在仿真中加一个偏振监测器验证实际入射光的椭圆率否则左旋入射变成了椭圆偏振两个焦点的响应会同时出现串扰指标完全失真。边界条件全部设为PMLx和y方向留出足够空气间隙z方向从衬底内部到空气传播域都要覆盖完整。监视器放在两个焦点平面位置以及经过透镜中心的x-z截面。计算内存方面整镜网格精度默认为2即网格尺寸约λ/10时内存占用已经不小如果透镜直径超过30μm普通16GB内存的机器就会吃力。我的做法是先用1/4象限对称性做仿真或者先用小口径透镜比如10μm验证整套设计逻辑确认无误后再放大到目标口径。4. 结果评估与优化迭代两个焦点之外还要看什么4.1 焦点位置的验证与指标量化仿真跑完第一步看x-z截面电场强度分布确认两个焦点分别出现在f1和f2附近并且两个焦点所在深度上没有明显的串扰焦点。但看到两个亮点只不过是最低标准真正要量化的是以下几个指标。第一个是焦斑尺寸。取焦点所在xy平面对光强分布做二维拟合计算焦斑的半高全宽FWHM。和理论衍射极限对比[ FWHM_{theory} \approx 0.515\frac{\lambda}{NA} ]对应前面参数直径30μm、f130μm时NA≈0.5衍射极限FWHM约652nmf245μm时NA≈0.33衍射极限约988nm。实测仿真FWHM如果和理论值偏差超过10%就需要检查相位离散误差或散焦像差是否太大。第二个是聚焦效率。定义上可以用焦平面焦点区域的功率积分除以入射到透镜口径的总功率。由于结构是偏振复用的效率还要区分信号效率和串扰功率。我统计时习惯只计算通过交叉偏振通道的能量因为共偏振通道不携带PB相位只能形成无用的背景光。单层介质超透镜的聚焦效率一般在30%~60%都是常见水平如果低于20%大概率是单元库透射率太低或者传播相位匹配误差太大。4.2 串扰是怎么评判的双焦点器件的串扰评判没有一个统一标准我的做法是对两种偏振分别做完整场仿真然后提取LCP入射时在f2平面的强度积分与f1平面的强度积分之比定义为串扰比。一个好的圆偏振复用器件串扰比应该低于-15dB。仿真中如果发现串扰偏大常见原因有三个入射圆偏振纯度不够这是仿真设置问题不是器件问题单元结构的几何相位幅度没有达到理想值交叉偏振转换效率低导致共偏振本底抬高双焦点所需的两套相位分布在某些位置自相容性差使得查表折中后每个焦点的波前误差都偏大。4.3 迭代优化的几个方向做完第一轮仿真我一般会根据结果做三方向优化。相位误差偏大时在单元库中把纳米柱宽度、长度扫描步长从10nm收细到5nm并重新匹配效率偏低时增加纳米柱高度或者换更高折射率材料如果特定焦点位置出现明显彗差检查相位公式——大NA情况下不能用paraxial近似代替精确双曲相位必须用完整的表达式。这些优化循环通常两三轮就能收敛。5. 仿真过程中的坑与排错清单5.1 网格收敛性测试FDTD仿真最容易被质疑的就是网格精度。不同软件默认网格策略差异很大我遇到过单元仿真切换软件后同一结构相位差了0.3rad的情况。科学的做法是做网格收敛测试在某一组代表性单元参数下把网格精度从默认值逐级加密记录相位和透射幅度随网格尺寸的变化直到两者波动小于1%取这个网格尺寸作为后续所有扫描的固定参数。这个测试虽然费时间但能避免整个单元库建在一个不准确的网格基准上那种错误后面几乎无法排查。5.2 相位参考平面不一致这是单元库和整镜仿真之间最容易出bug的地方。RCWA里默认的相位参考面可能在衬底表面、结构顶面或者某个自定义位置FDTD中监视器记录的相位同样依赖于参考平面。如果单元库和整镜仿真没有统一参考面整镜布局查表时所有传播相位会整体偏移一个常数。这个常数不会影响焦点位置相当于给整个透镜加了一个常数相位偏置但会严重影响焦斑幅值和焦平面上的零级背景。最稳妥的办法在单元仿真中就把参考平面定义为纳米柱顶面上方同一位置整镜验证时保持同样的监视器高度。5.3 旋转角映射到周期边界另外一个常见坑是旋转角和周期边界条件的矛盾。单元仿真时如果用了周期边界纳米柱旋转后理论上仍然满足周期性但当纳米柱宽度接近周期尺寸时旋转后结构可能超出原胞边界造成仿真结果突变。解决方法是把旋转角固定为某一组离散值通常每10°一个采样点每次旋转后独立重新构建几何模型而不是简单地平移坐标。很多新手直接把结构旋转后交给同一个网格文件去算结果相位谱出现毛刺这就是几何越界导致的。5.4 计算资源不够怎么办大口径超透镜的整镜FDTD仿真动辄需要几十GB内存。如果课题组没有高性能工作站我建议两条路。第一条是用Tidy3D这类云GPU求解器半小时就能跑完一个30μm口径的整镜成本不高适合学生党第二条是降维测试先用柱坐标旋转对称性把三维问题简化为二维或者只仿真一个扇形区域并配合对称/反对称边界条件虽然不能完全反映真实布局但能在设计定稿前快速筛查方案可行性。另外提一点经验之谈不要一上来就仿真大口径和高NA的完整器件。按单元库→小口径器件验证→放大口径→全口径仿真的顺序推进每一步都把预期结果先算好出了问题也能快速定位是单元库的问题还是布局算法的问题而不是在巨大的计算结果面前无从下手。结尾这套圆偏振复用双焦点超透镜的仿真流程我前前后后跑了大半年最大的感受是设计原理本身并不复杂难的是让仿真结果对得上设计预期。单元库质量、相位匹配算法、网格收敛性、偏振纯度这些环节任何一个有一点偏差最终的焦点质量都会明显恶化。特别是旋转角映射和相位参考平面这两个细节几乎每次重做都要提醒自己再检查一遍。如果你也在啃这个方向建议先从小口径、低NA起步把整个流程跑通再逐步放大指标这样能少走很多弯路。希望这篇记录对你的仿真研究有实际帮助。本文还有配套的精品资源点击获取
