Scaling ansatz for the jamming transition

Jamming is a first-order transition with quenched disorder in amorphous materials sheared by cyclic quasistatic deformations

Jamming is an athermal transition between flowing and rigid states in amorphous systems such as granular matter, colloidal suspensions, complex fluids and cells. The jamming transition seems to display mixed aspects of a first-order transition, evidenced by a discontinuity in the coordination number, and a second-order transition, indicated by power-law scalings and diverging lengths. Here we demonstrate that jamming is a first-order transition with quenched disorder in cyclically sheared systems with quasistatic deformations, in two and three dimensions. Based on scaling analyses, we show that fluctuations of the jamming density in finite-sized systems have important consequences on the finite-size effects of various quantities, resulting in a square relationship between disconnected and connected susceptibilities, a key signature of the first-order transition with quenched disorder. This study puts the jamming transition into the category of a broad class of transitions in disordered systems where sample-to-sample fluctuations dominate over thermal fluctuations, suggesting that the nature and behavior of the jamming transition might be better understood within the developed theoretical framework of the athermally driven random-field Ising model.

Field Theory for Mechanical Criticality in Disordered Fiber Networks

Strain-controlled criticality governs the elasticity of jamming and fiber networks. While the upper critical dimension of jamming is believed to be d_{u}=2, non-mean-field exponents are observed in numerical studies of 2D and 3D fiber networks. The origins of this remains unclear. In this study we propose a minimal mean-field model for strain-controlled criticality of fiber networks. We then extend this to a phenomenological field theory, in which non-mean-field behavior emerges as a result of the disorder in the network structure. We predict that the upper critical dimension for such systems is d_{u}=4 using a Gaussian approximation. Moreover, we identify an order parameter for the phase transition, which has been lacking for fiber networks to date.

Unified Understanding of Nonlinear Rheology near the Jamming Transition Point

When slowly sheared, jammed packings respond elastically before yielding. This linear elastic regime becomes progressively narrower as the jamming transition point is approached, and rich nonlinear rheologies such as shear softening and hardening or melting emerge. However, the physical mechanism of these nonlinear rheologies remains elusive. To clarify this, we numerically study jammed packings of athermal frictionless soft particles under quasistatic shear γ. We find the universal scaling behavior for the ratio of the shear stress σ and the pressure P, independent of the preparation protocol of the initial configurations. In particular, we reveal shear softening σ/P∼γ^{1/2} over an unprecedentedly wide range of strain up to the yielding point, which a simple scaling argument can rationalize.

Long-range correlations in elastic moduli and local stresses at the unjamming transition

We explore the behaviour of spatially heterogeneous elastic moduli as well as the correlations between local moduli in model solids with short-range repulsive potentials. We show through numerical simulations that local elastic moduli exhibit long-range correlations, similar to correlations in the local stresses. Specifically, the correlations in local shear moduli exhibit anisotropic behavior at large lengthscales characterized by pinch-point singularities in Fourier space, displaying a structural pattern akin to shear stress correlations. Focussing on two-dimensional jammed solids approaching the unjamming transition, we show that stress correlations exhibit universal properties, characterized by a quadratic p2 dependence of the correlations as the pressure p approaches zero, independent of the details of the model. In contrast, the modulus correlations exhibit a power-law dependence with different exponents depending on the specific interaction potential. Furthermore, we illustrate that while affine responses lack long-range correlations, the total modulus, which encompasses non-affine behavior, exhibits long-range correlations.

Critical-like slowdown in thermal soft-sphere glasses via energy minimization

Using hybrid molecular dynamics/SWAP Monte Carlo (MD/SMC) simulations, we show that while the terminal relaxation times τ(ϕ) for FIRE energy minimization of soft-sphere glasses can decrease by orders of magnitude as sample equilibration proceeds and the jamming density ϕ_{J} increases, they always scale as τ(ϕ)∼(ϕ_{J}-ϕ)^{-2}∼[Z_{iso}-Z_{ms}(τ)]^{-2}, where Z_{iso}=2d and Z_{ms}(τ) is the average coordination number of particles satisfying a minimal local mechanical stability criterion (Z≥d+1) at the top of the final potential-energy-landscape (PEL) sub-basin the system encounters. This scaling allows us to collapse τ datasets that look very different when plotted as a function of ϕ, and to address a closely related question: how does the character of the PEL basins that dense thermal glasses most typically occupy evolve as the glasses age at constant ϕ and T?

Jamming on convex deformable surfaces

Jamming is a fundamental transition that governs the behavior of particulate media, including sand, foams and dense suspensions. Upon compression, such media change from freely flowing to a disordered, marginally stable solid that exhibits non-Hookean elasticity. While the jamming process is well established for fixed geometries, the nature and dynamics of jamming for a diverse class of soft materials and deformable substrates, including emulsions and biological matter, remains unknown. Here we propose a new scenario, metric jamming, where rigidification occurs on a surface that has been deformed from its ground state. Unlike classical jamming processes that exhibit discrete mechanical transitions, surprisingly we find that metric jammed states possess mechanical properties continuously tunable between those of classically jammed and conventional elastic media. The compact and curved geometry significantly alters the vibrational spectra of the structures relative to jamming in flat Euclidean space, and metric jammed systems also possess new types of vibrational mode that couple particle and shape degrees of freedom. Our work provides a theoretical framework that unifies our understanding of solidification processes that take place on deformable media and lays the groundwork to exploit jamming for the control and stabilization of shape in self-assembly processes.

Designing the pressure-dependent shear modulus using tessellated granular metamaterials

Jammed packings of granular materials display complex mechanical response. For example, the ensemble-averaged shear modulus 〈G〉 increases as a power law in pressure p for static packings of soft spherical particles that can rearrange during compression. We seek to design granular materials with shear moduli that can either increase or decrease with pressure without particle rearrangements even in the large-system limit. To do this, we construct tessellated granular metamaterials by joining multiple particle-filled cells together. We focus on cells that contain a small number of bidisperse disks in two dimensions. We first study the mechanical properties of individual disk-filled cells with three types of boundaries: periodic boundary conditions (PBC), fixed-length walls (FXW), and flexible walls (FLW). Hypostatic jammed packings are found for cells with FLW, but not in cells with PBC and FXW, and they are stabilized by quartic modes of the dynamical matrix. The shear modulus of a single cell depends linearly on p. We find that the slope of the shear modulus with pressure λ_{c}<0 for all packings in single cells with PBC where the number of particles per cell N≥6. In contrast, single cells with FXW and FLW can possess λ_{c}>0, as well as λ_{c}<0, for N≤16. We show that we can force the mechanical properties of multicell granular metamaterials to possess those of single cells by constraining the end points of the outer walls and enforcing an affine shear response. These studies demonstrate that tessellated granular metamaterials provide a platform for the design of soft materials with specified mechanical properties.

Shear hardening in frictionless amorphous solids near the jamming transition

The jamming transition, generally manifested by a rapid increase of rigidity under compression (i.e. compression hardening), is ubiquitous in amorphous materials. Here we study shear hardening in deeply annealed frictionless packings generated by numerical simulations, reporting critical scalings absent in compression hardening. We demonstrate that hardening is a natural consequence of shear-induced memory destruction. Based on an elasticity theory, we reveal two independent microscopic origins of shear hardening: (i) the increase of the interaction bond number and (ii) the emergence of anisotropy and long-range correlations in the orientations of bonds-the latter highlights the essential difference between compression and shear hardening. Through the establishment of physical laws specific to anisotropy, our work completes the criticality and universality of jamming transition, and the elasticity theory of amorphous solids.

Hard-sphere jamming through the lens of linear optimization

The jamming transition is ubiquitous. It is present in granular matter, foams, colloids, structural glasses, and many other systems. Yet, it defines a critical point whose properties still need to be fully understood. Recently, a major breakthrough came about when the replica formalism was extended to build a mean-field theory that provides an exact description of the jamming transition of spherical particles in the infinite-dimensional limit. While such theory explains the jamming critical behavior of both soft and hard spheres, investigating the transition in finite-dimensional systems poses very difficult and different problems, in particular from the numerical point of view. Soft particles are modeled by continuous potentials; thus, their jamming point can be reached through efficient energy minimization algorithms. In contrast, the latter methods are inapplicable to hard-sphere (HS) systems since the interaction energy among the particles is always zero by construction. To overcome these difficulties, here we recast the jamming of hard spheres as a constrained optimization problem and introduce the CALiPPSO algorithm, capable of readily producing jammed HS packings without including any effective potential. This algorithm brings a HS configuration of arbitrary dimensions to its jamming point by solving a chain of linear optimization problems. We show that there is a strict correspondence between the force balance conditions of jammed packings and the properties of the optimal solutions of CALiPPSO, whence we prove analytically that our packings are always isostatic and in mechanical equilibrium. Furthermore, using extensive numerical simulations, we show that our algorithm is able to probe the complex structure of the free-energy landscape, finding qualitative agreement with mean-field predictions. We also characterize the algorithmic complexity of CALiPPSO and provide an open-source implementation of it.

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.

Jammed solids with pins: Thresholds, force networks, and elasticity

The role of fixed degrees of freedom in soft or granular matter systems has broad applicability and theoretical interest. Here we address questions of the geometrical role that a scaffolding of fixed particles plays in tuning the threshold volume fraction and force network in the vicinity of jamming. Our two-dimensional simulated system consists of soft particles and fixed "pins," both of which harmonically repel overlaps. On the one hand, we find that many of the critical scalings associated with jamming in the absence of pins continue to hold in the presence of even dense pin latices. On the other hand, the presence of pins lowers the jamming threshold in a universal way at low pin densities and a geometry-dependent manner at high pin densities, producing packings with lower densities and fewer contacts between particles. The onset of strong lattice dependence coincides with the development of bond-orientational order. Furthermore, the presence of pins dramatically modifies the network of forces, with both unusually weak and unusually strong forces becoming more abundant. The spatial organization of this force network depends on pin geometry and is described in detail. Using persistent homology, we demonstrate that pins modify the topology of the network. Finally, we observe clear signatures of this developing bond-orientational order and broad force distribution in the elastic moduli which characterize the linear response of these packings to strain.

Finite-size scaling of critical avalanches

We examine probability distribution for avalanche sizes observed in self-organized critical systems. While a power-law distribution with a cutoff because of finite system size is typical behavior, a systematic investigation reveals that it may also decrease with increasing the system size at a fixed avalanche size. We implement the scaling method and identify scaling functions. The data collapse ensures a correct estimation of the critical exponents and distinguishes two exponents related to avalanche size and system size. Our simple analysis provides striking implications. While the exact value for avalanches size exponent remains elusive for the prototype sandpile on a square lattice, we suggest the exponent should be 1. The simulation results represent that the distribution shows a logarithmic system size dependence, consistent with the normalization condition. We also argue that for the train or Oslo sandpile model with bulk drive, the avalanche size exponent is slightly less than 1, which differs significantly from the previous estimate of 1.11.

Criticality and marginal stability of the shear jamming transition of frictionless soft spheres

We study numerically the critical behavior and marginal stability of the shear jamming transition for frictionless soft spheres, observed to occur over a finite range of densities, associated with isotropic jamming for densities above the minimum jamming (J-point) density. Several quantities are shown to scale near the shear jamming point in the same way as the isotropic jamming point. We compute the exponents associated with the small force distribution and the interparticle gap distribution and show that the corresponding exponents are consistent with the marginal stability condition observed for isotropic jamming and with predictions of the mean-field theory of jamming in hard spheres.

Flow and arrest in stressed granular materials

Flowing granular materials often abruptly arrest if not driven by sufficient applied stresses. Such abrupt cessation of motion can be economically expensive in industrial materials handling and processing, and is significantly consequential in intermittent geophysical phenomena such as landslides and earthquakes. Using discrete element simulations, we calculate states of steady flow and arrest for granular materials under the conditions of constant applied pressure and shear stress, which are also most relevant in practice. Here the material can dilate or compact, and flow or arrest, in response to the applied stress. Our simulations highlight that under external stress, the intrinsic response of granular materials is characterized by uniquely-defined steady states of flow or arrest, which are highly sensitive to interparticle friction. While the flowing states can be equivalently characterized by volume fraction, coordination number or internal stress ratio, to characterize the states of shear arrest, one needs to also consider the structural anisotropy in the contact network. We highlight the role of dilation in the flow-arrest transition, and discuss our findings in the context of rheological transitions in granular materials.

Finite-size effects in the microscopic critical properties of jammed configurations: A comprehensive study of the effects of different types of disorder

Jamming criticality defines a universality class that includes systems as diverse as glasses, colloids, foams, amorphous solids, constraint satisfaction problems, neural networks, etc. A particularly interesting feature of this class is that small interparticle forces (f) and gaps (h) are distributed according to nontrivial power laws. A recently developed mean-field (MF) theory predicts the characteristic exponents of these distributions in the limit of very high spatial dimension, d→∞ and, remarkably, their values seemingly agree with numerical estimates in physically relevant dimensions, d=2 and 3. These exponents are further connected through a pair of inequalities derived from stability conditions, and both theoretical predictions and previous numerical investigations suggest that these inequalities are saturated. Systems at the jamming point are thus only marginally stable. Despite the key physical role played by these exponents, their systematic evaluation has yet to be attempted. Here, we carefully test their value by analyzing the finite-size scaling of the distributions of f and h for various particle-based models for jamming. Both dimension and the direction of approach to the jamming point are also considered. We show that, in all models, finite-size effects are much more pronounced in the distribution of h than in that of f. We thus conclude that gaps are correlated over considerably longer scales than forces. Additionally, remarkable agreement with MF predictions is obtained in all but one model, namely near-crystalline packings. Our results thus help to better delineate the domain of the jamming universality class. We furthermore uncover a secondary linear regime in the distribution tails of both f and h. This surprisingly robust feature is understood to follow from the (near) isostaticity of our configurations.

Shear response of granular packings compressed above jamming onset

We investigate the mechanical response of jammed packings of repulsive, frictionless spherical particles undergoing isotropic compression. Prior simulations of the soft-particle model, where the repulsive interactions scale as a power law in the interparticle overlap with exponent α, have found that the ensemble-averaged shear modulus 〈G(P)〉 increases with pressure P as ∼P^{(α-3/2)/(α-1)} at large pressures. 〈G〉 has two key contributions: (1) continuous variations as a function of pressure along geometrical families, for which the interparticle contact network does not change, and (2) discontinuous jumps during compression that arise from changes in the contact network. Using numerical simulations, we show that the form of the shear modulus G^{f} for jammed packings within near-isostatic geometrical families is largely determined by the affine response G^{f}∼G_{a}^{f}, where G_{a}^{f}/G_{a0}=(P/P_{0})^{(α-2)/(α-1)}-P/P_{0}, P_{0}∼N^{-2(α-1)} is the characteristic pressure at which G_{a}^{f}=0, G_{a0} is a constant that sets the scale of the shear modulus, and N is the number of particles. For near-isostatic geometrical families that persist to large pressures, deviations from this form are caused by significant nonaffine particle motion. We further show that the ensemble-averaged shear modulus 〈G(P)〉 is not simply a sum of two power laws, but 〈G(P)〉∼(P/P_{c})^{a}, where a≈(α-2)/(α-1) in the P→0 limit and 〈G(P)〉∼(P/P_{c})^{b}, where b≳(α-3/2)/(α-1), above a characteristic pressure that scales as P_{c}∼N^{-2(α-1)}.

Mean-Field Predictions of Scaling Prefactors Match Low-Dimensional Jammed Packings

No known analytic framework precisely explains all the phenomena observed in jamming. The replica theory for glasses and jamming is a mean-field theory which attempts to do so by working in the limit of infinite dimensions, such that correlations between neighbors are negligible. As such, results from this mean-field theory are not guaranteed to be observed in finite dimensions. However, many results in mean field for jamming have been shown to be exact or nearly exact in low dimensions. This suggests that the infinite dimensional limit is not necessary to obtain these results. In this Letter, we perform precision measurements of jamming scaling relationships between pressure, excess packing fraction, and number of excess contacts from dimensions 2-10 in order to extract the prefactors to these scalings. While these prefactors should be highly sensitive to finite dimensional corrections, we find the mean-field predictions for these prefactors to be exact in low dimensions. Thus the mean-field approximation is not necessary for deriving these prefactors. We present an exact, first-principles derivation for one, leaving the other as an open question. Our results suggest that mean-field theories of critical phenomena may compute more for d≥d_{u} than has been previously appreciated.

Rolling away from the Wall into Granular Matter

The sedimentation of solid objects into granular matter near boundaries is an almost virgin field of research. Here we describe in detail the penetration dynamics of a cylindrical object into a quasi-2D granular medium. By tracking the trajectory of the cylinder as it penetrates the granular bed, we characterize two distinct kinds of motion: its center of mass moves horizontally away from the lateral wall, and it rotates around its symmetry axis. While the repulsion is caused by the loading of force chains between the intruder and the wall, the rotation can be associated to the frictional forces between the grains and the intruder. Finally, we show the analogies between the sedimentation of twin intruders released far from any boundaries, and that of one intruder released near a vertical wall.

Finite size effects in critical fiber networks

Fibrous networks such as collagen are common in physiological systems. One important function of these networks is to provide mechanical stability for cells and tissues. At physiological levels of connectivity, such networks would be mechanically unstable with only central-force interactions. While networks can be stabilized by bending interactions, it has also been shown that they exhibit a critical transition from floppy to rigid as a function of applied strain. Beyond a certain strain threshold, it is predicted that underconstrained networks with only central-force interactions exhibit a discontinuity in the shear modulus. We study the finite-size scaling behavior of this transition and identify both the mechanical discontinuity and critical exponents in the thermodynamic limit. We find both non-mean-field behavior and evidence for a hyperscaling relation for the critical exponents, for which the network stiffness is analogous to the heat capacity for thermal phase transitions. Further evidence for this is also found in the self-averaging properties of fiber networks.

Direct measurement of force configurational entropy in jamming

Thermal fluctuations are not large enough to lead to state changes in granular materials. However, in many cases, such materials do achieve reproducible bulk properties, suggesting that they are controlled by an underlying statistical mechanics analogous to thermodynamics. While themodynamic descriptions of granular materials have been explored, they have not yet been concretely connected to their underlying statistical mechanics. We make this connection concrete by providing a first-principles derivation of the multiplicity and thus the entropy of the force networks in granular packings. We directly measure the multiplicity of force networks using a protocol based on the phase space volume of allowed force configurations. Analogous to Planck's constant, we find a scale factor, h_{f}, that discretizes this phase space volume into a multiplicity. To determine this scale factor, we measure angoricity over a wide range of pressures using the method of overlapping histograms and find that in the absence of a fundamental quantum scale it is set solely by the system size and dimensionality. This concretely links thermodynamic approaches of angoricity with the microscopic multiplicity of the force network ensemble.

Structured randomness: jamming of soft discs and pins

Simulations are used to find the zero temperature jamming threshold, φj, for soft, bidisperse disks in the presence of small fixed particles, or "pins", arranged in a lattice. The presence of pins leads, as one expects, to a decrease in φj. Structural properties of the system near the jamming threshold are calculated as a function of the pin density. While the correlation length exponent remains ν = 1/2 at low pin densities, the system is mechanically stable with more bonds, yet fewer contacts than the Maxwell criterion implies in the absence of pins. In addition, as pin density increases, novel bond orientational order and long-range spatial order appear, which are correlated with the square symmetry of the pin lattice.

Lindemann Unjamming of Emulsions

We study the bulk and shear elastic properties of barely-compressed, "athermal" emulsions and find that the rigidity of the jammed solid fails at remarkably large critical osmotic pressures. The minuscule yield strain and similarly small Brownian particle displacement of solid emulsions close to this transition suggests that this catastrophic failure corresponds to a plastic-entropic instability: the solid becomes too soft and weak to resist the thermal agitation of the droplets that compose it and fails. We propose a modified Lindemann stability criterion to describe this transition and derive a scaling law for the critical osmotic pressure that agrees quantitatively with experimental observations.

Pressure Dependent Shear Response of Jammed Packings of Frictionless Spherical Particles

The mechanical response of packings of purely repulsive, spherical particles to athermal, quasistatic simple shear near jamming onset is highly nonlinear. Previous studies have shown that, at small pressure p, the ensemble-averaged static shear modulus ⟨G-G_{0}⟩ scales with p^{α}, where α≈1, but above a characteristic pressure p^{**}, ⟨G-G_{0}⟩∼p^{β}, where β≈0.5. However, we find that the shear modulus G^{i} for an individual packing typically decreases linearly with p along a geometrical family where the contact network does not change. We resolve this discrepancy by showing that, while the shear modulus does decrease linearly within geometrical families, ⟨G⟩ also depends on a contribution from discontinuous jumps in ⟨G⟩ that occur at the transitions between geometrical families. For p>p^{**}, geometrical-family and rearrangement contributions to ⟨G⟩ are of opposite signs and remain comparable for all system sizes. ⟨G⟩ can be described by a scaling function that smoothly transitions between two power-law exponents α and β. We also demonstrate the phenomenon of compression unjamming, where a jammed packing unjams via isotropic compression.

Tuning and jamming reduced to their minima

Inspired by protein folding, we smooth out the complex cost function landscapes of two processes: the tuning of networks and the jamming of ideal spheres. In both processes, geometrical frustration plays a role-tuning pressure differences between pairs of target nodes far from the source in a flow network impedes tuning of nearby pairs more than the reverse process, while unjamming the system in one region can make it more difficult to unjam elsewhere. By modifying the cost functions to control the order in which functions are tuned or regions unjam, we smooth out local minima while leaving global minima unaffected, increasing the success rate for reaching global minima.

Characterizing nonaffinity upon decompression of soft-sphere packings

Athermal elastic moduli of soft-sphere packings are known to exhibit universal scaling properties near the unjamming point, most notably the vanishing of the shear-to-bulk moduli ratio G/B upon decompression. Interestingly, the smallness of G/B stems from the large nonaffinity of deformation-induced displacements under shear strains, compared to insignificant nonaffinity of displacements under compressive strains. In this work, we show using numerical simulations that the relative weights of the affine and nonaffine contributions to the bulk modulus, and their dependence on the proximity to the unjamming point, can differ qualitatively between different models that feature the same generic unjamming phenomenology. In canonical models of unjamming, we observe that the ratio of the nonaffine to total bulk moduli B_{na}/B approaches a constant upon decompression, while in other, less well-studied models, it vanishes. We show that the vanishing of B_{na}/B in noncanonical models stems from the emergence of an invariance of net (zero) forces on the constituent particles to compressive strains at the onset of unjamming. We provide a theoretical scaling analysis that fully explains our numerical observations, and allows us to predict the scaling behavior of B_{na}/B upon unjamming, given the functional form of the pairwise interaction potential.

Scaling Theory for Mechanical Critical Behavior in Fiber Networks

As a function of connectivity, spring networks exhibit a critical transition between floppy and rigid phases at an isostatic threshold. For connectivity below this threshold, fiber networks were recently shown theoretically to exhibit a rigidity transition with corresponding critical signatures as a function of strain. Experimental collagen networks were also shown to be consistent with these predictions. We develop a scaling theory for this strain-controlled transition. Using a real-space renormalization approach, we determine relations between the critical exponents governing the transition, which we verify for the strain-controlled transition using numerical simulations of both triangular lattice-based and packing-derived fiber networks.

The physics of jamming for granular materials: a review

Granular materials consist of macroscopic grains, interacting via contact forces, and unaffected by thermal fluctuations. They are one of a class systems that undergo jamming, i.e. a transition between fluid-like and disordered solid-like states. Roughly twenty years ago, proposals by Cates et al for the shear response of colloidal systems and by Liu and Nagel, for a universal jamming diagram in a parameter space of packing fraction, ϕ, shear stress, τ, and temperature, T raised key questions. Contemporaneously, experiments by Howell et al and numerical simulations by Radjai et al and by Luding et al helped provide a starting point to explore key insights into jamming for dry, cohesionless, granular materials. A recent experimental observation by Bi et al is that frictional granular materials have a a re-entrant region in their jamming diagram. In a range of ϕ, applying shear strain, γ, from an initially force/stress free state leads to fragile (in the sense of Cates et al), then anisotropic shear jammed states. Shear jamming at fixed ϕ is presumably conjugate to Reynolds dilatancy, involving dilation under shear against deformable boundaries. Numerical studies by Radjai and Roux showed that Reynolds dilatancy does not occur for frictionless systems. Recent numerical studies by several groups show that shear jamming occurs for finite, but not infinite, systems of frictionless grains. Shear jamming does not lead to known ordering in position space, but Sarkar et al showed that ordering occurs in a space of force tiles. Experimental studies seeking to understand random loose and random close packings (rlp and rcp) and dating back to Bernal have probed granular packings and their response to shear and intruder motion. These studies suggest that rlp's are anisotropic and shear-jammed-like, whereas rcp's are likely isotropically jammed states. Jammed states are inherently static, but the jamming diagram may provide a context for understanding rheology, i.e. dynamic shear in a variety of systems that include granular materials and suspensions.

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.

No unjamming transition in a Voronoi model of biological tissue

Vertex models are a popular approach to simulating the mechanical and dynamical properties of dense biological tissues, describing the tissue as a network of polygons. Recently a class of two-dimensional vertex models was shown to exhibit a disordered rigidity transition controlled by the preferred cellular geometry, which was subsequently echoed by experimental findings. An attractive variant of these models uses a Voronoi tessellation to describe the cells, which reduces the number of degrees of freedom as compared the original vertex model. The Voronoi model was also endowed with a non-equilibrium model of cellular motility, leading to rich, glassy behavior. This glassy behavior was suggested to be inextricably linked to an underlying jamming transition. We test this conjecture, exploring the low-effective-temperature limit of the 2D Voronoi model by studying cell trajectories from detailed dynamical simulations in combination with rigidity measurements of energy-minimized disordered cell configurations. We find that the zero-temperature limit of this model has no unjamming transition. We show that this absence of an unjamming transition is intimately linked to the marginality of the model, i.e. the fact that the constraints imposed on cell areas and perimeters precisely balance the number of degrees of freedom in the model. Our work suggests that constraint counting arguments are useful to understand rigidity in a broad class of models of dense biological tissues.

Softening and yielding of soft glassy materials

Solids deform and fluids flow, but soft glassy materials, such as emulsions, foams, suspensions, and pastes, exhibit an intricate mix of solid- and liquid-like behavior. While much progress has been made to understand their elastic (small strain) and flow (infinite strain) properties, such understanding is lacking for the softening and yielding phenomena that connect these asymptotic regimes. Here we present a comprehensive framework for softening and yielding of soft glassy materials, based on extensive numerical simulations of oscillatory rheological tests, and show that two distinct scenarios unfold depending on the material's packing density. For dense systems, there is a single, pressure-independent strain where the elastic modulus drops and the particle motion becomes diffusive. In contrast, for weakly jammed systems, a two-step process arises: at an intermediate softening strain, the elastic and loss moduli both drop down and then reach a new plateau value, whereas the particle motion becomes diffusive at the distinctly larger yield strain. We show that softening is associated with an extensive number of microscopic contact changes leading to a non-analytic rheological signature. Moreover, the scaling of the softening strain with pressure suggest the existence of a novel pressure scale above which softening and yielding coincide, and we verify the existence of this crossover scale numerically. Our findings thus evidence the existence of two distinct classes of soft glassy materials - jamming dominated and dense - and show how these can be distinguished by their rheological fingerprint.

Universal behaviour of the glass and the jamming transitions in finite dimensions for hard spheres

We investigate the glass and the jamming transitions of hard spheres in finite dimensions d, through a revised cell theory, that combines the free volume and the Random First Order Theory (RFOT). Recent results show that in infinite dimension the ideal glass transition and jamming transitions are distinct, while based on our theory we argue that they indeed coincide for finite d. As a consequence, jamming results into a percolation transition described by RFOT, with a static length diverging with exponent ν = 2/d, which we verify through finite size scaling, and standard critical exponents α = 0, β = 0 and γ = 2 independent on d.

Unjamming in models with analytic pairwise potentials

Canonical models for studying the unjamming scenario in systems of soft repulsive particles assume pairwise potentials with a sharp cutoff in the interaction range. The sharp cutoff renders the potential nonanalytic but makes it possible to describe many properties of the solid in terms of the coordination number z, which has an unambiguous definition in these cases. Pairwise potentials without a sharp cutoff in the interaction range have not been studied in this context, but should in fact be considered to understand the relevance of the unjamming phenomenology in systems where such a cutoff is not present. In this work we explore two systems with such interactions: an inverse power law and an exponentially decaying pairwise potential, with the control parameters being the exponent (of the inverse power law) for the former and the number density for the latter. Both systems are shown to exhibit the characteristic features of the unjamming transition, among which are the vanishing of the shear-to-bulk modulus ratio and the emergence of an excess of low-frequency vibrational modes. We establish a relation between the pressure-to-bulk modulus ratio and the distance to unjamming in each of our model systems. This allows us to predict the dependence of other key observables on the distance to unjamming. Our results provide the means for a quantitative estimation of the proximity of generic glass-forming models to the unjamming transition in the absence of a clear-cut definition of the coordination number and highlight the general irrelevance of nonaffine contributions to the bulk modulus.

Discontinuous change of shear modulus for frictional jammed granular materials

The shear modulus of jammed frictional granular materials with harmonic repulsive interaction under an oscillatory shear is numerically investigated. It is confirmed that the storage modulus, the real part of the shear modulus, for frictional grains with sufficiently small strain amplitude γ_{0} discontinuously emerges at the jamming transition point. The storage modulus for small γ_{0} differs from that of frictionless grains even in the zero friction limit, whereas they are almost identical with each other for sufficiently large γ_{0}, where the transition becomes continuous. The stress-strain curve exhibits a hysteresis loop even for a small strain, which connects a linear region for sufficiently small strain to another linear region for larger strain. We propose a scaling law to interpolate between the states of small and large γ_{0}.

Scaling Theory for the Frictionless Unjamming Transition

We develop a scaling theory of the unjamming transition of soft frictionless disks in two dimensions by defining local areas, which can be uniquely assigned to each contact. These serve to define local order parameters, whose distribution exhibits divergences as the unjamming transition is approached. We derive scaling forms for these divergences from a mean-field approach that treats the local areas as noninteracting entities, and demonstrate that these results agree remarkably well with numerical simulations. We find that the asymptotic behavior of the scaling functions arises from the geometrical structure of the packing while the overall scaling with the compression energy depends on the force law. We use the scaling forms of the distributions to determine the scaling of the total grain area A_{G} and the total number of contacts N_{C}.