1. 红外10.6μm几何相位超表面透镜设计解析
在红外光学系统中，超表面透镜因其轻薄、可定制化等优势正逐渐取代传统折射透镜。最近我在设计一个工作波长为10.6μm的红外超表面透镜时，遇到一个看似简单但颇具代表性的问题：设计的透镜有效孔径为50μm，但生成的GDSII版图中纳米柱基元却超出了这个范围。这个问题引发了我对超表面透镜设计细节的深入思考。
1.1 超表面透镜的基本原理
超表面透镜的核心在于通过亚波长结构的几何相位调控来实现波前整形。与传统折射透镜不同，它不依赖光程差，而是利用纳米结构的旋转角度来引入相位延迟。具体到我们这个设计：
- 工作波长：10.6μm（中红外波段）
- 透镜类型：几何相位（Pancharatnam-Berry相位）超表面
- 相位调控方式：通过旋转矩形纳米柱引入交叉极化分量相位延迟
这种设计的关键参数包括纳米柱的尺寸（长轴w1=3μm，短轴w2=1μm）、高度h=8μm，以及排列周期P=4μm。这些参数的选择不是随意的，而是基于严格的电磁仿真和相位调控需求。
关键点：几何相位超表面的相位调控能力取决于纳米结构的双折射特性。长宽比w1/w2=3确保了足够的双折射，而高度h=8μm则保证了π的相位延迟（对于折射率n≈3.4的红外硅材料）。
1.2 设计中的边界问题
原始设计设定的有效孔径为50μm（半径25μm），但在实际生成的GDSII版图中发现最外圈的纳米柱超出了这个范围。这主要是因为：
- 基元排列方式：采用极坐标网格排列，最外圈基元中心位于r=25μm处
- 基元尺寸：矩形纳米柱尺寸为3×1μm，对角线长度约3.16μm
- 实际占用空间：最外圈基元旋转后，其角点会超出r=25μm的范围约1.58μm（对角线一半）
因此，实际需要的衬底尺寸至少应为50 + 2×1.58 ≈ 53μm。实际操作中发现调整到60μm更为稳妥，这给了工艺一定的容错空间。
2. 代码实现细节与优化
2.1 核心参数设置
python复制# 1. 核心参数设置（修正重叠关键：grid_step=P）
lamda = 10.6 # 工作波长(μm)，红外10.6μm
f = 100.0 # 透镜焦距(μm)
D = 50.0 # 透镜有效孔径(μm)，直径
k0 = 2 * np.pi / lamda # 真空波数(1/μm)
P = 4.0 # 基元周期(μm)，<λ/2（10.6/2=5.3），满足亚波长
w1 = 3.0 # 矩形纳米柱长轴(μm)，强各向异性
w2 = 1.0 # 矩形纳米柱短轴(μm)，w1/w2=3>1.5
h = 8.0 # 矩形纳米柱高度(μm)，红外Si(n≈3.4)满足相位延迟差≈π
grid_step = P # 修正：径向栅格步长=周期P（原P/2，过小导致小半径环基元过密）
r_max = D / 2 # 透镜最大径向半径(μm)
这些参数的选择都有其物理意义：
- 周期P必须小于λ/2以避免衍射级次（对于10.6μm，P<5.3μm）
- w1/w2>1.5确保足够的双折射效应
- 高度h通过仿真确定能提供π的相位延迟
2.2 相位分布计算
超表面透镜的相位分布遵循抛物面相位轮廓：
python复制# 理想抛物面相位分布（小NA<0.3，无球差），模2π归一化
phi_target = -k0 * (r ** 2) / (2 * f) # 核心相位公式
phi_target = phi_target % (2 * np.pi) # 归一化到0~2π
这个公式实现了将平面波转换为会聚球面波的功能。相位模2π归一化是因为几何相位超表面只能提供0-2π的相位调控。
2.3 旋转角分布求解
几何相位超表面的核心原理是旋转各向异性纳米结构来引入相位延迟：
python复制# 核心公式：α=φ_target/2，旋转角范围0~π
alpha = phi_target / 2 # 基元旋转角(rad)
这个关系式表明，相位延迟是旋转角的两倍。因此，要实现2π的相位覆盖，只需要0-π的旋转角范围。
3. GDSII版图生成的关键技术
3.1 基元排列策略
为了避免基元重叠，代码中采用了以下优化措施：
- 径向步长等于周期P（grid_step=P），确保径向间距足够
- 周向基元数量由周长自然决定：N_theta = int(2 * np.pi * ri / P)
- 设置最小基元数量为4，避免极稀疏排列
python复制# 遍历每个径向半径，生成环形栅格的基元阵列（无重叠优化）
for i in range(len(r)):
ri = r[i]
alphai = alpha[i]
if ri == 0: # 透镜中心，单个基元（无重叠）
main_cell.add(gp.CellReference(rect_cell, rotation=alphai * 180 / np.pi))
else: # 非中心，环形排布基元
N_theta = int(2 * np.pi * ri / P)
N_theta = max(N_theta, 4) # 最小4个基元
theta = np.linspace(0, 2 * np.pi, N_theta, endpoint=False)
x = ri * np.cos(theta)
y = ri * np.sin(theta)
for j in range(N_theta):
ref = gp.CellReference(
rect_cell,
rotation=alphai * 180 / np.pi, # 弧度转角度
origin=(x[j], y[j]) # 平移到目标坐标
)
main_cell.add(ref)
3.2 版图边界问题解决方案
针对最初提到的边界超出问题，有以下几种解决方案：
-
调整衬底尺寸（最简单直接）：
- 将CST中的衬底尺寸从50μm调整为60μm
- 优点：无需修改设计代码，操作简单
- 缺点：实际有效孔径仍为50μm，外围区域不参与成像
-
修改代码中的有效孔径：
python复制D = 50.0 - 2*1.58 # 补偿边界溢出- 优点：保持衬底尺寸一致性
- 缺点：实际有效孔径会略小于设计值
-
优化基元排列算法：
- 对外圈基元进行特殊处理，适当内缩
- 优点：精确控制有效孔径
- 缺点：算法复杂度增加
实际建议：对于科研用途，方案1最为实用；对于需要精确控制孔径的工业应用，建议采用方案3。
4. 设计验证与性能评估
4.1 相位分布验证
通过代码生成的相位分布图应该呈现完美的抛物面轮廓（模2π）。我们可以通过以下指标评估设计质量：
- 相位连续性：相邻基元间的相位差应平滑变化
- 旋转角分布：应严格遵循α=φ/2的关系
- 对称性：对于圆形透镜，相位分布应具有旋转对称性
4.2 版图质量检查
在将GDSII文件导入CST或其他仿真软件前，应该检查：
- 基元密度：单位面积内的基元数量是否合理
- 边界完整性：最外圈基元是否完整
- 重叠检查：确保没有任何基元重叠
- 工艺兼容性：最小特征尺寸（1μm）是否符合工艺能力
5. 常见问题与解决方案
5.1 基元重叠问题
现象：在较小半径处出现基元重叠
原因：径向步长设置过小或周向基元数量强制下限过大
解决方案：
python复制grid_step = P # 确保径向步长不小于周期
N_theta = max(int(2 * np.pi * ri / P), 4) # 最小4个基元，避免过大
5.2 边界溢出问题
现象：最外圈基元超出设计孔径
原因：未考虑基元旋转后的实际占用空间
解决方案：
- 按前文所述三种方案调整
- 推荐补偿公式：D_compensated = D - 2*(w1/2)*√2
5.3 相位不连续
现象：相位图中出现突变
原因：
- 模2π操作不当
- 旋转角计算错误
检查点：
python复制phi_target = phi_target % (2 * np.pi) # 确保正确归一化
alpha = phi_target / 2 # 旋转角应为相位的一半
5.4 仿真与实测差异
可能原因：
- 材料折射率与实际不符
- 纳米柱高度误差
- 边缘效应未考虑
改进措施：
- 进行参数扫描分析敏感性
- 添加工艺误差补偿
- 考虑添加边缘过渡区
6. 设计扩展与优化建议
6.1 多波长设计
当前设计针对单一波长10.6μm优化。要实现多波长工作，可以考虑：
- 色散工程：设计不同区域对应不同波长
- 叠加设计：将多个波长设计融合在一个超表面上
- 宽带设计：优化纳米结构实现宽带响应
6.2 效率优化
几何相位超表面的主要效率限制来自：
- 交叉极化转换效率
- 材料吸收损耗
- 制造误差
优化方向：
- 采用更高折射率材料
- 优化纳米柱形状（如椭圆形）
- 添加抗反射结构
6.3 工艺适配建议
根据不同的制造工艺，可能需要调整：
- 电子束光刻：适应最小特征尺寸
- 反应离子刻蚀：优化侧壁垂直度
- 纳米压印：考虑脱模角度
在实际流片前，建议：
- 与工艺工程师沟通设计规则
- 添加测试结构和对准标记
- 考虑工艺误差的补偿设计
通过这个项目，我深刻体会到超表面透镜设计中细节的重要性。一个看似简单的边界溢出问题，背后涉及的是设计理念与实际制造的衔接。这种经验是教科书上难以学到的，只有在实际项目中才会遇到并需要解决。