Finite Temperature Stacking Fault Stability in Random and Locally Ordered CoCrNi beyond the Harmonic Approximation

本研究证明，尽管局部化学有序在 CoCrNi 中跨温度稳定了正层错能，但有限温度下的非谐效应未能热稳定随机固溶体 CoCrNi 所预测的负层错能，从而解决了先前基于谐近似密度泛函理论预测与实验观测之间的差异。

要点

想象一种名为CoCrNi（钴、铬和镍的混合物）的金属合金，将其视为一个巨大而拥挤的舞池。在这场舞蹈中，原子是舞者，它们通常以一种非常有序、重复的模式移动，这种模式被称为“面心立方”（FCC）结构。

有时，在一次舞蹈动作（变形）中，舞池的一部分会出现一点“故障”或滑移。在材料科学中，这被称为堆垛层错。你可以把它想象成一块稍微皱起或移位的地毯。

科学家们一直提出的核心问题是：这种“皱起”是保持微小且可控，还是会失控地扩散开来？

谜团：“负能量”问题

长期以来，计算机模拟（使用一种称为密度泛函理论 DFT 的方法）预测，在这些原子完全随机混合的情况下（称为随机固溶体或 RSS），这种“皱起”应该是不稳定的。

• 类比： 想象试图握住一根具有负张力的橡皮筋。它不会弹回，而是想要无限拉伸。
• 预测： 计算机指出，产生这种缺陷所需的“能量”是负的。这意味着原子会想要无限分离，形成一个巨大且无限的滑移。
• 现实： 现实世界的实验表明，滑移确实会发生，但它保持有限（它会在某个宽度处停止）。橡皮筋并没有无限拉伸，而是停住了。

科学家们提出了两种理论来解决计算机与现实之间的这种不匹配：

1. 理论 A（热量）： 也许房间的热量（温度）充当了稳定器，使橡皮筋停止拉伸。
2. 理论 B（有序）： 也许原子实际上并非完全随机。也许它们拥有小小的“朋友圈”或局部团簇（称为局部化学有序或 LCO），这些团簇自然地固定住了滑移。

本文做了什么

本文的作者旨在平息这场争论。他们使用了一个超精准的 AI 模型（“神经网络势函数”）来模拟原子，但有一个关键转折：他们不仅仅将原子视为轻微振动的僵硬小球（旧的“谐振动”方式）。他们将原子视为摇摆不定、混乱的舞者，彼此剧烈碰撞（“非谐振动”方式）。这更像现实生活，当原子受热时会变得混乱。

发现：实际上发生了什么？

1. “热稳定化”理论是错误的
作者首先测试了随机（RSS）混合。

• 旧观点： 他们认为加热会使“橡皮筋”停止拉伸。
• 新发现： 当他们考虑到热原子混乱、摇摆的振动时，他们发现了相反的情况。随着温度升高，“橡皮筋”实际上想要更拉伸。
• 结果： 在完全随机的混合中，堆垛层错并未被热量稳定。它仍然不稳定，想要无限扩展。那些声称“热量能解决它”的旧计算机模型，忽略了原子振动混乱的现实。

2. “局部有序”理论是英雄
接下来，他们观察了原子形成小小“朋友圈”（LCO）的混合情况。

• 发现： 即使在高温下，这些局部团簇也充当了安全网。它们产生了一种“恢复力”（像正常的橡皮筋一样），将滑移拉回到一个特定的、有限的尺寸。
• 结果： “皱起”保持微小且稳定，就像在真实实验中一样。正是局部化学有序阻止了滑移失控。

3. 位错舞蹈（证明）
为了绝对确定，他们运行了包含数百万个原子的巨大模拟，观察一条“位错”（缺陷线）如何在金属中移动。

• 在随机混合中： 位错分裂并持续扩散，直到撞到模拟盒的边缘。这是失控的混乱。
• 在有序混合中： 位错分裂，但随后停止了。它找到了一个舒适、稳定的宽度并停留在那里。

结论

本文得出结论：我们在 CoCrNi 合金中看到稳定、有限的“皱起”，并不是因为热量力挽狂澜。而是因为原子实际上并非完全随机。 它们拥有局部的有序区域，这些区域充当锚点，保持材料的稳定性。

简单来说：

• 随机混合： 就像一群陌生人互相推挤；如果一个人滑倒，整个群体可能会崩溃并无限扩散。
• 有序混合： 就像一群朋友手拉手组成小团体；如果一个人滑倒，团体会将他们拉回，将混乱控制在范围内。

这项研究证明，这些“朋友圈”（局部化学有序）是这种金属即使在受热时依然如此坚韧和稳定的真正原因。

技术摘要

技术摘要：超越谐近似的随机与局部有序 CoCrNi 中有限温度堆垛层错稳定性

问题陈述
关于等原子比 CoCrNi 中熵合金（MPEA）中的堆垛层错稳定性，理论预测与实验观测之间存在显著差异。实验研究一致观察到有限的堆垛层错宽度（SFW），这意味着存在正的本征堆垛层错能（ISFE）。然而，针对随机固溶体（RSS）模型的标准密度泛函理论（DFT）计算预测在 0 K 时平均 ISFE 值为负，这在理论上暗示了部分位错的连续、无界分离。

为了解决这一矛盾，提出了两种主要解释：

1. 有限温度稳定化：谐近似或准谐近似表明，ISFE 随温度升高而增加，可能在高温下变为正值。
2. 局部化学有序（LCO）：普遍存在于 MPEA 中的化学有序畴的存在，甚至可能在 0 K 时将 ISFE 移至正值。

然而，现有文献存在关键空白：以往针对 RSS 的有限温度计算往往忽略了显式的非谐振动效应，且 LCO 状态下 ISFE 的温度依赖性尚未被探索。

方法论
本研究采用多方位计算方法，探究 CoCrNi 在 RSS 和 LCO 构型中随温度变化的广义堆垛层错能（GSFE）：

• 原子间势：利用具有近量子精度的神经网络势（NNP）高效采样多样化的局部原子环境，克服了经验嵌入原子法（EAM）势的精度局限。
• 非谐自由能计算：研究采用完全非谐投影平均力积分器（PAFI）方法。与谐近似不同，PAFI 通过沿最小能量路径积分热平均平均力，显式地计入了非谐振动贡献。
• 系综采样：
  • RSS：由 50 个独立的 450 原子超胞组成的系综，其中 Co、Cr 和 Ni 呈随机等摩尔分布。
  • LCO：由 50 个独立的 2,520 原子超胞组成的系综，包含纳米尺度的化学有序畴，通过混合分子动力学/蒙特卡洛（MD/MC）协议生成。
• 验证：进行了大规模分子动力学（MD）模拟（约 167 万个原子），以模拟两种化学状态下的刃型位错结构，直接观测堆垛层错宽度的演化。
• 分析：计算了 Warren-Cowley（WC）参数以量化化学有序度，并采用间隔共同近邻分析（ICNA）技术进行缺陷识别。

关键结果

1. 谐近似在 RSS 中的失效：

   • 谐计算重现了 ISFE 随温度升高的趋势。
   • 相比之下，完全非谐 PAFI 计算揭示，RSS CoCrNi 的 ISFE 随温度升高而降低。
   • 因此，RSS ISFE 在所研究的温度范围（高达 1000 K）内保持负值，表明有限温度效应并未在热力学上稳定随机状态下的堆垛层错。

2. 通过局部化学有序（LCO）实现稳定化：

   • LCO 畴的存在将平均 ISFE 移至正值（0 K 时约为 38 mJ/m²）。
   • 虽然由于非谐振动，LCO ISFE 随温度单调下降，但在几乎整个范围（0–1000 K）内仍保持正值，仅在接近 1000 K 时趋近于零。
   • 这种正 ISFE 提供了必要的热力学恢复力，以平衡部分位错之间的排斥弹性相互作用，从而允许有限的 SFW 存在。

3. 位错行为：

   • RSS 状态：MD 模拟证实了无界位错解离，其中堆垛层错扩展至周期性边界，这与负 ISFE 一致。
   • LCO 状态：模拟显示位错核心被限制，在初始平衡后达到稳定的有限堆垛层错宽度。

4. 层错能的热响应：

   • 在 RSS 状态下，不稳定堆垛层错能（USFE）表现出比 ISFE 显著更强的热下降（约 2.2 倍），表明由于化学无序，这些能量发生了去耦合。
   • 在 LCO 状态下，该比率降至约 1.35，表明局部化学有序部分恢复了不稳定构型与本征层错构型之间更耦合的热响应。

意义与主张
本文得出结论：CoCrNi 在有限温度下出现有限堆垛层错宽度的根源是局部化学有序（LCO），而非随机合金状态的热稳定化。研究表明，以往对谐近似的依赖导致了对 RSS CoCrNi 热稳定化的错误预测。通过显式包含非谐效应，作者表明即使在高温下，RSS 堆垛层错在热力学上仍是不稳定的（负 ISFE）。

研究结果确立了 LCO 是提供热力学恢复力以稳定有限部分位错分离的关键机制。此外，该工作强调化学无序从根本上改变了 USFE 和 ISFE 的相对热响应，这是在模拟复杂 MPEA 力学行为时必须考虑的因素。结果得到了自由能计算和大规模位错模拟的双重证实，为这些合金中的堆垛层错稳定性提供了一致的物理图像。