✨ 要点🔬 技术摘要
在材料科学领域,理解电子如何运动和相互作用是开启新技术(从更好的电池到更快的计算机)的关键。几十年来,科学家们一直依赖强大的经典超级计算机来模拟这些微小粒子,但随着他们研究的系统变得更加复杂,计算变得异常困难。电子并非孤立运动;它们是深度纠缠的,这意味着一个电子的状态会瞬间影响其他电子,形成一个难以追踪的关系网络。为了解决这个问题,研究人员转向了量子计算机,这种机器旨在模仿自然界的自身规则。人们曾希望这些设备能够找到分子的最低能量状态,这个数值可以告诉我们一种材料的稳定性。然而,一个令人惊讶的转折出现了:对于仅仅寻找这一个静态能量值而言,巧妙的新型经典算法已经追平甚至在某些情况下超越了当前一代量子机器。现在的疑问不再是量子计算机能否计算一个静态数字,而是它们能否做一些经典机器难以胜任的事情:追踪材料如何随时间变化以及如何对能量做出响应,从而揭示其电子的动态行为。
一位研究人员通过将目标从静态能量转向动态运动,在这一方向上迈出了重要一步。研究人员不再要求量子计算机寻找分子的静止能量,而是要求它模拟电子的实时演化,并测量由此产生的模式。他们开发了一种方法,利用量子处理器生成一系列被称为“比特串”(bitstrings)的简单数据点,这些数据点代表了系统构型的快照。这些快照是从一个非常浅层、简单的电路中收集的,该电路不需要早期量子实验中困扰已久的复杂且易错的操作。一旦收集到这些快照,经典计算机就会接管工作,将它们组装成一张详细的材料谱函数图(spectral function)。这张图本质上是材料的指纹,展示了它在不同频率和动量下吸收和释放能量的情况。研究人员在一条原子链和十九种不同的分子上测试了这种方法,在所有案例中,重建的图谱都以高精度匹配了精确的理论结果。他们甚至在真实的量子处理器 IBM Heron 上运行了实验,成功重现了小规模原子链的谱函数,证明了该方法可以在实际硬件上运行,而无需依赖最先进的纠错技术。
这项研究还探讨了一个困扰该领域的根本问题:究竟是什么让一个问题对量子计算机来说变得困难?多年来,科学家们一直假设电子的一种特定属性——通常被称为“魔性”(magic)或非高斯性(non-Gaussianity)——是导致困难的主要驱动力。这种属性被认为衡量了一个系统偏离简单、无相互作用粒子的程度。研究人员在他们的分子系列中对这一假设进行了严格测试,发现这一假设是不正确的。他们发现,虽然这个“魔性”数值对于识别电子关系的复杂程度很有用,但它与运行模拟的实际成本之间并没有可靠的联系。一个系统可能具有很高的“魔性”,但却很容易模拟;或者“魔性”很低,却极其难以模拟。相反,真正的成本是由另一个因素决定的:即通过量子计算机特定的测量视角所看到的系统的纠缠结构。他们发现,计算机为了获得准确答案而必须采样的不同配置的数量,直接取决于“键维”(bond dimension)——这是一个衡量在机器使用的特定参考框架下,电子之间有多少种连接的度量。这一发现阐明了,困难并不在于一个单一的、易于测量的数字,而在于那个经典计算机难以压缩的复杂、高维的连接网络。
为了确保方法的稳健性,研究人员还探索了人工智能在过程中的作用。他们训练了一个机器学习模型来尝试预测数据的缺失部分,希望它能填补量子计算机有限采样留下的空白。结果是发人深省且富有启发性的:该学习模型的表现并未优于一种更简单的经典统计方法。事实上,经典方法的效果与之不相上下。这表明,量子计算机的角色不是被智能算法取代,而是提供没有任何经典方法能够独立生成的原始、高质量数据。量子机器充当了一个独特的传感器,捕捉模拟的原材料,而经典计算机和简单的统计工具则负责完成将这些材料组装成连贯故事的繁重工作。研究人员还发现,该方法对当前量子硬件中固有的噪声和误差具有惊人的韧性。通过使用一种自洽的恢复过程,他们可以过滤掉受损的数据点,即使在硬件出错的情况下,仍能得到准确的结果。这种噪声韧性是一个至关重要的特性,因为它意味着该方法可以在无需等待完美的、容错的未来量子计算机之前,就在当今不完美的机器上运行。
这项工作的意义不仅在于一种新算法,它还重新定义了我们对近期量子计算机的预期。计算动态谱函数的能力,为研究目前难以触及的现象打开了大门,例如在光或磁场影响下电子在复杂材料中的运动方式。这些正是实验物理学家在实验室中使用光电子能谱等工具进行的测量,能够模拟这些测量将彻底改变新材料的设计。该研究表明,前进的路径不是去追求基态能量(在这一领域经典方法已经占据主导地位),而是专注于这些动态的、频率分辨的观测值,因为在这些领域,量子优势最为明确。通过识别驱动成本的真正资源——即基于基底的纠缠而非轨道不变的“魔性”——研究人员为量子计算机最有用的方向提供了清晰的地图。他们表明,这些机器的真正力量在于它们能够采样那些经典计算机难以生成的复杂分布,并且只要提出正确的问题并使用正确的工具来解读答案,这种力量在今天就可以被用来解决重要的实际问题。
技术摘要:基于位串采样量子子空间的动力学谱函数
问题陈述 基于采样的量子对角化(SQD)和量子选择构型相互作用(QSCI)已成为电子结构领域领先的近即时量子算法,主要针对基态能量。然而,针对基态能量的经典前沿已经关闭;启发式经典方法(例如 matchgate 模拟器、自回归神经采样器、热浴构型相互作用)目前的性能已达到或超过了量子采样子空间在处理该特定可观测量时的表现。因此,关于基态化学具有通用指数级量子优势的说法正变得日益站不住脚。
本文识别了一个关键空白:目前尚无已发表的方法能从位串采样的量子子空间中重建频率解析的单体谱函数 A ( ω ) A(\omega) A ( ω ) 或动力学结构因子 S ( q , ω ) S(q, \omega) S ( q , ω ) ,尤其是在格林函数所需的带电(N ± 1 N \pm 1 N ± 1 )扇区中。此外,关于究竟是哪种资源决定了采样成本,目前仍缺乏明确解释。一个自然的假设是“一阶费米子魔力”(one-body fermionic magic,一种非高斯性的度量)驱动了计算难度,但这仍未得到证实,且可能具有误导性。
方法论 作者提出了一个统一的原语,旨在仅利用计算基测量,将 SQD/QSCI 工具箱从静态能量提升到连续的、频率解析的动力学领域。
采样原语: 该方法使用浅层实时演化电路,无需 Hadamard 测试、受控幺正算符或重叠估计。
种子(Seeds): 对于单体谱,系统在 ( N ± 1 ) (N \pm 1) ( N ± 1 ) 扇区中准备响应种子 ∣ ϕ ± ⟩ = c ^ p † ∣ 0 ⟩ |\phi^\pm\rangle = \hat{c}^\dagger_p |0\rangle ∣ ϕ ± ⟩ = c ^ p † ∣0 ⟩ 和 c ^ p ∣ 0 ⟩ \hat{c}_p |0\rangle c ^ p ∣0 ⟩ 。对于中性动力学结构因子,种子为守恒粒子数的密度或自旋算符 n ^ q ∣ 0 ⟩ \hat{n}_q |0\rangle n ^ q ∣0 ⟩ 或 S ^ q z ∣ 0 ⟩ \hat{S}^z_q |0\rangle S ^ q z ∣0 ⟩ 。
演化: 这些种子在时间演化(e − i H ^ t k e^{-i\hat{H}t_k} e − i H ^ t k )下进行,并通过计算基进行测量以生成位串。
子空间构建: 跨越不同演化时间的采样构型集合构成了一个行列式子空间 S \mathcal{S} S 。
经典重建: 量子处理器仅提供位串。哈密顿量被投影到 S \mathcal{S} S 上并在经典端进行对角化。利用莱曼表示(Lehmann representation)组装得到留退格林函数(retarded Green's function)和谱函数:G p q ( ω ) = ∑ m ⟨ 0 ∣ c ^ p ∣ m ~ + ⟩ ⟨ m ~ + ∣ c ^ q † ∣ 0 ⟩ ω − ( E ~ m + − E 0 ) + i η + … G_{pq}(\omega) = \sum_m \frac{\langle 0|\hat{c}_p|\tilde{m}^+\rangle \langle \tilde{m}^+|\hat{c}^\dagger_q|0\rangle}{\omega - (\tilde{E}^+_m - E_0) + i\eta} + \dots G pq ( ω ) = m ∑ ω − ( E ~ m + − E 0 ) + i η ⟨ 0∣ c ^ p ∣ m ~ + ⟩ ⟨ m ~ + ∣ c ^ q † ∣0 ⟩ + … 由此产生 A ( ω ) A(\omega) A ( ω ) 、动量解析的 动量解析的 动量解析的 A(k, \omega)以及动力学结构因子 以及动力学结构因子 以及动力学结构因子 S(q, \omega)和 和 和 S^{zz}(q, \omega)$。
噪声缓解: 该方法采用自洽构型恢复(S-CoRe)技术,将带有噪声的位串投影回正确的粒子数扇区,从而提高在现实设备噪声下的子空间质量。
核心贡献
动力学谱重建: 本文展示了首次从位串采样子空间中重建频率解析的单体谱函数和动力学结构因子。这使 SQD/QSCI 范式从离散能级扩展到了连续的动力学可观测量,并得到了 Hubbard 链和 19 个分子集的精确对角化验证。
魔力与成本的解耦: 作者严格证明了一阶费米子魔力 (通过费米子 AntiFlatness F 1 = 4 tr [ γ ( 1 − γ ) ] F_1 = 4 \text{tr}[\gamma(1-\gamma)] F 1 = 4 tr [ γ ( 1 − γ )] 量化)与采样成本是解耦 的。
F 1 F_1 F 1 是轨道旋转不变的(高斯不变性),而采样成本(行列式支撑 ∣ S ∣ |S| ∣ S ∣ )是依赖于基底的。
虽然 F 1 F_1 F 1 能忠实诊断静态多参考相关性(例如在键断裂期间上升),但它无法预测采样成本。例如,增加 Hubbard 链中的在位排斥力会提高 F 1 F_1 F 1 ,但由于局域化效应会降低 ∣ S ∣ |S| ∣ S ∣ 。
真正的成本驱动因素被确定为依赖于基底的纠缠 ,具体表现为最小键维数 χ \chi χ (张量网络键维数),它在经验上与 ∣ S ∣ |S| ∣ S ∣ 具有高度相关性(ρ ≈ 0.90 \rho \approx 0.90 ρ ≈ 0.90 )。
资源诊断: 本研究确立了真正的量子优势(如果存在的话)在于固定采样基底中连接的高阶多体累积量的秩(rank)和非高斯性,而非轨道不变的一阶魔力。
结果
验证: 该方法重构了 Hubbard 链和 19 个分子的精确谱函数(经全构型相互作用验证,误差 < 10 − 5 <10^{-5} < 1 0 − 5 Ha),其相对 L 1 L_1 L 1 误差低于 10 − 3 10^{-3} 1 0 − 3 。
硬件演示: 该算法在 IBM Heron 处理器(用于 L = 6 L=6 L = 6 Hubbard 模型的 i b m _ f e z ibm\_fez ibm _ f ez ;用于拉伸 N 2 N_2 N 2 的 i b m _ m a r r a k e s h ibm\_marrakesh ibm _ ma r r ak es h )上执行。
L = 6 L=6 L = 6 Hubbard 谱函数通过 50,000 次采样重建,在展宽分辨率内与精确结果相匹配。
对于拉伸的 N 2 N_2 N 2 ,噪声辅助构型恢复实现了 0.59 ± 0.14 0.59 \pm 0.14 0.59 ± 0.14 mHa 的基态能量误差(低于化学精度),优于相同浅层电路的无噪声模拟。
缩放与比较:
在固定精度下,采样成本(所需扇区的比例)随系统规模呈指数级下降,尽管绝对行列式支撑 ∣ S ∣ |S| ∣ S ∣ 随之指数增长。
在强耦合机制(U / t ≥ 12 U/t \ge 12 U / t ≥ 12 )下,一旦子空间超过约 120 个行列式,时间演化选择器在性能上优于确定性的经典 Krylov 和 CIPSI 选择器,但在可及规模下仍未超越热浴 CI。
AI/ML 限制: 作为“尾部补全器”来补充设备采样的学习型生成模型(GFlowNet)在统计上并未优于廉价的经典 CIPSI 基准,这表明当前的生成模型尚不能取代此类背景下的量子采样需求。
意义与主张
本文重新定义了电子结构领域的量子优势辩论:
目标转移: 作者认为,有力的量子优势不在于基态能量(因为经典方法已能与之抗衡),而在于动力学可观测量 (谱函数、结构因子),这些量对于材料发现和关联电子物理至 de 核心,但在维度 > 1 >1 > 1 时对经典张量网络而言依然困难。
资源澄清: 它纠正了“一阶魔力证明采样难度”的误解,并提供了一个精确的资源说明:一阶魔力是相关类型 (多参考特性)的诊断工具,而成本 则由依赖于基底的纠缠(键维数 χ \chi χ )决定,一阶魔力无法预测这一点。
实际效用: 该方法提供了一种硬件原生的方法,无需辅助比特或复杂的逻辑门,即可计算当前仅能处理能量的硬件无法提供的可观测量。它作为一个“分诊”工具:高效可测的费米子魔力可以标记出均场方法失效的强关联目标,但决定是否使用量子处理器的依据必须是对依赖于基底的成本(键维数)的估计,而一阶魔力无法预测此项。
作者总结道,虽然该方法并未声称在所测试的具体实例中具有“超越经典”的优势(因为对于一维和窄梯子,经典 DMRG 仍然有效),但它成功定义了量子采样变得必要的边界:即在固定采样基底的高秩、非高斯扇区中,那里超出了自由费米子估计和低键维数张量网络的覆盖范围。
每周获取最佳 quantum physics 论文。
受到斯坦福、剑桥和法国科学院研究人员的信赖。
请查收邮箱确认订阅。
出了点问题,再试一次?
无垃圾邮件,随时退订。