Probing the Ground State of the Antiferromagnetic Heisenberg Model on the Kagome Lattice using Geometrically Informed Variational Quantum Eigensolver

技术摘要：利用几何启发式变分量子特征值求解器探测卡格米（Kagome）晶格反铁磁海森堡模型的基态

问题陈述
本文旨在解决寻找卡格米晶格反铁磁海森堡模型基态的挑战，该系统以其几何挫折和潜在的量子自旋液体相而闻名。确定此类 $k$-局部哈密顿量的基态是一个 QMA-hard 问题，由于希尔伯特空间的指数级增长，使得大规模系统的经典模拟变得难以实现。虽然变分量子特征值求解器（VQE）提供了一种适用于含噪声中型量子（NISQ）设备的混合量子-经典方法，但它们面临着显著的障碍：需要深层电路、易受噪声（相干误差和非相干误差）的影响、梯度消失导致的贫瘠高原（barren plateaus）问题，以及设计既具表现力又符合硬件效率的拟设（ansatz）的困难。具体而言，作者旨在通过无需依赖重度特定问题简化或在优化循环中进行广泛误差缓解的方式，在真实的量子硬件上表征基态性质，如自旋-自旋相关性和静态结构因子。

方法论
作者提出了一种针对卡格米晶格定制的几何启发式 VQE 方法，重点关注两种基本的几何结构：单个三角形和一个卡格米星形结构（六个角共享的三角形）。

1. 几何启发式拟设：
   其核心创新在于一种定制的、浅层的、硬件高效的拟设。与过于深层的哈密顿量变分拟设（HVA）不同，该拟设通过一种由参数化单比特旋转（$R_y$）和受限于原生硬件连通性的纠缠门构成的模式来构建（避免使用 SWAP 门）。

   • 欧几里得参数空间： 该设计的一个关键特征是其 Fubini-Study 度规张量在解析上被证明是角对角且恒定的。这确保了参数空间自然呈现欧几里得特性，无奇异性且不存在过度参数化。因此，量子自然梯度等价于标准梯度，从而简化了优化景观。
   • 硬件兼容性： 电路深度被最小化（一层纠缠层），并通过交替放置 $R_y$ 门减少了参数数量。

2. 优化策略 (AQNGD/AGD)：
   作者采用了自适应梯度下降（AGD）优化器。由于该拟设产生了一个对角且恒定的 Fubini-Study 度规，因此自适应量子自然梯度下降（AQNGD）简化为 AGD。学习率使用 Armijo 回溯规则进行动态调整。这允许使用较大的初始步长以加速收敛，同时在优化器接近极小值时自动减小步长，从而避免在狭窄几何结构中停滞。将其与同时扰动随机逼近法（SPSA）进行了对比。

3. 误差缓解（后处理）：
   误差缓解是在 VQE 优化之后应用的，用于精炼最终的能量估计和结构性质，而不是在迭代训练期间进行。

   • 读取误差缓解 (REM)： 使用分块响应矩阵方法来反转状态准备与测量（SPAM）误差，对于稀疏测量而言，这种方法比全矩阵求逆更具扩展性。
   • 零噪声外推 (ZNE)： 作者使用二次外推（以及贝叶斯多项式回归）应用 ZNE，以将结果外推至零噪声极限。作者指出，尽管 ZNE 并不严格遵守 Rayleigh-Ritz 变分原理，但在此处被用于提高精度。

关键结果
研究在真实量子硬件（IBM Oslo, IBM Kyoto, IBM Torino, 和 IBM Algiers）上进行，并通过状态矢量模拟进行了验证。

• 收敛与优化：

  • 在 3-比特三角形系统中，经过 AGD 优化的 VQE 平滑地收敛至基态能量（$E \approx -3$），使用仅 19 次迭代即可达到 99.89% 的高保真度（在模拟中）。最终参数与理论精确值高度吻合。
  • 在 12-比特卡格米星形结构中，AGD 展示了优于 SPSA 的收敛性能，只需较少的迭代次数即可达到基态能量（$E \approx -18$）。在模拟中的最终保真度为 96.37%。
  • 通过 Armijo 规则进行的步长自适应被证明能有效导航代价景观，通过从大的步长开始并随着收敛接近而不断细化步长。

• 结构性质的噪声韧性：

  • 一个重要的发现是，自旋-自旋相关性和静态结构因子 ($S(q)$) 对噪声具有韧性。即使没有误差缓解，实验测得的 $S(q)$ 图谱仍保留了基态的定性特征（例如，零动量处的抑制表明不存在铁磁序）。
  • 噪声实验 $S(q)$ 与精确理论 $S(q)$ 之间的皮尔逊相关系数保持极高（>99.9%）。
  • 模拟实验揭示了静态结构因子对噪声水平的“阶梯式”依赖关系：相关结构在达到临界噪声阈值之前保持稳健，之后会迅速退化，类似于相变。这与能量随噪声单调衰减的行为形成对比。

• 误差缓解效能：

  • 应用 REM 和 ZNE 后处理显著改善了能量估计。对于卡格米星形结构，原始能量约为 $\approx -14.9$，仅使用 REM 可将其提升至 $\approx -17.3$，而 ZNE（二次）则使其更接近 $-18$ 的精确值。
  • 作者观察到，结合使用 REM 和 ZNE 有时会导致能量越过基态（违反变分原理），这一行为与 ZNE 的非变分性质一致。在测试配置中，ZNE 通常提供了最佳的准确度。

意义与主张
本文声称，一种几何启发式的拟设设计（确保了天然的欧几里得参数空间）使得在当前的 NISQ 硬件上进行高效的 VQE 训练成为可能，而无需复杂的特定问题简化。

• 稳健性： 本工作证明，可以从带有噪声的量子硬件中提取有意义的物理性质（特别是自旋-自旋相关性和静态结构因子），并获得足够的准确度来表征基态，即使在优化循环期间不进行误差缓解也是如此。
• 优化效率： 使用带有自适应步长的 AGD 被证明比 SPSA 对于此类硬件高效拟设更有效，在迭代次数方面实现了更快的收敛。
• 可扩展性： 通过避免深层电路和训练期间沉重的误差缓解，该方法为探测受挫量子磁体提供了一条可扩展的路径。作者建议该框架可以扩展到更大的卡格米晶格和其他受挫量子系统。
• 谦逊性： 作者承认，虽然其方法成功表征了小型基本单元（三角形和星形），但无限卡格米晶格的一般基态在文献中仍是一个开放性问题。他们并未声称解决了无限晶格问题，而是展示了一种在当前设备上探测此类系统的可行方法论。他们还指出，虽然 ZNE 提高了能量估计的精度，但它并不保持变分原理，且结合使用 REM 和 ZNE 需要谨慎解读。