Longitudinal magnon transport properties in the easy-axis XXZ Heisenberg ferromagnet on the face-centered cubic lattice

基本信息

• 论文ID: 2507.04082
• 标题: Longitudinal magnon transport properties in the easy-axis XXZ Heisenberg ferromagnet on the face-centered cubic lattice
• 作者: M. R. Parymuda (Institute for Condensed Matter Physics, National Academy of Sciences of Ukraine)
• 分类: cond-mat.str-el (凝聚态物理-强关联电子系统)
• 发表时间: 2025年10月24日 (arXiv v5版本)
• 论文链接: https://arxiv.org/abs/2507.04082

摘要

本文详细研究了面心立方(FCC)晶格上S=1/2铁磁XXZ模型的纵向磁热输运性质，该模型具有易轴交换各向异性(Δ>1)并由四个子晶格组成。研究采用线性自旋波理论(LSWT)和Kubo形式体系来计算纵向自旋电导率和热导率，两者在低温区域均表现出激活型温度依赖关系，并确定了它们的磁场依赖性。分析表明，磁振子能隙对于确保这些电导率的收敛性至关重要。此外，通过考察热导率与自旋电导率的比值，识别出低温下磁振子输运的Wiedemann-Franz定律类比。最后，证明这些结果可以推广到任意自旋的系统。

研究背景与动机

1. 研究问题

本文研究面心立方晶格上易轴各向异性铁磁体中磁振子(magnon)介导的自旋和热输运性质，特别关注：

• 低温区域的纵向自旋电导率σ和热导率κ的温度依赖关系
• 磁场对输运性质的调控机制
• 磁振子输运中Wiedemann-Franz定律的类比

2. 问题重要性

• 基础科学意义：自旋波输运是凝聚态物理的活跃研究领域，对理解磁有序系统中自旋激发的传播机制具有重要价值
• 技术应用前景：磁性绝缘体在磁振子自旋电子学中具有显著优势，因为缺少自由电荷载流子可减少能量耗散，使自旋波能够长距离传播且Gilbert阻尼参数较小
• 材料相关性：稀土硫族化合物如EuO和EuS是具有低Gilbert阻尼的FCC铁磁体，自旋-轨道耦合可产生交换各向异性，使XXZ模型成为合适的理论描述

3. 现有方法局限性

• 对各向同性Heisenberg模型(Δ=1)，色散关系中存在无能隙的声学模式，导致输运系数发散
• 二维系统中即使存在能隙和小的Gilbert阻尼，纵向输运系数仍会因Bose函数在零频率的红外奇异性而发散
• 缺乏对三维FCC晶格上各向异性铁磁体输运性质的系统研究

4. 研究动机

• 探索各向异性参数Δ、晶格几何结构和外磁场如何共同影响系统的输运性质
• 建立磁振子输运的Wiedemann-Franz定律，揭示自旋电导率和热导率之间的普适关系
• 为实际磁性材料的输运性质提供理论预测框架

核心贡献

1. 系统计算了FCC晶格上易轴XXZ铁磁体的输运系数：首次详细研究了包含四个子晶格的三维FCC结构中的纵向自旋电导率和热导率
2. 揭示了磁振子能隙的关键作用：证明了有限能隙(Δ>1)对于确保输运系数收敛的必要性，并给出了激活型温度依赖关系：
   • σ(T) ∝ T^(3/2)exp(-Δ_gp/k_BT)
   • κ(T) ∝ T^(5/2)exp(-Δ_gp/k_BT)
3. 建立了磁振子输运的Wiedemann-Franz定律：导出了普适的Lorentz数 L = 5/2(k_B/gμ_B)²，该定律独立于各向异性参数Δ和外磁场H
4. 分析了磁场对输运性质的调控：展示了外磁场如何通过调节能隙和磁振子密度来抑制输运系数
5. 提供了推广到任意自旋S的理论框架：证明了结果可通过适当修改Holstein-Primakoff变换推广到任意自旋系统

方法详解

任务定义

计算三维FCC晶格上S=1/2易轴XXZ铁磁体的纵向自旋电导率σ(T,H)和热导率κ(T,H)，其中：

• 输入：晶格结构、交换相互作用参数(J₁, J₂)、各向异性参数Δ、外磁场H、温度T
• 输出：自旋电导率σ、热导率κ及其温度和磁场依赖关系
• 约束：低温区域(k_BT ≪ Δ_gp)，小Gilbert阻尼参数α ≪ 1

模型架构

1. 晶格结构与哈密顿量

FCC晶格包含四个子晶格(α=1,2,3,4)，晶格位置定义为R_mα = R_m + r_α，其中R_m = m_xe_x + m_ye_y + m_ze_z。

XXZ交换哈密顿量：

H = J₁∑_{⟨mα;nβ⟩} S_mα·S_nβ + J₂∑_{⟨⟨mα;nβ⟩⟩} S_mα·S_nβ

其中J₁, J₂ < 0表示铁磁交换，各向异性标量积定义为：

S_mα·S_nβ = (1/2)(S⁺_mαS⁻_nβ + S⁻_mαS⁺_nβ) + ΔS^z_mαS^z_nβ

外磁场项：-gμ_BH∑_m∑_α S^z_mα

2. Holstein-Primakoff变换

将自旋算符映射到玻色子算符(领头阶近似)：

S⁺_mα = √(2S)a_mα
S⁻_mα = √(2S)a†_mα  
S^z_mα = S - a†_mαa_mα

满足对易关系：a_mα, a†_nβ = δ_mnδ_αβ

3. 动量空间对角化

引入Fourier变换：a_mα = (1/√N)∑_q e^(iq·R_m)a_qα

哈密顿量在动量空间写为4×4矩阵形式：

H = H₀ + ∑_q (a†_q1, a†_q2, a†_q3, a†_q4) M (a_q1, a_q2, a_q3, a_q4)^T

通过对角化矩阵M得到四支色散关系。

4. 磁振子色散关系

长波极限(|q|→0)下的色散关系：

• 两支平坦带：ω_q1 = ω_q2 = 2|J₁|Δ + gμ_BH
• 两支色散带：ω_q± = Δ^±_gp + D±|q|²

其中能隙：

Δ^±_gp = |J₂|Δ± + gμ_BH
Δ± = [Δ(3+8δ)-3]/2 ± (1/2)√C₁

刚度常数D±在附录A中给出。

技术创新点

1. 能隙的关键作用

创新：明确证明了磁振子能隙对输运系数收敛的必要性

• 对Δ=1(各向同性Heisenberg模型)，存在无能隙声学模式导致σ,κ→∞
• 对Δ>1，各向异性自动打开能隙，即使H=0也能确保输运系数有限

物理机制：

• 能隙抑制了低能磁振子的热激发
• Bose分布函数n_B ≃ exp(-ω_qα/k_BT)提供指数截断
• 磁振子密度N_mag ∝ T^(3/2)exp(-Δ_gp/k_BT)

2. Kubo形式体系计算

方法优势：

• 通过线性响应理论严格计算输运系数
• 引入Gilbert阻尼参数α确保Green函数的收敛性
• 利用解析延拓从Matsubara频率到实频率

输运系数表达式：

σ±(T) = L±₁₁
κ±(T) = (1/T)[L±₂₂ - (L±₁₂)²/L±₁₁]

其中L±_αβ是通过磁化流和热流关联函数定义的输运系数。

3. 低温渐近行为

关键技术：利用积分渐近展开

∫^∞_Δ dx (dn_B/dx)[(x-Δ)^(μ/2)/x] ≃ -Γ(μ/2+1)(k_BT)^(μ/2)e^(-Δ/k_BT)/Δ

当Δ/k_BT ≫ 1时成立。

得到：

σ(T) ∝ T^(3/2)e^(-Δ_gp/k_BT)
κ(T) ∝ T^(5/2)e^(-Δ_gp/k_BT)

4. Wiedemann-Franz定律

推导：

L = κ/(σT) = (1/T²)[L₂₂/L₁₁ - (L₁₂/L₁₁)²]
  = [Γ(9/2)/Γ(5/2) - (Γ(7/2)/Γ(5/2))²](k_B/gμ_B)²
  = (5/2)(k_B/gμ_B)²

普适性：Lorentz数L独立于Δ、H、温度T和材料参数，仅依赖于g因子。

实验设置

参数选择

本文为理论计算研究，使用的模型参数为：

• 交换相互作用：|J₁| = |J|(能量标度)，|J₂| = 0.8|J|，即δ = |J₁|/|J₂| = 1.25
• 各向异性参数：Δ = 1.0001, 1.00012, 1.00014, 1.00016, 1.00018(略大于1以打开小能隙)
• Gilbert阻尼：α = 0.01(小阻尼近似)
• g因子：g = 2(对应S=1/2)
• 温度范围：T ∈ 0, 2.5×10⁻⁴|J|/k_B(低温区域)
• 磁场范围：H ∈ 0, 7×10⁻⁴|J|/gμ_B

计算方法

1. 色散关系计算：对角化4×4动力学矩阵M得到四支色散关系
2. Brillouin区积分：将动量求和转换为积分(1/N)∑_q → (1/(2π)³)∫d³q
3. 频率积分：使用delta函数近似处理Gilbert阻尼引入的展宽
4. 数值积分：对温度和磁场相关的积分进行数值计算

理论框架有效性条件

• 低温条件：k_BT ≪ Δ_gp，确保磁振子密度很小
• 小阻尼近似：α ≪ 1，LSWT忽略磁振子-磁振子相互作用有效
• 长波近似：考虑|q|→0附近的低能激发

实验结果

主要结果

1. 温度依赖性(H=0)

自旋电导率σ(T)（图3(a)）：

• 在T=0时σ=0(无热激发磁振子)
• 低温区域呈激活型行为：σ ∝ T^(3/2)exp(-Δ_gp/k_BT)
• 随Δ增大(能隙增大)，σ显著减小
• 对Δ=1(Heisenberg模型)，σ发散

定量数据（从图3(a)读取，Δ=1.0001）：

• T=1.0×10⁻⁴: σ ≈ 0.004 (单位：1/(aμ²_B))
• T=2.0×10⁻⁴: σ ≈ 0.008
• 呈现近似线性增长(指数项主导)

热导率κ(T)（图3(b)）：

• 类似激活型行为：κ ∝ T^(5/2)exp(-Δ_gp/k_BT)
• 相比σ(T)有更高的幂次温度依赖
• 对Δ=1同样发散

定量数据（Δ=1.0001）：

• T=1.0×10⁻⁴: κ ≈ 0.5×10⁻⁷ (单位：1/(ak_B|J|))
• T=1.6×10⁻⁴: κ ≈ 1.0×10⁻⁷

2. 磁场依赖性(T=10⁻⁴固定)

自旋电导率σ(H)（图4(a)）：

• H=0时达到最大值
• 随H增加单调递减
• 高磁场下趋于零

定量数据（Δ=1.0001）：

• H=0: σ ≈ 0.9×10⁻⁴
• H=3×10⁻⁴: σ ≈ 0.45×10⁻⁴(约减半)
• H=7×10⁻⁴: σ ≈ 0.1×10⁻⁴

热导率κ(H)（图4(b)）：

• 与σ(H)表现出相同的定性行为
• 数值上相差因子(5/2)T

定量数据（Δ=1.0001）：

• H=0: κ ≈ 2.0×10⁸
• H=7×10⁻⁴: κ ≈ 0.2×10⁸(减少约10倍)

3. 色散关系（图2）

• 高对称路径Γ-M-A-R-Γ上的四支色散关系
• 上下两支色散带在Γ点能量分别约为4|J|和2|J|
• 两支平坦带(不贡献输运)

关键发现

1. 能隙的双重作用

• 必要性：确保输运系数收敛(Δ>1)
• 权衡：最优输运需要Δ_gp ≪ 1(能隙尽可能小但非零)

2. 低能支主导

在低温下，仅下支色散带(ω_q-)贡献输运，因为：

• Δ⁺_gp ≫ Δ⁻_gp
• 上支的Bose占据数被指数抑制

3. 平坦带不贡献

虽然平坦带有很高的态密度，但群速度v_g=0，因此对σ和κ无贡献。

4. Wiedemann-Franz定律

κ/(σT) = L = (5/2)(k_B/gμ_B)² ≈ 3.125×10⁻⁵ eV²/K² (对g=2)

该比值与Δ、H、T无关，具有普适性。

5. 磁场调控机制

外磁场通过两种方式抑制输运：

• 增大能隙：Δ_gp → Δ_gp + gμ_BH
• 减少磁振子密度：N_mag ∝ exp(-gμ_BH/k_BT)

与文献的比较

与Ref.25的一致性：

• 简单立方晶格上反铁磁和铁磁Heisenberg模型也得到L=5/2
• 证实了WF定律的普适性(与晶格结构无关)

与二维系统的差异：

• 二维系统即使有能隙也会因红外发散导致输运系数发散
• 三维系统的相空间积分∫q²dq收敛

相关工作

1. 磁振子自旋电子学

• YIG材料6-8：最著名的低阻尼磁性绝缘体，室温下自旋波可长距离传播
• 稀土化合物9,10：EuO和EuS是FCC铁磁体，与本文模型直接相关
• 其他材料11-13：MnFe₂O₄, Cr₂O₃, Cu₂OSeO₃等展示了丰富的磁振子输运现象

2. 理论方法

• 线性自旋波理论14-19：Holstein-Primakoff变换是研究低温磁性的标准工具
• Kubo形式体系15,16：线性响应理论计算输运系数的严格框架
• Boltzmann方程25：另一种推导WF定律的方法，与Kubo形式结果一致

3. 输运性质研究

• 频率依赖18：FCC晶格上自旋电导率的频率响应
• 横向输运34-36：破缺反演对称性的系统中的反常输运
• 二维系统34,35：棋盘晶格上的拓扑磁振子输运

4. Wiedemann-Franz定律

• 电子系统：经典WF定律L=(π²/3)(k_B/e)²
• 磁振子系统24-27：
  • Ref.24首次提出磁振子WF定律
  • Ref.25用Boltzmann方程推导L=5/2
  • Ref.27研究强非线性区域的破坏

5. 本文的定位

• 补充现有研究：提供FCC晶格上纵向输运的系统分析
• 推广到各向异性：研究Δ>1情况，而非仅限Δ=1
• 强调能隙作用：明确指出能隙对收敛性的必要性
• 三维vs二维：避免二维系统的红外发散问题

结论与讨论

主要结论

1. 激活型输运行为：在易轴各向异性(Δ>1)FCC铁磁体中，自旋电导率和热导率在低温下表现为σ∝T^(3/2)exp(-Δ_gp/k_BT)和κ∝T^(5/2)exp(-Δ_gp/k_BT)
2. 能隙的关键作用：磁振子能隙是确保输运系数收敛的必要条件；对各向同性Heisenberg模型(Δ=1)，输运系数发散
3. 普适的Wiedemann-Franz定律：建立了磁振子输运的WF定律，Lorentz数L=(5/2)(k_B/gμ_B)²，独立于Δ、H和T
4. 磁场调控：外磁场通过增大能隙和减少磁振子密度来抑制输运，提供了磁热效应的微观机制
5. 推广性：理论框架可通过修改HP变换推广到任意自旋S的系统

局限性

1. 理论近似的限制

• 低温限制：LSWT仅在k_BT≪Δ_gp时有效，当温度升高磁振子-磁振子相互作用变得重要
• 小阻尼假设：α≪1的条件在高温下可能不满足，实际上α应具有温度依赖性α(T)∝T²37,38
• 领头阶近似：忽略了HP变换的高阶项，无法描述磁振子相互作用效应

2. 模型简化

• 各向异性参数：假设Δ₁=Δ₂=Δ(最近邻和次近邻相同)，实际材料可能更复杂
• 交换相互作用：仅考虑最近邻和次近邻，忽略了更远程相互作用
• 拓扑平凡：考虑的是保持反演对称性的系统，不涉及拓扑非平凡相

3. 实验验证缺失

• 论文为纯理论工作，缺少与实验数据的直接比较
• 需要在EuO、EuS等实际材料中验证预测

4. 适用范围

• 晶格结构：特定于FCC晶格，其他晶格结构需要重新计算
• 磁有序类型：仅适用于铁磁体，反铁磁体需要不同处理
• 维度效应：三维结果不能直接应用于二维系统

未来方向

1. 理论扩展

• 包含磁振子相互作用：考虑HP变换的高阶项，研究有限温度下的输运性质
• 温度依赖的阻尼：引入α(T)以更准确描述高温行为
• 非线性输运：研究大温度梯度或大磁场梯度下的非线性效应

2. 材料应用

• 实际材料参数：使用EuO、EuS等材料的实验参数进行定量预测
• 与实验对比：与自旋Seebeck效应、磁热效应实验比较
• 材料设计：优化Δ和J以实现最佳输运性质

3. 拓展到其他系统

• 反铁磁体：研究反铁磁XXZ模型的输运性质
• 其他晶格：简单立方、体心立方、蜂窝晶格等
• 拓扑系统：破缺反演对称性引入拓扑非平凡相

4. 多体效应

• 磁振子-声子耦合：研究晶格振动对磁振子输运的影响
• 磁振子-磁振子散射：超越LSWT的多体计算
• 量子涨落：低维系统中的量子效应

深度评价

优点

1. 方法的严格性

• 理论框架完整：从微观哈密顿量出发，通过LSWT和Kubo形式体系系统推导输运系数
• 数学推导详细：附录中给出了完整的解析表达式和积分计算过程
• 近似条件明确：清楚说明了低温、小阻尼等适用条件

2. 物理洞察深刻

• 能隙作用的阐明：明确指出能隙对收敛性的必要性，并解释了物理机制
• WF定律的普适性：证明了Lorentz数独立于多个参数，揭示了深层的物理联系
• 低能支主导：解释了为何只有下支色散带贡献输运

3. 结果的系统性

• 多参数研究：考察了Δ、T、H的影响
• 定性与定量结合：既给出解析渐近行为，又提供数值结果
• 图表清晰：色散关系、温度依赖、磁场依赖均有清晰展示

4. 学术贡献

• 填补研究空白：首次系统研究FCC晶格上易轴XXZ模型的纵向输运
• 与文献一致：WF定律与其他方法(Boltzmann方程)的结果一致，增强可信度
• 推广性强：框架可应用于其他晶格和任意自旋

不足

1. 理论局限

• 温度范围窄：仅适用于极低温(T~10⁻⁴|J|/k_B)，实际应用受限
• 忽略相互作用：LSWT的领头阶近似在磁振子密度较高时失效
• 阻尼处理简化：常数α不能描述温度依赖的耗散机制

2. 实验联系薄弱

• 缺少材料参数：未使用具体材料(如EuO)的实验数据
• 无实验验证：纯理论预测，缺少与实验的比较
• 测量可行性：未讨论如何实验测量预测的输运系数

3. 分析深度

• 平坦带讨论不足：仅简单说明不贡献输运，未深入探讨其物理意义
• 拓扑性质缺失：未考虑可能的拓扑相和边界态
• 有限尺寸效应：未讨论实际样品的尺寸限制

4. 技术细节

• 数值方法未详述：未说明积分的具体数值算法
• 误差分析缺失：未给出数值结果的精度评估
• 参数选择依据：为何选择|J₂|=0.8|J|未充分说明

影响力

1. 学术贡献

• 理论框架：为FCC磁性材料的输运性质提供了标准理论描述
• 普适规律：WF定律的建立具有基础科学价值
• 方法示范：LSWT+Kubo形式的组合可应用于其他系统

2. 实用价值

• 材料筛选：理论预测可指导寻找最优输运性质的材料
• 器件设计：磁场调控机制可用于自旋电子学器件
• 磁热应用：为磁制冷技术提供微观理解

3. 可复现性

• 高可复现性：
  • 理论推导步骤清晰
  • 参数选择明确
  • 附录提供完整公式
• 代码实现：基于给出的公式可独立编程验证
• 推广容易：框架可直接应用于其他参数和晶格

4. 后续研究

• 引用潜力：作为FCC磁性材料输运的标准参考
• 研究方向：激发对其他晶格、反铁磁体、拓扑系统的研究
• 实验推动：理论预测可促进相关实验测量

适用场景

1. 材料类型

• 稀土硫族化合物：EuO, EuS等FCC铁磁绝缘体
• 低阻尼材料：Gilbert阻尼α<0.1的系统
• 易轴各向异性：Δ>1的材料(自旋-轨道耦合产生)

2. 温度条件

• 极低温：T < 0.1Δ_gp/k_B，通常需要T < 1K
• 稀磁振子气体：磁振子密度≪晶格位点数

3. 实验技术

• 自旋Seebeck效应测量：温度梯度驱动的自旋流
• 磁热效应：磁场调控的熵变和温度变化
• 非局域输运：长距离自旋流传输

4. 理论研究

• 基准比较：作为更复杂理论(包含相互作用)的基准
• 参数拟合：从实验数据反推交换参数和各向异性
• 推广研究：扩展到其他晶格、自旋、维度

5. 不适用场景

• 高温区域：T~J/k_B时磁振子相互作用重要
• 二维系统：存在红外发散
• 强阻尼材料：α>0.1时LSWT失效
• 拓扑非平凡相：破缺反演对称性的系统

参考文献（关键文献）

14 T. Holstein and H. Primakoff, Phys. Rev. 58, 1098 (1940) - HP变换的原始文献

15 K. Aoyama and H. Kawamura, Phys. Rev. B 100, 144416 (2019) - 方形晶格反铁磁体的输运性质

16 K. Aoyama, Phys. Rev. B 106, 224407 (2022) - 三维反铁磁体的临界行为

24 K. Nakata et al., Phys. Rev. B 92, 134425 (2015) - 首次提出磁振子WF定律

25 K. Nakata et al., Phys. Rev. B 96, 224414 (2017) - 用Boltzmann方程推导WF定律

37 S. Tyc and B. I. Halperin, Phys. Rev. B 42, 2096 (1990) - 二维反铁磁体的自旋波阻尼

总体评价：这是一篇严谨的理论物理论文，系统研究了FCC晶格上易轴XXZ铁磁体的磁振子输运性质。主要贡献是建立了完整的理论框架，揭示了能隙的关键作用，并证明了WF定律的普适性。方法上结合了LSWT和Kubo形式体系，推导详细，结果可靠。局限性在于仅适用于极低温和小阻尼情况，且缺少与实验的直接比较。论文为理解磁性绝缘体的输运机制提供了重要理论基础，对磁振子自旋电子学领域有潜在影响。