Determining the ensemble N-representability of Reduced Density Matrices

技术摘要：确定约化密度矩阵的系综 N-可表示性

问题陈述
N-可表示性问题是电子结构理论中的一个基础性挑战，涉及一个 $p$-体约化密度矩阵（$p$-RDM）必须满足哪些必要且充分的条件，才能由 $N$-电子波函数或混合（系综）态导出。尽管在纯态 N-可表示性方面已取得了显著进展，但确定 $p$-RDM 的系综 N-可表示性仍然是一个关键障碍，特别是在涉及热混合、开放量子系统和统计系综的应用中。验证给定的 $p$-RDM 是否对应于一个物理的 $N$-电子态在一般情况下在计算上是不可行的，属于量子 Merlin-Arthur 完全复杂度类（quantum Merlin-Arthur-complete）。现有的近似条件对于强关联系统往往不足，且对于 $p > 1$ 的 $p$-RDM，其完整的充分必要条件仍未可知。

方法论
本研究提出了一种通过将先前开发的幺正演化算法（纯态 ADAPT-VQA）扩展到系综情况，来确定 $p$-体矩阵系综 N-可表示性的实用框架。核心方法论依赖于一种纯化策略：

1. 纯化嵌入（Purification Embedding）： 根据 Schrödinger-Hughston-Jozsa-Wootters 定理，将定义在希尔伯特空间 $\mathcal{H}_s$ 上的系综态 $\rho$ 嵌入到定义在扩展希尔伯特空间 $\mathcal{H}_{sb} = \mathcal{H}_s \otimes \mathcal{H}_b$ 上的纯态 $|\Psi_{sb}\rangle$ 中，其中 $\mathcal{H}_b$ 是辅助空间。纯化过程的构建确保了纯化态的约化密度矩阵能够重现原始系综的约化密度矩阵，即 $\rho = \text{Tr}_b(|\Psi_{sb}\rangle\langle\Psi_{sb}|)$。
2. 变分幺正演化（Variational Unitary Evolution）： 该算法采用一种自适应导数组装伪 Trotter（ADAPT）变分量子算法（VQA）。它通过迭代地对初始纯化态（由初始 $N$-粒子态 $|\phi_0\rangle$ 及其纯化构造而成）应用变分幺正算符来进行演化。
3. 算符池（Operator Pool）： 演化由适用于费米子的反埃尔米特激发-去激发算符池驱动，并受到正则系综约束（系统或辅助空间内粒子数不变）。这些算符通过 Jordan-Wigner 变换映射到 Pauli 算符，使得该算法适用于量子计算机模拟。
4. 优化目标： 算法最小化演化后的纯化态的 $p$-RDM 与目标 $p$-RDM（$p\rho_{\text{target}}$）之间的希尔伯特-施密特距离（$D$）。
   • 如果算法收敛到距离 $D_{\min} \approx 0$，则认为目标是系综 N-可表示的。
   • 如果 $D_{\min} > 0$，则该距离的大小作为 N-可表示性违背程度的定量度量。
   • 演化后的状态提供了与缺陷目标最接近的、物理上可实现的系综可表示约化态。

主要贡献

• 系综 ADAPT-VQA 算法： 引入了一种通过纯化将纯态可表示性算法扩展到系综态的混合框架。这使得算法能够区分目标矩阵是源自纯量子态还是系综混合态。
• 误差修正与精炼： 该方法提供了一条密度矩阵精炼的路径。它可以修正有缺陷的系综约化密度矩阵（即违反 N-可表示性条件的矩阵），通过将其向最近的物理态进行演化来实现。
• 量子态重建： 该算法促进了量子计算机上的量子态重建，能够从目标 RDM 中恢复兼容的波函数或纯化态。
• 数值验证： 该框架通过在模型系统和分子系统（有限温度下）上的数值模拟进行了验证，展示了其处理 1-RDM 和 2-RDM 的能力。

结果
作者使用无噪声量子模拟器对算法进行了验证：

• 模型系统（零温）： 在具有 3 个和 4 个轨道的 4-电子系统上进行了测试。
  • 对于 1-RDM，算法成功区分了纯态和系综 N-可表示态。在 Klyachko 的纯态 N-可表示性条件（广义泡利条件）被违反（混合态）的情况下，纯态 ADAPT-VQA 无法收敛到零距离，而系综 ADAPT-VQA 则收敛至 $D_{\min} \approx 0$。
  • 对于 2-RDM（其完全纯态 N-可表示性条件尚不明确），算法正确识别了所有构建目标的系综 N-可表示性。它进一步展示了检测系综 2-RDM 是否能收缩为纯 1-RDM 的能力。
  • 噪声鲁棒性： 通过在目标矩阵中引入随机噪声，研究表明 $D_{\min}$ 随噪声强度（$\epsilon$）的变化而缩放，证实了该算法在量化 N-可表示性违背严重程度方面的效用。
• 分子系统（有限温度）： 算法被应用于有限温度下的 $H_2$ 和线性 $H_3$ 分子。目标矩阵由正则系综热态构造并添加了噪声。
  • 结果表明，收敛距离 $D_{\min}$ 对添加的噪声（$\epsilon$）比对底层无噪声态的内在相关特性更敏感。
  • 随着噪声增加，$D_{\min}$ 显著增加（例如，对于 $\epsilon=0.1$ 的 2-RDM，$D_{\min}$ 达到 $\sim 10^1$），证明了该算法在现实分子场景下检测和量化可表示性缺陷的能力。

意义
本文声称，这项工作为解决 N-可表示性问题提供了一个实用且稳健的框架，而这一问题在历史上在计算上是难以处理的。通过利用纯化和变分量子算法，所提出的方法提供了双重能力：它既可以对给定的 RDM 进行分类（纯态 vs 系综），也可以修正有缺陷的矩阵以强制执行物理约束。这对于提高量子化学计算的准确性具有重要意义，特别是对于强关联电子系统和有限温度应用。作者强调，该算法与量子硬件兼容，为未来的量子模拟中的密度矩阵精炼和量子态重建提供了路径。