Hydrodynamic correlation functions of a driven granular fluid in steady state

We study a homogeneously driven granular fluid of hard spheres at intermediate volume fractions and focus on time-delayed correlation functions in the stationary state. Inelastic collisions are modeled by incomplete normal restitution, allowing for efficient simulations with an event-driven algorithm. The incoherent scattering function Fincoh(q,t ) is seen to follow time-density superposition with a relaxation time that increases significantly as the volume fraction increases. The statistics of particle displacements is approximately Gaussian. For the coherent scattering function S(q,ω), we compare our results to the predictions of generalized fluctuating hydrodynamics, which takes into account that temperature fluctuations decay either diffusively or with a finite relaxation rate, depending on wave number and inelasticity. For sufficiently small wave number q we observe sound waves in the coherent scattering function S(q,ω) and the longitudinal current correlation function Cl(q,ω). We determine the speed of sound and the transport coefficients and compare them to the results of kinetic theory.

Temperature-Dependent Defect Dynamics in the Network Glass SiO2

We investigate the long time dynamics of a strong glass former, SiO2, below the glass transition temperature by averaging single-particle trajectories over time windows which comprise roughly 100 particle oscillations. The structure on this coarse-grained time scale is very well defined in terms of coordination numbers, allowing us to identify ill-coordinated atoms, which are called defects in the following. The most numerous defects are O-O neighbors, whose lifetimes are comparable to the equilibration time at low temperature. On the other hand, SiO and OSi defects are very rare and short lived. The lifetime of defects is found to be strongly temperature dependent, consistent with activated processes. Single-particle jumps give rise to local structural rearrangements. We show that in SiO2 these structural rearrangements are coupled to the creation or annihilation of defects, giving rise to very strong correlations of jumping atoms and defects.

Anomalous Velocity Distributions in Active Brownian Suspensions

Large-scale simulations and analytical theory have been combined to obtain the nonequilibrium velocity distribution, f(v), of randomly accelerated particles in suspension. The simulations are based on an event-driven algorithm, generalized to include friction. They reveal strongly anomalous but largely universal distributions, which are independent of volume fraction and collision processes, which suggests a one-particle model should capture all the essential features. We have formulated this one-particle model and solved it analytically in the limit of strong damping, where we find that f (v) decays as 1/v for multiple decades, eventually crossing over to a Gaussian decay for the largest velocities. Many particle simulations and numerical solution of the one-particle model agree for all values of the damping.

Strong Dynamical Heterogeneity and Universal Scaling in Driven Granular Fluids

Large-scale simulations of two-dimensional bidisperse granular fluids allow us to determine spatial correlations of slow particles via the four-point structure factor S-4 (q, t). Both cases, elastic (epsilon = 1) and inelastic (epsilon < 1) collisions, are studied. As the fluid approaches structural arrest, i.e., for packing fractions in the range 0.6 <= phi <= 0.805, scaling is shown to hold: S-4 (q, t)/chi(4)(t) = s(q xi(t)). Both the dynamic susceptibility chi(4)(tau(alpha)) and the dynamic correlation length xi(tau(alpha)) evaluated at the alpha relaxation time tau(alpha) can be fitted to a power law divergence at a critical packing fraction. The measured xi(tau(alpha)) widely exceeds the largest one previously observed for three-dimensional (3d) hard sphere fluids. The number of particles in a slow cluster and the correlation length are related by a robust power law, chi(4)(tau(alpha)) approximate to xi(d-p) (tau(alpha)), with an exponent d - p approximate to 1.6. This scaling is remarkably independent of epsilon, even though the strength of the dynamical heterogeneity at constant volume fraction depends strongly on epsilon.