Interplay of non-local transport and local scattering during electron thermalization and spatial equilibration in laser-excited metals

本文采用重新构建的玻尔兹曼输运方程来证明，虽然非局部输运通过移除非热载流子加速了辐照表面的表观热化过程，但它同时也延迟了整个电子系统的完全平衡，揭示了随位置和能量变化的输运与散射之间复杂的相互作用。

要点

想象一个位于金属块内部、拥挤不堪的舞池。突然，一道超快激光脉冲击中了表面，就像 DJ 投下了一个只有最前排舞者才能听到的重低音节拍。这些“舞者”（电子）变得兴奋起来，跳跃着并疯狂地移动，而房间后方的舞者仍然平静地坐着。

这篇论文研究的是接下来的过程：这些兴奋的舞者是如何冷静下来的，以及能量是如何在整个房间中传播的？

两种主要力量的作用

研究人员发现，有两种主要现象在同时发生，且它们往往相互制约：

1. “局部散射”（舞池碰撞）： 兴奋的舞者们互相碰撞，并撞向房间的墙壁（金属的原子结构）。这就像是一个混乱的冲撞区（mosh pit），最终每个人都会慢下来，开始跳起一种同步、冷静的节奏。这就是热化（thermalization）。
2. “非局部输运”（人群涌动）： 因为激光只击中了前方，所以前方的舞者既拥挤又充满活力，而后面则是空旷的。自然而然地，前方的这些充满能量的舞者开始向后方奔跑，以填补空位。这就是输运（transport）。

大惊喜：“假象”效应

这篇论文中最有趣的发现是一个视觉上的“诡计”。

如果你站在前门（表面）观察舞者，你可能会想：“哇，他们冷静下来的速度真快！”研究人员发现，输运实际上让表面看起来冷却得更快了。

为什么呢？因为兴奋的舞者正从前门跑向房间后方。他们并不一定是在前门处冷静下来，他们只是离开了前门。所以，如果你只观察表面，看起来系统似乎很快就达到了平和的平衡状态。

然而，论文指出，整个系统其实还没有真正平静下来。 舞者们仍在奔跑，试图填满房间的后部。前方的“平静”是一种由于人群移开而产生的错觉。只有当舞者们从前到后均匀分布后，整个系统才会真正变得平静。

“能量窗口”类比

研究人员还根据舞者有多“狂野”（能量水平）将其分成了特定的群体。

• “轻微兴奋”组（低能量）： 这些是仅仅有些坐立不安的舞者。他们的运动主要受**人群涌动（输运）**控制。他们大多只是从拥挤的前方移动到空旷的后方。
• “极度兴奋”组（高能量）： 这些是正在跳上桌子乱跳的舞者。他们的行为主要受**互相碰撞（散射）**控制。无论他们在房间的什么位置，他们都会通过迅速撞击他人来消耗掉狂野的能量。

核心结论

论文得出结论：你不能仅仅通过观察表面或假设所有事情都发生在同一个点上，来理解被激光击中的金属内部发生了什么。

• 在表面： 看起来冷却得很快，是因为“热”电子正在跑开。
• 在金属内部： 系统实际上仍然是混乱的，因为这些电子正在扩散，在向后方移动的过程中创造出一种新型的不平衡。

研究人员建立了一个新的数学模型（就像一个超级精确的舞池模拟器），能够同时追踪碰撞（散射）和奔跑（输运）。这有助于科学家理解，在厚金属中，“冷却”不仅仅是关于放慢脚步，还涉及到你站在房间里的什么位置。

技术摘要

技术摘要：激光激发金属中非局部输运与局部散射的相互作用

问题陈述
超快激光激发金属会诱导电子在能量分布上的局部非平衡态，以及在空间上的非局部非平衡态，特别是当样品厚度超过光学穿透深度时。随后的动力学过程受控于非局部输运与局部散射（电子-电子和电子-声子散射）之间的相互作用，然而，将这些过程结合起来的微观描述仍然非常稀缺。现有的理论方法通常依赖于简化假设：要么假设均匀加热以详细解析散射过程，要么保留输运过程但采用近似的散射模型（例如弛豫时间近似）。这些近似阻碍了人们区分测量到的载流子动力学中，哪些是由内在热化贡献的，哪些是由空间输运贡献的，特别是在能量分辨实验（如时间分辨光电子能谱，PES）中。

研究方法
作者提出了一个同时包含局部散射和非局部输运的空间分辨玻尔兹曼框架。该方法的核心是将玻尔兹曼输运方程在能量空间进行重构，并利用完整的微观碰撞积分。关键的方法论特征包括：

• 双分支描述： 为了在不进行全动量解析的情况下解析空间运动，作者将电子动量投影到表面法线方向（$z$ 方向）。能带被分为两个分支：向体相运动的电子（$v_z > 0$）和向表面运动的电子（$v_z < 0$）。
• 全碰撞积分 (FCI)： 模型采用了基于随机-$k$ 近似的电子-电子和电子-声子散射的全碰撞积分（FCI）。这避免了弛豫时间近似的不准确性，同时与全动量解析计算相比，保持了计算上的可行性。
• 散射与输运的耦合： 该框架明确考虑了可能改变运动方向（分支切换）的散射事件，并将输运视为这两个分支的空间演化。前表面和后表面的边界条件通过在分支间反射电子来处理。
• 模拟设置： 该模型应用于类铝自由电子金属（匹配了费米速度和吸收率的第一性原理参数），厚度为 50 nm。模拟涵盖了一系列激光通量范围（$0.05$至$15 \text{ J m}^{-2}$），并将非均匀 Beer-Lambert 激发与在相同吸收能量密度下的均匀激发进行了对比。

核心贡献
这项工作的首要贡献是开发了一个一致的微观框架，弥合了空间平衡化与能量空间散射之间的鸿沟。通过直接结合散射与电子运动，该模型隐含地覆盖了从弹道输运到扩散输运的整个范围。这使得直接量化输运过程如何修正表观热化动力学成为可能，由于以往研究依赖于简化模型，这种区分此前难以实现。

结果
研究得出了关于热电子动力学的几个关键见解：

1. 加速表面热化： 输运过程加速了在辐照表面观察到的“表观”热化。通过将非热载流子从表面区域移入体相，输运起到了一种额外的弛豫通道作用。这种效应在低通量下最为显著，因为此时内在的电子-电子散射较慢。
2. 延迟全局平衡化： 虽然表面热化看起来更快，但空间的能量重新分布延迟了整个电子系统的完全平衡。空间梯度的存在导致在整个样品深度内，非热载流子种群比均匀激发下持续存在更长时间。
3. 深度依赖的动力学： 驱动动力学的 dominant 机制随位置和能量的变化而显著不同：
   • 前表面： 动力学由电子-电子散射与输运（移除载流子）之间的竞争驱动。
   • 后表面： 动力学由输运过程主导，因为直接激光激发可以忽略不计。由于反射作用，非热电子会在后表面发生瞬态积累。
   • 能量依赖性： 在费米能级附近，输运是驱动谱密度变化的主要因素。在较高能量处（例如高于费米能级 1.0 eV），电子-电子散射主导了衰减，其发生的时间尺度与激光脉冲相当。
4. 谱密度演化： 对谱密度（模拟 PES 可观测物理量）的分析表明，对弛豫时间的解释高度依赖于所探测的能量窗口和深度。例如，在后表面，谱密度的变化几乎完全由输运而非局部散射驱动。

意义与主张
作者声称，其工作改进了对光学厚系统在激光驱动下的电子非平衡态的理解。他们强调，在这种系统中，电子热化不能被解释为一个纯粹的局部过程。一个关键发现是表面敏感的可观测物理量与全局电子状态之间的区别：表面测量可能会显示出快速的热化，因为输运过程将非热载流子从探测体积中移走，即便由于空间的重新分布，整个电子系统仍处于非平衡态。

论文得出结论，这种区别对于能量分辨的超快测量尤衷重要，因为主要的弛豫通道取决于探测的能量窗口和探测深度。所提出的框架提供了一条微观路径，用以连接输运、散射以及实验可及的热载流子动力学，为解释与未来涉及热载流子的电子应用相关的实验数据提供了更准确的基础。