强关联分子的量子模拟快速初态制备方法

对基态能量估计的量子算法研究通常假设初始态与真实基态完全重合;然而在实际中,初始态与真实基态之间往往存在不完美的重叠。论文Rapid Initial-State Preparation for the Quantum Simulation of Strongly Correlated Molecules从两个方面应对这一问题:一是更快速地制备矩阵积态(MPS)近似态,二是通过更高效的“过滤”技术从已制备的状态中提取基态能量。该论文的研究为在复杂化学系统中制备高重叠度初态提供了切实可行的路径示例。

论文作者为Dominic W. Berry , Yu Tong , Tanuj Khattar , Alec White , Tae In Kim , Guang Hao Low , Sergio Boixo , Zhiyan Ding , Lin Lin , Seunghoon Lee , Garnet Kin-Lic Chan, Ryan Babbush 和 Nicholas C. Rubin,来自Macquarie University, CalTech, Duke University, Google Quantum AI, Quantum Simulation Technologies Inc., Seoul National University, UC Berkeley。

一、研究背景与问题动机

在量子计算中,通过量子相位估计(Quantum Phase Estimation, QPE)求解化学和材料系统的基态能量是一项关键任务。传统分析常假设初始态已完美为基态,这大大简化了分析。但在现实中,实际制备出的初始量子态通常只能与基态有有限重叠,而这种不完美的初始态将直接影响QPE的成功概率和最终的能量估计精度。因此,为实际系统开发一种高效且具有良好重叠度的初态制备方案显得尤为重要。

此外,即便初态有一定重叠度,如何进一步“过滤”该状态以准确提取基态信息也是挑战之一。本论文从两方面入手解决此问题:一是设计更高效的矩阵积态(Matrix Product State, MPS)制备方案;二是提出更优化的“过滤”策略用于提取基态能量。

二、主要贡献与方法概述

本论文的贡献主要体现在两个方面:

  1. 更高效的MPS制备方法:借助改进的单元算子合成策略,显著降低了Toffoli门复杂度(非Clifford门,是评估容错量子计算资源的关键指标),在合成MPS时Toffoli门数量比以往方法减少约7倍。
  2. 两种过滤策略用于提升基态提取效果:
    • 基于采样的方法:通过QPE多次采样,取最小值作为基态能量估计。
    • 基于二分搜索与振幅估计的方法:利用量子振幅放大机制,对未知基态能量区间进行二分搜索,从而有效规避低重叠导致的低采样概率问题。

此外,研究团队还将方法应用于含铁硫簇系统(如FeMo辅基)的资源评估,发现即使在中等的MPS bond dimension(如4000)下,初态与真实基态的重叠度平方约为0.9,显示出较高的实用潜力。

三、MPS制备算法的技术创新

研究在MPS制备方面的关键技术创新是提出一种新的单元算子合成方法,基于矩阵对角相位变换(diagonal-phasing)与Hadamard变换交替执行,并结合QROM(Quantum Read-Only Memory)进行参数加载,从而极大降低了Toffoli门数量。具体技术包括:

  1. 将通用干涉仪分解为若干Givens旋转构成的二维块对角矩阵。
  2. 利用光学系统的“50:50分束器+相位”的组合思想将这些块简化为标准电路结构。
  3. 采用增量器与QROM输出机制实现高效控制与数据加载。

该算法还支持对MPS子空间的部分列合成(只合成需要的子矩阵列),进一步节省计算资源。

四、相位估计中的窗口函数与误差控制

为了更高效地从初态中提取基态能量,研究中引入了经典信号处理中的“窗口函数”理论用于QPE的相位误差控制:

  1. 传统QPE方法(如均匀幅度的|+⟩态)存在高方差与长尾分布,容易产生偏差。
  2. 本文采用Kaiser窗口和Slepian窗口,两者能以更小的尾概率抑制异常低估的能量值,提高估计稳定性。
  3. 对应的误差项被分类为:
    • 类型I误差:估计值小于真实基态能量的概率。
    • 类型II误差:估计值来自激发态但未能识别,从而导致高估。
      通过细致分析这两类误差的权衡关系,研究发现Kaiser窗口在实际中常优于Slepian窗口。

五、采样法与二分搜索法的比较

研究详细比较了直接采样法与基于振幅估计的二分搜索法在不同初态重叠度p下的资源消耗:

  1. 对于重叠度p ≥ 0.003,采样法资源更少,常规QPE多次测量即可满足目标置信区间。
  2. 对于极小的重叠度p < 0.003,振幅估计可带来平方根级别的加速(复杂度从1/p²降至1/√p),适用于极端稀疏重叠的初态。
  3. 本文还提供了误差项关于p、ε、置信度q的精确与渐进复杂度公式。

六、Fe-S簇系统模拟与资源估计

为验证该策略的实际应用潜力,研究团队将其应用于实际分子体系,包括[2Fe-2S]、[4Fe-4S]与FeMoco等难处理体系:

  1. 通过计算多个MPS初态与基态的重叠度,并引入一种基于双MPS态的外推方法,估计任意bond dimension对应的最终重叠度。
  2. 对FeMoco系统,在bond dimension为4000时能达到重叠度平方约0.9,这使得仅需2~3次QPE重复迭代即可获得足够置信的能量估计。
  3. 考虑所有步骤,总Toffoli门数量约为7.3×10¹⁰,仅为前人方法的2.3倍,说明此方案在实际容错量子计算中具有可行性。

七、结论与前景展望

本文提供了一套面向真实化学系统的高效初态制备与基态能量提取方法,克服了传统QPE中初始态完美重叠假设的局限,主要贡献如下:

  1. 提出新型高效的MPS合成方法,Toffoli复杂度大幅降低。
  2. 构建了两种基态提取过滤机制,在不同初态质量下具有良好适应性。
  3. 提供真实分子系统的资源估计与策略验证,为未来容错量子计算在化学模拟中的落地提供了可行路径。

发表评论

您的邮箱地址不会被公开。 必填项已用 * 标注