Thermal Stability of Two-Dimensional Crystals with Extended OH Hydrogen-Bonded Chains

Resumen Técnico: Estabilidad Térmica de Cristales Bidimensionales con Cadenas Extendidas de Puentes de Hidrógeno OH

Planteamiento del Problema
El desarrollo de membranas de intercambio de protones (PEMs) capaces de operar bajo temperaturas elevadas y condiciones anhidras sigue siendo un desafío crítico. Las PEMs actuales, que a menudo dependen del ácido fosfórico, sufren una degradación del rendimiento a altas temperaturas y baja humedad. Un cuello de botella principal es la falta de sistemas moleculares que puedan formar cadenas de puentes de hidrógeno de alta temperatura (específicamente cadenas O–H···O–H···O–H) necesarias para un transporte de protones eficiente. Mientras que los sistemas biológicos utilizan tales cadenas para la conducción de protones, el diseño de sistemas moleculares sintéticos que mantengan estas redes extendidas sobre un sustrato a altas temperaturas requiere un equilibrio entre una fuerte adsorción al sustrato y la formación de redes de puentes de hidrógeno continuas.

Metodología
El autor empleó simulaciones de dinámica molecular (MD) para evaluar la viabilidad de fabricar cristales 2D estables a partir de moléculas planas adsorbidas en una lámina de nitruro de boro hexagonal (h-BN). El estudio utilizó una aproximación de átomo unido, donde los grupos CH y CH₃ fueron tratados como sitios de interacción únicos.

• Campos de Fuerza: Las simulaciones utilizaron un campo de fuerza derivado del AMBER General Force Field (v2.1). Las interacciones intramoleculares se modelaron utilizando potenciales armónicos para el estiramiento de enlaces y deformaciones de ángulos de valencia, y potenciales periódicos para ángulos diedros. Las interacciones no enlazantes se describieron mediante potenciales de Lennard-Jones (LJ) e interacciones coulómbicas.
• Interacción con el Sustrato: La interacción entre las moléculas y el sustrato de h-BN se modeló utilizando un potencial de Lennard-Jones (m, l) con exponentes $l=10$ y $m=4.25$, representando un plano atractivo fijo.
• Protocolo de Simulación: El estudio minimizó la energía potencial para encontrar configuraciones de estado fundamental para monocapas de diversas moléculas. Para evaluar la estabilidad térmica, los sistemas se acoplaron a un termostato de Langevin. Las simulaciones rastrearon la capacidad calorífica adimensional ($c$), el número promedio de puentes de hidrógeno por molécula ($p_{hb}$) y la fracción de moléculas adsorbidas de forma plana sobre el sustrato ($p_s$) a medida que aumentaba la temperatura. Se modelaron tanto cristalitos finitos (con bordes libres) como cristales 2D infinitos (con condiciones de contorno periódicas) para observar las transiciones de fusión.

Moléculas Investigadas
El estudio se centró en moléculas planas que contienen anillos de benceno (para asegurar una fuerte adsorción al sustrato) y grupos hidroxilo (OH) o peptídicos (para formar puentes de hidrógeno). Las moléculas específicas analizadas fueron:

1. Fenol (C₆H₅OH)
2. Hidroquinona (C₆H₄(OH)₂)
3. 4-Fenilfenol (C₆H₅–C₆H₄OH)
4. 4-(4-Fenilfenil)fenol (C₆H₅–C₆H₄–C₆H₄OH)
5. Paracetamol (CH₃C(O)NHC₆H₄OH)
6. 4-Hidroxibenzanilida (C₆H₅C(O)NHC₆H₄OH)
7. 4,4-Dihidroxibenzanilida (C₆H₄OHC(O)NHC₆H₄OH)

Resultos Clave
Las simulaciones demostraron que estas moléculas pueden formar cristales bidimensionales estables caracterizados por cadenas lineales de puentes de hidrógeno.

• Formación Estructural: Moléculas como el paracetamol, las hidroxibenzanilidas y los fenoles formaron estructuras periódicas con cadenas moleculares paralelas. En la estructura de paracetamol homoquiral, se observaron cadenas lineales de grupos hidroxilo (Tipo 1) y grupos peptídicos (Tipo 12). En contraste, las estructuras racémicas a menudo formaron cadenas mixtas en zigzag, las cuales el autor señala como no aptas para la conducción de protones debido a que carecen de las cadenas continuas de OH Tipo 1.
• Estabilidad Térmica: Los cristales 2D exhibieron una alta estabilidad térmica, con la fusión ocurriendo a través de un rango de temperatura continuo en lugar de un punto único y agudo. Las temperaturas de inicio de fusión ($T_1$ para cristalitos finitos y $T_3$ para cristales infinitos) fueron significativamente más altas que las de sus contrapartes masivas en algunos casos, o comparables a ellas.
  • Fenol: Estable hasta ~320 K (cristal infinito).
  • Hidroquinona: Estable hasta ~460 K.
  • 4-Fenilfenol: Estable hasta ~400 K.
  • 4-(4-Fenilfenil)fenol: Estable hasta ~540 K.
  • Paracetamol: Estable hasta ~440 K (cristal infinito).
  • 4-Hidroxibenzanilida: Estable hasta ~580 K.
  • 4,4-Dihidroxibenzanilida: Estable hasta ~650 K.
• Correlación con Datos Masivos (Bulk): Las temperaturas de inicio de fusión simuladas para los cristales 2D mostraron un excelente acuerdo con las temperaturas de fusión experimentales de los correspondientes cristales 3D masivos (por ejemplo, el $T_0$ del paracetamol masivo $\approx 445$ K frente a la estabilidad del cristal 2D hasta 440 K).
• Mecanismo de Estabilidad: Se encontró que los anillos de benceno eran cruciales para una fuerte interacción de van der Waals con el sustrato plano de h-BN, evitando la desorción, mientras que la disposición específica de los grupos OH y peptídicos permitió la formación de redes extendidas de puentes de hidrógeno sin impedimento estérico.

Significancia y Afirmaciones
El artículo concluye que las moléculas planas con anillos de benceno y grupos hidroxilo/peptídicos pueden formar cristales 2D sobre láminas de h-BN que mantienen cadenas de puentes de hidrógeno lineales a temperaturas sustancialmente elevadas. El autor destaca que las temperaturas de inicio de fusión para estos cristales 2D oscilan entre 47 °C y 377 °C, superando con creces el punto de fusión del ácido fosfórico (42 °C), un electrolito común en las PEMs.

La principal significancia afirmada es que estas estructuras ofrecen una vía potencial para el desarrollo de nuevas membranas de intercambio de protones anhidras capaces de operar a altas temperaturas. Específicamente, el autor sugiere que las estructuras multicapa compuestas por láminas de h-BN y moléculas de hidroquinona, paracetamol o 4-hidroxibenzanilida representan plataformas prometedoras para tales aplicaciones. El estudio postula que la combinación de una fuerte adsorción al sustrato y redes extendidas de puentes de hidrógeno proporciona la estabilidad térmica necesaria para la conducción de protones en condiciones extremas, abordando una limitación clave en la tecnología actual de PEM.