Interatomic potentials for platinum

想象一下，将铂（Pt）视为一位极其 sophisticated、高性能的运动员。它强壮，即使在极端高温下也不易生锈，并被用于从催化转化器到医疗设备的各种领域。为了理解这位“运动员”在压力、热量或压强下的行为，科学家使用计算机模拟。但要运行这些模拟，他们需要一本规则手册——一套指导计算机如何描述铂的每一个原子与其邻近原子相互作用的指令集。这本规则手册被称为原子间势。

长期以来，可用的铂原子间势手册就像一些陈旧、磨损的地图。它们存在某些错误：预测的熔点温度不正确，或者预测某些内部化学键过于容易断裂。

在本文中，作者（Koju、Li 和 Mishin）决定为铂编写两本全新且高度精确的规则手册。以下是他们工作的简要分解：

1. “训练”（无需人工猜测）

通常，科学家在制作这些规则手册时，会参考现实世界的实验来验证其正确性。然而，该团队决定完全采用数字化方法。他们使用一种超精确的量子物理方法（称为 DFT）生成了一个庞大的“训练数据库”。

类比：想象训练一个机器人下棋。你不是向它展示人类进行的真实棋局，而是让机器人对抗一个基于完美数学的对手，进行数百万局对弈。机器人纯粹从数学中学习规则，而不是通过观察人类。
结果：他们利用这些纯数学数据训练了两个新模型。在训练阶段，他们没有使用任何实验测量数据。

2. 两本新规则手册

作者创建了两类不同风格的规则手册：

ADP 模型（“灵活”规则手册）：这是对一种较旧的标准方法的升级。旧方法就像一条规则，规定“原子只关心邻居离自己有多近”。而新的 ADP 版本增加了一个转折：“原子也关心邻居形成的角度。”这就像说，一个人不仅关心谁站在自己旁边，还关心谁站在自己的左边或右边。这使得该模型在预测金属如何弯曲和振动方面非常出色。
MT 模型（“适配”规则手册）：该模型最初是为金刚石或硅等材料（具有非常刚性、方向性强的化学键）设计的。作者将这个刚性模型“拉伸”以适应像铂这样的金属。
- 类比：想象一本为刚性木椅设计的规则手册。作者对其进行了修改，使其能够描述柔软、可挤压的金属枕头。令人惊讶的是，这种“拉伸”后的规则手册变得极其精确，有时甚至比 ADP 模型还要好。

3. 结果：谁胜出？

团队将这两个新规则手册与旧版本（“磨损的地图”）以及超精确的量子数学计算进行了测试对比。

熔点：旧规则手册预测的铂熔点温度低了几百度。新的 ADP 规则手册将熔点预测得几乎完全准确（误差在极小的度数范围内）。MT 规则手册也非常接近，只是略高。
断裂与弯曲：旧规则手册无法预测产生“缺陷”（缺失一个原子）所需的能量，也无法预测原子层相互滑移（就像洗牌一样）所需的能量。新模型修正了这些错误，更准确地预测了断裂或滑移金属所需的能量。
振动：当金属振动时（就像吉他弦一样），新模型预测的“音符”（频率）比旧模型好得多。

4. 权衡：速度 vs. 精度

这里有一个限制。

ADP 模型就像一辆快速的跑车。它非常精确，且运行模拟的速度很快。
MT 模型则像一辆高科技的重型坦克。它极其精确（有时甚至优于 ADP），但运行速度非常慢。使用 MT 模型运行模拟所需的时间是使用 ADP 模型的 100 多倍，因为它必须不断计算原子之间复杂的角度。

5. 为什么这很重要（根据论文所述）

作者指出，虽然 MT 模型对于纯铂来说速度较慢，但它可能是未来材料的“缺失环节”。

类比：想象你有一本关于水（液体）的规则手册和一本关于混凝土（固体）的规则手册。但如果你需要模拟一种一半是水、一半是混凝土的材料（比如湿水泥）呢？单独使用任何一本规则手册都无法很好地工作。
MT 模型的特殊之处在于，它能够使用相同的数学语言处理金属（如铂）和共价材料（如碳或硅）。
文中提到的具体应用：论文明确指出，这个新模型可用于模拟铂硅化物（用于微芯片）和基于铂的抗癌药物（其中铂与氮结合）。它使科学家能够模拟这些混合材料在原子层面的行为，而以前这非常困难。

总结

作者为铂原子构建了两本全新且高度精确的数字规则手册。他们利用纯数学而非实验数据对它们进行了训练。两者都比旧版本好得多，特别是在预测熔点和金属断裂方式方面。其中一个是快速的（ADP），另一个虽然慢但极其通用（MT）。那个通用的模型可能是模拟复杂材料的关键，这些材料混合了金属与其他元素，例如你手机中的芯片或特定的药物。

技术摘要：铂的原子间势函数

问题陈述
经典原子尺度模拟对于理解铂（Pt）的力学行为和微观结构演变至关重要，特别是在极端条件下。然而，这些模拟的可靠性取决于所用原子间势函数的准确性。现有的铂势函数主要采用嵌入原子法（EAM）格式，表现出显著的局限性。这些局限性包括堆垛层错能的低估，以及熔点与实验值偏差数百开尔文。此类缺陷可能导致定量精度降低，甚至在某些情况下引发非物理的模拟行为。此外，现有势函数难以捕捉控制缺陷构型和表面重构的微妙电子效应。

方法论
为了解决这些局限性，作者为铂开发了两种新的原子间势函数：一种采用角依赖势（ADP）格式，另一种采用修正 Tersoff（MT）格式。

训练数据库：两种势函数均完全基于第一性原理密度泛函理论（DFT）计算生成的参考数据库进行训练，计算使用了维也纳从头算模拟包（VASP）。训练过程中未使用任何实验数据。数据库包括：
- 各向同性变形下各种晶体结构（面心立方 FCC、体心立方 BCC、密排六方 HCP、简单立方 SC、A15 结构、金刚石立方）的状态方程（EOS）。
- 面心立方结构下的单轴拉伸/压缩和剪切变形。
- 通过 nudged elastic band（NEB）方法计算的空位迁移最小能量路径（MEP）。
- 表面能、孪晶界能和 $\gamma$ 面。
- 点缺陷形成能（空位和间隙原子）及二聚体能量。
- 声子色散关系和态密度。
势函数格式：
- ADP：EAM 的扩展形式，包含角依赖相互作用（偶极矢量和四极张量）以考虑非中心键合。它使用了 36 个拟合参数。
- MT：最初为共价元素（如 Si、C）设计，通过将截断半径扩展至包含多个配位壳层（高达 6.57 Å）而适配于金属铂。这使得该模型能够捕捉类似于 MEAM 和 ADP 的多体相互作用。它使用了 16 个拟合参数。
优化：两种势函数的参数均基于 Nelder-Mead 单纯形法的模拟退火算法进行优化，以最小化势函数预测值与 DFT 参考数据库之间的加权均方偏差。
验证：新势函数针对 DFT 数据库和实验数据进行了严格测试。为保持一致性，两项近期发表的 EAM 势函数（Zhou 等人的 EAM1 和 O'Brien 等人的 EAM2）也使用相同的计算协议进行了重新评估。

主要贡献与结果
本研究对新的 ADP 和 MT 势函数与现有 EAM 势函数、DFT 计算及实验数据进行了全面比较。

晶格与热力学性质：
- 两种新势函数均以高精度预测了平衡晶格常数和结合能。
- 熔点（ $T_m$ ）：ADP 势函数重现的实验熔点误差在 0.2% 以内（2041 K 对比实验值约 2045 K）。MT 势函数高估了 $T_m$ 7.3%（2195 K）。相比之下，EAM1 和 EAM2 势函数显著低估了 $T_m$ ，偏差达数百开尔文。
- 热膨胀：MT 势函数在高达 1200 K 的范围内与实验热膨胀数据表现出极好的一致性，优于 ADP 和 EAM 势函数。
- 声子：新势函数重现声子色散关系的能力显著优于 EAM 势函数，尽管所有势函数在高对称点处相对于实验值均略微低估了频率。
缺陷性质：
- 空位：新势函数正确预测了处于实验范围内的空位形成能，而 DFT 计算（PBE 泛函）严重低估了该值。ADP 和 MT 势函数均准确重现了空位迁移能垒，而 EAM 势函数则显著低估了该值。
- 间隙原子：DFT 计算确定八面体间隙原子是铂的基态，这与其他面心立方金属中典型的 $\langle 100 \rangle$ 哑铃构型偏好不同。虽然新势函数正确预测了 $\langle 100 \rangle$ 哑铃为基态（未能捕捉导致八面体偏好的微妙电子效应），但它们提供的能量值比 EAM 势函数更准确，后者系统性地低估了间隙原子能量。
- 表面：所有势函数均正确排序了表面能（ $\gamma_{111} < \gamma_{100} < \gamma_{110}$ ）。然而，四种势函数（ADP、MT、EAM1、EAM2）均未能正确预测与表面重构相关的能量降低（例如 (110) 面的 $(1 \times 2)$ 重构），这是一个已知的缺陷。
堆垛层错与孪晶界：
- 与 EAM 势函数相比，新势函数在预测堆垛层错能（SFE）和孪晶界能方面表现出显著改进。
- ADP 势函数重现了 DFT 计算的孪晶界能（142 mJ m $^{-2}$ ），一致性极佳。MT 势函数高估了约 50 mJ m $^{-2}$ ，而 EAM 势函数则严重低估（例如，EAM1 给出 46 mJ m $^{-2}$ ）。
- 同样，ADP 势函数在 SFE 方面与 DFT 最接近，而 EAM 势函数严重低估了这些值。
结构稳定性：
- 所有势函数均正确识别面心立方结构为基态。非平衡结构（BCC、HCP 等）相对能量的奇偶性图显示，ADP 和 MT 势函数相对于 DFT 表现出无偏散的分布，而 EAM 势函数则显示出低估能量的趋势。

意义与主张
该论文声称，新开发的 ADP 和 MT 势函数在铂的模拟中比现有 EAM 势函数提供了更高的精度，特别是在熔点、堆垛层错能和缺陷迁移能垒方面。

效率与精度：作者指出了计算成本上的权衡。MT 势函数虽然精度极高且所需参数较少（16 个对比 ADP 的 36 个），但由于在嵌套循环中显式计算键角，其计算成本高昂（在分子动力学模拟中比 ADP 慢 $10^2$ 倍以上）。ADP 势函数则在高精度和计算效率之间取得了平衡。
未来应用：作者强调了 MT 模型在弥合金属势函数与共价势函数差距方面的潜力。由于 MT 形式可以在单一分析框架内描述金属键和共价键，它被提议作为模拟混合键合系统的候选模型，例如铂硅化物（PtSi、Pt $_2$ Si）以及用于抗癌药物的铂 - 氮配位化合物。作者建议，尽管 MT 模型目前存在计算开销，但其处理混合键合化学的能力有望实现这些具有技术相关性材料的大规模模拟。

总之，这项工作提供了一套稳健的、基于第一性原理训练的势函数，纠正了文献中的特定缺陷，为研究人员模拟铂及其化合物的机械和物理行为提供了改进的工具。