Microscopic Mechanism of the Thermal Amorphization of ZIF-4 and Melting of ZIF-zni Revealed via Molecular Dynamics and Machine Learning Techniques

Resumen Técnico: Mecanismo Microscópico de la Amorfización Térmica de ZIF-4 y la Fusión de ZIF-zni

Planteamiento del Problema
Los Marcos Metal-Orgánicos (MOFs) amorfos ofrecen una combinación única de porosidad intrínseca y robustez mecánica; sin embargo, los mecanismos microscópicos que gobiernan su formación a partir de precursores cristalinos siguen estando insuficientemente comprendidos. Específicamente, las transiciones de fase ordenado–desordenado en los Marcos de Imidazolato Zeolítico (ZIFs), tales como la amorfización térmica de ZIF-4 y la fusión de ZIF-zni, presentan desafíos significativos para el modelado computacional. Si bien existen datos experimentales, las cuestiones relativas a los mecanismos específicos de amorfización, la naturaleza de los estados intermedios y el comportamiento más allá de las distancias de primer vecino permanecen sin respuesta. Los métodos estándar de dinámica molecular ab initio (AIMD) suelen estar limitados por el tamaño del sistema y las escalas de tiempo, mientras que los campos de fuerza tradicionales frecuentemente fallan al capturar la reactividad metal–ligando necesaria para modelar la ruptura y formación de enlaces durante estas transiciones.

Metodología
Los autores emplearon un enfoque híbrido que combina la dinámica molecular (MD) clásica con el aprendizaje automático (ML) para investigar las transiciones de fase de dos polimorfos de Zn(C3H3N2)2: el poroso ZIF-4 y el denso ZIF-zni.

1. Campo de Fuerza: Las simulaciones se realizaron utilizando el campo de fuerza nb-ZIF-FF (non-bonded ZIF-FF). Este campo de fuerza reactivo utiliza potenciales de Morse para modelar los enlaces metal–ligando, lo que permite la ruptura y formación de enlaces. Incorpora un modelo de "átomo ficticio" (dummy atom) para imponer la topología tetraédrica correcta alrededor de los centros de Zn²⁺, asegurando la reproducción precisa de las propiedades estructurales, mecánicas y termodinámicas para las fases cristalina, amorfa (ZIF_a) y líquida (ZIF_liq).
2. Generación de Estructuras: Las estructuras amorfas (ZIF_a) se generaron mediante recocido simulado (calentamiento a 900 K y enfriamiento) a partir de los precursores ZIF-4 y ZIF-zni. Las estructuras resultantes fueron validadas frente a la densidad experimental, el módulo de compresibilidad y las funciones de distribución radial (RDF).
3. Clasificación por Aprendizaje Automático: Para distinguir entre las cuatro fases (ZIF-4, ZIF-zni, ZIF_a, ZIF_liq) a nivel atómico local, los autores entrenaron una red neuronal.
   • Descriptores: Se computaron Funciones de Simetría de Behler-Parrinello (BPSF) para los entornos centrados en Zn²⁺, las cuales sirvieron como entradas no sesgadas.
   • Arquitectura: Una red neuronal con 12 nodos de entrada (funciones de simetría), una capa oculta (6 nodos) y 4 nodos de salida (probabilidades de fase).
   • Rendimiento: El modelo alcanzó una precisión del 90.3% en la clasificación de entornos locales, distinguiendo exitosamente entre las dos fases desordenadas (ZIF_a y ZIF_liq) basándose únicamente en la información estructural.
4. Protocolo de Simulación: Se realizaron simulaciones de MD NPT a gran escala utilizando LAMMPS. Las temperaturas de transición se determinaron identificando condiciones de equilibrio donde las fases coexisten sin un crecimiento neto de una sobre la otra.

Resultos Clave

• Termodinámica: El estudio determinó temperaturas de transición de equilibrio de 420 K para la amorfización de ZIF-4 y 625 K para la fusión de ZIF-zni. Aunque estos valores difieren de las observaciones experimentales (589 K y 863 K, respectivamente), el modelo reproduce correctamente las tendencias de estabilidad termodinámica, confirmando que el ZIF-4 es metaestable en condiciones ambientales. Las entalpías y entropías calculadas se alinean bien con los valores de referencia.
• Mecanismo de Amorfización de ZIF-4: La transición de ZIF-4 cristalino a ZIF_a amorfo es un proceso de dos pasos:
  1. Formación de una fase intermedia de baja densidad y tipo líquido (similar a ZIF_liq).
  2. Transformación de esta fase intermedia en la fase amorfa de alta densidad (ZIF_a).
     El proceso implica una disminución en la conectividad seguida de una reorganización. La formación de la fase amorfa es anisotrópica, propagándose más rápido a lo largo del eje x. Esta anisotropía se correlaciona con las distribuciones de densidad local y la labilidad de los enlaces de coordinación; los enlaces orientados en la dirección x son más propensos a romperse, facilitando el intercambio de conectividad.
• Mecanso de Fusión de ZIF-zni: La fusión de ZIF-zni también sigue una vía de dos pasos, pero en orden inverso de evolución de densidad comparado con ZIF-4:
  1. Formación de una fase de tipo amorfo a partir del cristal.
  2. Transición a la fase líquida (ZIF_liq).
     La fase líquida se nuclea dentro del seno de la fase desordenada intermedia.
• Disparadores Microscópicos: El inicio de la amorfización está vinculado a cambios específicos de conectividad más que a la simple ruptura de enlaces. Los defectos de Zn tricoordinados ocurren frecuentemente, pero son reversibles. El evento crítico que dispara la transición de fase es la formación de un nuevo enlace entre un centro de Zn²⁺ y una molécula de ligando distinta (cambiando la conectividad de la red), en lugar de una simple rotación o reformación del enlace original. Este evento viene acompañado de un aumento repentino en la población de entornos tipo líquido identificados por la red neuronal.
• Evolución Estructural: Se observa que la formación de anillos de 3 miembros (Zn–ligando–Zn) se correlaciona con el crecimiento de la fase desordenada, lo que sugiere que estas características topológicas son propias del estado amorfo, mientras que el cristal contiene principalmente anillos de 4, 6 y 8 miembros.

Significancia y Reivindicaciones
El artículo afirma proporcionar una comprensión microscópica detallada de la amorfización térmica y la fusión en los ZIFs, abordando las lagunas dejadas por estudios previos que ignoraron los estados intermedios o carecieron de un tamaño de sistema suficiente. Sus contribuciones principales son:

1. Validación de nb-ZIF-FF: Demostrar que este campo de fuerza puede reproducir con precisión las estructuras y propiedades de las fases cristalina, amorfa y líquida de los ZIF, incluyendo distribuciones angulares sutiles que otros campos de fuerza reactivos (como ReaxFF) no logran capturar.
2. Dinámica Aumentada por ML: Establecer una metodología donde una red neuronal, entrenada con funciones de simetría, puede rastrear la evolución temporal de las transiciones de fase con alta resolución espacial, distinguiendo entre fases desordenadas estructuralmente similares.
3. Perspectiva Mecanística: Identificar que tanto la amorfización como la fusión son procesos de dos pasos que involucran fases intermedias con densidades distintas, y que la amorfización de ZIF-4 es un proceso anisotrópico impulsado por la labilidad de enlaces de coordinación específicos.

Los autores concluyen que su trabajo ofrece un marco para estudiar procesos reactivos en MOFs, lo que podría ayudar en el diseño de MOFs amorfos para aplicaciones industriales como la purificación de agua y la catálisis, aunque no proponen nuevos protocolos experimentales ni aplicaciones futuras específicas más allá del alcance de la comprensión de los mecanismos fundamentales.