Jamming transition as probed by quasistatic shear flow

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.

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.

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?

Considerations on the relaxation time in shear-driven jamming

We study the jamming transition in a model of elastic particles under shear at zero temperature, with a focus on the relaxation time τ_{1}. This relaxation time is from two-step simulations where the first step is the ordinary shearing simulation and the second step is the relaxation of the energy after stopping the shearing. τ_{1} is determined from the final exponential decay of the energy. Such relaxations are done with many different starting configurations generated by a long shearing simulation in which the shear variable γ slowly increases. We study the correlations of both τ_{1}, determined from the decay, and the pressure, p_{1}, from the starting configurations as a function of the difference in γ. We find that the correlations of p_{1} are longer lived than the ones of τ_{1} and find that the reason for this is that the individual τ_{1} is controlled both by p_{1} of the starting configuration and a random contribution which depends on the relaxation path length-the average distance moved by the particles during the relaxation. We further conclude that it is γ_{τ}, determined from the correlations of τ_{1}, which is the relevant one when the aim is to generate data that may be used for determining the critical exponent that characterizes the jamming transition.

Fabric-based jamming phase diagram for frictional granular materials

A jamming phase diagram maps the phase states of granular materials to their intensive properties such as shear stress and density (or packing fraction). We investigate how different phases in a jamming phase diagram of granular materials are related to their fabric structure via three-dimensional discrete element method simulations. Constant-volume quasi-static simple shear tests ensuring uniform shear strain field are conducted on bi-disperse spherical frictional particles. Specimens with different initial solid fractions are sheared until reaching steady state at a large shear strain (200%). The jamming threshold in terms of stress, non-rattler fraction, and coordination numbers (Z's) of different contact networks is discussed. The evolution of fabric anisotropy (F) of each contact network during shearing is also examined. By plotting the fabric data in the F-Z space, a unique critical fabric surface (CFS) becomes apparent across all specimens, irrespective of their initial phase states. Through the correlation of this CFS with fabric signals corresponding to jamming transitions, we introduce a novel jamming phase diagram in the fabric F-Z space, offering a convenient approach to distinguish the various phases of granular materials solely through the direct observation of geometrical arrangements of particles. This jamming phase diagram underscores the importance of the microstructure underlying the conventional jamming phenomenon and introduces a novel standpoint for interpreting the phase transitions of granular materials that have been exposed to processes such as compaction, shearing, and other complex loading histories.

Shear flow of non-Brownian rod-sphere mixtures near jamming

We use the discrete element method, taking particle contact and hydrodynamic lubrication into account, to unveil the shear rheology of suspensions of frictionless non-Brownian rods in the dense packing fraction regime. We find that, analogously to the random close packing volume fraction, the shear-driven jamming point of this system varies in a nonmonotonic fashion as a function of the rod aspect ratio. The latter strongly influences how the addition of rodlike particles affects the rheological response of a suspension of frictionless non-Brownian spheres to an external shear flow. At fixed values of the total (rods plus spheres) packing fraction, the viscosity of the suspension is reduced by the addition of "short"(≤2) rods but is instead increased by the addition of "long"(≥2) rods. A mechanistic interpretation is provided in terms of packing and excluded-volume arguments.

Scaling relationships between viscosity and diffusivity in shear-thickening suspensions

Dense suspensions often exhibit a dramatic response to large external deformation. The recent body of work has related this behavior to transition from an unconstrained lubricated state to a constrained frictional state. Here, we use numerical simulations to study the flow behavior and shear-induced diffusion of frictional non-Brownian spheres in two dimensions under simple shear flow. We first show that both viscosity η and diffusivity D/ of the particles increase under characteristic shear stress, which is associated with lubrication to frictional transition. Subsequently, we propose a one-to-one relationship between viscosity and diffusivity using the length scale ξ associated with the size of collective motions (rigid clusters) of the particles. We demonstrate that η and D/ are controlled by ξ in two distinct flow regimes, i.e. in the frictionless and frictional states, where the one-to-one relationship is described as a crossover from D/ ∼ η (frictionless) to η1/3 (frictional). We also confirm that the proposed power laws are insensitive to the interparticle friction and system size.

Slow and fast particles in shear-driven jamming: Critical behavior

We do extensive simulations of a simple model of shear-driven jamming in two dimensions to determine and analyze the velocity distribution at different densities ϕ around the jamming density ϕ_{J} and at different low shear strain rates, γ[over ̇]. We then find that the velocity distribution is made up of two parts which are related to two different physical processes which we call the slow process and the fast process as they are dominated by the slower and the faster particles, respectively. Earlier scaling analyses have shown that the shear viscosity η, which diverges as the jamming density is approached from below, consists of two different terms, and we present strong evidence that these terms are related to the two different processes: the leading divergence is due to the fast process, whereas the correction-to-scaling term is due to the slow process. The analysis of the slow process is possible thanks to the observation that the velocity distribution for different γ[over ̇] and ϕ at and around the shear-driven jamming transition has a peak at low velocities and that the distribution has a constant shape up to and slightly above this peak. We then find that it is possible to express the contribution to the shear viscosity due to the slow process in terms of height and position of the peak in the velocity distribution and find that this contribution matches the correction-to-scaling term, determined through a standard critical scaling analysis. A further observation is that the collective particle motion is dominated by the slow process. In contrast to the usual picture in critical phenomena with a direct link between the diverging correlation length and a diverging order parameter, we find that correlations and shear viscosity decouple since they are controlled by different sets of particles and that shear-driven jamming is thus an unusual kind of critical phenomenon.

Local and global measures of the shear moduli of jammed disk packings

Strain-controlled isotropic compression gives rise to jammed packings of repulsive, frictionless disks with either positive or negative global shear moduli. We carry out computational studies to understand the contributions of the negative shear moduli to the mechanical response of jammed disk packings. We first decompose the ensemble-averaged, global shear modulus as 〈G〉=(1-F_{-})〈G_{+}〉+F_{-}〈G_{-}〉, where F_{-} is the fraction of jammed packings with negative shear moduli and 〈G_{+}〉 and 〈G_{-}〉 are the average values from packings with positive and negative moduli, respectively. We show that 〈G_{+}〉 and 〈|G_{-}|〉 obey different power-law scaling relations above and below pN^{2}∼1. For pN^{2}>1, both 〈G_{+}〉N and 〈|G_{-}|〉N∼(pN^{2})^{β}, where β∼0.5 for repulsive linear spring interactions. Despite this, 〈G〉N∼(pN^{2})^{β^{'}} with β^{'}≳0.5 due to the contributions from packings with negative shear moduli. We show further that the probability distribution of global shear moduli P(G) collapses at fixed pN^{2} and different values of p and N. We calculate analytically that P(G) is a Γ distribution in the pN^{2}≪1 limit. As pN^{2} increases, the skewness of P(G) decreases and P(G) becomes a skew-normal distribution with negative skewness in the pN^{2}≫1 limit. We also partition jammed disk packings into subsystems using Delaunay triangulation of the disk centers to calculate local shear moduli. We show that the local shear moduli defined from groups of adjacent triangles can be negative even when G>0. The spatial correlation function of local shear moduli C(r[over ⃗]) displays weak correlations for pn_{sub}^{2}<10^{-2}, where n_{sub} is the number of particles within each subsystem. However, C(r[over ⃗]) begins to develop long-ranged spatial correlations with fourfold angular symmetry for pn_{sub}^{2}≳10^{-2}.

Comparison of compression versus shearing near jamming, for a simple model of athermal frictionless disks in suspension

Using a simplified model for a non-Brownian suspension, we numerically study the response of athermal, overdamped, frictionless disks in two dimensions to isotropic and uniaxial compression, as well as to pure and simple shearing, all at finite constant strain rates ε[over ̇]. We show that isotropic and uniaxial compression result in the same jamming packing fraction ϕ_{J}, while pure-shear- and simple-shear-induced jamming occurs at a slightly higher ϕ_{J}^{*}, consistent with that found previously for simple shearing. A critical scaling analysis of pure shearing gives critical exponents consistent with those previously found for both isotropic compression and simple shearing. Using orientational order parameters for contact bond directions, we compare the anisotropy of the force and contact networks at both lowest nematic order, as well as higher 2n-fold order.

Orientational-induced strain hardening of axisymmetric grains

The rheological response of oriented axisymmetric grains has additional degrees of complexity associated with their microstructure orientation. These additional kinematic degrees of freedom that give rise to complex transient macroscale rheological responses are not well understood. In this Letter, we study the rheology of axisymmetric grains subjected to transient flow. We identify strong coupling between the microstructure rearrangement and strain hardening which, under certain conditions, can yield jamming. We identify the critical conditions corresponding to jamming and the dependency on the shape of the grains. It is shown that this is a particular form of jamming that is directional in nature, since unjamming occurs if the shear direction is reversed.

Relaxation times, rheology, and finite size effects for non-Brownian disks in two dimensions

We carry out overdamped simulations in a simple model of jamming-a collection of bidisperse soft core frictionless disks in two dimensions-with the aim to explore the finite size dependence of different quantities, both the relaxation time obtained from the relaxation of the energy and the pressure equivalent of the shear viscosity. The motivation for the paper is the observation [Nishikawa et al., J. Stat. Phys. 182, 37 (2021)0022-471510.1007/s10955-021-02710-8] that there are finite size effects in the relaxation time, τ, that give problems in the determination of the critical divergence, and the claim that this is due to a finite size dependence, τ∼lnN, which makes τ an ill-defined quantity. Beside analyses to determine the relaxation time for the whole system we determine particle relaxation times which allow us to determine both histograms of particle relaxation times and the average particle relaxation times-two quantities that are very useful for the analyses. The starting configurations for the relaxation simulations are of two different kinds-completely random or taken from steady shearing simulations-and we find that the difference between these two cases are bigger than previously noted and that the observed problems in the determination of the critical divergence obtained when starting from random configurations are not present when instead starting the relaxations from shearing configurations. We also argue that the the effect that causes the lnN dependence is not as problematic as asserted. When it comes to the finite size dependence of the pressure equivalent of the shear viscosity we find that our data don't give support for the claimed strong finite size dependence, but also that the finite size dependence is at odds with what one would normally expect for a system with a diverging correlation length, and that this calls for an alternative understanding of the phenomenon of shear-driven jamming.

Universality of stress-anisotropic and stress-isotropic jamming of frictionless spheres in three dimensions: Uniaxial versus isotropic compression

We numerically study a three-dimensional system of athermal, overdamped, frictionless spheres, using a simplified model for a non-Brownian suspension. We compute the bulk viscosity under both uniaxial and isotropic compression as a means to address the question of whether stress-anisotropic and stress-isotropic jamming are in the same critical universality class. Carrying out a critical scaling analysis of the system pressure p, shear stress σ, and macroscopic friction μ=σ/p, as functions of particle packing fraction ϕ and compression rate ε[over ̇], we find good agreement for all critical parameters comparing the isotropic and anisotropic cases. In particular, we determine that the bulk viscosity diverges as p/ε[over ̇]∼(ϕ_{J}-ϕ)^{-β}, with β=3.36±0.09, as jamming is approached from below. We further demonstrate that the average contact number per particle Z can also be written in a scaling form as a function of ϕ and ε[over ̇]. Once again, we find good agreement between the uniaxial and isotropic cases. We compare our results to prior simulations and theoretical predictions.

Moduli and modes in the Mikado model

We determine how low frequency vibrational modes control the elastic shear modulus of Mikado networks, a minimal mechanical model for semi-flexible fiber networks. From prior work it is known that when the fiber bending modulus is sufficiently small, (i) the shear modulus of 2D Mikado networks scales as a power law in the fiber line density, G ∼ ρα+1, and (ii) the networks also possess an anomalous abundance of soft (low-frequency) vibrational modes with a characteristic frequency ωκ ∼ ρβ/2. While it has been suggested that α and β are identical, the preponderance of evidence indicates that α is larger than theoretical predictions for β. We resolve this inconsistency by measuring the vibrational density of states in Mikado networks for the first time. Supported by these results, we then demonstrate analytically that α = β + 1. In so doing, we uncover new insights into the coupling between soft modes and shear, as well as the origin of the crossover from bending- to stretching-dominated response.

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.

Nonaffine displacements below jamming under athermal quasistatic compression

Critical properties of frictionless spherical particles below jamming are studied using extensive numerical simulations, paying particular attention to the nonaffine part of the displacements during the athermal quasistatic compression. It is shown that the squared norm of the nonaffine displacement exhibits a power-law divergence toward the jamming transition point. A possible connection between this critical exponent and that of the shear viscosity is discussed. The participation ratio of the displacements vanishes in the thermodynamic limit at the transition point, meaning that the nonaffine displacements are localized marginally with a fractal dimension. Furthermore, the distribution of the displacement is shown to have a power-law tail, the exponent of which is related to the fractal dimension.

Critical scaling of compression-driven jamming of athermal frictionless spheres in suspension

We study numerically a system of athermal, overdamped, frictionless spheres, as in a non-Brownian suspension, in two and three dimensions. Compressing the system isotropically at a fixed rate ε[over ̇], we investigate the critical behavior at the jamming transition. The finite compression rate introduces a control timescale, which allows one to probe the critical timescale associated with jamming. As was found previously for steady-state shear-driven jamming, we find for compression-driven jamming that pressure obeys a critical scaling relation as a function of packing fraction ϕ and compression rate ε[over ̇], and that the bulk viscosity p/ε[over ̇] diverges upon jamming. A scaling analysis determines the critical exponents associated with the compression-driven jamming transition. Our results suggest that stress-isotropic, compression-driven jamming may be in the same universality class as stress-anisotropic, shear-driven jamming.

Dilatancy, shear jamming, and a generalized jamming phase diagram of frictionless sphere packings

Granular packings display the remarkable phenomenon of dilatancy, wherein their volume increases upon shear deformation. Conventional wisdom and previous results suggest that dilatancy, also being the related phenomenon of shear-induced jamming, requires frictional interactions. Here, we show that the occurrence of isotropic jamming densities φ_j above the minimal density (or the J-point density) φ_J leads both to the emergence of shear-induced jamming and dilatancy in frictionless packings. Under constant pressure shear, the system evolves into a steady-state at sufficiently large strains, whose density only depends on the pressure and is insensitive to the initial jamming density φ_j. In the limit of vanishing pressure, the steady-state exhibits critical behavior at φ_J. While packings with different φ_j values display equivalent scaling properties under compression, they exhibit striking differences in rheological behaviour under shear. The yield stress under constant volume shear increases discontinuously with density when φ_j > φ_J, contrary to the continuous behaviour in generic packings that jam at φ_J. Our results thus lead to a more coherent, generalised picture of jamming in frictionless packings, which also have important implications on how dilatancy is understood in the context of frictional granular matter.

Dynamic length scales in athermal, shear-driven jamming of frictionless disks in two dimensions

We carry out numerical simulations of athermally sheared, bidisperse, frictionless disks in two dimensions. From an appropriately defined velocity correlation function, we determine that there are two diverging length scales, ξ and ℓ, as the jamming transition is approached. We analyze our results using a critical scaling ansatz for the correlation function and argue that the more divergent length ℓ is a consequence of a dangerous irrelevant scaling variable and that it is ξ, which is the correlation length that determines the divergence of the system viscosity as jamming is approached from below in the liquid phase. We find that ξ∼(ϕ_{J}-ϕ)^{-ν} diverges with the critical exponent ν=1. We provide evidence that ξ measures the length scale of fluctuations in the rotation of the particle velocity field, while ℓ measures the length scale of fluctuations in the divergence of the velocity field.

Translational and rotational velocities in shear-driven jamming of ellipsoidal particles

We study shear-driven jamming of ellipsoidal particles at zero temperature with a focus on the microscopic dynamics. We find that a change from spherical particles to ellipsoids with aspect ratio α=1.02 gives dramatic changes of the microscopic dynamics with much lower translational velocities and a new role for the rotations. Whereas the velocity difference at contacts-and thereby the dissipation-in collections of spheres is dominated by the translational velocities and reduced by the rotations, the same quantity is in collections of ellipsoids instead totally dominated by the rotational velocities. By also examining the effect of different aspect ratios we find that the examined quantities show either a peak or a change in slope at α≈1.2, which thus gives evidence for a crossover between different regions of low and high aspect ratio.

Packings of frictionless spherocylinders

We present simulation results on the properties of packings of frictionless spherocylindrical particles. Starting from a random distribution of particles in space, a packing is produced by minimizing the potential energy of interparticle contacts until a force-equilibrated state is reached. For different particle aspect ratios α=10⋯40, we calculate contacts z, pressure as well as bulk and shear modulus. Most important is the fraction f_{0}(α) of spherocylinders with contacts at both ends, as it governs the jamming threshold z_{c}(α)=8+2f_{0}(α). These results highlight the important role of the axial "sliding" degree of freedom of a spherocylinder, which is a zero-energy mode but only if no end contacts are present.

Stress Relaxation above and below the Jamming Transition

We numerically investigate stress relaxation in soft athermal disks to reveal critical slowing down when the system approaches the jamming point. The exponents describing the divergence of the relaxation time differ dramatically depending on whether the transition is approached from the jammed or unjammed phase. This contrasts sharply with conventional dynamic critical scaling scenarios, where a single exponent characterizes both sides. We explain this surprising difference in terms of the vibrational density of states, which is a key ingredient of linear viscoelastic theory. The vibrational density of states exhibits an extra slow mode that emerges below jamming, which we utilize to demonstrate the anomalous exponent below jamming.

Shear-driven flow of athermal, frictionless, spherocylinder suspensions in two dimensions: Particle rotations and orientational ordering

We use numerical simulations to study the flow of a bidisperse mixture of athermal, frictionless, soft-core two-dimensional spherocylinders driven by a uniform steady-state simple shear applied at a fixed volume and a fixed finite strain rate γ[over ̇]. Energy dissipation is via a viscous drag with respect to a uniformly sheared host fluid, giving a simple model for flow in a non-Brownian suspension with Newtonian rheology. Considering a range of packing fractions ϕ and particle asphericities α at small γ[over ̇], we study the angular rotation θ[over ̇]_{i} and the nematic orientational ordering S_{2} of the particles induced by the shear flow, finding a nonmonotonic behavior as the packing ϕ is varied. We interpret this nonmonotonic behavior as a crossover from dilute systems at small ϕ, where single-particle-like behavior occurs, to dense systems at large ϕ, where the geometry of the dense packing dominates and a random Poisson-like process for particle rotations results. We also argue that the finite nematic ordering S_{2} is a consequence of the shearing serving as an ordering field, rather than a result of long-range cooperative behavior among the particles. We arrive at these conclusions by consideration of (i) the distribution of waiting times for a particle to rotate by π, (ii) the behavior of the system under pure, as compared to simple, shearing, (iii) the relaxation of the nematic order parameter S_{2} when perturbed away from the steady state, and (iv) by construction, a numerical mean-field model for the rotational motion of a particle. Our results also help to explain the singular behavior observed when taking the α→0 limit approaching circular disks.

Rheology in dense assemblies of spherocylinders: Frictional vs. frictionless

Using molecular dynamics simulations, we study the steady shear flow of dense assemblies of anisotropic spherocylindrical particles of varying aspect ratios. Comparing frictionless and frictional particles we discuss the specific role of frictional inter-particle forces for the rheological properties of the system. In the frictional system we evidence a shear-thickening regime, similar to that for spherical particles. Furthermore, friction suppresses the alignment of the spherocylinders along the flow direction. Finally, the jamming density in frictional systems is rather insensitive to variations in aspect ratio, quite contrary to what is known from frictionless systems.

Soft-grain compression: Beyond the jamming point

We present the experimental studies of highly strained soft bidisperse granular systems made of hyperelastic and plastic particles. We explore the behavior of granular matter deep in the jammed state from local field measurement from the grain scale to the global scale. By means of a dedicated digital image correlation code and an accurate image recording method, we measure for each compression step the evolution of the particle geometries and their right Cauchy-Green strain tensor fields. We analyze the evolution of the usual macroscopic observables (stress, packing fraction, coordination, fraction of nonrattlers, etc.) along the compression process through the jamming point and far beyond. Analyzing the evolution of the local strain statistics, we evidence a crossover in the material behavior deep in the jammed state for both sorts of particles. We show that this crossover is due to a competition between material compression, dilation, and shear, so its position depends on the particle material. We argue that the strain field is a reliable observable to describe the evolution of a granular system through the jamming transition and deep in the dense packing state whatever the material behavior.

Shear-driven flow of athermal, frictionless, spherocylinder suspensions in two dimensions: Stress, jamming, and contacts

We use numerical simulations to study the flow of a bidisperse mixture of athermal, frictionless, soft-core two-dimensional spherocylinders driven by a uniform steady-state shear strain applied at a fixed finite rate. Energy dissipation occurs via a viscous drag with respect to a uniformly sheared host fluid, giving a simple model for flow in a non-Brownian suspension and resulting in a Newtonian rheology. We study the resulting pressure p and deviatoric shear stress σ of the interacting spherocylinders as a function of packing fraction ϕ, strain rate γ[over ̇], and a parameter α that measures the asphericity of the particles; α is varied to consider the range from nearly circular disks to elongated rods. We consider the direction of anisotropy of the stress tensor, the macroscopic friction μ=σ/p, and the divergence of the transport coefficient η_{p}=p/γ[over ̇] as ϕ is increased to the jamming transition ϕ_{J}. From a phenomenological analysis of Herschel-Bulkley rheology above jamming, we estimate ϕ_{J} as a function of asphericity α and show that the variation of ϕ_{J} with α is the main cause for differences in rheology as α is varied; when plotted as ϕ/ϕ_{J}, rheological curves for different α qualitatively agree. However, a detailed scaling analysis of the divergence of η_{p} for our most elongated particles suggests that the jamming transition of spherocylinders may be in a different universality class than that of circular disks. We also compute the number of contacts per particle Z in the system and show that the value at jamming Z_{J} is a nonmonotonic function of α that is always smaller than the isostatic value. We measure the probability distribution of contacts per unit surface length P(ϑ) at polar angle ϑ with respect to the spherocylinder spine and find that as α→0 this distribution seems to diverge at ϑ=π/2, giving a finite limiting probability for contacts on the vanishingly small flat sides of the spherocylinder. Finally, we consider the variation of the average contact force as a function of location on the particle surface.

Biomechanical Feedback Strengthens Jammed Cellular Packings

Growth in confined spaces can drive cellular populations through a jamming transition from a fluidlike state to a solidlike state. Experiments have found that jammed budding yeast populations can build up extreme compressive pressures (over 1 MPa), which in turn feed back onto cellular physiology by slowing or even stalling cell growth. Using numerical simulations, we investigate how this feedback impacts the mechanical properties of model jammed cell populations. We find that feedback directs growth toward poorly coordinated regions, resulting in an excess number of cell-cell contacts that rigidify cell packings. Cell packings possess anomalously large shear and bulk moduli that depend sensitively on the strength of feedback. These results demonstrate that mechanical feedback on the single-cell level is a simple mechanism by which living systems may tune their population-level mechanical properties.

Avalanche Interpretation of the Power-Law Energy Spectrum in Three-Dimensional Dense Granular Flow

Turbulence is ubiquitous in nonequilibrium systems, and it has been noted that even dense granular flows exhibit characteristics that are typical of turbulent flow, such as the power-law energy spectrum. However, studies on the turbulentlike behavior of granular flows are limited to two-dimensional (2D) flow. We demonstrate that the statistics in three-dimensional (3D) flow are qualitatively different from those in 2D flow. We also elucidate that avalanche dynamics can explain this dimensionality dependence. Moreover, we define clusters of collectively moving particles that are equivalent to vortex filaments. The clusters unveil complicated structures in 3D flows that are absent in 2D flows.

Dimensionality and Viscosity Exponent in Shear-driven Jamming

Collections of bidisperse frictionless particles at zero temperature in three dimensions are simulated with a shear-driven dynamics with the aim to compare with the behavior in two dimensions. Contrary to the prevailing picture, and in contrast to results from isotropic jamming from compression or quench, we find that the critical exponents in three dimensions are different from those in two dimensions and conclude that shear-driven jamming in two and three dimensions belong to different universality classes.

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.

Criticality of the zero-temperature jamming transition probed by self-propelled particles

We perform simulations of athermal systems of self-propelled particles (SPPs) interacting via harmonic repulsion in the vicinity of the zero-temperature jamming transition at point J. Every particle is propelled by a constant force f pointing to a randomly assigned and fixed direction. When f is smaller than the yield force f_y, the system is statically jammed. We find that f_y increases with packing fraction and exhibits finite size scaling, implying the criticality of point J. When f > f_y, SPPs flow forever and their velocities satisfy the k-Gamma distribution. Velocity distributions at various packing fractions and f collapse when the particle velocity is scaled by the average velocity v[combining macron], suggesting that v[combining macron] is a reasonable quantity to characterize the response to f. We thus define a response function R(ϕ,f) = v[combining macron](ϕ,f)/f. The function exhibits critical scaling nicely, implying again the criticality of point J. Our analysis and results indicate that systems of SPPs behave analogically to sheared systems, although their driving mechanisms are apparently distinct.

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.

On the apparent yield stress in non-Brownian magnetorheological fluids

We use simulations to probe the flow properties of dense two-dimensional magnetorheological fluids. Prior results from both experiments and simulations report that the shear stress σ scales with strain rate [small gamma, Greek, dot above] as σ ∼ [small gamma, Greek, dot above]^1-Δ, with values of the exponent ranging between 2/3 < Δ ≤ 1. However it remains unclear what properties of the system select the value of Δ, and in particular under what conditions the system displays a yield stress (Δ = 1). To address these questions, we perform simulations of a minimalistic model system in which particles interact via long ranged magnetic dipole forces, finite ranged elastic repulsion, and viscous damping. We find a surprising dependence of the apparent exponent Δ on the form of the viscous force law. For experimentally relevant values of the volume fraction ϕ and the dimensionless Mason number Mn (which quantifies the competition between viscous and magnetic stresses), models using a Stokes-like drag force show Δ ≈ 0.75 and no apparent yield stress. When dissipation occurs at the contact, however, a clear yield stress plateau is evident in the steady state flow curves. In either case, increasing ϕ towards the jamming transition suffices to induce a yield stress. We relate these qualitatively distinct flow curves to clustering mechanisms at the particle scale. For Stokes-like drag, the system builds up anisotropic, chain-like clusters as Mn tends to zero (vanishing strain rate and/or high field strength). For contact damping, by contrast, there is a second clustering mechanism due to inelastic collisions.

Anisotropic decay of the energy spectrum in two-dimensional dense granular flows

We study anisotropic collective motions of two-dimensional granular particles under simple shear deformations. Employing molecular-dynamics simulations of large system sizes, we find that anisotropic fluidized bands develop in the system yielding under quasistatic deformations, where the spectrum of nonaffine velocities, which is associated with the energy spectrum for turbulent flows, exhibits a quadrupole structure. To explain theoretically the anisotropic spectrum, we derive hydrodynamic modes from a continuum model of dense granular materials, where we find that fluidized bands are caused by long-lived hydrodynamic fluctuations characterized by compressibility.

Slow creep in soft granular packings

Transient creep mechanisms in soft granular packings are studied numerically using a constant pressure and constant stress simulation method. Rapid compression followed by slow dilation is predicted on the basis of a logarithmic creep phenomenon. Characteristic scales of creep strain and time exhibit a power-law dependence on jamming pressure, and they diverge at the jamming point. Microscopic analysis indicates the existence of a correlation between rheology and nonaffine fluctuations. Localized regions of large strain appear during creep and grow in magnitude and size at short times. At long times, the spatial structure of highly correlated local deformation becomes time-invariant. Finally, a microscale connection between local rheology and local fluctuations is demonstrated in the form of a linear scaling between granular fluidity and nonaffine velocity.

Shear banding, discontinuous shear thickening, and rheological phase transitions in athermally sheared frictionless disks

We report on numerical simulations of simple models of athermal, bidisperse, soft-core, massive disks in two dimensions, as a function of packing fraction ϕ, inelasticity of collisions as measured by a parameter Q, and applied uniform shear strain rate γ[over ̇]. Our particles have contact interactions consisting of normally directed elastic repulsion and viscous dissipation, as well as tangentially directed viscous dissipation, but no interparticle Coulombic friction. Mapping the phase diagram in the (ϕ,Q) plane for small γ[over ̇], we find a sharp first-order rheological phase transition from a region with Bagnoldian rheology to a region with Newtonian rheology, and show that the system is always Newtonian at jamming. We consider the rotational motion of particles and demonstrate the crucial importance that the coupling between rotational and translational degrees of freedom has on the phase structure at small Q (strongly inelastic collisions). At small Q, we show that, upon increasing γ[over ̇], the sharp Bagnoldian-to-Newtonian transition becomes a coexistence region of finite width in the (ϕ,γ[over ̇]) plane, with coexisting Bagnoldian and Newtonian shear bands. Crossing this coexistence region by increasing γ[over ̇] at fixed ϕ, we find that discontinuous shear thickening can result if γ[over ̇] is varied too rapidly for the system to relax to the shear-banded steady state corresponding to the instantaneous value of γ[over ̇].

Characterizing rare fluctuations in soft particulate flows

Soft particulate media include a wide range of systems involving athermal dissipative particles both in non-living and biological materials. Characterization of flows of particulate media is of great practical and theoretical importance. A fascinating feature of these systems is the existence of a critical rigidity transition in the dense regime dominated by highly intermittent fluctuations that severely affects the flow properties. Here, we unveil the underlying mechanisms of rare fluctuations in soft particulate flows. We find that rare fluctuations have different origins above and below the critical jamming density and become suppressed near the jamming transition. We then conjecture a time-independent local fluctuation relation, which we verify numerically, and that gives rise to an effective temperature. We discuss similarities and differences between our proposed effective temperature with the conventional kinetic temperature in the system by means of a universal scaling collapse.

Soft particulate flows such as granular media are prone to fluctuations like jamming and avalanches. Here Rahbari et al. consider the statistics of rare fluctuations to identify an effective temperature which, unlike previous ones, is valid for packing fractions both near and far from the jamming point.

Shear Yielding and Shear Jamming of Dense Hard Sphere Glasses

We investigate the response of dense hard sphere glasses to a shear strain in a wide range of pressures ranging from the glass transition to the infinite-pressure jamming point. The phase diagram in the density-shear strain plane is calculated analytically using the mean-field infinite-dimensional solution. We find that just above the glass transition, the glass generically yields at a finite shear strain. The yielding transition in the mean-field picture is a spinodal point in presence of disorder. At higher densities, instead, we find that the glass generically jams at a finite shear strain: the jamming transition prevents yielding. The shear yielding and shear jamming lines merge in a critical point, close to which the system yields at extremely large shear stress. Around this point, highly nontrivial yielding dynamics, characterized by system-spanning disordered fractures, is expected.

Enstrophy cascades in two-dimensional dense granular flows

Employing two-dimensional molecular dynamics simulations of dense granular materials under simple shear deformations, we investigate vortex structures of particle rearrangements. Introducing vorticity fields as a measure of local spinning motions of the particles, we observe their heterogeneous distributions, where statistics of vorticity fields exhibit the highly non-Gaussian behavior and typical domain sizes of vorticity fields significantly increase if the system is yielding under quasistatic deformations. In such dense granular flows, a power-law decay of vorticity spectra can be observed at mesoscopic scale, implying anomalous local structures of kinetic energy dissipation. We explain the power-law decay, or enstrophy cascades in dense granular materials, by a dimensional analysis, where the dependence of vorticity spectra not only on the wave number, but also on the shear rate, is well explained. From our dimensional analyses, the scaling of granular temperature and rotational kinetic energy is also predicted.

Self-Driven Jamming in Growing Microbial Populations

In natural settings, microbes tend to grow in dense populations [1-4] where they need to push against their surroundings to accommodate space for new cells. The associated contact forces play a critical role in a variety of population-level processes, including biofilm formation [5-7], the colonization of porous media [8, 9], and the invasion of biological tissues [10-12]. Although mechanical forces have been characterized at the single cell level [13-16], it remains elusive how collective pushing forces result from the combination of single cell forces. Here, we reveal a collective mechanism of confinement, which we call self-driven jamming, that promotes the build-up of large mechanical pressures in microbial populations. Microfluidic experiments on budding yeast populations in space-limited environments show that self-driven jamming arises from the gradual formation and sudden collapse of force chains driven by microbial proliferation, extending the framework of driven granular matter [17-20]. The resulting contact pressures can become large enough to slow down cell growth, to delay the cell cycle in the G1 phase, and to strain or even destroy the microenvironment through crack propagation. Our results suggest that self-driven jamming and build-up of large mechanical pressures is a natural tendency of microbes growing in confined spaces, contributing to microbial pathogenesis and biofouling [21-26].

Phase diagram for inertial granular flows

Flows of hard granular materials depend strongly on the interparticle friction coefficient μ_{p} and on the inertial number I, which characterizes proximity to the jamming transition where flow stops. Guided by numerical simulations, we derive the phase diagram of dense inertial flow of spherical particles, finding three regimes for 10^{-4}≲I≲10^{-1}: frictionless, frictional sliding, and rolling. These are distinguished by the dominant means of energy dissipation, changing from collisional to sliding friction, and back to collisional, as μ_{p} increases from zero at constant I. The three regimes differ in their kinetics and rheology; in particular, the velocity fluctuations and the stress ratio both display nonmonotonic behavior with μ_{p}, corresponding to transitions between the three regimes of flow. We rationalize the phase boundaries between these regimes, show that energy balance yields scaling relations between microscopic properties in each of them, and derive the strain scale at which particles lose memory of their velocity. For the frictional sliding regime most relevant experimentally, we find for I≥10^{-2.5} that the growth of the macroscopic friction μ(I) with I is induced by an increase of collisional dissipation. This implies in that range that μ(I)-μ(0)∼I^{1-2b}, where b≈0.2 is an exponent that characterizes both the dimensionless velocity fluctuations L∼I^{-b} and the density of sliding contacts χ∼I^{b}.

Reversibility and hysteresis of the sharp yielding transition of a colloidal glass under oscillatory shear

The mechanical response of glasses remains challenging to understand. Recent results indicate that the oscillatory rheology of soft glasses is accompanied by a sharp non-equilibrium transition in the microscopic dynamics. Here, we use simultaneous x-ray scattering and rheology to investigate the reversibility and hysteresis of the sharp symmetry change from anisotropic solid to isotropic liquid dynamics observed in the oscillatory shear of colloidal glasses (D. Denisov, M.T. Dang, B. Struth, A. Zaccone, P. Schall, Sci. Rep. 5 14359 (2015)). We use strain sweeps with increasing and decreasing strain amplitude to show that, in analogy with equilibrium transitions, this sharp symmetry change is reversible and exhibits systematic frequency-dependent hysteresis. Using the non-affine response formalism of amorphous solids, we show that these hysteresis effects arise from frequency-dependent non-affine structural cage rearrangements at large strain. These results consolidate the first-order-like nature of the oscillatory shear transition and quantify related hysteresis effects both via measurements and theoretical modelling.

Anomalous energy cascades in dense granular materials yielding under simple shear deformations

By using molecular dynamics (MD) simulations of dense granular particles in two dimensions, we study turbulent-like structures of their non-affine velocities under simple shear deformations. We find that the spectrum of non-affine velocities, introduced as an analog of the energy spectrum for turbulent flows, exhibits the power-law decay if the system is yielding in a quasi-static regime, where large-scale collective motions and inelastic interactions of granular particles are crucial for the anomalous cascade of kinetic energy. Based on hydrodynamic equations of dense granular materials, which include both kinetic and contact contributions in constitutive relations, we derive a theoretical expression for the spectrum, where a good agreement between the result of MD simulations and theoretical prediction is established over a wide range of length scales.

Viscoelastic properties of marginal networks in a solvent

Polymer networks at the margins of mechanical stability are known to be highly sensitive to applied forces and fields and to exhibit an anomalously large resistance to deformation. In this paper, we study the effects of hydrodynamic interactions on the behavior of marginal networks using a hybrid molecular dynamics and multiparticle collision dynamics simulation technique. We examine how the filament and solvent properties affect the response of marginal networks to shear. We find that the stiffening of the network shows a stronger dependence on the shear frequency when hydrodynamic interactions are present than when they are not. The network shear modulus scales as G'∼ω(α(c)), with a critical stiffening exponent α(c) that can be controlled by varying the relative concentrations of the network and the solvent. Our results show that this arises due to the solvent aiding the relaxation of the network and suppressing the network nonaffinity, with the system deforming more affinely when hydrodynamic interactions are maximized.

Length scales in glass-forming liquids and related systems: a review

The central problem in the study of glass-forming liquids and other glassy systems is the understanding of the complex structural relaxation and rapid growth of relaxation times seen on approaching the glass transition. A central conceptual question is whether one can identify one or more growing length scale(s) associated with this behavior. Given the diversity of molecular glass-formers and a vast body of experimental, computational and theoretical work addressing glassy behavior, a number of ideas and observations pertaining to growing length scales have been presented over the past few decades, but there is as yet no consensus view on this question. In this review, we will summarize the salient results and the state of our understanding of length scales associated with dynamical slow down. After a review of slow dynamics and the glass transition, pertinent theories of the glass transition will be summarized and a survey of ideas relating to length scales in glassy systems will be presented. A number of studies have focused on the emergence of preferred packing arrangements and discussed their role in glassy dynamics. More recently, a central object of attention has been the study of spatially correlated, heterogeneous dynamics and the associated length scale, studied in computer simulations and theoretical analysis such as inhomogeneous mode coupling theory. A number of static length scales have been proposed and studied recently, such as the mosaic length scale discussed in the random first-order transition theory and the related point-to-set correlation length. We will discuss these, elaborating on key results, along with a critical appraisal of the state of the art. Finally we will discuss length scales in driven soft matter, granular fluids and amorphous solids, and give a brief description of length scales in aging systems. Possible relations of these length scales with those in glass-forming liquids will be discussed.

Free volume under shear

Using an athermal quasistatic simulation protocol, we study the distribution of free volumes in sheared hard-particle packings close to, but below, the random-close packing threshold. We show that under shear, and independent of volume fraction, the free volumes develop features similar to close-packed systems - particles self-organize in a manner as to mimick the isotropically jammed state. We compare athermally sheared packings with thermalized packings and show that thermalization leads to an erasure of these structural features. The temporal evolution in particular the opening-up and the closing of free-volume patches is associated with the single-particle dynamics, showing a crossover from ballistic to diffusive behavior.

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.

Membrane Protein Structure, Function, and Dynamics: a Perspective from Experiments and Theory

Membrane proteins mediate processes that are fundamental for the flourishing of biological cells. Membrane-embedded transporters move ions and larger solutes across membranes; receptors mediate communication between the cell and its environment and membrane-embedded enzymes catalyze chemical reactions. Understanding these mechanisms of action requires knowledge of how the proteins couple to their fluid, hydrated lipid membrane environment. We present here current studies in computational and experimental membrane protein biophysics, and show how they address outstanding challenges in understanding the complex environmental effects on the structure, function, and dynamics of membrane proteins.

Softening of stressed granular packings with resonant sound waves

We perform numerical simulations of a two-dimensional bidisperse granular packing subjected to both a static confining pressure and a sinusoidal dynamic forcing applied by a wall on one edge of the packing. We measure the response experienced by a wall on the opposite edge of the packing and obtain the resonant frequency of the packing as the static or dynamic pressures are varied. Under increasing static pressure, the resonant frequency increases, indicating a velocity increase of elastic waves propagating through the packing. In contrast, when the dynamic amplitude is increased for fixed static pressure, the resonant frequency decreases, indicating a decrease in the wave velocity. This occurs both for compressional and for shear dynamic forcing and is in agreement with experimental results. We find that the average contact number Zc at the resonant frequency decreases with increasing dynamic amplitude, indicating that the elastic softening of the packing is associated with a reduced number of grain-grain contacts through which the elastic waves can travel. We image the excitations created in the packing and show that there are localized disturbances or soft spots that become more prevalent with increasing dynamic amplitude. Our results are in agreement with experiments on glass bead packings and earth materials such as sandstone and granite and may be relevant to the decrease in elastic wave velocities that has been observed to occur near fault zones after strong earthquakes, in surficial sediments during strong ground motion, and in structures during earthquake excitation.

Universal rescaling of flow curves for yield-stress fluids close to jamming

The experimental flow curves of four different yield-stress fluids with different interparticle interactions are studied near the jamming concentration. By appropriate scaling with the distance to jamming all rheology data can be collapsed onto master curves below and above jamming that meet in the shear-thinning regime and satisfy the Herschel-Bulkley and Cross equations, respectively. In spite of differing interactions in the different systems, master curves characterized by universal scaling exponents are found for the four systems. A two-state microscopic theory of heterogeneous dynamics is presented to rationalize the observed transition from Herschel-Bulkley to Cross behavior and to connect the rheological exponents to microscopic exponents for the divergence of the length and time scales of the heterogeneous dynamics. The experimental data and the microscopic theory are compared with much of the available literature data for yield-stress systems.