量子自旋链耗散基态制备实验解析

1. 量子自旋链耗散基态制备实验解析

在量子多体系统研究中，基态制备是理解材料性质、模拟量子系统的关键步骤。传统经典计算方法如精确对角化、量子蒙特卡洛等面临指数墙问题，而近期发展的量子计算方法为解决这一难题提供了新思路。RIKEN团队在Quantinuum的离子阱量子计算机"Reimei"上成功实现了19个自旋的横向场Ising模型耗散基态制备，为量子模拟开辟了新路径。

1.1 耗散基态制备原理

耗散基态制备的核心思想是通过设计特定的量子通道，使系统在与环境相互作用过程中自然弛豫到目标基态。与需要主动优化参数的变分量子本征求解器(VQE)不同，这种方法利用了开放量子系统的动力学特性。

关键创新点在于构造满足以下条件的量子通道Γ_K：

1. 基态是通道的唯一稳态：Γ_K[|E₀⟩⟨E₀|] = |E₀⟩⟨E₀|
2. 保真度单调性：F(Γ_K[ρ], |E₀⟩⟨E₀|) ≥ F(ρ, |E₀⟩⟨E₀|)

实验采用Kraus算子表示法实现这一通道： Γ_K[ρ] = M₀ρM₀† + M₁ρM₁† 其中M₀ = cos(√τK†K)，M₁ = -i√τK sinc(√τK†K)

1.2 离子阱量子计算机优势

选择Quantinuum的"Reimei"离子阱系统主要基于以下考虑：

• 高保真度门操作：单量子门误差约0.004%，双量子门误差约0.14%
• 原生支持参数化ZZ门：RZZ(α) = exp(-iπαZ_iZ_j/2)
• 中电路测量与重置(MCMR)能力：这对实现耗散通道至关重要
• 可扩展性：系统包含20个量子比特，满足19+1(ancilla)的需求

2. 实验设计与实现细节

2.1 横向场Ising模型设置

研究采用一维开边界横向场Ising模型： H = JΣZ_iZ_{i+1} + B_XΣX_i 参数设置为J = -1，B_X = -1.2，对应铁磁相。系统尺寸测试了N=4,6,19三种情况。

基态制备的挑战在于：

• 传统量子相位估计需要容错量子计算机
• VQE面临barren plateau问题
• 量子-经典混合算法测量开销大

2.2 耗散通道构造

实验采用的关键技术是算子傅里叶变换(OFT)实现的跳跃算子： K = ∫ds f(s)A(s) 其中A(s) = e^{iHs}Ae^{-iHs}，A选择为第一个自旋的Z算符(Z₀)

滤波器函数设计为： f̃(ω) = n_F(β(ω-b)) - n_F(β(ω-a)) 参数设置为β=8/Δ，b=Δ/4，a=-2|E₀|，Δ为能隙

这种设计确保：

• 只允许能量降低的跃迁(ω ≤ 0)
• 避免激发态污染基态
• 参数β控制滤波器锐度

2.3 量子电路实现

整体电路架构如图1所示，主要包含三个部分：

1. 相干演化部分(灰色框)：实现e^{-iHΔt}的Trotter分解
2. 耗散通道部分(蓝色框)：实现W(√τ)的Trotter分解
3. 测量部分(白色框)：基旋转后测量

具体实现时：

• 时间积分截断为S_s=4π/(b-a)，离散化为M_s=4点
• 采用二阶Trotter公式分解W(√τ)
• 每步演化包含57(N=4)或79(N=6)个原生RZZ门

3. 实验结果与噪声处理

3.1 能量收敛行为

图2展示了N=6系统的能量随演化步数m的变化：

1. 无噪声模拟(蓝色圆圈)：快速收敛到接近基态能量
2. 硬件结果(红色三角)：受噪声影响但仍保持收敛趋势
3. 模拟器结果(橙色倒三角)：与硬件行为高度一致

值得注意的是，即使m=20时电路包含1580个RZZ门，按门误差估算的保真度仅约11%，但实测能量远优于完全混合态的零能量预期，显示协议对噪声的固有鲁棒性。

3.2 零噪声外推(ZNE)技术应用

为克服噪声影响，实验采用门折叠零噪声外推：

1. 噪声缩放：通过门折叠产生G=1,3,5三个噪声水平 RZZ(θ) → [RZZ(θ)RZZ(-θ)]^{(G-1)/2}RZZ(θ)
2. 外推方法：
   • 线性拟合：E = ãG + b̃
   • 指数拟合：E = ãexp(b̃G)
3. 结果评估：
   • N=6系统在m=20时，相对误差从30%降至与无噪声模拟一致
   • N=19系统在m=30时(含4110个RZZ门)，仍能获得合理外推结果

3.3 大规模系统验证

在19自旋系统中观察到：

1. 无噪声模拟与精确基态能量存在明显差距，反映离散化误差
2. 硬件结果仍显示收敛趋势，验证协议可扩展性
3. ZNE后能量显著改善，证明方法在NISQ时代的实用性

关键数据点：

• m=30时电路包含4110(G=1)至20550(G=5)个RZZ门
• 仍能获得有意义的基态能量估计
• 验证了耗散方法对大系统的适用性

4. 技术挑战与解决方案

4.1 滤波器函数实现难点

实验中滤波器函数的精确实现面临两个主要挑战：

1. 时间截断引起的边缘展宽：约π/S_s
   • 解决方案：选择S_s = O(β)以保持滤波器锐度
   • 误差随e^{-πS_s/β}衰减
2. 离散化引起的混叠效应：周期ω_alias = 2π/Δs
   • 解决方案：确保ω_alias - (b-a) > 2||H||
   • 误差随e^{-β(ω_alias-2||H||)}衰减

4.2 资源消耗优化

实验中的资源消耗主要来自：

1. 量子门数量：随系统尺寸N和演化步数m线性增长
   • N=6时每步需79个RZZ门
   • 通过电路优化(如利用native RZZ门)减少开销
2. 测量次数：每个数据点100-1000次采样
   • 采用重要性采样等技术可进一步降低

4.3 误差来源分析

主要误差来源及其影响：

1. 硬件噪声：
   • 门误差、测量误差等
   • 通过ZNE有效抑制
2. 离散化误差：
   • 时间积分截断和Trotter分解引起
   • 可通过增加Ms和减小Δt改善
3. 滤波器非理想性：
   • 有限β导致的ω>0泄漏
   • 严格满足˜f(ω)=0 for ω>0需要β→∞

5. 应用前景与扩展方向

5.1 在NISQ时代的应用价值

耗散基态制备协议在含噪中等规模量子设备上具有独特优势：

1. 对噪声的固有鲁棒性：噪声态仍保持低能量特性
2. 避免参数优化：克服VQE的barren plateau问题
3. 资源效率：仅需1个辅助量子比特
4. 与错误缓解技术兼容：如ZNE、误差校正等

5.2 潜在扩展方向

基于当前成果，未来可能的发展包括：

1. 更复杂模型：如Hubbard模型、量子化学Hamiltonian
2. 混合协议：结合VQE优化耗散通道参数
3. 动态特性研究：非平衡态到基态的弛豫过程
4. 硬件改进：利用更高保真度门操作提升精度

5.3 与传统方法对比

与主流基态制备方法的比较：

实验结果表明，耗散方法在NISQ设备上实现了规模(19自旋)和精度(通过ZNE)的双重突破，为量子模拟提供了新范式。