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.
more »
« less
Scaling regimes and fluctuations of observables in computer glasses approaching the unjamming transition
Under decompression, disordered solids undergo an unjamming transition where they become under-coordinated and lose their structural rigidity. The mechanical and vibrational properties of these materials have been an object of theoretical, numerical, and experimental research for decades. In the study of low-coordination solids, understanding the behavior and physical interpretation of observables that diverge near the transition is of particular importance. Several such quantities are length scales (ξ or l) that characterize the size of excitations, the decay of spatial correlations, the response to perturbations, or the effect of physical constraints in the boundary or bulk of the material. Additionally, the spatial and sample-to-sample fluctuations of macroscopic observables such as contact statistics or elastic moduli diverge approaching unjamming. Here, we discuss important connections between all of these quantities and present numerical results that characterize the scaling properties of sample-to-sample contact and shear modulus fluctuations in ensembles of low-coordination disordered sphere packings and spring networks. Overall, we highlight three distinct scaling regimes and two crossovers in the disorder quantifiers χz and χμ as functions of system size N and proximity to unjamming δz. As we discuss, χX relates to the standard deviation σX of the sample-to-sample distribution of the quantity X (e.g., excess coordination δz or shear modulus μ) for an ensemble of systems. Importantly, χμ has been linked to experimentally accessible quantities that pertain to sound attenuation and the density of vibrational states in glasses. We investigate similarities and differences in the behaviors of χz and χμ near the transition and discuss the implications of our findings on current literature, unifying findings in previous studies.
more »
« less
- Award ID(s):
- 1951921
- PAR ID:
- 10518315
- Publisher / Repository:
- Journal of Chemical Physics
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 160
- Issue:
- 3
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We perform computational studies of jammed particle packings in two dimensions undergoing isotropic compression using the well-characterized soft particle (SP) model and deformable particle (DP) model that we developed for bubbles and emulsions. In the SP model, circular particles are allowed to overlap, generating purely repulsive forces. In the DP model, particles minimize their perimeter, while deforming at the fixed area to avoid overlap during compression. We compare the structural and mechanical properties of jammed packings generated using the SP and DP models as a function of the packing fraction ρ, instead of the reduced number density φ. We show that near jamming onset the excess contact number Δz=z-z J and shear modulus G scale as Δρ 0.5 in the large system limit for both models, where Δρ=ρ-ρ J and z J ≈4 and ρ J ≈0.842 are the values at jamming onset. Δz and G for the SP and DP models begin to differ for ρ≥0.88. In this regime, Δz∼G can be described by a sum of two power-laws in Δρ, i.e. Δz∼G∼C 0 Δρ 0.5 +C 1 Δρ 1.0 to lowest order. We show that the ratio C 1 /C 0 is much larger for the DP model compared to that for the SP model. We also characterize the void space in jammed packings as a function of ρ. We find that the DP model can describe the formation of Plateau borders as ρ→1.0. We further show that the results for z and the shape factor A versus ρ for the DP model agree with recent experimental studies of foams and emulsions.more » « less
-
Physical kinetic roughening processes are well-known to exhibit universal scaling of observables that fluctuate in space and time. Are there analogous dynamic scaling laws that are unique to the chemical reaction mechanisms available synthetically and occurring naturally? Here, we formulate an approach to the dynamic scaling of stochastic fluctuations in thermodynamic observables at and away from equilibrium. Both analytical expressions and numerical simulations confirm our dynamic scaling ansatz with associated scaling exponents, function, and law. A survey of common chemical mechanisms reveals classes that organize according to the molecularity of the reactions involved, the nature of the reaction vessel and external reservoirs, (non)equilibrium conditions, and the extent of autocatalysis in the reaction network. Varying experimental parameters, such as temperature, can cause coupled reactions capable of chemical feedback to transition between these classes. While path observables, such as the dynamical activity, have scaling exponents that are time-independent, the variance in the entropy production and flow can have time-dependent scaling exponents and self-averaging properties as a result of temporal correlations that emerge during thermodynamically irreversible processes. Altogether, these results establish dynamic universality classes in the nonequilibrium fluctuations of thermodynamic observables for well-mixed chemical reactions.more » « less
-
Discontinuous shear thickening (DST) is associated with a sharp rise in a suspension’s viscosity with increasing applied shear rate or stress. Key signatures of DST, highlighted in recent studies, are the very large fluctuations of the measured stress as the suspension thickens with increasing rate. A clear link between microstructural development and the dramatic increase in stress fluctuations has not been established yet. To identify the microstructural underpinnings of this behavior, we perform simulations of sheared dense suspensions. Through an analysis of the particle contact network, we identify a subset of constrained particles that contributes directly to the rapid rise in viscosity and large stress fluctuations. Indeed, both phenomena can be explained by the growth and percolation of constrained particle networks—in direct analogy to rigidity percolation. A finite size scaling analysis confirms this to be a percolation phenomenon and allows us to estimate the critical exponents. Our findings reveal the specific microstructural self-organization transition that underlies DST.more » « less
-
An approach is presented for simulating multipulse nuclear magnetic resonance (NMR) spectra of polycrystalline solids directly in the frequency domain. The approach integrates the symmetry pathway concept for multipulse NMR with efficient algorithms for calculating spinning sideband amplitudes and performing interpolated finite-element numerical integration over all crystallite orientations in a polycrystalline sample. The numerical efficiency is achieved through a set of assumptions used to approximate the evolution of a sparse density matrix through a pulse sequence as a set of individual transition pathway signals. The utility of this approach for simulating the spectra of complex materials, such as glasses and other structurally disordered materials, is demonstrated.more » « less
An official website of the United States government

