Computer Simulation of the Growth of a Metal-Organic Framework Proto-crystal at Constant Chemical Potential

Resumen Técnico: Simulación Computacional del Crecimiento de un Protocristal de Marco Organometálico a Potencial Químico Constante

Planteamiento del Problema
El diseño racional de los protocolos de síntesis de Marcos Organometálicos (MOF) se ve actualmente obstaculizado por la falta de una comprensión fundamental de los mecanismos de autoensamblaje a nivel molecular. Aunque trabajos computacionales y experimentales previos han investigado la nucleación de MOFs (específicamente el ZIF-8), la fase de crecimiento subsiguiente ha permanecido en gran medida inexplorada mediante simulaciones por computadora. Esta brecha existe debido a obstáculos computacionales significativos: modelar el crecimiento requiere un núcleo preformado (demasiado grande para métodos ab initio) mientras se mantiene simultáneamente una concentración constante de reactantes libres en la solución para imitar las condiciones de supersaturación experimental. Los ensambles estándar de dinámica molecular (MD) con números de partículas fijos no pueden capturar esto, ya que el agotamiento de los reactantes ocurre naturalmente a medida que el cristal crece, alterando el comportamiento del crecimiento.

Metodología
Para abordar estos desafíos, los autores emplearon un enfoque de simulación híbrido que combina la Dinámica Molecular de Potencial Químico Constante (C$\mu$MD) con una técnica de inserción de partículas.

• Configuración del Sistema: Se colocó una placa de protocristal de ZIF-8 preformada en el centro de una caja de simulación que contiene un solvente (DMSO) y reactantes (Zn$^{2+}$ y 2-metilimidazolato, MIm$^-$).
• Implementación de C$\mu$MD: La caja de simulación se dividió en una Región de Transición (TR) donde ocurre el crecimiento, Regiones de Control (CR) que flanquean el cristal donde las densidades iónicas se mantienen constantes mediante una fuerza externa, y un Reservorio (Res) que suministra iones.
• Inserción de Partículas: Para evitar el agotamiento del reservorio a lo largo del tiempo, se integró un esquema de Monte Carlo de Gran Canónico (GCMC). Cuando el recuento de iones en el reservorio caía por debajo de un umbral, se realizaban inserciones basadas en el criterio de Metropolis para restaurar las proporciones estequiométricas (1 Zn$^{2+}$ : 2 MIm$^-$) y la neutralidad de carga.
• Campo de Fuerza: Se utilizó el campo de fuerza nb-ZIF-FF, que presenta un potencial de Morse para las interacciones Zn–N para permitir la formación de enlaces reversibles.
• Condiciones: Las simulaciones se realizaron a dos temperaturas (298 K y 320 K) y tres concentraciones de reactantes ([Zn$^{2+}$] = 0.1 M, 0.2 M y 0.4 M).
• Análisis: El crecimiento se rastreó identificando el mayor clúster conectado de iones. El análisis estructural consistió en perfiles de densidad lineales, divergencia de Kullback–Leibler (KL) para medir la desviación de la red cristalina ideal y análisis de la distribución del tamaño de los anillos (de 3 a 8 miembros).

Resultados Clave

1. Mecanismo de Crecimiento y Desorden: El proceso de crecimiento implica la adhesión tanto de monómeros como de oligómeros. La formación de oligómeros es más prevalente a concentraciones más altas. Las capas resultantes no son perfectamente cristalinas, sino que exhiben una interfaz amorfa con un perfil de densidad que preserva la estructura de pseudo-poros de la semilla cristalina.
2. Efecto de las Condiciones Termodinámicas:
   • Concentración: El aumento de la concentración de reactantes conduce a tasas de crecimiento más altas y a un crecimiento significativamente más desordenado. A concentraciones bajas (0.1 M), el crecimiento es más ordenado y simétrico. A concentraciones altas (0.4 M), la superficie se vuelve irregular y dendrítica.
   • Temperatura: El aumento de la temperatura (de 298 K a 320 K) también aumenta la tasa de crecimiento y el grado de desorden estructural.
3. Evolución Estructural (Análisis de Anillos):
   • Tamaños de Anillo: La distribución de los tamaños de los anillos (formados por enlaces alternos metal-ligando) depende altamente de las condiciones. A concentraciones bajas, los anillos de 3 miembros son los más comunes. A medida que la concentración y la temperatura aumentan, la prevalencia de anillos más grandes (de 5, 7 y 8 miembros) aumenta, mientras que los anillos de 4 miembros siguen siendo dominantes en condiciones intermedias.
   • Comparación con la Nucleación: A diferencia de la fase de nucleación (estudiada previamente por el mismo grupo), la cual está dominada por anillos de 4 miembros, la fase de crecimiento favorece la formación de anillos más grandes, particularmente anillos de 6 miembros, lo que sugiere una progresión hacia la topología de sodalita del ZIF-8, aunque a través de intermediarios defectuosos transitorios.
4. Cinética:
   • Tasas de Estado Estacionario: Las tasas de crecimiento se computaron para el régimen de estado estacionario. Las tasas aumentan con la concentración y la temperatura.
   • Mecanismo de Control: La relación entre la tasa de crecimiento y la concentración en el seno (bulk) es no lineal. Esto sugiere que el crecimiento está controlado por la adsorción (específicamente el proceso activado de formación de enlaces Zn–ligando y disociación del solvente) en lugar de por la difusión.
   • Estequiometría: A pesar de la naturaleza amorfa de las capas en crecimiento, la superficie permanece eléctricamente neutra con una estricta relación 1:2 de Zn$^{2+}$:ligando.

Significancia y Reivindicaciones
Los autores afirman que este trabajo proporciona el primer estudio computacional a nivel atómico de la fase de crecimiento del autoensamblaje de MOFs, yendo más allá de la fase de nucleación estudiada previamente. Al implementar con éxito C$\mu$MD con inserción de partículas, demuestran un método viable para simular el crecimiento de cristales bajo un potencial químico constante, superando las limitaciones del agotamiento de reactantes en la MD estándar.

El estudio revela que el crecimiento del ZIF-8 procede a través de una fase intermedia amorfa caracterizada por distribuciones específicas de tamaño de anillo que difieren de la fase de nucleación. Los hallazgos sugieren que las condiciones de síntesis (concentración y temperatura) ajustan directamente el grado de cristalinidad frente a la de amorfo en el producto final. La dependencia no lineal de las tasas de crecimiento respecto a la concentración ofrece una visión mecánica de que la cinética de adsorción, en lugar de la difusión, gobierna el proceso bajo estas condiciones. Los autores postulan que este marco computacional puede aplicarse para investigar el crecimiento de otros MOFs bajo diversas condiciones de síntesis para dilucidar aún más los mecanismos de autoensamblaje.