1. 项目概述当神经网络遇见物理定律在光学和电磁学的前沿研究中我们常常面临一个两难困境一方面基于麦克斯韦方程组的全波电磁仿真如FDTD、FEM虽然精确但计算成本高昂一次仿真动辄数小时在需要遍历庞大设计空间的逆向设计任务中这几乎是一个无法承受的负担。另一方面深度神经网络DNN作为强大的通用函数逼近器理论上可以学习从器件几何参数到其电磁响应如散射矩阵S(ω)的复杂映射关系一旦训练完成其预测速度可比传统仿真快几个数量级。这听起来像是完美的解决方案对吧但现实往往骨感。我早期尝试用标准全连接网络FCN或卷积网络CNN直接预测超表面的透射谱时就踩过不少坑。最大的问题有两个数据饥渴和黑箱不可信。为了达到可接受的预测精度我们往往需要准备数万甚至数十万个高质量的仿真样本数据生成本身就是一场计算资源的消耗战。更令人头疼的是即便模型在测试集上表现良好其预测也可能在物理上不合理——比如出现能量不守恒透射率反射率吸收率1或违反因果律的频谱。这种模型就像一个记忆力超群但缺乏理解力的学生能复述例题却无法应对稍加变化的考题更无法解释其答案背后的“为什么”。物理信息机器学习Physics-Informed Machine Learning的出现为我们指明了另一条路。它的核心思想不是让神经网络从零开始“发明”物理而是将我们已经深刻理解的物理定律如能量守恒、因果律、对称性等作为“硬约束”或“软引导”整合到模型架构和训练过程中。这就像给这位学生一本写满定理和公式的教科书要求他的所有答案都必须符合这些基本原理。这样训练出的模型不仅数据效率更高需要的训练样本更少其预测结果也天然具备物理一致性并且学习到的中间参数往往具有明确的物理意义从而打开了模型的黑箱。今天要深入探讨的正是我们团队基于这一理念为电磁散射问题开发的一个通用框架基于准正规模Quasinormal Mode, QNM展开的物理信息神经网络我们称之为QNM-Net。它不是一个针对特定器件的定制化模型而是一个模块化的架构。其核心在于我们用神经网络去预测那些决定散射行为的、更本质的物理参数如谐振频率、模式振幅、背景散射然后通过一个严格基于QNM理论的、可微分的物理模型层将这些参数合成为最终的散射矩阵S(ω)。这样一来模型从“学习S(ω)的复杂图案”转变为“学习产生这些图案的物理原因”后者显然是一个更简单、更本质的任务。1.1 核心需求解析为什么是准正规模在深入架构之前我们必须先理解为什么选择准正规模作为物理模型的基石。在电磁散射理论中一个开放、有耗的谐振结构如光子晶体微腔、超表面单元的响应可以近似表示为一系列谐振模式的叠加。这些模式就是准正规模。每个QNM由一个复本征频率 $\tilde{\omega}m \omega_m i\gamma_m$ 描述其中实部 $\omega_m$ 是谐振频率虚部 $\gamma_m$通常为负代表模式的衰减率包括辐射损耗 $\gamma{r,m}$ 和非辐射损耗 $\gamma_{nr,m}$。QNM理论的美妙之处在于它为散射矩阵S(ω)与这些内在谐振模式之间建立了严格且普适的数学联系即QNM展开式。我们采用的是一种能保证能量守恒的近似形式$$ S(\omega) e^{i\omega \tau} \left[ C(\omega) D (i\omega - i\tilde{\Omega})^{-1} M^{-1} D^{\dagger} C(\omega) \right] e^{i\omega \tau} $$这个公式看起来复杂但我们可以拆解其物理含义$C(\omega)$: 一个缓变的背景散射矩阵代表非谐振的直接散射过程。$D$: 一个矩阵其列向量 $d_m$ 代表了第m个QNM在各个端口上的耦合振幅。这是连接内部模式与外部端口的关键。$\tilde{\Omega}$: 一个对角矩阵对角线元素就是各个QNM的复本征频率 $\tilde{\omega}_m$。$M$: 模式重叠矩阵用于保证展开的正交归一化。$\tau$: 端口延迟矩阵是一个实对角矩阵代表电磁波从端口到散射体中心的相位延迟。这个公式就是我们的“物理层”。它不是一个可训练的神经网络层而是一个确定的、符合物理定律的数学表达式。神经网络的职责变成了根据输入的器件设计参数去预测这个表达式中的所有物理参数$C(\omega)$, $\tau$, 以及每一个模式的 $\tilde{\omega}_m$ 和 $d_m$。这样做带来了几个根本性优势物理一致性保障只要QNM展开式本身满足能量守恒和因果律那么无论神经网络预测的参数值如何最终计算出的S(ω)也自动满足这些基本物理定律。这从根本上杜绝了非物理解的出现。数据效率飞跃神经网络不再需要记忆所有可能的、复杂的频谱形状只需要学习相对简单、平滑的物理参数随设计的变化关系。这极大地降低了学习任务的复杂度。可解释性突破模型的学习目标变成了具有明确物理意义的参数。例如我们可以直接检查网络预测的 $\tilde{\omega}_m$并与全波本征模求解器的结果进行对比从而验证模型是否真的“理解”了器件的物理本质。这为“AI for Science”中的知识发现提供了可能。逆向设计更直接在逆向设计中我们通常的目标是获得某个特定的频谱。在QNM-Net框架下这个目标可以转化为对特定谐振频率、耦合强度等物理参数的需求。优化过程可以直接在这些物理参数空间进行目标更明确路径更清晰。2. QNM-Net架构深度拆解模块化设计的艺术QNM-Net不是一个僵化的固定网络而是一个高度模块化的框架。这种设计哲学源于一个认识不同的电磁器件如对称的光子晶体板 vs. 非对称的自由形式超表面其物理约束和先验知识是不同的。一个好的物理信息模型应该能灵活地吸收这些知识。我们的架构如图2所示主要包含以下几个可定制的子模块2.1 特征提取器从设计参数到抽象特征输入是器件的设计参数化表示。对于参数化几何如光子晶体板的孔半径、周期这可能是一个向量对于自由形式设计如超表面的像素化掩模这可能是一张二维图像。特征提取器的任务是将这些原始设计映射到一个抽象的、高维的特征向量 $\phi$ 上。这个模块通常是一个标准的神经网络对于参数化输入可以使用全连接网络FCN。对于图像式输入卷积神经网络CNN或DenseNet等架构是更自然的选择因为它们能有效捕捉空间结构特征。这个模块是模型的主要可学习参数承载者负责捕捉设计中的复杂非线性关系。我们实践中发现为这个主干网络添加一个线性跳跃连接ResNet风格能显著提升训练的稳定性和收敛速度。2.2 模式模型预测每一个谐振的“指纹”这是QNM-Net的核心创新点之一。我们为预期中的每一个谐振模式都设置了一个独立的“模式模型”。每个模式模型是一个小型神经网络例如单层或双层全连接网络它以共享的特征向量 $\phi$ 为输入输出对应第m个QNM的物理参数复本征频率 $\tilde{\omega}_m$通常拆分为实部频率 $\omega_m$和虚部衰减率 $\gamma_m$分别预测。对于无耗材料我们可以通过设置 $\gamma_{nr,m}0$ 来施加约束。端口耦合振幅 $d_m$一个复数向量长度等于端口数N。其实部和虚部需要分别预测。关键技巧1模式数量的处理我们事先并不需要精确知道器件有多少个谐振模式。我们可以设置一个足够多的模式模型例如20个。如果某个设计在关心的频段内只有少数几个模式网络会学会将多余模式的谐振频率 $\omega_m$ 推到采样频段之外或者赋予其极大的衰减率 $\gamma_m$使得它们对S(ω)的贡献可以忽略不计。这赋予了模型处理未知模式数量的灵活性。关键技巧2对称性作为强约束这是体现“物理信息”威力的地方。例如对于具有镜面对称性的器件其QNM的模式对称性是确定的偶模的 $d_m$ 在两个对称端口上同号如(1, 1)奇模则反号如(1, -1)。如果我们从物理上知道主导模式是偶模就可以直接将 $d_1$ 固定为(1, 1)而不是让网络去学习。这极大地减少了需要学习的参数并强制模型遵循物理规律。2.3 背景模型与延迟模型刻画非谐振行为背景模型 $C(\omega)$负责预测非谐振的背景散射。它接收特征向量 $\phi$ 和频率 $\omega$ 作为输入。对于许多器件背景是缓变的。我们可以用一个小型网络预测背景参数如公式(4)中的相位角 $\alpha$在几个关键频率点的值然后进行插值来获得整个频段的 $C(\omega)$。同样可以利用物理约束简化它。对于无耗、互易的二端口系统$C(\omega)$ 必须是一个酉矩阵这可以转化为对输出参数的特定构造方式。延迟模型 $\tau$预测每个端口的相位延迟 $\tau_n$。这通常是一个与频率无关的实数。对于对称结构可以约束对称端口的延迟相等。2.4 物理模型层可微分的QNM合成器这是唯一没有可训练参数的“灰色盒子”。它接收所有子模型预测出的物理参数${ \tilde{\omega}_m, d_m }$, $C(\omega)$, $\tau$。然后严格按照前面给出的QNM展开公式计算出最终的散射矩阵 $S(\omega)$。这一层的关键在于它必须是完全可微分的这样才能允许梯度从损失函数比较预测的 $S(\omega)$ 和真实的 $S(\omega)$反向传播到前面所有神经网络子模块的参数中。我们利用现代深度学习框架如PyTorch、JAX的自动微分功能可以轻松实现这一点。训练流程简述输入一批设计参数和对应的频率点。特征提取器生成特征 $\phi$。各个模式模型、背景模型、延迟模型并行工作从 $\phi$ 预测出所有物理参数。物理模型层根据这些参数合成预测的 $S_{pred}(\omega)$。计算预测值与全波仿真真实值 $S_{true}(\omega)$ 之间的损失如S参数均方误差。通过自动微分计算梯度并更新所有神经网络子模块的权重。整个流程是端到端训练的目标是最小化最终的散射矩阵误差但学习的过程被物理定律严格地引导着。3. 实战应用与性能验证从光子晶体到自由形式超表面理论再优美也需要实验的验证。我们在两类典型且难度不同的电磁器件上测试了QNM-Net框架。3.1 案例一光子晶体平板——强物理约束下的高效学习系统描述我们研究一个在无损介质板上周期性排列空气孔的光子晶体平板。设计由五个参数控制板厚和四个环的孔半径。系统具有四重旋转对称性因此散射矩阵简化为2x2两个正交偏振态在正入射下解耦。频谱在关注频段内主要呈现一个明显的法诺谐振峰叠加在一个平滑的背景上。QNM-Net定制化对称性注入由于结构具有镜面对称性我们固定主导模式的端口振幅为 $d_1 (1, 1)$偶模。无耗约束材料无损因此设置非辐射衰减 $\gamma_{nr,1} 0$。背景矩阵参数化利用互易、无耗、对称性将背景矩阵 $C$ 参数化为一个仅由单个频率依赖相位角 $\alpha(\omega)$ 决定的酉矩阵如公式(4)所示。用一个小网络预测 $\alpha$ 在两个边界频率的值中间线性插值。模式数量仅使用一个模式模型因为频谱显示单谐振主导。数据与训练我们生成了10,000个随机设计及其在200个频率点上的散射矩阵作为数据集。QNM-Net仅用了约34k个可训练参数。结果令人振奋预测精度如图3(b)所示训练后的QNM-Net能极其精确地复现仿真得到的透射和反射谱测试集上的平均S参数均方误差S-MSE低至 $10^{-4}$ 量级。数据效率如图3(d)所示这是最具冲击力的结果。QNM-Net仅需160个训练样本占数据集的2%就能使验证损失低于 $10^{-3}$。而达到相同性能传统的全连接神经网络需要10倍以上的数据和更多的参数。这直接证明了物理先验知识带来的数据效率提升是数量级的。物理可解释性验证我们抽取了100个新的设计对比了QNM-Net预测的复本征频率 $\tilde{\omega}_1$ 与商用全波本征模求解器直接计算的结果。如图4所示两者高度吻合决定系数R² 0.999。这意味着QNM-Net不仅学会了拟合频谱曲线更真正学会了我们希望的底层物理——谐振频率。这打开了“知识发现”的大门如果我们事先不知道模式的对称性通过观察网络预测的 $d_1$ 向量我们可能会发现它总是接近(1, 1)或(1, -1)从而推断出结构的对称性。逆向设计演示我们以QNM-Net预测的谐振频率 $\tilde{\omega}_1$ 作为设计目标从随机初始设计出发利用自动微分计算梯度通过Adam优化器反向优化几何参数。在几秒钟内就成功设计出了一系列具有线性递增谐振频率和衰减率的光子晶体结构且全波仿真验证了设计结果的高度准确性图5。这展示了将QNM-Net作为可微分物理模型用于高效、精准逆向设计的潜力。3.2 案例二自由形式全介质超表面——弱物理先验下的挑战系统描述第二个案例更具挑战性一个自由形式的全介质超表面单元其结构由一个100x100的二值化位图定义。该系统设计空间巨大频谱包含多个重叠的谐振存在衬底导致结构不对称存在非辐射损耗且数据集只提供了从上表面入射的散射矩阵部分信息。QNM-Net定制化特征提取器升级由于输入是图像我们采用了卷积DenseNet作为特征提取器这是我们之前工作中验证有效的架构。背景简化对于如此复杂的结构没有明显的平滑背景。我们采用QNM展开中的常见做法设 $C -I$负单位矩阵让背景贡献由多个宽谐振来等效描述。模式数量设置了20个模式模型以容纳可能的多谐振特性。约束放松除了能量守恒和因果律由QNM公式保证外我们没有施加其他强约束如对称性。端口延 $\tau$ 设为固定值由端口到超表面的平均距离估算。结果分析数据效率优势依然明显如图6(d)所示在面对这个复杂得多的问题时QNM-Net要达到与最佳传统参考模型一个大型DenseNet后接分支网络相同的预测精仍然只需要后者约三分之一的训练数据。这证明了该框架的通用性和鲁棒性。预测精度与局限如图6(b)所示QNM-Net能够捕捉并复现频谱中的主要谐振峰但对于一些非常微弱的谐振特征预测存在偏差。我们认为这是由于自由形式超表面的设计空间过于庞大即使数万个训练样本也不足以让网络完全精确地学习所有细微的频谱特征。然而对于许多逆向设计应用抓住主要谐振特征已经足够。模式识别的智能性我们进一步对三个不同设计的预测模式进行了本征模验证图7。一个有趣的现象出现了QNM-Net预测的、在采样频段内且衰减较小的模式几乎都能在全波本征模求解结果中找到对应。然而本征模求解器还给出了许多QNM-Net没有预测到的模式。仔细分析发现这些“遗漏”的模式大多是与入射光耦合极弱的模式例如某些局域在单元内部特定区域、难以被外部端口激发的模式。这说明QNM-Net展现了一种“智能筛选”能力它自动识别并只学习那些对散射响应有显著贡献的模式而忽略了那些“沉默”的模式。相比之下传统的本征模求解器需要手动筛选大量模式QNM-Net的这一特性在实际应用中非常有用。4. 实现细节、调参心得与避坑指南将QNM-Net从理论框架转化为可运行的代码需要注意大量工程细节。以下是我在复现和实验过程中的一些核心经验。4.1 数据准备与预处理频率归一化输入的频率值 $\omega$ 需要进行归一化例如缩放到[0, 1]或[-1, 1]区间。这对于网络的稳定训练至关重要特别是当频率范围跨度较大时如从太赫兹到近红外。S参数处理散射矩阵 $S(\omega)$ 是复数。在训练时我们将其拆分为实部和虚部作为两个独立的输出通道进行回归。也可以考虑使用幅度和相位但对于谐振峰附近相位快速变化的情况实部/虚部表示通常更稳定。数据集划分务必确保训练集、验证集和测试集在参数空间分布上是一致的。对于几何参数化设计可以使用拉丁超立方采样来确保均匀覆盖。对于图像类设计随机划分通常即可但需检查划分后各集合的统计特性如像素密度分布是否相似。4.2 网络架构与初始化特征提取器深度与宽度这需要权衡。对于简单系统如光子晶体板一个3-5层的全连接网络可能就足够了。对于复杂系统如超表面则需要更深的卷积骨干网络。我们的经验是在参数量相近的情况下一个稍深但宽度较窄的网络往往比一个浅而宽的网络泛化能力更好。模式模型的大小每个模式模型是一个小型网络。我们发现对于大多数问题一个单隐藏层宽度8-32的MLP足以很好地预测[$\omega_m$, $\gamma_m$, Re($d_m$), Im($d_m$)]。过度参数化模式模型可能导致过拟合特别是当训练数据有限时。参数初始化谐振频率 $\omega_m$初始值应设定在目标频段内。可以均匀随机初始化。衰减率 $\gamma_m$必须初始化为正值因为 $\gamma_m -\Gamma_m/2$$\Gamma_m 0$ 是线宽。我们使用softplus激活函数来保证其输出恒为正gamma softplus(raw_output)。耦合振幅 $d_m$可以初始化为小的随机数。背景参数如果背景是缓变的初始化其网络输出接近零使得初始 $C(\omega)$ 接近单位矩阵或负单位矩阵这是一个合理的起点。4.3 损失函数与训练技巧损失函数选择我们使用S参数的均方误差MSE作为损失函数。对于复数即 $L \frac{1}{N_s N_\omega} \sum |S_{pred} - S_{true}|^2$。在某些强调带内性能的应用中可以引入频率加权。优化器与学习率Adam优化器是可靠的选择。学习率建议使用余弦退火或带热重启的余弦退火CosineAnnealingWarmRestarts。初始学习率可以设在1e-3到1e-4之间。梯度问题物理模型层涉及矩阵求逆(iωI - iΩ̃)^{-1}。当预测的谐振频率 $\omega_m$ 非常接近某个采样频率 $\omega$ 时会导致梯度爆炸。实践中我们从未遇到严重问题但一个稳健的做法是在训练初期对梯度进行裁剪gradient clipping。“模式竞争”与正则化当使用多个模式模型时可能会出现“模式竞争”现象几个模式试图解释同一个谐振峰。这通常不会影响最终的S(ω)预测精度但会使学习到的模式参数难以解释。可以加入轻微的正则化如鼓励不同模式的谐振频率 $\omega_m$ 彼此远离但这需要谨慎调整强度。4.4 常见问题排查训练损失震荡或不下降检查物理约束首先确认所有物理约束如对称性、无耗条件是否正确编码。一个错误的约束会导致网络在错误的假设下学习无法收敛。检查初始化糟糕的初始化可能导致网络一开始就陷入病态区域。尝试不同的随机种子或使用前面提到的针对性初始化策略。降低学习率这是最常用的手段。简化问题先用一个极简的、已知有解的设计如单个已知谐振频率的简单结构测试整个流程确保代码和物理公式正确无误。验证损失远高于训练损失过拟合增加数据这是最根本的但成本高。增强正则化对网络权重施加L2正则化或使用Dropout但在特征提取器的较低层使用需谨慎可能破坏特征连续性。减少网络容量减少特征提取器或模式模型的宽度/深度。利用数据增强对于几何参数化设计可以对参数施加微小扰动生成新样本对于超表面图像可以使用旋转、翻转等对称性操作如果物理允许来扩充数据。预测的频谱平滑但缺失尖锐谐振峰模式数量不足增加模式模型的数量。模式模型表达能力不足稍微增加模式模型小网络的容量。背景模型过于强势检查背景模型是否过于复杂以至于“吸收”了本该由谐振模式解释的特征。可以尝试简化背景模型如固定 $C$ 为常数矩阵。逆向设计优化不收敛或结果不物理检查梯度利用自动微分工具检查从目标物理参数如 $\tilde{\omega}_m$反向传播到设计参数的梯度是否存在非零。梯度消失可能意味着网络某部分不可微或出现了数值问题。优化目标可行性确认你设定的目标谐振频率等参数在物理上是可能实现的。要求一个超低损耗的谐振模式在一个本征损耗很大的材料中实现显然是不现实的。引入多目标或约束在逆向设计损失函数中除了主要目标如特定谐振频率可以加入对几何参数的约束如最小特征尺寸、面积约束以防止出现无法制造的极端设计。5. 框架的扩展潜力与未来方向QNM-Net的成功验证了将严格物理模型与深度学习结合的巨大潜力。这个框架本身是开放的有很多值得探索的扩展方向面向实验数据的训练传统神经网络对仿真数据的噪声和实验测量误差非常敏感。而QNM-Net由于内置了物理结构可能对噪声更具鲁棒性。我们的初步实验在训练数据中添加模拟噪声表明QNM-Net的性能下降确实比标准NN要小。这为直接利用验测量数据训练可解释的代理模型打开了大门对于校准困难或仿真-实验失配严重的系统尤其有价值。制造容差与鲁棒性设计在实际制造中几何尺寸偏差、材料不均匀性等会导致器件性能偏离设计值。可以将QNM-Net与表征制造误差的统计模型结合在逆向设计阶段就优化器件的鲁棒性。例如在优化过程中不仅优化目标性能还同时优化性能对关键几何参数的灵敏度。跨结构、跨频段的迁移学习QNM-Net的模块化特性非常适合迁移学习。例如在一个简单几何如纳米棒上预训练的特征提取器其学到的底层特征如边缘、角落、周期性的表征可能对更复杂的自由形式结构也有用。我们可以冻结特征提取器的部分层只微调模式模型从而用极少的新数据适配新器件。从频谱到本征模的联合训练目前QNM-Net仅从散射谱数据中学习模式参数。一个更激进的思路是进行多任务学习同时用散射谱数据和全波本征模求解器得到的精确QNM频率、场分布作为监督信号。这能进一步强化模型对物理本质的学习尤其对于频谱中难以分辨的尖锐谐振或准连续态BIC模式可能更有帮助。探索最优归纳偏置在本工作中我们手动注入了对称性、无耗等强约束。一个有趣的研究方向是系统性地量化每一种物理先验对称性、互易性、能量守恒、因果律对最终模型数据效率和泛化能力的贡献各有多大。这能指导我们在面对新问题时如何最有效地利用已知物理知识。在我个人的实践体会中QNM-Net这类物理信息模型最大的魅力在于它架起了“数据驱动”与“第一性原理”之间的桥梁。它没有抛弃我们数百年积累的物理智慧而是让神经网络成为执行这些物理定律的、更高效、更灵活的计算引擎。它解决的不仅仅是一个计算加速问题更是一个可信AI和可解释AI的问题。当你的模型不仅能告诉你“输出是什么”还能告诉你“为什么是这样的输出”因为谐振在这里耦合是那样的你对整个设计和优化过程的掌控感和洞察力会得到质的提升。当然这套方法对研究者也提出了更高要求你需要对你所研究系统的物理有深刻理解才能正确地将这些知识编码到模型架构中。这或许正是AI for Science走向深水区的必经之路——不是用AI替代物理学家而是让AI成为物理学家手中更强大的显微镜和计算尺。
