数值方法

RCWA 在周期开放结构与阈值分析中的应用

读完本章，你应能#

• 理解为什么在 第 5 章 建立能带概念、并在 第 13 章 用 PWEM 算出候选模目录之后，还必须引入 RCWA。

• 掌握 RCWA 的散射矩阵视角，明白它为什么特别适合处理层状周期开放结构。

• 学会区分“周期单胞被动共振结论”“光学阈值增益估计”和“有限器件级结论”这三个层次。

读前准备#

建议已掌握 第 13 章、第 10 章、第 7 章 中关于候选模筛选、耦合波图像和阈值条件的内容。

本章结构#

本章所处模型层级：L1→L2

第 5 章 已经说明怎样阅读能带图，第 13 章 又把它落实成了 PWEM 计算，于是到这里读者已经知道“候选模在哪里、怎样随几何变化、主导平面波分量大致是什么”。但这仍然不是面发射激光器的完整光学问题，因为 PWEM 还没有真正处理开放边界，也不能直接回答法向辐射如何离开器件。本章要做的事，就是把问题推进到“周期开放结构的散射与共振”层级。我们先说明为什么 RCWA 特别适合 PCSEL，再从分层和傅里叶展开建立散射矩阵框架，随后解释共振极点、被动谱和光学阈值估计之间的关系，最后明确它的边界：RCWA 适用范围宽，但它默认描述的仍是周期单胞，不是完整有限器件。

为什么在 PWEM 之后自然轮到 RCWA#

PWEM 解决的是被动本征问题，能告诉我们 Gamma 点附近有哪些候选模、主要由哪些平面波分量组成。但 PCSEL 的关键物理之一——模式如何向自由空间辐射——在 PWEM 的封闭本征框架里根本没有出现。

因此，接下来必须问的不再是”有哪些本征模”，而是”入射波和出射波如何经过层状周期结构发生散射”。这是 RCWA 的天然视角：它把周期结构看成散射系统，直接处理开放边界、入射通道与辐射耦合。对 PCSEL 来说，这一视角能直接回答三类问题：某个候选模是亮模还是暗模；某种几何扰动主要增强的是上方辐射还是下方辐射；以及在给定小信号增益下，开放极点何时靠近实轴。

与第 12 章的对比#

PWEM 擅长快速扫描晶格参数并找到候选模频段，但它默认封闭边界，无法区分向上辐射和向下辐射的贡献。RCWA 则以开放散射问题为出发点，可以直接计算$S_{11}$（反射）和$S_{21}$（透射），从而给出上下出光比例。

RCWA 的基本思路：平面内做傅里叶，纵向做分层#

RCWA 的名称里虽然有“耦合波”，但它和前面 CWT 的工作方式并不相同。CWT 的核心是主动保留少数主导波分量，用低维耦合模型换取解释性；RCWA 则更接近一种严格的半解析数值方法，它在每一层中对平面内周期方向做傅里叶展开，而在纵向把结构切成若干层逐层求解。

设器件在 $x$ 和 $y$ 方向周期，在 $z$ 方向分层。对每一层 $\ell$，把介电常数与场都在平面内写成傅里叶级数，就得到一个有限维模态系统。该层中的场可写成“向上”和“向下”传播的模叠加：

其中 $\bm{f}_\ell$ 是该层选取的场变量向量，$\mathbf{W}_\ell^{\pm}$ 是本征模基底矩阵，$\mathbf{\beta}_\ell$ 给出纵向传播常数，$\bm{a}_\ell^{\pm}$ 分别表示上、下行模振幅。 notebox 纵向传播常数矩阵 $\mathbf{\beta}_\ell$ 的定义。在周期层中，纵向传播常数应由层本征值问题给出：

这里 $\mathbf{M}_\ell$ 是由该层 Fourier 卷积矩阵构造的 RCWA 层算子，$\mathbf{W}_\ell$ 为本征模基底矩阵。仅在均匀层（无谐波耦合）时，才可退化为单谐波形式

当 $\beta_n$ 为纯虚数时，对应倏逝波；当为实数时，对应传播波。 notebox 这一步的物理意义非常明确：在每一层里，平面内周期导致的多谐波耦合已经被吸收进层本征模中；而真正的结构复杂性，被转化成"层与层之间如何通过边界条件连接"。于是，复杂三维问题被拆成两部分：层内求模，层间拼接。图 14.1直观展示了如何将连续的倾斜光栅剖面离散化为多层阶梯近似结构，这是RCWA处理任意形状周期结构的核心技巧。

图 14.1揭示的离散化思想是RCWA能够处理复杂周期结构的关键：任何连续的剖面都可以用足够多的薄层来逼近，而每一层内部的平面内周期性使得我们可以用傅里叶级数展开。这种"分层+展开"的策略，使 RCWA 可以在保持半解析离散结构的同时处理复杂剖面。

为什么散射矩阵是 RCWA 的真正核心#

如果只是逐层传递振幅，最直接的想法是用传输矩阵（transfer matrix）。但对高对比度、强倏逝波或厚层结构，这样做在数值上常常不稳定，因为指数增长和指数衰减模式会在矩阵乘法中导致严重病态。RCWA 在工程上较稳定的做法，是使用散射矩阵（scattering matrix）递推。

把“不稳定”说得更定量一点：若某层的纵向本征常数写成 $k_z=\ii\kappa$，则该层传播因子会同时含有

在简单传输矩阵连乘中，指数衰减项和指数增长项会被强行放进同一矩阵乘法链，结果条件数往往随层数和倏逝强度指数恶化。散射矩阵递推的价值就在于：它把入射、反射和透射通道分开追踪，避免把“物理上应当衰减的分量”和“数值上会炸掉的增长因子”直接混在同一个乘法序列里。

把整个结构视为一个多端口散射体，就有

其中 $\bm{a}_{\mathrm{up}}$ 与 $\bm{a}_{\mathrm{down}}$ 表示从上、下方入射的通道振幅，$\bm{b}_{\mathrm{up}}$ 与 $\bm{b}_{\mathrm{down}}$ 则是对应的出射振幅。$\mathbf{S}$ 编码了在给定频率和面内 Bloch 波矢下，整个周期单胞如何把入射通道重排为出射通道。

若把 式 (14.2) 写成块矩阵形式，

则 $\mathbf{S}_{11}$ 表示“从上方入射再向上方反射”的块，$\mathbf{S}_{21}$ 表示“从上方入射后向下方透射”的块；其余两块则对应从下方入射的情形。把这层物理含义紧贴在 式 (14.2) 后面，是为了让后面的 $S_{11}$、$S_{21}$ 不再像纯软件符号。

若只保留零级向上传播与向下传播通道，那么最常见的两个标量散射参量就是

其中下标 $0$ 表示零级衍射通道，通常要求端口振幅已经按单位功率通量归一化。这样处理后，$|S_{11}|^2$ 与 $|S_{21}|^2$ 才能直接对应反射率与透射率的主通道贡献。对存在多个衍射级次的周期结构，较稳妥的做法是分别统计各传播级次的功率通量，再把它们与同一个入射通量相除。明确这一步后，反射谱、透射谱和上下出光比例都可以对应到具体衍射通道的功率比，而不是只对应求解器输出的一条曲线。

对 PCSEL 来说，式 (14.2) 的威力在于它直接把“内部结构”和“外部辐射通道”连在一起了。于是，上下辐射不对称、偏振选择、法向出光窗口和反射谱中的共振尖峰，都能在同一个框架下被读取。不过，开放散射仍是被动层面的问题，它不会自动给出激射结论2。

Fourier factorization：由界面连续性决定的离散化规则#

RCWA 还有一个经常被忽略、却直接决定收敛质量的技术点：Fourier factorization。对具有介电常数跳变的周期结构，像 $\varepsilon E$ 这样的乘积项并不能一律用“傅里叶系数逐项卷积”来处理。Li 的 factorization 规则要求根据场分量在界面处的连续性，选择对 $\varepsilon$ 还是对 $\varepsilon^{-1}$ 构造 Toeplitz 卷积矩阵[Li, 1996]。记 $[\varepsilon]_{mn}=\hat{\varepsilon}_{m-n}$ 为对介电函数傅里叶系数构造的 Toeplitz 矩阵，$[\varepsilon^{-1}]$ 为对 $1/\varepsilon$ 同样构造的矩阵。若乘积中的场分量在界面处连续（如 TE 偏振的电场切向分量），可直接使用 Laurent 规则

若法向位移矢量更自然连续（TM 偏振下电场法向分量在界面跳跃），则应使用 inverse rule

用错规则的后果通常是伪收敛：Gibbs 震荡增强，收敛速度从指数级降为代数级，并出现系统性频移，而非发散。以 GaAs/空气二元光栅（折射率对比约 3.4:1）为例，错误的卷积规则可使 TM 偏振共振峰位产生 5--10 nm 的系统性偏移，且把傅里叶阶数从 $N=5$ 提高到 $N=21$ 后偏移仍不消失。因此，如果 RCWA 预测的共振波长与实验测量总有几纳米的固定偏差，在排除工艺误差后，应优先检查 factorization 规则是否与界面连续性一致。式 (14.4)、式 (14.5) 想让读者记住的只有一件事：Fourier factorization 是由界面连续性决定的离散化规则，而非可由提高截断阶数自动消除的数值细节。

被动谱、共振与极点：谱线为什么会尖起来#

很多读者第一次用 RCWA 时，最直观的输出是一条反射谱或透射谱。看到尖峰或窄谷之后，容易把它直接称为“高 $Q$ 模”或“低阈值候选模”。这种说法并非总错，但如果不说明其背后的极点结构，就会显得很经验化。

从散射矩阵观点看，开放共振模对应于复频率平面上散射矩阵的极点。通俗地说，如果系统内部存在一个能量停留时间较长的开放模，那么当实频扫描靠近该模的共振位置时，谱线就会出现尖锐变化。谱峰、谱谷、相位跃迁、Fano 线形，本质上都与极点及其和背景散射的干涉有关。

因此，RCWA 中看到尖谱线时，真正应该问的不是“峰高不高”，而是：

• 该特征是否随傅里叶阶数和分层数收敛；

• 它对应的是哪一支候选模，和 PWEM/CWT 的哪一支模族能对上；

• 其线宽对应的被动损耗量级是否与其他方法一致；

• 它主要通过哪个辐射通道向外耦合。

只有把这些问题答清楚，被动谱才能转化为器件设计信息。

随书实算：从 PWEM 频窗到 RCWA 反射谱#

图 14.2 是脚本 code/rcwa_square_hole/run_rcwa_square_hole.py 用 grcwa 实际计算的结果。模型取正方晶格圆形空气孔，$r/a=0.30$，高折射率层 $n=3.4$、厚度 $0.5a$，上方为空气、下方衬底 $n=1.5$；请求 49 个 Fourier 级次，求解器实际保留 45 个。被动扫描在 $\Omega=a/\lambda\approx0.43$ 附近出现明显的 Fano 型反射结构，且全频窗内 $|R+T-1|$ 的最大误差小于 $10^{-12}$。

这份算例把“PWEM 告诉 RCWA 去哪里看”变成了可执行接口：先在候选窗口做无源宽扫，再对谱线结构细扫，最后才改变增益代理。复现时至少应同时提高 Fourier 阶数与实空间采样网格；若峰位、线形或 $R+T$ 守恒随截断显著变化，就不能把该谱线用于模式排序。

RCWA 怎样进入光学阈值分析#

前面 第 7 章 已经明确区分了阈值增益和阈值电流。RCWA 在这里的作用，是估计周期开放结构层面上的光学阈值增益，而非直接给出 $I_{\mathrm{th}}$。最常见的做法，是在有源层引入小信号增益对应的等效虚部，例如写成

然后扫描增益参数 $g$，寻找使目标共振从净衰减变为恰好补偿损耗的条件。

如果只停在这一步，读者容易把“扫到线宽为零”误听成“已经建立了完整激光阈值理论”。教材写法必须把这件事分层说清。对 PCSEL 而言，至少应区分以下四个层次：

1. 冷腔开放共振层：先在无增益条件下找到复频率极点 $\tilde{\omega}_0=\omega_0-\ii\gamma_0$，它描述的是开放系统如何储能与泄能。

2. 线性小信号阈值代理层：在给定小信号增益模型后继续追踪极点 $\tilde{\omega}(g)$，并用 $\Im \tilde{\omega}(g_{\mathrm{th}})=0$ 估计光学阈值增益。

3. 阈值处自洽有源场层：真正的阈值问题要求场、极化和反转分布同时满足自洽条件，尤其当增益色散明显时，实部与虚部不能被拆成两个彼此独立的旋钮。

4. 阈上稳态激光层：一旦超过阈值，增益饱和、空间烧孔、模式竞争和注入非均匀性都会进入主导地位，这时需要更完整的有源激光形式化，例如基于 Maxwell--Bloch 方程的稳态或时域求解。

本章 RCWA 所做的，严格地说位于第二层。它有用，但它仍是线性小信号阈值代理，不是阈上稳态激光解。

在这个过程中，可以从几种等价角度理解“达到阈值”：

• 被动极点在复频率平面上向实轴靠近；

• 某个共振的净线宽收缩到零；

• 在小信号线性框架中，目标模的总损耗被增益恰好补偿。

这一步得到的是 $g_{\mathrm{th}}$ 这类光学量。它对设计极其有用，因为它告诉我们结构本身对目标模和竞争模提出了多高的增益要求。但必须再次强调，$g_{\mathrm{th}}$ 不是 $I_{\mathrm{th}}$。从前者到后者，还需要材料增益模型、漂移-扩散输运、复合和热反馈。

为什么这只是线性小信号阈值代理#

式 (14.6) 默认我们已经把有源区压缩成了一个“可调复介电常数”模型。只要这一模型仍停留在线性响应层面，那么我们实际上做的是：给开放冷腔模外加一个尚未饱和的增益补偿，并观察净衰减是否被抵消。这个判据与 第 7 章 的阈值条件相容，但它并没有告诉我们阈上光场会稳定在什么幅度、是否发生多模竞争、反转是否仍近似均匀，也没有告诉我们空间烧孔会不会把模式排序翻转。

增益色散为什么要求更谨慎#

在许多工程扫描中，人们把 $\varepsilon''_g$ 当成独立参数扫动，而把 $\varepsilon'$ 固定在冷腔值附近。这样做的目的是把问题简化成“先补偿损耗，再看哪一支模最早到零线宽”。作为前端筛选，这种近似经常足够。但若增益谱较窄、工作点离谱峰不远、或载流子引起的色散不可忽略，那么 $\varepsilon'(\omega)$ 与 $\varepsilon''(\omega)$ 由因果材料模型联系在一起，不能任意独立指定。实际有源介质是会同时改变频率位置、线宽和模场重叠的复响应系统，而非简单的“只加负虚部”。

阈值频率并不必然等于冷腔共振频率#

把阈值理解成“在冷腔频率上把线宽调到零”也过于粗糙。更准确的说法是：随着增益参数变化，极点本身会沿复频率平面移动，

阈值条件是 $\gamma(g_{\mathrm{th}})=0$，而真正的阈值频率应写成

它不一定与冷腔共振频率 $\omega_0$ 完全相同。造成这种偏移的原因包括：增益色散引起的频率拉动、场分布随增益微调而改变，以及不同辐射通道和材料色散共同作用下的复频率漂移。对 PCSEL 来说，这种偏移尤其重要，因为 $\Gamma$ 点附近的模式排序本来就可能很近，微小的频移就足以改变“哪一支模先达到阈值”的判断。

它与更完整 active lasing formulations 的关系#

从更宽的理论坐标系看，本章 RCWA 增益扫描与更完整的有源激光稳态理论之间是“前端代理”和“自洽主问题”的关系。前者把目标模的阈值判定简化成线性极点继续问题，因此很适合做候选模筛选和参数排序。后者则把 Maxwell 方程与极化、反转和饱和非线性联立起来，直接求解阈值处或阈上稳态激光解；典型代表包括 Maxwell--Bloch 方程的直接时域求解，以及稳态自洽形式如 steady-state ab initio laser theory（SALT）[Türeci et al., 2007] [Türeci et al., 2009]。本书在这里不展开 SALT 的完整推导，但读者必须知道：RCWA 扫 $\varepsilon''<0$ 到零线宽，是为真正的有源激光问题提供一个计算成本较低、可解释、可筛选的线性近似入口，而非激光理论本体。

为什么 RCWA 特别适合做“单胞层面的设计判断”#

一旦接受“RCWA 主要给光学阈值增益和辐射通道信息”，它在工作流中的位置就非常清晰了。它适合做以下几类判断：

第一，改变孔半径、占空比、刻蚀深度或双晶格扰动时，目标模与竞争模的被动损耗顺序怎样变化。

第二，上下辐射比例、不同偏振通道强弱和法向出光窗口如何随单胞几何变化。

第三，某种外延层微调是否主要通过改变纵向重叠影响辐射，而不是简单改频率。

第四，在后续有限器件仿真之前，哪些参数区间根本不值得继续深挖，因为开放周期单胞层面就已经显示目标模没有优势。

这也是 RCWA 的工程价值所在：它以远低于全尺寸三维时域仿真的代价，把大量“不可能成功”的单胞方案尽早筛掉。

RCWA 不能替你解决什么问题#

RCWA 适用范围宽，但它解决的是“周期单胞开放散射”问题，而不是“完整器件工作”问题。这个边界必须明确。

首先，它默认平面内无限周期。因此，有限孔阵边缘泄漏、横向模包络和大面积器件的边界损耗不在默认描述范围内。若这些效应开始主导，必须转向 第 15 章 的 FDTD/FEM 或更高阶包络模型。

其次，它通常在小信号线性框架里处理增益。于是，增益饱和、空间烧孔和工作点非线性不会自动出现在结果里。

进一步说，即便把小信号增益继续到“零线宽”，也仍未自动得到阈上稳态解。若问题已经进入强色散、强模竞争或反转明显不再静止的区域，就必须转向更完整的有源激光稳态理论，或至少用它们来校准 RCWA 前端筛选的可信边界。

再次，它不包含真实电流路径。均匀增益分布在 RCWA 中只是一个光学代理模型，而不是电注入器件的默认现实。

最后，它也不解决温升问题。很多结构在室温被动谱上看起来很好，一旦工作后温升导致折射率和增益谱漂移，目标模顺序就可能改变。这类问题必须交给热-电-光耦合章节处理。

2024--2026 前沿补充：modified RCWA 正在覆盖 full-PCSEL 结构#

近年一个值得单独标注的进展是 modified RCWA 对 full-PCSEL / VCSEL 结构的扩展：通过双向散射矩阵与分层拼接策略，模型不再局限于“仅从上/下单向入射”的传统流程，而可以更直接地处理“腔内源同时向上下传播”的器件场景[Xu et al., 2024]。这类方法在共振波长、发散角等关键量上与 FDTD 给出了可比精度，同时保留了更高扫描效率。

这对本章方法边界的影响是双面的：

• 正面：RCWA 可承担更大比例的前端筛选任务，尤其适合大参数空间下的快速候选淘汰；

• 边界：当问题进入强非线性增益、强有限尺寸边缘效应或多物理场强耦合时，仍需回到全波与器件级模型完成回写链。

因此，modified RCWA 应被视为“前端能力增强”，不是“对后端验证层的替代”。

收敛与数值卫生：哪些检查是必须交代的#

RCWA 对数值细节敏感。一个结果是否可信，不取决于图形是否平滑，而取决于它是否通过基本的数值检查。至少应当检查四件事。

第一，傅里叶阶数收敛。目标模的谱位、线宽、反射峰形和上下辐射分配，不能随着阶数增加继续大幅变化。

第二，纵向分层收敛。对刻蚀轮廓、锥形孔或多层外延而言，过粗的分层会把真实边界模糊掉。

第三，能量守恒或功率平衡检查。在无源条件下，反射、透射和吸收的总和应当满足相应守恒关系。

第四，结论对微小参数扰动的鲁棒性。如果某个“最好方案”只在一个极细的参数点上存在，而周围稍微变化就消失，那么它通常不是可靠设计。

这些检查之所以重要，不只是为了数值严谨，也是为了避免把求解器伪影误读成器件物理。

RCWA 在整条工作流中的准确位置#

现在可以把 RCWA 放回全书主线。第 13 章 的 PWEM 给出候选模目录；本章的 RCWA 则对这些候选模做开放周期结构层面的精修。接下来 第 15 章 的 FDTD/FEM 会把有限尺寸和更真实几何带回来，再由 第 17 章、第 18 章、第 19 章、第 20 章 把这些光学结果接进外延、增益、电注入和热耦合模型。

也就是说，RCWA 是一个关键的中间关口，而非终点。它回答的问题是：如果先只讨论单胞级开放光学，哪些方案值得继续投入更昂贵的器件级验证。把这个关口守好，后面的工作量会大幅下降；把它省掉，后面的全器件仿真往往会变成盲目试错。

回看本章#

RCWA 的核心是散射矩阵视角：把周期分层开放结构看成一个输入输出系统，直接分析内部模式如何与外界辐射通道交换能量，而非某个软件界面操作。对 PCSEL 而言，它特别适合处理法向辐射、偏振通道、被动共振和光学阈值增益估计。但本章给出的阈值分析属于线性小信号阈值代理：它通过继续追踪开放极点来估计哪一支模最先补偿净损耗，而不是直接求出阈上稳态激光解。只有把这一层级和更完整的有源激光理论分清，RCWA 才能被放在正确位置上使用：作为“开放单胞精修器”和“后续器件仿真的前置筛选器”，而不是最终器件结论的替代品。

练习题#

1. [基础] 用自己的话解释为什么 PWEM 之后自然需要 RCWA，而不是直接跳到器件级电热模型。

2. [基础] 说明散射矩阵方法为什么在数值上通常比简单传输矩阵更稳定。

3. [进阶] 讨论从 RCWA 得到光学阈值增益后，为什么仍然不能直接报告阈值电流。

4. [进阶] 解释为什么“在固定 $\varepsilon'$ 下扫描负虚部直到零线宽”属于线性小信号阈值代理，而不是完整激光阈值理论。

5. [综合] 设计一套你认为最基本的 RCWA 收敛检查清单，并说明每一项各自防止什么误判。

延伸阅读#

itemize 下一章将把视角从“周期单胞开放结构”推进到“有限尺寸真实器件”，也就是 FDTD 和 FEM 的领域。读完下一章后再回看本章，会更清楚两者的分工边界。

若想打牢开放周期结构的基础，可先对照 [Joannopoulos et al., 2008] [Sakoda, 2005] 中关于周期散射和开放通道的章节阅读。

若想补充 RCWA 的经典原始文献，可参考 [Moharam & Gaylord, 1981] [Li, 1996]。

若想理解本章线性增益扫描与更完整有源激光理论之间的关系，可进一步阅读 [Türeci et al., 2007] [Türeci et al., 2009]。 itemize