Jamming in finite systems: stability, anisotropy, fluctuations, and scaling

Scaling ansatz for the jamming transition

We propose a Widom-like scaling ansatz for the critical jamming transition. Our ansatz for the elastic energy shows that the scaling of the energy, compressive strain, shear strain, system size, pressure, shear stress, bulk modulus, and shear modulus are all related to each other via scaling relations, with only three independent scaling exponents. We extract the values of these exponents from already known numerical or theoretical results, and we numerically verify the resulting predictions of the scaling theory for the energy and residual shear stress. We also derive a scaling relation between pressure and residual shear stress that yields insight into why the shear and bulk moduli scale differently. Our theory shows that the jamming transition exhibits an emergent scale invariance, setting the stage for the potential development of a renormalization group theory for jamming.

Topological boundary modes in jammed matter

Granular matter at the jamming transition is poised on the brink of mechanical stability, and hence it is possible that these random systems have topologically protected surface phonons. Studying two model systems for jammed matter, we find states that exhibit distinct mechanical topological classes, protected surface modes, and ubiquitous Weyl points. The detailed statistics of the boundary modes shed surprising light on the properties of the jamming critical point and help inform a common theoretical description of the detailed features of the transition.

Beyond linear elasticity: jammed solids at finite shear strain and rate

The shear response of soft solids can be modeled with linear elasticity, provided the forcing is slow and weak. Both of these approximations must break down when the material loses rigidity, such as in foams and emulsions at their (un)jamming point - suggesting that the window of linear elastic response near jamming is exceedingly narrow. Yet precisely when and how this breakdown occurs remains unclear. To answer these questions, we perform computer simulations of stress relaxation and shear start-up tests in athermal soft sphere packings, the canonical model for jamming. By systematically varying the strain amplitude, strain rate, distance to jamming, and system size, we identify characteristic strain and time scales that quantify how and when the window of linear elasticity closes, and relate these scales to changes in the microscopic contact network.

Critical scaling of Bagnold rheology at the jamming transition of frictionless two-dimensional disks

We carry out constant volume simulations of steady-state shear-driven rheology in a simple model of bidisperse soft-core frictionless disks in two dimensions, using a dissipation law that gives rise to Bagnoldian rheology. We discuss in detail the critical scaling ansatz for the shear-driven jamming transition and carry out a detailed scaling analysis of our resulting data for pressure p and shear stress σ. Our analysis determines the critical exponent β that describes the algebraic divergence of the Bagnold transport coefficients lim_{γ[over ̇]→0}p/γ[over ̇]^{2},σ/γ[over ̇]^{2}∼(ϕ_{J}-ϕ)^{-β} as the jamming transition ϕ_{J} is approached from below. For the low strain rates considered in this work, we show that it is still necessary to consider the leading correction-to-scaling term in order to achieve a self-consistent analysis of our data, in which the critical parameters become independent of the size of the window of data used in the analysis. We compare our resulting value β≈5.0±0.4 against previous numerical results and competing theoretical models. Our results confirm that the shear-driven jamming transition in Bagnoldian systems is well described by a critical scaling theory and we relate this scaling theory to the phenomenological constituent laws for dilatancy and friction.

Fluctuations in flows near jamming

Bubbles, droplets or particles in flowing complex media such as foams, emulsions or suspensions follow highly complex paths, with the relative motion of the constituents setting the energy dissipation rate. What is their dynamics, and how is this connected to the global rheology? To address these questions, we probe the statistics and spatio-temporal organization of the local particle motion and energy dissipation in a model for sheared disordered materials. We find that the fluctuations in the local dissipation vary from nearly Gaussian and homogeneous at low densities and fast flows, to strongly intermittent for large densities and slow flows. The higher order moments of the relative particle velocities reveal strong evidence for a qualitative difference between two distinct regimes which are nevertheless connected by a smooth crossover. In the critical regime, the higher order moments are related by novel multiscaling relations. In the plastic regime the relations between these moments take on a different form, with higher moments diverging rapidly when the flow rate vanishes. As these velocity differences govern the energy dissipation, we can distinguish two qualitatively different types of flow: an intermediate density, critical regime related to jamming, and a large density, plastic regime.

Residual Defect Density in Random Disks Deposits

We investigate the residual distribution of structural defects in very tall packings of disks deposited randomly in large channels. By performing simulations involving the sedimentation of up to 50 × 10(9) particles we find all deposits to consistently show a non-zero residual density of defects obeying a characteristic power-law as a function of the channel width. This remarkable finding corrects the widespread belief that the density of defects should vanish algebraically with growing height. A non-zero residual density of defects implies a type of long-range spatial order in the packing, as opposed to only local ordering. In addition, we find deposits of particles to involve considerably less randomness than generally presumed.

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.

Rigidity loss in disordered systems: three scenarios

We reveal significant qualitative differences in the rigidity transition of three types of disordered network materials: randomly diluted spring networks, jammed sphere packings, and stress-relieved networks that are diluted using a protocol that avoids the appearance of floppy regions. The marginal state of jammed and stress-relieved networks are globally isostatic, while marginal randomly diluted networks show both overconstrained and underconstrained regions. When a single bond is added to or removed from these isostatic systems, jammed networks become globally overconstrained or floppy, whereas the effect on stress-relieved networks is more local and limited. These differences are also reflected in the linear elastic properties and point to the highly effective and unusual role of global self-organization in jammed sphere packings.

Shear modulus and dilatancy softening in granular packings above jamming

We investigate experimentally the mechanical response to shear of a monolayer of bidisperse frictional grains across the jamming transition. We inflate an intruder inside the packing and use photoelasticity and tracking techniques to measure the induced shear strain and stresses at the grain scale. We quantify experimentally the constitutive relations for strain amplitudes as low as 10(-3) and for a range of packing fractions within 2% variation around the jamming transition. At the transition strong nonlinear effects set in: both the shear modulus and the dilatancy shear soften at small strain until a critical strain is reached where effective linearity is recovered. The scaling of the critical strain and the associated critical stresses on the distance to jamming are extracted. We check that the constitutive laws, together with mechanical equilibrium, correctly predict to the observed stress and strain profiles. These profiles exhibit a spatial crossover between an effective linear regime close to the inflater and the truly nonlinear regime away from it. The crossover length diverges at the jamming transition.

Penrose tilings as jammed solids

Penrose tilings form lattices, exhibiting fivefold symmetry and isotropic elasticity, with inhomogeneous coordination much like that of the force networks in jammed systems. Under periodic boundary conditions, their average coordination is exactly four. We study the elastic and vibrational properties of rational approximants to these lattices as a function of unit-cell size N(S) and find that they have of order sqrt[N(S)] zero modes and states of self-stress and yet all their elastic moduli vanish. In their generic form, obtained by randomizing site positions, their elastic and vibrational properties are similar to those of particulate systems at jamming with a nonzero bulk modulus, vanishing shear modulus, and a flat density of states.

Contact changes near jamming

We probe the onset and effect of contact changes in soft harmonic particle packings which are sheared quasistatically. We find that the first contact changes are the creation or breaking of contacts on a single particle. We characterize the critical strain, statistics of breaking versus making a contact, and ratio of shear modulus before and after such events, and explain their finite size scaling relations. For large systems at finite pressure, the critical strain vanishes but the ratio of shear modulus before and after a contact change approaches one: linear response remains relevant in large systems. For finite systems close to jamming the critical strain also vanishes, but here linear response already breaks down after a single contact change.