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

Resumo Técnico: Estabilidade Térmica de Cristais Bidimensionais com Cadeias Estendidas de Pontes de Hidrogênio OH

Definição do Problema
O desenvolvimento de membranas de troca de prótons (PEMs) capazes de operar sob temperaturas elevadas e condições anidras permanece um desafio crítico. As PEMs atuais, frequentemente dependentes de ácido fosfórico, sofrem deterioração de desempenho em altas temperaturas e baixa umidade. Um gargalo primário é a falta de sistemas moleculares que possam formar cadeias de pontes de hidrogênio termicamente estáveis em alta temperatura (especificamente cadeias O–H···O–H···O–H) necessárias para o transporte eficiente de prótons. Embora sistemas biológicos utilizem tais cadeias para condução de prótons, projetar sistemas moleculares sintéticos que mantenham essas redes estendidas sobre um substrato em altas temperaturas requer um equilíbrio entre uma forte adsorção ao substrato e a formação de redes contínuas de pontes de hidrogênio.

Metodologia
O autor empregou simulações de dinâmica molecular (MD) para avaliar a viabilidade de fabricar cristais 2D estáveis a partir de moléculas planas adsorvidas em uma folha de nitreto de boro hexagonal (h-BN). O estudo utilizou uma aproximação de átomo unido, onde os grupos CH e CH₃ foram tratados como sítios de interação únicos.

• Campos de Força: As simulações utilizaram um campo de força derivado do AMBER General Force Field (v2.1). As interações intramoleculares foram modeladas usando potenciais harmônicos para deformações de estiramento de ligação e ângulos de valência, e potenciais periódicos para ângulos diedros. As interações não ligadas foram descritas por potenciais de Lennard-Jones (LJ) e interações coulombicas.
• Interação com o Substrato: A interação entre as moléculas e o substrato de h-BN foi modelada usando um potencial Lennard-Jones (m, l) com expoentes $l=10$ e $m=4,25$, representando um plano atrativo fixo.
• Protocolo de Simulação: O estudo minimizou a energia potencial para encontrar configurações de estado fundamental para monocamadas de várias moléculas. Para avaliar a estabilidade térmica, os sistemas foram acoplados a um termostato de Langevin. As simulações rastrearam a capacidade térmica adimensional ($c$), o número médio de pontes de hidrogênio por molécula ($p_{hb}$) e a fração de moléculas adsorvidas de forma plana no substrato ($p_s$) conforme a temperatura aumentava. Tanto cristalitos finitos (com bordas livres) quanto cristais 2D infinitos (com condições de contorno periódicas) foram modelados para observar transições de fusão.

Moléculas Investigadas
O estudo focou em moléculas planas contendo anéis de benzeno (para garantir uma forte adsorção ao substrato) e grupos hidroxila (OH) ou peptídeos (para formar pontes de hidrogênio). As moléculas específicas analisadas foram:

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)

Principais Resultados
As simulações demonstraram que estas moléculas podem formar cristais bidimensionais estáveis caracterizados por cadeias lineares de pontes de hidrogênio.

• Formação Estrutural: Moléculas como o paracetamol, hidroxibenzanilidas e fenóis formaram estruturas periódicas com cadeias moleculares paralelas. Na estrutura homoquiral do paracetamol, observaram-se cadeias lineares de grupos hidroxila (Tipo 1) e grupos peptídeos (Tipo 12). Em contraste, estruturas racêmicas frequentemente formaram cadeias em zigue-zague mistas, que o autor observa serem inadequadas para a condução de prótons por carecerem das cadeias contínuas de OH do Tipo 1.
• Estabilidade Térmica: Os cristais 2D exibiram alta estabilidade térmica, com a fusão ocorrendo ao longo de uma faixa de temperatura contínua, em vez de em um único ponto agudo. As temperaturas de início de fusão ($T_1$ para cristalitos finitos e $T_3$ para cristais infinitos) foram significativamente superiores às de seus equivalentes em massa (bulk) em alguns casos, ou comparáveis a eles.
  • Fenol: Estável até ~320 K (cristal infinito).
  • Hidroquinona: Estável até ~460 K.
  • 4-Fenilfenol: Estável até ~400 K.
  • 4-(4-Fenilfenil)fenol: Estável até ~540 K.
  • Paracetamol: Estável até ~440 K (cristal infinito).
  • 4-Hidroxibenzanilida: Estável até ~580 K.
  • 4,4-Dihidroxibenzanilida: Estável até ~650 K.
• Correlação com Dados de Massa (Bulk): As temperaturas de início de fusão simuladas para os cristais 2D mostraram excelente concordância com as temperaturas de fusão experimentais dos correspondentes cristais 3D em massa (ex: $T_0 \approx 445$ K do paracetamol em massa vs. estabilidade do cristal 2D até 440 K).
• Mecanismo de Estabilidade: Os anéis de benzeno foram considerados cruciais para a forte interação de van der Waals com o substrato plano de h-BN, impedindo a dessorção, enquanto o arranjo específico dos grupos OH e peptídeos permitiu a formação de redes estendidas de pontes de hidrogênio sem impedimento estérico.

Significância e Alegações
O artigo conclui que moléculas planas com anéis de benzeno e grupos hidroxila/peptídeos podem formar cristais 2D em folhas de h-BN que mantêm cadeias de pontes de hidrogênio lineares em temperaturas substancialmente elevadas. O autor destaca que as temperaturas de início de fusão para estes cristais 2D variam de 47 °C a 377 °C, excedendo em muito o ponto de fusão do ácido fosfórico (42 °C), um eletrólito comum em PEMs.

A principal significância alegada é que estas estruturas oferecem um caminho potencial para o desenvolvimento de novas membranas de troca de prótons anidras capazes de operar em altas temperaturas. Especificamente, o autor sugere que estruturas multicamadas compostas por folhas de h-BN e moléculas de hidroquinona, paracetamol ou 4-hidroxibenzanilida representam plataformas promissoras para tais aplicações. O estudo postula que a combinação de forte adsorção ao substrato e redes estendidas de pontes de hidrogênio fornece a estabilidade térmica necessária para a condução de prótons em condições extremas, abordando uma limitação fundamental na tecnologia atual de PEM.