Nearly logarithmic decay in the colloidal hard-sphere system

Anomalous diffusion, non-Gaussianity, nonergodicity, and confinement in stochastic-scaled Brownian motion with diffusing diffusivity dynamics

Scaled Brownian motions (SBMs) with power-law time-dependent diffusivity have been used to describe various types of anomalous diffusion yet Gaussian observed in granular gases kinetics, turbulent diffusion, and molecules mobility in cells, to name a few. However, some of these systems may exhibit non-Gaussian behavior which can be described by SBM with diffusing diffusivity (DD-SBM). Here, we numerically investigate both free and confined DD-SBM models characterized by fixed or stochastic scaling exponent of time-dependent diffusivity. The effects of distributed scaling exponent, random diffusivity, and confinement are considered. Different regimes of ultraslow diffusion, subdiffusion, normal diffusion, and superdiffusion are observed. In addition, weak ergodic and non-Gaussian behaviors are also detected. These results provide insights into diffusion in time-fluctuating diffusivity landscapes with potential applications to time-dependent temperature systems spreading in heterogeneous environments.

Impulse response function for Brownian motion

Motivated from the central role of the mean-square displacement and its second time-derivative - that is the velocity autocorrelation function in the description of Brownian motion and its implications to microrheology, we revisit the physical meaning of the first time-derivative of the mean-square displacement of Brownian particles. By employing a rheological analogue for Brownian motion, we show that the time-derivative of the mean-square displacement of Brownian microspheres with mass m and radius R immersed in any linear, isotropic viscoelastic material is identical to , where h(t) is the impulse response function (strain history γ(t), due to an impulse stress τ(t) = δ(t - 0)) of a rheological network that is a parallel connection of the linear viscoelastic material with an inerter with distributed inertance . The impulse response function of the viscoelastic material-inerter parallel connection derived in this paper at the stress-strain level of the rheological analogue is essentially the response function of the Brownian particles expressed at the force-displacement level by Nishi et al. after making use of the fluctuation-dissipation theorem. By employing the viscoelastic material-inerter rheological analogue we derive the mean-square displacement and its time-derivatives of Brownian particles immersed in a viscoelastic material described with a Maxwell element connected in parallel with a dashpot and we show that for Brownian motion of microparticles immersed in such fluid-like materials, the impulse response function h(t) maintains a finite constant value in the long term.

Theoretical Estimate of the Glass Transition Line of Yukawa One-Component Plasmas

The mode coupling theory of supercooled liquids is combined with advanced closures to the integral equation theory of liquids in order to estimate the glass transition line of Yukawa one-component plasmas from the unscreened Coulomb limit up to the strong screening regime. The present predictions constitute a major improvement over the current literature predictions. The calculations confirm the validity of an existing analytical parameterization of the glass transition line. It is verified that the glass transition line is an approximate isomorphic curve and the value of the corresponding reduced excess entropy is estimated. Capitalizing on the isomorphic nature of the glass transition line, two structural vitrification indicators are identified that allow a rough estimate of the glass transition point only through simple curve metrics of the static properties of supercooled liquids. The vitrification indicators are demonstrated to be quasi-universal by an investigation of hard sphere and inverse power law supercooled liquids. The straightforward extension of the present results to bi-Yukawa systems is also discussed.

Persistent Anti-Correlations in Brownian Dynamics Simulations of Dense Colloidal Suspensions Revealed by Noise Suppression

Transport properties of a hard-sphere colloidal fluid are investigated by Brownian dynamics simulations. We implement a novel algorithm for the time-dependent velocity-autocorrelation function (VACF) essentially eliminating the noise of the bare random motion. The measured VACF reveals persistent anti-correlations manifested by a negative algebraic power-law tail t^{-5/2} at all densities. At small packing fractions the simulations fully agree with the analytic low-density prediction, yet the amplitude of the tail becomes dramatically suppressed as the packing fraction is increased. The mode-coupling theory of the glass transition provides a qualitative explanation for the strong variation in terms of the static compressibility as well as the slowing down of the structural relaxation.

Observation of Strongly Heterogeneous Dynamics at the Depinning Transition in a Colloidal Glass

We study experimentally the origin of heterogeneous dynamics in strongly driven glass-forming systems. Thereto, we apply a well-defined force with a laser line trap on individual colloidal polystyrene probe particles seeded in an emulsion glass composed of droplets of the same size. Fluid and glass states can be probed. We monitor the trajectories of the probe and measure displacements and their distributions. Our experiments reveal intermittent dynamics around a depinning transition at a threshold force. For smaller forces, linear response connects mean displacement, and quiescent mean squared displacement. Mode coupling theory calculations rationalize the observations.

Underdamped scaled Brownian motion: (non-)existence of the overdamped limit in anomalous diffusion

It is quite generally assumed that the overdamped Langevin equation provides a quantitative description of the dynamics of a classical Brownian particle in the long time limit. We establish and investigate a paradigm anomalous diffusion process governed by an underdamped Langevin equation with an explicit time dependence of the system temperature and thus the diffusion and damping coefficients. We show that for this underdamped scaled Brownian motion (UDSBM) the overdamped limit fails to describe the long time behaviour of the system and may practically even not exist at all for a certain range of the parameter values. Thus persistent inertial effects play a non-negligible role even at significantly long times. From this study a general questions on the applicability of the overdamped limit to describe the long time motion of an anomalously diffusing particle arises, with profound consequences for the relevance of overdamped anomalous diffusion models. We elucidate our results in view of analytical and simulations results for the anomalous diffusion of particles in free cooling granular gases.

Glass-transition properties of Yukawa potentials: from charged point particles to hard spheres

The glass transition is investigated in three dimensions for single and double Yukawa potentials for the full range of control parameters. For vanishing screening parameter, the limit of the one-component plasma is obtained; for large screening parameters and high coupling strengths, the glass-transition properties cross over to the hard-sphere system. Between the two limits, the entire transition diagram can be described by analytical functions. Unlike other potentials, the glass-transition and melting lines for Yukawa potentials are found to follow shifted but otherwise identical curves in control-parameter space.

Probability Densities of a Forced Probe Particle in Glass: Results from Mode Coupling Theory and Simulations of Active Microrheology

We investigate the displacements of a probe particle inside a glass, when a strong external force is applied to the probe (active nonlinear microrheology). Calculations within mode coupling theory are presented for glasses of hard spheres and compared to Langevin and Brownian dynamics simulations. Under not too strong forces where the probe remains trapped, the probe density distribution becomes anisotropic. It is shifted towards the direction of the force, develops an enhanced tail in that direction (signalled by a positive skewness), and exhibits different variances along and perpendicular to the force direction. A simple model of an harmonically trapped probe rationalizes the low force limit, with strong strain softening setting in at forces of the order of a few thermal energies per particle radius.

Structural relaxation of polydisperse hard spheres: comparison of the mode-coupling theory to a Langevin dynamics simulation

We analyze the slow glassy structural relaxation as measured through collective and tagged-particle density correlation functions obtained from Brownian dynamics simulations for a polydisperse system of quasi-hard spheres in the framework of the mode-coupling theory (MCT) of the glass transition. Asymptotic analyses show good agreement for the collective dynamics when polydispersity effects are taken into account in a multicomponent calculation, but qualitative disagreement at small q when the system is treated as effectively monodisperse. The origin of the different small-q behavior is attributed to the interplay between interdiffusion processes and structural relaxation. Numerical solutions of the MCT equations are obtained taking properly binned partial static structure factors from the simulations as input. Accounting for a shift in the critical density, the collective density correlation functions are well described by the theory at all densities investigated in the simulations, with quantitative agreement best around the maxima of the static structure factor and worst around its minima. A parameter-free comparison of the tagged-particle dynamics however reveals large quantitative errors for small wave numbers that are connected to the well-known decoupling of self-diffusion from structural relaxation and to dynamical heterogeneities. While deviations from MCT behavior are clearly seen in the tagged-particle quantities for densities close to and on the liquid side of the MCT glass transition, no such deviations are seen in the collective dynamics.

Glass transition of hard spheres in high dimensions

We have investigated analytically and numerically the liquid-glass transition of hard spheres for dimensions d-->infinity in the framework of mode-coupling theory. The numerical results for the critical collective and self-nonergodicity parameters fc(k;d) and fc(s)(k;d) exhibit non-Gaussian k dependence even up to d=800.fc(s)(k;d) and fc(k;d) differ for k approximately d1/2, but become identical on a scale k approximately d, which is proven analytically. The critical packing fraction phic(d) approximately d(2)2(-d) is above the corresponding Kauzmann packing fraction phiK(d) derived by a small cage expansion. Its quadratic pre-exponential factor is different from the linear one found earlier. The numerical values for the exponent parameter and therefore the critical exponents a and b depend on d, even for the largest values of d.

Comparison of dynamic light scattering measurements and mode-coupling theory for the tagged particle dynamics of a hard-sphere suspension

The mean-squared displacement, velocity autocorrelation function, and the non-Gaussian parameter, obtained by dynamic light scattering on suspensions of particles with hard-sphere interactions, are compared with the results of the idealized version of mode-coupling theory. Both leading order asymptotic and full numerical solutions of the mode-coupling equations are considered. Experiment and the full numerical results of the theory expose similar qualitative changes at the volume fraction of the first order freezing transition. In particular, the emergence of negative algebraic decays in the velocity autocorrelation function of the undercooled suspension suggest the emergence of clusters in which particles are trapped. Consistency of experiment, computer simulation, and theory in this regard suggests that, at particular strengths of the delayed, nonlinear feedback, contained in mode coupling theory, the latter predicts not only structural arrest which, as already established, is symptomatic of a glass transition, but also a more subtle change in dynamics that signals the onset of the first order transition.

Dynamical heterogeneity and the freezing transition in hard-sphere suspensions: further analysis of the mean square displacement and the velocity autocorrelation function

The velocity autocorrelation function is derived from the mean-squared displacement measured on a colloidal suspension of particles with hard-sphere-like interactions. It decays to zero from below and follows a stretched exponential function of delay time for the thermodynamically stable suspension. For the metastable suspension a power-law decay emerges. The results are discussed in terms of the classical Lorentz gas and the model that describes diffusion confined to one dimension. With the aid of these models, the experimental results provide a characterization of the dynamical heterogeneities which are observed microscopically, and an explanation for the enhanced resistance to flow and diffusion usually found in undercooled fluids upon approaching the glass transition.

Depolarized light scattering versus optical Kerr effect. II. Insight into the dynamic susceptibility of molecular liquids

We have previously discussed [J. Chem. Phys. 125, 114502 (2006)] that optical Kerr effect (OKE) and depolarized light scattering (DLS) data of molecular liquids reveal, each in their native domain, the same characteristic signatures of the glass transition dynamics; in particular, the intermediate power law of OKE is equivalent with the excess wing of the frequency-domain data, long since known in dielectric spectroscopy. We now extend the discussion to show that the excess wing is an equally common feature in DLS. We further discuss the time-temperature superposition property of OKE data in relation to our DLS and literature dielectric-spectroscopic results, and the merits of their mode coupling theory analyses. Spectroscopic signatures of a liquid-crystal-forming system (nematogen) are discussed in the same frame.

Cole-Cole law for critical dynamics in glass-forming liquids

Within the mode-coupling theory (MCT) for glassy dynamics, the asymptotic low-frequency expansions for the dynamical susceptibilities at critical points are compared to the expansions for the dynamic moduli; this shows that the convergence properties of the two expansions can be quite different. In some parameter regions, the leading-order expansion formula for the modulus describes the solutions of the MCT equations of motion outside the transient regime successfully; at the same time, the leading- and next-to-leading-order expansion formulas for the susceptibility fail. In these cases, one can derive a Cole-Cole law for the susceptibilities; and this law accounts for the dynamics for frequencies below the band of microscopic excitations and above the high-frequency part of the alpha peak. It is shown that this scenario explains the optical-Kerr-effect data measured for salol and benzophenone (BZP). For BZP it is inferred that the depolarized light-scattering spectra exhibit a wing for the alpha peak within the Gigahertz band. This wing results from the crossover of the von Schweidler law part of the alpha peak to the high-frequency part of the Cole-Cole peak; and this crossover can be described quantitatively by the leading-order formulas of MCT for the modulus.

Dynamical perspective of the freezing transition of a suspension of hard spheres from the velocity autocorrelation function

The velocity autocorrelation function of concentrated colloidal fluids of hard-sphere particles, measured by dynamic light scattering, decays to the experimental noise floor from below. The decay follows a stretched exponential function of delay time for the colloidal fluid in thermodynamic equilibrium, and a power law for the nonequilibrium, undercooled colloidal fluid. Consideration of this difference between the two states points to a possible dynamical mechanism of freezing.

Random-walk analysis of displacement statistics of particles in concentrated suspensions of hard spheres

Mean-squared displacements (MSDs) of colloidal fluids of hard spheres are analyzed in terms of a random walk, an analysis which assumes that the process of structural relaxation among the particles can be described in terms of thermally driven memoryless encounters. For the colloidal fluid in thermodynamic equilibrium the magnitude of the stretching of the MSD is able to be reconciled by a bias in the walk. This description fails for the under-cooled colloidal fluid.