Comsol超表面折射率传感器。 电磁诱导透明EIT和典型连续体中的束缚态BIC。

超表面这玩意儿最近在光学传感领域火得不行，特别是玩折射率检测的实验室，谁家要没摆弄过几个超表面结构都不好意思说自己是搞光子器件的。今天咱们就扒一扒COMSOL里怎么折腾电磁诱导透明（EIT）和束缚态（BIC）这两个物理现象来做高灵敏度传感器。

先看个典型的超表面单元结构代码：

% 超表面单元参数
period = 400e-9;  % 周期400纳米
bar_length = 300e-9;
bar_width = 80e-9;
gap = 30e-9;       % 纳米级间隙

这个参数设置有个坑要注意——当gap小于20nm时，实际加工会出现量子隧穿效应，仿真结果和实测会差出十万八千里。建议先用30nm起步，等摸清规律再挑战极限。

在COMSOL里建模时，边界条件设置直接决定仿真成败。下面这段代码锁定了电磁场边界：

physics.create("emw", "ElectromagneticWaves", "geom1");
physics.feature("pc1").set("Physicssurface", "PerfectConductor");
physics.feature("pc1").selection.named("metalsurface");

这里把金属结构设为完美导体简化计算，但实际做损耗分析时需要改用Drude模型，特别是当研究EIT效应时，金属损耗会显著影响透明窗口的Q值。

当我们在频域扫描时，突然在透射谱里看到尖锐的透明峰——恭喜，EIT现象出现了。这时候赶紧运行参数化扫描：

study.create("param", "Parametric");
study.feature("param").set("plist", "linspace(1.33,1.38,10)");
study.feature("param").set("pname", "n_analyte");

通过改变环境折射率n_analyte，观察透明峰偏移量。有个反直觉的现象：当分析物折射率增大时，某些结构的共振峰会红移，但EIT窗口反而蓝移，这跟模式耦合强度有关。

BIC的情况更魔幻。原本应该无限Q值的理论束缚态，在实际仿真中会因为结构不对称出现漏光。我们通过在COMSOL里微调对称性参数来找准BIC点：

asymmetry = 0.01;  % 1%的结构不对称
theta = 0:0.1:10;   % 入射角扫描

当theta超过3度时，原本的BIC状态突然崩解，Q值断崖式下跌。这种非线性响应对应着超高的折射率灵敏度，实测能达到10^5 nm/RIU量级，比传统表面等离子体共振传感器高两个数量级。

最后给个实战技巧：在优化结构参数时，别直接用COMSOL内置的优化模块，先跑个粗糙的全局扫描锁定大致范围。比如用蒙特卡洛方法随机撒点：

import numpy as np
param_space = {
    'length': np.random.uniform(200e-9, 400e-9, 50),
    'width': np.random.uniform(50e-9, 150e-9, 50)
}

抓出透射谱线型最优的五个样本，再上梯度下降法精细调参。这么玩比全程自动优化省时得多，特别适合实验室那帮用着老旧工作站的苦逼研究生们。

搞超表面传感器这行的都知道，真正的核心技术往往藏在看似普通的参数设置里。就像那个把灵敏度做到2500nm/RIU的组，核心机密不过是把纳米棒阵列旋转了15度——物理有时候就是这么不讲道理又充满惊喜。