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

Cet article présente nb-UiO-FF, un nouveau champ de forces partiellement réactif qui modélise avec précision les propriétés structurelles, mécaniques et de défauts du réseau organométallique UiO-66 et permet des simulations de dynamique moléculaire de sa synthèse solvothermale et de ses mécanismes d'auto-assemblage.

L'essentiel

Imaginez un cadre métal-organique (MOF) comme le UiO-66 comme un gigantesque puzzle tridimensionnel microscopique. Il est construit à partir de « nœuds » métalliques (comme des clusters de zirconium) agissant comme les coins, et de « connecteurs » organiques (comme des cycles benzéniques) agissant comme les tiges de liaison. Les scientifiques adorent ces puzzles car ils sont incroyablement robustes et peuvent être ajustés pour piéger des gaz, délivrer des médicaments ou accélérer des réactions chimiques.

Cependant, il y a un problème : comprendre exactement comment ces puzzles s'assemblent au départ, c'est comme essayer de regarder un film joué à la vitesse de la lumière. Les liaisons chimiques qui maintiennent le métal et les connecteurs ensemble sont difficiles à simuler sur un ordinateur. La plupart des modèles informatiques traitent ces connexions comme une colle permanente ; ils ne peuvent pas montrer l'application de la colle, l'assemblage des pièces, ni même ce qui se passe si une pièce manque.

La Solution : Un champ de force « Colle Intelligente »
Dans cet article, les auteurs présentent un nouvel outil informatique appelé nb-UiO-FF. Imaginez cela comme un nouveau jeu de règles pour un jeu de simulation qui permet aux pièces du puzzle d'être « partiellement réactives ».

Voici comment ils l'ont rendu fonctionnel, en utilisant quelques analogies simples :

• Les Atomes « Fictifs » (Les Mains Invisibles) :
  Dans le monde réel, le nœud métallique de zirconium possède une charge électrique complexe qui attire les connecteurs dans des directions spécifiques. Les modèles informatiques standards peinent à imiter cela sans devenir désordonnés. Les auteurs ont résolu ce problème en attachant des atomes « fictifs » invisibles (comme de petits dummies magnétiques) aux nœuds métalliques. Ces dummies agissent comme des mains invisibles qui maintiennent les connecteurs dans la bonne forme et l'orientation, assurant que le puzzle se construit correctement sans nécessiter de calculs complexes et lourds.

• Le « Potentiel de Morse » (Le Ressort Extensible) :
  Habituellement, les modèles informatiques traitent les liaisons comme des bâtons rigides. Si vous les tirez, ils cassent instantanément. Les auteurs ont remplacé ces bâtons rigides par un potentiel de Morse, qui agit davantage comme un ressort extensible. Cela permet à la simulation de montrer le métal et le connecteur s'étirer, osciller, et même s'assembler ou se désassembler dynamiquement. Ceci est crucial pour observer la « naissance » du matériau.

Ce qu'ils ont testé
Les auteurs n'ont pas seulement construit l'outil ; ils l'ont soumis à un test de stress rigoureux pour s'assurer de sa fiabilité :

1. Le Puzzle Parfait : Ils ont vérifié si l'outil pouvait recréer la forme exacte d'un cristal UiO-66 parfait. Il correspondait presque parfaitement aux mesures réelles (à une infime fraction de pourcentage près).
2. Le Puzzle Imbibé : Ils ont testé l'outil avec le cristal imbibé de deux liquides différents utilisés pour le fabriquer (DMF et éthanol). Le modèle a montré que le cristal reste solide et ne se désintègre pas lorsqu'il est humide.
3. Le Puzzle Cassé : Les cristaux réels ont souvent des pièces manquantes (défauts). Les auteurs ont intentionnellement retiré des connecteurs ou des nœuds entiers dans la simulation. L'outil a montré avec succès que le cristal pouvait toujours maintenir sa forme même avec ces trous, tout comme le matériau réel.
4. Le Puzzle Rebondissant : Ils ont testé à quel point on pouvait comprimer le cristal avant qu'il ne se déforme. Les résultats correspondaient à des calculs de physique de haut niveau, prouvant que le modèle comprend la résistance du matériau.
5. Le Puzzle Cousin : Ils ont essayé l'outil sur une version légèrement plus grande du puzzle (UiO-67) et cela a fonctionné là aussi, prouvant que les règles sont flexibles.

Observer la Magie se Produire
La partie la plus excitante de l'article est l'utilisation de cet nouvel outil pour observer le processus d'auto-assemblage. Imaginez déposer toutes les pièces du puzzle (nœuds métalliques et connecteurs) et le solvant liquide dans une boîte et appuyer sur « lecture ».

La simulation a montré les pièces dérivant et commençant lentement à coller ensemble.

• Ils ont vu les nœuds métalliques et les connecteurs se trouver et former les blocs de construction initiaux.
• Ils ont remarqué que parfois, les pièces se coincent dans des positions « incorrectes » (pièges cinétiques), comme une pièce de puzzle qui s'adapte lâchement mais n'est pas tout à fait juste.
• Ils ont observé que le processus est lent ; les pièces sont lourdes et se déplacent lentement, de sorte que le puzzle complet ne s'assemble pas entièrement dans le temps qu'ils ont simulé.

La Conclusion
Cet article présente un nouveau modèle informatique hautement précis qui agit comme un « microscope intelligent » pour le matériau UiO-66. Il peut simuler la structure du matériau, sa résistance et sa capacité à gérer les défauts. Plus important encore, c'est le premier outil de ce genre capable de simuler de manière réaliste le processus dynamique de la construction du matériau à partir de zéro, aidant les scientifiques à comprendre comment ces matériaux étonnants naissent et comment contrôler leurs imperfections.

Résumé technique

Résumé technique : Champ de force partiellement réactif pour le cadre organométallique UiO-66

Énoncé du problème
Les cadres organométalliques (MOF), en particulier UiO-66, sont largement étudiés pour leur modulabilité structurelle et leur stabilité. Cependant, la compréhension de leurs mécanismes de synthèse et de la formation de défauts reste difficile en raison du manque de modèles computationnels capables de simuler des processus réactifs. Les champs de force existants (par exemple, UFF4MOF, BTW-FF, MOF-FF) reposent généralement sur une connectivité prédéfinie, ce qui les rend incapables de modéliser la rupture et la formation de liaisons, essentielles pour étudier l'auto-assemblage des MOF et l'ingénierie des défauts. Bien que des champs de force réactifs comme Reax-FF et des potentiels d'apprentissage automatique existent, ils font face à des défis concernant leur précision pour les MOF, leur coût computationnel, ou la difficulté de les étendre pour inclure des solvants et des modulateurs. De plus, des modèles spécifiques à UiO-66 peinent souvent à représenter avec précision la forme hydroxylée ou nécessitent des configurations d'atomes fictifs coûteuses en calcul (par exemple, 48 atomes fictifs par nœud).

Méthodologie
Les auteurs ont développé nb-UiO-FF, un champ de force partiellement réactif conçu pour modéliser le cadre UiO-66 et son analogue isoréticulaire, UiO-67. La méthodologie s'appuie sur une approche précédemment établie pour ZIF-8, utilisant une combinaison d'interactions non liées et d'atomes fictifs cationiques (CDA).

• Construction du champ de force : Le modèle remplace les potentiels de liaison rigides entre le nœud métallique et les ligands par un potentiel de Morse non lié afin de permettre la rupture et la formation dynamiques de liaisons.
• Atomes fictifs (CDA) : Pour reproduire la distribution de charge anisotrope des atomes de Zr et maintenir une géométrie de coordination correcte sans termes de liaison, quatre atomes fictifs (D) sont introduits par centre Zr. Ceux-ci sont connectés à l'atome Zr via des liaisons et des angles harmoniques forts. La charge de +1,968 de l'atome Zr est entièrement distribuée parmi ces quatre atomes fictifs (+0,492 chacun), rendant l'atome Zr électriquement neutre.
• Paramétrisation : Les paramètres du champ de force ont été dérivés à partir d'un champ de force non réactif de Yang et al., ajusté pour le nœud hydroxylé Zr6O4(OH)4. Les paramètres liés au cluster métallique ont été tirés de modèles publiés, tandis que les longueurs de liaison et les angles d'équilibre proviennent de résultats DFT de Valenzano et al. Les paramètres du potentiel de Morse ($D_0$ et $\alpha$) pour l'interaction Zr–Ligand ont été optimisés en explorant systématiquement l'espace des paramètres afin d'équilibrer la précision de la densité numérique et des longueurs de liaison locales (spécifiquement la fonction de distribution radiale Zr–Ligand) par rapport aux références DFT et expérimentales.
• Protocoles de simulation : Toutes les simulations ont été effectuées à l'aide de LAMMPS. La validation a impliqué une analyse structurelle, des tests de stabilité dans des solvants (DMF et éthanol) et des calculs de propriétés mécaniques. Les études d'auto-assemblage ont utilisé un système composé de 16 clusters Zr, de 96 ligands BDC et de 2136 molécules de DMF, simulé sur 6 ns dans l'ensemble NPT.

Résultats clés

• Précision structurelle : Le nb-UiO-FF reproduit le paramètre de maille expérimental de UiO-66 pur avec une déviation d'environ 0,046 % (20,71 Å contre 20,7004 Å). Les fonctions de distribution radiale (FDR) pour les paires d'atomes clés (Zr–$\mu_3$O, Zr–$\mu_3$OH, Zr–Zr, Zr–Ligand) montrent un accord raisonnable avec les données DFT et expérimentales, avec des déviations généralement inférieures à 10 %.
• Stabilité vis-à-vis des solvants et des défauts : Le champ de force maintient l'intégrité du cadre en présence de solvants DMF et éthanol, ainsi que dans des structures défectueuses contenant des ligands manquants et des nœuds manquants. Dans les modèles défectueux, la neutralité de charge a été préservée en coiffant les centres sous-coordonnés avec des ions chlorure ou des fragments BDC monodentés.
• Propriétés mécaniques : Les constantes élastiques calculées ($C_{11}$, $C_{12}$, $C_{44}$) et les propriétés mécaniques dérivées (module de compressibilité : 44 GPa, module de Young : 40 GPa) sont en bon accord avec les références DFT, bien que légèrement plus élevées, probablement en raison de la rigidité introduite par les fortes interactions liées entre Zr et les atomes fictifs.
• Transférabilité : Le modèle reproduit avec succès la structure du cadre isoréticulaire UiO-67, démontrant une transférabilité à travers différentes longueurs de ligands.
• Dynamique d'auto-assemblage : Dans les simulations dynamiques, le champ de force a capturé avec succès la formation progressive des liaisons de coordination entre les clusters métalliques et les ligands organiques. L'étude a identifié des motifs structuraux transitoires, notamment des unités de construction secondaires (SBU) initiales et des configurations de liaison non idéales (états piégés cinétiquement). Cependant, la convergence complète du nombre de coordination et la formation de grands agrégats ioniques n'ont pas été atteintes dans la fenêtre de simulation de 6 ns, suggérant la nécessité de méthodes d'échantillonnage amélioré pour des échelles de temps plus longues.

Signification et revendications
L'article revendique que nb-UiO-FF fournit une description cohérente et physiquement fiable de UiO-66 et UiO-67, comblant le fossé entre la modélisation structurelle statique et la dynamique réactive. En utilisant un ensemble réduit d'atomes fictifs (quatre par Zr) par rapport aux modèles précédents, les auteurs parviennent à un équilibre entre l'efficacité du modèle et la capacité à simuler des processus réactifs. Le champ de force permet l'étude des mécanismes d'auto-assemblage et de la formation de défauts, qui sont inaccessibles aux champs de force traditionnels à connectivité fixe. Les auteurs positionnent ce travail comme une base pour l'investigation des étapes précoces de la synthèse des MOF et de l'ingénierie des défauts ponctuels, tout en reconnaissant que des travaux supplémentaires avec des schémas d'échantillonnage amélioré (par exemple, la métadynamique) sont nécessaires pour surmonter pleinement les pièges cinétiques et modéliser l'auto-assemblage complet.