[論文レビュー] The fluid-fluid interface in a model colloid-polymer mixture: Application of grand canonical Monte Carlo to asymmetric binary mixtures

R. L. C. Vink, Jürgen Horbach|arXiv (Cornell University)|Oct 17, 2003

Phase Equilibria and Thermodynamics参考文献 3被引用数 77

ひとこと要約

本稿では、非対称コロイド-ポリマー混合系をシミュレートするため、集団的クラスターモーブとアンバレラサンプリングを組み合わせたグランドカノニカルモンテカルロ法を提示する。この手法により、アサクラ・オオサワ模型における流体-流体相分離と界面張力の高精度な計算が可能となり、3次元イジング臨界行動が確認された。また、密度汎関数理論（DFT）は界面張力を約30%過大評価していることが判明し、従来の定量的主張に疑問を呈する。

ABSTRACT

We present a Monte Carlo method to simulate asymmetric binary mixtures in the grand canonical ensemble. The method is used to study the colloid-polymer model of Asakura and Oosawa. We determine the phase diagram of the fluid-fluid unmixing transition and the interfacial tension, both at high polymer density and close to the critical point. We also present density profiles in the two-phase region. The results are compared to predictions of a recent density functional theory.

研究の動機と目的

• 標準的なグランドカノニカル・モーブが、大きな粒子の挿入に起因する大きな問題を引き起こす非対称二成分混合系（コロイド-ポリマー系など）のシミュレーションに適した、堅牢なモンテカルロ法の開発。
• 臨界点付近におけるアサクラ・オオサワ模型の流体-流体共存相図を高精度に特定すること。
• 二相領域における共存液体および気相反心間の界面張力を高精度に計算すること。
• アサクラ・オオサワ模型が、普遍性により予測されるように3次元イジング臨界行動を示すかどうかの検証。
• シミュレーション結果と最近の密度汎関数理論（DFT）予測を比較し、その定量的精度を評価すること。

提案手法

• 本手法は、小さな粒子（ポリマー）のクラスタを交換する集団的モンテカルロモーブを採用し、大きな粒子（コロイド）の重なりを避けて効率的に挿入を可能にする。
• これらのモーブに加え、共存領域全体でのサンプリングを向上させ、収束性を高めるためにアンバレラサンプリング技術を組み合わせる。
• 界面張力は、シミュレーションボックス内におけるコロイド体積分率の最大確率と最小確率の対数比を用いたバインダーの方法により計算する。
• 相図は、異なる化学ポテンシャルと密度におけるシミュレーションデータを再重み付けすることで得る。
• 臨界行動の評価のため、コロイド密度プロファイルを双曲正 tangent 関数にフィットさせ、界面幅を抽出する。
• 界面の安定化と有限サイズ効果の低減のため、アスペクト比が1:3の周期的境界条件を用いたボックスを採用する。

実験結果

リサーチクエスチョン

• RQ1 アサクラ・オオサワ模型は、流体-流体相分離転移において3次元イジング臨界行動を示すか？
• RQ2 最近の密度汎関数理論（DFT）が予測する界面張力は、直接的なシミュレーションと比較してどの程度正確か？
• RQ3 集団的モーブを用いたグランドカノニカルモンテカルロ法は、非対称二成分混合系において、標準的手法のサンプリング制限を克服できるか？
• RQ4 アサクラ・オオサワ模型の流体-流体転移における共存曲線（ビノード）の正確な形状は、臨界点付近でどのように変化するか？
• RQ5 界面幅は臨界点からの距離に従ってどのようにスケーリングされるか？また、3次元イジングのべき乗則に従うか？

主な発見

• シミュレーションにより3次元イジング臨界行動が確認され、界面幅が W ∝ (η_p^r - η_p,crit^r)^{-ν} とスケーリングし、ν ≈ 0.63 となることが判明。これはイジング普遍性クラスと整合的である。
• 無次元化された界面張力 γ* は、最近のDFT予測から約30%のずれを示し、DFTが界面張力を過大評価していることを示唆する。
• シミュレーションから得られたビノード曲線は、DFTの予測よりも平坦であり、β ≈ 0.325（イジングクラス）であるのに対し、DFTでは β = 0.5 となる。
• コロイド体積分率 η_c の確率分布 P(η_c) は、二つのピークの間で平坦な領域を示し、バインダーの方法による界面張力推定には不可欠である。
• 界面幅 W は、臨界点から離れていても系のサイズ L_z に強く依存しており、おそらくキャピラリー波効果によるものと推定される。
• 本手法により、二相領域における相図および界面張力の高精度な決定が可能となり、ギブズ集合とカノニカル集合の制限を克服した。