comsol二维光子晶体角态。 单胞能带，超胞能带，边界态以及角态计算。

二维光子晶体角态的研究就像玩拓扑积木——把周期性结构掰碎再重组总能发现神奇的光学现象。今天咱们用COMSOL实战演练从单胞到超胞的拓扑演变，手把手带你看边界态如何拐个弯就变成角态。

先打开COMSOL新建二维模型，选波动光学模块。画个六边形晶格当单胞，设置a=500nm的晶格常数。材料参数得讲究：

epsilon_air = 1.0^2;
epsilon_si = 3.5^2;  //硅介质柱
radius = 0.3*a;      //占空比控制

这里的关键是Floquet周期性边界条件的设置，k·p参数化扫描要覆盖整个布里渊区边界。扫频时注意设置足够密的采样点，否则狄拉克锥可能显示不全。计算完单胞能带后，重点检查Γ点附近的简并态——这玩意儿就是后续拓扑操作的种子。

当把单胞扩展成3×3超胞时，边界条件要改为开放边界。这时候需要微调模型尺寸：

supercell_size = 3;
model.component("comp1").geom("geom1").scale([supercell_size, supercell_size]);

超胞能带计算会蹦出新的边缘态，但别急着截图，先在结果里加个电场模分布图。如果看到光场像被502胶水粘在边界线附近，恭喜你捕获到边界态了。

角态的玄机藏在晶格转角里。试试把超胞顺时针旋转15度：

rotation_angle = 15;
model.component("comp1").geom("geom1").rotate(rotation_angle);

这时候边界态会发生拓扑相变，电场开始往拐角处聚集。为了验证这是真·角态而非普通局域模，可以做局域态密度(LDOS)分析：

study = model.study("std1");
study.feature("freq").set("plist", "linspace(0.45,0.55,50)");
model.result().dataset("dset1").run;

当LDOS在特定频率出现尖峰，且场分布呈现点状局域，基本实锤角态的存在。最后用参数扫描批量跑不同转角的数据，对比能带折叠情况，整套拓扑演化链条就完整了。

搞光子晶体的都知道，这种角态对微纳传感特别有用——光场被锁在拐角处时，周围介质折射率的微小变化会导致共振峰剧烈偏移。下次如果有人问怎么做高灵敏度光学传感器，直接把角态电场分布图甩他脸上准没错。