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

0. 执行摘要 (Executive Summary)

在现代凝聚态物理与量子多体理论中，对称性及其破缺机制是确立物质相与相变分类的基石。传统的朗道-金兹堡-威尔逊（Landau-Ginzburg-Wilson）范式通过局部序参量（Local Order Parameters）的真空期望值来表征自发对称性破缺（SSB）。然而，随着对称保护拓扑相（SPT）以及广义对称性（如非逆对称性，Non-invertible Symmetries）的发现，局域表征方法已显现出本质上的局限。SPT 相无法通过任何局域序参量探测，而必须依赖非局域的弦序参量（String Order Parameters, SOPs）。此外，当系统具备非逆对称性时，其电荷部门（Charge Sectors）不再由群表示论描述，而需诉诸更为宏大的数学工具——融合范畴论（Fusion Category Theory）。

本篇学术博文深度解析了最新研究成果《Algebras of order parameters in one-dimensional spin systems》。该工作利用**矩阵乘积态（MPS）和矩阵乘积密度算符（MPDO）等张量网络形式，在晶格模型中严格证明了弦序参量在一维量子自旋系统中是如何自组织为对称性融合范畴 $\mathcal{C}$ 的德林费尔德中心（Drinfel’d Center） $\mathcal{Z}(\mathcal{C})$ 中的拉格朗日代数（Lagrangian Algebra）的。该研究不仅解决了一直以来被忽视的拉格朗日代数“乘法规则”（Multiplication Rule）在红外（IR）极限下的物理佯谬，还清晰地表明：在自发对称性破缺中，基态所保留的对称性不应被理解为子范畴，而应被理解为融合范畴中的线算符代数（Algebra of Lines）。最后，研究将这一套优美而宏大的范畴学晶格算符表征方案推广到了对称混合态（Symmetric Mixed States）的凝聚代数（Condensable Algebras）**中，为开放量子系统的相分类提供了强有力的解析工具。

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

1.1 核心科学问题

本项研究聚焦于以下几个核心物理与数学问题：

1. 弦序参量的代数本质是什么？ 在给定的自旋链中，具备何种空间变换与对称变换特性的弦算符才能在特定的有间隙相（Gapped Phase）中获得非零的真空期望值？
2. 拉格朗日代数的乘法规则在红外重整化群（RG）流下扮演了什么物理角色？ 过去的研究多关注拉格朗日代数作为对象的分解，而忽略了其代数乘法。如何从格点微观算符的融合中自然诱导出此乘法？
3. 非逆对称性自发破缺的精确数学描述是什么？ 当群对称性 $G$ 被推广为非逆融合范畴 $\mathcal{C}$ 时，自发破缺后的基态空间保留了怎样的“残余对称性”？这些残余对称性在数学上如何互不相同却又在 Morita 等价意义下调和？
4. 混合态下的对称相与序参量如何推广？ 开放系统中的退相干会破坏纯态纠缠。在矩阵乘积密度算符（MPDO）框架下，如何定义强/弱对称性破缺相以及相应的非局域序参量？

1.2 理论基础

1.2.1 融合范畴（Fusion Category）与非逆对称性

在（1+1）维时空中，一个有限自发对称性（含不可逆线缺陷）可以用一个一阶半单、刚性的 $\mathbb{C}$-线性单子范畴 $\mathcal{C}$（即融合范畴）来描述。其简单对象 $X \in \text{Irr}(\mathcal{C})$ 代表拓扑对称线（Topological Symmetry Lines）。当对称性为常规群 $G$ 时，$\mathcal{C} = \text{Vec}_G$（群级联向量空间范畴）；若将其对偶化（如通过规范化/Gauging），则对称线由群表示范畴 $\mathcal{C} = \text{Rep}(G)$ 给出。此时，对称线的融合规则遵循群表示的张量积分解：

其中 $N_{V_1, V_2}^{V_3}$ 为非负整数，指示融合通道数。不可逆性体现在简单对象 $V$ 融合时其对偶 $V^*$ 不满足 $V \otimes V^* = \mathbf{1}$，而是包含多个通道。

1.2.2 德林费尔德中心 $\mathcal{Z}(\mathcal{C})$ 与管代数（Tube Algebra）

在物理上，弦算符的端点可以看作是半无穷对称线缺陷的终止端。这些端点承载着所谓的扭曲区局域算符（Twisted Sector Local Operators）。多重态（Multiplets）的扭曲区局域算符在拓扑缺陷线的“套索（Lasso）”作用下形成代数表示。这一代数即为管代数 $\text{Tube}(\mathcal{C})$。 管代数 $\text{Tube}(\mathcal{C})$ 的模块范畴 $\text{Mod}(\text{Tube}(\mathcal{C}))$ 与融合范畴 $\mathcal{C}$ 的德林费尔德中心 $\mathcal{Z}(\mathcal{C})$ 范畴等价：

$\mathcal{Z}(\mathcal{C})$ 是一个非退化的辫子融合范畴（Braided Fusion Category），它在物理上精确对应于相应的对称拓扑场论（Symmetry TFT, 简称 SymTFT）在体（Bulk）中的线缺陷算符。这就建立了一阶边界自旋链与二阶拓扑体（Bulk）之间的全息对应关系。

 2D 体 (Bulk SymTFT) 
 ----------------------------------------
 |  体缺陷线 (Bulk Line Operators)       |
 |  属于德林费尔德中心 Z(C)             |
 ----------------------------------------
 | (体缺陷线冷凝/Condense)                |
 v                                      v
======================================== (1D 边界自旋链 Gapped Phase)
弦序参量 SOPs (1D Boundary String Operators)
承载端点扭曲区算符，构成 Z(C) 中的拉格朗日代数 L

1.2.3 拉格朗日代数（Lagrangian Algebra）

在辫子融合范畴 $\mathcal{Z}(\mathcal{C})$ 中，拉格朗日代数 $L$ 是一个特殊的连接、可积（Etale）代数，其费罗贝尼乌斯-佩龙维度（Frobenius-Perron Dimension）满足：

拉格朗日代数在物理上描述了体（Bulk）拓扑线缺陷在边界上的冷凝（Condensation）。冷凝在边界上的算符在有间隙相的基态中获得非零真空期望值（VEV），从而成为该相的弦序参量；而那些不属于拉格朗日代数的算符则在边界上被限制（Confined），其期望值严格为零。

1.3 技术难点

将上述高度抽象的范畴学结构与晶格上的具体多体波函数及算符联系起来，面临几项严苛的技术难点：

1. 非局域弦算符的晶格表征：如何在晶格 Hilbert 空间中定义扭曲区局域算符？非逆对称算符表现为矩阵乘积算符（MPO），而非传统的张量积。如何在其端点拼接局域算符而不破坏融合范畴的 $F$-符号和乘法一致性？
2. 红外极限下的投影与代数乘法提取：在紫外（UV）格点上定义的弦算符往往非常繁杂。如何将其投射到重整化群（RG）不动点，并从中无损地提取出拉格朗日代数的乘法规则？这要求精确对角化包含缺陷线作用的混合转移矩阵（Mixed Transfer Matrix）。
3. 非逆对称自发破缺中的 Morita 等价性：非逆对称性破缺时，不同的基态（通常对应更高的多重度表示）所保留的自发破缺残余对称线并不相同（即它们对应的局域代数不同）。如何证明这些不同的代数在范畴学上是 Morita 等价的？

1.4 方法细节：张量网络形式与模范畴分解

为了克服上述难点，论文采用了**块单射矩阵乘积态（Block-Injective MPS）**作为有间隙基态的普适表示。

1.4.1 微观 Hilbert 空间与 MPO 对称性

考虑一维量子链，其物理格点承载自发对称性 $\mathcal{C} = \text{Vec}_G$。我们选取一个表示 $V \in \text{Rep}(G)$ 作为局域物理空间。拓扑缺陷线 $C_g \in \text{Vec}_G$ 在 Hilbert 空间上的作用由 MPO 给出。若对称性推广为非逆表示范畴 $\mathcal{C} = \text{Rep}(G)$，物理空间则由有限 $G$- graded 向量空间 $K$ 承载。对称性 MPO 此时通过群的逆作用来实现（如公式 4.1 所示）：

该 MPO 的局域张量满足 Clebsch-Gordan 系数融合（公式 4.5）和包含 $F$-符号的五边形一致性关系。

1.4.2 块单射 MPS 与基态分类

任意对称的有间隙基态总能表示为如下块对角形式的 MPS（公式 2.10）：

每一块 $\theta_a$（其虚拟空间为 $M_a$）在热力学极限下是单射的（Injective）且彼此正交。对称线 MPO 作用在基态上时，会置换这些块（公式 2.19）：

这表明简单基态块之间的置换构成一个左 $\mathcal{C}$-模范畴 $\mathcal{M}$。简单对象 $C_g \in \mathcal{C}$ 作用于简单模块 $M_a \in \mathcal{M}$ 产生新的基态块 $M_{g(a)}$，并由局域张量 $\phi^{g, a}$ 和 $\bar{\phi}^{g, a}$ 连接。这些连接算符在连续融合两次时，会产生模范畴结合子（Module Associator）$\alpha^{\triangleright}(g_1, g_2)(a)$（公式 2.24）：

这一结合子 $\alpha^{\triangleright}$ 满足五边形公理，标志着晶格 MPS 直接决定了半单模范畴 $\mathcal{M}$ 的核心数学输入。对于对称性为 $\text{Vec}_G$ 的情况，模范畴由子群 $H \le G$ 及其 2-上同调类 $\beta \in H^2(H, U(1))$ 完全分类，记作 $\mathcal{M}(H, \beta)$。

1.4.3 弦序参量（SOP）与转移矩阵分析

我们将弦算符定义为：在自旋链的 $i$ 和 $j$ 位点之间施加截断的对称性缺陷线 $C_x$，其端点分别附着局域扭曲区算符 $c_1$ 和 $c_2$（见公式 3.21）：

计算该算符在基态 $| \Psi(\theta_a) \rangle$ 下的期望值，本质上转化为在 $i$ 到 $j$ 之间插入缺陷转移矩阵（Defect Transfer Matrix）$\mathbb{E}^x_{\theta_a}$。在热力学极限 $|j - i| \to \infty$ 下，期望值收敛到最大特征值对应的左、右不动点的收缩（Contraction）：

为了使该式不为零，缺陷线 $C_x$ 必须包含在基态保留的局域对称性代数中，即 $x \in r_a H r_a^{-1}$。进一步地，通过精细分析管代数 $\text{Tube}(\mathcal{C})$ 作用在端点局域算符上的表现，论文证明：能够获得非零期望值的端点局域算符必须构成 $\text{Tube}(\mathcal{C})$ 模 $L(H, \beta)$ 的直接加数。 此模的辫子化（Half-braiding）正是由模范畴的结合子 $\alpha^{\triangleright}$ 所赋予的，这直接确立了 SOP 多重态与体拉格朗日代数 $L$ 在对象层面的一一对应。

2. 关键 Benchmark 体系、计算数据与物理案例解析

该研究提供了一套极其完备的解析与微观算符 Benchmark 体系。我们针对两个极具代表性的群自旋链体系进行微观重构，定量展示代数结构是如何物理显现的。

2.1 体系一：可逆克莱因四群对称性自旋链 $\mathbb{Z}_2 \oplus \mathbb{Z}_2$

2.1.1 物理背景与微观态

克莱因四群 $G = \mathbb{Z}_2 \oplus \mathbb{Z}_2 = \{ (), (12)(34), (13)(24), (14)(23) \}$ 是描述一维 Haldane 相和 AKLT 模型（自旋-1 链）的经典对称群。当边界对称性为 $\text{Vec}_{\mathbb{Z}_2 \oplus \mathbb{Z}_2}$ 时，存在非平凡的二阶上同调类 $H^2(\mathbb{Z}_2 \oplus \mathbb{Z}_2, U(1)) \cong \mathbb{Z}_2$。其非平凡元素 $\beta$ 对应非平凡的 SPT 相（即簇态/Cluster State 或 Haldane 拓扑相）。

2.1.2 拉格朗日代数分析

在该体系中，存在 6 个不同的有间隙相，分别对应 6 个拉格朗日代数。我们对非平凡 SPT 相 $L(\mathbb{Z}_2 \oplus \mathbb{Z}_2, \beta)$ 进行计算。选取 $\beta$ 分支，其在一维自旋-1/2 双量子比特表示下的 MPO 对称线 $C_{(12)(34)}$、$C_{(13)(24)}$ 和 $C_{(14)(23)}$ 分别作用为：

利用公式 (3.60)，拉格朗日代数 $L(\mathbb{Z}_2 \oplus \mathbb{Z}_2, \beta)$ 分解为德林费尔德中心 $\mathcal{Z}(\text{Vec}_{\mathbb{Z}_2 \oplus \mathbb{Z}_2})$ 中的简单对象和特征：

这一特征分解直接告诉我们，在红外极限下，该相中只有满足上述荷（Charge）配对的弦算符才能凝聚。具体地，其端点算符根据公式 (3.61) 作用在双量子比特物理空间上呈现为：

由此我们微观重构出在自旋-1/2 链（以双量子比特表示单格点）中对应的弦序参量（公式 3.62）：

2.1.3 自旋-1 经典 Haldane 弦序的复现

当我们将该微观局域表示通过等变同构映射到自旋-1 表示（AKLT 态）时，利用拓扑线在物理态上的对角作用：

将算符映射代入后，完美复现了 Nijs 和 Rommelse 经典的自旋-1 链弦序参量（公式 3.63）：

这证明了经典的物理弦序参量其本质上是由拉格朗日代数的特征荷结构决定的。

2.2 体系二：交错群对称性自旋链 $A_4$

2.2.1 物理背景与对称破缺

交错群 $A_4$ 阶数为 12，是典型的非阿贝尔群。它有四个共轭类：$[()]$、$[(12)(34)] = \{(12)(34), (13)(24), (14)(23)\}$、$[123]$、$[132]$。其不可约表示包含三个一维表示 $\mathbb{C}_1, \mathbb{C}_\omega, \mathbb{C}_{\bar{\omega}}$ 和一个三维不可约表示 $V_3$。

2.2.2 7种有间隙相与拉格朗日代数分解表

通过对 $A_4$ 子群及其 2-上同调类进行范畴学分类，我们可以确定系统具备 7 种对称有间隙相。论文中首次完整解析写出了这 7 个拉格朗日代数的德林费尔德中心对象分解，具体计算结果归纳为下表：

2.2.3 物理数据解析：利用 3 个序参量区分所有 7 个相

这是一个极具洞察力的物理应用。传统的表征方法需要极其繁多的多点关联函数，而通过上述拉格朗日代数的交叉交并关系，我们仅需测量 3 个特定序参量的期望值，即可将所有 7 个相完全区分开：

1. 局域序参量 $O_1$：属于表示 $([()], \mathbb{C}_\omega)$。其非零意味着系统的 $\mathbb{C}[\mathbb{Z}_3]$ 这一对称线缺陷代数被破坏。
2. 局域序参量 $O_2$：属于表示 $([()], V_3)$。其非零探测了整个 $\mathbb{C}[\mathbb{Z}_2 \oplus \mathbb{Z}_2]$ 缺陷线代数自发破缺。
3. 非局域弦序参量 $S_3$：属于扭曲区简单对象 $([(12)(34)], \mathbb{C}_{1,1})$。该弦算符的非零期望值表征了拓扑凝聚线在相中的凝聚特性。

通过计算这三个算符在各相基态下的期望值，可得到如下特征矩阵（1 代表非零期望值，0 代表期望值为零）：

注：$\mathbb{Z}_2$ 相与 $\mathbb{Z}_2 \oplus \mathbb{Z}_2$ 相虽然在上述三指标下相同，但可以通过进一步分析其基态简并度（分别为 6 和 3）或补充测量伴随的共轭类算符轻易拆分。这展示了拉格朗日代数分类法在相探测设计上的无与伦比的穿透力。

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

为了使量子化学和量子多体物理科研人员能够快速复现论文中的核心结论，本节提供了一套基于 Julia 语言和 Python 语言的多体数值模拟指南，用于构建对称 MPO、求解 MPS 基态以及数值提取弦序参量。

3.1 算法流程设计

弦序参量的数值计算和拉格朗日代数提取可以总结为以下七步算法流程：

 算法步骤:
 [Step 1] 定义对称群 G/范畴 C 及其局域表示 V 
   |--> 构造局域物理 Hilbert 空间 
 [Step 2] 构造对称性 MPO 算符 
   |--> 用于对系统施加全局或局部拓扑线缺陷作用 
 [Step 3] 运行 DMRG/iTEBD 算法 
   |--> 求解目标 Hamiltonian 得到基态的 block-injective MPS 
 [Step 4] 提取转移矩阵 E 的最大左/右特征向量 
   |--> \Xi_l 与 \Xi_r，作为不动点投影算符 
 [Step 5] 在特定格点截断对称性 MPO 
   |--> 在端点 i 和 j 拼接扭曲区局域算符 O(w_c) 
 [Step 6] 数值收缩整个张量网络 
   |--> 计算在投影不动点下的 SOP 期望值 
 [Step 7] 改变端点算符表示，扫描整个 
    --> 验证非零通道与拉格朗日代数 L 的分解完全一致

3.2 Julia 代码框架（基于 TensorKit.jl）

TensorKit.jl 是处理具备任意对称范畴结构张量网络的顶级 Julia 库。以下是复现可逆群自旋系统 SOP 凝聚的核心代码架构。

using TensorKit
using LinearAlgebra

# =====================================================================
# Step 1: 定义对称群及物理空间 (以 Z2 对称性自旋-1/2 链为例)
# =====================================================================
G = ProductGroup{Tuple{Z2Irrep}}()
# 物理格点承载平移与 Z2 对称性，虚拟空间维数为 D
物理空间 = RepresentationSpace{Z2Irrep}(0 => 1, 1 => 1) # 0为偶，1为奇
虚拟空间 = RepresentationSpace{Z2Irrep}(0 => 4, 1 => 4) # 键维数 D=8

# =====================================================================
# Step 2: 构造块单射 MPS 局域张量 A
# =====================================================================
A = TensorMap(rand, ComplexF64, 虚拟空间 * 物理空间, 虚拟空间)
# 对张量进行左正则化 (Left Orthonormalization)
A, R = leftorth(A)

# =====================================================================
# Step 3: 构建对称性转移矩阵 E (带有 Z2 缺陷作用)
# =====================================================================
# Z2 缺陷算符由简单对象 (1) 的一维表示作用实现
Z2_缺陷 = TensorMap(zeros, ComplexF64, 物理空间, 物理空间)
blocks(Z2_缺陷)[Z2Irrep(0)] .= 1.0
blocks(Z2_缺陷)[Z2Irrep(1)] .= -1.0 # 奇部门赋予负相因子

# 构造包含缺陷的转移矩阵作用子
function apply_transfer_matrix(v_in, A, Z2_缺陷)
    # v_in 是虚拟空间的左边界向量
    # 使用 TensorOperations 的爱因斯坦求和约定进行高效收缩
    @tensor v_out[left_out, right_out] := v_in[left_in, right_in] * 
                                         A[left_in, phys, center] * 
                                         Z2_缺陷[phys, phys_prime] * 
                                         conj(A[right_in, phys_prime, center_prime])
    return v_out
end

# =====================================================================
# Step 4: 幂指数法对角化转移矩阵，提取不动点投影算符
# =====================================================================
Ξ = Tensor(rand, ComplexF64, 虚拟空间 * 虚拟空间')
for iter in 1:100
    Ξ_new = apply_transfer_matrix(Ξ, A, Z2_缺陷)
    Ξ_new /= norm(Ξ_new)
    global Ξ = Ξ_new
end
println("转移矩阵最大特征值收敛，不动点张量范数: ", norm(Ξ))

# =====================================================================
# Step 5: 计算弦序参量 (SOP) 期望值
# =====================================================================
# 端点局域算符 O(w_c) 选取为带奇宇称的自旋算符
O_endpoint = TensorMap(rand, ComplexF64, 物理空间, 物理空间)
# 计算 10 个格点间距的 SOP
L = 10
expectation_val = 1.0
# 收缩端点和体缺陷
# (此处为演示性张量计算，实际中对各格点依次作用 apply_transfer_matrix 即可)

3.3 开源工具包与代码库推荐

对于研究非逆对称性与高维自旋链的学者，以下几个开源 Repo 是不可或缺的利器：

1. TensorKit.jl (Julia):
   • Link: https://github.com/Jutho/TensorKit.jl
   • 推荐理由: 能够完美原生支持任意辫子融合范畴（包括任意非阿贝尔群表示与超对称）的张量网络工具包。
2. MPSKit.jl (Julia):
   • Link: https://github.com/Jutho/MPSKit.jl
   • 推荐理由: 专为一维无限和有限 MPS 算法设计的上层工具箱，内置一维有间隙相基态求解的主流算法。
3. TeNPy (Tensor Network Python):
   • Link: https://github.com/tenpy/tenpy
   • 推荐理由: Python 生态下最流行、成熟的张量网络库。非常适合快速搭建自旋链、运行 DMRG 算法，并能方便地手动截断对称流计算非局域弦关联。
4. YASTN (Yet Another Symmetric Tensor Network library):
   • Link: https://github.com/yastn/yastn
   • 推荐理由: 专为包含任意阿贝尔对称性张量网络设计的 Python/PyTorch 库，支持利用 GPU 加速大键维数转移矩阵的求解。

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

4.1 关键参考文献

1. [dNR89] den Nijs, M. & Rommelse, K. Preroughening transitions in crystal surfaces and valence-bond phases in quantum spin chains. Phys. Rev. B 40 (1989) 4709.
   • 贡献: 首次引入非局域弦序参量，成功探测了自旋-1 链中的隐藏对称性破缺（Haldane 相）。
2. [CGW11] Chen, X., Gu, Z.-C. & Wen, X.-G. Classification of gapped symmetric phases in one-dimensional spin systems. Phys. Rev. B 83 (2011) 035107.
   • 贡献: 利用 MPS 的虚拟空间表示理论，对 1D 经典可逆对称性下的有间隙相进行了上同调分类。
3. [DNO13] Davydov, A., Nikshych, D. & Ostrik, V. On the structure of the Witt group of braided fusion categories. Selecta Mathematica 19 (2013) 237.
   • 贡献: 严格定义并分类了辫子融合范畴中的拉格朗日代数，奠定了 SymTFT 边界冷凝的数学根基。
4. [LOST22] Lin, Y.-H., Okada, M., Seifnashri, S. & Tachikawa, Y. Asymptotic density of states in 2d CFTs with non-invertible symmetries. JHEP 03 (2023) 094.
   • 贡献: 阐明了管代数（Tube Algebra）与扭曲区算符多重态在一维晶格系统中的微观对应。
5. [BBPSN23] Bhardwaj, L., Bottini, L. E., Pajer, D. & Schäfer-Nameki, S. Gapped phases with non-invertible symmetries: (1+1)d. SciPost Phys. 18 (2025) 032.
   • 贡献: 首次系统提出非逆对称性下有间隙相可以通过拉格朗日代数分类的理论假说。

4.2 局限性评论

虽然该论文在数学表征和一维格点复现上取得了里程碑式的进展，但作为面向量子化学与多体物理的研究者，我们必须指出其在实际应用和物理机制上的几点重要局限：

4.2.1 强烈的红外重整化群（RG）不动点依赖性

论文在证明“拉格朗日代数乘法规则”与“端点算符融合规则”的一致性时，关键性地使用了重整化群不动点近似（即转移矩阵 $\mathbb{E}$ 严格为投影算符，相关长度 $\xi = 0$）。然而，真实的自旋链（特别是具有长程相互作用或接近临界点的相，如量子化学中的多活性空间体系）偏离 RG 不动点极远，具有显著的非零相关长度 $\xi > 0$。在此情况下，非零期望值的 SOP 无法呈现为干净的常数，而是带有指数衰减项：

此时，代数乘法结构的数值提取会受到背景噪声和重正化因子的强烈干扰，使得微观算符乘法在非不动点处的物理定义变得模糊不清。

4.2.2 高维（2D 及以上）推广的数学瓶颈

一维系统的对称性对应 2D 体的线缺陷，由融合一阶范畴描述。若将该框架推广至二维自旋系统，对称性缺陷演化为 2D 空间中的面缺陷（Surface Defects），对称范畴必须升级为融合二阶范畴（Fusion 2-Categories）。其对应的德林费尔德中心和拉格朗日代数其数学结构呈指数级复杂化，目前在数学上尚未完全分类。晶格上的张量网络也必须从 MPS 升级为PEPS（投影纠缠配对态），其收缩（Contraction）本身就是 NP-hard 的。因此，该工作所依托的精细解析推导和数值验证在 2D/3D 体系中将遭遇极大的计算与数学瓶颈。

4.2.3 非逆对称 MPO 的构建与计算复杂度

对于较复杂的非逆对称性范畴（例如 Fibonacci 范畴或 Haagerup 范畴），其融合通道极多，对应的晶格 MPO 键维数（Bond Dimension）和物理张量极其复杂。在数值计算中，存储这些高度非阿贝尔的 Clebsch-Gordan 系数和 $F$-符号需要消耗巨量的内存。此外，量子多体算法（如 DMRG）中收缩此类 MPO 的计算复杂度极高。这导致该方法目前在很大程度上局限于解析推导以及极小群对偶（如 $\text{Rep}(A_4)$）的数值展示，距离解决分子体系和实际关联过渡金属氧化物中的非逆对称性还有很长的路要走。

5. 其他必要的补充：全息 SymTFT 视角与线代数自发破缺

为了给读者提供更加宏大、完整的理论图景，本节补充两个深度物理概念：对称拓扑场论（SymTFT）的全息夹心构造，以及为什么非逆对称性破缺中的基态保留对称性必须由“线代数”而非“子范畴”刻画。

5.1 全息夹心构造 (The Sandwich Construction)

在现代高能物理与凝聚态物理的交叉领域，Symmetry TFT (SymTFT) 提供了理解广义对称性的最强有力全息视角。考虑一个一维物理系统（边界自旋链），其对称性为 $\mathcal{C}$。我们可以在高一维（2D 空间）中构造一个无能隙的拓扑秩序（Topological Order），该拓扑秩序由德林费尔德中心 $\mathcal{Z}(\mathcal{C})$ 唯一确定。这个 2D 体（Bulk）被称为 SymTFT。

 2D Bulk (SymTFT: Topological Order Z(C)) 
=========================================
|                                       |
|               体线算符                |
|              (Bulk Lines)             |
|                                       |
=========================================
|                   |                   |
| (冷凝)             | (被禁闭)           |
| Condense          | Confined          |
|                   v                   v
-----------------------------------------
 1D Boundary (拉格朗日代数 L 决定边界相)

SymTFT 包含两个边界。一侧边界是对称边界（Symmetric Boundary），它是不动的，编码了整个对称性范畴 $\mathcal{C}$。另一侧边界是物理边界（Physical Boundary），它可以通过选择不同的拓扑边界条件（Topological Boundary Conditions）来模拟一维自旋链所处的所有可能对称相。这些拓扑边界条件在数学上与 $\mathcal{Z}(\mathcal{C})$ 中的拉格朗日代数 $L$ 一一对应。

• 弦算符的凝聚与禁闭： 当我们在 1D 物理边界上拉出一条弦序参量（SOP）时，它在 2D 体中对应于一条体缺陷线的一端终止于物理边界。如果这条体缺陷线属于拉格朗日代数 $L$，意味着该线可以在边界上自流冷凝（Condense），从而在边界上终止时不需要任何额外的能量代价。对应的 1D SOP 因而具有长程非零期望值。反之，如果体缺陷线不属于 $L$，它在边界上就会被“禁闭”（Confined），边界算符的期望值随距离呈指数衰减至零。因此，有间隙相中弦序参量的选择与计算，本质上是对全息体中哪些拓扑线可以在物理边界上冷凝的测量。

5.2 突破传统：为什么对称自发破缺属于线代数？

这是一个极其重要的概念革新。传统的凝聚态物理思维认为，当一个系统自发破缺一部分对称性时，其基态保留的对称性是一个子群（Subgroup）。例如，$G$ 破缺到 $H$。当推广到融合范畴对称性 $\mathcal{C}$ 时，直觉上会认为基态保留的对称性应该是一个融合子范畴（Fusion Subcategory）。

本工作用严谨的数学推导彻底打破了这一直觉！

以非逆对称性 $\text{Rep}(G)$ 自发破缺为例（其中 $G$ 为非阿贝尔群）。当我们考虑一维基态空间时，不同的简并基态通常对应于表示 $G$ 的不同不可约表示多重态。例如，$A_4$ 的三维不可约表示基态 $|\Psi(V_3)\rangle$。计算表明：

1. 基态 $|\Psi(V_3)\rangle$ 所保留（即作用在基态上不产生正交态）的拓扑缺陷线在 $\text{Rep}(A_4)$ 中是 $\mathcal{A} = \text{End}_{\mathbb{C}}(V_3)$。
2. 这一代数 $\mathcal{A}$ 作为范畴中的对象，其分解为： $$\mathcal{A} \cong \mathbb{C}_1 \oplus \mathbb{C}_\omega \oplus \mathbb{C}_{\bar{\omega}} \oplus 2 \cdot V_3$$
3. 显然，由于包含三维不可约表示 $V_3$ 的两次直接加数，这一简单对象的集合根本无法在融合范畴的张量积下封闭！也就是说，它们不能构成一个合法的子范畴。

这一结论揭示了非逆对称性自发破缺的核心物理本质：对称性自发破缺中，基态所保留的对称性不应在范畴级别（Category Level）上理解，而应在代数级别（Algebra Level）上理解。基态保留的是对称线构成的代数（Algebra of Lines），具体为简单对象的内部自同构代数（Internal Hom）。 虽然不同的简并基态保留的线代数并不相同（如同表 2.2.2 所示），但这些代数在数学上是 **Morita 等价（Morita Equivalent）**的。这种 Morita 等价性保证了它们指向同一个独一无二的有间隙物理相。这一深刻的物理洞察是该论文做出的最杰出的学术贡献之一，它将群的对称性破缺理论推向了范畴化代数的终极高度。