Non-equilibrium dynamics of glass-forming liquid mixtures

Cooperative activated hopping dynamics in binary glass-forming liquids: effects of the size ratio, composition, and interparticle interactions

Slow dynamics in supercooled and glassy liquids is an important research topic in soft matter physics. Compared to the traditionally focused one-component systems, glassy dynamics in mixture systems adds in a rich set of new complexities, which are fundamentally interesting and also relevant for many technological applications. In this paper, we apply the recently developed self-consistent cooperative hopping theory (SCCHT) to systematically investigate the effects of the size ratio, composition and interparticle interactions on the cooperative activated hopping dynamics of matrix (in larger size) and penetrant (in smaller size) particles in varied binary sphere mixture model systems, with a specific focus on ultrahigh mixture packing fractions that mimic the deeply supercooled glass transition conditions for molecular/polymeric mixture materials. Analysis shows that in these high activation barrier cases, the long-range elastic distortion associated with a matrix particle hopping over its cage confinement always generates an elastic barrier of a nonnegligible magnitude, although the ratio between the elastic barrier and local barrier contribution is sensitively dependent on all three mixture-specific system factors considered in this work. SCCHT predicts two general scenarios of penetrant-matrix cooperative activated hopping dynamics: matrix/penetrant co-hopping (regime 1) or the penetrant mean barrier hopping time shorter than that of the matrix (regime 2). Increasing the penetrant-to-matrix size ratio or the penetrant-matrix cross-attraction strength is found to universally enlarge the composition window of regime 1. Diverse dynamical properties characterising different aspects of the cooperative activated hopping process, including the penetrant and matrix transient localization lengths, penetrant and matrix hopping jump distances, different types of local and elastic activated barriers, and matrix long-time diffusivity, relaxation time and dynamic fragility are quantitatively studied against a wide range of variations over the three system factors. Of particular interest is the universal "anti-plasticization" phenomenon achievable for sufficiently strong cross-attractive interactions. The prospects this work opens for the exploration of a wide variety of polymer-based mixture materials are briefly discussed at the end.

Non-equilibrium view of the amorphous solidification of liquids with competing interactions

The interplay between short-range attractions and long-range repulsions (SALR) characterizes the so-called liquids with competing interactions, which are known to exhibit a variety of equilibrium and non-equilibrium phases. The theoretical description of the phenomenology associated with glassy or gel states in these systems has to take into account both the presence of thermodynamic instabilities (such as those defining the spinodal line and the so called λ line) and the limited capability to describe genuine non-equilibrium processes from first principles. Here, we report the first application of the non-equilibrium self-consistent generalized Langevin equation theory to the description of the dynamical arrest processes that occur in SALR systems after being instantaneously quenched into a state point in the regions of thermodynamic instability. The physical scenario predicted by this theory reveals an amazing interplay between the thermodynamically driven instabilities, favoring equilibrium macro- and micro-phase separation, and the kinetic arrest mechanisms, favoring non-equilibrium amorphous solidification of the liquid into an unexpected variety of glass and gel states.

Phase separation and dynamical arrest of protein solutions dominated by short-range attractions

The interplay of liquid-liquid phase separation (LLPS) and dynamical arrest can lead to the formation of gels and glasses, which is relevant for such diverse fields as condensed matter physics, materials science, food engineering, and the pharmaceutical industry. In this context, protein solutions exhibit remarkable equilibrium and non-equilibrium behaviors. In the regime where attractive and repulsive forces compete, it has been demonstrated, for example, that the location of the dynamical arrest line seems to be independent of ionic strength, so that the arrest lines at different ionic screening lengths overlap, in contrast to the LLPS coexistence curves, which strongly depend on the salt concentration. In this work, we show that the same phenomenology can also be observed when the electrostatic repulsions are largely screened, and the range and strength of the attractions are varied. In particular, using lysozyme in brine as a model system, the metastable gas-liquid binodal and the dynamical arrest line as well as the second virial coefficient have been determined for various solution conditions by cloud-point measurements, optical microscopy, centrifugation experiments, and light scattering. With the aim of understanding this new experimental phenomenology, we apply the non-equilibrium self-consistent generalized Langevin equation theory to a simple model system with only excluded volume plus short-range attractions, to study the dependence of the predicted arrest lines on the range of the attractive interaction. The theoretical predictions find a good qualitative agreement with experiments when the range of the attraction is not too small compared with the size of the protein.

Non-equilibrium relaxation and aging in the dynamics of a dipolar fluid quenched towards the glass transition

The recently developed non-equilibrium self-consistent generalized Langevin equation theory of the dynamics of liquids of non-spherically interacting particles [2016J. Phys. Chem. B1207975] is applied to the description of the irreversible relaxation of a thermally and mechanically quenched dipolar fluid. Specifically, we consider a dipolar hard-sphere liquid quenched (att_w= 0) from full equilibrium conditions towards different ergodic-non-ergodic transitions. Qualitatively different scenarios are predicted by the theory for the time evolution of the system after the quench (t_w> 0), that depend on both the kind of transition approached and the specific features of the protocol of preparation. Each of these scenarios is characterized by the kinetics displayed by a set of structural correlations, and also by the development of two characteristic times describing the relaxation of the translational and rotational dynamics, allowing us to highlight the crossover from equilibration to aging in the system and leading to the prediction of different underlying mechanisms and relaxation laws for the dynamics at each of the glass transitions explored.

Waiting-time dependent non-equilibrium phase diagram of simple glass- and gel-forming liquids

Under numerous circumstances, many soft and hard materials are present in a puzzling wealth of non-equilibrium amorphous states, whose properties are not stationary and depend on preparation. They are often summarized in unconventional "phase diagrams" that exhibit new "phases" and/or "transitions" in which time, however, is an essential variable. This work proposes a solution to the problem of theoretically defining and predicting these non-equilibrium phases and their time-evolving phase diagrams, given the underlying molecular interactions. We demonstrate that these non-equilibrium phases and the corresponding non-stationary (i.e., aging) phase diagrams can indeed be defined and predicted using the kinetic perspective of a novel non-equilibrium statistical mechanical theory of irreversible processes. This is illustrated with the theoretical description of the transient process of dynamic arrest into non-equilibrium amorphous solid phases of an instantaneously quenched simple model fluid involving repulsive hard-sphere plus attractive square well pair interactions.

Glassy dynamics in asymmetric binary mixtures of hard spheres

We perform a systematic and detailed study of the glass transition in highly asymmetric binary mixtures of colloidal hard spheres, combining differential dynamic microscopy experiments, event-driven molecular dynamics simulations, and theoretical calculations, exploring the whole state diagram and determining the self-dynamics and collective dynamics of both species. Two distinct glassy states involving different dynamical arrest transitions are consistently described, namely, a double glass with the simultaneous arrest of the self-dynamics and collective dynamics of both species, and a single glass of large particles in which the self-dynamics of the small species remains ergodic. In the single-glass scenario, spatial modulations in the collective dynamics of both species occur due to the structure of the large spheres, a feature not observed in the double-glass domain. The theoretical results, obtained within the self-consistent generalized Langevin equation formalism, are in agreement with both simulations and experimental data, thus providing a stringent validation of this theoretical framework in the description of dynamical arrest in highly asymmetric mixtures. Our findings are summarized in a state diagram that classifies the various amorphous states of highly asymmetric mixtures by their dynamical arrest mechanisms.

Physical aging and compressed exponential behaviors in a model soft colloidal system

Diffusing wave spectroscopy (DWS)-based micro-rheology has been used in different optical geometries (backscattering and transmission) as well as different sample thicknesses in order to probe system dynamics at different length scales [D. J. Pine, D. A. Weitz, J. X. Zhu, E. Herbolzheimer. J. Phys., 1990, 51(18), 2101-2127]. Previous study from this lab [Q. Li, X. Peng, G. B. McKenna. Soft Matter, 2017, 13(7), 1396-1404] indicates the DWS-based micro-rheology observes the system non-equilibrium behaviors differently from macro-rheology. The object of the present work was to further explore the non-equilibrium dynamics and to address the range of utility of DWS as a micro-rheological method. A thermo-sensitive core-shell colloidal system was investigated both during aging and subsequent to aging into a metastable equilibrium state using temperature-jump induced volume fraction-jump experiments. We find that in the non-equilibrium state, significant differences in the measured dynamics are observed for the different geometries and length scales. Compressed exponential relaxations for the autocorrelation function g2(t) were observed for large length scales. However, upon converting the g2(t) data to the mean square displacement (MSD), such differences with length scale diminished and the long-time MSD behavior was consistent with diffusive behavior. These observations in the non-equilibrium behaviors for different length scales leads to questioning of some interpretations in the current field of light scattering-based micro-rheology and provides a possibility to interrogate the aging mechanisms in colloidal glasses from a broader perspective than normally considered in measurements of g2(t) using DWS-based micro-rheology.

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.