Anisotropic Crystallization Kinetics and Interfacial Dynamics of Phase-Change Material Sb$_2$S$_3$ from Machine Learning Force Field Simulations

本研究利用机器学习力场揭示，Sb$_2$S$_3$受其准一维带状结构驱动表现出各向异性结晶，其界面控制生长动力学具有显著低于扩散的活化能，为优化其在数据存储和光子学应用中的性能提供了关键见解。

要点

想象你拥有一种名为**硫化锑（Sb₂S₃）**的神奇材料。这种材料就像计算机和光基技术的“变色龙”：它能瞬间在有序排列的固体晶体（如同整齐堆叠的图书馆）与杂乱无序的液体（如同散落的书堆）之间切换。这种来回切换的能力，使其成为存储数据和调控光线的理想选择。

然而，科学家很难在单个原子的层面上确切观察到这种切换是如何发生的。其过程太快、尺度太小，常规显微镜无法捕捉。为此，本文的研究人员构建了一个超级智能计算机大脑（称为机器学习力场），它充当高速、超精准的模拟引擎。这个“大脑”从复杂的物理计算中学习了原子相互作用的规则，使团队能够运行一段长达 40 纳秒的原子运动“大片”——这在原子世界中已是巨大的时间跨度。

以下是他们发现的成果，通过简单的类比进行解释：

1. “丝带”结构

请将这种材料的固体晶体形态想象成不是冰块，而是一束长而坚固的丝带。

• 快车道：原子沿着这些丝带的长度方向被紧密地粘合在一起（如同强共价键）。
• 慢车道：丝带之间的连接则弱得多，如同温柔的拥抱（范德华力）。

正因如此，材料在丝带方向上的生长速度最快。研究人员发现，晶体沿 [100] 方向（即丝带方向）的生长速度比其他方向快约4 倍。这就像拉上拉链：沿着齿牙方向能迅速扣合，但横向拉开布料则要困难得多。

2. 切换的“速度极限”

团队测量了发生两种不同过程所需的能量：

• 原子移动（扩散）：想象原子试图在拥挤的泳池中游泳。这是一项艰巨的工作。完成此过程所需的能量很高（约 1.16 至 1.56 eV）。
• 锁定到位（晶体生长）：想象原子到达晶体边缘并“咔哒”一声嵌入最终位置。这出乎意料地容易。所需能量要低得多（约 0.55 至 0.57 eV）。

重大发现：在许多其他类似材料中，“游泳”（原子移动）是缓慢且困难的部分，从而限制了速度。但对于 Sb₂S₃ 而言，“游泳”并非瓶颈。真正的瓶颈在于原子能多快附着到晶体边缘。该材料属于“界面控制型”。这就像一家工厂：工人（原子）可以非常快地跑向流水线，但机器（晶体边缘）只能以一定的速度将它们“咔哒”锁定到位。

3. “金发姑娘”温度

研究人员发现，该材料在极热或极冷时生长并非最快。

• 如果太热，原子过于躁动，难以结合。
• 如果太冷，原子过于迟钝，无法移动。
• 存在一个“甜蜜点”（约在熔点以下 100 度），此时生长效率最高。有趣的是，对于 Sb₂S₃ 而言，这个甜蜜点比其他常见材料更接近熔点，这意味着它只需更小的温度变化就能极快地切换状态。

4. “液态”记忆

即使材料熔化成液体，它也不会变成完全随机的汤。原子仍保留着对其丝带状结构的微弱记忆。它们保持了一些与固态形式相似的局部“舞步”（键角）。这就是为什么它重新凝固如此迅速且可靠的原因——原子无需学习新的舞蹈，只需记住它们原本就在做的舞步即可。

总结

简而言之，该论文利用强大的计算机模拟，观察了 Sb₂S₃ 从液态转变为固态的过程。他们发现：

1. 它沿其“丝带”方向生长最快。
2. 切换速度受限于原子在边缘锁定到位的速度，而非它们在液体中移动的速度。
3. 这使得它成为快速切换技术的极高效材料，因为它无需等待原子长途跋涉来形成晶体。

这项研究提供了该材料运作机制的清晰、原子级图谱，帮助工程师理解为何它能如此快速且可靠地切换状态。

技术摘要

技术摘要：基于机器学习力场模拟的相变材料 Sb2S3 各向异性结晶动力学与界面动力学

问题陈述
由于锑（Sb）在地壳中储量丰富、相稳定性强，且相较于 Ge2Sb2Te5 (GST) 等传统相变材料具有更优异的光学性能（高折射率对比度与低消光系数），硫化锑（Sb2S3）成为一种极具潜力的用于数据存储和光子应用的相变材料（PCM）。然而，实验技术在解析控制宽温域内结晶动力学的原子尺度机制、瞬态中间相以及局部成核事件方面存在固有局限。目前文献中尚缺乏对晶体至液相转变过程中结构演变的全面理解，特别是关于依赖于晶面的生长速率和界面动力学的研究。

方法论
为了克服从头算分子动力学（MD）在时间尺度和系统尺寸上的局限性，作者基于张量矩势（MTP）框架开发了一种机器学习力场（MLFF）。

• 力场开发：采用主动学习（自举）方法，起始于一组晶体结构和 20 组非晶体结构。训练集通过从 MD 模拟中采样构型并利用密度泛函理论（DFT）计算（使用 VASP 软件，采用 PBE-GGA 泛函及范德华修正）进行标记，从而迭代扩展。最终数据集包含 2,094 个结构，能量和力的均方根误差（RMSE）分别达到 7 meV/atom 和 0.17 eV/Å。
• 验证：MTP 模型针对实验 X 射线衍射（XRD）数据和从头算MD 结果进行了验证。该模型成功复现了晶体相和液相的晶格参数、径向分布函数（RDFs）以及原子间距离。
• 模拟设置：利用 LAMMPS-MLIP 接口，在包含 7,680 个原子的双相共存（DPC）模型上进行了大规模 MD 模拟。该设置允许在 600 K 至 1200 K 的温度范围内模拟长达 40 ns 的晶体生长过程。研究考察了 [100]、[010] 和 [001] 三个特定晶面，以分析各向异性生长。

主要贡献与结果

1. 结构表征：

   • 晶体相与液相：研究证实，晶体相具有独特的配位几何结构（三价 Sb 位于畸变的三角锥中，五价 Sb 位于四方锥中），并具有特定的键角。相比之下，液相表现出围绕 90°分布的宽键角分布（0–180°），且平均 Sb 配位数降低（从晶体中的 3/5 降至液相中的 2–4）。
   • 局部有序：尽管长程有序丧失，液相仍保留了类似畸变金字塔或准四面体 Sb–S 配位的局部基序，这促进了快速且可逆的相变。

2. 各向异性晶体生长动力学：

   • 晶面依赖性：生长速率（$V_g$）表现出高度的各向异性。[100] 晶面表现出最快的生长速度，其次是 [010] 和 [001]。在 875 K 时，[100] 的生长速率约为 [010] 的 1.6 倍，为 [001] 的 2.8 倍。
   • 机制：沿 [100] 方向的快速生长归因于 Sb2S3 晶体的准一维带状结构，其中强 Sb-S 共价键沿带轴连续形成，降低了表面能。带与带之间的侧向相互作用较弱（范德华力）。
   • 热力学：最大生长速率出现在 875 K，对应于相对于模拟熔点（975 K）约 100 K 的最佳过冷度（$\Delta T_{opt}$）。该$\Delta T_{opt}$显著低于 GST（270 K）和 GeTe（210 K），表明 Sb2S3 具有更好的热响应性。

3. 活化能与生长机制：

   • 生长与扩散：利用 Wilson-Frenkel 模型，计算得出各晶面的晶体生长活化能（$\Delta E_{a}^{WF}$）为 0.55–0.57 eV。相比之下，界面液层中原子扩散的活化能（$\Delta E_a$）显著更高，范围在 1.16 至 1.56 eV 之间。
   • 界面控制：发现$\Delta E_{a}^{WF} < \Delta E_a$，表明 Sb2S3 的结晶是界面控制而非扩散限制。固 - 液界面的原子附着在能量上优于长程原子输运。这与主要以扩散限制为主的 GST 和 GeTe 形成对比。

意义
本文证明，所开发的基于 MTP 的 MLFF 能够准确捕捉 Sb2S3 结构演变和结晶动力学的原子起源，其模拟尺度（7,680 个原子，40 ns）是从头算方法无法达到的。该研究首次提供了 Sb2S3 中依赖于晶面生长的详细机理视图，揭示了其结晶是由界面处的局部有序而非体相扩散驱动的。这些关于各向异性生长速率、低最佳过冷度以及界面控制动力学的见解，为优化 Sb2S3 的功能奠定了基础，特别是针对下一代光电器件中的开关速度、可靠性和能效。