来源论文: https://arxiv.org/abs/2604.02307v1 生成时间: Apr 03, 2026 11:56

0. 执行摘要

纠缠熵（Entanglement Entropy, EE）作为探测量子多体系统非平凡特性的强大工具，其重要性已在凝聚态物理与量子信息领域得到广泛公认。然而，传统的纠缠熵度量往往忽略了物理系统普遍存在的全局对称性。**对称性分辨纠缠（Symmetry-Resolved Entanglement, SRE）**的提出，旨在将总纠缠分解为来自不同对称性部门（Symmetry Sectors）的贡献，从而提供更精细的量子态结构探测。

尽管在(1+1)D共形场论（CFT）框架下，SRE 的理论研究已较为成熟，但在强相互作用的高维量子体系中，如何精确计算 SRE 仍是一个巨大的挑战。传统的精确对角化（ED）受限于系统尺寸，而张量网络（Tensor Networks）在处理高维体系时往往面临收敛与算力瓶颈。

本工作（来自复旦大学、香港中文大学等机构的研究团队）提出了一种基于量子蒙特卡洛（QMC）的通用框架。该方法通过在副本流形（Replica Manifolds）上测量“无序算子”（Disorder Operators），即对称性扭曲算子，成功实现了对强相互作用 1D 和 2D 体系对称性分辨瑞尼纠缠熵（SRRE）的精确提取。研究重点验证了纠缠均分（Entanglement Equipartition）原理在 (2+1)D 关键点的适用性，为研究超越一维的强关联量子物质提供了高效的数值工具。

1. 核心科学问题，理论基础，技术难点与方法细节

1.1 核心科学问题：为何需要对称性分辨？

对称性是量子力学的基础。当一个系统具有全局对称性（如 $U(1)$ 或 $Z_2$）时，其总电荷 $Q$ 是守恒的。对于一个双分区系统（子系统 $A$ 及其补集 $\bar{A}$），其电荷算符可分解为 $Q = Q_A + Q_{\bar{A}}$。由于纠缠谱与电荷算符在基底下是对角化的，约化密度矩阵 $\rho_A$ 会呈现出分块对角结构：

$$\rho_A = \bigoplus_q \tilde{\rho}_A(q)$$

其中 $q$ 是子系统电荷 $Q_A$ 的特征值。对称性分辨瑞尼熵（SRRE） $S_n(q)$ 定义为在特定电荷部门 $q$ 内的瑞尼熵。科学界的核心疑问在于：这些不同部门的纠缠是如何分布的？是否存在所谓的“纠缠均分”，即在热力学极限下，$S_n(q)$ 与电荷 $q$ 无关？

1.2 理论基础：从带电矩到 SRRE

计算 SRRE 的传统难点在于无法直接获取分块后的密度矩阵。本文采用了**带电矩（Charged Moments）**框架。定义第 $n$ 阶瑞尼瑞纠缠熵的带电矩为：

$$Z_n(\alpha) = \text{Tr}(\rho_A^n e^{i\alpha Q_A})$$

通过对 $\alpha$ 进行傅里叶变换，可以提取出对应于电荷 $q$ 的分量：

$$Z_n(q) = \text{Tr}(\Pi_q \rho_A^n) = \int_{-\pi}^{\pi} \frac{d\alpha}{2\pi} e^{-i\alpha q} Z_n(\alpha)$$

最终，SRRE 可表示为：

$$S_n(q) = \frac{1}{1-n} \ln \frac{Z_n(q)}{[Z_1(q)]^n}$$

其中 $Z_1(q) = P(q)$ 是在子系统 $A$ 中测量到电荷 $q$ 的概率分布。

1.3 技术难点：高维体系的采样瓶颈

在 QMC 模拟中，直接测量 $\text{Tr}(\rho_A^n)$ 本身就极其困难（即著名的“比率挑战”）。而要进一步测量带有相位因子 $e^{i\alpha Q_A}$ 的算符，则会引入严重的相位不确定性。此外，对于 (2+1)D 体系，纠缠熵遵循面积律，其修正项（如对数项）极其微小，需要极高的数值精度才能分辨。

1.4 方法细节：副本流形上的无序算子测量

本文的创新之处在于将 $Z_n(\alpha)$ 转化为副本流形上的物理观测值测量。对于 $n=2$（二阶瑞尼熵）：

  1. 单副本系综（Single-replica）：测量 $\langle e^{i\alpha Q_A} \rangle_Z$，这等同于在标准配分函数下的无序算子期望值。
  2. 双副本流形（Two-replica manifold）：构造两个相互独立的系统副本，在子系统 $A$ 区域进行“胶合”（即沿着虚时方向进行交换连接）。在此流形上测量 $\langle e^{i\alpha Q_A} \rangle_{Z_A^{(2)}}$。

通过这种方式，SRRE 与总瑞尼熵之差（减去后的 SRRE）可以表示为：

$$S_2(q) - S_A^{(2)} = \ln \frac{\mathcal{F}[\langle e^{i\alpha Q_A} \rangle_{Z_A^{(2)}}]}{\mathcal{F}[\langle e^{i\alpha Q_A} \rangle_Z]^2}$$

其中 $\mathcal{F}$ 代表傅里叶变换。该方法的优势在于，我们不需要单独计算绝对的纠缠熵值，而是直接通过无序算子的比率来提取对称性分辨的信息，极大地抵消了系统误差。

2. 关键 Benchmark 体系、计算数据与性能表现

2.1 横场伊辛模型 (TFIM) 的验证

研究者首先在 $Z_2$ 对称性的 TFIM 上测试了该算法。

  • 1D TFIM (量子临界点 $h_c/J=1$)：

    • 预期：CFT 预测无序算子 $\langle X \rangle$ 具有对数标度：$-\ln\langle X \rangle \sim \frac{1}{4n} \ln L$。
    • 结果：QMC 数据显示单副本斜率为 0.251(1)，双副本斜率为 0.1242(4)，与理论值 1/4 和 1/8 完美契合。
    • SRRE 表现：$S_2(0)$ 与 $S_2(1)$ 在 $L \to \infty$ 时收敛于同一常数，验证了 1D 情况下的纠缠均分，尽管收敛速度极慢（需 $L \sim 10^{20}$）。

  • 2D TFIM (正方晶格临界点 $h_c/J \approx 3.044$)：

    • 关键发现：在 2D 情况下，无序算子遵循面积律，并带有对数修正：$-\ln\langle X \rangle = aL + b \ln L + c$。
    • 均分证据：不同部门的纠缠熵差值在 $L \approx 500$ 时迅速趋于零。这是首次在 (2+1)D 相互作用体系中观察到纠缠均分的数值证据。

2.2 1D 海森堡链 (Heisenberg Chain)

针对具有连续 $U(1)$ 对称性的系统，研究者测量了磁化强度部门的纠缠分布。

  • 电荷分布 $P(q)$：呈现高斯分布，随着尺寸 $L$ 增加，展宽增加。
  • 电荷方差 $\text{Var}_q$：随 $\ln L$ 线性增长，斜率为 0.0547(1)，符合 Luttinger 液体理论。
  • 数值稳定性：通过离散余弦变换处理 $U(1)$ 算子，在 $L=768$ 的大规模系统下，依然保持了极高的信噪比。

2.3 性能数据总结

  • 系统尺寸：1D 可达 $L=1024$，2D 可达 $512 \times 512$（区域 $A$ 为 $256 \times 512$）。
  • 逆温度 $\beta$：取 $4L$ 或 $2L$ 以确保达到基态。
  • 精度：在 $Z_2$ 部门的探测中，统计误差控制在 $10^{-4}$ 量级，足以捕捉微弱的对数修正项。

3. 代码实现细节、复现指南与开源资源

3.1 核心算法：随机级数展开 (SSE)

该工作主要基于 SSE QMC 算法。复现该工作的技术栈建议如下：

  • 编程语言：C++ (推荐用于核心计算) 或 Julia。
  • 并行框架：MPI (用于多副本并行采样) + OpenMP。

3.2 复现关键步骤

  1. 构造副本流形：
    • 维护两个格点数组副本 lattice_1 和 lattice_2。
    • 在 SSE 的算符环更新（Operator Loop Update）中，对于属于区域 $A$ 的格点，其虚时边界条件需要交叉连接：副本1的末端 -> 副本2的开头，反之亦然。
  2. 测量无序算子：
    • 对于 $Z_2$ 对称性，测量 $X = \prod_{i \in A} \sigma_i^x$。在 SSE 中，这涉及到统计跨越区域 $A$ 的传播子特性。
    • 对于 $U(1)$ 对称性，测量 $e^{i\alpha \sum S_i^z}$。由于 $S^z$ 是 SSE 的自然基底，这只是一个简单的相位累积过程。
  3. 比率法 (Ratio Method)：
    • 由于直接计算 $Z_n(\alpha)$ 可能存在方差问题，建议采用增量法，将 $\alpha$ 从 $0$ 逐渐增加到 $\pi$，分段计算比率并相乘。

3.3 推荐软件包与链接

  • ALPSCore：提供了基础的 QMC 框架和格点定义。GitHub Link
  • SSE 基础代码参考：可参考 Anders Sandvik 的经典 SSE 实现（文章引用 [45, 46]）。
  • 项目相关 Repo：作者在文中提到了北京万源开天（PARATERA）的算力支持，具体的算法优化逻辑可参照 cond-mat/0210462 中关于纠缠熵测量的部分。

4. 关键引用文献与局限性评论

4.1 关键引用

  1. Calabrese & Cardy (2004) [1]：奠定了 1D 纠缠熵计算的基础。
  2. Goldstein & Sela (2018) [16]：首次正式提出对称性分辨纠缠（SRE）的理论框架。
  3. Xavier et al. (2018) [17]：对纠缠均分原理的早期探讨。
  4. Zhao, Yan, Cheng & Meng (2021) [36]：在 2D 体系中使用 QMC 测量无序算子的先驱工作。

4.2 工作局限性评论

尽管本工作取得了突破，但仍存在以下局限：

  • 符号问题（Sign Problem）：作为 QMC 方法，该框架依然无法直接应用于具有严重符号问题的费米子哈伯德模型或高度受挫的量子磁体。
  • 傅里叶变换的噪声：对于高电荷部门（电荷 $q$ 较大时），$P(q)$ 极小，傅里叶变换后的相对误差会指数级放大。这意味着该方法目前只能稳定探测低电荷部门的纠缠分布。
  • 边界效应：2D 体系中 smooth boundary 的定义在离散格点上仍有微小歧义，这可能会影响对数修正项 $b \ln L$ 的精确提取。

5. 其他必要补充：物理直觉与未来方向

5.1 物理直觉：为什么“均分”会发生？

纠缠均分的本质是系统的低能有效理论（如 CFT 或 Luttinger Liquid）中纠缠谱的普适性。在临界点，长程关联占主导，局部电荷的波动被平滑化，导致不同电荷部门感受到的有效希尔伯特空间维度在对数精度上是等价的。

5.2 未来研究方向：非阿贝尔对称性

目前的工作集中在阿贝尔对称性（$Z_2, U(1)$）。然而，许多量子物质（如非阿贝尔拓扑序、$SU(N)$ 磁体）具有更复杂的对称性。如何将该 QMC 框架扩展到非阿贝尔群的不可约表示（Irreducible Representations）分解，将是下一个竞争高地。

5.3 总结

本项研究不仅是一个算法的进步，更是对量子多体理论的一次深度验证。它证明了即便是在强相互作用的 2D 极限下，对称性依然以一种极其优雅且均等的方式约束着纠缠的分布。对于量子化学家而言，这种基于副本流形的 QMC 采样思路，或许能为分子体系中的电子相关能拆解提供新的灵感。