来源论文: https://arxiv.org/abs/2607.00762v2 生成时间: Jul 04, 2026 00:04

跨越禁区：量子蒙特卡洛揭示反铁磁绝缘体与d波超导体的解禁闭量子临界点（DQCP）深度解析

0. 执行摘要

在强关联电子体系和高温超导物理的研究中，反铁磁（AFM）Néel绝缘体与$d$波超导体（dSC）之间的相变行为一直是最核心的谜题之一。经典的Landau-Ginzburg-Wilson（LGW）对称破缺相变范式认为，两个具有不同对称性破缺的物相（反铁磁打破了$O(3)$自旋旋转对称性，超导打破了$U(1)$电磁规范对称性）之间的直接相变必须是强一阶的，除非通过极其精细的参数微调。然而，在铜氧化物（Cuprates）高温超导体中，实验观测表明随着空穴掺杂，系统在极低掺杂浓度下即经历了一个极其陡峭且直接的从反铁磁绝缘体到$d$波超导态的转变，几乎不存在显著的共存区。

为了解释这一“非LGW”的物理现象，理论物理学家提出了**解禁闭量子临界点（Deconfined Quantum Critical Point, DQCP）**的概念。在DQCP处，电子发生了电荷与自旋的分数化（Fractionalization），分裂为无质量的狄拉克自旋子（Spinon）和玻色型电荷子（Chargon），它们通过强耦合的涌现规范场相互作用。由于两边的相变序参量都是由这些分数化“部分子”（Partons）的复合算符所构成，因此可以实现两个非同源对称性破缺物相之间的二阶连续相变。

然而，由于强关联多体系统中的量子波动，特别是量子蒙特卡洛（QMC）在处理此类体系时必然遭遇的**“负符号问题”（Negative Sign Problem），这一非平庸的理论长期以来缺乏精确的、无近似的数值证据。近期，由陈闯（Chuang Chen）、Subir Sachdev以及孟子阳（Zi Yang Meng）组成的联合研究团队，在半满（Half-filling）正方晶格上设计了一种巧妙的电子分数化机制，利用自旋子-电荷子部分子表象**构建了一个自旋子和电荷子均在$\pi$-flux背景下运动并强耦合于动力学活跃的$SU(2)$晶格规范场的晶格模型。最关键的是，该模型在半满且具备粒子-空穴对称性的条件下，彻底规避了QMC的负符号问题。

通过大规模的行列式量子蒙特卡洛（DQMC）数值模拟，研究团队首次在热力学极限下给出了反铁磁Néel绝缘体到节点$d$波超导体（nodal d-wave superconductor）之间存在二阶连续解禁闭量子临界点的直接、铁证般的数值证据：

1. 相变连续性与临界点锁定：通过对反铁磁结构因子和超导配对结构因子的有限尺寸标度（Finite-Size Scaling）以及相应的关联比（Correlation Ratio）的系统测量，精确定位了量子临界点$\kappa_\tau^c = 2.2(2)$，并证实两者的序参量在相变点同时、连续地消失，同时伴随着非平庸的价键固体（VBS）通道中涨落的剧烈增强（呈现SO(5)涌现高度对称性特征）。
2. 分数化电子谱函数的直接谱学证据：利用随机解析延拓（Stochastic Analytic Continuation, SAC）技术，首次计算了规范不变的单粒子电子谱函数 $A(\mathbf{k}, \omega)$。结果表明，在$d$波超导相内，系统呈现完美的、无能隙的节点狄拉克色散（nodal Dirac dispersion），而在跨越临界点进入反铁磁相后，该色散连续地转变为全能隙的激发，完美重现了高低温超导谱学中准粒子行为的演化规律。

本篇博客将从物理图像、数学公式推导、DQMC算法细节、数值Benchmark结果以及局限性等维度，对该工作进行全景式的深度解析。

1. 核心科学问题、理论基础与方法细节

1.1 核心科学问题：超越Landau-Ginzburg-Wilson佯谬

凝聚态物理的核心基石之一是Landau提出的物相与相变理论。该理论基于局部序参量（Local Order Parameter）和自发对称性破缺（Spontaneous Symmetry Breaking）。按照LGW范式：

• 反铁磁Néel态由实空间自旋交错磁矩（staggered magnetization）$\mathbf{n} \sim (-1)^i \mathbf{S}_i$ 表征，其对称性破缺为 $SU(2)_\text{spin} \to U(1)$（或 $O(3) \to O(2)$）。
• $d$波超导态（dSC）由超导配对序参量 $\Delta_{ij} = \langle c_{i\alpha} \epsilon_{\alpha\beta} c_{j\beta} \rangle$ 表征，其打破了全局的 $U(1)_\text{electromagnetic}$ 电磁规范对称性。

在传统的LGW机制下，由于两个序参量属于完全不同的对称性群，它们的涨落是相互独立的。如果系统通过改变某个控制参数从AFM相变到dSC相，在相图上只能出现两种可能：

1. 一阶不连续相变（序参量突变）。
2. 共存相或中间无序相（两个序参量在一个区间内共存，或者中间隔着一个既无反铁磁又无超导的对称性恢复相，经历两次二阶相变）。

然而，实验中强关联电子系统（如空穴掺杂铜氧化物）的相图却极度逼近直接的连续转变。为了在不需要精细微调参数的前提下，实现AFM与dSC之间的直接二阶连续相变，必须引入解禁闭量子临界点（DQCP）。在DQCP处，凝聚态系统中的“基本粒子”——电子，由于强关联效应不再是不可分割的整体，而是发生分数化，分裂为携带自旋1/2的费米型自旋子（spinon）和携带电荷$e$的玻色型电荷子（chargon）。这两个部分子通过一个涌现出来的规范场（在本文中为$SU(2)$ gauge field）强耦合。AFM序和dSC序在物理上不过是自旋子和电荷子在不同物理真空下的凝聚或配对表现。因此，在规范场的主导下，两个序参量被内在锁定，从而能够实现直接的连续相变。

1.2 理论基础：部分子（Parton）表象与$SU(2)$规范场

为了在晶格模型上具体实现这一构想，研究团队采用了**自旋-电荷分离（Spin-Charge Separation）**的部分子构造。在正方晶格的每个格点 $i$ 上，电子算符 $c_{i\alpha}$（$\alpha = \uparrow, \downarrow$）被写为自旋子 $f_{i\alpha}$ 和玻色型电荷子（由二维复玻色子场 $B_i = (B_{1i}, B_{2i})^T$ 表示）的复合体。其具体的代数表示（满足规范不变性）为：

这里，$\epsilon = i\sigma_y$ 是反对称矩阵。这种表示引入了一个局域的 $SU(2)$ 规范自由度。为了保持物理希尔伯特空间的维数不变，必须满足局域约束。具体而言，自旋子 $f_i$ 和电荷子 $B_i$ 均在基本的 $SU(2)$ 规范场下变换：

规范变换作用为：$\Psi_i \to V_i \Psi_i$, $B_i \to V_i B_i$，其中 $V_i \in SU(2)$。显然，物理电子算符 (1.1) 在这一规范变换下是完全不变的（即规范不变量），因此代表了物理上可观测的带电激发。这一部分子表象允许自旋子携带自旋 $S=1/2$ 但不带电荷，电荷子携带电荷 $+e$ 但自旋为 $0$。它们之间唯一的纽带是强相互作用的 $SU(2)$ 晶格规范场 $U_{ij}$。

1.3 技术难点：量子蒙特卡洛中的“负符号问题”及其破局

对于大部分强关联费米子系统或具有几何挫败的磁体，在进行路径积分量子蒙特卡洛模拟时，玻色化后的行列式权重（fermi determinant weight）往往会出现负数甚至复数。这被称为负符号问题。一旦出现负符号，蒙特卡洛采样的相对误差将随系统尺寸 $L$ 和投影虚时 $\beta$ 呈指数级放大：

这使得大尺寸的极限计算变得完全不可能。为了克服这一限制，该工作做出了最关键的对称性设计。研究团队将自旋子和电荷子都置于 $\pi$-flux（$\pi$磁通）的背景中运动。由于两者感受到的 $\pi$-flux 相互抵消，物理电子在通过闭合路径时感受到的净磁通依然为零，从而维持了真实的晶格物理不变。

更妙的是，在半满（Half-filling）以及满足粒子-空穴对称性的条件下，通过引入反幺正时间反演对称性（Antiunitary Time-Reversal Symmetry），可以将自旋子行列式严格锁定为非负实数。具体步骤是对自旋向下（spin-down）的自旋子施加一个局域的粒子-空穴变换：

在此新基（$\tilde{f}$ 基）下，哈密顿量中所有涉及 $f^\dagger f^\dagger$ 和 $ff$ 的配对项都被转换为常规的二次型跳跃项 $\tilde{f}^\dagger \tilde{f}$。经过这一变换，自旋子部分的有效费米作用量可以用一个一维矢量 $\vec{\tilde{f}} = (\tilde{f}_{1\uparrow}, \tilde{f}_{1\downarrow}, \tilde{f}_{2\uparrow}, \tilde{f}_{2\downarrow}, \dots)^T$ 表示。相对应的跳跃矩阵 $K[U]$ 满足在如下定义的时间反演变换下不变：

其中 $\mathcal{K}$ 是复共轭算符，$\sigma_y$ 作用在自旋空间上。由于 $\mathcal{T}^2 = -1$，根据克拉默定理（Kramers Theorem），该二次型费米矩阵的所有本征值必定成对出现（共轭对）。因此，其费米决定子满足：

对于任意的 $SU(2)$ 规范场构型 $U_{ij}$，该决定子恒大于等于零。这在数学上严谨地证明了该物理模型完全不存在负符号问题。这也是国际上首次在保留完整 $SU(2)$ 规范自由度的同时，实现无符号问题的自旋子-电荷子强耦合系统的精确数值模拟。

1.4 晶格模型与拉格朗日量的严谨代数结构

系统整体的晶格拉格朗日量写为：

其中 $\mathcal{L}_s$ 是自旋子（spinon）与规范场耦合的部分，而 $\mathcal{L}_c$ 是电荷子（chargon）部分。具体展开如下：

1.4.1 自旋子与SU(2)规范拉格朗日量

• 自旋子动力学项（空间项）：

  $$\mathcal{L}_f = iJ \sum_{\langle i,j \rangle} e_{ij} \left( \Psi_i^\dagger U_{ij} \Psi_j - \Psi_j^\dagger U_{ji} \Psi_i \right) \quad (1.9)$$

  其中，$J=1$ 为能量单位。$e_{ij}$ 代表反平行的 $\pi$-flux 背景，满足 $e_{ij} = -e_{ji}$。在正方晶格上，我们选择：$e_{i, i+\hat{x}} = 1$，$e_{i, i+\hat{y}} = (-1)^{x_i}$。这一选择保证了每个基本格点（plaquette）包含一个 $\pi$ 规范磁通。$U_{ij} \in SU(2)$ 是连接格点 $i$ 与邻近格点 $j$ 的动力学规范场链（link variables）。

• 规范场空间Yang-Mills项：

  $$\mathcal{L}_U = \kappa \sum_{\Box} \left( 1 - \frac{1}{2} \text{ReTr} \prod_{ij \in \Box} U_{ij} \right) \quad (1.10)$$

  该项惩罚空间格片（plaquette）上的规范场涨落，$\kappa$ 为空间规范场耦合常数。

• 规范场虚时Yang-Mills项：

  $$\mathcal{L}_{U\tau} = \frac{1}{\kappa_\tau \Delta \tau^2} \sum_{ij, \tau} \left( 1 - \frac{1}{2} \text{ReTr} \left[ U_{ij}^\dagger(\tau) U_{ij}(\tau+1) \right] \right) \quad (1.11)$$

  $\kappa_\tau$ 控制时空方向上的规范波动，是微调系统由自旋解禁闭相进入 confinement 相的重要物理控制参数。$\Delta \tau = 0.1$ 为虚时步长。

1.4.2 电荷子拉格朗日量

• 电荷子二次型耦合项：

  $$\mathcal{E}_2(U, B) = (r + 2\sqrt{2}w) \sum_{i} B_i^\dagger B_i + iw \sum_{\langle ij \rangle} e_{ij} \left( B_i^\dagger U_{ij} B_j - B_j^\dagger U_{ji} B_i \right) \quad (1.13)$$

  其中 $r$ 是调控电荷子质量的重要参数。当 $r > 0$ 时，电荷子是有能隙的（Mott 绝缘体）；当 $r < 0$ 时，电荷子凝聚，从而诱导物理电子配对，系统进入超导态。$w$ 为电荷子的跳跃振幅，同样是在 $\pi$-flux 背景中运动。

• 电荷子四次自相互作用项（保持规范不变性与电荷守恒）：

  $$\mathcal{E}_4(U, B) = \frac{u}{2} \sum_{i} \rho_i^2 + V_1 \sum_{i} \rho_i (\rho_{i+\hat{x}} + \rho_{i+\hat{y}}) + g \sum_{\langle ij \rangle} |\Delta_{ij}|^2 + J_1 \sum_{\langle ij \rangle} Q_{ij}^2 + K_1 \sum_{\langle ij \rangle} J_{ij}^2 \quad (1.14)$$

  这里，定义了四种规范不变的局域物理量（复合算符）：

  1. 格点电荷密度：$\rho_i = B_i^\dagger B_i$
  2. 键能密度：$Q_{ij} = \text{Im}(B_i^\dagger e_{ij} U_{ij} B_j)$
  3. 键流密度：$J_{ij} = \text{Re}(B_i^\dagger e_{ij} U_{ij} B_j)$
  4. 超导配对算符：$\Delta_{ij} = \epsilon_{ab} B_{ia} e_{ij} U_{ij} B_{jb}$

• 虚时动能项（虚时协变微商）：

  $$\mathcal{L}_{B\tau} = \sum_{i} \frac{1}{m} |\mathcal{D}_\tau B_i(\tau)|^2 = \sum_{i} \frac{|B_i(\tau+\Delta \tau) - U_{i,\tau}(\tau) B_i(\tau)|^2}{m \Delta\tau^2} \quad (1.15)$$

  其中 $m=0.1$ 为电荷子的质量参数，时空规范链接选择时空规范固定的条件：$U_{i,\tau} = \mathbf{1}_2$。

2. 关键 Benchmark 体系、计算数据与性能展示

为了完整、精确地勾勒出该强关联模型的相图，并提取相变特征，研究团队对线性尺寸 $L = 4, 8, 12, 16$ 的正方晶格（相应格点数 $N = L^2$ 最高至 256）进行了极低温度下的模拟（投影温度设为 $\beta = 2L$，即对于最大尺寸 $L=16$，$\beta = 32$），从而确保能有效滤除有限温热涨落，准确直达多体基态。

2.1 纯自旋子体系的 Deconfined-to-AFM 转变 (无电荷子极限：$r = +\infty$)

第一步，研究团队验证了当没有电荷子参与（即 $r = +\infty$ 的极限情况，等效于单纯的自旋子晶格规范场系统）时系统的物相表现。该分析发表于相图（图 1）的左边界纵轴。通过固定空间规范耦合常数 $\kappa = 20$，并改变虚时耦合常数 $\kappa_\tau$，对反铁磁（Néel）结构因子进行精确计算：

其中自旋算符定义为 $S_i^+ = f_{i\uparrow}^\dagger f_{i\downarrow}$。当波矢处于反铁磁波矢 $\mathbf{q} = \mathbf{M} = (\pi, \pi)$ 时，$\chi_S(\mathbf{M})$ 随系统的演化如图 2(a) 所示。当 $\kappa_\tau$ 增加时，$\chi_S(\mathbf{M})$ 在所有尺寸上均显著增大。为了精确界定相变点，采用规范不变量的自旋关联比（Correlation Ratio）：

其中 $\delta \mathbf{q} = (2\pi/L, 0)$ 为布里渊区最小动量分辨率。如图 2(b) 所示，对于不同尺寸 $L = 4, 8, 12, 16$ 的 $r_\text{AFM}$ 曲线展示出了极其清晰且收敛的单一交叉点。交叉点精确地对应了自旋系统的量子临界点：

• 当 $\kappa_\tau < 1.0$ 时：$r_\text{AFM}$ 随尺寸 $L$ 的增大而减小，说明在热力学极限下无磁长程序，自旋子处于解禁闭（deconfined）的狄拉克自旋液体相（Dirac Spin Liquid）。
• 当 $\kappa_\tau > 1.0$ 时：$r_\text{AFM}$ 随尺寸 $L$ 递增并趋于 $1$，表明系统在虚时规范场的强烈束缚下发生了手征对称性破缺（Chiral Symmetry Breaking），自旋子被禁闭并产生了长程反铁磁 Néel 序。

2.2 dSC 到 AFM 量子相变的共轭测度 (有限电荷子：$r = -15$ 路径)

在加入活性的玻色电荷子后，研究团队沿着图 1 相图中右侧的红色扫描路径（固定 $r = -15$，变温/变耦合强度 $\kappa_\tau$ 扫描）探索 $d$ 波超导相与反铁磁相之间的直接相变。系统的相互作用四次项参数设为：

通过对 $d$ 波配对结构因子进行测量：

其中，$\Delta_i = \sum_{j \in \text{NN}} (-1)^{\eta(j)} \Delta_{ij}$，$\eta(j)$ 代表 $d$ 波形状因子 $\{1, -1, 1, -1\}$。在 $\mathbf{q} = \mathbf{\Gamma} = (0,0)$ 处的测量数据如图 3(a) 和图 3(d) 所示。当 $\kappa_\tau$ 较低时，超导关联非常强，且其关联比 $r_{\text{el}, d}$ 逼近 $1$。随着 $\kappa_\tau$ 升高，超导序连续下降。最震撼的发现来自自旋与超导两个独立通道关联比的协同演化：

1. 超导关联比 $r_{\text{el}, d}$（图 3d）：不同尺寸的曲线在 $\kappa_\tau^c \sim 2.2(2)$ 处交汇。
2. 反铁磁关联比 $r_\text{AFM}$（图 3e）：同样在 $\kappa_\tau^c \sim 2.2(2)$ 处展现出了完美的相变交叉行为。

这表明：$d$ 波超导序的消失与反铁磁 Néel 序的涌现是在同一个临界参数 $\kappa_{\tau, c} = 2.2(2)$ 处同时发生的。这一“不约而同”的交汇行为是直接二阶连续 DQCP 相变的教科书级数值铁证。

2.3 价键固体（VBS）通道的临界涨落与涌现SO(5)对称性

为了检验 DQCP 理论的另一个核心预测——即在相变点处，原本非同源的反铁磁序与价键固体（VBS）序会由于规范场的临界涨落融合成更高维的对称性（如SO(5)），研究团队测量了二聚体关联函数及其相应的结构因子 $\chi_\text{VBS}(\mathbf{q})$（定义二聚体算符 $D_i = \mathbf{S}_i \cdot \mathbf{S}_{i+\hat{x}}$，测量动量 $\mathbf{X} = (\pi, 0)$）。

• 无长程VBS序：如图 3(c) 所示，虽然 $\chi_\text{VBS}$ 随 $\kappa_\tau$ 有所抬升，但它完全不随系统尺寸 $L$ 呈发散行为，$\chi_\text{VBS}/L^2 \to 0$，证实了在相变两边均没有发生自发的空间平移对称性破缺。
• 临界涨落极大化：如图 3(f) 所示，VBS 的关联比 $r_\text{VBS}$ 在临界点 $\kappa_\tau^c = 2.2$ 处显现出了一个极其尖锐的峰值（Peak）。这一非平庸的峰值直接反映了临界区内 VBS 的相关长度显著被拉长。这与双 flavor 的狄拉克费米子场 DQCP（其有效连续拉格朗日量如 Eq. (20) 所示）中由于4-费米项（dangerously irrelevant $L_4$）导致的涌现更高对称性（SO(5) 涨落）完全契合。

2.4 SAC 谱学分析：Dirac 准粒子行为的演化

研究团队利用先进的**随机解析延拓（SAC）**方法，从虚时 Green 函数中重构出实频电子谱函数 $A(\mathbf{k}, \omega)$。图 4 展示了系统跨越量子临界区时，光谱行为发生的深刻重构：

这一连续的实空间光谱重构过程（图 4a-c 与 图 4d-f）不仅是数值计算的成功，更在物理机理上深刻地回应了角分辨光电子能谱（ARPES）及非弹性中子散射（INS）在高温超导体相变边界上观测到的谱学演化规律。

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

3.1 DQMC 核心算法流程与代数骨架

行列式量子蒙特卡洛（DQMC）的计算核心在于通过标度变换将多体费米相互作用转化为单粒子在时空随机辅助场中的运动，并通过费米子迹（Trace）的严格解析求解，将配分函数写为单粒子传播子乘积形式：

在代码实现中，计算流程可以概括为以下伪代码步骤：

# 伪代码：DQMC-SU(2) 耦合电荷子系统核心循环
for sweep in range(total_sweeps):
    # 1. 虚时切片逐级循环
    for tau in range(beta_slices):
        # 2. 遍历晶格上的所有规范场 link U_ij
        for link in lattice_links:
            # a. 提议一个局域四元数微扰 R
            R = propose_small_su2_quaternion(variance)
            U_prime = R * U[link, tau]
            
            # b. 计算由于 U 场改变导致的费米矩阵跳跃算符微扰 Delta-matrix (4x4)
            Delta = compute_delta_matrix(U[link, tau], U_prime, link)
            
            # c. 利用 Sherman-Morrison 公式计算费米决定子的比率 (Determinant Ratio)
            ratio_fermi = calculate_det_ratio_local(Green_function[tau], Delta)
            
            # d. 计算玻色作用量变化量 (Yang-Mills 项 + 电荷子二次项 + 动能项)
            dS_boson = calculate_boson_action_change(U_prime, U[link, tau], B_field, link, tau)
            
            # e. 蒙特卡洛接受度判定 (Metropolis-Hastings)
            acceptance_prob = min(1.0, ratio_fermi * np.exp(-dS_boson))
            if random() < acceptance_prob:
                # 接受更新，利用 Fast-Update 方案更新全局格点 Green 函数
                U[link, tau] = U_prime
                update_green_function_fast(Green_function[tau], Delta)
                
    # 3. 每隔若干个 Sweep 进行全局更新 (Global Overrelaxation & Block Update)
    if sweep % 10 == 0:
        perform_block_global_update_su2_gauge(U, Green_function)
        perform_global_update_chargon_B(B_field, U)
        
    # 4. 测量物理量 (Wick 收缩展开测量自旋、超导和VBS关联)
    if sweep > thermalization_sweeps:
        measure_observables(Green_function, B_field)

3.2 四元数（Quaternion）参数化方法及其实施细节

在数值采样中，直接对满足 $U \in SU(2)$ 的 $2 \times 2$ 复矩阵进行随机扰动十分困难（因为必须严格保持幺正性 $\det U = 1$）。因此，代码使用**四元数（Quaternion）**方法进行无偏参数化。任意的 $SU(2)$ 矩阵可以唯一写为四个实数组成的超球面坐标：

在代码中，propose_small_su2_quaternion 函数首先生成满足高斯分布的微小向量 $\delta \mathbf{a} = (\delta a_0, \delta a_1, \delta a_2, \delta a_3)$（方差极小，从而控制微扰幅度），然后将该向量投影回四维单位超球面上：

通过该四元数计算出的扰动 $R$ 被左乘到原链变量上 $U'_{ij} = R U_{ij}$，由此可天然且严格地保证新矩阵依然处于 $SU(2)$ 群内，完全避免了由于浮点数累计误差导致的群流失。

3.3 块全局更新（Block Global Update）

当空间耦合常数 $\kappa_\tau$ 非常小时，系统的规范场关联极强，局域的单点随机游走更新会遭遇极其严重的“临界慢化”（Critical Slowing Down），导致自关联时间无限拉长，采样效率降为零。为此，代码必须并入块全局更新（Block Global Update）：

• 时空链锁定：选定实空间某一固定的键 $\langle ij \rangle$，对该键在所有虚时切片 $\tau = 1, \dots, \beta/\Delta \tau$ 上的幺正矩阵实施完全相同的 $SU(2)$ 变换：

  $$U_{ij}(\tau) \to R U_{ij}(\tau) \quad (3.4)$$

• 计算开销：由于该更新会一次性改变所有时空切片上的场，无法使用局域的 Fast-Update 方案。其对应的费米决定子必须进行完全重构，代数计算开销为 $\mathcal{O}(N^3)$。但为了实现物理基态的有效 ergodic 采样，这一全局操作是绝对必不可少的。

3.4 测量中的 Wick 收缩：异常 Green 函数的代数灾难

由于在规避符号问题时，对 $\tilde{f}_{\downarrow}$ 基施加了局域粒子-空穴变换（Eq. 1.4），原本单轨道的常规 Green 函数在回到物理基下测量时，会无可避免地发生“混叠”。具体来说，除了正常的 Green 函数：

代码还必须计算并存储异常 Green 函数（Anomalous Green’s Functions）：

在测量自旋或价键固体等四费米子关联函数时，利用多体 Wick 定理展开，项数会呈爆炸式增长。以自旋测量为例，其收缩展开式为：

在处理更复杂的价键固体 $\langle (\mathbf{S}_i \cdot \mathbf{S}_{i+\hat{x}}) (\mathbf{S}_j \cdot \mathbf{S}_{j+\hat{x}}) \rangle$ 时，涉及的费米算符多达 8 个，Wick 收缩项高达数十项。该研究的代码中设计了一个自动化符号推导与收缩引擎（automated algebraic contraction engine），在编译期将收缩项解析并翻译为高度优化的密集矩阵乘法（BLAS/LAPACK 库调用），从而实现了超高性能的测量计算。

3.5 复现指南与开源链接

为了复现本工作的核心结论，推荐读者参考和基于以下开源量子蒙特卡洛框架进行定制化开发：

1. ALF (Algorithms for Lattice Fermions) 软件包：
   • GitHub 仓库：https://github.com/ALF-General/ALF
   • 物理定位：ALF 是一个功能极其强大的、通用的辅助场量子蒙特卡洛（Auxiliary-Field QMC）开源平台，支持任意复杂相互作用和具有多 Flavor 的晶格规范场系统。通过定制 Hamiltonian_mod.F90，可以非常便捷地复现本文中的 $SU(2)$ 晶格规范场与自旋子耦合体系。
2. Meng Group GPU-DQMC 专用库：
   • 学术主页与开源资源：孟子阳教授团队长期致力于开发基于 CUDA 加速的高性能行列式量子蒙特卡洛代码（可参见论文引用文献 [38]）。其内部库实现了超高性能的四元数 Metropolis 快速更新以及多 GPU 并行下的块全局更新。

4. 关键引用文献深度梳理与局限性评论

4.1 核心文献脉络

本工作建立在凝聚态理论过去四十年发展的数个重要里程碑基础之上：

1. Affleck & Marston (1988) [Ref 32]：提出了正方晶格上 Heisenberg 模型的 $\pi$-flux 鞍点解。这是本文自旋子在 $\pi$-flux 背景下运动的代数发端。
2. Senthil, Vishwanath, Balents, Sachdev, Fisher (2004) [Ref 11]：开创性地提出了解禁闭量子临界点（DQCP）的理论范式，首次证明了可以通过部分子凝聚绕过 LGW 理论限制，实现 AFM 到 VBS 的二阶连续相变。
3. Sachdev 团队系列工作 (2023-2026) [Refs 11, 12, 14, 15]：将传统的磁性 DQCP 推广到超导领域，提出了反铁磁绝缘体到 $d$ 波超导体的 DQCP 理论，并预言了由于分数化导致的 Fermi 弧（Fermi Arcs）以及非平庸的伪能隙（Pseudogap）物理。本篇论文正是对该系列预言给出的最强数值实证。
4. Hermele (2007) [Ref 40]：给出了强关联 Hubbard 模型向 $SU(2)$ 晶格规范场物理模型过渡的完整路径，为本文中电荷子玻色场的引进奠定了坚实的微观理论基础。

4.2 本文方案的局限性与物理代价

尽管本工作取得了极具里程碑意义的成果，但在对其进行学术评估时，我们必须理性看待其为了“免受负符号问题困扰”而在物理上付出的几处关键代价（局限性）：

1. 严格局限于半满（Half-filling, $p = 0$）与粒子-空穴对称性

• 局限性阐述：无符号问题的代数证明（Eq. 1.6）完全依赖于克拉默简并。这意味着系统必须处于半满状态，并且具有粒子-空穴对称性。然而，真正具有高 $T_c$ 超导性质的铜氧化物系统是在**掺杂（doped, $p > 0$）**条件下工作的。一旦引入空穴掺杂以偏离半满，时间反演对称性将立即被打破，负符号问题会卷土重来。
• 学术评论：虽然作者在文中辩称“非零掺杂下的系统物理本质仍将继承自半满处的 DQCP”，但这毕竟是一种外推。我们依然无法直接在数值上通过该模型计算真正掺杂铜氧化物的低温性质。

2. “部分子”表象的物理实在性争议

• 局限性阐述：在真实物质中，电子是一个不可分割的凝聚态基本单元。尽管在强关联下可以出现部分子分数化，但在低能下，规范场往往会走向禁闭（confinement），使得自旋子和电荷子重新组合回物理电子。本文中所模拟的解禁闭相，其参数区间是否对应真实铜氧化物的真实物理参数，仍有待进一步的实验检验。

3. 随机解析延拓（SAC）的固有病态本质

• 局限性阐述：从虚时 Green 函数 $G(\tau)$ 求解实频谱函数 $A(\omega)$ 本质上是一个反问题，数学上是一个极其病态的积分方程求解：

  $$G(\tau) = \int K(\tau, \omega) A(\omega) d\omega \quad (4.1)$$

  即使 SAC 算法在统计学上非常先进，它仍然高度依赖于噪声水平和人为设置的正则化参数。这导致所得出的电子色散谱（图 4）在细节和谱权重的绝对分布上，可能会存在由于平滑化处理导致的伪影，无法达到像角分辨光电子能谱（ARPES）那样绝对精准的能量分辨率。

5. 跨学科视角：量子化学与凝聚态物理的碰撞

5.1 量子化学视角下的强关联电子系统多体方法对比

对于量子化学（Quantum Chemistry）背景的研究者来说，分子和固体的电子结构计算通常围绕着密度泛函理论（DFT）、耦合簇理论（CCSD(T)）以及多配置自洽场（CASSCF）展开。然而，当面对如过渡金属氧化物（包括铜氧化物）这类的典型强关联电子系统时，传统的量子化学方法会遭遇极大的挑战：

• 静态关联（Static Correlation）灾难：在强关联下，单决定子近似完全失效，费米面附近的电子态高度简并。CASSCF 能够处理静态关联，但其活动空间（Active Space）的开销随轨道数呈指数级增长，通常最多只能处理约 18 个轨道（CAS(18,18)），这对于二维无穷大晶格系统来说无异于杯水车薪。
• 密度矩阵重整化群（DMRG）的维度局限：DMRG 能够精确处理一维或伪一维条带（Strips）系统。但在真正的二维晶格上，由于纠缠熵满足面积律（Area Law），DMRG 的键合维度（Bond Dimension）需要随系统宽度呈指数级增加，难以在大尺寸二维系统上提供完全收敛的数据。

本工作所展示的行列式量子蒙特卡洛（DQMC）在处理此类问题时展现了无与伦比的跨越性优势。只要设计得当（避免负符号问题），DQMC 的计算复杂度仅随格点数 $N$ 呈多项式标度（通常为 $\mathcal{O}(\beta N^3)$），这使得它能够直接在热力学极限下计算真正的二维体相材料。它为量子化学家提供了一条新路径：在面对复杂的晶格多体系统时，放弃从分子轨道出发的繁琐构建，转而借助规范场与分数化部分子的物理底色，在多体有效晶格模型上直接开展大规模无偏采样，以此直接获取材料基态的非平庸物相特征。

5.2 非零掺杂的 FL* 理论与费米面重构（Fermi Surface Reconstruction）

虽然本项工作在半满处探索了直接的二阶相变，但它对于掺杂下铜氧化物相图中的“伪能隙（Pseudogap）”和“费米弧（Fermi Arcs）”现象有着极其深刻的指导作用。在传统的 Luttinger 定理限制下，费米面的体积必须严格由系统中的总电子数决定：

其中 $p$ 为掺杂浓度。在空穴掺杂下，这应当对应一个大费米面。

然而，实验却在伪能隙区观测到了不满足 Luttinger 定理的小费米面（即费米弧）。本工作的理论基础指出，在半满 DQCP 附近的有限掺杂区，系统将诱导出一个被称为**分数化费米液体（Fractionalized Fermi Liquid, FL*）**的非平庸金属相。在 FL* 金属相中，存在一个背景自旋液体，其自旋子将多体纠缠锁定在后台，不直接贡献于电荷输运。结果，物理电子只需要在剩余的自由度上形成 Fermi pockets。这些 Fermi pockets 的有效面积仅正比于掺杂浓度 $p$（而非 $1-p$）：

这完美且自然地解释了高温超导体中困扰理论界数十年的小费米面之谜。本工作在半满处给出的确定性 QMC 临界数据，为这一大一小费米面转换的 FL* 理论大厦打下了最坚实的一块数值地基。

5.3 跨学科展望：凝聚态与量子化学的深度融合

未来的研究应当致力于打破领域藩篱。对于量子化学家而言，可以尝试将第一性原理计算（如通过 DFT 算得的 Wannier 轨道和局域库仑作用量 $U$）作为微观输入参数，映射到类似于本文所提出的无符号问题的自旋子-电荷子晶格哈密顿量中。通过 DQMC 算法直接计算出其低温物理行为（如超导转变温度 $T_c$ 趋势、自旋激发能谱），从而真正实现强关联超导材料的“自下而上”（bottom-up）免经验式精准预测和理性设计。这无疑将是未来凝聚态物理学与量子化学深度交叉融合的最激动人心的前沿方向之一。