工作流、例题与失败模式

一套实用的 PCSEL 仿真工作流

读完本章，你应能#

• 把前面二十余章建立的概念、方程和近似组织成一条可执行的研究工作流。

• 明确不同阶段该计算什么、向下一阶段交付什么，以及哪些结论在该阶段绝不能提前宣称。

• 在求解器无关的原则下，理解 Lumerical 与 COMSOL 分别更适合承担哪些计算，并掌握最小可执行建模过程。

读前准备#

建议已掌握 第 16 章、第 17 章、第 19 章、第 20 章、第 22 章、第 23 章、第 24 章 中关于方法分工、外延层、电学、热学、鲁棒性分析、动态稳定性和性能报告的内容。

本章结构#

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

到这里，全书已经把 PCSEL 的主要物理链条铺开了。剩下的问题不再是“某个公式怎么推”，而是“研究时应该先做什么、后做什么、哪些结果能信、哪些结果还只是候选”。本章把这些环节整理成可执行的工作流。我们先建立分层结论和数据账本，再逐阶段说明纵向模检查、周期单胞筛选、有限器件验证、电热建模、热-电-光回写以及鲁棒性分析的输入输出，随后专门讨论两条常见实现路径：以 Lumerical 为主的光学优先链，以及以 COMSOL 为主的多物理场优先链。最后给出停机判据和结果交付模板。 图 25.1 把指标、结构、能带、阈值、远场、热与容差串成七步闭环；其中 $\Delta T\lesssim\SI{50}{K}$ 和几何偏差 $\pm\SI{5}{nm}$ 仅是[本书算例]筛选门槛，发布级判据仍须绑定波段、封装散热与工艺统计重新定义[Moodie et al., 2024] [Gautam et al., 2024]。

本章如何与下一章例题和失败模式章配合使用#

如果只把本章当成“流程图章节”，读者容易知道名词，却不知道每一步真正长什么样；如果只看下一章例题，又可能误以为那只是某一个特定结构的演示，而不是一般工作流的落地版本。较稳妥的读法是把 第 25 章、第 26 章、第 27 章 连成一套使用：本章给出正向流程，第 26 章 把它落实成一次完整作答，而 第 27 章 则告诉你每一步最常见的失误会长成什么样。

因此，建议不要孤立阅读本章。较合适的顺序是：先读本章建立流程骨架，再在 第 26 章 中按表完成各阶段检查，最后用 第 27 章 反向审计整条链。

notebox FDTD 内存需求的快速估算。对三维 FDTD 仿真，内存占用可由网格数量粗略估算：

其中$N_x, N_y, N_z$是各方向的网格数，$N_\text{fields} \approx 6$--$12$（电场、磁场、辅助变量、PML 等）。对典型 PCSEL 仿真：

• 单胞周期结构：$100\times100\times50$网格 $\rightarrow$ 约 0.5--1 GB

• 有限阵列（$10\times10$晶格）：$500\times500\times100$网格 $\rightarrow$ 约 15--30 GB

• 大尺寸器件（$100\times100$晶格）：$5000\times5000\times100$网格 $\rightarrow$ 约 1.5--3 TB

因此，对有限阵列 PCSEL 的全波仿真，通常需要至少 32--64 GB 内存的工作站，或使用子网格、对称性、周期边界等技巧降低网格总量。 notebox

工作流的第一原则：先定义结论层级，再跑仿真#

仿真失败常常是因为在第一步就没有定义清楚自己打算得到哪一层结论，而非因为不会点界面操作。对 PCSEL 来说，最安全的做法是把结论分为四层：

1. L1 结构候选层：只谈候选模、对称性、被动辐射趋势和纵向重叠；

2. L2 光学阈值层：只谈开放结构中的阈值增益、偏振趋势、上下出光比例和有限尺寸模式排序；

3. L3 电学映射层：开始谈注入分布、J-V、串联电阻、阈值电流区间的来源；

4. L4 热-电-光工作层：开始谈工作窗口、稳定性、亮度、偏振稳定性与热漂移后的器件表现。

如果还停留在 L1 或 L2，就绝不能提前宣称阈值电流、线宽窗口或高温稳定工作。这个分层是为了防止求解器输出和物理结论错位，而非文风问题。

图 25.2 把这一纪律画成阶段门禁。每跨过一级，是必须增加一种新的物理证据，而非多跑一个软件名；缺少对应证据时，结论应停留在较低层级。

工作流的第二原则：先建数据账本，再建几何#

在任何求解器启动之前，先建立参数账本 $\mathcal{D}$。它至少应记录四类信息：几何参数、材料参数、外延层参数和电热参数。每个参数都至少带四个元信息：单位、来源、适用温度或频段、当前不确定度。没有元信息的参数，不应直接进入主结果模型。

这个步骤看似偏向项目管理，但直接影响后续结果的可复核性。因为 PCSEL 的仿真链跨越多个求解器，若不在一开始统一单位、坐标和参数来源，后面最容易发生的是参数错配，而非数值误差。例如，光学模型用的是冷态折射率，电热模型却默认了工作温度；或者 RCWA 的单胞孔径与 CHARGE 的电极开窗尺寸根本不是同一个版本。一个看似微小的版本错配，足以让整条工作流失去可复核性。

把“数据账本”落实成可执行的 manifest#

为了让这条原则不止停留在原则表述，建议把参数账本直接写成一个可版本控制的 manifest。最小字段至少应包括：参数唯一 ID、显示符号、文件字段名、单位、数值、来源、适用波长/温度、被哪些阶段使用，以及当前不确定度。一个最小 YAML 模板可以长成下面这样：

verbatim parameter: id: MAT-GAAS-N-001 name: n_GaAs symbol: n file_key: n_gaas unit: 1 value: 3.42 source: literature valid_wavelength_nm: [930, 990] valid_temperature_C: [20, 80] uncertainty: 0.02 used_in: - stage_A_waveguide - stage_B_unitcell - stage_C_finite_array version: v1.2 last_updated: 2026-03-01 verbatim

在项目实施中，这个 manifest 还应再加两列：propagates_to（它会影响哪些输出量）和 status（direct-measured / literature / fitted / inferred）。这样处理后，读者第一次独立搭建建模链时，就可以避免把文献冷态参数与实验反演的工作点参数混入同一张无版本表格。

最小命令行复现骨架：即使使用商业软件，也应把文件流做成终端原生#

本书不把 GUI 界面操作教程当主线，但这不等于 workflow 只能停在“软件名词表”。较稳妥的做法，是先把目录、参数文件、阶段脚本和输出制品固定成一条 terminal-native 骨架，再决定每一步内部到底由 Lumerical、COMSOL 还是其他求解器来完成。一个最小项目骨架可以写成：

verbatim pcsel_project/ params.yaml 01_layers/ 02_unitcell/ 03_finite_array/ 04_electrical_thermal/ 05_postprocess/ run_stage_A.py run_stage_B.py run_stage_C.py run_stage_D.py run_stage_E.py verbatim

其中每一步至少应有固定的输入、输出、sanity checks 与失败回扫项：

• Stage A：输入层栈与材料表，输出 $n_{\mathrm{eff}}$、$\Gamma_{\mathrm{QW}}$、$\Gamma_{\mathrm{loss}}$；sanity check 是模式主峰位置和量纲自洽。

• Stage B：输入单胞几何与层栈背景，输出候选模频段、被动损耗排序、上下辐射比例；sanity check 是 PWEM / RCWA 频段与对称性标签一致。

• Stage C：输入有限阵列几何，输出近场、远场、主模排序和尺寸收敛；sanity check 是尺寸/PML/网格继续加密后主结论不翻转。

• Stage D：输入外延、电极与接触参数，输出 $J(\rvec)$、$N(\rvec)$、$T(\rvec)$；sanity check 是边界条件、单位和功率守恒一致。

• Stage E：输入前四阶段结果，输出回写后的净模裕量、工作窗口和鲁棒性报告；sanity check 是结论层级不越权。

这套骨架的目的是强制读者把文件流、参数流和结论层级固定下来，而非承诺“整本书都提供开源全复现脚本”。只要这一步做对，后续无论内部调用商业软件还是自写脚本，研究链条都会清楚得多。

阶段 A：先做纵向模和层栈检查，而不是急着扫平面晶格#

PCSEL 的第一步不应直接扫单胞孔径，而应先确认外延层给出的纵向模式背景是否合理。也就是说，先解层栈中的纵向模、有效折射率、量子阱重叠和损耗重叠。若这一关没有通过，后续平面内晶格优化很可能建立在错误的纵向模型上。

在 Lumerical 路径里，这一阶段通常最适合用 MODE/FDE 之类的波导模求解器：建立二维横截面层栈，定义材料色散和损耗，求出目标偏振附近的导模，输出 $n_{\mathrm{eff}}$、场分布、$\Gamma_{\mathrm{QW}}$ 与与损耗层重叠积分。最小建模过程可以写成：

1. 画出外延层截面并赋予折射率、掺杂相关吸收代理；

2. 在工作波长附近做本征模求解；

3. 提取目标模的 $n_{\mathrm{eff}}$、纵向场峰值位置和重叠因子；

4. 对层厚与折射率做小扰动扫描，判断结果是否过于脆弱。

在 COMSOL 路径里，这一阶段通常用二维横截面的 Wave Optics 本征/模态分析来完成。最小过程是：建立二维层栈几何、指定材料、用适合导模的边界条件和本征求解得到传播常数，再做后处理积分得到重叠与损耗趋势。这里的核心是输出必须一致：都要产出一个可供下一阶段使用的纵向模式背景，而不仅仅是一张场图，而非软件界面差异。

阶段 B：周期单胞筛选，目标是排除错误结构而不是宣布最终器件#

有了可靠的纵向背景，才进入单胞筛选阶段。这里的任务是回答：在无限周期、开放结构、小信号框架下，哪些 Gamma 点候选模、哪些偏振状态和哪些辐射通道值得继续投入有限器件验证。

这一阶段最好再拆成两个子步骤。先按 第 5 章、第 13 章 的方法做一次被动无限周期候选模筛选：用 PWEM 建立候选频段、Gamma 点简并关系和主要平面波分量；再把已经缩小后的候选模窗口交给 RCWA 或周期频域模型，去检查开放辐射、上下出光比例和偏振趋势。这样做的好处是，RCWA 不必在过大的参数空间里盲扫，而是围绕已经识别出的候选模做精修。

Lumerical 路径中，若目标是快速扫描单胞参数和辐射分配，优先使用 RCWA；若需要更灵活地看时域响应或近场细节，也可用周期边界下的 FDTD，但成本通常更高。一个最小 RCWA 过程可以概括为：

1. 建立单胞几何、上/下介质和层栈；

2. 设定法向附近的入射波或端口条件，并扫描波长、孔径、刻蚀深度和微扰参数；

3. 读取反射、透射、衍射效率、局域场增强和上下辐射分配；

4. 用收敛的傅里叶阶数确认候选模排序与偏振趋势。

COMSOL 路径中，这一步通常用三维单胞频域模型完成：在横向使用相位周期边界，在上下使用端口或开放边界，并通过频率扫或参数扫读取散射量。若要说得再具体一点，就是“用 Wave Optics 的频域单胞模型来做 RCWA 所擅长的那类周期散射判断”，只不过实现框架是以有限元频域周期散射问题出现，而非以 RCWA 名义出现。

这一阶段的交付物应当是候选模目录，而不是“最终器件最低阈值”的宣称。目录中至少应包含：候选频段、PWEM 给出的 Gamma 点模族和主要傅里叶分量、目标/竞争模的被动损耗排序、偏振趋势、上下出光比例和对工艺扰动的敏感性。

阶段 C：有限器件全波验证，检查无限周期假设到底失效到什么程度#

通过单胞筛选的结构，下一步就要离开无限周期近似，进入有限孔阵模型。这里的任务不再是看单胞散射特征，而是看真实有限阵列中模式是否仍保持正确排序、边缘是否引入新的泄漏通道、远场是否还能维持目标主瓣和偏振状态。

在 Lumerical 路径里，这一步最常见的是 FDTD。最小过程通常是： enumerate 把候选单胞复制成有限阵列，并加入真实上下层栈； 在开放方向使用 PML，在必要时利用对称边界降低计算量； 用脉冲源、偶极源或与目标模匹配的激励源激发结构； 通过时域衰减、频谱、场分布和近远场变换提取复频率、$Q$、远场和主模排序； 对阵列尺寸、PML 厚度、网格和仿真时间窗做收敛检查。 enumerate

在 COMSOL 路径里，这一步常用三维 Wave Optics 的频域或本征频率研究来做。若重点是共振位置和模式形状，可用本征频率研究并在开放边界加 PML；若重点是给定频率下的响应和远场，则更适合频域研究。需要强调的是，COMSOL 的有限元框架在处理复杂材料依赖和后续多物理场联动时很灵活，但三维有限孔阵的计算代价往往不小，因此必须严格控制模型尺寸和收敛验证。

这一阶段的输出至少应包括：有限阵列主模与竞争模的排序、目标模近场与远场、边缘敏感性、以及“无限周期结论在有限阵列中是否仍成立”的判断。

阶段 D：电注入与热建模，开始把阈值增益映射成阈值电流区间#

如果研究目标是器件而不是纯光学结构，那么到这里就必须进入电学与热学。任务不再是找模，而是问：电流怎样进入器件、载流子密度怎样分布、哪里掉压、哪里发热、这些分布怎样改写模增益与模损耗。

在 Lumerical 路径里，常见的分工是用 CHARGE 解漂移-扩散和泊松方程，用 HEAT 解热传导。最小过程可以概括为： enumerate 依据真实外延层和电极定义半导体几何与掺杂； 在 CHARGE 中指定接触边界、复合模型和工作电压/电流，求得 $J(\mathbf{r})$ 与 $N(\mathbf{r})$； 依据 Joule 热和非辐射复合构造热源，送入 HEAT 得到 $T(\mathbf{r})$； 把 $N$ 和 $T$ 场转写成折射率修正、增益修正和吸收修正，回写给光学阶段。 enumerate

在 COMSOL 路径里，对应的是 Semiconductor 与 Heat Transfer 接口。几何通常可与光学模型共享，只是在网格、边界和研究类型上要做针对性修改。最小过程是： enumerate 在 Semiconductor 接口中设置接触、掺杂、复合与偏置，解出电势和载流子分布； 在 Heat Transfer 接口中加入 Joule 热与复合热源，求稳态或瞬态温度场； 将温度和载流子结果通过自定义材料关系回写到折射率和增益代理中； 重新运行 Wave Optics，检查目标模与竞争模的净裕量是否改变。 enumerate

这里必须再次提醒：除非你真的完成了这个映射链，否则不能把“阈值增益估计值”写成“阈值电流”。

阶段 E：热-电-光迭代，不追求一轮算完，而追求工作窗口收敛#

初次建模时，容易把多物理场耦合理解为“把所有接口同时打开，一次求解结束”。对 PCSEL 来说，这种理解往往既低效也不稳健。较实用的策略通常是分阶段迭代：光学给出模式和重叠，电学给出注入分布，热学给出温度场，再把这些结果回写到光学模型，重复直到目标指标收敛。

工作流应检查的是工作窗口是否稳定，而非某一次迭代里的最佳单点。一个结构若只有在第一次冷腔迭代中表现良好，而一旦写回温升就发生模排序翻转，它就不能算成功设计。相反，一个方案即使冷腔 $Q$ 不是最高，若热回写后仍保持目标模的净裕量为正，它在器件层面上更有价值。

阶段 F：鲁棒性验证，检查工作窗口是否同时具有工艺窗口#

若阶段 A 到阶段 E 用于判断“名义器件在当前模型下可行”，那么阶段 F 要回答的是“这是不是一个能加工出来、能批量重复出来的器件方案”。也就是说，在名义工作窗口已经建立之后，还应把 第 22 章 中的 tolerance analysis 显式纳入流程。最小动作包括：定义关键误差参数、先做局部灵敏度筛查、对高敏感参数做联合抽样，并用通过率而不是单一样本图来报告结果。

这个阶段之所以不能省略，是因为很多设计在名义值上成功，却在孔半径、刻蚀深度或孔径错位上异常脆弱。若不把阶段 F 写进工作流，所谓“可行设计”往往只是在数值环境里可行，而不是在工艺环境里可行。

阶段 G：实验回写链，不让模型停在“可运行”而要走到“可校准”#

可用于研究和器件开发的工作流，还必须再多走一步：把模型接回实验。最小回写链应写成“测量 $\rightarrow$ 参数提取 $\rightarrow$ 分层校准 $\rightarrow$ 预测 $\rightarrow$ 回验”。这里的关键是把它当成模型可用性的约束，而非把实验当成结果展示。

具体做法通常是：先用截面、几何、电学和热学的直接测量值约束几何与边界条件；再用 I--V、L--I、热阻或谱漂移数据分别校准电学和热学子模型；最后才把这些中间结果回写给光学与增益模型，去预测未参与拟合的阈值漂移、远场变化或偏振窗口。只有当“未参与校准的观测量”也能被解释时，模型才具备受检验的外推能力。

这一阶段的更细致写法，会在 第 29 章 中系统展开。本章只强调一个纪律：不要把所有实验量都同时拿来拟合同一套自由参数。那样做容易得到“处处都像、处处都不可信”的过参数化模型。

阶段 H：代理模型与 AI 加速（新增）#

2024--2026 的一个明显趋势是：PCSEL 设计正在把“快速近似模型”进一步扩展为“近似模型 + 代理模型 + 全波验收”的混合流程[Lang et al., 2025] [Qi et al., 2025] [Xu et al., 2024]。这一步可降低参数扫描成本，但可靠性取决于是否保留物理验收门槛。

本书建议的最小门禁是： enumerate 数据来源可追溯：训练集必须标注生成器（CWT/RCWA/FDTD）、参数范围与版本； 外推边界可声明：仅在训练分布内做优化建议，超出分布必须回退到物理求解； 关键结论可复核：所有准备进入发布级结论的点，必须用 RCWA/FDTD/FEM 复算； 器件级回写链不跳步：代理模型不能替代电热回写与实验校准阶段。 enumerate 代理模型应定位为筛选加速层，而不是真值层。它在流程中的合理位置是阶段 B/C 的前置筛选，而不是阶段 G 的替代。

Lumerical 与 COMSOL 的两条常见实现路径#

前面分阶段讲的是物理链，下面把它压缩成两条常用工程路径，供读者实际启动项目时参考。需要强调，这两条路径都是推荐流程，不是唯一流程。

路径一：以 Lumerical 为主的光学优先链#

这条路径适合“先把单胞和有限阵列光学迅速扫清，再把最有希望的结构接到器件级电热模型”这一工作风格。典型顺序是： enumerate MODE/FDE：检查纵向模、$n_{\mathrm{eff}}$、$\Gamma_{\mathrm{QW}}$ 与损耗重叠； RCWA：扫单胞、提取候选模、上下出光和偏振趋势； FDTD：有限阵列验证、提取远场和边缘敏感性； CHARGE：求注入分布和串联电阻来源； HEAT：求温升和热热点； 脚本化回写：把 $N(\mathbf{r})$ 与 $T(\mathbf{r})$ 回写给光学模型，更新模增益和损耗排序。 enumerate 这条路径的特点是光学前端筛选较快，但多物理场耦合依赖明确的接口管理和脚本化回写。

路径二：以 COMSOL 为主的多物理场优先链#

这条路径适合“在统一几何和 PDE 框架下做电-热-光之间的参数回写和灵敏度分析”这一工作风格。典型顺序是： enumerate 二维或简化三维 Wave Optics：建立纵向模背景； 三维单胞 Wave Optics：做周期单胞散射与偏振筛选； 三维有限阵列 Wave Optics：做有限器件近场、远场与本征模验证； Semiconductor：求注入和载流子分布； Heat Transfer：求温度场与热热点； 参数回写与分阶段迭代：更新材料参数、重跑光学、检查工作窗口。 enumerate 这条路径的优点是多物理场接口统一、参数回写自然，缺点是大尺寸三维光学问题的计算成本往往更高，对网格和求解设置更敏感。

若研究目标是快速筛选 50 组单胞设计，再从中挑 2 组做器件级深入验证，通常更适合先走 Lumerical 光学优先链；若研究目标是深入分析一个给定外延与电极方案中，温升如何改变模式竞争，则更适合走 COMSOL 多物理场优先链。这里的“更适合”是工程判断，不是排他性规则。

何时停止迭代：停机判据比“再跑一轮看看”更重要#

完整工作流容易陷入无止境试错，因此必须在开工前就设定停机判据。最有用的停机条件通常有三类： enumerate 数值收敛停机：网格、时间窗、PML、傅里叶阶数继续加密时，关键输出变化已低于预设阈值； 物理一致性停机：功率守恒、参数量纲、层级结论和趋势解释均已自洽； 目标窗口停机：在目标工作区间内，目标模对竞争模的净裕量、远场和温升都满足要求。 enumerate 如果这三条没有满足，即使结果图平滑，也不应进入最终结论。

本章到此给出的是正向流程骨架与停机判据。当流程走到需要把仿真数字与实验测量放在同一张桌上时，完整的仿真--实验对标协议（层级对齐五项前置检查、三级对标体系、关键物理量技术规范、不符时的诊断树与报告伦理）已收入 第 29 章 的证据链部分（第 29.12 节），建议在进入 第 26 章 的阶段 G 之前通读一遍。

回看本章#

一套实用的 PCSEL 仿真工作流必须以结论分层和参数账本为起点，然后依次经历纵向模检查、周期单胞筛选、有限器件验证、电学与热学建模，以及热-电-光回写。Lumerical 与 COMSOL 都能承载这条链，但常见工程分工不同：前者通常更擅长高效光学筛选，后者通常更擅长统一多物理场耦合。无论选哪条路线，应追求的都是一个通过收敛和物理一致性检查的工作窗口，而非“某一轮最优图片”。若把本章和后面两章一起看，本章给出的是正向流程骨架，第 26 章 给出的是一次完整走法，第 27 章 给出的则是反向诊断顺序。

练习题#

enumerate [基础] 为一个你感兴趣的 PCSEL 方案写出 L1 到 L4 的结论层级，并说明每一级允许说什么、不允许说什么。

[基础] 分别给出阶段 A 到阶段 E 的输入与输出变量，特别指出哪些量需要在不同求解器之间传递。

[进阶] 比较以 Lumerical 为主和以 COMSOL 为主的两条工作流，说明它们各自最可能在哪一步最耗时、最容易出错。

[综合] 为你的项目设计三条停机判据，并解释为什么这些判据能防止你过早下结论。 enumerate

延伸阅读#

itemize 下一章将把本章工作流落实到一个完整示例，给出参数、建模顺序、求解设置模板和可信度检查。阅读时建议对照 表 25.1，逐步核对每一阶段到底在例题中对应哪一步。

结合 第 27 章 阅读，可以把本章工作流直接转化为项目级排错流程；最有效的用法，是先按本章正向走，再按失败模式章反向复盘。 itemize