宽带任意阶 贝塞尔光束 超表面 模型 fdtd仿真 复现论文2017年Light ScienceApplicationsGeneration of wavelength-independent subwavelength Bessel beams using metasurfaces 论文介绍介质超表面实现宽带任意阶贝塞尔光束的产生贝塞尔光束是无衍射光束的一种可以在较长的传播距离内保持很好的横向分布特性广泛应用于例子操控、成像等领域 案例内容主要包括文章的两个贝塞尔光束模型一个零阶贝塞尔光束一个一阶贝塞尔光束采用二氧化钛介质单元执行几何相位来构建 案例包括fdtd模型、fdtd设计脚本、Matlab计算代码和复现结果以及一份word教程附带计算超表面的远场光场分布的脚本可以得到任意位置的贝塞尔光场。1. 系统概述与架构设计本代码库实现了一个完整的超表面光学仿真系统专门用于生成和分析任意阶数的贝塞尔光束。系统采用模块化设计主要由三个功能模块组成MATLAB相位计算模块、FDTD结构构建模块和远场分析模块。这三个模块相互协作完成了从理论计算到物理仿真再到结果分析的全流程。宽带任意阶 贝塞尔光束 超表面 模型 fdtd仿真 复现论文2017年Light ScienceApplicationsGeneration of wavelength-independent subwavelength Bessel beams using metasurfaces 论文介绍介质超表面实现宽带任意阶贝塞尔光束的产生贝塞尔光束是无衍射光束的一种可以在较长的传播距离内保持很好的横向分布特性广泛应用于例子操控、成像等领域 案例内容主要包括文章的两个贝塞尔光束模型一个零阶贝塞尔光束一个一阶贝塞尔光束采用二氧化钛介质单元执行几何相位来构建 案例包括fdtd模型、fdtd设计脚本、Matlab计算代码和复现结果以及一份word教程附带计算超表面的远场光场分布的脚本可以得到任意位置的贝塞尔光场。从代码组织结构来看系统清晰地分离了不同阶段的任务MATLAB代码负责前期的数学计算和相位分布生成FDTD的structurelens脚本负责将数学模型转化为物理结构而metalensfarfield系列脚本则负责后处理和数据可视化。这种架构体现了良好的工程实践使得每个模块可以独立开发和测试。2. MATLAB相位计算模块深度分析2.1 核心算法实现phase_Bessel.m文件是整个系统的理论基础其核心功能是基于贝塞尔光束的物理原理计算超表面所需的相位分布。代码采用了标准的贝塞尔光束相位公式phi_Bessel0_405(i,j) 2*pi-2*pi/lamda1*sqrt(x1(i)^2y1(j)^2)*NA1l1*atan2(y1(j),x1(i));这个公式包含三个关键部分基础相位项2*pi提供相位基准径向相位项2pi/lamda1sqrt(x1(i)^2y1(j)^2)*NA1实现光束的锥形波前这是贝塞尔光束无衍射特性的来源角向相位项l1*atan2(y1(j),x1(i))引入轨道角动量决定光束的拓扑荷数2.2 参数处理流程代码中定义了两组工作参数分别对应405nm和532nm波长NA1 0.5; NA2 0.3; lamda1 0.405e-6; lamda2 0.532e-6; period1 0.155e-6; period2 0.25e-6;参数处理体现了实际光学系统的约束数值孔径NA影响光束的锥角和传播特性单元周期必须小于工作波长以满足采样定理拓扑荷数l直接控制贝塞尔光束的阶数2.3 数据后处理与输出计算得到的相位数据经过标准化处理phi_Bessel0_405 wrapTo2Pi(phi_Bessel0_405)/2/pi; target_angle1_405 round(phi_Bessel0_405/2*360);这个处理流程将理论相位转换为实际可用的纳米柱旋转角度其中wrapTo2Pi确保相位值在[0,2π]范围内除以2π进行归一化乘以360转换为角度制round四舍五入到整数角度便于制造3. FDTD结构构建模块功能详解3.1 超表面物理结构创建structure_lens.lsf脚本实现了从相位数据到物理结构的转换。代码首先加载MATLAB计算的角度数据matlabload(target_angle3_532.mat);然后通过双重循环在指定区域内放置纳米柱for(i1:np) { for (j1:np) { r sqrt(x_mask(i)^2y_mask(j)^2); if (r lens_radius) { addrect; set(rotation 1,target_angle3_532(i,j)); } } }这个循环逻辑确保只在圆形区域内放置结构r lens_radius每个位置根据预计算的角度旋转纳米柱保持超表面的周期性排列3.2 材料与几何参数设置纳米柱的具体尺寸参数体现了实际制造约束set(x span,0.21e-6); # 纳米柱X方向长度 set(y span,0.065e-6); # 纳米柱Y方向长度 set(z min,0); set(z max,0.6e-6); # 纳米柱高度各向异性的设计0.21μm vs 0.065μm是实现几何相位调制的关键这种尺寸差异使得纳米柱对不同偏振的光产生不同的相位响应。3.3 光源系统设计代码中设置了两个正交的TFSTF光源# 第一个光源 - X偏振 addtfsf; set(name,source_x); # 第二个光源 - Y偏振带90度相位延迟 addtfsf; set(name,source_y); set(polarization angle,90); set(phase,90);这种配置实际上产生了圆偏振光入射是激发几何相位的标准方法。两个光源的联合作用确保了超表面能够正确响应并产生所需的贝塞尔光束。4. 远场分析模块功能解析4.1 计算区域定义metalens_farfield系列脚本定义了详细的仿真分析区域L 25e-6; x linspace(-L/2,L/2,200); y linspace(-L/2,L/2,200); z linspace(1e-6,1.5*d,200);参数选择考虑了计算效率和物理需求的平衡25μm计算区域足够覆盖光束的主要能量区域200个采样点在精度和计算成本间取得平衡1.5倍焦距范围充分观察光束的传播特性4.2 远场变换算法核心的远场计算使用FDTD的内置函数E_xz farfieldexact3d(monitor,x,0,z); E_xz pinch(sum(abs(E_xz)^2,4));farfielddexact3d函数执行严格的矢量衍射计算将近场结果变换到远场。随后的处理步骤abs(E_xz)^2计算电场强度sum(...,4)对偏振维度求和得到总强度pinch移除单例维度准备可视化4.3 多平面分析策略代码分别在两个关键平面进行分析# X-Z平面传播截面 E_xz farfieldexact3d(monitor,x,0,z); # X-Y平面焦平面 E_xy farfieldexact3d(monitor,x,y,28e-6);这种多平面分析策略能够全面表征贝塞尔光束的特性X-Z平面展示光束的无衍射传播范围X-Y平面显示焦平面的强度分布和相位特征5. 代码的工程实践价值5.1 参数化设计思想整个代码库体现了良好的参数化设计关键参数都在文件开头明确定义便于修改和优化。例如波长、数值孔径、拓扑荷数等关键参数都可以直接调整而不需要深入理解算法细节。5.2 模块间数据接口系统设计了清晰的数据接口MATLAB输出.mat格式的角度数据FDTD脚本直接加载这些数据文件统一的坐标系统和单位制全部使用米制5.3 误差控制机制代码中多处体现了数值误差控制角度值的四舍五入处理相位值的周期性包裹计算区域的合理截断6. 物理效应与代码实现的对应关系通过分析代码我们可以清楚地看到物理原理如何转化为具体的算法实现几何相位原理→ 纳米柱旋转角度数组贝塞尔光束方程→ 相位分布计算函数衍射传播理论→ 远场变换算法偏振光学→ 正交光源设置这种从物理理论到代码实现的直接映射使得整个系统具有很好的可解释性和可扩展性。7. 总结这个代码库提供了一个完整、实用的超表面光学仿真解决方案。通过深入分析每个模块的功能实现我们不仅理解了如何用代码生成贝塞尔光束更重要的是学习了如何将复杂的光学理论转化为可执行的计算流程。代码中体现的模块化设计、参数化思想和物理准确性为类似的光子器件仿真项目提供了有价值的参考模板。
宽带任意阶贝塞尔光束生成：超表面技术与FDTD仿真实践
宽带任意阶 贝塞尔光束 超表面 模型 fdtd仿真 复现论文2017年Light ScienceApplicationsGeneration of wavelength-independent subwavelength Bessel beams using metasurfaces 论文介绍介质超表面实现宽带任意阶贝塞尔光束的产生贝塞尔光束是无衍射光束的一种可以在较长的传播距离内保持很好的横向分布特性广泛应用于例子操控、成像等领域 案例内容主要包括文章的两个贝塞尔光束模型一个零阶贝塞尔光束一个一阶贝塞尔光束采用二氧化钛介质单元执行几何相位来构建 案例包括fdtd模型、fdtd设计脚本、Matlab计算代码和复现结果以及一份word教程附带计算超表面的远场光场分布的脚本可以得到任意位置的贝塞尔光场。1. 系统概述与架构设计本代码库实现了一个完整的超表面光学仿真系统专门用于生成和分析任意阶数的贝塞尔光束。系统采用模块化设计主要由三个功能模块组成MATLAB相位计算模块、FDTD结构构建模块和远场分析模块。这三个模块相互协作完成了从理论计算到物理仿真再到结果分析的全流程。宽带任意阶 贝塞尔光束 超表面 模型 fdtd仿真 复现论文2017年Light ScienceApplicationsGeneration of wavelength-independent subwavelength Bessel beams using metasurfaces 论文介绍介质超表面实现宽带任意阶贝塞尔光束的产生贝塞尔光束是无衍射光束的一种可以在较长的传播距离内保持很好的横向分布特性广泛应用于例子操控、成像等领域 案例内容主要包括文章的两个贝塞尔光束模型一个零阶贝塞尔光束一个一阶贝塞尔光束采用二氧化钛介质单元执行几何相位来构建 案例包括fdtd模型、fdtd设计脚本、Matlab计算代码和复现结果以及一份word教程附带计算超表面的远场光场分布的脚本可以得到任意位置的贝塞尔光场。从代码组织结构来看系统清晰地分离了不同阶段的任务MATLAB代码负责前期的数学计算和相位分布生成FDTD的structurelens脚本负责将数学模型转化为物理结构而metalensfarfield系列脚本则负责后处理和数据可视化。这种架构体现了良好的工程实践使得每个模块可以独立开发和测试。2. MATLAB相位计算模块深度分析2.1 核心算法实现phase_Bessel.m文件是整个系统的理论基础其核心功能是基于贝塞尔光束的物理原理计算超表面所需的相位分布。代码采用了标准的贝塞尔光束相位公式phi_Bessel0_405(i,j) 2*pi-2*pi/lamda1*sqrt(x1(i)^2y1(j)^2)*NA1l1*atan2(y1(j),x1(i));这个公式包含三个关键部分基础相位项2*pi提供相位基准径向相位项2pi/lamda1sqrt(x1(i)^2y1(j)^2)*NA1实现光束的锥形波前这是贝塞尔光束无衍射特性的来源角向相位项l1*atan2(y1(j),x1(i))引入轨道角动量决定光束的拓扑荷数2.2 参数处理流程代码中定义了两组工作参数分别对应405nm和532nm波长NA1 0.5; NA2 0.3; lamda1 0.405e-6; lamda2 0.532e-6; period1 0.155e-6; period2 0.25e-6;参数处理体现了实际光学系统的约束数值孔径NA影响光束的锥角和传播特性单元周期必须小于工作波长以满足采样定理拓扑荷数l直接控制贝塞尔光束的阶数2.3 数据后处理与输出计算得到的相位数据经过标准化处理phi_Bessel0_405 wrapTo2Pi(phi_Bessel0_405)/2/pi; target_angle1_405 round(phi_Bessel0_405/2*360);这个处理流程将理论相位转换为实际可用的纳米柱旋转角度其中wrapTo2Pi确保相位值在[0,2π]范围内除以2π进行归一化乘以360转换为角度制round四舍五入到整数角度便于制造3. FDTD结构构建模块功能详解3.1 超表面物理结构创建structure_lens.lsf脚本实现了从相位数据到物理结构的转换。代码首先加载MATLAB计算的角度数据matlabload(target_angle3_532.mat);然后通过双重循环在指定区域内放置纳米柱for(i1:np) { for (j1:np) { r sqrt(x_mask(i)^2y_mask(j)^2); if (r lens_radius) { addrect; set(rotation 1,target_angle3_532(i,j)); } } }这个循环逻辑确保只在圆形区域内放置结构r lens_radius每个位置根据预计算的角度旋转纳米柱保持超表面的周期性排列3.2 材料与几何参数设置纳米柱的具体尺寸参数体现了实际制造约束set(x span,0.21e-6); # 纳米柱X方向长度 set(y span,0.065e-6); # 纳米柱Y方向长度 set(z min,0); set(z max,0.6e-6); # 纳米柱高度各向异性的设计0.21μm vs 0.065μm是实现几何相位调制的关键这种尺寸差异使得纳米柱对不同偏振的光产生不同的相位响应。3.3 光源系统设计代码中设置了两个正交的TFSTF光源# 第一个光源 - X偏振 addtfsf; set(name,source_x); # 第二个光源 - Y偏振带90度相位延迟 addtfsf; set(name,source_y); set(polarization angle,90); set(phase,90);这种配置实际上产生了圆偏振光入射是激发几何相位的标准方法。两个光源的联合作用确保了超表面能够正确响应并产生所需的贝塞尔光束。4. 远场分析模块功能解析4.1 计算区域定义metalens_farfield系列脚本定义了详细的仿真分析区域L 25e-6; x linspace(-L/2,L/2,200); y linspace(-L/2,L/2,200); z linspace(1e-6,1.5*d,200);参数选择考虑了计算效率和物理需求的平衡25μm计算区域足够覆盖光束的主要能量区域200个采样点在精度和计算成本间取得平衡1.5倍焦距范围充分观察光束的传播特性4.2 远场变换算法核心的远场计算使用FDTD的内置函数E_xz farfieldexact3d(monitor,x,0,z); E_xz pinch(sum(abs(E_xz)^2,4));farfielddexact3d函数执行严格的矢量衍射计算将近场结果变换到远场。随后的处理步骤abs(E_xz)^2计算电场强度sum(...,4)对偏振维度求和得到总强度pinch移除单例维度准备可视化4.3 多平面分析策略代码分别在两个关键平面进行分析# X-Z平面传播截面 E_xz farfieldexact3d(monitor,x,0,z); # X-Y平面焦平面 E_xy farfieldexact3d(monitor,x,y,28e-6);这种多平面分析策略能够全面表征贝塞尔光束的特性X-Z平面展示光束的无衍射传播范围X-Y平面显示焦平面的强度分布和相位特征5. 代码的工程实践价值5.1 参数化设计思想整个代码库体现了良好的参数化设计关键参数都在文件开头明确定义便于修改和优化。例如波长、数值孔径、拓扑荷数等关键参数都可以直接调整而不需要深入理解算法细节。5.2 模块间数据接口系统设计了清晰的数据接口MATLAB输出.mat格式的角度数据FDTD脚本直接加载这些数据文件统一的坐标系统和单位制全部使用米制5.3 误差控制机制代码中多处体现了数值误差控制角度值的四舍五入处理相位值的周期性包裹计算区域的合理截断6. 物理效应与代码实现的对应关系通过分析代码我们可以清楚地看到物理原理如何转化为具体的算法实现几何相位原理→ 纳米柱旋转角度数组贝塞尔光束方程→ 相位分布计算函数衍射传播理论→ 远场变换算法偏振光学→ 正交光源设置这种从物理理论到代码实现的直接映射使得整个系统具有很好的可解释性和可扩展性。7. 总结这个代码库提供了一个完整、实用的超表面光学仿真解决方案。通过深入分析每个模块的功能实现我们不仅理解了如何用代码生成贝塞尔光束更重要的是学习了如何将复杂的光学理论转化为可执行的计算流程。代码中体现的模块化设计、参数化思想和物理准确性为类似的光子器件仿真项目提供了有价值的参考模板。
