轨道角动量OAM超表面与自旋-轨道角动量耦合结构设计的FDTD仿真案例
轨道角动量 OAM 超表面 自旋-轨道角动量耦合结构设计 fdtd仿真 复现2017 OE：Spin-to-orbital angular momentum conversion in dieletric metasurface 介绍：自旋-轨道角动量转换超表面，入射自旋角动量光束，即左右旋圆偏振光时，经过介质超表面的调制，出射光为交叉偏振态的涡旋光束，携带轨道角动量，相位调控为几何相位 案例包括fdtd模型，建模脚本，远场涡旋光束的强度和相位计算脚本，空间复用超表面设计脚本，二氧化钛纳米结构单元仿真等，及一份word教程；
自旋-轨道角动量耦合超构表面仿真平台
—— 从单元设计到远场涡旋的全链路功能拆解
----------------------------------------------------------------
- 引言
在经典光学链路里，「相位」通常由光程差决定；而在亚波长尺度的介电超构表面（dielectric metasurface）上，相位可以由纳米结构的「朝向角」一次性赋予——这就是几何相（Pancharatnam–Berry phase）的核心思想。本文档围绕一套已验证的 Lumerical-FDTD + MATLAB 脚本集，系统梳理“自旋-轨道角动量（spin-to-orbital angular momentum, STOAM）转换”仿真平台的数据流、控制流与业务边界，帮助后续开发者快速理解“哪些参数可以动、哪些结果可信、哪些坑不能踩”。文中所有关键变量、文件名将用等宽字体标注，但不会暴露完整核心代码，仅保留足以定位功能的伪代码与片段。
----------------------------------------------------------------
- 业务全景图
用户目标：
“给定圆偏振泵浦光，让超构表面输出携带指定拓扑荷l的交叉圆偏振涡旋光束，并在远场验证其相位奇点。”
平台用例拆解：
- 单元级：验证单根 TiO₂ 纳米柱的交叉偏振转换效率与相位随旋转角θ的线性度。
- 阵列级：按l=2q规律排布朝向角，生成实际可打印的 GDS 等价模型。
- 场级：FDTD 模拟近→远场变换，提取强度与相位涡旋。
- 理论级：MATLAB 解析公式生成理想涡旋及干涉图样，用于“对标”仿真。
----------------------------------------------------------------
- 数据流与文件拓扑
以下用“节点→节点”描述数据依赖，括号内为典型文件或变量名。
单元 sweep 数据流
unit_TiO2.fsp (参数化模型)
└─ theta 参数表 (0°–180°，步长 10°)
└─ 监视器Exz_y(交叉偏振分量)
└─unitsweep.lsf→ 绘制LT(θ)与L_phase(θ)曲线
└─ 预期：相位斜率 ≈ 2（PB 规律）
阵列生成数据流
OAM_polar.lsf
└─ 输入：拓扑荷乘子q(例：5)
└─ 业务规则：
for 每一环i
for 每一柱j
alpha = qphi_j360/(2π)
end
end
轨道角动量 OAM 超表面 自旋-轨道角动量耦合结构设计 fdtd仿真 复现2017 OE：Spin-to-orbital angular momentum conversion in dieletric metasurface 介绍：自旋-轨道角动量转换超表面，入射自旋角动量光束，即左右旋圆偏振光时，经过介质超表面的调制，出射光为交叉偏振态的涡旋光束，携带轨道角动量，相位调控为几何相位 案例包括fdtd模型，建模脚本，远场涡旋光束的强度和相位计算脚本，空间复用超表面设计脚本，二氧化钛纳米结构单元仿真等，及一份word教程；
└─ 输出：结构组structure（材料 TiO₂，高度 600 nm）
远场验证数据流
farfield.lsf
└─ 输入：监视器monitor（近场复振幅E）
└─ 处理：
E_xy = farfieldexact3d(..., z = f)
ELCP = Ex – 1i*E_y
ERCP = Ex + 1i*E_y
└─ 输出：强度ILCP、相位phaseLCP（100×100 采样）
理论对照数据流
Vortex_beam.m
└─ 函数E(rho,z,theta,l)
└─ 可视化：|E|^2与干涉条纹（球面波、平面波、MZ 干涉仪三种）
----------------------------------------------------------------
- 关键控制参数与边界
| 参数 | 物理含义 | 推荐范围 | 越界风险 |
|---|---|---|---|
q | 拓扑荷乘子 | 0.5–5 | 过大 (>6) 时相位梯度 > 100°/μm， fabrication 不可控 |
period | 径向周期 | 325 nm | 小于 280 nm 会出现高阶衍射，大于 400 nm 欠采样 |
nanofin x_span | 柱宽 | 250 nm | 宽 < 200 nm 效率骤降；宽 > 300 nm 支持多模，相位非 2θ |
z_max结构高 | 柱高 | 600 nm | 高 < 500 nm 相位滞后不足；高 > 700 nm 刻蚀深宽比 > 10:1，工艺难 |
waist radius w0 | 高斯源束腰 | 5 μm | 束腰 > 15 μm 时远场涡旋环半径过大，采样窗不足 |
----------------------------------------------------------------
- 单元验证：确保“2θ”线性度
目标：证明单柱的交叉偏振相位 = 2×几何旋转角。
操作序列：
- 打开
unit_TiO2.fsp，确认材料TiO2 (Palik)已加载。 - 在“Optimizations & Sweeps”里运行
thetasweep（19 点）。 - 运行
unit_sweep.lsf，自动绘制：
- 左轴：归一化透射
L_T（应 > 0.85）； - 右轴：相位斜率
Lphase/max(Lphase)（应 ≈ 直线，斜率 1 归一化后）。
通过条件：决定系数R²> 0.98，且截距 < 5°。
未通过：检查柱宽/周期是否落在波导谐振区，或材料折射率是否漂移。
----------------------------------------------------------------
- 阵列生成：环状排布算法
输入：lens_radius = 13 µm，q = 5
算法要点（伪代码）：
np = floor(lens_radius / period) for i = 1:np nphi = 6*i # 每环柱数 ∝ 周长 for j = 1:nphi phi = 2π*(j-1)/nphi alpha = q*phi*180/π # 几何相位 = 2*alpha place_rect('x', r*cos(phi), 'y', r*sin(phi), 'rotation', alpha) end end输出：结构组structure，可直接runFDTD。
内存提示：q=5 时总柱数 ≈ 6×(1+2+…+40)= 4 920，网格数 ≈ 1800×1800×80，需 GPU ≥ 32 GB。
----------------------------------------------------------------
- 远场涡旋诊断
运行farfield.lsf后，自动弹出四张图：
I_LCP：应呈“甜甜圈”状，暗斑对比度 > 20 dB。phase_LCP：绕奇点一周相位变化 2πl，计数器顺时针递增。I_RCP：应几乎为均匀背景（< 0.05×环峰值），证明交叉偏振消光比合格。- 相位展开图：无跳变、无分叉，确认网格采样足够。
诊断脚本隐藏技巧：
- 用
unwrap(phase,[],2)前先mask = ILCP > max(ILCP)*0.1，排除噪声区。 - 若暗斑偏移中心 > 1 µm，说明高斯源未对准，需回退检查
waist position。
----------------------------------------------------------------
- 理论对照与误差预算
MATLAB 端提供三种干涉模型：
- 球面波干涉：
I = 2E₀²[1+cos(lθ + k(x²+y²)/2f)] - 平面波干涉：
I = E₁²+E₂²+2E₁E₂cos(lθ – kx sinα) - MZ 干涉仪：
I = 4I₀(r²/σ²)^L exp(–2r²/σ²)cos²[arg((x+iy)^L)]
误差定义：
ε = |lmeasured – ltarget| / l_target
经验值：
- 仿真–理论：ε < 3 %（采样 256×256）
- 实验–仿真：若加工 RMS 粗糙度 < 3 nm，ε 可控制在 8 % 以内。
----------------------------------------------------------------
- 常见故障速查表
| 现象 | 可能根因 | 快速验证 |
|---|---|---|
| 远场环断裂 | 柱高误差导致 2π 跳变 | 单柱 sweep 看相位是否仍线性 |
| 暗斑对比 < 10 dB | 源偏振纯度不足 | 加polarizer监视器，看 ER |
| 相位图多奇点 | 周期 > 0.4 µm，出现栅瓣 | 降周期至 0.32 µm 以下 |
| 运行时报“out of memory” | 结构组未加addtogroup | 检查selectpartial后是否及时清理 |
----------------------------------------------------------------
- 扩展路线（不修改核心代码）
- 双荷复用：已提供
OAMpolar2.lsf，奇偶环分别赋q1=2.5 / q2=5，可生成l=±5双涡旋。 - 波长调谐：保持柱宽/高不变，仅改
wavelength，看相位斜率漂移，可评估带宽。 - 逆向设计：把
thetasweep 换成optimization节点，以max(L_T)为目标，可自动搜最优柱尺寸。
----------------------------------------------------------------
- 小结
整套平台把“单元验证→阵列生成→近场仿真→远场诊断→理论对照”五个环节封装成可脚本化、可定量比对的流水线。开发者只要守住“PB 相位=2×旋转角”这一金线，即可在任意波长、任意拓扑荷场景下快速复现 STOAM 效应，而无需深入 FDTD 内核。文中给出的参数边界、误差预算与故障速查表，可显著降低二次开发时的试错成本。