来源论文: https://arxiv.org/abs/2605.26208v1 生成时间: May 28, 2026 06:35

从虚拟走向现实：基于张量网络将缠绕场映射至物理局部算符，攻克量子多体纠缠测量难题

0. 执行摘要

在现代量子化学与凝聚态物理的交汇处，量子纠缠已成为诊断强关联电子结构、拓扑物相以及量子临界现象的核心工具。特别是在多中心过渡金属催化剂、键断裂过程以及激发态非绝热耦合等强关联分子体系中，传统的密度泛函理论（DFT）往往失效，而基于张量网络方法（如矩阵乘积态 MPS 和密度矩阵重整化群 DMRG）的纠缠熵分析（如单轨道和双轨道熵）则展现出极强的解释力。然而，理论界长期面临一个痛点：表征多体系统纠缠特性的 Rényi 纠缠熵，其理论计算通常高度依赖于“缠绕场”（Twist Fields）和“虚拟交换算符”（Virtual Swap Operators）。这些算符作用在张量网络的辅助键维度（Virtual Bond Dimension）上。对于理论计算，这是一种优雅的数学构筑；但对于量子模拟硬件（如里德堡原子阵列、超导比特、离子阱）以及真实的物理测量而言，这些辅助维度是根本无法直接访问的“虚拟存在”。

近日，Andrea Bulgarelli、Marco Panero、Paolo Stornati 和 Luca Tagliacozzo 撰写的重磅论文 “Mapping twist fields to local operators via tensor networks” 彻底打破了这一壁垒。该研究巧妙地利用了张量网络状态中的单射性（Injectivity）性质，提出了一种通用的、严密的数学框架，成功将作用在虚拟层面的缠绕算符“提升”（Lift）并精确映射为直接作用在物理希尔伯特空间上的局部物理算符。这意味着，原本必须在虚拟键上进行的抽象交换操作，现在可以等价地写成物理可观测量的线性组合（如 Pauli 字符串）。这一成果不仅为量子模拟平台提供了一种高可扩展、无需全状态层析（Tomography）的 Rényi 纠缠熵测量新范式，也为经典张量网络量子化学计算中快速计算非连续轨道、多不相交子系统的纠缠性质提供了革命性的工具。

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

1.1 核心科学问题：为何虚拟缠绕场无法在实验中直接测量？

在量子信息与量子多体物理中，给定一个子系统 $A$ 及其补空间 $B$，其第 $n$ 阶 Rényi 纠缠熵定义为：

$$S_n(A) = \frac{1}{1-n} \log \text{Tr} \rho_A^n, \quad \rho_A = \text{Tr}_B \rho$$

当 $n \to 1$ 时，它收敛到经典的 von Neumann 纠缠熵。在共形场论（CFT）中，计算 $S_n(A)$ 的经典理论方法是“复制技巧”（Replica Trick），即在 $n$ 个相互独立的系统副本（Replicas）上构筑黎曼面，并通过在子系统 $A$ 的边界插入分支点缠绕场（Branch-point Twist Fields）来连接不同的副本。在张量网络（特别是 MPS）的语言中，这一理论被巧妙地实现为：在虚拟层引入一个循环交换算符（Cyclic Swap Operator），直接作用在张量网络的辅助键索引上。对于 $n=2$ 的情况，该交换算符通过在虚拟辅助空间将两个副本的键索引相互交叉，即可高效算得 $\text{Tr}\rho_A^2$。

然而，这带来了一个致命的科学问题：辅助键维度（Virtual Degrees of Freedom）是经典计算中为了压缩波函数信息而引入的数学辅助量，在真实的量子实验或物理系统中不存在物理对应物。 实验物理学家无法将一个激光脉冲或微波脉冲作用在“虚拟键”上，因此，缠绕场的虚拟算符公式在实验上是不可直接操作的。若要在实验中测量 Rényi 熵，传统的做法是使用作用在全系统物理希尔伯特空间上的全局物理交换算符（Physical Swap Operator）。然而，这种全局物理交换算符的支撑集覆盖了整个子系统 $A$。当子系统 $A$ 尺度很大、或者是多个空间不相交的碎片化区间时，制备多副本并实施全局交换操作的实验复杂度呈指数级上升，极其脆弱且难以扩展。

因此，该论文的核心科学问题是：能否找到一个作用在物理希尔伯特空间上的“局部”算符，其支撑集仅局限于子系统的边界（而不是填满整个子系统），且其物理期望值能够完美再现虚拟缠绕场的作用？

1.2 理论基础：张量网络的单射性（Injectivity）

解决这一问题的关键在于**单射矩阵乘积态（Injective MPS）**的深层数学性质。一个 MPS 在某个阻塞尺度（Blocking Scale）$l$ 上是单射的，意味着由 $l$ 个连续的局部 MPS 张量拼接而成的映射：

$$\mathcal{A}^{(l)} : \mathbb{C}^{\chi} \otimes \mathbb{C}^{\chi} \to \mathbb{C}^{d^l}$$

（其中 $\chi$ 为虚拟键维度，$d$ 为单轴物理维度）是一个单射映射。换言之，该映射将虚拟希尔伯特空间（维度为 $\chi^2$）无损地嵌入到了一个局域物理希尔伯特空间（维度为 $d^l$）中。只要物理维度足够大，满足 $d^l \ge \chi^2$，单射性在物理上几乎是普适成立的。对于任何具有有限关联长度 $\xi$ 的量子态，随着阻塞长度 $l$ 相比于 $\xi$ 按指数级增长，系统将迅速达到单射极限（Injectivity Limit）。

一旦单射性确立，意味着虚拟空间上的任何线性操作（包括虚拟交换算符）都可以被唯一且忠实地“提升”（Lift）并表示为物理空间上的局域操作。这就是本工作最核心的理论基石。

1.3 技术难点与方法细节：如何将算符从虚拟层“提升”到物理层？

1.3.1 矩阵乘积态（MPS）中的精确映射

考虑 $n=2$ 副本的情况。假设每个 MPS 局部张量为 $A$，其物理维度为 $d^l$（已阻塞 $l$ 个站点），虚拟维度为 $\chi$。我们将 $A$ 视为一个形状为 $\chi^2 \times d^l$ 的矩阵，其行索引由左右虚拟索引 $(\mu, \nu)$ 构成，列索引由物理索引 $i$ 构成。

由于单射性，我们可以利用奇异值分解（SVD）来显式构造 $A$ 的右伪逆（Right Pseudo-inverse，记为 $\bar{A}$）。具体过程如下： 对张量 $A$ 进行 SVD：

$$A = U S V^\dagger$$

其中 $U$ 是 $\chi^2 \times \chi^2$ 的斜正交矩阵，$S$ 是 $\chi^2 \times \chi^2$ 的奇异值对角矩阵（在单射极限下其奇异值全部非零），$V^\dagger$ 是 $\chi^2 \times d^l$ 的半正交矩阵（满足 $V^\dagger V = I_{\chi^2}$）。则 $A$ 的广义逆矩阵 $\bar{A}$（形状为 $d^l \times \chi^2$）定义为：

$$\bar{A} = V S^{-1} U^\dagger$$

这种构造保证了物理收缩关系：

$$\sum_{i=1}^{d^l} \bar{A}_{i, \mu\nu} A_{i, \alpha\beta} = \delta_{\mu\alpha} \delta_{\nu\beta}$$

即 $A$ 与 $\bar{A}$ 在物理索引上的收缩直接退化为虚拟空间上的单位算符（Identity）。这一神奇的性质允许我们将任何作用在虚拟索引上的算符，包裹在 $A$ 和 $\bar{A}$ 之间，从而转化为作用在物理索引 $i$ 上的等价算符。

基于此，作者定义了前向物理缠绕算符（Forward Twist Operator） $T^{(f)}_{ijkn}(x)$。对于两个副本，它作用于位置 $x$ 处的物理位点（每个副本各占 $l$ 个站点，总共 $2l$ 个物理站点）：

$$T^{(f)}_{ijkn}(x) = \bar{A}^i_{\mu\nu}(x) A^j_{\mu\sigma}(x) \bar{A}^k_{\rho\sigma}(x) A^n_{\rho\nu}(x)$$

通过精妙的张量图形收缩可以证明（见论文图 1(b) 与图 4），这一算符的期望值在单射极限下严格等价于在虚拟层插入的交换算符，从而完美给出：

$$\langle T^{(f)} \rangle = \text{Tr}\rho_A^2$$

对于反向缠绕算符 $T^{(b)}(x+1)$，可以通过对称地交换左侧虚拟键来实现，其物理构建逻辑完全一致。

1.3.2 2D 投影纠缠配对态（PEPS）中的变分推广

在二维空间（PEPS）中，单射性和正则形式（Canonical Form）的性质远没有一维（MPS）那样强大且易于解析处理。在 2D 系统中，若要计算一个区域 $A$ 的纠缠熵，虚拟层上的缠绕算符不再是一个点算符，而是沿着区域边界 $\partial A$ 延伸的弦状算符（String-like Operator）。

直接进行张量求逆在 2D 情况下会遭遇维数灾难。例如，即使对于最小的 PEPS 键维度 $D=2$，为了达到单射性，通常需要对一个 $4 \times 4$ 的斑块（16个自旋）进行阻塞。这导致需要求逆的物理算符维度高达 $2^{16} \times 2^{16}$，其对应的 Pauli 字符串展开将包含高达 $10^{18}$ 项，在经典计算和实验上都完全不可行。

为了解决这一难点，作者提出了一种基于变分最小化（Variational Minimization）的局部近似策略： 将问题定义在一个半无限圆柱体上，其边界状态由一维边界矩阵乘积算符（MPO）$\rho_L$ 和 $\rho_R$ 刻画（见论文图 7）。边界上的 Rényi 熵由作用在虚拟键上的局部交换算符 $S_D$ 表达。为了寻找一个仅作用在两个副本物理层上的局域算符 $O_T$，作者将其表述为一个在局部物理算符子空间中的最小二乘变分问题：

$$O_T = \text{argmin}_O \| O |L_T\rangle - (I \otimes S_D)|L_T\rangle \|^2$$

其中 $|L_T\rangle$ 是局部张量网络在应用虚拟交换之前的状态向量，$(I \otimes S_D)|L_T\rangle$ 是应用虚拟交换后的目标向量。通过将 $O_T$ 在局域物理算符基组（如 Pauli 基组）$P_\alpha$ 上展开：

$$O_T = \sum_{\alpha} a_\alpha P_\alpha$$

该变分问题直接转化为一组线性方程（法方程，Normal Equations）：

$$\sum_{\beta} G_{\alpha\beta} a_\beta = b_\alpha$$

其中 $G_{\alpha\beta} = \langle v_\alpha | v_\beta \rangle$ 为 Gram 矩阵，$b_\alpha = \langle v_\alpha | w \rangle$，并且定义了辅助态 $|v_\alpha\rangle = P_\alpha |L_T\rangle$ 和 $|w\rangle = (I \otimes S_D)|L_T\rangle$。若 Gram 矩阵奇异或病态，作者引入了 Tikhonov 正则化：

$$a = (G + \lambda I)^{-1} b$$

这一变分方法极其漂亮地绕过了直接求逆的维数灾难，在保证高精度的同时，极大地压缩了所需物理算符的支撑集和项数。

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

为了严谨地评估该映射方法的精确度、收敛性以及跨体系的可转移性（Transferability），作者选取了经典的 一维横场伊辛模型（1D TFIM） 作为 Benchmark 体系。其哈密顿量为：

$$H = -\sum_{x=1}^{L-1} \sigma^z_x \sigma^z_{x+1} - g \sum_{x=1}^L \sigma^x_x$$

其中，当耦合常数 $g_c = 1$ 时，系统处于量子临界点。以下是核心测试数据及物理分析：

2.1 单射极限下的指数级快速收敛

作者在链长 $L=120$、不同耦合强度 $g = 4.0$（远离临界区，短关联长度）、$g = 2.0$ 以及 $g=1.5$（接近临界区，长关联长度）下，测试了相对误差：

$$\text{Relative Error} = \frac{|\text{Exact } S_2 - \text{Approximated } S_2|}{\text{Exact } S_2}$$

其中，近似值由物理缠绕算符的期望值 $\langle T^{(f)} \rangle$ 算出。测试结果（见原论文图 2）表明：

  • 对于 $g=4.0$：由于系统关联长度极短，仅需阻塞长度 $l=4$，相对误差便瞬间坠落至机器精度下限（$\sim 10^{-11}$）。
  • 对于 $g=2.0$：随着关联长度增加，单射性需要稍大的阻塞尺度。在 $l=5$ 时，全链上除边界极小区域外，相对误差均降至 $10^{-10}$ 以下。
  • 对于接近临界点的 $g=1.5$：系统关联长度显著增大。在 $l=1$ 到 $l=5$ 时，仍有一定系统误差；但当 $l=6$ 时，整个物理链上的相对误差也完美收敛至 $10^{-11}$ 级别的机器精度。

这强有力地证明了：只要阻塞长度 $l$ 达到或超越系统关联长度尺度的单射极值，通过张量求逆构造的物理缠绕算符就是绝对精确的。

2.2 复合不相交区间的普适计算

利用 CFT 理论，通过多次插入构造的物理缠绕算符，可以极为简便地计算多个不相交（Disjoint）区间的纠缠。作者展示了双点关联函数与四点关联函数的计算结果（见原论文图 3）：

  • 图 3(a) 绘制了两个缠绕算符的关联函数作为其距离 $r$ 的函数，计算得到的第二阶 Rényi 熵 $S_2(r)$ 与精确 DMRG 结果的误差在整个 $r/L$ 范围内低于 $10^{-14}$。
  • 图 3(b) 展示了四个缠绕算符插入下，两个独立区间的纠缠熵随交叉比（Cross Ratio）$\eta = (x_{12}x_{34})/(x_{13}x_{24})$ 的变动曲线，同样达到了 $10^{-14}$ 的极致精度。

2.3 正交中心（Center of Orthogonality）下的免求逆捷径

论文指出一个极具实用价值的推论：如果我们将局部张量 $A$ 精确选择在张量网络的正交中心（Center of Orthogonality）上，那么即使完全不满足单射条件（即 $l=1$），公式（3）也严格成立！ 这是因为正交中心自带半正交收缩性质（见论文图 4），无需通过 SVD 进行伪逆计算。这一发现极大地简化了经典数值计算中提取特定单点/双点纠缠谱的复杂度。

2.4 颠覆性的发现：可转移性（Transferability）与标度律

对于量子模拟实验而言，在庞大的靶系统（Target System，尺寸为 $L_{\text{tar}}$）上直接校准或计算缠绕算符是不现实的。作者提出一个天才的猜想：能否在较小尺寸的参考系统（Reference System，尺寸为 $L$）上提取缠绕算符，然后直接“套用”到大系统上？

作者在 $g=4.0$ 和 $g=1.5$ 下进行了深入测试，将小系统（$L < L_{\text{tar}}$）提取的算符 $T_L$ 作用于 $L_{\text{tar}} = 100$ 的大系统上，观察其相对误差（见原论文图 5）：

  • 结果令人震惊：只要小系统的尺寸 $L$ 超越一个临界阈值 $L_c$，其在大系统上的估算误差便会以指数级衰减至零！ 且这一衰减速率完全独立于靶系统尺寸 $L_{\text{tar}}$。
  • 作者进一步对临界尺度 $L_c$ 随系统关联长度 $\xi$ 的变化进行了精确拟合，得到了精美的幂律关系式（式 5）： $$L_c(\xi) = A \xi^\omega + k$$ 经自助法（Bootstrap）拟合得到的物理参数为： $$A = 37(4), \quad \omega = 0.68(5), \quad k = 8(4)$$ 其中，指数 $\omega$ 被认为与缠绕场的共形标度维度（Scaling Dimension）有深层物理关联。

这一可转移性不仅极具理论美感，更具备巨大的实用价值：这意味着研究人员只需在经典超级计算机上对一个微小的分子或自旋链片段（$L \sim L_c$）进行 DMRG 计算并提取物理缠绕算符的 Pauli 系数，就能直接拿到超大尺寸的量子模拟器上测量纠缠熵。

3. 代码实现细节、复现指南及开源工具链

为了便于科研人员复现此工作，以下提供基于 Python 的完整实现逻辑与核心算法指南。我们将结合广受欢迎的张量网络库 TensorNetwork 或 quimb 来阐述关键步骤。

3.1 核心算法复现步骤

  1. 基态求解：利用 DMRG 算法求解 1D TFIM 的基态，得到其 MPS 表示，确保其处于右正则（Right-canonical）形式。
  2. 张量阻塞（Tensor Blocking）：选择一个中心位点 $x$，将左右各 $l/2$ 个局部三阶张量 $A$ 收缩为一个大张量 $A^{(l)}$，其形状为 $(\chi, \chi, d^l)$。
  3. 矩阵重塑与 SVD 求逆：
    • 将 $A^{(l)}$ 重塑为矩阵 $\mathbf{M}$，形状为 $\chi^2 \times d^l$（其中行索引合并了左右两个虚拟维度）。
    • 对 $\mathbf{M}$ 进行奇异值分解：$\mathbf{M} = U S V^\dagger$。
    • 设定阈值 $\epsilon \sim 10^{-12}$，过滤接近于零的奇异值，构造对角阵 $S^{-1}$。
    • 构造伪逆矩阵 $\mathbf{\bar{M}} = V S^{-1} U^\dagger$（形状为 $d^l \times \chi^2$）。
    • 将 $\mathbf{\bar{M}}$ 重新塑形回张量 $\bar{A}^{(l)}$，其形状为 $(d^l, \chi, \chi)$。
  4. 算符收缩（物理缠绕算符构建）：根据公式（2），在虚拟键上循环收缩两个 $\bar{A}^{(l)}$ 与两个 $A^{(l)}$，最终得到一个形状为 $(d^l, d^l, d^l, d^l)$ 的四阶张量，此即物理缠绕算符 $T^{(f)}$。
  5. Pauli 字符串分解：利用 Kronecker 积，将四阶张量 $T^{(f)}$ 投影到 $2l$ 阶的 Pauli 基组上，筛选出系数显著非零的局部 Pauli 字符串。

3.2 关键代码复现片段（Python 伪代码说明）

import numpy as np
from scipy.linalg import svd

def construct_physical_twist_operator(A_tensor, l, chi, d, threshold=1e-12):
    """
    A_tensor: 输入的已阻塞连续 l 个位点的 MPS 张量，形状为 (chi, chi, d**l)
    """
    # 1. 重塑为矩阵形式 (chi*chi, d**l)
    M = A_tensor.reshape(chi * chi, d**l)
    
    # 2. 进行 SVD 得到 U, S, Vh
    U, S, Vh = svd(M, full_matrices=False)
    
    # 3. 计算奇异值的逆，低于 threshold 的设为 0
    S_inv = np.zeros_like(S)
    non_zero = S > threshold
    S_inv[non_zero] = 1.0 / S[non_zero]
    
    # 4. 构造广义逆矩阵 M_bar = V * S_inv * U^dagger
    # V 是 Vh 的共轭转置
    V = Vh.T.conj()
    M_bar = V @ np.diag(S_inv) @ U.T.conj() # 形状为 (d**l, chi*chi)
    
    # 5. 重塑回逆张量 A_bar，形状为 (d**l, chi, chi)
    A_bar = M_bar.reshape(d**l, chi, chi)
    
    # 6. 构建物理缠绕算符 T^(f) 
    # T_(i,j,k,n) = sum_{mu, nu, sigma, rho} A_bar[i, mu, nu] * A[j, mu, sigma] * A_bar[k, rho, sigma] * A[n, rho, nu]
    # 我们使用 numpy.einsum 进行极其直观的高效张量收缩
    T_f = np.einsum('imn, jms, krs, rpn -> ijkn', A_bar, A_tensor, A_bar, A_tensor)
    
    return T_f

# 示例参数
chi = 10
l = 4
d = 2 # 自旋-1/2 体系
# 随机生成一个模拟的单射张量进行测试
np.random.seed(42)
mock_A = np.random.randn(chi, chi, d**l)
T_operator = construct_physical_twist_operator(mock_A, l, chi, d)
print("构建的物理缠绕算符维度为:", T_operator.shape)

3.3 推荐开源软件与 Repo 链接

  1. TeNPy (Tensor Network Python): 论文作者在数值计算中使用的核心库。它不仅对 1D 系统有极为完美的自适应 DMRG 支持，还提供了详尽的矩阵重整化群和有限温计算接口。
  2. quimb: 一个非常现代化、高度优化的张量网络与量子信息计算 Python 库。其对大尺度张量收缩（使用 opt_einsum）和 PEPS 的支持非常友好，特别适合用于实现本文中的 2D 变分最小化方案。
  3. ITensors: 如果你更倾向于 Julia 的极致性能与直观的物理索引设计，ITensors 是复现张量求逆与正交中心收缩的最佳选择。

4. 关键引用文献与局限性批判

4.1 关键引用文献分析

本工作的成功立足于前人多项里程碑式研究的奠基：

  1. Calabrese & Cardy (2004) [5]：奠定了利用共形场论（CFT）和复制技巧计算一维系统纠缠谱的理论大厦，指明了缠绕场与纠缠熵之间的对数标度律关系。
  2. Coser, Tagliacozzo, Tonni (2014) [25]：首次将不相交区间的 Rényi 熵计算推广到 MPS 框架中，指出了虚拟交换算符作为转移矩阵（Transfer Matrix）的核心机制。本论文正是对该工作“无法在物理层直接实验测量”这一核心痛点的完美终结。
  3. Pérez-García et al. (2007) [19, 26]：系统阐明了 MPS 表示的数学结构，特别是**单射性（Injectivity）**和正则形式的严格数学证明，为本文的伪逆构造提供了坚不可摧的数学保障。

4.2 对本工作局限性的客观评论

尽管这是一项兼具理论深度与实验指导意义的杰出工作，但在量子化学及大规模应用场景中，它依然存在以下不容忽视的局限性：

  1. 2D PEPS 的计算高墙依然坚固：虽然作者变分地提出了利用边界 MPO 近似构建 2D 物理缠绕算符的框架，但实际计算成本依然惊人。Gram 矩阵的构建和变分求解本身极其脆弱，若关联长度 $\xi$ 较大，所需的算符支撑集边界将变宽，导致 Pauli 字符串的数量暴增。在实际的 2D 强关联分子盘面模型中，该方法的实用性仍需打上问号。
  2. 高度的“状态依赖性”（State-dependency）：这是该方法最显著的物理局限。构造的物理缠绕算符 $T$ 的系数 $a_\alpha$ 是完全依赖于特定态（State-dependent）的。这意味着，对于分子体系，如果你改变了键长（比如拉伸化学键），或者改变了外部磁场，你必须重新计算和校准一套全新的 Pauli 系数。它不像全局物理交换算符那样具有普适性（Universal）。虽然“可转移性”在一定程度上缓解了这一问题，但在拓扑相变点附近或高度动态的非平衡态演化过程中，如何实时、稳健地更新算符仍是个巨大挑战。
  3. 实验测量中的“噪声放大器效应”：尽管该方法将物理测量的支撑集从“子系统全域”成功压缩到了“边界局域”，极大减少了所需的测量点数。但是，伪逆 $\bar{A} = V S^{-1} U^\dagger$ 中包含奇异值倒数 $S^{-1}$。当系统接近单射极限边缘时，较小的奇异值会导致 $S^{-1}$ 中的某些项数值极大。在真实的量子模拟器中，这会成倍地放大物理测量中的系统噪声与统计涨落，导致最终估算的 Rényi 熵具有极大的方差（Variance），亟需引入更先进的误差缓解（Error Mitigation）技术。

5. 补充探讨：对强关联分子体系与量子模拟的启示

5.1 对强关联量子化学轨道选择的深远启示

在强关联量子化学（如过渡金属催化剂、多核锰簇等）中，构建活动空间（Active Space）是所有高级电子相关计算（如 CASSCF, DMRG-CASSCF）的核心生死线。传统上，化学家依赖于经验或轨道的自然占用数（Natural Orbital Occupancies）来选择活动空间。然而，现代量子化学理论表明，单轨道熵（Single-orbital Entropy）和双轨道互信息（Two-orbital Mutual Information）才是筛选活性轨道、揭示化学键断裂机制最科学的物理判据。

$$\mathcal{I}_{p,q} = S_1(p) + S_1(q) - S_2(p,q)$$

传统的 DMRG 计算在提取空间不相邻、碎片化的轨道对 $p, q$ 的双轨道纠缠熵 $S_2(p,q)$ 时，往往面临繁琐的基底重组。本论文提出的“物理缠绕算符映射方法”提供了一条全新的思路：由于分子轨道可以被视为一维 MPS 链上的不同位点，我们可以直接利用局部物理算符 $T$，在一维轨道链上通过极其简单、局域的物理期望值测量，直接提取任意轨道碎片（即使它们空间上相隔极远）之间的 Rényi 纠缠熵。 这将极大加速 DMRG-CASSCF 计算中活性轨道的自动筛选流程。

5.2 费米子体系的映射延展

化学体系本质上是费米子系统。本文的代码与公式主要基于自旋/比特（Spin/Qubit）体系。在推广到真实的费米子波函数时，必须额外小心：

  • 费米子算符在非局域交换时会产生 Jordan-Wigner 弦因子（Sign String）。
  • 因此，在费米子张量网络中构造伪逆和物理缠绕算符时，必须将费米子奇偶性（Parity）和费米统计规律严密地编织进张量索引中。
  • 令人欣慰的是，由于本工作提出的物理缠绕算符 $T$ 本身是局部支持的（仅作用于 $2l$ 个站点），Jordan-Wigner 弦在算符支撑集外部会自动相互抵消，这预示着该方法在费米子体系（如 Fermi-Hubbard 模拟器和真实的分子体系）中同样具有极高的可行性和无与伦比的优雅度。

5.3 结语：弥合经典算法与现代量子硬件的终极桥梁

本工作最令人兴奋的愿景在于它为经典计算与现代模拟硬件之间架起了一座精美的桥梁。在经典计算中，我们拥有高精度的张量网络波函数，但受限于算力无法模拟极大的体系；在量子模拟器（如里德堡原子、超导芯片）中，我们拥有庞大的、能演化真实多体动力学的物理系统，但却极难测量纠缠等非线性可观测物理量。本研究通过“经典算符预校准 + 硬件局部测量”的共生模式，将两者的优势融为一体：在小参考系统上由经典超级计算机算出高精度的局部物理算符，然后直接部署在数万个原子规模的量子模拟器上进行局部的 Pauli 测量。这一范式必将深刻影响未来强关联化学系统模拟与量子多体纠缠探测的发展轨迹。