Huge ultrafast spin Seebeck effect mediated by laser-excited superdiffusive magnon currents

本文引入了一种基于量子玻尔兹曼方程的从头算参数化微观框架，以揭示激光激发的体心立方铁薄膜中的超扩散输运机制和巨大的超快自旋塞贝克效应，从而克服了传统扩散模型在描述非热磁振子动力学方面的局限性。

要点

以下是用通俗语言和日常类比对该论文的解读。

宏观图景：激光引发的“自旋风暴”

想象一块铁（一种铁磁性金属）就像一个拥挤的舞池，每个人都手拉手朝同一个方向旋转。这种同步旋转就是我们所说的磁性。

1996 年，科学家发现，如果你用超快、超亮的激光脉冲击中这个舞池，舞者们几乎会瞬间停止同步旋转。磁性在不到一万亿分之一秒的时间内消失。这被称为超快退磁。

几十年来，科学家们一直在争论这究竟是如何发生的。旧理论就像试图通过假设每个人只是均匀地慢慢变热和变冷，来描述一个混乱的冲撞舞池（mosh pit）。但本文认为，现实要暴力且迅速得多：这是一种能量混乱的奔涌，其速度之快，简单的热扩散理论无法解释。

新理论：“量子玻尔兹曼”交通模型

作者（来自乌普萨拉大学）建立了一个更详细的计算机模型，以模拟激光击中铁时发生的情况。

1. 旧方法（三温模型）：
想象一个房间里有三组人：电子（快速移动者）、声子（振动的地板）和磁振子（旋转的舞者）。旧模型假设，当激光击中时，电子变热，它们将热量传递给地板，地板再将热量传递给舞者。最终，每个人都会达到“热平衡”（所有人的温度相同）。它将舞者视为只是慢慢变暖。

2. 新方法（非局域模型）：
作者指出，在前几皮秒（万亿分之一秒）内，这是错误的。激光产生的不是缓慢的升温，而是一股冲击波。

• 类比： 想象将一大桶水倒在干燥的海绵上。旧模型说水会慢慢渗入。新模型说水会以高速喷射出来，在来得及渗入之前四处飞溅。
• 机制： 激光激发电子，电子随后猛烈地“踢”向“旋转者”（磁振子）。这些旋转者不会只是坐在那里；它们像子弹一样从表面射出，将其自旋能量带入材料深处。

关键发现：“超扩散”奔涌

该论文确定了一种特定的机制，称为超扩散输运。用通俗英语解释如下：

• 弹道相（子弹）： 激光击中后，激发的磁振子立即沿直线运动，就像从枪里射出的子弹。它们尚未与任何物体碰撞。它们移动得极快（约每皮秒 35–50 纳米）。
• 扩散相（人群）： 几皮秒后，它们开始与其他粒子碰撞，减速并随机扩散，就像一群人在走廊里徘徊。
• “超”的部分： “子弹”阶段和“人群”阶段之间的过渡，就是作者所称的“超扩散”。它比正常的扩散更快、更高效。

“自旋塞贝克”效应：自旋海啸

该论文声称，这一过程会产生巨大的超快自旋塞贝克效应。

• 类比： 通常，自旋塞贝克效应就像一条缓慢的自旋河流，从高温区流向低温区。
• 论文的声称： 在这种超快情境下，它不是一条河流，而是一场海啸。激光在表面瞬间产生巨大的温差。这触发了一股自旋电流的“爆发”，其强度是正常稳态加热所产生强度的1000 倍。
• 为何重要： 这种爆发如此强劲且迅速，作者计算出，从理论上讲，它足以在短短 10 皮秒内翻转一层薄铁的磁化方向。这是制造超快计算机内存的“圣杯”。

将理论与现实连接：“克尔角”

我们如何知道这正在发生？我们无法直接看到磁振子。相反，科学家使用一种称为**磁光克尔效应（MOKE）**的工具。

• 类比： 想象用手电筒照射镜子。如果镜子具有磁性，反射回来的光会带有略微不同的偏转（偏振）。
• 论文的贡献： 作者利用他们的模型，精确预测了随着磁性在铁的不同深度消失和重现，这种光的“偏转”将如何随时间变化。他们发现，光信号以一种非常具体且反直觉的方式表现（有时即使磁性变弱，信号反而变强），这是因为“自旋海啸”正在向材料深处移动。

他们声称的总结

1. 旧模型太慢： 它们忽略了粒子初始的“子弹般”速度。
2. 新模型准确： 通过使用“量子玻尔兹曼”方程，他们可以追踪这些快速移动的粒子如何从表面射向深处。
3. 巨大的自旋电流： 激光产生了一股巨大且超快的自旋（磁振子）流，其强度远超稳态实验中看到的任何情况。
4. 两阶段退磁： 首先，表面瞬间失去磁性。然后，随着自旋电流的到达，一股“退磁波”向材料深处传播。
5. 实验证明： 他们计算了激光实验将“看到”什么（克尔信号），并表明这些信号包含了这种超快、深穿自旋电流的“指纹”。

简而言之： 该论文声称，当你用激光照射铁时，你不仅仅是在加热它；你正在发射一股高速、超强磁能波，其传播的深度和速度都超过了任何人此前认为的可能。

技术摘要

技术摘要：激光激发的超扩散磁振子电流介导的巨大超快自旋塞贝克效应

问题陈述
铁磁金属中超快激光诱导退磁的微观起源仍是一个争论的焦点，特别是关于横向与纵向自旋激发以及非局域输运机制的作用。虽然广泛使用的三温模型（3TM）成功复现了许多实验结果，但它假设子系统瞬时热化，无法捕捉角动量转移或亚皮秒时间尺度上存在的强非平衡动力学。之前的扩展，如$(N+2)$温模型（$(N+2)$TM），解决了非热磁振子布居的问题，但仅限于局域（空间均匀）动力学。因此，超快时间尺度上磁振子输运的理论描述主要局限于扩散模型，而这些模型难以捕捉电子自旋输运中观察到的弹道和超扩散机制。亟需一种框架，将非热磁振子散射与空间输运相结合，以解释任意厚度薄膜中深度分辨的退磁和自旋流。

方法论
作者提出了一个基于第一性原理参数化的微观框架，将局域$(N+2)$TM扩展为非局域版本，能够描述空间动力学。该方法的核心是利用量子玻尔兹曼方程（QBE）来模拟磁振子输运，超越了早期研究中使用的扩散近似。

• 理论框架：该模型耦合了电子、声子和磁振子的动力学。电子和声子通过热扩散方程处理（包含带有源/汇项的傅里叶定律），而磁振子则由源自 QBE 的漂移 - 扩散方程控制。该方程包含一个代表弹道传播的漂移项（$v_q \partial n_q / \partial z$）和一个源自$sp$-$d$哈密顿量和费米黄金定则的散射项（$s_q$）。
• 散射机制：散射项考虑了驱动磁振子趋向与电子热平衡的电子 - 磁振子散射事件（一阶），以及重新分配磁振子布居的磁振子数守恒散射事件（二阶）。
• 数值实现：该系统在体心立方（bcc）Fe 薄膜的空间离散域上进行求解。电子和声子方程采用隐式 Crank-Nicolson 格式求解，而磁振子平流方程则采用显式总变差减小（TVD）格式求解。
• 参数：研究利用了参考文献 [43] 中计算的第一性原理磁振子频率和散射率参数。模拟针对厚度为 100 nm 和 1 $\mu$m 的 Fe 薄膜进行，由 800 nm 激光脉冲（60 fs 脉宽）激发，注量范围为 0.01 至 4 mJ/cm$^2$。
• 可观测量：为了将理论与实验联系起来，作者计算了深度分辨的磁化分布，并使用广义传输矩阵法（pyGTM）模拟了时间分辨的磁光克尔效应（MOKE）信号。

主要贡献与结果

1. 超快自旋塞贝克效应：该模型预测了一种由强非均匀温度梯度驱动的“巨大”超快自旋塞贝克效应。激光激发后立即在表面附近产生非热磁振子的瞬态爆发。这些磁振子向体内传播，携带自旋流，其振幅达到约$5 \hbar \text{ nm}^{-2} \text{ ps}^{-1}$。该电流比线性响应理论预测的稳态温度梯度下的电流高出几个数量级。
2. 超扩散输运机制：模拟揭示了磁振子输运从初始的弹道机制向后期扩散机制的交叉。通过分析均方位移（MSD），作者识别出一个“超扩散”机制，其中输运指数$\alpha$从 2（弹道）过渡到 1（扩散）。该机制持续长达 10 ps，距离达 100 nm，这是纯扩散模型失效的时间尺度和长度尺度。
3. 深度分辨退磁动力学：非局域模型预测了取决于薄膜厚度和深度的不同退磁分布。
   • 表面：表现出快速退磁，随后是迅速的局部再磁化。
   • 体内：更深层的区域显示出延迟的退磁，因为热量和磁振子从表面传播而来。
   • 厚度依赖性：在较厚的薄膜（如 1 $\mu$m）中，系统需要数十皮秒才能达到空间均匀的磁化状态，而超薄薄膜（<20 nm）几乎瞬间均匀化，这验证了局域$(N+2)$TM 适用于薄膜。
4. MOKE 信号模拟：计算得到的深度依赖磁化分布被用于模拟 L-MOKE 信号。结果表明，在超快时间尺度上，克尔旋转和椭圆率并不线性地跟踪平均磁化强度。值得注意的是，从泵浦表面的对侧进行探测揭示了由弹道磁振子爆发到达引起的延迟的第二退磁阶段。这一阶段可能导致反直觉的行为，例如克尔分量幅度的增加或椭圆率符号的变化，而在没有输运的模拟中这些现象并不存在。

意义与主张
本文主张其框架提供了一条超越扩散模型描述超快非热磁振子输运的途径。通过结合量子玻尔兹曼方程，该模型成功捕捉了从弹道到扩散输运的过渡，这一特征对于理解厚度大于激光穿透深度的薄膜中的动力学至关重要。

作者断言，计算出的磁振子自旋流具有技术相关性，表明激光诱导的磁振子输运有可能在 10 ps 内诱导薄膜中磁化的显著重定向或反转。此外，MOKE 信号的模拟提供了一种直接方法，能够在无需复杂深度敏感重建的情况下弥合理论预测与实验可观测量之间的差距。对侧探测中出现的延迟退磁阶段被强调为磁振子介导输运的潜在实验特征，支持将超快退磁归因于横向自旋激发。该工作被呈现为一种工具，有助于时间分辨自旋输运实验的设计与解释。