Searching for the Gardner Transition in Glassy Glycerol.

We search for a Gardner transition in glassy glycerol, a standard molecular glass, measuring the third harmonics cubic susceptibility χ_{3}^{(3)} from slightly below the usual glass transition temperature down to 10 K. According to the mean-field picture, if local motion within the glass were becoming highly correlated due to the emergence of a Gardner phase then χ_{3}^{(3)}, which is analogous to the dynamical spin-glass susceptibility, should increase and diverge at the Gardner transition temperature T_{G}. We find instead that upon cooling |χ_{3}^{(3)}| decreases by several orders of magnitude and becomes roughly constant in the regime 100-10 K. We rationalize our findings by assuming that the low temperature physics is described by localized excitations weakly interacting via a spin-glass dipolar pairwise interaction in a random magnetic field. Our quantitative estimations show that the spin-glass interaction is twenty to fifty times smaller than the local random field contribution, thus rationalizing the absence of the spin-glass Gardner phase. This hints at the fact that a Gardner phase may be suppressed in standard molecular glasses, but it also suggests ways to favor its existence in other amorphous solids and by changing the preparation protocol.

Universality of jamming of nonspherical particles.

Amorphous packings of nonspherical particles such as ellipsoids and spherocylinders are known to be hypostatic: The number of mechanical contacts between particles is smaller than the number of degrees of freedom, thus violating Maxwell's mechanical stability criterion. In this work, we propose a general theory of hypostatic amorphous packings and the associated jamming transition. First, we show that many systems fall into a same universality class. As an example, we explicitly map ellipsoids into a system of "breathing" particles. We show by using a marginal stability argument that in both cases jammed packings are hypostatic and that the critical exponents related to the contact number and the vibrational density of states are the same. Furthermore, we introduce a generalized perceptron model which can be solved analytically by the replica method. The analytical solution predicts critical exponents in the same hypostatic jamming universality class. Our analysis further reveals that the force and gap distributions of hypostatic jamming do not show power-law behavior, in marked contrast to the isostatic jamming of spherical particles. Finally, we confirm our theoretical predictions by numerical simulations.

Critical Jammed Phase of the Linear Perceptron.

Criticality in statistical physics naturally emerges at isolated points in the phase diagram. Jamming of spheres is not an exception: varying density, it is the critical point that separates the unjammed phase where spheres do not overlap and the jammed phase where they cannot be arranged without overlaps. The same remains true in more general constraint satisfaction problems with continuous variables where jamming coincides with the (protocol dependent) satisfiability transition point. In this work we show that by carefully choosing the cost function to be minimized, the region of criticality extends to occupy a whole region of the jammed phase. As a working example, we consider the spherical perceptron with a linear cost function in the unsatisfiable jammed phase and we perform numerical simulations which show critical power laws emerging in the configurations obtained minimizing the linear cost function. We develop a scaling theory to compute the emerging critical exponents.

Jamming in Multilayer Supervised Learning Models.

Critical jamming transitions are characterized by an astonishing degree of universality. Analytic and numerical evidence points to the existence of a large universality class that encompasses finite and infinite dimensional spheres and continuous constraint satisfaction problems (CCSP) such as the nonconvex perceptron and related models. In this Letter we investigate multilayer neural networks (MLNN) learning random associations as models for CCSP that could potentially define different jamming universality classes. As opposed to simple perceptrons and infinite dimensional spheres, which are described by a single effective field in terms of which the constraints appear to be one dimensional, the description of MLNN involves multiple fields, and the constraints acquire a multidimensional character. We first study the models numerically and show that similarly to the perceptron, whenever jamming is isostatic, the sphere universality class is recovered, we then write the exact mean-field equations for the models and identify a dimensional reduction mechanism that leads to a scaling regime identical to the one of spheres.

Can a Large Packing be Assembled from Smaller Ones?

We consider zero temperature packings of soft spheres that undergo a jamming to unjamming transition as a function of packing fraction. We compare differences in the structure, as measured from the contact statistics, of a finite subsystem of a large packing to a whole packing with periodic boundaries of an equivalent size and pressure. We find that the fluctuations of the ensemble of whole packings are smaller than those of the ensemble of subsystems. Convergence of these two quantities appears to occur at very large systems, which are usually not attainable in numerical simulations. Finding differences between packings in two dimensions and three dimensions, we also consider four dimensions and mean-field models, and find that they show similar system size dependence. Mean-field critical exponents appear to be consistent with the 3D and 4D packings, suggesting they are above the upper critical dimension. We also find that the convergence as a function of system size to the thermodynamic limit is characterized by two different length scales. We argue that this is the result of the system being above the upper critical dimension.

Microscopic Theory of Two-Step Yielding in Attractive Colloids.

Attractive colloids display two distinct amorphous solid phases: the attractive glass due to particle bonding and the repulsive glass due to the hard-core repulsion. By means of a microscopic mean field approach, we analyze their response to a quasistatic shear strain. We find that the presence of two distinct interaction length scales may result in a sharp two-step yielding process, which can be associated with a hysteretic stress response or with a reversible but nonmonotonic stress-strain curve. We derive a generic phase diagram characterized by two distinct yielding lines, an inverse yielding and a critical point separating the hysteretic and reversible regimes. Our results should be applicable to a large class of glassy materials characterized by two distinct interaction length scales.

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.

Shear Yielding and Shear Jamming of Dense Hard Sphere Glasses.

We investigate the response of dense hard sphere glasses to a shear strain in a wide range of pressures ranging from the glass transition to the infinite-pressure jamming point. The phase diagram in the density-shear strain plane is calculated analytically using the mean-field infinite-dimensional solution. We find that just above the glass transition, the glass generically yields at a finite shear strain. The yielding transition in the mean-field picture is a spinodal point in presence of disorder. At higher densities, instead, we find that the glass generically jams at a finite shear strain: the jamming transition prevents yielding. The shear yielding and shear jamming lines merge in a critical point, close to which the system yields at extremely large shear stress. Around this point, highly nontrivial yielding dynamics, characterized by system-spanning disordered fractures, is expected.

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.

Quantitative field theory of the glass transition.

We develop a full microscopic replica field theory of the dynamical transition in glasses. By studying the soft modes that appear at the dynamical temperature, we obtain an effective theory for the critical fluctuations. This analysis leads to several results: we give expressions for the mean field critical exponents, and we analytically study the critical behavior of a set of four-points correlation functions, from which we can extract the dynamical correlation length. Finally, we can obtain a Ginzburg criterion that states the range of validity of our analysis. We compute all these quantities within the hypernetted chain approximation for the Gibbs free energy, and we find results that are consistent with numerical simulations.

A note on weakly discontinuous dynamical transitions.

We analyze mode coupling discontinuous transition in the limit of vanishing discontinuity, approaching the so called "A(3)" point. In these conditions structural relaxation and fluctuations appear to have universal form independent from the details of the system. The analysis of this limiting case suggests new ways for looking at the mode coupling equations in the general case.

Static replica approach to critical correlations in glassy systems.

We discuss the slow relaxation phenomenon in glassy systems by means of replicas by constructing a static field theory approach to the problem. At the mean field level we study how criticality in the four point correlation functions arises because of the presence of soft modes and we derive an effective replica field theory for these critical fluctuations. By using this at the gaussian level we obtain many physical quantities: the correlation length, the exponent parameter that controls the mode-coupling dynamical exponents for the two-point correlation functions, and the prefactor of the critical part of the four point correlation functions. Moreover, we perform a one-loop computation in order to identify the region in which the mean field gaussian approximation is valid. The result is a Ginzburg criterion for the glass transition. We define and compute in this way a proper Ginzburg number. Finally, we present numerical values of all these quantities obtained from the hypernetted chain approximation for the replicated liquid theory.

High-dimensional optimization under nonconvex excluded volume constraints.

We consider high-dimensional random optimization problems where the dynamical variables are subjected to nonconvex excluded volume constraints. We focus on the case in which the cost function is a simple quadratic cost and the excluded volume constraints are modeled by a perceptron constraint satisfaction problem. We show that depending on the density of constraints, one can have different situations. If the number of constraints is small, one typically has a phase where the ground state of the cost function is unique and sits on the boundary of the island of configurations allowed by the constraints. In this case, there is a hypostatic number of marginally satisfied constraints. If the number of constraints is increased one enters a glassy phase where the cost function has many local minima sitting again on the boundary of the regions of allowed configurations. At the phase transition point, the total number of marginally satisfied constraints becomes equal to the number of degrees of freedom in the problem and therefore we say that these minima are isostatic. We conjecture that by increasing further the constraints the system stays isostatic up to the point where the volume of available phase space shrinks to zero. We derive our results using the replica method and we also analyze a dynamical algorithm, the Karush-Kuhn-Tucker algorithm, through dynamical mean-field theory and we show how to recover the results of the replica approach in the replica symmetric phase.

A stability-reversibility map unifies elasticity, plasticity, yielding, and jamming in hard sphere glasses.

Amorphous solids, such as glasses, have complex responses to deformations, with substantial consequences in material design and applications. In this respect, two intertwined aspects are important: stability and reversibility. It is crucial to understand, on the one hand, how a glass may become unstable due to increased plasticity under shear deformations, and, on the other hand, to what extent the response is reversible, meaning how much a system is able to recover the original configuration once the perturbation is released. Here, we focus on assemblies of hard spheres as the simplest model of amorphous solids such as colloidal glasses and granular matter. We prepare glass states quenched from equilibrium supercooled liquid states, which are obtained by using the swap Monte Carlo algorithm and correspond to a wide range of structural relaxation time scales. We exhaustively map out their stability and reversibility under volume and shear strains using extensive numerical simulations. The region on the volume-shear strain phase diagram where the original glass state remains solid is bounded by the shear yielding and the shear jamming lines that meet at a yielding-jamming crossover point. This solid phase can be further divided into two subphases: the stable glass phase, where the system deforms purely elastically and is totally reversible, and the marginal glass phase, where it experiences stochastic plastic deformations at mesoscopic scales and is partially irreversible. The details of the stability-reversibility map depend strongly on the quality of annealing of the glass. This study provides a unified framework for understanding elasticity, plasticity, yielding, and jamming in amorphous solids.

Fractal free energy landscapes in structural glasses.

Glasses are amorphous solids whose constituent particles are caged by their neighbours and thus cannot flow. This sluggishness is often ascribed to the free energy landscape containing multiple minima (basins) separated by high barriers. Here we show, using theory and numerical simulation, that the landscape is much rougher than is classically assumed. Deep in the glass, it undergoes a 'roughness transition' to fractal basins, which brings about isostaticity and marginal stability on approaching jamming. Critical exponents for the basin width, the weak force distribution and the spatial spread of quasi-contacts near jamming can be analytically determined. Their value is found to be compatible with numerical observations. This advance incorporates the jamming transition of granular materials into the framework of glass theory. Because temperature and pressure control what features of the landscape are experienced, glass mechanics and transport are expected to reflect the features of the topology we discuss here.

Exact theory of dense amorphous hard spheres in high dimension. II. The high density regime and the Gardner transition.

We consider the theory of the glass phase and jamming of hard spheres in the large space dimension limit. Building upon the exact expression for the free-energy functional obtained previously, we find that the random first order transition (RFOT) scenario is realized here with two thermodynamic transitions: the usual Kauzmann point associated with entropy crisis and a further transition at higher pressures in which a glassy structure of microstates is developed within each amorphous state. This kind of glass-glass transition into a phase dominating the higher densities was described years ago by Elisabeth Gardner, and may well be a generic feature of RFOT. Microstates that are small excitations of an amorphous matrix-separated by low entropic or energetic barriers-thus emerge naturally, and modify the high pressure (or low temperature) limit of the thermodynamic functions.