工作流、例题与失败模式

从单胞到器件：PCSEL 仿真路径

读前准备#

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

若一开始就建立毫米尺度三维全波模型，算力首先消耗在尚未确认的层栈和单胞参数上；若只算周期单胞，又看不到有限边缘、电流扩展和热回写。更可靠的路径是让问题随证据逐级变复杂：纵向模先固定光学背景，周期单胞筛出候选模，有限阵列检查边缘与远场，最后才加入电注入、温升和容差。图 25.1给出这条路径；图中的 $\Delta T\lesssim\SI{50}{K}$ 与几何偏差 $\pm\SI{5}{nm}$ 只是[本书算例]门槛，不能直接移作其他器件的验收标准[Moodie et al., 2024] [Gautam et al., 2024]。

同一结构的正向求解与反向诊断#

表 25.1把本章的求解步骤与例题、故障症状对应起来。它的用途是定位“哪一层物理尚未进入模型”，而不是规定唯一的软件流程。

正向计算时按 A--G 逐步增加物理；结果异常时从症状所在行反查边界、输入和模型层级。

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

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

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

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

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

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

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

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

建模前先核对参数语义#

跨求解器最常见的差异并非离散误差，而是输入根本没有描述同一结构：冷态折射率被传入热态模型、孔半径与直径混用、RCWA 单胞和有限阵列来自不同版几何，或坐标原点在电学模型中发生偏移。开始求解前，应为几何、材料、层栈和电热参数同时注明单位、来源、适用温度/波段与不确定度；缺少其中任一项的参数只能用于探索计算。

参数表应由各模型共同读取或由自动检查证明等价。真正需要保存的是物理输入及其来源，而不是一套固定目录名称；商业软件工程文件、脚本和导出数据可按具体复现环境组织。

阶段 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。最小过程通常是：

1. 把候选单胞复制成有限阵列，并加入真实上下层栈；

2. 在开放方向使用 PML，在必要时利用对称边界降低计算量；

3. 用脉冲源、偶极源或与目标模匹配的激励源激发结构；

4. 通过时域衰减、频谱、场分布和近远场变换提取复频率、$Q$、远场和主模排序；

5. 对阵列尺寸、PML 厚度、网格和仿真时间窗做收敛检查。

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

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

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

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

在 Lumerical 路径里，常见的分工是用 CHARGE 解漂移-扩散和泊松方程，用 HEAT 解热传导。最小过程可以概括为：

1. 依据真实外延层和电极定义半导体几何与掺杂；

2. 在 CHARGE 中指定接触边界、复合模型和工作电压/电流，求得 $J(\mathbf{r})$ 与 $N(\mathbf{r})$；

3. 依据 Joule 热和非辐射复合构造热源，送入 HEAT 得到 $T(\mathbf{r})$；

4. 把 $N$ 和 $T$ 场转写成折射率修正、增益修正和吸收修正，回写给光学阶段。

在 COMSOL 路径里，对应的是 Semiconductor 与 Heat Transfer 接口。几何通常可与光学模型共享，只是在网格、边界和研究类型上要做针对性修改。最小过程是：

1. 在 Semiconductor 接口中设置接触、掺杂、复合与偏置，解出电势和载流子分布；

2. 在 Heat Transfer 接口中加入 Joule 热与复合热源，求稳态或瞬态温度场；

3. 将温度和载流子结果通过自定义材料关系回写到折射率和增益代理中；

4. 重新运行 Wave Optics，检查目标模与竞争模的净裕量是否改变。

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

阶段 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