Comsol经典小案例：铌酸锂LNOI体系中的法诺Fano共振

Comsol经典小案例。 铌酸锂LNOI体系法诺Fano共振。

法诺共振这玩意儿在光学微腔里特别有意思——它那种不对称的线型总让我联想到被狗啃过的面包。铌酸锂（LNOI）作为新一代光子学材料，凭借其超高的二阶非线性系数，在调制器、频率转换器里混得风生水起。今天咱们用Comsol折腾个LNOI微盘腔，看看怎么玩转这个奇葩的共振现象。

打开Comsol先搞几何建模。这里有个骚操作：用参数化曲线画微盘边界。代码区敲入：

theta = linspace(0,2*pi);
R = 10e-6;  % 半径10μm
ecc = 0.05; % 椭圆度
x = R*(1 + ecc*cos(3*theta)).*cos(theta);
y = R*(1 + ecc*cos(3*theta)).*sin(theta);

这可不是普通的圆盘，故意加了3θ的扰动项。椭圆度参数ecc控制着模式耦合强度，就像往平静的湖面扔石子，让原本简并的模式开始眉来眼去。

材料设置阶段要特别注意铌酸锂的各向异性。在材料库找到LiNbO3后，手动覆盖其折射率张量：

epsilon = [n_o^2 0 0;
           0 n_o^2 0;
           0 0 n_e^2];

这里no=2.211，ne=2.138可不是随便填的数值，它们决定了TE和TM模式的群速度差异。记得把z轴对准晶体的光轴方向，不然仿真结果会变成玄学。

Comsol经典小案例。 铌酸锂LNOI体系法诺Fano共振。

边界条件配置最容易被新手搞砸。这里要用到完美匹配层（PML），但得调整其拉伸函数：

pml.set('coord_system', 'cartesian');
pml.set('stretching', [1+1i, 1+1i, 1]);

复数拉伸参数相当于给电磁场加了缓刑，让它们在边界处自然死亡而不是暴力反射。要是看到结果里有奇怪的驻波，八成是这里的实部虚部比例没调好。

扫频计算时有个坑爹细节——必须用自适应频点扫描。在求解器设置里勾选"自适应细化"，系统会自动在共振峰附近加密采样。看这段日志输出：

发现Q因子>1e4，在246.7-247.2THz区间增加5个采样点

这功能简直救命，手动扫频的话要么错过尖峰，要么算到电脑冒烟。记得设置最大细化次数别超过5次，否则容易陷入无限循环。

最后可视化阶段，在结果里截取电场分布的动画帧。用MATLAB后处理时注意：

mphplot(model, 'pg2', 'frame', 'all', 'savepath', 'frames');
imwrite(imread('frames_01.png'), 'resonance.gif', 'DelayTime',0.1);

生成动态GIF能清晰看到模式跳变过程。当基模和高阶模开始勾肩搭背时，传输谱上就会冒出标志性的不对称峰——这货就是法诺共振本尊了。搞完这波操作，记得喝杯咖啡犒劳自己，毕竟和Comsol斗智斗勇堪比马拉松。