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.

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.

Anomalous dynamic arrest of non-interacting spheres ("polymer") diluted in a hard-sphere ("colloid") liquid.

Upon compression, the equilibrium hard-sphere liquid [pair potential uHS(r)] freezes at a packing fraction ϕf = 0.494 or, if crystallization is prevented, becomes metastable up to its glass transition at ϕg ≈ 0.58. Throughout the fluid regime (ϕ < ϕg), we are, thus, certain that this model liquid does not exhibit any form of kinetic arrest. If, however, a small portion of these spheres (packing fraction ϕ2 ≪ ϕ) happen to ignore each other [u22(r) = 0] but do not ignore the remaining "normal" hard spheres [u12(r) = u21(r) = u11(r) = uHS(r)], whose packing fraction is thus ϕ1 = ϕ - ϕ2, they run the risk of becoming dynamically arrested before they demix from the "normal" particles. This unexpected and counterintuitive scenario was first theoretically predicted and then confirmed by simulations.

Brownian motion in non-equilibrium systems and the Ornstein-Uhlenbeck stochastic process.

The Ornstein-Uhlenbeck stochastic process is an exact mathematical model providing accurate representations of many real dynamic processes in systems in a stationary state. When applied to the description of random motion of particles such as that of Brownian particles, it provides exact predictions coinciding with those of the Langevin equation but not restricted to systems in thermal equilibrium but only conditioned to be stationary. Here, we investigate experimentally single particle motion in a two-dimensional granular system in a stationary state, consisting of 1 mm stainless balls on a plane circular surface. The motion of the particles is produced by an alternating magnetic field applied perpendicular to the surface of the container. The mean square displacement of the particles is measured for a range of low concentrations and it is found that following an appropriate scaling of length and time, the short-time experimental curves conform a master curve covering the range of particle motion from ballistic to diffusive in accordance with the description of the Ornstein-Uhlenbeck model.

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.

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.

Transition from diffusive to subdiffusive motion in colloidal liquids.

In this work we report experimental and theoretical results for the motion of single colloidal particles embedded in complex fluids with different interparticle interactions. The motion of particles is found to follow a similar behavior for the different systems. In particular, the transition from the short-time diffusive motion to the subdiffusive intermediate-time motion is found to occur when the square root of its mean squared displacement is in the order of 1 tenth of the neighbors' interparticle distance, thus following a quantitative criterion similar to Lindemann's criterion for melting.

Generalized mean spherical approximation for the multicomponent restricted primitive model.

The generalized mean spherical approximation of Stell and Sun [J. Chem. Phys. 63, 5333 (1975)] for the binary charge-symmetric restricted primitive model (electroneutral mixture of equally sized hard spheres) is extended to charge-asymmetric binary electrolytes and to the generally multicomponent, but still restricted (i.e., equally sized) primitive model.

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.

Non-equilibrium relaxation and near-arrest dynamics in colloidal suspensions.

In this work we propose a theory to describe the irreversible diffusive relaxation of the local concentration of a colloidal dispersion that proceeds toward its stable thermodynamic equilibrium state, but which may in the process be trapped in metastable or dynamically arrested states. The central assumption of this theory is that the irreversible relaxation of the macroscopically observed mean value [Formula: see text] of the local concentration of colloidal particles is described by a diffusion equation involving a local mobility b(*)(r,t) that depends not only on the mean value [Formula: see text] but also on the covariance [Formula: see text] of the fluctuations [Formula: see text]. This diffusion equation must hence be solved simultaneously with the relaxation equation for the covariance σ(r,r';t), and here we also derive the corresponding relaxation equation. The dependence of the local mobility b(*)(r,t) on the mean value and the covariance is determined by a self-consistent set of equations involving now the spatially and temporally non-local time-dependent correlation functions, which in a uniform system in equilibrium reduces to the self-consistent generalized Langevin equation (SCGLE) theory of colloid dynamics. The resulting general theory considers the possibility that these relaxation processes occur under the influence of external fields, such as gravitational forces acting in the process of sedimentation. In this paper, however, we describe a simpler application, in which the system remains spatially uniform during the irreversible relaxation process, and discuss the general features of the glass transition scenario predicted by this non-equilibrium theory.

Alternative view of dynamic arrest in colloid-polymer mixtures.

Colloid-polymer mixtures are frequently viewed as an effective one-component fluid (the colloid) with polymer-mediated depletion interactions. This view, together with conventional mode coupling theory, constitutes the current description of the reentrant glass transition experimentally observed in these systems. A more fundamental view is to consider these systems as what they actually are, namely, genuine highly size-asymmetric binary colloidal mixtures. In this Letter we demonstrate that the recently developed multicomponent self-consistent generalized Langevin equation theory of dynamic arrest correctly predicts the observed reentrance in excellent quantitative agreement with the experimental glass transition line of a colloid-polymer mixture. In this scenario the polymer plays a much more active dynamic role than in the conventional one-component description.

Theory of dynamic arrest in colloidal mixtures.

We present a first-principles theory of dynamic arrest in colloidal mixtures based on the multicomponent self-consistent generalized Langevin equation theory of colloid dynamics [M. A. Chávez-Rojo and M. Medina-Noyola, Phys. Rev. E 72, 031107 (2005); M. A. Chávez-Rojo and M. Medina-Noyola, Phys. Rev. E76, 039902 (2007)]. We illustrate its application with a description of dynamic arrest in two simple model colloidal mixtures: namely, hard-sphere and repulsive Yukawa binary mixtures. Our results include observation of the two patterns of dynamic arrest, one in which both species become simultaneously arrested and the other involving the sequential arrest of the two species. The latter case gives rise to mixed states in which one species is arrested while the other species remains mobile. We also derive the ("bifurcation" or fixed-point") equations for the nonergodic parameters of the system, which takes the surprisingly simple form of a system of coupled equations for the localization length of the particles of each species. The solution of this system of equations indicates unambiguously which species is arrested (finite localization length) and which species remains ergodic (infinite localization length). As a result, we are able to draw the entire ergodic-nonergodic phase diagram of the binary hard-sphere mixture.

Diffusion of colloidal fluids in random porous media.

The diffusive relaxation of a colloidal fluid adsorbed in a porous medium depends on many factors, including the concentration and composition of the adsorbed colloidal fluid, the average structure of the porous matrix, and the nature of the colloid-colloid and colloid-substrate interactions. A simple manner to describe these effects is to model the porous medium as a set of spherical particles fixed in space at random positions with prescribed statistical structural properties. Within this model one may describe the relaxation of concentration fluctuations of the adsorbed fluid by simply setting to zero the short-time mobility of one species (the porous matrix) in a theory of the dynamics of equilibrium colloidal mixtures, or by extending such dynamic theory to explicitly consider the porous matrix as a random external field, as recently done in the framework of mode coupling theory [V. Krakoviack, Phys. Rev. Lett. 94, 065703 (2005)]. Here we consider the first approach and employ the self-consistent generalized Langevin equation (SCGLE) theory of the dynamics of equilibrium colloidal mixtures, to describe the dynamics of the mobile component. We focus on the short- and intermediate-time regimes, which we compare with Brownian dynamics simulations involving a binary mixture with screened Coulomb interactions for two models of the average static structure of the matrix: a porous matrix constructed by quenching configurations of an equilibrium mixture in which both species were first equilibrated together, and a preexisting matrix with prescribed average structure, in which we later add the mobile species. We conclude that in both cases, if the correct static structure factors are provided as input, the SCGLE theory correctly predicts the main features of the dynamics of the permeating fluid.

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.

Dynamic arrest within the self-consistent generalized Langevin equation of colloid dynamics.

This paper presents a recently developed theory of colloid dynamics as an alternative approach to the description of phenomena of dynamic arrest in monodisperse colloidal systems. Such theory, referred to as the self-consistent generalized Langevin equation (SCGLE) theory, was devised to describe the tracer and collective diffusion properties of colloidal dispersions in the short- and intermediate-time regimes. Its self-consistent character, however, introduces a nonlinear dynamic feedback, leading to the prediction of dynamic arrest in these systems, similar to that exhibited by the well-established mode coupling theory of the ideal glass transition. The full numerical solution of this self-consistent theory provides in principle a route to the location of the fluid-glass transition in the space of macroscopic parameters of the system, given the interparticle forces (i.e., a nonequilibrium analog of the statistical-thermodynamic prediction of an equilibrium phase diagram). In this paper we focus on the derivation from the same self-consistent theory of the more straightforward route to the location of the fluid-glass transition boundary, consisting of the equation for the nonergodic parameters, whose nonzero values are the signature of the glass state. This allows us to decide if a system, at given macroscopic conditions, is in an ergodic or in a dynamically arrested state, given the microscopic interactions, which enter only through the static structure factor. We present a selection of results that illustrate the concrete application of our theory to model colloidal systems. This involves the comparison of the predictions of our theory with available experimental data for the nonergodic parameters of model dispersions with hard-sphere and with screened Coulomb interactions.

State diagram for the electrostatic adsorption of charged colloids on confining walls: simulation and theory.

We study the structure of charged colloidal suspensions under confinement and determine a state diagram for the occurrence of electrostatic adsorption onto the confining walls, an effect that results in the accumulation of particles on the bounding surfaces and that could be relevant in experiments. We use Monte Carlo simulations to quantify this structural transition and perform theoretical calculations based on integral equations. Overall, our results provide a guide for experimentalists dealing with charged colloidal systems to determine the relevance of this purely electrostatic effect.

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.

Fundamental considerations on the mechanisms of silver cementation onto zinc particles in the Merril-Crowe process.

Studies on the Merrill-Crowe process as applied to silver recovery have shown that one half of the used zinc powder is wasted in water reduction at high cyanide concentrations, while the other half reduces silver ions from the cyanide solution. However, the cementation mechanisms as an electrochemical process taking place on the zinc surface do not explain the split of the electric current resulting from the anodic dissolution of zinc into two equal values. This study demonstrates that the mechanism for silver precipitation at high and low cyanide concentrations differs considerably. At low cyanide concentrations cementation is essentially an electrochemically-controlled process following a shrinking-core behavior. At high cyanide concentrations, the process seems not to be electrochemically controlled. The areas for zinc dissolution and silver deposition are not connected by an electrical-conducting medium and reduction of silver-cyano complex ions takes place by hydrogen adsorbed onto silver growing outward from the cementing zinc particles. The results are based on scanning electron microscopy of solids recovered from cementations in stirred reactors and in situ observations by optical microscopy of the cementation process on the edge of thin zinc disks in cyanide solutions.

Self-consistent generalized Langevin equation for colloid dynamics.

We present a general self-consistent theory of colloid dynamics which, for a system without hydrodynamic interactions, allows us to calculate F(k,t), and its self-diffusion counterpart F(S)(k,t), given the effective interaction pair potential u(r) between colloidal particles, and the corresponding equilibrium static structural properties. This theory is build upon the exact results for F(k,t) and F(S)(k,t) in terms of a hierarchy of memory functions, derived from the application of the generalized Langevin equation formalism, plus the proposal of Vineyard-like connections between F(k,t) and F(S)(k,t) through their respective memory functions, and a closure relation between these memory functions and the time-dependent friction function Delta zeta(t). As an illustrative application, we present and analyze a selection of numerical results of this theory in the short- and intermediate-time regimes, as applied to a two-dimensional repulsive Yukawa Brownian fluid. For this system, we find that our theory accurately describes the dynamic properties contained in F(k,t) in a wide range of conditions, including strongly correlated systems, at the longest times available from our computer simulations.

Long-time tracer diffusion of nonspherical brownian particles

The long-time tracer-diffusion properties of a nonspherical Brownian particle that interacts with a suspension of spherical particles are studied in terms of an idealized but nontrivial two-dimensional model system. For this system, the predictions of the generalized Langevin equation approach to tracer diffusion can be calculated, and compared with the (extrapolated) results of a computer simulation experiment. In the model, the nonspherical particle is represented by a rigid linear array of N(T) (=2 or 3) spherical particles with nearest-neighbor separation DeltaL. We calculate the long-time rotational and (transverse and longitudinal) translational diffusion coefficients. The theory is found to reproduce qualitatively and quantitatively the main features of the extrapolated results. Finally, we also present theoretical results that derive from still simpler approximate theoretical schemes.