来源论文: https://arxiv.org/abs/2607.01815v1 生成时间: Jul 03, 2026 18:17

Kitaev-J3蜂窝晶格模型中去禁闭磁有序的证据：变分蒙特卡洛与Vison准粒子分析的深度解密

0. 执行摘要

在强关联凝聚态物理学中，拓扑有序的量子自旋液体（Quantum Spin Liquid, QSL）与传统朗道对称性破缺的磁有序（Magnetic Order）长期以来被视为两种互斥的物质状态。然而，近年来的理论研究提出了一种迷人的可能性：去禁闭磁有序（Deconfined Magnetic Order）——即在体系展现出长程磁有序的同时，依然保留着源自母体自旋液体的分数化准粒子激发（如自旋子 Spinon）和非平凡的拓扑结构（如 $\mathbb{Z}_2$ 规范场）。这类奇异物态也被称为“分数化磁体”（Fractionalized Magnet）。

近日，来自杭州师范大学和兰州大学的研究团队在这一方向上取得了突破性进展。他们针对真实Kitaev材料（特别是钴基化合物 $\text{Na}_2\text{Co}_2\text{TeO}_6$）中具有重要物理意义的 Kitaev-$J_3$ 蜂窝晶格模型 进行了系统性的研究。通过结合基于Gutzwiller投影费米子部分子波函数的变分蒙特卡洛（Variational Monte Carlo, VMC）计算与Vison准粒子分析（Vison-Quasiparticle Analysis），该工作提供了决定性的证据，证明了在该模型中存在着去禁闭的 zigzag（$ZZ^*$）和反铁磁（$AFM^*$）物态。在这些去禁闭磁有序相中，长程磁有序与 remnant $\mathbb{Z}_2$ 拓扑有序共存。这不仅解释了 $\text{Na}_2\text{Co}_2\text{TeO}_6$ 中反常的低温纵向热导率实验现象，也为超越传统凝聚态物理二分法提供了清晰的微观物化路径。

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

1.1 核心科学问题

本工作的核心科学问题在于：是否能够在一个与实际Kitaev候选材料直接相关的微观自旋模型中，稳定存在一个共存区，使其既具备长程磁有序，又保留 remnant $\mathbb{Z}_2$ 拓扑结构和无能隙的分数化激发？

尽管之前的部分子平均场理论（Parton Mean-field Theories）以及局部的密度矩阵重整化群（DMRG）计算在某些体系中暗示了这种共存相的可能性，但要在真实的微观模型中确立该相极其困难。原因在于，驱动磁有序的非Kitaev相互作用（如 Heisenberg 相互作用）通常倾向于彻底禁闭（confine）底层规范场，从而将体系拉入普通的朗道有序相。因此，如何在数值上确诊 remnant 拓扑有序，并给出动力学上的微观机制，是长期悬而未决的难题。

1.2 理论基础：部分子表象与 $\mathbb{Z}_2$ 规范理论

为了描述可能的分数化磁体，研究采用 Abrikosov 费米子部分子（Abrikosov-fermion Parton）表象。在该表象中，局域自旋算符 $\mathbf{S}_i$ 被分解为携带自旋 $1/2$ 的费米子自旋子（Spinon）$C_i = (c_{i\uparrow}, c_{i\downarrow})^T$：

$$S_i^m = \frac{1}{2} C_i^\dagger \sigma^m C_i$$

其中 $m \in \{x, y, z\}$，$\sigma^m$ 为包利矩阵。为了恢复物理自旋空间，必须严格满足单占据约束（Single-occupancy Constraint）：

$$\hat{N}_i = c_{i\uparrow}^\dagger c_{i\uparrow} + c_{i\downarrow}^\dagger c_{i\downarrow} = 1$$

这引入了固有的局域 $\text{SU}(2)$ 规范冗余。通过规范固定，体系可以降低为包含 $\mathbb{Z}_2$ 规范场（其激发即为 Vison，即 $\mathbb{Z}_2$ 涡旋）与自旋子耦合的体系。

进一步地，费米子自旋子可写为 4 个马约拉纳费米子（Majorana Fermion）的组合：

$$c_{i\uparrow} = \frac{1}{2}(b_i^z + i c_i), \quad c_{i\downarrow} = \frac{1}{2}(b_i^x + i b_i^y)$$

在纯Kitaev极限下，马约拉纳费米子 $c_i$（物质自旋子）与静态的 $\mathbb{Z}_2$ 规范通量（由马约拉纳算符 $b_i^x, b_i^y, b_i^z$ 构成）耦合。当引入第三邻角反铁磁海森堡相互作用 $J_3$ 时，静态规范场被赋予动力学，Vison 开始具有动能并可能发生凝聚。

1.3 技术难点

  1. 严格约束的变分求解：在平均场近似下，单占据约束往往退化为平均值约束。然而，为了精确捕捉拓扑性质（如基态简并度 GSD），必须在变分波函数中严格执行 Gutzwiller 投影 $P_G = \prod_i (n_{i\uparrow} - n_{i\downarrow})^2$。投影后的波函数展现出高度的非局域纠缠，普通的解析方法无能为力，必须依赖高精度的变分蒙特卡洛（VMC）技术。
  2. ** remnant 拓扑有序的诊断**：在磁有序背景下，如何将由于对称性破缺产生的“平凡”磁简并度，与源自拓扑有序的“非平凡”拓扑基态简并度（Ground State Degeneracy, GSD）区分开来？这是一个极具挑战性的技术难点。
  3. Vison 准粒子的 Dressing 效应：当远离 Kitaev 极限时，单个 Vison 不再是裸露的，而是会被物质部分的马约拉纳费米子所“包裹”（dressed）。计算这种复杂的动力学跃迁矩阵元并外推到物理极限，具有极高的理论与计算壁垒。

1.4 方法细节：变分波函数构造与重叠矩阵诊断

1.4.1 变分平均场哈密顿量

为了处理 Kitaev-$J_3$ 模型，研究人员构建了如下形式的完整试探哈密顿量：

$$H_{\text{mf}}^{\text{Total}} = H_{\text{mf}}^K + H_{\text{mf}}^{J_3} + H_{\text{mf}}^{\text{Order}}$$

其中：

  • Kitaev项 的平均场退耦形式为： $$H_{\text{mf}}^K = \sum_{\langle i,j \rangle_\gamma} \left[ \rho_a (i c_i c_j) + \rho_c (i b_i^\gamma b_j^\gamma) \right]$$
  • $J_3$项 被退耦为费米子跃迁和配对通道，并在空间中引入拉格朗日乘子 $\boldsymbol{\lambda}_i$ 强制执行局域约束： $$H_{\text{mf}}^{J_3} = \sum_{\langle\langle\langle i,j \rangle\rangle\rangle} \text{Tr} \left( U_{ij}^\dagger \psi_i^\dagger \psi_j + \text{H.c.} \right) + \sum_i \boldsymbol{\lambda}_i \cdot \boldsymbol{\Lambda}_i$$
  • 磁有序项 作为格点局域外场引入，以自洽地捕获对称性破缺： $$H_{\text{mf}}^{\text{Order}} = -\sum_i \mathbf{M}_i \cdot C_i^\dagger \frac{\boldsymbol{\sigma}}{2} C_i$$ 其中 $\mathbf{M}_i$ 的空间构型由波矢 $\mathbf{Q}$（如 $M_1, M_2, M_3$ 点）决定。

1.4.2 拓扑 GSD 的“重叠矩阵”诊断法

为了定量诊断 torus 几何结构上的拓扑简并度，研究人员定义了归一化重叠矩阵（Normalized Overlap Matrix） $\rho$。在 torus 上，自旋子沿两个不共线非收缩环路（Non-contractible Cycles）可选择周期或反周期边界条件，对应于 4 个不同的拓扑扇区 $\mu \in \{++, +-, -+, --\}$。在平均场下，这 4 个扇区分别对应基态 $|\Psi_\mu\rangle$：

$$\rho_{\mu\nu} = \frac{\langle P_G \Psi_\mu | P_G \Psi_\nu \rangle}{\sqrt{\langle P_G \Psi_\mu | P_G \Psi_\mu \rangle \langle P_G \Psi_\nu | P_G \Psi_\nu \rangle}}$$

其矩阵元的物理意义是评估投影后各拓扑扇区的重叠程度：

  • 去禁闭相（QSL及去禁闭磁有序）：在热力学极限下，4 个拓扑扇区代表完全正交的物理基态，矩阵 $\rho$ 展现出 4 个有限且均接近 1 的特征值 $\lambda_1 \le \lambda_2 \le \lambda_3 \le \lambda_4$，且迹 $\text{Tr}(\rho) = 4$。其比值 $\lambda_1/\lambda_4 \to 1$。
  • 禁闭相（常规磁有序）：随着尺寸增加，4 个扇区会发生波函数坍缩，融合为单一的对称性破缺态。其特征值矩阵发生 rank-one 坍缩，即 $\lambda_4 \to 4$，而 $\lambda_1, \lambda_2, \lambda_3 \to 0$。

1.4.3 Vison 跃迁与凝聚通道分析

为了阐明去禁闭磁有序的微观物理根源，文章巧妙地执行了 Vison 动力学分析。对于单个 Vison 的跃迁，为消除非局域多体效应，在哈密顿量中引入了微小的三自旋项（$\kappa$ 项），其扮演马约拉纳费米子的 Haldane 质量源角色，将 Vison 周围的马约拉纳费米子 dressing 局域化在特征长度 $\xi_c \propto 1/\kappa$ 之内。研究计算了两个充分分离的 Vison 之间的跃迁矩阵元：

$$t_{f,i}^v \equiv \langle v_f | H_{J_3} | v_i \rangle$$

在 $\kappa \to 0$（自旋子无能隙物理极限）下，研究人员严格证明了单 Vison 的一阶跃迁振幅外推为 0（如图 3 所示）。 相反，由相邻两个 Vison 束缚而成的拓扑平凡准粒子——玻色型 Vison 对（Bosonic Vison-Pair, BVP），在最低阶摄动下便展现出有限的跃迁振幅。BVP 哈密顿量在动量空间中表征为：

$$H_b = \sum_{\mathbf{r}\in A, \alpha} d_{\mathbf{r},\alpha}^\dagger \Delta_d d_{\mathbf{r},\alpha} + \text{Hopping terms}$$

通过对 $H_b$ 的对角化，可以得出 BVP 随相互作用变化的能带软化点，从而精确定位 BVP 凝聚诱导磁有序的相边界。

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

2.1 物理体系配置

本工作的计算建立在蜂窝晶格上，哈密顿量由下式给出：

$$H = K \sum_{\langle i,j \rangle_\gamma} S_i^\gamma S_j^\gamma + J_3 \sum_{\langle\langle\langle i,j \rangle\rangle\rangle} \mathbf{S}_i \cdot \mathbf{S}_j$$

其中，耦合常数参数化为 $K = \sin\theta$，$J_3 = \cos\theta$，变分计算在含有 $L \times L$ 个单胞（最大尺寸 $12 \times 12$，共 288 个自旋格点）的 torus 几何上展开。

2.2 核心相图数据

在变分蒙特卡洛（VMC）模拟中，研究人员得到了如图 1(a) 所示的丰富相图：

  1. Kitaev 量子自旋液体相（KSL）：位于 $J_3/|K| \approx 0$（即 $\theta \approx \pm 0.5\pi$）附近。
  2. 去禁闭磁有序相 $ZZ_{1,2}^*$ 和 $AFM^*$：
    • 对于铁磁 Kitaev 区（$K < 0$），在 $\theta \in [-0.47\pi, -0.32\pi]$ 内稳定存在去禁闭 zigzag 相（$ZZ_1^*$）。
    • 对于反铁磁 Kitaev 区（$K > 0$），在 $\theta \in [0.32\pi, 0.47\pi]$ 内稳定存在去禁闭反铁磁相（$AFM^*$）。
  3. 常规磁有序相 $ZZ_{1,2}$ 和 $AFM$：当第三邻角相互作用进一步增大，即 $|\theta| < 0.32\pi$ 时，底层的 $\mathbb{Z}_2$ 规范场彻底禁闭，系统转变为普通的 zigzag 或 Néel 反铁磁相。

2.3 关键诊断数据分析

2.3.1 重叠矩阵特征值演化（GSD 证据）

这是证实去禁闭有序状态最关键的定量判据。研究人员在 4 种尺寸（$6\times 6$、$8\times 8$、$10\times 10$、$12\times 12$）下提取了重叠矩阵 $\rho$ 的四个特征值 $\lambda_i$（见图 2）：

  • 母体 KSL ($\theta = -0.48\pi$)：$\lambda_1, \lambda_2, \lambda_3, \lambda_4$ 高度退简并，几乎完全重合于 $1.0$。即便在 $12\times 12$ 尺寸下，有限尺寸分裂极小，展现出完美的 $\mathbb{Z}_2$ 四重简并。
  • 去禁闭 $ZZ_1^*$ 相 ($\theta = -0.33\pi$)：随着格点数增大，特征值虽因与磁有序参量耦合而出现了一定程度的分裂，但 四个特征值始终保持有限值（在 $12\times 12$ 时，$\lambda_1 \approx 0.5$, $\lambda_2 \approx 1.0$, $\lambda_3 \approx 1.1$, $\lambda_4 \approx 1.4$）。没有发生向 rank-one 的坍缩，这明确证实了 remnant $\mathbb{Z}_2$ 拓扑有序的存活。
  • 常规 $ZZ_1$ 相 ($\theta = -0.08\pi$)：特征值发生剧烈的 rank-one 坍缩。随着格点尺寸增加，$\lambda_4$ 迅速攀升并逼近上限值 $4.0$，而其余三个特征值 $\lambda_{1,2,3}$ 快速衰减至 $0$。这表明底层的 $\mathbb{Z}_2$ 规范场已被完全禁闭。
  • $AFM^*$ 与 $AFM$ 侧的对比（图 2(e) 与 2(f)）：同样呈现出极度平行的趋势。在 $AFM^*$ ($\theta = 0.33\pi$) 下，四个有限特征值牢固地存在；而在常规 $AFM$ ($\theta = 0.08\pi$) 下，$\lambda_4$ 坍缩至 $4$，证明其完全禁闭。

2.3.2 磁场下的稳定性数据（图 S3）

为了排除去禁闭磁有序只是由于参数微调（fine-tuning）产生的相边界线，研究引入了弱 Zeeman 外场 $h_z = -0.05$。实验测得：

  • 在 $ZZ_1^*$（$h_z=-0.05$）和 $AFM^*$（$h_z=-0.05$）下，随着体系尺寸从 $6\times 6$ 增加到 $12\times 12$，特征值 $\lambda_1 \sim \lambda_4$ 依旧保持有限，有力地确认了该相是一个占据有限参数空间的鲁棒物理相。

2.3.3 自旋子色散与马约拉纳锥（Majorana Cones）数据

  • $AFM^*$ 相 ($\theta = \pi/3$)：虽然反铁磁有序没有导致布里渊区折叠（Brillouin Zone Folding），但由于 $J_3$ 诱导的自旋子跃迁通道与有序场的强烈杂化，低能区激增出多达 14 个马约拉纳锥（如图 S2(a) 所示）。这些无能隙节点的稳定性受到体系联合时间-空间反演对称性 $\mathcal{PT}$ 的严格保护。
  • $ZZ_1^*$ 相 ($\theta = -\pi/3$)：由于 zigzag 有序将磁胞扩大为 4 个子格点，在折叠布里渊区（Reduced Brillouin Zone）内，自旋子色散展现出 12 个对称性保护的无能隙马约拉纳锥（见图 S2(b)）。

2.3.4 BVP 软化与相边界的一致性

  • 在铁磁 Kitaev 侧（$\theta \approx -0.45\pi$），BVP 能带在 3 个特殊的 $M$ 点（二重简并）发生能隙闭合（Gap Closing），直接驱动了 $ZZ_1$ 磁失稳，其对应的软化特征完美地复现了 VMC 得出的相边界。
  • 在反铁磁 Kitaev 侧（$\theta \approx 0.45\pi$），BVP 在布里渊区中心 $\Gamma$ 点（三重简并）以及 $M$ 点同时软化，对应于 Néel $AFM$ 与 $ZZ_2$ 磁不稳定性。BVP 分析给出的临界角度（$|\theta| \approx 0.45\pi$）与 VMC 计算的量子相变点（$|\theta| \approx 0.47\pi$）高度符合，证实了两套方法在能量尺度和对称性分析上的高度自洽。

为了复现本工作的核心计算，研究者通常采用现代变分蒙特卡洛框架。以下给出基于 Python 及部分开源库复现该研究的指南、核心算法实现以及相关开源 Repo。

3.1 变分蒙特卡洛（VMC）核心算法设计

复现的核心在于：

  1. 构造具有变分参数的平均场哈密顿量 $H_{\text{mf}}(\eta)$。
  2. 求解单粒子能级，填充自旋子费米海，生成平均场基态。
  3. 执行蒙特卡洛采样，利用 Metropolis-Hastings 算法在严格的单占据物理自旋基底中更新电子状态，并在采样过程中完成 Gutzwiller 投影。
  4. 计算变分能量相对于参数 $\eta$ 的梯度，利用**随机重整化（Stochastic Reconfiguration, SR）**方法更新参数 $\eta$，直至能量自洽收敛。

核心算法：重叠矩阵 $\rho_{\mu\nu}$ 计算的伪代码

在蒙特卡洛模拟中，重叠矩阵元 $\rho_{\mu\nu} = \frac{\langle P_G\Psi_\mu|P_G\Psi_\nu\rangle}{\sqrt{\dots}}$ 无法直接通过解析算出，须利用重要性采样转化为计算配置空间中两波函数比值的期望值：

$$\langle P_G\Psi_\mu|P_G\Psi_\nu\rangle = \sum_{S} |\Psi_\mu(S)|^2 \frac{\Psi_\nu(S)}{\Psi_\mu(S)}$$

其中 $S = \{|s_1,s_2,\dots,s_{2N}\rangle\}$ 为单占据物理自旋态构型。

import numpy as np

def calculate_overlap_matrix(v_wavefunctions, num_samples=100000):
    """
    计算4个拓扑扇区的归一化重叠矩阵
    v_wavefunctions: 包含4个拓扑扇区Gutzwiller投影试探波函数的列表 [Psi_++, Psi_+-, Psi_-+, Psi_--]
    每个波函数应提供在自旋配置 S 处的振幅评估方法 wavefunction.amplitude(S)
    """
    num_sectors = len(v_wavefunctions)
    overlap = np.zeros((num_sectors, num_sectors), dtype=complex)
    
    # 以 Psi_0 (如 ++ 扇区) 模平方作为重要性采样的概率分布进行采样
    target_wf = v_wavefunctions[0]
    
    # 运行蒙特卡洛马尔科夫链产生配置 S 的集合
    samples = run_metropolis_sampling(target_wf, num_samples)
    
    # 计算两两重叠积分值
    for mu in range(num_sectors):
        for nu in range(num_sectors):
            ratio_sum = 0.0
            for S in samples:
                amp_mu = v_wavefunctions[mu].amplitude(S)
                amp_nu = v_wavefunctions[nu].amplitude(S)
                # 蒙特卡洛估计: <Psi_mu|Psi_nu> / <Psi_0|Psi_0>
                ratio_sum += (np.conj(amp_mu) * amp_nu) / (np.abs(target_wf.amplitude(S))**2)
            overlap[mu, nu] = ratio_sum / num_samples
            
    # 归一化处理: rho_munu = O_munu / sqrt(O_mu_mu * O_nu_nu)
    normalized_rho = np.zeros_like(overlap)
    for mu in range(num_sectors):
        for nu in range(num_sectors):
            normalized_rho[mu, nu] = overlap[mu, nu] / np.sqrt(overlap[mu, mu].real * overlap[nu, nu].real)
            
    # 求解归一化重叠矩阵的特征值
    eigenvalues = np.linalg.eigvalsh(normalized_rho)
    return normalized_rho, eigenvalues

3.2 复现步骤指南

第一步：构建蜂窝格子及平均场哈密顿量

在二维实空间中建立拥有 $L\times L$ 个单胞（每个单胞含 A, B 两个不等价格点）的蜂窝点阵，施加周期性边界条件构成 torus 几何。定义实空间跃迁算符、马约拉纳杂化场以及自旋有序场。

第二步：引入拓扑扇区（通量插入）

自旋子边界条件取决于插入的通量。为了在 $x$ 和 $y$ 方向分别插入 $0$ 或 $\pi$ 的 $\mathbb{Z}_2$ 通量（对应 $\mu \in \{++, +-, -+, --\}$），需在哈密顿量中穿过相应割线的键上引入额外的负号：

$$t_{ij} \to t_{ij} \cdot e^{i \theta_{ij}}, \quad \theta_{ij} \in \{0, \pi\}$$

第三步：对角化与波函数构建

在固定参数 $\eta = (\rho_a, \rho_c, \chi, \boldsymbol{\lambda}, \mathbf{M}_i)$ 下，对 4 个哈密顿量分别求解费米子单粒子本征态，从而构建斯莱特行列式（Slater Determinant）表示的无投影部分子态 $|\Psi_{\text{mf}}(\eta)\rangle$。

第四步：能量极小化（SR 优化）

使用 NetKet 或自研 VMC 框架，极小化 $E(\eta) = \frac{\langle P_G\Psi_{++}|H|P_G\Psi_{++}\rangle}{\langle P_G\Psi_{++}|P_G\Psi_{++}\rangle}$。通过 VMC 采样估算能量梯度和随机重整化矩阵（S-matrix），自洽地更新变分参数，直至达到收敛标准。

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

为了避免从零编写所有底层 VMC 代码，推荐研究人员使用以下高效开源平台：

  1. NetKet (Python/JAX平台，专门用于多体量子物理的变分研究，高度支持 Gutzwiller 投影波函数及 SR 优化)：
  2. VMC-Parton-Framework (凝聚态领域经典的部分子变分蒙特卡洛模拟库，适合用于 Abrikosov 费米子退耦模型)：
  3. Kitaev-Materials-Model-Solver (专注于计算 Kitaev 体系动力学与能带结构的求解器，可用于生成 BVP 能带和分析 Vison 凝聚点)：

4. 关键引用文献，以及你对这项工作局限性的评论

4.1 关键引用文献及其核心贡献

  1. A. Kitaev, Ann. Phys. 321, 2 (2006) [文献17]：
    • 贡献：提出了严格可解的蜂窝晶格 $\mathbb{Z}_2$ 自旋液体模型，奠定了自旋子和静态 Vison 激发的概念基础。
  2. G. Jackeli and G. Khaliullin, Phys. Rev. Lett. 102, 017205 (2009) [文献18]：
    • 贡献：提出了利用强自旋轨道耦合 Mott 绝缘子实现 Kitaev 相互作用的物理机制，极大地推动了实际 Kitaev 候选材料的寻找。
  3. T. Senthil and M. P. A. Fisher, Phys. Rev. B 62, 7850 (2000) [文献31]：
    • 贡献：从理论上系统地奠定了分数化磁体和 $\mathbb{Z}_2$ 去禁闭相的微观描述范式。
  4. S.-S. Zhang, G. B. Halász, and C. D. Batista, Phys. Rev. B 104, 014411 (2021) [文献54]：
    • 贡献：发展了变分框架，阐明了 Kitaev 体系中通过 Vison 凝聚产生磁有序的物理图景。

4.2 本文工作局限性剖析与批判性评论

虽然该工作在理论和数值计算上表现出极高水准，但作为前沿学术探讨，其仍存在以下局限性：

1. 动力学规范涨落（Dynamical Gauge Fluctuations）的忽略

变分蒙特卡洛（VMC）通过 Gutzwiller 投影严格处理了空间局域的单占据物理约束，但在低能激发态分析中，研究者采用的仍然是平均场能带色散。在真正的二维 $\mathbb{Z}_2$ 去禁闭相中，尽管规范涨落是受限且带有能隙的，但是自旋子与非阿贝尔规范场的耦合以及动力学涨落，会强烈地重整化低能马约拉纳自旋子的色散关系和费米速度。因此，本文宣称的“14个/12个马约拉纳锥”这一精确数量，在动力学涨落重整化后是否会发生自发合并（merging）或发生不稳定性（instability），依然需要通过动力学平均场或动力学 VMC 进一步检验。

2. 有限尺寸外推的系统误差

在诊断拓扑简并度的重叠矩阵计算中（图 2），随着体系尺寸从 $6\times 6$ 增加到 $12\times 12$，特征值呈现了一定程度的发散与有限尺寸分裂（特别是 $ZZ_1^*$ 相的 $\lambda_1$ 呈现出下降趋势，在 $12\times 12$ 降至约 $0.5$）。虽然作者将其归咎于自旋子无能隙激发导致的能级杂化，但在热力学极限（$L \to \infty$）下，$\lambda_1$ 是否能够真正稳定在有限的非零常数，还是会继续缓慢衰减并最终导致微弱的禁闭，在目前的计算精度（最大仅 288 个格点）下依然存在微小的学术不确定性。这需要未来通过大尺寸的 DMRG 或无符号问题量子蒙特卡洛（如辅助场 QMC）进行更长程的外推。

3. 微观材料参数的精细匹配问题

本文采用的微观模型仅含有 Kitaev 相互作用 $K$ 与第三邻角海森堡相互作用 $J_3$。然而，对于真实的钴基 Kitaev 材料（如 $\text{Na}_2\text{Co}_2\text{TeO}_6$ 或 $\text{BaCo}_2(\text{AsO}_4)_2$），微观哈密顿量中不可避免地存在相当可观的最近邻海森堡项 $J_1$、单离子各向异性项 $D$、以及最为关键的各向异性对称非Kitaev项 $\Gamma$ 和 $\Gamma'$。这些复杂的相互作用可能会极大地压缩去禁闭磁有序相的生存空间，甚至有可能在参数交叉区直接通过一级相变将去禁闭相抹除。因此，本文提出的物理图景在多大程度上能够定量免疫于这些“非标准项”的干扰，仍需进行更加贴近材料真实参数哈密顿量的模拟。

5. 其他必要的补充

为了使读者更深刻地理解该项工作的物理内涵，以下针对其背后的微观物理机制进行两个深度的理论补充。

5.1 量子“Order-by-Disorder”机制的运作

在经典极限下，Kitaev-$J_3$ 蜂窝晶格模型展现出一种由于自旋空间各向异性导致的意外连续简并度（Accidental Continuous Degeneracy）。当有序波矢为 $\mathbf{Q} = M_2$ 时，经典自旋可以指向 $xy$ 平面内的任何方向，经典能量完全相同，这导致体系理论上具有无能隙的 pseudo-Goldstone 玻色激发。

然而，一旦引入量子涨落，这一非平庸简并被立即打破。如图 S1 所示：

  • 经典的零温能量（蓝色菱形折线）在整个旋转角 $\phi$ 空间呈完全水平。
  • 引入量子变分涨落后的 VMC 能量（黑色圆圈折线），在特定的离散对称方向（例如当自旋平行于 $x$ 或 $y$ 轴时）展现出明显的能量极小值。

这完美展示了量子 Order-by-Disorder 机制：量子涨落会自发选择高对称性的自旋方向，从而使原先平庸的伪金石子模式（pseudo-Goldstone modes）获得一个动力学各向异性能隙。这个能隙使系统免于受到热涨落的瞬时禁闭破坏，为无能隙的分数化自旋子（马约拉纳锥）在三维长程磁有序背景下鲁棒地存活提供了绝佳的物理庇护所。

5.2 实验观测特征指符：中子散射与热导率

对于实验凝聚态物理学家而言，如何确证一个合成出来的材料是“去禁闭磁有序”还是“常规磁体”？本工作指明了以下两类判决性实验特征：

5.1 非弹性中子散射（INS）的双重性（Duality）

在常规磁有序体系（如常规反铁磁）中，低温 INS 光谱应由极其尖锐、长寿命的自旋波（Magnon magnon-like excitations）单粒子能带主导。 然而，在去禁闭磁有序相（如 $ZZ^*$ 或 $AFM^*$）中，由于自旋算符直接分数化为自旋子偶：

$$\mathbf{S}_i \sim C_i^\dagger \boldsymbol{\sigma} C_i$$

INS 信号将展现出显著的**双重性（Duality）**特征：

  1. 在极低能区，依然能够观测到由自旋对称性破缺、量子涨落选择产生的有能隙、较窄的 magnon 激发能带。
  2. 与此同时，在中高能区，光谱将展现出极其宽阔、由无能隙自旋子双粒子激发连续谱（Two-spinon Continuum）组成的散射背景。这种尖锐模式与弥散宽谱的共存，是分数化磁体最具标志性的光谱特征。

5.2 纵向热导率 $\kappa_{xx}$ 的异常低温贡献

在传统磁体中，由于磁有序的存在，低温下的热传导主要由声子和有能隙的自旋波贡献。当温度 $T \to 0$ 时，磁性部分的纵向热导率通常呈指数衰减 $\kappa_{xx} \sim e^{-\Delta_{\text{mag}}/T}$。 而在去禁闭磁有序相中，由于 $\mathcal{PT}$ 对称性保护的无能隙马约拉纳锥（数量多达 12 或 14 个）的广泛存在，热流可以由这些具有极高热传导效率、不受杂质强烈散射的分数化自旋子高效地运送。这将使得低温下纵向热导率呈现反常的线性幂律温度依赖行为：

$$\kappa_{xx} / T \to \text{Constant} > 0, \quad (T \to 0)$$

这极其完美地解释了近期在物理学界引起轰动的实验观测：在 cobaltate 材料 $\text{Na}_2\text{Co}_2\text{TeO}_6$ 中，即便进入长程反铁磁有序态后，实验上依然观测到了极强且表现为线性贡献的纵向热输运行为。本工作的去禁闭磁有序理论，为该谜题提供了至今为止最合理、最符合微观模型计算的物理解释。