General nonequilibrium theory of colloid dynamics

A simplified longitudinal model for the development of Type 2 Diabetes Mellitus

Obesity and diabetes are a progressively more and more deleterious hallmark of modern, well fed societies. In order to study the potential impact of strategies designed to obviate the pathological consequences of detrimental lifestyles, a model for the development of Type 2 diabetes geared towards large population simulations would be useful. The present work introduces such a model, representing in simplified fashion the interplay between average glycemia, average insulinemia and functional beta-cell mass, and incorporating the effects of excess food intake or, conversely, of physical activity levels. Qualitative properties of the model are formally established and simulations are shown as examples of its use.

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 T_f 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 (T_f > T^a), 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 t^eq(T_f) that becomes increasingly longer as T_f approaches T^a. For quenches inside the glass (T_f ≤ T^a), the growth rate of the structural relaxation time becomes progressively larger as T_f 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 T_f = 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.

Perspective: New directions in dynamical density functional theory

Classical dynamical density functional theory (DDFT) has become one of the central modeling approaches in nonequilibrium soft matter physics. Recent years have seen the emergence of novel and interesting fields of application for DDFT. In particular, there has been a remarkable growth in the amount of work related to chemistry. Moreover, DDFT has stimulated research on other theories such as phase field crystal models and power functional theory. In this perspective, we summarize the latest developments in the field of DDFT and discuss a variety of possible directions for future research.

From Subaging to Hyperaging in Structural Glasses

We demonstrate nonequilibrium scaling laws for the aging and equilibration dynamics in glass formers that emerge from combining a relaxation equation for the static structure with the equilibrium scaling laws of glassy dynamics. Different scaling regimes are predicted for the evolution of the structural relaxation time τ with age (waiting time t_{w}), depending on the depth of the quench from the liquid into the glass: "simple" aging (τ∼t_{w}) applies for quenches close to the critical point of mode-coupling theory (MCT) and implies "subaging" (τ≈t_{w}^{δ} with δ<1) as a broad equilibration crossover for quenches to nearly arrested equilibrium states; "hyperaging" (or superaging, τ∼t_{w}^{δ^{'}} with δ^{'}>1) emerges for quenches deep into the glass. The latter is cut off by non-mean-field fluctuations that we account for within a recent extension of MCT, the stochastic β-relaxation theory (SBR). We exemplify the scaling laws with a schematic model that quantitatively fits simulation data.

"Inner clocks" of glass-forming liquids

Providing a physically sound explanation of aging phenomena in non-equilibrium amorphous materialsis a challenging problem in modern statistical thermodynamics. The slow evolution of physical propertiesafter quenches of control parameters is empirically well interpreted via the concept of material time (orinternal clock), based on the Tool-Narayanaswamy-Moynihan (TNM) model. Yet, the fundamental reasonsof its striking success remain unclear. We propose a microscopic rationale behind the material time onthe basis of the linear laws of irreversible thermodynamics and its extension that treats the correspondingkinetic coefficients as state functions of a slowly evolving material state. Our interpretation is based onthe recognition that the same mathematical structure governs both the Tool model and the recently devel-oped non-equilibrium extension of the self-consistent generalized Langevin equation theory (NE-SCGLE),guided by the universal principles of Onsager's theory of irreversible processes. This identification opensthe way for a generalization of the material-time concept to aging systems where several relaxation modeswith very different equilibration processes must be considered, and partially frozen glasses manifest theappearance 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 (φvsT) plane.

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.

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 (T_F) of the glassy structure. We solve for the kinetics of T_F(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 T_F(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.

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.

Vitrification is a spontaneous non-equilibrium transition driven by osmotic pressure

Persistent dynamics in colloidal glasses suggest the existence of a non-equilibrium driving force for structural relaxation during glassy aging. But the implicit assumption in the literature that colloidal glasses form within the metastable state bypasses the search for a driving force for vitrification and glassy aging and its connection with a metastable state. The natural relation of osmotic pressure to number-density gradients motivates us to investigate the osmotic pressure as this driving force. We use dynamic simulation to quench a polydisperse hard-sphere colloidal liquid into the putative glass region while monitoring structural relaxation and osmotic pressure. Following quenches to various depths in volume fractionϕ(whereϕ_RCP≈ 0.678 for 7% polydispersity), the osmotic pressure overshoots its metastable value, then decreases with age toward the metastable pressure, driving redistribution of coordination number and interparticle voids that smooths structural heterogeneity with age. For quenches to 0.56 ⩽ϕ⩽ 0.58, accessible post-quench volume redistributes with age, allowing the glass to relax into a strong supercooled liquid and easily reach a metastable state. At higher volume fractions, 0.59 ⩽ϕ< 0.64, this redistribution encounters a barrier that is subsequently overcome by osmotic pressure, allowing the system to relax toward the metastable state. But forϕ⩾ 0.64, the overshoot is small compared to the high metastable pressure; redistribution of volume stops as particles acquire contacts and get stuck, freezing the system far from the metastable state. Overall, the osmotic pressure drives structural rearrangements responsible for both vitrification and glassy age-relaxation. The connection of energy, pressure, and structure identifies the glass transition, 0.63 <ϕ_g⩽ 0.64. We leverage the connection of osmotic pressure to energy density to put forth the mechanistic view that relaxation of structural heterogeneity in colloidal glasses occurs via individual particle motion driven by osmotic pressure, and is a spontaneous energy minimization process that drives the glass off and back to the metastable state.

"Dense diffusion" in colloidal glasses: short-ranged long-time self-diffusion as a mechanistic model for relaxation dynamics

Despite decades of exploration of the colloidal glass transition, mechanistic explanation of glassy relaxation processes has remained murky. State-of-the-art theoretical models of the colloidal glass transition such as random first order transition theory, active barrier hopping theory, and non-equilibrium self-consistent generalized Langevin theory assert that relaxation reported at volume fractions above the ideal mode coupling theory prediction φg,MCT requires some sort of activated process, and that cooperative motion plays a central role. However, discrepancies between predicted and measured values of φg and ambiguity in the role of cooperative dynamics persist. Underlying both issues is the challenge of conducting deep concentration quenches without flow and the difficulty in accessing particle-scale dynamics. These two challenges have led to widespread use of fitting methods to identify divergence, but most a priori assume divergent behavior; and without access to detailed particle dynamics, it is challenging to produce evidence of collective dynamics. We address these limitations by conducting dynamic simulations accompanied by experiments to quench a colloidal liquid into the putative glass by triggering an increase in particle size, and thus volume fraction, at constant particle number density. Quenches are performed from the liquid to final volume fractions 0.56 ≤ φ ≤ 0.63. The glass is allowed to age for long times, and relaxation dynamics are monitored throughout the simulation. Overall, correlated motion acts to release dynamics from the glassy plateau - but only over length scales much smaller than a particle size - allowing self-diffusion to re-emerge; self-diffusion then relaxes the glass into an intransient diffusive state, which persists for φ < 0.60. We observe similar relaxation dynamics up to φ = 0.63 before achieving the intransient state. We find that this long-time self-diffusion is short-ranged: analysis of mean-square displacement reveals a glassy cage size a fraction of a particle size that shrinks with quench depth, i.e. increasing volume fraction. Thus the equivalence between cage size and particle size found in the liquid breaks down in the glass, which we confirm by examining the self-intermediate scattering function over a range of wave numbers. The colloidal glass transition can hence be viewed mechanistically as a shift in the long-time self-diffusion from long-ranged to short-ranged exploration of configurations. This shift takes place without diverging dynamics: there is a smooth transition as particle mobility decreases dramatically with concomitant emergence of a dense local configuration space that permits sampling of many configurations via local particle motion.

Active glasses

Active glassy matter has recently emerged as a novel class of non-equilibrium soft matter, combining energy-driven, active particle movement with dense and disordered glass-like behavior. Here we review the state-of-the-art in this field from an experimental, numerical, and theoretical perspective. We consider both non-living and living active glassy systems, and discuss how several hallmarks of glassy dynamics (dynamical slowdown, fragility, dynamical heterogeneity, violation of the Stokes-Einstein relation, and aging) are manifested in such materials. We start by reviewing the recent experimental evidence in this area of research, followed by an overview of the main numerical simulation studies and physical theories of active glassy matter. We conclude by outlining several open questions and possible directions for future work.

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.

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.

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.

Different routes into the glass state for soft thermo-sensitive colloids

We report an experimental and theoretical investigation of glass formation in soft thermo-sensitive colloids following two different routes: a gradual increase of the particle number density at constant temperature and an increase of the radius in a fixed volume at constant particle number density. Confocal microscopy experiments and the non-equilibrium self-consistent generalized Langevin equation (NE-SCGLE) theory consistently show that the two routes lead to a dynamically comparable state at sufficiently long aging times. However, experiments reveal the presence of moderate but persistent structural differences. Successive cycles of radius decrease and increase lead instead to a reproducible glass state, indicating a suitable route to obtain rejuvenation without using shear fields.

Strong attractions and repulsions mediated by monovalent salts

Controlling interactions between proteins and nanoparticles in electrolyte solutions is crucial for advancing biological sciences and biotechnology. The assembly of charged nanoparticles (NPs) and proteins in aqueous solutions can be directed by modifying the salt concentration. High concentrations of monovalent salt can induce the solubilization or crystallization of NPs and proteins. By using a multiscale coarse-grained molecular dynamics approach, we show that, due to ionic correlations in the electrolyte, NPs pairs at high monovalent salt concentrations interact via remarkably strong long-range attractions or repulsions, which can be split into three regimes depending on the surface charge densities of the NPs. NPs with zero-to-low surface charge densities interact via a long-range attraction that is stronger and has a similar range to the depletion attraction induced by polymers with radius of gyrations comparable to the NP diameter. On the other hand, moderately charged NPs with smooth surfaces as well as DNA-functionalized NPs with no possibility of hybridization between them interact via a strong repulsion of range and strength larger than the repulsion predicted by models that neglect ionic correlations, including the Derjaguin-Landau-Vervey-Overbeek (DLVO) model. Interactions between strongly charged NPs (>2 e/nm²), both types smooth and DNA-functionalized NPs, show an attractive potential well at intermediate-to-high salt concentrations, which demonstrates that electrolytes can induce aggregation of strongly charged NPs. Our work provides an improved understanding of the role of ionic correlations in NP assembly and design rules to utilize the salting-out process to crystallize NPs.

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.

Renormalized dynamics of the Dean-Kawasaki model

We study the model of a supercooled liquid for which the equation of motion for the coarse-grained density ρ(x,t) is the nonlinear diffusion equation originally proposed by Dean and Kawasaki, respectively, for Brownian and Newtonian dynamics of fluid particles. Using a Martin-Siggia-Rose (MSR) field theory we study the renormalization of the dynamics in a self-consistent form in terms of the so-called self-energy matrix Σ. The appropriate model for the renormalized dynamics involves an extended set of field variables {ρ,θ}, linked through a nonlinear constraint. The latter incorporates, in a nonperturbative manner, the effects of an infinite number of density nonlinearities in the dynamics. We show that the contributing element of Σ which renormalizes the bare diffusion constant D(0) to D(R) is same as that proposed by Kawasaki and Miyazima [Z. Phys. B Condens. Matter 103, 423 (1997)]. D(R) sharply decreases with increasing density. We consider the likelihood of a ergodic-nonergodic (ENE) transition in the model beyond a critical point. The transition is characterized by the long-time limit of the density correlation freezing at a nonzero value. From our analysis we identify an element of Σ which arises from the above-mentioned nonlinear constraint and is key to the viability of the ENE transition. If this self-energy would be zero, then the model supports a sharp ENE transition with D(R)=0 as predicted by Kawasaki and Miyazima. With the full model having nonzero value for this self-energy, the density autocorrelation function decays to zero in the long-time limit. Hence the ENE transition is not supported in the model.

Bond lifetime and diffusion coefficient in colloids with short-range interactions

We use molecular dynamics simulations to study the influence of short-range structures in the interaction potential between hard-sphere-like colloidal particles. Starting from model potentials and effective potentials in binary mixtures computed from the Ornstein-Zernike equations, we investigate the influence of the range and strength of a possible tail beyond the usual core repulsion or the presence of repulsive barriers. The diffusion coefficient and mean "bond" lifetimes are used as indicators of the effect of this structure on the dynamics. The existence of correlations between the variations of these quantities with the physical parameters is discussed to assess the interpretation of dynamics slowing down in terms of long-lived bonds. We also discuss the question of a universal behaviour determined by the second virial coefficient B ((2)) and the interplay of attraction and repulsion. While the diffusion coefficient follows the B ((2)) law for purely attractive tails, this is no longer true in the presence of repulsive barriers. Furthermore, the bond lifetime shows a dependence on the physical parameters that differs from that of the diffusion coefficient. This raises the question of the precise role of bonds on the dynamics slowing down in colloidal gels.

Self-consistent generalized Langevin-equation theory for liquids of nonspherically interacting particles

A self-consistent generalized Langevin-equation theory is proposed to describe the self- and collective dynamics of a liquid of linear Brownian particles. The equations of motion for the spherical harmonics projections of the collective and self-intermediate-scattering functions, F_{lm,lm}(k,t) and F_{lm,lm}^{S}(k,t), are derived as a contraction of the description involving the stochastic equations of the corresponding tensorial one-particle density n_{lm}(k,t) and the translational (α=T) and rotational (α=R) current densities j_{lm}^{α}(k,t). Similar to the spherical case, these dynamic equations require as an external input the equilibrium structural properties of the system contained in the projections of the static structure factor, denoted by S_{lm,lm}(k). Complementing these exact equations with simple (Vineyard-like) approximate relations for the collective and the self-memory functions we propose a closed self-consistent set of equations for the dynamic properties involved. In the long-time asymptotic limit, these equations become the so-called bifurcation equations, whose solutions (the nonergodicity parameters) can be written, extending the spherical case, in terms of one translational and one orientational scalar dynamic order parameter, γ_{T} and γ_{R}, which characterize the possible dynamical arrest transitions of the system. As a concrete illustrative application of this theory we determine the dynamic arrest diagram of the dipolar hard-sphere fluid. In qualitative agreement with mode coupling theory, the present self-consistent equations also predict three different regions in the state space spanned by the macroscopic control parameters η (volume fraction) and T* (scaled temperature): a region of fully ergodic states, a region of mixed states, in which the translational degrees of freedom become arrested while the orientational degrees of freedom remain ergodic, and a region of fully nonergodic states.

Dynamic equivalences in the hard-sphere dynamic universality class

We perform systematic simulation experiments on model systems with soft-sphere repulsive interactions to test the predicted dynamic equivalence between soft-sphere liquids with similar static structure. For this we compare the simulated dynamics (mean squared displacement, intermediate scattering function, α-relaxation time, etc.) of different soft-sphere systems, between them and with the hard-sphere liquid. We then show that the referred dynamic equivalence does not depend on the (Newtonian or Brownian) nature of the microscopic laws of motion of the constituent particles, and hence, applies independently to colloidal and to atomic simple liquids. Finally, we verify another more recently proposed dynamic equivalence, this time between the long-time dynamics of an atomic liquid and its corresponding Brownian fluid (i.e., the Brownian system with the same interaction potential).

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)).

Tunable soft structure in charged fluids confined by dielectric interfaces

Fluids of charged particles act as the supporting medium for chemical reactions and physical, dynamical, and biological processes. The local structure in an electrolytic background is deformed by micro- and nanoscopic polarizable objects. Vice versa, the forces between the objects are regulated by the cohesive properties of the background. We study here the range and strength of these forces and the microscopic origin from which they emerge. We find the forces to be sensitively dependent on the material properties of the charged fluid and the immersed solutes. The induced interactions can be varied over decades, offering high tunability and aided by accurate theory, control in experiments and applications. To distinguish correlational effects from simple ionic screening, we describe electrolyte-induced forces between neutral objects. The interplay of thermal motion, short-range repulsions, and electrostatic forces is responsible for a soft structure in the fluid. This structure changes near polarizable interfaces and causes diverse attractions between confining walls that seem well-exploited by microbiological systems. For parameters that correspond to monovalent electrolytes in biologically and technologically relevant aqueous environments, we find induced forces between nanoscopic areas of the order of piconewtons over a few nanometers.

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.

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.

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.