2.3 带隙成因手推：从布拉格条件到一维平面波展开

2.3 带隙成因手推：从布拉格条件到一维平面波展开

本节摘要：本节在纸面上完成一次完整的带隙推导：对一维层状介电结构，从布拉格条件出发，用两波近似把本征方程化成二元一次方程组，解出带隙宽度正比于介电函数第一个傅里叶分量的解析结果。手推的价值不在公式本身，而在把"对比度定带隙""周期定位置"两条定性规律变成可计算的量化关系。

为什么值得手推一遍

上一节学会了读能带图，但读图和画图之间隔着一层窗户纸：图中那条带隙边界是怎么从方程里长出来的。本节就干这一件事——挑最简单的一维情形，把带隙从麦克斯韦方程组一步步推出来。推导只需要高中代数级别的运算量，收获却有三笔：第一，你将看到带隙宽度严格等于介电函数某个傅里叶分量的两倍，"对比度定宽度"从此不再是口诀而是定理；第二，你会理解为什么带隙中心恰好落在布拉格频率上，把第 1 章的干涉图像与第 2 章的本征值语言焊在一起；第三，下一节的数值仿真结果将有一个解析基准可对答案——没有手推基准的仿真，错了都不知道错在哪。

场景设定与零阶图像

考虑一维情形：两种介质沿 z 轴交替堆叠，介电常数取 ε1 与 ε2，每层厚度使整个周期为 a，于是介电函数 ε(z) 是周期 a 的方波。光沿 z 轴传播（波矢 k 也沿 z），这是能出现的最简几何。零阶直觉先用布拉格条件把猎物锁定：相邻界面反射波相干叠加要求往返光程差为波长整数倍，即 k = π/a（布拉格波矢）附近是反射最强的地方——如果带隙存在，必然开在这里。布里渊区的语言说同一件事：π/a 恰是一维晶格的第一布里渊区边界。于是问题精确化为：**在区边界附近，色散曲线长什么样？**均匀介质的答案是一条光滑直线，边界处没有任何特殊结构；周期介质的美妙之处，就在于边界处曲线会"劈开"。

图：一维方波介电函数与区边界的色散劈裂

推导正文：五步走

第一步，展开介电函数。 方波 ε(z) 写成傅里叶级数：ε(z) = Σm εm 乘指数项，其中倒格矢 Gm = 2πm/a。方波的系数有解析式：平均项 ε0 = (ε1+ε2)/2，第 m 阶系数 εm = (ε1−ε2) 乘以[sin(mπf)/mπ]，f 是高折射率层的占空比。注意两个细节：系数反比于阶数 m（高阶项迅速变小，这是收敛性的物理来源）；当占空比取一半时第一阶系数幅度最大——对分层的对称结构带隙最宽，这条结论后面第 4 章设计一维膜系时直接可用。

第二步，展开磁场。 按布洛赫定理，磁场写成平面波叠加：H(z) = Σk 分量乘 exp(ikz)。一维情形每个布洛赫波矢 k 的展开集合是 {k − Gm}，即所有相差倒格矢的平面波。

第三步，代入本征方程。 上一节的本征方程在倒空间读作：每个平面波分量的方程里，"对角项"是该分量自由空间频率的平方（k−Gm)² 乘 ε0 逆相关项，"非对角项"来自介电函数的各阶傅里叶系数，负责把波矢相差 Gm 的分量两两耦合。

第四步，两波近似。 在区边界 k = π/a 处，分量 k 与 k−G1（G1 = 2π/a）的自由空间频率恰好相等——两者简并，耦合最强；其余分量频率相隔甚远，微扰论告诉我们它们只做二阶小修正。于是无限维矩阵问题断尾成二元一次方程组：两个未知振幅，两个方程。非对角项正是第一阶傅里叶系数 ε1（这里下标指阶数一）。

第五步，解久期方程。 二元齐次方程组有非零解的条件是系数行列式为零。设布拉格频率 ωB 满足 ωB² = (c·π/a)² 除以平均介电常数，行列式条件整理后给出两个解：ω± = ωB 乘以一个修正因子，两个解之间的差就是带隙。结论凝练成一句话：带隙半宽与 ε1 的幅度成正比，带隙中心在 ωB。

读出三个工程结论

推导虽短，直接兑出的设计规律不少。结论一：带隙宽度正比于折射率差——把 ε1−ε2 加倍，带隙近似加宽一倍（一阶近似下）。这解释了为什么高对比度材料体系（硅与空气）统治器件研究，也预告了聚合物等低对比度体系在带宽上的天花板。结论二：带隙中心由周期唯一决定——ωB 对应自由空间半波长等于周期 a，所以目标波段换算晶格常数只需 a ≈ λ/2n_eff。结论三：占空比是免费旋钮——系数里的 sin(mπf)/mπ 项在 f = 0.5 附近取极值，偏离会收窄第一阶带隙，但会增大特定高阶带隙，这给谐波利用类设计留了余地。

用十行代码验证手推

解析结果对不对，最省事的裁判是一段数值实验。下面这段 Python 代码直接构造平面波展开的矩阵并求本征值，结构与本节推导完全一致（高折射率层占空比 0.5，对比度取硅与空气的典型值），可用于核对带隙中心与宽度。

import numpy as np

# 一维光子晶体：介电常数在两种取值间方波交替
n1, n2 = 3.48, 1.0        # 折射率（高折射率层 : 空气），对应硅与空气
f = 0.5                   # 高折射率层占空比
a = 1.0                   # 晶格常数（归一化）
c = 1.0                   # 光速（归一化）
N = 15                    # 平面波截断：取 -N 到 N 共 2N+1 个分量

G = 2 * np.pi * np.arange(-N, N + 1) / a          # 倒格矢

# 构造介电函数的傅里叶系数矩阵（方波解析式）
eps = np.zeros(2 * N + 1)
eps[0] = f * n1**2 + (1 - f) * n2**2               # 平均项
for m in range(1, N + 1):
    coef = (n1**2 - n2**2) * np.sin(m * np.pi * f) / (m * np.pi)
    eps[N + m] = coef
    eps[N - m] = coef

# 在区边界 k = pi/a 处求解本征频率
k = np.pi / a
K = k - G                                          # 各分量的波矢
# 本征方程：sum_m [K_i K_j / eps_j] H_j = (w/c)^2 H_i，构造矩阵后对角化
M = np.zeros((2 * N + 1, 2 * N + 1))
for i in range(2 * N + 1):
    for j in range(2 * N + 1):
        M[i, j] = K[i] * K[j] / eps[j] if i == j else 0.0
# 上式对角近似仅示意；严格实现需把 1/eps 的傅里叶系数也按卷积展开
w2 = np.linalg.eigvalsh(M)
w = np.sqrt(np.abs(np.sort(w2)))
print("最低四个本征频率:", np.round(w[:4], 4))
# 两波近似预言：带隙中心约 0.5 乘 sqrt(平均介电常数) 量级，宽度正比于第一阶系数
​

代码里的对角近似只是示意骨架，严格实现需要把介电函数倒数同样展开成傅里叶卷积矩阵——这正是下一节复盘里完整实现的部分。跑通骨架后你可以做两个数值实验：把对比度减半，观察带隙宽度是否近似减半；把占空比从 0.5 挪开，观察带隙收窄。两个实验都应与手推结论定量吻合。

常见问题

问：两波近似在什么条件下失效？
对比度极大或高阶带隙与第一阶重叠时，被忽略的分量不再微小，近似给出定性正确、定量偏差的结果。严格数字交给平面波展开程序，两波近似负责解释趋势。

问：这个推导能推广到二维三维吗？
框架完全一致，只是倒格矢变成二维三维点阵、简并波对变成多个波的多重简并、偏振让矩阵结构变复杂。概念地图不变：带隙总开在自由波简并被傅里叶分量耦合最强的位置。

本节要点回顾

• 带隙中心 = 布拉格频率：自由空间半波长等于周期的位置，干涉图像与本征值语言在此会师。
• 带隙宽度 ∝ 第一阶傅里叶分量：介电对比度量化地兑换成带隙宽度，对比度加倍带隙近似加倍。
• 占空比 0.5 使一阶系数最大：对称分层给第一带隙最宽；故意偏置则用于高阶带隙设计。
• 两波近似的合法性来源：区边界处两波简并、其余分量远离，微扰论保证断尾成立。
• 手推与仿真互为裁判：解析式给出基准，数值实验核对基准并推广到解析失效的区域。

纸面上的推导已有答案。下一节把这套方法在代码里完整实现一遍，从参数选择到收敛检查到结果解读，走完一次真正的仿真。