Ultrafast dynamics of moments in bulk ferromagnets

基本信息

• 论文ID: 2502.07375
• 标题: Ultrafast dynamics of moments in bulk ferromagnets
• 作者: Mouad Fattouhi, Pascal Thibaudeau, Liliana D. Buda-Prejbeanu
• 分类: cond-mat.mes-hall cond-mat.mtrl-sci physics.app-ph
• 发表时间: May 22, 2025 (arXiv预印本)
• 论文链接: https://arxiv.org/abs/2502.07375

摘要

本文创建了一个稳健高效的模型来研究激光脉冲激发的磁性材料的超快动力学，该模型整合了动态Landau-Lifshitz-Bloch方程、量子热浴和双温度模型。该模型已成功应用于文献中的三种典型材料：镍、钴和铁。此外，对超快动态磁化率张量的分析表明，非对角元素根据是否存在连续外磁场而显示特定特征。

研究背景与动机

解决的核心问题

1. 精确控制热涨落：在小型磁性器件中，精确控制相关物理量的热涨落以实现可靠操作仍然是一个挑战
2. 超快磁化动力学建模：由于涉及不同长度和时间尺度的复杂过程，目前尚无单一的超快磁化动力学模型
3. 计算效率问题：现有的随机Landau-Lifshitz-Gilbert (sLLG)方程需要多次随机实现和平均化，计算效率低

研究重要性

1. 超快自旋电子学发展：飞秒激光激发方案能够快速产生自旋流，进入超快自旋电子学领域
2. THz自旋电子器件：为THz自旋电子器件的发展铺平道路，这是超快自旋物理学的首个实际应用
3. 理论理解的迫切需求：超快磁性器件的优化和可靠控制只有通过理论理解才能实现

现有方法局限性

1. 经典LLB方程：需要预先了解精确的热力学函数，包括交换刚度、自发磁化和阻尼参数的温度依赖性
2. 随机方法计算成本：sLLG方程需要大量随机实现才能获得可靠的平均值
3. 热连接问题：电子系统温度与原子自旋动力学的随机或确定性形式的耦合方式仍然是一个问题

核心贡献

1. 动态LLB (dLLB)模型：提出了一个基于统计平均程序的dLLB方程组，能够描述自旋变量的一阶和二阶矩动力学
2. 量子热浴集成：将量子涨落-耗散定理与dLLB模型结合，使用Bose-Einstein统计替代经典Maxwell-Boltzmann统计
3. 高效确定性方法：相比sLLG方程，dLLB方法的计算速度提高了100倍以上
4. 超快动态磁化率分析：首次系统分析了超快动态磁化率张量的对角和非对角元素行为

方法详解

任务定义

研究激光脉冲激发下铁磁材料的超快去磁和再磁化过程，包括：

• 输入：激光脉冲参数（功率密度、脉宽、波长）、材料参数
• 输出：磁化强度随时间演化、磁化率张量动力学
• 约束条件：涨落-耗散定理、能量守恒

模型架构

1. dLLB方程组

对于每个原子位点i，dLLB模型简化为以下方程组：

一阶矩动力学方程：

(1 + α²) ds_i/dt = ω_i × s_i - α(Σ_i - Tr(Σ_i)1)ω_i - 2D/(1 + α²) s_i

二阶矩动力学方程：

(1 + α²) dΣ_i/dt = 2ΩΣ_i - 2α([M^A_i, Σ_i] + Tr(M_i)(Σ_i - 2Γ_i) - (Tr(Σ_i) - 2Tr(Γ_i))M^S_i) - 2D/(1 + α²)(3Σ_i - 1Tr(Σ_i))

其中：

• s_i ≡ ⟨S_i⟩：平均磁矩方向
• Σ_i ≡ ⟨S_i ⊗ S_i⟩：二阶矩张量
• α：横向Gilbert阻尼参数
• D：噪声幅度

2. 量子热浴模型

经典涨落-耗散关系 (CFDR)：

D = (α/ℏ)k_B T

量子涨落-耗散关系 (QFDR)：

D = (α/ℏ) ∫₀^∞ [ℏω/(exp(ℏω/k_B T) - 1)] g_m(ω, T) dω

其中g_m(ω, T)是磁振子态密度。

3. 双温度模型 (2TM)

C_e dT_e/dt = -G_{e-ph}(T_{ph} - T_e) + P(t)
C_{ph} dT_{ph}/dt = G_{e-ph}(T_{ph} - T_e)

激光激发项：

P(t) = (μF)/(ℓ_p τ) exp(-(t-t₀)²/(τ²/4log(2)))

技术创新点

1. 统计平均方法：基于泛函微积分的严格统计平均程序，避免了预先提供热力学函数的需求
2. 量子修正：使用四次磁振子色散关系和Van Hove奇点参数的磁振子态密度
3. 自洽闭合：通过二阶矩动力学方程实现层次方程的自洽闭合
4. 确定性求解：避免了随机方法需要多次实现的问题

实验设置

材料参数

研究了三种典型3d铁磁金属：

模拟设置

• Ni和Co：5×5×5单胞，N=500个原子自旋
• Fe：7×7×7单胞，N=686个原子自旋
• 边界条件：周期性边界条件
• 数值方法：4阶Runge-Kutta方案

评价指标

1. 居里温度预测精度：与实验值的比较
2. 去磁化时间尺度：激光脉冲后磁化强度恢复的特征时间
3. 计算效率：相对于sLLG方法的速度提升

实验结果

主要结果

1. 居里温度预测

• CFDR方法：系统性高估居里温度
• QFDR方法：与实验数据excellent agreement
• 改进机制：量子热浴捕获了低温下的磁振子物理

2. 超快去磁化

对三种材料在不同激光功率密度下的测试：

• Ni: 22-50 J/m²范围内与实验数据良好一致
• Fe: 15-60 J/m²范围内准确重现去磁和再磁化过程
• Co: 22-50 J/m²范围内捕获主要动力学特征

3. 计算效率提升

• 速度提升：比sLLG方法快100倍以上
• 原因：确定性方法避免了多次随机实现的需求

动态磁化率分析

无外磁场情况

• 对角元素：χ_, χ_, χ_在顺磁相趋近于1/3（各向同性）
• 非对角元素：χ_, χ_, χ_保持为零

有外磁场情况 (B = 1.5T along z)

• 对角元素：行为与无场情况类似
• 非对角元素：出现峰值和弛豫行为，反映空间取向间的关联

脉宽依赖性

不同激光脉宽(50 fs, 200 fs, 800 fs)下：

• 对角元素：对脉宽变化不敏感
• 非对角元素：脉宽增加时峰值幅度减小，宽度增大

相关工作

超快磁学理论模型分类

1. 直接相干相互作用：激光脉冲与自旋自由度的直接相互作用
2. 局域自旋动力学：激光加热或激发触发的局域自旋动力学
3. 非局域角动量转移：如超扩散自旋输运

现有热自旋模型

• 随机方法：直接模拟原子自旋的随机形式
• 确定性方法：如热宏自旋模型和LLB方程
• 原子尺度方法：基于已知哈密顿算符的原子自旋动力学

本文优势

1. 预测性：相比经典LLB方程，不需要预先知识
2. 效率性：相比sLLG方法，计算速度显著提升
3. 准确性：量子修正提供更准确的低温和临界行为

结论与讨论

主要结论

1. dLLB-QFDR模型：成功结合量子热浴的dLLB方法能够准确预测铁磁材料的热效应
2. 居里温度预测：量子修正显著改善了居里温度和低温磁化曲线的预测
3. 超快动力学：模型准确捕获了3d铁磁元素的超快去磁化主要特征
4. 计算优势：确定性方法相比随机方法有显著的计算效率优势

局限性

1. 平均场近似：忽略了空间关联，可能影响临界行为的精确描述
2. 简化2TM：假设电子和自旋温度相同，可能不够准确
3. 材料范围：目前仅验证了简单3d铁磁金属

未来方向

1. 扩展应用：方法可扩展到铁磁体、反铁磁体和多层结构
2. 改进2TM：考虑更复杂的温度依赖性和非线性效应
3. 空间关联：包含更高阶涨落以改进临界行为描述

深度评价

优点

1. 理论创新：dLLB方程的统计平均推导严格且系统
2. 量子修正：QFDR的引入显著改善了模型的物理准确性
3. 计算效率：确定性方法带来的效率提升具有重要实用价值
4. 系统验证：在多种材料和条件下进行了充分的验证

不足

1. 理论近似：平均场近似可能限制了模型在临界区域的精度
2. 参数依赖：模型需要多个材料参数，部分参数的实验确定较困难
3. 适用范围：目前主要针对简单铁磁金属，对复杂磁性材料的适用性有待验证

影响力

1. 学术贡献：为超快磁学提供了新的理论工具和计算方法
2. 实用价值：高效的计算方法有利于器件设计和优化
3. 可复现性：方法和参数描述详细，便于其他研究者复现和应用

适用场景

1. 超快自旋电子学器件设计：THz频段磁性器件的理论设计
2. 激光诱导磁化动力学研究：飞秒激光与磁性材料相互作用的基础研究
3. 磁性材料热效应分析：有限温度下磁性材料性质的理论预测

参考文献

论文引用了73篇相关文献，涵盖了超快磁学、自旋动力学、统计物理等多个领域的重要工作，为该研究提供了坚实的理论基础和实验对比数据。