Jamming transition as probed by quasistatic shear flow

Relaxation times and rheology in dense athermal suspensions

We study the jamming transition in a model of elastic particles under shear at zero temperature. The key quantity is the relaxation time τ which is obtained by stopping the shearing and letting energy and pressure decay to zero. At many different densities and initial shear rates we do several such relaxations to determine the average τ. We establish that τ diverges with the same exponent as the viscosity and determine another exponent from the relation between τ and the coordination number. Though most of the simulations are done for the model with dissipation due to the motion of particles relative to an affinely shearing substrate, we also examine a model, where the dissipation is instead due to velocity differences of disks in contact, and confirm that the above-mentioned exponent is the same for these two models. We also consider finite size effects on both τ and the coordination number.

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.

Universal and non-universal features in coarse-grained models of flow in disordered solids

We study the two-dimensional (2D) shear flow of amorphous solids within variants of an elastoplastic model, paying particular attention to spatial correlations and time fluctuations of, e.g., local stresses. The model is based on the local alternation between an elastic regime and plastic events during which the local stress is redistributed. The importance of a fully tensorial description of the stress and of the inclusion of (coarse-grained) convection in the model is investigated; scalar and tensorial models yield similar results, while convection enhances fluctuations and breaks the spurious symmetry between the flow and velocity gradient directions, for instance when shear localisation is observed. Besides, correlation lengths measured with diverse protocols are discussed. One class of such correlation lengths simply scale with the spacing between homogeneously distributed, simultaneous plastic events. This leads to a scaling of the correlation length with the shear rate as γ̇(-1/2) in 2D in the athermal regime, regardless of the details of the model. The radius of the cooperative disk, defined as the near-field region in which plastic events induce a stress redistribution that is not amenable to a mean-field treatment, notably follows this scaling. On the other hand, the cooperative volume measured from the four-point stress susceptibility and its dependence on the system size and the shear rate are model-dependent.

Measuring the configurational temperature of a binary disc packing

Jammed packings of granular materials differ from systems normally described by statistical mechanics in that they are athermal. In recent years a statistical mechanics of static granular media has emerged where the thermodynamic temperature is replaced by a configurational temperature X which describes how the number of mechanically stable configurations depends on the volume. Four different methods have been suggested to measure X. Three of them are computed from properties of the Voronoi volume distribution, the fourth takes into account the contact number and the global volume fraction. This paper answers two questions using experimental binary disc packings: first we test if the four methods to measure compactivity provide identical results when applied to the same dataset. We find that only two of the methods agree quantitatively. This implies that at least two of the four methods are wrong. Secondly, we test if X is indeed an intensive variable; this becomes true only for samples larger than roughly 200 particles. This result is shown to be due to recently measured correlations between the particle volumes [Zhao et al., Europhys. Lett., 2012, 97, 34004].

Aspects of jamming in two-dimensional athermal frictionless systems

In this work we provide an overview of jamming transitions in two dimensional systems focusing on the limit of frictionless particle interactions in the absence of thermal fluctuations. We first discuss jamming in systems with short range repulsive interactions, where the onset of jamming occurs at a critical packing density and where certain quantities show a divergence indicative of critical behavior. We describe how aspects of the dynamics change as the jamming density is approached and how these dynamics can be explored using externally driven probes. Different particle shapes can produce jamming densities much lower than those observed for disk-shaped particles, and we show how jamming exhibits fragility for some shapes while for other shapes this is absent. Next we describe the effects of long range interactions and jamming behavior in systems such as charged colloids, vortices in type-II superconductors, and dislocations. We consider the effect of adding obstacles to frictionless jamming systems and discuss connections between this type of jamming and systems that exhibit depinning transitions. Finally, we discuss open questions such as whether the jamming transition in all these different systems can be described by the same or a small subset of universal behaviors, as well as future directions for studies of jamming transitions in two dimensional systems, such as jamming in self-driven or active matter systems.

Jamming of frictional particles: a nonequilibrium first-order phase transition

We propose a phase diagram for the shear flow of dry granular particles in two dimensions based on simulations and a phenomenological Landau theory for a nonequilibrium first-order phase transition. Our approach incorporates both frictional as well as frictionless particles. The most important feature of the frictional phase diagram is reentrant flow and a critical jamming point at finite stress. In the frictionless limit the regime of reentrance vanishes and the jamming transition is continuous with a critical point at zero stress. The jamming phase diagrams derived from the model agree with the experiments of Bi et al. [Nature (London) 480, 355 (2011)] and brings together previously conflicting numerical results.

Rheology near jamming: the influence of lubrication forces

We study, by computer simulations, the roles of different dissipation forces in the rheological properties of highly dense particle-laden flows. In particular, we are interested in the close-packing limit (jamming) and the question of whether "universal" observables can be identified that do not depend on the details of the dissipation model. To this end, we define a simplified lubrication force and systematically vary the range h(c) of this interaction. For fixed h(c) a crossover is seen from a Newtonian flow regime at small strain rates to inertia-dominated flow at larger strain rates. The same crossover is observed as a function of the lubrication range h(c). At the same time, but only at high densities close to jamming, single-particle velocities as well as local density distributions are unaffected by changes in the lubrication range--they are candidates for universal behavior. At densities away from jamming, this invariance is lost: short-range lubrication forces lead to pronounced particle clustering, while longer-ranged lubrication does not. These findings highlight the importance of "geometric" packing constraints for particle motion--independent of the specific dissipation model. With the free volume vanishing at random close packing, particle motion is more and more constrained by the ever smaller amount of free space. On the other hand, macroscopic rheological observables as well as higher-order correlation functions retain the variability of the underlying dissipation model.

Size and density avalanche scaling near jamming

The current microscopic picture of plasticity in amorphous materials assumes local failure events to produce displacement fields complying with linear elasticity. Indeed, the flow properties of nonaffine systems, such as foams, emulsions and granular materials close to jamming, that produce a fluctuating displacement field when failing, are still controversial. Here we show, via a thorough numerical investigation of jammed materials, that nonaffinity induces a critical scaling of the flow properties dictated by the distance to the jamming point. We rationalize this critical behavior by introducing a new universal jamming exponent and hyperscaling relationships, and we use these results to describe the volume fraction dependence of the friction coefficient.

Dynamics of sheared ellipses and circular disks: effects of particle shape

Much recent effort has focused on glassy and jamming properties of spherical particles. Very little is known about such phenomena for nonspherical particles, and we take a first step by studying ellipses. We find important differences between the dynamical and structural properties of disks and two-dimensional ellipses subject to continuous Couette shear. In particular, ellipses show slow dynamical evolution, without a counterpart in disks, in the mean velocity, local density, orientational order, and local stress. Starting from an unjammed state, ellipses can first jam under shear, and then slowly unjam. The slow unjamming process is understood as a result of gradual changes in their orientations, leading to a denser packing. For disks, the rotation of particles only contributes to the relaxation of frictional forces, and hence, does not significantly cause structural changes. For the shear-jammed states, the global building up and relaxation of stress, which occurs in the form of stress avalanches, is qualitatively different for disks and ellipses, and is manifested by different forms of rate dependence for ellipses versus disks. Unlike the weak rate dependence typical for many granular systems, ellipses show power-law dependence on the shearing rate Ω.

Finite size analysis of zero-temperature jamming transition under applied shear stress by minimizing a thermodynamic-like potential

By finding local minima of a thermodynamic-like potential, we generate jammed packings of frictionless spheres under constant shear stress σ and obtain the yield stress σy by sampling the potential energy landscape. For three-dimensional systems with harmonic repulsion, σy satisfies the finite size scaling with the limiting scaling relation σy∼ϕ-ϕc,∞, where ϕc,∞ is the critical volume fraction of the jamming transition at σ=0 in the thermodynamic limit. The finite size scaling implies a length ξ∼(ϕ-ϕc,∞)-ν with ν=0.81±0.05, which turns out to be a robust and universal length scale exhibited as well in the finite size scaling of multiple quantities measured without shear and independent of particle interaction. Moreover, comparison between our new approach and quasistatic shear reveals that quasistatic shear tends to explore low-energy states.

Length scales and self-organization in dense suspension flows

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

Pressure distribution and critical exponent in statically jammed and shear-driven frictionless disks

We numerically study the distributions of global pressure that are found in ensembles of statically jammed and quasistatically sheared systems of bidisperse, frictionless disks at fixed packing fraction ϕ in two dimensions. We use these distributions to address the question of how pressure increases as ϕ increases above the jamming point ϕ(J), p ∼ |ϕ-ϕ(J)(y). For statically jammed ensembles, our results are consistent with the exponent y being simply related to the power law of the interparticle soft-core interaction. For sheared systems, however, the value of y is consistent with a nontrivial value, as found previously in rheological simulations.

Shear thickening in granular suspensions: interparticle friction and dynamically correlated clusters

We consider the shear rheology of concentrated suspensions of non-Brownian frictional particles. The key result of our study is the emergence of a pronounced shear-thickening regime, where frictionless particles would normally undergo shear thinning. We can clarify that shear thickening in our simulations is due to enhanced energy dissipation via frictional interparticle forces. Moreover, we evidence the formation of dynamically correlated particle clusters of size ξ, which contribute to shear thickening via an increase in viscous dissipation. A scaling argument gives for the associated viscosity η(v)~ξ(2), which is in very good agreement with the data.

Rheology and shear band suppression in particle and chain mixtures

Using numerical simulations we consider an amorphous particle mixture which exhibits shear localization, and find that the addition of even a small fraction of chains strongly enhances the material strength, creating pronounced overshoot features in the stress-strain curves. The strengthening occurs in the case where the chains are initially perpendicular to the shear direction, leading to a suppression of the shear band. This also leads to stiffening effects that are typical of polymeric systems. For large strain, the chains migrate to the region where a shear band forms, resulting in a stress drop. For chains larger than the linear system size we find oscillatory behavior, which does not resemble polymeric systems since the second stress peak is larger than the first. Our results are also useful for providing insights into methods of controlling and strengthening granular materials against failure.

Jamming in systems with quenched disorder

We numerically study the effect of adding quenched disorder in the form of randomly placed pinning sites on jamming transitions in a disk packing that jams at a well-defined point J in the clean limit. Quenched disorder decreases the jamming density and introduces a depinning threshold. The onset of a finite threshold coincides with point J at the lowest pinning densities, but for higher pinning densities there is always a finite depinning threshold even well below jamming. We find that proximity to point J strongly affects the transport curves and noise fluctuations, and we observe a change from plastic behavior below jamming, where the system is highly heterogeneous, to elastic depinning above jamming. Many of the general features we find are related to other systems containing quenched disorder, including the peak effect observed in vortex systems.

Shear flow of non-Brownian suspensions close to jamming

The dynamical mechanisms controlling the rheology of dense suspensions close to jamming are investigated numerically, using simplified models for the relevant dissipative forces. We show that the velocity fluctuations control the dissipation rate and therefore the effective viscosity of the suspension. These fluctuations are similar in quasi-static simulations and for finite strain rate calculations with various damping schemes. We conclude that the statistical properties of grain trajectories-in particular the critical exponent of velocity fluctuations with respect to volume fraction φ-only weakly depend on the dissipation mechanism. Rather they are determined by steric effects, which are the main driving forces in the quasistatic simulations. The critical exponent of the suspension viscosity with respect to φ can then be deduced, and is consistent with experimental data.

Herschel-Bulkley shearing rheology near the athermal jamming transition

We consider the rheology of soft-core frictionless disks in two dimensions in the neighborhood of the athermal jamming transition. From numerical simulations of bidisperse, overdamped particles, we argue that the divergence of the viscosity below jamming is characteristic of the hard-core limit, independent of the particular soft-core interaction. We develop a mapping from soft-core to hard-core particles that recovers all the critical behavior found in earlier scaling analyses. Using this mapping we derive a relation that gives the exponent of the nonlinear Herschel-Bulkley rheology above jamming in terms of the exponent of the diverging viscosity below jamming.

Microscale rheology of a soft glassy material close to yielding

Using confocal microscopy, we study the flow of a model soft glassy material: a concentrated emulsion. We demonstrate the micro-macro link between in situ measured movements of droplets during the flow and the macroscopic rheological response of a concentrated emulsion, in the form of scaling relationships connecting the rheological "fluidity" with local standard deviation of the strain-rate tensor. Furthermore, we measure correlations between these local fluctuations, thereby extracting a correlation length which increases while approaching the yielding transition, in accordance with recent theoretical predictions.

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

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

Inhomogeneous shear flows in soft jammed materials with tunable attractive forces

We perform molecular dynamics simulations to characterize the occurrence of inhomogeneous shear flows in soft jammed materials. We use rough walls to impose a simple shear flow and study the athermal motion of jammed assemblies of soft particles in two spatial dimensions, both for purely repulsive interactions and in the presence of an additional short-range attraction of varying strength. In steady state, pronounced flow inhomogeneities emerge for all systems when the shear rate becomes small. Deviations from linear flow are stronger in magnitude and become very long lived when the strength of the attraction increases, but differ from permanent shear bands. Flow inhomogeneities occur in a stress window bounded by the dynamic and static yield stress values. Attractive forces enhance the flow heterogeneities because they accelerate stress relaxation, thus effectively moving the system closer to the yield stress regime where inhomogeneities are most pronounced. The present scenario for understanding the effect of particle adhesion on shear localization, which is based on detailed molecular dynamics simulations with realistic particle interactions, differs qualitatively from previous qualitative explanations and ad hoc theoretical modeling.

Contact percolation transition in athermal particulate systems

Typical quasistatic compression algorithms for generating jammed packings of purely repulsive, frictionless particles begin with dilute configurations and then apply successive compressions with the relaxation of the elastic energy allowed between each compression step. It is well known that during isotropic compression these systems undergo a first-order-like jamming transition at packing fraction φ(J) from an unjammed state with zero pressure and no force-bearing contacts to a jammed, rigid state with nonzero pressure, a percolating network of force-bearing contacts, and contact number z=2d, where d is the spatial dimension. Using computer simulations of two-dimensional systems with monodisperse and bidisperse particle size distributions, we investigate the second-order-like contact percolation transition, which precedes the jamming transition with φ(P)<φ(J) and signals the formation of a system-spanning cluster of non-force-bearing contacts between particles. By measuring the number of nonfloppy modes of the dynamical matrix, the displacement field between successive compression steps, and the overlap between the adjacency matrix, which represents the network of contacting grains, at φ and φ(J), we find that the contact percolation transition also signals the onset of a nontrivial mechanical response to applied stress. Our results show that cooperative particle motion occurs in unjammed systems significantly below the jamming transition for φ(P)<φ<φ(J), not only for jammed systems with φ>φ(J).

Similarities between protein folding and granular jamming

Grains and glasses, widely different materials, arrest their motions upon decreasing temperature and external load, respectively, in common ways, leading to a universal jamming phase diagram conjecture. However, unified theories are lacking, mainly because of the disparate nature of the particle interactions. Here we demonstrate that folded proteins exhibit signatures common to both glassiness and jamming by using temperature- and force-unfolding molecular dynamics simulations. Upon folding, proteins develop a peak in the interatomic force distributions that falls on a universal curve with experimentally measured forces on jammed grains and droplets. Dynamical signatures are found as a dramatic slowdown of stress relaxation upon folding. Together with granular similarities, folding is tied not just to the jamming transition, but a more nuanced picture of anisotropy, preparation protocol and internal interactions emerges. Results have implications for designing stable polymers and can open avenues to link protein folding to jamming theory.

Local anisotropy in globally isotropic granular packings

We report on two-dimensional computer simulations of frictionless granular packings at various area fractions φ above the jamming point φ(c). We measure the anisotropy in coarse-grained stress ε(s) and shear modulus ε(m) as functions of coarse-graining scale, R. ε(s) can be collapsed onto a master curve after rescaling R by a characteristic length scale ξ and ε(s) by an anisotropy magnitude A. Both A and ξ accelerate as φ→φ(c) from above, consistent with a divergence at φ(c). ε(m) shows no characteristic length scale and has a nontrivial power-law form, ε(m)~R(-0.62), over almost the entire range of R at all φ. These results suggest that the force chains present in the spatial structure of the quenched stress may be governed by different physics than the anomalous elastic response near jamming.

Jamming by shear

A broad class of disordered materials including foams, glassy molecular systems, colloids and granular materials can form jammed states. A jammed system can resist small stresses without deforming irreversibly, whereas unjammed systems flow under any applied stresses. The broad applicability of the Liu-Nagel jamming concept has attracted intensive theoretical and modelling interest but has prompted less experimental effort. In the Liu-Nagel framework, jammed states of athermal systems exist only above a certain critical density. Although numerical simulations for particles that do not experience friction broadly support this idea, the nature of the jamming transition for frictional grains is less clear. Here we show that jamming of frictional, disk-shaped grains can be induced by the application of shear stress at densities lower than the critical value, at which isotropic (shear-free) jamming occurs. These jammed states have a much richer phenomenology than the isotropic jammed states: for small applied shear stresses, the states are fragile, with a strong force network that percolates only in one direction. A minimum shear stress is needed to create robust, shear-jammed states with a strong force network percolating in all directions. The transitions from unjammed to fragile states and from fragile to shear-jammed states are controlled by the fraction of force-bearing grains. The fractions at which these transitions occur are statistically independent of the density. Jammed states with densities lower than the critical value have an anisotropic fabric (contact network). The minimum anisotropy of shear-jammed states vanishes as the density approaches the critical value from below, in a manner reminiscent of an order-disorder transition.

Microscopic theory of the jamming transition of harmonic spheres

We develop a microscopic theory to analyze the phase behavior and compute correlation functions of dense assemblies of soft repulsive particles both at finite temperature, as in colloidal materials, and at vanishing temperature, a situation relevant for granular materials and emulsions. We use a mean-field statistical mechanical approach which combines elements of liquid state theory to replica calculations to obtain quantitative predictions for the location of phase boundaries, macroscopic thermodynamic properties, and microstructure of the system. We focus, in particular, on the derivation of scaling properties emerging in the vicinity of the jamming transition occurring at large density and zero temperature. The new predictions we obtain for pair correlation functions near contact are tested using computer simulations. Our work also clarifies the conceptual nature of the jamming transition and its relation to the phenomenon of the glass transition observed in atomic liquids.

Jamming phase diagram for frictional particles

We investigate the jamming transition of frictional particulate systems via discrete element simulations, reporting the existence of new regimes, which are conveniently described in a jamming phase diagram with axes density, shear stress, and friction coefficient. The resulting jammed states are characterized by different mechanical and structural properties and appear not to be "fragile" as speculated. In particular, we find a regime, characterized by extremely long processes, with a diverging time scale, whereby a suspension first flows but then suddenly jams. We link this sudden jamming transition to the presence of impeded dilatancy.

Jamming in granular polymers

We examine the jamming transition in a two-dimensional granular polymer system using compressional simulations. The jamming density Φ(c) decreases with increasing length of the granular chain due to the formation of loop structures, in excellent agreement with recent experiments. The jamming density can be further reduced in mixtures of granular chains and granular rings, also as observed in experiment. We show that the nature of the jamming in granular polymer systems has pronounced differences from the jamming behavior observed for polydisperse two-dimensional disk systems at point J. This result provides further evidence that there is more than one type of jamming transition.

Critical scaling near jamming transition for frictional granular particles

The critical rheology of sheared frictional granular materials near jamming transition is numerically investigated. It is confirmed that there exists a true critical density which characterizes the onset of the yield stress and two fictitious critical densities which characterize the scaling laws of rheological properties. We find the existence of a hysteresis loop between two of the critical densities for each friction coefficient. It is noteworthy that the critical scaling law for frictionless jamming transition seems to be still valid even for frictional jamming despite using fictitious critical density values.

Critical scaling of shearing rheology at the jamming transition of soft-core frictionless disks

We perform numerical simulations to determine the shear stress and pressure of steady-state shear flow in a soft-disk model in two dimensions at zero temperature in the vicinity of the jamming transition ϕ{J}. We use critical point scaling analyses to determine the critical behavior at jamming, and we find that it is crucial to include corrections to scaling for a reliable analysis. We find that the relative size of these corrections are much smaller for pressure than for shear stress. We furthermore find a superlinear behavior for pressure and shear stress above ϕ{J}, both from the scaling analysis and from a direct analysis of pressure data extrapolated to the limit of vanishing shear rate.

Finite-size scaling at the jamming transition: corrections to scaling and the correlation-length critical exponent

We carry out a finite-size scaling analysis of the jamming transition in frictionless bidisperse soft core disks in two dimensions. We consider two different jamming protocols: (i) quench from random initial positions and (ii) quasistatic shearing. By considering the fraction of jammed states as a function of packing fraction for systems with different numbers of particles, we determine the spatial correlation length critical exponent ν ≈ 1 and show that corrections to scaling are crucial for analyzing the data. We show that earlier numerical results yielding ν < 1 are due to the improper neglect of these corrections.

Glassiness, rigidity, and jamming of frictionless soft core disks

The jamming of bidisperse soft core disks is considered, using a variety of different protocols to produce the jammed state. In agreement with other works, we find that cooling and compression can lead to a broad range of jamming packing fractions ϕ{J}, depending on cooling rate and initial configuration; the larger the degree of big particle clustering in the initial configuration, the larger will be the value of ϕ{J}. In contrast, we find that shearing disrupts particle clustering, leading to a much narrower range of ϕ{J} as the shear strain rate varies. In the limit of vanishingly small shear strain rate, we find a unique nontrivial value for the jamming density that is independent of the initial system configuration. We conclude that shear driven jamming is a unique and well-defined critical point in the space of shear driven steady states. We clarify the relation between glassy behavior, rigidity, and jamming in such systems and relate our results to recent experiments.

Length scales of dynamic heterogeneities in a network of fluctuating mechanical constraints

We describe a random bond network with fluctuating mechanical constraints for which structural and dynamical length scales are explicitly defined and can be compared with the observable properties of the system. The constraints give well-defined individual structures, describing rigid clusters of particles, and there are distinct structural changes between different configurations. In addition, unambiguous unconstrained motions, including extended collective motions, can be defined within each configuration of constraints. By considering fluctuations of bonds in the network, we can analyze the relaxation of a simple model amorphous material subject to randomly arranged mechanical constraints. The lengths of the collective motions are obtained directly and compared with both the relaxation kinetics and the susceptibility χ(4). Above percolation, the latter quantity is shown to decrease with increasing degree of constraint, a behavior linked, in this paper, to the presence of overconstraint in the network at bond densities beyond the percolation point. Our model system, despite a propensity for highly collective motions, shows that an increasing timescale for relaxation can be accompanied by a decreasing χ(4).

Fluidization of wet granulates under shear

Small amounts of a wetting liquid render sand a stiff and moldable material. The cohesive forces between the sand grains are caused by capillary bridges at the points of contact. Due to the finite strength of these bridges wet sand undergoes a transition from an arrested (i.e., solidified) to a fluidized state under an externally applied shear force. The transition between these two dynamic states is studied in a MD-type simulation of a two-dimensional assembly of bidisperse frictionless disks under the action of a cosine force profile. In addition to soft core repulsion the disks interact through a hysteretic and short ranged attractive force modeling the effect of the capillary bridges. In this model the transition between the fluidized and the arrested state is discontinuous and hysteretic. The parameter dependence of the critical force for solidification is modeled by combining theoretical arguments with a detailed numerical exploration of the transition. We address a range of densities from slightly below close packing until slightly above densities where the system approaches a shear-banded state. Differences and similarities of the transition in wet granulates to the jamming transition are also addressed.

Fluctuations, jamming, and yielding for a driven probe particle in disordered disk assemblies

Using numerical simulations we examine the velocity fluctuations and velocity-force curve characteristics of a probe particle driven with constant force through a two-dimensional disordered assembly of disks which has a well-defined jamming point J at a density of ϕJ=0.843. As ϕ increases toward ϕJ, the average velocity of the probe particle decreases and the velocity fluctuations show an increasingly intermittent or avalanchelike behavior. When ϕ is within a few percent of the jamming density, the velocity distributions are exponential, while when ϕ is less than 1% away from jamming, the velocity distributions have a power-law character with exponents in agreement with recent experiments. The velocity power spectra exhibit a crossover from a Lorentzian form to a 1/f shape near jamming. We extract a correlation length exponent ν which is in good agreement with recent shear simulations. For ϕ>ϕJ, there is a critical threshold force F(c) that must be applied for the probe particle to move through the sample which increases with increasing ϕ. The velocity-force curves are linear below jamming, while at jamming they have a power-law form. The onset of the probe motion above ϕJ occurs via a local yielding of the particles around the probe particle which we term a local shear banding effect.

Focused force transmission through an aqueous suspension of granules

We investigate force transmission through a layer of shear-thickening fluid, here a concentrated aqueous cornstarch suspension. When a solid body is pushed through this complex fluid and approaches its containing wall, a hardened volume of the suspension is observed that adds to the leading side of the body. This volume leads to an imprint on the wall which is made of molding clay. By studying the geometry of the hardened volume, inferred by the imprint shapes, we find that its geometry is determined by the size and speed of the body. By characterizing the response of the clay to deformation we show that the force transmitted through the suspension to the wall is localized. We also study other aspects of this dynamical hardening of the suspension, such as the effect of the substrate and body shape, and its relaxation as the imposed straining is stopped.

Mesoscopic length scale controls the rheology of dense suspensions

From the flow properties of dense granular suspensions on an inclined plane, we identify a mesoscopic length scale strongly increasing with volume fraction. When the flowing layer height is larger than this length scale, a diverging Newtonian viscosity is determined. However, when the flowing layer height drops below this scale, we evidence a nonlocal effective viscosity, decreasing as a power law of the flow height. We establish a scaling relation between this mesoscopic length scale and the suspension viscosity. These results support recent theoretical and numerical results implying collective and clustered granular motion when the jamming point is approached from below.

Asymmetric velocity correlations in shearing media

A model of soft frictionless disks in two dimensions at zero temperature is simulated with a shearing dynamics to study various kinds of asymmetries in sheared systems. We examine both single particle properties, the spatial velocity correlation function, and a correlation function designed to separate clockwise and counterclockwise rotational fields from one another. Among the rich and interesting behaviors we find that the velocity correlation along the two different diagonals corresponding to compression and dilation, respectively, are almost identical and, furthermore, that a feature in one of the correlation functions is directly related to irreversible plastic events.

Model for the scaling of stresses and fluctuations in flows near jamming

We probe flows of soft, viscous spheres near the jamming point, which acts as a critical point for static soft spheres. Starting from energy considerations, we find nontrivial scaling of velocity fluctuations with strain rate. Combining this scaling with insights from jamming, we arrive at an analytical model that predicts four distinct regimes of flow, each characterized by rational-valued scaling exponents. Both the number of regimes and the values of the exponents depart from prior results. We validate predictions of the model with simulations.

Unjamming dynamics: the micromechanics of a seismic fault model

The unjamming transition of granular systems is investigated in a seismic fault model via three dimensional molecular dynamics simulations. A two-time force-force correlation function, and a susceptibility related to the system response to pressure changes, allow us to characterize the stick-slip dynamics, consisting in large slips and microslips leading to creep motion. The correlation function unveils the micromechanical changes occurring both during microslips and slips. The susceptibility encodes the magnitude of the incoming microslip.

Jamming transitions in amorphous packings of frictionless spheres occur over a continuous range of volume fractions

We numerically produce fully amorphous assemblies of frictionless spheres in three dimensions and study the jamming transition these packings undergo at large volume fractions. We specify four protocols yielding a critical value for the jamming volume fraction which is sharply defined in the limit of large system size, but is different for each protocol. Thus, we directly establish the existence of a continuous range of volume fractions where nonequilibrium jamming transitions occur. However, these jamming transitions share the same critical behavior. Our results suggest that, even in the absence of partial crystalline ordering, a unique location of a random close packing does not exist.

Diffusion and velocity autocorrelation at the jamming transition

We perform numerical simulations to examine particle diffusion at steady shear in a soft-disk model in two dimensions and zero temperature around the jamming density. We find that the diffusion constant depends on shear rate as D approximately gamma below jamming and as D approximately gammaqD with qD<1 at the transition and set out to relate this to properties of the velocity autocorrelation function. It is found that this correlation function is governed by two processes with different time scales. The first time scale, the inverse of the externally applied shear rate, controls an exponential decay of the correlations whereas the second time scale, equal to the inverse shear stress, governs an algebraic decay with time. The obtained value of qD is related to these properties of the correlation function.

Scaling of the glassy dynamics of soft repulsive particles: a mode-coupling approach

We combine the hypernetted chain approximation of liquid state theory with the mode-coupling theory of the glass transition to analyze the structure and dynamics of soft spheres interacting via harmonic repulsion. We determine the locus of the fluid-glass dynamic transition in a temperature--volume fraction phase diagram. The zero-temperature (hard-sphere) glass transition influences the dynamics at finite temperatures in its vicinity. This directly implies a form of dynamic scaling for both the average relaxation time and dynamic susceptibilities quantifying dynamic heterogeneity. We discuss several qualitative disagreements between theory and existing simulations at equilibrium. Our theoretical results are, however, very similar to numerical results for the driven athermal dynamics of repulsive spheres, suggesting that "mean-field" mode-coupling approaches might be good starting points to describe these nonequilibrium dynamics.

Jamming at zero temperature, zero friction, and finite applied shear stress

Shear stress, as temperature, induces particle motion, and affects the jamming transition from a fluid to a disordered solid state. Here, we show that at finite shear stress, the jamming transition is characterized by the presence of hysteresis, as in a given range of control parameters a flowing or a jammed state can be found, depending on whether the system is prepared coming from the fluid or the jammed phase. At small shear stress, where the hysteresis is negligible, the jamming transition has a mixed first-order second-order character close to that found at the glass transition of thermal systems, with discontinuities in the asymptotic values of two time quantities such as the self-intermediate scattering function.

Soft Dynamics simulation. 2. Elastic spheres undergoing a T(1) process in a viscous fluid

Robust empirical constitutive laws for granular materials in air or in a viscous fluid have been expressed in terms of timescales based on the dynamics of a single particle. However, some behaviours such as viscosity bifurcation or shear localization, observed also in foams, emulsions, and block copolymer cubic phases, seem to involve other micro-timescales which may be related to the dynamics of local particle reorganizations. In the present work, we consider a T(1) process as an example of a rearrangement. Using the Soft Dynamics simulation method introduced in the first paper of this series, we describe theoretically and numerically the motion of four elastic spheres in a viscous fluid. Hydrodynamic interactions are described at the level of lubrication (Poiseuille squeezing and Couette shear flow) and the elastic deflection of the particle surface is modeled as Hertzian. The duration of the simulated T(1) process can vary substantially as a consequence of minute changes in the initial separations, consistently with predictions. For the first time, a collective behaviour is thus found to depend on a parameter other than the typical volume fraction of particles.

Random organization and plastic depinning

We provide evidence that the general phenomenon of plastic depinning can be described as an absorbing phase transition, and shows the same features as the random organization which was recently studied in periodically driven particle systems [L. Corte, Nature Phys. 4, 420 (2008)]. In the plastic flow system, the pinned regime corresponds to the absorbing state and the moving state corresponds to the fluctuating state. When an external force is suddenly applied, the system eventually organizes into one of these two states with a time scale that diverges as a power law at a nonequilibrium transition. We propose a simple experiment to test for this transition in systems with random disorder.