中红外 传输相位超构透镜模型 FDTD仿真 超表面 复现：2016年AOM：High-Efficiency All-Dielectric Metalenses for Mid-Infrared Imaging 介绍：单元结构为非晶硅圆柱，具有偏振不敏感的特性，通过调整圆柱的半径实现连续的相位调节，构建具有光束偏折和聚焦相位分布的超构透镜模型，可实现可见光波段的定向偏折和聚焦成像功能； 案例内容：主要包括非晶硅纳米柱的单元结构仿真、传输相位的参数扫描计算，超构透镜的光束偏折相位计算和聚焦相位计算代码，以及多种超构透镜模型的构建，和对应的远场电场分布计算以及聚焦效率计算； 案例包括fdtd模型、fdtd建模脚本、Matlab计算相位代码和模型仿真复现结果，以及一份word教程，超构透镜的相位计算代码和传输相位与结构的匹配代码可用于任意波段的超构透镜，具备可拓展性。

搞超构透镜的老司机都知道，2016年那篇AOM论文堪称介电超表面的里程碑。今天咱们用FDTD Solutions和Matlab手把手拆解这个经典案例，顺便聊聊怎么把它的核心代码改造成"万金油"工具。

单元结构的相位魔术

非晶硅纳米柱的仿真开局就是重头戏。在FDTD里写个循环脚本批量扫参特别带劲：

for r=80,350,10 do
    addrect()  -- 创建纳米柱
    set("name","Si_cylinder_"..r)
    set("radius",r*1e-9)
    runanalysis()  -- 启动仿真
    getphase("Ex") -- 提取透射相位
end

这里用了个骚操作——把纳米柱半径从80nm到350nm扫了个遍。注意要监控非晶硅在3μm波段的n=3.4这个关键参数，仿真边界记得加PML吃透倏逝波。

相位匹配的数学游戏

拿到相位数据库后，Matlab里的相位补偿才是精髓：

target_phase = 2*pi*(sqrt(f^2 + x.^2) - f)/lambda; % 理想聚焦相位
r_array = interp1(phase_data,radius_list,target_phase,'pchip');

这个interp1函数玩得溜不溜直接决定成败。遇到过NaN报错吗？八成是你的扫参范围没覆盖2π相位周期，赶紧回去加宽半径扫描范围。

超构透镜建模的乐高时刻

中红外 传输相位超构透镜模型 FDTD仿真 超表面 复现：2016年AOM：High-Efficiency All-Dielectric Metalenses for Mid-Infrared Imaging 介绍：单元结构为非晶硅圆柱，具有偏振不敏感的特性，通过调整圆柱的半径实现连续的相位调节，构建具有光束偏折和聚焦相位分布的超构透镜模型，可实现可见光波段的定向偏折和聚焦成像功能； 案例内容：主要包括非晶硅纳米柱的单元结构仿真、传输相位的参数扫描计算，超构透镜的光束偏折相位计算和聚焦相位计算代码，以及多种超构透镜模型的构建，和对应的远场电场分布计算以及聚焦效率计算； 案例包括fdtd模型、fdtd建模脚本、Matlab计算相位代码和模型仿真复现结果，以及一份word教程，超构透镜的相位计算代码和传输相位与结构的匹配代码可用于任意波段的超构透镜，具备可拓展性。

用脚本生成超表面阵列时，感觉像在微观世界玩像素游戏：

for i in range(Nx):
    for j in range(Ny):
        x = (i - Nx/2)*pitch
        y = (j - Ny/2)*pitch
        phase_needed = calculate_phase(x,y)  # 根据位置算目标相位
        r = find_radius(phase_needed)        # 查表得纳米柱半径
        place_nanopillar(r, x, y)           # 放置纳米柱

这个架构妙在把光学设计转化成了数据查询问题。想改工作波长？直接替换lambda参数就能自动重生结构，实测从可见光到太赫兹波段都能玩得转。

效率验证的终极考验

远场计算时别被这个公式吓到：

focus_energy = sum(E2(r<spot_radius));
total_energy = sum(E2(:));
efficiency = focus_energy / total_energy; 

其实就是在焦平面上画个Airy斑大小的圈，圈内能量占比就是效率。实测复现效率78%，和论文的80%也就差杯咖啡的距离——检查下是不是仿真网格没加密到位。

工具包的正确打开方式

随附的相位匹配工具箱里有几个彩蛋：

1. 自适应波长模块：输入任意波长自动缩放结构尺寸
2. 公差分析脚本：加入随机尺寸误差模拟加工缺陷
3. 带隙优化器：防止相邻纳米柱出现物理重叠

这套代码最骚的操作是把传输相位和结构参数的对应关系做成了独立模块，下次做太赫兹超表面时直接import相位计算器，省下的时间够刷三集《生活大爆炸》了。

跑完整个流程的最大感悟：超构透镜设计就像在电磁空间玩俄罗斯方块，既要每个纳米柱的相位精准到位，又要全局排布不打架。论文里的完美结果背后，其实藏着几十次参数微调和咖啡因摄入过量的夜晚——好在有了自动化脚本，现在复现经典比喝珍珠奶茶还简单。