1. 量子基态制备从理论到离子阱实验的突破量子基态制备一直是量子多体物理和量子计算领域的核心挑战。传统方法如变分量子本征求解器(VQE)虽然被广泛使用但面临着参数优化困难、测量开销大等问题。最近RIKEN量子计算研究中心团队在Quantinuum的离子阱量子计算机Reimei上成功实现了一种基于耗散动力学的量子基态制备新方法。这项工作的创新之处在于首次在离子阱系统中验证了Ding等人提出的耗散基态制备协议通过Kraus表示扩展了Lindblad动力学框架使其适用于任意时间离散化步骤在含19个自旋的横场Ising模型中实现了稳定收敛即使量子电路包含多达4110个纠缠门结合零噪声外推(ZNE)技术显著提升了能量期望值的精度这项研究为NISQ(含噪声中等规模量子)设备上的量子-经典混合算法提供了新思路在量子化学模拟和强关联系统研究中展现出重要应用潜力。1.1 传统基态制备方法的局限在深入探讨这项突破性工作之前我们需要理解传统量子基态制备方法面临的挑战变分量子本征求解器(VQE)的瓶颈需要优化大量变分参数随着电路深度增加会出现贫瘠高原现象能量期望值的测量需要大量采样导致资源开销急剧增长对噪声敏感参数优化过程容易陷入局部最优量子相位估计(QPE)的限制需要容错量子计算机才能处理经典难解问题当前NISQ设备无法满足其相干时间和门保真度要求量子选择构型相互作用(QSCI)方法避免了VQE的参数优化问题但仍面临与经典选择CI方法比较优势不明确的问题这些限制促使研究者寻找更鲁棒、更适合NISQ设备的基态制备方案而耗散动力学方法正是这样一种有前景的替代方案。2. 耗散基态制备的理论框架2.1 核心思想与数学表述耗散基态制备协议的核心思想是设计一个完全正定且保迹(CPTP)的量子通道Γ_K使得目标哈密顿量Ĥ的基态|E₀⟩成为该通道的唯一稳态Γ_K[|E₀⟩⟨E₀|] |E₀⟩⟨E₀|为了实现这一点研究团队引入了一个辅助比特(ancilla)和跳跃算符K̂构建了如下的酉扩张Ŵ(√τ) exp(-iK̂√τ)其中K̂是扩展的跳跃算符K̂ ≡ [ 0 K̂† ] [ K̂ 0 ]关键设计在于跳跃算符K̂需要满足K̂|E₀⟩0即基态被K̂湮灭。通过这种方式基态自然成为量子通道的稳态。2.2 Kraus表示与保真度单调性研究团队的一个重要理论贡献是推导出了该耗散通道的Kraus表示Γ_K[ρ̂] M̂₀ρ̂M̂₀† M̂₁ρ̂M̂₁†其中M̂₀ cos(√(τK̂†K̂))M̂₁ -i√τ K̂ sinc(√(τK̂†K̂))这种表示的优势在于它适用于任意时间步长τ而不仅限于Lindblad动力学的小τ极限。更关键的是该协议保证了基态保真度的单调非减性。对于任意初始态ρ̂经过m次通道应用后F(Γ_K^{m1}[ρ̂], |E₀⟩⟨E₀|) ≥ F(Γ_K^m[ρ̂], |E₀⟩⟨E₀|)这意味着随着迭代次数的增加系统会越来越接近目标基态。注意虽然保真度单调增加但能量期望值不一定单调下降这与虚时间演化不同。这是耗散动力学的一个特点在实际应用中需要特别注意。2.3 跳跃算符与滤波函数设计跳跃算符的具体形式对协议性能至关重要。研究团队采用了算子傅里叶变换(OFT)表示K̂ ∫_{-∞}^∞ ds f(s)Â(s)其中Â(s) e^{iĤs}Âe^{-iĤs}是海森堡绘景下的算符演化。滤波函数f(s)的设计需要满足关键条件其傅里叶变换f̃(ω)在ω0时为零确保不引入能量增加的跃迁。研究中采用了费米-狄rac分布形式的滤波函数f̃(ω) n_F(β(ω-b)) - n_F(β(ω-a))其中a b 0β 0。当β→∞时这个函数趋近于理想的矩形窗函数。3. 离子阱实验实现细节3.1 量子硬件与模型系统实验在Quantinuum的20比特离子阱量子计算机Reimei上进行该系统具有以下特点单比特门平均保真度99.996%两比特RZZ门平均保真度99.86%状态制备和测量(SPAM)误差约0.35%研究团队选择了一维横场Ising模型作为测试平台Ĥ JΣ_{i0}^{N-2}Ẑ_iẐ_{i1} B_XΣ_{i0}^{N-1}Ŷ_i参数设置为J -1B_X -1.2研究系统尺寸N4,6,19。3.2 量子电路实现实验电路的核心是将连续的时间积分离散化处理。具体步骤包括截断积分区间为[-S_s, S_s]离散化为2M_s1个点使用二阶Trotter公式近似酉演化通过ancilla比特实现耗散通道的部分迹操作电路中的一个关键优化是移除了端点处的相干演化操作这在不影响结果的情况下显著减少了门数量。3.3 零噪声外推(ZNE)技术为了缓解硬件噪声的影响研究团队采用了零噪声外推技术通过门折叠方法增加噪声水平将每个RZZ(θ)门替换为[RZZ(θ)RZZ(-θ)]^{(G-1)/2}RZZ(θ)取G1,3,5三个噪声水平进行测量使用线性和指数两种拟合方式外推到G→0极限在N6的系统中最大电路包含1580个两比特门(G1,m20)通过ZNE技术成功将相对误差从约30%降低到与无噪声模拟一致的水平。4. 实验结果与性能分析4.1 能量收敛行为实验观察到了几个重要现象即使存在显著噪声系统仍能收敛到一个远离最大混合态的低能态噪声使稳态能量高于理想情况但远非完全混合态的零能量无噪声模拟也存在系统误差主要来自时间离散化近似图3展示了N4和N6系统的能量收敛曲线其中几个关键点值得注意硬件结果与噪声模拟器(Reimei-E)结果吻合良好收敛后的能量明显低于初始态能量(初始态为Y方向极化态E(0)0)即使电路包含上千个门协议仍表现出对噪声的鲁棒性4.2 大规模系统验证在N19的系统中(图5)研究团队实现了以下突破成功运行包含多达20550个两比特门的量子电路(G5,m30)观测到能量收敛行为与较小系统一致ZNE技术有效改善了能量估计值这一结果特别令人印象深刻因为按照两比特门保真度(0.9986)估算G1,m30的电路(4110个门)理论保真度仅为(0.9986)^4110≈0.003但实际获得的能量信号远好于这一预期表明协议具有内在的噪声鲁棒性。4.3 误差来源分析实验中的误差主要来自三个方面硬件噪声包括门误差、退相干和测量误差离散化误差时间积分截断和Trotter分解引入的系统误差滤波函数不完美有限的β值导致滤波函数边缘展宽可能引入能量增加的跃迁研究团队通过以下方式缓解这些误差采用ZNE技术校正硬件噪声优化离散化参数平衡精度与资源开销选择适当的β值确保滤波函数质量5. 应用前景与未来方向这项工作的意义不仅在于理论框架的提出更在于其实验实现展示了在现有量子硬件上解决实际问题的可行性。潜在应用包括量子化学模拟分子基态能量计算电子结构问题求解化学反应路径研究强关联系统研究Hubbard模型基态性质高温超导机制探索量子自旋液体研究算法发展结合变分量子本征求解器开发混合量子-经典优化协议探索早期容错量子算法未来研究方向可能包括扩展到更高维系统和更复杂相互作用开发更高效的滤波函数设计方法与其他误差缓解技术结合在超导量子处理器等其他平台实现这项研究为NISQ时代的量子计算应用开辟了新途径展示了即使在不完美的量子硬件上通过巧妙的算法设计也能获得有物理意义的结果。随着量子处理器性能的持续提升这类耗散协议有望在量子模拟和优化问题中发挥更大作用。