轨道角动量 OAM 超表面 自旋-轨道角动量耦合结构设计 fdtd仿真 复现2017 OESpin-to-orbital angular momentum conversion in dieletric metasurface 介绍自旋-轨道角动量转换超表面入射自旋角动量光束即左右旋圆偏振光时经过介质超表面的调制出射光为交叉偏振态的涡旋光束携带轨道角动量相位调控为几何相位 案例包括fdtd模型建模脚本远场涡旋光束的强度和相位计算脚本空间复用超表面设计脚本二氧化钛纳米结构单元仿真等及一份word教程一、项目概述本项目围绕自旋-轨道角动量转换超构表面展开通过FDTD时域有限差分法仿真与Matlab理论计算实现了超构表面结构建模、涡旋光束生成及干涉场分析。项目核心目标是验证介质超表面对入射光自旋角动量的调控能力使其转换为特定拓扑荷数的轨道角动量为光通信、量子光学等领域的涡旋光应用提供技术支撑。轨道角动量 OAM 超表面 自旋-轨道角动量耦合结构设计 fdtd仿真 复现2017 OESpin-to-orbital angular momentum conversion in dieletric metasurface 介绍自旋-轨道角动量转换超表面入射自旋角动量光束即左右旋圆偏振光时经过介质超表面的调制出射光为交叉偏振态的涡旋光束携带轨道角动量相位调控为几何相位 案例包括fdtd模型建模脚本远场涡旋光束的强度和相位计算脚本空间复用超表面设计脚本二氧化钛纳米结构单元仿真等及一份word教程项目文件分为三大核心模块Nanofins单元纳米柱仿真模块、FDTD Models超表面全场仿真模块、Matlab Models涡旋光理论计算模块各模块协同完成“单元验证-结构建模-性能仿真-理论验证”的完整技术流程。二、各模块功能解析2.1 Nanofins单元纳米柱仿真模块2.1.1 核心定位单元纳米柱是超构表面的基本功能单元本模块通过仿真验证单个纳米柱的偏振调控特性确保其满足自旋-轨道角动量转换所需的相位调控能力为后续超构表面阵列设计提供参数依据。2.1.2 关键功能与流程参数配置与扫参设置- 基础参数采用TiO₂作为纳米柱材料高折射率介质利于相位调控SiO₂玻璃作为基底根据文献设置纳米柱周期、高度、横截面尺寸等关键参数。- 扫参逻辑通过sweep功能对纳米柱的旋转角度θ进行扫描0°-180°监测正交偏振态下的相位变化与透射率验证“相位-旋转角度”的2倍关系几何相位调控核心原理。数据处理与结果可视化- 运行unitsweep.lsf脚本读取扫参数据旋转角度θ、正交偏振相位Lphase、透射率L_T并生成两类曲线- 透射率曲线反映不同旋转角度下纳米柱的光透射效率确保核心工作角度范围内透射率稳定避免能量损失影响后续仿真。- 归一化相位曲线验证相位变化与旋转角度的线性关系且斜率为2几何相位调控的关键特征直接决定后续轨道角动量转换精度。场分布监测通过Exz方向电场监视器观察纳米柱内部及近场的电场分布分析电磁场与纳米柱的相互作用排除结构设计缺陷如边缘场畸变、模式杂散等。2.2 FDTD Models超构表面全场仿真模块2.2.1 核心定位基于单元纳米柱的验证结果构建超构表面阵列仿真完整的自旋-轨道角动量转换过程生成携带特定拓扑荷数的涡旋光束并通过远场监测获取光强与相位分布。2.2.2 模块文件构成与功能分工本模块包含两类核心脚本超构表面建模脚本OAMpolar.lsf、OAMpolar_2.lsf与远场分析脚本farfield.lsf覆盖“单一拓扑荷转换”与“空间复用双拓扑荷转换”两类应用场景。2.2.3 单一拓扑荷转换模型OAM_polar.lsf核心功能构建均匀超构表面阵列实现“入射自旋角动量→固定轨道角动量”的转换支持不同拓扑荷数l的定制化生成。关键流程与原理- 结构建模逻辑基底创建通过addcircle函数生成SiO₂圆形基底尺寸略大于超构表面避免边缘效应影响厚度设为2倍入射波长确保基底无透射损耗。纳米柱阵列生成采用双层循环构建六边形密排阵列兼顾高填充率与均匀性- 外层循环i控制阵列半径方向的环数根据“超构表面半径/纳米柱周期”计算最大环数np。- 内层循环j在每一环上生成6i个纳米柱六边形排列特征通过极坐标r周期×环数φ角度间隔计算每个纳米柱的中心坐标xmask, ymask。相位调控通过set(rotation 1, alpha)设置纳米柱旋转角度其中αq×(φ/2π×360)α₀q为调控参数φ为纳米柱极角α₀为初始相位偏移。根据文献原理轨道角动量拓扑荷数l2q通过修改q值可生成l1、2、5等不同拓扑荷的涡旋光。- 光源与监测设置- 光源采用两个正交偏振的高斯光源“x plane”偏振角0°、相位0°“y plane”偏振角90°、相位90°组合生成圆偏振入射光左旋/右旋可通过相位差调整光源腰斑半径与工作距离匹配超构表面尺寸确保均匀照明。- 监测器通过addpower设置2D Z-normal功率监测器位于z1e-6处超构表面上方记录近场光强分布为后续远场分析提供原始数据。拓扑荷数定制方法修改脚本中q的取值重新运行即可生成对应拓扑荷的超构表面。例如q2.5时l5q0.5时l1满足不同应用场景下的涡旋光需求。2.2.4 空间复用双拓扑荷转换模型OAM_polar_2.lsf核心功能突破单一拓扑荷限制通过空间复用技术不同区域采用不同旋转角度规律实现“入射自旋方向→不同拓扑荷”的差异化转换即入射右旋圆偏振光时出射l5的左旋圆偏振光入射左旋圆偏振光时出射l-5的右旋圆偏振光。关键创新点空间复用逻辑- 环数分区控制通过mod(i,2)1奇数环与mod(i,2)0偶数环判断环数类型对不同环的纳米柱采用不同旋转角度规律- 奇数环采用q12.5计算旋转角度α调控逻辑与单一拓扑荷模型一致。- 偶数环采用q25计算旋转角度并额外添加负号set(rotation 1, -alpha)通过旋转方向的反向实现拓扑荷的符号调控。- 结果匹配性该设计与文献Figure5的结构一致仿真结果可直接验证空间复用技术在多拓扑荷转换中的可行性为高密度光通信如多路涡旋光复用提供技术参考。2.2.5 远场分析脚本farfield.lsf核心功能将监测器获取的近场数据转换为远场光强与相位分布直观验证涡旋光束的拓扑荷特性如相位奇点、螺旋光强分布。关键流程- 数据转换调用farfieldexact3d函数基于近场监测数据E场计算远场分布支持x-z平面传播方向与x-y平面垂直传播方向的场分布分析。- 偏振分离与分析- 左旋圆偏振光LCP通过Ex - iEy计算提取光强|ELCP|²与相位angle(ELCP)验证其拓扑荷数相位呈螺旋状分布奇点位于中心。- 右旋圆偏振光RCP通过Ex iEy计算同理分析其场分布验证交叉偏振转换效果入射RCP时仅LCP远场有信号且携带目标拓扑荷。2.3 Matlab Models涡旋光理论计算模块2.3.1 核心定位通过理论公式计算涡旋光束的传播特性及干涉场分布与FDTD仿真结果进行对比验证仿真的准确性同时为超构表面设计提供理论预研支持如拓扑荷数与干涉条纹的对应关系。2.3.2 核心功能与理论公式涡旋光束传播模式计算Vortex_beam.m- 核心公式基于高斯涡旋光场理论考虑光束在自由空间中的传播损耗与相位演变公式如下\[E(\rho,z,\theta,l) -i^{(l1)} \cdot \left(\frac{r(z)}{r}\right)^{|l|} \cdot \exp\left(ikz \frac{i2z\rho^2}{kr^2r(z)^2}\right) \cdot 10 \cdot \left(\frac{\rho^2}{r(z)}\right)^{|l|} \cdot \exp\left(-\frac{\rho^2}{r(z)^2}\right) \cdot \exp(i|l|\theta)\]其中l为拓扑荷数ρ为径向距离θ为极角z为传播距离k为波数r(z)为传播后的光束腰斑半径。- 结果输出通过surf函数绘制x-y平面的光强分布|E|²直观展示涡旋光的“甜甜圈”状光强中心为相位奇点光强为0支持不同l值如l2、5的场分布对比。干涉场计算- 球面波干涉模拟涡旋光与球面波的干涉公式为\[I 2E0^2 \left[1 \cos\left(l\theta k\frac{x^2y^2}{2}\right)\right]\]干涉条纹呈螺旋状条纹数量与拓扑荷数l正相关可用于快速识别涡旋光的拓扑荷。- 平面波干涉模拟涡旋光与倾斜平面波的干涉公式为\[I E1^2 E2^2 2E1E2\cos\left(l\theta - kx\sin\alpha\right)\]其中α为平面波入射角干涉条纹呈周期性偏移的螺旋状可用于分析涡旋光的相位梯度。- MZ干涉仪干涉模拟涡旋光在马赫-曾德尔干涉仪中的干涉公式为\[I 4I0 \left(\frac{x^2y^2}{\delta^2}\right)^L \cdot \exp\left(-\frac{2(x^2y^2)}{\delta^2}\right) \cdot \cos^2\left(\text{angle}((xiy)^L)\right)\]其中δ为光束腰斑相关参数L为拓扑荷数干涉结果可用于涡旋光的定量检测如拓扑荷数标定。三、整体工作流程与应用场景3.1 标准工作流程单元验证阶段运行Nanofins模块确认纳米柱的“相位-旋转角度”关系与透射率满足要求输出合格的单元参数如周期、旋转角度范围。超构表面建模与仿真阶段- 单一拓扑荷需求修改OAMpolar.lsf的q值生成对应超构表面运行FDTD仿真获取近场数据。- 双拓扑荷需求运行OAMpolar_2.lsf生成空间复用超构表面执行仿真。远场分析阶段运行farfield.lsf将近场数据转换为远场光强与相位图验证涡旋光拓扑荷的正确性。理论验证阶段运行Matlab Models模块计算相同参数下的理论场分布与FDTD仿真结果对比确保仿真精度如光强分布形态、相位奇点位置一致。3.2 典型应用场景光通信领域利用不同拓扑荷的涡旋光如l1、3、5作为信息载体实现多路信号的并行传输提升通信带宽。量子光学领域通过空间复用超构表面生成l±5的涡旋光用于量子比特的编码与操控。显微成像领域利用涡旋光的相位奇点特性实现超分辨成像如涡旋光显微镜突破衍射极限。四、关键技术特点与注意事项4.1 技术特点高精度相位调控基于几何相位原理通过纳米柱旋转角度实现2倍相位调控转换精度高拓扑荷数误差0.1。高兼容性支持不同拓扑荷数的快速定制仅需修改q值且兼容左旋/右旋圆偏振入射光的交叉转换。理论-仿真闭环Matlab理论计算与FDTD仿真相互验证确保结果的可靠性与准确性。4.2 注意事项参数一致性FDTD仿真中的波长默认0.532μm、材料参数TiO₂、SiO₂的折射率需与实验或文献保持一致避免因参数偏差导致转换效率下降。边界条件设置FDTD区域的x/y/z范围需覆盖超构表面及近场区域避免边界反射对场分布的干扰。计算资源配置超构表面阵列规模较大如半径13μm、周期0.325μm约千级纳米柱需合理配置计算资源如网格精度、时间步长平衡仿真精度与计算效率。五、总结本项目通过模块化设计实现了自旋-轨道角动量转换超构表面的全流程仿真与理论验证从单元纳米柱的特性验证到超构表面的阵列建模再到远场涡旋光的分析形成了完整的技术链条。代码的可扩展性强支持拓扑荷数定制、空间复用设计理论与仿真的结合确保了结果的可靠性可为相关领域的科研与工程应用提供直接的技术参考。