数值方法
CWT、RCWA、FDTD、FEM 的比较与互证
读前准备
建议已掌握 第 10 章、第 11 章、第 14 章、第 15 章 中关于半解析模型、开放周期结构和全波有限器件求解的内容。
取同一个方格孔阵：CWT 把它压缩成少数耦合通道，RCWA 求无限周期单胞的端口散射，FDTD 与 FEM 则把有限边缘和复杂层栈带回模型。四种计算即使都正确，输出也不会逐项相同，因为它们解的不是同一个边值问题。本章按待回答的问题比较方法，并给出可交叉核对的公共量，而不按软件名称排高低。
为什么单一方法几乎注定不够
PCSEL 同时是三种对象：它是周期开放电磁系统，是半导体有源腔，也是带真实注入与热分布的器件。任何一种方法，只要天然擅长其中一层，就通常会在另外两层变弱。
CWT 擅长把复杂物理压缩成可解释的耦合通道，因此适合建立直觉、做灵敏度分析和前期筛选；但它的代价是近似较强，不能天然承担高保真几何与开放边界的全部负担。
RCWA 擅长周期开放单胞，特别适合研究辐射通道、偏振和开放共振；但它默认的是无限周期单胞，不天然处理有限器件边缘。
FDTD 和 FEM 擅长真实有限器件全波分析，能把边缘、复杂层栈和真实几何带回来；但它们计算代价更高，也通常不像 CWT 那样自动给出清晰的机制分解。
于是，单一方法不够是因为 PCSEL 的问题本来就不是单层问题，而非因为“每种方法都不行”。方法必须分层合作，才和物理对象的层次相匹配。
先按“问题类型”而不是“方法名字”来比较
最有用的比较方式，是先问“我现在究竟在回答什么问题”，而非问“哪个方法更准”。如果问题是候选模目录、主导耦合通道和几何灵敏度，那么 CWT 最合适；如果问题是开放周期单胞的辐射、反射谱和光学阈值增益趋势，RCWA 最合适；如果问题是有限器件边缘、近远场、真实 和复杂几何，FDTD/FEM 更合适；如果问题进一步变成阈值电流、热稳定性和模式切换，则必须把这些光学结果交给器件模型。
也就是说，方法比较的第一原则是“按问题选择方法”，而不是“按方法反推问题”。这条原则一旦建立，许多常见争论其实可被重新表述。比如，“RCWA 准还是 FDTD 准”这种问题本身就不完整，因为两者默认解决的问题层级并不一样。
| 方法 | 典型自由度 | 典型对象尺度 | 单次求解量级 |
|---|---|---|---|
| PWEM | -- | 二维单胞、被动无限周期 | 秒级到分钟级 |
| CWT | -- | 少数主导波与参数扫描 | 秒级 |
| RCWA | -- | 三维层状周期单胞 | 秒级到分钟级 |
| 3D FDTD | -- | 数十周期有限孔阵、宽频瞬态 | 小时级到天级 |
| 3D FEM | -- | 有限器件、本征模或窄频响应 | 小时级到天级 |
表 16.1以 16 核 CPU、64 GB 内存和达到本章误差门槛的最小离散规模为参考。硬件、求解器实现与收敛要求变化时，耗时可跨越一个数量级；表中数值只说明“候选筛选”和“三维有限结构验证”的成本差异[Moharam et al., 1995] [Li, 1996]。
图 16.1是定性索引，多边形不代表实测评分。跨软件定量比较需要固定几何、边界、误差门槛、版本和硬件。
四种方法各自最擅长什么
为了把分工说得更清楚，可以用一句话概括四种方法的适用场景。
CWT 最擅长回答：“为什么这类几何扰动会主要改这个耦合通道，而不是另一个？”
RCWA 最擅长回答：“这个周期单胞在开放边界下会怎样把能量耦合到各个辐射通道？”
FDTD 最擅长回答：“这个有限器件在宽频和瞬态条件下，实际会怎样建立、衰减和辐射？”
FEM 最擅长回答：“这个复杂几何、复杂层栈的有限器件在某个频率或某支模式附近，精确场分布和复本征值是什么？”
如果把四者的适用场景混在一起用，就会出现常见误判。例如：用 CWT 试图给出精确远场角宽，用 RCWA 试图给出有限器件边缘模泄漏，用 FDTD 的一组冷腔结果直接断言器件阈值电流。这些都是问题与方法没有对齐，而非方法本身的问题。
| 方法 | 不应单独承担的结论 | 原因 |
|---|---|---|
| CWT | 精确远场角宽、有限阵列边缘泄漏、器件级阈值电流 | 基底维数过低，且有限尺寸与电热回写未显式进入 |
| PWEM | 开放辐射损耗、真实有限器件排序、阈值结论 | 典型形式是无损被动无限周期本征问题 |
| RCWA | 有限孔阵边缘主导问题、横向包络与边缘吸收主导问题 | 单胞周期模型本身不含真实有限尺寸终止 |
| FDTD/FEM 冷腔 | 电注入阈值电流、阈上稳态、热稳定工作窗口 | 若无电热回写，只能给冷腔或小信号光学层级结论 |
比较方法时，应比较哪些量
多方法互证并不意味着所有数值都要逐位相同。不同方法的基底、边界和离散方式不同，过分追求“完全一致”往往会掩盖更重要的信息。更有意义的是比较那些应当在方法变化下保持稳定的物理量或趋势。
第一类是候选模的排序。例如目标频段内究竟是哪几支模式最值得关注，这种排序不应因方法切换而完全翻转。
第二类是对称性与偏振标签。虽然数值细节可能不同，但目标模属于哪类对称性表示、主要偏向哪种偏振，应该能在方法间对得上。
第三类是参数微扰趋势。若增大孔半径在 CWT 中预测会增强某种辐射通道，那么 RCWA 与全波方法至少应在趋势上给出相容结论。
第四类是损耗来源分解。即使绝对 不完全一致，也应能判断损耗主要来自向上辐射、向下辐射、边缘泄漏还是掺杂吸收重叠。
换句话说，互证的重点是“物理解释链条在不同分辨率下是否一致”，而非“每个小数点都一样”。
跨方法比较时，最危险的是层级混淆，而非少量数值偏差。若一种方法给的是无限周期被动结果，而另一种方法给的是有限器件冷腔结果，再拿二者的绝对阈值优劣直接比较，本身就没有严格可比性。先对齐问题层级，再比较结果，才有意义。
数值误差从哪来：互证前的最后一项功课
每种方法都有自己的主导误差源与收敛行为：PWEM 的平面波截断与有效折射率近似、CWT 的模型降阶、RCWA 的 Fourier 收敛与 Gibbs 现象、FDTD 的时域误差累积、FEM 的网格与弱形式实现。逐方法的误差源清单、三层收敛门槛口径（探索级/本书算例/发布级）、不确定度传递与报告规范，统一收入 第 30 章 的数值不确定度部分（第 30.12.1 小节）。互证时只需带着一条纪律：任何进入比较的数字，都必须附有收敛证据和误差量级。
为什么互证必须按“由粗到细”的顺序组织
如果一开始就直接上三维全器件全波扫描，当然不是做不到，但通常代价极高，而且一旦结果不理想，很难知道问题究竟出在晶格、外延还是边界设置。相反，若先由粗到细组织工作流，每一层都能承担一个清楚的筛选任务，后续搜索空间会大幅缩小。
最自然的顺序通常是：先用 第 5 章 的能带阅读框架配合 第 13 章 的 PWEM 找到候选模，再用 CWT 识别主导耦合通道与参数方向，再用 RCWA 精修单胞开放散射与光学阈值趋势，再用 FDTD/FEM 把有限尺寸和真实外延层栈带回来，最后把这些结果交给电注入和热耦合模型。这条顺序并非人为规定，而是与问题层级天然匹配。
这样做的另一个好处是，任何阶段一旦出问题，都能更容易定位。例如，如果 CWT 和 RCWA 已经都显示某方案目标模并无明显优势，那就没必要花大代价做全器件验证；反过来，如果 RCWA 很乐观、但全波有限器件结果很差，问题通常就出在边缘、有限尺寸或纵向层栈细节，而不是单胞本身。
一条适合 PCSEL 的推荐工作流
把前面的思路压缩成一条实际可执行的工作流，可以概括为五步。
第一步，建立候选模地图。使用 第 5 章 的能带概念和 第 13 章 的 PWEM 计算，在目标波长附近找出 Gamma 点模族、简并关系和对称性标签。
第二步，做机理级预设计。使用 CWT 或更高阶 CWT，把注意力集中到少数关键耦合通道，判断哪些几何扰动最可能改变模式分裂、偏振选择和辐射窗口。
第三步，做单胞开放精修。使用 RCWA，围绕前两步锁定的候选模区间，提取被动谱、上下辐射分配和光学阈值增益趋势。
第四步，做有限器件全波验证。使用 FDTD/FEM，把有限孔阵边缘、真实层栈、窗口和边界条件带进模型，检查真实近场、远场、有限尺寸 和模式排序。
第五步，做器件级回写链。把第四步输出的模式重叠、损耗和导数信息交给增益模型、电学输运与热耦合模型，得到阈值电流、工作点漂移和高功率稳定性判断。
这条工作流最大的价值，不在于“把所有工具都用一遍”，而在于每一步都有清楚的淘汰与交接标准。
| 当前层 | 主要中间量 | 下一层如何使用 |
|---|---|---|
| PWEM | 候选频段、Gamma 点模目录、简并与对称性标签 | 交给 CWT 选择最小耦合子空间，交给 RCWA 锁定扫描频段 |
| CWT | 主导耦合通道、模式分裂方向、参数灵敏度符号 | 交给 RCWA 选重点参数和候选模，减少盲扫 |
| RCWA | 单胞开放谱、上下辐射分配、光学阈值代理、偏振趋势 | 交给 FDTD/FEM 构造有限阵列验证目标，交给器件模型提供辐射损耗输入 |
| FDTD/FEM | 有限器件近场/远场、有限尺寸 、场重叠与边缘泄漏 | 交给增益、电学和热模型作为器件级输入 |
| 电热光模型 | 阈值电流、温升、工作点漂移、稳定工作窗口 | 交给实验回写链做回验与设计更新 |
| 层级 | 继续前进前至少要确认的事 | 若未通过应如何处理 |
|---|---|---|
| PWEM | 候选模频段、简并标签和带图读取无歧义 | 回到单胞几何与归一化定义 |
| CWT | 参数变化方向与主导耦合通道已识别 | 先缩小参数空间，不要急于全波 |
| RCWA | 候选模开放谱与辐射排序稳定，功率守恒通过 | 回查傅里叶截断、材料色散和边界设定 |
| FDTD/FEM | 阵列尺寸、PML、监视面和网格收敛，模式身份稳定 | 回查有限尺寸与开放边界，不要越级宣称器件优劣 |
| 电热光 | 注入、温升和工作点回写后目标模仍占优 | 不能输出最终器件窗口，应回到光学或注入分布层重设方案 |
2024--2026 前沿补充：代理模型与 AI 设计器如何接入互证链
前沿工作显示，PCSEL 的快速设计工具正在从“低阶解析近似”扩展到“解析模型 + 数据驱动代理”的混合范式：一方面，快速工具评估研究给出不同近似器在精度/速度上的可操作对照[Lang et al., 2025]；另一方面，基于 Transformer 的代理模型已能在复杂晶格参数空间中实现高吞吐预测[Qi et al., 2025]。 这并不改变本章的主结论，反而强化了“层级互证”的必要性：代理模型可以提前压缩搜索空间，但不能绕过物理保真层验收。
本书建议把代理模型接入点固定在“阶段 B 与阶段 C 之间”：
先用代理模型做粗筛（候选区间、非可行区排除）；
再用 RCWA/FDTD/FEM 做物理验收（模式身份、边界条件、收敛与误差）；
最后用器件级电热模型验证工作窗口是否仍成立。
若缺少第二步和第三步，代理模型给出的“最优设计”只能算假设生成器，不应直接进入发布级结论。
互证不是礼仪，而是发现错误的主要手段
研究里常见的一种误区，是把多方法互证理解成“为了让论文更好看，多放几张图”。其实互证的主要价值是发现错误。原因很简单：不同方法对不同误差源敏感。若 RCWA 和 FDTD 对同一趋势意见一致，通常说明单胞开放物理和有限器件物理在该问题上相容；若它们结论冲突，就说明至少有一层物理被漏掉或被误设。
举例来说，若 CWT 预测某种椭圆孔扰动会增强偏振分裂，而 RCWA 也显示对应偏振通道的辐射损耗排序拉开，但 FEM 在真实有限外延层栈里却显示两支模式仍接近简并，那么问题很可能不在面内晶格，而在纵向层栈和边界。反之，若 RCWA 与 FEM 一致乐观，而器件级热电模型仍显示工作点下模式切换频繁，那么问题多半不在冷腔光学，而在注入与热反馈。
互证的作用不只是“确认”，还包括“定位”。
假设某方案在 RCWA 中显示目标模的光学阈值增益最低，于是看起来很理想；但 FEM 有限器件分析发现该模式在边缘附近场很强，导致边缘吸收一并增强；再进一步，热-电-光耦合分析又显示电流拥挤把增益集中在另一片区域。最终器件先起振的模式可能并非 RCWA 中最优的那支。这是它只负责单胞开放层级，后面两层信息把排序改写了，而非 RCWA 错。
这个例子在于多方法链存在的意义：每一层都在回答自己的问题，任何一层都不应擅自越权宣布最终器件结论。
什么时候该停在某一层，不再继续加方法
并非每个课题都需要把所有方法全部用满。若你的目标只是理解某种晶格扰动的基本耦合机制，CWT 加 RCWA 可能已经足够；若你的目标是论文级有限器件远场和偏振验证，通常还需要 FEM 或 FDTD；若你的目标是器件级阈值电流与高功率稳定工作区，那就必须继续进入后续的外延、电学和热学模型。
判断是否需要继续升级的一个实用标准是：当前层输出的信息，是否足以支持下一步设计决策。如果还不能，就必须升级；如果已经足够支持当前问题，不应为了“形式上更完整”而无条件加重模型。
本章结论
CWT 给出耦合机制，RCWA 给出开放周期单胞的端口响应，FDTD/FEM 给出有限结构的场、复频率与远场。互证应比较同一物理量的模式身份、排序和参数趋势；阈值电流与热稳定性还需把这些光学输出传入外延、增益、输运和热模型。
练习题
[基础] 解释为什么“哪种方法更准”这个问题在没有说明问题层级时往往没有意义。
[基础] 设计一条适用于大面积单模 PCSEL 的方法工作流，并说明每一步的淘汰标准。
[进阶] 给出一个 RCWA 与 FDTD 结论可能冲突的场景，并说明冲突背后最可能缺失的物理是什么。
[综合] 说明为什么互证时更应比较模式排序、对称性标签和参数趋势，而不是只比较单一绝对数值。
延伸阅读
下一章开始正式进入外延层的光学结构。后续章节将说明，前面光学方法的输出量如何被转换成外延设计中的纵向重叠、损耗重叠和上下辐射不对称。
若想补充方法基础，可回看 [Joannopoulos et al., 2008] [Taflove & Hagness, 2005] 中关于本征法、散射法与时域法的相关章节。
若想检验自己是否理解了“方法层输出如何喂给后续模块”，可带着本章的分工框架回读 第 9 章、第 23 章：前者依赖你正确分辨开放辐射通道，后者依赖你正确继承模式损耗和阈值排序。