Nontrivial Critical Fixed Point for Replica-Symmetry-Breaking Transitions

Evidence of a second-order phase transition in the six-dimensional Ising spin glass in a field

The very existence of a phase transition for spin glasses in an external magnetic field is controversial, even in high dimensions. We carry out massive simulations of the Ising spin-glass in a field, in six dimensions (which, according to classical-but not generally accepted-field-theoretical studies, is the upper critical dimension). We obtain results compatible with a second-order phase transition and estimate its critical exponents for the simulated lattice sizes. The detailed analysis performed by other authors of the replica symmetric Hamiltonian, under the hypothesis of critical behavior, predicts that the ratio of the renormalized coupling constants remain bounded as the correlation length grows. Our numerical results are in agreement with this expectation.

Dynamic Gardner cross-over in a simple glass

The criticality of the jamming transition responsible for amorphous solidification has been theoretically linked to the marginal stability of a thermodynamic Gardner phase. While the critical exponents of jamming appear independent of the preparation history, the pertinence of Gardner physics far from equilibrium is an open question. To fill this gap, we numerically study the nonequilibrium dynamics of hard disks compressed toward the jamming transition using a broad variety of protocols. We show that dynamic signatures of Gardner physics can be disentangled from the aging relaxation dynamics. We thus define a generic dynamic Gardner cross-over regardless of the history. Our results show that the jamming transition is always accessed by exploring increasingly complex landscape, resulting in anomalous microscopic relaxation dynamics that remains to be understood theoretically.

Machine learning glass caging order parameters with an artificial nested neural network

Around a glass transition, the dynamics of a supercooled liquid dramatically slow down, exhibited by caging of particles, while the structural changes remain subtle. Alternative to recent machine learning studies searching for structural predictors of glassy dynamics, here we propose to learn directly particle caging features defined purely according to dynamics. We focus on three transitions in a simulated hard sphere glass model, the melting of ultra-stable glasses, the Gardner transition and the liquid to ordinary glass transition. Implementing the machine learning algorithm based on a two-level nested neural network, we attain not only appropriate caging order parameters for all three transitions, but also a phase classification for input samples. A finite-size scaling analysis of the phase classification results identifies the order of melting (first) and Gardner (second) transitions. A false positive is avoided, as the liquid to glass transition is indicated as a crossover, rather than a phase transition with a well-defined transition point. This study paves the way to a generic approach for learning dynamical features in glassy systems, with a minimum requirement of system-specific knowledge.

Theory and experiments for disordered elastic manifolds, depinning, avalanches, and sandpiles

Domain walls in magnets, vortex lattices in superconductors, contact lines at depinning, and many other systems can be modeled as an elastic system subject to quenched disorder. The ensuing field theory possesses a well-controlled perturbative expansion around its upper critical dimension. Contrary to standard field theory, the renormalization group (RG) flow involves a function, the disorder correlator Δ(w), and is therefore termed the functional RG. Δ(w) is a physical observable, the auto-correlation function of the center of mass of the elastic manifold. In this review, we give a pedagogical introduction into its phenomenology and techniques. This allows us to treat both equilibrium (statics), and depinning (dynamics). Building on these techniques, avalanche observables are accessible: distributions of size, duration, and velocity, as well as the spatial and temporal shape. Various equivalences between disordered elastic manifolds, and sandpile models exist: an elastic string driven at a point and the Oslo model; disordered elastic manifolds and Manna sandpiles; charge density waves and Abelian sandpiles or loop-erased random walks. Each of the mappings between these systems requires specific techniques, which we develop, including modeling of discrete stochastic systems via coarse-grained stochastic equations of motion, super-symmetry techniques, and cellular automata. Stronger than quadratic nearest-neighbor interactions lead to directed percolation, and non-linear surface growth with additional Kardar-Parisi-Zhang (KPZ) terms. On the other hand, KPZ without disorder can be mapped back to disordered elastic manifolds, either on the directed polymer for its steady state, or a single particle for its decay. Other topics covered are the relation between functional RG and replica symmetry breaking, and random-field magnets. Emphasis is given to numerical and experimental tests of the theory.

Probing Gardner Physics in an Active Quasithermal Pressure-Controlled Granular System of Noncircular Particles

To search for experimental signals of the Gardner crossover, an active quasithermal granular glass is constructed using a monolayer of air-fluidized star-shaped particles. The pressure of the system is controlled by adjusting the tension exerted on an enclosing boundary. Velocity distributions of the internal particles and the scaling of the pressure, density, effective temperature, and relaxation time are examined, demonstrating that the system has key features of a thermal system. Using a pressure-based quenching protocol that brings the system into deeper glassy states, signals of the Gardner crossover are detected via cage size and separation order parameters for both particle positions and orientations, offering experimental evidence of Gardner physics for a system of anisotropic quasithermal particles in a low spatial dimension.

Numerical test of the replica-symmetric Hamiltonian for correlations of the critical state of spin glasses in a field

A growing body of evidence indicates that the sluggish low-temperature dynamics of glass formers (e.g., supercooled liquids, colloids, or spin glasses) is due to a growing correlation length. Which is the effective field theory that describes these correlations? The natural field theory was drastically simplified by Bray and Roberts in 1980. More than 40 years later, we confirm the tenets of Bray and Roberts's theory by studying the Ising spin glass in an externally applied magnetic field, both in four spatial dimensions (data obtained from the Janus collaboration) and on the Bethe lattice.

Unexpected Upper Critical Dimension for Spin Glass Models in a Field Predicted by the Loop Expansion around the Bethe Solution at Zero Temperature

The spin-glass transition in a field in finite dimension is analyzed directly at zero temperature using a perturbative loop expansion around the Bethe lattice solution. The loop expansion is generated by the M-layer construction whose first diagrams are evaluated numerically and analytically. The generalized Ginzburg criterion reveals that the upper critical dimension below which mean-field theory fails is D_{U}≥8, at variance with the classical result D_{U}=6 yielded by finite-temperature replica field theory. Our expansion around the Bethe lattice has two crucial differences with respect to the classical one. The finite connectivity z of the lattice is directly included from the beginning in the Bethe lattice, while in the classical computation the finite connectivity is obtained through an expansion in 1/z. Moreover, if one is interested in the zero temperature (T=0) transition, one can directly expand around the T=0 Bethe transition. The expansion directly at T=0 is not possible in the classical framework because the fully connected spin glass does not have a transition at T=0, being in the broken phase for any value of the external field.

Determining the nonequilibrium criticality of a Gardner transition via a hybrid study of molecular simulations and machine learning

Apparent critical phenomena, typically indicated by growing correlation lengths and dynamical slowing down, are ubiquitous in nonequilibrium systems such as supercooled liquids, amorphous solids, active matter, and spin glasses. It is often challenging to determine if such observations are related to a true second-order phase transition as in the equilibrium case or simply a crossover and even more so to measure the associated critical exponents. Here we show that the simulation results of a hard-sphere glass in three dimensions are consistent with the recent theoretical prediction of a Gardner transition, a continuous nonequilibrium phase transition. Using a hybrid molecular simulation-machine learning approach, we obtain scaling laws for both finite-size and aging effects and determine the critical exponents that traditional methods fail to estimate. Our study provides an approach that is useful to understand the nature of glass transitions and can be generalized to analyze other nonequilibrium phase transitions.

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.

Exploratory study of the glassy landscape near jamming

We present a study of the landscape structure of hard and soft spheres as a function of the packing fraction and of the energy. We find that, on approaching the jamming transition, the number of different configurations available to the system increases steeply and a hierarchical organization of the landscape emerges. We use the knowledge of the structure of the landscape to predict the values of thermodynamic observables on the edge of the transition.

One-step replica-symmetry-breaking phase below the de Almeida-Thouless line in low-dimensional spin glasses

The de Almeida-Thouless (AT) line is the phase boundary in the temperature-magnetic field plane of an Ising spin glass at which a continuous (i.e., second-order) transition from a paramagnet to a replica-symmetry-breaking (RSB) phase occurs, according to mean-field theory. Here, using field-theoretic perturbative renormalization group methods on the Bray-Roberts reduced Landau-Ginzburg-type theory for a short-range Ising spin glass in space of dimension d, we show that at nonzero magnetic field the nature of the corresponding transition is modified as follows: (a) For d-6 small and positive, with increasing field on the AT line, first, the ordered phase just below the transition becomes the so-called one-step RSB, instead of the full RSB that occurs in mean-field theory; the transition on the AT line remains continuous with a diverging correlation length. Then at a higher field, a tricritical point separates the latter transition from a quasi-first-order one, that is one at which the correlation length does not diverge, and there is a jump in part of the order parameter, but no latent heat. The location of the tricritical point tends to zero as d→6^{+}. (b) For d≤6, we argue that the quasi-first-order transition could persist down to arbitrarily small nonzero fields, with a transition to full RSB still expected at lower temperature. Whenever the quasi-first-order transition occurs, it is at a higher temperature than the AT transition would be for the same field, preempting it as the temperature is lowered. These results may explain the reported absence of a diverging correlation length in the presence of a magnetic field in low-dimensional spin glasses in some simulations and in high-temperature series expansions. We also draw attention to the similarity of the "dynamically frozen" state, which occurs at temperatures just above the quasi-first-order transition, and the "metastate-average state" of the one-step RSB phase, and discuss the issue of the number of pure states in either.

Gardner physics in amorphous solids and beyond

One of the most remarkable predictions to emerge out of the exact infinite-dimensional solution of the glass problem is the Gardner transition. Although this transition was first theoretically proposed a generation ago for certain mean-field spin glass models, its materials relevance was only realized when a systematic effort to relate glass formation and jamming was undertaken. A number of nontrivial physical signatures associated with the Gardner transition have since been considered in various areas, from models of structural glasses to constraint satisfaction problems. This perspective surveys these recent advances and discusses the novel research opportunities that arise from them.

Rejuvenation and Memory Effects in a Structural Glass

We show numerically that a three-dimensional model for structural glass displays aging, rejuvenation, and memory effects when subjected to a temperature cycle. These effects indicate that the free energy landscape of structural glasses may possess the complex hierarchical structure that characterizes materials such as spin and polymer glasses. We use the theoretical concept of marginal stability to interpret our results, and explain in which physical conditions a complex aging dynamics can emerge in dense supercooled liquids, paving the way for future experimental studies of complex aging dynamics in colloidal and granular glasses.

Impact of jamming criticality on low-temperature anomalies in structural glasses

We present a mechanism for the anomalous behavior of the specific heat in low-temperature amorphous solids. The analytic solution of a mean-field model belonging to the same universality class as high-dimensional glasses, the spherical perceptron, suggests that there exists a cross-over temperature above which the specific heat scales linearly with temperature, while below it, a cubic scaling is displayed. This relies on two crucial features of the phase diagram: (i) the marginal stability of the free-energy landscape, which induces a gapless phase responsible for the emergence of a power-law scaling; and (ii) the vicinity of the classical jamming critical point, as the cross-over temperature gets lowered when approaching it. This scenario arises from a direct study of the thermodynamics of the system in the quantum regime, where we show that, contrary to crystals, the Debye approximation does not hold.

Morphology of renormalization-group flow for the de Almeida-Thouless-Gardner universality class

A replica-symmetry-breaking phase transition is predicted in a host of disordered media. The criticality of the transition has, however, long been questioned below its upper critical dimension, six, due to the absence of a critical fixed point in the renormalization-group flows at one-loop order. A recent two-loop analysis revealed a possible strong-coupling fixed point, but given the uncontrolled nature of perturbative analysis in the strong-coupling regime, debate persists. Here we examine the nature of the transition as a function of spatial dimension and show that the strong-coupling fixed point can go through a Hopf bifurcation, resulting in a critical limit cycle and a concomitant discrete scale invariance. We further investigate a different renormalization scheme and argue that the basin of attraction of the strong-coupling fixed point (or limit cycle) may stay finite for all dimensions.

Anomalous nonlinear damping in metallic glasses: Signature of elasticity breakdown

Solid materials, whether crystalline or glasses, are characterized by their elasticity. Generally, elasticity is independent of the probing strain if it is not exceeding the yielding point. Here, by contrast, we experimentally capture a pronounced strain-dependent elasticity in metallic glasses, as manifested by nonlinear mechanical damping in the apparent elastic deformation regime (∼1/100 of the yielding strain). Normal damping behaviors recover at higher temperatures but still below the glass transition. Atomistic simulations reproduce these features and reveal that they could be related to avalanche-like local structural instabilities. Our findings demonstrate that the standard elasticity is not held for metallic glasses at low temperatures and plastic events can be triggered at small perturbations. These results are consistent with previous simulations of model glasses and a scenario of hierarchical free-energy landscape of mean-field theory.

Marginally stable phases in mean-field structural glasses

A novel form of amorphous matter characterized by marginal stability was recently discovered in the mean-field theory of structural glasses. Using this approach, we provide complete phase diagrams delimiting the location of the marginally stable glass phase for a large variety of pair interactions and physical conditions, extensively exploring physical regimes relevant to granular matter, foams, emulsions, hard and soft colloids, and molecular glasses. We find that all types of glasses may become marginally stable, but the extent of the marginally stable phase highly depends on the preparation protocol. Our results suggest that marginal phases should be observable for colloidal and non-Brownian particles near jamming and for poorly annealed glasses. For well-annealed glasses, two distinct marginal phases are predicted. Our study unifies previous results on marginal stability in mean-field models and will be useful to guide numerical simulations and experiments aimed at detecting marginal stability in finite-dimensional amorphous materials.

Glassy, Gardner-like phenomenology in minimally polydisperse crystalline systems

We report on a nonequilibrium phase of matter, the minimally disordered crystal phase, which we find exists between the maximally amorphous glasses and the ideal crystal. Even though these near crystals appear highly ordered, they display glassy and jamming features akin to those observed in amorphous solids. Structurally, they exhibit a power-law scaling in their probability distribution of weak forces and small interparticle gaps as well as a flat density of vibrational states. Dynamically, they display anomalous aging above a characteristic pressure. Quantitatively, this disordered crystal phase has much in common with the Gardner-like phase seen in maximally disordered solids. Near crystals should be amenable to experimental realizations in commercially available particulate systems and are to be indispensable in verifying the theory of amorphous materials.

Identity crisis in alchemical space drives the entropic colloidal glass transition

A universally accepted explanation for why liquids sometimes vitrify rather than crystallize remains hotly pursued, despite the ubiquity of glass in our everyday lives, the utilization of the glass transition in innumerable modern technologies, and nearly a century of theoretical and experimental investigation. Among the most compelling hypothesized mechanisms underlying glass formation is the development in the fluid phase of local structures that somehow prevent crystallization. Here, we explore that mechanism in the case of hard particle glasses by examining the glass transition in an extended alchemical (here, shape) space; that is, a space where particle shape is treated as a thermodynamic variable. We investigate simple systems of hard polyhedra, with no interactions aside from volume exclusion, and show via Monte Carlo simulation that glass formation in these systems arises from a multiplicity of competing local motifs, each of which is prevalent in-and predictable from-nearby ordered structures in alchemical space.

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.

Spin-glass-like aging in colloidal and granular glasses

Motivated by the mean field prediction of a Gardner phase transition between a "normal glass" and a "marginally stable glass", we investigate the off-equilibrium dynamics of three-dimensional polydisperse hard spheres, used as a model for colloidal or granular glasses. Deep inside the glass phase, we find that a sharp crossover pressure PG separates two distinct dynamical regimes. For pressure P < PG, the glass behaves as a normal solid, displaying fast dynamics that quickly equilibrate within the glass free energy basin. For P > PG, instead, the dynamics become strongly anomalous, displaying very large equilibration timescales, aging, and a constantly increasing dynamical susceptibility. The crossover at PG is strongly reminiscent of the one observed in three-dimensional spin-glasses in an external field, suggesting that the two systems could be in the same universality class, consistent with theoretical expectations.

Multicritical Point on the de Almeida-Thouless Line in Spin Glasses in d>6 Dimensions

The de Almeida-Thouless (AT) line in Ising spin glasses is the phase boundary in the temperature T and magnetic field h plane below which replica symmetry is broken. Using perturbative renormalization group (RG) methods, we show that, when the dimension d of space is just above six, there is a multicritical point (MCP) on the AT line, which separates a low-field regime, in which the critical exponents have mean-field values, from a high-field regime, where the RG flows run away to infinite coupling strength; as d approaches six from above, the MCP approaches the zero-field critical point exponentially in 1/(d-6). Thus, on the AT line, perturbation theory for the critical properties breaks down at a sufficiently large magnetic field even above 6 dimensions, as well as for all nonzero fields when d≤6, as was known previously. We calculate the exponents at the MCP to first order in ϵ=d-6>0. The fate of the MCP as d increases from just above six to infinity is not known.

Cycle-expansion method for the Lyapunov exponent, susceptibility, and higher moments

Lyapunov exponents characterize the chaotic nature of dynamical systems by quantifying the growth rate of uncertainty associated with the imperfect measurement of initial conditions. Finite-time estimates of the exponent, however, experience fluctuations due to both the initial condition and the stochastic nature of the dynamical path. The scale of these fluctuations is governed by the Lyapunov susceptibility, the finiteness of which typically provides a sufficient condition for the law of large numbers to apply. Here, we obtain a formally exact expression for this susceptibility in terms of the Ruelle dynamical ζ function for one-dimensional systems. We further show that, for systems governed by sequences of random matrices, the cycle expansion of the ζ function enables systematic computations of the Lyapunov susceptibility and of its higher-moment generalizations. The method is here applied to a class of dynamical models that maps to static disordered spin chains with interactions stretching over a varying distance and is tested against Monte Carlo simulations.

Absence of Marginal Stability in a Structural Glass

Marginally stable solids have peculiar physical properties that were first analyzed in the context of the jamming transition. We theoretically investigate the existence of marginal stability in a prototypical model for structural glass formers, combining analytical calculations in infinite dimensions to computer simulations in three dimensions. While mean-field theory predicts the existence of a Gardner phase transition towards a marginally stable glass phase at low temperatures, simulations show no hint of diverging time scales or length scales, but reveal instead the presence of sparse localized defects. Our results suggest that the Gardner transition is deeply affected by finite dimensional fluctuations, and raise issues about the relevance of marginal stability in structural glasses far away from jamming.