[論文レビュー] Freezing of Polydisperse Hard Spheres

N. G. Almarza, E. Enciso|ArXiv.org|May 19, 1999

Textile materials and evaluations被引用数 31

ひとこと要約

本研究は、コンピュータシミュレーションおよび高精度な状態方程式を用いて、多分散性硬球系における流体-結晶相平衡を検討し、結晶相の組成を特定するための新しい手法を開発した。先行研究とは対照的に、流体相に終末的多分散性が存在しないことが判明した。以前の研究で観察された限界は、物理的に不適切な化学ポテンシャルモデルの結果であり、真の熱力学的制約ではない。

ABSTRACT

The fluid - crystal equilibria of polydisperse mixtures of hard spheres have been studied by computer simulation of the solid phase and using an accurate equation of state for the fluid. A new scheme has been developed to evaluate the composition of crystalline phases in equilibrium with a given polydisperse fluid. Some common assumptions in theoretical approaches and their results are discussed on the light of the simulation results. Finally, no evidence of the existence of a terminal polydispersity in the fluid phase is found for polydisperse hard spheres, the disagreement of this finding with previous molecular simulation results is explained in terms of the inherent limitations of some ways of modeling the chemical potential as a function of the particle size.

研究の動機と目的

• 正確なシミュレーションおよび状態方程式法を用いて、多分散性硬球系における流体-結晶相平衡を再検討すること。
• 特定の流体多分散性と平衡をとる結晶相の組成を特定するための新しい計算手法を開発すること。
• シミュレーションデータを用いて、相の組成に関する一般的な理論的仮定の妥当性を批判的に評価すること。
• 先行の分子シミュレーション研究が流体相に終末的多分散性を主張する点における矛盾を解消すること。
• 過去の研究で観察された見かけの相境界の物理的起源を明確にすること、特に化学ポテンシャルモデルの役割を特定すること。

提案手法

• 流体相の熱力学的性質を高精度に表現するため、粒子径分布のモーメントでパrameter化されたBMCSL状態方程式を用いる。
• 固定圧力、全粒子数、および基準成分に対する化学ポテンシャル差を固定する半グランドアンサンブル法を適用する。
• 圧力および多分散性に依存する強制化学ポテンシャル分布（ICPD）を導入し、目的の流体径分布が達成されるようにする。
• 共存線を(βp, λ)平面で追跡するために、Clausius-Clapeyron統合スキームを用いる。ここでλは多分散性を定量化するパラメータである。
• 共存路を進めるために、化学ポテンシャル係数の圧力および多分散性に関する微分係数を数値的に計算する。
• より小さな粒子を豊富に含む面心立方結晶を用いた制御シミュレーションを実施し、相の競合や転移の有無を検証する。

実験結果

リサーチクエスチョン

• RQ1 多分散性硬球系の流体相に、先行研究が示唆するように終末的多分散性が実際に存在するのか？
• RQ2 強制化学ポテンシャル分布の形が、相共存および見かけの限界の予測にどのように影響を与えるか？
• RQ3 一般的に用いられる理論的仮定（例えば、流体相と固体相の多分散性が等しい）が、現実的なシミュレーションにおいてどの程度成立するのか？
• RQ4 凍結過程において、複数の結晶相が安定性をめぐって競合する可能性はあり、これが相律の制約を変えるのか？
• RQ5 なぜ過去のシミュレーションでは終末的多分散性が観察されたのか？これは物理的限界なのか、モデル化の誤りによるものなのか？

主な発見

• 本研究で検討された条件下では、多分散性硬球系の流体相に終末的多分散性の証拠は得られなかった。
• 先行研究で観察された見かけの終末的多分散性は、特定の強制化学ポテンシャル分布の形（式6）に起因し、最大達成可能な流体充填率が人工的に制限されているためである。
• 参照文献[4]で用いられた化学ポテンシャルモデルは、最大充填率が流体共存曲線と交差するように設計されており、真の熱力学的限界を反映するものではなく、誤った終点を生じさせている。
• 本研究では、化学ポテンシャル展開の係数a₃が正であり、硬球反発を相殺するために大きな粒子が好まれる傾向を示す。これは、過去のICPD形式が引き起こす物理的に不適切な振る舞いとは対照的である。
• 小さな粒子を豊富に含む結晶から開始したシミュレーションは、いずれも同じ固体相に再平衡化するか、融解するだけであり、異なる結晶構造への相転移は本研究範囲内で観察されなかった。
• 同じ相対的多分散度を示す結晶相でも、圧力に応じて異なる流体径分布に対応する。高圧状態ではより負の歪度とより大きなm₁値を示す。