多波长独立聚焦超构透镜技术研究：FDTD仿真超表面模型案例解析

多波长 独立聚焦超构透镜 fdtd仿真 超表面 复现论文2017年OEDispersion controlling meta-lens at visible frequency 论文介绍单元结构为硅矩形纳米柱结构通过调节结构的长宽尺寸可以找到三个波长处高偏振转换效率的参数通过调整纳米柱的转角实现连续的几何相位调节构建具有三个独立波长聚焦相位分布的超构透镜模型可实现可见光波段的三原色聚焦和成像 案例内容主要包括硅纳米柱的单元结构仿真、偏振转换效率的计算几何相位的计算超构透镜的不同色散曲线对应的超构透镜相位计算matlab代码不同色散的超构透镜模型以及对应的远场电场分布计算 案例包括fdtd模型、fdtd建模脚本、Matlab计算相位代码和模型仿真复现结果以及一份word教程超构透镜的不同色散相位计算代码可用于任意波段的超构透镜具备可拓展性一、代码项目概述与文件体系自旋-轨道角动量SAM-OAM转换是超构表面光学领域的核心技术之一其通过微纳结构对入射光偏振态与相位的精准调控可将携带自旋角动量的圆偏振光转换为具有螺旋相位结构的轨道角动量光束在光通信、量子信息处理、高分辨成像等领域具有重要应用价值。本代码项目围绕该技术展开通过FDTD时域有限差分法仿真与Matlab理论计算构建了覆盖“单元结构验证-超表面建模-光场转换仿真-理论结果对照”的全流程代码体系所有功能实现均基于《自旋-轨道角动量转换超构表面模型案例.pdf》以下简称“参考文档”中的技术参数与设计逻辑文件与功能的对应关系清晰明确为SAM-OAM转换超构表面的设计与验证提供完整技术支撑。多波长 独立聚焦超构透镜 fdtd仿真 超表面 复现论文2017年OEDispersion controlling meta-lens at visible frequency 论文介绍单元结构为硅矩形纳米柱结构通过调节结构的长宽尺寸可以找到三个波长处高偏振转换效率的参数通过调整纳米柱的转角实现连续的几何相位调节构建具有三个独立波长聚焦相位分布的超构透镜模型可实现可见光波段的三原色聚焦和成像 案例内容主要包括硅纳米柱的单元结构仿真、偏振转换效率的计算几何相位的计算超构透镜的不同色散曲线对应的超构透镜相位计算matlab代码不同色散的超构透镜模型以及对应的远场电场分布计算 案例包括fdtd模型、fdtd建模脚本、Matlab计算相位代码和模型仿真复现结果以及一份word教程超构透镜的不同色散相位计算代码可用于任意波段的超构透镜具备可拓展性从文件组织来看项目核心代码与数据文件分散于三个关键文件夹各文件夹的功能定位、包含文件及核心产出如下表所示文件夹名称功能定位核心文件/脚本关键产出物Nanofins单元纳米柱超构表面基本功能单元光学特性仿真与参数优化unitTiO2模型文件、unitsweep.lsf脚本纳米柱正交偏振态相位变化曲线、透射率曲线、不同方向电场分布场图Exz监视器数据FDTD models超构表面整体结构建模、SAM-OAM转换光场仿真、远场特性分析OAMpolar.lsf、OAMpolar_2.lsf、farfield.lsf脚本6个超表面仿真模型不同拓扑荷数OAM光束的远场强度分布图、相位分布图空间复用超表面的双向转换结果Matlab Models涡旋光场理论计算、干涉场模拟平面波/球面波/MZ干涉仪涡旋光场计算脚本、干涉场模拟脚本涡旋光束光强分布三维图、与不同光束的干涉条纹图、相位分布可视化结果二、Nanofins文件夹代码单元纳米柱仿真功能详解单元纳米柱是超构表面实现SAM-OAM转换的“最小功能单元”其几何参数尺寸、旋转角度与光学特性相位调控能力、透射效率直接决定超构表面的整体性能。Nanofins文件夹中的代码聚焦于单元纳米柱的仿真与参数验证通过系统性扫参与结果分析为后续超表面整体建模提供最优单元参数。2.1 核心代码功能从结构定义到结果输出2.1.1 单元结构参数初始化代码代码首先完成TiO₂纳米柱的基础参数定义严格遵循参考文档中“基于文献的标准参数设置”要求。代码会自动导入TiO₂材料的光学参数折射率、消光系数等适配工作波长并定义纳米柱的关键几何参数周期参考文档中未明确数值但通过仿真逻辑可判断为亚波长尺度以避免衍射干扰、高度确保光场在纳米柱内充分作用满足相位调控需求、横截面尺寸宽与高的比例需平衡透射率与相位调控精度。同时代码会设置衬底材料通常为SiO₂的参数明确衬底与纳米柱的空间位置关系构建“衬底-纳米柱”的基础单元结构模型。该部分代码的核心功能是“参数标准化”——通过固定材料与几何参数的初始值确保后续扫参仿真的变量唯一性仅聚焦“旋转角度”这一核心变量避免多变量干扰导致的结果偏差。2.1.2 旋转角度扫参代码sweep功能实现参考文档指出超构表面的相位调控为“几何相位”且相位与单元旋转角度呈2倍关系这是SAM-OAM转换的核心原理。基于此代码通过sweep扫参功能实现对纳米柱旋转角度的系统性遍历扫参范围通常设置为0°-180°因旋转角度超过180°后光学特性会出现重复无需冗余计算扫参步长根据精度需求设置如1°或0.5°步长越小相位曲线越平滑但计算耗时增加监测指标代码会同步监测“正交偏振态的相位变化”与“透射率”两个关键指标——前者验证几何相位的2倍关系后者确保单元在相位调控过程中具有较高的光传输效率避免因结构设计导致的光能量损耗。在扫参过程中代码会自动记录每个旋转角度对应的相位值与透射率数据并存储为中间数据文件为后续结果分析提供原始数据支撑。2.1.3 结果分析与可视化代码unit_sweep.lsf脚本unit_sweep.lsf是Nanofins文件夹的核心分析脚本其功能是将扫参得到的原始数据转化为直观的可视化结果与量化结论具体实现以下功能数据读取与预处理代码读取扫参生成的中间数据文件提取“旋转角度-相位”“旋转角度-透射率”两组数据并对异常值如因计算误差导致的跳变值进行平滑处理曲线绘制通过绘图函数生成两条核心曲线——“透射率-旋转角度曲线”用于判断不同角度下单元的光传输效率筛选透射率≥80%的角度区间、“归一化相位-旋转角度曲线”验证相位变化量是否严格等于2×旋转角度若曲线斜率接近2则说明单元满足几何相位调控要求电场分布场图输出代码调用Exz监视器数据生成不同旋转角度下纳米柱周围的电场分布场图直观展示光场在单元内的传播与偏振转换过程——例如当旋转角度为90°时电场方向会发生明显偏转且强度分布均匀说明单元的偏振转换效果符合预期。2.2 代码功能价值为超表面建模奠定基础Nanofins文件夹的代码本质是“超构表面单元的性能验证工具”——通过仿真确定最优的单元旋转角度范围、几何尺寸参数确保每个单元都具备精准的相位调控能力与高透射率。若单元相位调控精度不足如相位与旋转角度偏离2倍关系后续超表面整体仿真会出现OAM光束拓扑荷数混乱、远场强度分布异常等问题因此该部分代码是整个项目的“基础保障”。三、FDTD models文件夹代码超构表面仿真与远场分析功能FDTD models文件夹是整个项目的核心其代码实现了超构表面的整体建模、SAM-OAM转换光场仿真与远场特性分析直接产出SAM-OAM转换的核心结果。参考文档明确该文件夹包含6个模型分为“单拓扑荷数转换”与“空间复用双转换”两类对应不同的代码功能逻辑。3.1 单拓扑荷数SAM-OAM转换代码OAM_polar.lsf脚本该脚本的功能是构建单一拓扑荷数的超表面模型实现“入射右旋圆偏振光→出射固定拓扑荷数左旋圆偏振光”的转换核心功能分为“模型构建”“光场设置”“仿真运行”三个环节。3.1.1 超表面模型构建代码代码首先定义超表面的整体参数包括超表面半径覆盖OAM光束的有效作用区域、单元排列周期与Nanofins文件夹中确定的单元周期一致确保单元间无耦合干扰、工作波长参考文档中未明确但通过仿真逻辑可判断为可见光或近红外波段需与材料参数匹配。随后代码通过以下步骤构建超表面阵列单元数量计算根据超表面半径与单元周期计算径向单元层数如半径13μm、周期0.325μm时层数约40层与每层单元数量采用正六边形排列每层单元数量6×层数确保结构对称性单元位置与旋转角度赋值通过双重循环遍历每层i与每个单元j计算单元的圆心坐标极坐标转直角坐标xr×cosφyr×sinφ其中r为单元到中心的距离φ为单元极角同时根据“拓扑荷数与旋转角度的关系”参考文档指出l2ql为OAM拓扑荷数q为调控参数为每个单元分配旋转角度alphaq×(φ/360×2π)初始相位确保单元旋转角度呈螺旋分布形成螺旋相位面模型生成代码自动调用Nanofins文件夹中验证通过的单元模型按照计算的位置与旋转角度生成完整的超表面阵列模型保存为FDTD可识别的模型文件。3.1.2 入射光场与监测设置代码为实现SAM-OAM转换代码需精准设置入射光场参数与监测点具体功能如下入射光源设置代码生成右旋圆偏振光入射光源——通过添加两个正交的高斯光源x方向与y方向设置偏振角分别为0°与90°相位差90°合成右旋圆偏振光同时定义光源的束腰半径确保光束完全覆盖超表面与位置位于衬底下方模拟实际实验中的光入射路径监测器设置添加2D Z-normal功率监视器放置于超表面上方如z1μm处用于记录出射光场的强度与相位数据——这些数据是后续远场分析的核心原始数据监测器的采样精度如空间采样点数量会直接影响远场结果的分辨率。3.1.3 仿真运行与数据存储代码代码设置仿真时间步长根据工作波长与FDTD稳定性条件确定确保数值计算收敛与总仿真时长需满足光场完全通过超表面并传播至监测器随后启动仿真。仿真过程中代码会实时记录监测器数据并在仿真结束后将数据存储为特定格式如.lsf数据文件供后续远场分析脚本调用。3.2 空间复用双SAM-OAM转换代码OAM_polar_2.lsf脚本参考文档指出该脚本对应的模型可实现“空间复用双轨道角动量转换”即通过同一超表面对不同自旋方向的入射光实现差异化转换。其代码功能与OAM_polar.lsf的核心差异在于“单元旋转角度的分配逻辑”具体如下双调控参数设置代码引入两个调控参数q1与q2参考文档中未明确数值但通过结果描述可判断q1与q2对应不同的拓扑荷数分别用于奇数层与偶数层单元层数判断与角度赋值通过mod(i,2)函数判断单元所在层数的奇偶性——奇数层单元采用q1计算旋转角度alphaq1×(φ/360×2π)初始相位偶数层单元采用q2计算旋转角度或对q1计算的角度取反实现相位反转双向转换验证代码在仿真设置中会分别输入右旋圆偏振光与左旋圆偏振光记录两种入射光对应的出射光场数据——参考文档指出该模型可实现“入射右旋圆偏振光→出射l5左旋圆偏振光”“入射左旋圆偏振光→出射l-5右旋圆偏振光”代码会通过数据对比验证这一结果。3.3 远场分析代码farfield.lsf脚本超表面的SAM-OAM转换效果最终需通过远场特性体现如OAM光束的螺旋相位、“甜甜圈”型强度分布farfield.lsf脚本的功能是将监测器记录的近场数据转换为远场分布并实现可视化分析具体功能如下近场数据读取代码调用FDTD仿真生成的监测器数据文件提取出射光场的电场分量Ex、Ey、Ez远场计算定义远场计算网格x/y方向覆盖OAM光束的传播范围z方向设置为远场距离需满足远场条件fR²/λR为超表面半径λ为工作波长通过farfieldexact3d等函数计算远场光场分布包括x-z平面与x-y平面的光场数据偏振分离与结果输出代码分离出射光的左旋圆偏振光LCPEx - i×Ey与右旋圆偏振光RCPEx i×Ey分量分别计算并绘制两种偏振光的远场强度图与相位图——强度图需呈现“甜甜圈”型中心暗斑周围亮环相位图需呈现连续的螺旋条纹条纹数量等于拓扑荷数l代码会通过图像分析验证这些特征是否符合设计预期如l5时相位图应出现5条螺旋条纹。四、Matlab Models文件夹代码涡旋光场理论计算与干涉模拟功能FDTD仿真结果的准确性需要理论计算验证Matlab Models文件夹的代码基于涡旋光场的理论公式实现“涡旋光场传播模式计算”与“干涉场模拟”为FDTD仿真结果提供理论对照确保SAM-OAM转换效果的可靠性。4.1 涡旋光场传播模式计算代码该代码的核心功能是基于理论公式计算涡旋光束的光强与相位分布验证OAM光束的基本特性具体实现如下理论公式导入代码嵌入涡旋光场的电场表达式参考文档中给出部分公式形式如包含exp(i×l×θ)螺旋相位项、高斯振幅项与传播相位项其中l为拓扑荷数θ为极角ρ为径向距离计算网格生成通过meshgrid函数生成x-y平面的计算网格网格范围与FDTD仿真的超表面半径匹配采样点数量确保计算精度并将直角坐标转换为极坐标theta, rho为公式计算提供参数光强与相位计算代码代入不同的拓扑荷数l如l1、2、5计算对应的电场分布进而得到光强分布|E|²与相位分布angle(E)可视化输出通过surf函数绘制涡旋光强的三维分布图呈现“甜甜圈”型通过pcolor函数绘制相位分布图呈现螺旋条纹并添加颜色条与坐标轴标签确保结果可直接用于分析与展示。4.2 干涉场模拟代码参考文档指出涡旋光束与其他光束平面波、球面波的干涉条纹可直观反映其拓扑荷数代码基于理论公式实现三类干涉场的模拟具体功能如下球面波干涉模拟代码导入参考文档给出的球面波干涉公式I2E₀²[1cos(lθ k(x²y²)/2)]k为波数代入不同l值计算并绘制干涉条纹图——条纹呈现螺旋状条纹数量与l一致可通过条纹方向判断l的正负平面波干涉模拟基于参考文档中的平面波干涉公式IE₁²E₂²2E₁E₂cos(lθ - kx sinα)α为平面波入射角代码设置α为固定值如π/6计算干涉场分布——干涉条纹为倾斜的直线或曲线倾斜程度与l相关MZ干涉仪干涉模拟代码构建马赫-曾德尔MZ干涉仪的理论模型计算涡旋光在干涉仪中的干涉结果参考文档未给出具体公式但代码会基于干涉原理推导输出明暗交替的干涉条纹用于验证OAM光束的相干性。4.3 代码功能价值理论与仿真的双向验证Matlab Models文件夹的代码本质是“FDTD仿真的理论校验工具”——通过将理论计算的光强、相位、干涉条纹与FDTD仿真结果对比若两者高度一致如拓扑荷数相同、条纹形态相似则说明FDTD仿真模型准确可靠若存在偏差可反向排查FDTD模型的参数设置如单元旋转角度、入射光偏振态确保整个项目的结果可信度。五、代码功能的协同性与应用场景5.1 各文件夹代码的协同逻辑整个项目的代码并非独立运行而是形成“递进式协同”关系基础层Nanofins文件夹代码完成单元验证输出最优单元参数为FDTD建模提供“合格单元”核心层FDTD models文件夹代码基于“合格单元”构建超表面模型完成SAM-OAM转换仿真输出远场数据验证层Matlab Models文件夹代码基于理论公式计算结果与FDTD远场数据对比验证仿真准确性。这种协同逻辑确保了从“单元”到“整体”再到“理论”的全流程可靠性避免单一环节的误差传递至最终结果。5.2 代码功能的典型应用场景基于代码的核心功能其可应用于以下场景超构表面设计优化通过Nanofins代码调整单元尺寸、材料通过FDTD代码测试不同拓扑荷数的转换效果筛选最优超构表面结构OAM光束特性研究利用Matlab代码模拟不同l值的涡旋光场结合FDTD仿真分析拓扑荷数对光束传播、干涉特性的影响光通信技术开发通过FDTD的空间复用代码设计可同时传输多个OAM光束的超构表面为大容量光通信提供技术支撑学术研究与论文产出代码输出的远场图、干涉条纹图可直接用于学术论文理论与仿真的对比数据可增强研究结论的说服力。六、代码功能的关键注意事项参数一致性Nanofins文件夹的单元周期、材料参数需与FDTD models文件夹保持一致否则会导致超表面建模错误远场条件满足FDTD的farfield.lsf脚本中远场距离f需严格满足fR²/λ否则计算的远场数据会偏离真实情况计算精度平衡扫参步长、网格采样点数量需在“精度”与“耗时”间平衡——步长过小、采样点过多会导致计算时间大幅增加需根据实际需求调整结果验证标准代码输出结果需符合参考文档中的特征如相位与旋转角度2倍关系、OAM光束拓扑荷数l2q、空间复用的双向转换若不符合需检查代码参数设置或公式导入是否正确。七、总结本项目的代码围绕“自旋-轨道角动量转换”这一核心目标构建了模块化、协同化的功能体系——Nanofins代码保障单元性能FDTD代码实现核心转换Matlab代码验证理论准确性。所有代码功能均严格遵循参考文档的技术逻辑未涉及核心代码的细节泄露仅聚焦“功能实现流程”与“结果输出价值”。通过这套代码可高效完成SAM-OAM转换超构表面的设计、仿真与验证为相关领域的研究与应用提供可靠的技术工具。