Elasticity of frictionless particles near jamming

Scaling regimes and fluctuations of observables in computer glasses approaching the unjamming transition

Under decompression, disordered solids undergo an unjamming transition where they become under-coordinated and lose their structural rigidity. The mechanical and vibrational properties of these materials have been an object of theoretical, numerical, and experimental research for decades. In the study of low-coordination solids, understanding the behavior and physical interpretation of observables that diverge near the transition is of particular importance. Several such quantities are length scales (ξ or l) that characterize the size of excitations, the decay of spatial correlations, the response to perturbations, or the effect of physical constraints in the boundary or bulk of the material. Additionally, the spatial and sample-to-sample fluctuations of macroscopic observables such as contact statistics or elastic moduli diverge approaching unjamming. Here, we discuss important connections between all of these quantities and present numerical results that characterize the scaling properties of sample-to-sample contact and shear modulus fluctuations in ensembles of low-coordination disordered sphere packings and spring networks. Overall, we highlight three distinct scaling regimes and two crossovers in the disorder quantifiers χz and χμ as functions of system size N and proximity to unjamming δz. As we discuss, χX relates to the standard deviation σX of the sample-to-sample distribution of the quantity X (e.g., excess coordination δz or shear modulus μ) for an ensemble of systems. Importantly, χμ has been linked to experimentally accessible quantities that pertain to sound attenuation and the density of vibrational states in glasses. We investigate similarities and differences in the behaviors of χz and χμ near the transition and discuss the implications of our findings on current literature, unifying findings in previous studies.

Microrheology near jamming

The jamming transition is a nonequilibrium critical phenomenon, which governs characteristic mechanical properties of jammed soft materials, such as pastes, emulsions, and granular matters. Both experiments and theory of jammed soft materials have revealed that the complex modulus measured by conventional macrorheology exhibits a characteristic frequency dependence. Microrheology is a new type of method to obtain the complex modulus, which transforms the microscopic motion of probes to the complex modulus through the generalized Stokes relation (GSR). Although microrheology has been applied to jammed soft materials, its theoretical understanding is limited. In particular, the validity of the GSR near the jamming transition is far from obvious since there is a diverging length scale lc, which characterizes the heterogeneous response of jammed particles. Here, we study the microrheology of jammed particles by theory and numerical simulation. First, we develop a linear response formalism to calculate the response function of the probe particle, which is transformed to the complex modulus via the GSR. Then, we apply our formalism to a numerical model of jammed particles and find that the storage and loss modulus follow characteristic scaling laws near the jamming transition. Importantly, the observed scaling law coincides with that in macrorheology, which indicates that the GSR holds even near the jamming transition. We rationalize this equivalence by asymptotic analysis of the obtained formalism and numerical analysis on the displacement field of jammed particles under a local perturbation.

Finite amplitude waves in jammed matter

Here we use simulations and theory to show that, close to the jamming point, an arbitrary initial distortion of a granular media induces the formation of forward and backward non-linear finite amplitude waves. There are two regimes in the evolution of these waves (near field and far field). Initially, non-linear interactions between forward and backward waves dominate the propagation, leading to complex early evolution (near field). At longer times, forward and backwards waves cease interacting in the far field, and the propagation enters a new regime. Here the waves acquire a triangular-like profile, and evolve in a self-similar fashion characterized by a power law attenuation, whose exponent is weakly dependent on the initial pressure of the system. The finite amplitude waves gradually become linear waves when the amplitude of the initial distortion decreases, or the confining pressure on the system increases.

Universal scaling for disordered viscoelastic matter near the onset of rigidity

The onset of rigidity in interacting liquids, as they undergo a transition to a disordered solid, is associated with a rearrangement of the low-frequency vibrational spectrum. In this Letter, we derive scaling forms for the singular dynamical response of disordered viscoelastic networks near both jamming and rigidity percolation. Using effective-medium theory, we extract critical exponents, invariant scaling combinations, and analytical formulas for universal scaling functions near these transitions. Our scaling forms describe the behavior in space and time near the various onsets of rigidity, for rigid and floppy phases and the crossover region, including diverging length scales and timescales at the transitions.

Structural, mechanical, and vibrational properties of particulate physical gels

Our lives are surrounded by a rich assortment of disordered materials. In particular, glasses are well known as dense, amorphous materials, whereas gels exist in low-density, disordered states. Recent progress has provided a significant step forward in understanding the material properties of glasses, such as mechanical, vibrational, and transport properties. In contrast, our understanding of particulate physical gels is still highly limited. Here, using molecular dynamics simulations, we study a simple model of particulate physical gels, the Lennard-Jones (LJ) gels, and provide a comprehensive understanding of their structural, mechanical, and vibrational properties, all of which are markedly different from those of LJ glasses. First, the LJ gels show sparse, heterogeneous structures, and the length scale ξ_s of the structures grows as the density is lowered. Second, the LJ gels are extremely soft, with both shear G and bulk K moduli being orders of magnitude smaller than those of LJ glasses. Third, many low-frequency vibrational modes are excited, which form a characteristic plateau with the onset frequency ω_* in the vibrational density of states. Structural, mechanical, and vibrational properties, characterized by ξ_s, G, K, and ω_*, respectively, show power-law scaling behaviors with the density, which establishes a close relationship between them. Throughout this work, we also reveal that LJ gels are multiscale, solid-state materials: (i) homogeneous elastic bodies at long lengths, (ii) heterogeneous elastic bodies with fractal structures at intermediate lengths, and (iii) amorphous structural bodies at short lengths.

The role of friction in statistics and scaling laws of avalanches

We investigate statistics and scaling laws of avalanches in two-dimensional frictional particles by numerical simulations. We find that the critical exponent for avalanche size distributions is governed by microscopic friction between the particles in contact, where the exponent is larger and closer to mean-field predictions if the friction coefficient is finite. We reveal that microscopic "slips" between frictional particles induce numerous small avalanches which increase the slope, as well as the power-law exponent, of avalanche size distributions. We also analyze statistics and scaling laws of the avalanche duration and maximum stress drop rates, and examine power spectra of stress drop rates. Our numerical results suggest that the microscopic friction is a key ingredient of mean-field descriptions and plays a crucial role in avalanches observed in real materials.

Universal Relaxation Dynamics of Sphere Packings below Jamming

We show that non-Brownian suspensions of repulsive spheres below jamming display a slow relaxational dynamics with a characteristic timescale that diverges at jamming. This slow timescale is fully encoded in the structure of the unjammed packing and can be readily measured via the vibrational density of states. We show that the corresponding dynamic critical exponent is the same for randomly generated and sheared packings. Our results show that a wide variety of physical situations, from suspension rheology to algorithmic studies of the jamming transition are controlled by a unique diverging timescale, with a universal critical exponent.

Boson peak, elasticity, and glass transition temperature in polymer glasses: Effects of the rigidity of chain bending

The excess low-frequency vibrational spectrum, called boson peak, and non-affine elastic response are the most important particularities of glasses. Herein, the vibrational and mechanical properties of polymeric glasses are examined by using coarse-grained molecular dynamics simulations, with particular attention to the effects of the bending rigidity of the polymer chains. As the rigidity increases, the system undergoes a glass transition at a higher temperature (under a constant pressure), which decreases the density of the glass phase. The elastic moduli, which are controlled by the decrease of the density and the increase of the rigidity, show a non-monotonic dependence on the rigidity of the polymer chain that arises from the non-affine component. Moreover, a clear boson peak is observed in the vibrational density of states, which depends on the macroscopic shear modulus G. In particular, the boson peak frequency ω_BP is proportional to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{G}$$\end{document}G. These results provide a positive correlation between the boson peak, shear elasticity, and the glass transition temperature.

Mechanical response of dense pedestrian crowds to the crossing of intruders

The increasing number of mass events involving large crowds calls for a better understanding of the dynamics of dense crowds. Inquiring into the possibility of a mechanical description of these dynamics, we experimentally study the crossing of dense static crowds by a cylindrical intruder, a mechanical test which is classical for granular matter. The analysis of our experiments reveals robust features in the crowds' response, comprising both similarities and discrepancies with the response of granular media. Common features include the presence of a depleted region behind the intruder and the short-range character of the perturbation. On the other hand, unlike grains, pedestrians anticipate the intruder's passage by moving much before contact and their displacements are mostly lateral, hence not aligned with the forces exerted by the intruder. Similar conclusions are reached when the intruder is not a cylinder, but a single crossing pedestrian. Thus, our work shows that pedestrian interactions even at high densities (3 to 6 ped/m2) do not reduce to mechanical ones. More generally, the avoidance strategies evidenced by our findings question the incautious use of force models for dense crowds.

Two Diverging Length Scales in the Structure of Jammed Packings

At densities higher than the jamming transition for athermal, frictionless repulsive spheres we find two distinct length scales, both of which diverge as a power law as the transition is approached. The first, ξ_{Z}, is associated with the two-point correlation function for the number of contacts on two particles as a function of the particle separation. The second, ξ_{f}, is associated with contact-number fluctuations in subsystems of different sizes. On scales below ξ_{f}, the fluctuations are highly suppressed, similar to the phenomenon of hyperuniformity usually associated with density fluctuations. The exponents for the divergence of ξ_{Z} and ξ_{f} are different and appear to be different in two and three dimensions.

Quasilocalized states of self stress in packing-derived networks

States of self stress (SSS) are assignments of forces on the edges of a network that satisfy mechanical equilibrium in the absence of external forces. In this work we show that a particular class of quasilocalized SSS in packing-derived networks, first introduced by D.M. Sussman, C.P. Goodrich, A.J. Liu (Soft Matter 12, 3982 (2016)), are characterized by a decay length that diverges as [Formula: see text] , where [Formula: see text] is the mean connectivity of the network, and [Formula: see text] is the Maxwell threshold in two dimensions, at odds with previous claims. Our results verify the previously proposed analogy between quasilocalized SSS and the mechanical response to a local dipolar force in random networks of relaxed Hookean springs. We show that the normalization factor that distinguishes between quasilocalized SSS and the response to a local dipole constitutes a measure of the mechanical coupling of the forced spring to the elastic network in which it is embedded. We further demonstrate that the lengthscale that characterizes quasilocalized SSS does not depend on its associated degree of mechanical coupling, but instead only on the network connectivity.

Correlation and shear bands in a plastically deformed granular medium

Recent experiments (Le Bouil et al., Phys. Rev. Lett., 2014, 112, 246001) have analyzed the statistics of local deformation in a granular solid undergoing plastic deformation. Experiments report strongly anisotropic correlation between events, with a characteristic angle that was interpreted using elasticity theory and the concept of Eshelby transformations with dilation; interestingly, the shear bands that characterize macroscopic failure occur at an angle that is different from the one observed in microscopic correlations. Here, we interpret this behavior using a mesoscale elastoplastic model of solid flow that incorporates a local Mohr-Coulomb failure criterion. This differs from the interpretation of Le Bouil et al., which is based on purely elastic considerations ignoring the potential role of local friction on deformation patterns. We show that the angle observed in the microscopic correlations can be understood by combining the elastic interactions associated with Eshelby transformation with the local failure criterion. At large strains, we also induce permanent shear bands at an angle that is different from the one observed in the correlation pattern. We interpret this angle as the one that leads to the maximal instability of slip lines.

Viscous forces and bulk viscoelasticity near jamming

When weakly jammed packings of soft, viscous, non-Brownian spheres are probed mechanically, they respond with a complex admixture of elastic and viscous effects. While many of these effects are understood for specific, approximate models of the particles' interactions, there are a number of proposed force laws in the literature, especially for viscous interactions. We numerically measure the complex shear modulus G* of jammed packings for various viscous force laws that damp relative velocities between pairs of contacting particles or between a particle and the continuous fluid phase. We find a surprising sensitive dependence of G* on the viscous force law: the system may or may not display dynamic critical scaling, and the exponents describing how G* scales with frequency can change. We show that this sensitivity is closely linked to manner in which viscous damping couples to floppy-like, non-affine motion, which is prominent near jamming.

Nonlocal Elasticity near Jamming in Frictionless Soft Spheres

We use simulations of frictionless soft sphere packings to identify novel constitutive relations for linear elasticity near the jamming transition. By forcing packings at varying wavelengths, we directly access their transverse and longitudinal compliances. These are found to be wavelength dependent, in violation of conventional (local) linear elasticity. Crossovers in the compliances select characteristic length scales, which signify the appearance of nonlocal effects. Two of these length scales diverge as the pressure vanishes, indicating that critical effects near jamming control the breakdown of local elasticity. We expect these nonlocal constitutive relations to be applicable to a wide range of weakly jammed solids, including emulsions, foams, and granulates.

Elastic moduli and vibrational modes in jammed particulate packings

When we elastically impose a homogeneous, affine deformation on amorphous solids, they also undergo an inhomogeneous, nonaffine deformation, which can have a crucial impact on the overall elastic response. To correctly understand the elastic modulus M, it is therefore necessary to take into account not only the affine modulus M_{A}, but also the nonaffine modulus M_{N} that arises from the nonaffine deformation. In the present work, we study the bulk (M=K) and shear (M=G) moduli in static jammed particulate packings over a range of packing fractions φ. The affine M_{A} is determined essentially by the static structural arrangement of particles, whereas the nonaffine M_{N} is related to the vibrational eigenmodes. We elucidate the contribution of each vibrational mode to the nonaffine M_{N} through a modal decomposition of the displacement and force fields. In the vicinity of the (un)jamming transition φ_{c}, the vibrational density of states g(ω) shows a plateau in the intermediate-frequency regime above a characteristic frequency ω^{*}. We illustrate that this unusual feature apparent in g(ω) is reflected in the behavior of M_{N}: As φ→φ_{c}, where ω^{*}→0, those modes for ω<ω^{*} contribute less and less, while contributions from those for ω>ω^{*} approach a constant value which results in M_{N} to approach a critical value M_{Nc}, as M_{N}-M_{Nc}∼ω^{*}. At φ_{c} itself, the bulk modulus attains a finite value K_{c}=K_{Ac}-K_{Nc}>0, such that K_{Nc} has a value that remains below K_{Ac}. In contrast, for the critical shear modulus G_{c}, G_{Nc} and G_{Ac} approach the same value so that the total value becomes exactly zero, G_{c}=G_{Ac}-G_{Nc}=0. We explore what features of the configurational and vibrational properties cause such a distinction between K and G, allowing us to validate analytical expressions for their critical values.

Spatial structure of states of self stress in jammed systems

States of self stress, organizations of internal forces in many-body systems that are in equilibrium with an absence of external forces, can be thought of as the constitutive building blocks of the elastic response of a material. In overconstrained disordered packings they have a natural mathematical correspondence with the zero-energy vibrational modes in underconstrained systems. While substantial attention in the literature has been paid to diverging length scales associated with zero- and finite-energy vibrational modes in jammed systems, less is known about the spatial structure of the states of self stress. In this work we define a natural way in which a unique state of self stress can be associated with each bond in a disordered spring network derived from a jammed packing, and then investigate the spatial structure of these bond-localized states of self stress. This allows for an understanding of how the elastic properties of a system would change upon changing the strength or even existence of any bond in the system.

Spatial Distributions of Local Elastic Moduli Near the Jamming Transition

Recent progress on studies of the nanoscale mechanical responses in disordered systems has highlighted a strong degree of heterogeneity in the elastic moduli. In this contribution, using computer simulations, we study the elastic heterogeneities in athermal amorphous solids--composed of isotropic static sphere packings--near the jamming transition. We employ techniques based on linear response methods that are amenable to experimentation. We find that the local elastic moduli are randomly distributed in space and are described by Gaussian probability distributions, thereby lacking any significant spatial correlations, that persist all the way down to the transition point. However, the shear modulus fluctuations grow as the jamming threshold is approached, which is characterized by a new power-law scaling. Through this diverging behavior we are able to identify a characteristic length scale, associated with shear modulus heterogeneities, that distinguishes between bulk and local elastic responses.