来源论文: https://arxiv.org/abs/2606.05837v1 生成时间: Jun 06, 2026 15:50

拓扑激子凝聚与奇宇称磁振子：Haldane-Hubbard 模型的集体激发态深度理论解析

0. 执行摘要

自 1988 年 Haldane 提出无朗道能级的量子反常霍尔效应模型以来，拓扑物态的设计与微观物理解析一直是凝聚态物理和量子化学的核心前沿。随着实验技术（如 Floquet 光学驱动、冷原子光晶格）的发展，研究视线已从单纯的单粒子拓扑能带结构，拓展至强关联体系下的集体激发态。关联拓扑体系中的一个基础课题是：当电子关联与本征拓扑相互竞争、交织时，系统的集体激发（如激子、磁振子）展现出何种非平凡的拓扑物性？

本工作针对经典的 Haldane-Hubbard (HH) 模型，首次在全动量空间及整个相图范围内，系统地研究了其在顺磁（Paramagnetic）相和共线 Néel 磁有序相中的集体激发态动力学。利用自共轭时协 Hartree-Fock (TDHF) 近似与自洽随机相位近似 (RPA) 理论：

1. 拓扑激子的发现：在顺磁相中，由于电子能带本身的非平凡拓扑（Chern 数 $C = \pm 1$），体系在粒子-空穴连续区下方存在尖锐的、空间离散的拓扑激子束缚态，其 Chern 数满足两带差值：$C_{\text{exc}}^{\sigma\to\bar{\sigma}} = C_{c\bar{\sigma}} - C_{v\sigma} = -2$。
2. 激子凝聚机制：通过调节 Hubbard 斥力 $U$ 与次近邻跃迁 $t_2$，发现在相边界处，该拓扑激子在 $\Gamma$ 点（零动量）发生能量崩塌，其**凝聚（Condensation）**直接驱动了向共线 Néel 反铁磁相的二级量子相变，建立了拓扑激子凝聚与共线磁有序的直接对应关系。
3. 奇宇称磁振子与 $f$ 波分裂：在磁有序相中，凝聚后的激子转化为线性色散的 Type-I 金斯特（Goldstone）玻色子——磁振子。由于次近邻复跃迁（对应自发环形电流）打破了时间反演对称性 $\mathcal{T}$，但保留了空间反演与自旋旋转的联合对称性，磁振子能带展现出独特的非相对论性 $f$ 波型奇宇称自旋分裂。
4. 电子带隙闭合迫使的磁振子拓扑相变：伴随着次近邻参数 $t_2$ 的改变，电子能带在 $K/\bar{K}$ 点发生带隙闭合与重开（拓扑向非拓扑 AFM 相变），这必然迫使上方的粒子-空穴连续区向下塌缩，拉动并迫使磁振子带隙在同一点闭合，引发磁振子拓扑数的突变（$C = -2 \to 0$）。这一机制超越了常规局域自旋模型的描述范畴，凸显了巡游电子背景对磁振子拓扑的决定性反馈作用。

下文将面向具有量子化学、关联电子系统和拓扑磁学背景的科研人员，对该工作的核心物理图像、数学公式、数值基准、复现路径及理论局限性进行深度的技术解析。

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

1.1 核心科学问题

常规的反铁磁体（如偶宇称 altermagnets）其磁振子自旋分裂通常表现出偶宇称（如 $d$ 波或 $g$ 波）特征，而奇宇称自旋分裂（如 $p$ 波或 $f$ 波）传统上被认为只存在于非共线磁结构（如螺旋磁体、非共面非共线织构）中。然而，近期研究表明，在存在本征**电子环形电流（Loop Currents）**的共线磁体中，即使不引入相对论性自旋-轨道耦合（SOC），也能实现非相对论性的奇波自旋分裂。

本工作回答的核心科学问题是：

• 在具有自发环形电流的拓扑晶格体系（以 Haldane-Hubbard 模型为代表）中，共线磁有序态下的集体激发态（磁振子）的微观形成机制是什么？
• 顺磁相中的拓扑激子是如何通过凝聚过渡为奇宇称自旋分裂磁振子的？
• 巡游电子的能带拓扑演化（带隙闭合）对这些集体激发的拓扑和色散结构有何深层次的制约与反馈？

1.2 Haldane-Hubbard 模型的微观哈密顿量

考虑二维蜂窝晶格（如石墨烯结构，包含 $A, B$ 双子晶格），Haldane-Hubbard 模型哈密顿量表示为：

其中：

• $t_1$ 是最近邻（NN）实跃迁系数。
• $t_2$ 是次近邻（NNN）复跃迁系数，引入相位因子 $e^{i \nu_{r, r'} \phi}$。本工作设置相位 $\phi = \pi/2$，使得 NNN 跃迁纯为虚数。$\nu_{r, r'} = \pm 1$ 决定了手性（如图 1(a) 所示），在每个六角形斑块（Plaquette）中产生大小为 $6\phi$ 的净磁通，从而打破了时间反演对称性 $\mathcal{T}$，并自发感生内部的电子环形电流。
• $U$ 是位置实空间的 Hubbard 排斥力，控制着体系自旋涨落与局域化程度。本研究侧重半满（Half-filling）情况，重点探讨中等关联强度 $U/t_1 = 4$ 的物理区间。

1.3 理论工具之一：自洽 Mean-Field Hartree-Fock (HF)

为了建立集体激发的研究基准，必须首先自洽求解哈密顿量的平均场基态。将 Hubbard 项进行常规的 Hartree-Fock 平均场解耦：

由于目标相为顺磁（PM）和共线 Néel 反铁磁（AFM）态，忽略含有非共线自旋翻转的 Fock 项。在动量空间中，平均场哈密顿量可写为分块矩阵形式：

其中，基底为 $\hat{\mathbf{c}}_{\mathbf{k}\sigma} = (c_{\mathbf{k}, A, \sigma}, c_{\mathbf{k}, B, \sigma})^{T}$，各项定义为：

• 最近邻色散项：$\gamma(\mathbf{k}) = -t_1 \sum_{j=1}^3 e^{i \mathbf{k} \cdot \mathbf{\delta}_j} = -t_1 (1 + e^{-i \mathbf{k}\cdot\mathbf{a}_1} + e^{-i \mathbf{k}\cdot\mathbf{a}_2})$
• 次近邻拓扑项：$\iota(\mathbf{k}) = -2t_2 \sum_{i=1}^3 \cos(\mathbf{k} \cdot \mathbf{a}_i + \phi)$

通过对该矩阵自洽对角化，获取占有轨道（价带 $v$）与空轨道（导带 $c$）的单粒子本征能量 $\epsilon_{\mathbf{k}, v/c, \sigma}$ 及其 Bloch 波函数 $|\mathbf{k}, v/c, \sigma\rangle$。自洽相图如图 1(c) 所示：当 $t_2/t_1 \lesssim 0.4$ 时，体系随着 $U$ 的增加经历：

• Phase A (拓扑顺磁相)：自旋对称性保持，电子能带具有本征拓扑 Chern 数 $C_{v\uparrow} = C_{v\downarrow} = +1$。
• Phase B (拓扑共线 AFM 相)：自旋对称性破缺产生 Néel 序，但能带仍保持非平凡拓扑（$C_{v\uparrow} = C_{v\downarrow} = +1$）。
• Phase C (普通/非拓扑共线 AFM 相)：自旋对称性破缺，能带拓扑 Chern 数转变为 $C = 0$。

1.4 核心技术：时协 Hartree-Fock (TDHF) 理论与准玻色子近似

为了超越单粒子物理，描述体系中的多体关联束缚态——激子（Excitons）以及磁性相中的玻色子激发——磁振子（Magnons），必须在双粒子格林函数或运动方程层面进行求解。本工作采用 TDHF 理论（在核物理中常称为准粒子随机相位近似 QRPA）。

定义自旋翻转通道（以自旋向下至向上 $\sigma \to \bar{\sigma}$ 为例，这里指 $\uparrow \to \downarrow$）的激子产生算符 $Q^{\dagger}_{\mathbf{q}, \nu, \uparrow \to \downarrow}$。在 Bogoliubov-de Gennes (BdG) 形式下，它由两部分线性组合而成：

其中，$\alpha \in [0, 1]$ 是 Berry 联络规范相关的物理位置分配参数。经过严格数值验证，$\alpha=0$ 与 $\alpha=1$ 在极限下收敛一致，后续计算均采用 $\alpha=1$ 规范。

激子的本征频率 $\omega_{\mathbf{q}, \nu}$ 通过求解如下运动方程获得：

在 TDHF 近似下，该算符方程投影至多体平均场基态 $|0\rangle$，演化为如下非厄米（更准确地说是拟辛/拟埃尔米特，Paraunitary）特征值问题：

矩阵元具体形式由多体微扰论给出：

这里 $\sigma^{\pm}_{\alpha}$ 是作用于蜂窝晶格子格空间 $\alpha = \{A, B\}$ 上的自旋阶梯算符：$\sigma^{\pm}_{\alpha} = \frac{1}{2}(\sigma^x_{\alpha} \pm i\sigma^y_{\alpha})$。

1.5 技术难点及应对策略

• 难点一：拟辛矩阵的数值不稳定性。上式的对角化矩阵形式为 $\mathcal{M} = \sigma_z \mathcal{K}$。由于 $\mathcal{M}$ 是非厄米矩阵，常规的 eigen 求解器极易在激子能量临近零点或体系发生相变时，由于舍入误差产生虚数特征值（物理上意味着动力学不稳定性，但在数值上需严格区分是物理失稳还是数值失稳）。
  • 对策：利用 Colpa 对角化算法。在系统处于稳定相区时，确保双线性形式矩阵 $\mathcal{K} = \begin{pmatrix} A & B \\ B^* & A^* \end{pmatrix}$ 严格正定。执行 Cholesky 分解 $\mathcal{K} = L L^{\dagger}$，将原拟辛矩阵特征值问题转化为一个标准的厄米矩阵特征值求解问题，极大提高了计算稳定性。
• 难点二：大动量网格上的格点效应与奇异性。当计算激子的 Berry 曲率和 Chern 数时，在狄拉克点 $K/\bar{K}$ 和高对称点 $M$ 附近，波函数的规范不确定性会导致传统求导法失效。
  • 对策：引入 Fukui-Hatsugai-Suzuki (FHS) 离散格点 Wilson 圈方法，在离散的 $30 \times 30$ $\mathbf{k}$ 空间网格上定义离散 Berry 联络与流，严格确保导出的激子和磁振子 Chern 数为精确整数。

2. 关键 Benchmark 体系、计算所得数据与物理性能分析

2.1 顺磁拓扑激子相 (Phase A)

设置参数 $U/t_1 = 4.0$，以 $t_2/t_1 = 0.3$ 为典型基准点进行计算。计算所得动态自旋结构因子 $S(\mathbf{q}, \omega) = S_{+-}(\mathbf{q}, \omega) + S_{-+}(\mathbf{q}, \omega)$ 谱图（结合自洽 RPA 与 TDHF）如图 2(a) 所示。

物理分析：激子的 Chern 数不仅可以由多体波函数 Wilson 圈方法直接求得，也可以通过公式 $C_{\text{exc}}^{\sigma\to\bar{\sigma}} = C_{c\bar{\sigma}} - C_{v\sigma}$ 完美验证。在本区，由于 $C_{c\downarrow} = -1$， $C_{v\uparrow} = +1$，因此 $C_{\text{exc}} = -1 - (+1) = -2$。图 3 所示的激子 Berry 曲率分布清晰地表明，曲率的极大值并不分布在电子带隙极小值点 $K/\bar{K}$ 处，而是高度汇聚在三个等价的 $M$ 点。这充分揭示了激子作为粒子-空穴束缚态的多体关联本征效应，并非简单的单粒子能带行为的简单线性组合。

2.2 激子凝聚与向拓扑共线 AFM 相 (Phase B) 的转变

当减小次近邻虚跃迁参数 $t_2/t_1$ 时，电子能隙结构改变。图 2(b) 展示了相变临界点 $t_2/t_1 = 0.215$ 处的激发谱。

• 计算数据点：在 $t_2/t_1 = 0.215$ 时，$\Gamma$ 点的激子能量 $\hbar\omega_{\text{exc}}(\Gamma)$ 骤降至 $\approx 0.08 t_1$。
• 临界行为：随着 $t_2/t_1$ 进一步逼近 $0.212$，$\Gamma$ 点激子能隙完全闭合。激子凝聚算符的期望值 $\langle Q^{\dagger}_{\mathbf{q}=0} \rangle \neq 0$ 建立，标志着拓扑激子发生凝聚。这一凝聚过程自发打破了自旋空间的 $SO(3)$ 旋转对称性，将系统推向具有交错磁化强度的 Néel 反铁磁相。值得注意的是，TDHF 预测的激子凝聚临界线（图 1(c) 绿线）与单粒子自洽 Hartree-Fock 的磁性相变边界完美重合。

2.3 拓扑共线 AFM 相 (Phase B) 中的奇宇称自旋分裂磁振子

在相变后的 Phase B（以 $t_2/t_1 = 0.210$ 为基准点，图 2(c), (d)）中：

在除高对称点外的任意动量 $\mathbf{q}$ 处，磁振子谱发生分裂。具体分裂形态满足三重旋转对称性约束：

$f$ 波自旋分裂数据特征：如图 2(d) 插入图的等能面轮廓线所示，该分裂呈现典型的六角瓣状（$f$ 波对称性）。在长波极限（$\mathbf{q} \to 0$）下，分裂绝对值：

此处 $\theta_{\mathbf{q}}$ 为动量空间极角。这一关系在极低能区保证了 Goldstone 玻色子的线性耗散特征（Type-I 磁振子，声速不分裂），而在中高能区展现强烈的 $f$ 波自旋极化。

2.4 电子带隙闭合引起的磁振子拓扑相变 (Phase B $\to$ Phase C)

本工作的一项重大科学突破在于揭示了巡游电子带隙闭合对磁振子拓扑的强约束机制。调节次近邻参数跨越 $t_2/t_1 = 0.200$：

• 电子能带端（图 4(a, d, g)）：当 $t_2/t_1 = 0.200$ 时，自旋向上与向下的电子能带在 $K/\bar{K}$ 点同时发生狄拉克锥接触（能带闭合）。越过此点（如 $t_2/t_1 = 0.100$），体系进入非拓扑 AFM 相 C（电子 Chern 数变为 $0$）。
• 集体激发端（图 4(b, e, h) & 图 4(c, f, i)）：随着电子能带在 $K/\bar{K}$ 的闭合，由电子-空穴跃迁构成的连续体下边缘必然塌塌至零能。因为磁振子是束缚态，其能级处于连续体下方，这就不可避免地迫使磁振子能带在 $K/\bar{K}$ 点也发生带隙闭合（图 4(f)）。
• 拓扑数演化：
  • Phase B ($t_2/t_1 = 0.210$): 磁振子 Chern 数 $C_{\text{mag}} = -2$
  • 临界点 ($t_2/t_1 = 0.200$): 磁振子能隙在 $K/\bar{K}$ 闭合
  • Phase C ($t_2/t_1 = 0.100$): 磁振子能隙重开，Chern 数突变为 $C_{\text{mag}} = 0$

这一完整的拓扑演化路径由图 4 清晰展现，成功建立了“电子拓扑相变 $\to$ 连续体边缘塌陷 $\to$ 磁振子拓扑相变”这一明确的物理因果链条。

3. 代码实现细节、复现指南与开源工具推荐

为了使科研人员和量子化学计算人员能够高效复现本文的研究成果，下面给出基于自洽平均场自旋流以及 TDHF 拟辛特征值求解的算法架构及复现步骤。

3.1 核心复现算法架构设计

复现的核心分为两大模块：

1. 自洽 Hartree-Fock 模块：获得确定动量格点上的平均场哈密顿量、能级和 Bloch 波函数矩阵。
2. TDHF 拟辛矩阵构造与对角化模块：利用得到的自洽波函数构造大矩阵 $A, B$ 并求解。

模块一：自洽 HF 求解算法（Julia 伪代码伪算法）

using LinearAlgebra

# 参数定义
const t1 = 1.0
const U = 4.0
const phi = pi / 2
const L = 30 # 动量空间 L x L 的网格

# 倒空间格点构建
function build_k_mesh(L)
    k_mesh = []
    # 生成蜂窝晶格对应的第一布里渊区动量点
    # ...
    return k_mesh
end

# 自洽循环主程序
function solve_self_consistent_hf(t2, U, L)
    # 1. 初始猜测子格占有数 n_A_up, n_A_dn, n_B_up, n_B_dn
    n_A_up = 0.8; n_A_dn = 0.2; n_B_up = 0.2; n_B_dn = 0.8
    tol = 1e-8
    diff = 1.0
    
    # 动量格点
    k_mesh = build_k_mesh(L)
    N_k = length(k_mesh)
    
    while diff > tol
        n_A_up_new, n_A_dn_new = 0.0, 0.0
        n_B_up_new, n_B_dn_new = 0.0, 0.0
        
        for k in k_mesh
            # 分别构造自旋向上和向下平均场哈密顿 H_up, H_dn
            # H_up = [ iota(k) + U*n_B_dn,  gamma(k);
            #          gamma*(k),         -iota(k) + U*n_A_dn ]
            # 对角化 H_up 和 H_dn，取出占有态波函数并累加到电荷密度
            # ...
        end
        
        # 归一化并更新占有数，计算 diff
        # ...
    end
    return converged_wavefunctions_and_energies
end

模块二：TDHF 矩阵构造与 Colpa 对角化

对于给定的激发动量 $\mathbf{q}$，构造维度为 $2N_k \times 2N_k$ 的大矩阵（假定系统有 $N_k$ 个单粒子动量点，每个点有价带到导带的单粒子激发）：

function solve_tdhf_exciton(q, hf_data, U, N_k)
    # hf_data 包含所有 k 点的自洽能级 ε(k, v/c, σ) 和波函数 |k, v/c, σ>
    # 1. 初始化 A 和 B 矩阵，维度均为 (N_k, N_k)
    A = zeros(ComplexF64, N_k, N_k)
    B = zeros(ComplexF64, N_k, N_k)
    
    # 2. 填充矩阵元（依据公式 5 和 6）
    # 需要高效计算子格自旋算符的泡利矩阵元：<k, v| σ_α | k+q, c>
    # ...
    
    # 3. 构建大厄米双线性矩阵 K
    # K = [ A    B;
    #       B^H  A^* ]
    K = [A B; B' conj(A)]
    
    # 4. 采用 Colpa 对角化方法（确保 K 正定）
    L_mat = cholesky(Hermitian(K)).L
    
    # 构造辅助厄米矩阵 H_tilde = L^H * J * L ，其中 J = diag(I, -I)
    J = diagonal_matrix_J(N_k)
    H_tilde = L_mat' * J * L_mat
    
    # 对 H_tilde 进行标准厄米对角化
    evals, evecs = eigen(Hermitian(H_tilde))
    
    # 物理激发能量 ω = J * evals
    # 转换波函数 envelope functions (X, Y)
    # ...
    return exciton_energies, exciton_wavefunctions
end

3.2 物理量计算：Chern 数的高效格点计算法

在 FHS 方法中，激子的 Chern 数通过计算第一布里渊区内每个离散网格小方块（由 $k, k+\mu, k+\mu+\nu, k+\nu$ 围成）的 U(1) 规范场磁通量之和给出。对于最低激子态波函数 $|\Psi_{\mathbf{q}, \nu}\rangle$：

1. 在网格方向 $\mu = \mathbf{e}_x, \mathbf{e}_y$ 上计算重叠矩阵元 $U_\mu(\mathbf{q}) = \langle\Psi_{\mathbf{q}, \nu}|\Psi_{\mathbf{q}+\mu, \nu}\rangle / |\langle\Psi_{\mathbf{q}, \nu}|\Psi_{\mathbf{q}+\mu, \nu}\rangle|$。
2. 计算单个小方块的格点曲率：$F_{xy}(\mathbf{q}) = \ln \left[ U_x(\mathbf{q}) U_y(\mathbf{q}+\mathbf{e}_x) U_x(\mathbf{q}+\mathbf{e}_y)^{-1} U_y(\mathbf{q})^{-1} \right]$。
3. 格点共轭 Chern 数即为：$C_{\text{exc}} = \frac{1}{2\pi i} \sum_{\mathbf{q}} F_{xy}(\mathbf{q})$。

3.3 推荐使用的开源工具与物理软件库

• 自洽 Tight-Binding 框架：推荐使用 Julia 的 PyBinding 接口或 Tequila、Kwant 开源库，它们能快速搭建二维六角晶格并自动计算复杂格点跃迁和自洽平均场场量。
• 大规模高精度矩阵对角化：由于对角化涉及到极高维度的稀疏/稠密非厄米特征值问题，Julia 的 IterativeSolvers.jl 或 C++ 的 SLEPc / PETSc 可以提供高效的高对称性特征值求解支持。

4. 关键引用文献与局限性评述

4.1 核心引用文献分析

本工作立足于以下数个关键研究支柱之上：

• [1] F. D. M. Haldane, Phys. Rev. Lett. 61, 2015 (1988): 奠定了无朗道能级的反常霍尔效应以及非平庸拓扑能带的基石。本工作在 NNN 跃迁中完全继承了 Haldane 的虚数跃迁拓扑相因子设计。
• [6] V. S. Arun et al., Phys. Rev. B 93, 115110 (2016): 系统构建了半满 Haldane-Hubbard 模型的 Hartree-Fock 平均场相图，明确指出了 Phase A, B, C 的单粒子物性。本工作的激发态计算正是建立在其基态的基础之上。
• [15, 16] L. Šmejkal et al., Phys. Rev. X 12, 031042 (2022): 确立了 altermagnetism（自旋分裂偶宇称磁体）的物理范式，证明无需强自旋-轨道耦合，仅依靠晶格对称性破缺即可实现非平凡的自旋极化色散。本工作则是将其进一步拓展至“奇宇称共线磁体”与非相对论性磁振子自旋分裂的研究领域。
• [50] Y. H. Kwan et al., Phys. Rev. Lett. 126, 137601 (2021): 开创性地将激子能带拓扑和 Chern 数的概念引入莫尔（Moiré）超晶格体系中，构成了本工作求解激子 Chern 数的基础数学方法学。
• [59] Y. Hwang et al., arXiv:2604.21643: 提出了通过激子波函数零点（Stable Wave-function Zeros）来定义和保护激子拓扑的物理假说。本工作在第 4 页第 2 栏公式 10 附近，完美验证了该假说在拟辛 BdG 激子系统中的精确适用性。

4.2 局限性深度评述

尽管该工作在关联拓扑激发态领域做出了里程碑式的理论发现，但在量子物理与材料学层面上，仍存在如下不可忽视的理论局限性：

1. 平均场（TDHF/RPA）近似对强关联涨落的低估： TDHF 和自洽 RPA 本质上仍属于微扰论在弱-中等关联区（$U/t_1 \le 4$）的推广。在更强关联或存在阻挫的区段，动力学平均场理论（DMFT）或变分蒙特卡洛（VMC）已经证实 Haldane-Hubbard 模型会演化出非共面非共线非平庸磁有序（如四面体磁结构，Tetrahedral Order）甚至手性自旋液体（Chiral Spin Liquid）相。在这些区域，共线 Néel 态的假设不再成立，基于简并粒子-空穴激子凝聚的图像将需要重构。
2. 忽略了强激子-激子散射效应： TDHF 在计算激子凝聚时，采用的是二次准玻色子近似（Quadratic Bosonic Approximation），即将激子算符视为完全尊从对易关系的玻色子。实际上，在临近凝聚点处，激子密度极大，激子间的 Pauli 阻塞效应（Pauli Blocking）以及非简谐激子-激子相互作用会显著改变相变动力学，可能导致二级量子相变演化为一级的物理行为。
3. 实验体系中的非平衡态（Floquet）热化与寿命挑战： 论文作者指出，该模型可利用圆偏振光驱动（Floquet 工程）在石墨烯等蜂窝状材料中实现。但在真实的激光泵浦-探测（Pump-Probe）实验中，高频强激光不可避免地导致强烈的电子雪崩热化（Heating Effect），使系统退相干。在这种强非平衡态背景下，拓扑激子和 $f$ 波磁振子激发的寿命可能极短，如何实现在热化前的超快光谱探测，仍是实验物理学亟待突破的瓶颈。

5. 补充物理直觉与对称性深度解析

为了给科研人员提供更为透彻的物理图像，本节对该体系的对称性保护机制、自旋分裂的物理直觉以及实验探测路径进行补充剖析。

5.1 为什么是 $f$ 波？——基于时空对称性群论的严谨物理解释

在无相对论 SOC 的二维晶格中，电子自旋空间与实空间是完全解耦的。为什么能在自旋共线反铁磁态中，实现高度各向异性的奇宇称 $f$ 波分裂？这可以通过体系的空间-时间联合群操作来深刻理解。

• 时间反演对称性 $\mathcal{T}$：在顺磁相 A 中，$\mathcal{T}$ 已经被虚跃迁（次近邻磁通）显式打破。进入 Néel 反铁磁相 B 后，自旋向上和向下子格交错，自旋反转操作 $\mathcal{T}$ 无法单独保持系统哈密顿量不变。
• 空间反演对称性 $\mathcal{P}$：蜂窝晶格中，空间反演将 $A$ 子格映射至 $B$ 子格。而在 Néel 反铁磁态下，由于 $A, B$ 子格磁化强度相反（$m_A = -m_B$），普通的空间反演 $\mathcal{P}$ 不对称。
• 联合对称性 $\mathcal{P}\mathcal{T}$：由于空间反演交换子格，时间反演翻转自旋，两者的联合操作 $\mathcal{P}\mathcal{T}$ 是哈密顿量的系统对称性。
• 自旋旋转对称性 $C_{2\perp}\mathcal{T}$ 与 $C_{3\parallel}$：在共线磁结构中，垂直于自旋取向（设为 $z$ 轴）的 $x$ 轴或 $y$ 轴方向旋转 $\pi$ 的操作 $C_{2\perp}$ 会翻转 $z$ 方向的自旋。因此，自旋空间旋转 $C_{2\perp}$ 与时间反演 $\mathcal{T}$ 组成的联合对称性 $C_{2\perp}\mathcal{T}$ 被打破。但是，由于蜂窝晶格保留了绕垂直平面轴的 3 阶旋转对称性 $C_{3\parallel}$，磁振子能带必须保持 3 阶旋转对称。

由于 $\mathcal{P}\mathcal{T}$ 对称性的保持，磁振子在动量相反点 $\mathbf{q}$ 和 $-\mathbf{q}$ 处的两支能带满足自旋简并解除下的轨道反对称：

结合 3 阶旋转对称性 $C_{3\parallel}$ 的要求，自旋劈裂函数 $\Delta\omega(\mathbf{q}) = \omega_{\mathbf{q}, +-} - \omega_{\mathbf{q}, -+}$ 必须在动量空间中展现出反对称（$\Delta\omega(\mathbf{q}) = -\Delta\omega(-\mathbf{q})$）且具有 3 倍周期性。在二维动量平面上，最符合这一群论要求的最低阶谐波展开即为：

这正是经典的 $f$ 波物理形式。这一群论分析彻底撇开了具体的模型细节，在对称性层面上宣告了 Haldane-Hubbard 模型反铁磁相中必然存在 $f$ 波磁振子自旋分裂。

5.2 避开衰减（Avoided Decay）——巡游电子层面的能级排斥

在局域自旋模型（如海森堡模型）中，磁振子是无法感知到巡游电子粒子-空穴连续体的。而在真实的巡游电子体系（Haldane-Hubbard 模型）中，TDHF 的计算显示了超越局域模型的现象：避开衰减（Avoided Decay）机制。

在 $K/\bar{K}$ 附近，磁振子能谱极度逼近上方的粒子-空穴连续体边缘（见图 4(e, h)）。在强多体排斥力 $U$ 产生的能级排斥（Level Repulsion）作用下，磁振子束缚态能级被强行“压低”，使其在不掉入连续体内部（一旦掉入则发生朗道衰减）的前提下，最大程度地弯曲色散。这种能级避开效应使得磁振子的自旋分裂差值相比于简单局域自旋模型的预测结果，被放大了数倍之多，体现了巡游性对奇波自旋劈裂的极强放大作用。这一结论对寻找具有大自旋分裂的磁振子器件具有重要的指导意义。

5.3 实验观测方法学建议

如何在真实或人工量子系统中探测该工作预言的物理效应？本指南给出以下具体观测路径：

+--------------------------------------------------------------+
|                       实验观测路径设计                        |
+--------------------------------------------------------------+
                               |
      +------------------------+------------------------+
      |                                                 |
      v                                                 v
[方案一：固态量子材料层]                        [方案二：超冷原子气人工模拟]
- 候选材料：双层石墨烯、过渡金属硫族化合物          - 晶格：光晶格（Optical Lattice）人工构建
- 机制：Floquet 偏振激光驱动感生虚跃迁              - 环形电流：通过拉曼辅助跃迁（Raman Assisted
- 探测 1：偏振非弹性中子散射（INS）                  Hopping）调控次近邻相位
- 探测 2：共振非弹性 X 射线散射（RIXS）            - 探测：利用原位量子气体显微镜（Quantum
- 信号：观测动量各向异性的磁振子谱分裂               Gas Microscope）配合布拉格光谱仪直接测
                                                     量拓扑激子在动量空间的相干凝聚信号

通过上述多维度的理论解析、数值 benchmark 及实验复现方案，本工作不仅为拓扑凝聚态物理开辟了“利用拓扑激子凝聚创制新型奇宇称自旋分裂磁体”的崭新道路，也为量子化学界在更广泛的分子-晶格关联激发态物理计算中，提供了极佳的研究范例。