来源论文: https://arxiv.org/abs/2606.28521v1 生成时间: Jul 04, 2026 12:17

磁场中交替混合自旋 (1/2, 1) 海森堡梯子的量子相图与相变：基于 DMRG 与线性自旋波理论的深度解析

0. 执行摘要

低维量子磁性系统（Low-dimensional quantum magnetic systems）由于强烈的量子涨落和丰富的拓扑激发，一直是凝聚态物理和量子化学领域的研究热点。在这些系统中，磁化高原（Magnetization Plateaus）的出现不仅标志着具有有限能隙的非平凡量子相的形成，还对应着独特的拓扑序和对称性破缺机制。特别是交替混合自旋（Mixed-spin）链或梯子系统，由于兼具不同自旋大小（如 $s=1/2$ 与 $S=1$）的局部自由度，展现出远比同同自旋系统（Homogeneous-spin systems）更为复杂的相图。

本文针对 D. S. Almeida 等人的最新研究工作《Mixed-spin Heisenberg ladders in a magnetic field》进行深度技术解析。该工作系统地研究了处于外加磁场 $h$ 下的交替混合自旋 $(s, S) = (1/2, 1)$ 海森堡梯子。作者联合运用了高精度密度矩阵重整化群（Density Matrix Renormalization Group, DMRG）方法和基于完全极化真空（Fully Polarized Vacuum）的线性自旋波理论（Linear Spin-Wave Theory, LSWT），绘制了完整的磁场 $h$ 对比层间（rung）耦合强度 $J_\perp$ 的量子相图。

这项工作最显著的物理发现和技术突破包括：

1. 对称性简化与非平庸相图：通过发现和利用系统的滑移反射对称性（Glide Reflection Symmetry），成功地将具有 4 个自旋的原始胞（4-spin unit cell）简化为仅含 2 个自旋的简化胞（2-spin unit cell），大大简化了解析解析和自旋波对角化过程，从而在展开布里渊区（Unfolded Brillouin Zone, UBZ）中计算出了精确的激发谱。
2. 1/3 磁化高原的稳定性区间：发现并证实了在反铁磁层间耦合（$J_\perp > 0$）以及有限的铁磁层间耦合（$-1.32 < J_\perp < 0$）区域内，系统存在一个极其稳定的 $1/3$ 磁化高原（1/3-plateau）。这一发现严格满足了低维磁性中的 Oshikawa-Yamanaka-Affleck (OYA) 判据。
3. Kosterlitz-Thouless (KT) 相变点的精准确定：在铁磁层间耦合区域，随着 $J_\perp$ 的减小，$1/3$ 高原的能隙逐渐闭合。作者通过计算不同系统尺寸（最高至 $L=256$ 聚体）下的横向自旋关联函数，并进行有限尺寸标度（Finite-Size Scaling）分析，确定了闭合点处的 Luttinger 液体参数 $K = 2$，从而极其精准地定位了 KT 临界点，得出 $J_{\perp, KT} = -1.32 \pm 0.02$。
4. 磁结构的多样性：揭示了系统从非重构的铁磁/反铁磁关联态到由磁场驱动的 gapless Luttinger 液体相的转变过程，为分子磁体（Molecular Magnets）的实验合成与表征提供了坚实的理论基础。

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

1.1 核心科学问题与 OYA 判据

在强关联电子系统和低维量子磁性中，最基本的物理效应之一就是外加磁场下的磁化高原。对于一个一维自旋链，若其磁化强度 $m$（每个单元的平均磁化）在某一磁场范围内保持恒定，这意味着系统在该状态下存在一个有限的电荷或自旋能隙。这种量子磁化高原的存在必须满足由 Oshikawa, Yamanaka 和 Affleck 提出的 OYA 判据：

其中 $S_{\text{cell}}$ 是系统磁性基态最小重复单元（晶胞）内的总自旋大小，$p$ 是物理磁性基态空间调制的周期（以原始晶胞为单位）。 对于本文所研究的 $(s, S) = (1/2, 1)$ 交替梯子系统，若采用图 1(a) 所示的 4 自旋胞（包含两个 $s=1/2$ 和两个 $S=1$ 自旋），其 $S_{\text{cell}} = 1/2 + 1/2 + 1 + 1 = 3$。当周期 $p=1$（即空间平移对称性不发生自发破缺）时，为了形成磁化高原，其磁化强度 $m$ 必须满足：

然而，如果引入系统的滑移反射对称性（图 1(b)），我们可以将晶胞缩小为仅含一个 $s=1/2$ 和一个 $S=1$ 的 2 自旋胞，此时 $S_{\text{cell}} = 1/2 + 1 = 3/2$。在 $p=1$ 的情况下，OYA 判据变为：

这里的 $m=1/2$ 刚好对应于饱和磁化强度 $m_s = 3/2$ 的 $1/3$（即 $1/3$-plateau）。因此，滑移反射对称性的存在构成了该系统展现稳定的 $1/3$ 磁化高原的拓扑与对称性基础。本工作致力于解决的核心科学问题是：这一 $1/3$ 高原在强烈的量子涨落和竞争性的层间耦合 $J_\perp$ 下如何保持稳定？它在什么临界条件下会通过量子相变闭合？

1.2 理论模型与滑移反射对称性分析

交替混合自旋海森堡梯子的哈密顿量可以写为两种等价形式。

1.2.1 原始 4 自旋胞表象（Folded Brillouin Zone, FBZ）

在 4 自旋胞表示下（如图 1(a)），系统的哈密顿量为：

其中 $j$ 是晶胞索引，$N_c$ 是总晶胞数。上标 $(1), (2)$ 分别代表梯子的上、下两条腿。$\mathbf{s}_j^{(\alpha)}$ 为自旋-1/2 算符，$\mathbf{S}_j^{(\alpha)}$ 为自旋-1 算符。$J$ 设定为能量单位（$J=1$），表示梯子腿方向的内层耦合；$J_\perp$ 表示两条腿之间的耦合。$S^z$ 是系统沿着磁场方向的总自旋投影：

1.2.2 2 自旋胞表象（Unfolded Brillouin Zone, UBZ）

显然，哈密顿量 (1) 在以下操作下是不变的：将晶格沿腿方向平移距离 $u$（这会将上腿的 $s$ 变换到下腿的 $S$ 的水平位置），然后交换上下腿的标记 $1 \leftrightarrow 2$。这种复合对称性被称为滑移反射对称性。

通过对自旋索引进行重新排列（如图 1(b) 所示），将奇数二聚体和偶数二聚体交叉放置，我们可以构造一个周期为 $u$ 的等价单链哈密顿量：

其中 $L = 2N_c$ 是系统中二聚体 $(\mathbf{s}_l, \mathbf{S}_l)$ 的总数。这种表示法将原先含有 4 个磁自旋的物理晶胞减少为仅包含 2 个自旋（一个自旋-1/2 和一个自旋-1）的简化胞，从而在自旋波计算中将 4 个色散支（折叠布里渊区）简化为 2 个色散支（展开布里渊区，如图 1(c)），大幅降低了解析和数值计算的复杂性。

1.3 线性自旋波理论（LSWT）的数学推导与技术细节

为了确定系统进入完全极化（Fully Polarized, FP）状态的精确临界磁场 $h_{FP}$，我们可以将完全极化态作为经典真空，并引入 Holstein-Primakoff (HP) 变换。

对于自旋 $s = 1/2$ 场源 $\mathbf{s}_l$ 和自旋 $S = 1$ 场源 $\mathbf{S}_l$，其 $z$ 分量写为：

其中 $a_l (a_l^\dagger)$ 和 $b_l (b_l^\dagger)$ 均为玻色子消灭（产生）算符。其阶梯算符在最低阶近似（线性近似）下为：

将这些算符代入简化哈密顿量 (2) 中，忽略常数项及高于二次的玻色子相互作用项，得到二次型自旋波哈密顿量。通过傅里叶变换：

我们得到动量空间中的二次型哈密顿量：

其中耦合常数与局部自旋势能分别为：

通过对系数矩阵进行标准的对角化（由于哈密顿量在 HP 变换下不含双产生/双消灭项，只需普通幺正矩阵对角化即可），可以求得两个自由自旋波激发的色散关系：

将 $(s, S) = (1/2, 1)$ 代入，最终色散关系化简为：

在外加磁场 $h$ 存在时，磁激发的单粒子谱线被整体移频：$\omega_h^{(\pm)}(k) = \omega^{(\pm)}(k) + h$。完全极化态保持稳定的条件是最低能激发的自旋波能量必须非负，即 $\min_k [\omega_h^{(-)}(k)] \ge 0$。根据系统的色散关系，极小值点的位置取决于 $J_\perp$ 的符号：

1. 对于 $J_\perp > 0$，极小值发生在布里渊区中心 $k = 0$ 处，此时临界磁场为：

   $$h_{FP}(J_\perp) = -\omega^{(-)}(0) = \frac{3}{2}(J_\perp + 2) \quad (13)$$

2. 对于 $J_\perp < 0$，极小值发生在布里渊区边界 $k = \pi/u$ 处，此时临界磁场为：

   $$h_{FP}(J_\perp) = -\omega^{(-)}\left(\frac{\pi}{u}\right) = \frac{3}{4}(J_\perp + 2) + \frac{1}{4}\sqrt{(J_\perp + 2)^2 + 8(J_\perp - 2)^2} \quad (14)$$

这两条解析临界线与 DMRG 数值计算所得结果在外推至热力学极限时展现了完美的符合（误差在 $10^{-5}$ 量级以内）。

1.4 数值重整化群（DMRG）与有限尺寸标度分析的技术难点

尽管线性自旋波理论能极其精准地预测完全极化区（饱和区）的临界边界，但由于其忽略了玻色子之间的相互作用（即 magnon-magnon interactions），它无法定量描述低磁场下的 gapless 液体相，更无法准确预测 $1/3$ 磁化高原的下边界 $h_-$ 及其在铁磁耦合区的 KT 闭合点。为此，必须引入高精度的 DMRG 数值方法。

技术难点一：开边界条件（OBC）引起的边界效应

DMRG 计算通常在开边界条件（Open Boundary Conditions, OBC）下执行，因为 OBC 下的矩阵乘积态（MPS）变分优化效率远高于周期边界条件（PBC）。然而，OBC 会在链的两端引入显著的边界效应（Boundary Effects），导致局部磁化强度 $\langle s_i^z \rangle$ 和关联函数在空间上产生强烈的非均匀振荡。为了提取真实的热力学极限物理量，必须采用系统尺寸标度（Finite-Size Scaling）方法。本工作通过计算 $L=100, 128, 192, 256$ 聚体等一系列尺寸，并使用多项式和指数外推技术来消除边界效应影响。

技术难点二：KT 相变点附近的能隙指数闭合

Kosterlitz-Thouless (KT) 相变是一个无限阶相变。在相变点附近，激发能隙（Gap）以指数形式闭合：$\Delta \sim \exp(-C/\sqrt{J_\perp - J_{\perp, KT}})$。由于这种极慢的闭合特性，直接通过能量能隙（能级差）的外推来定位 KT 临界点会产生极大的数值误差。为此，必须绕开能隙计算，转而通过高精度的**横向自旋关联函数（Transverse Spin Correlation Functions）**的长程渐近行为来确定 Luttinger 液体参数 $K$。而在无隙 Luttinger 液体中，DMRG 的截断误差（Truncation Error）和纠缠熵（Entanglement Entropy）会随着系统尺寸对数发散，这对 MPO 的键维（Bond Dimension）提出了极高要求。本研究将最大保留键维设为 $m=700$，截断权重（Discarded Weight）严格控制在 $10^{-10}$ 以下，确保了关联函数计算的保真度。

2. 关键 Benchmark 体系、计算所得数据与物理图像解析

为了建立清晰的物理图像并定量验证理论，本工作选择交替混合自旋 $(1/2, 1)$ 海森堡梯子作为 Benchmark 体系，系统地改变层间耦合 $J_\perp \in [-1.5, 1.0]$ 和磁场 $h \in [0.0, 4.5]$。

2.1 整体量子相图（磁场 $h$ 对比层间耦合 $J_\perp$）

DMRG 模拟所得的完整相图绘制于图 2 中。相图展现出四个截然不同的量子相：

1. 完全极化（Fully Polarized, FP）高原：位于高磁场区。系统所有的自旋都沿磁场方向排列，总磁化强度达到饱和值 $m = 1.5$（每个二聚体）。
2. 1/3 磁化高原区：在 $J_\perp > 0$（反铁磁耦合）的整个区域以及 $-1.32 < J_\perp < 0$ 的铁磁耦合区域内，系统展示了一个宽阔且稳定的 $1/3$ 高原，磁化强度平盘保持在 $m = 1/2$。该区域是一个自旋能隙相，基态具有非平庸的自旋排布。
3. 无隙 Luttinger 液体相 (Gapless Phases)：夹在 FP 高原和 1/3 高原之间，以及 1/3 高原与零磁场（或低磁场）基态之间。在这两个无隙相中，自旋波激发没有能隙，横向自旋关联函数表现为幂律衰减（Power-law decay）。
4. KT 相变点（Kosterlitz-Thouless Transition Point）：位于 $J_{\perp, KT} = -1.32 \pm 0.02, h_{KT} = 1.02$。在此临界点，1/3 磁化高原的能隙完全闭合，系统直接从高磁场 gapless 相过渡到低磁场 gapless 相。

2.2 局部磁化强度与磁有序分析

通过对 $J_\perp = 0.25$（反铁磁层间耦合）和 $J_\perp = -0.25$（铁磁层间耦合）进行磁场 $h$ 扫描（h-scans），计算得到的子格平均磁化强度如图 4(a, b) 所示：

物理图像分析

• 反铁磁耦合下 ($J_\perp = 0.25$)：在 1/3 高原处，由于反铁磁性的层间相互作用，自旋-1/2 的 $\langle s^z \rangle$ 被强烈地压制为负值（约为 $-0.19$），与顺磁场排列的自旋-1（$\langle S^z \rangle \approx 0.69$）呈亚铁磁共线排列（Ferrimagnetic alignment）。
• 铁磁耦合下 ($J_\perp = -0.25$)：即便层间相互作用是铁磁性的（倾向于让 $s$ 和 $S$ 平行排列），在 1/3 高原处，量子涨落和内层反铁磁耦合 $J$ 依然强迫自旋-1/2 的 $\langle s^z \rangle$ 保持负值（约为 $-0.12$），这也解释了为什么 1/3 高原在一定强度的铁磁耦合范围内仍能稳定存在。

2.3 横向自旋关联函数与 KT 临界行为的定量标度

为了精确测定 KT 相变点，作者深入研究了沿着梯子腿方向的横向自旋关联函数：

在无隙（Gapless）Luttinger 液体相中，该关联函数在远距离下遵循幂律衰减：

其中 $K$ 即为 Luttinger 参数。在 KT 相变点，1/3 高原对应的能隙闭合，根据 Berezinskii-Kosterlitz-Thouless 理论，相变点处的 Luttinger 参数具有普适值：

作者计算了尺寸为 $L = 128, 192, 256$ 在不同 $J_\perp$ 下的 $\Gamma(r)$，并选取三个空间窗口区间进行数据拟合以消除有限尺寸效应（图 7(c, d)）：

• 拟合区间 1: $r \in [8, 16]$
• 拟合区间 2: $r \in [16, 32]$
• 拟合区间 3: $r \in [32, 48]$

对每个尺寸 $L$ 的拟合值 $K(L)$ 绘制成关于 $1/L$ 的函数，并线性外推至 $1/L \rightarrow 0$（无边界效应的热力学极限）。最终得到热力学极限下的 $K$ 随 $J_\perp$ 的变化曲线（图 7(e)）：

• 当 $J_\perp = -1.4$ 时，外推得到 $K \approx 1.7 < 2$（系统处于有隙的高原相，能隙未闭合）。
• 当 $J_\perp = -1.2$ 时，外推得到 $K \approx 2.15 > 2$（系统进入了无隙相）。
• 曲线与 $K = 2$ 的交点精确定位在：

这一结果代表了目前对该类低维混合自旋系统最精准的临界相变数值界定。

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

本节提供一套基于 Julia 语言和著名张量网络开源库 ITensors.jl 的完整复现方案，用于计算该交替 $(1/2, 1)$ 海森堡链的基态磁化强度和关联函数。

3.1 物理建模与 ITensors 自定义 SiteType 映射

由于哈密顿量 (2) 是通过交替二聚体 $(\mathbf{s}_l, \mathbf{S}_l)$ 构造的单链，物理系统中的奇数格点（odd sites）全为自旋-1/2，偶数格点（even sites）全为自旋-1。在 ITensors 中，我们可以通过显式构建一个非均匀的 SiteSet 来实现这一混合自旋结构。

3.2 完整 Julia 复现代码

using ITensors
using Printf

"""
构建交替混合自旋 (1/2, 1) 哈密顿量的 MPO
N: 总格点数 (必须为偶数，包含 N/2 个 SpinHalf 和 N/2 个 SpinOne)
J_leg: 腿方向内层耦合 (对应公式2中的1.0)
J_rung: 层间耦合 (对应公式2中的 J_perp)
h_field: 外加磁场
"""
function build_mixed_spin_hamiltonian(sites, J_leg::Float64, J_rung::Float64, h_field::Float64)
    N = length(sites)
    os = OpSum()
    
    for j in 1:N
        # 1. 磁场项 -h * S^z_j
        os += -h_field, "Sz", j
        
        # 2. 相互作用项
        if j < N
            if j % 2 == 1
                # 奇数格点到偶数格点：层间耦合 J_perp * (s_j . S_{j+1})
                os += J_rung, "Sz", j, "Sz", j+1
                os += J_rung * 0.5, "S+", j, "S-", j+1
                os += J_rung * 0.5, "S-", j, "S+", j+1
            else
                # 偶数格点到下一个奇数格点：内层腿耦合 J * (S_j . s_{j+1})
                os += J_leg, "Sz", j, "Sz", j+1
                os += J_leg * 0.5, "S+", j, "S-", j+1
                os += J_leg * 0.5, "S-", j, "S+", j+1
            end
        end
    end
    
    return MPO(os, sites)
end

function run_dmrg_simulation(L_dimers::Int, J_perp::Float64, h::Float64)
    # 总物理格点数 N = 2 * L_dimers
    N = 2 * L_dimers
    
    # 创建混合自旋 SiteSet
    # 奇数位置是 S=1/2, 偶数位置是 S=1
    sites_array = Vector{Index}(undef, N)
    for j in 1:N
        if j % 2 == 1
            sites_array[j] = siteind("S=1/2"; conserve_qns=true)
        else
            sites_array[j] = siteind("S=1"; conserve_qns=true)
        end
    end
    sites = SiteSet(sites_array)
    
    # 构建哈密顿量
    H = build_mixed_spin_hamiltonian(sites, 1.0, J_perp, h)
    
    # 设定初始状态为交替量子态，确保总 Sz 为我们期望的扇区
    # 为了探究 1/3 高原，此处可设定总 Sz 量子数
    state = [isodd(n) ? "Dn" : "Up" for n in 1:N]
    psi0 = randomMPS(sites, state)
    
    # 设定 DMRG 扫描参数 (根据论文要求，最大键维为 700，截断阈值为 1e-10)
    sweeps = Sweeps(8)
    setmaxdim!(sweeps, 50, 100, 200, 400, 700, 700)
    setcutoff!(sweeps, 1e-10)
    setnoise!(sweeps, 1e-6, 1e-7, 1e-8, 1e-10, 0.0)
    
    # 执行 DMRG 变分计算
    energy, psi = dmrg(H, psi0, sweeps; silent=true)
    
    # 计算子格磁化强度
    sz_values = expect(psi, "Sz")
    
    # 区分 spin-1/2 和 spin-1 的局部平均磁化
    mag_half = mean(sz_values[1:2:end])
    mag_one  = mean(sz_values[2:2:end])
    total_m  = mag_half + mag_one
    
    @printf("h = %.3f | Energy = %.8f | <s^z> = %.5f | <S^z> = %.5f | m_dimer = %.5f\n", 
            h, energy, mag_half, mag_one, total_m)
            
    return psi, sz_values, total_m
end

# 示例运行：计算 J_perp = 0.25, h = 1.25（处于 1/3 高原区）时的基态
L_dimers = 64 # 对应 128 个自旋格点
psi, sz_local, m_dimer = run_dmrg_simulation(L_dimers, 0.25, 1.25)

3.3 数据复现步骤与性能优化建议

1. 安装环境：确保安装了 Julia 1.8+ 以及最新版的 ITensors 库（using Pkg; Pkg.add("ITensors")）。
2. 多线程加速：ITensors 底层的张量收缩高度依赖 BLAS 库。在运行脚本前，可以通过设置环境变量 export JULIA_NUM_THREADS=4 和 using LinearAlgebra; BLAS.set_num_threads(4) 启用多线程收缩，大幅降低大尺寸（$L=256$）下的计算延迟。
3. 量子数守恒（QN Conservation）：在上面的代码中，我们显式开启了 conserve_qns=true。这是极为关键的性能优化步骤。由于哈密顿量守恒总 $S^z$ 算符，开启 QN 守恒可以将密集张量收缩转化为分块对角矩阵收缩，不仅将计算速度提升 5-10 倍，还能确保 DMRG 严格收敛在特定的磁化扇区内。

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

4.1 关键引用文献

本项研究建立在低维磁性领域几十年来数个里程碑式工作的基础之上：

1. Oshikawa, Yamanaka, & Affleck (1997) [Phys. Rev. Lett. 78, 1984]：提出了低维自旋系统产生磁化高原的普适拓扑判据（OYA 判据），这是本文物理图像的理论基石。
2. Lieb & Mattis (1962) [J. Math. Phys. 3, 749]：证明了反铁磁交替磁性链中基态自旋有序性定理，为解释 $J_\perp > 0$ 时 1/3 高原的亚铁磁共线排列提供了坚实的数学证明。
3. Steven R. White (1992) [Phys. Rev. Lett. 69, 2863]：发明了密度矩阵重整化群（DMRG）算法，使一维强关联体系的精确数值求解成为现实。
4. Montenegro-Filho et al. (2020) [Phys. Rev. B 102, 035137]：详细讨论了一维开链中多体边缘态与 1/3 磁化高原、Luttinger 液体的拓扑关联，为本研究提供了直接的研究方法参考。

4.2 对这项工作的局限性与改进空间的学术评论

尽管该工作在相图精度和临界点的测定上达到了极高的水准，但从现代强关联物理和实验凝聚态化学的角度来看，仍存在以下局限性：

1. 经典 LSWT 理论无法自洽处理强量子涨落下的多 magnon 相互作用

论文中第四节尝试使用自由玻色子气体近似（硬核玻色子填充）来解释 1/3 高原的产生。然而，线性自旋波理论（LSWT）本质上是将自旋算符大 $S$ 展开的最低阶近似，在接近无隙区或强量子起伏区（如 $J_\perp < 0$）时，LSWT 给出的高原边界与 DMRG 精确结果存在显著偏差（LSWT 预测的闭合点在 $J_\perp = -2.0, h=0.0$，而精确 DMRG 发现闭合在 $J_{\perp} = -1.32, h=1.02$）。这表明该体系中的**多 magnon 相互作用和束缚态（Bound States）**极其显著，单粒子自旋波图像在定量层面上是失效的。未来的理论改进需要引入高阶的相互作用自旋波理论（Interacting Spin-Wave Theory）或者动力学平均场方法。

2. 有限尺寸外推中对指数闭合规律的拟合敏感性

在第六节确定 KT 转变点时，作者选用了 $K(L)$ 对 $1/L$ 的多项式线性外推。然而，KT 相变点附近的能隙闭合极其缓慢，伴随着强烈的对角修正项（logarithmic corrections）。在有限的尺寸范围（最高 $L=256$ 聚体，即 $512$ 个自旋格点）内，这些对数修正项可能尚未完全衰减，导致外推得到的 $K(L \rightarrow \infty)$ 存在系统性的微小偏差。采用**无限系统尺寸 DMRG（iDMRG）**直接在热力学极限下计算纠缠熵的标度行为（Entanglement Entropy Scaling），或许能够提供不依赖于人工选取拟合窗口的、更具说服力的临界点验证。

3. 缺乏对实验真实材料中“三维耦合”和“晶格畸变”的讨论

现实中能模拟该模型的物理材料（如过渡金属配位聚合物或双金属链 $[\text{MnCu}(\text{pba})(\text{H}_2\text{O})_3] \cdot 2\text{H}_2\text{O}$）不可避免地存在梯子之间的弱三维反铁磁耦合以及由自旋-声子耦合引起的 Peierls 晶格畸变（Spin-Peierls Distortion）。晶格畸变会自发破缺空间平移对称性，进而极大地改变磁化高原的稳定区间。未来的理论研究应当将系统与一维声子晶格的耦合纳入哈密顿量中。

5. 补充理论：线性自旋波对角化的详细解析推导

为了帮助读者彻底掌握线性自旋波理论的技术精髓，本节详细给出哈密顿量矩阵 $\tau_k$ 的对角化过程。这是论文中 Appendix A 的核心数学内容，但在论文中被高度浓缩。以下给出完整的中间推导步骤。

在动量空间中，系统的无场二次型哈密顿量可以表示为：

其中，玻色子基底向量定义为 $\Psi_k = \begin{pmatrix} a_k \\ b_k \end{pmatrix}$，$\Psi_k^\dagger = \begin{pmatrix} a_k^\dagger & b_k^\dagger \end{pmatrix}$。对应的 $2 \times 2$ 系数矩阵 $\mathbf{M}_k$ 为：

其中，由公式 (9)-(10) 可知：

为了求得矩阵 $\mathbf{M}_k$ 的本征值，我们解其特征方程：

整理展开得到关于 $\omega$ 的一元二次方程：

利用求根公式，直接解得：

我们将具体的物理参数值代入上式中的各项：

1. 均值项：

   $$\frac{\epsilon_a + \epsilon_b}{2} = \frac{-(J_\perp + 2) - \frac{1}{2}(J_\perp + 2)}{2} = -\frac{3}{4}(J_\perp + 2)$$

2. 差值项的平方：

   $$\left(\frac{\epsilon_a - \epsilon_b}{2}\right)^2 = \left( \frac{-(J_\perp + 2) + \frac{1}{2}(J_\perp + 2)}{2} \right)^2 = \left( -\frac{1}{4}(J_\perp + 2) \right)^2 = \frac{1}{16}(J_\perp + 2)^2$$

3. 交叉项的平方：

   $$t_{kk}^2 = \left( \frac{\sqrt{2}}{2}(J_\perp + 2 \cos ku) \right)^2 = \frac{1}{2}(J_\perp + 2 \cos ku)^2$$

将上述 1, 2, 3 项代入特征值通式，得到：

为了与论文中的标准表达一致，将根号内的分母提出来：

这便是论文中第四节公式 (11)-(12) 的完整、严谨的推导过程。自旋波谱的下支 $\omega^{(-)}(k)$ 在外加磁场 $h$ 作用下上移。通过分析 $\omega^{(-)}(k)$ 动量空间极小值的位置，我们可以精确无误地描绘出完全极化态的相边界，这充分体现了线性自旋波理论在研究饱和磁场动力学时的强大威力。