Point-group symmetry analysis of many-electron wavefunctions on a quantum computer

本文提出并验证了一种无需辅助比特的混合量子方法，用于计算多电子波函数的点群对称性投影权重，并在苯和二茂铁等分子的数值模拟与真实硬件上验证了其准确性。

要点

想象一下，你正试图聆听一首复杂的交响乐，但音乐如此混乱，以至于你无法分辨哪些乐器在演奏，也听不出它们敲击的是哪些音符。在量子化学的世界里，这首“交响乐”就是多电子波函数——一种描述分子中原子周围电子如何“舞蹈”的数学描述。

问题在于，这些电子遵循严格的对称性规则，就像编排好的舞蹈中的舞者一样。如果你知道对称性规则（即“点群”），你就能预测分子的行为、反应方式以及其能级。然而，在当前的量子计算机上，很难检查你对这些电子的数字模拟是否真的遵循了这些对称性规则。

本文介绍了一种新的、实用的工具，无需额外昂贵的设备即可检查这种“舞蹈编排”。以下是他们所做工作的简要分解：

1. 问题：“幽灵”量子比特

传统上，为了检查量子态是否具有正确的对称性，科学家使用一种需要“辅助”量子比特的方法。可以将这些辅助量子比特想象为幽灵助手。你需要引入这些额外的帮手来执行检查，但它们会占用量子计算机上宝贵的空间，并引入更多噪声（误差）。这就像试图用一台需要第二台重型天平来平衡的秤，去测量一根羽毛的重量。

2. 解决方案：“镜像”技巧

作者提出了一种巧妙的、无需辅助量子比特的方法。他们不使用幽灵助手，而是采用一种“镜像”技术。

• 类比：想象你有一个旋转的陀螺（即电子态）。为了观察它是否完美对称地旋转，你并不需要第二个陀螺。相反，你旋转对陀螺的视角（应用数学旋转），然后测量它与原始状态有多相似。
• 工作原理：他们取量子态，利用从分子形状推导出的特定数学规则对其进行旋转，然后测量旋转后的版本与原始版本之间的重叠度。这能确切地告诉他们，该状态有多少部分属于特定的对称性“家族”（称为不可约表示）。

3. 试驾：苯和二茂铁

为了证明该方法有效，他们在两种分子上进行了模拟：

• 苯：一个碳原子环（像六边形的蜂巢）。
• 二茂铁：一个铁原子夹在两个碳环之间（像分子三明治）。

他们在两类“舞者”上测试了该方法：

• 独舞者（斯莱特行列式）：描述电子的简单、单一编排。
• 复杂剧团（关联波函数）：更混乱、更真实的描述，其中电子之间相互作用强烈。

结果：他们的方法成功识别了简单和复杂情况下的对称性“权重”（即该状态有多少部分属于哪个对称性家族）。他们发现，有时即使模拟在能量表现上看起来不错，其对称性“风味”可能是错误的，而他们的工具能立即捕捉到这一点。

4. 现实世界演示：“嘈杂”的量子计算机

最激动人心的是，他们不仅是在完美的计算机模拟中完成这项工作，还在IBM 的"ibm kawasaki"量子设备上运行了该实验。

• 挑战：真实的量子计算机充满噪声。这就像试图在摇滚音乐会上听清耳语。信号会变得混乱。
• 对策：他们使用了先进的“降噪”技术（称为误差缓解）。这就像使用一种高科技麦克风，过滤掉人群的嘈杂声，从而清晰地听到耳语。
• 成果：利用多达 32 个量子比特（这对当前技术而言数量众多），他们成功测量了苯基态及其第一激发态的对称性权重。即使存在噪声，他们的“降噪”方法也使他们能够以极高的精度（误差在几个百分点以内）复现正确结果。

为什么这很重要

本文并不声称能一夜之间治愈疾病或制造新材料。相反，它为在当今不完美的量子计算机上工作的科学家提供了一套实用工具包。

这就像给机械师提供了一台新的、简单的诊断扫描仪，即使是在一辆老旧生锈的汽车上也能工作。在此之前，检查复杂量子模拟的对称性既困难又需要额外资源。现在，科学家们可以：

1. 检查工作：查看他们的量子模拟是否真正遵守了对称性定律。
2. 基准测试设备：利用此方法测试量子计算机的性能表现，以及其纠错工具的有效性。

简而言之，他们建立了一种方法，即使量子计算机有些嘈杂，也能清晰地聆听电子的“对称性音乐”，而无需任何额外的“幽灵”助手。

技术摘要

技术摘要：量子计算机上多电子波函数的点群对称性分析

问题陈述
点群代表了分子系统中的空间对称操作，是分析分子轨道和光谱学的化学基础工具。尽管已提出利用这些对称性的量子算法（例如用于减少量子比特数量或简化电路拟设），但点群对称性操作的实际实现以及对真实多电子波函数的详细对称性分析在很大程度上仍未被充分探索。具体而言，缺乏在近期量子硬件上具体实现对称性操作并分析真实关联多电子态对称性性质（例如不可约表示权重）的方法。

方法论
作者提出了一种无需辅助量子比特（ancilla-free）的量子 - 经典混合方法，用于分析多电子态的点群对称性，该方法适用于阿贝尔群和非阿贝尔群。该方法的核心在于计算给定波函数 $|\Psi\rangle$ 中不可约表示（$\Gamma$）的“权重”（$w_\Gamma$）。该权重定义为 $|\Psi\rangle$ 投影到 $\Gamma$ 子空间上的系数的平方。

理论框架依赖于投影算符公式：
$$w_\Gamma = \frac{d_\Gamma}{|G|} \sum_C \chi^*_\Gamma(C) \sum_{g \in C} \langle\Psi|\hat{g}|\Psi\rangle$$
其中 $d_\Gamma$ 是表示的维数，$|G|$ 是群的阶，$\chi_\Gamma(C)$ 是共轭类 $C$ 的特征标，$\hat{g}$ 是对称操作。

为了在量子计算机上无需辅助量子比特地评估 $\langle\Psi|\hat{g}|\Psi\rangle$，作者利用了表示矩阵的对角化。对于对称性适配轨道（例如 Hartree-Fock 轨道），对称操作 $U(g)$ 被分解为轨道旋转。通过对简并轨道集的表示矩阵 $\bar{D}_\gamma(g)$ 进行对角化，幺正算符 $U(g)$ 被重写为基旋转 $\tilde{U}(g)$ 与对角相位算符的乘积。随后，通过制备状态 $|\tilde{\Psi}(g)\rangle = \tilde{U}(g)|\Psi\rangle$，在计算基下测量，并将结果与适当的相位因子结合，从而获得期望值。这种方法避免了哈达玛测试（Hadamard tests）或线性组合幺正算符（LCU）技术所需的深层电路和辅助量子比特要求。

本研究采用了两种主要的验证策略：

1. 经典数值模拟：利用 ffsim 和 pyscf 库，将该方法应用于苯（$D_{6h}$）和交错二茂铁（$D_{5d}$）。作者分析了单 Slater 行列式（Hartree-Fock 基态和单激发态）以及关联波函数（Unitary Cluster Jastrow 和 Local Unitary Cluster Jastrow 函数）。
2. 硬件演示：该方法在 IBM 的 ibm_kawasaki 超导量子设备（最多 32 个量子比特）上执行。目标系统为具有 $D_{2h}$ 对称性（$D_{6h}$ 的子群）的苯分子。试验波函数（基态和第一激发态）通过密度矩阵重整化群（DMRG）计算制备，经张量网络技术压缩为砖墙拟设（brick-wall ansatz），并映射到量子电路。为应对噪声，比较了两种误差缓解技术：通过门折叠实现的零噪声外推（ZNE）和准概率误差缓解方法（QESEM）。

关键结果

• 经典模拟：
  • 对于苯和二茂铁的 Hartree-Fock 基态，该方法正确识别出全对称不可约表示（$A_{1g}$）的权重接近 1.0，这与闭壳层构型一致。
  • 对于单激发（SE）态，该方法成功将波函数分解为多个不可约表示。例如，苯中的一个特定 SE 态被分解为 $\approx 0.005 B_{1u} + 0.495 B_{2u} + 0.5 E_{1u}$。对所有独立 SE 构型的权重求和得出了与群论约化公式一致的不变结果。
  • 在关联 LUCJ 波函数的分析中，该方法揭示出最小化能量误差（无正则化）的优化程序可能会显著降低对称性权重（偏离正确的 $A_{1g}$ 权重），而正则化有助于保持对称性性质，尽管会带来不同的能量权衡。这表明该指标可用于评估试验波函数的质量。
• 硬件演示：
  • 在 ibm_kawasaki 设备上，原始测量结果显示，随着量子比特数量的增加，对称性权重显著下降（例如，对于 32 个量子比特，基态的 $A_g$ 权重在无缓解情况下降至约 0.54）。
  • 应用误差缓解至关重要。虽然 ZNE 改善了结果，但准概率 QESEM 方法产生了最高的保真度，在百分之几的误差范围内复现了正确的不可约表示权重（例如，32 个量子比特上基态的结果为 $0.968 \pm 0.010$，而精确值为 $0.999$）。
  • 实验证实，所提出的无需辅助量子比特的方法，结合张量网络态制备和先进的误差缓解技术，能够在当前硬件上忠实复现对称性性质。

意义与主张
本文主张填补理论对称性投影与量子硬件实际实现之间的空白。所提出的方法作为实用工具，可用于：

1. 对称性分析：在近期和早期容错量子计算机上，分析真实材料模拟中多电子波函数的对称性性质。
2. 波函数评估：通过量化试验波函数（例如在变分算法中）对特定对称性的遵循程度来评估其质量。
3. 基准测试：为真实量子设备和误差缓解技术提供基准，因为该方法无需辅助量子比特，且对于对称性适配轨道的开销极小。

作者强调，虽然硬件演示利用了阿贝尔子群（$D_{2h}$）和基于泡利算符的评估，但底层方法是通用的，适用于非阿贝尔群和任意基函数。这项工作表明，对称性分析是未来复杂材料量子模拟中可行且必要的组成部分。