2. 多物理场仿真基础
异质集成光子芯片的仿真难点，不在于单个物理场算不准，而在于多个物理场彼此“牵动”。例如，一个硅基微环调制器在电学上靠载流子注入改变折射率，注入电流同时产生焦耳热，热又进一步改变折射率与谐振波长；而温度升高带来的热膨胀还会引入应力，应力通过光弹效应再次修改折射率。若只做电磁仿真，得到的谐振峰位置往往与实测相差数纳米，这正是多物理场耦合需要解决的问题。
一、多物理场耦合的基本概念
多物理场耦合，指两个或多个物理场通过共享的物理量相互影响。按耦合方向可分为单向耦合与双向耦合：单向耦合中，场A影响场B，但场B对场A的反馈可忽略，例如“热→光”在弱光功率下可视为单向；双向耦合则需迭代求解，例如高功率下的“光→热→光”自热效应。按耦合方式又可分为域耦合（同一求解域内多个场同时求解）与边界耦合（通过界面传递通量或位移）。
在异质集成光子中，最常见的耦合对包括：
- 电–光耦合：载流子色散效应（等离子体色散），用于调制器、开关。
- 热–光耦合：热光系数（硅约1.8×10⁻⁴ /K），用于热光调谐与热失谐分析。
- 力–光耦合：光弹效应，应力改变折射率张量，影响偏振与相位。
- 热–力耦合：热膨胀产生应力，异质键合界面尤为敏感。
- 量子–光耦合：量子点、色心与腔模的耦合，涉及自发辐射率与Purcell因子。
理解耦合的关键是明确“谁是驱动量、谁是响应量、耦合系数是什么”。例如热光耦合中，驱动量是温度变化ΔT，响应量是有效折射率变化Δn_eff，耦合系数即热光系数∂n/∂T。只有把耦合系数、边界条件和时间尺度写清楚，后续数值方法才有意义。
二、多物理场仿真流程
一个可复现的多物理场仿真流程通常包含六个环节，但并非线性一次通过，而是带反馈的迭代：
- 几何与材料定义：建立波导、电极、键合层、衬底的三维或二维截面模型，赋予各区域材料参数（折射率、热导率、杨氏模量、泊松比、热膨胀系数、光弹系数）。异质集成中尤其要注意III-V与Si之间的键合层厚度，它常是热阻与应力的集中区。
- 物理场选择与耦合设置：确定求解哪些场、耦合方向、耦合项来源。例如热–光耦合中，把热求解得到的T插值到光学网格上，再通过∂n/∂T更新折射率。
- 网格划分：光学场要求网格小于λ/(4n)以捕捉模式，热场网格可较粗，力场在界面附近需加密。多物理场常采用不同网格，通过插值算子传递数据。
- 边界条件与激励：光学端口设模式输入，电学设电压或电流，热学设热沉温度与对流系数，力学设固定约束与自由表面。
- 求解与迭代：单向耦合按顺序求解；双向耦合需在每步迭代中检查残差，直到Δn_eff或ΔT收敛。
- 后处理与验证：提取有效折射率、谐振波长、温度分布、应力张量，并与解析估计或实验数据对比。
三、常用数值方法概述
不同物理场适合的数值方法不同，选择时需权衡精度、计算量与边界处理能力。
| 方法 | 适用物理场 | 核心思想 | 典型工具/场景 |
|---|---|---|---|
| FEM（有限元法） | 电磁、热、力学 | 将求解域离散为单元，用变分或加权残差建立代数方程组 | COMSOL、ANSYS；适合复杂几何与多物理场耦合 |
| FDTD（时域有限差分） | 电磁 | 在Yee网格上交替更新E和H，显式时间推进 | Lumerical FDTD；适合宽带、色散、非线性光学 |
| BPM（光束传播法） | 电磁（缓变包络） | 沿传播方向步进求解傍轴波动方程 | 长波导、 taper、定向耦合器快速仿真 |
| 热传导有限差分/有限元 | 热 | 离散傅里叶导热方程，求解温度场 | 稳态/瞬态热分析，热沉与界面热阻 |
| 弹性力学有限元 | 力学 | 离散平衡方程与几何方程，求解位移与应力 | 热应力、键合应力、光弹分析 |
FEM的强项是处理不规则几何与多场耦合，但网格色散误差需注意；FDTD适合宽带瞬态，但窄带长距离仿真效率低；BPM计算极快，但忽略反射，不适合高折射率对比的尖锐结构。实际工程中常采用“FDTD/FEM算模式与局部场，BPM算长距离传输，热/力FEM算慢变场”的分工策略。
四、耦合求解的实践要点
多物理场仿真的常见陷阱是“网格不匹配导致插值误差”和“时间尺度差异导致迭代不收敛”。光学时间尺度为飞秒，热时间尺度为微秒到毫秒，若做瞬态耦合，需采用多时间步或准静态近似。另一个要点是材料参数的温度依赖性：硅的热光系数并非常数，在高温区会变化，若仿真温升超过50 K，应使用非线性热光模型。
最后，验证是仿真可信度的关键。建议至少做三件事：一是网格收敛性检查，将网格加密一倍，观察关键结果变化是否小于1%；二是能量守恒检查，热场中输入功率应等于传导、对流与辐射之和；三是与解析解或已发表实验数据对比。只有通过验证的仿真结果，才能用于指导异质集成光子器件的设计与优化。