Abstract
In a recent paper (Zhao et al., Phys Rev X, 2022, 12: 031,021), we reported experimental observations of “ultrastable” states in a shear-jammed granular system subjected to small-amplitude cyclic shear. In such states, all the particle positions and contact forces are reproduced after each shear cycle so that a strobed image of the stresses and particle positions appears static. In the present work, we report further analyses of data from those experiments to characterize both global and local responses of ultrastable states within a shear cycle, not just the strobed dynamics. We find that ultrastable states follow a power-law relation between shear modulus and pressure with an exponent β ≈ 0.5, reminiscent of critical scaling laws near jamming. We also examine the evolution of contact forces measured using photoelasticimetry. We find that there are two types of contacts: non-persistent contacts that reversibly open and close; and persistent contacts that never open and display no measurable sliding. We show that the non-persistent contacts make a non-negligible contribution to the emergent shear modulus. We also analyze the spatial correlations of the stress tensor and compare them to the predictions of a recent theory of the emergent elasticity of granular solids, the Vector Charge Theory of Granular mechanics and dynamics (VCTG) (Nampoothiri et al., Phys Rev Lett, 2020, 125: 118,002). We show that our experimental results can be fit well by VCTG, assuming uniaxial symmetry of the contact networks. The fits reveal that the response of the ultrastable states to additional applied stress is substantially more isotropic than that of the original shear-jammed states. Our results provide important insight into the mechanical properties of frictional granular solids created by shear.
1 Introduction
Granular materials are athermal collections of particles that interact with each other only when they form direct, frictional contacts. These materials can jam into solid packings that statically resist applied stresses [–]. Shear-induced jamming occurs in a variety of disordered, complex systems, including granular suspensions [–] and dry granular materials with [–] or without friction [–]. The stability of shear-jammed states, however, remains at best partially understood. In a recent experiment [], the stability of shear-jammed states in a frictional granular system was systematically examined by monitoring their evolution under small-amplitude, volume-conserving cyclic shear. Many shear-jammed packings relaxed to a stress-free, diffusive steady state under cyclic strain amplitude as small as 1%. However, in some cases, the shear-jammed system relaxed into an unexpected state in which all microscopic degrees of freedom, including particle positions, orientations, and contact forces, remain the same for thousands of shear cycles. These states were termed “ultrastable” to distinguish them from originally formed shear-jammed states that would deform plastically under a single shear cycle with same strain amplitude. They emerged in athermal, frictional granular packings and are thus qualitatively different from the states of glasses obtained using vapor deposition that have also been termed ultrastable []. Nevertheless, the two systems are similar in that the ultrastable shear-jammed granular packings have smaller pressure and behave more like an elastic ordinary solid than other shear-jammed packings, and ultrastable glasses have lower energy and are more stable against shearing than ordinary glasses [, ].
In our earlier work [], we found that a reversibility transition and a jamming/unjamming transition coincide at the phase boundary between the two types of nonequilibrium steady states induced by cyclic shearing: the ultrastable states that return to the same microscopic configuration after each cycle and the fluid-like unjammed states in which particles undergo diffusive displacements. Without changing the volume fraction, the different types of steady states can be realized by changing either the shear strain γI used to form an original shear-jammed state or the cyclic strain amplitude δγ. A stability diagram is given in Ref. []. Notably, ultrastable states formed by larger γI survive under larger δγ. The transition from ultrastable states to unjammed states with increasing δγ or decreasing γI may be viewed as a yielding transition. This transition is similar to the oscillatory yielding of amorphous solids [–], which is accompanied by a microscopic reversibility transition that can be classified as an absorbing-state transition [, , ]. The ultrastable states reported in Ref. [] are both reversible and mechanically stable and are thus similar to the absorbing states of amorphous solids under cyclic shear with a strain amplitude below the threshold for oscillatory yielding [], but different from other absorbing states that do not have a mechanically stable structure as in the case of dilute suspensions []. The global stress-strain curves for the ultrastable states reported in Ref. [] appear to be highly elastic, but the internal deformation occurring within individual shear cycles was not examined for evidence of reversible plastic events [] or loops in particle trajectories [–].
In addition to characterizing the grain scale deformations occurring in the ultrastable state, we study the relation between the global elasticity features and the local stresses at the grain scale. A recent theory, termed the Vector Charge Theory of Granular mechanics and dynamics (VCTG) [, ], suggests a promising approach for relating the global elastic behaviour to features of the contact forces between individual particles. VCTG is a stress-only framework for amorphous solids; it does not rely on a unique reference structure to define strain. This theory maps the mechanical response of granular solids to the static, dielectric response of a tensorial electromagnetism with the electric polarizability of the medium mapping to emergent elastic moduli. VCTG relates the spatial correlations of the stress-tensor to these elastic moduli, which emerge from the underlying contact and force network. While previous experiments confirmed some features of the stress correlation functions predicted by the theory [], there has not been a direct comparison of the elastic constants obtained from fitting the stress correlations to the elastic constants measured from stress-strain curves in experiments. It is thus of interest to examine our data in the framework of such a theory.
In the present work, we report a detailed analysis of the elastic properties of the previously reported ultrastable shear-jammed states []. We find that the emergent shear modulus follows a power-law relation with pressure with an exponent consistent with some numerical models. However, the shear response contains a special non-linear feature: there are many contacts that reversibly open and close under low amplitude cyclic shear. These non-persistent contacts contribute a non-negligible portion to the global effective shear modulus. We also examine the relation between the global elastic constants and internal stress correlations predicted by VCTG. The analysis leads to intriguing scalings of the emergent elastic properties of the system and uncovers a feature that reflects how cyclic shear modifies the elastic properties of a jammed packing. Our results bring new insights to the elasticity of frictional granular materials near jamming.
2 Materials and methods
The analyses in the present paper are performed on the same set of experiments reported in Ref. []. In this work, we focus on the evolution of the system within several shear cycles after an ultrastable state is formed while Ref. [] focused on the strobed states. The materials and experimental protocols are briefly summarized here. More details can be found in Ref. [].
Our model granular system consists of a bidisperse layer of photoelastic disc with same height, 6.8 mm, but different diameters: db = 15.9 mm and ds = 12.7 mm. The static friction coefficient between the particles is μs = 0.87 ± 0.03. Under static diametric loading, the normal contact force law is roughly Hertzian. For a small disc squeezed between to rigid surfaces, we measurewhere fn is the normal contact force, δ/ds is the diametric strain, rs is the radius of the small disc, and ϵs = 2.73 N⋅m. (This expression slightly overestimates the weak forces. Details on contact force law calibration are given in Appendix A of Ref. [].) The discs are placed in a simple shear box with a parallelogram boundary and a multi-slat base that promotes homogeneous shear when the angle between the boundaries changes. A schematic of the shear apparatus is shown in Figure 1A, and more details can be found in Refs. [, ]. The number of particles is fixed at 1040 for all experiments. The area of the shear box is also kept constant throughout. The packing fraction, defined as the total area of particles divided by the area of the shear box, is ϕ = 0.816 for all experiments, which is below the frictionless isotropic jamming point ϕJ ≈ 0.835 estimated using the same apparatus [].
FIGURE 1
At the beginning of each experiment, the particles are randomly placed in a zero-stress, unjammed configuration. When the boundary walls impose a volume-conserving simple shear deformation, the parallel bottom slats move accordingly to impose a uniform internal shear strain field. Static friction causes the particles to move with the slats in an unjammed configuration. Such a substrate-assisted shear protocol avoids boundary-induced density heterogeneity and leads to homogeneous shear-jammed states []. The frictional forces between particles and the slats are, however, much smaller in magnitude than the mean contact force in the jammed states.
Starting from the unjammed initial state, each shear experiment consists of two stages: (i) an initial shear that forms a shear-jammed state; and (ii) a number of consecutive shear cycles that cause the shear-jammed state to transform. The two stages are sketched in Figure 1B. In the initial shear stage, we apply a shear strain γI that transforms an initial chosen parallelogram into a rectangle. Note that the jammed states formed by different γI have different features of the contact and force networks [, , ]. In the cyclic shear stage, we apply a series of N small-amplitude shear cycles with strain amplitude δγ ≪ γI. The value of N is at least 1500, and the largest one we used is 4800. In this stage, we only monitor states before and after each complete shear cycle, as sketched by the small circles in Figure 1B. In the last two or three cycles, where an ultrastable state has been formed, we also record images of the system within the shear cycles. The present work focuses on these cycles, as highlighted by the light blue region in Figure 1B.
At all stages, the imposed shear can be considered quasistatic; the system reaches mechanical equilibrium much faster than the overall shearing rate []. A high-resolution camera is used to take images of the system, and physical quantities of interest are measured using image processing techniques. For each state of interest, we measure all the vector contact forces between individual particles using a nonlinear fitting algorithm [], details of which can be found in Appendix C of Ref. []. The stress tensor is calculated from the measured contact forces as in Refs. [, , ].where fij is the vector force applied to partilce i particle by the particle j, rij is the displacement of the contact point of particles i and j from the center of particle i, and the summation runs over particle indices i and j from 1 to the number of particles Np. We exclude the particles that are in direct contact with the boundary walls, and S in Eq. 2 is the total Voronoi area for the particles that do not belong to the boundary layer.
Ultrastable states are formed when γI is large and δγ is small. For smaller γI or larger δγ, the system relaxes to an unjammed, fluid-like state. Figure 1 shows snapshots illustrating the different behaviors. Figures 1C–F shows example images obtained from typical runs with different γI. All of the images are taken through a polariscope so that only the discs supporting finite stress are visible. The second row (labeled n = 0) shows the stress state after the initial shear γI, and the third and fourth rows show the states after 1000 and 2000 shear cycles with δγ = 0.95%. The nearly blank images in columns (C) and (D) indicate that the system has relaxed to a steady state with nearly zero pressure. The close similarity between images in the third and fourth rows of columns (E) and (F) indicate that ultrastable states are reached within 1000 cycles. Our focus in this paper is on the elastic properties of these ultrastable states.
3 Results
3.1 Emergent shear modulus
We first examine the global shear modulus G for the ultrastable states under cyclic shear. The insert panel of Figure 2A plots an example σxy evolution in a shear cycle. The filled red circle is the ultrastable state being considered. The behavior appears similar to a viscoelastic material in a highly elastic regime. The finite area enclosed by the curve suggests that there is measurable energy dissipation inside the system, although this hysteresis in the stress-strain curve is much smaller than that of shear-jammed states formed by initial shear alone []. The microscopic mechanisms responsible for this small dissipation could be the sliding of particles over the base and the confining walls or the sliding at inter-particle contacts. We have examined the distribution of the tangential to normal contact force ratio and find that most force-carrying contacts are far from the Coulomb threshold for sliding. However, we cannot exclude the existence of reversible sliding at weak contacts. In addition, the viscoelasticity of the polyurethane photoelastic discs leads to small but measurable hysteresis in the force-displacement curve for a single particle under cyclic diametric loading, as shown in Appendix A in Ref. []. This material effect may also contribute to the global hysteresis in the stress-strain curve.
FIGURE 2
From numerical simulations, it is known that both the elastic constants and the stresses of jammed granular materials follow scaling laws in the vicinity of the jamming point [–]. While the exponents associated with stresses and contact numbers have been examined in experiments [], the scaling of elastic moduli remains largely unexplored, especially for frictional systems. Previous experiments measured the scaling indirectly through acoustic propagation []. Our experimental system allows us to study the scaling of the shear modulus of the ultrastable states.
We define the shear modulus G as the slope of the curve in the vicinity of the ultrastable state. In practice, we fit a straight line to the rising branch of the curve and obtain the fitted slope. The jamming point has zero pressure. Thus, the pressure p serves as the measure of distance to the jamming point. While the excess contact number may be a more fundamental quantity, pressure is measured with higher accuracy than contact number in photoelastic experiments. Thus, in this paper, we focus on the relation between G and p. We also note that this relation is of great engineering interest [], as in real applications it is easier to control the pressure in a packing than the average contact number. Figure 2A shows the measured values of G plotted as a function of pressure for various small values of δγ. Without any rescaling, the data falls on a single curve, suggesting that all the ultrastable states are governed by a universal scaling relation despite the rather special protocol used to generate them. The independence of δγ also suggests that the ultrastable states are below the onset of softening []. The solid curve shows a power-law fit of the formwith the fit parameters G0 = (95 ± 15)N/m and β = 0.50 ± 0.06. A log-log plot of same data is shown in Figure 2B. Interestingly, the shear modulus measured using sound propagation exhibits similar dependence on pressure []. Numerical simulations with frictionless spherical particles near jamming interacting through linear spring force laws also show β = 0.5 for both isotropic [, , ] and shear-jammed systems []. For frictionless Hertzian contact models simulations give β = 2/3 over the range of dimensionless pressures studied in our system [, ]. Simulations of 2D packings of frictional spheres with Hertz-Mindlin forces and with friction coefficients similar to ours show β between 1/2 and 2/3 []. As mentioned above, the contact force law between our discs is roughly Hertzian, and our measured β appears similar to those obtained in numerical simulations. In addition, the range of dimensionless pressures in our experiments falls in the range studied in simulations, suggesting that our system is indeed close enough to the jamming point to be in the scaling regime. We note that a recent simulation of 2D frictional particles under oscillatory shear with linear-dashpot contact model also found β = 1/2 for the small-strain plateau shear modulus []. Another recent simulation shows that β can be different for deformable particles whose shape is controlled by surface tension rather than internal bulk stresses [].
3.2 Persistent and non-persistent contacts
The internal deformation of the system exhibits non-trivial features. Particle displacements are non-affine, and many force-bearing contacts are activated and deactivated reversibly during a shear cycle. These contacts contribute to the emergent elastic moduli through a nonlinear process that can not be predicted by analyzing a single contact network.
We here report experimental characterizations of the two types of contacts that contribute to the emergent elastic modulus of the packing in a shear cycle. The non-persistent contacts are those that break reversibly during a shear cycle, while the persistent contacts never break once the ultrastable state is reached.
3.2.1 Non-persistent contacts
We first demonstrate the existence of non-persistent contacts in an example ultrastable state formed by initial shear γI = 14.7% and cyclic shear amplitude δγ = 0.95%. After the system has settled in the ultrastable state for thousands of cycles, we examine the response of the system during the next three shear cycles. Figure 3A plots the global shear strain for these three shear cycles.
FIGURE 3
We show that the contact between particles 672 and 827 shown in Figure 3C is a non-persistent contact. The magnitude of the normal component of the contact force, fn, on this contact is plotted in Figure 3B. It reversibly drops to zero. To further show that the contact actually opens, we show snapshots of the system in five typical states A to E in Figure 3C. In state D, as shown by the snapshot in the third column of Figure 3C, the contact is clearly opened. In state B, the clearly visible photoelastic fringes confirm that the contact is carrying finite forces. Thus, this contact opens and closes reversibly in a shear cycle, and is called a non-persistent contact in this paper. See Supplementary Video S1 for a video of this process. Note that the evolution of fn is consistent with the global shear. The branch vector pointing to contact 672–827 from the center of particle 672 is roughly parallel to the y′ direction (see Figure 1A), along which the system is compressed from state A to state B and is stretched from B to D. Accordingly, fn on contact 672–827 grows from A to B and drops from B to D, during which it opens.
To quantitatively classify the contacts and characterize the behavior of the non-persistent contacts, we consider two characteristic quantities: (1) the fraction of time that the contact is open, ρopen, defined as the number of steps in one complete shear cycle for which fn < 0.01 N divided by the total number of steps in the cycle; and (2) the maximum value of fn for this contact over the whole shear cycle, fn,max. Figure 3D shows a scatter plot of fn,max and ρopen for all contacts measured over three consecutive shear cycles.
By definition, a non-persistent contact should have a non-zero ρopen and a non-zero fn,max. For present purposes, we detect the non-persistent contacts with ρopen > 0.1 and fn,max > 0.06 N (the points in the purple region in Figure 3D. We intentionally choose a large value for ρopen to ensure that the contact actually opens and a large value for fn,max to ensure that the contact actually closes. We emphasize that the goal here is to demonstrate the existence of non-persistent contacts and their relevance in contributing to the elastic responses of the packing. This conservative classification method ensures that a positive result is meaningful.
The threshold for fn,max is chosen based on the uncertainty of our force solving algorithm. Figure 3F plots the probability density function of the difference between action and reaction contact forces determined by our fitting algorithm for all contacts detected in 61 jammed states over the three shear cycles. The width of this distribution is an estimation of the uncertainty of our force measurements because Newton’s third law ensures that these differences must actually be zero. A Gaussian fit gives a width around 0.03 N. Thus, a contact is convincingly closed at least once in a shear cycle if fn,max > 0.06 N.
The threshold for ρopen is chosen based on our sampling frequency. For a shear cycle with strain amplitude δγ = 0.95%, there are 20 quasi-static data collection steps per cycle, as shown in Figure 3A. Thus, the resolution of ρopen is 1/20. Therefore, we expect that a threshold value of 1/10 probes contacts that actually opens during a shear cycle. In Figure 3D there appear to be some data points with ρopen between 0 and 1/20 because they are averaged values over three shear cycles.
The detected non-persistent contacts for the example ultrastable state are plotted on top of the photoelastic fringes in Figure 3G. These contacts are scattered in space and do not form a percolating network. As a measure of the prevalence of non-persistent contacts, we calculate fnpc, the number of non-persistent contacts divided by the total number of contacts that were ever closed during a shear cycle. Figure 5A plots fnpc as a function of the initial shear strain γI for ultrastable states formed using same δγ = 0.95%. We find that fnpc is larger for smaller γI and can be as large as about 10%.
Notably, reversible plastic events observed in frictionless systems (e.g. in two-dimensional foams []) also involve reversibly activated inter-particle contacts. A distinction in our frictional system is that there is no obvious T1 event as observed in Ref. []. For example, the reversible activation of contact 672–827 in Figure 3C is not accompanied by a neighbor switching event for the four particle 672, 827, 381 and 584. In other words, there is no obvious local plastic event triggered by an opening of a non-persistent contact. Thus, the particles with non-persistent contacts are not equivalent to bucklers in isostatic frictionless packings [], where breaking a contact will immediately induce the formation of a new contact.
3.2.2 Persistent contacts
For an ultrastable state, most of the contacts are persistent; they remain closed throughout the shear cycle. One example is the contact between particles 672 and 381, shown in Figure 3C. It always bears a finite normal force component fn under cyclic shear, as shown in Figure 3B. In practice, we classify a contact as persistent if the minimum normal force during a whole cycle, fn,min, is greater than 0.06 N. The threshold is chosen as twice our force measurement uncertainty to make sure that fn,min is convincingly larger than 0. Figure 3E plots all the contacts for the example ultrastable state in the fn,min and fn,max plane; all the contacts in the light blue region are classified as persistent contacts.
Unlike the scattered non-persistent contacts, the persistent contacts form a percolating network. Figure 3G shows the network formed by the persistent contacts (light blue lines). Notably, the persistent contact network formed shown in Figure 5G features large holes reminiscent of the sponge-like structures revealed by rigidity analysis [].
We find that there is no measurable sliding at the persistent contacts. The ultrastable states show the same strobed state over thousands of shear cycles, and the particles do not rotate from cycle to cycle, as can be seen from the supplementary videos of Ref. []. The lack of rotation contrasts with recent experiments where particles do not move much but rotate significantly under cyclic loading, displaying contact sliding that leads to energy dissipation [, ]. For an ultrastable state, there should be only two possible cases for a given persistent contact: (1) there is no sliding at contact; or (2) the two particles slide against each other during a shear cycle but return to the same position and stress states after a complete cycle. In the latter case, the tangential to normal force ratio μ = ft/fn should reach both + μs and − μs in a cycle, where μs = 0.87 is the static friction coefficient of the particles. The evolution of μ on an example contact is shown in Figure 4A. We measure the magnitude of variation of μ, denoted as Δμ in Figure 4A. The example contact plotted in Figure 4A has the largest Δμ among all persistent contacts in the ultrastable state shown in Figure 4. In an ultrastable state, a contact that slides must have Δμ ≈ 2μs. (If sliding occurs in only one direction during the cycle, it would necessarily produce relative rotations of the two particles, which is not observed.) Figure 4C plots the number statistics of Δμ for all persistent contacts in the example ultrastable state shown in Figure 3. Clearly, almost no persistent contact has Δμ near 2μs, suggesting that these is no reversible sliding. We note that rolling without sliding is allowed and was observed at persistent contacts, but a complete characterization of rolling is beyond the scope of this paper. For completeness, we show the evolution of some example persistent contacts in the (fn, ft) space in Figure 4B. We did not perform the same analysis to non-persistent contacts because the measured μ becomes unreliable for weak forces due to the resolution limit of photoelastic force measurements. Extremely slow accumulation of plasticity induced by ratchet-like sliding on weak contacts, as found in numerical simulations [, ], cannot be completely ruled out.
FIGURE 4
The conditions that we used to identify persistent contacts and non-persistent contacts give high true positive ratio and a small true negative ratio. There are many contacts in our system that do not satisfy either criterion. Whether these contacts are persistent or non-persistent could conceivably be resolved in future experiments with higher force, distance, and time resolutions. These unclassified contacts are not plotted in Figure 3G.
3.2.3 Contribution to the global elastic modulus
We now show that the non-persistent contacts contribute a non-negligible amount to the emergent global elastic modulus. We calculate the shear stress contributed from the non-persistent contacts (npc) aswhere the summation is only over all non-persistent contacts. Figure 5B plots the total σxy and for an example ultrastable state under cyclic shear. The contribution from non-persistent contacts to the shear modulus, Gnpc, is the slope of , as sketched in Figure 5D. Figure 5C plots the ratio Gnpc/G for ultrastable states formed under δγ = 0.95% but different γI. We see that the contribution to the shear modulus from the non-persistent contacts can be as large as 10% for γI near the onset value for creating ultrastable states, and the actual contribution could be even larger because some unclassified contacts are likely non-persistent ones. Thus, the non-persistent contacts make an appreciable contribution to the mechanical response of our ultrastable packings. In addition, the importance of the non-persistent contacts suggests that grains that might be identified as rattlers during some portion of the cycle may actually contribute to the elastic behavior observed for finite amplitude shear deformations. It also worth mentioning here that the scaling relation of Eq. 3 is a global relation that contains contributions from all contacts. A detailed discussion of the separate contributions of persistent and non-persistent contacts to G(p) is beyond the scope of the present paper.
FIGURE 5
3.3 Particle center trajectories
Characterizing the particle motion within a shear cycle gives valuable insights into the nature of the mechanical responses of the packing. If an ultrastable packing deforms like a linear elastic continuum, all the particle displacements should define an affine deformation field. However, in our ultrastable states, the particles display clearly detectable, spatially correlated, non-affine displacements. In addition, some particle trajectories form loops with measurable enclosed area.
Figure 6A shows all the particle center trajectories in a shear cycle for the example ultrastable state shown in Figure 3. The trajectories are nearly vertical lines that appear roughly consistent with an affine simple shear deformation field, suggesting that both the non-affine displacements and the enclosed loop areas are small. Figure 6B plots the non-affine displacements of particles measured in the strain interval between state B to state D in Figure 3A. Here, the non-affine displacement of particle i, δrna,i, is defined aswhere Θ(x) is the Heaviside step function, xi is the x coordinate of particle i, δri is the real displacement of particle i, rs is the radius of the small particle, and Np is the total number of particles. In Figure 6B, the arrows are colored according to the magnitude of the non-affine displacements |δrna|, and the lengths of the arrows are 20 times |δrna|. While the particles go back exactly to the same position after a full cycle, it is clear that their displacements within a shear cycle often contain significant non-affine components. In the future, it could be interesting to compare our results to a recent theory which considers non-affine deformations while assuming no sliding at frictional contacts [].
FIGURE 6
We further show that there are measurable loops formed by particle center trajectories, and also by the non-affine center trajectories. Figure 6 shows center trajectories (C-E) and their corresponding non-affine center trajectories (F-H) calculated from Eq. 5 for three example particles over three consecutive shear cycles. Note that the non-affine displacements appear noisier because they are near the accuracy of our particle center detection (about 0.01ds). The trajectory in (C) clearly is a loop and the one in (E) does not show measurable area. For the non-affine trajectories, only (G) shows a noticeable loop above the noise level. In numerical simulations of frictional granular systems, loops in particle trajectories [] and in non-affine trajectories [] were observed. In particular, the areas of these loops were found to obey a scaling relation with the elastic moduli of the system []. While relating these loops to global elastic responses is beyond the scope of the present paper, our work establishes their existence in this experimental frictional granular system.
3.4 Stress correlations and emergent properties
To obtain more insight on the stress responses of the ultrastable states, we consider the Vector Charge Theory of Granular mechanics and dynamics (VCTG) [, ]. VCTG relates features of the stress correlations in the continuum limit to emergent elastic properties of the packing. We show that such a theory predicts forms of stress correlations that reasonably match our data. Notably, the results obtained from fitting the data to VCTG predictions uncover a feature that distinguishes the ultrastable states from the original shear-jammed states formed by initial shear alone.
3.4.1 Defining ensembles with similar stress states
To obtain the correlation functions, we group ultrastable states with similar stress fields to form ensembles and calculate the averaged correlation functions over these ensembles. As shown in Figures 7A,B, we show that all ultrastable states fall roughly on a same curve when plotting the non-rattler contact number Znrversus pressure p or when plotting the shear stress σxyversusp. This observation suggests that we group states according to p, and the states with similar p will have similar σxy and Znr. Specifically, we group states with a pressure interval of 3 N/m, and the averaged state variables for ultrastable states in these intervals are plotted using purple circles in Figures 7A,B. The error bars mark the standard deviations.
FIGURE 7
For completeness, we also plot data from the original shear-jammed states that are formed by initial shear only in Figures 7A,B. Comparing ultrastable states and original states provides additional insights into how cyclic shear modifies the mechanical properties of a jammed granular packing. Notably, for packings with similar p, ultrastable states usually contain more contacts and exhibit lower shear stress, suggesting that they are more stable and less anisotropic. We also group the original states according to intervals of pressure and apply the same stress correlation analysis below.
3.4.2 Stress correlation functions
We calculate the correlation functions between components of the stress tensor following the procedure detailed in Ref. [] and Ref. []. We note that it is more convenient to examine the stress correlation functions in a reference frame x′y′ that is rotated 45° from the original reference frame xy (see Figure 1). After the rotation, x′ is the principal dilation direction and y′ is the principal compression direction of the initial shear. The force chains in Figure 3G mostly align with direction y′. We consider below correlation functions and stress tensor expressions in this rotated frame.
The correlation functions in Fourier space are calculated as follows [].where δσi′j′ = σi′j′ − ⟨σi′j′⟩ and ⟨σi′j′⟩ is the spatially averaged value in a packing and the ⟨⟩ in Eq. 6 refers to average over different packings in an ensemble, and
We have used primed indices to emphasize that all the calculations are done in the rotated frame x′y′. More details on the calculation of the stress correlation functions are provided in the Supplementary Material S1.
As an example, Figure 8 shows the six stress-stress correlation functions in Fourier space obtained from an ensemble that contains 6 packings with a averaged pressure p = 8.6 ± 0.7 N/m. We note that the general features of the correlation functions are consistent with those reported in Refs. [, ], including the pinch-point singularities at |q| → 0 and the obvious radial variation for wavelengths shorter than about 4ds, where the continuuum theory is affected by the granularity of the medium. In Figure 8 all correlation functions are normalized by , a parameter in the VCTG fitting form, which is presented below (Eq. 17). Correlation functions for ultrastable states with other stress states share similar features.
FIGURE 8
3.4.3 Elastic moduli from stress correlations
The emergent elastic moduli appear in the VCTG predictions of stress correlations, and we compute these by fitting the data. We first extract the angular dependence of the correlation functions in the long-wavelength limit and compare them to the VCTG predictions. Specifically, for each correlation function, we average the data in a radial range between 2π/6ds and 2π/16ds and plot the radially averaged data versus the azimuthal angle θ. Note that 16ds is the size of the region of interest that we used to calculate these correlation functions, and 4ds is the length below which the correlation functions clearly deviates from the values at smaller |q|. The radially averaged correlation functions are plotted in Figure 9 for ultrastable states with different stress states as labeled by color. Note that each curve is averaged over several independent experimental realizations. We believe the scattering of data originates from that the number of packings used in the averaging process is rather small, and our system is not large enough. Nonetheless, we find that these curves can be reasonably characterized by the VCTG predictions.
FIGURE 9
A key prediction from the VCTG theory is the form of the stress correlation functions in the long-wavelength limit [].whereand
Einstein notation applies to Eqs. 8–9, where ϵi′j′ and δi′j′ denote the Levi-Civita symbol and Kronecker delta. Detailed derivations of Eq. 8 can be found in Refs. [, ]. Here, the only unknown variables are the elements of the 4-rank tensor Λ. These elements will be obtained by fitting the experimentally calculated stress correlation functions to Eq. 8. In the VCTG framework, Λ maps to the inverse elastic constant tensor. Using the symmetries of the elastic constant tensor Λijkl = Λjikl and Λijkl = Λijlk, there iswhere
For subscripts of Λ we have used 1 and 2 to represent x′ and y′ for simplicity. Considerations of additional symmetries of the system may lead to particular forms of the elastic constant tensor that simplify the analysis [, ]. In shear-jammed systems, it is reasonable to assume that the elastic moduli have uniaxial symmetry []. Note that in these jammed states created by external stresses, the elastic moduli are determined by the geometry and topology of the force-bearing network that emerges from the jamming process []. In the Voigt form [], such an elastic modulus tensor reads:where Ex′ and Ey′ are the Young’s moduli along x′, y′ directions, while νx′ and νy′ stand for the Poisson’s ratios along x′ and y′ directions. G′ is a shear modulus such that in the x′y′ coordinate system there is σx′y′ = G′ɛx′y′. G′ is not the shear modulus G which follows σxy = Gγ. In addition, our imposed simple shear strain field has ɛx′y′ = 0. Thus, G′ can not be extracted from our measured stress-strain curves.
The inverse elastic tensor then reads
Note that this tensor maps to the Λ tensor in the VCTG framework as following,
Comparing Eq. 14 and Eq. 15 we get Λ1211 = Λ1222 = Λ1112 = Λ2212 = 0. Thus, according to Eq. 12, we expect and
Thus, we fit the correlation functions Ci′j′k′l′(θ) to the following form.
Note that, we fit all six correlation functions together with three parameters , , and . In Figure 9 we plot both the raw data and the fitted curves. All data and fitted curves are normalized by the fitting parameter which depends on p. Notably, the fitted curves matches with the experimental data reasonably well. In addition, it also appears that most of the data collapse after the normalization by , suggesting roughly constant and for ultrastable states with different p.
3.4.4 Emergent elastic response
The emergent elastic moduli, determined from fitting the measured stress-stress correlations to the VCTG predictions, depend on preparation protocols and the average stress state of a shaer-jammed solid. Here, we analyze the dependence of the three fitting parameters, , , and , on the pressure p of the ultrastable states. We note that as the system unjams at p = 0, what we are considering is scaling of the emergent properties near jamming. However, according to Eq. 16, while and can be directly related to the two Young’s moduli, we can not extract the shear modulus and the Poisson’s ratios from just the three fitting parameters. There are 5 unknowns but only three equations in Eq. 16. Reference [] demonstrates that additional equations can be obtained by considering material responses to additionally applied forces. In the present work, we report only scalings for the fitting parameters, and we leave the actual solutions for the shear modulus and Poisson’s ratios to future investigations.
We plot and as functions of p in Figure 10A. While for all p, they both decay with increasing p. Notably, they appear to follow power laws with same exponent ∼ − 1.6. The divergence of these parameters suggest the Young’s moduli vanish at the jamming point. We plot the ratios and in Figure 10B. Both ratios appear to be roughly constant. Notably, , meaning that the system is always stiffer along the y′ direction. It is interesting that . We do not yet have a clear understanding of this feature.
FIGURE 10
3.4.5 A feature that distinguish ultrastable states from original shear-jammed states
We show that the ratio can be used as a indicator to distinguish ultrastable states and the original shear-jammed states that are formed by initial shear alone. Figures 10A,B also show fit results obtained by performing same analysis on ensembles of the original shear-jammed states. Interestingly and follow power law scaling versusp with same exponent as the ultrastable states. We find, however, that the ratio for original shear-jammed states, which is much larger than for the ultrastable states. This means the original states have much more anisotropic elastic properties compared to the ultrastable states. We thus show that small-amplitude cyclic shearing changes the elastic response of a jammed packing. We emphasize that this observation is not equivalent to the changes of stress states as can be evidenced in Figure 7B. Instead, it highlights that elasticity of these shear-jammed solids is truly an emergent phenomenon reflecting a rigidity that emerges from the complex interplay of local and global force and torque balance contstraints [, ].
4 Concluding discussion
In summary, we report a set of analyses on both global and local features of ultrastable shear-jammed granular materials in response to cyclic shear. We present three major findings.
First, we show that the emergent shear modulus G for ultrastable states formed by different γI and δγ falls on a single curve when plotted versus pressure p. A critical scaling near jamming between G and p is examined extensively in numerical simulations [, , , , , ], and is of key interests in the scaling theories of jamming []. Notably, the ultrastable states follow G ∼ pβ with β ≈ 0.5, consistent with a numerical simulation with particles having similar friction coefficient and contact force law []. To our knowledge, the range of boundary strain within which a frictional system behaves elastically is usually very small because boundary strain may induce sliding at contacts [–]. Thus the shear modulus has typically been determined from measurements of sound speeds [] rather than from stress-strain curves. Consistent with this picture, the original shear-jammed states in our experiments were observed to deform plastically under any given cyclic strain amplitude []. As previously reported, highly elastic ultrastable states emerged under cyclic shearing, which induces changes in the distribution of friction forces at contacts []. We have now found that there is no measurable sliding for the persistent contacts that carry the majority of the forces in these ultrastable states.
Second, we find nontrivial grain scale motions within a shear cycle in an ultrastable state. A measurable fraction of contacts open and close reversibly during a cycle, and these contacts make a non-negligible contribution to the emergent elastic modulus. Thus predictions based on a contact network with a fixed geometry presumably cannot completely account for the macroscopic elasticity of these states. It is known that the distribution of small inter-particle gaps and weak contact forces are intimately connected to packing stability [, , ]. Our work demonstrates that reversible activation of these gaps may lead to non-trivial dynamical phases in a frictional, shear-jammed system. We also observe non-affine particle displacements, with some particles moving around loops with finite enclosed area. It would be interesting to compare the observed particle displacement fields to the low-frequency vibrational modes that can be calculated from the experimental data [], where one may find analogies to features found in model glasses, such as string-like dynamical defects [].
Third, we examine the relation between the spatial stress fluctuations and the emergent elastic constants of the ultrastable states from the perspective of the Vector Charge Theory of Granular mechanics and dynamics (VCTG) [, ]. In the long-wavelength limit, the stress-stress correlation functions measured from ultrastable and original shear-jammed states matches well with the predictions by VCTG for an anisotropic system with uniaxial symmetry. Fitting our data to the theory, we extract the values of three parameters. Two of these are the Young’s moduli Ex′ and Ey′, and the third is a linear combination of the Poisson ratios νx′, νy′ and a shear modulus G′. Note that x′ and y′ are the principal dilation and compression directions of the initial shear. We find that, for both original shear-jammed states and the ultrastable states, Ex′ and Ey′ scale as power-laws with pressure p, sharing same exponent α ≈ 1.6. The vanishing of the Young’s moduli as p → 0 is consistent qualitatively with the vanishing of the shear modulus G measured independently from the stress-strain curves. The relation between the exponent α and the exponent β that links G and p is an interesting topic for further investigation. We note that the elastic moduli obtained from VCTG fittings are linear elastic constants, while stress-strain curves contain contributions from non-linear features like the non-persistent contacts. The ratio Ex′/Ey′ is always at least twice as large for original shear-jammed states as for the ultrastable states, suggesting that small-amplitude cyclic shearing significantly alters the elastic properties of a jammed packing. In addition, the ratio Ex′/Ey′ does not approach 1 as p → 0, suggesting that the system remains anisotropic at jamming point, which is a feature of the shear jamming transition [, , , ]. Additional experiments and analysis probing the system response to a point force [] may help to determine the Poisson ratios and the shear modulus G′.
Statements
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Author contributions
YiZ, BC, and JS contributed to conception and design of the study and collaborated on the interpretation of the experimental results. YiZ, YuZ, DW, and HZ designed the experimental procedure and calibrated the apparatus. YiZ and YuZ conducted the experiments. YiZ performed the data analysis. YiZ, BC, and JS wrote the first draft of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.
Funding
This work was primarily supported by NSF grant DMR-1809762. BC was supported by NSF grants CBET-1916877, and CMMT-2026834, and BSF-2016188.
Acknowledgments
YiZ thanks Yinqiao Wang for helpful discussions on calculating the correlation functions. YiZ thanks Peter K. Morse, Shuai Zhang, Yuliang Jin, and Deng Pan for helpful discussions about the scaling of shear modulus.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphy.2022.1048683/full#supplementary-material
References
1.
ZhaoYZhaoYWangDZhengHChakrabortyBSocolarJES. Ultrastable shear-jammed granular material. Phys Rev X (2022) 12:031021. 10.1103/PhysRevX.12.031021
2.
NampoothiriJNWangYRamolaKZhangJBhattacharjeeSChakrabortyB. Emergent elasticity in amorphous solids. Phys Rev Lett (2020) 125:118002. 10.1103/PhysRevLett.125.118002
3.
LiuAJNagelSR. Jamming is not just cool any more. Nature (1998) 396:21–2. 10.1038/23819
4.
MajmudarTSSperlMLudingSBehringerRP. Jamming transition in granular systems. Phys Rev Lett (2007) 98:058001. 10.1103/PhysRevLett.98.058001
5.
van HeckeM. Jamming of soft particles: Geometry, mechanics, scaling and isostaticity. J Phys: Condens Matter (2009) 22:033101. 10.1088/0953-8984/22/3/033101
6.
SarkarSBiDZhangJBehringerRPChakrabortyB. Origin of rigidity in dry granular solids. Phys Rev Lett (2013) 111:068301. 10.1103/PhysRevLett.111.068301
7.
LudingS. So much for the jamming point. Nat Phys (2016) 12:531–2. 10.1038/nphys3680
8.
BehringerRPChakrabortyB. The physics of jamming for granular materials: A review. Rep Prog Phys (2018) 82:012601. 10.1088/1361-6633/aadc3c
9.
PetersIRMajumdarSJaegerHM. Direct observation of dynamic shear jamming in dense suspensions. Nature (2016) 532:214–7. 10.1038/nature17167
10.
HanEJamesNMJaegerHM. Stress controlled rheology of dense suspensions using transient flows. Phys Rev Lett (2019) 123:248002. 10.1103/PhysRevLett.123.248002
11.
MorrisJF. Shear thickening of concentrated suspensions: Recent developments and relation to other phenomena. Annu Rev Fluid Mech (2020) 52:121–44. 10.1146/annurev-fluid-010816-060128
12.
BiDZhangJChakrabortyBBehringerRP. Jamming by shear. Nature (2011) 480:355–8. 10.1038/nature10667
13.
VinuthaHASastryS. Disentangling the role of structure and friction in shear jamming. Nat Phys (2016) 12:578–83. 10.1038/nphys3658
14.
WangDRenJDijksmanJAZhengHBehringerRP. Microscopic origins of shear jamming for 2d frictional grains. Phys Rev Lett (2018) 120:208004. 10.1103/PhysRevLett.120.208004
15.
ZhaoYBarésJZhengHSocolarJESBehringerRP. Shear-jammed, fragile, and steady states in homogeneously strained granular materials. Phys Rev Lett (2019) 123:158001. 10.1103/PhysRevLett.123.158001
16.
VinuthaHSastryS. Force networks and jamming in shear-deformed sphere packings. Phys Rev E (2019) 99:012123. 10.1103/physreve.99.012123
17.
OtsukiMHayakawaH. Shear jamming, discontinuous shear thickening, and fragile states in dry granular materials under oscillatory shear. Phys Rev E (2020) 101:032905. 10.1103/PhysRevE.101.032905
18.
KumarNLudingS. Memory of jamming–multiscale models for soft and granular matter. Granular Matter (2016) 18:58. 10.1007/s10035-016-0624-2
19.
Baity-JesiMGoodrichCPLiuAJNagelSRSethnaJP. Emergent so(3) symmetry of the frictionless shear jamming transition. J Stat Phys (2017) 167:735–48. 10.1007/s10955-016-1703-9
20.
BabuVPanDJinYChakrabortyBSastryS. Dilatancy, shear jamming, and a generalized jamming phase diagram of frictionless sphere packings. Soft Matter (2021) 17:3121–7. 10.1039/D0SM02186E
21.
EdigerMD. Perspective: Highly stable vapor-deposited glasses. J Chem Phys (2017) 147:210901. 10.1063/1.5006265
22.
KimSHilgenfeldtS. Structural measures as guides to ultrastable states in overjammed packings. Phys Rev Lett (2022) 129:168001. 10.1103/PhysRevLett.129.168001
23.
RegevILookmanTReichhardtC. Onset of irreversibility and chaos in amorphous solids under periodic shear. Phys Rev E (2013) 88:062401. 10.1103/PhysRevE.88.062401
24.
FioccoDFoffiGSastryS. Oscillatory athermal quasistatic deformation of a model glass. Phys Rev E (2013) 88:020301. 10.1103/physreve.88.020301
25.
Hima NagamanasaKGokhaleSSoodAKGanapathyR. Experimental signatures of a nonequilibrium phase transition governing the yielding of a soft glass. Phys Rev E (2014) 89:062308. 10.1103/PhysRevE.89.062308
26.
KeimNCArratiaPE. Mechanical and microscopic properties of the reversible plastic regime in a 2d jammed material. Phys Rev Lett (2014) 112:028302. 10.1103/PhysRevLett.112.028302
27.
KawasakiTBerthierL. Macroscopic yielding in jammed solids is accompanied by a nonequilibrium first-order transition in particle trajectories. Phys Rev E (2016) 94:022615. 10.1103/PhysRevE.94.022615
28.
NessCCatesME. Absorbing-state transitions in granular materials close to jamming. Phys Rev Lett (2020) 124:088004. 10.1103/PhysRevLett.124.088004
29.
ReichhardtCRegevIDahmenKOkumaSReichhardtCJO. Perspective on reversible to irreversible transitions in periodic driven many body systems and future directions for classical and quantum systems. arXiv preprint arXiv:2211.03775 (2022).
30.
CortéLChaikinPMGollubJPPineDJ. Random organization in periodically driven systems. Nat Phys (2008) 4:420–4. 10.1038/nphys891
31.
RoyerJRChaikinPM. Precisely cyclic sand: Self-organization of periodically sheared frictional grains. Proc Natl Acad Sci U S A (2015) 112:49–53. 10.1073/pnas.1413468112
32.
OtsukiMHayakawaH. Shear modulus and reversible particle trajectories of frictional granular materials under oscillatory shear. Eur Phys J E (2021) 44:70. 10.1140/epje/s10189-021-00075-0
33.
NagasawaKMiyazakiKKawasakiT. Classification of the reversible–irreversible transitions in particle trajectories across the jamming transition point. Soft matter (2019) 15:7557–66. 10.1039/c9sm01488h
34.
NampoothiriJND’EonMRamolaKChakrabortyBBhattacharjeeS. Tensor electromagnetism and emergent elasticity in jammed solids. arXiv preprint arXiv:2204.11811 (2022).
35.
RenJDijksmanJABehringerRP. Reynolds pressure and relaxation in a sheared granular system. Phys Rev Lett (2013) 110:018302. 10.1103/PhysRevLett.110.018302
36.
RenJ. Nonlinear dynamics and network properties in granular materials under shear. Ph.D. thesis. Durham, NC: Duke University (2013).
37.
SarkarSBiDZhangJRenJBehringerRPChakrabortyB. Shear-induced rigidity of frictional particles: Analysis of emergent order in stress space. Phys Rev E (2016) 93:042901. 10.1103/PhysRevE.93.042901
38.
MajmudarTSBehringerRP. Contact force measurements and stress-induced anisotropy in granular materials. Nature (2005) 435:1079–82. 10.1038/nature03805
39.
ChristoffersenJMehrabadiMMNemat-NasserS. A micromechanical description of granular material behavior. J Appl Mech (1981) 48:339–44. 10.1115/1.3157619
40.
RadjaiFWolfDEJeanMMoreauJJ. Bimodal character of stress transmission in granular packings. Phys Rev Lett (1998) 80:61–4. 10.1103/PhysRevLett.80.61
41.
O’HernCSSilbertLELiuAJNagelSR. Jamming at zero temperature and zero applied stress: The epitome of disorder. Phys Rev E (2003) 68:011306. 10.1103/PhysRevE.68.011306
42.
SomfaiEvan HeckeMEllenbroekWGShundyakKvan SaarloosW. Critical and noncritical jamming of frictional grains. Phys Rev E (2007) 75:020301. 10.1103/PhysRevE.75.020301
43.
GoodrichCPLiuAJSethnaJP. Scaling ansatz for the jamming transition. Proc Natl Acad Sci U S A (2016) 113:9745–50. 10.1073/pnas.1601858113
44.
MakseHAGlandNJohnsonDLSchwartzL. Granular packings: Nonlinear elasticity, sound propagation, and collective relaxation dynamics. Phys Rev E (2004) 70:061302. 10.1103/PhysRevE.70.061302
45.
IshimaDHayakawaH. Scaling laws for frictional granular materials confined by constant pressure under oscillatory shear. Phys Rev E (2020) 101:042902. 10.1103/PhysRevE.101.042902
46.
VanderWerfKBoromandAShattuckMDO’HernCS. Pressure dependent shear response of jammed packings of frictionless spherical particles. Phys Rev Lett (2020) 124:038004. 10.1103/PhysRevLett.124.038004
47.
WangPZhangSTuckmanPOuelletteNTShattuckMDO’HernCS. Shear response of granular packings compressed above jamming onset. Phys Rev E (2021) 103:022902. 10.1103/PhysRevE.103.022902
48.
TreadoJDWangDBoromandAMurrellMPShattuckMDO’HernCS. Bridging particle deformability and collective response in soft solids. Phys Rev Mater (2021) 5:055605. 10.1103/PhysRevMaterials.5.055605
49.
LundbergMKrishanKXuNO’HernCSDenninM. Reversible plastic events in amorphous materials. Phys Rev E (2008) 77:041505. 10.1103/PhysRevE.77.041505
50.
CharbonneauPCorwinEIParisiGZamponiF. Jamming criticality revealed by removing localized buckling excitations. Phys Rev Lett (2015) 114:125504. 10.1103/PhysRevLett.114.125504
51.
LiuKKollmerJEDanielsKESchwarzJMHenkesS. Spongelike rigid structures in frictional granular packings. Phys Rev Lett (2021) 126:088002. 10.1103/PhysRevLett.126.088002
52.
PeshkovAGirvanMRichardsonDCLosertW. Reversibility of granular rotations and translations. Phys Rev E (2019) 100:042905. 10.1103/PhysRevE.100.042905
53.
BensonZAPeshkovAYunger HalpernNRichardsonDCLosertW. Experimentally measuring rolling and sliding in three-dimensional dense granular packings. Phys Rev Lett (2022) 129:048001. 10.1103/PhysRevLett.129.048001
54.
Alonso-MarroquínFHerrmannHJ. Ratcheting of granular materials. Phys Rev Lett (2004) 92:054301. 10.1103/PhysRevLett.92.054301
55.
McNamaraSGarcía-RojoRHerrmannHJ. Microscopic origin of granular ratcheting. Phys Rev E (2008) 77:031304. 10.1103/PhysRevE.77.031304
56.
IshimaDSaitohKOtsukiMHayakawaH. Theory of rigidity and density of states of two-dimensional amorphous solids of dispersed frictional grains in the linear response regime. arXiv preprint arXiv:2207.06632 (2022).
57.
OttoMBouchaudJPClaudinPSocolarJES. Anisotropy in granular media: Classical elasticity and directed-force chain network. Phys Rev E (2003) 67:031302. 10.1103/PhysRevE.67.031302
58.
MorsePKRoySAgoritsasEStaniferECorwinEIManningML. A direct link between active matter and sheared granular systems. Proc Natl Acad Sci U S A (2021) 118:e2019909118. 10.1073/pnas.2019909118
59.
PanDMengFJinY. Shear hardening in frictionless amorphous solids near the jamming transition. arXiv preprint arXiv:2208.08793 (2022).
60.
AssimakiDKauselEWhittleA. Model for dynamic shear modulus and damping for granular soils. J Geotech Geoenviron Eng (2000) 126126:85910–69. 10.1061/(asce)1090-0241(2000)126:10(859)
61.
Alonso-MarroquínFLudingSHerrmannHJVardoulakisI. Role of anisotropy in the elastoplastic response of a polygonal packing. Phys Rev E (2005) 71:051304. 10.1103/PhysRevE.71.051304
62.
OtsukiMHayakawaH. Discontinuous change of shear modulus for frictional jammed granular materials. Phys Rev E (2017) 95:062902. 10.1103/PhysRevE.95.062902
63.
SunAWangYChenYShangJZhengJYuSet alTurbulent-like velocity fluctuations in two-dimensional granular materials subject to cyclic shear. Soft Matter (2022). 10.1039/D1SM01516H
64.
BabuVSastryS. Criticality and marginal stability of the shear jamming transition of frictionless soft spheres. Phys Rev E (2022) 105:L042901. 10.1103/PhysRevE.105.L042901
65.
WangYShangJJinYZhangJ. Experimental observations of marginal criticality in granular materials. Proc Natl Acad Sci U S A (2022) 119:e2204879119. 10.1073/pnas.2204879119
66.
ZhangLZhengJWangYZhangLJinZHongLet alExperimental studies of vibrational modes in a two-dimensional amorphous solid. Nat Commun (2017) 8:67. 10.1038/s41467-017-00106-5
67.
HuYCTanakaH. Origin of the boson peak in amorphous solids. Nat Phys (2022) 18:669–77. 10.1038/s41567-022-01628-6
68.
ChenSBertrandTJinWShattuckMDO’HernCS. Stress anisotropy in shear-jammed packings of frictionless disks. Phys Rev E (2018) 98:042906. 10.1103/PhysRevE.98.042906
69.
XiongFWangPClarkAHBertrandTOuelletteNTShattuckMDet alComparison of shear and compression jammed packings of frictional disks. Granular Matter (2019) 21:109. 10.1007/s10035-019-0964-9
Summary
Keywords
reversibility, jamming, elasticity, friction, granular microstructure, cyclic shear, stress correlation, yielding
Citation
Zhao Y, Zhao Y, Wang D, Zheng H, Chakraborty B and Socolar JES (2022) Microscopic reversibility and emergent elasticity in ultrastable granular systems. Front. Phys. 10:1048683. doi: 10.3389/fphy.2022.1048683
Received
19 September 2022
Accepted
14 November 2022
Published
28 November 2022
Volume
10 - 2022
Edited by
Yang Jiao, Arizona State University, United States
Reviewed by
Yuanchao Hu, Arizona State University, United States
Qinyi Liao, Institute of Theoretical Physics (CAS), China
Updates
Copyright
© 2022 Zhao, Zhao, Wang, Zheng, Chakraborty and Socolar.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yiqiu Zhao, yiqiuzhao@ust.hk; Joshua E. S. Socolar, socolar@duke.edu
This article was submitted to Soft Matter Physics, a section of the journal Frontiers in Physics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.