Explicit pore pre-conditioning and inverse design of amorphous nanoporous carbon

技術要約：アモルファスナノ多孔質炭素の明示的な細孔事前調製および逆設計

問題提起
無秩序なナノ多孔質炭素は、エネルギー貯蔵、ガス分離、水浄化において極めて重要であり、その性能は比表面積だけでなく、空隙ネットワークの特定の幾何学的形状、特に細孔径分布（PSD）によって決定される。既存のモデリング手法には、以下のような重大な限界がある：

• 理想化された幾何学的モデル： 硬いスリット細孔や単純な形状を使用しており、実際の無秩序な炭素が持つ曲率、非六角形環、および連結性が欠如しているため、キャパシタンスなどの特性予測が不正確になる。
• 再構成法（例：逆モンテカルロ法）： 実測データに基づいて原子配置を適合させるが、既存のサンプルを必要とし、非一意性の問題がある。また、これらは局所的な秩序を制約するものの、グローバルな細孔ネットワークを制約するものではない。
• 模倣分子動力学（MD）： 物理的な合成プロセス（例：クエンチングやアニーリング）をシミュレートして、緩和された構造を生成する。しかし、細孔の幾何学はシミュレーションの「出力」であり、密度や熱スケジュールによって決定される。ユーザーは目標とするPSDを「入力」として指定することができず、特定の分布に一致させるには、プロトコルのパラメータに対する困難な探索が必要となる。また、任意の形状（例：バイモーダル分布）を生成できる模倣プロトコルはこれまで示されていない。

手法
著者らは、標準的なワークフローを反転させた生成フレームワークを提案している。すなわち、細孔ネットワークを定義として入力し、その周囲に炭素構造を成長させる。プロセスは以下の5つのステップで構成される：

1. 明示的な細孔事前調機（構築）：

   • 細孔配置： 目標とするPSD（任意の形状）から描画された球体（空隙を表す）を、目標とする細孔体積分率（$\phi$）に達するまで周期セル内に挿入する。
   • 環ごとの成長： これらの細孔の周囲に、炭素ネットワークを確率的に成長させる。環は、幾何学的基準（原子が細孔に入らないこと）およびエネルギー的基準（エネルギーが閾値以下であること）に基づいて受理される。
   • 配置ポテンシャル： 2種類の相互作用モデルをテストする：
     • 調和結合モデル（Harmonic Bonded Model）： 結合、角度、および二面角の項を用いて、親シートに対する環の幾何学を設定する。
     • Tersoff接触モデル（Tersoff-Contact Model）： Tersoffポテンシャルを使用するが、局所的な結合項（アンカーおよび近傍原子）を除外し、新しい環とネットワークの遠方部分との間の接触に対してのみ作用する。
       パラメータ： 構築は、探索パラメータ（細孔体積分率 $\phi$、細孔壁の深さ $A$ および範囲 $B$、二面角剛性 $k_\phi$、およびマージ係数 $m$）によって制御される。

2. 特性評価：

   • 生成された「プロキシ構造」は、ターゲットと同じルーチン（PSD、有効表面積、動径分布関数など）を用いて特性評価される。
   • カスタムPythonライブラリ nanocarbon が、これらの計算を効率的に実行する。

3. スコアリング：

   • プロキシ構造の特性とターゲット特性との間の無次元偏差の総和として、コスト関数が計算される。

4. 逆設計探索（最適化）：

   • 粒子群最適化（PSO）を用いて、コスト関数を最小化するように構築パラメータ空間を探索する。
   • 群（スウォーム）はパラメータセットを提案し、二次サロゲートモデルが候補の提案を補助する。探索の履歴は、どのパラメータが実際にターゲットによって制約されているか（識別性分析）を判断するために分析される。

5. アニーリング：

   • 最適化されたプロキシ構造に対し、細孔を反発壁として保持した状態で、短い分子動力学アニーリング（EDIPポテンシャルを用い、2500 Kで行う）を行う。
   • このステップにより、局所的なトポロジー（配位欠陥の修正）を修復し、設定されたグローバルな細孔幾何学を維持しながら、ネットワークを緻密化させる。

主な結果

• グローバル幾何学の回復： 探索は、クエンチ分子動力学（QMD）によって生成された参照構造のPSD、密度、および表面積を、約5%以内の誤差で正常に回収した。
• バイモーダル分布によるストレス・テスト： 調和配置ポテンシャルを用いることで、合成されたバイモーダルPSD（モードが7 Åと16 Å）の再構築に成功した。これは、模倣MDではこれまで達成できなかった形状である。一方、Tersoff接触構築では、この特定のターゲットに対して細孔を配置することに失敗した。これは、手法が配置ポテンシャルの選択に依存していることを浮き彫りにした。
• 識別性分析： 探索履歴の分析により、細孔体積分率（$\phi$）が、成功した実行においてターゲットによって制約される主要なパラメータであることが明らかになった。細孔壁の範囲（$B$）は一部の実行で制約されるが、他のパラメータ（二面角剛性、壁の深さ、マージ係数）は「鈍感（insensitive）」であった（コストはそれらに無関心であった）。この分析は、完全な構築を行う前に、失敗した探索（バイモーダルターゲットに対するTersoffの実行など）をフラグ立てすることも可能である。
• 局所的秩序とアニーリング：
  • プロキシ構造は、初期段階では不完全なトポロジー（2配位原子の高割合、4配位原子の低割合）を持つ。
  • 短いアニール（約2500 Kで15 ps）により、配位数が参照構造（約80%が3配位、12%が2配位）に一致するように回復し、骨格密度が増加する。
  • アニーリングにより、ネットワークの緻密化に伴って平均細孔径が大きな値の方へわずかにシフトする（0.8–1.5 Å）が、PSDの全体的な形状は維持される。
  • 最終的な構造は、アニーリング後に、細孔の開放性を維持するために反発壁に依存しなくなる（壁が保持するエネルギーは無視できる程度である）。

意義と主張
本論文は、細孔幾何学がシミュレーションの結果ではなく、設計仕様となる手法を確立したと主張している。

• 設計変数： PSDを「入力」として扱うことで、この手法は、既知の合成経路で作られる構造に限定されることなく、特定の用途（例：特定の電解質や輸送要件）に応じた炭素構造の設計を可能にする。
• 効率性： 細孔径分布の特性評価コストは構造のサイズに対して線形に増加するため、完全な模倣MDのトラジェトリと比較して、逆設計探索は計算コストを抑えて実行可能である。
• 限定的な範囲： 著者らは以下の制限事項を挙げている：本手法は現在、元素炭素に限定されていること、および最終的なアニール構造における環統計（特に5/7員環と6員環の比率）が、参照となるQMD構造よりも低いままであること。局所的な化学組成は依然として分子動力学によって提供されており、構築パラメータは必ずしも物理的な値ではなく、ジェネレータのためのチューニングノブである。
• 将来の展望： 本研究は、吸着測定から直接得られたターゲットや、特定の電解質のための電極設計など、設計プロセスを「イン・シリコ（in silico）」へと移行させる可能性を示唆している。