来源论文: https://arxiv.org/abs/2605.26219v1 生成时间: May 28, 2026 00:45

经典非平衡相变与量子纠缠图样的桥梁：基于定向渗透与等距张量网络状态（isoTNS）的深度解析

0. 执行摘要

量子相变通常在平衡态多体系统的基态中展开，其特征是纠缠结构的长程关联与拓扑序的涌现。然而，非平衡态动力学中的相变（例如经典随机过程中的吸收相变）是否能为理解凝聚态物理和量子化学中的非平凡量子纠缠提供全新的视角？这正是本博文深度解析的学术前沿成果——由 Julian Boesl、Frank Pollmann 和 Michael Knap 撰写的最新论文《Entanglement Pattern Transition of Quantum States from Directed Percolation》所回答的核心科学问题。

该研究构建了一种精确的映射机制，将经典 $(1+1)$ 维随机元胞自动机——Domany-Kinzel (DK) 自动机——映射到二维等距张量网络状态（isoTNS）中。通过这一映射，经典的定向渗透（Directed Percolation, DP）非平衡相变在量子世界中转化为一种全新的纠缠图样转变：从吸收相中具有系统尺度的、类似于 $W$ 态的非平凡两体长程纠缠，转变到活性相中高度局域化的、平凡的短程纠缠。在临界点处，该量子态展示了在所有空间方向上的代数（幂律）关联，这在传统的等距张量网络研究中是极为罕见的。这一发现不仅突破了稳定平衡态物质相的传统范式，而且为量子化学在开放系统非平衡态动力学模拟、多体疤痕态构建以及量子计算的主动误差校正领域提供了极具启发性的理论框架。

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

1.1 核心科学问题

本研究的核心科学问题在于：非平衡经典动力学中的临界现象如何映射并决定量子多体波函数的纠缠结构？

在量子信息与凝聚态物理中，人们通常利用张量网络状态（Tensor Network States, TNS）来分类和表示不同的量子物相。例如，拓扑序对应于弦网模型（String-Net Models），对称性破缺态对应于 GHZ 态或簇态。然而，这些相通常是处于热力学平衡态的稳定相，其基态简并度在局部微扰下是稳健的。与此相反，存在一类更为“脆弱”的非平凡纠缠态（如 $W$ 态或 Dicke 态），它们不构成宏观热力学极限下的稳定相，但其独特的纠缠特性在量子信息传输和量子计算中具有至关重要的价值。如何在一个统一的、可调参数的框架中系统地产生、控制并研究这些脆弱的非平凡纠缠态，并揭示其背后的相变机制，一直是物理学与量子化学领域的一大技术难题。本文通过引入经典随机自动机到等距张量网络的局部映射，成功实现了这一目标。

1.2 理论基础

1.2.1 定向渗透（Directed Percolation, DP）与吸收相变

定向渗透是非平衡统计物理中最普遍、最经典的临界现象之一。它描述了一个反应-扩散系统在“活跃”（Active）和“吸收”（Absorbing）状态之间的转变。吸收态的定义是：一旦系统进入该状态，由于局部动力学规则的限制，系统将永远无法自发脱离。在低生存概率下，系统不可避免地跌入全局空无的吸收态；而在高生存概率下，系统能够维持一个具有有限活性密度的稳态。在两者的交界处，系统展现出分形几何特征和无标度关联，其临界行为属于 DP 普适类。

1.2.2 等距张量网络状态（isoTNS）

张量网络是表示高维希尔伯特空间中量子态的高效工具。等距张量网络（isoTNS）是 TNS 的一个特殊子集，其局部张量满足等距条件（Isometry Condition）：

$$\sum_{j,k,c,d} (T)_{abcd}^{jk} (T)_{a'b'cd}^{jk*} = \delta_{a,a'}\delta_{b,b'}$$

这一等距性保证了张量网络在进行部分收缩时无需复杂的变分优化即可高效进行，同时它也直接对应于一维量子电路在空间或时间方向上的幺正演化。这使得 isoTNS 能够被极高效率地制备，并允许研究者深入探讨其临界行为。

1.3 技术难点

将经典非平衡相变映射到量子态中存在三大核心理论与技术瓶颈：

  1. 非局域哈密顿量避免：传统的经典-量子映射往往会导致所得的量子母哈密顿量（Parent Hamiltonian）具有长程、非局域的相互作用。在物理和化学上，我们迫切需要局部（Local）且无挫折（Frustration-Free）的哈密顿量，以便能真正模拟真实物理系统。
  2. 等距性约束的维持：确保经典转移矩阵在量子化映射后能够完全保留 isoTNS 的等距性，从而避免在收缩和测量计算中出现指数级复杂度。
  3. 临界方向的限制：先前的 isoTNS 模型通常只能在单一方向（如时间轴）上支持临界代数关联，而在其他空间方向上则表现为平凡的有限关联长度。如何在 isoTNS 中实现全方向的二维代数关联是一个悬而未决的难题。

1.4 方法细节与映射机制

1.4.1 自动机到 isoTNS 的精细映射

本工作始于一个 $(1+1)$ 维的 Domany-Kinzel (DK) 元胞自动机。系统在一个斜置的方格点阵上交替更新偶数和奇数格点。格点 $i$ 的状态（$0$ 表示空无或吸收，$1$ 表示活跃）取决于其前一时刻相邻格点 $j$ 和 $k$ 的状态。定义条件概率如下：

  • $P(1|00) = p_1$
  • $P(1|01) = P(1|10) = p_2$
  • $P(1|11) = p_3$
  • 且 $P(0|jk) = 1 - P(1|jk)$

当 $p_1 = 0$ 时，全空无态 $|00\cdots0\rangle$ 成为严格的吸收态。该经典概率流可以映射为一个二维量子态：

$$|\Psi\rangle = \sum_{\alpha} \sqrt{p_{\alpha}} |\alpha\rangle$$

其中，$\alpha$ 代表所有可能的演化轨迹（时空路径配置），$p_{\alpha}$ 为该轨迹在经典自动机中的发生概率。利用这一映射，作者将 DK 自动机的局部更新规则嵌入到一个 bond 维度 $\chi = 2$ 的二维张量 $T_{DK}$ 中：

$$(T_{DK})_{abcd}^{jk} = \sqrt{P(i|jk)} \delta_{a,j} \delta_{b,k} \delta_{c,i} \delta_{d,i}$$

这里，虚拟索引 $a, b$ 锁定了前一层的格点状态，而虚拟索引 $c, d$ 与物理索引 $i$ 锁定。通过这种构造，概率守恒性 $\sum_i P(i|jk) = 1$ 自然且完美地转化为张量的局部等距条件。这就构成了二维等距张量网络态 $|DK\rangle$。

1.4.2 母哈密顿量的构造

为了在物理上寻找该量子态的自洽物理背景，作者推导了其局部无挫折母哈密顿量 $H_{DK}$。它表示为局部投影算符之和：

$$H_{DK} = \sum_{\vee} B_{\vee}$$

每个投影算符 $B_{\vee}$ 作用于晶格顶点周围的 $8$ 个量子比特上。在 $p_1 = 0$ 的面上，哈密顿量具备一个极为关键的对称性——缺陷数 $N_D$ 的 $U(1)$ 守恒律。由于当 $p_1 = 0$ 时，禁止从双零状态自发产生活跃态（即禁止在零背景上创生粒子），这形成了一个局部缺陷算符：

$$n_{D,\vee} = \frac{1}{16}(1 - Z_{i_1})(1 - Z_{i_2})(1 + Z_{j_1})(1 + Z_{j_2})$$

整个系统的总缺陷数 $N_D = \sum_{\vee} n_{D,\vee}$ 与哈密顿量对易，即 $[H_{DK}, N_D] = 0$。因此，希尔伯特空间被精确地划分为不同的缺陷数扇区，所有的基态和量子转变都在 $N_D = 0$ 的扇区内发生，这为吸收相与活性相之间的量子相变奠定了数学物理基石。

2. 关键 Benchmark 体系，计算所得数据，性能数据

为了验证上述理论映射的有效性，并阐明经典非平衡相变在量子态中的体现，作者使用了一系列具有严格边界条件（Open Boundary Conditions, OBC）以及周期性边界条件（Periodic Boundary Conditions, PBC）的二维晶格作为 Benchmark 系统，进行了大规模的数值收缩与解析计算。

2.1 临界点处的全向代数关联（OBC Benchmark）

作者聚焦于 $p_2 = p_3 = p$ 的参数线（此处属于典型的场地向渗透普适类），其经典临界点已被高精度确定为 $p_c \approx 0.7055$。利用经典采样技术（对于张量网络的对角算符测量），作者在大小为 $3000 \times 3000$ 的庞大系统上测量了归一化关联函数：

$$C^i_{norm}(j) = \frac{\langle(1 - Z_i)(1 - Z_j)\rangle}{2(1 - \langle Z_i \rangle)} = \frac{\overline{n_i n_j}}{\overline{n_i}}$$

计算所得的数据揭示了令人震惊的物理结果：

  • 在时空临界点 $p = p_c$ 处：

    • 在时间方向 $y$（对应自动机的演化方向），关联函数呈现完美的代数幂律衰减：$C^i_{norm}(j) \sim |i-j|^{-\beta/\nu_{\parallel}}$，数值拟合得到的临界指数 $\beta/\nu_{\parallel} \approx 0.1595$。
    • 在空间方向 $x$ 及其对角方向上，关联函数同样呈现清晰的代数衰减：$C^i_{norm}(j) \sim |i-j|^{-\beta/\nu_{\perp}}$，拟合指数 $\beta/\nu_{\perp} \approx 0.2521$。
    • 这一数据无可辩驳地证明了：该 isoTNS 在所有空间方向上均支持临界代数关联，彻底打破了以往 isoTNS 只能在单方向上临界的局限。

  • 在吸收相内 $p = 0.55 < p_c$（测试规模 $L=100$）：

    • 关联在所有方向上均呈指数衰减，关联长度为有限值，表现为经典的吸收态行为。

  • 在活性相内 $p = 0.80 > p_c$（测试规模 $L=3000$）：

    • 关联函数在长程极限下饱和至一个非零的常数，对应于稳态下活跃格点的均匀有限密度。

2.2 周期性边界条件下的纠缠负度（PBC Benchmark）

当系统被赋予双向周期性（即环面拓扑）时，母哈密顿量 $H_{DK}$ 的基态简并度被打开。简并基态包含完全空无的真空态 $|vac\rangle$ 以及与其正交的非平凡基态 $|GS\rangle$。

为了量化 $|GS\rangle$ 在吸收相深处（$p_1 = 0, p_2 = p_3 = \epsilon \ll 1$）的纠缠结构，作者引入了子系统纠缠负度（Subsystem Negativity, $N(\rho_{AB})$）作为计算指标。选取两个在 $y$ 方向上贯穿整个系统、在 $x$ 方向上宽度为 $\ell$、彼此相距 $d$ 的条带区域 $A$ 和 $B$（保证 $\ell > L_y$）。

通过对约化密度矩阵 $\rho_{AB}$ 的部分转置进行严格的谱分析，作者推导并计算了负度的下界：

$$N(\rho_{AB}) \ge \frac{1}{L_x(L_x - 2)} + O\left(\frac{1}{L_x^4}\right)$$

这一极具理论价值的数据表明：

  1. 距离无关性：负度的下界完全不依赖于子系统 $A$ 和 $B$ 之间的空间距离 $d$。即使在热力学极限下将两区域无限拉远，它们之间依然保留常数大小的纠缠。
  2. $W$ 态等价性：标准的 $L$-比特 $W$ 态： $$|W\rangle = \frac{1}{\sqrt{L}} \sum_i X_i |00\cdots0\rangle$$ 其两体纠缠负度为 $N(\rho_{AB}) = \frac{1}{L(L-2)}$。本工作计算的负度在形式和标度上与 $W$ 态精确契合。这表明在吸收相中，量子态本质上是由单根跨越系统闭合路径的“活跃链”构成的超叠加态，呈现出高度稳健的、非平凡的 W 态两体长程纠缠图样。而一旦跃迁到活性相，由于活跃点的扩散和凝聚，纠缠图样迅速坍缩为普通的、平凡的短程局域纠缠。

3.1 算法架构与物理量计算逻辑

要复现本论文的成果，核心需要实现两个层面：一是经典马尔可夫链蒙特卡洛（MCMC）采样，用于获取对角可观测量的关联函数；二是基于张量网络算符重正化群（TRG）或边界矩阵乘积态（bMPS）的纠缠度量计算，用于计算非对角密度矩阵的偏转置和负度。

对于对角可观测量的计算，由于张量网络的等距性质，对角算符的真空期望值可以直接等价于经典 Domany-Kinzel 自动机的统计平局。这意味着，读者不需要直接对庞大的二维张量网络进行昂贵的变分收缩，而是可以通过极低计算成本的经典蒙特卡洛算法来模拟时空动力学，从而以极高的精度复现图 2 中的幂律关联曲线。而对于非对角纠缠物理量（如负度 $N$），则需要显式收缩张量网络来构建约化密度矩阵 $\rho_{AB}$。

3.2 逐步复现指南（Julia 语言实现示例）

Julia 凭借其无缝的多维数组操作和出色的运算性能，是复现本项工作最理想的编程语言。以下提供一个基于 Julia 进行经典动力学采样以及利用 ITensors.jl 构建 isoTNS 张量网络的逻辑框架复现指南：

步骤 1：经典 Domany-Kinzel 自动机的临界演化采样（用于复现图 2）

using Random
using Statistics

# 经典自动机演化参数设置
const Lx = 3000
const Ly = 3000
const p2 = 0.7055  # 临界概率
const p3 = 0.7055

function evolve_dk_automaton(Lx, Ly, p2, p3)
    # 初始化状态矩阵，1表示活跃，0表示空无
    state = zeros(Int8, Lx, Ly)
    # 边界初始化：边缘全部设置为活跃（对应论文中 1 1 ... 1 的边界条件）
    state[:, 1] .= 1 
    
    for t in 2:Ly
        for i in 1:Lx
            # 计算前驱节点索引（带周期性边界条件）
            left = i == 1 ? Lx : i - 1
            right = i
            
            s_left = state[left, t-1]
            s_right = state[right, t-1]
            
            # 转移概率规则
            if s_left == 0 && s_right == 0
                prob = 0.0
            elseif s_left == 1 && s_right == 1
                prob = p3
            else
                prob = p2
            end
            
            state[i, t] = rand() < prob ? 1 : 0
        end
    end
    return state
end

步骤 2：对角关联函数 $C^i_{norm}(j)$ 的高精度计算

function calculate_correlations(state, mid_x, mid_y, max_dist)
    correlations = zeros(max_dist)
    ni_mean = state[mid_x, mid_y]
    
    if ni_mean == 0
        return correlations # 若中心点未激活，则本次采样不贡献归一化关联
    end
    
    for r in 1:max_dist
        # y 方向关联
        nj = state[mid_x, mid_y + r]
        correlations[r] = (ni_mean * nj) / ni_mean
    end
    return correlations
end

步骤 3：利用 ITensors.jl 构建 isoTNS 的局部张量并求取母哈密顿量

在量子物理复现中，我们需要利用张量库来精确定义投影算符 $B_{\vee}$。利用 Julia 的张量网络核心库：

using ITensors

# 定义 DK 张量的物理索引与虚拟索引
function build_dk_tensor(p1, p2, p3)
    # 定义 bond 维度为 2 的索引
    a = Index(2, "Link,a")
    b = Index(2, "Link,b")
    c = Index(2, "Link,c")
    d = Index(2, "Link,d")
    i_phys = Index(2, "Site,i")
    
    T = ITensor(a, b, c, d, i_phys)
    
    # 填充张量元素（遵循映射方程）
    for s_a in 1:2, s_b in 1:2, s_c in 1:2, s_d in 1:2, s_i in 1:2
        # 注意 Julia 索引从 1 开始，对应状态 0 和 1
        jk_state = (s_a-1, s_b-1)
        i_state = s_i - 1
        
        # 只有在物理索引和后向虚拟索引相同时非零（克罗内克 δ 约束）
        if s_c - 1 == i_state && s_d - 1 == i_state
            prob = get_probability(i_state, jk_state, p1, p2, p3)
            T[a=>s_a, b=>s_b, c=>s_c, d=>s_d, i_phys=>s_i] = sqrt(prob)
        else
            T[a=>s_a, b=>s_b, c=>s_c, d=>s_d, i_phys=>s_i] = 0.0
        end
    end
    return T
end

3.3 推荐开源软件包与资源链接

  1. ITensors (Julia/C++)：量子多体物理中功能最强大、使用最广泛的张量网络开源库，极度适合构建本工作中的局部等距张量，进行投影算符的对角化和张量网络收缩。
  2. TensorKit.jl (Julia)：专门针对具有对称性（如同构、U(1) 守恒）的张量网络设计的工具包，极度适合处理本论文中缺陷数 $N_D$ 守恒的扇区计算。
  3. Zenodo 开源数据与代码库：作者在论文中明确指出了其核心计算代码已托管于 Zenodo，读者可直接访问进行数据对比。

4. 关键引用文献及局限性深度点评

4.1 关键引用文献分析

本研究的理论支柱深深植根于以下里程碑式的学术成果之中，对这些文献的梳理有助于全面把握本工作的学术脉络：

  1. 定向渗透与吸收相变经典理论：
    • Hinrichsen, Adv. Phys. 49, 815 (2000)：该综述是非平衡态吸收相变领域的圣经，系统定义了 DP 普适类、临界指数和标度行为，为本文经典自动机的物理解析提供了坚实的统计物理基础。
    • Domany & Kinzel, Phys. Rev. Lett. 53, 311 (1984)：提出了著名的 Domany-Kinzel 自动机，将经典渗透过程与一维统计模型（如 Ising 模型）联系起来，是本工作最直接的物理模型来源。
  2. 等距张量网络状态（isoTNS）：
    • Zaletel & Pollmann, Phys. Rev. Lett. 124, 037201 (2020)：首次系统阐明了二维等距张量网络状态（isoTNS）的规范形式、高效收缩算法和物理制备协议。本论文正是该技术在非平衡相变研究中的一次极其重大的应用拓展。
  3. W 态与母哈密顿量不稳定性：
    • Gioia & Thorngren (2024), arXiv:2310.10716 ； Gioia, Moudgalya, & Motrunich (2025), arXiv:2510.24713：这两个最新工作深入探讨了 $W$ 态无法作为一个稳定相在热力学极限下独立存在的“脆弱性”，并指出其必须依附于特定对称性哈密顿量的退化流形。本论文的研究结果强力印证并丰富了这一理论预言。

4.2 工作的局限性与严厉点评

尽管该工作在统计物理与量子信息交叉领域取得了极具创见性的突破，但站在严谨的学术视角，本工作依然存在以下不可忽视的理论与应用局限性：

  1. 基态简并度的极端脆弱性（Fragility）： 论文中展现出的 W 态到平凡态的量子相变，完全依赖于缺陷数 $N_D$ 的严格 $U(1)$ 守恒。在真实物理系统（如冷原子晶格、超导量子比特或分子量子晶格）中，任何极其微小的局部微扰（例如杂散磁场引起的单点 $Z$ 场微扰，或者 $p_1 \neq 0$ 的非零噪声）都将不可避免地彻底破坏该 $U(1)$ 对称性。这会导致简并的基态流形瞬间解除，真空态 $|vac\rangle$ 成为唯一的平凡基态。因此，这种纠缠相变并不对应于热力学极限下稳定的拓扑相或热力学相，其在实验上的观测需要对系统进行极为苛刻的参数调控。
  2. 母哈密顿量的多体相互作用高昂代价（Frustration-Free Parent Hamiltonian）： 为了支撑这一精细的 isoTNS，推导出的无挫折哈密顿量 $H_{DK}$ 是一个由 8 qubit 投影算符 构成的极度复杂的算符。在实际的量子模拟硬件（如量子芯片或里德堡原子阵列）中，直接实现 8 体相互作用是极其困难的，通常需要引入大量的辅助比特和多步门操作，这大大降低了该量子态作为工程制备纠缠源的可行性。
  3. 非对角可观测量的计算复杂度瓶颈： 虽然对角算符可以通过经典的 MCMC 算法以 $\mathcal{O}(L)$ 的时间复杂度极速收缩并高精度采样，但在计算诸如量子相干性、非局域纠缠熵、拓扑纠缠负度等非对角性质时，该映射无法直接应用经典随机过程。研究者仍不得不面对张量网络全面收缩的指数墙，这限制了该框架在超大体系非平衡动力学全物理量模拟中的应用。

5. 物理化学延伸讨论与未来展望

5.1 量子化学与非平衡态分子动力学的交汇

作为量子化学科研人员，本项工作看似抽象的物理模型，实际上为复杂的分子体系和化学动力学演化提供了极其宝贵的新思路：

  • 激子输运与光合作用人工捕光系统中的相干相变： 在光合作用捕光复合物（如 FMO 复合物）中，激子（Excitons）在叶绿素分子间的输运过程本质上是一个开放量子系统中的非平衡反应-扩散过程。在强耗散限制下，激子的物理行为高度类似于定向渗透（DP）过程。本项工作提供的“经典自动机到量子 TNS”映射，允许化学家将激子在空间分子网络中的经典渗透动力学，完美改写为基于 isoTNS 的量子高效算法，这不仅能高精度捕获分子局部的激子捕获效率，还能量化激子输运通道中非平凡两体纠缠的相干寿命。
  • 催化反应表面相变与多点协同催化： 传统的化学催化剂表面反应（如 CO 在金属表面的氧化过程，即著名的 ZGB 模型）也展现出丰富的吸收相变（催化剂中毒态）。若能利用本文提出的张量映射机制，将催化剂表面的经典动力学态映射为具备特定局部守恒律的量子哈密顿量基态，化学家就可以从纠缠图样的角度，全新理解多点协同催化中的量子纠缠辅助效应，进而利用量子模拟器设计出在临界点处具备代数关联、具有极高协同效率的新一代量子纳米催化材料。

5.2 经典模型对比：为何 Kasteleyn 模型无法实现该相变？

在补充材料中，作者给出了一个绝佳的物理对比体系——Kasteleyn 二聚体模型。在 honeycomb 晶格上，紧密堆积的二聚体（Dimers）可以通过经典玻尔兹曼权重映射为量子态 $|\psi_K(g)\rangle$。该态同样会在临界参数 $g = 1/\sqrt{2}$ 处发生波函数的宏观突变。然而，这两者背后的量子纠缠物理却有着天壤之别：

物理特性Domany-Kinzel 自动机量子态 (本工作)Kasteleyn 二聚体模型量子态
经典对应过程1+1D 随机更新，存在粒子创生、消灭与分支（Branching）经典紧密堆积二聚体，满足严格的单点覆盖硬核约束
守恒律特征仅全局缺陷数 $N_D$ 守恒，子通道允许分支与合并极度严格的弦数 $N_{strings}$ 守恒，不允许任何弦分支
相变性质真正的非平衡临界相变，临界点处全向代数关联参数 $g$ 只是平滑调节各基态权重的参数化流形，无动力学相变
纠缠转变非平凡 W 类长程纠缠 $\rightarrow$ 平凡局域纠缠弦结构被拓扑锁定，无法通过局部变宽实现纠缠拓扑结构改变

这一精细的对比深刻揭示了：粒子分支（Branching）机制与弱守恒律相结合，是促成量子相干纠缠图样相变的关键物理诱因，而过度苛刻的硬核约束（如二聚体模型）反而会锁死系统的纠缠流形，使其失去非平衡临界相变的活性。这为分子量子剪裁（Quantum Molecular Cropping）和高维纠缠分子链的设计提供了极具价值的理论指导。

5.3 未来展望

Boesl 等人的这项前沿工作仅仅掀起了经典非平衡相变与量子张量网络交叉研究的冰山一角。未来的三大核心演进方向非常明确且充满吸引力：

  1. 高维通用非平衡映射的建立：将 2D 键定向渗透映射到 3D 复杂 isoTNS 架构，探索多自由度下的量子临界点。这对于三维高分子晶格中的量子纠缠传输具有直接指导意义。
  2. 量子主动误差校正的全新设计：定向渗透的吸收态天然具备抵御局部噪声的自修复能力。利用本文哈密顿量，研究者可望设计出基于非平衡动力学相变的新型“耗散量子纠错码”，在主动纠错领域展现威力。
  3. 多体疤痕态（Many-body Scars）的系统构建：通过对 $U(1)$ 守恒扇区的精细切割，该模型中的 $W$ 态流形可以被自然地嵌入到高能激发态中，这为探索化学反应中非遍历性热化通道提供了极好的数学建模工具。

本项研究所展现的深邃物理洞察力与优雅的数学框架，必将持续激发凝聚态物理学家与理论化学家在探索希尔伯特空间那片广袤非平衡疆域时的无限灵感。