Accurate Solvation Properties in supercritical CO$_2$ with Molecular Density Functional Theory

技術要約：分子密度汎関数理論を用いた超臨界CO2における正確な溶媒和特性

問題提起
超臨界二酸化炭素（scCO2）は、抽出からナノ粒子合成に至るまで、産業用途において非常に効率的かつ環境に優しい溶媒である。これらのプロセスを最適化するには、様々な熱力学的条件下での溶媒和特性、特に溶媒和自由エネルギーと溶媒和構造を正確に予測することが求められる。分子動力学（MD）シミュレーションは、高い予測精度と微視的な洞察を提供するが、高度なサンプリング戦略（ウィドム挿入法や熱力学的積分法など）を用いた熱力学的性質の計算においては、膨大な計算コストに直面する。一方、既存のパラメトリック手法（状態方程式、密度ベース、または溶媒和モデル）は、計算コストは低いものの、転移性に乏しく、新しい溶質やシステムの変更に対して再パラメータ化が必要になることが多い。したがって、MDの精度と、大規模なパラメトリック研究に必要な計算効率を兼ね備えた手法が必要とされている。

手法
本研究では、分子流体に適合させた古典密度汎関数理論（cDFT）の形式である、分子密度汎関数理論（MDFT）を採用している。このアプローチでは、平衡分子密度 $\rho(\mathbf{r}, \omega)$ と溶媒和自由エネルギーを決定するために、溶媒和汎関数 $\Delta\Omega_{solv}[\rho]$ を最小化する。この汎関数は、以下の3つの成分に分解される：

1. 理想的寄与 ($\Delta\Omega_{id}$): 溶媒密度のエントロピーを考慮する。
2. 外部寄与 ($F_{ext}$): 標準的な力場（CO2にはEPM2、分子溶質にはOPLS-AA、水にはSPC/E）から導かれる溶質－溶媒相互作用を表す。
3. 過剰寄与 ($\Delta\Omega_{exc}$): 溶媒－溶媒相互作用を表す。この項は、**均一参照流体（HRF）**近似を用いて近似される。これは、溶質近傍の密度揺らぎが均一なバルクの状態と同様であると仮定するものである。HRF項は、事前のMDシミュレーションおよび積分方程式論から得られた純粋バルク溶媒の直接対相関関数（DCF）に依存する、二次展開として近似される。

HRF近似の限界（高次の揺らぎを無視している点）に対処するため、著者らはブリッジ汎関数（$\Delta\Omega_{bridge}$：厳密な過剰汎関数とHRF項との差）に関する3つの近似レベルを調査した：

• ブリッジなし: $\Delta\Omega_{bridge} = 0$。
• 硬球ブリッジ (HSB): バルクの自由エネルギーに較正された硬球過剰汎関数（FMT）を用いてブリッジを近似する。
• 粗視化ブリッジ (CGB): バルクの状態方程式の多項式展開から導出された、粗視化密度場の汎関数としてブリッジをモデル化する。

MDFT計算は、空間格子 $0.3$Å、最大角度指数$n_{max}=4$の立方体ボックス（$L=30$ Å）で行われた。これらの結果は、4000個のCO2分子を用い、ベネット受容比法による自由エネルギー計算を用いたMDベンチマークと比較された。

主な貢献

• 形式論の適用: 本研究は、空間位置と分子配向の両方に依存する分子cDFT形式をscCO2に適用することに成功した。
• ブリッジ汎関数のパラメータ化: 著者らは、scCO2に対して2つの異なるブリッジ汎関数（HSBおよびCGB）を構築し、パラメータ化した。CGB汎関数が、粗視化長さ $\sigma_g = 1.45$ Å を用いて、状態方程式から導かれる均一ブリッジ寄与を正確に再現することを実証した。
• 計算効率: 本研究はMDFTの計算上の利点を定量化しており、単一の計算がノートパソコン上で1分未満（約0.01 CPU・時間）で完了するのに対し、同等のMDシミュレーションには約640 CPU・時間を要することを指摘している。

結果
MDFFフレームワークは、2種類の溶質クラス（球状のアルゴン様溶質、および水、エタノール、トルエン、CO2自身を含む9種類の分子溶質）に対して、準臨界状態点（$n_0 = 0.8 n_c$, $T = 1.05 T_c$）でMDベンチマークと比較検証された。

• 構造的特性: 全ての溶質について、MDFTによって予測された動径分布関数（RDF）は、MDの結果とほぼ完全に重なった。この一致は、使用したブリッジ近似（HRF、HSB、またはCGB）に関わらず維持された。
• 溶媒和自由エネルギー:
  • 球状溶質については、HRF近似（ブリッジなし）は大きな溶質（$\sigma > 7$ Å）に対して溶媒和自由エネルギーを系統的に過大評価した。しかし、HSBおよびCGBによる補正は、幅広いサイズおよび相互作用エネルギーにおいて、MDベンチマークとほぼ区別がつかない結果を与えた。
  • 分子溶質については、3つのMDFTバリアントすべて（HRF、HSB、CGB）が、平均相対誤差約1%でMDの溶媒和自由エネルギーを再現した。
  • 特筆すべき点として、scCO2中の分子溶質において、ブリッジ汎関数の寄与は無視できるほど小さいことが判明した。つまり、HRF近似単独で熱力学を捉えるのに十分であった。これは、液体の水におけるMDFT研究とは対照的であり、scCO2には強い方向性を持つ水素結合ネットワークが存在しないことに起因する。
• 精度とコスト: MDFTアプローチは、構造的特性および熱力学的特性の両方において、MDと比較して計算量を約60,000分の1に削減しながら、量子的な精度に近い精度を達成した。

意義と主張
本論文は、MDFTがscCO2における溶媒和をモデリングするための、堅牢かつ計算効率の高いMDに代わる選択肢であることを主張している。主な意義は、バルク溶媒の相関が正確に知られていれば、均一参照流体（HRF）近似が、scCO2の溶媒和熱力学を回復するのに十分であることを示した点にある。この知見は、他の流体（水など）で必要とされる複雑なブリッジ補正が、scCO2においてはその特有の流体特性により不要である可能性を示唆している。

著者らは、現在の研究は相図全体にわたる転移性を確立するものではなく、代表的な準臨界状態点における手法の妥当性を検証するものであることを明示している。結論として、このフレームワークは、界面や閉じ込めといったより複雑なシナリオへMDFTを拡張するための基礎を提供するものであり、今後の発展においては、熱力学的条件間の転移性を達成するために機械学習ベースの汎関数を活用できる可能性があるとしている。特定された直接的な有用性は、迅速かつ正確な溶媒和特性のスクリーニングを通じて、scCO2ベースの産業プロセスの合理的設計および最適化を促進できる可能性である。