🔬 materials science

Thermal diffuse scattering in TEM: complex absorptive potentials compared to the frozen phonon model

本文将广泛使用的复吸收势方法与更为复杂的冻结声子模型（基于相关原子运动和爱因斯坦模型）进行了对比，用于模拟透射电子显微镜中的热弥散散射和非弹性吸收效应。

原作者： Martin Hájek (Institute of Physical Engineering, Brno University of Technology), Ján Rusz (Department of Physics and Astronomy, Uppsala University)

发布于 2026-08-26

📖 1 分钟阅读☕ 轻松阅读

CC BY 4.0

原作者： Martin Hájek (Institute of Physical Engineering, Brno University of Technology), Ján Rusz (Department of Physics and Astronomy, Uppsala University)

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

当科学家们窥探原子世界时，他们经常使用一束穿过薄片材料的电子束，就像用手电筒照射彩色玻璃窗一样。在完美且静止的晶体中，这些电子会以可预测的模式从原子上反弹，从而创造出材料结构的清晰、锐利的图像。然而，原子从未真正静止；由于热量，它们不断地振动，这种现象被称为热运动。当电子与这些振动的原子碰撞时，它们会发生散射，从而使图像变得模糊，并产生一层朦胧的光背景。这种散射被称为热弥散散射，是研究人员必须予以考虑的主要噪声来源，以便观察材料的真实细节。为了在计算机上模拟这一过程，科学家们开发了不同的数学工具来预测电子束的行为。一种方法将能量损失视为一种简单的、稳定的损耗；而另一种方法则试图模仿原子原地抖动的、如同被冻结在时间中的混乱快照。了解哪种工具能提供最准确的图像，对于解读揭示新材料奥秘的高分辨率图像至关重要。

在最近的一项研究中，研究人员马丁·哈杰克（Martin Hájek）和扬·鲁斯（Ján Rusz）着手测试这些不同模拟工具的可靠性。他们专注于比较一种被称为复吸收势模型（complex absorptive potential model）的广泛使用的、计算效率高的模型，与一种被称为冻结声子模型（frozen phonon model）的更详细的方法。高效的模型通过假设原子独立振动且能量损失可以被视为一种局部效应来简化问题，本质上平滑了复杂的相互作用。相比之下，更复杂的冻结声子模型则生成了数千个特定的“快照”，在这些快照中，原子被冻结在略微不同的位置，从而模拟它们实际的相关运动。该团队在两种截然不同的材料上运行了这些模拟：一种是由轻碳原子组成的金刚石，另一种是含有锶和钛的较重晶体——钛酸锶，其振动方式更加复杂且不均匀。他们的目标是观察这种更简单、更快速的方法是否真的可以在不同条件下替代更严谨、更耗时的模型。

研究人员发现，虽然这两种方法通常一致，但当材料变重或样品变厚时，较简单的模型就开始失效。对于金刚石样品中的轻碳原子，快速模型与详细模拟之间的差异很小，尽管仍可测量。然而，当他们转向钛酸锶晶体时，这种差距显著扩大。在这种较重的材料中，简化模型忽略了关于原子如何协同振动以及如何吸收能量的重要细节，导致随着电子束向晶体深处传播，误差也随之增大。研究表明，这种差异不仅仅是一个微小的故障，而是一种系统性误差，在散射角较高处（即最常发现详细结构信息的区域）变得尤为明显。研究人员发现，简化模型无法捕捉到钛酸锶中的重原子彼此之间运动的具体方式，而详细的模拟则能自然地捕捉到这一细微差别。

该团队还调查了用于捕捉电子信号的探测器尺寸如何影响结果。他们发现，如果探测器足够大，能够收集广泛的散射电子，那么两种模型之间的差异就会显著缩小。在这种情况下，更简单、更快速的方法可以成为更费力模型的可靠替代品。然而，对于涉及重元素的精密工作或厚样品，研究表明，详细的冻结声子法仍然是更优的选择。作者总结道，虽然高效模型是许多应用中的宝贵工具，但科学家在使用它处理复杂的重元素系统时必须保持谨慎，因为它可能会引入微妙但显著的误差，从而导致对材料真实原子结构的误读。

技术摘要：透射电子显微镜（TEM）中的热弥散散射：复吸收势与冻结声子模型的比较

问题陈述
在高分辨率透射电子显微镜（HRTEM）和电子衍射中，通过晶格振动引起的非弹性散射（即热弥散散射，TDS）显著影响图像形成和衍射图样。虽然 TDS 表现为弥散背景，但对其进行建模对于准确模拟诸如高角环形暗场（HAADF）成像技术中的弹性通道损耗及强度分布至关重要。目前存在两种主要的建模方法：复吸收势（CAP）法，将非弹性散射视为弹性通道中的一个吸收项；以及冻结声子模型（FPM），通过模拟电子在代表原子热运动的静态晶格构型（“快照”）中的传播来进行模拟。虽然 CAP 方法计算效率较高，但其有效性条件（特别是关于原子相关性和各向异性方面）需要通过更具物理完整性的 FPM 进行严格验证。

方法论
作者使用多层切片模拟（multislice simulations）在 300 kV 加平行照明条件下，对三种特定模型进行了系统比较：

复吸收势 (CAP)： 利用参数化的复形式因子（Weickenmeier 和 Kohl；Peng 等）引入虚部势分量，从而使弹性波函数发生衰减。
具有相关运动的冻结声子模型 (MD-FPM)： 使用分子动力学（MD）模拟（利用 LAMMPS 软件及 Tersoff [针对金刚石] 或 deepMD [针对 STO] 原子间势）生成原子快照。这捕捉了原子的相关运动。
具有独立运动的冻协议声子模型 (Einstein FPM)： 基于爱因斯坦模型生成快照，假设原子进行独立的谐振运动。

研究通过两种不同的晶体系统对模型极限进行了测试：

金刚石 (C)： 一种轻元素立方晶体（3.5 nm × 3.5 nm × 50 nm），作为 TDS 吸收较小的参考系。
钛酸锶 (STO)： 一种含有重原子（Sr, Ti）且具有各向异性原子振动的重元素钙钛矿结构（3.1 nm × 3.1 nm × 30 nm），代表了更复杂的案例。

研究分析了总强度、TDS 分布以及弹性通道强度。通过对 MD-FPM 采用 100 个快照、Einstein FPM 采用 50 个快照进行平均，确保了收敛性。作者根据布拉格斑点强度的标准差建立了统计阈值，以区分显著的模型差异与数值噪声。此外，通过将衍射图样与不同直径（27, 55, 110 mrad）的圆形孔径进行卷积，模拟了探测器收集面积的影响。

关键结果

TDS 分布： MD-FPM 和 Einstein FPM 产生了定性相似的 TDS 图样。然而，由于存在原子声子相关性，MD-FPM 在布拉格斑点附近显示出更宽的强度分布，而这是独立的 Einstein 模型所不具备的特征。CAP 模型不产生空间 TDS 分布，仅产生一个弹性通道。
弹性通道差异：
- 金刚石： 三种模型在定性上表现出广泛的一致性。然而，MD-FPM 与其他两种模型之间的绝对差异在中心束处达到 $\approx 10^{-3} I_0$ ，并持续延伸至高阶劳厄区（HOLZ）区域（ $\theta \approx 110$ mrad）。
- STO： 差异显著增大。在低散射角下，绝对差异达到 $10^{-1.5} I_0$ 至 $10^{-2} I_0$ （百分比级误差）。即使在高散射角下，相关 MD-FPM 与 CAP/Einstein 模型之间的分歧仍然显著。
厚度依赖性： 对于 STO，相关 MD-FPM 与较简单模型（CAP 和 Einstein）之间的分歧随样品厚度的增加而逐渐增大，尤其是在高散射角下。这表明 CAP 和 Einstein 模型中的近似会在传播深度上累积误差。
探测器收集面积： 与大孔径探测器（110 mrad）进行卷积显著降低了模型间的差异。对于金刚石，最大差异降至 $\approx 10^{-4.5} I_0$ ，证实了较大的收集面积可以抑制模型间的分歧。

意义与主张
本文主张，虽然 CAP 模型是一种计算高效的替代方案，但它在全散射角范围内会引入不可忽视的系统误差。作者确定了 CAP 模型不再作为可靠近似的具体条件：

重元素： TDS 吸收对原子序数的强依赖性放大了重元素系统（如 STO）中 CAP 与 FPM 物理机制之间的失配。
各向异性振动： CAP 和 Einstein 模型依赖于各向同性的德拜-沃勒因子，无法捕捉复杂结构（如 STO）中存在的各向异性振动，而相关的 MD-FPM 则能自然地捕捉到这一点。
样品厚度： CAP 模型的误差会随着样品厚度的增加而累积，这在重元素样品中尤为明显。

作者得出结论，具有相关原子运动的 FPM (MD-FPM) 是多层切片框架内对 TDS 最完整的物理描述。然而，他们也确立了 CAP 模型仅在特定条件下才能作为 MD-FPM 的可靠替代品：轻元素系统、中等样品厚度以及较大的探测器收集角。这些发现为定量 TEM 模拟中选择吸收模型提供了实践指导。

技术摘要：透射电子显微镜（TEM）中的热弥散散射：复吸收势与冻结声子模型的比较

类似论文