[論文レビュー] Finite-size effects and thermodynamic limit in one-dimensional Janus fluids

Riccardo Fantoni, M A G Maestre+1 more|arXiv (Cornell University)|Jul 29, 2021

Pickering emulsions and particle stabilization参考文献 48被引用数 5

ひとこと要約

本稿では、1次元のジャンスフリクイドに対して、クエンチドおよびアンネールド系の有限サイズにおける正確なギブズ自由エネルギーを導出し、両者が同じ熱力学的極限式に収束することを証明している。正確な統計力学とモンテカルロシミュレーションを用いて、熱力学的極限の頑健性を確認し、特にハードロッドおよびスクェアーベル・インタラクションを有する Kern–Frenkel モデルにおいて、クエンチドとアンネールド系の間の等価性が有限Nでも成立することを検証している。

ABSTRACT

The equilibrium properties of a Janus fluid made of two-face particles confined to a one-dimensional channel are revisited. The exact Gibbs free energy for a finite number of particles $N$ is exactly derived for both quenched and annealed realizations. It is proved that the results for both classes of systems tend in the thermodynamic limit ($N o\infty$) to a common expression recently derived (Maestre M A G and Santos A 2020 J Stat Mech 063217). The theoretical finite-size results are particularized to the Kern--Frenkel model and confirmed by Monte Carlo simulations for quenched and (both biased and unbiased) annealed systems.

研究の動機と目的

• クエンチドおよびアンネールドの1次元ジャンスフリクイドの有限Nにおける正確なギブズ自由エネルギーを導出すること。
• 有限系サイズにおいてもクエンチド-アンネールド等価性が成立するかを調査すること。
• バイアスありおよびバイアスなしのアンネールド系の両方において、理論的予測をモンテカルロシミュレーションと照合すること。
• アンネールド系における面の向き確率q₁の関数としての平均モル分率⟨x₁⟩の頑健性を検討すること。

提案手法

• 正確な配置積分を用いて、クエンチド二成分混合系の等温等圧系分配関数を導出する。
• 変数変換および母関数の応用により、有限Nにおける正確なギブズ自由エネルギーを計算する。
• 固定確率q₁およびq₂におけるスピン配置の平均化を用いて、アンネールド系への形式の拡張を行う。
• 理論的結果の具体化のため、ハードロッドおよびスクェアーベル・ポテンシャルを有する Kern–Frenkel モデルを適用する。
• クエンチド（MCq）、アンネールド（MCa）、および代替クエンチド-アンネールド（MCaq）プロトコルにおける等温等圧系で古典的モンテカルロシミュレーションを実装する。
• 観測量⟨x₁⟩、v、およびu_exの収束を保証するため、相関時間τOおよび統計誤差を推定する。

実験結果

リサーチクエスチョン

• RQ1 有限Nにおけるクエンチドおよびアンネールドの1次元ジャンスフリクイドの正確なギブズ自由エネルギーを導出可能か？
• RQ2 有限Nにおいてクエンチド-アンネールド等価性は破綻するか、それとも両系が同じ熱力学的極限に収束するか？
• RQ3 クエンチドおよびアンネールド両設定において、有限サイズ自由エネルギーに関する理論的予測がモンテカルロシミュレーションで確認されるか？
• RQ4 特に有限Nにおいて、アンネールド系における面の向き確率q₁の変化に対して、平均モル分率⟨x₁⟩はどの程度頑健か？

主な発見

• 正確な統計力学および配置積分を用いて、クエンチドおよびアンネールドの両方の1次元ジャンスフリクイドについて、有限Nにおける正確なギブズ自由エネルギーが導出された。
• クエンチドおよびアンネールド両系がN → ∞のとき、同じ熱力学的極限式に収束することを確認し、熱力学的挙動の普遍性を裏付けた。
• Kern–Frenkel モデル（ハードロッドおよびスクェアーベル・インタラクション）に対する理論的予測は、クエンチド、バイアスなしアンネールド、バイアスありアンネールドのすべての系においてモンテカルロシミュレーションで妥当性を確認された。
• 有限Nに対しても、アンネールド系における平均モル分率⟨x₁⟩は頑健であり、理論的予測と良好に一致した。
• モンテカルロシミュレーションにおける相関時間τOは、誤差推定に不可欠であり、収束を保証するためには長大な平衡化時間が必要であった（T=0.2で最大2×10⁷ MCステップ）。
• 大Nにおける漸近的展開および重み関数を用いて、バイアスありアンネールド系における⟨x₁⟩のヒューリスティックな近似式が導出され、経験的パラメータを調整することでシミュレーションデータと良好な一致を示した。