Equilibration of Concentrated Hard Sphere Fluids

TL;DR

This study uses molecular dynamics to analyze how concentrated hard-sphere fluids approach equilibrium near the glass transition, revealing that equilibration times grow rapidly and become practically unattainable at high densities.

Contribution

It provides a systematic analysis of the equilibration process in hard-sphere fluids near the glass transition, highlighting the growth of equilibration times and the saturation of the crossover volume fraction.

Findings

Equilibration times grow faster than relaxation times as density increases.

The crossover volume fraction saturates around 0.582.

Equilibrium properties are practically unmeasurable above this saturation point.

Abstract

We report a systematic molecular dynamics study of the isochoric equilibration of hard-sphere fluids in their metastable regime close to the glass transition. The thermalization process starts with the system prepared in a non-equilibrium state with the desired final volume fraction {\phi} but with a prescribed non-equilibrium static structure factor S_0(k; {\phi}). The evolution of the {\alpha}- relaxation time {\tau}{\alpha} (k) and long-time self-diffusion coefficient DL as a function of the evolution time tw is then monitored for an array of volume fractions. For a given waiting time the plot of {\tau}{\alpha} (k; {\phi}, tw) as a function of {\phi} exhibits two regimes corresponding to samples that have fully equilibrated within this waiting time ({\phi} \leq {\phi}(c) (tw)), and to samples for which equilibration is not yet complete ({\phi} \geq {\phi}(c) (tw)). The crossover volume fraction {\phi}(c) (tw) increases with tw but seems to saturate to a value {\phi}(a) \equiv {\phi}(c) (tw \rightarrow \infty) \approx 0.582. We also find that the waiting time t^(eq)_w({\phi}) required to equilibrate a system grows faster than the corresponding equilibrium relaxation time, t^(eq)({\phi}) \approx 0.27 \times [{\tau}{\alpha} (k; {\phi})]^1.43, and that both characteristic times increase strongly as {\phi} approaches {\phi}^(a), thus suggesting that the measurement of equilibrium properties at and above {\phi}(a) is experimentally impossible.