Partially reactive force field for the UiO-66 metal-organic framework

Este artigo apresenta o nb-UiO-FF, um novo campo de forças parcialmente reativo que modela com precisão as propriedades estruturais, mecânicas e de defeitos da estrutura metal-orgânica UiO-66 e permite simulações de dinâmica molecular de sua síntese solvotermal e mecanismos de auto montagem.

O essencial

Imagine uma estrutura metal-orgânica (MOF) como a UiO-66 como um quebra-cabeça tridimensional microscópico gigante. Ela é construída a partir de "nós" metálicos (como aglomerados de Zircônio) que atuam como os cantos, e "conectores" orgânicos (como anéis de benzeno) que atuam como as hastes de conexão. Os cientistas adoram esses quebra-cabeças porque são incrivelmente resistentes e podem ser ajustados para aprisionar gases, administrar medicamentos ou acelerar reações químicas.

No entanto, há um problema: descobrir exatamente como esses quebra-cabeças se encaixam desde o início é como tentar assistir a um filme reproduzido na velocidade da luz. As ligações químicas que mantêm o metal e os conectores unidos são difíceis de simular em um computador. A maioria dos modelos computacionais trata essas conexões como cola permanente; eles não conseguem mostrar a cola sendo aplicada, as peças se encaixando ou mesmo o que acontece se uma peça estiver faltando.

A Solução: Um Campo de Força de "Cola Inteligente"
Neste artigo, os autores apresentam uma nova ferramenta computacional chamada nb-UiO-FF. Pense nisso como um novo conjunto de regras para um jogo de simulação que permite que as peças do quebra-cabeça sejam "parcialmente reativas".

Veja como eles fizeram isso funcionar, usando algumas analogias simples:

• Os Átomos "Fictícios" (As Mãos Invisíveis):
  No mundo real, o nó metálico de Zircônio possui uma carga elétrica complexa que puxa os conectores em direções específicas. Modelos computacionais padrão lutam para imitar isso sem se tornar confusos. Os autores resolveram isso ao anexar átomos "fictícios" invisíveis (como pequenos manequins magnéticos) aos nós metálicos. Esses manequins atuam como mãos invisíveis que seguram os conectores na forma e orientação corretas, garantindo que o quebra-cabeça seja construído da maneira certa sem a necessidade de cálculos complexos e pesados.

• O "Potencial de Morse" (A Mola Esticável):
  Geralmente, os modelos computacionais tratam as ligações como hastes rígidas. Se você as puxar, elas se quebram instantaneamente. Os autores substituíram essas hastes rígidas por um potencial de Morse, que atua mais como uma mola esticável. Isso permite que a simulação mostre o metal e o conector esticando, oscilando e até se encaixando ou se separando dinamicamente. Isso é crucial para observar o "nascimento" do material.

O Que Eles Testaram
Os autores não apenas construíram a ferramenta; eles a submeteram a um teste de estresse rigoroso para garantir que fosse confiável:

1. O Quebra-Cabeça Perfeito: Eles verificaram se a ferramenta podia recriar a forma exata de um cristal perfeito de UiO-66. Ela combinou com medições do mundo real quase perfeitamente (dentro de uma fração minúscula de um por cento).
2. O Quebra-Cabeça Encharcado: Eles testaram a ferramenta com o cristal encharcado em dois líquidos diferentes usados para produzi-lo (DMF e etanol). O modelo mostrou que o cristal permanece forte e não se desintegra quando molhado.
3. O Quebra-Cabeça Quebrado: Cristais do mundo real frequentemente têm peças faltando (defeitos). Os autores removeram intencionalmente conectores ou nós inteiros na simulação. A ferramenta mostrou com sucesso que o cristal ainda podia manter sua forma mesmo com esses buracos, assim como o material real faz.
4. O Quebra-Cabeça Saltitante: Eles testaram o quão forte você pode apertar o cristal antes que ele se deforme. Os resultados combinaram com cálculos de física de alto nível, provando que o modelo entende a resistência do material.
5. O Quebra-Cabeça Primo: Eles tentaram usar a ferramenta em uma versão ligeiramente maior do quebra-cabeça (UiO-67) e ela funcionou lá também, provando que as regras são flexíveis.

Observando a Mágica Acontecer
A parte mais emocionante do artigo é usar essa nova ferramenta para observar o processo de auto-montagem. Imagine jogar todas as peças do quebra-cabeça (nós metálicos e conectores) e o solvente líquido em uma caixa e apertar "play".

A simulação mostrou as peças flutuando e começando lentamente a grudar umas nas outras.

• Eles viram os nós metálicos e os conectores encontrando-se e formando os blocos de construção iniciais.
• Eles notaram que, às vezes, as peças ficam presas em posições "erradas" (armadilhas cinéticas), como uma peça de quebra-cabeça que se encaixa frouxamente, mas não está totalmente certa.
• Eles observaram que o processo é lento; as peças são pesadas e se movem com dificuldade, então o quebra-cabeça completo não se monta totalmente no tempo que eles simularam.

A Conclusão
Este artigo apresenta um novo modelo computacional altamente preciso que atua como um "microscópio inteligente" para o material UiO-66. Ele pode simular a estrutura do material, sua resistência e sua capacidade de lidar com defeitos. Mais importante ainda, é a primeira ferramenta do seu tipo que pode simular realisticamente o processo dinâmico do material construindo-se a partir do zero, ajudando os cientistas a entender como esses materiais incríveis nascem e como controlar suas imperfeições.

Resumo técnico

Resumo Técnico: Campo de Força Parcialmente Reativo para a Estrutura Metal-Orgânica UiO-66

Declaração do Problema
As estruturas metal-orgânicas (MOFs), particularmente a UiO-66, são amplamente estudadas devido à sua sintonizabilidade estrutural e estabilidade. No entanto, compreender seus mecanismos de síntese e formação de defeitos permanece desafiador devido à falta de modelos computacionais capazes de simular processos reativos. Os campos de força existentes (por exemplo, UFF4MOF, BTW-FF, MOF-FF) geralmente dependem de conectividade pré-definida, tornando-os incapazes de modelar a quebra e formação de ligações, o que é essencial para estudar a auto-organização de MOFs e a engenharia de defeitos. Embora existam campos de força reativos como o Reax-FF e potenciais de aprendizado de máquina, eles enfrentam desafios relacionados à precisão para MOFs, custo computacional ou à dificuldade de expandi-los para incluir solventes e moduladores. Além disso, modelos específicos para a UiO-66 frequentemente lutam para representar com precisão a forma hidroxilada ou exigem configurações computacionalmente caras de átomos fictícios (por exemplo, 48 átomos fictícios por nó).

Metodologia
Os autores desenvolveram o nb-UiO-FF, um campo de força parcialmente reativo projetado para modelar a estrutura UiO-66 e seu análogo isoreticular, a UiO-67. A metodologia baseia-se em uma abordagem previamente estabelecida para a ZIF-8, utilizando uma combinação de interações não ligadas e átomos fictícios catiônicos (CDA).

• Construção do Campo de Força: O modelo substitui potenciais rígidos ligados entre o nó metálico e os ligantes por um potencial de Morse não ligado para permitir a quebra e formação dinâmica de ligações.
• Átomos Fictícios (CDA): Para reproduzir a distribuição anisotrópica de carga dos átomos de Zr e manter a geometria de coordenação correta sem termos ligados, quatro átomos fictícios (D) são introduzidos por centro de Zr. Estes são conectados ao átomo de Zr por meio de ligações e ângulos harmônicos fortes. A carga de +1,968 do átomo de Zr é totalmente distribuída entre esses quatro átomos fictícios (+0,492 cada), tornando o átomo de Zr eletricamente neutro.
• Parametrização: Os parâmetros do campo de força foram derivados a partir de um campo de força não reativo de Yang et al., ajustado para o nó hidroxilado Zr6O4(OH)4. Os parâmetros ligados para o cluster metálico foram retirados de modelos relatados, enquanto os comprimentos de ligação e ângulos de equilíbrio foram obtidos a partir de resultados de DFT de Valenzano et al. Os parâmetros do potencial de Morse ($D_0$ e $\alpha$) para a interação Zr–Oligante foram otimizados explorando sistematicamente o espaço de parâmetros para equilibrar a precisão na densidade numérica e nos comprimentos de ligação locais (especificamente a função de distribuição radial Zr–Oligante) contra referências de DFT e experimentais.
• Protocolos de Simulação: Todas as simulações foram realizadas usando LAMMPS. A validação envolveu análise estrutural, testes de estabilidade em solventes (DMF e etanol) e cálculos de propriedades mecânicas. Estudos de auto-organização utilizaram um sistema de 16 clusters de Zr, 96 ligantes BDC e 2136 moléculas de DMF, simulados por 6 ns no ensemble NPT.

Principais Resultados

• Precisão Estrutural: O nb-UiO-FF reproduz o parâmetro de rede experimental da UiO-66 pura com um desvio de aproximadamente 0,046% (20,71 Å vs. 20,7004 Å). As funções de distribuição radial (FDRs) para pares atômicos-chave (Zr–$\mu_3$O, Zr–$\mu_3$OH, Zr–Zr, Zr–Oligante) mostram concordância razoável com dados de DFT e experimentais, com desvios geralmente dentro de 10%.
• Estabilidade de Solvente e Defeitos: O campo de força mantém a integridade da estrutura na presença de solventes DMF e etanol, bem como em estruturas defeituosas contendo ligantes e nós ausentes. Em modelos defeituosos, a neutralidade de carga foi preservada ao tampar centros subcoordenados com íons cloreto ou fragmentos BDC monodentados.
• Propriedades Mecânicas: As constantes elásticas calculadas ($C_{11}$, $C_{12}$, $C_{44}$) e as propriedades mecânicas derivadas (Módulo de Bulk: 44 GPa, Módulo de Young: 40 GPa) estão em boa concordância com os benchmarks de DFT, embora ligeiramente mais altas, provavelmente devido à rigidez introduzida pelas fortes interações ligadas entre Zr e átomos fictícios.
• Transferibilidade: O modelo reproduz com sucesso a estrutura da estrutura isoreticular UiO-67, demonstrando transferibilidade através de diferentes comprimentos de ligante.
• Dinâmica de Auto-organização: Em simulações dinâmicas, o campo de força capturou com sucesso a formação progressiva de ligações de coordenação entre clusters metálicos e ligantes orgânicos. O estudo identificou motivos estruturais transitórios, incluindo unidades construtoras secundárias (SBUs) iniciais e configurações de ligação não ideais (estados presos cineticamente). No entanto, a convergência completa do número de coordenação e a formação de grandes agregados iônicos não foram alcançadas dentro da janela de simulação de 6 ns, sugerindo a necessidade de métodos de amostragem aprimorada para escalas de tempo mais longas.

Significado e Alegações
O artigo alega que o nb-UiO-FF fornece uma descrição consistente e fisicamente confiável da UiO-66 e da UiO-67, preenchendo a lacuna entre modelagem estrutural estática e dinâmica reativa. Ao utilizar um conjunto reduzido de átomos fictícios (quatro por Zr) em comparação com modelos anteriores, os autores alcançam um equilíbrio entre eficiência do modelo e a capacidade de simular processos reativos. O campo de força permite o estudo de mecanismos de auto-organização e formação de defeitos, que são inacessíveis a campos de força tradicionais de conectividade fixa. Os autores posicionam este trabalho como uma fundação para investigar os estágios iniciais da síntese de MOFs e a engenharia de defeitos pontuais, ao mesmo tempo em que reconhecem que trabalhos adicionais com esquemas de amostragem aprimorada (por exemplo, metadinâmica) são necessários para superar completamente as armadilhas cinéticas e modelar a auto-organização completa.