1. 复合波导光栅中的准BIC与古斯汉森位移效应解析
在光子学器件设计中，复合波导光栅结构因其独特的电磁场调控能力而备受关注。其中，准束缚态连续体（quasi-Bound state in the Continuum, 准BIC）与古斯汉森（Goos-Hänchen, GH）位移的耦合效应，为光场操控提供了新的物理机制。这种效应本质上是通过精心设计的周期性结构，在特定波长和入射角度下，使反射光束产生远超常规情况的横向位移。
准BIC是一种特殊的共振态，其品质因数（Q值）理论上可以趋近于无限大。在实际结构中，由于工艺限制和材料损耗，我们通常只能实现"准"BIC状态。当光栅参数调整至准BIC状态时，电磁场会被强烈局域在波导层内，形成极高的能量密度。这种异常的场增强效应会显著改变光与物质的相互作用，进而影响反射光束的GH位移。
GH位移是光束在界面反射时产生的横向位移现象，传统情况下位移量通常在波长量级。但在准BIC条件下，由于场梯度剧烈变化，位移量可提升两个数量级以上。这种增强效应在光学传感、集成光子器件和光通信等领域具有重要应用价值。
2. COMSOL仿真模型构建详解
2.1 几何建模技巧
构建复合波导光栅模型时，采用分层设计思路能有效提高建模效率。基础结构由下至上依次为：
- 硅基底（折射率~3.5）
- 二氧化硅波导层（厚度200nm，折射率~1.45）
- 硅光栅层（厚度150nm）
使用参数化曲面生成光栅结构是提升建模效率的关键。在COMSOL中，FourierSeries特征可以快速生成锯齿波形：
matlab复制geom.feature('wp1').geom.create('f1', 'FourierSeries').set({
'0', 'h_sio2', 'period/2', 'h_si',
'0', 'period', '0', 'h_sio2',
'sawtooth', '2'});
注意：单位制必须统一使用纳米(nm)量级，COMSOL对单位一致性要求极为严格。建议在模型开头就通过
model.param.set定义所有参数变量，避免后续混淆。
2.2 材料参数设置
材料属性设置需要考虑色散效应，特别是在近红外波段（1500-1600nm）：
- 硅的折射率采用Sellmeier方程：
matlab复制n_Si = sqrt(1 + 10.6684293*λ^2/(λ^2-0.301516485^2) + 0.003043475*λ^2/(λ^2-1.13475115^2) + 1.54133408*λ^2/(λ^2-1104.0^2)) - 二氧化硅采用固定折射率1.45
- 背景介质设为空气（n=1.0）
3. 物理场设置与边界条件
3.1 电磁波频域求解器配置
选择"电磁波，频域"物理场接口，设置波动方程：
code复制∇×(μ_r^-1 ∇×E) - k0^2 ε_r E = 0
其中k0为自由空间波数，ε_r和μ_r分别为相对介电常数和磁导率。
边界条件设置要点：
- 上下边界设为完美匹配层(PML)，厚度设为半波长
- 左右边界设为周期性边界条件：
matlab复制physics.feature('prt1').set('PortType', 'PeriodicPort'); physics.feature('prt1').set('PeriodicityParameter', {'period', '0'}); - 端口激励采用TE偏振模式（电场垂直于入射面）
3.2 模式分析参数设置
传播常数β扫描是关键步骤，其与共振频率ω的关系为：
code复制β(ω) = n_eff(ω)·ω/c
其中n_eff为有效折射率，c为光速。在准BIC点附近，需要设置密集的扫描步长（建议Δβ<0.001μm^-1）。
4. 参数扫描与结果分析
4.1 双参数扫描策略
采用两阶段扫描方法提高效率：
- 粗扫阶段：硅层厚度100-200nm，步长20nm；二氧化硅层150-250nm，步长25nm
- 精扫阶段：在Q值峰值区域，步长缩小至5nm
扫描脚本示例：
matlab复制study.feature('param').set('plist', {
'100:20:200', // 硅层扫描范围
'150:25:250' // 二氧化硅层扫描范围
});
4.2 品质因数Q的计算
Q值通过共振线宽Δλ计算：
code复制Q = λ_resonance / Δλ
在COMSOL中可通过以下步骤获取：
- 在共振频率附近进行频域扫描
- 提取反射谱的共振峰
- 拟合洛伦兹线型获取半高宽
4.3 古斯汉森位移计算
GH位移量D_GH由场相位梯度决定：
code复制D_GH = - (λ/2π) · (dφ/dθ)
其中φ为反射相位，θ为入射角。在COMSOL中可通过以下步骤实现：
- 在端口边界设置场监视器
- 计算反射系数的相位变化
- 对入射角求数值微分
5. 关键结果可视化技巧
5.1 场分布图优化
准BIC状态下的场分布具有以下特征：
- 能量高度局域在波导层
- 电场强度增强10-100倍
- 明显的驻波模式形成
可视化建议：
- 使用对数尺度显示场强（log10|E|^2）
- 叠加等值线显示介质边界
- 添加箭头图表示能流方向
5.2 参数依赖关系图
制作二维参数扫描图时：
- X轴：硅层厚度（100-200nm）
- Y轴：二氧化硅层厚度（150-250nm）
- 颜色映射：Q值或GH位移量
- 添加等高线标识关键区域
6. 常见问题与解决方案
6.1 收敛性问题
现象：求解器报错或结果不收敛
解决方法：
- 检查网格尺寸是否足够小（建议λ/10）
- 调整PML厚度和位置
- 尝试不同的初始猜测值
6.2 假共振峰识别
现象：反射谱中出现非物理的尖峰
解决方法：
- 检查网格是否足够精细
- 验证材料参数是否正确
- 改变扫描步长重新计算
6.3 计算资源优化
大型参数扫描时：
- 使用批处理模式运行
- 采用分布式计算（如有条件）
- 先进行低精度计算定位关键区域
7. 实际应用中的设计考量
7.1 工艺容差分析
准BIC对结构参数极其敏感，实际制造时需要考虑：
- 刻蚀侧壁角度影响（80°-90°）
- 层厚均匀性（±5nm）
- 表面粗糙度（<2nm RMS）
7.2 温度稳定性
温度变化会导致：
- 硅的热光系数：1.86×10^-4/K
- 二氧化硅热膨胀系数：0.55×10^-6/K
建议工作温度控制在±0.1K以内
7.3 耦合效率优化
提高光栅耦合效率的方法：
- 优化光栅齿形（梯形vs矩形）
- 添加抗反射涂层
- 采用渐变周期设计
在完成整套仿真流程后，建议将关键参数（如峰值Q值、最大GH位移量）整理成表格，方便后续论文写作和实验对比。同时保存完整的模型文件和脚本，便于结果复现和参数调整。