来源论文: https://arxiv.org/abs/2603.00536v1 生成时间: Mar 09, 2026 09:46

0. 执行摘要

本研究论文深入探讨了在3x3几何阻挫方格（具有周期性边界条件）上的Hubbard模型中，电子关联如何驱动交错磁性（Altermagnetism）的涌现。交错磁性是一种新颖的共线磁性相，其特征是具有显著的非相对论性自旋劈裂，但净磁化为零，并保持了独特的组合对称性。与传统依赖单粒子各向异性或特定亚晶格几何结构的模型不同，这项工作通过精确对角化（Exact Diagonalization, ED）方法，系统地分析了在位库仑排斥U和近邻（NN）库仑相互作用V与不同填充因子（N=8、9、10）下这些参数的相互作用及其对磁性的影响。研究发现，在半填充（N=9）时，几何阻挫会抑制交错磁性，但在载流子掺杂（空穴掺杂N=8或电子掺杂N=10）的情况下，却能稳定形成鲁棒的、旋转对称的交错磁性关联。尤其是在电子掺杂体系中，当U处于中等耦合强度（U=4）时，NN相互作用V存在一个临界阈值，在该阈值处交错磁性态会发生类似一级相变的急剧变化。而在强耦合（U=10）区域，局部磁矩形成良好，并能保持其固有的自旋各向异性纹理，同时相变点被推迟到更高的V值。这项研究为在对称几何阻挫晶格上实现交错磁性提供了一条波动介导的途径，并确定了载流子浓度和NN库仑相互作用作为调控强关联体系中磁各向异性的关键可调参数，为未来新材料的设计和实验研究提供了宝贵的基础。本文的核心贡献在于，通过严谨的精确对角化计算，证实了即使在没有自旋轨道耦合或单粒子各向异性等外部因素的情况下，纯粹的电子关联作用也能在阻挫体系中诱导出稳定的交错磁性，这极大地拓展了我们对这种新型磁性相起源的理解，并为强关联交错磁性研究开辟了新途径。

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

核心科学问题

本研究的核心科学问题聚焦于在没有单粒子各向异性或自旋轨道耦合（Spin-Orbit Coupling, SOC）这一通常被认为是交错磁性先决条件的机制时，纯粹的电子关联作用如何在几何阻挫的晶格上驱动交错磁性的涌现？更进一步，这种关联驱动的交错磁性是否能够保持其鲁棒性和重要的旋转对称性？传统的交错磁性理论往往依赖于特定的晶体结构对称性，或者通过引入自旋轨道耦合来产生动量空间中的自旋劈裂。然而，在诸多强关联材料中，几何阻挫是普遍存在的，它通常会抑制常规磁序的形成，为更复杂的量子相提供温床。因此，理解在这些高度竞争且复杂的相互作用背景下，交错磁性能否作为一种内禀的电子现象自发出现，对于拓展现有磁性分类，并为新型自旋电子学材料的设计提供理论指导具有关键意义。

理论基础

  1. 交错磁性 (Altermagnetism)：

    • 定义与性质：交错磁性是一种独特的新型共线磁性相，它在宏观上表现为零净磁化强度，与反铁磁类似。然而，其核心特征在于，它在动量空间中存在巨大的、非相对论性的自旋劈裂，类似于铁磁。与传统反铁磁不同，交错磁性的自旋劈裂并非由时间反演对称性与平移对称性或反演对称性的组合所保护。相反，它通过其独特的晶体对称性与自旋对称性的结合，产生了一种动量依赖的自旋劈裂，这种劈裂模式可以类比为d波、g波或i波配对对称性。这种独特的性质组合使得交错磁性在下一代自旋电子学器件、太赫兹振荡器以及拓扑量子计算等前沿领域具有巨大的潜在应用价值。
    • 微观起源：交错磁性的微观起源与其内部局域磁矩的形成以及这些磁矩之间自旋关联的空间组织模式紧密相关。本研究旨在探索，在缺乏外部对称性破缺源（如SOC）的情况下，纯粹的电子-电子相互作用能否自发地产生这种特殊的、具有交错磁性特征的自旋关联模式。

  2. Hubbard 模型及其扩展 (Hubbard Model and its Extensions)：

    • Hubbard 模型：作为描述强关联电子系统最简洁的微观模型，Hubbard模型成功捕捉了电子动能（t）与在位库仑排斥（U）之间的基本竞争。动能倾向于电子的离域化，而U则促进电子在格点上的局域化，进而形成局域磁矩。通过调节U/t比值，该模型可以揭示系统从巡游金属态到Mott绝缘体态的转变，并能模拟多种量子相，如非常规超导体和各种奇异磁性态的涌现。在强耦合极限下，Hubbard模型可以约化为有效的Heisenberg自旋哈密顿量，而在中等耦合区域，它保留了电荷和自旋涨落的相互作用，使其对载流子掺杂和几何阻挫特别敏感。
      • 哈密顿量 (Simple Hubbard Model)：H = -t Σ_⟨i,j⟩,σ (c†_i,σ c_j,σ + h.c.) + U Σ_i n_i↑ n_i↓
        • 第一项代表电子在近邻格点间的跳跃，强度为t。c†_i,σ 和 c_j,σ 分别是格点i和j处自旋为σ的电子的产生和湮灭算符。
        • 第二项描述在同一格点上两个自旋相反的电子之间的在位库仑排斥，强度为U。n_i↑ 和 n_i↓ 分别是格点i处自旋向上和向下的电子数算符。
    • 扩展 Hubbard 模型 (Extended Hubbard Model)：为进一步探讨非局域关联的影响，该模型在简单Hubbard模型的基础上引入了近邻（NN）库仑相互作用V。这种扩展使得模型能够捕捉更丰富的物理现象，例如电荷序、相分离以及它们与磁性序的竞争。
      • 哈密顿量 (Extended Hubbard Model)：H = H_Hubbard + V Σ_⟨i,j⟩ n_i n_j
        • 额外项 V Σ_⟨i,j⟩ n_i n_j 描述了近邻格点i和j上的电子总数n_i和n_j之间的库仑相互作用，强度为V。n_i = n_i↑ + n_i↓。

  3. 几何阻挫 (Geometric Frustration)：

    • 核心作用：几何阻挫是指晶格的几何结构导致竞争性的相互作用无法同时满足所有最低能量的构型，从而阻止了简单长程磁序（如Néel反铁磁序）的形成。在这种情况下，系统倾向于形成高度简并的基态流形、自旋液体态或复杂的空间调制自旋纹理等非常规关联态。阻挫体系的能量景观通常非常平坦，存在许多简并或准简并的基态，使得系统对外部扰动或微小参数变化异常敏感。
    • 3x3 方格作为研究平台：本研究特别选择了具有周期性边界条件的3x3方格作为几何阻挫晶格（如图1所示）。这个最小的环形几何结构不仅保留了完整的量子相干性，更重要的是，它内在包含了几何阻挫和非平凡的环路连接。作为精确对角化可处理的最大尺寸之一，该体系足够小以进行精确计算，同时又足够复杂以展示阻挫效应。它通过抑制经典的Néel序倾向，增强了对对称性选择性自旋响应的敏感性，因此成为识别交错磁性关联是内禀电子现象还是自发对称性破缺产物的理想受控平台。

技术难点

  1. 希尔伯特空间维度爆炸 (Exponential Scaling of Hilbert Space)：
    • 精确对角化方法的核心挑战在于其计算成本随系统尺寸呈指数增长。对于包含 N_sites 个格点，每个格点有自旋向上和自旋向下两种占据态的Hubbard模型，理论上的希尔伯特空间维度是 4^N_sites。即使利用总粒子数N和总S_z（自旋的z分量）守恒等对称性进行分块，对于9个格点的系统，维度依然相当庞大。例如，在半填充（N=9）且 S_z=1/2 的情况下，希尔伯特空间维度约为15876。如何高效构建如此大规模的哈密顿量矩阵，并准确求解其最低本征态（基态），是计算上的核心难点。这需要精密的算法设计，如稀疏矩阵存储和迭代求解器。
  2. 阻挫体系的基态复杂性：
    • 几何阻挫常常导致系统基态的高度简并或准简并。在这种情况下，单个基态的平均物理量可能不足以完全描述系统的真实性质，因为系统可能在多个简并态之间存在量子涨落。对这些复杂基态流形的分析需要仔细区分真正的相变与由于简并性导致的平滑演化，并可能需要评估基态的稳定性。
  3. 交错磁性关联的精确量化：
    • 交错磁性并非简单的长程磁序，它通过特定的空间和动量特征来定义。因此，如何设计和计算合适的观测量，以准确地捕捉其存在的微观证据，并能够清晰地区分其与传统反铁磁等其他磁性序的不同之处，是一个重要的挑战。本研究通过引入“平均交错磁性自旋响应量度” (Aspin)、“平均局域磁矩” (m) 和“各向异性参数” (A) 等多个精心定义的观测量来解决这一问题，但这些观测量本身在数值计算和物理诠释上需要严谨对待。
  4. 多参数空间的探索：
    • 本研究需要在由在位库仑排斥U、近邻库仑相互作用V、总电子数N（或填充因子）以及总自旋S（或其z分量S_z）构成的多维参数空间中进行系统探索。在如此广阔的参数空间中，高效地识别关键的物理区域（如相变点、稳定相的边界）并精确计算相应的物理量，需要庞大的计算资源和细致的数据分析策略。

方法细节

  1. 精确对角化 (Exact Diagonalization, ED)：

    • ED是一种强大的数值方法，它通过构建并对角化有限尺寸量子系统的哈密顿量全矩阵，以获得其精确的本征态和本征值。作为一种无近似方法，ED被认为是研究强关联系统基态性质和低激发态的黄金标准。其主要限制是计算成本随系统尺寸呈指数增长。
    • 实现策略：
      • Lanczos 算法：为高效地提取大型稀疏哈密顿量矩阵的最低（或最高）几个本征值及其对应的本征向量（即基态），通常采用Lanczos算法。该算法将原矩阵投影到一个小得多的三对角矩阵上，从而显著降低了计算复杂性。
      • 对称性利用：为有效降低希尔伯特空间的维度，本研究利用了守恒量进行分块对角化。具体包括：
        • 总电子数N：通过固定系统中的电子总数来缩小计算空间。
        • 总自旋的z分量S_z：通过固定自旋向上和自旋向下电子的数量差来分块哈密顿量。论文中针对半填充（N=9）设定总自旋S=1/2（这通常在S_z=1/2扇区计算），以及对掺杂系统（N=8和N=10）设定总自旋S=0（这通常在S_z=0扇区计算）。这种设定意味着基态是特定总自旋S的本征态，通常通过在S_z守恒的扇区中找到能量最低的本征态并验证其S值来确定。
        • 总动量K：对于周期性边界条件下的均匀晶格，总动量K也是一个守恒量，可以用于进一步分块哈密顿量矩阵，但论文中未明确提及利用K对称性。

  2. 晶格与填充因子：

    • 3x3 方格：研究对象是一个包含9个格点的方格晶格，通过施加周期性边界条件（PBC）形成一个环面几何结构（如图1所示）。这种结构确保了所有格点等价，并引入了固有几何阻挫。
    • 填充因子：研究了三种具有代表性的电子填充情况：
      • N=9：半填充情况，即每个格点平均一个电子。这通常是Mott绝缘体相出现的典型区域。基态总自旋设定为S=1/2。
      • N=8：空穴掺杂情况，比半填充少一个电子。基态总自旋设定为S=0。
      • N=10：电子掺杂情况，比半填充多一个电子。基态总自旋设定为S=0。

  3. 观测量定义：

    • 平均交错磁性自旋响应量度 (Aspin)：
      • 定义为 Aspin = (1/Nb) Σ_⟨i,j⟩ Ω_ij，其中 Nb 是近邻键的总数。Ω_ij = ⟨S_i ⋅ S_j⟩ η_ij，⟨S_i ⋅ S_j⟩ 是格点i和j之间的自旋关联函数，η_ij 是一个精心设计的交错符号因子。这个符号因子η_ij被构造用来捕捉交错磁性特有的C(2,2,1,2’)对称性，其具体形式为(-1)^(r_ix+r_iy+r_jx+r_jy)，其中(r_ix, r_iy) 和 (r_jx, r_jy) 分别是格点i和j的晶格坐标。这个量度通过平均近邻自旋关联的加权和来量化交错磁性关联的强度和空间组织模式，其非零值直接指示交错磁性的存在。
    • 平均局域磁矩 (m)：
      • 定义为 m = (1/N_sites) Σ_i ⟨(n_i↑ - n_i↓)^2⟩。N_sites是格点总数（这里是9）。这个量度直接量化了电子在格点上的局域化程度和局域磁极化的强度。当m值接近1时，表明每个格点上形成了稳定的局域磁矩，系统处于Mott绝缘体行为；当m值接近0时，表明电子处于巡游态。
    • 各向异性参数 (A)：
      • 定义为 A = (1/M) Σ_i |Δ^x_i - Δ^y_i|，其中M是格点总数。Δ^x_i 和 Δ^y_i 分别是针对格点i，沿x方向和y方向计算的平均交错磁性自旋响应。更具体地说，Δ^x_i 可以理解为格点i与所有沿x方向（包括直接近邻和通过PBC连接的更远格点）的键的Ω_ij的平均贡献，类似地定义Δ^y_i。这个参数量化了交错磁性关联在x和y晶格方向上的方向依赖性（即旋转对称性破缺程度）。A的有限值表示存在方向性各向异性，而A接近零则表明体系保持了旋转对称性（C4对称性）。

  4. 参数探索范围：

    • 在位库仑排斥U：从0t到12t的宽泛范围，覆盖了从弱关联金属态到强耦合Mott绝缘体态。
    • 近邻库仑相互作用V：在0到U/2的范围内进行探索，主要关注其在中等耦合（U=4）和强耦合（U=10）情况下的影响。
    • 通过系统地改变这些参数并计算上述观测量，研究人员能够全面揭示交错磁性相在不同物理条件下的涌现、稳定性、相变行为以及对对称性的影响。

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

关键 Benchmark 体系

本研究的核心 benchmark 体系是 3x3 几何阻挫方格上的Hubbard模型。这个体系的选择是经过深思熟虑的，具有以下重要原因：

  • 尺寸与精确性：9个格点（3x3）的系统是当前精确对角化方法能够在不引入近似的情况下进行计算的最大尺寸之一。这确保了计算结果的精确性，避免了平均场或其他近似方法可能引入的偏差，从而为理解强关联系统提供了无偏见的基础数据。
  • 几何阻挫特性：3x3方格结合周期性边界条件（PBC）形成一个环形结构，自然引入了几何阻挫。这种阻挫效应对于抑制常规磁序（如Néel反铁磁）并促进非常规量子相（如交错磁性）的形成至关重要。作为最小的阻挫体系，它能够捕捉到阻挫对磁性关联的深刻影响。
  • 对称性保留：该体系具有中心反演对称性。在这样的对称晶格中，交错磁性的出现不能归因于单粒子各向异性或晶格对称性破缺。因此，它提供了一个理想的平台来验证交错磁性是否可以纯粹由电子关联驱动，是内禀的电子现象。
  • 载流子敏感性：小尺寸体系对载流子浓度（掺杂）的微小变化非常敏感，这使得研究掺杂效应如何调控磁性相变得可行。

计算所得数据及分析

论文通过对多个关键观测量（Aspin、m、A）在不同参数U、V和填充因子N下的行为进行系统分析，揭示了交错磁性的涌现机制和稳定性。以下是主要结果的详细分析：

  1. 简单 Hubbard 模型 (V=0) 中的结果：

    • Aspin 随 U 的变化 (图2)：
      • 半填充 (N=9, S=1/2)：Aspin 随 U 的增加而单调衰减，并在强耦合极限（U > t）下趋近于零。这表明，在半填充状态下，尽管形成了局域磁矩，但几何阻挫效应阻止了这些磁矩组织成具有交错磁性特征的特定各向异性模式。相反，系统可能进入一种自旋液体状或高度无序的自旋态，导致交错磁性关联被有效抑制。
      • 掺杂系统 (N=8, S=0 和 N=10, S=0)：与半填充情况形成鲜明对比，空穴掺杂 (N=8) 和电子掺杂 (N=10) 系统都表现出 Aspin 随 U 增加而显著增强的趋势。这种增强对于电子掺杂 (N=10) 尤其明显，Aspin 迅速上升并在大U极限下饱和。这强烈表明，载流子掺杂是诱导和稳定交错磁性关联的关键机制。移动载流子通过优先选择特定的各向异性跳跃路径来降低动能，并通过在位U相互作用将这种动能耦合到自旋自由度，从而稳定了交错磁性特有的交错、补偿自旋模式。
    • m 随 U 的变化 (图3)：
      • 所有填充情况 (N=8, 9, 10)：m 随 U 的增加而单调递增。这反映了系统从弱关联的巡游金属态向局域磁矩明确的Mott绝缘态的普遍交叉行为。U越大，电子占据双位的能量代价越高，电子越倾向于局域化形成局域磁矩。
      • 对比 Aspin 和 m：半填充 (N=9) 尽管表现出最强的局域磁矩 (m 增长最快且达到最高值)，但 Aspin 却最小。这再次强调了在阻挫体系中，局域磁矩的形成本身不足以产生交错磁性，还需要特定的空间组织。而掺杂系统 (N=8, 10) 中，局域磁矩的增长与 Aspin 的显著增强同步发生，表明移动载流子在促进局域磁矩形成的同时，也促进了它们的有序组织。
    • 实空间热图 Ω_ij (U=4, V=0.2) (图4)：
      • 半填充 (N=9)：实空间自旋关联 Ω_ij 普遍较弱且高度碎片化，没有清晰的d波或交错对称性。这支持了半填充下关联混乱且无序的结论。
      • 掺杂系统 (N=8, N=10)：显示出清晰、结构化的各向异性图案，表明存在自旋关联的特定空间组织。电子掺杂 (N=10) 的图案尤其强烈和清晰，并呈现出显著的粒子-空穴不对称性，这可能与晶格能带色散和几何阻挫的复杂相互作用有关。
    • 各向异性参数 A 随 U 的变化 (图5)：
      • 半填充 (N=9)：在小U时 A 值有限（约 10^-3 量级），但随U增加而逐渐减小，并在强耦合极限下趋近于零。这表明几何阻挫在半填充下抑制了相干的对称相形成，系统倾向于各向同性关联。
      • 掺杂系统 (N=8, N=10)：A 值在整个U范围内几乎为零。这证实了掺杂系统中的交错磁性关联是旋转对称的，即使在强关联区域也保持了C4晶格对称性，表明形成了一个鲁棒的、方向性平衡的交错磁性基态。

  2. 扩展 Hubbard 模型 (V>0) 中的结果 (U=4)：

    • Aspin 随 V 的变化 (图6)：
      • 半填充 (N=9)：Aspin 始终保持在低水平，表明NN相互作用V并未能克服几何阻挫，反而可能进一步抑制了本已微弱的交错磁性关联。
      • 空穴掺杂 (N=8)：Aspin 随V的增加而平稳下降，表明交错磁性态在此中间V区域内保持稳定，但强度有所减弱。
      • 电子掺杂 (N=10)：Aspin 初始平稳下降，但在 V ≈ 1.8 附近出现急剧的、类似一级相变的骤降。这强烈表明，当NN相互作用V达到临界值时，电子态发生快速重构，可能向电荷序或非磁性相转变，竞争性构型变得更具能量优势。
    • m 随 V 的变化 (U=4) (图7)：
      • 所有填充情况 (N=8, 9, 10)：m 随 V 的增加而普遍下降，这反映了NN库仑排斥倾向于抑制局域磁矩的形成或增强电荷涨落。对于电子掺杂 (N=10) 系统，在 V ≈ 1.8 处同样出现急剧下降，与 Aspin 的变化同步，再次确认了相变的存在。局域磁矩在此处大幅减少但未完全消失，表明系统并非转向完全非磁性态，而是向一种保留有限局域磁矩的竞争性构型转变。
    • 实空间热图 Ω_ij (U=4, V=0.2, 1, 2) (图8)：
      • 半填充 (N=9)：关联依然弱且碎片化，V的增加未能诱导任何相干有序，证实了阻挫的主导作用。
      • 空穴掺杂 (N=8)：特征性的各向异性图案依然清晰可见，但强度随V增加而逐渐减弱，结构保持完整。
      • 电子掺杂 (N=10)：在 V=0.2 和 V=1 时，显示清晰的交错磁性图案；在 V=2 时（相变后），图案变得模糊且强度显著降低，证实了相变导致空间关联的重组。
    • 各向异性参数 A 随 V 的变化 (U=4) (图9)：
      • 半填充 (N=9)：A 值小而有限，随V增加而下降。未发生明显的对称性破缺。
      • 空穴掺杂 (N=8)：A 始终接近零，表明 C4 对称性在 NN 库仑排斥存在下依然保持，交错磁性态具有鲁棒的旋转对称性。
      • 电子掺杂 (N=10)：A 在 V ≈ 1.8 附近出现尖锐的峰值，与 Aspin 和 m 的骤降同步。这明确指示了在此相变点处，旋转对称性被瞬时打破或发生剧烈重构。

  3. 强耦合下的结果 (U=10)：

    • Aspin 随 V 的变化 (图10)：
      • 半填充 (N=9)：Aspin 仍然较低，但相比U=4略有增加。强耦合并未完全消除阻挫效应。
      • 掺杂系统 (N=8, N=10)：Aspin 显著高于U=4时的情况，并随V增加而平稳下降。电子掺杂 (N=10) 的相变点被推迟到 V ≈ 4.9，表明强耦合U增强了交错磁性相的稳定性，使其对NN相互作用的鲁棒性更强。
    • m 随 V 的变化 (图11)：
      • 所有填充情况 (N=8, 9, 10)：m 随 V 增加而缓慢下降，但始终保持较高的局域磁矩。这表明强在位排斥确保了即使有NN排斥，电子也倾向于保持局域化。电子掺杂 (N=10) 的相变点同样发生在 V ≈ 4.9，m值在此处急剧下降，但仍保持有限值，再次证明相变不是完全非磁性转变。
    • 实空间热图 Ω_ij (U=10, V=1, 3, 5) (图12)：
      • 掺杂系统 (N=8, N=10)：即使在较高的V值下，特征性的交错磁性图案依然保持结构完整，尽管强度有所降低。这证明了强在位排斥U对交错磁性纹理的保护作用，使其结构免受NN相互作用的剧烈破坏。
    • 各向异性参数 A 随 V 的变化 (图13)：
      • 半填充 (N=9)：A 依然很小且有限，但其涨落反映了阻挫Mott态的复杂性。
      • 掺杂系统 (N=8, N=10)：A 几乎可以忽略不计。即使在电子掺杂 (N=10) 的相变点附近，也只有很小的增强，没有U=4时的尖锐峰值。这表明，在强耦合条件下，强在位排斥U有效地抑制了对称性破缺不稳定性，确保了交错磁性相在更广阔的参数范围内保持旋转对称性。

性能数据

论文中并未直接提供具体的计算性能数据，如CPU/GPU时间、内存消耗或并行效率。然而，基于精确对角化（ED）方法的特点和计算物理领域的实践，我们可以对其性能进行合理推断：

  • 计算复杂度：ED的计算成本主要由其希尔伯特空间的维度决定，并呈指数增长。对于一个包含 N_sites 个格点的Hubbard模型，在利用粒子数N和S_z对称性后，希尔伯特空间的维度D仍然可以很大。对于本研究的9个格点系统：
    • N=9 (半填充，S_z=1/2)：希尔伯特空间维度为 C(9,5) * C(9,4) = 126 * 126 = 15876。
    • N=8 或 N=10 (掺杂，S_z=0)：以 N=8 为例，S_z=0对应 N_up = 4, N_down = 4，维度为 C(9,4) * C(9,4) = 126 * 126 = 15876。这些维度对于现代计算资源而言是完全可管理的。
  • 求解算法：Lanczos算法是求解大型稀疏矩阵基态的常用方法，其每次迭代主要涉及稀疏矩阵与向量的乘法运算，计算成本与非零矩阵元素的数量（通常与 D * N_sites 成正比）和迭代次数有关。对于维度在 10^4 级别的矩阵，Lanczos算法通常能在几十到几百次迭代内收敛到基态，在标准工作站或小型计算集群上，单次计算通常只需数分钟到数小时。
  • 内存需求：存储稀疏哈密顿量矩阵和基态波函数所需的内存，对于 D ~ 1.6 x 10^4 的维度，通常在几十GB到几百GB的范围内，这在现代高性能计算节点上是标准配置。
  • 并行化：矩阵的构建和Lanczos迭代中的矩阵向量乘法都可以有效地并行化，利用多核CPU来加速计算。尽管论文没有明确说明，但高效的ED实现通常会利用OpenMP或MPI进行并行计算。

总的来说，虽然论文没有提供具体的性能数字，但从其方法和系统规模来看，本研究所涉及的计算在计算物理界是可行的，并且通过精确对角化获得的无近似结果，为理解强关联交错磁性提供了宝贵的基准。

代码实现细节

论文中并未提供具体的代码实现细节，也没有公开任何代码仓库。然而，从其方法描述和引用的作者团队过往文献（如[53, 54, 55, 56]）来看，可以推断出以下可能的实现策略和技术：

  1. 编程语言与库：精确对角化通常使用高性能的编译型语言（如C++或Fortran）实现，以最大限度地提高计算效率。常用的科学计算库可能包括：

    • 线性代数库：LAPACK、BLAS（用于底层的稠密矩阵运算，但更常用的是作为稀疏矩阵库的底层加速器），或者更专业的稀疏矩阵库。例如，如果使用Python作为接口，SciPy.sparse模块是常见的选择，其底层是C/Fortran。
    • 并行计算：考虑到计算量，OpenMP（共享内存并行）或MPI（分布式内存并行）可能被用于加速矩阵构建和Lanczos迭代。对于9个格点的系统，通常OpenMP在单个多核CPU节点上已足够高效。

  2. 基态构建与对称性利用：

    • Fock 态表示：最直接的方法是使用Fock态（即每个格点占据态的二进制表示）作为计算基。每个格点i可以有四种状态：空、自旋向上占据、自旋向下占据或双占据。使用二进制位操作可以高效地表示和操纵这些态。
    • 守恒量分块：为降低希尔伯特空间的维度，作者利用了守恒量进行分块对角化：
      • 总粒子数 (N)：只考虑具有特定总电子数N的Fock态。这直接缩小了基空间。
      • 总自旋的z分量 (S_z)：通过固定自旋向上和自旋向下电子的数量（N_up, N_down）来进一步分块希尔伯特空间。例如，对于N=9，S=1/2，通常会在N_up=5, N_down=4或N_up=4, N_down=5的扇区进行计算。
      • 总自旋 (S)：论文明确指出了一些填充情况下的总自旋S（例如N=9时S=1/2，N=8和N=10时S=0）。这暗示可能使用了更复杂的自旋适配基（如Slater行列式组成的线性组合，或通过Clebsch-Gordan系数构建的自旋本征态）。这种方法比单纯的S_z分块更为复杂，但能直接得到特定总S值的基态，并进一步降低维度。作者团队在参考文献[55, 56]中曾提及“配置空间截断方案”和“Rumer基方法”，这些可能与高效构建或处理自旋适配基的原理相关，即便这些工作主要是针对更大系统的截断近似。
      • 总动量 (K)：对于周期性边界条件下的均匀晶格，总动量K也是一个守恒量，可用于进一步分块哈密顿量矩阵，但论文中未明确提及是否利用此对称性。

  3. 哈密顿量矩阵的构建：

    • 稀疏矩阵：Hubbard哈密顿量，尤其是对于只有近邻跳跃和在位/近邻相互作用的模型，通常是高度稀疏的。因此，采用稀疏矩阵存储格式（如Compressed Sparse Row, CSR 或 Coordinate List, COO）是必需的，以节省内存并提高计算效率。
    • 矩阵元计算：实现一个函数，根据哈密顿量中的每一项（跳跃项、在位库仑U项、近邻库仑V项），计算其在Fock态基下的矩阵元 ⟨α|H|β⟩。这通常涉及遍历希尔伯特空间中的所有基态 |β⟩，并对每个|β⟩应用哈密顿量算符，将其映射到其他基态 |α⟩。例如，跳跃项会改变两个近邻格点的占据情况和自旋；U项检查双占据情况；V项检查近邻占据情况。

  4. 基态求解：

    • Lanczos 算法：这是求解大型稀疏矩阵极端本征值（如基态能量）和对应本征向量（基态）的标准迭代算法。它通过将矩阵投影到一个Krylov子空间，生成一个三对角矩阵，然后对这个小得多的三对角矩阵进行对角化。算法的关键在于高效实现稀疏矩阵-向量乘法。

复现指南

由于论文没有提供具体的代码，以下是一个通用的、基于标准ED技术复现该研究的指南：

  1. 环境设置：

    • 编程语言：推荐使用C++（搭配Eigen、Intel MKL或OpenBLAS/LAPACK等线性代数库）或Python（搭配NumPy和SciPy）。Python由于其简洁性和丰富的库生态系统，在原型开发和数据分析中非常受欢迎，但对于大规模矩阵操作，C++通常提供更高的性能。
    • 计算资源：建议在一个配备多核CPU的计算集群上运行，以利用并行计算加速矩阵构建和Lanczos迭代。

  2. 定义晶格和哈密顿量：

    • 晶格拓扑：首先，需要实现3x3方格的拓扑结构，包括周期性边界条件。可以构建一个邻接列表或邻接矩阵来表示格点之间的连接关系，并实现一个函数来获取每个格点的近邻列表。
    • 哈密顿量算符：根据简单Hubbard模型和扩展Hubbard模型的表达式，定义跳跃项、在位相互作用项和近邻相互作用项的矩阵元计算规则。这需要理解电子的产生和湮灭算符如何作用于Fock态，以及如何处理费米子的反对易关系（如引入Jordan-Wigner变换或直接处理Slater行列式）。

  3. 构建希尔伯特空间：

    • 态基生成：编写代码生成满足特定总粒子数N和总S_z（或S）的Fock态。这些态可以通过二进制位串或整数表示。例如，对于9个格点，可以表示为18位的二进制数（每格点2位，分别代表自旋向上和向下占据）。
    • 态到索引的映射：实现一个高效的双向映射（如哈希表或排序列表），将每个有效的Fock态映射到唯一的整数索引，反之亦然。这对于在矩阵中定位元素和访问波函数系数至关重要。

  4. 构建哈密顿量矩阵：

    • 选择一种稀疏矩阵格式（如CSR或COO）。
    • 遍历希尔伯特空间中的所有基态 |α⟩。对于每个 |α⟩，应用哈密顿量中的每一个项（跳跃、U、V），计算其作用结果 H|α⟩ = Σ_β C_βα |β⟩。其中 C_βα 就是矩阵元 H_βα。
    • 将计算出的非零矩阵元填充到稀疏矩阵中。此步骤通常是最耗时的部分，需要高度优化。

  5. 求解基态：

    • 利用Lanczos算法（可以自己实现，或使用现有库函数，如SciPy.sparse.linalg.eigs (Python), ARPACK for C/Fortran）计算哈密顿量矩阵的最低本征值（基态能量）和对应的本征向量（基态波函数）。
    • 确保Lanczos迭代收敛到所需的能量和波函数精度。

  6. 计算观测量：

    • 自旋关联：基于基态波函数 |ψ_0⟩，计算近邻自旋关联 ⟨S_i ⋅ S_j⟩。自旋算符 S_i 可以通过电子数算符 n_iσ 和电子产生/湮灭算符 c†_iσ, c_iσ 表示。例如，S_i ⋅ S_j = S_ix S_jx + S_iy S_jy + S_iz S_jz，其中S_iz = (n_i↑ - n_i↓)/2，S_i+ = c†_i↑ c_i↓，S_i- = c†_i↓ c_i↑。
    • η_ij 因子：根据格点坐标和交错模式计算 η_ij 符号因子。
    • Aspin, m, A：根据论文中给出的公式，利用计算出的 ⟨S_i ⋅ S_j⟩ 期望值和基态波函数来计算这些观测量。

  7. 数据分析与可视化：

    • 将计算结果（能量、Aspin、m、A值）存储为结构化数据（如CSV、JSON或HDF5文件）。
    • 使用Matplotlib (Python)、Gnuplot或其他绘图工具生成与论文中相似的图表，以便进行比较和验证。

所用的软件包及开源 Repo 链接

论文中 未明确提及 所使用的具体软件包和开源代码仓库链接。在计算物理领域，精确对角化研究常常依赖于研究团队自己开发和维护的定制化代码，以实现对模型、对称性利用和算法优化的高度控制。考虑到该论文作者团队在精确对角化方面有丰富的发表历史（参考文献[53, 54, 55, 56]），很可能使用的是他们内部长期迭代和优化的ED框架。

因此，目前没有可供直接访问的开源代码仓库链接。

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

关键引用文献

本研究工作建立在强关联物理和磁性材料的深厚理论基础之上，并积极吸收了交错磁性领域的最新发展。以下是其引用文献中对本研究具有奠基性或直接相关的关键部分：

  1. 交错磁性理论与综述：

    • [27] L. Šmejkal, J. Sinova, and T. Jungwirth, “Emerging research landscape of altermagnetism,” Phys. Rev. X 12, 040501 (2022)。这篇具有里程碑意义的综述为交错磁性领域描绘了广阔的研究图景，强调了其作为一种新兴磁性相的独特地位和巨大潜力。
    • [28] L. Šmejkal, J. Sinova, and T. Jungwirth, “Beyond conventional ferromagnetism and antiferromagnetism: A phase with nonrelativistic spin and crystal rotation symmetry,” Phys. Rev. X 12, 031042 (2022)。此文从严格的对称性角度详细定义了交错磁性，并将其与传统磁性序（铁磁、反铁磁）进行明确区分。
    • [29] I. Mazin (The PRX Editors), “Editorial: Altermagnetism—a new punch line of fundamental magnetism,” Phys. Rev. X 12, 040002 (2022)。《Physical Review X》编辑部这篇社论概括了交错磁性作为“基础磁学的新亮点”的地位，突出了其在科学界引起的热烈关注。
    • [30] T. Jungwirth, R. M. Fernandes, E. Fradkin, A. H. MacDonald, J. Sinova, and L. Šmejkal, “Altermagnetism: An unconventional spin-ordered phase of matter,” Newton 1, 100162 (2025)。这是一篇近期发表的关于交错磁性的全面综述，可能深入探讨了其微观起源、对称性分类以及在各种材料中的应用前景。
    • [41]-[44] (Yuan, Isidori, Giuli, Cheong) 等研究则侧重于关联驱动的交错磁性，为本研究提出纯电子关联起源的假设提供了重要的理论支撑和先行探索。

  2. Hubbard 模型与强关联理论：

    • [19] J. Hubbard, “Electron correlations in narrow energy bands,” Proceedings of the Royal Society A 276, 238-257 (1963)。Hubbard模型创始性论文，为研究强关联电子系统奠定了理论基础。
    • [1] E. Dagotto, “Correlated electrons in high-temperature superconductors,” Rev. Mod. Phys. 66, 763 (1994)。这篇综述概述了强关联电子系统中的基本物理现象，包括Mott绝缘体和高温超导，为理解本研究背景提供了广阔的视角。

  3. 精确对角化方法与作者团队工作：

    • [53] M. F. Equbal, S. R. Hassan, and M. A. H. Ahsan, “Exact diagonalization study of ground state properties and level statistics in simple and extended hubbard lattice cluster,” Physica Scripta 100, 035964 (2025)。这篇论文是本研究作者团队近期发表的精确对角化工作，显示了他们在ED方法上的专业性，可能与本研究的ED实现细节和计算框架紧密相关。
    • [54] M. F. Equbal, S. R. Hassan, and M. A. H. Ahsan, “Principal component analysis of competing correlations in quarter-filled hubbard models,” arXiv preprint arXiv:2511.12551 (2025)。同样是作者团队的近期ED工作，反映了他们对关联函数的深入分析能力。
    • [55] M. A. H. Ahsan and C. R. Sarma, “Configuration space truncation scheme for hubbard hamiltonian,” Journal of Chemical Sciences 106, 379 (1994) 和 [56] C. R. Sarma and M. A. H. Ahsan, “Electron correlation studies: Rumer basis approach,” International Journal of Quantum Chemistry 60, 147 (1996)。这些早期工作虽然主要关注截断近似，但其中可能包含了作者团队在构建和处理自旋对称基函数方面的经验和技术，这些原理可能被应用于本研究的精确对角化实现。

  4. 几何阻挫物理：

    • [47] M. Alfonso-Moro, V. Guisset, P. David, B. Canals, J. Coraux, and N. Rougemaille, “Geometrical frustration, correlated disorder, and emerging order in a corrugated c60 monolayer,” Phys. Rev. Lett. 131, 186201 (2023)。这篇论文讨论了几何阻挫如何导致关联无序以及非常规序的涌现，为理解本研究中阻挫的作用提供了背景。

对这项工作局限性的评论

尽管本研究通过精确对角化为关联驱动的交错磁性提供了深刻而严格的见解，但作为一项理论计算工作，它也存在一些固有的局限性：

  1. 有限尺寸效应 (Finite Size Effects)：

    • 体系尺寸限制：3x3（9个格点）晶格是精确对角化可处理的最大尺寸之一。然而，它仍然是一个高度受限的有限尺寸系统。在这样的尺度下，系统可能无法完全展现出热力学极限下（即无限大系统）的物理行为。例如，在无限大晶格中可能存在的长程磁序，在小尺寸体系中可能表现为短程关联或涨落增强。相变点的精确位置也可能受到尺寸效应的影响。
    • 外推至热力学极限：将小尺寸系统的精确结果可靠地外推到热力学极限通常具有挑战性，需要进行细致的有限尺寸标度分析，但这并非总能直接适用或完全准确。

  2. 模型简化 (Model Simplification)：

    • 无自旋轨道耦合 (SOC)：论文明确指出，本研究是在“没有单粒子各向异性或自旋轨道耦合”的条件下进行的。虽然这有助于突出纯粹电子关联的作用，但SOC在许多真实材料中对于产生交错磁性是至关重要的。忽略SOC可能限制了将这些发现直接应用于那些SOC效应显著的实际材料。在真实世界中，SOC、晶体场效应、电子关联通常是复杂交织的。
    • 相互作用范围：Hubbard模型及其扩展版仅考虑了在位U和近邻V相互作用。然而，在真实材料中，更长程的库仑相互作用、电子-声子耦合、晶格畸变以及更复杂的能带结构都可能扮演重要角色。这些简化虽然有助于聚焦于核心物理机制，但也可能忽略了真实材料系统中的一些重要细节和耦合效应。

  3. 仅限于基态分析：

    • 有限温度效应：精确对角化主要提供系统的基态性质（T=0K）。然而，许多有趣的物理现象和相变通常发生在有限温度下。本研究未探讨有限温度下的行为，这限制了其对实验观测的直接比较，因为大多数实验都是在非零温度下进行的。
    • 动态性质缺失：研究主要关注静态的自旋关联和局域磁矩，并未深入探讨系统的激发态或动态相关函数。而材料的输运性质、光学响应以及在自旋电子学中至关重要的超快自旋动力学（例如太赫兹自旋振荡）都与动态响应密切相关。

  4. 无代码开源：

    • 复现性挑战：论文未提供任何公开的代码仓库或详细的实现细节，这使得其他研究人员难以直接复现其结果或在此基础上进行扩展研究。在计算物理和计算化学领域，代码的可复现性和可访问性对于推动科学进步和成果验证至关重要。缺乏透明度可能会影响研究成果的广泛采纳和信任。

  5. 缺乏直接实验对比或材料指导：

    • 虽然论文讨论了其发现对材料设计的潜在意义，但并未直接与现有实验数据进行详细对比，也没有明确指出哪些具体材料可能展现出此处发现的关联驱动交错磁性。更具体的材料候选建议和与实验的直接联系将有助于加速理论与实验的相互验证。

尽管存在这些局限性，本研究依然在理论上取得了重要进展，通过精确对角化为理解关联驱动的交错磁性提供了严格的微观证据，并为未来的研究奠定了坚实的基础。

5. 其他你认为必要的补充

工作的深远意义与更广泛的科学背景

这项研究的意义远超出了Hubbard模型的特定数值计算范畴，它对凝聚态物理领域，特别是对新型磁性材料的探索和理解，具有深远的影响：

  1. 拓展磁性物理的分类边界：交错磁性作为一种新颖的磁性相，挑战并拓展了传统的铁磁和反铁磁分类范式。它既具有共线的自旋结构（与反铁磁类似），又表现出在动量空间中巨大的非相对论性自旋劈裂（通常与铁磁或某些非共线磁体相关），同时保持零净磁化。本研究的突破在于，它严谨地证实了在没有自旋轨道耦合或单粒子各向异性这些外部条件的情况下，纯粹的电子关联作用也能在几何阻挫的对称晶格中自发地诱导出鲁棒的交错磁性。这为理解磁性起源提供了一个“纯关联”的全新途径，极大地丰富了我们对磁性序的微观机制的认识，并为凝聚态物理学界提供了一个令人兴奋的新研究方向。

  2. 为下一代自旋电子学提供新平台：交错磁性独一无二的性质——零净磁化与动量空间自旋劈裂的结合，使其成为设计下一代自旋电子器件的理想候选材料。它可能用于开发具有超低能耗、超快速响应的新型自旋器件，例如无杂散场的磁存储器、高效的自旋电流产生器，甚至在太赫兹频率下工作的自旋振荡器。本研究通过揭示载流子掺杂和近邻相互作用（V）作为调控交错磁性关键参数的作用，为实验上实现和优化交错磁性材料提供了重要的理论指导，有助于加速从基础研究到实际应用的转化。

  3. 深化对强关联和几何阻挫物理的理解：几何阻挫和强电子关联是凝聚态物理中两个高度复杂且相互作用的概念，它们常常导致复杂的基态和非常规的量子现象。本研究通过精确对角化这一严格的数值方法，在这些复杂因素的交织下，精确捕捉到了交错磁性关联的形成过程。这本身就是对强关联物理理解的重大贡献，它展示了波动介导的机制如何在阻挫体系中选择性地稳定特定的自旋模式，并克服了阻挫导致的简并性。

  4. 提供高精度基准数据：精确对角化（ED）虽然受限于系统尺寸，但其结果是无近似的“精确”解。因此，本研究为未来开发和验证更适用于大尺寸系统的近似理论方法（如动力学平均场理论、张量网络态方法如DMRG和PEPS、以及量子蒙特卡洛模拟）提供了宝贵的基准数据。这些精确结果可以帮助校准和评估近似方法的准确性和适用范围，推动计算物理方法学的发展。

粒子-空穴不对称性及其重要性

论文中一个特别引人注目的发现是，电子掺杂（N=10）和空穴掺杂（N=8）系统在交错磁性响应上表现出显著的粒子-空穴不对称性。这种不对称性体现在：

  • 响应强度差异：在简单Hubbard模型中，电子掺杂系统（N=10）的Aspin增强幅度明显大于空穴掺杂系统（N=8），表明电子掺杂更能有效地稳定交错磁性关联。
  • 相变行为差异：在扩展Hubbard模型中，电子掺杂系统（N=10）表现出在近邻相互作用V达到临界值时发生急剧的、类似一级相变的复杂行为，而空穴掺杂系统（N=8）则显示出相对平稳连续的演化。

这种粒子-空穴不对称性可能源于晶格能带结构的内在不对称性，以及电子和空穴在几何阻挫晶格中的不同动力学行为。例如，电子和空穴的有效质量、散射机制以及它们与局域磁矩的耦合强度可能存在差异。理解这种不对称性在铜氧化物高温超导体等强关联体系中也普遍存在，这提示了其在更广泛的强关联物理背景下的普遍性和重要性。对于未来精确设计和调控掺杂型交错磁性材料而言，深入理解并利用这种粒子-空穴不对称性至关重要。

卓越的可调控性与实验潜力

本研究识别并强调了载流子浓度（掺杂）和近邻库仑相互作用V作为调控交错磁性的关键参数，这为实验研究和材料设计提供了明确的方向：

  • 载流子掺杂：研究结果清晰地表明，通过增加或减少电子（即N=10或N=8），可以有效地将系统从半填充的、几何阻挫导致的无序态转变为具有鲁棒交错磁性关联的态。这意味着在实验上，可以通过化学掺杂（例如通过替换元素改变价态）或电场诱导的门电压调控（如在场效应晶体管结构中）等手段，来精确地诱导和控制材料中的交错磁性。
  • 近邻库仑相互作用V：V的引入揭示了Hubbard模型的丰富相图和填充依赖的相变行为。特别是电子掺杂系统中的一级相变，表明V可以作为一个临界调控参数，将系统从一个稳定的交错磁性态推向一个竞争性的电荷序或非磁性态。在实验中，V可以通过改变材料的介电常数、施加外部压力（影响键长和电子屏蔽）或通过构建异质结界面来间接调控，为探索和控制交错磁性提供了多维度的方法。此外，在超冷原子实验等量子模拟平台中，光格中的相互作用参数可以被精确控制，这些理论预测有望在这些受控环境中得到直接验证。

未来研究方向的展望

本研究为未来的强关联和交错磁性研究指明了几个有前景和重要的方向：

  1. 更大尺寸体系的探索：尽管ED在小系统上提供了精确结果，但为了接近热力学极限，未来的工作可以结合ED与其他更具扩展性的数值方法，如张量网络态方法（如密度矩阵重正化群DMRG、投影纠缠对态PEPS）或结合动力学平均场理论（DMFT）。这些方法能够处理更大尺度的系统，但需要谨慎处理阻挫引起的符号问题（对于量子蒙特卡洛方法而言）。
  2. 系统引入自旋轨道耦合：鉴于SOC在许多真实材料中的重要性，未来的研究应系统地引入SOC，探索其与关联效应、几何阻挫的协同作用或竞争关系，以更全面地理解和预测真实材料中的交错磁性。
  3. 有限温度和动态性质：利用有限温度ED、量子蒙特卡洛或有限温度DMFT来研究交错磁性相在非零温度下的稳定性、相变行为以及动态响应（如自旋波谱、光学导电性、非线性磁响应），这对于实验观测和实际应用至关重要。
  4. 材料识别与实验验证：基于本研究的理论发现，可以与材料科学专家合作，寻找具有类似晶格结构和电子关联特征的候选材料（例如某些过渡金属氧化物或硫化物），并通过中子散射、角分辨光电子能谱（ARPES）或X射线磁圆二色性（XMCD）等实验技术进行验证。
  5. 探索非平衡态动力学：研究强光场或其他外部驱动如何诱导、激发和调控交错磁性，这对于设计超快自旋电子器件和探索新奇量子现象具有重要意义。
  6. 更复杂的几何结构和拓扑：探索其他几何阻挫晶格（如kagome、三角格、Pyrochlore等）以及非平凡拓扑结构对关联驱动交错磁性的影响，以揭示晶格几何普适性地影响关联驱动磁性的机制。

总而言之，这项工作不仅深化了我们对交错磁性起源的理论理解，也为未来在强关联物理、材料科学和量子技术领域的交叉研究提供了坚实的基础和丰富的研究思路。