COMSOL六边形光子晶体能带模型

先看晶格搭建。六边形的精髓在于空间周期性排列，这里用参数化构建更灵活：

double a = 1e-6; // 晶格常数1微米
double r = 0.4*a; // 介质柱半径
// 生成六边形基准单元
model.geom.create("hex", "Hexagon");
model.geom("hex").set("type", "regular");
model.geom("hex").set("size", a);

这段代码在几何中创建了基础六边形，size参数控制蜂窝结构的间距。注意这里的regular类型会自动生成正六边形，比手动画坐标省事得多。不过有个隐藏细节——COMSOL的正六边形默认边长为size参数除以√3，需要自己验证尺寸是否准确。

介质柱的排列采用晶格副本生成：

% 周期性阵列扩展
model.component("comp1").geom.create("array1", "Array");
model.component("comp1").geom("array1").set("size", {'3*a' '3*a' '0'});
model.component("comp1").geom("array1").set("displ", {'-a' '-a' '0'});

这里故意设置3x3阵列并位移-a，是为了消除边界效应。实际操作中发现，当扫描波矢到布里渊区边界时，太小的阵列会导致模式泄露，5x5可能更稳但计算量暴增，需要自己权衡。

重点来了——Floquet边界条件设置：

// 波矢扫描参数
double kx_max = pi/a;
double ky_points = 21;
// 创建参数化扫描
ModelUtil.study("std1").feature.create("param", "Parametric");
ModelUtil.study("std1").feature("param").set("pname", "kx");
ModelUtil.study("std1").feature("param").set("plist", range(0,kx_max/20,kx_max));

这里偷了个懒，只扫描Γ到M点方向。真正的行家会走完整个不可约布里渊区边界，但新手建议先搞定单方向。注意ky_points的数量别设太大，否则特征频率求解器可能报特征值丢失——尤其是当介质柱占空比过高时。

求解后的后处理藏着玄机：

# 提取特征频率
results = model.study("std1").getSolution()
freqs = results.getEigenvalues()
# 过滤物理模式
valid_bands = [f for f in freqs if f.imag < 1e-3]

这个过滤条件专治各种"幽灵模式"。COMSOL在计算时会吐出复数解，但物理模式其实应该接近纯实数。遇到过有人把虚部全当损耗分析，结果发了篇漏洞百出的论文，这事咱得避免。

最后秀张典型的TE极化能带图（想象此处有图）：低频区域出现明显的带隙，在0.3-0.4c/a频率范围内模式密度骤降。但要注意，这里假设介质柱是完美硅材料（ε=12），实际材料色散会让计算复杂度翻倍。

调试时常见翻车现场：1）网格太粗导致狄拉克锥消失 2）周期性边界条件方向设反 3）特征值搜索范围没覆盖实际频率。下次碰到能带曲线出现诡异震荡，先查这三项准没错。