来源论文: https://arxiv.org/abs/2605.18492v1 生成时间: May 24, 2026 06:54

0. 执行摘要

在现代微电子技术中，多层印刷电路板（PCB）的电磁建模是设计优化与信号完整性分析的核心。本文解析的科研工作提出了一种基于2.5D矩量法（Method of Moments, MoM）与散射矩阵（S-matrix）框架的通用仿真方案。该方案的核心贡献在于：

  1. 格林函数构建的创新：利用S矩阵形式化方法导出了层状屏蔽介质中的完整双向格林函数，有效解决了传统传输线法在处理纵向电流（过孔）时的复杂性。
  2. 数值稳定性控制：通过S矩阵的递推公式与指数因子处理，确保了即使在处理极薄电介质层或高频Evanescent模态时，计算依然保持高精度且不发散。
  3. 计算效率优化：结合了FFT加速算法与块Toeplitz/Hankel矩阵技术，将稠密矩阵的填充与求解复杂度显著降低。
  4. 混合基函数建模：利用面屋顶（Rooftop）函数、体积脉冲（Pulse）函数及线性函数，实现了对平面布线与竖直过孔的统一建模。

本研究不仅为PCB辅助设计提供了高效工具，其在格林函数解析构造上的思路对于量子化学中处理准二维层状体系的电子结构计算亦有跨学科的借鉴意义。

1. 核心科学问题，理论基础，技术难点与方法细节

1.1 核心科学问题：复杂层状屏蔽环境下的场耦合

在屏蔽式多层PCB中，电磁场被限制在一个金属腔体（波导）内。传统的3D全波仿真（如有限元法FEM）需要对整个空间进行体网格剖分，这在处理厚度极小的金属薄片与高电导率材料时极其低效。矩量法（MoM）通过只对导体表面电流建模，极大地减少了未知数数量，但其核心难点在于构建能够精确反映层状波导边界条件的电场积分方程（EFIE）内核——即格林函数（Green’s Function）。

1.2 理论基础：2.5D矩量法与EFIE

研究基于频率响应的麦克斯韦方程组，采用时间谐波因子 $exp(-i\omega t)$。对于任意位置 $r$ 的总切向电场，其满足积分方程：

$$E^{exc}(r) = -ikZ \int_V \mathbf{G}(r, r') \mathbf{J}(r') dV' + Z_{met} \mathbf{J}(r)$$

其中，$\mathbf{G}(r, r')$ 是双向格林函数。2.5D的含义在于：虽然几何结构是三维层状的，但电流分布被假定在特定的金属平面（横向）或贯穿层的过孔（纵向）中。

1.3 技术难点：层状格林函数的数值不稳定性

在多层介质中，格林函数通常表达为无穷双重级数。传统的传输线类比法（Transmission-line approach）在计算不同层间的场耦合时，往往涉及双曲函数或大指数项的相乘，这容易导致浮点溢出或精度损失（特别是处理衰减模态时）。此外，处理包含横向电流（$J_x, J_y$）和纵向电流（$J_z$）的完整张量格林函数需要复杂的代数推导。

1.4 方法细节：S矩阵形式化构建格林函数

本工作的精髓在于放弃了直接推导场分量，而是借鉴了量子散射理论中的**Redheffer星乘（Star Product）**概念：

  1. 模式分解：将波导内的电磁场分解为TE（$h$）和TM（$e$）模式。在均匀层内，场振幅随 $z$ 轴以 $exp(i\beta z)$ 变化。

  2. S矩阵递推：定义每一层介质的传输矩阵和每一个界面的反射/透射矩阵。通过星乘运算 $S_{total} = S_1 \star S_2$，构建出从源层到观测层的整体散射响应。

  3. 自洽振幅计算：利用S矩阵组件定义“向上”和“向下”的传播因子。对于源点 $z_s$，计算其激发的自洽模式振幅，这包括了由于多层界面反复反射形成的几何级数效应。最终，格林函数分量被表达为：

    $$Z_{mn}^{\alpha\alpha'} (z_v, z_s) = F(S_{above}, S_{below}, S_{pair}, e^{i\beta z})$$

    这种表达方式通过提取正指数因子，确保了数值计算的绝对稳定性。

  4. 基函数选择：

    • 平面金属：采用平面Rooftop基函数，能够自然地满足电流连续性方程。
    • 垂直过孔：采用体积Pulse基函数模拟恒定电流，或线性基函数模拟沿过孔长度变化的电流。

2. 关键 Benchmark 体系、计算数据与性能数据

2.1 Benchmark 体系 1：三层六带状线滤波器

该模型参照了经典文献中的屏蔽滤波器结构。其物理参数如下：

  • 层结构：总共三层介质。层1（$h=7.0$ mm, $\epsilon_r=1.0$），层2（$h=0.51$ mm, $\epsilon_r=2.33$），层3（$h=7.0$ mm, $\epsilon_r=1.0$）。
  • 金属层：C1层包含4个金属条，C2层包含2个金属条。所有金属均为理想导体（PEC）且厚度忽略不计。
  • 仿真设置：频率范围 1-4 GHz，网格步长 0.1 mm，对应 $300 \times 250$ 的宏观网格。

计算结果：

  • S参数精度：S11（反射系数）与S21（传输系数）的曲线与已知文献（Reference [82]）高度吻合。在2-3 GHz的通带内，损耗极低，准确捕捉到了谐振峰位置。
  • 收敛性：由于采用了S矩阵法，格林函数级数的项数即使增加到数千项，也不会出现数值振荡。

2.2 Benchmark 体系 2：带过孔的复杂互连结构

在体系1的基础上，通过在C1和C2层之间增加两个垂直过孔（Vias）来测试纵向电流处理能力。过孔被建模为跨越介质层的体积电流元素。

关键数据对比：

  • 引入过孔的影响：S21曲线在 2.5 GHz 附近出现了明显的频移和带宽变化，这反映了过孔引入的寄生电感效应。模型能够精确计算 $Z^{xz}, Z^{zx}, Z^{zz}$ 等交叉分量，这些分量在纯2D仿真中是被忽略的。
  • 性能数据：使用FFT加速后，矩阵填充时间从 $O(N^2)$ 降低到近乎 $O(N \log N)$。对于上述网格规模，单频点求解时间保持在秒级。

2.3 性能优势总结

  • 内存占用：由于采用规则网格，MoM矩阵呈现块Toeplitz-Hankel结构。只需存储矩阵的第一行和第一列，内存开销降低了约 2 个数量级。
  • 数值一致性：通过公式 (78) 引入的重缩放（Rescaling）技术，确保了横向和纵向电流在不同基函数量纲下的数值稳定性，避免了条件数病态问题。

3. 代码实现细节与复现指南

3.1 软件包与开发环境

  • 核心语言：C++ (C++17标准)。
  • 数学库：
    • FFTW3：用于执行高效的二维快速傅里叶变换。
    • Eigen/OpenBLAS：处理稠密矩阵运算与线性方程组求解（如BiCGSTAB或GMRES迭代法）。
  • 开源参考：虽然本文代码为专有实现，但其算法逻辑高度遵循 Sonnet Software 的理论架构。开发者可参考开源矩量法项目如 OpenEMS 的底层逻辑，结合本文提供的格林函数公式进行修改。

3.2 关键算法实现步骤

  1. S矩阵预计算模块：

    • 输入：各层厚度 $h_l$、介电常数 $\epsilon_l$、频率 $\omega$、网格尺寸 $a_x, a_y$。
    • 计算：利用公式 (13)-(17) 递归生成每一层的 $S^{(L)}$ 和界面 $S^{(I)}$。
    • 输出：存储三类聚合S矩阵（$S^A, S^B, S^P$）。

  2. 格林函数加速填充模块：

    • 实现公式 (85)-(89) 中的核心加速逻辑。
    • 注意：在执行二维FFT前，需要对矩阵进行零填充（Padding），以模拟非循环卷积效果。
    • 难点提示：必须根据公式 (80)-(83) 正确映射波导模态索引 $(m, n)$ 到FFT频率索引。

  3. 重缩放处理：

    • 在组装全局 $Z$ 矩阵前，应用公式 (78) 的缩放因子，将电压/电流单位统一为欧姆/伏特/安培。这是复现成功、避免求解器不收敛的关键。

3.3 复现检查点

  • 单层波导校验：先实现自由空间波导格林函数，验证其收敛性。
  • TE/TM 匹配：检查 $Z^{xy}$ 和 $Z^{yx}$ 的对称性。
  • 奇异性处理：在计算重叠积分时，对于 $r=r'$ 的自相互作用项，应采用解析积分（公式 59-61）以消除奇异性。

4. 关键引用文献与局限性评论

4.1 关键引用文献

  • [8] W. C. Gibson, “The Method of Moments in Electromagnetics”：矩量法的基础教材，提供了EFIE的基础框架。
  • [71] Sonnet Software：工业界屏蔽层状介质仿真的标杆，本工作在思路上是对其商业算法的一种学术透明化实现。
  • [80] R. Redheffer, “On the relation of transmission-line theory to scattering and transfer”：散射矩阵星乘运算的理论源头，是本文数值稳定性的保证。
  • [31] B. J. Rautio et al.：统一FFT加速算法（Unified-FFT）的提出者，本文的加速方案深受其启发。

4.2 工作局限性评论

  1. 几何约束：该方法高度依赖于波导壁的侧向屏蔽条件（PEC壁）。虽然可以通过匹配阻抗边界模拟开放空间，但在处理非规则腔体时灵活性较差。
  2. 过孔电流假设：文章假设过孔长度远小于波长，从而将电流视为恒定。对于极高频（太赫兹级别）或超厚PCB，需要引入更高级的多模态纵向基函数。
  3. 金属厚度：目前的2.5D模型忽略了平面金属的物理厚度。在趋肤深度与金属厚度可比拟时，需要引入表面阻抗修正或厚金属模型。
  4. 计算资源：尽管有FFT加速，但当层数超过50层且布线极其复杂时，S矩阵的递归次数和格林函数级数的项数仍会导致显著的计算开销。

5. 补充论述：格林函数在计算物理中的普适性

虽然本文讨论的是PCB仿真，但其处理层状体系的思想在量子化学和凝聚态物理中具有深远的普适性。

5.1 从电磁波到电子波函数

在处理二维范德华异质结（如石墨烯/氮化硼叠层）的电子输运时，其薛定谔方程或狄拉克方程的解同样可以表达为层状格林函数的形式。本文中使用的S矩阵递推法，在物理上等价于电子结构计算中的递归格林函数法（RGF）。两者的核心目标一致：通过消去内部层自由度，仅保留界面或源位置的响应。

5.2 离散化与基函数的跨学科映射

  • MoM vs LCAO：本文的矩量法将连续电流投影到基函数上，这与量子化学中将波函数投影到原子轨道线性组合（LCAO）具有数学上的对等性。本文提供的解析重叠积分计算（公式 59-69）对于提高量子化学计算中库仑积分的精度具有参考意义。
  • FFT加速：本文对块Toeplitz矩阵的FFT处理，可以平移到周期性体系的哈密顿矩阵快速求逆中。

5.3 数值稳定性的启示

本研究特别强调了“正长度指数因子”对稳定性的重要性。在计算重费米子体系或具有复杂势垒的隧穿效应时，传统的格林函数传播子常因指数衰减项过小而导致数值下溢。本文的S矩阵框架通过“双向分解+星乘合并”，优雅地规避了这一问题，为开发高稳定性的量子输运模拟器提供了清晰的路径。

5.4 未来方向：非局域响应与材料色散

未来的研究可以进一步引入金属的非局域电磁响应模型，或者结合量子流体力学模型，将此PCB仿真框架扩展到等离激元光子学芯片的模拟中。这种跨越宏观电路与微观量子效应的建模方式，将是未来算力芯片设计的关键技术。