多cell超导腔瞬态行为研究
在大多数对多单元超导加速结构瞬态行为的仿真中，腔体往往被简化表示为一个由离散电阻 $R$、电感 $L$ 和电容 $C$ 组成的单一谐振器。然而在实际中，构成腔体的各个单元应当被视为彼此耦合但相互独立的谐振器，每个单元都有其特定的谐振频率失谐、$Q$ 值和并联阻抗 $R$。本文对多单元腔进行了全面分析，并讨论了腔体制造缺陷对其时间相关行为的影响。
1. 理论模型
我们构建了一个基于并联 RLC 谐振器电路的多单元腔模型。各个单元中依次流过电阻、电感和电容的电流之和为零。在此等效下，相邻的单元通过电容 $\hat C$ 相互耦合，并由输入阻抗为 $Z_0$ 的发生器驱动：
图1：基于并联电路 RLC 谐振器的多单元腔，通过电容 $\hat C$ 相互耦合，并由发生器驱动。
具有内阻 $Z_0$ 的发生器提供电流，并将其送入腔体。输入功率通过一个输入耦合器注入第一个单元中。由于电流边界在谐振时可近似解为慢变幅度 $V_j=\tilde V_je^{i\omega t}$ 的微分组合，这为后续的理论推导奠定了分析前提。第一和最后一个单元还受到与束流管的耦合影响，用电容 $C_a$ 和 $C_b$ 来表示。
2. 稳态分析
在讨论完整的时变问题之前，我们先观察无驱动条件下的裸多单元腔的稳态情况，从而确定束流管电容大小以及腔体本征模频率。
图2：未加驱动的多单元腔，保留与束流管耦合的电容 $C_a$ 和 $C_b$。
假设系统以频率 $\omega$ 振荡，可以用阻抗 $Z_C=1/i\omega C$ 和 $Z_L=i\omega L$ 分别表征电容和电感。对于处于中部的单元节点 $V_j$ 进行电流平衡列式可得：
$$
0 = I_j - I_{j-1} + I_{LC} = i\omega\hat C(V_j - V_{j-1}) - i\omega\hat C(V_{j+1} - V_j) + \left(i\omega C + \frac{1}{i\omega L}\right) V_j
$$
整理后，引入缩写 $\kappa=C_a/C$ 以及 $\hat\omega^2=1/LC$，可以得到类似于如下形式的本征值方程（以 $n=5$ 单元为例）：
$$
\left(
\begin{array}{ccccc}
1+3\kappa & -\kappa & 0 & 0 & 0 \
-\kappa & 1+ 2\kappa & -\kappa & 0 & 0 \
0 &-\kappa & 1+ 2\kappa & -\kappa & 0 \
0 &0 &-\kappa & 1+ 2\kappa & -\kappa \
0 & 0 & 0 & -\kappa & 1+3\kappa
\end{array}
\right)
\left(\begin{array}{c} V_1\ V_2\ V_3\ V_4\ V_5\end{array}\right)
= \frac{\hat \omega^2}{\omega^2}
\left(\begin{array}{c} V_1\ V_2\ V_3\ V_4\ V_5\end{array}\right)
$$
2.1 模态分布
在 $\kappa=0.01$ 时，这五个本征模的电压分布和相应的数值本征值如图3所示。其中最大的本征值对应着 $\pi$ 模，这里即对于 $\pi$ 模，满足要求条件 $\omega_\pi = \frac{\hat\omega}{\sqrt{1+4\kappa}}$。
图3：$\kappa=0.01$ 时五单元腔的模态（左）与相应本征值（右）。
如果用各自的频率激发不同模式，组合后的激励分布就可以被控制为在不同的特定单元中达到峰值，例如可以在特定单元中局部激发等离子体来进行清洗。
2.2 对单元耦合误差的敏感度
研究某一个单元耦合误差对理解稳态失谐有很大启发。下图展示了在第二与第三单元之间的耦合误差 $\Delta\kappa_2$ 对本征值及 $\pi$ 模相对频率偏移的影响。
极微小的耦合变动也会极大地影响超导腔极为狭窄的带宽，即使 $\Delta\kappa_2=10^{-4}$ 的误差都可能将频率偏移出腔体的固有带宽。
3. 瞬态现象理论建模与推导
为了将静态方程扩展至时变系统，我们需要代入慢变幅度，建立完整的矩阵微分方程。假设由电流产生并激励出特定主谐振分布，推导所得时域下对应的矩阵方程为：
$$
D\frac{d\vec V}{dt}\Delta t = -(A-iB)\vec V +\frac{\hat\omega\Delta t R_1}{Q_1}\vec I
$$
其中矩阵 $D$ 用于耦合与调整各项谐振偏移，矩阵 $A$ 包含阻尼有关（与 $Q$ 相关的常数），矩阵 $B$ 由单元间的电容耦合系数 $\kappa$ 构造得出。为实际数值解法，进一步将其重组为包含实部、虚部的双倍维度实矩阵。当此矩阵经过本征分解（如利用 eig() 计算其时间步长为 $\Delta t$ 时对应的矩阵指数 $\exp\left[-\tilde A s\right]$）后，我们可以通过解常微分方程组求出发生器开启后不同标度时间 $s$ 的演化规律。
当驱动频率设定为 $\pi$ 模激励时，仿真前几微秒的启动行为呈现出特征：
可以观察到，极短时间内各单元场依次上升出现滞后；但在 0.2 μs 后就达到了整体一致上升趋势。腔体在短时间尺度上实现同步驱动。
长脉冲仿真中发电机在 $t=0$ 开启并在 $t=3.5$ ms 时关闭。此时各单元电压差异微乎其微。发射停止时，能观察到明显的相控反射尖峰。
为更清晰地揭示腔体的能量吸收表现，可以考察上述时间演化矩阵的 10 个本征值：复数本征值的实部决定了对应宏观模态的时域放电（或填充）常数，虚部则主导在脉冲建立沿上的振荡高低。编号 9 与 10 的本征值（对应的正是 $\pi$ 模）具有最小的虚部与适中的实部，因此它们不表现为强振荡，而是极有效地吸收发生器能量。填充这种多单元的本质特征就是时间常数会被放大（例如 5 单元需要多花费 5 倍时间），带宽也将按比例缩小。
对比多单元模型与单单元腔的电压解析解，在时间常数匹配（即 $Q$ 放大约五倍以弥补级联延时）的前提下两者包络高度吻合，这充分证明了在大多数场景下采用单单元等效模型近似多单元正常腔室是合理的。
*图7：单单元腔中的电压（品红）叠加在五单元腔中的电压（黑色）之上。两者呈现极高的一致性。*
4. 腔体缺陷与公差影响分析
在稳健的矩阵微分方程基础上，我们进一步从理论中引入反映各类制造缺陷或运行问题的参数（$\Delta\kappa_i$, $\Delta\omega_i$, $Q$ 改变等），来探究这些失常现象在时域上的外在表现。这对于测定真实腔体各类公差指标具有非常高的工程利用价值。
4.1 单元间耦合误差的影响
若加工过程中出现误差，使得单元之间具体的偏耦合项 $\kappa_i$ 偏离设计值 $\kappa$。当把修正带有 $\Delta\kappa$ 的偏置量代入 $B$ 和 $D$ 矩阵之后，系统本征值会带有一定的虚部，引发额外的腔内电压透射和反射振铃（Ringing）现象。
当改变中间单元的耦合 $\Delta\kappa_2$ 时，可以看到 $\pi$ 模出现了类共振衰减，此时共振曲线较窄；而若是改变前后端与束流管耦合的误差 $\Delta\kappa_a$，共振响应曲线则明显变宽（大约宽3倍）。这说明系统对腔体内部单元间的耦合加工公差更为敏感，而两端的制造公差则可以适度放宽。
4.2 本征品质因数 $Q$ 衰减与失谐（Detuning）
当特定单元的电阻耗散突然增加，即导致了特定部分 $Q$ 值下降或“热失控”。如果单一单元的品质因数只是温和降低，则对整体电压的影响往往不明显。
单元失谐影响：我们还针对谐振频率的失配进行了分析，引入频率偏离项 $\delta_j = 2\Delta\omega_j / \hat\omega$ 作用于矩阵方程式。当所有单元以相同量失谐时，系统表现出极高的敏感性（频率响应骤降）；而若是仅有某个单一单元失谐（如单元4失谐量为 $\delta_4$），对系统整体透过率和反射率波形的影响则相对缓和很多，只有轻微的振幅虚部产生。因为在正常驱动中，$\pi$ 模作为一个相干整体在运转，“平均”并掩盖了微弱的个体失谐特性。
4.3 腔室淬灭（Quench）
如果因耗散致使发生腔室局部发热与电阻损耗剧增发生极大的淬灭，导致 $Q$ 发生极大下降：
只要热量损耗仍在这个单元可以通过输入耦合弥补的阈值（对应于耦合比 $\beta$）内，该淬灭是可以被稳住的。然而，如果因淬灭造成等效电阻极剧减小（例如因子达 $x=10^5 \gg \beta$），将无法补充耗散能量，最终直接导致 $\pi$ 模甚至整个建立的超导腔电场近乎完全崩溃。
5. 结论
通过电路模型和微分矩阵方程对多单元并联耦合等效结构的推导，我们可以成功重建超导腔内的不同模式稳态，并分析高精度的时间相关充放电过程。
- 瞬态主导与整体近似：除了高频瞬态开启部分会涉及多个本征杂模，$\pi$ 模主导了工作时的电磁场动力学作为相干整体振荡，因此只要时间常数匹配，单单元谐振器非常适合用作物理评估时的宏观模拟替代工具。
- 制造公差依赖性探究：$\pi$ 模的整体性导致多数缺陷对输出性能展现出“平均”效用。基于此方法推导，诸如针对单元内耦合系数及个别 $Q$ 值的衰退，可以明确不同位置元件针对失配的允许误差范围，为超导制造和诊断检测（截取兆赫兹高频采样）提供了极佳的理论支撑。
该计算框架有助于更深入理解超导腔在实际加工和工况环境中的物理特性，对于后续研究并分析热淬灭和制程不均匀性具有重要的指导意义。