超导量子比特贝尔测试中的准备非平稳性漏洞解析

1. 超导量子比特贝尔测试中的准备非平稳性漏洞解析

量子计算领域中的贝尔测试是验证量子非局域性的重要工具，其核心原理基于CHSH不等式在经典与量子体系中的不同表现。在理想情况下，量子系统可以违反经典CHSH不等式（|S|≤2），达到理论最大值2√2≈2.828。然而，在实际的超导量子处理器上实施贝尔测试时，硬件的不完美性会引入各种系统性误差，其中准备过程的非平稳性（Preparation Nonstationarity）是一个长期被忽视的关键问题。

传统贝尔测试分析都基于一个核心假设：重复电路执行采样于单一平稳准备系综（single stationary preparation ensemble）。这意味着每次量子态制备都遵循完全相同的物理过程，可以用同一个完全正定保迹（CPTP）映射来描述。但最新实验研究表明，在当代超导量子硬件上，这个假设经常被违反。控制参数的慢时变（slow temporal drift）会导致准备过程产生上下文依赖的演化，即使保持测量独立性和局域性，也会影响观测到的贝尔违规值。

关键发现：当存在准备非平稳性时，经典CHSH边界应松弛为|S|≤2+6δens，其中δens量化了不同测量上下文采样的准备系综之间的最大统计差异。实验数据显示，在IBM量子处理器上观测到的贝尔违规值对应的δens与Hall型测量依赖模型所需的量级相当，但这里完全源于准备漂移与实验调度的结合。

2. 准备非平稳性的物理机制与量化框架

2.1 超导量子处理器中的非平稳性来源

在超导量子硬件中，准备非平稳性主要源于以下几类物理效应：

1. 低频控制参数波动：包括微波幅度/相位漂移、混频器失衡、直流偏置漂移、通量噪声和残余量子比特失谐等。这些波动的时间尺度通常在秒到小时量级，远长于单次电路执行时间（微秒量级），但短于完整实验持续时间。

2. 共享控制串扰：在多量子比特架构中，校准漂移和共享控制线之间的串扰会导致上下文依赖性——名义相同的操作会因电路结构或测量基的选择而表现出系统性偏差。

3. 环境参数漂移：制冷机温度波动、磁通噪声等环境因素会随时间影响量子比特的相干性和门保真度。

这些效应导致重复执行实际上是从一系列缓慢演化的准备映射中采样，而非单一的平稳过程。图1展示了这种非平稳性在时域上的表现，其中不同时间窗口内的量子态制备统计特性存在明显差异。

2.2 系综发散的理论量化

为了严格量化准备非平稳性，研究团队引入了系综发散参数δens，定义为不同测量上下文(a,b)和(a',b')对应的准备系综πab(λ)和πa'b'(λ)之间的最大总变差距离（total variation distance）：

δens = max_{(a,b),(a',b')} dTV(πab, πa'b')

其中总变差距离dTV(μ,ν) = (1/2)∫dλ|μ(λ)-ν(λ)|。基于此定义，可以推导出松弛的CHSH不等式：

|S| ≤ 2 + 6δens

这个结果表明，即使保持局域性和测量独立性，仅准备非平稳性就足以放宽经典贝尔边界。值得注意的是，δens与Hall提出的测量依赖参数M有本质区别——前者源于准备过程的时变，后者则涉及隐藏变量与测量设置之间的相关性。

2.3 操作漂移参数δop的测量方案

由于直接观测πab(λ)不可行，实验中采用操作漂移参数δop作为δens的可观测代理。对于固定测量上下文(a,b)，定义：

δop(ab) = max_{i<j} dTV(p_{ab}^{(i)}, p_{ab}^{(j)}) δop^global = max_{(a,b)} δop(ab)

其中p_{ab}^{(i)}表示第i个时间仓内的测量结果统计分布。根据总变差距离的收缩性，有δop ≤ δens，因此非零的δop即可证明准备非平稳性的存在。

实验设计的关键创新在于使用固定泡利轴测量（X、Y、Z）而非CHSH最优设置来评估δop。这是因为：

1. 泡利测量定义了正交、时间无关的经典通道，其结果可直接解释且对小误差稳定
2. CHSH最优设置涉及连续校准的模拟旋转，其漂移会与准备非平稳性混淆
3. IBM的校准协议将泡利测量作为固定参考基准持续监控

3. 实验设计与结果分析

3.1 硬件实验配置

实验在IBM的ibm_fez和ibm_torino超导量子处理器上进行，主要技术参数包括：

实验选用特定量子比特对（如ibm_torino上的(40,41)和ibm_fez上的(0,1)），通过标准Hadamard-CNOT序列制备单重态|Ψ⁻⟩=(|01⟩-|10⟩)/√2。采用两种执行调度方案：

1. 轮询调度：每个时间仓内执行所有四种测量上下文
2. 非平衡调度：每种上下文在连续时间仓中执行

每种配置使用B=6个时间仓，每个仓1024次测量，总测量次数达24,576次/设备。

3.2 操作漂移的观测结果

图2展示了从泡利轴测量提取的δop实验结果，关键发现包括：

1. 统计显著性：在轮询调度下，两个设备的δop均显著超过蒙特卡洛零模型预期（p<0.05至p<10^-3）
2. 调度依赖性：非平衡调度显著降低了ibm_fez的δop，但ibm_torino仍保持较高值，表明后者存在更强的固有非平稳性
3. 读数误差排除：经过双量子比特读数误差缓解后，δop仍显著存在，排除了经典测量误差主导的可能性

具体数值对比如下：

3.3 CHSH最优测量的对比研究

作为对照，实验也使用CHSH最优测量设置（一个量子比特沿0/π/2轴，另一个沿±π/4轴）评估贝尔违规值。关键发现：

1. 原始数据显示明显的操作漂移（δop≈0.15-0.20）
2. 但经过读数误差缓解后，δop降至与零模型一致的水平
3. 这表明CHSH最优设置不适合诊断准备非平稳性，因其对测量轴漂移过于敏感

图3对比了两种测量方案在误差缓解前后的δop变化，验证了泡利轴测量的必要性。

4. 对量子认证的影响与解决方案

4.1 准备非平稳性的实际影响

实验观测到的δop≈0.06-0.09对应松弛的CHSH边界2.36-2.54，这与实际测量的贝尔违规值范围（2.2-2.6）相当。这意味着：

1. 观测到的贝尔违规可能部分源于准备非平稳性，而非纯量子非局域性
2. 在NISQ设备上，传统CHSH测试的解释需要重新评估
3. 量子认证协议必须考虑这种准备依赖的漏洞

4.2 漂移感知的测试协议

基于此研究，我们建议采用以下改进方案：

1. 双重测量策略：

   • 使用泡利轴测量诊断准备非平稳性（δop）
   • 使用CHSH最优测量评估贝尔违规值
   • 联合分析确保|S|-2 > 6δop才有意义

2. 动态校准协议：

while True: # 监控泡利轴漂移 δop = measure_pauli_drift() if δop > threshold: recalibrate() # 触发重新校准 # 执行贝尔测试 S = measure_CHSH() # 验证是否满足严格量子违规 if S > 2 + 6*δop + margin: return "Genuine violation"

3. 时序优化方案：
   • 采用交错执行模式，最小化不同设置间的时间间隔
   • 实现"时间平衡"的CHSH测试，消除调度相关偏差
   • 结合实时漂移补偿算法

4.3 未来研究方向

这一发现开辟了几个重要研究方向：

1. 硬件层面：开发更稳定的控制电子学，降低低频噪声
2. 协议层面：设计对准备非平稳性鲁棒的贝尔测试变体
3. 理论层面：建立更完备的噪声模型，量化δens与门误差、相干时间等参数的关系

在实际操作中，我们建议研究团队定期监控泡利轴漂移作为硬件稳定性的诊断工具，并在解释贝尔测试结果时明确报告δop值。对于关键认证实验，应采用多次重复和交叉验证来区分真正的量子违规与准备伪影。

这项研究揭示了超导量子处理器上一个微妙但重要的系统误差来源，为提升量子基准测试的可靠性提供了新的方法论框架。随着量子硬件的不断发展，此类基础性研究将继续在确保量子优势声明的严谨性方面发挥关键作用。