COMSOL仿真弯月型结构BIC：几何建模与斜入射设置

1. 弯月型结构BIC仿真概述

在光子学研究中，弯月型结构因其独特的几何特性常被用于实现连续域束缚态（Bound states in the continuum, BIC）。这种特殊的光学模式能将电磁能量局域在特定区域，同时保持与辐射连续态的耦合。通过COMSOL Multiphysics进行仿真时，斜入射条件下的线偏振设置尤为关键，直接影响BIC模式的激发效率。

弯月型结构本质上是一种非对称谐振腔，其几何参数（如圆弧半径、圆心偏移量）决定了电磁场的局域特性。当入射光满足特定条件时，结构会支持高品质因子（Q值）的共振模式，这正是BIC的典型特征。理解这一点对后续的仿真设置至关重要。

2. 几何建模与参数化设计

2.1 弯月型结构创建

在COMSOL中构建弯月型结构，最有效的方法是使用参数化曲线。通过两个圆弧的布尔运算（差集）可以精确控制月牙形状的几何特征。以下是详细的建模步骤：

1. 创建第一个圆弧（arc1）：

   • 圆心坐标设为(0,0)
   • 起始角度0°，终止角度180°
   • 半径R1为主要几何参数

2. 创建第二个圆弧（arc2）：

   • 圆心坐标设为(d,0)，其中d为圆心偏移量
   • 起始角度180°，终止角度0°
   • 半径R2通常与R1相同或略有差异

3. 执行布尔差集运算：

   matlab复制% COMSOL LiveLink脚本示例
   arc1 = mpharc(model, 'arc1', 'base', 'circle', 'point1', [0,0], 'angle1', 0, 'angle2', 180, 'radius', R1);
   arc2 = mpharc(model, 'arc2', 'base', 'circle', 'point1', [d,0], 'angle1', 180, 'angle2', 0, 'radius', R2);
   diff = mphgeom(model, 'diff1', 'setdiff', {'arc1','arc2'});
   

关键提示：圆心偏移量d与半径比(d/R)决定了结构的非对称程度，这直接影响BIC的产生条件。建议初始设置d/R≈0.2-0.3，后续通过参数扫描优化。

2.2 几何参数优化策略

在实际操作中，建议采用以下参数扫描方法：

1. 固定R1=1μm（参考值），扫描d从0.1到0.5μm
2. 记录每个d值对应的共振频率和Q值
3. 特别注意Q值突然升高的参数区域，这可能是BIC出现的信号

通过这种方法可以系统性地研究几何参数对BIC模式的影响，避免盲目尝试。记得保存每个参数组合的场分布图，便于后续分析模式特性。

3. 斜入射条件设置详解

3.1 波矢量分解原理

斜入射与正入射的最大区别在于波矢量的空间分量分布。设入射角为θ，波数k0=2π/λ，则波矢量的各分量为：

• x分量：kx = k0·sinθ
• z分量：kz = k0·cosθ
• y分量：ky = 0（二维情况）

在COMSOL中实现时，需要特别注意软件默认的坐标系定义。以下是具体设置步骤：

1. 定义入射角度变量：

   matlab复制theta = 30; % 入射角度(度)
   theta_rad = theta*pi/180; % 转换为弧度
   

2. 计算波矢量分量：

   matlab复制lambda = 1.55e-6; % 工作波长(m)
   k0 = 2*pi/lambda;
   kx = k0*sin(theta_rad);
   kz = k0*cos(theta_rad);
   

3. 在散射边界条件中正确输入这些分量值

3.2 坐标系投影要点

COMSOL的默认坐标系以z轴为主要传播方向，斜入射时需要特别注意：

1. 确保计算域足够大，使斜入射波能完整通过
2. 在端口设置中选择"用户定义"波矢量
3. 输入计算得到的kx和kz值
4. 检查波前方向是否与预期一致

常见错误是忘记考虑投影效应，导致实际入射角度与设定不符。建议先用均匀介质验证波矢量设置是否正确，再应用到复杂结构中。

4. 偏振状态精确控制

4.1 线偏振方向设置

对于斜入射情况，偏振方向的定义需要格外小心。COMSOL中偏振方向是相对于波矢量的局部坐标系定义的：

1. TM偏振（p偏振）：

   • 电场平行于入射平面
   • 设置方法：

     matlab复制polarization_TM = [cos(theta_rad), 0, -sin(theta_rad)];
     

2. TE偏振（s偏振）：

   • 电场垂直于入射平面
   • 设置方法：

     matlab复制polarization_TE = [0, 1, 0]; % y方向
     

4.2 偏振纯度验证技巧

在实际仿真中，可能会遇到偏振不纯的问题，影响BIC的观测。以下是验证和解决方法：

1. 添加场监视器检查各电场分量
2. 在均匀介质区域验证偏振方向是否正确
3. 如果发现偏振退化，检查：
   • 边界条件是否对称
   • 网格是否足够精细
   • 端口设置是否正确

特别提醒：弯月型结构的非对称性会导致偏振转换效应，这是正常现象，但需要与设置错误区分开来。

5. 边界条件与网格划分

5.1 散射边界条件设置

散射边界条件(Scattering Boundary Condition, SBC)用于模拟开放空间，防止非物理反射：

1. 选择所有外部边界
2. 应用散射边界条件
3. 输入之前计算的波矢量分量
4. 设置适当的阻尼系数（通常0.1-0.5）

重要提示：在斜入射情况下，必须正确设置入射波方向，否则会导致虚假反射。建议先用平面波测试边界条件的吸收效果。

5.2 网格划分策略

弯月型结构的网格划分直接影响结果精度，特别是BIC这种对场分布敏感的模式：

具体操作步骤：

1. 对弯月形边缘添加边界层网格
2. 设置尖端区域网格加密
3. 使用曲率自适应网格细化
4. 整体最大单元尺寸控制在λ/8以下

常见错误是只在全局设置细网格，而忽略了局部加密，导致计算资源浪费且尖端分辨率不足。

6. 后处理与BIC识别

6.1 远场散射分析

BIC的一个重要特征是远场散射的对称性破缺。设置监视器的关键步骤：

1. 添加远场计算域
2. 定义散射截面计算
3. 设置角度分辨率为1°或更小
4. 特别注意90°和270°方向的散射差异

6.2 近场特征识别

真正的BIC模式在近场会表现出独特分布：

1. 创建电场幅值切片图
2. 观察涡旋状场分布
3. 检查相位变化（应呈现拓扑特征）
4. 寻找类似"太极图"的场结构

典型BIC近场特征包括：

• 明显的场局域化
• 清晰的相位涡旋
• 高Q值共振（>10^4）

6.3 Q值计算方法

品质因子Q是识别BIC的关键指标：

1. 进行频率扫描（范围要足够窄）
2. 记录共振峰处的频率f0和半高宽Δf
3. 计算Q = f0/Δf
4. 寻找Q值突变的频率点

建议使用对数坐标绘制Q值曲线，便于观察数量级变化。

7. 常见问题排查指南

7.1 仿真结果异常排查表

7.2 高级调试技巧

1. 使用场探针实时监测关键点场值
2. 分阶段验证：先简后繁
3. 保存中间结果，便于回溯
4. 尝试不同的求解器设置

对于特别复杂的情况，建议：

1. 先从正入射开始调试
2. 使用更简单的几何验证设置
3. 逐步增加复杂度

8. 参数优化与性能提升

8.1 自动化参数扫描

利用COMSOL的批处理功能实现高效优化：

1. 创建参数化扫描研究
2. 定义关键变量范围（d, θ等）
3. 设置并行计算选项
4. 使用MATLAB LiveLink进行后处理

示例脚本框架：

matlab复制for d = 0.1:0.05:0.3
    model.param.set('d', num2str(d));
    model.study('std1').run;
    % 提取Q值等结果
end

8.2 计算资源优化

大型仿真可以采取以下加速策略：

1. 使用对称性简化模型（如可行）
2. 采用扫频法替代瞬态仿真
3. 合理使用自适应网格
4. 利用高性能计算集群

对于频域求解器，建议：

1. 选择合适的频段范围
2. 使用稀疏矩阵求解器
3. 合理设置迭代容差

9. 实际案例演示

9.1 典型参数设置参考

以下是一组经过验证的参数组合，可作为起点：

9.2 预期结果分析

成功仿真应观察到：

1. 共振频率约193THz附近
2. Q值>10^4的尖锐峰
3. 近场涡旋分布
4. 远场不对称散射

建议保存以下关键结果图：

1. 电场幅值分布（切片）
2. 相位分布
3. 散射方向图
4. Q值随参数变化曲线

10. 进阶技巧与扩展应用

10.1 多物理场耦合考虑

在实际器件中，可能需要考虑：

1. 热效应（高Q导致局部加热）
2. 机械应力（光力效应）
3. 材料非线性（高功率时）

这些效应可以通过COMSOL的多物理场接口实现耦合分析。

10.2 实验验证准备

仿真结果指导实验时需注意：

1. 考虑制造公差的影响
2. 添加适当的损耗机制
3. 模拟实际测量条件
4. 准备参数容差分析

建议进行蒙特卡洛分析，评估工艺波动对性能的影响。