Self-consistent Generalized Langevin Equation Research Articles

Arrested spinodal decomposition of the screened symmetric restricted primitive model.

Amorphous solids, such as glasses and gels, arise as the asymptotic limit of non-equilibrium and irreversible relaxation aging processes. These amorphous solids form when the system is suddenly and deeply quenched in the dynamic arrest region. We use the non-equilibrium self-consistent generalized Langevin equation (NE-SCGLE) theory to investigate the formation of such structures via arrested spinodal decomposition in the screened symmetric restricted primitive model. We propose a direct correlation function that allows us to derive an expression for the functional derivative of the chemical potential, which serves as the necessary input in the NE-SCGLE theory. By analyzing the asymptotic localization length and the asymptotic characteristic length, we identify different classes of dynamic arrest states as a function of the density and the final quench temperature. The system features simultaneously attractive and repulsive interactions, resulting in different arrested regions in the non-equilibrium phase diagram for a given screened parameter: (i) ionic glasses, (ii) electrostatic gels, and (iii) attractive glasses. Finally, by collapsing the asymptotic effective structure factors at a typical fractal dimension of df = 3, we confirm the formation of gels slightly above the glass-gel transition in the arrested phase diagram.

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.

Structural relaxation, dynamical arrest, and aging in soft-sphere liquids.

We investigate the structural relaxation of a soft-sphere liquid quenched isochorically (ϕ = 0.7) and instantaneously to different temperatures Tf above and below the glass transition. For this, we combine extensive Brownian dynamics simulations and theoretical calculations based on the non-equilibrium self-consistent generalized Langevin equation (NE-SCGLE) theory. The response of the liquid to a quench generally consists of a sub-linear increase of the α-relaxation time with system's age. Approaching the ideal glass-transition temperature from above (Tf > Ta), sub-aging appears as a transient process describing a broad equilibration crossover for quenches to nearly arrested states. This allows us to empirically determine an equilibration timescale teq(Tf) that becomes increasingly longer as Tf approaches Ta. For quenches inside the glass (Tf ≤ Ta), the growth rate of the structural relaxation time becomes progressively larger as Tf decreases and, unlike the equilibration scenario, τα remains evolving within the whole observation time-window. These features are consistently found in theory and simulations with remarkable semi-quantitative agreement and coincide with those revealed in a previous and complementary study [P. Mendoza-Méndez et al., Phys. Rev. 96, 022608 (2017)] that considered a sequence of quenches with fixed final temperature Tf = 0 but increasing ϕ toward the hard-sphere dynamical arrest volume fraction ϕHS a=0.582. The NE-SCGLE analysis, however, unveils various fundamental aspects of the glass transition, involving the abrupt passage from the ordinary equilibration scenario to the persistent aging effects that are characteristic of glass-forming liquids. The theory also explains that, within the time window of any experimental observation, this can only be observed as a continuous crossover.

General framework for the study of dynamical properties and arrested states of ionic liquids

In this work, we present a first-principles theoretical framework for the description of structural and dynamical behaviors of ionic liquids. This framework applies the self-consistent generalized Langevin equation (SCGLE) theory, which predicts the dynamically arrested states of several physical systems, including charge and size asymmetric ionic liquids. In its current form, the SCGLE is a particularly useful framework for predictions of arrested and partially arrested states in both classical and room-temperature molten salts, electrolytic solutions, and superionic conductors. Some illustrative examples of the predictive power of our framework are presented in this work, showing its potential application in design and development of novel conducting materials.

"Inner clocks" of glass-forming liquids.

Providing a physically sound explanation of aging phenomena in non-equilibrium amorphous materials is a challenging problem in modern statistical thermodynamics. The slow evolution of physical properties after quenches of control parameters is empirically well interpreted via the concept of material time (or internal clock) based on the Tool-Narayanaswamy-Moynihan model. Yet, the fundamental reasons of its striking success remain unclear. We propose a microscopic rationale behind the material time on the basis of the linear laws of irreversible thermodynamics and its extension that treats the corresponding kinetic coefficients as state functions of a slowly evolving material state. Our interpretation is based on the recognition that the same mathematical structure governs both the Tool model and the recently developed non-equilibrium extension of the self-consistent generalized Langevin equation theory, guided by the universal principles of Onsager's theory of irreversible processes. This identification opens the way for a generalization of the material-time concept to aging systems where several relaxation modes with very different equilibration processes must be considered, and partially frozen glasses manifest the appearance of partial ergodicity breaking and, hence, materials with multiple very distinct inner clocks.

On the dynamically arrested states of equilibrium and non-equilibrium gels: a comprehensive Brownian dynamics study

In this work a systematic study over a wide number of final thermodynamic states for two gel-forming liquids was performed. Such two kind of gel formers are distinguished by their specific interparticle interaction potential. We explored several thermodynamic states determining the thermodynamic, structural and dynamic properties of both liquids after a sudden temperature change. The thermodynamic analysis allows to identify that the liquid with short range attraction and long range repulsion lacks of a stable gas–liquid phase separation liquid, in contrast with the liquid with short range attractions. Thus, although for some thermodynamic states the structural behavior, measured by the static structure factor, is similar to and characteristic of the gel phase, for the short range attractive fluid the gel phase is a consequence of a spinodal decomposition process. In contrast, gelation in the short range attraction and long range repulsion liquid is not due to a phase separation. We also analyze the similarities and differences of the dynamic behavior of both systems through the analysis of the mean square displacement, the self part of the intermediate scattering function, the diffusion coefficient and the α relaxation time. Finally, using one of the main results of the non-equilibrium self-consistent generalized Langevin equation theory (NE-SCGLE), we determine the dynamic arrest phase diagram in the volume fraction and temperature (φ vs T) plane.

Spatially heterogeneous dynamics and locally arrested density fluctuations from first principles

We present a first-principles formalism for studying dynamical heterogeneities in glass-forming liquids. Based on the non-equilibrium self-consistent generalized Langevin equation theory, we were able to describe the time-dependent local density profile during the particle interchange among small regions of the fluid. The final form of the diffusion equation contains both the contribution of the chemical potential gradient written in terms of a coarse-grained density and a collective diffusion coefficient as well as the effect of a history-dependent mobility factor. With this diffusion equation, we captured interesting phenomena in glass-forming liquids such as the cases when a strong density gradient is accompanied by a very low mobility factor attributable to the denser part: in such circumstances, the density profile falls into an arrested state even in the presence of a density gradient. On the other hand, we also show that above a certain critical temperature, which depends on the volume fraction, any density heterogeneity relaxes to a uniform state in a finite time, known as equilibration time. We further show that such equilibration time varies little with the temperature in diluted systems but can change drastically with temperature in concentrated systems.

Ultra-slow and arrested density-fluctuations as precursor of spatial heterogeneity

Dynamical heterogeneities in glass-forming liquids subjected to cooling processes are studied by a theoretical framework based on the non-equilibrium self-consistent generalized Langevin equation theory. This theory predicts that slow cooling rates permit the relaxation to the equilibrium state distinguished by a homogeneous local density. In contrast, fast cooling rates provoke dynamically arrested density-fluctuations and the establishment of permanent spatial heterogeneities even in the presence of density gradients. We further show that the dynamics toward the arrested state has two steps: a truncated relaxation followed by a second relaxation of the diluted part of the system.

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 [2016 J. Phys. Chem. B 120 7975] 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 (at t 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.

On a fundamental description of the Kovacs' kinetic signatures in glass-forming systems.

The time-evolution equation for the time-dependent static structure factor of the non-equilibrium self-consistent generalized Langevin equation (NE-SCGLE) theory was used to investigate the kinetics of glass-forming systems under isochoric conditions. The kinetics are studied within the framework of the fictive temperature (TF) of the glassy structure. We solve for the kinetics of TF(t) and the time-dependent structure factor and find that they are different but closely related by a function that depends only on temperature. Furthermore, we are able to solve for the evolution of TF(t) in a set of temperature-jump histories referred to as the Kovacs' signatures. We demonstrate that the NE-SCGLE theory reproduces all the Kovacs' signatures, namely, intrinsic isotherm, asymmetry of approach, and memory effect. In addition, we extend the theory into largely unexplored, deep glassy state, regions that are below the notionally "ideal" glass temperature.

First-principles prediction of multiple stationary states in glass-forming liquids.

In the present work, the Non-Equilibrium Self-Consistent Generalized Langevin Equation (NESCGLE) theory is used to predict the final state of glass-forming liquids subjected to different cooling processes. We show that the NESCGLE theory correctly describes two essential features of the glass transition. Such features are the structural recovery and the dependence of the final state with the cooling rate. We demonstrate that below a particular temperature Tc, the system is unable to equilibrate, independently of the cooling rate. We show that the equilibrium state is only reached for the quasistatic process. Additionally, we show how, from the NESCGLE theory, it is possible to deduce a relaxation model of structural recovery, for which we obtain molecular expressions of the parameters.

Interference between the glass, gel, and gas-liquid transitions

Recent experiments and computer simulations have revealed intriguing phenomenological fingerprints of the interference between the ordinary equilibrium gas-liquid phase transition and the non-equilibrium glass and gel transitions. We thus now know, for example, that the liquid-gas spinodal line and the glass transition loci intersect at a finite temperature and density, that when the gel and the glass transitions meet, mechanisms for multistep relaxation emerge, and that the formation of gels exhibits puzzling latency effects. In this work we demonstrate that the kinetic perspective of the non-equilibrium self-consistent generalized Langevin equation (NE-SCGLE) theory of irreversible processes in liquids provides a unifying first-principles microscopic theoretical framework to describe these and other phenomena associated with spinodal decomposition, gelation, glass transition, and their combinations. The resulting scenario is in reality the competition between two kinetically limiting behaviors, associated with the two distinct dynamic arrest transitions in which the liquid-glass line is predicted to bifurcate at low densities, below its intersection with the spinodal line.

Glass-transition asymptotics in two theories of glassy dynamics: Self-consistent generalized Langevin equation and mode-coupling theory.

We contrast the generic features of structural relaxation close to the idealized glass transition that are predicted by the self-consistent generalized Langevin equation theory (SCGLE) against those that are predicted by the mode-coupling theory of the glass transition (MCT). We present an asymptotic solution close to conditions of kinetic arrest that is valid for both theories, despite the different starting points that are adopted in deriving them. This in particular provides the same level of understanding of the asymptotic dynamics in the SCGLE as was previously done only for MCT. We discuss similarities and different predictions of the two theories for kinetic arrest in standard glass-forming models, as exemplified through the hard-sphere system. Qualitative differences are found for models where a decoupling of relaxation modes is predicted, such as the generalized Gaussian core model, or binary hard-sphere mixtures of particles with very disparate sizes. These differences, which arise in the distinct treatment of the memory kernels associated to self- and collective motion of particles, lead to distinct scenarios that are predicted by each theory for partially arrested states and in the vicinity of higher-order glass-transition singularities.

The Subtle Kinetics of Arrested Spinodal Decomposition: Colloidal Gels and Porous Glasses

ABSTRACTThe non-equilibrium self-consistent generalized Langevin equation (NE-SCGLE) theory of irreversible processes in liquids has been proposed as a theoretical framework capable of predicting the age- and preparation-dependent properties of highly ubiquitous non-equilibrium amorphous solids, such as like glasses and gels. By this formalism, we discuss the main kinetic features of the irreversible relaxation of simple liquids involved in the arrested spinodal decomposition of suddenly and deeply quenched. At some lower temperature we identify, by means of a latency time within which particles retain a finite apparently stationary mobility, the crossover from full phase separation to arrested spinodal decomposition which leads to recognize the onset of gelation.

Nonequilibrium kinetics of the transformation of liquids into physical gels

A major stumbling block for statistical physics and materials science has been the lack of a uni- versal principle that allows us to understand and predict elementary structural, morphological, and dynamical properties of non-equilibrium amorphous states of matter. The recently-developed non- equilibrium self-consistent generalized Langevin equation (NE-SCGLE) theory, however, has been shown to provide a fundamental tool for the understanding of the most essential features of the transformation of liquids into amorphous solids, such as their aging kinetics or their dependence on the protocol of fabrication. In this work we focus on the predicted kinetics of one of the main fingerprints of the formation of gels by arrested spinodal decomposition of suddenly and deeply quenched simple liquids, namely, the arrest of structural parameters associated with the morpho- logical evolution from the initially uniform fluid, to the dynamically arrested sponge-like amorphous material. The comparison of the theoretical predictions (based on a simple specific model system), with simulation and experimental data measured on similar but more complex materials, suggests the universality of the predicted scenario.

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.

Self-consistent generalized Langevin equation theory of the dynamics of multicomponent atomic liquids

A fundamental challenge of the theory of liquids is to understand the similarities and differences in the macroscopic dynamics of both colloidal and atomic liquids, which originate in the (Newtonian or Brownian) nature of the microscopic motion of their constituents. Starting from the recently discovered long-time dynamic equivalence between a colloidal and an atomic liquid that share the same interparticle pair potential, in this work we develop a self-consistent generalized Langevin equation theory for the dynamics of equilibrium multicomponent atomic liquids, applicable as an approximate but quantitative theory describing the long-time diffusive dynamical properties of simple equilibrium atomic liquids. When complemented with a Gaussian-like approximation, this theory is also able to provide a reasonable representation of the passage from a ballistic to diffusive behavior. We illustrate the applicability of the resulting theory with three particular examples, namely, a monodisperse and a polydisperse monocomponent hard-sphere liquid and a highly size-asymmetric binary hard-sphere mixture. To assess the quantitative accuracy of our results, we perform event-driven molecular dynamics simulations, which corroborate the general features of the theoretical predictions.

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.

Equilibration and Aging of Liquids of Non-Spherically Interacting Particles.

The nonequilibrium self-consistent generalized Langevin equation theory of irreversible processes in liquids is extended to describe the positional and orientational thermal fluctuations of the instantaneous local concentration profile n(r,Ω,t) of a suddenly quenched colloidal liquid of particles interacting through nonspherically symmetric pairwise interactions, whose mean value n(r,Ω,t) is constrained to remain uniform and isotropic, n (r,Ω, t) = n (t). Such self-consistent theory is cast in terms of the time-evolution equation of the covariance [Formula: see text] of the fluctuations [Formula: see text] of the spherical harmonics projections nlm(k;t) of the Fourier transform of n(r,Ω,t). The resulting theory describes the nonequilibrium evolution after a sudden temperature quench of both, the static structure factor projections Slm(k,t) and the two-time correlation function [Formula: see text], where τ is the correlation delay time and t is the evolution or waiting time after the quench. As a concrete and illustrative application we use the resulting self-consistent equations to describe the irreversible processes of equilibration or aging of the orientational degrees of freedom of a system of strongly interacting classical dipoles with quenched positional disorder.