当前位置：首页 > news >正文

COMSOL光子晶体中的角态与边界态

news 2026/5/12 11:19:23

comsol光子晶体中的角态与边界态。

当光子晶体开始玩“角落杀”

光子晶体这玩意儿，搞光的人应该不陌生。但最近发现它还能在角落里搞事情——比如“角态”（Corner States）和“边界态”（Edge States）。简单来说，角态就是光被锁在晶格的角落蹦迪，而边界态则是沿着两种不同结构的交界处溜达。听起来像不像物理版的密室逃脱？今天咱们用COMSOL来扒一扒这些态是怎么形成的，顺便搞点代码实操。

先搭个蜂窝棋盘

想玩角态，得先有个周期性的晶格结构。六方晶格（蜂窝结构）是经典选择。COMSOL里可以用脚本快速建模，比如用参数化几何画蜂窝晶格：

% COMSOL MATLAB API示例：生成六方晶格 model = ModelUtil.create('Model'); geom = model.geom.create('geom', 3); cyl = geom.create('cyl', 'Cylinder'); cyl.set('pos', [0 0 0]); cyl.set('r', 'r0'); % 柱体半径 cyl.set('h', 'h0'); % 高度 pattern = geom.create('pattern', 'Array'); pattern.set('type', 'hexagonal'); pattern.set('spacing', 'a'); % 晶格常数 pattern.set('fullsize', [10 10 0]); % 扩展10x10阵列 geom.runAll;

这段代码干了啥？先建了个圆柱体（模拟介质柱），然后用六方阵列无限复制，最后生成蜂窝状光子晶体。参数a控制晶格间距，r0决定柱子粗细。调这两个值，能改变光子带隙的位置——这直接关系到角态能不能出现。

边界条件：给光子上锁

comsol光子晶体中的角态与边界态。

想让光在边界或角落憋着，得设置合适的边界条件。比如在晶格边缘用Floquet周期边界模拟无限大结构，而在角落区域用散射边界条件或者完美电导体（PEC）。这里有个坑：如果边界条件设反了，光直接跑路，角态就没了。

% 设置Floquet边界 model.physics('emw').prop('BoundaryCondition').set('Floquet', 'on'); model.physics('emw').prop('FloquetX').set('kx', 'k0'); % Bloch波矢 model.physics('emw').prop('FloquetY').set('ky', '0');

Bloch波矢k0对应动量空间中的位置，扫频时从0到π/a遍历布里渊区。这时候你可能会想：为什么边界态只在特定k0出现？因为能带结构在这里劈叉了（带隙开了），光只能在带隙边缘的动量区域存活。

扫频：找光的藏身之处

跑完仿真，重点看频带图。角态通常出现在带隙中间，而边界态可能出现在带隙边缘。COMSOL的后处理功能可以提取本征频率和场分布：

% 提取本征频率 study = model.study('eig'); study.set('paramvalues', {'k0'}, linspace(0, pi/a, 50)); % 扫k0 solutions = study.getSolverSequences(); freqs = solutions.getEigenfrequencies(); % 绘制频带图 plot(real(k0_values), real(freqs), 'LineWidth', 1.5); xlabel('Bloch wavevector k'); ylabel('Frequency (Hz)');

如果看到带隙中有孤立的频率点（比如0.42c/a附近），恭喜，角态出现了！这时候切到场分布图，应该能看到光场挤在角落瑟瑟发抖。

翻车预警：参数别乱调