🔬 材料科学

Coherent and Incoherent Time-Resolved Diffuse Scattering Intensity from First Principles

本文提出了一种基于密度矩阵形式理论的第一性原理公式，该公式将时间分辨弥散散射强度中的相干与非相干贡献统一起来，成功重现了光激发锗中的实验观测结果，从而能够对超快晶格动力学进行系统性的解释。

原作者： Matteo Furci, Giovanni Marini, Stefano Mocatti, Matteo Calandra

发布于 2026-10-16

📖 1 分钟阅读☕ 轻松阅读

原作者： Matteo Furci, Giovanni Marini, Stefano Mocatti, Matteo Calandra

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

为了理解固体对象的行为——从其导热方式到其相变过程——科学家必须深入内部观察其原子的运动。在晶体中，原子并非冻结在原地；它们不断地振动，像连接着弹簧的小质量块一样，围绕其平衡位置上下跳动。这些被称为声子（phonons）的振动携带能量并定义了材料的特性。几十年来，研究人员一直使用 X 射线和电子束来研究这些振动。当一束光线击中晶体时，大部分光线会以可预测的模式反射，但有一小部分微弱的光会以弥散、模糊的方式散射。这种弥散散射是原子之间如何相互运动的直接映射。通过测量这种微弱的光芒，科学家可以重建材料的振动谱。

挑战始终在于速度。在正常的稳定晶体中，这些振动是随机且混沌的，会产生一个稳定的背景信号。但当材料被超短脉冲光击中时，原子会被撞离平衡状态。它们开始以新的、复杂的方式运动，有时同步振动，有时则发生随机加热。捕捉这种瞬时的演化需要一种足够快到可以“冻结”运动的探测手段，以及一个足够先进的理论来解释所见之物。直到现在，解释这些快速变化通常依赖于简化的模型，这些模型将有序、同步的原子运动与随机的加热运动分开处理。这种分离使得在激发后的最初几万亿分之一秒内，难以获得关于能量如何在固体中流动的完整图像。

由伦托大学（University of Trento）的一组研究人员开发了一种全新的、全面的方法，可以从第一性原理出发计算这些时间分辨信号。他们没有依赖简化的假设，而是构建了一个数学框架，将原子的同步运动和随机运动置于同等地位进行处理。他们首先考虑探测脉冲（例如来自自由电子激光器的 X 射线）如何与刚刚被激光脉冲击中的晶体发生相互作用。其核心洞察是：假设探测脉冲足够长，以至于可以平均掉快速的电子变化，但又足够短，足以捕捉到较重的原子核的较慢运动。利用这种方法，他们推导出了追踪原子位置随时间如何相互关联的方程。

研究人员发现，实验中观察到的信号实际上是两种不同行为的混合物。第一种是不相干背景，它在几皮秒内缓慢增长。这代表了晶格随着能量从受激电子流向振动的原子而升温的过程，导致随机的跳动增加。第二种是相干振荡，即信号的一种有节奏的脉动，其频率是原子振动频率的两倍。这种有节奏的部分产生于激光脉冲突然改变了原子键的刚度，导致原子在一种同步的、被挤压的状态下进行振荡。新理论表明，这种同步运动是由两方面驱动的：一是激光引起的振动频率的突然偏移，二是系统处于非平衡态时电子与振动原子之间的散射。

为了测试他们的理论，团队将其应用于一种常见的半导体——锗。他们使用了一个名为 EPIq 的计算机代码进行模拟，将有关该材料电子结构及其原子如何相互作用的数据输入其中。模拟设置旨在模仿一个真实的实验，即用一个 800 纳米的激光脉冲轰击锗晶体，并用 X 射线进行探测。研究人员计算了动量空间中五个不同点的弥散散射强度，追踪了信号在十皮秒窗口内的变化。他们发现，他们的计算能够以极高的准确度重现实验数据。模拟正确预测了信号发生振荡的具体频率，匹配了实验室中观察到的双频节奏。它还捕捉到了不相干背景的稳步上升，这对应于晶格的逐渐加热。

研究表明，随着动量转移的增加，振荡幅度会减小，这一细节符合实验观察结果，并反映了原子在不同方向上运动的特定方式。此外，研究人员能够观察到，随着受激电子密度的增加，原子的振动频率如何发生软化。这种软化是因为激光脉冲将电子从成键态提升到反键态，从而有效地削弱了将原子连接在一起的化学键。团队指出，虽然他们的模型在处理原子振动基本保持独立的锗时表现良好，但若要处理振动相互混合或晶体结构发生对称性变化的材料，则需要对其进行扩展。

目前这项工作的一个局限性在于，研究人员必须假设电子和原子如何随时间升温，而不是从头开始计算这个加热过程。他们将温度升高建模为一个平滑的指数过程，这足以匹配实验数据，但并不是一个完全微观的描述。作者建议，未来版本的代码将整合电子和原子的完整动力学，从而消除这些假设的需求，并提供一个更加完整的超快材料行为图景。目前，这一新框架为解释超快衍射实验中观察到的复杂信号提供了一个强大的工具，使科学家能够以以往无法达到的清晰度，区分原子的有序运动与混沌运动。

技术摘要：基于第一性原理的相干与非相干时间分辨弥散散射强度

问题陈述
来自 X 射线自由电子激光器（XFEL）和超快电子衍射（UED）的时间分辨弥散散射为研究光激发固体中的非平衡晶格动力学提供了途径。实验信号表现出两种截然不同的时间调制：一是由于晶格加热而在皮秒量级内上升的缓慢非相干背景；二是频率为两倍声子频率的快速相干振荡分量。现有的理论解释通常依赖于简化的模型，将这两个分量分开处理。非相干背景通常由多温度模型（例如 Allen 的双温度模型）来描述，而相干振荡则常通过假设瞬时频率偏移的“骤降淬火”（sudden-quench）近似进行建模。目前缺乏一个统一的第一性原理框架，能够在同等条件下推导出这两个贡献，特别是要考虑到光诱导频率偏移与非平衡载流子分布产生的电子-声子散射之间的相互作用。

方法论
作者利用密度矩阵形式法，提出了时间相关一阶弥散强度的第一性原理表述。该方法假设探测脉冲的相干时间长于电子时间尺度，但短于离子运动时间尺度。

理论框架：
- 散射截面根据原子位移进行展开。时间相关的一阶弥散强度被确定为与位移-位移相关函数 $\langle u_{q\nu} u_{-q\nu'} \rangle_t$ 成正比。
- 系统由电子-声子哈密顿量描述，其中声子简正模在与随时间变化的、处于非平衡态的 Born-Oppenheimer 势能相关的动力学矩阵瞬时对角化。
- 针对位移-位移相关函数导出运动方程（EOM）。通过四点相关函数的簇展开并结合马尔可夫近似来封闭该层级。
- 该形式区分了非相干贡献（受声子占据数 $n_{q\nu}$ 控制）和相干贡献（受反常声子相关函数 $\langle \hat{a}_{q\nu} \hat{a}_{-q\nu} \rangle$ 控制）。
关键方程：
- 反常相关函数的运动方程（方程 11）表明，相干动力学由两个来源驱动：
  1. 一个“挤压驱动”（squeezing drive）项，与 $\dot{\omega}_{q\nu}/\omega_{q\nu}$ 成正比，源于光诱导的声子频率偏移。
  2. 一个电子-声子源项，源于与非平衡载流子的散射，这可以在没有晶格不稳定性（instability）的情况下产生反常相关性。
- 总弥散强度（方程 12）表示为非相干项（ $2n_{q\nu} + 1$ ）与相干项（ $2\text{Re}\langle \hat{a}_{q\nu} \hat{a}_{-q\nu} \rangle$ ）之和。
实现：
- 该方案在 EPIq 代码中实现，利用 Wannier 插值高效评估电子-声子矩阵元的倒空间积分。
- 进行了两次独立的计算（基态和光激发态），以获得初始和最终的声子谱及变形势。
- 电子和声子占据数的时间演化采用唯象方法进行建模：电子占据数在 ~200 fs 内向准平衡态转变，而声子频率和占据数分别在可比的以及皮秒量级的尺度上演化，使用了平滑切换函数和随时间变化的晶格温度。

主要结果
该框架应用于光激发的锗（Ge），并将计算结果与飞秒弥散散射实验（文献 7）进行了对比。

重现实验特征： 计算得到的弥散强度成功重现了：
- 振荡频率，其对应于光激发声子支的 $2\omega_{q\nu}$ 。
- 振荡幅度的波矢依赖性。
- 上升的非相干背景，该背景随着能量从载流子转移到晶格导致声子占据数增加而自然产生。这一特征在标准的骤降淬火模型中是不存在的。
声子软化： 计算出的横向声学声子支随激发载流子密度的增加呈现单调软化，这与共价键的削弱一致。软化的幅度具有波矢依赖性，在靠近布里渊区边界处最大。
定量一致性： 通过假设激发载流子密度为 $n_{exc} = 0.0075$ 每原胞，且电子温度为 1000 K，该模型实现了与实验时间迹的良好吻合，包括相干振荡的衰减（通过引入唯象阻尼项来解释未被完全捕捉到的缺陷和非简谐性）。

意义与主张
本文声称提供了一种系统的、基于第一性原理的途径，用于解释时间分辨弥散散射，通过统一处理相干与非相干信号来实现。

统一性： 本文证明了相干与非相干贡献可以自然地从同一个密度矩阵形式中产生，既恢复了骤降淬火图景作为其极限情况，又将其扩展到了包含电子-声子散射效应的范畴。
机制识别： 文中明确指出，非平凡的相干动力学是由光诱导的声子频率偏移以及与非平衡载流子的电子-声离子散射共同诱导的。
局限性与未来工作： 作者谦逊地指出，目前的实现依赖于对载流子和声子群体的唯象处理（使用切换函数和共同的晶格温度），而非完全微观的多体传播。他们指出，将目前的方案与 EPIq 中已实现的基于第一性原理的载流子和声子动力学（碰撞积分）相结合，是消除这些唯象近似并改进背景增长及非简谐衰减描述的必要步骤。关于“无模式混合”（no-mode-mixing）近似的说明是：该近似对于 Ge 是有效的，但预计在经历对称性降低结构转变（如 VO $_2$ ）的系统中会失效，届需要使用更通用的模式混合方程。

技术摘要：基于第一性原理的相干与非相干时间分辨弥散散射强度

类似论文