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

Resumo Técnico: Mecanismo Microscópico da Amorfização Térmica de ZIF-4 e Fusão de ZIF-zni

Definição do Problema
Estruturas Metal-Orgânicas Amorfas (MOFs) oferecem uma combinação única de porosidade intrínseca e robustez mecânica, mas os mecanismos microscópicos que regem sua formação a partir de precursores cristalinos permanecem mal compreendidos. Especificamente, as transições de fase ordenado–desordenado em Estruturas Zeolíticas de Imidazolato (ZIFs), como a amorfização térmica de ZIF-4 e a fusão de ZIF-zni, apresentam desafios significativos para a modelagem computacional. Embora existam dados experimentais, questões relativas aos mecanismos específicos de amorfização, à natureza dos estados intermediários e ao comportamento além das distâncias de primeiro vizinho permanecem sem resposta. As dinâmicas moleculares ab initio (AIMD) padrão são frequentemente limitadas pelo tamanho do sistema e pelas escalas de tempo, enquanto campos de força tradicionais frequentemente falham em capturar a reatividade metal–ligante necessária para modelar a quebra e a formação de ligações durante essas transições.

Metodologia
Os autores empregaram uma abordagem híbrida combinando dinâmica molecular (MD) clássica com aprendizado de máquina (ML) para investigar as transições de fase de dois polimorfos de Zn(C3H3N2)2: o poroso ZIF-4 e o denso ZIF-zni.

1. Campo de Força: As simulações foram conduzidas utilizando o campo de força nb-ZIF-FF (non-bonded ZIF-FF). Este campo de força reativo utiliza potenciais de Morse para modelar ligações metal–ligante, permitindo a quebra e a formação de ligações. Ele incorpora um modelo de "átomo fantasma" (dummy atom) para impor a topologia tetraédrica correta ao redor dos centros de Zn²⁺, garantindo a reprodução precisa das propriedades estruturais, mecânicas e termodinâmicas para as fases cristalina, amorfa (ZIF_a) e líquida (ZIF_liq).
2. Geração de Estrutura: Estruturas amorfas (ZIF_a) foram geradas via recozimento simulado (aquecimento a 900 K e resfriamento) a partir de precursores de ZIF-4 e ZIF-zni. As estruturas resultantes foram validadas em relação à densidade experimental, módulo de bulk e funções de distribuição radial (RDFs).
3. Classificação por Aprendizado de Máquina: Para distinguir entre as quatro fases (ZIF-4, ZIF-zni, ZIF_a, ZIF_liq) no nível atômico local, os autores treinaram uma rede neural.
   • Descritores: Funções de Simetria de Behler-Parrinello (BPSF) foram computadas para ambientes centrados em Zn²⁺, servindo como entradas não enviesadas.
   • Arquitetura: Uma rede neural com 12 nós de entrada (funções de simetria), uma camada oculta (6 nós) e 4 nós de saída (probabilidades de fase).
   • Desempenho: O modelo alcançou 90,3% de precisão na classificação de ambientes locais, distinguindo com sucesso entre os dois estados desordenados (ZIF_a e ZIF_liq) baseando-se apenas em informações estruturais.
4. Protocolo de Simulação: Simulações de MD NPT de larga escala foram realizadas usando o LAMMPS. As temperaturas de transição foram determinadas identificando condições de equilíbrio onde as fases coexistem sem o crescimento líquido de uma sobre a outra.

Principais Resultados

• Termodinâmica: O estudo determinou temperaturas de transição de equilíbrio de 420 K para a amorfização de ZIF-4 e 625 K para a fusão de ZIF-zni. Embora esses valores difiram das observações experimentais (589 K e 863 K, respectivamente), o modelo reproduz corretamente as tendências de estabilidade termodinâmica, confirmando que o ZIF-4 é metaestável em condições ambientes. Entalpias e entropias calculadas alinham-se bem com os valores de referência.
• Mecanismo de Amorfização de ZIF-4: A transição do cristal ZIF-4 para o amorfo ZIF_a é um processo de duas etapas:
  1. Formação de uma fase intermediária de baixa densidade, do tipo líquido (ZIF_liq-like).
  2. Transformação desta fase intermediária na fase amorfa de alta densidade (ZIF_a).
     O processo envolve uma diminuição na conectividade seguida por uma reorganização. A formação da fase amorfa é anisotrópica, propagando-se mais rapidamente ao longo do eixo x. Esta anisotropia correlaciona-se com as distribuições de densidade local e a labilidade das ligações de coordenação; as ligações orientadas na direção x são mais propensas a quebrar, facilitando a troca de conectividade.
• Mecanismo de Fusão de ZIF-zni: A fusão do ZIF-zni também segue um caminho de duas etapas, mas em ordem inversa de evolução de densidade em comparação ao ZIF-4:
  1. Formação de uma fase do tipo amorfa a partir do cristal.
  2. Transição para a fase líquida (ZIF_liq).
     A fase líquida nucleia dentro do volume da fase desordenada intermediária.
• Gatilhos Microscópicos: A iniciação da amorfização está ligada a mudanças específicas de conectividade, e não apenas à simples quebra de ligações. Defeitos de Zn tricoordenados ocorrem frequentemente, mas são reversíveis. O evento crítico que desencadeia a transição de fase é a formação de uma nova ligação entre um centro de Zn²⁺ e uma molécula de ligante diferente (mudando a conectividade da rede), em vez de uma simples rotação ou reformação da ligação original. Este evento é acompanhado por um aumento súbito na população de ambientes do tipo líquido identificados pela rede neural.
• Evolução Estrutural: A formação de anéis de 3 membros (Zn–ligante–Zn) é observada como correlacionada ao crescimento da fase desordenada, sugerindo que estas características topológicas são características do estado amorfo, enquanto o cristal contém primariamente anéis de 4, 6 e 8 membros.

Significância e Alegações
O artigo afirma fornecer uma compreensão microscópica detalhada da amorfização térmica e da fusão em ZIFs, abordando lacunas deixadas por estudos anteriores que ignoraram estados intermediários ou careciam de tamanho de sistema suficiente. As principais contribuições são:

1. Validação do nb-ZIF-FF: Demonstrar que este campo de força pode reproduzir com precisão as estruturas e propriedades das fases cristalina, amorfa e líquida de ZIF, incluindo distribuições angulares sutis que outros campos de força reativos (como o ReaxFF) não capturam.
2. Dinâmica Aumentada por ML: Estabelecer uma metodologia onde uma rede neural, treinada em funções de simetria, pode rastrear a evolução temporal de transições de fase com alta resolução espacial, distinguindo entre fases desordenadas estruturalmente semelhantes.
3. Insight Mecanístico: Identificar que tanto a amorfização quanto a fusão são processos de duas etapas envolvendo fases intermediárias com densidades distintas e que a amorfização de ZIF-4 é um processo anisotrópico impulsionado pela labilidade de ligações de coordenação específicas.

Os autores concluem que seu trabalho oferece um arcabouço para estudar processos reativos em MOFs, podendo auxiliar no design de MOFs amorfas para aplicações industriais, como purificação de água e catálise, embora não proponham novos protocolos experimentais ou aplicações futuras específicas além do escopo da compreensão dos mecanismos fundamentais.