Validation of an Ab Initio-Informed Electronic Stopping Model for Large-Scale Atomistic Simulations

当高能粒子撞击固体材料时，它们并不仅仅像台球一样弹开；它们会撕裂原子晶格，撞开原子，并引发一场混乱的碰撞级联。这种被称为辐射损伤的过程，是核反应堆和航天器中材料面临的无声敌人，因为强烈的辐射会逐渐降低结构的完整性。为了预测一种材料能维持多久，科学家必须精确理解这些高速运动的粒子在穿透固体时是如何损失能量的。几十年来，模拟这一过程的标准方法是将原子视为在邻近原子影响下运动的经典粒子，同时忽略围绕它们的电子。这种简化对于许多情况都适用，但当粒子移动得足够快，足以从这些电子中夺走能量时（这种现象被称为电子阻尼），它就会失效。如果没有一种精确的方法来计算这种能量损失，关于粒子穿透深度或其造成的损伤程度的预测将变得不可靠，往往依赖于武断的猜测而非基础物理学。

芬兰阿尔托大学的一个研究小组通过利用现实世界的数据测试一种更先进的模型，在填补这一空白方面迈出了重要一步。他们专注于一种被称为“统一双温度模型”的具体方法，该模型试图描述能量如何在运动的原子与电子云之间流动，且无需使用任何可调节的、人为设定的数值。该模型并非通过猜测离子感受到的摩擦力，而是直接从量子力学的基本定律出发，通过观察离子移动时电子密度的变化来推导其行为。研究人员想要知道，这种理论上完美的模型是否真的能预测现实实验中的情况，或者在面对物理材料的复杂现实时是否会崩溃。为了验证这一点，他们模拟了一束携带100千电子伏特能量的硅离子，射穿一片仅53纳米厚的硅晶体薄片。随后，他们将计算机生成的模拟结果与实验室中的实际测量值进行了对比，在实验中，离子穿过完全相同的箔片，并在29厘米外的探测器上被捕捉。

实验设计得对模拟细节极其敏感。通过旋转硅晶体，研究人员可以迫使离子沿着不同的路径旅行：有些离子会沿着原子行之间的开放通道平滑滑动，而另一些则会随机撞击原子晶格。这种变化至关重要，因为离子损失能量的多少在很大程度上取决于其轨迹；通过晶体通道的平滑路径造成的能量损失比通过原子的混乱路径要小。以往的模型通常将所有路径视为相同，或者使用基于速度发生突变的简单规则，从而忽略了这些微妙的方向性效应。然而，新模型考虑到了电子云如何根据离子相对于其邻居的位置和移动速度做出不同的响应。在模拟过程中，研究人员还必须考虑自然形成在硅箔表面的薄氧化层。他们发现，在箔片的进入侧和退出侧各包含一层22埃厚的氧化层，对于重现实验中观察到的散射模式至关重要，因为这一层轻微地引导了离子，使其分布范围恰好能与观测数据相匹配。

当团队将模拟结果与实验测量值进行比较时，一致性令人惊叹。新模型不仅成功预测了离子损失能量的平均量，还预测了能量分布的具体形状。对于沿晶体通道直线行驶的离子，模型预测了一个窄能量范围；而对于采取随机路径的离子，则显示出更宽的分布。最重要的是，模拟捕捉到了在特定角度下数据中呈现出的独特的双峰结构，这一特征源于保持平滑路径的离子与被散射离子的混合。旧的、更简单的模型无法重现这种复杂性，往往预测的能量损失要么过高，要么过低，或者完全缺失了数据的分布范围。研究人员还检查了模型在离子向随机方向运动（即晶体结构不提供平滑通道）时的表现。即使在这种情况下，该模型依然表现稳健，无需针对射束的具体角度进行调整或重新校准，便能与实验数据相吻合。

这项研究证实了统一双温度模型是预测辐射损伤材料的一种可靠工具。与以往依赖任意阈值或经验调整的方法不同，这种方法使用的电子摩擦描述植根于第一性原理计算，这意味着它是从材料本身的基本属性中推导出来的。研究人员证明，该方法可以准确追踪高速离子穿过晶体时的能量耗散，捕捉到平滑路径与混乱路径之间的细微差别。这种精确度对于设计能够承受太空和核技术严酷环境的材料至关重要。通过证明该模型在电子阻尼机制（即离子将能量损失给电子的特定阶段）中有效，该团队为未来的大规模模拟奠定了坚实的基础。这些模拟现在可以用来以以往难以企及的信心度来预测辐射损伤，为理解材料在持续不断的辐射轰击下如何演变提供了更清晰的路径。