Gustavo Düring

A unified framework for non-Brownian suspension flows and soft amorphous solids

While the rheology of non-Brownian suspensions in the dilute regime is well understood, their behavior in the dense limit remains mystifying. As the packing fraction of particles increases, particle motion becomes more collective, leading to a growing length scale and scaling properties in the rheology as the material approaches the jamming transition. There is no accepted microscopic description of this phenomenon. However, in recent years it has been understood that the elasticity of simple amorphous solids is governed by a critical point, the unjamming transition where the pressure vanishes, and where elastic properties display scaling and a diverging length scale. The correspondence between these two transitions is at present unclear. Here we show that for a simple model of dense flow, which we argue captures the essential physics near the jamming threshold, a formal analogy can be made between the rheology of the flow and the elasticity of simple networks. This analogy leads to a new conceptual framework to relate microscopic structure to rheology. It enables us to define and compute numerically normal modes and a density of states. We find striking similarities between the density of states in flow, and that of amorphous solids near unjamming: both display a plateau above some frequency scale ω∗ ∼ |zc - z|, where z is the coordination of the network of particle in contact, zc = 2D where D is the spatial dimension. However, a spectacular difference appears: the density of states in flow displays a single mode at another frequency scale ωmin ≪ ω∗ governing the divergence of the viscosity.

Breakdown of continuum elasticity in amorphous solids

We show numerically that the response of simple amorphous solids (elastic networks and particle packings) to a local force dipole is characterized by a lengthscale lc that diverges as unjamming is approached as lc ∼ (z - 2d)(-1/2), where z ≥ 2d is the mean coordination, and d is the spatial dimension, at odds with previous numerical claims. We also show how the magnitude of the lengthscale lc is amplified by the presence of internal stresses in the disordered solid. Our data suggests a divergence of lc ∼ (pc - p)(-1/4) with proximity to a critical internal stress pc at which soft elastic modes become unstable.

Low-energy non-linear excitations in sphere packings

We study theoretically and numerically how hard frictionless particles in random packings can rearrange. We demonstrate the existence of two distinct unstable non-linear modes of rearrangement, both associated with the opening and the closing of contacts. The first mode, whose density is characterized by some exponent θ′, corresponds to motions of particles extending throughout the entire system. The second mode, whose density is characterized by an exponent θ ≠ θ′, corresponds to the local buckling of a few particles. Extended modes are shown to yield at a much higher rate than local modes when a stress is applied. We show that the distribution of contact forces follows P(f) ∼ fmin(θ′,θ), and that imposing the restriction that the packing cannot be densified further leads to the bounds and , where γ characterizes the singularity of the pair distribution function g(r) at contact. These results extend the theoretical analysis of [Wyart, Phys. Rev. Lett., 2012, 109, 125502] where the existence of local modes was not considered. We perform numerics that support that these bounds are saturated with γ ≈ 0.38, θ ≈ 0.17 and θ′ ≈ 0.44. We measure systematically the stability of all such modes in packings, and confirm their marginal stability. The principle of marginal stability thus allows us to make clearcut predictions on the ensemble of configurations visited in these out-of-equilibrium systems, and on the contact forces and pair distribution functions. It also reveals the excitations that need to be included in a description of plasticity or flow near jamming, and suggests a new path to study two-level systems and soft spots in simple amorphous solids of repulsive particles.

Why glass elasticity affects the thermodynamics and fragility of supercooled liquids

Supercooled liquids are characterized by their fragility: The slowing down of the dynamics under cooling is more sudden and the jump of specific heat at the glass transition is generally larger in fragile liquids than in strong ones. Despite the importance of this quantity in classifying liquids, explaining what aspects of the microscopic structure controls fragility remains a challenge. Surprisingly, experiments indicate that the linear elasticity of the glass—a purely local property of the free energy landscape—is a good predictor of fragility. In particular, materials presenting a large excess of soft elastic modes, the so-called boson peak, are strong. This is also the case for network liquids near the rigidity percolation, known to affect elasticity. Here we introduce a model of the glass transition based on the assumption that particles can organize locally into distinct configurations that are coupled spatially via elasticity. The model captures the mentioned observations connecting elasticity and fragility. We find that materials presenting an abundance of soft elastic modes have little elastic frustration: Energy is insensitive to most directions in phase space, leading to a small jump of specific heat. In this framework strong liquids turn out to lie the closest to a critical point associated with a rigidity or jamming transition, and their thermodynamic properties are related to the problem of number partitioning and to Hopfield nets in the limit of small memory.

Wave turbulence in vibrating plates: The effect of damping

The effect of damping in the wave turbulence regime for thin vibrating plates is studied. An experimental method, allowing measurements of dissipation in the system at all scales, is first introduced. Practical experimental devices for increasing the dissipation are used. The main observable consequence of increasing the damping is a significant modification in the slope of the power spectral density, so that the observed power laws are not in a pure inertial regime. However, the system still displays a turbulent behavior with a cut-off frequency that is determined by the injected power which does not depend on damping. By using the measured damping power-law in numerical simulations, similar conclusions are drawn out.

Statistics and Properties of Low-Frequency Vibrational Modes in Structural Glasses

Low-frequency vibrational modes play a central role in determining various basic properties of glasses, yet their statistical and mechanical properties are not fully understood. Using extensive numerical simulations of several model glasses in three dimensions, we show that in systems of linear size L sufficiently smaller than a crossover size L_{D}, the low-frequency tail of the density of states follows D(ω)∼ω^{4} up to the vicinity of the lowest Goldstone mode frequency. We find that the sample-to-sample statistics of the minimal vibrational frequency in systems of size L<L_{D} is Weibullian, with scaling exponents in excellent agreement with the ω^{4} law. We further show that the lowest-frequency modes are spatially quasilocalized and that their localization and associated quartic anharmonicity are largely frequency independent. The effect of preparation protocols on the low-frequency modes is elucidated, and a number of glassy length scales are briefly discussed.

Unified theory of inertial granular flows and non-Brownian suspensions

Rheological properties of dense flows of hard particles are singular as one approaches the jamming threshold where flow ceases both for aerial granular flows dominated by inertia and for over-damped suspensions. Concomitantly, the length scale characterizing velocity correlations appears to diverge at jamming. Here we introduce a theoretical framework that proposes a tentative, but potentially complete, scaling description of stationary flows. Our analysis, which focuses on frictionless particles, applies both to suspensions and inertial flows of hard particles. We compare our predictions with the empirical literature, as well as with novel numerical data. Overall, we find a very good agreement between theory and observations, except for frictional inertial flows whose scaling properties clearly differ from frictionless systems. For overdamped flows, more observations are needed to decide if friction is a relevant perturbation. Our analysis makes several new predictions on microscopic dynamical quantities that should be accessible experimentally.

Length scales and self-organization in dense suspension flows.

Dense non-Brownian suspension flows of hard particles display mystifying properties: As the jamming threshold is approached, the viscosity diverges, as well as a length scale that can be identified from velocity correlations. To unravel the microscopic mechanism governing dissipation and its connection to the observed correlation length, we develop an analogy between suspension flows and the rigidity transition occurring when floppy networks are pulled, a transition believed to be associated with the stress stiffening of certain gels. After deriving the critical properties near the rigidity transition, we show numerically that suspension flows lie close to it. We find that this proximity causes a decoupling between viscosity and the correlation length of velocities ξ, which scales as the length l(c) characterizing the response to a local perturbation, previously predicted to follow l(c)∼ 1/sqrt[z(c)-z] ∼ p(0.18), where p is the dimensionless particle pressure, z is the coordination of the contact network made by the particles, and z(c) is twice the spatial dimension. We confirm these predictions numerically and predict the existence of a larger length scale l(r)∼sqrt[p] with mild effects on velocity correlation and of a vanishing strain scale δγ ∼ 1/p that characterizes decorrelation in flow.

Simulations of driven overdamped frictionless hard spheres

Abstract We introduce an event-driven simulation scheme for overdamped dynamics of frictionless hard spheres subjected to external forces, neglecting hydrodynamic interactions. Our event-driven approach is based on an exact equation of motion which relates the driving force to the resulting velocities through the geometric information characterizing the underlying network of contacts between the hard spheres. Our method allows for a robust extraction of the instantaneous coordination of the particles as well as contact force statistics and dynamics, under any chosen driving force, in addition to shear flow and compression. It can also be used for generating high-precision jammed packings under shear, compression, or both. We present a number of additional applications of our method.

Symmetry Induced Four-Wave Capillary Wave Turbulence

We report theoretical and experimental results on 4-wave capillary wave turbulence. A system consisting of two immiscible and incompressible fluids of the same density can be written in a Hamiltonian way for the conjugated pair (eta, Psi). Adding the symmetry z --> -z, the set of capillary waves display a Kolmogorov-Zakharov spectrum k(-4) in wave vector space and f(-8/3) in the frequency domain. The wave system is studied experimentally with two immiscible fluids of almost equal densities (water and silicon oil) where the capillary surface waves are excited by a low-frequency random forcing. The probability density function of the local wave amplitude shows a quasi-Gaussian behavior and the power spectral density is shows a power-law behavior in frequency with a slope of -2.75. Theoretical and experimental results are in fairly good agreement with each other.