Ergodicity breaking in matrix-product-state effective Hamiltonians

这是一份关于论文《Ergodicity breaking in matrix-product-state effective Hamiltonians》（矩阵乘积态有效哈密顿量中的遍历性破缺）的详细技术总结。

1. 研究背景与问题 (Problem)

• 核心挑战：在相互作用的量子多体系统中，热化（Thermalization）及其失效（如多体局域化 MBL 和量子多体疤痕 QMBS）主要由中谱本征态（mid-spectrum eigenstates）决定。然而，这些本征态通常具有体积律纠缠（volume-law entanglement），导致传统的数值方法（如精确对角化 ED）受限于系统尺寸（通常 $N < 30$），难以研究大尺度下的遍历性破缺现象。
• 现有方法的局限：
  • 矩阵乘积态（MPS）和密度矩阵重整化群（DMRG）擅长处理低纠缠的基态或低激发态，但传统上被认为无法准确描述高纠缠的热相。
  • 虽然已有 DMRG-X（针对 MBL）和 DMRG-S（针对 QMBS）等变分方法可以靶向特定激发态，但它们通常只关注单个本征态，缺乏对整个能谱统计特性（如能级间距分布、谱形因子）的系统性探测能力。
• 科学问题：能否利用 routinely used 的 DMRG 有效哈密顿量（Effective Hamiltonian, $H_{\text{eff}}$）来提取远离平衡态的量子多体系统的谱信息，特别是用于探测遍历性破缺（强局域化和弱遍历性破缺）？

2. 方法论 (Methodology)

作者提出了一种利用 DMRG 有效哈密顿量 ($H_{\text{eff}}$) 作为谱探针的新框架：

• $H_{\text{eff}}$ 的构建：
  • 基于一个参考 MPS 态 $|\psi(A)\rangle$（通常通过 DMRG-S 算法获得的近似本征态）。
  • 将总哈密顿量 $H$ 表示为矩阵乘积算符（MPO）。
  • 将 MPO 与 MPS 在除 $M$ 个中心格点外的所有位置进行缩并（contraction），得到一个作用在中心格点及其左右虚拟指标上的有效哈密顿量 $H_{\text{eff}}$。
  • $H_{\text{eff}}$ 的维度为 $\chi^{M} d^M$（其中 $\chi$ 是最大键维，$d$ 是物理维）。
• 谱分析流程：
  1. 构造参考态：使用 DMRG-S 算法在特定能量 $E_{\text{target}}$ 附近寻找近似本征态 $|\Psi_0\rangle$。
  2. 对角化 $H_{\text{eff}}$：对构建好的 $H_{\text{eff}}$ 进行精确对角化，获得一组本征值 $E_j$ 和本征矢 $v_j$。
  3. 重构物理态：将 $v_j$ 重塑为 MPS 张量，结合固定的左右 Schmidt 基，构造出整个希尔伯特空间中的正交态 $|\psi(A, v_j)\rangle$。
  4. 物理量计算：计算这些重构态的纠缠熵（$S_E$）、能级间距比（$\langle r \rangle$）、谱形因子（SFF）以及能量方差（用于筛选高质量本征态）。
• 关键创新点：尽管 $|\psi(A, v_j)\rangle$ 并非原哈密顿量 $H$ 的精确本征态，但作者证明 $H_{\text{eff}}$ 的谱统计特性能够编码原系统的遍历性破缺特征，即使在全局波函数近似不完美时，局部关联和谱统计依然保持准确。

3. 主要研究结果 (Key Results)

作者分别在两个经典模型中验证了该方法的有效性：

A. 随机场 XXZ 链（多体局域化 MBL）

• 模型：一维自旋 1/2 XXZ 链，带有随机磁场 $h_j \in [-W, W]$。
• 能级统计：
  • 计算了平均能级间距比 $\langle r \rangle$。结果显示，随着无序度 $W$ 增加，$\langle r \rangle$ 从高斯正交系综（GOE, $\approx 0.536$）平滑过渡到泊松分布（Poisson, $\approx 0.386$），准确捕捉了从热化相到 MBL 相的转变。
  • 谱形因子（SFF）$K(\tau)$ 在弱无序下呈现典型的“凹陷 - 斜坡 - 平台”（dip-ramp-plateau）结构，而在强无序下凹陷消失，符合 MBL 特征。
• 纠缠熵：
  • 在热化相，中谱态的纠缠熵遵循体积律（$S_E \propto L_{\text{eff}}$）；在 MBL 相，遵循面积律（$S_E \propto \text{const}$）。
  • 通过有限尺寸标度分析（Finite-size scaling），利用纠缠熵方差 $\sigma(S_E)$ 的峰值确定了临界无序度 $W_c \approx 4.4$ 和临界指数 $\nu \approx 0.89$，与 ED 结果高度一致。
• 遍历性气泡（Ergodic Bubbles）：
  • 在 End Matter 中，作者模拟了嵌入在局域化背景中的低无序区域（“热 grains"）。
  • 结果显示，$H_{\text{eff}}$ 能够空间分辨地探测到热 grains 对周围的影响范围。在接近相变点时，热 grains 能显著增强周围区域的纠缠和谱关联（预示雪崩不稳定性）；而在深 MBL 相中，这种影响被局域化抑制。

B. PXP 模型（量子多体疤痕 QMBS）

• 模型：描述里德堡原子阵列的 PXP 模型，具有动力学约束（相邻自旋不能同时向上）。该模型是混沌的，但存在违反 ETH 的“疤痕”态。
• 疤痕塔结构：
  • 通过靶向特定能量 $E_{\text{target}}(n)$ 构建 $H_{\text{eff}}$，成功复现了 QMBS 的“塔”状能级结构（能量间隔 $\Delta E \approx 1.331$）。
  • 重构的本征态与初始态 $|Z_2\rangle$（$|\uparrow\downarrow\uparrow\downarrow\dots\rangle$）具有高重叠度，且表现出低纠缠特征。
• 不对称性发现：
  • 发现了一个有趣的现象：靶向能量 $E_{\text{target}}(n)$ 构建的 $H_{\text{eff}}$ 能准确捕捉能量更低的疤痕塔（$n' < n$），但对能量更高的塔（$n' > n$）捕捉效果较差。
  • 作者将此归因于 QMBS 塔近似满足的 $su(2)$ 代数结构。
• 弱遍历性破缺：该方法成功区分了热化背景中的低纠缠非热态，证明了其在探测弱遍历性破缺方面的能力。

4. 关键贡献 (Key Contributions)

1. 重新定义 DMRG 有效哈密顿量的用途：将通常用于变分基态优化的 $H_{\text{eff}}$ 转化为一种强大的谱探针，用于研究远离基态的量子多体物理。
2. 突破系统尺寸限制：该方法能够处理 $N=50$ 甚至更大的系统（取决于 $\chi$），远超精确对角化（ED）的能力范围，同时保留了探测遍历性破缺所需的谱统计信息。
3. 统一框架：提供了一个统一的框架，既能探测强遍历性破缺（MBL，完全抑制热化），也能探测弱遍历性破缺（QMBS，部分违反 ETH）。
4. 空间分辨能力：展示了该方法在探测 MBL 相中罕见的“热 grains"及其引发的雪崩不稳定性方面的独特优势，这是传统全局方法难以做到的。
5. 理论验证：通过解析可解的自旋 -1 XY 模型，从理论上证明了 $H_{\text{eff}}$ 能够重构精确的疤痕塔结构，为 PXP 模型的结果提供了坚实的理论支撑。

5. 意义与展望 (Significance)

• 方法论突破：这项工作解决了 MPS 方法难以处理高纠缠热态的长期痛点。它表明，即使全局波函数近似不完美，$H_{\text{eff}}$ 的局部谱信息依然能忠实反映系统的物理本质（如能级排斥、纠缠标度）。
• 实验关联：该方法适用于当前实验平台（如冷原子、超导电路、里德堡原子阵列）中观测到的 MBL 和 QMBS 现象，为解释实验数据提供了新的数值工具。
• 未来方向：
  • 结合机器学习技术识别遍历性破缺相。
  • 应用于其他类型的遍历性破缺（如 Stark MBL、希尔伯特空间碎片化、无 Disorder 局域化）。
  • 在量子处理器上利用变分激发态制备技术直接提取谱信息。

总结：该论文证明了 DMRG 有效哈密顿量不仅是基态研究的工具，更是探索量子多体系统非平衡动力学、遍历性破缺及相变机制的通用且强大的谱学探针，极大地扩展了张量网络方法在量子多体物理中的应用边界。