\textit{\textbf{First-principles}} description of pumped inelastic X-ray scattering: example of K-edge RIXS in graphite

本文提出了一种结合 Bethe-Salpeter 方程与实时时间依赖密度泛函理论的从头算框架，用于预测光泵浦材料中的时间分辨共振非弹性 X 射线散射（RIXS），并通过成功模拟不同延迟时间下石墨的角相关 K 边 RIXS 光谱，证明了其准确性。

要点

想象一下你正在试图理解一台复杂机器的工作原理，比如一架大钢琴。你可以通过听它演奏一首曲子来了解它（这就像普通的谱系学/光谱学），但**共振非弹性X射线散射（RIXS）**则像是用一把由光组成的锤子敲击特定的琴键，聆听它发出的声音，并分析内部琴弦和琴槌是如何响应并振动的。它不仅能告诉你这台机器是由什么组成的，还能告诉你它的部件是如何运动和相互作用的。

这篇论文介绍了一种全新的、超高精度的计算机程序，它可以精确预测这种“光锤”实验看起来是什么样的，即使是在材料被第二波更快的脉冲光（类似于“泵浦”）剧烈摇晃时。

以下是使用日常类比对他们工作的拆解：

1. 问题所在：预测不可预测之物

科学家们长期以来一直能够利用X射线对材料进行“快照”拍摄。然而，要精确预测这些快照看起来是什么样的——尤其是当材料正被激光“泵浦”以将其从睡眠中唤醒时——是非常困难的。

• 旧方法： 以前的计算机模型就像是在观察一群人，并假设每个人都在原地不动且各自行动。它们忽略了人们（电子）实际上是如何手拉手并协同运动的（这种现象被称为“激子效应”）。
• 新方法： 作者构建了一个新的框架，它就像一个高速的3D电影模拟器。它不仅观察个体，还观察整个群体是如何共同起舞的，并考虑了他们之间是如何互相牵引的。

2. 方法论：两步舞曲

研究人员结合了两种强大的工具来创建他们的模拟：

• 第一步（“泵浦”）： 他们使用了一种叫做 RT-TDDFT 的工具来模拟激光撞击材料时发生的情况。想象一下向蹦床照亮手电筒：这个工具可以计算在光线撞击后，蹦床如何弹跳以及上面的人们如何移动重心。这为他们提供了在激光脉冲之后电子所处的“非平衡态”地图。
• 第二步（“探测”）： 然后他们使用了 Bethe-Salpeter 方程 (BSE)。你可以把它看作是一本关于X射线如何与那个跳动的蹦床相互作用的超精确规则手册。它计算了被踢出的电子与它留下的“空穴”（空位）之间复杂的舞蹈。

通过结合这两者，他们可以预测对于任何入射光角度和任何出射光角度的“回声”（散射X射线）。

3. 测试案例：石墨（铅笔芯）

为了证明其方法的有效性，他们用石墨（铅笔芯的成分）进行了测试。

• 为什么选择石墨？ 它就像一叠纸张。每一层纸内部的原子被紧紧粘在一起（就像强力胶），但层与层之间却只松散地堆叠在一起（就像一叠松散的纸）。这使得这种材料具有很强的“各向异性”，意味着取决于你是从侧面还是从顶部观察，它的表现会非常不同。
• 结果： 计算机模拟成功预测了石墨会演奏出的两种截然不同的“音符”：
  • $\pi$ (Pi) 音符： 这些来自于电子在层与层之间移动（松散的纸）。
  • $\sigma$ (Sigma) 音符： 这些来自于电子在层内部紧密运动（强力胶）。
    模拟显示，如果你从侧面照射，你会听到主要是“胶水”的音符；如果你从顶部照射，你会听到“纸张”的音符。这与现实世界的实验完美契合。

4. “泵浦”实验：摇晃桌子

这篇论文最令人兴奋的部分是，在用X射线撞击石墨之前，先用激光对其进行“泵浦”。

• 类比： 想象石墨是一个平静的池塘。激光泵浦就像向池塘里扔了一块石头，产生了涟漪。X射线则是发出的声呐脉冲，用来观察涟漪如何改变了水面。
• 发现： 当石墨被“泵浦”时，模拟显示“音符”发生了细微的变化。在低能区出现了新的、微弱的声音，且现有声音的音量发生了偏移。
• 结论： 计算机预测，即使是极短的激光脉冲也会改变材料的电子“情绪”，创造出一种不同于其静止状态的临时状态。模拟结果与实验数据高度吻合，甚至能观察到这些细微的偏移，证明了该方法适用于“时间分辨”（逐帧电影式）的研究。

总结

简单来说，这篇论文是在说：“我们构建了一个全新的、高精度的计算机模型，它可以精确预测材料在受到激光剧烈震动时，对X射线的反应。”

他们用石墨进行了测试，计算机的“预测”与现实中的“实验”完美匹配，正确识别了材料的内部结构（紧密的层与松散的层）如何响应不同角度和不同时间的照射。这为科学家提供了一个强大的新工具，让他们无需为每一次猜测都进行昂贵的实验，就能理解材料在实时状态下的行为。

技术摘要

技术摘要：泵浦非弹性 X 射线散射的第一性原理描述

问题与背景
共振非弹性 X 射线散射（RIXS）是一种强大的光谱技术，用于探测具有元素、轨道和动量分辨率的电子激发。虽然 RIXS 的理论框架已从独立粒子近似发展到利用 Bethe-Salpeter 方程（BSE）的多体微扰理论（MBPT），以捕捉激子效应，但在描述时间分辨（泵浦）RIXS 方面仍存在显著差距。现有的理论模型通常缺乏同时兼顾以下三点的能力：

1. 核心态与价带激发中的激子效应，这对于准确描述光谱至关重要。
2. 由光学泵浦诱导的非平衡电荷载流子分布。
3. 入射与出射光子的任意极化几何结构，这对于各向异性材料至关重要。

时间分辨 RIXS（tr-RIXS）实验正在兴起，但仍处于起步阶段；理论建模需要一个能够整合电子-空穴相关性与非平衡动力学的框架。

方法论
作者提出了一种结合**实时时间依赖密度泛函理论（RT-TDDFT）**与 **Bethe-Salpeter 方程（BSE）**的从头算框架，用于预测极化和时间分辨的 RIXS 光谱。该方法流程如下：

1. 理论形式： 该方法基于用于计算双微分截面（DDCS）的 Kramers-Heisenberg 公式。RIXS 过程被描述为一个两步机制：

   • 吸收： 入射 X 射线光子（$\omega_1$）产生一个内壳层空穴中间态（$|n\rangle$）。
   • 发射： 价电子填充内壳层空穴，发射一个光子（$\omega_2$）并留下一个最终的电子-空穴对（$|f\rangle$）。
     DDCS 使用由 BSE 本征向量和本征值导出的跃迁算符进行计算。

2. 集成非平衡态： 为了模拟泵浦情况，该框架采用 RT-TDDFT（具体在 exciting 代码中）在光学泵浦激光下传播 Kohn-Sham 系统。这产生了随时间变化的、非平衡的电子占据数。

3. 工作流：

   • 基态： 通过 DFT 计算提供起点。
   • BSE 计算： 进行两次独立的 BSE 计算：
     • 一次针对内壳层激发（以获得内壳层振子强度 $t^{(1)}$）。
     • 一次针对价带/光学激发（以获得最终态路径 $t^{(2)}$）。
   • 非平衡扩展： 对于 tr-RIXS，使用来自 RT-TDDFT 的非平衡占据数作为 BSE 计算的约束条件。
   • 极化泛化： 该实现被扩展以处理入射和出射光的任意、独立的极化矢量（$\mathbf{e}_1$ 和 $\mathbf{e}_2$），超越了以往 $\mathbf{e}_1 = \mathbf{e}_2$ 的限制。

4. 计算细节： 计算使用全电子全势能软件包 exciting 进行。研究采用了带有范德华修正的 GGA-PBE 泛函来处理石墨的碳 K 吸收边。BSE 包括用于光学光谱的 8 个占据带和 15 个非占据带，以及用于 X 射线光谱的 40 个导带。

主要贡献

• 统一框架： 本文建立了一个完全无参数的从头算工作流，结合了 RT-TDDFT 和 BSE 来模拟 tr-RIXS，捕捉了多体激子效应和泵浦诱导的非平衡动力学。
• 极化泛化： 作者将 RIXS 实现进行了泛化，以支持入射和散射光子的任意极化方向，从而能够研究复杂的散射几何结构和混合极化。
• 在 exciting 中的实现： 该方法实现在于全势能代码 exciting 中，并利用 pyBRIXS 程序包进行截面计算。

结果：石墨 K 吸收边 RIXS
该框架被应用于石墨，这是一种在面内（$\sigma$）和面外（$\pi$）状态之间具有强电子各向异性的半金属。

• 平衡态光谱：

  • 计算的吸收光谱（包括光学和内壳层吸收）重现了实验特征，证实了利用 BSE 捕捉激子效应（特别是 $\pi \to \sigma^*$ 和 $\sigma \to \pi^*$ 跃迁）的必要性。
  • 各向异性： 平衡态 RIXS 光谱显示出强烈的极化依赖性。面内极化（$\mathbf{e} \parallel x$）突出了高能 $\sigma$ 衍生激发（$>8$ eV），而面外极化（$\mathbf{e} \parallel z$）则强调了低能 $\pi$ 衍生激发。
  • 散射几何： 对匹配实验设置的 $30^\circ$ 入射角的模拟揭示了 $\pi$ 和 $\sigma$ 贡献的非线性混合，显示出干涉效应并增强了低能量损失处的 $\pi \to \pi^*$ 跃迁。

• 泵浦（非平衡）光谱：

  • 模拟针对由 266 nm 激光（持续 45 fs）泵浦且延迟时间为 150 fs 的系统进行。
  • 光谱变化： 泵浦光谱与平衡态情况非常相似，表明发生了适度的修改。然而，对于两种极化，都出现了一个明显的额外的低能量损失特征（几个 eV），这归因于费米能级附近电子分布的改变。
  • 混合极化： 在混合极化设置中，非平衡光谱与平衡态相比表现出光谱权重的重新分布和特征展宽。
  • 与实验对比： 理论结果很好地重现了实验峰位和一般特征。微小的偏差（例如峰宽和特定位移）归因于未包含在当前静态框架中的电子-声子耦合以及实验中潜在的自吸收效应。

意义与主张
作者声称，其工作为计算复杂材料（从金属到大带隙半导体）的 RIXS 光谱提供了一条通用且完全无参数的途径，适用于任意散射几何结构。通过成功重现石墨的特征各向异性及其在光学泵浦下的演化，该方法证明了其具备以下能力：

1. 解析截然不同的 $\pi$ 和 $\sigma$ 衍生特征。
2. 捕捉电子结构与极化之间的相互作用。
3. 模拟由光学激发引起的随时间变化的现象。

论文结论指出，这种方法为解释 tr-RIXS 实验和理解固体中的非平衡现象提供了强大的工具，特别强调了基本激发的方向性依赖。作者对实验中的差异保持谦逊，将其归因于已知局限性，如排除了电子-声子耦合和自吸收效应。