Reactive Force Field for P/Sn/I System: Atomistic Insight into the Early Stage of Black Phosphorus and Phosphorene Synthesis Process

技术摘要：P/Sn/I 体系的反应力场及黑磷合成的原子级洞察

问题陈述
黑磷（BP）及其二维对应物磷烯（phosphorene）由于其高载流子迁移率和可调控的带隙，被视为下一代电子器件中极具前景的材料。然而，通过化学气相输运法（CVT）合成高质量、大面积黑磷薄膜仍面临挑战。该过程通常使用红磷并添加锡（Sn）和碘（I₂）作为添加剂，但其在原子水平上的机制尚未得到充分理解。关于添加剂的具体作用、稳定气相物种的身份、磷输运机制以及前成核中间体的性质，仍存在关键的不确定性。由于高温和高反应性环境，实验上难以观察到这些早期阶段的动力学过程；而标准的密度泛函理论（DFT）受限于小规模体系和静态 0K 配置，无法捕捉模拟生长机制所需的动态、有限温度反应网络。

方法论
为了解决这些局限性，作者开发了一种专门针对 P/Sn/I 体系的新型反应力场，命名为 P/Sn/I-2512。该力场通过使用 DFT（采用 PBE 泛函、TZ2P 基组及 Grimme 的 D3 色散校正）生成的全面量子力学（QM）数据集进行训练。训练集包括：

• Sn–Sn、Sn–I、Sn–P、I–P 和 I–I 键的键解离能。
• 各种分子构型的键角和扭转剖面。
• 关键物种（SnI₂、SnI₄、PI₃）的生成热。
• 凝聚态锡（$\alpha$-Sn 和 $\beta$-Sn）的状态方程以及锡簇（Sn$_x$, $x=2-6$）的形成能。

力场参数使用 Amsterdam Modeling Suite 中的协方差矩阵自适应进化策略（CMA-ES）进行了优化。通过将 ReaxFF 结果与 DFT 数据在键解离、状态方程和反应能方面进行对比进行验证，结果显示出高度的一致性。

随后，利用 LAMMPS 进行了反应分子动力学（ReaxFF-MD）模拟，以研究 923 K 下的早期前成核过程。研究采用了带有 Nosé-Hoover 热浴的 NVT 系综。模拟通过隔离特定变量来开展：

1. 系统验证：对不同密度（有效压力分别为 0.01、0.1 和 1 atm）下的 Sn/I₂ 和 P/I₂ 系统进行 MD 模拟，以验证关键物种（如 SnI₂、SnI₄、PI₂ 和 PI₃）的形成。
2. 碘的作用：对比了三个系统（仅 P；P + Sn；P + Sn + I₂）以分离碘对磷重组的影响。
3. 密度效应：在三种密度区间（有效压力分别为 5、10 和 100 atm）下对 P/Sn/I₂ 系统进行模拟，以模拟从稀薄蒸气到凝聚相的转变，模拟实验中观察到黑磷形成的条件。

关键结果

• 力场验证：P/Sn/I-2512 力场准确重现了 DFT 参考数据中的键解离、键角、扭转势垒以及锡晶体和团簇的状态方程。它成功预测了高温下 SnI₂、SnI₄ 和 PI₃ 的形成。
• Sn/I₂ 和 P/I₂ 中的密度依赖性：在 Sn/I₂ 系统中，增加密度使组分从元素物种转向复杂的配位环境，促进了 SnI₄ 和 Sn₂I₄（SnI₂ 的二聚体）连同 SnI₂ 一起的形成。在 P/I₂ 系统中，更高的密度促进了 PI₃ 和聚合 PI$_x$ 物种的形成，这与实验中难以分离纯 PI₃ 的现象一致。
• 碘在重组中的作用：模拟显示，在缺失碘的情况下，磷原子主要保持原子状态或形成极小的团簇。加入 I₂ 显著增加了 P–P 重组速率，导致了二元 P$_x$I$_y$ 基团和瞬态三元 Sn$_x$P$_y$I$_z$ 化合物的形成。这证实了碘在溶解磷并促进输运方面的作用。
• 高密度下的前成核路径：在低密度（5 atm）下，系统主要由原子磷和小团簇（P₂, P₃）组成。在中等密度（10 atm）下，较大的团簇（P₄, P₅）开始形成。在高密度（100 atm）下，系统迅速凝聚成高阶 P$_x$ 团簇（x = 8–12）和大体积三元 Sn$_x$P$_y$I$_z$ 化合物。
• 中间体的识别：在最高密度下，模拟产生了位于 Sn$_x$P$_y$I$_z$ 团簇边缘的凝聚态 Hittorf 型磷结构。此外，还观察到了由碘稳定的类黑磷（BP-like）六角形结构，这些结构被认为是潜在的成核种子。

意义与主张
本文声称提供了一个用于研究黑磷/磷烯蒸气合成中生长条件和添加剂效应的实用原子级模型。其主要贡献包括：

1. 开发了经过验证的力场：P/Sn/I-2512 力场提供了一个可靠工具，用于模拟涉及磷、锡和碘的复杂反应系统，克服了 DFT 在尺寸和时间尺度上的限制。
2. 机制洞察：研究阐明了碘和系统密度是控制磷重组的关键因素。研究表明，碘不仅是输运剂，还通过形成中间体 P$_x$I$_y$ 和三元 Sn$_x$P$_y$I$_z$ 物种，积极促进了 P–P 重组。
3. 前成核模型：作者提出了一个特定的前成核路径，即高密度条件结合碘，促进了磷向 Hittorf 型中间体和类黑磷六角形结构的凝聚。这些物种被认为是随后在冷却阶段生长为黑磷的前驱体。

该工作并非提出新的实验方案，而是提供了一个理论框架，用以解释现有的实验观察，特别是黑磷 CVT 合成中高压（密度）的必要性以及碘的关键作用。