来源论文: https://arxiv.org/abs/2501.13371 生成时间: Feb 24, 2026 14:21

0. 执行摘要

在 NISQ（有噪声的中等规模量子）时代，如何在有限的量子比特数和较浅的电路深度下准确描述电子强关联系统，是量子化学领域的核心挑战。传统的变分量子特征值求解器（VQE）受限于硬件高效模拟算符（HEA）的表达能力，往往难以捕捉复杂的电子相关性能；而基于化学启发式的酉耦合簇（UCC）方法虽然准确，但电路深度过大，超出了当前硬件的相干时间。

近日，来自芝加哥大学（University of Chicago）、阿贡国家实验室（Argonne National Laboratory）及 NVIDIA 的研究团队提出了一种名为 LAS-nuVQE 的新框架。该工作通过以下三大支柱实现了突破：

  1. 理论创新：结合局部活跃空间自洽场（LASSCF）作为初始态，利用 Jastrow 启发的非酉（Non-unitary）算符在不增加电路深度的前提下补偿电子相关能。
  2. 物理约束：引入自旋惩罚项，解决了 HEA 在对称性保持方面的固有缺陷，确保了计算结果的物理含义。
  3. 工程优化：通过 Pauli 分组和“吝啬采样”（Shot-frugal sampling）技术，将测量开销降低了两个数量级，使得 100 量子比特规模的计算在墙钟时间内变得可行。

该研究在 $H_4$ 链和环丁二烯（Cyclobutadiene）系统上进行了基准测试，结果证明 LAS-nuVQE 能够以不到 70 个量子门达到化学精度，显著优于传统的 HEA-VQE 方法。

1. 核心科学问题，理论基础，技术难点与方法细节

1.1 核心科学问题：精度与深度的权衡

量子化学模拟的核心在于求解分子哈密顿量的基态能量。对于强关联系统（如过渡金属配合物、断键过程），单参考态方法（如 Hartree-Fock）失效，必须采用多参考态方法。然而，在量子计算机上实现高精度的多参考态模拟（如 CASSCF）面临指数级增长的资源需求。LASSCF 方法通过将活跃空间分割为相互独立或弱相互作用的碎片（Fragments），在经典计算机上实现了线性缩放，但在处理碎片间的相关性时精度受损。如何通过量子算法以低廉的代价补足这部分能量，是本论文解决的首要问题。

1.2 理论基础：LASSCF 与 nuVQE 的协同

1.2.1 局部活跃空间自洽场（LASSCF）

LASSCF（Localized Active Space Self-Consistent Field）是 CASSCF 的近似方案。它将总活跃空间划分为若干个局部碎片活跃空间。其波函数定义为各个碎片波函数的反对称乘积（Wedge Product）：

$$ |\text{LAS}\rangle = \bigwedge_K (\Psi_{A_K}) \wedge \Phi_D $$

其中 $\Psi_{A_K}$ 表示第 $K$ 个碎片的波函数。LASSCF 的优势在于它忽略了碎片间的高阶关联，从而大幅降低了计算量。在 LAS-nuVQE 中，研究者首先在经典计算机上完成 LASSCF 计算，将其结果作为量子电路的初始参考态。

1.2.2 非酉 VQE (nuVQE) 的数学原理

传统的 VQE 依赖于酉算符 $U(\theta)$。nuVQE 则引入了一个基于 Jastrow 因子的非酉算符 $\hat{O}(\vec{\lambda})$。其试探波函数表示为：

$$ |\Psi^o(\vec{\theta})\rangle = \hat{O}(\vec{\lambda})\hat{U}(\vec{\theta})|\Psi_0\rangle $$

由于 $\hat{O}$ 是非酉的，波函数不再归一化。因此，能量期待值的计算公式变为：

$$ E = \frac{\langle \Psi(\vec{\theta}) | \hat{O}^\dagger \hat{H} \hat{O} | \Psi(\vec{\theta}) \rangle}{\langle \Psi(\vec{\theta}) | \hat{O}^\dagger \hat{O} | \Psi(\vec{\theta}) \rangle} $$

通过这种构造，我们不需要在量子电路上真正实现非酉操作（这通常需要辅助比特），而是通过在经典后处理中测量 $\hat{O}^\dagger \hat{H} \hat{O}$ 和 $\hat{O}^\dagger \hat{O}$ 两个算符的期待值来获得能量。这种方法的神奇之处在于，它通过改变测量目标，在不增加量子电路深度的情况下，人为地增加了波函数的参数化空间。

1.3 技术难点：算符的线性化与对称性缺失

  • 算符复杂度：Jastrow 算符通常具有指数形式 $\exp(\sum Z_i Z_j)$。直接测量该算符极其困难。本文采用了线性化近似： $$ J(\vec{\alpha}, \vec{\lambda}) = 1 - \sum_i \alpha_i Z_i - \sum_{i
  • 对称性破缺：HEA（硬件高效拟合）的一个致命伤是它不保证自旋和粒子数守恒。对于环丁二烯这种对自旋极度敏感的分子，HEA 往往会产生物理上错误的态。研究团队引入了自旋惩罚项 $\mu_C \langle (\hat{S}^2 - c)^2 \rangle$ 到代价函数中，强行将优化过程导向正确的自旋多重度空间。

1.4 方法细节：LAS-nuVQE 算法流程

  1. 经典预处理：运行 LASSCF 获得碎片轨道和碎片相互作用的有效哈密顿量。
  2. 电路初始化：在量子计算机上利用直接初始化（DI）加载 LASSCF 态，并构建 HEA 电路层。
  3. 算符映射：将算符 $\hat{J}^\dagger \hat{H} \hat{J}$ 映射为大量 Pauli 算符的加权和。
  4. 混合优化：在经典优化器（如 L-BFGS-B）中，同时调节 HEA 的转动角 $\theta$ 和 Jastrow 参数 $\alpha, \lambda$。

2. 关键 Benchmark 体系与数据分析

2.1 体系选择：从 H4 到环丁二烯

研究者选择了两个具有代表性的体系来验证算法：

  • $H_4$ 链：通过调节几何构型，模拟碎片间关联从弱到强的演变。构型 1 具有较强的碎片间关联（LASSCF 表现差），构型 2 关联较弱。
  • 正方形环丁二烯 ($C_4H_4$)：典型的抗芳香性分子，具有两个简并的单电子轨道，是单参考态方法的“噩梦”。

2.2 核心性能数据

2.2.1 $H_4$ 系统的精度突破

在 $H_4$ 构型 1 中，LASSCF 的原始误差超出了化学精度（1.6 mEh）。使用常规的 LAS-VQE（只有酉操作），在 2 层电路深度下，能量误差依然显著。然而，LAS-nuVQE 在 2 层电路下即达到了化学精度，总门数少于 50 个。这证明了非酉算符能够有效捕获 LASSCF 丢失的碎片间电子相关能。

2.2.2 环丁二烯的相关能回收

对于 $C_4H_4$（8 量子比特，cc-pvdz 基组），研究测量了“碎片间相关能回收率”：

  • LAS-VQE：在 4 层电路时仅能回收约 92.5% 的相关能。
  • LAS-nuVQE：在 2 层电路时回收率即达到 95% 以上，4 层以上时回收率超过 99.9%。 这意味着对于复杂的强关联体系，nuVQE 具有极高的参数效率。

2.2.3 自旋约束的效果（Table 1 & 2）

未经约束的 HEA 计算显示，其自旋期待值 $\langle S^2 \rangle$ 偏差巨大（从 0.67 到 1.52 不等，理想值为 0）。引入惩罚项后，$\langle S^2 \rangle$ 降至 $10^{-4}$ 到 $10^{-9}$ 数量级，且能量依然保持在化学精度范围内。这解决了 HEA 在实际应用中的可靠性问题。

3. 代码实现细节与复现指南

3.1 软件包生态

该研究建立在成熟的开源生态之上：

  • 量子模拟器：Qiskit (version 0.46) 及其 Aer 模块（statevector_simulator 和 qasm_simulator）。
  • 经典化学计算：
    • PySCF：用于 Hartree-Fock 和基础轨道计算。
    • mrh：由 Matthew Otten 开发的专门用于 LASSCF 和多参考态计算的工具包。
  • 优化算法：SciPy 中的 L-BFGS-B，适合处理带约束的大规模参数优化。

3.2 关键实现逻辑

  1. 哈密顿量转换：利用 Jordan-Wigner 变换将费米子算符转为 Pauli 算符。
  2. Jastrow 构造：Jastrow 算符的线性化形式需要在经典侧维护一个 Pauli 算符池。由于 $\hat{J}$ 只包含 $Z$ 和 $I$，$\hat{J}^\dagger \hat{J}$ 依然是全对角的，这大大简化了分母的测量。
  3. Pauli 分组逻辑：为了缓解 $O(N^8)$ 的测量压力，必须实现 QWC（Qubit-Wise Commutativity）或 FC（Full Commutativity）分组算法。通过寻找图的最大团（Clique Cover），将数千个待测项合并为几十个测量任务。

3.3 开源资源

研究团队秉持开放科学的精神，已将核心代码托管至 GitHub：

  • 仓库链接：https://github.com/joannaqw/LAS-nuVQE
  • 复现建议：建议从 $H_2$ 分子开始，使用 STO-3G 基组熟悉 nuVQE 的非归一化测量逻辑，再逐步扩展到 LAS 碎片化场景。

4. 关键引用文献与局限性评论

4.1 关键参考文献

  1. nuVQE 起源：Benfenati et al. (J. Chem. Theory Comput. 2021, 17, 3946) 首次提出了在噪声设备上利用非酉 VQE 提升精度的想法。
  2. LASSCF 理论：Hermes and Gagliardi (J. Chem. Theory Comput. 2019, 15, 972) 奠定了局部活跃空间方法的基础。
  3. 测量优化：Verteletskyi et al. (J. Chem. Phys. 2020, 152, 124114) 提供的最小团覆盖理论是本文测量减缓技术的物理依据。

4.2 局限性分析（技术评论）

尽管 LAS-nuVQE 表现优异，但在迈向生产级应用前仍有几点局限：

  • 经典优化器压力：nuVQE 引入了额外的参数（$N_q$ 个一体项和 $N_q(N_q-1)/2$ 个二体项）。随着比特数增加，经典优化器在非凸能量景观中寻找全局最优解的难度会剧增，可能会遭遇“贫瘠高原”（Barren Plateaus）问题。
  • 非线性近似的局限：线性化 Jastrow 算符虽然方便测量，但在极强关联区域（如多重键断裂中心），可能需要更高阶的展开项才能保持精度，这会重新拉高测量成本。
  • 测量统计误差：尽管有吝啬采样技术，但在实际硬件上，非归一化波函数的能量估算对分母（归一化因子）的精度极度敏感。如果分母由于噪声偏小，会导致能量值出现剧烈波动。

5. 补充讨论：算法的可扩展性与未来展望

5.1 墙钟时间预测：通往 100 量子比特之路

论文最引人注目的部分之一是对大规模系统的资源预测（Figure 9 & 10）。

  • 原始 nuVQE：对于 100 量子比特，单次能量查询需要惊人的 3000 万年（在离子阱平台上）。
  • 经过缓解的 nuVQE：通过 FC 分组和吝啬采样，超导平台的计算时间可缩短至 28 小时以内。这意味着该算法在现有的“准容错”或高质量 NISQ 硬件上已经具备了开展大规模真实分子研究的可能性。

5.2 对比 HEA 与 UCC

LAS-nuVQE 实际上开辟了第三条道路：它既不像 HEA 那样盲目搜索（通过 LASSCF 提供了物理良好的起点），也不像 UCC 那样受困于电路深度（通过非酉后处理转移了计算负担）。这是一种典型的“以空间换深度”或“以测量换门控”的策略。

5.3 未来研究方向

  1. 噪声鲁棒性研究：本文主要基于模拟器。在真实量子芯片上，退相干对非酉算符测量精度的影响需要进一步量化。
  2. LUCJ 算符的应用：将局部酉耦合簇 Jastrow（LUCJ）与 LAS 结合，可能会进一步压缩参数量。
  3. 自动活跃空间选择：结合机器学习自动划分 LAS 碎片，将使该方法从专家级工具变为通用化软件。

总结：LAS-nuVQE 的出现标志着变分量子算法正在从“概念验证”走向“高效实用”。通过深入理解电子结构的局部性，并利用非酉数学技巧规避硬件深度的限制，我们离在量子计算机上解决真正的工业级化学难题又近了一步。