Search

We study numerically the Hessian of low-lying minima of vector spin glass models defined on random regular graphs. We consider the two-component (XY) and three-component (Heisenberg) spin glasses at zero temperature, subjected to the action of a randomly oriented external field. Varying the intensity of the external field, these models undergo a zero temperature phase transition from a paramagnet at high field to a spin glass at low field. We study how the spectral properties of the Hessian depend on the magnetic field. In particular, we study the shape of the spectrum at low frequency and the localization properties of low energy eigenvectors across the transition. We find that in both phases the edge of the spectral density behaves as $\lambda^{3/2}$: such a behavior rules out the presence of a diverging spin-glass susceptibility $\chi_{SG}=\langle 1/\lambda^2 \rangle$. As to low energy eigenvectors, we find that the softest eigenmodes are always localized in both phases of the two models. However, by studying in detail the geometry of low energy eigenmodes across different energy scales close to the lower edge of the spectrum, we find a different behavior for the two models at the transition: in the XY case, low energy modes are typically localized; at variance, in the Heisenberg case low-energy eigenmodes with a multi-modal structure (sort of ``delocalization'') appear at an energy scale that vanishes in the infinite size limit. These geometrically non-trivial excitations, which we call Concentrated and Delocalised Low Energy Modes (CDLEM), coexist with trivially localised excitations: we interpret their existence as a sign of critical behavior related to the onset of the spin glass phase., Comment: 15 pages, 8 figures

We study the chaotic behavior of the Gibbs state of spin-glasses under the application of an external magnetic field, in the crossover region where the field intensity scales proportional to $1/\sqrt{N}$, being $N$ the system size. We show that Replica Symmetry Breaking (RSB) theory provides universal predictions for chaotic behavior: they depend only on the zero-field overlap probability function $P(q)$ and are independent of other features of the system. Using solely $P(q)$ as input we can analytically predict quantitatively the statistics of the states in a small field. In the infinite volume limit, each spin-glass sample is characterized by an infinite number of states that have a tree-like structure. We generate the corresponding probability distribution through efficient sampling using a representation based on the Bolthausen-Sznitman coalescent. In this way, we can compute quantitatively properties in the presence of a magnetic field in the crossover region, the overlap probability distribution in the presence of a small field and the degree of decorrelation as the field is increased. To test our computations, we have simulated the Bethe lattice spin glass and the 4D Edwards-Anderson model, finding in both cases excellent agreement with the universal predictions., Comment: New version with modifications included after the review process and publication. 17 pages, 13 figures

We release a set of GPU programs for the study of the Quantum ($S=1/2$) Spin Glass on a square lattice, with binary couplings. The library contains two main codes: MCQSG (that carries out Monte Carlo simulations using both the Metropolis and the Parallel Tempering algorithms, for the problem formulated in the Trotter-Suzuki approximation), and EDQSG (that obtains the extremal eigenvalues of the Transfer Matrix using the Lanczos algorithm). EDQSG has allowed us to diagonalize transfer matrices with size up to $2^{36}\times2^{36}$. From its side, MCQSG running on four NVIDIA A100 cards delivers a sub-picosecond time per spin-update, a performance that is competitive with dedicated hardware. We include as well in our library GPU programs for the analysis of the spin configurations generated by MCQSG. Finally, we provide two auxiliary codes: the first generates the lookup tables employed by the random number generator of MCQSG; the second one simplifies the execution of multiple runs using different input data., Comment: Code available in the Computer Physics Communications repository

Quantum annealers are commercial devices aiming to solve very hard computational problems named spin glasses. Just like in metallurgic annealing one slowly cools a ferrous metal, quantum annealers seek good solutions by slowly removing the transverse magnetic field at the lowest possible temperature. The field removal diminishes quantum fluctuations but forces the system to traverse the critical point that separates the disordered phase (at large fields) from the spin-glass phase (at small fields). A full understanding of this phase transition is still missing. A debated, crucial question regards the closing of the energy gap separating the ground state from the first excited state. All hopes of achieving an exponential speed-up, as compared to classical computers, rest on the assumption that the gap will close algebraically with the number of qspins, but renormalization group calculations predict that the closing will be instead exponential. Here we solve this debate through extreme-scale numerical simulations, finding that both parties grasped parts of the truth. While the closing of the gap at the critical point is indeed super-algebraic, it remains algebraic if one restricts the symmetry of possible excitations. Since this symmetry restriction is experimentally achievable (at least nominally), there is still hope for the Quantum Annealing paradigm.

We have studied numerically the appearance of multiscaling behavior in the three-dimensional ferromagnetic Ising site diluted model, in the form of a multifractal distribution of the decay exponents for the spatial correlation functions at the critical temperature. We have computed the exponents of the long-distance decay of higher moments of the correlation function, up to the 10th power, by studying three different quantities: global susceptibilities, local susceptibilities and correlation functions. We have found very clear evidences for multiscaling behavior., Comment: 19 pages and 5 figures. Final version of the paper

Rejuvenation and memory, long considered the distinguishing features of spin glasses, have recently been proven to result from the growth of multiple length scales. This insight, enabled by simulations on the Janus~II supercomputer, has opened the door to a quantitative analysis. We combine numerical simulations with comparable experiments to introduce two coefficients that quantify memory. A third coefficient has been recently presented by Freedberg et al. We show that these coefficients are physically equivalent by studying their temperature and waiting-time dependence., Comment: 10 pages, 8 figures

Finite-size scaling above the upper critical dimension is a long-standing puzzle in the field of Statistical Physics. Even for pure systems various scaling theories have been suggested, partially corroborated by numerical simulations. In the present manuscript we address this problem in the even more complicated case of disordered systems. In particular, we investigate the scaling behavior of the random-field Ising model at dimension $D = 7$, i.e., above its upper critical dimension $D_{\rm u} = 6$, by employing extensive ground-state numerical simulations. Our results confirm the hypothesis that at dimensions $D > D_{\rm u}$, linear length scale $L$ should be replaced in finite-size scaling expressions by the effective scale $L_{\rm eff} = L^{D / D_{\rm u}}$. Via a fitted version of the quotients method that takes this modification, but also subleading scaling corrections into account, we compute the critical point of the transition for Gaussian random fields and provide estimates for the full set of critical exponents. Thus, our analysis indicates that this modified version of finite-size scaling is successful also in the context of the random-field problem., Comment: 10 pages, 5 figures, Appendix included, version published in Phys. Rev. E

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., Comment: 15 pages, 15 figures

This is an extended version of my Nobel Lecture, delivered on December $8^{th}$ 2021. I will recall the genesis of the concept of multiple equilibria in natural sciences. I will then describe my contribution to the development of this concept in the framework of statistical mechanics. Finally, I will briefly mention the cornucopia of applications of these ideas both in physics and in other disciplines., Comment: 18 pages, 7 figures, a few typos corrected and very minor corrections added

By means of enhanced Monte Carlo numerical simulations parallelized on GPU's we study the critical properties of the spin-glass-like model for the mode-locked glassy random laser, a $4$-spin model with complex spins with a global spherical constraint and quenched random interactions. Using two different boundary conditions for the mode frequencies we identify the critical points and the critical indices of the random lasing phase transition using , with finite size scaling techniques. The outcome of the scaling analysis is that the mode-locked random laser is in a mean-field universality class, though different from the mean-field class of the Random Energy Model and the glassy random laser in the narrow band approximation, that is, the fully connected version of the present model. The low temperature (high pumping) phase is finally characterized by means of the overlap distribution and evidence for the onset of replica symmetry breaking in the lasing regime is provided., Comment: 17 pages, 11 figures

The extended principle of superposition has been a touchstone of spin glass dynamics for almost thirty years. The Uppsala group has demonstrated its validity for the metallic spin glass, CuMn, for magnetic fields $H$ up to 10 Oe at the reduced temperature $T_\mathrm{r}=T/T_\mathrm{g} = 0.95$, where $T_\mathrm{g}$ is the spin glass condensation temperature. For $H > 10$ Oe, they observe a departure from linear response which they ascribe to the development of non-linear dynamics. The thrust of this paper is to develop a microscopic origin for this behavior by focusing on the time development of the spin glass correlation length, $\xi(t,t_\mathrm{w};H)$. Here, $t$ is the time after $H$ changes, and $t_\mathrm{w}$ is the time from the quench for $T>T_\mathrm{g}$ to the working temperature $T$ until $H$ changes. We connect the growth of $\xi(t,t_\mathrm{w};H)$ to the barrier heights $\Delta(t_\mathrm{w})$ that set the dynamics. The effect of $H$ on the magnitude of $\Delta(t_\mathrm{w})$ is responsible for affecting differently the two dynamical protocols associated with turning $H$ off (TRM, or thermoremanent magnetization) or on (ZFC, or zero field-cooled magnetization). In this paper, we display the difference between the zero-field cooled $\xi_{\text {ZFC}}(t,t_\mathrm{w};H)$ and the thermoremanent magnetization $\xi_{\text {TRM}}(t,t_\mathrm{w};H)$ correlation lengths as $H$ increases, both experimentally and through numerical simulations, corresponding to the violation of the extended principle of superposition in line with the finding of the Uppsala Group., Comment: 23 pages and 18 figures

Memory and rejuvenation effects in the magnetic response of off-equilibrium spin glasses have been widely regarded as the doorway into the experimental exploration of ultrametricity and temperature chaos (maybe the most exotic features in glassy free-energy landscapes). Unfortunately, despite more than twenty years of theoretical efforts following the experimental discovery of memory and rejuvenation, these effects have thus far been impossible to simulate reliably. Yet, three recent developments convinced us to accept this challenge: first, the custom-built Janus II supercomputer makes it possible to carry out "numerical experiments" in which the very same quantities that can be measured in single crystals of CuMn are computed from the simulation, allowing for parallel analysis of the simulation/experiment data. Second, Janus II simulations have taught us how numerical and experimental length scales should be compared. Third, we have recently understood how temperature chaos materializes in aging dynamics. All three aspects have proved crucial for reliably reproducing rejuvenation and memory effects on the computer. Our analysis shows that (at least) three different length scales play a key role in aging dynamics, while essentially all theoretical analyses of the aging dynamics emphasize the presence and the crucial role of a single glassy correlation length., Comment: Main text:13 pages, 4 figures SM: 7 pages, 8 figures

We investigate the replica symmetry broken (RSB) phase of spin glass (SG) models in a random field defined on Bethe lattices at zero temperature. From the properties of the RSB solution we deduce a closed equation for the extreme values of the cavity fields. This equation turns out not to depend on the parameters defining the RSB, and it predicts that the spontaneous RSB does not take place homogeneously on the whole system. Indeed, there exist spins having the same effective local field in all local ground states, exactly as in the replica symmetric (RS) phase, while the spontaneous RSB manifests only on the remaining spins, whose fraction vanishes at criticality. The characterization in terms of spins having fixed or fluctuating local fields can be extended also to the random field Ising model (RFIM), in which case the fluctuating spins are the only responsible for the spontaneous magnetization in the ferromagnetic phase. Close to criticality we are able to connect the statistics of the local fields acting on the spins in the RSB phase with the correlation functions measured in the paramagnetic phase. Identifying the two types of spins on given instances of SG and RFIM, we show that they participate very differently to avalanches produced by flipping a single spin. From the scaling of the number of spins inducing RSB effects close to the critical point and using the $M$-layer expansion we estimate the upper critical dimension $D_U \geq 8$ for SG., Comment: 8 pages, 6 figures

The jamming transition is ubiquitous. It is present in granular matter, colloids, glasses, and many other systems. Yet, it defines a critical point whose properties still need to be fully understood. A major breakthrough came about when the replica formalism was extended to build a mean-field theory that provides an exact description of the jamming transition of spherical particles in the infinite-dimensional limit. While such theory explains the jamming critical behavior of both soft and hard spheres, investigating the transition in finite-dimensional systems poses very difficult and different problems, in particular from the numerical point of view. Soft particles are modeled by continuous potentials; thus, their jamming point can be reached through efficient energy minimization algorithms. In contrast, the latter methods are inapplicable to hard-sphere (HS) systems since the interaction energy among the particles is always zero by construction. To overcome these difficulties, here we recast the jamming of hard spheres as a constrained optimization problem and introduce the CALiPPSO algorithm, capable of readily producing jammed HS packings without including any effective potential. This algorithm brings a HS configuration of arbitrary dimensions to its jamming point by solving a chain of linear optimization problems. We show that there is a strict correspondence between the force balance conditions of jammed packings and the properties of the optimal solutions of CALiPPSO, whence we prove analytically that our packings are always isostatic and in mechanical equilibrium. Furthermore, using extensive numerical simulations, we show that our algorithm is able to probe the complex structure of the free-energy landscape, finding qualitative agreement with mean-field predictions. We also characterize the algorithmic complexity of CALiPPSO and provide an open-source implementation of it., Comment: Changed title to match the version published in PRE. Find our implementation of the CALiPPSO algorithm here: https://github.com/rdhr/CALiPPSO

Recently, the large-order behavior of correlation functions of the $O(N)$-anharmonic oscillator has been analyzed by us in [L. T. Giorgini et el., Phys. Rev. D 101, 125001 (2020)]. Two-loop corrections about the instanton configurations were obtained for the partition function, and the two-point and four-point functions, and the derivative of the two-point function at zero momentum transfer. Here, we attempt to verify the obtained analytic results against numerical calculations of higher-order coefficients for the $O(1)$, $O(2)$, and $O(3)$ oscillators, and demonstrate the drastic improvement of the agreement of the large-order asymptotic estimates and perturbation theory upon the inclusion of the two-loop corrections to the large-order behavior., Comment: 18 pages, 4 figures

We study the energy minima of the fully-connected $m$-components vector spin glass model at zero temperature in an external magnetic field for $m\ge 3$. The model has a zero temperature transition from a paramagnetic phase at high field to a spin glass phase at low field. We study the eigenvalues and eigenvectors of the Hessian in the minima of the Hamiltonian. The spectrum is gapless both in the paramagnetic and in the spin glass phase, with a pseudo-gap behaving as $\lambda^{m-1}$ in the paramagnetic phase and as $\sqrt{\lambda}$ in the spin glass phase. Despite the long-range nature of the model, the eigenstates close to the edge of the spectrum display quasi-localization properties. We show that the paramagnetic to spin glass transition corresponds to delocalization of the edge eigenvectors. We solve the model by the cavity method in the thermodynamic limit. We also perform numerical minimization of the Hamiltonian for $N\le 2048$ and compute the spectral properties, that show very strong corrections to the asymptotic scaling approaching the critical point.

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 forty years later, we confirm the tenets of Bray and Roberts 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., Comment: 8 pages, 4 figures, accepted in PRE

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 \le 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., Comment: 14 pages, 4 figures

A recent 3-XORSAT challenge required to minimize a very complex and rough energy function, typical of glassy models with a random first order transition and a golf course like energy landscape. We present the ideas beyond the quasi-greedy algorithm and its very efficient implementation on GPUs that are allowing us to rank first in such a competition. We suggest a better protocol to compare algorithmic performances and we also provide analytical predictions about the exponential growth of the times to find the solution in terms of free-energy barriers., Comment: 7 pages, 7 figure, EPL format + SM (2 pages)

The synergy between experiment, theory, and simulations enables a microscopic analysis of spin-glass dynamics in a magnetic field in the vicinity of and below the spin-glass transition temperature $T_\mathrm{g}$. The spin-glass correlation length, $\xi(t,t_\mathrm{w};T)$, is analysed both in experiments and in simulations in terms of the waiting time $t_\mathrm{w}$ after the spin glass has been cooled down to a stabilised measuring temperature $T

In critical systems, the effect of a localized perturbation affects points that are arbitrarily far from the perturbation location. In this paper, we study the effect of localized perturbations on the solution of the random dimer problem in $2D$. By means of an accurate numerical analysis, we show that a local perturbation of the optimal covering induces an excitation whose size is extensive with finite probability. We compute the fractal dimension of the excitations and scaling exponents. In particular, excitations in random dimer problems on non-bipartite lattices have the same statistical properties of domain walls in the $2D$ spin glass. Excitations produced in bipartite lattices, instead, are compatible with a loop-erased self-avoiding random walk process. In both cases, we find evidence of conformal invariance of the excitations that is compatible with $\mathrm{SLE}_\kappa$ with parameter $\kappa$ depending on the bipartiteness of the underlying lattice only., Comment: 8 pages

Jamming criticality defines a universality class that includes systems as diverse as glasses, colloids, foams, amorphous solids, constraint satisfaction problems, neural networks, etc. A particularly interesting feature of this class is that small interparticle forces ($f$) and gaps ($h$) are distributed according to nontrivial power laws. A recently developed mean-field (MF) theory predicts the characteristic exponents of these distributions in the limit of very high spatial dimension, $d\rightarrow\infty$ and, remarkably, their values seemingly agree with numerical estimates in physically relevant dimensions, $d=2$ and $3$. These exponents are further connected through a pair of inequalities derived from stability conditions, and both theoretical predictions and previous numerical investigations suggest that these inequalities are saturated. Systems at the jamming point are thus only marginally stable. Despite the key physical role played by these exponents, their systematic evaluation has yet to be attempted. Here, we carefully test their value by analyzing the finite-size scaling of the distributions of $f$ and $h$ for various particle-based models for jamming. Both dimension and the direction of approach to the jamming point are also considered. We show that, in all models, finite-size effects are much more pronounced in the distribution of $h$ than in that of $f$. We thus conclude that gaps are correlated over considerably longer scales than forces. Additionally, remarkable agreement with MF predictions is obtained in all but one model, namely near-crystalline packings. Our results thus help to better delineate the domain of the jamming universality class. We furthermore uncover a secondary linear regime in the distribution tails of both $f$ and $h$. This surprisingly robust feature is understood to follow from the (near) isostaticity of our configurations., Comment: Version accepted for publication

We find a dynamic effect in the non-equilibrium dynamics of a spin glass that closely parallels equilibrium temperature chaos. This effect, that we name dynamic temperature chaos, is spatially heterogeneous to a large degree. The key controlling quantity is the time-growing spin-glass coherence length. Our detailed characterization of dynamic temperature chaos paves the way for the analysis of recent and forthcoming experiments. This work has been made possible thanks to the most massive simulation to date of non-equilibrium dynamics, carried out on the Janus~II custom-built supercomputer., Comment: Version accepted for publication in Communication Physics 10 pages, 9 figures

The correlation length $\xi$, a key quantity in glassy dynamics, can now be precisely measured for spin glasses both in experiments and in simulations. However, known analysis methods lead to discrepancies either for large external fields or close to the glass temperature. We solve this problem by introducing a scaling law that takes into account both the magnetic field and the time-dependent spin-glass correlation length. The scaling law is successfully tested against experimental measurements in a CuMn single crystal and against large-scale simulations on the Janus II dedicated computer., Comment: Revised version, including supplemental material

Spin Glasses (SG) are paradigmatic models for physical, computer science, biological and social systems. The problem of studying the dynamics for SG models is NP hard, i.e., no algorithm solves it in polynomial time. Here we implement the optical simulation of a SG, exploiting the N segments of a wavefront shaping device to play the role of the spin variables, combining the interference at downstream of a scattering material to implement the random couplings between the spins (the J ij matrix) and measuring the light intensity on a number P of targets to retrieve the energy of the system. By implementing a plain Metropolis algorithm, we are able to simulate the spin model dynamics, while the degree of complexity of the potential energy landscape and the region of phase diagram explored is user-defined acting on the ratio the P/N = \alpha. We study experimentally, numerically and analytically this peculiar system displaying a paramagnetic, a ferromagnetic and a SG phase, and we demonstrate that the transition temperature T g to the glassy phase from the paramagnetic phase grows with \alpha. With respect to standard in silico approach, in the optical SG interaction terms are realized simultaneously when the independent light rays interferes at the target screen, enabling inherently parallel measurements of the energy, rather than computations scaling with N as in purely in silico simulations.

In this Letter, we analyze the quantum dynamics of the perceptron model: a particle is constrained on a $N$-dimensional sphere, with $N\to \infty$, and subjected to a set of randomly placed hard-wall potentials. This model has several applications, ranging from learning protocols to the effective description of the dynamics of an ensemble of infinite-dimensional hard spheres in Euclidean space. We find that the jamming transition with quantum dynamics shows critical exponents different from the classical case. We also find that the quantum jamming transition, unlike the typical quantum critical points, is not confined to the zero-temperature axis, and the classical results are recovered only at $T=\infty$. Our findings have implications for the theory of glasses at ultra-low temperatures and for the study of quantum machine-learning algorithms.

We discuss a phase transition in spin glass models which 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 huge 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)

In this paper, we investigate the finite connectivity spin-glass problem. Our work is focused on the expansion around the point of infinite connectivity of the free energy of a spin glass on a graph with Poissonian distributed connectivity: we are interested to study the first-order correction to the infinite connectivity result for large values or the connectivity $z$. The same calculations for one and two replica symmetry breakings were done in previous works; the result for the first-order correction was divergent in the limit of zero temperature and it was suggested that it was an artifact for having a finite number of replica symmetry breakings. In this paper we are able to calculate the expansion for an infinite number of replica symmetry breakings: in the zero-temperature limit, we obtain a well defined free energy. We have shown that cancellations of divergent terms occur in the case of an infinite number of replica symmetry breakings and that the pathological behavior of the expansion was due only to the finite number of replica symmetry breakings., Comment: 16 pages, 4 figures

We consider the storage properties of temporal patterns, i.e. cycles of finite lengths, in neural networks represented by (generally asymmetric) spin glasses defined on random graphs. Inspired by the observation that dynamics on sparse systems have more basins of attractions than the dynamics of densely connected ones, we consider the attractors of a greedy dynamics in sparse topologies, considered as proxy for the stored memories. We enumerate them using numerical simulation and extend the analysis to large systems sizes using belief propagation. We find that the logarithm of the number of such cycles is a non monotonic function of the mean connectivity and we discuss the similarities with biological neural networks describing the memory capacity of the hippocampus., Comment: 10 pages, 11 figures

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., Comment: Same as published version; better figures placement

In the paper [Angelini M C, Parisi G, and Ricci-Tersenghi F, Ensemble renormalization group for disordered systems, Phys. Rev. B 87 134201 (2013)] we introduced a real-space renormalization group called Ensemble Renormalization Group (ERG) and we applied it to the Edwards-Anderson model, obtaining estimates for the critical exponents in good agreement with those from Monte Carlo simulations. Recently the paper [Castellana M, Real-space renormalization-group methods for hierarchical spin glasses, J. Phys. A: Math. Theor. 52 445002 (2019)] re-examined the ERG method from a different perspective, concluding that the previous results were wrong, and claiming that the ERG method predicts trivially wrong critical exponents. In this comment we explain why the conclusions reached by Castellana are wrong, as they are based on a misinterpretation of finite-size effects. We conclude that the ERG method remains a good RG method to obtain critical exponents in strongly disordered models (if properly used)., Comment: Comment on arXiv:1910.04521

We study the random-link matching problem on random regular graphs, alongside with two relaxed versions of the problem, namely the fractional matching and the so-called "loopy" fractional matching. We estimated the asymptotic average optimal cost using the cavity method. Moreover, we also study the finite-size corrections due to rare topological structures appearing in the graph at large sizes. We estimate these contributions using the cavity approach, and we compare our results with the output of numerical simulations. The analysis also clarifies the meaning of the finite-size contributions appearing in the fully-connected version of the problem, that has been already analyzed in the literature.

We present the study of the landscape structure of athermal soft spheres both 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 has a steep increase and that 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., Comment: 24 pages, 12 figures

We present a complementary estimation of the critical exponent $\alpha$ of the specific heat of the 5D random-field Ising model from zero-temperature numerical simulations. Our result $\alpha = 0.12(2)$ is consistent with the estimation coming from the modified hyperscaling relation and provides additional evidence in favor of the recently proposed restoration of dimensional reduction in the random-field Ising model at $D = 5$., Comment: 11 pages, 7 figures, final version with minor corrections

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 GPUs, we have been able to run the latter model for more than $2^{64}$ 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 $T

We apply to the Random Field Ising Model at zero temperature (T= 0) the perturbative loop expansion around the Bethe solution. A comparison with the standard epsilon-expansion is made, highlighting the key differences that make the new expansion much more appropriate to correctly describe strongly disordered systems, especially those controlled by a T = 0 RG fixed point. This new loop expansion produces an effective theory with cubic vertices. We compute the one-loop corrections due to cubic vertices, finding new terms that are absent in the epsilon-expansion. However, these new terms are subdominant with respect to the standard, supersymmetric ones, therefore dimensional reduction is still valid at this order of the loop expansion., Comment: 6 pages + Supporting Information. Submitted to PNAS

Using the replica approach and the cavity method, we study the fluctuations of the optimal cost in the random-link matching problem. By means of replica arguments, we derive the exact expression of its variance. Moreover, we study the large deviation function, deriving its expression in two different ways, namely using both the replica method and the cavity method., Comment: 9 pages, 3 figures

The ferromagnetic XY model on sparse random graphs in a randomly oriented field is analyzed via the belief propagation algorithm. At variance with the fully connected case and with the random field Ising model on the same topology, we find strong evidences of a tiny region with Replica Symmetry Breaking (RSB) in the limit of very low temperatures. This RSB phase is robust against different choices of the external field direction, while it rapidly vanishes when increasing the graph mean degree, the temperature or the directional bias in the external field. The crucial ingredients to have such a RSB phase seem to be the continuous nature of vector spins, mostly preserved by the O(2)-invariant random field, and the strong spatial heterogeneity, due to graph sparsity. We also uncover that the ferromagnetic phase can be marginally stable despite the presence of the random field. Finally, we study the proper correlation functions approaching the critical points to identify the ones that become more critical., Comment: 15 pages, 10 figures. v3: appendix improved

Recent numerical studies on glassy systems provide evidences for a population of non-Goldstone modes (NGMs) in the low-frequency spectrum of the vibrational density of states $D(\omega)$. 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(\omega)$. In particular, the phonon attenuation is reflected in a parental temperature dependency of the exponent $s(T^*)$ in the low-frequency power law $D(\omega) \sim \omega^{s(T^*)}$, with $2 \leq s(T^*) \leq 4 $. Secondly, by comparing $s(T^*)$ with $s(p)$, i. e., the same quantity obtained by pinning \mttp{a} $p$ particle fraction, we suggest that $s(T^*)$ reflects the presence of dynamical heterogeneous regions of size $\xi^3 \propto p$. Finally, we provide an estimate of $\xi$ as a function of $T^*$, finding a mild power law divergence, $\xi \sim (T^* - T_d)^{-\alpha/3}$, with $T_d$ the dynamical crossover temperature and $\alpha$ falling in the range $\alpha \in [0.8,1.0]$.

We provide a non-trivial 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., Comment: 8 pages, 5 figures, final version with minor modifications to be published in Physical Review Letters

We study finite-size corrections to the free energy of the Sherrington-Kirkpatrick spin glass in the low temperature phase. We investigate the role of longitudinal fluctuations in these corrections, neglecting the transverse contribution. In particular, we are interested in the exponent $\alpha$ defined by the relation $f-f_\infty\sim N^{-\alpha}$. We perform both an analytical and numerical estimate of the analytical result for $\alpha$. From both the approaches we get the result: $\alpha=0.8$., Comment: 12 pages

We study the Anderson model on the Bethe lattice by working directly with propagators at real energies $E$. We introduce a novel criterion for the localization-delocalization transition based on the stability of the population of the propagators, and show that it is consistent with the one obtained through the study of the imaginary part of the self-energy. We present an accurate numerical estimate of the transition point, as well as a concise proof of the asymptotic formula for the critical disorder on lattices of large connectivity, as given in [P.W. Anderson 1958]. We discuss how the forward approximation used in analytic treatments of localization problems fits into this scenario and how one can interpolate between it and the correct asymptotic analysis., Comment: Close to published version

We present a novel mechanism for the anomalous behaviour 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 crossover 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 (ii) The vicinity of the classical jamming critical point, as the crossover 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., Comment: 7 pages + 38 pages SI, 5 figures

We study the problem of glassy relaxations in the presence of an external field in the highly controlled context of a spin-glass simulation. We consider a small spin glass in three dimensions (specifically, a lattice of size L=8, small enough to be equilibrated through a Parallel Tempering simulations at low temperatures, deep in the spin glass phase). After equilibrating the sample, an external field is switched on, and the subsequent dynamics is studied. The field turns out to reduce the relaxation time, but huge statistical fluctuations are found when different samples are compared. After taking care of these fluctuations we find that the expected linear regime is very narrow. Nevertheless, when regarded as a purely numerical method, we find that the external field is extremely effective in reducing the relaxation times., Comment: 22 pages, 10 figures; Published version

In this work we want to enhance the calculation performed by Franz, Parisi and Virasoro (FPV) to estimate the free energy cost of interfaces in spin glasses and evaluate the lower critical dimension at which replica symmetry is restored. In particular we evaluate the free energy cost for a general class of effective Hamiltonians showing full replica symmetry breaking, and study the dependence of this cost on the order parameter and on the temperature. We confirm the findings of the FPV papers for the scaling of the free energy, recovering a value for the lower critical dimension of $D_{lc} = 2.5$. In addition to their results we find a non-trivial dependence of the free energy density cost on the order parameter and the temperature. Apart from the case of a restricted class of effective Hamiltonians this dependence cannot be expressed in terms of functions with a clear physical interpretation, as is the case in hierarchical models. In addition we connect the results on the lower critical dimension with recent simulations.

The comprehension of the mechanisms at the basis of the functioning of complexly interconnected networks represents one of the main goals of neuroscience. In this work, we investigate how the structure of recurrent connectivity influences the ability of a network to have storable patterns and in particular limit cycles, by modeling a recurrent neural network with McCulloch-Pitts neurons as a content-addressable memory system. A key role in such models is played by the connectivity matrix, which, for neural networks, corresponds to a schematic representation of the "connectome": the set of chemical synapses and electrical junctions among neurons. The shape of the recurrent connectivity matrix plays a crucial role in the process of storing memories. This relation has already been exposed by the work of Tanaka and Edwards, which presents a theoretical approach to evaluate the mean number of fixed points in a fully connected model at thermodynamic limit. Interestingly, further studies on the same kind of model but with a finite number of nodes have shown how the symmetry parameter influences the types of attractors featured in the system. Our study extends the work of Tanaka and Edwards by providing a theoretical evaluation of the mean number of attractors of any given length $L$ for different degrees of symmetry in the connectivity matrices., Comment: 35 pages, 12 figures

In thermal glasses at temperatures sufficiently lower than the glass transition, the constituent particles are trapped in their cages for sufficiently long time such that their {\em time-averaged positions} can be determined before diffusion and structural relaxation takes place. The effective forces are those that hold these average positions in place. In numerical simulations the effective forces $\B F_{ij}$ between any pair of particles can be measured as a time average of the {\em bare} forces $\B f_{ij}(\B r_{ij}(t))$. In general even if the bare forces come from two-body interactions, thermal dynamics dress the effective forces to contain many-body interactions. Here we develop the effective theory for systems with generic interactions, where the effective forces are derivable from an effective potential and in turn they give rise to an effective Hessian whose eigenvalues are all positive when the system is stable. In this Letter we offer analytic expressions for the effective theory, and demonstrate the usefulness and the predictive power of the approach., Comment: 2 figures

We consider the spatial correlation function of the two-dimensional Ising spin glass under out-equilibrium conditions. We pay special attention to the scaling limit reached upon approaching zero temperature. The field-theory of a non-interacting field makes a surprisingly good job at describing the spatial shape of the correlation function of the out-equilibrium Edwards-Anderson Ising model in two dimensions., Comment: 20 pages + 5 Figures

Recent high precision experimental results on spin-glass films ask for a detailed understanding of the domain-growth dynamics of two-dimensional spin glasses. To achieve this goal, we numerically simulate the out-equilibrium dynamics of the Ising spin glass for a time that spans close to twelve orders of magnitude (from picoseconds to order of a second), in systems large enough to avoid finite-size effects. We find that the time-growth of the size of the glassy domains is excellently described by a single scaling function. A single time-scale $\tau(T)$ controls the dynamics. $\tau(T)$ diverges upon approaching the $T=0$ critical point. The divergence of $\tau(T\to 0)$ is Arrhenius-like, with a barrier height that depends very mildly on temperature. The growth of this barrier-height is best described by critical dynamics. As a side product we obtain an impressive confirmation of universality of the equilibrium behavior of two-dimensional spin-glasses., Comment: 21 pages, 9 figures. Updated references. Added DOI and Journal ref

Bethe Lattice Spin Glasses (BLSG) are models with finite connectivity which undergo a Replica Symmetry Breaking (RSB) phase transition in field, even at zero temperature. The order parameter describing this phase is a function related to the non-trivial organization of the states. We compute numerically the spin glass order parameter at zero temperature in BLSG with coordination number $z=3$ and couplings $J=\pm1$, near the transition. The method is based on a universal formula which relates the order parameter to the joint probability distribution of the energy gaps of the lowest lying states. Such distribution is derived numerically by introducing a convenient bulk perturbation to the Hamiltonian., Comment: 10 pages, 4 figures