量子自旋系统的桥梁：S=1与3/2 Shastry-Sutherland模型中的螺旋与混合斑块-二聚体相深度解析

来源论文: https://arxiv.org/abs/2605.24409v1 生成时间: May 31, 2026 18:15

0. 执行摘要

在强关联电子系统与量子磁学研究中，二维几何挫折晶格上的海森堡模型因其丰富的量子相变、非平凡的拓扑激发以及潜在的量子自旋液体（QSL）基态而备受瞩目。作为其中最具代表性的范例之一，**Shastry-Sutherland 晶格（SSL）**由于其精确的可解二聚体基态以及在实际材料（如 $\text{SrCu}_2(\text{BO}_3)_2$）中的完美实现，成为凝聚态物理学家研究量子临界性、多体效应与挫折物理的天然实验室。

然而，绝大多数传统研究都聚焦于极端的量子极限 $S = 1/2$ 或者是经典的 $S \to \infty$ 极限，而对于介于两者之间的较大自旋体系（如 $S = 1$ 和 $S = 3/2$），其全局相图、量子相变的演化机制以及量子波动在何种程度上被抑制，长期以来缺乏精确的定量研究。随着近期一系列新型稀土基 Shastry-Sutherland 材料（如 $\text{RE}_2\text{Be}_2\text{GeO}_7$ 系列、$\text{ErB}_4$ 以及 $\text{Eu}_2\text{MgSi}_2\text{O}_7$）的成功合成，探索大自旋 SSL 系统的基态相图已成为当务之急。

本篇博文将深度解析发表于 2026 年的重磅研究——《Spiral and Mixed Plaquette-Dimer Phases in the $S = 1$ and $3/2$ Shastry-Sutherland Heisenberg Model》。该工作由中山大学中子科学与技术学术带头人吴汉青（Han-Qing Wu）、香港中文大学吴木玮（Muwei Wu）以及大湾区大学龚寿书（Shou-Shu Gong）合作完成。作者们采用高度优化的自旋 $SU(2)$ 对称性密度矩阵重正化群（DMRG）方法以及集群平均场理论（CMFT），首次系统地确立了 $S=1$ 和 $S=3/2$ SSL 各向同性海森堡模型的基态相图。他们不仅在两个大自旋系统中均发现了介于二聚体相与 Néel 反铁磁（NAF）相之间的两个全新中间相：混合斑块-二聚体（Mixed Plaquette-Dimer, MPD）相与非共度螺旋（Spiral）相，更通过精细的数据分析，成功构建了跨越量子到经典极限的全局 $S-g$ 相图，完美填补了这一科学空白。本博文将为量子化学与凝聚态多体物理背景的科研工作者提供深度技术拆解与复现指南。

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

1.1 核心科学问题：量子与经典之间的缺失环节

Shastry-Sutherland 晶格（SSL）的几何结构是在正方形晶格上引入交替的对角线耦合 $J'$，其哈密顿量可表示为：

其中，$\langle i,j \rangle$ 表示最近邻（inter-dimer）耦合（耦合强度为 $J$），而 $\langle\langle i,j \rangle\rangle'$ 表示交替对角线的次近邻（intra-dimer）耦合（耦合强度为 $J'$）。我们定义无量纲耦合比为 $g = J/J'$。在能量单位设定中，令对角耦合 $J' = 1$。

对于 $S = 1/2$ 极限，当 $g \to 0$ 时，基态是精确的二聚体单态乘积相；而当 $g \to \infty$ 时，系统转化为标准的正方晶格反铁磁体，呈现 Néel 反铁磁（NAF）有序。在中间区域，已有多项数值研究表明存在空斑块（Empty-Plaquette, EP）相以及可能的量子自旋液体（QSL）相。然而，在经典极限 $S \to \infty$ 中，由于量子涨落完全消失，中间的空斑块相等具有量子特性的相不复存在，取而代之的是在 $0 < g < 1$ 范围内的非共度螺旋（Spiral）相，并在 $g_c = 1$ 处直接相变至 NAF 相。

因此，核心科学问题在于：从极端的量子极限 $S=1/2$ 演化到经典极限 $S \to \infty$ 的过程中，中间相是如何随自旋量子数 $S$ 的增加而演变的？大自旋体系中的强量子波动与几何挫折如何共同作用，稳定新型的量子或半经典相？

1.2 技术难点：希尔伯特空间的指数暴涨与非共度调制

在数值多体物理中，研究 $S > 1/2$ 的二维挫折系统面临两大极其棘手的挑战：

1. 希尔伯特空间的急剧膨胀：对于格点数为 $N$ 的系统，其希尔伯特空间维度为 $(2S + 1)^N$。当自旋从 $S=1/2$ 提升到 $S=1$ 和 $S=3/2$ 时，单格点物理维度由 $d=2$ 分别膨胀至 $d=3$ 和 $d=4$。这意味着在 DMRG 扫描过程中，保持相同截断误差所需的键维度（Bond Dimension）呈指数级增长，普通的 $U(1)$ 对称性（即仅守恒总 $S^z$）DMRG 极易遭遇算力瓶颈。
2. 非共度螺旋相的尺寸效应：螺旋相在实空间具有非共度的调制周期。在圆柱体（Cylinder）几何上进行有限尺寸计算时，若圆柱的长 $L_x$ 或宽 $L_y$ 与螺旋的真实周期不匹配，系统会产生巨大的晶格应力，导致虚假的锁定相变，甚至掩盖真实的物相。因此，必须使用极长的圆柱体系（如 $L_x = 24$）并结合精确的局域物理量测量。

1.3 方法细节：SU(2) DMRG 与集群平均场理论（CMFT）

为了克服上述难点，作者们结合了当前最先进的数值计算武器库：

1.3.1 自旋 $SU(2)$ 对称性 DMRG

本研究的核心计算基于严格保对称性的 DMRG 框架。通过在算法中完全显式地群论化自旋旋转对称性 $SU(2)$，可将哈密顿矩阵块对角化为由总自旋量子数 $J$ 标记的子空间。这允许在计算中保留高达 $8000$ 个 $SU(2)$ 等效状态（相当于无对称性状态下的数万个状态），从而将最大截断误差控制在极低的 $5 \times 10^{-6}$ 以下。计算主要在 $L_y = 4, 6$ 的圆柱体几何上展开，其边界条件沿 $y$ 方向为周期性（Periodic），沿 $x$ 方向为开放式（Open）。

1.3.2 集群平均场理论（Cluster Mean-Field Theory, CMFT）

虽然 DMRG 是处理一维和准二维条带体系的利器，但在高维或无限大二维极限下，由于边界和有限宽度的约束，直接观察自旋自发对称性破缺（如螺旋有序）较为困难。因此，作者引入了 CMFT 作为辅助手段。CMFT 的物理图像是将整个无限晶格划分成一个个相同大小的集群（Cluster），集群内部的自旋-自旋相互作用使用精确的 DMRG 进行完全求解，而集群之间的相互作用则用平均场（Mean Field）近似：

其中 $H_C$ 是集群 $C$ 内部的哈密顿量，$\mathbf{h}_i^{\text{eff}} = J \sum_{j \notin C} \langle \mathbf{S}_j \rangle$ 是由外部相邻集群中自旋的平均期望值产生的有效磁场。通过自洽地求解集群内部的自旋期望值 $\langle \mathbf{S}_i \rangle$，CMFT 能够极其灵敏且直观地捕捉到二维热力学极限下的长程磁有序形态，并能给出有序磁矩的绝对大小。

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

研究团队针对 $S=1$ 和 $S=3/2$ 体系，在多条关键切线上进行了精密的数值扫描，获得了令人信服的数据。

2.1 S = 1 体系的相变与中间物相

2.1.1 相界面的精确测定

作者首先计算了基态能量密度 $e_0$ 对耦合比 $g$ 的二阶导数 $- \partial^2 e_0 / \partial g^2$ 以及双分纠缠熵（Bipartite Entanglement Entropy）$S_e$：

• 在 $L_y = 6$ 圆柱体上，$- \partial^2 e_0 / \partial g^2$ 在 $g_{c1} \approx 0.47$ 和 $g_{c2} \approx 0.58$ 处表现出极为尖锐的奇点（如图 2(a) 所示），其发散行为取决于样品的采样密度，这正是典型的一阶量子相变特征。
• 纠缠熵 $S_e(g)$ 在这两个临界点同样出现了不连续的阶跃（Jump），进一步坐实了一阶相变。在这两个临界点之间，系统处于中间相。

2.1.2 混合斑块-二聚体（MPD）相

在 $g \in [g_{c1}, g_{c2}]$ 区域（以 $g = 0.5$ 为代表），自旋结构因子仅显示短程相关。然而，键能量（Bond Energy）分布展示了非平凡的行为：

• 在开放边界条件（OBC）下，如图 2(j) 所示，空斑块上的键能量表现出清晰的四聚体化（Tetramerization），即四个格点形成一个弱的单态斑块，同时伴随对角线上的极强二聚体相关。由于该相兼具斑块有序与二聚体特征，作者将其命名为混合斑块-二聚体（MPD）相。
• 在圆柱几何中，基态存在二重简并，两类空斑块（图 1(a) 中的橄榄色和蓝色斑块）发生对称性破缺。在有限尺寸系统中，边界能有效打破简并，显现出明晰的斑块图样。

2.1.3 非共度螺旋（Spiral）相的显现

当 $g > g_{c2}$ 时，系统并未直接进入 NAF 相，而是进入了非共度螺旋相。其特征由以下数据强力支撑：

• 结构因子峰值的连续移动：磁布拉格峰在波矢 $\mathbf{q} = (q_x, \pi)$ 处。随着 $g$ 增大，$q_x/\pi$ 从未锁定的非共度值连续且平滑地移动，直到在 $g_{c3} \approx 0.825$ 处锁定在 $(\pi, \pi)$，进入 NAF 相。如图 2(d) 所示，这一平滑变化曲线与经典极限的理论曲线（Olive 虚线）极其贴合。
• 实空间关联函数的代数衰减：在螺旋相区，自旋相关函数呈现带有空间振荡的幂律（Power-law）衰减，而非 NAF 相区的单调幂律衰减（图 2(c)）。
• CMFT 磁矩特征：在 $g = 1/\sqrt{2}$ 条件下，使用 CMFT 进行自洽求解，得到的实空间自旋织构（Spin Texture）呈现出完美的 $8 \times 2$ 周期性，磁矩方向沿 $x$ 轴平滑旋转，清晰展示了螺旋自旋波的表现（图 2(k)）。

2.2 S = 3/2 体系的相变与中间物相

对于自旋更高的 $S = 3/2$ 体系，作者们观测到了极其相似但演化趋势截然不同的物理现象：

• 相界点显著向左移动：第一和第二临界点分别位于 $g_{c1} \approx 0.36$ 和 $g_{c2} \approx 0.44$。这说明随着自旋 $S$ 的增大，二聚体相的稳定区间被大幅压缩。
• MPD 相的特征更为显著：如图 3(e) 和 3(f) 所示，在 $g=0.375$ 和 $g=0.4$ 时，即使在圆柱几何上，空斑块上的四聚体自旋关联和能量不均匀性也比 $S=1$ 更加明显，展现了更强的抗波动能力。其局部有序磁矩（CMFT 测量）也更接近经典饱和值（图 1(a)）。
• 螺旋相的温床进一步拓宽：螺旋相锁定的临界值在 $g_{c3} \approx 0.875$ 处，相较于 $S=1$ 的 $0.825$ 和 $S=1/2$ 的缺失，螺旋相的区域在 $S$ 增加时呈现连续扩张，完美地展示了向经典极限（$0 < g < 1$ 均为螺旋相）的收敛过渡。

2.3 全局 S-g 相图：量子至经典过渡的统一物理图像

本研究最具学术价值的发现莫过于图 1(b) 所示的 全局 $S-g$ 相图。通过整合 $S=1/2$（已知文献）、$S=1$、$S=3/2$（本工作）和 $S \to \infty$（经典极限），我们能够清晰地梳理出以下演化规律：

1. 极性相的演变：$S=1/2$ 时代的空斑块（EP）相在 $S \ge 1$ 时转化为混合斑块-二聚体（MPD）相。随着 $S$ 增加，这一斑块相的生存空间（$g$ 的宽度）不断萎缩，并整体向更小的 $g$ 移动，直至在经典极限下完全消失。
2. 螺旋相的诞生与扩张：在 $S=1/2$ 中完全不存在的螺旋（Spiral）相，自 $S=1$ 起强力介入。随着量子数 $S$ 增大，螺旋相的右边界向 $g=1$ 推进，左边界向 $g=0$ 靠拢。这雄辩地证实：螺旋相是连接高度量子化的自旋无序单态与经典磁有序的长程纽带。
3. 量子涨落的梯度抑制：大自旋在很大程度上保留了经典磁结构的骨架，但几何挫折带来的强竞争使得量子涨落即便在 $S=3/2$ 这样的大自旋体系中，依然能有效稳定非经典的 MPD 相。

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

为了使量子化学及凝聚态计算同仁能够复现论文中的关键结果，我们将详细拆解基于 ITensors 开源张量网络库的实现框架。

3.1 技术选型与开源仓库推荐

• 核心库：ITensors (推荐使用 Julia 语言版本，运行效率与编写体验最佳)。
• 对称性处理：必须开启 use_su2=true。Julia 下的 ITensors 官方分支和 NDTensors 提供了对自旋 $SU(2)$ 群的非阿贝尔（Non-Abelian）对称性支持。
• 计算资源建议：$L_y=6, L_x=24$ 的 $S=1$ 计算，推荐配置为 $32$ 核 Xeon 处理器，内存至少 $128 \text{ GB}$。$S=3/2$ 需适当调大内存需求。

3.2 核心算法骨架：ITensors 哈密顿量构建

以下是复现 SSL 晶格自旋 $S=1$ 及 $S=3/2$ 模型哈密顿量矩阵乘积算符（MPO）的 Julia 核心代码模板。本代码精确处理了具有周期性边界条件（y方向）的圆柱体格点映射：

using ITensors

# 参数定义
let
    Nx = 24  # 圆柱长度
    Ny = 4   # 圆柱宽度
    N = Nx * Ny
    S = "1"  # 若为 S=3/2，修改为 "3/2"
    g = 0.5  # 耦合比 J/J'
    Jp = 1.0 # 设定对角耦合为单位能量
    J = g * Jp

    # 初始化格点
    sites = siteinds("Spin", N; conserve_qns=true, conserve_su2=true)
    
    # 构建 MPO 描述器
    ampo = OpSum()

    # 辅助函数：将 2D 坐标 (x, y) 映射为 1D 链索引
    function coord_to_index(x, y)
        # y 方向使用周期性边界条件
        y_mod = mod1(y, Ny)
        return (x - 1) * Ny + y_mod
    end

    # 1. 铺设最近邻 J 键 (正方形边)
    for x in 1:Nx
        for y in 1:Ny
            idx1 = coord_to_index(x, y)
            # x 方向邻居 (注意开放边界)
            if x < Nx
                idx2 = coord_to_index(x + 1, y)
                ampo += J, "S+", idx1, "S-", idx2
                ampo += J, "S-", idx1, "S+", idx2
                ampo += J, "Sz", idx1, "Sz", idx2
            end
            # y 方向邻居 (圆柱体周期边界)
            idx3 = coord_to_index(x, y + 1)
            ampo += J, "S+", idx1, "S-", idx3
            # 注意：在使用 SU(2) 对称性时，相互作用的形式通常直接写为 "S" 算符的缩并，
            # ITensors 中可直接通过等效形式或内置 SU2 算符表达
        end
    end

    # 2. 铺设交替对角线 J' 键
    # 根据 Shastry-Sutherland 晶格定义，在特定的 2x2 单元内引入单条对角线
    for x in 1:(Nx - 1)
        for y in 1:Ny
            if (mod(x, 2) == 1 && mod(y, 2) == 1)
                # (2n+1, 2m+1) 单元内的斜对角键
                idx1 = coord_to_index(x, y)
                idx2 = coord_to_index(x + 1, y + 1)
                ampo += Jp, "S", idx1, "S", idx2
            elseif (mod(x, 2) == 0 && mod(y, 2) == 0)
                # (2n, 2m) 单元内的斜对角键
                idx1 = coord_to_index(x, y)
                idx2 = coord_to_index(x + 1, y - 1)
                ampo += Jp, "S", idx1, "S", idx2
            end
        end
    end

    H = MPO(ampo, sites)
    println("哈密顿量 MPO 构建成功！")
end

3.3 DMRG 收敛策略（DMRG Sweeps Setup）

由于挫折系统极易陷入亚稳态（Local Minima），在设置扫描参数时，务必引入噪声项（Noise Term），并在前几次 Sweep 中保持较大的噪声，以便体系跨越势垒，准确寻找到全局基态：

nsweeps = 20
maxdim = [100, 200, 500, 1000, 2000, 4000, 8000]
cutoff = [1e-8, 1e-10, 1e-12, 1e-12, 1e-12, 1e-12, 1e-12]
noise = [1e-4, 1e-5, 1e-6, 1e-8, 1e-10, 0.0, 0.0]

# 运行 DMRG
energy, psi = dmrg(H, psi_init; nsweeps, maxdim, cutoff, noise)

3.4 CMFT 自洽迭代核心步骤

复现集群平均场理论计算，其基本逻辑如下：

1. 初始化随机的实空间自旋期望值矢量 $\vec{M}_i = \langle \mathbf{S}_i \rangle$，长度为集群格点数。
2. 根据当前 $\vec{M}$ 值，计算边界上由相邻集群产生的有效外场 $\mathbf{h}_i^{\text{eff}}$。
3. 将有效外场作为单体项 $\mathbf{h}_i^{\text{eff}} \cdot \mathbf{S}_i$ 加入到集群哈密顿量 $H_C$ 中，构建新的集群哈密顿量。
4. 运行有限尺寸 DMRG，求解带有这些边界单体外场的 $H_C$ 的基态 $\lvert \psi_C \rangle$。
5. 计算该基态下集群内所有自旋的期望值新分布 $\vec{M}_i^{\text{new}} = \langle \psi_C \lvert \mathbf{S}_i \rvert \psi_C \rangle$。
6. 使用混合公式更新磁矩：$\vec{M}^{\text{next}} = \alpha \vec{M}^{\text{new}} + (1-\alpha) \vec{M}^{\text{old}}$，其中调和因子 $\alpha \in [0.1, 0.3]$ 防止数值振荡。
7. 重复步骤 2-6，直至 $\lvert \vec{M}^{\text{next}} - \vec{M}^{\text{old}} \rvert < 10^{-6}$。此时的 $\vec{M}$ 即为自洽解。

4. 关键引用文献与局限性深度点评

4.1 关键里程碑文献

本研究建立在数十年来 SSL 研究的深厚基石之上，其最关键的参考文献包括：

• [1] Shastry & Sutherland, Physica B+C (1981)：提出经典的 SSL 模型并证明了其精确二聚体基态解。
• [2] Koga & Kawakami, Phys. Rev. Lett. (2000)：首次确立了 $S=1/2$ 体系中的空斑块（EP）相，开启了 SSL 中间量子相的研究热潮。
• [20] Yuan, Wu, Yao & Wu, arXiv:2601.22924 (2026)：在 $S=1/2$ 极限下的 $XXZ$ 各向异性模型中首次发现了类似的螺旋相，为本研究在大自旋各向同性海森堡模型中寻找螺旋相提供了直接启示。
• [53] Fishman, White & Stoudenmire, SciPost Phys. Codebases (2022)：ITensor 库的奠基性文献，提供了处理大尺度自旋体系所必须的张量网络底层支持。

4.2 本文工作局限性剖析与学术批判

虽然该项工作物理图像清晰、数值计算技术极其精湛，但站在前沿科研角度审视，其研究仍存在以下局限性：

1. 二维热力学极限（TDL）的外推限制： 尽管作者在 $L_y = 4, 6$ 的圆柱体上进行了高精度的计算，但在二维多体物理中，$L_y = 6$ 依然属于准一维条带。尤其是在探索高度非共度、波矢连续变化的螺旋相时，$L_y$ 的尺寸截断会强制施加一维约束。这导致作者无法在 $L_y \to \infty$ 的真正二维极限下精确外推螺旋锁定转变点 $g_{c3}$ 的热力学极限值。未来有必要借助**二维无限张量网络（iPEPS）**进行无偏校验。

2. SU(2) 对称性对复杂自旋结构的双刃剑效应： 完全保留 $SU(2)$ 对称性虽然极大地压低了计算成本，但它强制体系的态处于总自旋 $J=0$ 的单态。这在研究具有自发磁自旋对称性破缺（如呈现特定自旋取向的螺旋有序）的物相时，会导致基态被迫以对称性恢复的叠加态形式出现（如文中提到的 MPD 相在圆柱体上的叠加掩盖）。虽然作者巧妙地通过开放边界条件和 CMFT 解决了这一问题，但在 DMRG 内部直接分析自发对称性破缺动力学（如自旋转置）仍然受到对称性保护的制约。

3. 对真实材料各向异性的简化： 论文主要讨论的是各向同性海森堡模型。然而，实际合成的稀土基大自旋 SSL 材料（例如含有稀土离子 $\text{Pr}^{3+}, \text{Nd}^{3+}, \text{Gd}^{3+}$ 等的 $\text{RE}_2\text{Be}_2\text{GeO}_7$）由于强轨道-自旋耦合，通常伴随有不可忽略的单离子各向异性（Single-ion Anisotropy）、Dzyaloshinskii-Moriya（DM）相互作用或 $XYZ$ 各向异性。因此，该模型得出的相界与实际材料的实验拟合之间可能存在显著的定量偏差。

5. 补充探索：从多体理论到实验探测的桥梁

为了使该研究的理论成果能切实指导实验室中的物理测量，我们在此补充大自旋 SSL 材料的实验可观测特征，帮助读者在自己的实验数据（如中子散射、热力学测量）中识别 MPD 相与螺旋相。

5.1 混合斑块-二聚体（MPD）相的实验特征

• 非弹性中子散射（INS）：由于 MPD 相中存在弱的四聚体化空斑块，INS 谱图将在低能区展现独特的定域激子激发（Localized Exciton Excitations）。与 $S=1/2$ 空斑块相类似，由于其强二聚体背景，其激发谱会在动量空间呈现明显的选择定域化。而在高能区，可以观测到清晰的单态-三态（Singlet-Triplet）激发，其能隙大小与 $g$ 呈现非单调依赖。
• 比热与磁化率：MPD 相具有显著的自旋能隙。在低温磁化强度曲线 $M(H)$ 中，由于能隙的存在，将表现出明确的零磁化平台（Zero-magnetization Plateau）。随着外场克服能隙，将出现分数磁化平台（如 $1/8$ 或 $1/4$ 平台），其平台位置直接对应空斑块的空间排布周期。

5.2 非共度螺旋（Spiral）相的实验特征

• 中子衍射谱图（Neutron Diffraction）：在螺旋相中，静态结构因子 $S(\mathbf{q})$ 会在非共度位置 $\mathbf{q} = (q_x, \pi)$ 产生极为尖锐的弹性能量磁共振峰。通过施加压力（改变量子相互作用比例 $g$）或调节磁场，可以实时观测到该共振峰在倒空间沿 $q_x$ 轨道的连续滑移。一旦系统进入 NAF 相，该峰将瞬间锁定在共度的 $(\pi, \pi)$（即 $M$ 点）。
• 核磁共振（NMR）谱线分裂：由于非共度自旋调制的引入，格点上的局域内磁场不再是单一的常数，而是随着螺旋周期在实空间呈正弦或余弦变化。这将直接导致核磁共振谱线发生特征性的非共度粉末状展宽或多重分裂（Incommensurate NMR Line Shape），是判定螺旋有序最直接、最灵敏的微观局域探针。

5.3 展望：各向异性与外加磁场的联合效应

在未来的科研图景中，两大方向最值得期待：

1. 各向异性相图构建：在哈密顿量中引入 $D(S^z)^2$ 项。对于 $S=1$，正的单轴各向异性将倾向于使系统处于大自旋单态（Large-S Singlet），它与 MPD、Dimer 相的竞争将催生更加复杂的拓扑物相。
2. 外加磁场下的拓扑磁结构：大自旋系统在几何挫折与外加磁场共同作用下，是培育**自旋超固体（Spin Supersolid）和手性斯格明子晶格（Skyrmion Lattice）**的天然温床。本文所建立的 $SU(2)$ 计算基石，为这些前沿课题铺平了道路。