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

本論文は、UiO-66 金属有機構造体の構造、機械的、および欠陥特性を正確にモデル化し、その溶媒熱合成および自己集合メカニズムの分子動力学シミュレーションを可能にする新規の部分的に反応性のある力場 nb-UiO-FF を導入する。

要約

金属有機構造体（MOF）であるUiO-66を、巨大で微細な 3 次元パズルだと想像してみてください。これは、金属の「ノード」（ジルコニウムクラスターなど）が角となり、有機の「リンカー」（ベンゼン環など）が接続棒となるように組み立てられています。科学者たちは、これらのパズルが非常に頑丈であり、ガスを閉じ込めたり、薬物を送達したり、化学反応を加速させたりするように調整できるため、これらを愛用しています。

しかし、問題があります。これらのパズルが最初にどのように組み合わさるのかを正確に理解することは、光速で再生される映画を観ようとするようなものです。金属とリンカーを結びつけている化学結合は、コンピュータ上でシミュレーションするのが難しいのです。ほとんどのコンピュータモデルは、これらの結合を永久的な接着剤として扱います。接着剤が塗布される様子や、部品がカチッと嵌まる様子、あるいは部品が欠けた場合に何が起こるかを示すことはできません。

解決策：「スマートな接着剤」力場
この論文では、著者らはnb-UiO-FFと呼ばれる新しいコンピュータツールを紹介しています。これは、パズルの部品が「部分的に反応性を持つ」ことを可能にする、シミュレーションゲームの新しいルールセットだと考えてください。

彼らがどのように機能させたか、いくつかの単純な比喩を用いて説明します。

• 「ダミー」原子（見えない手）：
  現実世界では、ジルコニウム金属ノードは複雑な電荷を持っており、リンカーを特定の方向に引き寄せます。標準的なコンピュータモデルは、これを messy（ごちゃごちゃした）状態にならずに模倣するのが困難です。著者らは、金属ノードに目に見えない「ダミー」原子（小さな磁気ダミーのようなもの）を結合させることでこれを解決しました。これらのダミーは、リンカーを正しい形状と向きに保持する「見えない手」として機能し、複雑で重たい計算を必要とせずに、パズルが正しく組み立てられるようにします。

• 「モーゼポテンシャル」（伸縮するバネ）：
  通常、コンピュータモデルは結合を剛直な棒として扱います。引っ張れば、瞬時に壊れてしまいます。著者らは、これらの剛直な棒を、より伸縮するバネのように働くモーゼポテンシャルに置き換えました。これにより、シミュレーションは金属とリンカーが伸びたり、揺らめいたり、動的にカチッと嵌まったり、バラバラになったりする様子を示すことができます。これは、物質の「誕生」を観察する上で不可欠です。

彼らがテストしたこと
著者らは単にツールを構築しただけでなく、それが信頼できることを確認するために厳格なストレステストを行いました。

1. 完璧なパズル： ツールが完璧な UiO-66 結晶の正確な形状を再現できるか確認しました。それは実世界の測定値とほぼ完璧に一致しました（ごくわずかなパーセントの範囲内）。
2. 水に浸ったパズル： ツールを、製造に使用される 2 種類の液体（DMF とエタノール）に浸した結晶でテストしました。モデルは、結晶が濡れても強く保たれ、崩壊しないことを示しました。
3. 壊れたパズル： 現実世界の結晶には、欠けた部品（欠陥）がよくあります。著者らはシミュレーションで意図的にリンカーやノード全体を除去しました。ツールは、これらの穴があいていても、実際の物質がそうであるように、結晶が形状を保ち続けることを正常に示しました。
4. 跳ねるパズル： 結晶が変形するまで、どれほど強く押しつぶせるかをテストしました。結果は高レベルの物理計算と一致し、モデルが物質の強度を理解していることを証明しました。
5. いとこのパズル： 彼らはツールを、パズルの少し大きなバージョン（UiO-67）で試しましたが、そこでも機能し、ルールが柔軟であることを証明しました。

魔法が起こるのを見る
この論文の最もエキサイティングな部分は、この新しいツールを使って自己組織化過程を観察することです。すべてのパズルの部品（金属ノードとリンカー）と液体溶媒を箱に落とし、「再生」ボタンを押すことを想像してください。

シミュレーションは、部品が漂い、ゆっくりと貼り付き始める様子を示しました。

• 彼らは、金属ノードとリンカーが互いを見つけ、初期の構築ブロックを形成するのを目撃しました。
• 彼らは、部品が時折「間違った」位置に引っかかる（運動論的トラップ）ことに気づきました。これは、パズルの部品が緩くはまるが、完全に正しいわけではないような状態です。
• 彼らは、このプロセスが遅いことを観察しました。部品は重く、鈍重に動くため、彼らがシミュレーションした時間内では、完全なパズルは完全に組み立てられませんでした。

結論
この論文は、UiO-66 物質のための「スマートな顕微鏡」として機能する、新しい高精度のコンピュータモデルを提示しています。これは、物質の構造、強度、欠陥への耐性をシミュレートできます。最も重要なのは、これが物質がゼロから自ら組み立てられる動的過程を現実的にシミュレートできる、同種の最初のツールであることです。これにより、科学者たちはこれらの驚くべき物質がどのように生まれるか、そしてその欠陥をどのように制御するかを理解する手助けをされています。

技術概要

技術的概要：UiO-66 金属有機構造体に対する部分的反応性力場

問題提起
金属有機構造体（MOF）、特に UiO-66 は、その構造の調整可能性と安定性から広く研究されています。しかし、反応過程をシミュレートできる計算モデルの欠如により、その合成メカニズムや欠陥形成の理解は依然として困難です。既存の力場（UFF4MOF、BTW-FF、MOF-FF など）は通常、事前に定義された結合性に依存しており、MOF の自己集合や欠陥工学の研究に不可欠な結合の切断と形成をモデル化することができません。Reax-FF や機械学習ポテンシャルなどの反応性力場は存在しますが、MOF に対する精度、計算コスト、あるいは溶媒やモジュレーターを含めることの難しさといった課題に直面しています。さらに、UiO-66 固有のモデルは、水酸化された形態を正確に表現する際に苦労するか、計算的に高価なダミー原子配置（例：ノードあたり 48 個のダミー原子）を必要とすることが多いです。

手法
著者らは、UiO-66 骨格とその同系物である UiO-67 をモデル化するために設計された部分的反応性力場「nb-UiO-FF」を開発しました。この手法は、以前 ZIF-8 に対して確立されたアプローチに基づき、非結合相互作用と陽イオン性ダミー原子（CDA）の組み合わせを利用しています。

• 力場の構築: 金属ノードとリガンド間の硬い結合ポテンシャルを、動的な結合の切断と形成を可能にする非結合モーセポテンシャルに置き換えました。
• ダミー原子（CDA）: 結合項なしで Zr 原子の異方性電荷分布を再現し、正しい配位幾何構造を維持するために、Zr 中心あたり 4 つのダミー原子（D）を導入しました。これらは強い調和結合と角度を介して Zr 原子に結合されています。Zr 原子の +1.968 の電荷は、これら 4 つのダミー原子（それぞれ +0.492）に完全に分配され、Zr 原子は電気的に中性となります。
• パラメータ化: 力場パラメータは、Yang らによる非反応性力場から導出され、水酸化された Zr6O4(OH)4 ノードに合わせて調整されました。金属クラスターの結合パラメータは報告されたモデルから採取され、平衡結合長と角度は Valenzano らによる DFT 結果から引用されました。Zr-リガンド相互作用のモーセポテンシャルパラメータ（$D_0$および$\alpha$）は、DFT および実験参照値に対する数密度と局所結合長（特に Zr-リガンドの動径分布関数）の精度をバランスさせるために、パラメータ空間を体系的に探索して最適化されました。
• シミュレーションプロトコル: すべてのシミュレーションは LAMMPS を使用して行われました。検証には、構造解析、溶媒（DMF およびエタノール）中での安定性テスト、および機械的特性の計算が含まれました。自己集合研究では、16 の Zr クラスター、96 の BDC リガンド、および 2136 の DMF 分子からなる系を用い、NPT アンサンブルで 6 ns 間シミュレートしました。

主要な結果

• 構造精度: nb-UiO-FF は、純粋な UiO-66 の実験格子定数を約 0.046% の偏差（20.71 Å 対 20.7004 Å）で再現します。主要な原子対（Zr–$\mu_3$O、Zr–$\mu_3$OH、Zr–Zr、Zr–リガンド）の動径分布関数（RDF）は、DFT および実験データと妥当な一致を示し、偏差は一般的に 10% 以内です。
• 溶媒および欠陥安定性: この力場は、DMF およびエタノール溶媒の存在下、ならびに欠けたリンカーや欠けたノードを含む欠陥構造においても、骨格の完全性を維持します。欠陥モデルでは、電荷中性を維持するために、配位不足の中心を塩化物イオンまたは単座 BDC 断片でキャップしました。
• 機械的特性: 計算された弾性定数（$C_{11}$、$C_{12}$、$C_{44}$）および導出された機械的特性（体積弾性率：44 GPa、ヤング率：40 GPa）は、DFT ベンチマークと良好な一致を示しますが、Zr とダミー原子間の強い結合相互作用によって導入された剛性のため、わずかに高い値となっています。
• 転移性: このモデルは、異なるリンカー長にわたる転移性を示し、同系物である UiO-67 骨格の構造を成功裡に再現しました。
• 自己集合ダイナミクス: 動的シミュレーションにおいて、この力場は金属クラスターと有機リンカー間の配位結合の段階的な形成を成功裡に捉えました。本研究では、初期の二次構築単位（SBU）や非理想的な結合配置（運動学的に捕捉された状態）を含む一時的な構造モチーフを同定しました。しかし、配位数の完全な収束と大型のイオン凝集体の形成は、6 ns のシミュレーションウィンドウ内では達成されず、より長い時間スケールに対する強化サンプリング手法の必要性を示唆しています。

意義と主張
本論文は、nb-UiO-FF が静的構造モデルと反応的ダイナミクスの間のギャップを埋める、UiO-66 および UiO-67 の一貫性があり物理的に信頼できる記述を提供すると主張しています。以前のモデルと比較して削減されたダミー原子セット（Zr あたり 4 つ）を利用することで、著者らはモデルの効率性と反応過程のシミュレーション能力の間のバランスを達成しています。この力場は、従来の固定結合性力場ではアクセスできない自己集合メカニズムや欠陥形成の研究を可能にします。著者らは、この研究を MOF 合成の初期段階や点欠陥の工学を調査するための基盤として位置付けていますが、運動学的トラップを完全に克服し、完全な自己集合をモデル化するには、メタダイナミクスなどの強化サンプリング手法を用いたさらなる作業が必要であると認識しています。