自旋-轨道力矩驱动的交错磁体磁动力学

自旋-轨道力矩驱动的交错磁体磁动力学

 

物理学报 · 2026 · 应用磁学

自旋-轨道力矩驱动的交错磁体磁动力学

Magnetization Dynamics of Altermagnet Driven by Spin-Orbit Torque

导读 导读：本文基于原子自旋模拟方法（VAMPIRE），以 d 波交错磁体 RuO₂ 为原型，系统研究了 SOT 作用下交错磁体的磁动力学演化及其在撤去 SOT 后的弛豫过程，并与反铁磁体进行了系统比较。研究发现：SOT 驱动阶段两者行为相似，但交错磁体翻转时间更短；弛豫阶段交错磁体展现出混沌特征（λ_max = 0.0117 > 0, D₂ = 3.582），磁矩最终稳定于非易轴方向；多循环 SOT 脉冲揭示多稳定态与随机弛豫特性。该发现为理解交错磁体的奇异物性提供了新的视角，也为基于交错磁体的存算融合功能型磁性器件设计提供了理论支持。

一、前言背景

交错磁体：第三类基本磁性相

交错磁体（Altermagnet, AM）是继铁磁体和反铁磁体之后发现的第三类基本磁性相。其实空间呈现共线且补偿的磁矩排布——净磁矩严格为零，因而无杂散场；而动量空间则表现为自旋能带劈裂，兼具铁磁体中自旋极化输运的优势。

交错磁体的核心特征是"交错自旋劈裂效应"——一种不依赖于自旋-轨道耦合的非相对论性效应，源于晶体对称性所关联的磁亚晶格所处各向异性晶体场环境的差异。

Šmejkal 等人于 2022 年建立了交错磁体的自旋群对称性分类框架 [Phys. Rev. X 12, 031042 & 040501]，将交错磁体按自旋劈裂的动量空间对称性分为 d 波、g 波、i 波等类型，其中 RuO₂ 是典型的 d 波交错磁体。

交错磁体同时具备反铁磁体的零杂散场 + 太赫兹本征磁动力学频率，以及铁磁体的自旋劈裂能带特征——被视为实现高密度、高速度且易读出的磁随机存储器（MRAM）的关键候选材料。

交错磁体的非平庸磁电输运性质

交错自旋劈裂能带结构赋予交错磁体一系列非平庸的磁电输运性质，包括由自旋劈裂能带控制的反常横向输运效应 [Feng et al., Nat. Electron. 5, 735 (2022)]、自旋流-电荷流相互转化 [Bai et al., Phys. Rev. Lett. 130, 216701 (2023)]、巨磁阻效应和隧穿磁电阻效应 [Šmejkal et al., Phys. Rev. X 12, 011028 (2022)] 等。

由于同时具备反铁磁体无杂散场和高磁动力学本征频率的优势，以及铁磁体自旋能带劈裂的特征，交错磁体被视为实现高密度、高速度且易读出的磁随机存储器的关键候选材料。

SOT-MRAM 与交错磁体磁动力学

磁性隧道结（MTJ）是 MRAM 的基本功能单元——通过控制自由层和参考层磁序的平行或反平行排列，实现高低阻态的切换，即隧穿磁电阻效应。

自旋-轨道力矩（SOT）被视为下一代 MRAM 的关键写入技术——在保持高器件耐久度的同时，可进一步降低写入电流密度。

SOT 驱动交错磁体磁序翻转已在 Mn₅Si₃ [Han et al., Sci. Adv. 2024]、CrSb [Zhou et al., Nature 638, 645 (2025)] 和 RuO₂ [Zhang et al., Nat. Commun. 16, 5646 (2025)] 等体系中得到实验验证。

然而，关于 SOT 引起的交错磁体磁动力学行为的理论研究尚不充足，尤其是其与反铁磁体的差异仍有待揭示。

本文研究目标与创新点

基于原子自旋模拟方法（VAMPIRE 软件包），以 d 波交错磁体 RuO₂ 为原型构建原子自旋模型，系统研究 SOT 作用下交错磁体的磁动力学演化及其在撤去 SOT 作用后的弛豫过程。

与相同条件下的反铁磁体磁动力学行为进行系统比较，揭示各向异性交换相互作用对磁动力学的本质影响。

创新点：首次揭示交错磁体在 SOT 撤除后弛豫过程中的混沌动力学特征——这是传统反铁磁体中从未被报道的新现象。

二、原子自旋模型

以 RuO₂ 为原型的原子自旋模型

本文采用的交错磁体原子自旋模型以 d 波交错磁体 RuO₂ 为原型进行构建。蓝色与红色球体表示磁性 Ru 原子，具有反平行排列的磁化方向，构成两套磁矩补偿的磁亚晶格。

非磁原子（O）在磁性 Ru 原子周围构成八面体结构——围绕磁性原子的 8 个 O 原子到中心 Ru 原子的距离不完全相同，磁性原子所处晶体场呈现各向异性。

晶体场的各向异性导致同一磁亚晶格内沿不同晶向的交换相互作用不相等：沿 [110] 晶向的 J₅（原子 1 与原子 2 之间）与沿 [1̄10] 晶向的 J₄（原子 3 与原子 4 之间）存在显著差异。

在传统反铁磁体中，晶体场对称性使这两类交换相互作用通常相等；交错磁体中的这种不等性——即各向异性交换相互作用——是区分两类磁体的核心特征。

两套磁亚晶格周围的八面体结构由旋转对称操作相互联系，各向异性交换相互作用也可通过相同的旋转对称操作互相映射。

图 1 | 交错磁体的原子自旋模型。(a) 交错磁体晶格结构和交换相互作用，蓝/红球体表示反平行磁矩的两套磁亚晶格；(b) 磁亚晶格内的交换相互作用 J₄ 和 J₅；(c) 磁亚晶格间的交换相互作用映射关系。

模型参数与物理合理性

交换相互作用参数取自第一性原理计算 [Šmejkal et al., Phys. Rev. Lett. 131, 256703 (2023)]：J₁=3×10⁻²² J/原子对，J₂=−8×10⁻²³ J/原子对，J₃=3×10⁻²³ J/原子对，J₄=1×10⁻²² J/原子对，J₅=−2×10⁻²² J/原子对。

由于第一性原理计算得到的 J₄ 和 J₅ 差异较小，其对磁矩排布造成的影响易被亚晶格间的交换相互作用 J₂ 掩盖。为揭示交错磁体中可能存在的新奇磁动力学行为，在保证体系无外场、无外力矩条件下仍稳定为反平行排列的前提下，对 J₂ 的数值在合理范围内进行优化调整。

磁晶各向异性能 K_u=1.6×10⁻²³ J/atom，原子磁矩 2.83 μ_B，磁易轴沿 [001] 晶向，晶胞常数 a=b=4.493 Å, c=3.106 Å。模型边界采用周期性边界条件。

包含 SOT 作用的 LLG 方程

SOT 器件的典型三层结构：下层重金属自旋源层（Pt, W, Ta 等）→ 中间交错磁体层 → 上层覆盖层。电流在自旋源层中流动，通过自旋霍尔效应产生沿膜面法向流动的自旋流。

注入的自旋流携带的自旋极化方向 σ̂ 与电荷流方向和自旋流方向均相互垂直。注入的自旋流主要以类阻尼力矩（DL-SOT）和类场力矩（FL-SOT）两种形式作用于两套亚晶格磁矩。

考虑到非绝缘性交错磁体层中电流分流产生的层内自旋流所携带的自旋极化方向与交错磁体的 Néel 矢量平行，根据自旋力矩中叉乘形式的关系，其在亚晶格磁矩上产生的力矩为零，不会对磁动力学行为产生实质性影响。

公式 (1) — LLG 方程：m₁,₂ 为两套亚晶格磁矩，γ 为旋磁比，α 为 Gilbert 阻尼因子，σ̂ 为注入自旋流的极化方向单位矢量，ζ_FL 和 ζ_DL 分别为类场和类阻尼有效场强度系数。

系统哈密顿量：包含磁性原子间的交换相互作用（J_ij 项）和磁晶各向异性能（K_u 项），求和遍历所有磁性原子对。

有效场 H_eff 由哈密顿量对磁矩的变分给出：包含磁晶各向异性场和交换场两部分，是驱动磁矩进动的核心物理量。

图 2 | SOT 激发交错磁体磁动力学行为示意图。下层自旋源层（重金属）中电荷流通过自旋霍尔效应转化为自旋流注入交错磁体层，SOT 作用于磁矩激发磁动力学响应。

VAMPIRE 模拟方法与数值细节

采用由约克大学 Evans 等人开发的 VAMPIRE 原子自旋模拟软件包开展计算 [J. Phys. Condens. Mater. 26, 103202 (2014)]，采用 Heun 数值积分格式求解 LLG 方程。

模拟超胞 4.493 nm × 4.493 nm × 3.106 nm，时间步长 Δt = 2×10⁻¹⁷ s，总模拟时间 200 ps。初始条件：两套亚晶格磁矩沿 [001] 易轴方向反平行排列，注入自旋极化方向 σ̂ ∥ [110]，模拟温度 T = 0 K。

SOT 脉冲时序：直流 SOT 持续 50 ps 后撤除，弛豫阶段持续 150 ps。多循环研究中，10/100 次 SOT 脉冲间隔足够长以确保每次脉冲前磁矩已充分弛豫。

三、SOT 驱动磁动力学

SOT 作用下的磁动力学行为

初始条件设定：两套亚晶格磁矩沿 [001] 易轴方向反平行排列，注入自旋极化方向 σ̂ ∥ [110]。模拟表明初始磁矩排列偏离易轴方向不会对磁动力学行为产生定性影响。

对类阻尼有效场强度 μ₀H_DL 进行参数扫描，发现临界阈值 μ₀H_DL^crit = 0.45 T——该阈值两侧交错磁体展示出截然不同的磁动力学行为。

低场区域（μ₀H_DL < 0.45 T）：类谐振子振荡

当类阻尼有效场低于临界阈值时，SOT 强度不足以驱动磁矩跨越 (001) 面——磁矩在类阻尼力矩驱动下偏离易轴方向，产生类谐振子振荡。

物理机制：磁矩偏离易轴方向时，磁晶各向异性场随偏离角增大而增强，形成类似于弹簧振子中回复力的作用。振荡频率约 0.2 THz，与交错磁体理论预期的本征磁动力学频率量级一致。

在 0.45 T 以下的低场范围内，振荡频率随类阻尼有效场的变化无明显偏移——这是因为材料固有的磁晶各向异性场在振荡过程中占据主导地位。

在 Gilbert 阻尼的持续作用下，振荡幅度逐渐衰减，亚晶格磁矩最终稳定在由磁晶各向异性与类阻尼有效场共同决定的倾斜方向。

高场区域（μ₀H_DL > 0.45 T）：磁矩跨越翻转

当类阻尼有效场超过临界阈值时，SOT 能量足以克服磁晶各向异性能垒——亚晶格磁矩跨越 (001) 面，在单位球面上进行大振幅圆周振荡，最终翻转至与注入自旋极化平行的 [110] 晶向。

两阶段动力学特征：0-15 ps 区间内磁矩在 SOT 驱动下快速偏离初始方向，进动幅度和频率均较大；15-50 ps 区间内磁矩已接近翻转终点，进动幅度快速衰减。

FFT 分析显示，两个时间区间内的进动频率均随类阻尼有效场的增大而明显升高——这表明在磁矩翻转过程中，外加类阻尼有效场占据主导作用。

图 3 | SOT 作用下交错磁体和反铁磁体的磁动力学行为。(a) 交错磁体磁动力学频率随 SOT 有效场的变化；(b) 交错磁体亚晶格磁矩在 0.4 T 和 0.5 T SOT 有效场作用下的时域谱线；(c) 反铁磁体磁动力学频率随 SOT 有效场的变化；(d) 反铁磁体亚晶格磁矩在 0.4 T 和 0.5 T SOT 有效场作用下的时域谱线。

AM vs AFM：SOT 驱动阶段的定量对比

将 J₄ 和 J₅ 设定为等同正值后，交错磁体模型转化为反铁磁体模型——在相同模拟条件下进行对比计算。从频域信息来看，两者在 SOT 驱动下呈现出相似的动力学响应特征。

关键差异：当类阻尼有效场超过翻转阈值后，交错磁体表现出显著更短的翻转时间——各向异性交换相互作用对亚晶格磁矩施加了额外的力矩作用，显著加快了翻转过程。

反铁磁体在 15 ps 之后保持进动幅度的缓慢衰减（Gilbert 阻尼主导），而交错磁体在相同时段内进动幅度快速衰减。

这一发现表明：基于交错磁体的 SOT 存储器件有望实现较反铁磁体器件更快的写入速度——这是交错磁体在器件应用中的一个重要优势。

四、弛豫过程混沌特征

SOT 作用后的弛豫行为：混沌特征的发现

实验设计：选取能充分翻转磁矩的类阻尼有效场强度 μ₀H_DL = 1.5 T 与作用时间 t = 50 ps，得到交错磁体与反铁磁体在相同条件下的弛豫行为。

反铁磁体的弛豫（期望行为）：亚晶格磁矩在磁晶各向异性场的作用下首先围绕 [110] 晶向进行球面圆周进动，进动轴逐渐由 [110] 转变为易轴 [001] 方向，最终在 Gilbert 阻尼作用下稳定至易轴方向——全过程平滑、可预测。

交错磁体的弛豫（意外发现）：短暂围绕 [110] 晶向进行圆周进动之后，磁矩运动路径迅速变得紊乱——时域谱线和球面轨迹均呈现明显的非周期性与不可预测性。

这一发现是本工作的核心创新点：交错磁体中各向异性交换相互作用使 SOT 撤除后的弛豫过程从简单的阻尼振荡转变为复杂的混沌动力学行为。

图 4 | SOT 作用后弛豫过程反铁磁体和交错磁体的磁动力学行为。(a) 反铁磁体亚晶格磁矩在 SOT 作用后弛豫过程的时域谱线；(b) 反铁磁体亚晶格磁矩的弛豫轨迹；(c) 交错磁体亚晶格磁矩在 SOT 作用后弛豫过程的时域谱线；(d) 交错磁体亚晶格磁矩的弛豫轨迹。

混沌定量诊断：λ_max 与 D₂

对 100-200 ps 区间内的磁矩分量 m₁z(t)（共 1000 个数据采样点）进行混沌动力学定量分析。嵌入维数通过假近邻法确定，时间延迟使用平均互信息法进行优化。

最大李雅普诺夫指数 λ_max = 0.0117：采用 Rosenstein 方法计算。λ_max > 0 是混沌动力学的严格判据——相邻轨线指数发散，系统对初始条件极端敏感。

关联维数 D₂ = 3.582：采用 Grassberger-Procaccia 方法计算。D₂ 为非整数且有限——表明系统具有低维混沌特征，其动力学行为由少数几个自由度决定，但具有不可预测性。

λ_max > 0 与 D₂ 为非整数这两个指标共同确认了弛豫过程的混沌本质——这是首次在磁体系中通过各向异性交换相互作用系统揭示的混沌动力学行为。

在经历一段进动之后，亚晶格磁矩最终收敛至球面坐标 (0.24, −0.90, 0.35)，明显偏离易轴方向——暗示了交错磁体中各向异性的交换相互作用可能导致不同于易轴方向的稳态或亚稳态的出现。

最大李雅普诺夫指数 λ_max 的定义：δm(t) 为初始间距 δm(0) 的两条相邻轨线在 t 时刻的间距。λ_max > 0 表征轨道指数发散，是混沌系统的严格判据。

弛豫行为的影响因素：参数扫描

为系统研究混沌弛豫过程的影响因素，在更多参数组合下进行弛豫行为的数值计算。三个关键参数：各向异性交换相互作用（J₄, J₅）、Gilbert 阻尼因子 α、磁晶各向异性能 K_u。

参数扫描策略：保持 J₄+J₅ 不变以确保其余交换项与磁晶各向异性项在哈密顿量中的相对贡献不变，仅改变 J₄ 与 J₅ 的比例（−J₄/J₅ 和 −J₅/J₄，从 1.2 开始）。

各向异性交换相互作用的影响：在参数扫描范围内改变 J₄ 和 J₅ 比例后，交错磁体在弛豫阶段仍表现出正的 λ_max，即混沌特征具有一定鲁棒性。λ_max 随比例增大而逐渐减小，混沌程度降低。

Gilbert 阻尼因子 α 的影响：当 α < 0.01（低阻尼）时，磁矩进动的能量耗散较弱，弛豫过程中保持较高的 λ_max，对应更显著的混沌特征；随着 α 增大，阻尼增强，混沌特征逐渐消失（λ_max 由正转负）。

磁晶各向异性 K_u 的影响：在不同 K_u 条件下，交错磁体在弛豫过程中均表现出明显的混沌轨迹特征。K_u 增大时，进动幅度显著增大，同时最终稳定状态中 z 分量的绝对值也随之增大。

图 5 | 交错磁体亚晶格磁矩弛豫过程 λ_max 的影响因素。(a) λ_max 与各向异性交换相互作用 J₄, J₅ 的关系；(b) λ_max 与 Gilbert 阻尼因子 α 的关系。

图 6 | 不同磁晶各向异性场条件下 SOT 作用后交错磁体弛豫过程的动力学行为。(a) K_u = 8×10⁻²⁴ J/atom 时的交错磁体亚晶格磁矩弛豫过程；(b) K_u = 3.2×10⁻²³ J/atom 时的交错磁体亚晶格磁矩弛豫过程。

参数调控的物理意义与实验前景

上述规律揭示了各向异性交换相互作用和 Gilbert 阻尼因子对混沌程度的调控作用——交换相互作用对温度和晶格应力敏感 [Kunes et al., J. Phys. Soc. Jpn. 74, 1408 (2005)]，Gilbert 阻尼因子亦可通过材料生长工艺进行调控 [Gilmore et al., Phys. Rev. Lett. 99, 027204 (2007)]。

交错磁体最终的稳定取向由磁晶各向异性能与磁性原子间交换相互作用能之间的竞争共同决定。

通过合理选择材料体系，使其磁相关能量尺度显著高于室温下的热涨落能量，仍有望在室温条件下保持稳定的磁有序结构，从而观测到本文所预测的交错磁体动力学特性。

五、脉冲响应与多稳定态

多循环 SOT 脉冲实验设计

对交错磁体和反铁磁体施加周期性 SOT 脉冲，比较两者在多次写入-弛豫循环中的动力学行为差异。

反铁磁体的行为：10 个脉冲循环中，亚晶格磁矩的动力学行为保持高度一致——弛豫终态由磁晶各向异性决定且始终位于易轴 [001] 方向，表现出完美的可重复性。

交错磁体的行为：每一次脉冲作用后均弛豫至与前一循环不同的位置——表现出对初始条件的敏感依赖性，这正是混沌系统的典型特征之一。

图 7 | 交错磁体和反铁磁体对周期性 SOT 脉冲的响应。(a) 交错磁体和反铁磁体亚晶格磁矩对 10 个 SOT 脉冲响应的时域谱线；(b) 交错磁体亚晶格磁矩在 100 个 SOT 脉冲作用下的弛豫最终态。

100 次脉冲统计：结构化的多稳定态分布

进一步将脉冲作用次数扩展至 100 次，并统计每个循环结束后亚晶格磁矩的稳定位置。

尽管交错磁体的稳定位置在不同循环间存在差异，但并非在球面上均匀随机分布，而是呈现出特定方向的聚集。

这一现象为前文提出的观点提供了更加明确的支持：交错磁体中各向异性的交换相互作用可能为体系引入多个不同于易轴方向的稳态或亚稳态，使其弛豫终态呈现结构化的多重分布。

这些稳态或亚稳态的产生源于亚晶格内部磁性原子之间各向异性交换相互作用的存在——与常规反铁磁体不同，交错磁体中同一亚晶格内的磁矩不再严格保持平行排列，从而引入了由交换相互作用产生的附加力矩项。

多稳定态的应用前景：存算融合器件

交错磁体在 SOT 脉冲作用后可在多个稳定态之间"随机"弛豫——这正是随机性或概率性磁性元件所期望的关键特征。

传统存储功能：通过调控 SOT 脉冲参数使系统在规则区域工作，实现确定性翻转，完成"0"和"1"的信息存储。

概率计算：利用多稳定态之间的随机弛豫特性，实现 p-bit（概率比特）——其输出为概率分布而非确定的二进制值，是概率计算和贝叶斯推理的基本单元。

神经形态计算与随机逻辑：多稳定态可对应神经网络的多个权重状态，随机弛豫可模拟神经元的随机发放行为。

由此可见，基于交错磁体的器件在实现传统存储功能的同时，还具备开展概率计算、神经形态计算与随机逻辑等新型功能的潜力，有望成为兼具计算与存储能力的功能型存算融合磁性器件。

六、总结展望

核心结论一：SOT 驱动下的动力学差异

交错磁体与反铁磁体在 SOT 作用下呈现出类似的动力学响应特征。然而，当类阻尼有效场达到磁矩翻转阈值（μ₀H_DL^crit = 0.45 T）后，交错磁体表现出显著更短的翻转时间。这表明基于交错磁体的自旋轨道力矩存储器件有望实现较反铁磁体器件更快的写入速度。

核心结论二：具有混沌特征的弛豫行为

在 SOT 脉冲撤去后的弛豫过程中，交错磁体展现出混沌特征——最大李雅普诺夫指数 λ_max = 0.0117 > 0，关联维数 D₂ = 3.582（非整数）。λ_max 可通过调控各向异性交换相互作用参数及 Gilbert 阻尼因子进行有效调节。该可调控性为实验条件下克服不可避免的噪声干扰、提升混沌特征的可探测性提供了理论支撑。

核心结论三：多稳定态与随机弛豫特性

多次循环 SOT 写入的模拟结果显示，交错磁体在每次脉冲作用后的动力学行为呈现出多稳定态与随机弛豫的特征。这表明交错磁体在传统存储功能之外，还具备实现概率计算、神经形态计算和随机逻辑等应用的潜在优势，在构建存算融合功能磁性器件方面展示出可能性。

展望与未来方向

实验验证：通过时间分辨磁光克尔效应（TR-MOKE）、铁磁共振（FMR）和自旋泵浦等时间分辨磁动力学测量技术，在室温条件下直接观测交错磁体的 SOT 驱动动力学和混沌弛豫行为。

材料优化：通过应变工程、界面效应和化学掺杂来调控各向异性交换、Gilbert 阻尼和磁晶各向异性，优化器件性能。

温度效应：将有限温度引入 VAMPIRE 模拟，通过朗之万热场描述热涨落，系统研究温度对混沌弛豫的阈值、λ_max 和多稳定态分布的影响。

新型计算范式：基于交错磁体混沌弛豫和多稳定态特性，探索其在储层计算、伊辛机和随机计算等新型计算范式中的应用。

主要参考文献

核心引用文献

[1] 曹彦彰, 韩磊, 潘峰, 王钰言, 宋成. 物理学报 75, 060709 (2026) — 本工作：SOT 驱动的交错磁体磁动力学。

[2] Šmejkal L, Sinova J, Jungwirth T. Phys. Rev. X 12, 031042 (2022) — 交错磁体的自旋群对称性分类。

[3] Šmejkal L, Sinova J, Jungwirth T. Phys. Rev. X 12, 040501 (2022) — 超越铁磁与反铁磁的第三类磁性相。

[4] Šmejkal L, Marmodoro A, et al. Phys. Rev. Lett. 131, 256703 (2023) — RuO₂ 交换相互作用的第一性原理计算。

[5] Evans R F L, Fan W J, et al. J. Phys. Condens. Mater. 26, 103202 (2014) — VAMPIRE 原子自旋模拟软件包。

[6] Han L, Fu X Z, Peng R, et al. Sci. Adv. 10, eadn0479 (2024) — Mn₅Si₃ 中 SOT 驱动的 Néel 矢量翻转。

[7] Zhou Z Y, Cheng X K, Hu M L, et al. Nature 638, 645 (2025) — CrSb 中 SOT 驱动的磁序翻转。

[8] Zhang Y C, Bai H, Dai J K, et al. Nat. Commun. 16, 5646 (2025) — RuO₂ 中 SOT 驱动的磁序翻转。

[9] Song C, Bai H, et al. Nat. Rev. Mater. 10, 473 (2025) — 交错磁体自旋电子学综述。

[10] Wolf A, Swift J B, et al. Physica D 16, 285 (1985) — 最大李雅普诺夫指数的计算方法。

[11] Grassberger P, Procaccia I. Phys. Rev. Lett. 50, 346 (1983) — 关联维数的计算方法。

[12] Roy K, Sengupta A, Shim Y. J. Appl. Phys. 123, 210901 (2018) — 概率计算与神经形态计算中的自旋电子学器件。

DFT 计算 Tips

【DFT Tip 1】交换耦合参数 J_ij 的 DFT 提取

从 DFT 提取交换耦合参数是原子自旋模拟的第一步，也是最容易出错的步骤。

四态法（four-state method）：对不同磁构型（铁磁、反铁磁、交错磁等）计算总能量，解线性方程组得到 J_ij。

关键：至少需要 4 种磁构型（铁磁 + 3 种反铁磁排列），推荐 6-8 种以降低误差。

常见错误：只用 2 种磁构型（铁磁和反铁磁）→ 无法区分 J₁、J₂、J₃ 的贡献 → 模型参数误差可达 50%。

建议：提取的 J 值必须能再现 DFT 计算的磁基态能量和磁矩方向，否则不可靠。

【DFT Tip 2】磁晶各向异性能 K_u 的计算

磁晶各向异性能 K_u 对 LLG 动力学的弛豫终态有决定性影响——计算误差直接影响磁矩稳定方向。

VASP：LNONCOLLINEAR=.TRUE. + LSORBIT=.TRUE.，计算不同磁化方向（如 [001], [100], [110]）的总能差。

关键：SOC 计算需要比 Non-SOC 更密的 k 网格（至少 2 倍），否则 K_u 的误差可达 50%。

常见错误：用 Non-SOC 的 k 网格直接跑 SOC → K_u 值被系统性低估 → 磁矩弛豫方向错误。

建议：K_u 的符号决定了易轴方向——正值为易轴，负值为易面。交错磁中 K_u 的面内各向异性同样重要。

【DFT Tip 3】Gilbert 阻尼因子 α 的理论估算

Gilbert 阻尼因子 α 是 LLG 方程中决定磁矩进动衰减速度的关键参数——直接从 DFT 计算 α 是当前研究的前沿。

方法：基于线性响应理论，从自旋分辨的态密度和自旋-轨道耦合矩阵元计算 α——VASP 中通过 LDAU 和 LORBIT 辅助。

实践：α 的理论值通常比实验值小 1-2 个数量级——因为 DFT 无法描述 magnon-magnon 散射和电子-声子散射等额外阻尼通道。

常见错误：直接将 DFT 计算的 α 用于 LLG 模拟 → 阻尼太弱 → 磁矩进动永不衰减。

建议：以实验 α 为参考（如 RuO₂ 的 α~0.01-0.05），理论 α 可作为下限，参数扫描 α 以评估阻尼对动力学的影响。

【DFT Tip 4】SOT 有效场参数 ζ_DL 和 ζ_FL 的提取

SOT 有效场强度 ζ_DL 和 ζ_FL 是原子自旋模拟中 SOT 力矩的输入参数——直接决定磁矩翻转的临界电流密度。

DFT 提取方法：计算交错磁体/重金属异质结中的自旋霍尔电导 σ_SH，结合界面自旋透明度 T_int 得到 ζ = (ħ/2e) σ_SH T_int。

关键：界面自旋透明度 T_int 对界面原子结构极其敏感——界面弛豫 vs 不弛豫，T_int 可差 3-5 倍。

常见错误：用体相自旋霍尔角代替界面有效场 → ζ 被高估或低估 → 翻转阈值预测错误。

建议：必须显式计算异质结（至少 4 层 AM + 6 层 HM），弛豫界面结构后再计算非平衡自旋密度。

【DFT Tip 5】VAMPIRE 原子自旋模拟的参数映射

VAMPIRE 需要将 DFT 的交换耦合参数（通常以 J/magnetic atom 为单位）映射到模拟参数中。

关键：注意 J 的符号——正 J 为铁磁耦合，负 J 为反铁磁耦合。交错磁中 J₄ 和 J₅ 的符号/大小不等是区分 AM 和 AFM 的关键。

若 DFT 计算的 J₄/J₅ 差异太小（< 5%），在模拟中可能被 J₂ 掩盖——此时可适当放大 J₄/J₅ 差异以揭示交错磁特有的动力学效应。

常见错误：将 DFT 的 J 值直接代入，不检查磁基态是否与 DFT 一致 → 模拟的基态可能与 DFT 完全不同。

建议：先在 VAMPIRE 中弛豫无外场条件下的磁矩，确认基态能量最低的磁序与 DFT 一致。

【DFT Tip 6】李雅普诺夫指数的计算陷阱

最大李雅普诺夫指数 λ_max > 0 是混沌动力学的严格判据，但计算中容易产生伪正值的假阳性。

关键参数：(1) 嵌入维数 m 需通过假近邻法确定——m 太小导致轨线折叠，伪混沌；m 太大导致噪声放大。(2) 时间延迟 τ 需通过平均互信息法确定。

常见错误：直接用原始数据不经过状态空间重构 → λ_max 无物理意义；时间序列太短（< 500 点）→ λ_max 的统计误差大。

建议：至少用 1000 个数据点，用 Rosenstein 和 Kantz 两种方法交叉验证 λ_max，用替代数据法（surrogate data）排除线性随机过程的假阳性。

【DFT Tip 7】混沌 vs 随机：关联维数 D₂ 的判别

关联维数 D₂ 的非整数值是混沌系统的特征——但 D₂ 的计算对噪声极为敏感。

实践：Grassberger-Procaccia 方法计算 D₂ 时，ln C(r) vs ln r 曲线的线性区域必须跨越至少 1 个数量级。

陷阱：如果线性区域太窄（< 0.5 个数量级），D₂ 的拟合值不可靠——可能是噪声而非混沌。

建议：D₂ 的计算应结合替代数据检验——如果替代数据的 D₂ 与原数据 D₂ 在误差范围内重合，则原数据可能是线性随机过程而非混沌。

【DFT Tip 8】交错磁体对称性识别：ISYM 与自旋群

VASP 默认使用磁空间群对称性（ISYM=2），但交错磁的对称性需要自旋群描述——两者在某些情况下不一致。

实践：在 INCAR 中设置 ISYM=0 可禁用对称性（代价是计算量增加），确保 VASP 不会错误地施加反铁磁对称性。

但 ISYM=0 会增加 k 点约减失效——建议在 Relax 阶段用 ISYM=0，Static 阶段用 ISYM=2（如果对称性一致）。

验证：用 pymatgen 或 FINDSYM 检查 VASP 输出的对称性，确保 CONTCAR 的磁序与期望一致。

建议：如果 VASP 自动将 AM 结构识别为 AFM 并施加了反演对称性，能带劈裂会消失——这是最常见的"假阴性"。

【DFT Tip 9】多稳定态搜索：从 DFT 到 VAMPIRE

SOT 作用后交错磁体弛豫到非易轴方向的稳定态——这些稳定态是 DFT 基态搜索中通常不会出现的。

实践：在 VAMPIRE 中设置不同的初始磁化方向，弛豫后记录所有不同的终态磁矩方向——这些就是多稳定态。

DFT 验证：选取 VAMPIRE 预测的稳定态磁矩方向，在 DFT 中进行 constrained magnetic moment 计算，验证这些态是否真的是（亚）稳态。

建议：多稳定态的能量差通常在 1-10 meV/atom 量级——DFT 收敛精度必须足够高（EDIFF < 1e-6 eV）才能可靠区分。

【DFT Tip 10】从计算到实验：可验证的预测

理论计算的最终目标是提出可被实验验证的预测。SOT 驱动混沌动力学研究中，最重要的可验证预测是：

(1) 翻转时间：交错磁比反铁磁快 30-50%，可通过时间分辨电输运测量验证。

(2) 弛豫终态：非易轴方向的稳定态，可通过 XMCD-PEEM 直接成像。

(3) 多稳定态：100 次脉冲后的终态分布，可通过统计 TR-MOKE 信号验证。

建议：在论文中明确列出"可验证的预测"——这是连接理论计算和实验验证的桥梁，也是审稿人最看重的。

九、知识扩展

【知识扩展 1】原子自旋模拟（Atomistic Spin Modeling）

【理论解释】原子自旋模拟是介观尺度的磁动力学方法——将每个磁性原子视为一个经典磁矩矢量，通过 LLG 方程描述其时间演化。体系哈密顿量包含交换相互作用、磁晶各向异性、偶极相互作用和外场项。与微磁学（连续介质）不同，原子自旋模拟能分辨原子尺度的磁矩非共线性和局部涨落。

【方法比较】DFT（第一性原理）：精度最高，但限于 ~1000 原子/ps 时间尺度；原子自旋模拟：精度受限于输入参数（J, K_u, α），但可处理 ~10⁶ 原子/ns 时间尺度；微磁学：速度最快，但无法描述原子尺度的磁矩非共线性。交错磁中各向异性交换相互作用是原子尺度效应——必须用原子自旋模拟而非微磁学。

【经典参考】Evans et al., J. Phys. Condens. Mater. 26, 103202 (2014) — VAMPIRE 软件包；Skubic et al., J. Phys. D 41, 134002 (2008) — 原子自旋模拟综述；Antropov et al., Phys. Rev. B 54, 1019 (1996) — 原子自旋动力学的基础理论。

【迁移能力】原子自旋模拟可用于任何磁有序体系的动力学研究：包括铁磁/反铁磁共振、skyrmion 动力学、畴壁运动、自旋波传播和热辅助磁化翻转。特别适合研究交错磁中由各向异性交换相互作用引起的非平庸动力学行为。

【知识扩展 2】SOT（自旋-轨道力矩）驱动的磁化翻转

【理论解释】SOT 是自旋-轨道耦合在异质结界面产生的力矩效应。当电流通过重金属层（如 Pt, W, Ta）时，自旋霍尔效应产生垂直于电流方向的自旋流注入磁性层，注入的自旋极化电子在磁性原子磁矩上施加力矩。SOT 包含类阻尼力矩（DL-SOT ∝ m×(σ×m)）和类场力矩（FL-SOT ∝ m×σ）两个分量。

【方法比较】STT（自旋转移力矩）vs SOT：STT 需要电流通过磁性层，写入电流密度高（~10⁷ A/cm²）；SOT 的写入电流不通过磁性层，写入电流密度可降低 1-2 个数量级。SOT 的三端器件结构将读写路径分离，显著提升器件耐久度。

【经典参考】Manchon et al., Rev. Mod. Phys. 91, 035004 (2019) — SOT 综述；Miron et al., Nature 476, 189 (2011) — SOT 磁化翻转实验；Liu et al., Science 336, 555 (2012) — SOT 在 HT-MTJ 中的应用。

【迁移能力】SOT 不仅用于 MRAM 写入，还可用于：自旋波激发、skyrmion 驱动、自旋神经形态计算和自旋逻辑器件。在交错磁体中，SOT 的研究为理解自旋-轨道耦合在补偿磁体系中的角色提供了独特视角。

十、科研经验

【科研经验 1】为什么 DFT 计算的交换耦合参数误差大，导致磁动力学预测不可靠？

【问题】用 DFT 提取的 J_ij 参数跑 VAMPIRE，发现磁矩翻转阈值、弛豫终态与预期完全不同——参数误差的放大效应。

【原因】(1) DFT 总能量差（~1-10 meV/atom）在 J_ij 提取中被放大——k 点不够密或 ENCUT 不够高都会导致 J_ij 误差；(2) 四态法中磁构型选择不当——某些构型能量接近简并，导致线性方程组病态；(3) 远邻 J 的截断误差——只取到最近邻 J 而忽略远邻，误差累积。

【解决方案】(1) 提高 DFT 计算精度至 EDIFF < 1e-7 eV，ENCUT ≥ 1.3×ENMAX；(2) 使用更多磁构型（6-8 种），通过最小二乘法而非直接求解；(3) 验证截断半径——逐步增加 J 的范围，直到模拟的磁子谱与 DFT 声子谱一致。

【建议】交换参数提取是"垃圾进垃圾出"的典型——如果 DFT 精度不够，后续所有原子自旋模拟结果都是不可靠的。在开始 VAMPIRE 之前，至少花 50% 的时间验证 J_ij 的收敛性。

【科研经验 2】为什么原子自旋模拟的弛豫过程与实验预期不符？

【问题】VAMPIRE 模拟中磁矩弛豫到非易轴方向——但实验测量（如 XMCD）显示磁矩沿易轴方向。为什么不一致？

【原因】可能的原因包括：(1) 模拟温度为 T=0 K，而实验在有限温度下热涨落可能使磁矩逃离亚稳态回到易轴方向；(2) 模拟中 Gilbert 阻尼 α 与实际值不符——α 太小导致亚稳态不易逃离；(3) 实验的时间尺度（秒）远大于模拟（ns），在更长时间尺度上磁矩可能通过热激活回到易轴方向。

【解决方案】(1) 引入朗之万热场（T > 0 K）跑更长时间的模拟，看亚稳态是否在有限温度下不稳定；(2) 参数扫描 α 和 T，给出亚稳态寿命的相图；(3) 计算亚稳态和基态之间的能垒，通过 Arrhenius 公式估算室温下的寿命。

【建议】不要因为模拟与实验不一致就否定模拟结果——亚稳态在有限温度下可能存在 nm 到 μs 量级的寿命，在模拟时间尺度（ns）内表现为"稳定"，但在实验时间尺度（s）内可能已弛豫到基态。

 

【进一步计算 1】有限温度磁动力学：朗之万热场

为什么值得算：当前模拟为 T=0 K，实际器件在室温（300 K）工作——热涨落对混沌动力学和亚稳态寿命的影响至关重要。

能回答的问题：室温下混沌特征是否仍然存在？亚稳态在室温下的寿命是多长？λ_max 在什么温度下由正转负？

适合体系：RuO₂（T_N > 300 K）和 MnTe（T_N=307 K）等室温磁有序交错磁体。

输入：VAMPIRE + 朗之万热场（温度参数 T）+ 多次独立模拟（至少 50 次取平均）。

 

【进一步计算 2】更多交错磁体材料体系

为什么值得算：当前研究以 RuO₂ 为原型——但 RuO₂ 的交错磁性近期受到争议。需要扩展到更多公认的交错磁体。

能回答的问题：MnTe、CrSb、Mn₅Si₃ 等体系中混沌动力学是否普遍存在？不同自旋劈裂对称性（d/g/i 波）是否导致不同的混沌特征？

适合体系：MnTe（d 波）、CrSb（g 波）、Mn₅Si₃（d 波）、Fe₂WTe₄（d 波半导体）。

输入：每种材料的 DFT 交换参数 + VAMPIRE 模拟 + 混沌分析。

 

【进一步计算 3】自旋波激发与传播

为什么值得算：SOT 不仅驱动磁矩翻转，还激发自旋波——自旋波的传播和干涉可能影响翻转的可靠性和写入速度。

能回答的问题：SOT 激发的自旋波模式是什么？自旋波如何在交错磁体中传播？手性自旋波是否具有方向选择性？

适合体系：所有磁有序交错磁体，优先选择自旋波谱已知的体系。

输入：DFT 交换参数 + 自旋波色散计算 + 微磁学/原子自旋模拟中的自旋波激发。

 

【进一步计算 4】拓扑磁结构：skyrmion 与反 skyrmion

为什么值得算：交错磁中的各向异性交换相互作用可能稳定非平庸的拓扑磁结构——skyrmion 是未来自旋电子学器件的信息载体。

能回答的问题：交错磁中是否存在 skyrmion？skyrmion 的尺寸和稳定性？SOT 能否驱动 skyrmion 运动？

适合体系：具有 DMI 的交错磁异质结（如 AM/Pt 或应变调控的 AM）。

输入：DFT 交换参数 + DMI 参数 + 微磁学模拟（MuMax3 或 OOMMF）。

 

【进一步计算 5】器件级别模拟：SOT-MRAM 写入/读取

为什么值得算：将原子自旋模拟的结果外推到器件级别——评估交错磁 SOT-MRAM 的实际性能指标。

能回答的问题：写入电流密度和功耗？写入速度和可靠性？读取信号的信噪比？器件 endurance？

适合体系：RuO₂/Pt SOT 器件、Mn₅Si₃/Pt 器件——已有实验验证的体系。

输入：原子自旋模拟参数 + 微磁学模拟 + 电输运模型（自旋扩散方程）。

 

【进一步计算 6】机器学习力场：加速 DFT 精度的原子自旋模拟

为什么值得算：当前原子自旋模拟使用经典 Heisenberg 模型——但 DFT 精度的交换相互作用是磁矩方向的函数，经典模型无法捕捉。

能回答的问题：机器学习力场能否以 DFT 精度描述非共线磁矩排列的能量？能否在 ML 力场中捕捉到混沌动力学？

适合体系：RuO₂ 或 MnTe——简单二元化合物，适合作为概念验证。

输入：DFT 计算的非共线磁构型数据库（~1000 种构型）→ 机器学习训练 → VAMPIRE 集成。

 

曹彦彰, 韩磊, 潘峰, 王钰言, 宋成 | 物理学报 75, 060709 (2026) | DOI: 10.7498/aps.75.20251628 | 交错磁体 · SOT · 磁动力学 · 混沌 · 存算融合