第四章：RCWA方法入门——严格耦合波分析原理、傅里叶展开、特征值求解、收敛性分析
各位好，我是老张。做MicroLED光学仿真这些年，RCWA（严格耦合波分析）是我用得最多的频域方法之一。说实话，刚入行那会儿我也被这个名字唬住了——什么「严格耦合波」，听着就头大。但用久了你会发现，它本质上就是个处理周期性结构电磁场问题的数学工具。今天咱们就把它拆开揉碎了讲清楚。
4.1 为什么需要RCWA？
做MicroLED的人都知道，芯片表面有周期性光栅结构——比如PSS（图形化蓝宝石衬底）、光子晶体、或者纳米柱阵列。这些结构的周期往往和波长在同一量级，传统的几何光学（比如光线追迹）根本算不准。而FDTD虽然能算，但计算量太大，尤其当你需要扫描参数时，跑一次就要半天。
RCWA的优势在哪？说白了，它把空间域的问题转换到频域，用傅里叶级数展开电磁场，然后通过求解特征值问题得到衍射效率。对于周期性结构，它比FDTD快得多——我做过对比，同样一个光栅结构，RCWA比FDTD快至少一个数量级。
核心思想：将电磁场和介电常数在周期方向上做傅里叶展开，把麦克斯韦方程组转化为一组代数特征值问题。你不需要关心每个网格点的场值，只需要关心各个衍射级次的振幅。
4.2 傅里叶展开——把连续问题离散化
RCWA的第一步，就是把介电常数和电磁场展开成傅里叶级数。假设你的光栅沿x方向周期为Λ，那么介电常数ε(x)可以写成：
ε(x) = Σ εₙ · exp(j·2π·n·x/Λ)这里n从负无穷到正无穷。当然实际计算时我们只能取有限项——这就是截断阶数，通常用N表示。N越大，精度越高，但计算量也越大。
我记得第一次做RCWA仿真时，取了个N=5，结果算出来的衍射效率和实验对不上。后来发现是傅里叶展开的项数不够——对于高折射率对比度的结构（比如GaN和空气界面），需要更多的傅里叶项才能准确描述介电常数的突变。
我的经验：对于MicroLED常见的GaN/空气光栅，建议N至少取20以上。如果结构中有金属（比如Ag反射镜），N可能需要50甚至更高。你可以先跑一个粗网格看看趋势，再逐步增加N验证收敛性。
4.3 特征值求解——RCWA的核心步骤
傅里叶展开之后，麦克斯韦方程组就变成了一个矩阵特征值问题。具体来说，我们需要求解：
∂²U/∂z² = A · U其中U是电场或磁场的切向分量，A是一个由傅里叶系数构成的矩阵。求解这个方程，本质上就是求矩阵A的特征值和特征向量。
嗯，这里要注意：A矩阵的维度是(2N+1)×(2N+1)。如果你取N=20，矩阵就是41×41，求解特征值对现代计算机来说小菜一碟。但如果你取N=200，矩阵就是401×401，计算量就上来了。
我个人习惯用MATLAB或Python的numpy.linalg.eig来求解。不过要提醒一句：当结构中有损耗材料（比如金属）时，特征值可能是复数，对应的传播常数也有虚部——这代表光在传播过程中有衰减。
避坑指南：我曾经在计算金属光栅时，发现特征值求解结果不稳定，衍射效率出现非物理的负值。后来排查发现是矩阵条件数太大——介电常数实部和虚部相差好几个数量级。解决办法是使用双精度浮点数，或者改用更稳定的广义特征值求解器。
4.4 收敛性分析——你的结果可信吗？
RCWA的收敛性是个大问题。说白了，你取的傅里叶项数N不同，结果可能天差地别。怎么判断你的结果收敛了？
我的做法是：
- 逐步增加N：从N=10开始，每次翻倍，直到N=80或100。观察衍射效率的变化。
- 看相对误差：当N增加一倍时，衍射效率的变化小于0.1%，我就认为收敛了。
- 检查能量守恒：所有衍射级次的反射率和透射率之和应该等于1（无损耗材料）。如果偏差超过1%，说明N不够。
举个例子，我之前做MicroLED的PSS结构仿真，N=20时能量守恒误差约3%，N=40时降到0.5%，N=80时只有0.05%。最后我选了N=60作为折中——精度够了，计算时间也能接受。
收敛性判断标准：
| N值 | 能量守恒误差 | 建议 |
|---|---|---|
| 10-20 | > 2% | 仅用于趋势分析 |
| 30-50 | 0.5% - 2% | 可用于定性分析 |
| 60-100 | < 0.5% | 推荐用于定量设计 |
| > 100 | < 0.1% | 高精度研究使用 |
4.5 知识体系框架
下面这张图是我自己整理的RCWA方法知识结构，你可以把它当作学习路线图：
4.6 一个简单的RCWA计算流程
最后，我给大家梳理一下RCWA的完整计算流程。你写代码时可以照着这个框架来：
- 定义结构参数：周期Λ、光栅厚度d、材料折射率n1/n2、入射角θ、波长λ
- 设置截断阶数N：决定傅里叶展开的项数
- 构建介电常数傅里叶矩阵：计算εₙ系数，组装成Toeplitz矩阵
- 构建特征值矩阵A：根据TE/TM偏振，组装对应的矩阵
- 求解特征值：得到传播常数γₘ和对应的特征向量
- 匹配边界条件：在光栅上下界面，用电磁场连续条件求解未知系数
- 计算衍射效率：得到各衍射级次的反射率和透射率
- 收敛性检查：增加N，重复计算，直到结果稳定
嗯，这个流程看起来步骤多，但实际写代码也就一两百行。我最早用MATLAB写RCWA时，花了三天调试边界条件匹配那部分——总是算出来能量不守恒。后来发现是界面处的场匹配方程写反了符号。这种错误，你光看公式看不出来，必须跑个简单例子验证。
验证小技巧：先用一个最简单的例子——均匀介质平板（没有光栅，就是一层膜）。RCWA应该能退化成薄膜干涉的解析解。如果这个例子都算不对，那肯定是代码有bug。
好了，关于RCWA的入门内容就讲到这里。记住，这个方法的核心就三件事：傅里叶展开、特征值求解、收敛性分析。把这三件事搞明白，你就能处理大部分MicroLED的周期性结构仿真了。