Partially reactive force field for the UiO-66 metal-organic framework

Este artículo presenta nb-UiO-FF, un nuevo campo de fuerzas parcialmente reactivo que modela con precisión las propiedades estructurales, mecánicas y de defectos del marco metal-orgánico UiO-66 y permite simulaciones de dinámica molecular de sus mecanismos de síntesis solvotermal y autoensamblaje.

Lo esencial

Imagina un marco metal-orgánico (MOF) como UiO-66 como un rompecabezas 3D gigante y microscópico. Está construido a partir de "nodos" metálicos (como cúmulos de zirconio) que actúan como las esquinas, y "enlaces" orgánicos (como anillos de benceno) que actúan como las varillas de conexión. A los científicos les encantan estos rompecabezas porque son increíblemente resistentes y pueden ajustarse para atrapar gases, administrar fármacos o acelerar reacciones químicas.

Sin embargo, hay un problema: averiguar exactamente cómo se ensamblan estos rompecabezas desde el principio es como intentar ver una película reproducida a la velocidad de la luz. Los enlaces químicos que mantienen unidos al metal y a los enlaces son difíciles de simular en una computadora. La mayoría de los modelos informáticos tratan estas conexiones como un pegamento permanente; no pueden mostrar la aplicación del pegamento, el encaje de las piezas, ni siquiera lo que sucede si falta una pieza.

La Solución: Un Campo de Fuerza de "Pegamento Inteligente"
En este artículo, los autores introducen una nueva herramienta informática llamada nb-UiO-FF. Piensa en esto como un nuevo conjunto de reglas para un juego de simulación que permite que las piezas del rompecabezas sean "parcialmente reactivas".

Así es como lo hicieron funcionar, utilizando algunas analogías simples:

• Los átomos "ficticios" (Las manos invisibles):
  En el mundo real, el nodo metálico de zirconio tiene una carga eléctrica compleja que atrae a los enlaces en direcciones específicas. Los modelos informáticos estándar luchan por imitar esto sin volverse desordenados. Los autores resolvieron esto adjuntando átomos "ficticios" invisibles (como pequeños títeres magnéticos) a los nodos metálicos. Estos títeres actúan como manos invisibles que mantienen los enlaces en la forma y orientación correctas, asegurando que el rompecabezas se construya de la manera correcta sin necesidad de cálculos complejos y pesados.

• El "Potencial de Morse" (El resorte elástico):
  Por lo general, los modelos informáticos tratan los enlaces como palos rígidos. Si los tiras, se rompen instantáneamente. Los autores reemplazaron estos palos rígidos con un potencial de Morse, que actúa más como un resorte elástico. Esto permite que la simulación muestre cómo el metal y el enlace se estiran, oscilan e incluso se unen o se separan dinámicamente. Esto es crucial para observar el "nacimiento" del material.

Lo que probaron
Los autores no solo construyeron la herramienta; la sometieron a una prueba de estrés rigurosa para asegurarse de que fuera confiable:

1. El rompecabezas perfecto: Verificaron si la herramienta podía recrear la forma exacta de un cristal UiO-66 perfecto. Coincidió con las mediciones del mundo real casi perfectamente (dentro de una fracción diminuta de un porcentaje).
2. El rompecabezas empapado: Probaron la herramienta con el cristal empapado en dos líquidos diferentes utilizados para fabricarlo (DMF y etanol). El modelo mostró que el cristal se mantiene fuerte y no se desmorona cuando está húmedo.
3. El rompecabezas roto: Los cristales del mundo real a menudo tienen piezas faltantes (defectos). Los autores eliminaron intencionalmente enlaces o nodos completos en la simulación. La herramienta mostró con éxito que el cristal aún podía mantener su forma incluso con estos agujeros, al igual que lo hace el material real.
4. El rompecabezas rebotando: Probaron cuánto podían apretar el cristal antes de que se deformara. Los resultados coincidieron con cálculos de física de alto nivel, demostrando que el modelo comprende la resistencia del material.
5. El rompecabezas primo: Probaron la herramienta en una versión ligeramente más grande del rompecabezas (UiO-67) y también funcionó allí, demostrando que las reglas son flexibles.

Observando la magia suceder
La parte más emocionante del artículo es utilizar esta nueva herramienta para observar el proceso de autoensamblaje. Imagina soltar todas las piezas del rompecabezas (nodos metálicos y enlaces) y el disolvente líquido en una caja y presionar "reproducir".

La simulación mostró cómo las piezas se deslizaban y comenzaban lentamente a pegarse.

• Vieron cómo los nodos metálicos y los enlaces se encontraban y formaban los bloques de construcción iniciales.
• Notaron que a veces las piezas se quedan atascadas en posiciones "incorrectas" (trampas cinéticas), como una pieza de rompecabezas que encaja holgadamente pero no es del todo correcta.
• Observaron que el proceso es lento; las piezas son pesadas y se mueven con lentitud, por lo que el rompecabezas completo no se ensambla completamente en el tiempo que simularon.

La conclusión
Este artículo presenta un nuevo modelo informático altamente preciso que actúa como un "microscopio inteligente" para el material UiO-66. Puede simular la estructura del material, su resistencia y su capacidad para manejar defectos. Lo más importante es que es la primera herramienta de su tipo que puede simular de manera realista el proceso dinámico de construcción del material desde cero, ayudando a los científicos a entender cómo nacen estos materiales asombrosos y cómo controlar sus imperfecciones.

Resumen técnico

Resumen Técnico: Campo de Fuerza Parcialmente Reactivo para el Marco Metal-Orgánico UiO-66

Planteamiento del Problema
Los marcos metal-orgánicos (MOF), particularmente el UiO-66, son ampliamente estudiados por su capacidad de ajuste estructural y su estabilidad. Sin embargo, comprender sus mecanismos de síntesis y la formación de defectos sigue siendo un desafío debido a la falta de modelos computacionales capaces de simular procesos reactivos. Los campos de fuerza existentes (por ejemplo, UFF4MOF, BTW-FF, MOF-FF) dependen típicamente de una conectividad predefinida, lo que los hace incapaces de modelar la ruptura y formación de enlaces, algo esencial para estudiar el autoensamblaje de MOF y la ingeniería de defectos. Aunque existen campos de fuerza reactivos como Reax-FF y potenciales de aprendizaje automático, enfrentan desafíos relacionados con la precisión para MOF, el costo computacional o la dificultad de expandirlos para incluir disolventes y moduladores. Además, los modelos específicos para UiO-66 a menudo luchan por representar con precisión la forma hidroxilada o requieren configuraciones de átomos ficticios computacionalmente costosas (por ejemplo, 48 átomos ficticios por nodo).

Metodología
Los autores desarrollaron nb-UiO-FF, un campo de fuerza parcialmente reactivo diseñado para modelar el marco UiO-66 y su análogo isoreticular, UiO-67. La metodología se basa en un enfoque previamente establecido para ZIF-8, utilizando una combinación de interacciones no enlazantes y átomos ficticios catiónicos (CDA).

• Construcción del Campo de Fuerza: El modelo reemplaza los potenciales enlazantes rígidos entre el nodo metálico y los ligandos con un potencial de Morse no enlazante para permitir la ruptura y formación dinámica de enlaces.
• Átomos Ficticios (CDA): Para reproducir la distribución de carga anisotrópica de los átomos de Zr y mantener la geometría de coordinación correcta sin términos enlazantes, se introducen cuatro átomos ficticios (D) por centro de Zr. Estos están conectados al átomo de Zr mediante enlaces y ángulos armónicos fuertes. La carga de +1.968 del átomo de Zr se distribuye completamente entre estos cuatro átomos ficticios (+0.492 cada uno), haciendo que el átomo de Zr sea eléctricamente neutro.
• Parametrización: Los parámetros del campo de fuerza se derivaron partiendo de un campo de fuerza no reactivo de Yang et al., ajustado para el nodo hidroxilado Zr6O4(OH)4. Los parámetros enlazantes para el cluster metálico se tomaron de modelos reportados, mientras que las longitudes de enlace y los ángulos de equilibrio se obtuvieron de resultados de DFT de Valenzano et al. Los parámetros del potencial de Morse ($D_0$ y $\alpha$) para la interacción Zr–Oligando se optimizaron explorando sistemáticamente el espacio de parámetros para equilibrar la precisión en la densidad numérica y las longitudes de enlace locales (específicamente la función de distribución radial Zr–Oligando) frente a referencias de DFT y experimentales.
• Protocolos de Simulación: Todas las simulaciones se realizaron utilizando LAMMPS. La validación incluyó análisis estructural, pruebas de estabilidad en disolventes (DMF y etanol) y cálculos de propiedades mecánicas. Los estudios de autoensamblaje utilizaron un sistema de 16 clusters de Zr, 96 ligandos BDC y 2136 moléculas de DMF, simulados durante 6 ns en el ensemble NPT.

Resultados Clave

• Precisión Estructural: El nb-UiO-FF reproduce el parámetro de red experimental del UiO-66 prístino con una desviación de aproximadamente 0.046% (20.71 Å frente a 20.7004 Å). Las funciones de distribución radial (FDR) para pares atómicos clave (Zr–$\mu_3$O, Zr–$\mu_3$OH, Zr–Zr, Zr–Oligando) muestran un acuerdo razonable con datos de DFT y experimentales, con desviaciones generalmente dentro del 10%.
• Estabilidad de Disolvente y Defectos: El campo de fuerza mantiene la integridad del marco en presencia de disolventes DMF y etanol, así como en estructuras defectuosas que contienen ligandos faltantes y nodos faltantes. En modelos defectuosos, la neutralidad de carga se preservó encapsulando centros subcoordinados con iones cloruro o fragmentos BDC monodentados.
• Propiedades Mecánicas: Las constantes elásticas calculadas ($C_{11}$, $C_{12}$, $C_{44}$) y las propiedades mecánicas derivadas (Módulo de compresibilidad: 44 GPa, Módulo de Young: 40 GPa) están en buen acuerdo con las referencias de DFT, aunque ligeramente más altas, probablemente debido a la rigidez introducida por las fuertes interacciones enlazantes entre Zr y los átomos ficticios.
• Transferibilidad: El modelo reproduce con éxito la estructura del marco isoreticular UiO-67, demostrando transferibilidad a través de diferentes longitudes de ligandos.
• Dinámica de Autoensamblaje: En simulaciones dinámicas, el campo de fuerza capturó con éxito la formación progresiva de enlaces de coordinación entre clusters metálicos y ligandos orgánicos. El estudio identificó motivos estructurales transitorios, incluidas unidades de construcción secundarias (SBUs) iniciales y configuraciones de unión no ideales (estados atrapados cinéticamente). Sin embargo, la convergencia completa del número de coordinación y la formación de grandes agregados iónicos no se lograron dentro de la ventana de simulación de 6 ns, lo que sugiere la necesidad de métodos de muestreo mejorado para escalas de tiempo más largas.

Significado y Afirmaciones
El artículo afirma que nb-UiO-FF proporciona una descripción consistente y físicamente fiable de UiO-66 y UiO-67, cerrando la brecha entre el modelado estructural estático y la dinámica reactiva. Al utilizar un conjunto reducido de átomos ficticios (cuatro por Zr) en comparación con modelos anteriores, los autores logran un equilibrio entre la eficiencia del modelo y la capacidad de simular procesos reactivos. El campo de fuerza permite el estudio de mecanismos de autoensamblaje y formación de defectos, que son inaccesibles para los campos de fuerza tradicionales de conectividad fija. Los autores posicionan este trabajo como una base para investigar las etapas tempranas de la síntesis de MOF y la ingeniería de defectos puntuales, reconociendo al mismo tiempo que se requiere más trabajo con esquemas de muestreo mejorado (por ejemplo, metadinámica) para superar completamente las trampas cinéticas y modelar el autoensamblaje completo.