轨道角动量OAM超表面设计：自旋到轨道角动量转换与几何相位调控的FDTD仿真研究
轨道角动量 OAM 超表面 自旋-轨道角动量耦合结构设计 fdtd仿真 复现2017 OE：Spin-to-orbital angular momentum conversion in dieletric metasurface 介绍：自旋-轨道角动量转换超表面，入射自旋角动量光束，即左右旋圆偏振光时，经过介质超表面的调制，出射光为交叉偏振态的涡旋光束，携带轨道角动量，相位调控为几何相位 案例包括fdtd模型，建模脚本，远场涡旋光束的强度和相位计算脚本，空间复用超表面设计脚本，二氧化钛纳米结构单元仿真等，及一份word教程；
一、项目概述
本项目围绕自旋-轨道角动量转换超构表面展开，通过FDTD（时域有限差分法）仿真与Matlab理论计算，实现了超构表面结构建模、涡旋光束生成及干涉场分析。项目核心目标是验证介质超表面对入射光自旋角动量的调控能力，使其转换为特定拓扑荷数的轨道角动量，为光通信、量子光学等领域的涡旋光应用提供技术支撑。
轨道角动量 OAM 超表面 自旋-轨道角动量耦合结构设计 fdtd仿真 复现2017 OE：Spin-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为调控参数，φ为纳米柱极角，α₀为初始相位偏移）。根据文献原理，轨道角动量拓扑荷数l=2q，通过修改q值可生成l=1、2、5等不同拓扑荷的涡旋光。
- 光源与监测设置
- 光源：采用两个正交偏振的高斯光源（“x plane”偏振角0°、相位0°；“y plane”偏振角90°、相位90°），组合生成圆偏振入射光（左旋/右旋可通过相位差调整），光源腰斑半径与工作距离匹配超构表面尺寸，确保均匀照明。
- 监测器：通过addpower设置2D Z-normal功率监测器（位于z=1e-6处，超构表面上方），记录近场光强分布，为后续远场分析提供原始数据。 - 拓扑荷数定制方法：修改脚本中q的取值，重新运行即可生成对应拓扑荷的超构表面。例如：q=2.5时，l=5；q=0.5时，l=1，满足不同应用场景下的涡旋光需求。
2.2.4 空间复用双拓扑荷转换模型（OAM_polar_2.lsf）
- 核心功能：突破单一拓扑荷限制，通过空间复用技术（不同区域采用不同旋转角度规律），实现“入射自旋方向→不同拓扑荷”的差异化转换，即入射右旋圆偏振光时出射l=5的左旋圆偏振光，入射左旋圆偏振光时出射l=-5的右旋圆偏振光。
- 关键创新点：空间复用逻辑
- 环数分区控制：通过mod(i,2)==1（奇数环）与mod(i,2)==0（偶数环）判断环数类型，对不同环的纳米柱采用不同旋转角度规律：
- 奇数环：采用q1=2.5计算旋转角度α，调控逻辑与单一拓扑荷模型一致。
- 偶数环：采用q2=5计算旋转角度，并额外添加负号（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^{(l+1)} \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值（如l=2、5）的场分布对比。 - 干涉场计算
- 球面波干涉：模拟涡旋光与球面波的干涉，公式为：
\[
I = 2E0^2 \left[1 + \cos\left(l\theta + k\frac{x^2+y^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^2+y^2}{\delta^2}\right)^L \cdot \exp\left(-\frac{2(x^2+y^2)}{\delta^2}\right) \cdot \cos^2\left(\text{angle}((x+iy)^L)\right)
\]
其中δ为光束腰斑相关参数，L为拓扑荷数，干涉结果可用于涡旋光的定量检测（如拓扑荷数标定）。
三、整体工作流程与应用场景
3.1 标准工作流程
- 单元验证阶段：运行
Nanofins模块，确认纳米柱的“相位-旋转角度”关系与透射率满足要求，输出合格的单元参数（如周期、旋转角度范围）。 - 超构表面建模与仿真阶段：
- 单一拓扑荷需求：修改OAMpolar.lsf的q值，生成对应超构表面，运行FDTD仿真，获取近场数据。
- 双拓扑荷需求：运行OAMpolar_2.lsf，生成空间复用超构表面，执行仿真。 - 远场分析阶段：运行
farfield.lsf，将近场数据转换为远场光强与相位图，验证涡旋光拓扑荷的正确性。 - 理论验证阶段：运行
Matlab Models模块，计算相同参数下的理论场分布，与FDTD仿真结果对比，确保仿真精度（如光强分布形态、相位奇点位置一致）。
3.2 典型应用场景
- 光通信领域：利用不同拓扑荷的涡旋光（如l=1、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，约千级纳米柱），需合理配置计算资源（如网格精度、时间步长），平衡仿真精度与计算效率。
五、总结
本项目通过模块化设计，实现了自旋-轨道角动量转换超构表面的全流程仿真与理论验证，从单元纳米柱的特性验证到超构表面的阵列建模，再到远场涡旋光的分析，形成了完整的技术链条。代码的可扩展性强（支持拓扑荷数定制、空间复用设计），理论与仿真的结合确保了结果的可靠性，可为相关领域的科研与工程应用提供直接的技术参考。