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 φ for which we can obtain a well-defined non-equilibrium static structure factor S0(k; φ). The evolution of the α-relaxation time τα(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 τα(k; φ, tw) as a function of φ exhibits two regimes corresponding to samples that have fully equilibrated within this waiting time (φ ≤ φ (c) (tw)), and to samples for which equilibration is not yet complete (φ ≥ φ (c) (tw)). The crossover volume fraction φ (c) (tw) increases with tw but seems to saturate to a value φ (a) ≡ φ (c) (tw → ∞) ≈ 0.582. We also find that the waiting time t eq w (φ) required to equilibrate a system grows faster than the corresponding equilibrium relaxation time,1.43 , and that both characteristic times increase strongly as φ approaches φ (a) , thus suggesting that the measurement of equilibrium properties at and above φ (a) is experimentally impossible. exp ≈ 0.58 [6,7], although a number of intrinsic experimental uncertainties render the precise determination of φ (a) exp a topic of recurrent scientific discussion [4,[6][7][8][9][10].

The binary hard-sphere mixture is one of the simplest representations of a many-body system with competing time and length scales. This model is relevant to fundamentally understand both the structural and dynamical properties of materials, such as metallic melts, colloids, polymers and bio-based composites. It also allows us to study how different scales influence the physical behavior of a multicomponent glass-forming liquid; a question that still awaits a unified description. In this contribution, we report on distinct dynamical arrest transitions in highly asymmetric binary colloidal mixtures, namely, a single glass of big particles, in which the small species remains ergodic, and a double glass with the simultaneous arrest of both components. When the mixture approaches any glass transition, the relaxation of the collective dynamics of both species becomes coupled. In the single glass domain, spatial modulations occur due to the structure of the large spheres, a feature not observed in the two-glass domain. The relaxation of the self dynamics of small and large particles, in contrast, become decoupled at the boundaries of both transitions; the large species always displays dynamical arrest, whereas the small ones appear arrested only in the double glass. Thus, in order to obtain a complete picture of the distinct glassy states, one needs to take into account the dynamics of both species.

Asymmetric binary mixtures of hard-spheres exhibit several interesting thermodynamic phenomena, such as multiple kinds of glassy states. When the degrees of freedom of the small spheres are integrated out from the description, their effects are incorporated into an effective pair interaction between large spheres known as the depletion potential. The latter has been widely used to study both the phase behavior and dynamic arrest of the big particles. Depletion forces can be accounted for by a contraction of the description in the multicomponent Ornstein-Zernike equation [R. Castañeda-Priego, A. Rodríguez-López, and J. M. Méndez-Alcaraz, Phys. Rev. E 73, 051404 (2006)]. Within this theoretical scheme, an approximation for the difference between the effective and bare bridge functions is needed. In the limit of infinite dilution, this difference is irrelevant and the typical Asakura-Osawa depletion potential is recovered. At higher particle concentrations, however, this difference becomes important, especially where the shell of first neighbors is formed, and, as shown here, cannot be simply neglected. In this work, we use a variant of the Verlet expression for the bridge functions to highlight their importance in the calculation of the depletion potential at high densities and close to the spinodal decomposition. We demonstrate that the modified Verlet closure predicts demixing in binary mixtures of hard spheres for different size ratios and compare its predictions with both liquid state and density functional theories, computer simulations, and experiments. We also show that it provides accurate correlation functions even near the thermodynamic instability; this is explicitly corroborated with results of molecular dynamics simulations of the whole mixture. Particularly, our findings point toward a possible universal behavior of the depletion potential around the spinodal line.

We report the combined results of molecular dynamics simulations and theoretical calculations concerning various dynamical arrest transitions in a model system representing a dipolar fluid, namely, N (soft core) rigid spheres interacting through a truncated dipole-dipole potential. By exploring different regimes of concentration and temperature, we find three distinct scenarios for the slowing down of the dynamics of the translational and orientational degrees of freedom: At low (η = 0.2) and intermediate (η = 0.4) volume fractions, both dynamics are strongly coupled and become simultaneously arrested upon cooling. At high concentrations (η ≥ 0.6), the translational dynamics shows the features of an ordinary glass transition, either by compressing or cooling down the system, but with the orientations remaining ergodic, thus indicating the existence of partially arrested states. In this density regime, but at lower temperatures, the relaxation of the orientational dynamics also freezes. The physical scenario provided by the simulations is discussed and compared against results obtained with the self-consistent generalized Langevin equation theory, and both provide a consistent description of the dynamical arrest transitions in the system. Our results are summarized in an arrested states diagram which qualitatively organizes the simulation data and provides a generic picture of the glass transitions of a dipolar fluid. PACS numbers: xxxxxxxx yyyyyyyy