Vibrations and diverging length scales near the unjamming transition

Jamming is a first-order transition with quenched disorder in amorphous materials sheared by cyclic quasistatic deformations

Jamming is an athermal transition between flowing and rigid states in amorphous systems such as granular matter, colloidal suspensions, complex fluids and cells. The jamming transition seems to display mixed aspects of a first-order transition, evidenced by a discontinuity in the coordination number, and a second-order transition, indicated by power-law scalings and diverging lengths. Here we demonstrate that jamming is a first-order transition with quenched disorder in cyclically sheared systems with quasistatic deformations, in two and three dimensions. Based on scaling analyses, we show that fluctuations of the jamming density in finite-sized systems have important consequences on the finite-size effects of various quantities, resulting in a square relationship between disconnected and connected susceptibilities, a key signature of the first-order transition with quenched disorder. This study puts the jamming transition into the category of a broad class of transitions in disordered systems where sample-to-sample fluctuations dominate over thermal fluctuations, suggesting that the nature and behavior of the jamming transition might be better understood within the developed theoretical framework of the athermally driven random-field Ising model.

Experimental measurements of the granular density of modes via impact

The jamming transition is an important feature of granular materials, with prior work showing an excess of low-frequency modes in the granular analog to the density of states, the granular density of modes. In this work, we present an experimental method for acoustically measuring the granular density of modes using a single impact event to excite vibrational modes in an experimental, three-dimensional, granular material. We test three different granular materials, all of which are composed of spherical beads. The first two systems are monodisperse collections of either 6 mm or 8 mm diameter beads. The third system is a bidisperse mixture of the previous two bead sizes. During data collection, the particles are confined to a box; on top of this box, and resting on the granular material, is a light, rigid sheet onto which pressure can be applied to the system. To excite the material, a steel impactor ball is dropped on top of the system. The response of the granular material to the impact pulse is recorded by piezoelectric sensors buried throughout the material, and the density of modes is computed from the spectrum of the velocity autocorrelation of these sensors. Our measurements of the density of modes show more low-frequency modes at low pressure, consistent with previous experimental and numerical results, as well as several low-frequency peaks in the density of modes that shift with applied pressure. Our method represents an experimentally simple technique for investigating the granular density of modes and may increase the accessibility and number of such measurements.

Long-range correlations in elastic moduli and local stresses at the unjamming transition

We explore the behaviour of spatially heterogeneous elastic moduli as well as the correlations between local moduli in model solids with short-range repulsive potentials. We show through numerical simulations that local elastic moduli exhibit long-range correlations, similar to correlations in the local stresses. Specifically, the correlations in local shear moduli exhibit anisotropic behavior at large lengthscales characterized by pinch-point singularities in Fourier space, displaying a structural pattern akin to shear stress correlations. Focussing on two-dimensional jammed solids approaching the unjamming transition, we show that stress correlations exhibit universal properties, characterized by a quadratic p2 dependence of the correlations as the pressure p approaches zero, independent of the details of the model. In contrast, the modulus correlations exhibit a power-law dependence with different exponents depending on the specific interaction potential. Furthermore, we illustrate that while affine responses lack long-range correlations, the total modulus, which encompasses non-affine behavior, exhibits long-range correlations.

Critical-like slowdown in thermal soft-sphere glasses via energy minimization

Using hybrid molecular dynamics/SWAP Monte Carlo (MD/SMC) simulations, we show that while the terminal relaxation times τ(ϕ) for FIRE energy minimization of soft-sphere glasses can decrease by orders of magnitude as sample equilibration proceeds and the jamming density ϕ_{J} increases, they always scale as τ(ϕ)∼(ϕ_{J}-ϕ)^{-2}∼[Z_{iso}-Z_{ms}(τ)]^{-2}, where Z_{iso}=2d and Z_{ms}(τ) is the average coordination number of particles satisfying a minimal local mechanical stability criterion (Z≥d+1) at the top of the final potential-energy-landscape (PEL) sub-basin the system encounters. This scaling allows us to collapse τ datasets that look very different when plotted as a function of ϕ, and to address a closely related question: how does the character of the PEL basins that dense thermal glasses most typically occupy evolve as the glasses age at constant ϕ and T?

An explicit granular-mechanics approach to marine sediment acoustics

Here, we theoretically and computationally study the frequency dependence of phase speed and attenuation for marine sediments from the perspective of granular mechanics. We leverage recent theoretical insights from the granular physics community as well as discrete-element method simulations, where the granular material is treated as a packing of discrete objects that interact via pairwise forces. These pairwise forces include both repulsive contact forces as well as dissipative terms, which may include losses from the fluid as well as losses from inelasticity at grain-grain contacts. We show that the structure of disordered granular packings leads to anomalous scaling laws for frequency-dependent phase speed and attenuation that do not follow from a continuum treatment. Our results demonstrate that granular packing structure, which is not explicitly considered in existing models, may play a crucial role in a complete theory of sediment acoustics. While this simple approach does not explicitly treat sound propagation or inertial effects in the interstitial fluid, it provides a starting point for future models that include these and other more complex features.

Effects of coordination and stiffness scale separation in disordered elastic networks

Many fibrous materials are modeled as elastic networks featuring a substantial separation between the stiffness scales that characterize different microscopic deformation modes of the network's constituents. This scale separation has been shown to give rise to emergent complexity in these systems' linear and nonlinear mechanical response. Here we study numerically a simple model featuring said stiffness scale separation in two-dimensions and show that its mechanical response is governed by the competition between the characteristic stiffness of collective nonphononic soft modes of the stiff subsystem, and the characteristic stiffness of the soft interactions. We present and rationalize the behavior of the shear modulus of our complex networks across the unjamming transition at which the stiff subsystem alone loses its macroscopic mechanical rigidity. We further establish a relation in the soft-interaction-dominated regime between the shear modulus, the characteristic frequency of nonphononic vibrational modes, and the mesoscopic correlation length that marks the crossover from a disorder-dominated response to local mechanical perturbations in the near field, to a linear, continuumlike response in the far field. The effects of spatial dimension on the observed scaling behavior are discussed, in addition to the interplay between stiffness scales in strain-stiffened networks, which is relevant to understanding the nonlinear mechanics of non-Brownian fibrous biomatter.

Connecting polymer collapse and the onset of jamming

Previous studies have shown that the interiors of proteins are densely packed, reaching packing fractions that are as large as those found for static packings of individual amino-acid-shaped particles. How can the interiors of proteins take on such high packing fractions given that amino acids are connected by peptide bonds and many amino acids are hydrophobic with attractive interactions? We investigate this question by comparing the structural and mechanical properties of collapsed attractive disk-shaped bead-spring polymers to those of three reference systems: static packings of repulsive disks, of attractive disks, and of repulsive disk-shaped bead-spring polymers. We show that the attractive systems quenched to temperatures below the glass transition T≪T_{g} and static packings of both repulsive disks and bead-spring polymers possess similar interior packing fractions. Previous studies have shown that static packings of repulsive disks are isostatic at jamming onset, i.e., the number of interparticle contacts N_{c} matches the number of degrees of freedom, which strongly influences their mechanical properties. We find that repulsive polymer packings are hypostatic at jamming onset (i.e., with fewer contacts than degrees of freedom) but are effectively isostatic when including stabilizing quartic modes, which give rise to quartic scaling of the potential energy with displacements along these modes. While attractive disk and polymer packings are often considered hyperstatic with excess contacts over the isostatic number, we identify a definition for interparticle contacts for which they can also be considered as effectively isostatic. As a result, we show that the mechanical properties (e.g., scaling of the potential energy with excess contact number and low-frequency contribution to the density of vibrational modes) of weakly attractive disk and polymer packings are similar to those of isostatic repulsive disk and polymer packings. Our results demonstrate that static packings generated via attractive collapse or compression of repulsive particles possess similar structural and mechanical properties.

Jamming on convex deformable surfaces

Jamming is a fundamental transition that governs the behavior of particulate media, including sand, foams and dense suspensions. Upon compression, such media change from freely flowing to a disordered, marginally stable solid that exhibits non-Hookean elasticity. While the jamming process is well established for fixed geometries, the nature and dynamics of jamming for a diverse class of soft materials and deformable substrates, including emulsions and biological matter, remains unknown. Here we propose a new scenario, metric jamming, where rigidification occurs on a surface that has been deformed from its ground state. Unlike classical jamming processes that exhibit discrete mechanical transitions, surprisingly we find that metric jammed states possess mechanical properties continuously tunable between those of classically jammed and conventional elastic media. The compact and curved geometry significantly alters the vibrational spectra of the structures relative to jamming in flat Euclidean space, and metric jammed systems also possess new types of vibrational mode that couple particle and shape degrees of freedom. Our work provides a theoretical framework that unifies our understanding of solidification processes that take place on deformable media and lays the groundwork to exploit jamming for the control and stabilization of shape in self-assembly processes.

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.

Shape-driven confluent rigidity transition in curved biological tissues

Collective cell motions underlie structure formation during embryonic development. Tissues exhibit emergent multicellular characteristics such as jamming, rigidity transitions, and glassy dynamics, but there remain questions about how those tissue-scale dynamics derive from local cell-level properties. Specifically, there has been little consideration of the interplay between local tissue geometry and cellular properties influencing larger-scale tissue behaviors. Here, we consider a simple two-dimensional computational vertex model for confluent tissue monolayers, which exhibits a rigidity phase transition controlled by the shape index (ratio of perimeter to square root area) of cells, on surfaces of constant curvature. We show that the critical point for the rigidity transition is a function of curvature such that positively curved systems are likely to be in a less rigid, more fluid, phase. Likewise, negatively curved systems (saddles) are likely to be in a more rigid, less fluid, phase. A phase diagram we generate for the curvature and shape index constitutes a testable prediction from the model. The curvature dependence is interesting because it suggests a natural explanation for more dynamic tissue remodeling and facile growth in regions of higher surface curvature. Conversely, we would predict stability at the base of saddle-shaped budding structures without invoking the need for biochemical or other physical differences. This concept has potential ramifications for our understanding of morphogenesis of budding and branching structures.

Emergence of rigidity percolation in flowing granular systems

Jammed granular media and glasses exhibit spatial long-range correlations as a result of mechanical equilibrium. However, the existence of such correlations in the flowing matter, where the mechanical equilibrium is unattainable, has remained elusive. Here, we investigate this problem in the context of the percolation of interparticle forces in flowing granular media. We find that the flow rate introduces an effective long-range correlation, which plays the role of a relevant perturbation giving rise to a spectrum of varying exponents on a critical line as a function of the flow rate. Our numerical simulations along with analytical arguments predict a crossover flow rate [Formula: see text] below which the effect of induced disorder is weak and the universality of the force chain structure is shown to be given by the standard rigidity percolation. We also find a power-law behavior for the critical exponents with the flow rate [Formula: see text].

Designing the pressure-dependent shear modulus using tessellated granular metamaterials

Jammed packings of granular materials display complex mechanical response. For example, the ensemble-averaged shear modulus 〈G〉 increases as a power law in pressure p for static packings of soft spherical particles that can rearrange during compression. We seek to design granular materials with shear moduli that can either increase or decrease with pressure without particle rearrangements even in the large-system limit. To do this, we construct tessellated granular metamaterials by joining multiple particle-filled cells together. We focus on cells that contain a small number of bidisperse disks in two dimensions. We first study the mechanical properties of individual disk-filled cells with three types of boundaries: periodic boundary conditions (PBC), fixed-length walls (FXW), and flexible walls (FLW). Hypostatic jammed packings are found for cells with FLW, but not in cells with PBC and FXW, and they are stabilized by quartic modes of the dynamical matrix. The shear modulus of a single cell depends linearly on p. We find that the slope of the shear modulus with pressure λ_{c}<0 for all packings in single cells with PBC where the number of particles per cell N≥6. In contrast, single cells with FXW and FLW can possess λ_{c}>0, as well as λ_{c}<0, for N≤16. We show that we can force the mechanical properties of multicell granular metamaterials to possess those of single cells by constraining the end points of the outer walls and enforcing an affine shear response. These studies demonstrate that tessellated granular metamaterials provide a platform for the design of soft materials with specified mechanical properties.

Microrheology near jamming

The jamming transition is a nonequilibrium critical phenomenon, which governs characteristic mechanical properties of jammed soft materials, such as pastes, emulsions, and granular matters. Both experiments and theory of jammed soft materials have revealed that the complex modulus measured by conventional macrorheology exhibits a characteristic frequency dependence. Microrheology is a new type of method to obtain the complex modulus, which transforms the microscopic motion of probes to the complex modulus through the generalized Stokes relation (GSR). Although microrheology has been applied to jammed soft materials, its theoretical understanding is limited. In particular, the validity of the GSR near the jamming transition is far from obvious since there is a diverging length scale lc, which characterizes the heterogeneous response of jammed particles. Here, we study the microrheology of jammed particles by theory and numerical simulation. First, we develop a linear response formalism to calculate the response function of the probe particle, which is transformed to the complex modulus via the GSR. Then, we apply our formalism to a numerical model of jammed particles and find that the storage and loss modulus follow characteristic scaling laws near the jamming transition. Importantly, the observed scaling law coincides with that in macrorheology, which indicates that the GSR holds even near the jamming transition. We rationalize this equivalence by asymptotic analysis of the obtained formalism and numerical analysis on the displacement field of jammed particles under a local perturbation.

Jamming and unusual charge density fluctuations of strange metals

The strange metallic regime across a number of high-temperature superconducting materials presents numerous challenges to the classic theory of Fermi liquid metals. Recent measurements of the dynamical charge response of strange metals, including optimally doped cuprates, have revealed a broad, featureless continuum of excitations, extending over much of the Brillouin zone. The collective density oscillations of this strange metal decay into the continuum in a manner that is at odds with the expectations of Fermi liquid theory. Inspired by these observations, we investigate the phenomenology of bosonic collective modes and the particle-hole excitations in a class of strange metals by making an analogy to the phonons of classical lattices falling apart across an unconventional jamming-like transition associated with the onset of rigidity. By making comparisons to the experimentally measured dynamical response functions, we reproduce many of the qualitative features using the above framework. We conjecture that the dynamics of electronic charge density over an intermediate range of energy scales in a class of strongly correlated metals can be at the brink of a jamming-like transition.

Finite amplitude waves in jammed matter

Here we use simulations and theory to show that, close to the jamming point, an arbitrary initial distortion of a granular media induces the formation of forward and backward non-linear finite amplitude waves. There are two regimes in the evolution of these waves (near field and far field). Initially, non-linear interactions between forward and backward waves dominate the propagation, leading to complex early evolution (near field). At longer times, forward and backwards waves cease interacting in the far field, and the propagation enters a new regime. Here the waves acquire a triangular-like profile, and evolve in a self-similar fashion characterized by a power law attenuation, whose exponent is weakly dependent on the initial pressure of the system. The finite amplitude waves gradually become linear waves when the amplitude of the initial distortion decreases, or the confining pressure on the system increases.

Anomalous linear elasticity of disordered networks

Continuum elasticity is a powerful tool applicable in a broad range of physical systems and phenomena. Yet, understanding how and on what scales material disorder may lead to the breakdown of continuum elasticity is not fully understood. We show, based on recent theoretical developments and extensive numerical computations, that disordered elastic networks near a critical rigidity transition, such as strain-stiffened fibrous biopolymer networks that are abundant in living systems, reveal an anomalous long-range linear elastic response below a correlation length. This emergent anomalous elasticity, which is non-affine in nature, is shown to feature a qualitatively different multipole expansion structure compared to ordinary continuum elasticity, and a slower spatial decay of perturbations. The potential degree of universality of these results, their implications (e.g. for cell-cell communication through biological extracellular matrices) and open questions are briefly discussed.

Universal scaling for disordered viscoelastic matter near the onset of rigidity

The onset of rigidity in interacting liquids, as they undergo a transition to a disordered solid, is associated with a rearrangement of the low-frequency vibrational spectrum. In this Letter, we derive scaling forms for the singular dynamical response of disordered viscoelastic networks near both jamming and rigidity percolation. Using effective-medium theory, we extract critical exponents, invariant scaling combinations, and analytical formulas for universal scaling functions near these transitions. Our scaling forms describe the behavior in space and time near the various onsets of rigidity, for rigid and floppy phases and the crossover region, including diverging length scales and timescales at the transitions.

An acoustic signature of extreme failure on model granular materials

Unexpectedly, granular materials can fail, the structure even destroyed, spontaneously in simple isotropic compression with stick-slip-like frictional behaviour. This extreme behaviour is conceptually impossible for saturated two-phase assembly in classical granular physics. Furthermore, the triggering mechanisms of these laboratory events remain mysterious, as in natural earthquakes. Here, we report a new interpretation of these failures in under-explored isotropic compression using the time-frequency analysis of Cauchy continuous wavelet transform of acoustic emissions and multiphysics numerical simulations. Wavelet transformation techniques can give insights into the temporal evolution of the state of granular materials en route to failure and offer a plausible explanation of the distinctive hearing sound of the stick-slip phenomenon. We also extend the traditional statistical seismic Gutenberg-Richter power-law behaviour for hypothetical biggest earthquakes based on the mechanisms of stick-slip frictional instability, using very large artificial isotropic labquakes and the ultimate unpredictable liquefaction failure.

Orientational-induced strain hardening of axisymmetric grains

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

Phonon transport properties of particulate physical gels

Particulate physical gels are sparse, low-density amorphous materials in which clusters of glasses are connected to form a heterogeneous network structure. This structure is characterized by two length scales, ξ s and ξ G: ξ s measures the length of heterogeneities in the network structure and ξ G is the size of glassy clusters. Accordingly, the vibrational states (eigenmodes) of such a material also exhibit a multiscale nature with two characteristic frequencies, [Formula: see text] and ω G, which are associated with ξ s and ξ G, respectively: (i) phonon-like vibrations in the homogeneous medium at [Formula: see text], (ii) phonon-like vibrations in the heterogeneous medium at [Formula: see text], and (iii) disordered vibrations in the glassy clusters at ω > ω G. Here, we demonstrate that the multiscale characteristics seen in the static structures and vibrational states also extend to the phonon transport properties. Phonon transport exhibits two distinct crossovers at frequencies ω* and ω G (or at wavenumbers of [Formula: see text] and [Formula: see text]). In particular, both transverse and longitudinal phonons cross over between Rayleigh scattering at [Formula: see text] and diffusive damping at [Formula: see text]. Remarkably, the Ioffe–Regel limit is located at the very low frequency of ω*. Thus, phonon transport is localized above ω*, even where phonon-like vibrational states persist. This markedly strong scattering behavior is caused by the sparse, porous structure of the gel.

Transient learning degrees of freedom for introducing function in materials

SignificanceMany protocols used in material design and training have a common theme: they introduce new degrees of freedom, often by relaxing away existing constraints, and then evolve these degrees of freedom based on a rule that leads the material to a desired state at which point these new degrees of freedom are frozen out. By creating a unifying framework for these protocols, we can now understand that some protocols work better than others because the choice of new degrees of freedom matters. For instance, introducing particle sizes as degrees of freedom to the minimization of a jammed particle packing can lead to a highly stable state, whereas particle stiffnesses do not have nearly the same impact.

Structural, mechanical, and vibrational properties of particulate physical gels

Our lives are surrounded by a rich assortment of disordered materials. In particular, glasses are well known as dense, amorphous materials, whereas gels exist in low-density, disordered states. Recent progress has provided a significant step forward in understanding the material properties of glasses, such as mechanical, vibrational, and transport properties. In contrast, our understanding of particulate physical gels is still highly limited. Here, using molecular dynamics simulations, we study a simple model of particulate physical gels, the Lennard-Jones (LJ) gels, and provide a comprehensive understanding of their structural, mechanical, and vibrational properties, all of which are markedly different from those of LJ glasses. First, the LJ gels show sparse, heterogeneous structures, and the length scale ξ_s of the structures grows as the density is lowered. Second, the LJ gels are extremely soft, with both shear G and bulk K moduli being orders of magnitude smaller than those of LJ glasses. Third, many low-frequency vibrational modes are excited, which form a characteristic plateau with the onset frequency ω_* in the vibrational density of states. Structural, mechanical, and vibrational properties, characterized by ξ_s, G, K, and ω_*, respectively, show power-law scaling behaviors with the density, which establishes a close relationship between them. Throughout this work, we also reveal that LJ gels are multiscale, solid-state materials: (i) homogeneous elastic bodies at long lengths, (ii) heterogeneous elastic bodies with fractal structures at intermediate lengths, and (iii) amorphous structural bodies at short lengths.

Achieving adjustable elasticity with non-affine to affine transition

For various engineering and industrial applications it is desirable to realize mechanical systems with broadly adjustable elasticity to respond flexibly to the external environment. Here we discover a topology-correlated transition between affine and non-affine regimes in elasticity in both two- and three-dimensional packing-derived networks. Based on this transition, we numerically design and experimentally realize multifunctional systems with adjustable elasticity. Within one system, we achieve solid-like affine response, liquid-like non-affine response and a continuous tunability in between. Moreover, the system also exhibits a broadly tunable Poisson's ratio from positive to negative values, which is of practical interest for energy absorption and for fracture-resistant materials. Our study reveals a fundamental connection between elasticity and network topology, and demonstrates its practical potential for designing mechanical systems and metamaterials.

Low-energy quasilocalized excitations in structural glasses

Glassy solids exhibit a wide variety of generic thermomechanical properties, ranging from universal anomalous specific heat at cryogenic temperatures to nonlinear plastic yielding and failure under external driving forces, which qualitatively differ from their crystalline counterparts. For a long time, it has been believed that many of these properties are intimately related to nonphononic, low-energy quasilocalized excitations (QLEs) in glasses. Indeed, recent computer simulations have conclusively revealed that the self-organization of glasses during vitrification upon cooling from a melt leads to the emergence of such QLEs. In this Perspective, we review developments over the past three decades toward understanding the emergence of QLEs in structural glasses and the degree of universality in their statistical and structural properties. We discuss the challenges and difficulties that hindered progress in achieving these goals and review the frameworks put forward to overcome them. We conclude with an outlook on future research directions and open questions.

Unified quantifier of mechanical disorder in solids

Mechanical disorder in solids, which is generated by a broad range of physical processes and controls various material properties, appears in a wide variety of forms. Defining unified and measurable dimensionless quantifiers, allowing quantitative comparison of mechanical disorder across widely different physical systems, is therefore an important goal. Two such coarse-grained dimensionless quantifiers (among others) appear in the literature: one is related to the spectral broadening of discrete phononic bands in finite-size systems (accessible through computer simulations) and the other is related to the spatial fluctuations of the shear modulus in macroscopically large systems. The latter has been recently shown to determine the amplitude of wave attenuation rates in the low-frequency limit (accessible through laboratory experiments). Here, using two alternative and complementary theoretical approaches linked to the vibrational spectra of solids, we derive a basic scaling relation between the two dimensionless quantifiers. This scaling relation, which is supported by simulational data, shows that the two apparently distinct quantifiers are in fact intrinsically related, giving rise to a unified quantifier of mechanical disorder in solids. We further discuss the obtained results in the context of the unjamming transition taking place in soft sphere packings at low confining pressures, in addition to their implications for our understanding of the low-frequency vibrational spectra of disordered solids in general, and in particular those of glassy systems.

Structural and short-time vibrational properties of colloidal glasses and supercooled liquids in the vicinity of the re-entrant glass transition

We investigate the short-time vibrational properties and structure of two-dimensional, bidisperse, colloidal glasses and supercooled liquids in the vicinity of the re-entrant glass transition, as a function of interparticle depletion attraction strength. The long-time spatiotemporal dynamics of the samples are measured to be non-monotonic, confirming that the suspensions evolve from repulsive glass to supercooled liquid to attractive glass with increasing depletion attraction. Here, we search for vibrational signatures of the re-entrant behavior in the short-time spatiotemporal dynamics, i.e., dynamics associated with particle motion inside its nearest-neighbor cage. Interestingly, we observe that the anharmonicity of these in-cage vibrations varies non-monotonically with increasing attraction strength, consistent with the non-monotonic long-time structural relaxation dynamics of the re-entrant glass. We also extract effective spring constants between neighboring particles; we find that spring stiffness involving small particles also varies non-monotonically with increasing attraction strength, while stiffness between large particles increases monotonically. Last, from study of depletion-dependent local structure and vibration participation fractions, we gain microscopic insight into the particle-size-dependent contributions to short-time vibrational modes in the glass and supercooled liquid states.

Nonlocal Effects Reflect the Jamming Criticality in Frictionless Granular Flows Down Inclines

The jamming transition is accompanied by a rich phenomenology such as hysteresis or nonlocal effects that is still not well understood. Here, we experimentally investigate a model frictionless granular layer flowing down an inclined plane as a way to disentangle generic collective effects from those arising from frictional interactions. We find that thin frictionless granular layers are devoid of hysteresis of the avalanche angle, yet the layer stability increases as it gets thinner. Steady rheological laws obtained for different layer thicknesses can be collapsed into a unique master curve, supporting the idea that nonlocal effects are the consequence of the usual finite-size effects associated with the presence of a critical point. This collapse indicates that the so-called isostatic length l^{*}, the scale on which pinning a boundary freezes all remaining floppy modes, governs the effect of boundaries on flow and rules out other propositions made in the past.

Sound damping in soft particle packings: the interplay between configurational disorder and inelasticity

We numerically investigate sound damping in disordered two-dimensional soft particle packings. We simulate evolution of standing waves of particle displacements and analyze time correlation functions of particle velocities and power spectra. We control the strength of inelastic interactions between the particles in contact to show how the inelasticity affects anomalous sound characteristics of disordered systems: Increasing the strength of inelastic interactions, we find that (i) sound softening vanishes and (ii) attenuation coefficients exhibit a transition from the Rayleigh law to quadratic growth. We also report (iii) how the Ioffe-Regel limit frequencies depend on the strength of inelasticity as useful information for experiments and applications of the sound in disordered media. Our findings suggest that sound damping in soft particle packings is determined by the interplay between elastic heterogeneities and inelasticity.

Spongelike Rigid Structures in Frictional Granular Packings

We show how rigidity emerges in experiments on sheared two-dimensional frictional granular materials by using generalizations of two methods for identifying rigid structures. Both approaches, the force-based dynamical matrix and the topology-based rigidity percolation, agree with each other and identify similar rigid structures. As the system becomes jammed, at a critical contact number z_{c}=2.4±0.1, a rigid backbone interspersed with floppy, particle-filled holes of a broad range of sizes emerges, creating a spongelike morphology. While the pressure within rigid structures always exceeds the pressure outside the rigid structures, they are not identified with the force chains of shear jamming. These findings highlight the need to focus on mechanical stability arising through arch structures and hinges at the mesoscale.

A structural approach to vibrational properties ranging from crystals to disordered systems

Many scientists generally attribute the vibrational anomalies of disordered solids to the structural disorder, which, however, is still under intense debate. Here we conduct simulations on two-dimensional packings with a finite temperature, whose structure is tuned from a crystalline configuration to an amorphous one, then the amorphous from very dense state to a relatively loose state. By measuring the vibrational density of states and the reduced density of states, we clearly observe the evolution of the boson peak with the change of the disorder and volume fractions. Meanwhile, to understand the structural origin of this anomaly, we identify the soft regimes of all systems with a novel machine-learning method, where the "softness", a local structural quantity, is defined. Interestingly, we find a strong monotonic relationship between the shape of the boson peak and the softness as well as its spatial heterogeneity, suggesting that the softness of a system may be a new structural approach to the anomalous vibrational properties of amorphous solids.

Mechanical disorder of sticky-sphere glasses. I. Effect of attractive interactions

Recent literature indicates that attractive interactions between particles of a dense liquid play a secondary role in determining its bulk mechanical properties. Here we show that, in contrast with their apparent unimportance to the bulk mechanics of dense liquids, attractive interactions can have a major effect on macro- and microscopic elastic properties of glassy solids. We study several broadly applicable dimensionless measures of stability and mechanical disorder in simple computer glasses, in which the relative strength of attractive interactions-referred to as "glass stickiness"-can be readily tuned. We show that increasing glass stickiness can result in the decrease of various quantifiers of mechanical disorder, on both macro- and microscopic scales, with a pair of intriguing exceptions to this rule. Interestingly, in some cases strong attractions can lead to a reduction of the number density of soft, quasilocalized modes, by up to an order of magnitude, and to a substantial decrease in their core size, similar to the effects of thermal annealing on elasticity observed in recent works. Contrary to the behavior of canonical glass models, we provide compelling evidence indicating that the stabilization mechanism in our sticky-sphere glasses stems predominantly from the self-organized depletion of interactions featuring large, negative stiffnesses. Finally, we establish a fundamental link between macroscopic and microscopic quantifiers of mechanical disorder, which we motivate via scaling arguments. Future research directions are discussed.

Size distribution dependence of collective relaxation dynamics in a two-dimensional wet foam

Foams can be ubiquitously observed in nature and in industrial products. Despite the relevance of their properties to deformation, fluidity, and collapse, all of which are essential for applications, there are few experimental studies of collective relaxation dynamics in a wet foam. Here, we directly observe how the relaxation dynamics changes with increasing liquid fraction in both monodisperse and polydisperse two-dimensional foams. As we increase the liquid fraction, we quantitatively characterize the slowing-down of the relaxation, and the increase of the correlation length. We also find two different relaxation modes which depend on the size distribution of the bubbles. It suggests that the bubbles which are simply near to each other play an important role in large rearrangements, not just those in direct contact. Finally, we confirm the generality of our experimental findings by a numerical simulation for the relaxation process of wet foams.

Role of Anharmonic Interactions for Vibration Density of States in α-Cristobalite

The vibrational density of states (VDOS) of solids in the low-energy regime controls the thermal and transport properties of materials, such as heat capacity, heat conduction, free energy and entropy. In α-Cristobalite, the low-frequency part of vibration density of states (VDOS) has many common features with the Boson peak in silica glass of matched densities. Recent theoretical work reported that anharmonic phonon–phonon interactions were critical for the low-frequency part of VDOS in α-Cristobalite. Therefore, it is urgent to identify the role of different anharmonic interactions from first principles. In this paper, we focus on the main peak of the low-frequency part of VDOS in α-Cristobalite. Calculated by our own developed codes and first principles, we find that the quartic anharmonic interaction can increase the frequency of the peak, while the cubic anharmonic can reduce the frequency and change the shape of the peak. Meanwhile, the anharmonic interactions are critical for the temperature effect. Therefore, we calculated the temperature-dependent property of the peak. We find that the frequency of the peak is directly proportional to the temperature. The atomic displacement patterns of different temperatures also confirm the above conclusion. All our calculations converged well. Moreover, our basic results agree well with other published results. Finally, we highlight that our codes offer a general and reliable way to calculate the VDOS with temperature.

Novel elastic instability of amorphous solids in finite spatial dimensions

Recently, progress has been made in the understanding of anomalous vibrational excitations in amorphous solids. In the lowest-frequency region, the vibrational spectrum follows a non-Debye quartic law, which persists up to zero frequency without any frequency gap. This gapless vibrational density of states (vDOS) suggests that glasses are on the verge of instability. This feature of marginal stability is now highlighted as a key concept in the theories of glasses. In particular, the elasticity theory based on marginal stability predicts the gapless vDOS. However, this theory yields a quadratic law and not the quartic law. To address this inconsistency, we presented a new type of instability, which is different from the conventional one, and proposed that amorphous solids are marginally stable considering the new instability in the preceding study [M. Shimada, H. Mizuno and A. Ikeda, Soft Matter, 2020, 16, 7279]. In this study, we further extend and detail the results for these instabilities. By analyzing various examples of disorder, we demonstrate that real glasses in finite spatial dimensions can be marginally stable by the proposed novel instability.

Critical point of jamming transition in two-dimensional monodisperse systems

The existence of amorphous packings in two-dimensional monodisperse system is a classical unsolved problem. We get the energy minimum state by the energy minimization method of enthalpy under constant pressure conditions. Firstly, we find that there are two peaks in the experiment, which demonstrate the interesting features of the coexistence of crystals and amorphous crystals. And then, we confirm the critical point of jamming transition of the two-dimensional monodisperse is [Formula: see text]. Finally, we prove that the jamming scaling is still satisfied in two-dimensional monodispersed system: [Formula: see text] and vanishes as [Formula: see text], and the boson peak shifts to lower frequencies for less compressed systems.

Vibrational density of states of amorphous solids with long-ranged power-law-correlated disorder in elasticity

A theory of vibrational excitations based on power-law spatial correlations in the elastic constants (or equivalently in the internal stress) is derived, in order to determine the vibrational density of states D([Formula: see text]) of disordered solids. The results provide the first prediction of a boson peak in amorphous materials where spatial correlations in the internal stresses (or elastic constants) are of power-law form, as is often the case in experimental systems, leading to a logarithmic enhancement of (Rayleigh) phonon attenuation. A logarithmic correction of the form [Formula: see text] is predicted to occur in the plot of the reduced excess DOS for frequencies around the boson peak in 3D. Moreover, the theory provides scaling laws of the density of states in the low-frequency region, including a [Formula: see text] regime in 3D, and provides information about how the boson peak intensity depends on the strength of power-law decay of fluctuations in elastic constants or internal stress. Analytical expressions are also derived for the dynamic structure factor for longitudinal excitations, which include a logarithmic correction factor, and numerical calculations are presented supporting the assumptions used in the theory.

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

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

Propagation and attenuation of mechanical signals in ultrasoft 2D solids

The propagation of elastic waves in soft materials plays a crucial role in the spatiotemporal transmission of mechanical signals, e.g., in biological mechanotransduction or in the failure of marginal solids. At high Reynolds numbers Re ≫ 1, inertia dominates and wave propagation is readily observed. However, mechanical cues in soft and biological materials often occur at low Re, where waves are overdamped. Overdamped waves are not only difficult to observe experimentally, also theoretically their description remains incomplete. Here, we present direct measurements of the propagation and attenuation of mechanical signals in colloidal soft solids, induced by an optical trap. We derive an analytical theory for low Re wave propagation and damping, which is in excellent agreement with the experiments. Our results present both a previously unexplored method to characterize damped waves in soft solids and a theoretical framework showing how localized mechanical signals can provoke a remote and delayed response.

Power-Law Scaling in Granular Rheology across Flow Geometries

Based on discrete element method simulations, we propose a new form of the constitutive equation for granular flows independent of packing fraction. Rescaling the stress ratio μ by a power of dimensionless temperature Θ makes the data from a wide set of flow geometries collapse to a master curve depending only on the inertial number I. The basic power-law structure appears robust to varying particle properties (e.g., surface friction) in both 2D and 3D systems. We show how this rheology fits and extends frameworks such as kinetic theory and the nonlocal granular fluidity model.

Two-Scale Scenario of Rigidity Percolation of Sticky Particles

In the presence of attraction, the jamming transition of packings of frictionless particles corresponds to the rigidity percolation. When the range of attraction is long, the distribution of the size of rigid clusters, P(s), is continuous and shows a power-law decay. For systems with short-range attractions, however, P(s) appears discontinuous. There is a power-law decay for small cluster sizes, followed by a low probability gap and a peak near the system size. We find that this appearing "discontinuity" does not mean that the transition is discontinuous. In fact, it signifies the coexistence of two distinct length scales, associated with the largest cluster and smaller ones, respectively. The comparison between the largest and second largest clusters indicates that their growth rates with system size are rather different. However, both cluster sizes tend to diverge in the large system size limit, suggesting that the jamming transition of systems with short-range attractions is still continuous. In the framework of the two-scale scenario, we also derive a generalized hyperscaling relation. With robust evidence, our work challenges the former single-scale view of the rigidity percolation.

Vibrational Properties of Hard and Soft Spheres Are Unified at Jamming

The unconventional thermal properties of jammed amorphous solids are directly related to their density of vibrational states. While the vibrational spectrum of jammed soft sphere solids has been fully described, the vibrational spectrum of hard spheres, a model glass former often related to physical colloidal glasses, is still unknown due to the difficulty of treating the nonanalytic interaction potential. We bypass this difficulty using the recently described effective interaction potential for the free energy of thermal hard spheres. By minimizing this effective free energy, we mimic the rapid compression of hard spheres and produce typical configurations of the thermal system. We measure the resulting vibrational spectrum and characterize its evolution toward the jamming point where configurations of hard and soft spheres are trivially unified. For densities approaching jamming from below, we observe low-frequency modes which agree with those found in numerical simulations of jammed soft spheres. Our measurements of the vibrational structure demonstrate that the jamming universality extends away from jamming: hard sphere thermal systems below jamming exhibit the same vibrational spectra as thermal and athermal soft sphere systems above the transition.

Structured randomness: jamming of soft discs and pins

Simulations are used to find the zero temperature jamming threshold, φj, for soft, bidisperse disks in the presence of small fixed particles, or "pins", arranged in a lattice. The presence of pins leads, as one expects, to a decrease in φj. Structural properties of the system near the jamming threshold are calculated as a function of the pin density. While the correlation length exponent remains ν = 1/2 at low pin densities, the system is mechanically stable with more bonds, yet fewer contacts than the Maxwell criterion implies in the absence of pins. In addition, as pin density increases, novel bond orientational order and long-range spatial order appear, which are correlated with the square symmetry of the pin lattice.

Low-frequency vibrations of jammed packings in large spatial dimensions

Amorphous packings prepared in the vicinity of the jamming transition play a central role in theoretical studies of the vibrational spectrum of glasses. Two mean-field theories predict that the vibrational density of states g(ω) obeys a characteristic power law, g(ω)∼ω^{2}, called the non-Debye scaling in the low-frequency region. Numerical studies have, however, reported that this scaling breaks down at low frequencies, due to finite-dimensional effects. In this study, we prepare amorphous packings of up to 128000 particles in spatial dimensions from d=3 to d=9 to characterize the range of validity of the non-Debye scaling. Our numerical results suggest that the non-Debye scaling is obeyed down to a frequency that gradually decreases as d increases, and possibly vanishes for large d, in agreement with mean-field predictions. We also show that the prestress is an efficient control parameter to quantitatively compare packings across different spatial dimensions.

Stress Relaxation above and below the Jamming Transition

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

Universal Relaxation Dynamics of Sphere Packings below Jamming

We show that non-Brownian suspensions of repulsive spheres below jamming display a slow relaxational dynamics with a characteristic timescale that diverges at jamming. This slow timescale is fully encoded in the structure of the unjammed packing and can be readily measured via the vibrational density of states. We show that the corresponding dynamic critical exponent is the same for randomly generated and sheared packings. Our results show that a wide variety of physical situations, from suspension rheology to algorithmic studies of the jamming transition are controlled by a unique diverging timescale, with a universal critical exponent.

Boson peak, elasticity, and glass transition temperature in polymer glasses: Effects of the rigidity of chain bending

The excess low-frequency vibrational spectrum, called boson peak, and non-affine elastic response are the most important particularities of glasses. Herein, the vibrational and mechanical properties of polymeric glasses are examined by using coarse-grained molecular dynamics simulations, with particular attention to the effects of the bending rigidity of the polymer chains. As the rigidity increases, the system undergoes a glass transition at a higher temperature (under a constant pressure), which decreases the density of the glass phase. The elastic moduli, which are controlled by the decrease of the density and the increase of the rigidity, show a non-monotonic dependence on the rigidity of the polymer chain that arises from the non-affine component. Moreover, a clear boson peak is observed in the vibrational density of states, which depends on the macroscopic shear modulus G. In particular, the boson peak frequency ω_BP is proportional to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{G}$$\end{document}G. These results provide a positive correlation between the boson peak, shear elasticity, and the glass transition temperature.

Jammed packings of 3D superellipsoids with tunable packing fraction, coordination number, and ordering

We carry out numerical studies of static packings of frictionless superellipsoidal particles in three spatial dimensions. We consider more than 200 different particle shapes by varying the three shape parameters that define superellipsoids. We characterize the structural and mechanical properties of both disordered and ordered packings using two packing-generation protocols. We perform athermal quasi-static compression simulations starting from either random, dilute configurations (Protocol 1) or thermalized, dense configurations (Protocol 2), which allows us to tune the orientational order of the packings. In general, we find that superellipsoid packings are hypostatic, with coordination number zJ < ziso, where ziso = 2df and df = 5 or 6 depending on whether the particles are axi-symmetric or not. Over the full range of orientational order, we find that the number of quartic modes of the dynamical matrix for the packings always matches the number of missing contacts relative to the isostatic value. This result suggests that there are no mechanically redundant contacts for ordered, yet hypostatic packings of superellipsoidal particles. Additionally, we find that the packing fraction at jamming onset for disordered packings of superellipsoidal depends on at least two particle shape parameters, e.g. the asphericity A and reduced aspect ratio β of the particles.

Vibrational density of states and specific heat in glasses from random matrix theory

The low-temperature properties of glasses present important differences with respect to crystalline matter. In particular, models such as the Debye model of solids, which assume the existence of an underlying regular lattice, predict that the specific heat of solids varies with the cube of temperature at low temperatures. Since the 1970s at least, it is a well-established experimental fact that the specific heat of glasses is instead just linear in T at T∼1 K and presents a pronounced peak when normalized by T^{3}, known as the boson peak. Here we present an approach which suggests that the vibrational and thermal properties of amorphous solids are affected by the random-matrix part of the vibrational spectrum. The model is also able to reproduce, for the first time, the experimentally observed inverse proportionality between the boson peak in the specific heat and the shear modulus.

Vibrational properties of two-dimensional dimer packings near the jamming transition

Jammed particulate systems composed of various shapes of particles undergo the jamming transition as they are compressed or decompressed. To date, sphere packings have been extensively studied in many previous works, where isostaticity at the transition and scaling laws with the pressure of various quantities, including the contact number and the vibrational density of states, have been established. Additionally, much attention has been paid to nonspherical packings, and particularly recent work has made progress in understanding ellipsoidal packings. In this work, we study the dimer packings in two dimensions, which have been much less understood than systems of spheres and ellipsoids. We first study the contact number of dimers near the jamming transition. It turns out that packings of dimers have "rotational rattlers," each of which still has a free rotational motion. After correcting this effect, we show that dimers become isostatic at the jamming, and the excess contact number obeys the same critical law and finite-size scaling law as those of spheres. We next study the vibrational properties of dimers near the transition. We find that the vibrational density of states of dimers exhibits two characteristic plateaus that are separated by a peak. The high-frequency plateau is dominated by the translational degree of freedom, while the low-frequency plateau is dominated by the rotational degree of freedom. We establish the critical scaling laws of the characteristic frequencies of the plateaus and the peak near the transition. In addition, we present detailed characterizations of the real space displacement fields of vibrational modes in the translational and rotational plateaus.