Vibrations and diverging length scales near the unjamming transition

Elastic moduli and vibrational modes in jammed particulate packings

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

Spatial structure of states of self stress in jammed systems

States of self stress, organizations of internal forces in many-body systems that are in equilibrium with an absence of external forces, can be thought of as the constitutive building blocks of the elastic response of a material. In overconstrained disordered packings they have a natural mathematical correspondence with the zero-energy vibrational modes in underconstrained systems. While substantial attention in the literature has been paid to diverging length scales associated with zero- and finite-energy vibrational modes in jammed systems, less is known about the spatial structure of the states of self stress. In this work we define a natural way in which a unique state of self stress can be associated with each bond in a disordered spring network derived from a jammed packing, and then investigate the spatial structure of these bond-localized states of self stress. This allows for an understanding of how the elastic properties of a system would change upon changing the strength or even existence of any bond in the system.

Heterogeneity in the Small-Scale Deformation Behavior of Disordered Nanoparticle Packings

Atomic force microscopy-based nanoindentation is used to image and probe the local mechanical properties of thin disordered nanoparticle packings. The probed region is limited to the size of a few particles, and an individual particle can be loaded and displaced to a fraction of a single particle radius. The results demonstrate heterogeneous mechanical response that is location-dependent. The weak locations may be analogous to the "soft spots" previously predicted in glasses and other disordered packings.

Grains unchained: local fluidization of a granular packing by focused ultrasound

We report experimental results on the dynamics of a granular packing submitted to high-intensity focused ultrasound. Acoustic radiation pressure is shown to remotely induce local rearrangements within a pile as well as global motion around the focal spot in an initially jammed system. We demonstrate that this fluidization process is intermittent for a range of acoustic pressures and hysteretic when the pressure is cycled. Such a first-order-like unjamming transition is reproduced in numerical simulations in which the acoustic pressure field is modeled by a localized external force. Further analysis of the simulated packings suggests that in the intermittent regime unjamming is not associated with any noticeable prior structural signature. A simple two-state model based on effective temperatures is proposed to account for these findings.

Length scales in glass-forming liquids and related systems: a review

The central problem in the study of glass-forming liquids and other glassy systems is the understanding of the complex structural relaxation and rapid growth of relaxation times seen on approaching the glass transition. A central conceptual question is whether one can identify one or more growing length scale(s) associated with this behavior. Given the diversity of molecular glass-formers and a vast body of experimental, computational and theoretical work addressing glassy behavior, a number of ideas and observations pertaining to growing length scales have been presented over the past few decades, but there is as yet no consensus view on this question. In this review, we will summarize the salient results and the state of our understanding of length scales associated with dynamical slow down. After a review of slow dynamics and the glass transition, pertinent theories of the glass transition will be summarized and a survey of ideas relating to length scales in glassy systems will be presented. A number of studies have focused on the emergence of preferred packing arrangements and discussed their role in glassy dynamics. More recently, a central object of attention has been the study of spatially correlated, heterogeneous dynamics and the associated length scale, studied in computer simulations and theoretical analysis such as inhomogeneous mode coupling theory. A number of static length scales have been proposed and studied recently, such as the mosaic length scale discussed in the random first-order transition theory and the related point-to-set correlation length. We will discuss these, elaborating on key results, along with a critical appraisal of the state of the art. Finally we will discuss length scales in driven soft matter, granular fluids and amorphous solids, and give a brief description of length scales in aging systems. Possible relations of these length scales with those in glass-forming liquids will be discussed.

Universal spectrum of normal modes in low-temperature glasses

We report an analytical study of the vibrational spectrum of the simplest model of jamming, the soft perceptron. We identify two distinct classes of soft modes. The first kind of modes are related to isostaticity and appear only in the close vicinity of the jamming transition. The second kind of modes instead are present everywhere in the glass phase and are related to the hierarchical structure of the potential energy landscape. Our results highlight the universality of the spectrum of normal modes in disordered systems, and open the way toward a detailed analytical understanding of the vibrational spectrum of low-temperature glasses.

Fluctuations in flows near jamming

Bubbles, droplets or particles in flowing complex media such as foams, emulsions or suspensions follow highly complex paths, with the relative motion of the constituents setting the energy dissipation rate. What is their dynamics, and how is this connected to the global rheology? To address these questions, we probe the statistics and spatio-temporal organization of the local particle motion and energy dissipation in a model for sheared disordered materials. We find that the fluctuations in the local dissipation vary from nearly Gaussian and homogeneous at low densities and fast flows, to strongly intermittent for large densities and slow flows. The higher order moments of the relative particle velocities reveal strong evidence for a qualitative difference between two distinct regimes which are nevertheless connected by a smooth crossover. In the critical regime, the higher order moments are related by novel multiscaling relations. In the plastic regime the relations between these moments take on a different form, with higher moments diverging rapidly when the flow rate vanishes. As these velocity differences govern the energy dissipation, we can distinguish two qualitatively different types of flow: an intermediate density, critical regime related to jamming, and a large density, plastic regime.

Diffusion in Jammed Particle Packs

Using random walk simulations we explore diffusive transport through monodisperse sphere packings over a range of packing fractions ϕ in the vicinity of the jamming transition at ϕ(c). Various diffusion properties are computed over several orders of magnitude in both time and packing pressure. Two well-separated regimes of normal "Fickian" diffusion, where the mean squared displacement is linear in time, are observed. The first corresponds to diffusion inside individual spheres, while the latter is the long-time bulk diffusion. The intermediate anomalous diffusion regime and the long-time value of the diffusion coefficient are both shown to be controlled by particle contacts, which in turn depend on proximity to ϕ(c). The time required to recover normal diffusion t* scales as (ϕ-ϕ(c))(-0.5) and the long-time diffusivity D(∞)∼(ϕ-ϕ(c))0.5, or D(∞)∼1/t*. It is shown that the distribution of mean first passage times associated with the escape of random walkers between neighboring particles controls both t* and D(∞) in the limit ϕ→ϕ(c).

Strain fluctuations and elastic moduli in disordered solids

Recently there has been a surge in interest in using video-microscopy techniques to infer the local mechanical properties of disordered solids. One common approach is to minimize the difference between particle vibrational displacements in a local coarse-graining volume and the displacements that would result from a best-fit affine deformation. Effective moduli are then inferred under the assumption that the components of this best-fit affine deformation tensor have a Boltzmann distribution. In this paper, we combine theoretical arguments with experimental and simulation data to demonstrate that the above does not reveal information about the true elastic moduli of jammed packings and colloidal glasses.

Elasticity of frictionless particles near jamming

We study the linear elastic response of harmonic disk packings near jamming via three types of probes: (i) point forcing, (ii) constrained homogeneous deformation of subregions of large systems, and (iii) unconstrained deformation of the full system subject to periodic boundary conditions. For the point forcing, our results indicate that the transverse component of the response is governed by a lengthscale ξT, which scales with the confining pressure, p, as ξT∼p-0.25, while the longitudinal component is governed by ξL, which scales as ξL∼p-0.4. The former scaling is precisely the transverse lengthscale, which has been invoked to explain the structure of normal modes near the density of states anomaly in sphere packings, while the latter is much closer to the rigidity length, l*∼p-0.5, which has been invoked to describe the jamming scenario. For the case of constrained homogeneous deformation, we find that μ(R), the value of the shear modulus measured in boxes of size R, gives a value much higher than the continuum result for small boxes and recedes to its continuum limit only for boxes bigger than a characteristic length, which scales like p-0.5, precisely the same way as l*. Finally, for the case of unconstrained homogeneous deformation, we find displacement fields with power spectra, which are consistent with independent, uncorrelated Eshelby transformations. The transverse sector is amazingly invariant with respect to p and very similar to what is seen in Lennard-Jones glasses. The longitudinal piece, however, is sensitive to p. It develops a plateau at long wavelength, the start of which occurs at a length that grows in the p→0 limit. Strikingly, the same behavior is observed both for applied shear and dilation.

Softening of stressed granular packings with resonant sound waves

We perform numerical simulations of a two-dimensional bidisperse granular packing subjected to both a static confining pressure and a sinusoidal dynamic forcing applied by a wall on one edge of the packing. We measure the response experienced by a wall on the opposite edge of the packing and obtain the resonant frequency of the packing as the static or dynamic pressures are varied. Under increasing static pressure, the resonant frequency increases, indicating a velocity increase of elastic waves propagating through the packing. In contrast, when the dynamic amplitude is increased for fixed static pressure, the resonant frequency decreases, indicating a decrease in the wave velocity. This occurs both for compressional and for shear dynamic forcing and is in agreement with experimental results. We find that the average contact number Zc at the resonant frequency decreases with increasing dynamic amplitude, indicating that the elastic softening of the packing is associated with a reduced number of grain-grain contacts through which the elastic waves can travel. We image the excitations created in the packing and show that there are localized disturbances or soft spots that become more prevalent with increasing dynamic amplitude. Our results are in agreement with experiments on glass bead packings and earth materials such as sandstone and granite and may be relevant to the decrease in elastic wave velocities that has been observed to occur near fault zones after strong earthquakes, in surficial sediments during strong ground motion, and in structures during earthquake excitation.

Phonons and elasticity in critically coordinated lattices

Much of our understanding of vibrational excitations and elasticity is based upon analysis of frames consisting of sites connected by bonds occupied by central-force springs, the stability of which depends on the average number of neighbors per site z. When z  <  zc  ≈  2d, where d is the spatial dimension, frames are unstable with respect to internal deformations. This pedagogical review focuses on the properties of frames with z at or near zc, which model systems like randomly packed spheres near jamming and network glasses. Using an index theorem, N0  -NS  =  dN  -NB relating the number of sites, N, and number of bonds, NB, to the number, N0, of modes of zero energy and the number, NS, of states of self stress, in which springs can be under positive or negative tension while forces on sites remain zero, it explores the properties of periodic square, kagome, and related lattices for which z  =  zc and the relation between states of self stress and zero modes in periodic lattices to the surface zero modes of finite free lattices (with free boundary conditions). It shows how modifications to the periodic kagome lattice can eliminate all but trivial translational zero modes and create topologically distinct classes, analogous to those of topological insulators, with protected zero modes at free boundaries and at interfaces between different topological classes.

Thermal fluctuations, mechanical response, and hyperuniformity in jammed solids

Jamming is a geometric phase transition occurring in dense particle systems in the absence of temperature. We use computer simulations to analyze the effect of thermal fluctuations on several signatures of the transition. We show that scaling laws for bulk and shear moduli only become relevant when thermal fluctuations are extremely small, and propose their relative ratio as a quantitative signature of jamming criticality. Despite the nonequilibrium nature of the transition, we find that thermally induced fluctuations and mechanical responses obey equilibrium fluctuation-dissipation relations near jamming, provided the appropriate fluctuating component of the particle displacements is analyzed. This shows that mechanical moduli can be directly measured from particle positions in mechanically unperturbed packings, and suggests that the definition of a "nonequilibrium index" is unnecessary for amorphous materials. We find that fluctuations of particle displacements are spatially correlated, and define a transverse and a longitudinal correlation length scale which both diverge as the jamming transition is approached. We analyze the frozen component of density fluctuations and find that it displays signatures of nearly hyperuniform behavior at large length scales. This demonstrates that hyperuniformity in jammed packings is unrelated to a vanishing compressibility and explains why it appears remarkably robust against temperature and density variations. Differently from jamming criticality, obstacles preventing the observation of hyperuniformity in colloidal systems do not originate from thermal fluctuations.

Jamming and unjamming transition of oil-in-water emulsions under continuous temperature change

To analyze the jamming and unjamming transition of oil-in-water emulsions under continuous temperature change, we simulated an emulsion system whose critical volume fraction was 0.3, which was validated with experimental results under oscillatory shear stress. In addition, we calculated the elastic modulus using the phase lag between strain and stress. Through heating and cooling, the emulsion experienced unjamming and jamming. A phenomenon-which is when the elastic modulus does not reach the expected value at the isothermal system-occurred when the emulsion system was cooled. We determined that this phenomenon was caused by the frequency being faster than the relaxation of the deformed droplets. We justified the relation between the frequency and relaxation by simulating the frequency dependency of the difference between the elastic modulus when cooled and the expected value at the same temperature.

Disordered surface vibrations in jammed sphere packings

We study the vibrational properties near a free surface of disordered spring networks derived from jammed sphere packings. In bulk systems, without surfaces, it is well understood that such systems have a plateau in the density of vibrational modes extending down to a frequency scale ω*. This frequency is controlled by ΔZ = 〈Z〉 - 2d, the difference between the average coordination of the spheres and twice the spatial dimension, d, of the system, which vanishes at the jamming transition. In the presence of a free surface we find that there is a density of disordered vibrational modes associated with the surface that extends far below ω*. The total number of these low-frequency surface modes is controlled by ΔZ, and the profile of their decay into the bulk has two characteristic length scales, which diverge as ΔZ(-1/2) and ΔZ(-1) as the jamming transition is approached.

Rigidity percolation by next-nearest-neighbor bonds on generic and regular isostatic lattices

We study rigidity percolation transitions in two-dimensional central-force isostatic lattices, including the square and the kagome lattices, as next-nearest-neighbor bonds ("braces") are randomly added to the system. In particular, we focus on the differences between regular lattices, which are perfectly periodic, and generic lattices with the same topology of bonds but whose sites are at random positions in space. We find that the regular square and kagome lattices exhibit a rigidity percolation transition when the number of braces is ∼LlnL, where L is the linear size of the lattice. This transition exhibits features of both first-order and second-order transitions: The whole lattice becomes rigid at the transition, and a diverging length scale also exists. In contrast, we find that the rigidity percolation transition in the generic lattices occur when the number of braces is very close to the number obtained from Maxwell's law for floppy modes, which is ∼L. The transition in generic lattices is a very sharp first-order-like transition, at which the addition of one brace connects all small rigid regions in the bulk of the lattice, leaving only floppy modes on the edge. We characterize these transitions using numerical simulations and develop analytic theories capturing each transition. Our results relate to other interesting problems, including jamming and bootstrap percolation.

Disordered solids without well-defined transverse phonons: the nature of hard-sphere glasses

We probe the Ioffe-Regel limits of glasses with repulsions near the zero-temperature jamming transition by calculating the dynamical structure factors. The Ioffe-Regel limit (frequency) is reached when the phonon wavelength is comparable to the mean free path, beyond which phonons are no longer well defined. At zero temperature, the transverse Ioffe-Regel frequency vanishes at the jamming transition with a diverging length, but the longitudinal one does not, which excludes the existence of a diverging length associated with the longitudinal excitations. At low temperatures, the transverse and longitudinal Ioffe-Regel frequencies approach zero at the jamminglike transition and glass transition, respectively. As a consequence, glasses between the glass transition and the jamminglike transition, which are hard-sphere glasses in the low temperature limit, can only carry well-defined longitudinal phonons and have an opposite pressure dependence of the ratio of the shear modulus to the bulk modulus from glasses beyond the jamminglike transition.

Vibrational and structural signatures of the crossover between dense glassy and sparse gel-like attractive colloidal packings

We investigate the vibrational modes of quasi-two-dimensional disordered colloidal packings of hard colloidal spheres with short-range attractions as a function of packing fraction. Certain properties of the vibrational density of states (vDOS) are shown to correlate with the density and structure of the samples (i.e., in sparsely versus densely packed samples). Specifically, a crossover from dense glassy to sparse gel-like states is suggested by an excess of phonon modes at low frequency and by a variation in the slope of the vDOS with frequency at low frequency. This change in phonon mode distribution is demonstrated to arise largely from localized vibrations that involve individual and/or small clusters of particles with few local bonds. Conventional order parameters and void statistics did not exhibit obvious gel-glass signatures as a function of volume fraction. These mode behaviors and accompanying structural insights offer a potentially new set of indicators for identification of glass-gel transitions and for assignment of gel-like versus glass-like character to a disordered solid material.

Evolution of covalent networks under cooling: contrasting the rigidity window and jamming scenarios

We study the evolution of structural disorder under cooling in supercooled liquids, focusing on covalent networks. We introduce a model for the energy of networks that incorporates weak noncovalent interactions. We show that at low temperature these interactions considerably affect the network topology near the rigidity transition that occurs as the coordination increases. As a result, this transition becomes mean field and does not present a line of critical points previously argued for, the "rigidity window." Vibrational modes are then not fractons but instead are similar to the anomalous modes observed in packings of particles near jamming. These results suggest an alternative interpretation for the intermediate phase observed in chalcogenides.

Force distribution affects vibrational properties in hard-sphere glasses

We theoretically and numerically study the elastic properties of hard-sphere glasses and provide a real-space description of their mechanical stability. In contrast to repulsive particles at zero temperature, we argue that the presence of certain pairs of particles interacting with a small force f soften elastic properties. This softening affects the exponents characterizing elasticity at high pressure, leading to experimentally testable predictions. Denoting P(f) ~ f(θ(e)), the force distribution of such pairs and ϕ(c) the packing fraction at which pressure diverges, we predict that (i) the density of states has a low-frequency peak at a scale ω*, rising up to it as D(ω) ~ ω(2+a), and decaying above ω* as D(ω) ~ ω(-a) where a = (1 - θ(e))/(3 + θ(e)) and ω is the frequency, (ii) shear modulus and mean-squared displacement are inversely proportional with ⟨δR²⟩ ~ 1/μ ~ (ϕ(c) - ϕ)(κ), where κ = 2 - 2/(3 + θ(e)), and (iii) continuum elasticity breaks down on a scale ℓ(c) ~ 1/√(δz) ~ (ϕ(c) - ϕ)(-b), where b = (1 + θ(e))/(6 + 2θ(e)) and δz = z - 2d, where z is the coordination and d the spatial dimension. We numerically test (i) and provide data supporting that θ(e) ≈ 0.41 in our bidisperse system, independently of system preparation in two and three dimensions, leading to κ ≈ 1.41, a ≈ 0.17, and b ≈ 0.21. Our results for the mean-square displacement are consistent with a recent exact replica computation for d = ∞, whereas some observations differ, as rationalized by the present approach.

Penrose tilings as jammed solids

Penrose tilings form lattices, exhibiting fivefold symmetry and isotropic elasticity, with inhomogeneous coordination much like that of the force networks in jammed systems. Under periodic boundary conditions, their average coordination is exactly four. We study the elastic and vibrational properties of rational approximants to these lattices as a function of unit-cell size N(S) and find that they have of order sqrt[N(S)] zero modes and states of self-stress and yet all their elastic moduli vanish. In their generic form, obtained by randomizing site positions, their elastic and vibrational properties are similar to those of particulate systems at jamming with a nonzero bulk modulus, vanishing shear modulus, and a flat density of states.

Effects of coordination and pressure on sound attenuation, boson peak and elasticity in amorphous solids

Connectedness and applied stress strongly affect elasticity in solids. In various amorphous materials, mechanical stability can be lost either by reducing connectedness or by increasing pressure. We present an effective medium theory of elasticity that extends previous approaches by incorporating the effect of compression, of amplitude e, allowing one to describe quantitative features of sound propagation, transport, the boson peak, and elastic moduli near the elastic instability occurring at a compression ec. The theory disentangles several frequencies characterizing the vibrational spectrum: the onset frequency where strongly-scattered modes appear in the vibrational spectrum, the pressure-independent frequency ω* where the density of states displays a plateau, the boson peak frequency ωBP found to scale as , and the Ioffe-Regel frequency ωIR where scattering length and wavelength become equal. We predict that sound attenuation crosses over from ω(4) to ω(2) behaviour at ω0, consistent with observations in glasses. We predict that a frequency-dependent length scale ls(ω) and speed of sound ν(ω) characterize vibrational modes, and could be extracted from scattering data. One key result is the prediction of a flat diffusivity above ω0, in agreement with previously unexplained observations. We find that the shear modulus does not vanish at the elastic instability, but drops by a factor of 2. We check our predictions in packings of soft particles and study the case of covalent networks and silica, for which we predict ωIR ≈ ωBP. Overall, our approach unifies sound attenuation, transport and length scales entering elasticity in a single framework where disorder is not the main parameter controlling the boson peak, in agreement with observations. This framework leads to a phase diagram where various glasses can be placed, connecting microscopic structure to vibrational properties.

Contact changes near jamming

We probe the onset and effect of contact changes in soft harmonic particle packings which are sheared quasistatically. We find that the first contact changes are the creation or breaking of contacts on a single particle. We characterize the critical strain, statistics of breaking versus making a contact, and ratio of shear modulus before and after such events, and explain their finite size scaling relations. For large systems at finite pressure, the critical strain vanishes but the ratio of shear modulus before and after a contact change approaches one: linear response remains relevant in large systems. For finite systems close to jamming the critical strain also vanishes, but here linear response already breaks down after a single contact change.

Jamming in finite systems: stability, anisotropy, fluctuations, and scaling

Athermal packings of soft repulsive spheres exhibit a sharp jamming transition in the thermodynamic limit. Upon further compression, various structural and mechanical properties display clean power-law behavior over many decades in pressure. As with any phase transition, the rounding of such behavior in finite systems close to the transition plays an important role in understanding the nature of the transition itself. The situation for jamming is surprisingly rich: the assumption that jammed packings are isotropic is only strictly true in the large-size limit, and finite-size has a profound effect on the very meaning of jamming. Here, we provide a comprehensive numerical study of finite-size effects in sphere packings above the jamming transition, focusing on stability as well as the scaling of the contact number and the elastic response.

Contact nonlinearities and linear response in jammed particulate packings

Packings of frictionless athermal particles that interact only when they overlap experience a jamming transition as a function of packing density. Such packings provide the foundation for the theory of jamming. This theory rests on the observation that, despite the multitude of disordered configurations, the mechanical response to linear order depends only on the distance to the transition. We investigate the validity and utility of such measurements that invoke the harmonic approximation and show that, despite particles coming in and out of contact, there is a well-defined linear regime in the thermodynamic limit.

Breakdown of continuum elasticity in amorphous solids

We show numerically that the response of simple amorphous solids (elastic networks and particle packings) to a local force dipole is characterized by a lengthscale lc that diverges as unjamming is approached as lc ∼ (z - 2d)(-1/2), where z ≥ 2d is the mean coordination, and d is the spatial dimension, at odds with previous numerical claims. We also show how the magnitude of the lengthscale lc is amplified by the presence of internal stresses in the disordered solid. Our data suggests a divergence of lc ∼ (pc - p)(-1/4) with proximity to a critical internal stress pc at which soft elastic modes become unstable.

Stress-dependent normal-mode frequencies from the effective mass of granular matter

A zero-temperature critical point has been invoked to control the anomalous behavior of granular matter as it approaches jamming or mechanical arrest. Criticality manifests itself in an anomalous spectrum of low-frequency normal modes and scaling behavior near the jamming transition. The critical point may explain the peculiar mechanical properties of dissimilar systems such as glasses and granular materials. Here we study the critical scenario via an experimental measurement of the normal modes frequencies of granular matter under stress from a pole decomposition analysis of the effective mass. We extract a complex-valued characteristic frequency which displays scaling |ω (σ)| ∼ σΩ' with vanishing stress σ for a variety of granular systems. The critical exponent is smaller than that predicted by mean-field theory opening new challenges to explain the exponent for frictional and dissipative granular matter. Our results shed light on the anomalous behavior of stress-dependent acoustics and attenuation in granular materials near the jamming transition.

Aspects of jamming in two-dimensional athermal frictionless systems

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

Physics in ordered and disordered colloidal matter composed of poly(N-isopropylacrylamide) microgel particles

This review collects and describes experiments that employ colloidal suspensions to probe physics in ordered and disordered solids and related complex fluids. The unifying feature of this body of work is its clever usage of poly(N-isopropylacrylamide) (PNIPAM) microgel particles. These temperature-sensitive colloidal particles provide experimenters with a 'knob' for in situ control of particle size, particle interaction and particle packing fraction that, in turn, influence the structural and dynamical behavior of the complex fluids and solids. A brief summary of PNIPAM particle synthesis and properties is given, followed by a synopsis of current activity in the field. The latter discussion describes a variety of soft matter investigations including those that explore formation and melting of crystals and clusters, and those that probe structure, rearrangement and rheology of disordered (jammed/glassy) and partially ordered matter. The review, therefore, provides a snapshot of a broad range of physics phenomenology which benefits from the unique properties of responsive microgel particles.

High-velocity drag friction in granular media near the jamming point

Drag friction that acts on a disk in a two-dimensional granular medium is studied at high packing fractions. We concentrate on a high-velocity region, in which the dynamic component of the force, obtained as an average of a strongly fluctuating force, clearly scales with velocity squared. We find that the total force composed of dynamic and static components, as well as its fluctuation, diverges with practically the same exponent as the packing fraction approaches the jamming point. To explain the critical behavior, we propose a simple theory equipped with a diverging length scale, which agrees well with the data and elucidates physical pictures for the divergence.

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

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

Anharmonic damping of terahertz acoustic waves in a network glass and its effect on the density of vibrational states

We report the observation, by means of high-resolution inelastic x-ray scattering, of an unusually large temperature dependence of the sound attenuation of a network glass at terahertz frequency, an unprecedentedly observed phenomenon. The anharmonicity can be ascribed to the interaction between the propagating acoustic wave and the bath of thermal vibrations. At low temperatures the sound attenuation follows a Rayleigh-Gans scattering law. As the temperature is increased the anharmonic process sets in, resulting in an almost quadratic frequency dependence of the damping in the entire frequency range. We show that the temperature variation of the sound damping accounts quantitatively for the temperature dependence of the density of vibrational states.

Role of disorder in the thermodynamics and atomic dynamics of glasses

We measured the density of vibrational states (DOS) and the specific heat of various glassy and crystalline polymorphs of SiO2. The typical (ambient) glass shows a well-known excess of specific heat relative to the typical crystal (α-quartz). This, however, holds when comparing a lower-density glass to a higher-density crystal. For glassy and crystalline polymorphs with matched densities, the DOS of the glass appears as the smoothed counterpart of the DOS of the corresponding crystal; it reveals the same number of the excess states relative to the Debye model, the same number of all states in the low-energy region, and it provides the same specific heat. This shows that glasses have higher specific heat than crystals not due to disorder, but because the typical glass has lower density than the typical crystal.

Phonon dispersion and elastic moduli of two-dimensional disordered colloidal packings of soft particles with frictional interactions

Particle tracking and displacement covariance matrix techniques are employed to investigate the phonon dispersion relations of two-dimensional colloidal glasses composed of soft, thermoresponsive microgel particles whose temperature-sensitive size permits in situ variation of particle packing fraction. Bulk, B, and shear, G, moduli of the colloidal glasses are extracted from the dispersion relations as a function of packing fraction, and variation of the ratio G/B with packing fraction is found to agree quantitatively with predictions for jammed packings of frictional soft particles. In addition, G and B individually agree with numerical predictions for frictional particles. This remarkable level of agreement enabled us to extract an energy scale for the interparticle interaction from the individual elastic constants and to derive an approximate estimate for the interparticle friction coefficient.

Shear shocks in fragile networks

A minimal model for studying the mechanical properties of amorphous solids is a disordered network of point masses connected by unbreakable springs. At a critical value of its mean connectivity, such a network becomes fragile: it undergoes a rigidity transition signaled by a vanishing shear modulus and transverse sound speed. We investigate analytically and numerically the linear and nonlinear visco-elastic response of these fragile solids by probing how shear fronts propagate through them. Our approach, which we tentatively label shear front rheology, provides an alternative route to standard oscillatory rheology. In the linear regime, we observe at late times a diffusive broadening of the fronts controlled by an effective shear viscosity that diverges at the critical point. No matter how small the microscopic coefficient of dissipation, strongly disordered networks behave as if they were overdamped because energy is irreversibly leaked into diverging nonaffine fluctuations. Close to the transition, the regime of linear response becomes vanishingly small: the tiniest shear strains generate strongly nonlinear shear shock waves qualitatively different from their compressional counterparts in granular media. The inherent nonlinearities trigger an energy cascade from low to high frequency components that keep the network away from attaining the quasi-static limit. This mechanism, reminiscent of acoustic turbulence, causes a superdiffusive broadening of the shock width.

Using shape anisotropy to toughen disordered nanoparticle assemblies

Assemblies of disordered nanoparticles constitute an important class of materials that have numerous applications in energy conversion and storage, electronics, photonics, and sensing. One major roadblock that limits the widespread utilization of disordered nanoparticle assemblies (DNAs) is their poor damage tolerance; they fracture under small loads and, thus, have low toughness. The absence of fundamental understanding on the mechanical behavior and failure mechanism of disordered nanoparticle assemblies makes it even more challenging to develop new strategies to toughen these structures without compromising their mechanical strength. Here we show the formation of shear bands, highly localized regions of mechanical strain that prelude fracture, in disordered assemblies of spherical nanoparticles, which bear striking resemblance to the deformation mechanism of a different class of disordered materials, metallic glasses. We demonstrate that anisotropic nanoparticles greatly suppress shear band formation and toughen nanoparticle packings without sacrificing their strength, implying that tuning constituent anisotropy can be used to enhance toughness in disordered packings of nanoparticles.

Lattice model of correlated forces in granular solids near jamming

We have devised a lattice model to study force correlations in jamming granular solids in d=2 dimensions. We perform biased Monte Carlo simulations, favoring configurations with more bonds that bear no force, to "starve" the network of bonds and thereby control the distance from the isostatic point J. Increasingly long-ranged correlations are visible as point J is approached, not in the structure of the network of force-bearing bonds but in the spatial extent of perturbations of the force magnitudes consistent with a given starved network. The correlation length so defined diverges as the isostatic point is approached as a power law with an exponent of about ξ~δZ(-5). This divergence is much stronger than for the length scale of "soft modes" observed in jammed systems approaching point J from above.

Effective-medium theory of a filamentous triangular lattice

We present an effective-medium theory that includes bending as well as stretching forces, and we use it to calculate the mechanical response of a diluted filamentous triangular lattice. In this lattice, bonds are central-force springs, and there are bending forces between neighboring bonds on the same filament. We investigate the diluted lattice in which each bond is present with a probability p. We find a rigidity threshold p(b) which has the same value for all positive bending rigidity and a crossover characterizing bending, stretching, and bend-stretch coupled elastic regimes controlled by the central-force rigidity percolation point at p(CF)=/~2/3 of the lattice when fiber bending rigidity vanishes.

Dynamic critical response in damped random spring networks

The isostatic state plays a central role in organizing the response of many amorphous materials. We construct a diverging length scale in nearly isostatic spring networks that is defined both above and below isostaticity and at finite frequencies and relate the length scale to viscoelastic response. Numerical measurements verify that proximity to isostaticity controls the viscosity, shear modulus, and creep of random networks.

Uniform shock waves in disordered granular matter

The confining pressure P is perhaps the most important parameter controlling the properties of granular matter. Strongly compressed granular media are, in many respects, simple solids in which elastic perturbations travel as ordinary phonons. However, the speed of sound in granular aggregates continuously decreases as the confining pressure decreases, completely vanishing at the jamming-unjamming transition. This anomalous behavior suggests that the transport of energy at low pressures should not be dominated by phonons. In this work we use simulations and theory to show how the response of granular systems becomes increasingly nonlinear as pressure decreases. In the low-pressure regime the elastic energy is found to be mainly transported through nonlinear waves and shocks. We numerically characterize the propagation speed, shape, and stability of these shocks and model the dependence of the shock speed on pressure and impact intensity by a simple analytical approach.

Marginal stability constrains force and pair distributions at random close packing

The requirement that packings of frictionless hard spheres, arguably the simplest structural glass, cannot be compressed by rearranging their network of contacts is shown to yield a new constraint on their microscopic structure. This constraint takes the form a bound between the distribution of contact forces P(f) and the pair distribution function g(r): if P(f)∼f(θ) and g(r)∼(r-σ(0))(-γ), where σ(0) is the particle diameter, one finds that γ ≥ 1/(2 + θ). This bound plays a role similar to those found in some glassy materials with long-range interactions, such as the Coulomb gap in Anderson insulators or the distribution of local fields in mean-field spin glasses. There are grounds to believe that this bound is saturated, yielding a mechanism to explain the avalanches of rearrangements with power-law statistics that govern plastic flow in packings.

Heterogeneous Structure, Heterogeneous Dynamics, and Complex Behavior in Two-Dimensional Liquids

Analysis of the metrical and topological features of the local structure in a freezing two-dimensional Lennard-Jones system found that in a narrow strip [Formula: see text] of thermodynamic states close to the melting line, the liquid becomes a complex liquid characterized by a super-Arrhenius increase of relaxation times, stretched-exponential decay of correlations in time, and a power-law distribution of waiting times for changes in the local order. In [Formula: see text], the structure of the liquid and its dynamics are spatially heterogeneous; the sizes of ordered clusters are power-law distributed. Those features are governed by local structure evolution between solid-like and liquid-like (disordered) patterns. The liquid inside the strip [Formula: see text] gives a unique opportunity to study how heterogeneous structure, dynamics and complexity are intertwined with each other on a microscopic level.

Soft-sphere packings at finite pressure but unstable to shear

When are athermal soft-sphere packings jammed? Any experimentally relevant definition must, at the very least, require a jammed packing to resist shear. We demonstrate that widely used (numerical) protocols, in which particles are compressed together, can and do produce packings that are unstable to shear-and that the probability of generating such packings reaches one near jamming. We introduce a new protocol which, by allowing the system to explore different box shapes as it equilibrates, generates truly jammed packings with strictly positive shear moduli G. For these packings, the scaling of the average of G is consistent with earlier results, while the probability distribution P(G) exhibits novel and rich scalings.

Finite-size scaling at the jamming transition

We present an analysis of finite-size effects in jammed packings of N soft, frictionless spheres at zero temperature. There is a 1/N correction to the discrete jump in the contact number at the transition so that jammed packings exist only above isostaticity. As a result, the canonical power-law scalings of the contact number and elastic moduli break down at low pressure. These quantities exhibit scaling collapse with a nontrivial scaling function, demonstrating that the jamming transition can be considered a phase transition. Scaling is achieved as a function of N in both two and three dimensions, indicating an upper critical dimension of 2.

Force propagation in isostatic granular packs

We investigate how forces spread through frictionless granular packs at the jamming transition. Previous work has indicated that such packs are isostatic, and thus obey a null stress law which, independent of the packing history, causes rays of stress to propagate away from a point force at oblique angles. Prior verifications of the null stress law have used a sequential packing method which yields packs with anisotropic packing histories. We create packs without this anisotropy, and then later break the symmetry by adding a boundary. Our isotropic packs are very sensitive, and their responses to point forces diverge wildly, indicating that they cannot be described by any continuum stress model. We stabilize the packs by supplying an additional boundary, which makes the response much more regular. The response of the stabilized packs resembles what one would expect in a hyperstatic pack, despite the isostatic bulk. The expected stress rays characteristic of null stress behavior are not present. This suggests that isostatic packs do not need to obey a null stress condition. We argue that the rays may arise instead from more simple geometric considerations, such as preferred contact angles between beads.