Pedro E. Ramírez-González

Aging of a homogeneously quenched colloidal glass-forming liquid.

The nonequilibrium self-consistent generalized Langevin equation theory of colloid dynamics is used to describe the nonstationary aging processes occurring in a suddenly quenched model colloidal liquid with hard-sphere plus short-ranged attractive interactions, whose static structure factor and van Hove function evolve irreversibly from the initial conditions before the quench to a final dynamically arrested state. The comparison of our numerical results with available simulation data are highly encouraging.

Simplified self-consistent theory of colloid dynamics.

One of the main elements of the self-consistent generalized Langevin equation (SCGLE) theory of colloid dynamics [Phys. Rev. E 62, 3382 (2000); 72, 031107 (2005)] is the introduction of exact short-time moment conditions in its formulation. The need to previously calculate these exact short-time properties constitutes a practical barrier for its application. In this Brief Report, we report that a simplified version of this theory, in which this short-time information is eliminated, leads to the same results in the intermediate and long-time regimes. Deviations are only observed at short times, and are not qualitatively or quantitatively important. This is illustrated by comparing the two versions of the theory for representative model systems.

Density-temperature-softness scaling of the dynamics of glass-forming soft-sphere liquids.

We employ the principle of dynamic equivalence between soft-sphere and hard-sphere fluids [Phys. Rev. E 68, 011405 (2003)] to describe the interplay of the effects of varying the density n, the temperature T, and the softness (characterized by a softness parameter ν(-1)) on the dynamics of glass-forming soft-sphere liquids in terms of simple scaling rules. The main prediction is the existence of a dynamic universality class associated with the hard-sphere fluid, constituted by the soft-sphere systems whose dynamic parameters depend on n, T, and ν only through the reduced density n*≡nσ(HS)(T*,ν). A number of scaling properties observed in recent experiments and simulations involving glass-forming fluids with repulsive short-range interactions are found to be a direct manifestation of this general dynamic equivalence principle.

Glass transition in soft-sphere dispersions

The concept of dynamic equivalence among mono-disperse soft-sphere fluids is employed in the framework of the self-consistent generalized Langevin equation (SCGLE) theory of colloid dynamics to calculate the ideal glass transition phase diagram of model soft-sphere colloidal dispersions in the softness-concentration state space. The slow dynamics predicted by this theory near the glass transition is compared with available experimental data for the decay of the intermediate scattering function of colloidal dispersions of soft-microgel particles. Increasing deviations from this simple scheme occur for increasingly softer potentials, and this is studied here using the Rogers-Young static structure factor of the soft-sphere systems as the input of the SCGLE theory, without assuming a priori the validity of the equivalence principle above.

Equilibration of concentrated hard-sphere fluids.

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 nonequilibrium state with the desired final volume fraction ϕ for which we can obtain a well-defined nonequilibrium static structure factor S(0)(k;ϕ). The evolution of the α-relaxation time τ(α)(k) and long-time self-diffusion coefficient D(L) as a function of the evolution time t(w) is then monitored for an array of volume fractions. For a given waiting time the plot of τ(α)(k;ϕ,t(w)) as a function of ϕ exhibits two regimes corresponding to samples that have fully equilibrated within this waiting time [ϕ≤ϕ(c)(t(w))] and to samples for which equilibration is not yet complete [ϕ≥ϕ(c)(t(w))]. The crossover volume fraction ϕ(c)(t(w)) increases with t(w) but seems to saturate to a value ϕ(a)≡ϕ(c)(t(w)→∞)≈0.582. We also find that the waiting time t(w)(eq)(ϕ) required to equilibrate a system grows faster than the corresponding equilibrium relaxation time, t(w)(eq)(ϕ)≈0.27[τ(α)(eq)(k;ϕ)](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.

Equilibration and aging of dense soft-sphere glass-forming liquids.

The recently developed nonequilibrium extension of the self-consistent generalized Langevin equation theory of irreversible relaxation [Ramírez-González and Medina-Noyola, Phys. Rev. E 82, 061503 (2010); Ramírez-González and Medina-Noyola, Phys. Rev. E 82, 061504 (2010)] is applied to the description of the irreversible process of equilibration and aging of a glass-forming soft-sphere liquid that follows a sudden temperature quench, within the constraint that the local mean particle density remains uniform and constant. For these particular conditions, this theory describes the nonequilibrium evolution of the static structure factor S(k;t) and of the dynamic properties, such as the self-intermediate scattering function F(S)(k,τ;t), where τ is the correlation delay time and t is the evolution or waiting time after the quench. Specific predictions are presented for the deepest quench (to zero temperature). The predicted evolution of the α-relaxation time τ(α)(t) as a function of t allows us to define the equilibration time t(eq)(φ), as the time after which τ(α)(t) has attained its equilibrium value τ(α)(eq)(φ). It is predicted that both, t(eq)(φ) and τ(α)(eq)(φ), diverge as φ→φ((a)), where φ((a)) is the hard-sphere dynamic-arrest volume fraction φ((a))(≈0.582), thus suggesting that the measurement of equilibrium properties at and above φ((a)) is experimentally impossible. The theory also predicts that for fixed finite waiting times t, the plot of τ(α)(t;φ) as a function of φexhibits two regimes, corresponding to samples that have fully equilibrated within this waiting time (φ≤φ((c))(t)), and to samples for which equilibration is not yet complete (φ≥φ((c))(t)). The crossover volume fraction φ((c))(t) increases with t but saturates to the value φ((a)).

Crossover from equilibration to aging: Nonequilibrium theory versus simulations

Understanding glasses and the glass transition requires comprehending the nature of the crossover from the ergodic (or equilibrium) regime, in which the stationary properties of the system have no history dependence, to the mysterious glass transition region, where the measured properties are nonstationary and depend on the protocol of preparation. In this work we use nonequilibrium molecular dynamics simulations to test the main features of the crossover predicted by the molecular version of the recently developed multicomponent nonequilibrium self-consistent generalized Langevin equation theory. According to this theory, the glass transition involves the abrupt passage from the ordinary pattern of full equilibration to the aging scenario characteristic of glass-forming liquids. The same theory explains that this abrupt transition will always be observed as a blurred crossover due to the unavoidable finiteness of the time window of any experimental observation. We find that within their finite waiting-time window, the simulations confirm the general trends predicted by the theory.

Communication: Probing the existence of partially arrested states in ionic liquids.

The recent predictions of the self-consistent generalized Langevin equation theory, describing the existence of unusual partially arrested states in the context of ionic liquids, were probed using all-atom molecular dynamics simulations of a room-temperature ionic liquid. We have found a slower diffusion of the smaller anions compared with the large cations for a wide range of temperatures. The arrest mechanism consists on the formation of a strongly repulsive glass by the anions, stabilized by the long range electrostatic potential. The diffusion of the less repulsive cations occurs through the holes left by the small particles. All of our observations in the simulated system coincide with the theoretical picture.

Glass–liquid–glass reentrance in mono-component colloidal dispersions

The self-consistent generalized Langevin equation (SCGLE) theory of colloid dynamics is employed to describe the ergodic-non-ergodic transition in model mono-disperse colloidal dispersions whose particles interact through hard-sphere plus short-ranged attractive forces. The ergodic-non-ergodic phase diagram in the temperature-concentration state space is determined for the hard-sphere plus attractive Yukawa model within the mean spherical approximation for the static structure factor by solving a remarkably simple equation for the localization length of the colloidal particles. Finite real values of this property signals non-ergodicity and determines the non-ergodic parameters f(k) and f(s)(k). The resulting phase diagram for this system, which involves the existence of reentrant (repulsive and attractive) glass states, is compared with the corresponding prediction of mode coupling theory. Although both theories coincide in the general features of this phase diagram, there are also clear qualitative differences. One of the most relevant is the SCGLE prediction that the ergodic-attractive glass transition does not preempt the gas-liquid phase transition, but always intersects the corresponding spinodal curve on its high-concentration side. We also calculate the ergodic-non-ergodic phase diagram for the sticky hard-sphere model to illustrate the dependence of the predicted SCGLE dynamic phase diagram on the choice of one important constituent element of the SCGLE theory.