Pinching a glass reveals key properties of its soft spots.

It is now well established that glasses feature quasilocalized nonphononic excitations-coined "soft spots"-, which follow a universal [Formula: see text] density of states in the limit of low frequencies ω. All glass-specific properties, such as the dependence on the preparation protocol or composition, are encapsulated in the nonuniversal prefactor of the universal [Formula: see text] law. The prefactor, however, is a composite quantity that incorporates information both about the number of quasilocalized nonphononic excitations and their characteristic stiffness, in an apparently inseparable manner. We show that by pinching a glass-i.e., by probing its response to force dipoles-one can disentangle and independently extract these two fundamental pieces of physical information. This analysis reveals that the number of quasilocalized nonphononic excitations follows a Boltzmann-like law in terms of the parent temperature from which the glass is quenched. The latter, sometimes termed the fictive (or effective) temperature, plays important roles in nonequilibrium thermodynamic approaches to the relaxation, flow, and deformation of glasses. The analysis also shows that the characteristic stiffness of quasilocalized nonphononic excitations can be related to their characteristic size, a long sought-for length scale. These results show that important physical information, which is relevant for various key questions in glass physics, can be obtained through pinching a glass.

Statistical mechanics of local force dipole responses in computer glasses.

Soft quasilocalized modes (QLMs) are universally featured by structural glasses quenched from a melt, and are involved in several glassy anomalies such as the low-temperature scaling of their thermal conductivity and specific heat, and sound attenuation at intermediate frequencies. In computer glasses, QLMs may assume the form of harmonic vibrational modes under a narrow set of circumstances; however, direct access to their full distribution over frequency is hindered by hybridizations of QLMs with other low-frequency modes (e.g., phonons). Previous studies to overcome this issue have demonstrated that the response of a glass to local force dipoles serves as a good proxy for its QLMs; we, therefore, study here the statistical-mechanical properties of these responses in computer glasses, over a large range of glass stabilities and in various spatial dimensions, with the goal of revealing properties of the yet-inaccessible full distribution of QLMs' frequencies. We find that as opposed to the spatial-dimension-independent universal distribution of QLMs' frequencies ω (and, consequently, also of their stiffness κ = ω2), the distribution of stiffnesses associated with responses to local force dipoles features a (weak) dependence on spatial dimension. We rationalize this dependence by introducing a lattice model that incorporates both the real-space profiles of QLMs-associated with dimension-dependent long-range elastic fields-and the universal statistical properties of their frequencies. Based on our findings, we propose a conjecture about the form of the full distribution of QLMs' frequencies and its protocol-dependence. Finally, we discuss possible connections of our findings to basic aspects of glass formation and deformation.

Particle pinning suppresses spinodal criticality in the shear-banding instability.

Strained amorphous solids often fail mechanically by creating a shear band. It had been understood that the shear-banding instability is usefully described as crossing a spinodal point (with disorder) in an appropriate thermodynamic description. It remained contested, however, whether the spinodal is critical (with divergent correlation length) or not. Here we offer evidence for critical spinodal by using particle pinning. For a finite concentration of pinned particles the correlation length is bounded by the average distance between pinned particles, but without pinning it is bounded by the system size.

Wave attenuation in glasses: Rayleigh and generalized-Rayleigh scattering scaling.

The attenuation of long-wavelength phonons (waves) by glassy disorder plays a central role in various glass anomalies, yet it is neither fully characterized nor fully understood. Of particular importance is the scaling of the attenuation rate Γ(k) with small wavenumbers k → 0 in the thermodynamic limit of macroscopic glasses. Here, we use a combination of theory and extensive computer simulations to show that the macroscopic low-frequency behavior emerges at intermediate frequencies in finite-size glasses, above a recently identified crossover wavenumber k†, where phonons are no longer quantized into bands. For k < k†, finite-size effects dominate Γ(k), which is quantitatively described by a theory of disordered phonon bands. For k > k†, we find that Γ(k) is affected by the number of quasilocalized nonphononic excitations, a generic signature of glasses that feature a universal density of states. In particular, we show that in a frequency range in which this number is small, Γ(k) follows a Rayleigh scattering scaling ∼k¯d+1 (¯d is the spatial dimension) and that in a frequency range in which this number is sufficiently large, the recently observed generalized-Rayleigh scaling of the form ∼k¯d+1 log(k0/k) emerges (k0 > k† is a characteristic wavenumber). Our results suggest that macroscopic glasses-and, in particular, glasses generated by conventional laboratory quenches that are known to strongly suppress quasilocalized nonphononic excitations-exhibit Rayleigh scaling at the lowest wavenumbers k and a crossover to generalized-Rayleigh scaling at higher k. Some supporting experimental evidence from recent literature is presented.

Robustness of mean field theory for hard sphere models.

In very recent work the mean field theory of the jamming transition in infinite-dimensional hard sphere models was presented. Surprisingly, this theory predicts quantitatively the numerically determined characteristics of jamming in two and three dimensions. This is a rare and unusual finding. Here we argue that this agreement is nongeneric: only for hard sphere models does it happen that sufficiently close to the jamming density (at any temperature) the effective interactions are binary, in agreement with mean field theory, justifying the truncation of many-body interactions (which is the exact protocol in infinite dimensions). Any softening of the bare hard sphere interactions results in many-body effective interactions that are not mean field at any density, making the d=∞ results not applicable.

Shear bands as manifestation of a criticality in yielding amorphous solids.

Amorphous solids increase their stress as a function of an applied strain until a mechanical yield point whereupon the stress cannot increase anymore, afterward exhibiting a steady state with a constant mean stress. In stress-controlled experiments, the system simply breaks when pushed beyond this mean stress. The ubiquity of this phenomenon over a huge variety of amorphous solids calls for a generic theory that is free of microscopic details. Here, we offer such a theory: The mechanical yield is a thermodynamic phase transition, where yield occurs as a spinodal phenomenon. At the spinodal point, there exists a divergent correlation length that is associated with the system-spanning instabilities (also known as shear bands), which are typical to the mechanical yield. The theory, the order parameter used, and the correlation functions that exhibit the divergent correlation length are universal in nature and can be applied to any amorphous solids that undergo mechanical yield.

Mechanical failure in amorphous solids: Scale-free spinodal criticality.

The mechanical failure of amorphous media is a ubiquitous phenomenon from material engineering to geology. It has been noticed for a long time that the phenomenon is "scale-free," indicating some type of criticality. In spite of attempts to invoke "Self-Organized Criticality," the physical origin of this criticality, and also its universal nature, being quite insensitive to the nature of microscopic interactions, remained elusive. Recently we proposed that the precise nature of this critical behavior is manifested by a spinodal point of a thermodynamic phase transition. Demonstrating this requires the introduction of an "order parameter" that is suitable for distinguishing between disordered amorphous systems. At the spinodal point there exists a divergent correlation length which is associated with the system-spanning instabilities (known also as shear bands) which are typical to the mechanical yield. The theory, the order parameter used and the correlation functions which exhibit the divergent correlation length are universal in nature and can be applied to any amorphous solid that undergoes mechanical yield. The phenomenon is seen at its sharpest in athermal systems, as is explained below; in this paper we extend the discussion also to thermal systems, showing that at sufficiently high temperatures the spinodal phenomenon is destroyed by thermal fluctuations.

Emergent interparticle interactions in thermal amorphous solids.

Amorphous media at finite temperatures, be them liquids, colloids, or glasses, are made of interacting particles that move chaotically due to thermal energy, continuously colliding and scattering off each other. When the average configuration in these systems relaxes only at long times, one can introduce effective interactions that keep the mean positions in mechanical equilibrium. We introduce a framework to determine the effective force laws that define an effective Hessian that can be employed to discuss stability properties and the density of states of the amorphous system. We exemplify the approach with a thermal glass of hard spheres; these experience zero forces when not in contact and infinite forces when they touch. Close to jamming we recapture the effective interactions that at temperature T depend on the gap h between spheres as T/h [C. Brito and M. Wyart, Europhys. Lett. 76, 149 (2006)EULEEJ0295-507510.1209/epl/i2006-10238-x]. For hard spheres at lower densities or for systems whose binary bare interactions are longer ranged (at any density), the emergent force laws include ternary, quaternary, and generally higher-order many-body terms, leading to a temperature-dependent effective Hessian.

Breakdown of nonlinear elasticity in amorphous solids at finite temperatures.

It is known [H. G. E. Hentschel et al., Phys. Rev. E 83, 061101 (2011)PLEEE81539-375510.1103/PhysRevE.83.061101] that amorphous solids at zero temperature do not possess a nonlinear elasticity theory: besides the shear modulus, which exists, none of the higher order coefficients exist in the thermodynamic limit. Here we show that the same phenomenon persists up to temperatures comparable to that of the glass transition. The zero-temperature mechanism due to the prevalence of dangerous plastic modes of the Hessian matrix is replaced by anomalous stress fluctuations that lead to the divergence of the variances of the higher order elastic coefficients. The conclusion is that in amorphous solids elasticity can never be decoupled from plasticity: the nonlinear response is very substantially plastic.

Mechanical Yield in Amorphous Solids: A First-Order Phase Transition.

Amorphous solids yield at a critical value of the strain (in strain-controlled experiments); for larger strains, the average stress can no longer increase-the system displays an elastoplastic steady state. A long-standing riddle in the materials community is what the difference is between the microscopic states of the material before and after yield. Explanations in the literature are material specific, but the universality of the phenomenon begs a universal answer. We argue here that there is no fundamental difference in the states of matter before and after yield, but the yield is a bona fide first-order phase transition between a highly restricted set of possible configurations residing in a small region of phase space to a vastly rich set of configurations which include many marginally stable ones. To show this, we employ an order parameter of universal applicability, independent of the microscopic interactions, that is successful in quantifying the transition in an unambiguous manner.

Numerical detection of the Gardner transition in a mean-field glass former.

Recent theoretical advances predict the existence, deep into the glass phase, of a novel phase transition, the so-called Gardner transition. This transition is associated with the emergence of a complex free energy landscape composed of many marginally stable sub-basins within a glass metabasin. In this study, we explore several methods to detect numerically the Gardner transition in a simple structural glass former, the infinite-range Mari-Kurchan model. The transition point is robustly located from three independent approaches: (i) the divergence of the characteristic relaxation time, (ii) the divergence of the caging susceptibility, and (iii) the abnormal tail in the probability distribution function of cage order parameters. We show that the numerical results are fully consistent with the theoretical expectation. The methods we propose may also be generalized to more realistic numerical models as well as to experimental systems.

Calorimetric glass transition in a mean-field theory approach.

The study of the properties of glass-forming liquids is difficult for many reasons. Analytic solutions of mean-field models are usually available only for systems embedded in a space with an unphysically high number of spatial dimensions; on the experimental and numerical side, the study of the properties of metastable glassy states requires thermalizing the system in the supercooled liquid phase, where the thermalization time may be extremely large. We consider here a hard-sphere mean-field model that is solvable in any number of spatial dimensions; moreover, we easily obtain thermalized configurations even in the glass phase. We study the 3D version of this model and we perform Monte Carlo simulations that mimic heating and cooling experiments performed on ultrastable glasses. The numerical findings are in good agreement with the analytical results and qualitatively capture the features of ultrastable glasses observed in experiments.

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

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

Dynamical arrest with zero complexity: The unusual behavior of the spherical Blume-Emery-Griffiths disordered model.

The short- and long-time dynamics of model systems undergoing a glass transition with apparent inversion of Kauzmann and dynamical arrest glass transition lines is investigated. These models belong to the class of the spherical mean-field approximation of a spin-1 model with p-body quenched disordered interaction, with p>2, termed spherical Blume-Emery-Griffiths models. Depending on temperature and chemical potential the system is found in a paramagnetic or in a glassy phase and the transition between these phases can be of a different nature. In specific regions of the phase diagram coexistence of low-density and high-density paramagnets can occur, as well as the coexistence of spin-glass and paramagnetic phases. The exact static solution for the glassy phase is known to be obtained by the one-step replica symmetry breaking ansatz. Different scenarios arise for both the dynamic and the thermodynamic transitions. These include: (i) the usual random first-order transition (Kauzmann-like) for mean-field glasses preceded by a dynamic transition, (ii) a thermodynamic first-order transition with phase coexistence and latent heat, and (iii) a regime of apparent inversion of static transition line and dynamic transition lines, the latter defined as a nonzero complexity line. The latter inversion, though, turns out to be preceded by a dynamical arrest line at higher temperature. Crossover between different regimes is analyzed by solving mode-coupling-theory equations near the boundaries of paramagnetic solutions and the relationship with the underlying statics is discussed.