超表面设计实现宽带任意阶贝塞尔光束-代码聚汇网

超表面设计实现宽带任意阶贝塞尔光束

采璇

1. 项目概述：宽带任意阶贝塞尔光束超表面设计与仿真

贝塞尔光束作为一种特殊的无衍射光束，在光学微操控、显微成像等领域具有重要应用价值。传统方法生成贝塞尔光束通常需要复杂的光学系统，而基于超表面的解决方案则能大幅简化系统结构。本项目复现了2017年发表在Light: Science & Applications上的研究成果，通过二氧化钛介质超表面实现宽带任意阶贝塞尔光束的生成。

超表面是一种由亚波长结构单元组成的二维平面光学元件，通过精心设计每个单元的几何参数，可以实现对光波前相位、振幅和偏振态的精确调控。本项目采用几何相位（Pancharatnam-Berry相位）原理，通过旋转矩形纳米柱的角度来实现所需的相位调制。这种方法的独特优势在于其宽带特性——由于几何相位与波长无关，设计的超表面可以在宽光谱范围内工作。

2. 核心原理与技术路线

2.1 贝塞尔光束的基本特性

贝塞尔光束最显著的特征是其无衍射性，即在传播过程中光束剖面保持不变的特性。数学上，n阶贝塞尔光束的电场分布可以表示为：

E(r,φ,z) = A₀Jₙ(kᵣr)exp(inφ)exp(ikₓz)

其中Jₙ是n阶贝塞尔函数，kᵣ和kₓ分别是径向和轴向波矢分量，φ是方位角。零阶贝塞尔光束(n=0)具有中心亮斑的强度分布，而高阶贝塞尔光束则呈现环形结构，并携带轨道角动量。

注意：实际应用中，理想的无限能量贝塞尔光束无法实现，我们通常通过有限孔径近似产生准贝塞尔光束。

2.2 几何相位超表面设计原理

几何相位，又称Pancharatnam-Berry相位，源于光波在偏振态空间中的几何演化。当圆偏振光通过一个光学轴方向旋转θ的波片时，出射光会获得±2θ的附加相位（符号取决于入射光的旋向性）。

本项目采用矩形二氧化钛纳米柱作为基本单元，通过旋转纳米柱的角度来调控几何相位。每个纳米柱相当于一个半波片，当左旋圆偏振光入射时，右旋圆偏振光分量会获得2θ的相位延迟，其中θ是纳米柱相对于参考方向的旋转角度。

3. 仿真模型构建与实现

3.1 单元结构设计与验证

首先需要验证单个纳米柱的相位调制能力。我们选择二氧化钛(TiO₂)作为纳米柱材料，其在可见光波段具有高折射率(n≈2.4)和低吸收特性。纳米柱的尺寸设计遵循以下原则：

高度h≈λ₀/(n-nₛ)，其中λ₀是中心波长，nₛ是衬底折射率
长宽比≠1，以确保偏振相关响应
周期P<λ₀/NA，避免高阶衍射

通过FDTD仿真扫描不同旋转角度下的相位延迟，可以验证相位与转角是否满足φ=±2θ的关系。典型的纳米柱参数可能为：长200nm，宽100nm，高600nm，周期400nm。

3.2 超表面相位分布计算

贝塞尔光束所需的相位分布包含两部分：

径向相位：φ_r(r) = -2πrNA/λ
方位角相位：φ_φ(φ) = nφ

其中r=√(x²+y²)是径向坐标，φ是方位角，n是拓扑荷数（决定贝塞尔光束的阶数），NA是数值孔径。总相位分布为这两部分的叠加：

φ(x,y) = mod(2π - 2πrNA/λ + nφ, 2π)

在MATLAB中实现这一计算的代码如下：

matlab复制% 参数设置
lambda = 532e-9;   % 波长(m)
NA = 0.7;          % 数值孔径
n = 1;             % 拓扑荷数
N = 500;           % 单元数
P = 400e-9;        % 周期(m)

% 坐标网格
x = linspace(-N/2,N/2,N)*P;
y = linspace(-N/2,N/2,N)*P;
[X,Y] = meshgrid(x,y);
[phi,r] = cart2pol(X,Y);

% 相位计算
phase = mod(2*pi - 2*pi*r*NA/lambda + n*phi, 2*pi);

% 转角转换
theta = phase/2;   % 几何相位关系

3.3 FDTD仿真模型搭建

在FDTD Solutions中搭建超表面仿真模型的关键步骤包括：

创建SiO₂衬底
根据MATLAB生成的转角数据阵列放置TiO₂纳米柱
设置圆偏振平面波光源
添加近场和远场监视器

一个典型的LSF脚本片段如下：

lua复制# 创建纳米柱阵列
for(i=1:Nx) {
    for(j=1:Ny) {
        theta = theta_matrix(i,j);  # 从MATLAB数据读取转角
        addrect;
        set("name",sprintf("pillar_%d_%d",i,j));
        set("x",x_pos(i));
        set("y",y_pos(j));
        set("z span",height);
        set("x span",length);
        set("y span",width);
        set("material","TiO2");
        set("rotation 1",theta);
    }
}

# 设置光源
addplane;
set("injection axis","z");
set("direction","forward");
set("wavelength start",lambda*1e9-50);
set("wavelength stop",lambda*1e9+50);
set("polarization angle",45);
set("phase",90);  # 圆偏振光

4. 结果分析与验证

4.1 零阶贝塞尔光束特性

对于n=0的零阶贝塞尔光束，仿真结果应显示：

近场相位呈同心圆环分布
远场强度分布为中心亮斑+多个同心圆环
轴向强度分布显示无衍射特性（一定范围内强度保持）

典型的仿真结果参数：

中心斑直径(FWHM)≈0.36λ/NA
无衍射距离Z_max≈R/NA，其中R是超表面半径

4.2 一阶贝塞尔光束特性

n=1的一阶贝塞尔光束具有以下特征：

近场相位呈螺旋形分布（0到2π变化）
远场强度分布为暗中心环形
携带轨道角动量，可用于光学微操控

通过改变拓扑荷数n，可以生成不同阶数的贝塞尔光束。高阶光束的环形半径随n增加而增大。

5. 宽带特性验证

由于几何相位与波长无关，同一超表面可以在不同波长下工作。验证宽带特性的方法：

保持超表面结构不变
改变光源波长（如450nm, 532nm, 635nm）
比较各波长下的远场强度分布

需要注意的是，虽然相位关系与波长无关，但纳米柱的尺寸效率会随波长变化。因此，要实现最优性能，可能需要针对不同波长优化纳米柱尺寸。

6. 实际应用中的注意事项

制造公差影响：
- 纳米柱尺寸误差应控制在±10nm以内
- 转角误差应小于5°
- 边缘粗糙度会影响散射损耗
材料选择考量：
- TiO₂在可见光波段透明且高折射
- 紫外波段可考虑Al₂O₃或HfO₂
- 近红外可用Si或GaN
性能优化方向：
- 采用渐变周期设计提高效率
- 优化纳米柱形状（如椭圆形）减少偏振相关性
- 多层结构扩展工作带宽
常见问题排查：
- 如果远场图案不对称：检查纳米柱排列和转角数据
- 如果效率过低：验证纳米柱高度和材料折射率
- 如果无衍射距离短：增加超表面尺寸或降低NA

7. 扩展应用与进阶设计

基于相同的设计框架，还可以实现以下功能：

多焦点贝塞尔光束：在相位分布中叠加多个径向相位项
非对称贝塞尔光束：引入倾斜相位项
偏振复用：设计对不同偏振态产生不同阶数光束
动态调控：结合液晶等可调材料实现可重构超表面

一个进阶示例是产生携带不同轨道角动量的多路贝塞尔光束，其相位分布可表示为：

φ(x,y) = mod(∑[aₙφₙ(x,y)], 2π)

其中aₙ是加权系数，φₙ是对应不同n值的相位分布。

在实际操作中，我发现超表面边缘的纳米柱对整体性能影响较大。通过逐渐减小边缘区域的单元密度（"切趾"技术），可以有效抑制旁瓣，提高光束质量。此外，对于高NA设计，需要考虑矢量衍射效应，简单的标量相位近似可能不够精确。