COMSOL光学模型中的回音壁腔PM强耦合模型（涉及耦合强度与反交叉曲线）

COMSOL光学模型：回音壁腔PM强耦合模型（耦合强度，反交叉曲线）

打开COMSOL的瞬间，电脑风扇开始狂转——这玩意儿吃硬件是真不客气。今天要折腾的是回音壁光学微腔的PM强耦合模型，说白了就是看两个光学模式怎么在特定条件下互相"撕扯"。先别急着点求解，咱们得把参数调对路子。

耦合间距这个参数贼关键，手动画网格太蠢，直接上参数扫描：

for d = linspace(100e-9, 500e-9, 20) model.param.set('gap', num2str(d)); model.study('std1').run; freq(:,k) = mphglobal(model, 'ewfd.freq'); k = k+1; end

这段脚本让间距从100纳米扫到500纳米，每次跑完把共振频率抠出来塞数组里。注意别把纳米写成米，否则你会看到频率曲线直接表演自由落体。

COMSOL光学模型：回音壁腔PM强耦合模型（耦合强度，反交叉曲线）

跑完数据别急着收工，重点在反交叉曲线的特征。上python做个暴力可视化：

plt.scatter(gap_list, freq[0,:], c='b', label='TE模式') plt.scatter(gap_list, freq[1,:], c='r', label='TM模式') plt.xlabel('耦合间距 (nm)') plt.ylabel('频率 (THz)') plt.legend()

这时候图上应该出现明显的X形交叉——如果没出现，八成是边界条件设成PEC了，赶紧回去检查材料设置。正常情况你会看到两个模式在某个临界间距突然互相回避，就像地铁里陌生人刻意保持的距离。

耦合强度计算不能照搬教科书公式，实测用洛伦兹拟合更靠谱：

[fitresult, gof] = fit( xData, yData, 'lorentz1' ); FWHM = fitresult.a1/(pi*fitresult.b1); g = abs(freq1 - freq2)/2; kappa = FWHM; C = g^2/(kappa^2); // 耦合强度指标

这里有个坑：当间距小于200纳米时，模式劈裂会突然增大，这时候得改用双洛伦兹拟合。曾经有哥们在这卡了三星期，最后发现是Q值设太高导致数值发散。

调通模型后记得存个.mph备份，COMSOL崩溃起来比女朋友翻脸还快。最后在论文里放图时，记得把颜色调成色盲友好系——审稿人指不定是红绿色弱呢。