A route to the thermodynamics of colloid-polymer mixtures from structural information

技術要約：構造情報からコロイド・ポリマー混合物の熱力学への経路

問題提起
液相理論は、微視的な構造（例えば、動径分布関数 $g(r)$ や静的構造因子 $S(k)$）とマクロな熱力学的性質との間の厳密な接続に依存している。しかし、コロイド・ポリマー混合物のような複雑な流体の熱力学や相挙動を予測することは、大きな課題である。従来のアプローチでは、通常、基礎となる相互作用ポテンシャルの明示的な知識を必要とする。多くの実験系や複雑なソフトマター系においては、これらの有効相互作用は未知であったり、状態点に依存したり、あるいは多体効果の影響を強く受けているため、単純な対ポテンシャル近似（ポテンシャルの平均力など）では不正確になる。さらに、正確な熱力学的性質を得るためには、広い範囲の状態点にわたって熱力学的積分（TI）を行う必要があり、これは計算コストが高い。逆問題的手法（Inverse Boltzmann Iteration: IBIなど）は、構造データからポテンシャルを再構成できるが、得られるポテンシャルは特定の熱力学的状態への依存性を保持してしまうことが多く、より広い相挙動に対する予測能力が制限される。

手法
著者らは、明示的な相互作用ポテンシャルを必要とせず、単一の熱力学的状態点から得られた構造相関を用いて、熱力学と相挙動を予測するフレームワークを開発した。このアプローチは、動径分布関数 ($g(r)$) または静的構造因子 ($S(k)$) を入力として自由エネルギー汎関数を構築するものである。

本定式化では、2種類の相互作用を区別している：

1. 排除体積相互作用： 硬いコアによる反発を持つ系（$r < \sigma_u$ で $g(r)=0$ となるもの）。
2. ソフトな相互作用： 粒子が互いに重なり合える系（$g(r)$ が有限の値を持つもの）。

排除体積相互作用を持つ系については、全相互作用ポテンシャルを、短距離の反発コア ($\phi_{ref}$) と残りの裾部分 ($\phi_{tail}$) に分解する。コアは、有効硬球参照系へと写像される。本研究の新規な側面は、有効硬球直径 ($\sigma_{opt}$) の決定にある。既知のポテンシャルに対して Barker-Henderson 法を用いる標準的な処方ではなく、$\sigma_{opt}$ は、IBI によって再構成されたポテンシャル（これは反発寄与を分離するための中間ステップとしてのみ使用される）から導かれる参照硬球の RDF との差を最小化することによって決定される。

その後、過剰自由エネルギー ($F^{exc}$) は、参照寄与 ($F^{exc}_{ref}$) と裾部分の寄与 ($F^{exc}_{tail}$) に分割される：

• 参照系： $F^{exc}_{ref}$ は、有効硬球流体に対する基本測定理論 (FMT) または Carnahan-Starling 状態方程式を用いて計算される。
• 裾部分の寄与： 著者らは、裾部分を表す積分相互作用強度 ($A$) を決定するために、2つの補完的な経路を提案している：
  1. 熱力学的積分 (TI) ルート： 構造相関が排除体積効果によって支配されていると仮定し、全系の $g(r)$ を用いて結合平均化された対相関関数を近似する。
  2. 直接相関関数 (DCF) ルート： Ornstein-Zernike (OZ) 関係式を通じて得られる、対直接相関関数 $c^{(2)}(r)$ を利用する。過剰自由エネルギーは、全系と参照系の間の DCF の差 ($\Delta c^{(2)}$) を用いて表現される。

ソフトな相互作用（例：ポリマー）については、ランダム位相近似 (RPA) の関係を用い、積分強さを DCF または第2ビリアル係数に直接結びつける。

本フレームワークは、すべての対の種類に関する積分相互作用強度の行列 ($A_{ij}$) を構築することにより、多成分混合物（コロイド・ポリマー系）へと拡張される。参照自由エネルギーは、硬球相関、デプレッション効果（自由体積理論による）、およびソフトなポリマー間相互作用を考慮する。そして、裾部分の寄与が摂動的に加えられる。

主な貢献

• 単一状態点での予測： 主要な貢献は、広範な熱力学的積分や相互作用ポテンシャルの事前知識を必要とせず、単一の状態点における構造データから状態方程式と相挙動を推論する手法を確立したことである。
• ポテンシャルを用いない逆問題フレームワーク： 状態依存的な有効ポテンシャルの再構成を回避する。代わりに、構造相関から自由エネルギー汎関数を直接構築する。ここでは、参照系を定義するための暫定的なツールとしてのみ IBI を使用する。
• ルートの堅牢性： DCF ルートと TI ルートが積分相互作用強度について一致した結果を与えることを示し、自由エネルギー計算のための構造カーネルとしての $c^{(2)}(r)$ または $g(r)$ の妥当性を検証した。
• 混合物への適用性： 非理想的な（ソフトな）ポリマー間相互作用を含む系において、自己相関と相互相関の両方を考慮した、二成分コロイド・ポリマー混合物の扱いを成功させた。

結果
本フレームワークは、以下の分子シミュレーションデータを用いて検証された：

1. 一成分系： Lennard-Jones (12-6) および一般化された Mie (24-12) ポテンシャルのテストにおいて、TI ルートと DCF ルートは、低密度から高密度にかけてシミュレーションデータと良好に一致する状態方程式を予測した。希薄極限近似に基づくビリアルルートは、高密度では性能が低下した。12-6 系における予測された共存曲線はシミュレーションの傾向と一致したが、臨界点は、おそらくポテンシャルの打ち切りによる違いにより、わずかに過小評価された。
2. コロイド・ポリマー混合物： ガウス型ポリマーポテンシャルと Mie コロイドポテンシャルを用いた二成分混合物に適用された。本フレームワークは、固定されたコロイド分率における状態方程式（圧力 vs 密度）を正確に予測した。
3. 相挙動の予測： コロイド・ポリマー混合物の流体－流体共存（ビノダル）を予測することに成功した。単一の入力状態からの構造相関のみから導出された相図は、ポリマー間相互作用が非理想的であっても、ポリマー豊富／コロイド乏しい相と、コロイド豊富／ポリマー乏しい相への偏析を正しく特定した。

意義と主張
本論文は、微視的な相互作用が事前に未知であっても、実験的に測定された構造情報から相挙容を直接推論するための実行可能な経路を確立したと主張している。これは、有効相互作用の決定が困難であり、かつ系の状態によって変化し得る実験的なソフトマター系において特に重要である。

著者らは、本定式化が、明示的な相互作用ポテンシャルや系固有の熱力学モデルなしに、自由エネルギー汎関数を構築するための一般的な経路を提供することを強調している。また、現在の定式化は等方的な相互作用に限定されており、配向の自由度や、対相関だけでは不十分な可能性がある強く相関した系（例：高度に電荷を持ったコロイド）を明示的に考慮していないことを謙虚に述べている。しかし、構造相関が等方成分によってよく記述される場合には、異方的な系にも適用可能であると彼らは考えている。本研究は、単一の状態点における構造情報が、複雑な流体のより広い熱力学的応答を予測するために十分な情報を含んでいることを示唆しており、計算コストの高いシミュレーションや不確実なポテンシャルモデルに代わる実用的な選択肢を提示している。