RCWA 严格耦合波分析原理：周期性结构光学响应的高效计算

引言：当”光栅方程”不够用时

站在纱帘后面看一盏小夜灯，会看到灯周围出现一圈彩色光带——纱帘的网格其实是一张天然的衍射光栅，周期只有一两毫米，却足以把白光按波长分开。把周期缩小到几百纳米，同样的现象就发生在光刻对准标记、光谱仪的核心元件、太阳能电池的表面纹理和超表面的纳米柱阵列上：一个周期性的微纳结构，会把一束入射光”拆”成许多不同方向的衍射级。

图 1：纱帘网格作为天然衍射光栅时，小夜灯发出的光被分为不同颜色的衍射级，图中标出了各级的位置。

周期结构的衍射方向可以用光栅方程 $d(\sin\theta_i + \sin\theta_m) = m\lambda$ 快速算出：周期 $d$、入射角 $\theta_i$、波长 $\lambda$ 和整数 $m$ 一起决定了第 $m$ 级衍射的方向。但方向只是”光往哪里走”，工程师更常问的是”每个方向有多少光”——也就是各级衍射效率。效率取决于结构内部的电磁场如何被散射、吸收和耦合，这在周期接近波长甚至小于波长时变得非常复杂：纳米柱之间会发生多次散射，金属结构会激发表面等离激元，电场在材料边界上剧烈变化。这时候，把结构当成”理想薄屏”的标量近似、以及只算几何光程的光栅方程，都不再可靠。

这正是严格耦合波分析（Rigorous Coupled-Wave Analysis，RCWA）登场的场景。RCWA 是一种半解析的频域方法：它不把整个空间切成时域网格，而是把周期结构内部的电磁场用一组”空间谐波”展开，再把麦克斯韦方程组变成可以高效求解的矩阵问题。它由 Moharam 和 Gaylord 在 1981 年系统提出，如今是光栅、光子晶体和超表面仿真中最常用的严格方法之一，被大量商业与开源仿真工具实现。本文从”为什么要严格求解”出发，逐步讲清 RCWA 的展开、分层、边界匹配三步，再解释稳定实现与收敛性这两个工程上最关键的细节。

二、核心思想：把周期结构”翻译”成傅里叶空间的模式

RCWA 的出发点是一个数学定理：Floquet 定理。周期为 $\Lambda$ 的结构里，电磁场不会以任意形式变化，而一定可以写成周期函数与平面波因子的乘积；换句话说，场可以展开成一组离散的空间谐波。用公式表达，在 $x$ 方向周期为 $\Lambda$ 的结构中，场可以写作

$$
E(x,z) = \sum_{m=-\infty}^{+\infty} E_m(z) e^{i k_{xm} x},
\qquad k_{xm} = k_0 \sin\theta_i + \frac{2\pi m}{\Lambda},
$$

其中 $k_0$ 是真空波矢大小，$\theta_i$ 是入射角。每一项对应一个”衍射级”，它的横向波矢 $k_{xm}$ 正是光栅方程给出的方向。也就是说，光栅方程里的衍射级 $m$，在 RCWA 里天然就是傅里叶展开的阶数——物理图像和数学展开在这里完全对应。

图 2：透射衍射光栅的几何示意：周期为 d 的周期结构把入射光散射到不同衍射角 θ 的方向，相邻周期的散射光在特定角度相长干涉。

图 3：衍射级的方向由入射波矢与周期共同决定，各级波矢对应 Floquet 展开中的不同空间谐波。

现实中我们不可能保留无穷多项。RCWA 的做法是截断：只保留 $m = -N$ 到 $+N$ 的谐波，把无穷级数变成一个有限的向量和矩阵问题。截断阶数 $N$ 越大，近似越精确，但矩阵尺寸随之增大，计算量通常按阶数的立方量级增长。于是”算得准”和”算得快”之间的平衡，就落在截断阶数的选择上——这是后文收敛性讨论的起点。

确定了展开方式之后，RCWA 还需要处理结构在垂直方向上的变化。光栅和超表面通常不是均匀的一整块，而是由衬底、薄膜、纳米柱等叠成的多层结构，每一层的折射率分布还可能随位置变化。RCWA 的处理办法是：把结构沿垂直方向（z 方向）切成若干薄层，每层内部的介电常数在 z 方向近似均匀；斜侧壁或连续曲面就用阶梯状薄层去逼近（称为 staircase approximation）。于是，”任意形状的周期结构”被翻译成”一串沿 z 方向均匀、在 x-y 平面内周期的薄层”。

三、分层与模式：RCWA 的数学结构

有了薄层划分，剩下的问题就是：对每一层，求出其中的电磁模式，再把层与层连接起来。

先看单层。层内介电常数只随 x-y 变化、不随 z 变化，而场已经用空间谐波展开。把这个展开代回麦克斯韦方程组，偏微分方程就化成一个以 z 为自变量的矩阵常微分方程：每一阶谐波的振幅随 z 如何演化，由该层材料的介电常数矩阵决定。求解这个矩阵问题（一个特征值问题），得到的就是该层内能够稳定存在的模式——每个模式对应一个特征向量和一个传播常数，物理上就是这一层”允许”的电磁场本征形态。这一步与波导分析中求解模式非常相似，只是这里的模式是一组离散的衍射波。

再来看层与层的连接。每一层的模式解给出层内的场分布，但要得到整个结构的响应，还必须让电场和磁场的切向分量在每个界面上连续。最直接的做法是逐层传递：从入射层开始，把场”传播”到出射层。这个方法思路简单，却在层数多或折射率对比大时数值上很不稳定——误差会在逐层传递中指数放大，导致结果完全失真。RCWA 的现代实现因此普遍改用散射矩阵（scattering matrix，S 矩阵）：每次只处理一层，把该层与周围层当成一个”散射问题”，用矩阵记录从各个方向入射的波如何被反射和透射，再把相邻层的 S 矩阵组合起来。这样一层一层”拼接”下去，最终得到整个多层结构的总反射和总透射。S 矩阵方法在数值上稳定得多，也因此成为 RCWA 乃至整个傅里叶模态法家族的标准做法。

四、稳定与收敛：决定成败的两个细节

RCWA 的概念在 1981 年就已提出，但真正成为工程上可靠的工具，靠的是两个后来才被解决的问题：数值稳定性和截断收敛速度。

稳定性问题来自边界连接方式。如前所述，早期实现用透射矩阵（T 矩阵）逐层传递，在层厚、层数或折射率对比度增大时，矩阵中的数值会剧烈增长，产生所谓的”数值不稳定”。1995 年，Moharam 等人给出了一套面向二元光栅的稳定高效实现，核心就是用增强透射矩阵（enhanced transmittance matrix，ETM）或等效的散射矩阵方式组织计算，把误差限制在逐层的局部范围内。如今各类 RCWA 代码几乎都以这类”稳定公式”为骨架，S 矩阵也成为标准组织方式。

收敛问题则更微妙。RCWA 是傅里叶空间方法，对介电常数突变的结构，截断后的傅里叶级数会在跳变点附近出现振荡——这就是 Gibbs 现象。对于横磁（TM）偏振或高折射率对比的结构，直接对介电常数做傅里叶展开再卷积，收敛会非常慢，甚至出现错误结果。1997 年，Li Lifeng 系统研究了这一问题：当两个具有共同跳变点的不连续函数相乘时，它们的傅里叶级数卷积必须遵守特定的”因子化”规则（即逆规则，inverse rule），才能保证截断后仍然收敛到正确解。把这条规则嵌入 RCWA 的介电常数处理，收敛速度显著改善。后续发展的快速傅里叶因子化（fast Fourier factorization，FFF）和矢量傅里叶模态法（vector FMM）进一步推广了这类思想：在单元内每个点上建立局部的法向/切向坐标系，把场的法向与切向分量区别对待，从而在任意形状、任意对比度的结构上都获得稳定快速的收敛。

五、工程中的 RCWA：能算什么，不能算什么

RCWA 在工程中的定位很明确：它最适合”平面波入射 + 周期性分层结构”这一类问题。超表面单元（一个周期内的纳米柱阵列）、衍射光栅、光子晶体板、抗反射纹理、半导体测量中的周期性沟槽，都是它的主场。它算得又快又稳，因此特别适合需要扫描大量参数的场景，比如超表面单元的相位-半径曲线扫描、光栅槽深与占空比的工艺窗口搜索。开源与商业实现都很成熟：Ansys/Lumerical 的 RCWA 求解器、斯坦福的 S4、Facebook Research 的 fmmax 等，都是实际可用的工具；fmmax 用 JAX 实现，支持 GPU 加速和自动微分，后者意味着可以把 RCWA 直接嵌入逆向设计流程，用梯度优化超表面结构。

RCWA 与平台此前介绍过的 FDTD 是互补的关系，它们的取舍可以总结为一张表：

表 1：RCWA 与 FDTD 的典型取舍（针对周期结构仿真场景）

对比项	RCWA（严格耦合波分析）	FDTD（时域有限差分）
求解域	频域，单波长单角度一次求解	时域，一次宽带激励覆盖一段频谱
基本思想	空间谐波展开 + 层内模式 + S 矩阵	空间网格上直接离散麦克斯韦方程组
典型对象	周期薄层结构：光栅、超表面单元、光子晶体板	任意形状：波导、谐振腔、非周期器件、瞬态问题
单点计算成本	通常较低，适合参数扫描	通常较高，但一次可得宽带响应
难点	周期边界与收敛阶数选择	网格分辨率、Courant 稳定性条件、边界吸收
适用边界	结构需可分层、激励为平面波；非周期问题需借助 PML 等扩展	大尺寸或极高 Q 值结构耗时显著

可以看到，RCWA 的”不能”其实也清晰：它天然面向周期性分层结构，对非周期器件需要额外技巧（例如在边界加上完美匹配层 PML 来模拟开放区域）；它对任意三维复杂形貌依赖阶梯近似，过于陡峭或精细的轮廓需要很多薄层；它的激励一般是平面波，局域光源或脉冲源不是它的强项。这些边界并不是缺陷，而是方法定位的一部分——正因为如此，实际工程中常见 RCWA 与 FDTD 配合使用：先用 RCWA 快速扫参数、缩小设计空间，再用 FDTD 对最终结构做更一般的验证。

六、意义与展望

RCWA 是连接”严格麦克斯韦求解”与”工程可用效率”的一座桥。它把周期结构的光学响应变成一个可以反复快速计算的矩阵问题，让光栅设计、超表面逆向优化和半导体光学量测从”靠经验试错”走向”靠仿真筛选”。从原理验证走向更广的实际应用，还有几个方向值得关注：更快的矢量公式与更好的预处理，让任意形状结构也能快速收敛；GPU 与自动微分结合的实现，使 RCWA 成为可微仿真器、直接服务于伴随法逆向设计；以及结合 PML、近场光源建模等扩展，把方法的适用对象从理想周期结构推向真实器件。可以预期，随着仿真工具的开源化和可微化的推进，RCWA 会在光电子设计与制造工艺优化中扮演越来越基础的角色。

参考资料

1. Moharam, M. G., & Gaylord, T. K. (1981). Rigorous coupled-wave analysis of planar-grating diffraction. Journal of the Optical Society of America, 71(7), 811–818. https://doi.org/10.1364/JOSA.71.000811
2. Moharam, M. G., Grann, E. B., Pommet, D. A., & Gaylord, T. K. (1995). Formulation for stable and efficient implementation of the rigorous coupled-wave analysis of binary gratings. Journal of the Optical Society of America A, 12(5), 1068–1076. https://doi.org/10.1364/JOSAA.12.001068
3. Li, L. (1997). New formulation of the Fourier modal method for crossed surface-relief gratings. Journal of the Optical Society of America A, 14(10), 2758–2767. https://doi.org/10.1364/JOSAA.14.002758
4. Liu, V., & Fan, S. (2012). S4: A free electromagnetic solver for layered periodic structures. Computer Physics Communications, 183(10), 2233–2244. https://doi.org/10.1016/j.cpc.2011.12.005
5. Rigorous coupled-wave analysis. Wikipedia. https://en.wikipedia.org/wiki/Rigorous_coupled-wave_analysis
6. Ansys Optics. RCWA Solver Introduction. https://optics.ansys.com/hc/en-us/articles/4414575008787-RCWA-Solver-Introduction
7. Facebook Research. FMMAX: Fourier Modal Method with JAX（README）. https://raw.githubusercontent.com/facebookresearch/fmmax/main/README.md
8. Liu, V., & Fan, S. S4: Stanford Stratified Structure Solver（GitHub）. https://github.com/morganblevins/S4
9. 逍遥科技. RCWA与FDTD：pMaxwell在PIC Studio平台上支持光电子器件仿真（讯石光通讯网）. http://www.iccsz.com/site/cn/News/2024/07/03/20240703011224601367.htm