COMSOL光学模型：魔角光子晶体激光器（能带，模式，参数化几何建模，激光，图中只放了一个例子）
在魔角光子晶体激光器的设计里，参数化建模是逃不开的套路。每次打开COMSOL，总得先搞明白周期性结构怎么摆弄——比如六边形蜂窝晶格，用参数控制晶格常数和空气孔半径，比手动拖拽坐标强多了。这里有个小技巧：用Java方法批量生成空气孔阵列，代码虽然看着像外星符号，但能省下两小时重复劳动。
// 参数化生成空气孔 for (int i = -5; i <= 5; i++) { for (int j = -5; j <= 5; j++) { double x = i*a + (j%2)*0.5*a; double y = j*a*Math.sqrt(3)/2; model.geom().create("circ"+i+j, "Circle").set("r", r); model.geom("circ"+i+j).set("pos", new double[]{x, y}); } }这段代码其实在玩空间对称性的魔术。当晶格常数a=300nm，空气孔半径r=0.3a时，系统会突然进入光子带隙状态——就像乐高积木找到了完美卡扣角度。这时候在能带图上会看到明显的平带区域，对应着光场被锁在特定位置的魔角效应。
仿真模式场的时候总得和边界条件较劲。PML层厚度设成工作波长的1.5倍是玄学经验值，但有时候偷懒减到1.2倍反而收敛更快。观察TE模的电场分布时，会发现光强在空气孔边缘形成环形驻波，这玩意儿就是激光的种子选手。调Q因子时要注意材料损耗的设置，硅的消光系数给到1e-4刚好能让模式线宽收窄到0.1nm量级。
激光阈值计算是场持久战。得在频域求解里嵌套自定义方程，把增益介质模型写成随粒子数反转变化的复折射率。有一次我把载流子寿命参数多打了个零，结果阈值电流直接跌到物理定律禁止的区域——这种数值陷阱在参数扫描时经常遇到，得像扫雷游戏一样小心。
最后放张能带图镇楼（假装这里有图）：横轴是布里渊区高对称点路径，纵轴归一化频率。那条突兀的水平线就是魔角对应的平带，模式局域在七个空气孔构成的超胞里。这种结构最骚的操作在于，稍微旋转晶格5度，Q值就能飙升三个数量级——所以说参数化建模不仅是偷懒工具，更是探索物理奇点的钥匙。
