Multifractality in spin glasses.

We unveil the multifractal behavior of Ising spin glasses in their low-temperature phase. Using the Janus II custom-built supercomputer, the spin-glass correlation function is studied locally. Dramatic fluctuations are found when pairs of sites at the same distance are compared. The scaling of these fluctuations, as the spin-glass coherence length grows with time, is characterized through the computation of the singularity spectrum and its corresponding Legendre transform. A comparatively small number of site pairs controls the average correlation that governs the response to a magnetic field. We explain how this scenario of dramatic fluctuations (at length scales smaller than the coherence length) can be reconciled with the smooth, self-averaging behavior that has long been considered to describe spin-glass dynamics.

Strong ergodicity breaking in aging of mean-field spin glasses.

Out-of-equilibrium relaxation processes show aging if they become slower as time passes. Aging processes are ubiquitous and play a fundamental role in the physics of glasses and spin glasses and in other applications (e.g., in algorithms minimizing complex cost/loss functions). The theory of aging in the out-of-equilibrium dynamics of mean-field spin glass models has achieved a fundamental role, thanks to the asymptotic analytic solution found by Cugliandolo and Kurchan. However, this solution is based on assumptions (e.g., the weak ergodicity breaking hypothesis) which have never been put under a strong test until now. In the present work, we present the results of an extraordinary large set of numerical simulations of the prototypical mean-field spin glass models, namely the Sherrington-Kirkpatrick and the Viana-Bray models. Thanks to a very intensive use of graphics processing units (GPUs), we have been able to run the latter model for more than [Formula: see text] spin updates and thus safely extrapolate the numerical data both in the thermodynamical limit and in the large times limit. The measurements of the two-times correlation functions in isothermal aging after a quench from a random initial configuration to a temperature [Formula: see text] provides clear evidence that, at large times, such correlations do not decay to zero as expected by assuming weak ergodicity breaking. We conclude that strong ergodicity breaking takes place in mean-field spin glasses aging dynamics which, asymptotically, takes place in a confined configurational space. Theoretical models for the aging dynamics need to be revised accordingly.

Loop expansion around the Bethe solution for the random magnetic field Ising ferromagnets at zero temperature.

We apply to the random-field Ising model at zero temperature ([Formula: see text]) the perturbative loop expansion around the Bethe solution. A comparison with the standard ϵ expansion is made, highlighting the key differences that make the expansion around the Bethe solution much more appropriate to correctly describe strongly disordered systems, especially those controlled by a [Formula: see text] renormalization group (RG) fixed point. The latter loop expansion produces an effective theory with cubic vertices. We compute the one-loop corrections due to cubic vertices, finding additional terms that are absent in the ϵ expansion. However, these additional terms are subdominant with respect to the standard, supersymmetric ones; therefore, dimensional reduction is still valid at this order of the loop expansion.

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.

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.

The Mpemba effect in spin glasses is a persistent memory effect.

The Mpemba effect occurs when a hot system cools faster than an initially colder one, when both are refrigerated in the same thermal reservoir. Using the custom-built supercomputer Janus II, we study the Mpemba effect in spin glasses and show that it is a nonequilibrium process, governed by the coherence length ξ of the system. The effect occurs when the bath temperature lies in the glassy phase, but it is not necessary for the thermal protocol to cross the critical temperature. In fact, the Mpemba effect follows from a strong relationship between the internal energy and ξ that turns out to be a sure-tell sign of being in the glassy phase. Thus, the Mpemba effect presents itself as an intriguing avenue for the experimental study of the coherence length in supercooled liquids and other glass formers.

Long-Range Anomalous Decay of the Correlation in Jammed Packings.

We numerically study the structure of the interactions occurring in three-dimensional systems of hard spheres at jamming, focusing on the large-scale behavior. Given the fundamental role in the configuration of jammed packings, we analyze the propagation through the system of the weak forces and of the variation of the coordination number with respect to the isostaticity condition, ΔZ. We show that these correlations can be successfully probed by introducing a correlation function weighted on the density-density fluctuations. The results of this analysis can be further improved by introducing a representation of the system based on the contact points between particles. In particular, we find evidence that the weak forces and the ΔZ fluctuations support the hypothesis of randomly jammed packings of spherical particles being hyperuniform by exhibiting an anomalous long-range decay. Moreover, we find that the large-scale structure of the density-density correlation exhibits a complex behavior due to the superimposition of two exponentially damped oscillating signals propagating with linearly depending frequencies.

Inferring the particle-wise dynamics of amorphous solids from the local structure at the jamming point.

Jamming is a phenomenon shared by a wide variety of systems, such as granular materials, foams, and glasses in their high density regime. This has motivated the development of a theoretical framework capable of explaining many of their static critical properties with a unified approach. However, the dynamics occurring in the vicinity of the jamming point has received little attention and the problem of finding a connection with the local structure of the configuration remains unexplored. Here we address this issue by constructing physically well defined structural variables using the information contained in the network of contacts of jammed configurations, and then showing that such variables yield a resilient statistical description of the particle-wise dynamics near this critical point. Our results are based on extensive numerical simulations of systems of spherical particles that allow us to statistically characterize the trajectories of individual particles in terms of their first two moments. We first demonstrate that, besides displaying a broad distribution of mobilities, particles may also have preferential directions of motion. Next, we associate each of these features with a structural variable computed uniquely in terms of the contact vectors at jamming, obtaining considerably high statistical correlations. The robustness of our approach is confirmed by testing two types of dynamical protocols, namely molecular dynamics and Monte Carlo, with different types of interaction. We also provide evidence that the dynamical regime we study here is dominated by anharmonic effects and therefore it cannot be described properly in terms of vibrational modes. Finally, we show that correlations decay slowly and in an interaction-independent fashion, suggesting a universal rate of information loss.

Support for the value 5/2 for the spin glass lower critical dimension at zero magnetic field.

We study numerically various properties of the free energy barriers in the Edwards-Anderson model of spin glasses in the low-temperature region in both three and four spatial dimensions. In particular, we investigated the dependence of height of free energy barriers on system size and on the distance between the initial and final states (i.e., the overlap distance). A related quantity is the distribution of large local fluctuations of the overlap in large 3D samples at equilibrium. Our results for both quantities (barriers and large deviations) are in agreement with the prediction obtained in the framework of mean-field theory. In addition, our result supports [Formula: see text] as the lower critical dimension of the model.

Probing the non-Debye low-frequency excitations in glasses through random pinning.

We investigate the properties of the low-frequency spectrum in the density of states [Formula: see text] of a 3D model glass former. To magnify the non-Debye sector of the spectrum, we introduce a random pinning field that freezes a finite particle fraction to break the translational invariance and shifts all of the vibrational frequencies of the extended modes toward higher frequencies. We show that non-Debye soft localized modes progressively emerge as the fraction p of pinned particles increases. Moreover, the low-frequency tail of [Formula: see text] goes to zero as a power law [Formula: see text], with [Formula: see text] and [Formula: see text] above a threshold fraction [Formula: see text].

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.

A statics-dynamics equivalence through the fluctuation-dissipation ratio provides a window into the spin-glass phase from nonequilibrium measurements.

We have performed a very accurate computation of the nonequilibrium fluctuation-dissipation ratio for the 3D Edwards-Anderson Ising spin glass, by means of large-scale simulations on the special-purpose computers Janus and Janus II. This ratio (computed for finite times on very large, effectively infinite, systems) is compared with the equilibrium probability distribution of the spin overlap for finite sizes. Our main result is a quantitative statics-dynamics dictionary, which could allow the experimental exploration of important features of the spin-glass phase without requiring uncontrollable extrapolations to infinite times or system sizes.

Spin Glasses in a Field Show a Phase Transition Varying the Distance among Real Replicas (And How to Exploit It to Find the Critical Line in a Field).

We discuss a phase transition in spin glass models that have been rarely considered in the past, namely, the phase transition that may take place when two real replicas are forced to be at a larger distance (i.e., at a smaller overlap) than the typical one. In the first part of the work, by solving analytically the Sherrington-Kirkpatrick model in a field close to its critical point, we show that, even in a paramagnetic phase, the forcing of two real replicas to an overlap small enough leads the model to a phase transition where the symmetry between replicas is spontaneously broken. More importantly, this phase transition is related to the de Almeida-Thouless (dAT) critical line. In the second part of the work, we exploit the phase transition in the overlap between two real replicas to identify the critical line in a field in finite dimensional spin glasses. This is a notoriously difficult computational problem, because of considerable finite size corrections. We introduce a new method of analysis of Monte Carlo data for disordered systems, where the overlap between two real replicas is used as a conditioning variate. We apply this analysis to equilibrium measurements collected in the paramagnetic phase in a field, h > 0 and T c ( h ) < T < T c ( h = 0 ) , of the d = 1 spin glass model with long range interactions decaying fast enough to be outside the regime of validity of the mean field theory. We thus provide very reliable estimates for the thermodynamic critical temperature in a field.

Relation between Heterogeneous Frozen Regions in Supercooled Liquids and Non-Debye Spectrum in the Corresponding Glasses.

Recent numerical studies on glassy systems provide evidence for a population of non-Goldstone modes (NGMs) in the low-frequency spectrum of the vibrational density of states D(ω). Similarly to Goldstone modes (GMs), i.e., phonons in solids, NGMs are soft low-energy excitations. However, differently from GMs, NGMs are localized excitations. Here we first show that the parental temperature T^{*} modifies the GM/NGM ratio in D(ω). In particular, the phonon attenuation is reflected in a parental temperature dependency of the exponent s(T^{*}) in the low-frequency power law D(ω)∼ω^{s(T^{*})}, with 2≤s(T^{*})≤4. Second, by comparing s(T^{*}) with s(p), i.e., the same quantity obtained by pinning a p particle fraction, we suggest that s(T^{*}) reflects the presence of dynamical heterogeneous regions of size ξ^{3}∝p. Finally, we provide an estimate of ξ as a function of T^{*}, finding a mild power law divergence, ξ∼(T^{*}-T_{d})^{-α/3}, with T_{d} the dynamical crossover temperature and α falling in the range α∈[0.8,1.0].

Evidence for Supersymmetry in the Random-Field Ising Model at D=5.

We provide a nontrivial test of supersymmetry in the random-field Ising model at five spatial dimensions, by means of extensive zero-temperature numerical simulations. Indeed, supersymmetry relates correlation functions in a D-dimensional disordered system with some other correlation functions in a D-2 clean system. We first show how to check these relationships in a finite-size scaling calculation and then perform a high-accuracy test. While the supersymmetric predictions are satisfied even to our high accuracy at D=5, they fail to describe our results at D=4.

Growing timescales and lengthscales characterizing vibrations of amorphous solids.

Low-temperature properties of crystalline solids can be understood using harmonic perturbations around a perfect lattice, as in Debye's theory. Low-temperature properties of amorphous solids, however, strongly depart from such descriptions, displaying enhanced transport, activated slow dynamics across energy barriers, excess vibrational modes with respect to Debye's theory (i.e., a boson peak), and complex irreversible responses to small mechanical deformations. These experimental observations indirectly suggest that the dynamics of amorphous solids becomes anomalous at low temperatures. Here, we present direct numerical evidence that vibrations change nature at a well-defined location deep inside the glass phase of a simple glass former. We provide a real-space description of this transition and of the rapidly growing time- and lengthscales that accompany it. Our results provide the seed for a universal understanding of low-temperature glass anomalies within the theoretical framework of the recently discovered Gardner phase transition.

Configurational entropy of polydisperse supercooled liquids.

We propose a computational method to measure the configurational entropy in generic polydisperse glass-formers. In particular, our method resolves issues related to the diverging mixing entropy term due to a continuous polydispersity. The configurational entropy is measured as the difference between the well-defined fluid entropy and a more problematic glass entropy. We show that the glass entropy can be computed by a simple generalisation of the Frenkel-Ladd thermodynamic integration method, which takes into account permutations of the particle diameters. This approach automatically provides a physically meaningful mixing entropy for the glass entropy and includes contributions that are not purely vibrational. The proposed configurational entropy is thus devoid of conceptual and technical difficulties due to continuous polydispersity, while being conceptually closer, but technically simpler, than alternative free energy approaches.

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.

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.