2.2 能带结构计算与分析

2.2 能带结构计算与分析

本节摘要：能带图是光子晶体设计的通用语言：横轴是布里渊区高对称路径，纵轴是归一化频率，每条曲线是一支布洛赫模式；带隙是所有方向曲线上方的"空白带"。本节训练能带图的三步读法，并建立平面波展开、时域有限差分、有限元三大数值方法的选型直觉——频域算带隙用前者，时域看传输与腔模用后两者。

能带图怎么读

上一节把求解问题折叠成了"沿高对称路径逐 k 解本征值"，本节学习它的输出格式。读能带图只需要三步，但每一步都有容易翻车的地方。第一步看横轴：路径上的每个刻度是一个高对称点（区中心、边中点、角点），图上每一个水平位置对应一个特定的布洛赫波矢；横轴不均匀是正常的，各段长度只代表路径走向，不代表真实距离比例。第二步看纵轴与曲线：纵轴是归一化频率（常用无量纲的频率乘周期再除光速，即以晶格常数为尺度的"光程"度量）；每条曲线是一支模式，k 处曲线的高度就是该波矢允许的频率。曲线的斜率有直接物理意义——正比于群速度，平坦的带意味着慢光，这一点在第 5 章器件章会反复用到。第三步找空白：把路径上所有 k 的曲线整体看一遍，如果存在一个频率区间在整条路径上都没有任何曲线穿过，这就是沿该路径的带隙；只有当补充计算确认区间外所有方向的 k 都无模式时，才能宣布完全带隙。

图：一张能带图的解剖：三步读法标注

归一化频率：为什么单位消失了

能带图的纵轴几乎从不用真实频率单位，而是用归一化量——频率乘晶格常数再除以真空光速，或者等价地写成波长与晶格常数之比的倒数。这么做的理由上一节已经埋下：缩放律。同一张能带图，晶格常数取半微米就是通信波段器件，放大到毫米级就是太赫兹器件。归一化让"设计"与"定制"彻底解耦：设计阶段只关心归一化频率上的带隙位置与宽度，落地阶段选一个晶格常数把带隙中心平移到目标波长即可。反过来看，这也会让初学者犯一个典型错误——看到文献里带隙在归一化频率 0.3 附近，就直接拿真实波长去比。先换算：归一化频率乘以晶格常数除以光速，才得到真实频率。

三大数值方法：各自的地盘

能带计算的方法论目前三分天下，选型的关键问题是"你要答案的哪个侧面"。

平面波展开法是频域方法，直接继承上一节的本征值框架：把磁场展开成平面波基组，介电函数展开成傅里叶级数，本征方程化成矩阵对角化。它的优点是概念干净、与布洛赫定理严丝合缝，输出就是标准能带图；缺点是截断收敛慢（尤其介电对比度尖锐时），且难以处理有限尺寸与缺陷结构。一句话：算无限周期结构的带隙，它是默认选择。

时域有限差分法是时域方法，把空间切成网格、时间离散步进，直接模拟电磁波的传播。它不算本征值，而是"放一炮看波怎么走"：输入一个脉冲，接收端记录时间信号，做傅里叶变换得到透射谱——带隙直接表现为透射谱的深谷。它天生擅长有限结构、缺陷、非线性与复杂材料，能一次算出宽频段响应；代价是网格色散误差与长时间仿真的计算量。验证器件的实际传输行为，它是主力。

有限元方法是空间离散的频域方法，用自适应网格贴合复杂几何，在处理弯曲边界、倾斜侧壁、异形基元时精度优势明显，工业软件多以此为核心。缺点是网格生成复杂、大体积计算内存吃紧。几何复杂、精度苛刻的单个器件（如高品质因子微腔）常交给它。

三者的关系不是竞争而是接力：平面波展开法快速筛出"有带隙的候选结构"，有限元精细核对候选结构的真实模式与品质因子，时域有限差分做整机传输验证。下一节的仿真复盘走的就是第一棒。

实践中的三个陷阱

第一，收敛性检验不可省。平面波展开法的截断数不足时，能带位置系统性偏低、带隙宽度虚高——文献里两条"同结构不同结果"的能带曲线，多半是截断数不同。惯例是把截断数翻倍重算一遍，能带移动小于目标精度才收工。第二，路径覆盖不全。高对称路径只保证极值大多落在路径上，"大多"不等于"全部"；宣布完全带隙前，严谨的做法是在整个布里渊区做网格扫描，确认空档对区面上每一个 k 都成立。第三，有效介电函数的处理方式影响精度。平面波展开法需要对介电函数或其倒数做傅里叶展开，两种做法收敛速度差别明显；高对比度结构上，对倒数做展开（即以介电函数倒数插值）通常显著更快。这些坑在第 4 节的仿真复盘里都会实际踩一遍。

读图实战：三道自测题

读图能力要用考题验收。以下三道题都来自真实场景，建议先自行判断再看答案。题一：一张二维光子晶体能带图上，某支曲线在整个布里渊区内几乎水平，它对器件设计是好事还是坏事？答案：看用途——是慢光的原料（群速度低、光与物质作用增强，做延迟线与非线性增强是好事），也是色散的祸根（对结构涨落极度敏感、带宽极窄，做传输波导是坏事）；脱离用途谈曲线形状没有意义。题二：文献报告"带隙相对宽度百分之三十二"，你复算只有百分之二十四，最可能的原因是什么？答案：先核对带隙宽度百分比的定义（宽度除以中心频率，有作者用上边界做分母，数字会偏大）；再核对计算口径（是否含完全性判定、截断数是否一致）；最后才怀疑结构参数抄录误差——先怀疑口径、再怀疑参数，是复算失败的标准排查顺序。题三：一张能带图画出了六条带，第二与第三带之间在整条路径上有空档，作者宣称"完全带隙"，你接受吗？答案：不接受——路径空档只证明高对称路径上无模式，完全带隙需要对布里渊区全空间扫描确认；同时还要确认偏振声明（该空档是否只对单一偏振成立）。三道题背后的同一句话：能带图是结论的载体而非结论本身，读图永远要追问口径与适用范围。

常见问题

问：能带图上的带隙宽度百分比怎么算？
带隙宽度除以带隙中心频率，乘以百分之百。注意分母是中心而非上边界，读文献数据时留意作者口径。

问：二维计算的结果能用于三维平板吗？
不能直接用。二维平板的"有效折射率"近似会移动带隙位置，平板修正后的带隙通常比纯二维结果窄。二维结果当定向参考，平板设计必须按三维（或等效二维加修正）重算。

三种方法的机时预算

选型之外补一组机时数量级，帮助建立计划感。同一个小体积二维原胞问题：平面波展开秒级到分钟级，时域有限差分单次分钟级但一次覆盖宽谱，有限元分钟级。换成含缺陷的超胞或三维腔：平面波展开进入小时到天级，时域有限差分小时级（内存成为瓶颈），有限元小时级但网格自适应能省一半以上。经验结论：设计迭代的前期用便宜方法大量试错，后期才动用昂贵方法做终审——把有限元用在粗扫阶段是最常见的机时浪费，把平面波展开用在最终验收是最常见的精度浪费。

本节要点回顾

• 三步读法：看横轴（高对称点）、看曲线（斜率即群速度）、找空白（路径共同空档为带隙，全方向确认才是完全带隙）。
• 归一化频率：设计只对归一化频率负责，落地时选晶格常数做缩放；文献数据先换算再比较。
• 方法选型口诀：无限周期算带隙用平面波展开，复杂几何精算用有限元，器件传输验证用时域有限差分；三者接力而非互斥。
• 三大陷阱：截断数不收敛、高对称路径覆盖不全、介电函数插值方式选错；宣布带隙前逐一排查。

读图的能力就位。但光会读别人的图还不够——下一节我们拿起纸笔，亲手把最简单的一维带隙从方程推到公式。