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

Resumen Técnico: Campo de Fuerza Reactivo para el Sistema P/Sn/I e Información Atómica sobre la Síntesis de Fósforo Negro

Planteamiento del Problema
El fósforo negro (BP) y su contraparte bidimensional, el fosforeno, son materiales prometedores para la electrónica de próxima generación debido a su alta movilidad de portadores y brechas de banda sintonizables. Sin embargo, la síntesis de películas de BP de alta calidad y gran área mediante el Transporte Químico de Vapor (CVT) sigue siendo un desafío. El proceso, que típicamente utiliza fósforo rojo con aditivos de estaño (Sn) yodio (I₂), no se comprende completamente a nivel atómico. Persisten incertidumbres clave respecto a los roles específicos de los aditivos, la identidad de las especies estables en fase gaseosa, los mecanismos de transporte de fósforo y la naturaleza de los intermediarios de pre-nucleación. La observación experimental de estas dinámicas tempranas es difícil debido a las altas temperaturas y los entornos reactivos, mientras que la Teoría del Funcional de la Densidad (DFT) estándar está limitada a sistemas pequeños y configuraciones estáticas a 0 K, fallando al capturar las redes de reacción dinámicas y de temperatura finita requeridas para modelar los mecanismos de crecimiento.

Metodología
Para abordar estas limitaciones, los autores desarrollaron un nuevo campo de fuerza reactivo (ReaxFF) específicamente para el sistema P/Sn/I, designado como P/Sn/I-2512. Este campo de fuerza fue entrenado contra un conjunto de datos de mecánica cuántica (QM) exhaustivo generado mediante DFT (funcional PBE con base TZ2P y corrección de dispersión de Grimme D3). El conjunto de entrenamiento incluyó:

• Energías de disociación de enlaces para Sn–Sn, Sn–I, Sn–P, I–P e I–I.
• Perfiles de ángulos de valencia y torsión para diversas configuraciones moleculares.
• Entalpías de formación para especies clave (SnI₂, SnI₄, PI₃).
• Ecuaciones de estado para el estaño en fase condensada ($\alpha$-Sn y $\beta$-Sn) y energías de formación de cúmulos de estaño (Sn$_x$, $x=2-6$).

Los parámetros del campo de fuerza fueron optimizados utilizando la Estrategia Evolutiva de Adaptación de la Matriz de Covarianza (CMA-ES) dentro de la Amsterdam Modeling Suite. La validación se realizó comparando los resultados de ReaxFF contra los datos de DFT para la disociación de enlaces, ecuaciones de estado y energías de reacción, mostrando una fuerte concordancia.

Posteriormente, se realizaron simulaciones de Dinámica Molecular Reactiva (ReaxFF-MD) utilizando LAMMPS para investigar los procesos de pre-nucleación en etapas tempranas a 923 K. El estudio empleó el ensamble NVT con un termostato de Nosé-Hoover. Las simulaciones se estructuraron para aislar variables específicas:

1. Validación del Sistema: Simulaciones de MD de los sistemas Sn/I₂ y P/I₂ a diferentes densidades (presiones efectivas de 0.01, 0.1 y 1 atm) para verificar la formación de especies clave como SnI₂, SnI₄, PI₂ y PI₃.
2. Rol del Yodo: Se compararon tres sistemas (P solo; P + Sn; P + Sn + I₂) para aislar el efecto del yodo en la recombinación de fósforo.
3. Efectos de Densidad: El sistema P/Sn/I₂ se simuló bajo tres regímenes de densidad (presiones efectivas de 5, 10 y 100 atm) para modelar la transición de vapor diluido a fases condensadas, imitando las condiciones donde se observa experimentalmente que ocurre la formación de BP.

Resultos Clave

• Validación del Campo de Fuerza: El campo de fuerza P/Sn/I-2512 reprodujo con precisión los datos de referencia de DFT para la disociación de enlaces, ángulos de valencia, barreras torsionales y las ecuaciones de estado para los cristales de estaño y sus cúmulos. Predijo con éxito la formación de SnI₂, SnI₄ y PI₃ a altas temperaturas.
• Dependencia de la Densidad en Sn/I₂ y P/I₂: En el sistema Sn/I₂, el aumento de la densidad desplazó la composición desde especies elementales hacia entornos de coordinación complejos, promoviendo la formación de SnI₄ y Sn₂I₄ (dímeros de SnI₂) junto con SnI₂. En el sistema P/I₂, las densidades más altas facilitaron la formación de PI₃ y especies poliméricas de PI$_x$, lo cual es consistente con las dificultades experimentales para aislar PI₃ puro.
• Rol del Yodo en la Recombinación: Las simulaciones revelaron que, en ausencia de yodo, los átomos de fósforo permanecían mayoritariamente atómicos o formaban cúmulos mínimos. La adición de I₂ aumentó significamente la tasa de recombinación de P–P, llevando a la formación de motivos binarios de P$_x$I$_y$ y compuestos ternarios transitorios de Sn$_x$P$_y$I$_z$. Esto confirma el rol del yodo en la solubilización del fósforo y la facilitación del transporte.
• Vía de Pre-nucleación a Alta Densidad: A bajas densidades (5 atm), el sistema consistía principalmente de fósforo atómico y pequeños cúmulos (P₂, P₃). A densidad media (10 atm), comenzaron a formarse cúmulos más grandes (P₄, P₅). A alta densidad (100 atm), el sistema se condensó rápidamente en cúmulos de alto orden de P$_x$ (x = 8–12) y grandes compuestos ternarios de Sn$_x$P$_y$I$_z$.
• Identificación de Intermediarios: En la densidad más alta, las simulaciones produjeron estructuras de fósforo tipo Hittorf condensadas en los bordes de los cúmulos de Sn$_x$P$_y$I$_z$. Además, se observaron motivos hexagonales de tipo BP estabilizados por yodo, identificados como posibles semillas de nucleación.

Significancia y Reivindicaciones
El artículo afirma proporcionar un modelo atómico práctico para estudiar las condiciones de crecimiento y los efectos de los aditivos en la síntesis de vapor de BP/fosforeno. Las principales contribuciones son:

1. Desarrollo de un Campo de Fuerza Validado: El campo de fuerza P/Sn/I-2512 ofrece una herramienta fiable para simular sistemas reactivos complejos que involucran fósforo, estaño e yodo, superando las limitaciones de tamaño y escala temporal de la DFT.
2. Información Mecanística: El estudio elucida que el yodo y la densidad del sistema son los factores de control para la recombinación de fósforo. Demuestra que el yodo no es meramente un agente de transporte, sino que promueve activamente la recombinación de P–P mediante la formación de intermediarios de tipo P$_x$I$_y$ y compuestos ternarios de Sn$_x$P$_y$I$_z$.
3. Modelo de Pre-nucleación: Los autores proponen una vía de pre-nucleación específica donde las condiciones de alta densidad, combinadas con el yodo, facilitan la condensación del fósforo en intermediarios tipo Hittorf y hexágonos tipo BP. Estos componentes se sugieren como precursores para el crecimiento posterior del fósforo negro durante la etapa de enfriamiento.

El trabajo no propone nuevos protocolos experimentales, sino que ofrece un marco teórico para interpretar observaciones experimentales existentes, específicamente la necesidad de alta presión (densidad) y el papel crítico del yodo en la síntesis de CVT de fósforo negro.