基于多线直通反射线校准的325 GHz互连表征技术详解

1. 项目概述：为什么325 GHz互连表征是个“硬骨头”？

在毫米波乃至太赫兹频段搞电路设计，工程师们最头疼的问题之一，可能就是“我的信号到底在互连上损耗了多少？相位变了多少？” 这听起来像是个基础问题，但在频率飙升至325 GHz甚至更高时，它却成了一个极具挑战性的测量难题。想象一下，你设计了一个精妙的异构集成芯片，不同功能的电路层通过垂直互连（比如硅通孔TSV、倒装焊球或者文中提到的宽带侧耦合共面波导）连接。你想精确知道这个互连本身的性能，而不是连带着它前后一大段传输线的混合结果。传统方法要求你在互连的两侧都能直接放置探针进行测量，这在高密度、多层堆叠的先进封装中，往往是个“不可能完成的任务”——下层电路可能被完全覆盖，根本没有给你下针的空间。

这就引出了我们这次要深入探讨的核心：基于多线直通反射线校准的互连表征技术。这项技术的目标非常明确——在不直接接触互连端口的情况下，依然能把它“揪”出来，单独分析它的散射参数。其核心思想颇具巧思：与其费力去直接测量一个“孤立”的互连，不如主动把它“藏”进我们整个测量系统的误差模型里，然后再通过精密的数学手段把它“剥离”出来。具体来说，就是设计一套特殊的校准结构，让待测互连成为校准误差盒的一部分，然后通过多级校准流程，最终反演出互连本身的特性。这项工作的价值，对于从事高频、高速电路（如5G/6G前端、雷达芯片、太赫兹成像系统）研发和封装的工程师而言，是实实在在的。它提供了一套可复现、高精度的“标尺”，让你能在设计阶段就准确预测互连性能，或在生产测试中有效进行质量管控，避免因互连性能不达标而导致整个系统失效。

2. 核心原理拆解：从“黑箱”到“透明盒”的数学与物理

要理解这套方法，我们得先抛开复杂的公式，从两个基本概念说起：误差盒和T矩阵（传输矩阵）。

2.1 误差盒与T矩阵：系统误差的“集装箱”

任何矢量网络分析仪的测量，都无法避免地包含了测试系统本身（如电缆、探针、适配器）引入的误差。我们可以把这些误差效应抽象成两个“黑箱”，分别放在被测器件（DUT）的左边和右边，这就是误差盒（Error Box）。VNA测到的数据，实际上是“左误差盒 + DUT + 右误差盒”这个整体级联的结果。校准的目的，就是通过测量一系列已知特性的标准件（如直通、短路、负载、延迟线），来解算出这两个误差盒的特性，从而在后续测量中将其“扣除”，得到DUT真实的S参数。

而T矩阵是描述双端口网络的一种非常有力的数学工具。它的最大优点在于，对于级联的网络，总的T矩阵等于各个网络T矩阵的乘积。这完美契合了“误差盒-DUT-误差盒”的级联模型。如果我们将左误差盒记为X，右误差盒记为Ȳ，DUT的T矩阵记为T_DUT，那么VNA测量到的整体T矩阵M就可以表示为：M = X * T_DUT * Ȳ。校准，就是通过测量已知T_DUT的标准件，来求解X和Ȳ。

2.2 传统直接法 vs. 本文嵌入法

传统的互连表征思路是直接法，如图1(a)所示。它要求互连两侧的传输线层都能被直接探测。这样，我们可以先在左侧和右侧分别进行独立的校准（例如SOLT或LRRM），得到针对各自端口的误差盒X和Ȳ。然后，将互连作为DUT进行测量，得到M = X * T_IC * Ȳ。由于X和Ȳ已知，通过矩阵运算就能直接解出互连的T矩阵T_IC。这个方法直观、准确，但前提苛刻：你必须能接触到互连的每一侧。

当无法直接探测互连的某一侧时（这在多层异构集成中非常常见），直接法就失效了。此时，报告的S参数实际上混杂了互连和一侧甚至两侧传输线的效应，导致插入损耗被高估，相位信息失真，无法用于精准建模。

本文提出的嵌入法则另辟蹊径，其核心思想如图1(b)所示。既然无法在互连两侧单独校准，那就把互连本身设计进校准结构里。我们制作一对相同的互连（T_IC1和T_IC2），让它们分别成为左、右误差盒的一部分。也就是说，新的左误差盒变成了X = X‘ * T_IC1，新的右误差盒变成了Ȳ = T_IC2 * Ȳ‘。这里的X‘和Ȳ‘是排除了互连之后、我们可以通过其他方式校准的“剩余”误差盒。

此时，我们测量一个校准标准件C_i（比如一段直通传输线），得到的测量模型是：M_i = (X‘ * T_IC1) * C_i * (T_IC2 * Ȳ‘)。如果我们能通过另一套独立的校准手段确定X‘和Ȳ‘，那么从M_i中就能同时解出T_IC1和T_IC2。这就像把互连“溶解”在了校准液中，再通过结晶把它单独提取出来。

2.3 为什么是mTRL？超越“直通唯一”法的局限

最简单的嵌入法可能是“直通唯一”法，它只用一个直通标准件。但这存在根本性缺陷：一个测量方程无法唯一确定两个未知矩阵（T_IC1和T_IC2）。它通常假设两个互连完全对称（T_IC1 = T_IC2），且要求直通标准件的长度非常短以减少串扰，但这又可能引入新的不确定性。

多线直通反射线校准的强大之处在于，它引入了多个不同长度的传输线作为标准件。这提供了丰富的相位信息，使得方程组足够多，从而能够唯一地求解出误差盒，进而也就能唯一地二分出两个互连的特性，即使它们不完全对称。mTRL校准本身是业界公认的高精度方法，尤其适用于非50欧姆环境和非同轴介质，它通过多个线标准有效剥离了系统误差，并能够提取传输线的传播常数和特性阻抗。本文将mTRL与嵌入思想结合，相当于用最精密的“尺子”去丈量被隐藏的“目标”，从而实现了在325 GHz高频下的精确表征。

实操心得：理解“参考面平移”这是整个技术的灵魂。校准的本质就是确定一个“参考面”，在此参考面上，VNA的测量值等于DUT的真实值。本文方法的高明之处在于，它通过嵌入和去嵌入，成功地将测量参考面从探针尖端，“平移”到了互连的边界，甚至更远。这意味着，即使互连深埋在封装内部，我们依然能在“纸面上”将虚拟的探针伸到其端口进行“测量”，这对于复杂系统的仿真与测试对标至关重要。

3. 实验实现：从理论到晶圆的每一步

纸上谈兵终觉浅，我们来看NIST团队是如何将这一套精妙的数学物理模型，落地为可重复、可测量的硅上现实的。

3.1 校准结构设计与分区策略

为了实现多级去嵌入，他们巧妙地设计了三个校准区域，这是一种“层层剥笋”的思路：

• 区域1：最基础的校准。仅包含从探针尖端到共面波导传输线（CPW）的过渡（T_P1, T_P2）以及第一段缓冲传输线（T_1）。这个区域的校准结果X1, Ȳ1主要用于验证从仿真中提取的传输线参数（如特性阻抗Z1）的准确性，为后续步骤奠定可信基础。
• 区域2：在区域1的基础上，增加了第二段传输线T_2。此时误差盒变为X2 = X1 * T_1 * T_2和Ȳ2 = T_2 * T_1 * Ȳ1。这个校准的目的是为了获取T_2段的传播常数γ2和特性阻抗Z2，并得到排除了互连的“干净”误差盒X‘和Ȳ‘（在此模型中，X‘ = X1 * T_1,Ȳ‘ = T_1 * Ȳ1）。
• 区域3：核心区域。在区域2的结构中，插入了我们待测的互连T_IC1和T_IC2，以及一段用于模式净化的缓冲线T_3。此时误差盒包含了互连：X3 = X‘ * T_2 * T_IC1,Ȳ3 = T_IC2 * T_2 * Ȳ‘。

设计的关键细节：

1. 缓冲线：每个区域之间、以及互连与标准件之间，都设计了0.21 mm的“半直通”缓冲线（如T_3）。它的作用至关重要——让可能被激发的非传播模式（凋落模）充分衰减，确保到达校准标准件（如直通、延迟线）端口的只有纯净的CPW主模。这是保证校准精度的物理基础。
2. 阻抗匹配：尽管三个区域的CPW几何尺寸因层叠结构不同而有所差异（参见原文表III），但设计时都努力将其特性阻抗控制在50欧姆附近。这最小化了因阻抗失配引起的反射，简化了分析模型。
3. 传输线长度优化：mTRL标准件中多个传输线的长度不是随意选的。他们通过优化算法，在10 GHz到1 THz的频带内最小化校准的归一化标准偏差，从而确定了各区域最优的线长组合（原文表II）。这确保了在整个宽频带内校准的稳定性。

3.2 制造工艺：实现可重复性的基石

所有校准结构都在同一片3英寸熔融石英晶圆上共面制造。共面制造是验证该方法有效性的关键前提，因为它确保了出现在不同区域校准结构中的相同模块（如T_P1, T_1, T_2）具有完全一致的物理和电学特性。如果这些模块是分开制造的，工艺波动会引入无法区分的差异，从而破坏去嵌入公式成立的基础假设。

工艺流主要包括五个层次（原文图5）：

1. PdAu电阻层：用于制作校准中的集总电阻标准件和晶圆对准标记。
2. SU-8基底层：光刻形成665 nm厚的SU-8介质层，定义了CPW传输线的“沟槽”区域。
3. CPW 1/2层：沉积665 nm金，形成底部两层CPW信号线，它们被嵌入SU-8沟槽中，实现了与介质的共面性。
4. SU-8间隔层：旋涂并图形化一层2.5 µm厚的SU-8，这层介质决定了顶部CPW与底部CPW之间的垂直间距，是形成宽带侧耦合的关键。
5. CPW 3层：最后沉积500 nm金，形成顶部的CPW信号线。

这种单片集成工艺保证了极高的结构一致性和对准精度，为高频率测量提供了可靠的硬件基础。

3.3 测量与仿真：双管齐下的验证

测量方面：工作覆盖了DC-110 GHz（同轴）、140-220 GHz（WR-5）和220-325 GHz（WR-3）三个频段，使用50 µm间距的接地-信号-接地波导探针。每个校准套件测量10个结构：1个直通、8个不同长度的线、1个短路器。对于区域1，额外测量了串联电阻和电容标准件，用于提取单位长度电容和设置参考阻抗。

仿真方面：采用了2D和3D仿真结合的策略。

• 2D仿真：使用ANSYS 2D Extractor分析各区域CPW的横截面，提取单位长度的分布电阻、电感、电容和电导参数。这些参数与从测量中得到的传播常数γ结合，通过公式Z = (R + jωL)/γ计算各区域的复特性阻抗，这是去嵌入公式（2）和（3）所必需的。
• 3D全波仿真：建立了包含所有校准结构的完整3D模型，模拟了整个测量过程。这相当于在计算机里完美复现了一次实验，用于验证实验设计和数据处理流程的正确性。仿真中采用了分频段自适应网格剖分等技术来保证高达325 GHz的计算精度。

注意事项：材料参数的重要性高频仿真，尤其是涉及SU-8这类高分子介质的仿真，材料模型的准确性至关重要。文中提到，在220 GHz以上，测量与仿真的衰减常数出现偏差，推测原因可能是SU-8的Cole-Cole介电模型在110 GHz以上数据缺失。这提醒我们，在进行极高频率设计时，必须尽可能获取或测量所用材料在目标频段的准确电磁参数，否则仿真结果会失去指导意义。

4. 多级去嵌入操作流程详解

有了测量数据，接下来就是施展“数学魔法”，一步步将互连的S参数剥离出来。这个过程清晰体现了“分层剥离”的思想。

4.1 第一步：执行各级mTRL校准

使用测量得到的原始S参数数据，分别对区域1、区域2、区域3执行标准的mTRL校准算法。这一步可以通过开源工具如Python的scikit-rf库方便地完成。每个校准会输出两样关键结果：

1. 误差盒矩阵：对于区域2，得到X2和Ȳ2；对于区域3，得到X3和Ȳ3。
2. 传输线参数：各区域的传播常数γ_k和特性阻抗Z_k。

4.2 第二步：计算中间传输线矩阵

利用第一步得到的γ_k和Z_k，我们可以根据传输线理论，计算出特定长度传输线段（如T_1, T_2）的T矩阵。计算公式的核心是考虑阻抗变换的传输矩阵：T_line = Q(50->Z) * diag(exp(-γ*l), exp(γ*l)) * Q(Z->50)其中，Q是阻抗变换矩阵，l是传输线长度，diag是对角矩阵。这个计算将传输线矩阵统一转换到50欧姆参考阻抗下，便于后续运算。

4.3 第三步：核心去嵌入计算

这是提取互连参数的关键步骤。根据模型，我们有：

• X3 = X‘ * T_2 * T_IC1
• Ȳ3 = T_IC2 * T_2 * Ȳ‘
• X2 = X‘ * T_2
• Ȳ2 = T_2 * Ȳ‘

通过矩阵运算，可以解出：

• T_IC1 = inv( Q(50->Z2) * inv(X2) * X3 * Q(Z3->50) )
• T_IC2 = inv( Q(50->Z3) * Ȳ3 * inv(Ȳ2) * Q(Z2->50) )

这里的inv表示矩阵求逆。公式（2）和（3）原文给出了考虑阻抗变换的完整形式。实际操作中，利用scikit-rf等工具的矩阵运算功能，可以相对简洁地实现这一过程。

4.4 第四步：结果验证与对比

提取出T_IC1和T_IC2后，将其转换为更常用的S参数形式S_IC，即可得到互连本身的散射参数。为了验证结果的正确性，本文采用了三重验证：

1. 实验去嵌入结果：即上述流程从实际测量数据中得到的结果。
2. 全波仿真去嵌入结果：对3D仿真模型的数据，采用完全相同的流程进行处理。
3. 解析模型结果：基于分布式传输线理论建立的三模解析模型（源自团队先前工作），该模型无任何拟合参数，仅基于几何尺寸和材料属性。

5. 结果分析与工程启示

将三种方法得到的互连S参数进行对比（对应原文图7），是验证该方法有效性的最终环节。

5.1 一致性分析

对于500 µm和250 µm两种耦合长度的宽带侧耦合CPW互连，在DC-325 GHz的整个频段内，实验测量、全波仿真和解析模型三者展现出了极好的一致性。无论是传输参数还是反射参数，三条曲线都基本重合。这强有力地证明了：

1. 嵌入与去嵌入方法的正确性：数学模型和实验流程是可靠的。
2. 制造工艺的可重复性：共面制造保证了结构的一致性。
3. 仿真模型的准确性：2D/3D仿真能有效预测高频行为。

5.2 差异与误差来源讨论

尽管一致性很好，细微的差异仍然存在，分析这些差异来源对工程实践更有价值：

• 测量中的纹波：实验曲线在部分频点存在微小纹波，而仿真和模型曲线则较平滑。这很可能源于辐射损耗或探针激励起的杂散模式。在325 GHz这样的高频下，任何微小的不连续性都可能激励起高次模或导致能量辐射，这些效应在理想的解析模型和边界条件设置完善的仿真中可能被抑制，但在实际测量中无法完全避免。
• 仿真中的纹波：全波仿真结果在反射系数较低（< -10 dB）时也存在细微纹波。这通常与网格剖分精度、仿真边界条件设置（如是否完美匹配层）以及频带内是否使用单一网格有关。自适应网格在不同频点可能略有差异，导致结果出现微小波动。
• 高频损耗偏差：在220 GHz以上，区域2和区域3的测量衰减常数高于仿真值。如之前所述，这指向了SU-8介质损耗模型在高频下的不确定性。这是当前太赫兹器件设计中的一个普遍挑战，即缺乏宽频带、高精度的介质材料数据库。

5.3 方法优势与适用范围总结

基于mTRL校准的互连表征技术，其核心优势在于：

• 解决不可达性问题：为无法直接探测端口的多层异构集成互连提供了唯一的精确表征途径。
• 高精度：得益于mTRL校准自身的高精度和多重验证，结果可信度高。
• 通用性强：方法原理适用于各种类型的互连（接触式、非接触式）、传输线（微带线、带状线、CPW）和频率范围。
• 提供完整模型：不仅能得到S参数，还能通过过程获取传输线的传播常数和特性阻抗，有利于建立更精准的分布参数模型。

它的“代价”是需要额外的芯片面积来制作包含互连的复杂校准结构，并且数据处理流程相对复杂。因此，它更适用于技术开发、模型提取和标准建立阶段。一旦通过该方法建立了某种互连的精准紧凑模型，该模型就可以用于后续大规模电路设计，而无需对每个互连都进行如此复杂的测量。

6. 常见问题与实操陷阱规避

在实际尝试复现或应用该方法时，可能会遇到一些典型问题，以下是一些排查思路和经验之谈。

6.1 校准失败或结果不收敛

• 问题现象：mTRL校准算法报错，或提取的误差盒、传播常数出现非物理值（如衰减为负）。
• 排查步骤：
  1. 检查标准件模型：确保在算法中为每个线标准件输入的物理长度和电气模型（是否为无损传输线）绝对准确。一个单位的输入错误会导致全盘皆输。
  2. 检查测量数据质量：观察原始S参数。直通件的插入损耗是否在合理范围？反射是否足够小？延迟线的相位变化是否平滑？在325 GHz，探针接触的重复性、电缆的轻微移动都可能引入误差。确保测量环境稳定，并多次测量取平均。
  3. 验证阻抗定义：mTRL校准可以提取出传输线的特性阻抗。将其与2D仿真或理论计算值对比。如果差异巨大，可能是标准件设计或测量有问题。
  4. 简化问题：先尝试对一个简单的、可直接探测的CPW传输线进行标准mTRL校准，验证整个测量和数据处理链路是否正常。

6.2 去嵌入后互连S参数不被动或不互易

• 问题现象：提取出的互连S参数不满足无源网络的条件（对于无源互连，S矩阵应满足无源性：I - S^H * S为正定矩阵；通常也满足互易性：S21 = S12）。
• 可能原因及解决：
  1. 误差累积：多级校准和矩阵运算会放大测量噪声。确保每一级校准（区域1，2，3）都尽可能精确。增加每个校准套件中线标准的数量（多于8个）可能有助于提高统计稳定性。
  2. 模型不匹配：在去嵌入公式中，我们假设不同区域中相同的传输线段（如T_2）的T矩阵完全相等。如果由于制造偏差导致它们实际上有细微差别，就会引入误差。共面制造是减少此误差的关键。在无法共面制造时，需通过严格的工艺控制或额外的表征来评估此偏差。
  3. 缓冲线长度不足：如果缓冲线（如T_3）太短，非主模能量未能充分衰减，会导致校准标准件端口处的模式不纯，从而使校准基准失效。确保缓冲线长度足以让高次模衰减到可忽略水平（通常需要几个趋肤深度的长度）。

6.3 高频段（>200 GHz）结果抖动大

• 问题现象：在频率高端，提取的S参数曲线出现明显抖动或跳变。
• 应对策略：
  1. 探针与校准：确保使用高质量、经过厂商校准的毫米波/太赫兹探针。在如此高频段，探针本身的性能及其校准状态（使用阻抗标准基板ISS）至关重要。考虑使用更精密的探针台和更稳定的连接。
  2. 仿真辅助诊断：进行全波仿真，并在仿真中引入微小的几何变化（如边缘粗糙度）或材料参数变化，观察是否会出现类似抖动。这有助于判断是测量系统问题还是结构本身的高频谐振效应。
  3. 数据平滑与滤波：在最终分析前，可对原始测量数据施加轻微的频率域平滑（如Savitzky-Golay滤波器），但需谨慎，避免掩盖真实的物理现象。更好的做法是增加测量点数，提高数据密度。

6.4 与紧凑模型对接困难

• 问题现象：提取的宽频带S参数数据量庞大，难以直接用于电路级仿真或转化为SPICE兼容的紧凑模型。
• 解决方案：
  1. 矢量拟合：使用矢量拟合技术，将宽频S参数用一组有理函数（极点、留数）来逼近。这能生成一个频域连续、因果且被动的模型，可以轻松转换为时域卷积模型或等效电路。
  2. 建立参数化模型：如果测试了不同尺寸（如不同耦合长度、宽度）的互连，可以利用提取的S参数库，通过机器学习或插值方法，建立一个几何参数到S参数的映射模型。这样，设计师只需输入尺寸，即可快速获得模型。
  3. 与EDA工具集成：将矢量拟合后的模型或S参数Touchstone文件，直接导入ADS、HFSS Circuit等EDA工具，用于系统级链路仿真。

在我个人的工程实践中，这套方法最深刻的体会是“规划优于补救”。成功的关键在于最初的实验设计：校准结构的布局、缓冲线的长度、阻抗的粗略匹配、以及最重要的——所有相关结构的共面制造。数据处理流程虽然复杂，但一旦脚本化后就变得可重复。真正的挑战往往在于前期的工艺实现和测量设置，任何一个环节的疏忽都可能导致后续精巧的数学处理失去意义。因此，对于希望应用此方法的团队，我建议投入足够资源在工艺流片设计和测量平台搭建上，并与仿真紧密协同，用仿真来指导设计和预判问题，这样才能在高达325 GHz的征途上，真正驾驭互连的“微观世界”。