Cooperative diffusion in colloidal mixtures

Rescaled mode-coupling scheme for the quantitative description of experimentally observed colloid dynamics

We describe experimentally observed collective dynamics in colloidal suspensions of model hard-sphere particles using a modified mode coupling theory (MCT). This rescaled MCT is capable of describing quantitatively the wave-vector and time-dependent diffusion in these systems. Intermediate scattering functions of liquidlike structured dispersions are determined by means of static and dynamic light-scattering experiments. The structure and short-time dynamics of the systems can be described quantitatively employing a multicomponent Percus-Yevick ansatz for the partial structure factors and an effective, one-component description of hydrodynamic interactions based on the semianalytical δγ expansion. Combined with a recently proposed empirical modification of MCT in which memory functions are calculated using effective structure factors at rescaled number densities, the scheme is able to model the collective dynamics over the entire accessible time and wave-vector range and predicts the volume-fraction-dependence of long-time self-diffusion coefficients and the zero-shear viscosity quantitatively. This highlights the potential of MCT as a practical tool for the quantitative analysis and prediction of experimental observations.

Microscopic theory for nonequilibrium correlation functions in dense active fluids

One of the key hallmarks of dense active matter in the liquid, supercooled, and solid phases is the so-called equal-time velocity correlations. Crucially, these correlations can emerge spontaneously, i.e., they require no explicit alignment interactions, and therefore represent a generic feature of dense active matter. This indicates that for a meaningful comparison or possible mapping between active and passive liquids one not only needs to understand their structural properties, but also the impact of these velocity correlations. This has already prompted several simulation and theoretical studies, though they are mostly focused on athermal systems and thus overlook the effect of translational diffusion. Here, we present a fully microscopic method to calculate nonequilibrium correlations in two-dimensional systems of thermal active Brownian particles (ABPs). We use the integration through transients formalism together with (active) mode-coupling theory and analytically calculate qualitatively consistent static structure factors and active velocity correlations. We complement our theoretical results with simulations of both thermal and athermal ABPs which exemplify the disruptive role that thermal noise has on velocity correlations.

A deep learning approach to the measurement of long-lived memory kernels from generalized Langevin dynamics

Memory effects are ubiquitous in a wide variety of complex physical phenomena, ranging from glassy dynamics and metamaterials to climate models. The Generalized Langevin Equation (GLE) provides a rigorous way to describe memory effects via the so-called memory kernel in an integro-differential equation. However, the memory kernel is often unknown, and accurately predicting or measuring it via, e.g., a numerical inverse Laplace transform remains a herculean task. Here, we describe a novel method using deep neural networks (DNNs) to measure memory kernels from dynamical data. As a proof-of-principle, we focus on the notoriously long-lived memory effects of glass-forming systems, which have proved a major challenge to existing methods. In particular, we learn the operator mapping dynamics to memory kernels from a training set generated with the Mode-Coupling Theory (MCT) of hard spheres. Our DNNs are remarkably robust against noise, in contrast to conventional techniques. Furthermore, we demonstrate that a network trained on data generated from analytic theory (hard-sphere MCT) generalizes well to data from simulations of a different system (Brownian Weeks-Chandler-Andersen particles). Finally, we train a network on a set of phenomenological kernels and demonstrate its effectiveness in generalizing to both unseen phenomenological examples and supercooled hard-sphere MCT data. We provide a general pipeline, KernelLearner, for training networks to extract memory kernels from any non-Markovian system described by a GLE. The success of our DNN method applied to noisy glassy systems suggests that deep learning can play an important role in the study of dynamical systems with memory.

Active glassy dynamics is unaffected by the microscopic details of self-propulsion

Recent years have seen a rapid increase of interest in dense active materials, which, in the disordered state, share striking similarities with the conventional passive glass-forming matter. For such passive glassy materials, it is well established (at least in three dimensions) that the details of the microscopic dynamics, e.g., Newtonian or Brownian, do not influence the long-time glassy behavior. Here, we investigate whether this still holds true in the non-equilibrium active case by considering two simple and widely used active particle models, i.e., active Ornstein-Uhlenbeck particles (AOUPs) and active Brownian particles (ABPs). In particular, we seek to gain more insight into the role of the self-propulsion mechanism on the glassy dynamics by deriving a mode-coupling theory (MCT) for thermal AOUPs, which can be directly compared to a recently developed MCT for ABPs. Both theories explicitly take into account the active degrees of freedom. We solve the AOUP- and ABP-MCT equations in two dimensions and demonstrate that both models give almost identical results for the intermediate scattering function over a large variety of control parameters (packing fractions, active speeds, and persistence times). We also confirm this theoretical equivalence between the different self-propulsion mechanisms numerically via simulations of a polydisperse mixture of active quasi-hard spheres, thereby establishing that, at least for these model systems, the microscopic details of self-propulsion do not alter the active glassy behavior.

Generalized mode-coupling theory for mixtures of Brownian particles

Generalized mode-coupling theory (GMCT) has recently emerged as a promising first-principles theory to study the poorly understood dynamics of glass-forming materials. Formulated as a hierarchical extension of standard mode-coupling theory (MCT), it is able to systematically improve its predictions by including the exact dynamics of higher-order correlation functions into its hierarchy. However, in contrast to Newtonian dynamics, a fully generalized version of the theory based on Brownian dynamics is still lacking. To close this gap, we provide a detailed derivation of GMCT for colloidal mixtures obeying a many-body Smoluchowski equation. We demonstrate that a hierarchy of coupled equations can again be established and show that these, consistent with standard MCT, are identical to the ones obtained from Newtonian GMCT when taking the overdamped limit. Consequently, the nontrivial similarity between Brownian and Newtonian MCT is maintained for our multicomponent GMCT. As a proof of principle, we also solve the generalized mode-coupling equations for the binary Kob-Andersen Lennard-Jones mixture undergoing Brownian dynamics and confirm the improved predictive power of the theory upon using more levels of the GMCT hierarchy of equations.

Multi-component generalized mode-coupling theory: predicting dynamics from structure in glassy mixtures

The emergence of glassy dynamics and the glass transition in dense disordered systems is still not fully understood theoretically. Mode-coupling theory (MCT) has shown to be effective in describing some of the non-trivial features of glass formation, but it cannot explain the full glassy phenomenology due to the strong approximations on which it is based. Generalized mode-coupling theory (GMCT) is a hierarchical extension of the theory, which is able to outclass MCT by carefully describing the dynamics of higher-order correlations in its generalized framework. Unfortunately, the theory has so far only been developed for single-component systems and as a result works poorly for highly polydisperse materials. In this paper, we solve this problem by developing GMCT for multi-component systems. We use it to predict the glassy dynamics of the binary Kob-Andersen Lennard-Jones mixture, as well as its purely repulsive Weeks-Chandler-Andersen analogue. Our results show that each additional level of the GMCT hierarchy gradually improves the predictive power of GMCT beyond its previous limit. This implies that our theory is able to harvest more information from the static correlations, thus being able to better understand the role of attraction in supercooled liquids from a first-principles perspective.

Assessment of the Wolf method using the Stillinger-Lovett sum rules: From strong electrolytes to weakly charged colloidal dispersions

The Ewald method has been the cornerstone in molecular simulations for modeling electrostatic interactions of charge-stabilized many-body systems. In the late 1990s, Wolf and collaborators developed an alternative route to describe the long-range nature of electrostatic interactions; from a computational perspective, this method provides a more efficient and straightforward way to implement long-range electrostatic interactions than the Ewald method. Despite these advantages, the validity of the Wolf potential to account for the electrostatic contribution in charged fluids remains controversial. To alleviate this situation, in this contribution, we implement the Wolf summation method to both electrolyte solutions and charged colloids with moderate size and charge asymmetries in order to assess the accuracy and validity of the method. To this end, we verify that the proper selection of parameters within the Wolf method leads to results that are in good agreement with those obtained through the standard Ewald method and the theory of integral equations of simple liquids within the so-called hypernetted chain approximation. Furthermore, we show that the results obtained with the original Wolf method do satisfy the moment conditions described by the Stillinger-Lovett sum rules, which are directly related to the local electroneutrality condition and the electrostatic screening in the Debye-Hückel regime. Hence, the fact that the solution provided by the Wolf method satisfies the first and second moments of Stillinger-Lovett proves, for the first time, the reliability of the method to correctly incorporate the electrostatic contribution in charge-stabilized fluids. This makes the Wolf method a powerful alternative compared to more demanding computational approaches.

Glass transition and reversible gelation in asymmetric binary mixtures: A study by mode coupling theory and molecular dynamics

The glass transition and the binodals of asymmetric binary mixtures are investigated from the effective fluid approach in the mode coupling theory and by molecular dynamics. Motivated by previous theoretical predictions, the hard-sphere mixture and the Asakura-Oosawa models are used to analyze experimental results from the literature, relative to polystyrene spheres mixed either with linear polymers or with dense microgel particles. In agreement with the experimental observations, the specificity of the depletant particles is shown to favor lower density gels. It further favors equilibrium gelation by reducing also the tendency of the system to phase separate. These results are confirmed by a phenomenological modification of the mode coupling theory in which the vertex functions are computed at an effective density lower than the actual one. A model effective potential in asymmetric mixtures of hard particles is used to further check this phenomenological modification against molecular dynamics simulation.

α-relaxation, shear viscosity, and elastic moduli of hard-particle fluids from a mode-coupling theory with a retarded vertex

The recently proposed modification of the mode-coupling theory (MCT) in which the static structure used in the vertex is computed at a lower density than the actual one is tested on several dynamics-related properties. The predictions from this modified version of MCT calibrated on the one-component hard-sphere fluid are found in very good agreement with simulation data for one-component and binary hard-sphere fluids. They are also relevant for the stress moduli for models with attractive tails beyond the hard core. The clear improvement observed on several properties should give a new impetus to the use of MCT as a quantitative tool.

Role of the Pair Correlation Function in the Dynamical Transition Predicted by Mode Coupling Theory

In a recent study, we have found that for a large number of systems the configurational entropy at the pair level S_{c2}, which is primarily determined by the pair correlation function, vanishes at the dynamical transition temperature T_{c}. Thus, it appears that the information of the transition temperature is embedded in the structure of the liquid. In order to investigate this, we describe the dynamics of the system at the mean field level and, using the concepts of the dynamical density functional theory, show that the dynamical transition temperature depends only on the pair correlation function. Thus, this theory is similar in spirit to the microscopic mode coupling theory (MCT). However, unlike microscopic MCT, which predicts a very high transition temperature, the present theory predicts a transition temperature that is similar to T_{c}. This implies that the information of the dynamical transition temperature is embedded in the pair correlation function.

Stratification and Size Segregation of Ternary and Polydisperse Colloidal Suspensions during Drying

We investigate the drying process of three-component and polydisperse colloidal suspensions using Brownian dynamics simulations. We have previously reported (Phys. Rev. Lett. 2016, 116, 118301) on the drying of binary mixtures. For binary mixtures, we found that a gradient of colloidal osmotic pressure develops during drying and that this leads to the final film being stratified with a layer of smaller particles on top of a layer of larger particles. Here, we find that stratification by size is very general and also occurs in ternary and polydisperse mixtures. We name the segregation effect colloidal diffusiophoresis. In particular, we show that by changing the composition of a ternary mixture, different stratification morphologies can be achieved and hence the film properties can be tuned. In polydisperse spheres, colloidal diffusiophoresis leads to enrichment in the large particles at the bottom part of the film, whereas the top part is enriched with smaller particles. This segregation means that in the final film, the particle size distribution depends on height. Thus, the properties of the film will then depend on height. We propose a model that predicts a power-law dependence of the phoretic velocity on particle size. Results from the model and simulation show a good agreement.

Glass transition in hard-core fluids and beyond, using an effective static structure in the mode coupling theory

The dynamical arrest in classical fluids is studied using a simple modification of the mode coupling theory (MCT) aimed at correcting its overestimation of the tendency to glass formation while preserving its overall structure. As in previous attempts, the modification is based on the idea of tempering the static pair correlations used as input. It is implemented in this work by computing the static structure at a different state point than the one used to solve the MCT equation for the intermediate scattering function, using the pure hard-sphere glass for calibration. The location of the glass transition predicted from this modification is found to agree with simulations data for a variety of systems --pure fluids and mixtures with either purely repulsive interaction potentials or ones with attractive contributions. Besides improving the predictions in the long-time limit, and so reducing the non-ergodicity domain, the same modification works as well for the time-dependent correlators.

Viscosity of electrolyte solutions: a mode-coupling theory

We present a versatile theoretical method for calculating the steady-state viscosity and shear relaxation function of strong electrolyte solutions. In this method, the ions are described on a primitive model level as charged Brownian spheres, and the essential ion-ion hydrodynamic interactions (HIs) are accounted for in the shear relaxation effect of the ionic atmosphere. The method combines a many-component mode-coupling theory (MCT) approach by Nägele et al (1998 J. Chem. Phys. 108 9893) with a simplified solution scheme, leading to an analytic expression for the shear relaxation contribution to the viscosity. This expression accounts for both the excluded volumes of the ions and their HIs. We show that the limiting law results for the viscosity of electrolyte mixtures by Falkenhagen and by Onsager and Fuoss are recovered at very low concentrations, and we discuss HIs corrections appearing at higher concentrations. Our numerical results for a 1:1 electrolyte reveal a strong enlargement of the viscosity caused by the HIs. The high-frequency viscosity gives the largest contribution to the total viscosity at higher concentrations.

Non-ergodicity transition and multiple glasses in binary mixtures: on the accuracy of the input static structure in the mode coupling theory

We examine the question of the accuracy of the static correlation functions used as input in the mode coupling theory (MCT) of non-ergodic states in binary mixtures. We first consider hard-sphere mixtures and compute the static pair structure from the Ornstein-Zernike equations with the Percus-Yevick closure and more accurate ones that use bridge functions deduced from Rosenfeld's fundamental measures functional. The corresponding MCT predictions for the non-ergodicity lines and the transitions between multiple glassy states are determined from the long-time limit of the density autocorrelation functions. We find that while the non-ergodicity transition line is not very sensitive to the input static structure, up to diameter ratios D(2)/D(1) = 10, quantitative differences exist for the transitions between different glasses. The discrepancies with the more accurate closures become even qualitative for sufficiently asymmetric mixtures. They are correlated with the incorrect behavior of the PY structure at high size asymmetry. From the example of ultra-soft potential it is argued that this issue is of general relevance beyond the hard-sphere model.

Mode-coupling theory for the glass transition: test of the convolution approximation for short-range interactions

We reexamine the convolution approximation commonly used in the mode-coupling theory (MCT) of nonergodic states of classical fluids. This approximation concerns the static correlation functions used as input in the MCT treatment of the dynamics. Besides the hard-sphere model, we consider interaction potentials that present a short-range tail, either attractive or repulsive, beyond the hard core. By using accurate static correlation functions obtained from the fundamental measures functional for hard spheres, we show that the role of three-body direct correlations can be more significant than what is inferred from previous simple ansatzs for pure hard spheres. This may in particular impact the location of the glass transition line and the nonergodicity parameter.

Testing "microscopic" theories of glass-forming liquids

We assess the validity of "microscopic" approaches of glass-forming liquids based on the sole knowledge of the static pair density correlations. To do so, we apply them to a benchmark provided by two liquid models that share very similar static pair density correlation functions while displaying distinct temperature evolutions of their relaxation times. We find that the approaches are unsuccessful in describing the difference in the dynamical behavior of the two models. Our study is not exhaustive, and we have not tested the effect of adding corrections by including, for instance, three-body density correlations. Yet, our results appear strong enough to challenge the claim that the slowdown of relaxation in glass-forming liquids, for which it is well established that the changes of the static structure factor with temperature are small, can be explained by "microscopic" approaches only requiring the static pair density correlations as nontrivial input.

Self-assembly of nanoparticles into heterogeneous structures with gradient material properties

We present a mechanism to form self-assembled functional gradient superlattice structures by subjecting binary nanoparticles in an electric field. The interaction among different dipoles leads to the controllable formation of diverse structures, including particle columns with gradient material properties from inside to outside and various hierarchical layered or three-dimensional particle chain networks. We elucidate how permittivity, volume fraction, particle size, and the frequency of the electric field can be utilized to control the morphology of the induced structures, which would enable designed nanofabrication.

Statistical mechanics of homogeneous partly pinned fluid systems

The homogeneous partly pinned fluid systems are simple models of a fluid confined in a disordered porous matrix obtained by arresting randomly chosen particles in a one-component bulk fluid or one of the two components of a binary mixture. In this paper, their configurational properties are investigated. It is shown that a peculiar complementarity exists between the mobile and immobile phases, which originates from the fact that the solid is prepared in presence of and in equilibrium with the adsorbed fluid. Simple identities follow, which connect different types of configurational averages, either relative to the fluid-matrix system or to the bulk fluid from which it is prepared. Crucial simplifications result for the computation of important structural quantities, both in computer simulations and in theoretical approaches. Finally, possible applications of the model in the field of dynamics in confinement or in strongly asymmetric mixtures are suggested.

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.

Dynamical density functional theory for colloidal dispersions including hydrodynamic interactions

A dynamical density functional theory (DDFT) for translational Brownian dynamics is derived which includes hydrodynamic interactions. The theory reduces to the simple Brownian DDFT proposed by Marconi and Tarazona (U. Marini Bettolo Marconi and P. Tarazona, J. Chem. Phys. 110, 8032 (1999); J. Phys.: Condens. Matter 12, A413 (2000)) when hydrodynamic interactions are neglected. The derivation is based on Smoluchowski's equation for the time evolution of the probability density with pairwise hydrodynamic interactions. The theory is applied to hard-sphere colloids in an oscillating spherical optical trap which switches periodically in time from a stable confining to an unstable potential. Rosenfeld's fundamental measure theory for the equilibrium density functional is used and hydrodynamics are incorporated on the Rotne-Prager level. The results for the time-dependent density profiles are compared to extensive Brownian dynamics simulations which are performed on the same Rotne-Prager level and excellent agreement is obtained. It is further found that hydrodynamic interactions damp and slow the dynamics of the confined colloid cluster in comparison to the same situation with neglected hydrodynamic interactions.

Dynamical density functional theory with hydrodynamic interactions and colloids in unstable traps

A density functional theory for colloidal dynamics is presented which includes hydrodynamic interactions between the colloidal particles. The theory is applied to the dynamics of colloidal particles in an optical trap which switches periodically in time from a stable to an unstable confining potential. In the absence of hydrodynamic interactions, the resulting density breathing mode exhibits huge single peaked oscillations in the trap center which become double peaked and damped by hydrodynamic interactions. The predicted dynamical density fields are in good agreement with Brownian dynamics computer simulations.

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.

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.

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.

Long-time self-diffusion of charged colloidal particles: electrokinetic and hydrodynamic interaction effects

The authors analyze the long-time self-diffusion of charge-stabilized colloidal macroions in nondilute suspensions using a mode-coupling scheme developed for multicomponent suspensions of interacting Brownian spheres. In this scheme, all ionic species, including counterions and electrolyte ions, are treated on an equal footing as charged hard spheres undergoing overdamped Brownian motion. Hydrodynamic interactions between all ions are accounted for on the far-field level. We show that the influence on the colloidal long-time self-diffusion coefficient arising from the relaxation of the microionic atmosphere surrounding the colloids, the so-called electrolyte friction effect, is usually insignificant in comparison with the friction contributions arising from direct and hydrodynamic interactions between the colloidal particles. This finding is true even for small colloid concentrations unless the mobility difference between colloidal particles and microions is not large. Furthermore, we observe an interesting nonmonotonic density dependence of the colloidal long-time self-diffusion coefficient in suspensions with low amount of added salt. We show that this unusual density dependence is due to colloid-colloid hydrodynamic interactions.

Irreducible memory function and slow dynamics in disordered systems

We show how the irreducible memory function can be obtained in a rather straightforward way, and that it can be expressed in terms of two contributions representing two parallel decay channels. This representation should be useful for treating systems with a slow time dependence and where eventually some internal degrees of freedom enters in the relaxation process, and cuts off an underlying ideal ergodic to nonergodic transition. We also show how the irreducible memory function under certain mild conditions defines a regenerative stochastic process, or a two level stochastic system. This leads to a picture with dynamical heterogeneities, where the statistical properties asymptotically are ruled by limit processes. This can explain the universal behavior observed in many glass-forming systems.

Self-consistent generalized Langevin equation for colloidal mixtures

A self-consistent theory of collective and tracer diffusion in colloidal mixtures is presented. This theory is based on exact results for the partial intermediate scattering functions derived within the framework of the generalized Langevin equation formalism, plus a number of conceptually simple and sensible approximations. The first of these consists of a Vineyard-like approximation between collective and tracer diffusion, which writes the collective dynamics in terms of the memory function related to tracer diffusion. The second consists of interpolating this only unknown memory function between its two exact limits at small and large wave vectors; for this, a phenomenologically determined, but not arbitrary, interpolating function is introduced: a Lorentzian with its inflection point located at the first minimum of the partial static structure factor. The small wave-vector exact limit involves a time-dependent friction function, for which we take a general approximate result, previously derived within the generalized Langevin equation formalism. This general result expresses the time-dependent friction function in terms of the partial intermediate scattering functions, thus closing the system of equations into a fully self-consistent scheme. This extends to mixtures a recently proposed self-consistent theory developed for monodisperse suspensions [Yeomans-Reyna and Medina-Noyola, Phys. Rev. E 64, 066114 (2001)]. As an illustration of its quantitative accuracy, its application to a simple model of a binary dispersion in the absence of hydrodynamic interactions is reported.

Relaxation in a glassy binary mixture: comparison of the mode-coupling theory to a Brownian dynamics simulation

We solved the mode-coupling equations for the Kob-Andersen binary mixture using structure factors calculated from Brownian dynamics simulations of the same system. We found, as was previously observed, that the mode-coupling temperature T(c) inferred from simulations is about two times greater than that predicted by the theory. However, we find that many time-dependent quantities agree reasonably well with the predictions of the mode-coupling theory if they are compared at the same reduced temperature epsilon = (T - T(c))/T(c), and if epsilon is not too small. Specifically, the simulation results for the incoherent intermediate scattering function, the mean square displacement, the relaxation time, and the self-diffusion coefficient agree reasonably well with the predictions of the mode-coupling theory. We find that there are substantial differences for the non-Gaussian parameter. At small reduced temperatures the probabilities of the logarithm of single particle displacements demonstrate that there is hopping-like motion present in the simulations, and this motion is not predicted by the mode-coupling theory. The wave-vector-dependent relaxation time is shown to be qualitatively different from the predictions of the mode-coupling theory for temperatures where hopping-like motion is present.

Self-diffusion in a system of interacting Langevin particles

The behavior of the self-diffusion constant of Langevin particles interacting via a pairwise interaction is considered. The diffusion constant is calculated approximately within a perturbation theory in the potential strength about the bare diffusion constant. It is shown how this expansion leads to a systematic double expansion in the inverse temperature beta and the particle density rho. The one-loop diagrams in this expansion can be summed exactly and we show that this result is exact in the limit of small beta and rhobeta constants. The one-loop result can also be resummed using a semiphenomenological renormalization group method which has proved useful in the study of diffusion in random media. In certain cases the renormalization group calculation predicts the existence of a diverging relaxation time signaled by the vanishing of the diffusion constant, possible forms of divergence coming from this approximation are discussed. Finally, at a more quantitative level, the results are compared with numerical simulations, in two dimensions, of particles interacting via a soft potential recently used to model the interaction between coiled polymers.

Self-consistent theory of collective Brownian dynamics: theory versus simulation

A recently developed theory of collective diffusion in colloidal suspensions is tested regarding the quantitative accuracy of its description of the dynamics of monodisperse model colloidal systems without hydrodynamic interactions. The idea is to exhibit the isolated effects of the direct interactions, which constitute the main microscopic relaxation mechanism, in the absence of other effects, such as hydrodynamic interactions. Here we compare the numerical solution of the fully self-consistent theory with the results of Brownian dynamics simulation of the van Hove function G(r,t) and/or the intermediate scattering function F(k,t) of four simple model systems. Two of them are representative of short-ranged soft-core repulsive interactions [(sigma/r)(mu), with mu>>1], in two and in three dimensions. The other two involve long-ranged repulsive forces in two (dipolar, r(-3) potential) and in three (screened Coulomb, or repulsive Yukawa interactions) dimensions. We find that the theory, without any sort of adjustable parameters or rescaling prescriptions, provides an excellent approximate description of the collective dynamics of these model systems, particularly in the short- and intermediate-time regimes. We also compare our results with those of the single exponential approximation and with the competing mode-mode coupling theory.

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.

Overdamped van hove function of colloidal suspensions

The generalized-hydrodynamic theory for collective diffusion of a monodisperse colloidal suspension is developed in the framework of the Onsager-Machlup theory of time-dependent fluctuations. The time evolution of the intermediate scattering function F(k,t) is derived as a contraction of the description involving the instantaneous particle number concentration, the particle current, and the stress tensor of the Brownian fluid as state variables. We show that the proper overdamped limit of this equation requires the explicit separation of the stress tensor in its mutually orthogonal kinetic and configurational contributions. Analogous results also follow for the self-intermediate scattering function F(s)(k,t). We show that neglecting the non-Markovian part of the configurational stress tensor memory, one recovers the single exponential memory approximation (based on sum rules derived from the Smoluchowski equation) for both F(s)(k,t) and F(k,t). We suggest simple approximate manners to relate the collective and the self-memory functions, leading to Vineyard-like approximate relations between F(s)(k,t) and F(k,t).