在量子物理的微观世界中，被称为费米子（如电子）的粒子可以以两种截然不同的方式进行排列。有时，它们的行为就像一种平静、可预测的流体，每一部分都以一种简单的线性方式与其他部分完美相关。物理学家称这种状态为“高斯”（Gaussian）态，由于其规则非常直观，即使是最强大的经典计算机也能轻松模拟这些系统。但当这些粒子以复杂且混沌的方式相互作用时，它们便进入了一个不同的领域。在这里，相关性变得错综复杂且具有层次感，难以进行简单描述。这种“非高斯”（non-Gaussian）行为是让量子计算机潜在地比经典计算机更强大的秘密成分，但它也是最难测量的东西。多年来，科学家们一直致力于量化大型系统中的这种复杂性，却经常撞上一堵墙：计算成本增长得过快，以至于变得不再实用，这实际上将我们的视野限制在了极小规模的粒子集合上。
东京大学的一个研究小组现在开发出了一种新技术，能够窥探过这堵墙。他们创造了一种技术，可以计算一种被称为“魔术雷尼熵”（magic Rényi entropy）的特定量子复杂性度量，该方法可用于包含多达 128 个格点的系统。这是一个显著的飞跃，因为以往的方法只能处理足够小到能放在标准台式电脑上的系统。研究人员通过设计一种巧妙的采样策略实现了这一目标，该策略直接作用于一种被称为“矩阵乘积态”（matrix product state）的量子态数学表示。他们并没有试图一次性计算整个系统——这项任务需要无法想象的内存——而是生成了系统行为的单个快照。通过仔细控制这些快照中的统计涨落，他们能够以高精度重建系统的整体复杂性图景。
该团队在一种由相互作用电子组成的链状模型上测试了其方法，这种系统被称为 XXZ 链。他们将模拟推向了包含多达 128 个格点的链，这一规模在以前是无法通过此类特定计算实现的。在系统行为的一个极端——即粒子以特定的磁有序排列时——他们的结果与已知的理论预测完美吻合，证实了其方法是正确的。更重要的是，在粒子相互作用更强的其他区域，他们的方法揭示了一种新现象：他们观察到，随着链条变长，量子复杂性的度量停止增长并趋于平缓，即达到饱和。这种饱和发生的规模对于旧有的精确计算方法来说已经太大了，无法触及。
这一发现至关重要，因为它表明系统的复杂性不仅仅是其各部分的简单加总。研究人员发现，这种趋于平缓的过程与粒子在链条上相关性的特定模式形成相吻合，这种模式看起来像是一道分隔两个不同磁性区域的墙。传统的工具如果只关注粒子之间简单的两点连接，将会完全忽略这一点，因为它们无法观测到这些高阶的多粒子关系。然而，这种新方法捕捉到了这些复杂的连接，证明了系统的量子复杂性与这些大规模结构有着深刻的联系。
这项工作的意义不仅限于这一个模型。该方法为科学家研究量子资源如何在大型系统中演化（无论是处于稳定的基态还是随时间变化）提供了一个实用的工具。由于他们计算的这种度量也可以通过使用多个量子态副本在现实世界的实验中进行观测，因此这些数值结果可以直接与来自未来量子模拟器的实验数据进行对比。这架起了理论预测与实验现实之间的桥梁，为基准测试和理解构建通用量子计算机所必需的非高斯资源提供了方法。通过使量化这些复杂相关性的规模达到现实水平成为可能，研究人员为观察多体费米子系统的行为开启了一扇新的窗口，将一个曾经难以处理的问题转化为了一个可控的计算过程。
技术摘要：多体费米子非高斯性的高效采样
问题陈述
费米子非高斯性是通用量子计算的关键资源，它表征了相互作用电子系统中超越平均场或高斯描述的高阶相关性。虽然基于协方差矩阵的度量提供了计算效率高的诊断方法，但它们无法捕捉到在强相关机制中至关重要的高阶相关性。一种更严谨的资源度量——基于卷积操作的魔术 Rényi 熵（Magic Rényi Entropy, MRE）已被提出。然而，由于高阶相关性的数量随系统规模呈指数级增长，且通过矩阵乘积态（MPS）进行的直接计算涉及极高的成本（例如，收缩维度为 O(D8) 的转移矩阵），因此 MRE 的精确评估仅限于较小的系统规模。
方法论
作者开发了一种完美采样法，用于从矩阵乘积态（MPS）中计算二阶 MRE (M2)，而无需显式构建卷积态或复制张量网络。该方法的核心包含两项关键创新：
通过 Rao-Blackwellization 实现的有界估计量：
纯度 ν2（从中导出 M2=−lnν2）被重新表述为输出态在 Fock 基底上的期望值。作者定义了一个与交换输出构型相关的振幅比 r(X)。由于无偏估计量 r(X) 是无界的且容易产生剧烈波动，作者应用了 Rao-Blackwellization 技术。通过在具有条件权重的交换对上对 r(X) 进行平均，他们构造了一个有界的估计量 h(X)，满足 −1≤h(X)≤1。这消除了由微小分母引起的巨大波动，同时保持了均值并降低了方差。
来自 MPS 的递归完美采样：
为了从 Born 分布 ∣Φ∣2 中进行采样，并在不构建完整输出 MPS 的情况下评估必要的振幅，作者采用了一种递归收缩方案：
- 他们将输入态表示为一个归一化的右正则（right-canonical）MPS。
- 他们引入了两个解析了宇称（parity-resolved）的边界矩阵，用以编码采样的前缀。
- 通过使用扫掠算法（Algorithm 1），他们通过收缩局部 MPS 张量、分束器幺正算符以及边界矩阵，来计算下一位点输出构型的条件概率。
- 该过程生成了精确的 Born 样本，并能以 O(nD3) 的计算成本（其中 n 是模式数，D 是 MPS 键维数）评估给定的振幅（包括估计量所需的“交叉”振幅）。
- 该方法可以通过利用粒子数守恒（U(1) 对称性）进一步优化，将计算限制在允许的电荷扇区内。
关键结果
该方法在包含多达 128 个位点的自旋化费米子 XXZ 链上进行了基准测试：
- 铁磁点 (Δ=−1)： 采样得到的 M2 在所有直至 L=128 的系统规模下，均与解析推导出的 Dicke 态渐近公式完美一致。
- 超越铁磁点 (Δ<−1)： 该方法揭示了 MRE 在精确对角化无法触及的系统规模（L>32）下的饱和现象。这种饱和现象与自旋-自旋相关函数中出现的畴壁（domain-wall）轮廓相吻合。
- 与协方差度量的比较： 研究表明，虽然基于协方差矩阵的度量（如费米子反平坦度/antiflatness）对于整个族系的地态都给出相同的数值（由于自旋翻转对称性和粒子数守恒），但 MRE 能够成功区分这些状态，并捕捉到向畴壁形成的转变。
- 计算效率： 采样成本随 O(NnD3) 缩放，其中 N 是样本数。对于广延性 MRE，所需样本数随 eM2 缩放，虽然这相对于资源值呈指数级增长，但在处理高达 ∼100 个位点的系统时仍然可行。对于 L=128 且 Δ=−1 的情况，约 1.4×107 个样本即可达到 0.01 的标准误差。
意义与主张
本文声称提供了一种实用的计算工具，通过考虑高阶相关性，定量地评估多体费米子的非高斯性。这项工作的意义在于：
- 可扩展性： 使研究高达 128 个位点的系统中的非高斯性成为可能，这在以前是 MRE 评估无法企及的规模。
- 物理洞察： 证明了 MRE 可以解析复杂的许多体相关性（例如畴壁形成），而这些相关性在基于高斯的诊断方法中是不可见的。
- 实验相关性： 由于 MRE 是通过纯度和通过复制件的干涉（如 SWAP 测试）进行定义的，其数值结果为费米子量子模拟器提供了直接的基准。
- 广泛适用性： 该方法适用于由 MPS 表示的平衡态和非平衡态，允许追踪量子猝灭（quantum quench）后的非高斯资源演化。
作者强调，他们的方法避免了显式构建分束器输出 MPS，从而使采样成本保持在关于系统规模和键维数的低阶多项式水平，使其成为大规模费米子系统的一种极具竞争力的替代方案。
每周获取最佳 quantum physics 论文。
受到斯坦福、剑桥和法国科学院研究人员的信赖。
请查收邮箱确认订阅。
出了点问题，再试一次？
无垃圾邮件，随时退订。