4.1 一维光子晶体

4.1 一维光子晶体

本节摘要：一维光子晶体是介电常数仅沿一个方向周期变化的多层膜结构，是整个家族里最简单、也唯一实现大规模商业化的形态。本节给出特征矩阵法的完整工具链（可手算可编程），拆解布拉格镜、窄带滤光片、宽带减反三种典型膜系，并用传递矩阵代码现场演示带隙的计算与设计。

一维：最简单的维度，最广泛的落地

维度阶梯的最底层是本章主线的起点，选它开局有双重理由。物理上，一维情形简单到可以把第 2 章的所有概念全部走一遍数值流程——每个界面、每次反射都看得见摸得着；工程上，它是不折不扣的商业冠军：增透膜、高反镜、滤光片、防晒涂层里的红外屏蔽层，镀膜工业每年量产的"一维光子晶体"以面积计可能是天文数字，只是行业里叫它多层膜而非光子晶体。第 1 章讲过的历史也在提醒我们：周期多层膜的选择性反射现象十九世纪就被瑞利研究过——这个维度上的"光子晶体"，人类其实已经造了一百多年。

结构与物理图像

一维光子晶体是两种介质沿 z 方向交替堆叠的层状结构，周期 a 由两层厚度决定，光垂直或斜入射。它的带隙是方向性的：沿堆叠方向，布拉格条件满足的频段被禁行——表现为高反射；垂直于堆叠方向的传播完全不受影响。物理图像就是第 1 章的合唱模型：每个界面的菲涅尔反射很小，但数十层按四分之一波长厚度堆叠后，反射波全部同相叠加，反射率轻松超过 99.9%。带隙的中心波长由"光学厚度"决定——每层的光学厚度（几何厚度乘折射率）取目标波长的四分之一时，第一带隙恰好覆盖目标波长，这就是镀膜行业祖传的"四分之一波长堆叠"公式。

特征矩阵法：一维的屠龙刀

一维多层膜有一套优雅到近乎奢侈的解析工具：特征矩阵法（也称传递矩阵法）。每层膜用一个 2 乘 2 矩阵描述，矩阵元素由该层折射率、厚度与光波在层内的相位厚度决定；整堆膜的矩阵就是各层矩阵的连乘；最后一个矩阵元素直接给出反射率与透射率。它的美妙之处有三：可手算——三四层膜的矩阵乘法一张草稿纸搞定；可编程——十行代码算任何层数的任意堆叠；可反演——给定目标反射谱，可以数值优化出膜系结构，这正是镀膜工程师的日常工具。相比第 2 章的平面波展开法，特征矩阵法不需要周期性假设，天然处理有限层数、任意层厚、含缺陷层的结构——一维世界里它比布洛赫方法更顺手。

图：四分之一波长堆叠：反射波的同相合唱

三种典型膜系：一维缺陷工程的三个段位

膜系设计的世界里，第 2 章的缺陷态理论以最朴素的方式呈现。段位一：纯堆叠布拉格镜——L 与 H 交替各十余层，目标是带隙内高反。设计自由度只有层数与两种材料，稳定性极高，是激光谐振腔高反镜、传感器窗口的标准答案。段位二：插层窄带滤光片——在堆叠正中插入一层光学厚度为半波长的"腔层"，带隙中央劈出一条窄透射峰：结构上这正是第 2 章讲的一维缺陷腔，物理上是法布里珀罗干涉仪的固体版。透射峰的宽度由两侧镜子的反射率决定——镜子越"高"，峰越窄，滤波选择性越强。段位三：非整数层优化——放弃四分之一波长的整齐阵列，让数值优化器在几十层的厚度空间里搜索最优解，逼近任意目标谱形。增透膜的完美曲线、陡边截止滤光片都是这么来的。三段位合成一条学习路线：先手算段位一建立直觉，再用代码玩段位二理解缺陷，最后敬畏段位三里优化器与工艺容差的博弈。

现场演算：二十行代码算一座布拉格镜

下面这段完整代码用特征矩阵法计算一个实际膜系在通信波段的反射谱，可以直接运行，改材料与层数即可变成自己的设计工具。

import numpy as np

# 光学常量：两种膜层材料与基底（近红外常用组合）
nH, nL = 2.10, 1.45          # 高折射率层（五氧化二铌）与低折射率层（二氧化硅）
n_sub = 1.52                 # 玻璃基底
pairs = 12                   # HL 交替的对数
wl0 = 1550.0                 # 设计中心波长（纳米）
dH, dL = wl0 / (4 * nH), wl0 / (4 * nL)   # 四分之一波长厚度

def char_matrix(n, d, wl):
    """单层特征矩阵：相位厚度 delta，导纳 n（垂直入射）"""
    delta = 2 * np.pi * n * d / wl
    return np.array([[np.cos(delta), 1j * np.sin(delta) / n],
                     [1j * n * np.sin(delta), np.cos(delta)]])

def reflectance(wl):
    """整堆膜的特征矩阵连乘 → 反射率"""
    M = np.eye(2)
    for _ in range(pairs):
        M = char_matrix(nH, dH, wl) @ M
        M = char_matrix(nL, dL, wl) @ M
    # 特征导纳 → 反射系数
    Y = (M[1, 0] + M[1, 1] * n_sub) / (M[0, 0] + M[0, 1] * n_sub)
    r = (1 - Y) / (1 + Y)
    return abs(r) ** 2

# 扫描 1200 到 2000 纳米，输出带隙边缘与峰值
for wl in range(1200, 2001, 100):
    print(f"波长 {wl} 纳米: 反射率 {reflectance(wl) * 100:6.2f}%")
# 预期输出形态：1300-1900 纳米区间反射率逼近 100%，
# 带隙中心 1550 纳米，边缘处反射率陡降到个位数
​

运行后可以立刻做三个实验：把对数从 12 减到 6，看带隙边缘反射率的斜坡变缓（镜子变"矮"）；把 nH 与 nL 的差缩小，看带隙变窄（对比度定律的一维版）；在堆叠正中插入一层半波厚的高折射率层，看反射谱中央裂出一条透射峰——三分钟做完三个定理的验证，这就是一维工具链的爽快感。

常见问题

问：一维带隙能挡住斜入射的光吗？
部分能。斜入射时带隙向短波移动且 TE 与 TM 的带隙边缘劈开，s 与 p 偏振的行为不同；全角度高反需要更复杂的共轭膜系设计，这也是"全向反射镜"曾是研究热点的原因。

问：为什么镀膜工程师不叫它光子晶体？
历史与实用双重原因：多层膜工业比光子晶体概念早几十年，特征矩阵工具链自成体系。概念上两者完全同源，理解了一维光子晶体就理解了镀膜工业的物理底盘。

膜系设计的两个实用技巧

收尾前补两个镀膜行业沉淀下来的实用技巧，它们在一维设计里出现频率极高。技巧一：等效层思想。 斜入射时一层膜可以按等效折射率与等效厚度拆算，利用这个性质，用两种材料就能近似出"第三种折射率"的层——材料受限的实验室里，这是凑出任意膜系的常用魔术。技巧二：监控方式决定误差分布。 膜厚监控用石英晶振（绝对厚度）还是光学极值法（相对光学厚度），决定了厚度误差在层间是独立还是相关；光学极值法让相邻层的误差自动抵消一部分，是窄带滤光片行业的偏好。两个技巧共同指向一个认知：一维设计的功夫一半在纸上、一半在制造与监控的选择里——这也是它作为"唯一大规模商业化形态"的深层原因：一百年的工艺沉淀，让每一个设计自由度都有配套的制造手段。带着这份认知进入二维，你会更清楚平板工艺还欠哪些账。

本节要点回顾

• 一维是最老牌的光子晶体：瑞利时代就有物理，镀膜工业已量产百年；带隙方向性限制它只管一个方向。
• 四分之一波长堆叠公式：光学厚度定带隙中心，折射率比定带隙宽度，层数定反射率峰值。
• 特征矩阵法三连：可手算、可编程、可反演；一维世界里比布洛赫方法更贴合工程。
• 三段位膜系：纯堆叠镜、插腔窄带滤光（一维缺陷腔）、数值优化任意谱；学习路线按段位递进。
• 二十行代码是本节的通行证：跑通并完成三个验证实验，一维维度即告毕业。

一维的纵深讲完，维度阶梯升到二级。下一节看二维平板如何用"面内带隙加面外折射率"的双保险成为产业主力。