4. 共面波导传输线特征阻抗与有效介电常数提取

为什么共面波导参数必须“先提取、再设计”

共面波导（CPW）由中心导带、两侧接地平面和它们之间的两条缝隙构成，所有导体位于衬底同一表面。对超导量子芯片而言，CPW 既是传输线，也是谐振腔、耦合器和布线的基础结构。它的两个核心电学参数是特征阻抗 Z₀ 与有效介电常数 ε_eff：前者决定与 50 Ω 测量链路的匹配，后者决定波在芯片上的相速度与谐振频率。二者都由横截面几何（中心导带宽度 W、缝隙宽度 S、衬底厚度 h、金属厚度 t）和材料（衬底相对介电常数 ε_r）共同决定。

本案例的任务是：给定一组 CPW 截面参数，用保角映射解析公式提取 Z₀ 与 ε_eff，并定量分析衬底厚度 h 与金属间隙 S 对微波传输特性的影响。解析公式适用于准静态近似（横截面尺寸远小于波长），在微波频段对常规 CPW 几何具有足够精度，可作为设计初值和趋势分析的快速工具；最终定版仍需全波电磁仿真校核。

提取模型与解析公式

采用无金属厚度修正的保角映射结果。定义两个几何比值：

k  = W / (W + 2S)
k' = sqrt(1 - k²)

K(k)  : 第一类完全椭圆积分
K(k') : 互补模的完全椭圆积分

当衬底厚度有限时，用 h 构造第二个模，用于描述电场穿透衬底的程度：

k1 = sinh(πW / (4h)) / sinh(π(W + 2S) / (4h))
k1' = sqrt(1 - k1²)

于是有效介电常数与特征阻抗为：

ε_eff = 1 + (ε_r - 1)/2 · [K(k')/K(k)] · [K(k1)/K(k1')]

Z₀ = (30π / √ε_eff) · K(k') / K(k)      [Ω]

注意两个极限行为：当 h → ∞ 时 k1 → k，方括号项趋于 1，退化为半无限厚衬底公式；当 h 很小时电场主要落在衬底内，ε_eff 趋近于 (1+ε_r)/2 附近。椭圆积分的比值用算术几何平均（AGM）迭代计算，收敛极快，适合写成短函数。

用 AGM 实现椭圆积分比

完全椭圆积分的比值 K(k')/K(k) 不需要单独求积分，用 AGM 迭代即可：

import math

def agm(a, b, tol=1e-15):
    while abs(a - b) > tol:
        a, b = 0.5 * (a + b), math.sqrt(a * b)
    return a

def K_over_Kp(k):
    """返回 K(k)/K(k')，k 为模（0<k<1）"""
    kp = math.sqrt(1.0 - k * k)
    return agm(1.0, kp) / agm(1.0, k)

def cpw_params(W, S, h, eps_r):
    """W, S, h 单位 μm；返回 (Z0 Ω, eps_eff)"""
    k  = W / (W + 2.0 * S)
    k1 = math.sinh(math.pi * W / (4.0 * h)) / math.sinh(math.pi * (W + 2.0 * S) / (4.0 * h))
    ratio_k  = K_over_Kp(k)     # K(k)/K(k')
    ratio_k1 = K_over_Kp(k1)    # K(k1)/K(k1')
    eps_eff = 1.0 + (eps_r - 1.0) / 2.0 * (1.0 / ratio_k) * ratio_k1
    Z0 = (30.0 * math.pi / math.sqrt(eps_eff)) * (1.0 / ratio_k)
    return Z0, eps_eff

注意公式中 Z₀ 用的是 K(k')/K(k)，即 1/ratio_k；ε_eff 中的 [K(k')/K(k)] 同样取倒数，而 [K(k1)/K(k1')] 直接取 ratio_k1。这两处方向容易写反，是提取结果偏差的常见来源。

基准算例：高阻硅上的 50 Ω 设计

取高阻硅衬底 ε_r = 11.9，衬底厚度 h = 500 μm，目标中心导带 W = 10 μm，缝隙 S = 6 μm。代入函数：

Z0, ee = cpw_params(W=10.0, S=6.0, h=500.0, eps_r=11.9)
print(f"Z0 = {Z0:.2f} Ω, eps_eff = {ee:.3f}")

按上述公式计算，结果约为 Z₀ ≈ 50.6 Ω、ε_eff ≈ 6.44（示例计算值，由给定输入直接推出）。ε_eff 落在空气（1）与硅（11.9）之间且明显偏向衬底一侧，说明大部分电场能量集中在衬底中；相速度 v_p = c/√ε_eff ≈ 1.18×10⁸ m/s，约为真空光速的 39%。若该 CPW 用作 λ/4 谐振腔，基频 f = v_p/(4L)，长度 L = 5 mm 时对应约 5.9 GHz，这一量级与超导量子芯片常用的读取/比特频率区间吻合。

衬底厚度 h 的敏感性

固定 W = 10 μm、S = 6 μm、ε_r = 11.9，只改变 h：

（表中为按公式计算的示例值。）趋势很清楚：h 从 50 μm 增到 500 μm，Z₀ 上升约 5.7 Ω、ε_eff 下降约 1.1；而 h 从 500 μm 再增到 1000 μm，Z₀ 只再升约 0.4 Ω。也就是说，当 h 大于约 3~4 倍 (W+2S) 后，电场已基本“看不到”衬底背面，参数趋于半无限厚衬底的渐近值。工程含义是：减薄衬底（如为倒装焊或通孔工艺减薄到 100 μm 以下）会显著降低 Z₀ 并抬高 ε_eff，必须重新提取而不能沿用原设计值。

金属间隙 S 的敏感性

固定 W = 10 μm、h = 500 μm、ε_r = 11.9，只改变 S：

（示例计算值。）S 对 Z₀ 的调节非常灵敏：S 从 6 μm 增到 10 μm，Z₀ 上升近 10 Ω，而 ε_eff 只变化约 0.12。这说明调阻抗主要靠改缝隙，调相速度/谐振频率主要靠改衬底与 ε_r。同时要注意，S 越小边缘电场越集中，电流密度在导带边缘的峰值越高，对超导薄膜的临界电流密度和微波损耗更不利；因此不能为了压低 Z₀ 一味缩小 S。

结果判读与验收条件

• 阻抗判据：提取的 Z₀ 与设计目标偏差在 ±2 Ω 内视为可接受初值；若用于与 50 Ω 链路直接对接，偏差超过 ±5 Ω 需调整 S。
• 介电常数判据：ε_eff 必须满足 1 < ε_eff < ε_r；若算出 ε_eff 接近 1，说明 k1 或椭圆积分比方向写错。
• 趋势自检：增大 S 应使 Z₀ 单调上升；增大 ε_r 应使 ε_eff 上升、Z₀ 下降；h 增大时参数应单调趋近半无限厚渐近值。任一趋势不符即为公式实现错误。
• 一致性判据：同一截面用解析公式与全波仿真提取的 Z₀ 差异通常应在几个百分点内，超出则检查金属厚度、接地平面宽度与辐射边界设置。