Avalanche contribution to shear modulus of granular materials

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.

An exact expression of three-body system for the complex shear modulus of frictional granular materials

We propose a simple model comprising three particles to study the nonlinear mechanical response of jammed frictional granular materials under oscillatory shear. Owing to the introduction of the simple model, we obtain an exact analytical expression of the complex shear modulus for a system including many monodispersed disks, which satisfies a scaling law in the vicinity of the jamming point. These expressions perfectly reproduce the shear modulus of the many-body system with low strain amplitudes and friction coefficients. Even for disordered many-body systems, the model reproduces the results by introducing a single fitting parameter.

Softening and Residual Loss Modulus of Jammed Grains under Oscillatory Shear in an Absorbing State

From a theoretical study of the mechanical response of jammed materials comprising frictionless and overdamped particles under oscillatory shear, we find that the material becomes soft, and the loss modulus remains nonzero even in an absorbing state where any irreversible plastic deformation does not exist. The trajectories of the particles in this region exhibit hysteresis loops. We succeed in clarifying the origin of the softening of the material and the residual loss modulus with the aid of Fourier analysis. We also clarify the roles of the yielding point in the softening to distinguish the plastic deformation from reversible deformation in the absorbing state.

Geometrical properties of mechanically annealed systems near the jamming transition

Geometrical properties of two-dimensional mixtures near the jamming transition point are numerically investigated using harmonic particles under mechanical training. The configurations generated by the quasi-static compression and oscillatory shear deformations exhibit anomalous suppression of the density fluctuations, known as hyperuniformity, below and above the jamming transition. For the jammed system trained by compression above the transition point, the hyperuniformity exponent increases. For the system below the transition point under oscillatory shear, the hyperuniformity exponent also increases until the shear amplitude reaches the threshold value. The threshold value matches with the transition point from the point-reversible phase where the particles experience no collision to the loop-reversible phase where the particles' displacements are non-affine during a shear cycle before coming back to an original position. The results demonstrated in this paper are explained in terms of neither of universal criticality of the jamming transition nor the nonequilibrium phase transitions.

Unified view of avalanche criticality in sheared glasses

Plastic events in sheared glasses are considered an example of so-called avalanches, whose sizes obey a power-law probability distribution with the avalanche critical exponent τ. Although the so-called mean-field depinning (MFD) theory predicts a universal value of this exponent, τ_{MFD}=1.5, such a simplification is now known to connote qualitative disagreement with realistic systems. Numerically and experimentally, different values of τ have been reported depending on the literature. Moreover, in the elastic regime, it has been noted that the critical exponent can be different from that in the steady state, and even criticality itself is a matter of debate. Because these confusingly varying results have been reported under different setups, our knowledge of avalanche criticality in sheared glasses is greatly limited. To gain a unified understanding, in this work, we conduct a comprehensive numerical investigation of avalanches in Lennard-Jones glasses under athermal quasistatic shear. In particular, by excluding the ambiguity and arbitrariness that has crept into the conventional measurement schemes, we achieve high-precision measurement and demonstrate that the exponent τ in the steady state follows the prediction of MFD theory, τ_{MFD}=1.5. Our results also suggest that there are two qualitatively different avalanche events. This binariness leads to the nonuniversal behavior of the avalanche size distribution and is likely to be the cause of the varying values of τ reported thus far. To investigate the dependence of criticality and universality on applied shear, we further study the statistics of avalanches in the elastic regime and the ensemble of the first avalanche event in different samples, which provide information about the unperturbed system. We show that while the unperturbed system is indeed off-critical, criticality gradually develops as shear is applied. The degree of criticality is encoded in the fractal dimension of the avalanches, which starts from zero in the off-critical unperturbed state and saturates in the steady state. Moreover, the critical exponent τ is consistent with the prediction of the MFD τ_{MFD} universally, regardless of the amount of applied shear, once the system becomes critical.

Shear modulus and reversible particle trajectories of frictional granular materials under oscillatory shear

In this study, we numerically investigated the mechanical responses and trajectories of frictional granular particles under oscillatory shear in the reversible phase where particle trajectories form closed loops below the yielding point. When the friction coefficient is small, the storage modulus exhibits softening, and the loss modulus remains finite in the quasi-static limit. As the friction coefficient increases, the softening and residual loss modulus are suppressed. The storage and loss moduli satisfy scaling laws if they are plotted as functions of the areas of the loop trajectories divided by the strain amplitude and diameter of grains, at least for small values of the areas.

Adhesion as a trigger of droplet polarization in flowing emulsions

Tissues are subjected to large external forces and undergo global deformations during morphogenesis. We use synthetic analogues of tissues to study the impact of cell-cell adhesion on the response of cohesive cellular assemblies under such stresses. In particular, we use biomimetic emulsions in which the droplets are functionalized in order to exhibit specific droplet-droplet adhesion. We flow these emulsions in microfluidic constrictions and study their response to this forced deformation via confocal microscopy. We find that the distributions of avalanche sizes are conserved between repulsive and adhesive droplets. However, adhesion locally impairs the rupture of droplet-droplet contacts, which in turn pulls on the rearranging droplets. As a result, adhesive droplets are a lot more deformed along the axis of elongation in the constriction. This finding could shed light on the origin of polarization processes during morphogenesis.

Scaling laws for frictional granular materials confined by constant pressure under oscillatory shear

Herein we numerically study the rheology of a two-dimensional frictional granular system confined by constant pressure under oscillatory shear. Several scaling laws for the storage and loss moduli against the scaled strain amplitude have been found. The scaling laws in plastic regime for large strain amplitude can be understood by the angular distributions of the contact force. The scaling exponents are estimated by considering the physical mechanism.

Shear jamming, discontinuous shear thickening, and fragile states in dry granular materials under oscillatory shear

We numerically study the linear response of two-dimensional frictional granular materials under oscillatory shear. The storage modulus G^{'} and the loss modulus G^{''} in the zero strain rate limit depend on the initial strain amplitude of the oscillatory shear before measurement. The shear jammed state (satisfying G^{'}>0) can be observed at an amplitude greater than a critical initial strain amplitude. The fragile state is defined by the emergence of liquid-like and solid-like states depending on the form of the initial shear. In this state, the observed G^{'} after the reduction of the strain amplitude depends on the phase of the external shear strain. The loss modulus G^{''} exhibits a discontinuous jump corresponding to discontinuous shear thickening in the fragile state.

Classification of the reversible-irreversible transitions in particle trajectories across the jamming transition point

The reversible-irreversible (RI) transition of particle trajectories in athermal colloidal suspensions under cyclic shear deformation is an archetypal nonequilibrium phase transition which has attracted much attention recently. Most studies of the RI transitions have focused on either dilute limit or very high densities well above the jamming transition point. The transition between the two limiting cases is largely unexplored. In this paper, we study the RI transition of athermal frictionless colloidal particles over a wide range of densities, with emphasis on the region below φ_J, by using oscillatory sheared molecular dynamics simulation. We reveal that the nature of the RI transitions in the intermediate densities is very rich. As demonstrated by the previous work by Schreck et al. [Phys. Rev. E: Stat., Nonlinear, Soft Matter Phys., 2013, 88, 052205], there exist the point-reversible and the loop-reversible phases depending on the density and the shear strain amplitude. We find that, between the two reversible phases, a quasi-irreversible phase where the particles' trajectories are highly non-affine and diffusive. The averaged number of contacts of particles is found to characterize the phase boundaries. We also find that the system undergoes the yielding transition below but in the vicinity of φ_J when the strain with a small but finite strain rate is applied. This yielding transition line matches with the RI transition line separating the loop-reversible from the irreversible phases. Surprisingly, the nonlinear rheological response called "softening" has been observed even below φ_J. These findings imply that geometrical properties encoded in the sheared configurations control the dynamical transitions.

Gardner physics in amorphous solids and beyond

One of the most remarkable predictions to emerge out of the exact infinite-dimensional solution of the glass problem is the Gardner transition. Although this transition was first theoretically proposed a generation ago for certain mean-field spin glass models, its materials relevance was only realized when a systematic effort to relate glass formation and jamming was undertaken. A number of nontrivial physical signatures associated with the Gardner transition have since been considered in various areas, from models of structural glasses to constraint satisfaction problems. This perspective surveys these recent advances and discusses the novel research opportunities that arise from them.

Drastic slowdown of the Rayleigh-like wave in unjammed granular suspensions

We present an experimental investigation of Rayleigh-like wave propagation along the surface of a dense granular suspension. Using an ultrafast ultrasound scanner, we monitor the softening of the shear modulus via the Rayleigh-like wave velocity slowdown in the optically opaque medium as the driving amplitude increases. For such nonlinear behavior two regimes are found when increasingthe driving amplitude progressively: First, we observe a significant shear modulus weakening due to the microslip on the contact level without macroscopic rearrangements of grains. Second, there is a clear macroscopic plastic rearrangement accompanied by a modulus decrease up to 88%. A friction model is proposed to describe the interplay between nonlinear elasticity and plasticity, which highlights the crucial effect of contact slipping before contact breaking or loss. Investigation of this nonlinear Rayleigh-like wave may bridge the gap between two disjoint approaches for describing the dynamics near unjamming: linear elastic soft modes and nonlinear collisional shock.

Seismicity in sheared granular matter

We report on experiments investigating the dynamics of a slider that is pulled by a spring across a granular medium consisting of a vertical layer of photoelastic disks. The motion proceeds through a sequence of discrete events, analogous to seismic shocks, in which elastic energy stored in the spring is rapidly released. We measure the statistics of several properties of the individual events: the energy loss in the spring, the duration of the movement, and the temporal profile of the slider motion. We also study certain conditional probabilities and the statistics of mainshock-aftershock sequences. At low driving rates, we observe crackling with Omori-Utsu, Båth, and waiting time laws similar to those observed in seismic dynamics. At higher driving rates, where the sequence of events shows strong periodicity, we observe scaling laws and asymmetrical event shapes that are clearly distinguishable from those in the crackling regime.

Coarsening and mechanics in the bubble model for wet foams

Aqueous foams are an important model system that displays coarsening dynamics. Coarsening in dispersions and foams is well understood in the dilute and dry limits, where the gas fraction tends to zero and one, respectively. However, foams are known to undergo a jamming transition from a fluidlike to a solidlike state at an intermediate gas fraction ϕ_{c}. Much less is known about coarsening dynamics in wet foams near jamming, and the link to mechanical response, if any, remains poorly understood. Here we probe coarsening and mechanical response using numerical simulations of a variant of the Durian bubble model for wet foams. As in other coarsening systems we find a steady state scaling regime with an associated particle size distribution. We relate the time rate of evolution of the coarsening process to the wetness of the foam and identify a characteristic coarsening time that diverges approaching jamming. We further probe mechanical response of the system to strain while undergoing coarsening. There are two competing timescales, namely the coarsening time and the mechanical relaxation time. We relate these to the evolution of the elastic response and the mechanical structure.

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.

Local and global avalanches in a two-dimensional sheared granular medium

We present the experimental and numerical studies of a two-dimensional sheared amorphous material composed of bidisperse photoelastic disks. We analyze the statistics of avalanches during shear including the local and global fluctuations in energy and changes in particle positions and orientations. We find scale-free distributions for these global and local avalanches denoted by power laws whose cutoffs vary with interparticle friction and packing fraction. Different exponents are found for these power laws depending on the quantity from which variations are extracted. An asymmetry in time of the avalanche shapes is evidenced along with the fact that avalanches are mainly triggered by the shear bands. A simple relation independent of the intensity is found between the number of local avalanches and the global avalanches they form. We also compare these experimental and numerical results for both local and global fluctuations to predictions from mean-field and depinning theories.

Stress relaxation in viscous soft spheres

We report the results of molecular dynamics simulations of stress relaxation tests in athermal viscous soft sphere packings close to their unjamming transition. By systematically and simultaneously varying both the amplitude of the applied strain step and the pressure of the initial condition, we access both linear and nonlinear response regimes and control the distance to jamming. Stress relaxation in viscoelastic solids is characterized by a relaxation time τ* that separates short time scales, where viscous loss is substantial, from long time scales, where elastic storage dominates and the response is essentially quasistatic. We identify two distinct plateaus in the strain dependence of the relaxation time, one each in the linear and nonlinear regimes. The height of both plateaus scales as an inverse power law with the distance to jamming. By probing the time evolution of particle velocities during relaxation, we further identify a correlation between mechanical relaxation in the bulk and the degree of non-affinity in the particle velocities on the micro scale.

Discontinuous change of shear modulus for frictional jammed granular materials

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

Exploring the complex free-energy landscape of the simplest glass by rheology

For amorphous solids, it has been intensely debated whether the traditional view on solids, in terms of the ground state and harmonic low energy excitations on top of it, such as phonons, is still valid. Recent theoretical developments of amorphous solids revealed the possibility of unexpectedly complex free-energy landscapes where the simple harmonic picture breaks down. Here we demonstrate that standard rheological techniques can be used as powerful tools to examine nontrivial consequences of such complex free-energy landscapes. By extensive numerical simulations on a hard sphere glass under quasistatic shear at finite temperatures, we show that above the so-called Gardner transition density, the elasticity breaks down, the stress relaxation exhibits slow, and ageing dynamics and the apparent shear modulus becomes protocol-dependent. Being designed to be reproducible in laboratories, our approach may trigger explorations of the complex free-energy landscapes of a large variety of amorphous materials.

Contact changes of sheared systems: Scaling, correlations, and mechanisms

We probe the onset and effect of contact changes in two-dimensional soft harmonic particle packings which are sheared quasistatically under controlled strain. First, we show that, in the majority of cases, the first contact changes correspond to the creation or breaking of contacts on a single particle, with contact breaking overwhelmingly likely for low pressures and/or small systems, and contact making and breaking equally likely for large pressures and in the thermodynamic limit. The statistics of the corresponding strains are near-Poissonian, in particular for large-enough systems. The mean characteristic strains exhibit scaling with the number of particles N and pressure P and reveal the existence of finite-size effects akin to those seen for linear response quantities [C. P. Goodrich et al., Phys. Rev. Lett. 109, 095704 (2012)PRLTAO0031-900710.1103/PhysRevLett.109.095704; C. P. Goodrich et al., Phys. Rev. E 90, 022138 (2014)].PLEEE81539-375510.1103/PhysRevE.90.022138 Second, we show that linear response accurately predicts the strains of the first contact changes, which allows us to accurately study the scaling of the characteristic strains of making and breaking contacts separately. Both of these show finite-size scaling, and we formulate scaling arguments that are consistent with the observed behavior. Third, we probe the effect of the first contact change on the shear modulus G and show in detail how the variation of G remains smooth and bounded in the large-system-size limit: Even though contact changes occur then at vanishingly small strains, their cumulative effect, even at a fixed value of the strain, are limited, so, effectively, linear response remains well defined. Fourth, we explore multiple contact changes under shear and find strong and surprising correlations between alternating making and breaking events. Fifth, we show that by making a link with extremal statistics, our data are consistent with a very slow crossover to self-averaging with system size, so the thermodynamic limit is reached much more slowly than expected based on finite-size scaling of elastic quantities or contact breaking strains.

Experimental Evidence of the Gardner Phase in a Granular Glass

Analyzing the dynamics of a vibrated bidimensional packing of bidisperse granular disks below jamming, we provide evidence of a Gardner phase deep into the glass phase. To do so, we perform several compression cycles within a given realization of the same glass and show that the particles select different average vibrational positions at each cycle, while the neighborhood structure remains unchanged. The separation between the cages obtained for different compression cycles plateaus with an increasing packing fraction, while the mean square displacement steadily decreases. This phenomenology is strikingly similar to that reported in recent numerical observations when entering the Gardner phase, for a mean-field model of glass as well as for hard spheres in finite dimension. We also characterize the distribution of the cage order parameters. Here we note several differences from the numerical results, which could be attributed to activated processes and cage heterogeneities.

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

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

Elastic moduli and vibrational modes in jammed particulate packings

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

Spontaneous instabilities and stick-slip motion in a generalized Hébraud-Lequeux model

We revisit the Hébraud-Lequeux (HL) model for the rheology of jammed materials and argue that a possibly important time scale is missing from HL's initial specification. We show that our generalization of the HL model undergoes interesting oscillating instabilities for a wide range of parameters, which lead to intermittent, stick-slip flows under constant shear rate. The instability we find is akin to the synchronization transition of coupled elements that arises in many different contexts (neurons, fireflies, financial bankruptcies, etc.). We hope that our scenario could shed light on the commonly observed intermittent, serrated flows of glassy materials under shear.

Divergence of Viscosity in Jammed Granular Materials: A Theoretical Approach

A theory for jammed granular materials is developed with the aid of a nonequilibrium steady-state distribution function. The approximate nonequilibrium steady-state distribution function is explicitly given in the weak dissipation regime by means of the relaxation time. The theory quantitatively agrees with the results of the molecular dynamics simulation on the critical behavior of the viscosity below the jamming point without introducing any fitting parameter.

Following the evolution of hard sphere glasses in infinite dimensions under external perturbations: compression and shear strain

We consider the adiabatic evolution of glassy states under external perturbations. The formalism we use is very general. Here we use it for infinite-dimensional hard spheres where an exact analysis is possible. We consider perturbations of the boundary, i.e., compression or (volume preserving) shear strain, and we compute the response of glassy states to such perturbations: pressure and shear stress. We find that both quantities overshoot before the glass state becomes unstable at a spinodal point where it melts into a liquid (or yields). We also estimate the yield stress of the glass. Finally, we study the stability of the glass basins towards breaking into sub-basins, corresponding to a Gardner transition. We find that close to the dynamical transition, glasses undergo a Gardner transition after an infinitesimal perturbation.