Comsol连续体中的束缚态BIC。 涉及能带计算与Q因子计算，包含一维光栅和二维光子晶体板。 注: 不包含拓扑荷计算。

在光子晶体领域，束缚态的光子晶体模式（Bound States in the Continuum, BICs）是一种特殊的光学模式，这些模式虽然存在于连续态的频率范围内，但仍然表现出极高的质量因子（Q factor），从而在光学集成电路中具有重要应用价值。本文将通过Comsol软件，分别对一维光栅和二维光子晶体板中的束缚态BIC进行建模与分析，涵盖能带计算和Q因子计算。

一、建模流程概述

在Comsol中建模BIC之前，需要明确几个关键点：

1. 选择合适的物理场：由于我们需要研究光子晶体的光学特性，选择“波的传播”模块，即“Electromagnetic Waves, Frequency Domain”。
2. 设置计算域：根据具体结构选择合适的几何形状，并确保计算域足够大以避免边界反射对结果的影响。
3. 材料和边界条件：定义材料的折射率分布，并设置适当的边界条件（如Perfectly Matched Layer, PML）以吸收边界处的反射波。

以下是Comsol中设置的基本代码结构：

% 创建一个新的模型
model = Model('Model');
% 添加波传播物理场
em = model.component('comp1').physics.create('emw');
% 设置计算域
em.create('dom', 'Domin');
% 设置材料
material_tags = {'Air', 'Dielectric'};
model.material.create('mat1').name('Air');
model.material('mat1').property('relpermittivity').set('1');
model.material('mat1').property('relpermeability').set('1');
model.material.create('mat2').name('Dielectric');
model.material('mat2').property('relpermittivity').set('n^2'); % 假设折射率为n
% 设置边界条件
em.create('bnd1', 'Boundary');
em.feature('bnd1').name('PML');
em.feature('bnd1').physics('emw').create('PML', 'PML');

二、一维光栅的能带计算

一维光栅是光子晶体的最简单形式，通常由周期性排列的介质柱或沟槽构成。通过计算其能带结构，我们可以找到潜在的BIC模式。

1. 几何结构：假设一维光栅由周期性排列的高折射率介质柱组成，周期为a，柱宽为d，介质折射率为n。
2. 能带计算：通过Comsol的“Eigenfrequency” study，计算不同波矢k点处的固有频率。

以下是计算能带时的代码片段：

% 设置扫描的k点（比如Γ到X点）
k_points = linspace(0, pi/a, 100);
% 初始化存储能带数据的数组
frequencies = zeros(size(k_points));

for i = 1:length(k_points)
    % 设置k矢
    em.feature('eig1').set('kx', k_points(i));
    % 求解
    model.solve();
    % 提取固有频率
    frequencies(i) = model.result('r1').feature('eig').result('freq').data(1);
end

通过上述计算，我们可以得到一维光栅的能带结构，通常在Γ点附近会出现明显的带隙，而在带隙内部可能存在BIC模式。

三、二维光子晶体板的能带计算

二维光子晶体板的计算相对复杂，但其设计自由度更高，能够实现更高的Q因子。

1. 几何结构：假设二维光子晶体板由平面上的周期性孔阵构成，周期为a，孔径为d，板厚度为t。
2. 能带计算：同样通过“Eigenfrequency” study，计算不同k点的固有频率。

以下是二维结构的代码片段：

% 设置二维模型
model.component('comp1').geom('geom1').create('r1', 'Rectangle');
model.component('comp1').geom('geom1').feature('r1').set('size', [a, a]);
model.component('comp1').geom('geom1').create('c1', 'Circle');
model.component('comp1').geom('geom1').feature('c1').set('radius', d/2);
model.component('comp1').geom('geom1').run('dif');
% 设置网格
model.component('comp1').mesh('mesh1').create('ftri');
model.component('comp1').mesh('mesh1').run();

通过计算二维光子晶体板的能带结构，我们可以找到BIC模式对应的频率点，并进一步计算其Q因子。

四、Q因子计算

Q因子是衡量BIC模式质量的重要指标，通常通过传输谱（transmission spectrum）计算。具体步骤如下：

1. 设置传输谱计算：将光子晶体板放置在波导中，通过计算波导输入端和输出端的传输率，得到传输谱。
2. 计算Q因子：通过传输谱的对称性和宽度估算Q因子。

以下是计算传输谱的代码片段：

% 设置波导
model.component('comp1').physics('emw').create('port1', 'Port');
model.component('comp1').physics('emw').feature('port1').set('location', 'left');
model.component('comp1').physics('emw').feature('port1').set('direction', 'normal');
% 设置传输谱计算
model.study('std1').feature('eig1').set('freq', 'freq');
model.study('std1').run();
% 计算传输率
T = model.result('r1').feature('port1').result('T').data';
Q = freq./ (fwhm/2); % 假设fwhm为峰宽

五、总结

通过Comsol软件，我们成功实现了对一维光栅和二维光子晶体板中BIC模式的建模与分析。能带计算帮助我们找到潜在的BIC模式，而Q因子计算则为我们提供了这些模式的质量评估。在实际应用中，可以通过优化几何参数（如周期、孔径等）来提高Q因子，从而实现更高性能的光学集成电路。未来，我们可以进一步研究这些BIC模式在光学传感、光学交换等领域的实际应用。

Comsol连续体中的束缚态BIC。 涉及能带计算与Q因子计算，包含一维光栅和二维光子晶体板。 注: 不包含拓扑荷计算。