MOSP: A User-interface Package for Simulating Metal Nanoparticle Structure and Reactivity under Operando Conditions

Resumen Técnico de MOSP: Un paquete de interfaz de usuario para la simulación de la estructura y reactividad de nanopartículas metálicas bajo condiciones operando

Planteamiento del Problema
Las nanopartículas (NPs) metálicas son fundamentales en la catálisis heterogénea, donde su reactividad catalítica está intrínsecamente ligada a su estructura atómica y al ensamblaje de los sitios activos. Mientras que las simulaciones cinéticas tradicionales, como el modelado microcinético de campo medio (MKM) y el Monte Carlo cinético (KMC), han predicho con éxito el comportamiento catalítico basándose en estructuras de modelos fijos, las observaciones experimentales recientes in situ revelan que las estructuras de las NPs se reconstruyen dinámicamente en respuesta a los cambios en las condiciones de reacción. Los enfoques tradicionales que dependen de estructuras estáticas no logran contabilizar estos cambios de forma inducidos por las condiciones de reacción, creando una brecha en la comprensión de las relaciones estructura-propiedad dinámicas, esenciales para el diseño racional de catalizadores.

Metodología
Para abordar esta limitación, los autores presentan MOSP (Multi-scale Operando Simulation Package), una interfaz gráfica de usuario (GUI) de código abierto basada en Python que integra dos modelos computacionales distintos:

1. Modelo de Reconstrucción Estructural Multiescala (MSR):

   • Basado en la teoría de Wulff, este módulo predice la forma de equilibrio de las NPs metálicas bajo atmósferas de reacción específicas.
   • Calcula la tensión interfacial ($\gamma_{hkl}^{int}$) de cada faceta cristalina modificando la tensión superficial de vacío mediante términos que contabilizan la cobertura de gases, las energías de adsorción y las interacciones laterales entre adsorbatos.
   • Las coberturas de adsorbato se determinan resolviendo la isoterma de adsorción de Fowler-Guggenheim, incorporando presiones parciales, temperaturas y cambios de entropía.
   • El modelo predice cuantitativamente las estructuras de equilibrio para partículas mayores de 3 nm, considerando efectos de tamaño finito.

2. Modelo de Monte Carlo Cinético (KMC):

   • Este módulo simula la dinámica de reacción sobre las estructuras específicas de las NPs generadas por el módulo MSR (o importando estructuras existentes).
   • El estado del sistema se define por una configuración superficial de sitios discretos, cada uno caracterizado por elemento, posición y valor de ocupación.
   • La simulación emplea el Método Directo (DM) del algoritmo de Gillespie para resolver la ecuación maestra de la evolución temporal.
   • Los eventos incluyen adsorción, desorción, difusión y reacciones superficiales. Las constantes de velocidad se derivan de la teoría de colisiones (para adsorción/desorción) y de la teoría del estado de transición (para reacciones/difusión), utilizando relaciones de Brønsted-Evans-Polanyi (BEP) y números de coordinación generalizados (GCN) para mapear los perfiles de energía específicos del sitio.

Contribuciones Clave

• Flujo de Trabajo Integrado: MOSP proporciona un marco unificado para acoplar la predicción estructural termodinámica (MSR) con la simulación de reacción cinética (KMC), permitiendo la exploración de las relaciones estructura-actividad dinámicas.
• Accesibilidad para el Usuario: El paquete cuenta con una GUI intuitiva construida con wxpython, que permite a los usuarios introducir parámetros complejos (por ejemplo, matrices de interacción lateral, relaciones de escala, redes de reacción) sin necesidad de escribir código. Incluye paneles de visualización interactiva para las estructuras de las NPs (usando pyOpenGL) y los resultados de la simulación (usando pandas y matplotlib).
• Modularidad: Aunque está diseñado para una operación acoplada, los módulos MSR y KMC pueden funcionar de forma independiente. El MSR puede exportar estructuras de equilibrio y el KMC puede cargar estructuras preexistentes en formato XYZ.
• Código Abierto: El código fuente está disponible públicamente, lo que facilita la adopción y extensión por parte de la comunidad.

Resultados y Demostraciones
El artículo valida MOSP a través de dos ejemplos ilustrativos:

1. Oxidación de CO sobre nanopartículas de Pt:

   • El módulo MSR reconstruyó con éxito las NPs de Pt (radio de 2.5 nm) bajo presiones variables (0–2000 Pa) a 800 K.
   • Los resultados mostraron una evolución estructural dependiente de la presión: a medida que la presión aumentaba, la fracción de facetas (111) disminuía mientras que las (110) aumentaban, desplazando la forma de la NP de un octaedro truncado a un dodecaedro casi rómbico.
   • Las simulaciones KMC posteriores (10 millones de pasos) sobre la estructura de 500 Pa revelaron que los bordes y las esquinas eran los principales contribuyentes a la actividad, con coberturas calculadas de CO (0.01) y O (0.11) y una frecuencia de rotación (TOF) de $8.33 \times 10^5$ sitio$^{-1}$s$^{-1}$.
   • Las pruebas de rendimiento indicaron que el tiempo de ejecución de MSR escala con el número de átomos ($N_{átomos}^{2.4}$), mientras que el KMC escala de forma casi lineal ($N_{átomos}^{0.93}$).

2. Reacción de Desplazamiento de Gas de Agua (WGS) sobre nanopartículas de Cu:

   • Se realizó una simulación KMC sobre una NP de Cu de 2190 átomos con forma de octaedro truncado a 650 K y 4000 Pa.
   • El modelo incluyó cinco adsorbatos y catorce eventos siguiendo el mecanismo de carboxilo.
   • Los resultados mostraron al CO como el adsorbato dominante (cobertura 0.67).
   • Crucialmente, la simulación destacó el acoplamiento cinético a nivel atómico: se identificaron los sitios de baja coordinación como favorables para la formación de CO$_2$, mientras que los sitios de alta coordinación favorecieron la formación de H$_2$, a pesar de tener TOFs generales similares para ambos productos.

Significado
El artículo posiciona a MOSP como una herramienta eficiente para que los investigadores predigan la estructura y las propiedades catalíticas de las NPs metálicas bajo condiciones operando. Al cerrar la brecha entre los supuestos estructurales estáticos y las realidades experimentales dinámicas, MOSP ayuda en el cribado de las condiciones de reacción óptimas y el diseño racional de catalizadores. Los autores enfatizan que el paquete está diseñado para ejecutarse en portátiles estándar, haciendo que la simulación operando avanzada sea accesible sin depender de clústeres de computación de alto rendimiento para el cribado inicial.