Source author record

Giorgio Parisi

Giorgio Parisi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

110works
22topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

110 published item(s)

preprint2026arXiv

Low energy excitations in a long prism geometry: computing the lower critical dimension of the Ising spin glass

We propose a general method for studying systems that display excitations with arbitrarily low energy in their low-temperature phase. We argue that in a rectangular right prism geometry, with longitudinal size much larger than the transverse size, correlations decay exponentially (at all temperatures) along the longitudinal dimension, but the scaling of the correlation length with the transverse size carries crucial information from which the lower critical dimension can be inferred. The method is applied in the particularly demanding context of Ising spin glasses at zero magnetic field. The lower critical dimension and the multifractal spectrum for the correlation function are computed from large-scale numerical simulations. Several technical novelties (such as the unexpectedly crucial performance of Houdayer's cluster method or the convenience of using open - rather than periodic - boundary conditions) allow us to study three-dimensional prisms with transverse dimensions up to $L=24$ and effectively infinite longitudinal dimensions down to low temperatures. The value that we find for the lower critical dimension turns out to be in agreement with expectations from both the Replica Symmetry Breaking theory and the Droplet model for spin glasses. We argue that our novel setting holds promise in clarifying which of the two competing theories more accurately describes three-dimensional spin glasses.

preprint2022arXiv

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 \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.

preprint2021arXiv

Correlation Functions of the Anharmonic Oscillator: Numerical Verification of Two-Loop Corrections to the Large-Order Behavior

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.

preprint2021arXiv

Delocalization transition in low energy excitation modes of vector spin glasses

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 $λ^{m-1}$ in the paramagnetic phase and as $\sqrtλ$ 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.

preprint2020arXiv

An exploratory study of the glassy landscape near jamming

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.

preprint2020arXiv

Free energy expansion of the spin glass with finite connectivity for $\infty$ RSB

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.

preprint2020arXiv

Random-link matching problems on random regular graphs

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.

preprint2020arXiv

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 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)<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.

preprint2019arXiv

Comment on "Real-space renormalization-group methods for hierarchical spin glasses"

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).

preprint2019arXiv

New loop expansion for the Random Magnetic Field Ising Ferromagnets at zero temperature

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.

preprint2019arXiv

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 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<T_c$ 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.

preprint2018arXiv

New analysis of the free energy cost of interfaces in spin glasses

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.

preprint2016arXiv

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 length scales 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.

preprint2016arXiv

Numerical estimate of the Kardar Parisi Zhang universality class in (2 + 1) dimensions

We study the Restricted Solid on Solid model for surface growth in spatial dimension $d=2$ by means of a multi-surface coding technique that allows to produce a large number of samples of samples in the stationary regime in a reasonable computational time. Thanks to: (i) a careful finite-size scaling analysis of the critical exponents, (ii) the accurate estimate of the first three moments of the height fluctuations, we can quantify the wandering exponent with unprecedented precision: $χ_{d=2} = 0.3869(4)$. This figure is incompatible with the long-standing conjecture due to Kim and Koesterlitz that hypothesized $χ_{d=2}=2/5$.

preprint2016arXiv

The jamming transition in high dimension: an analytical study of the TAP equations and the effective thermodynamic potential

We present a parallel derivation of the Thouless-Anderson-Palmer (TAP) equations and of an effective potential for the negative perceptron and soft sphere models in high dimension. Both models are continuous constrained satisfaction problems with a critical jamming transition characterized by the same exponents. Our analysis reveals that a power expansion of the potential up to the second order represents a successful framework to approach the jamming line from the SAT phase (the region of the phase diagram where at least one configuration verifies all the constraints), where the ground-state energy is zero. An interesting outcome is that close to jamming the effective thermodynamic potential has a logarithmic contribution, which turns out to be dominant in a proper scaling regime. Our approach is quite general and can be directly applied to other interesting models. Finally, we study the spectrum of small harmonic fluctuations in the SAT phase recovering the typical scaling $D(ω) \sim ω^2$ below the cutoff frequency but a different behavior characterized by a non-trivial exponent above it.

preprint2015arXiv

Composite operators in Cubic Field Theories and Link Overlap Fluctuations in Spin-Glass Models

We present a complete characterization of the fluctuations and correlations of the squared overlap in the Edwards-Anderson Spin-Glass model in zero field. The analysis reveals that the energy-energy correlations (and thus the specific heat) have a different critical behavior than the fluctuations of the link overlap in spite of the fact that the average energy and average link overlap have the same critical properties. More precisely the link-overlap fluctuations are larger than the specific heat according to a computation at first order in the $6-ε$ expansion. An unexpected outcome is that the link-overlap fluctuations have a subdominant power-law contribution characterized by an anomalous logarithmic prefactor which is missing in the specific heat. In order to compute the $ε$ expansion we consider the problem of the renormalization of quadratic composite operators in a generic multicomponent cubic field theory: the results obtained have a range of applicability beyond spin-glass theory.

preprint2015arXiv

Cross correlations of the American baby names

The quantitative description of cultural evolution is a challenging task. The most difficult part of the problem is probably to find the appropriate measurable quantities that can make more quantitative such evasive concepts as, for example, dynamics of cultural movements, behavior patterns and traditions of the people. A strategy to tackle this issue is to observe particular features of human activities, i.e. cultural traits, such as names given to newborns. We study the names of babies born in the United States of America from 1910 to 2012. Our analysis shows that groups of different correlated states naturally emerge in different epochs, and we are able to follow and decrypt their evolution. While these groups of states are stable across many decades, a sudden reorganization occurs in the last part of the twentieth century. We think that this kind of quantitative analysis can be possibly extended to other cultural traits: although databases covering more than one century (as the one we used) are rare, the cultural evolution on shorter time scales can be studied thanks to the fact that many human activities are usually recorded in the present digital era.

preprint2015arXiv

Diluted Mean-Field Spin-Glass Models at Criticality

We present a method derived by cavity arguments to compute the spin-glass and higher-order susceptibilities in diluted mean-field spin-glass models. The divergence of the spin-glass susceptibility is associated to the existence of a non-zero solution of a homogeneous linear integral equation. Higher order susceptibilities, relevant for critical dynamics through the parameter exponent $λ$, can be expressed at criticality as integrals involving the critical eigenvector. The numerical evaluation of the corresponding analytic expressions is discussed. The method is illustrated in the context of the de Almeida-Thouless line for a spin-glass on a Bethe lattice but can be generalized straightforwardly to more complex situations.

preprint2015arXiv

Language discrimination and clustering via a neural network approach

We classify twenty-one Indo-European languages starting from written text. We use neural networks in order to define a distance among different languages, construct a dendrogram and analyze the ultrametric structure that emerges. Four or five subgroups of languages are identified, according to the "cut" of the dendrogram, drawn with an entropic criterion. The results and the method are discussed.

preprint2015arXiv

Non-perturbative effects in spin glasses

We present a numerical study of an Ising spin glass with hierarchical interactions - the hierarchical Edwards-Anderson model with an external magnetic field (HEA). We study the model with Monte Carlo (MC) simulations in the mean-field (MF) and non-mean-field (NMF) regions corresponding to $d\geq4$ and $d<4$ for the $d$-dimensional ferromagnetic Ising model respectively. We compare the MC results with those of a renormalization-group (RG) study where the critical fixed point is treated as a perturbation of the MF one, along the same lines as in the $ε$-expansion for the Ising model. The MC and the RG method agree in the MF region, predicting the existence of a transition and compatible values of the critical exponents. Conversely, the two approaches markedly disagree in the NMF case, where the MC data indicates a transition, while the RG analysis predicts that no perturbative critical fixed point exists. Also, the MC estimate of the critical exponent $ν$ in the NMF region is about twice as large as its classical value, even if the analog of the system dimension is within only $\sim 2\%$ from its upper-critical-dimension value. Taken together, these results indicate that the transition in the NMF region is governed by strong non-perturbative effects.

preprint2015arXiv

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

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

preprint2015arXiv

Off-equilibrium finite-size method for critical behavior analyses

We present a new dynamic off-equilibrium method for the study of continuous transitions, which represents a dynamic generalization of the usual equilibrium cumulant method. Its main advantage is that critical parameters are derived from numerical data obtained much before equilibrium has been attained. Therefore, the method is particularly useful for systems with long equilibration times, like spin glasses. We apply it to the three-dimensional Ising spin-glass model, obtaining accurate estimates of the critical exponents and of the critical temperature with a limited computational effort.

preprint2015arXiv

Quasi equilibrium construction for the long time limit of glassy dynamics

In this paper we review a recent proposal to understand the long time limit of glassy dynamics in terms of an appropriate Markov Chain. [1]. The advantages of the resulting construction are many. The first one is that it gives a quasi equilibrium description on how glassy systems explore the phase space in the slow relaxation part of their dynamics. The second one is that it gives an alternative way to obtain dynamical equations starting from a dynamical rule that is static in spirit. This provides a way to overcome the difficulties encountered in the short time part of the dynamics where current conservation must be enforced. We study this approach in detail in a prototypical mean field disordered spin system, namely the p-spin spherical model, showing how we can obtain the well known equations that describes its dynamics. Then we apply the same approach to structural glasses. We first derive a set of dynamical Ornstein-Zernike equations which are very general in nature. Finally we consider two possible closure schemes for them, namely the Hypernetted Chain approximation of liquid theory and a closure of the BBGKY hierarchy that has been recently introduced by G. Szamel. From both approaches we finally find a set of dynamical Mode-Coupling like equations that are supposed to describe the system in the long time/slow dynamics regime.

preprint2015arXiv

Quasi-equilibrium in glassy dynamics: a liquid theory approach

We introduce a quasi-equilibrium formalism in the theory of liquids in order to obtain a set of coarse grained long time dynamical equations for the two point density correlation functions. Our scheme allows to use typical approximations devised for equilibrium to study long time glassy dynamics. We study the Hypernetted Chain (HNC) approximation and a recent closure scheme by Szamel. In both cases we get dynamical equations that have the structure of the Mode-Coupling (MCT) equations in the long time region. The HNC approach, that was so far used to get equilibrium quantities is thus generalized to a fully consistent scheme where long-time dynamic quantities can also be computed. In the context of this approximation we get an asymptotic description of both equilibrium glassy dynamics at high temperature and of aging dynamics at low temperature. The Szamel approximation on the other hand is shown to lead to the exact Mode Coupling equation of Götze for equilibrium dynamics. We clarify the way phase space is sampled according to MCT during dynamical relaxation.

preprint2015arXiv

Soft Modes, Localization and Two-Level Systems in Spin Glasses

In the three-dimensional Heisenberg spin glass in a random field we study the properties of the inherent structures that are obtained by an instantaneous cooling from infinite temperature. For not too large field the density of states $g(ω)$ develops localized soft {\em plastic} modes and reaches zero as $ω^4$ (for large fiels a gap appears). When we perturb the system adding a force along the softest mode one reaches very similar minima of the energy, separated by small barriers, that appear to be good candidates for classical two-level systems.

preprint2015arXiv

The large connectivity limit of bootstrap percolation

Bootstrap percolation provides an emblematic instance of phase behavior characterised by an abrupt transition with diverging critical fluctuations. This unusual hybrid situation generally occurs in particle systems in which the occupation probability of a site depends on the state of its neighbours through a certain threshold parameter. In this paper we investigate the phase behavior of the bootstrap percolation on the regular random graph in the limit in which the threshold parameter and lattice connectivity become both increasingly large while their ratio $α$ is held constant. We find that the mixed phase behavior is preserved in this limit, and that multiple transitions and higher-order bifurcation singularities occur when $α$ becomes a random variable.

preprint2015arXiv

The simplest model of jamming

We study a well known machine learning model -the perceptron- as a simple model of jamming of hard objects. We exhibit two regimes: 1) a convex optimisation regime where jamming is hypostatic and non-critical. 2) a non convex optimisation regime where jamming is isostatic and critical. We characterise the critical jamming phase through exponents describing the distributions law of forces and gaps. Surprisingly we find that these exponents coincide with the corresponding ones recently computed in high dimensional hard spheres. In addition, modifying the perceptron to a random linear programming problem, we show that isostaticity is not a sufficient condition for singular force and gap distributions. For that, fragmentation of the space of solutions (replica symmetry breaking) appears to be a crucial ingredient. We hypothesise universality for a large class of non-convex constrained satisfaction problems with continuous variables.

preprint2015arXiv

Universal non-Debye scaling in the density of states of amorphous solids

At the jamming transition, amorphous packings are known to display anomalous vibrational modes with a density of states (DOS) that remains constant at low frequency. The scaling of the DOS at higher densities remains, however, unclear. One might expect to find simple Debye scaling, but recent results from effective medium theory and the exact solution of mean-field models both predict an anomalous, non-Debye scaling. Being mean-field solutions, however, these solutions are only strictly applicable to the limit of infinite spatial dimension, and it is unclear what value they have for finite-dimensional systems. Here, we study packings of soft spheres in dimensions 3 through 7 and find, far from jamming, a universal non-Debye scaling of the DOS that is consistent with the mean-field predictions. We also consider how the soft mode participation ratio converges to the mean-field prediction as dimension increases.

preprint2015arXiv

Universal Spectrum of Normal Modes in Low-Temperature Glasses: an Exact Solution

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 towards a detailed analytical understanding of the vibrational spectrum of low-temperature glasses.

preprint2014arXiv

Anomalous finite size corrections in random field models

The presence of a random magnetic field in ferromagnetic systems leads, in the broken phase, to an anomalous $O(\sqrt{1/N})$ convergence of some thermodynamic quantities to their asymptotic limits. Here we show a general method, based on the replica trick, to compute analytically the $O(\sqrt{1/N})$ finite size correction to the average free energy. We apply this method to two mean field Ising models, fully connected and random regular graphs, and compare the results to exact numerical algorithms. We argue that this behaviour is present in finite dimensional models as well.

preprint2014arXiv

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 to thermalize the system in the supercooled liquid phase, where the thermalization time may be extremely large. We consider here an hard-sphere mean field model which is solvable in any number of spatial dimensions; moreover we easily obtain thermalized configurations even in the glass phase. We study the three dimensional version of this model and we perform Monte Carlo simulations which mimic heating and cooling experiments performed on ultra-stable glasses. The numerical findings are in good agreement with the analytical results and qualitatively capture the features of ultra-stable glasses observed in experiments.

preprint2014arXiv

Ensemble renormalization group for the random field hierarchical model

The Renormalization Group (RG) methods are still far from being completely understood in quenched disordered systems. In order to gain insight into the nature of the phase transition of these systems, it is common to investigate simple models. In this work we study a real-space RG transformation on the Dyson hierarchical lattice with a random field, which led to a reconstruction of the RG flow and to an evaluation of the critical exponents of the model at T = 0. We show that this method gives very accurate estimations of the critical exponents, by comparing our results with the ones obtained by some of us using an independent method.

preprint2014arXiv

Exact theory of dense amorphous hard spheres in high dimension. III. The full RSB solution

In the first part of this paper, we derive the general replica equations that describe infinite-dimensional hard spheres at any level of replica symmetry breaking (RSB) and in particular in the fullRSB scheme. We show that these equations are formally very similar to the ones that have been derived for spin glass models, thus showing that the analogy between spin glasses and structural glasses conjectured by Kirkpatrick, Thirumalai, and Wolynes is realized in a strong sense in the mean field limit. We also suggest how the computation could be generalized in an approximate way to finite dimensional hard spheres. In the second part of the paper, we discuss the solution of these equations and we derive from it a number of physical predictions. We show that, below the Gardner transition where the 1RSB solution becomes unstable, a fullRSB phase exists and we locate the boundary of the fullRSB phase. Most importantly, we show that the fullRSB solution predicts correctly that jammed packings are isostatic, and allows one to compute analytically the critical exponents associated with the jamming transition, which are missed by the 1RSB solution. We show that these predictions compare very well with numerical results.

preprint2014arXiv

Finite size correction to the spectrum of regular random graphs: an analytical solution

We develop a thorough analytical study of the $O(1/N)$ correction to the spectrum of regular random graphs with $N \rightarrow \infty$ nodes. The finite size fluctuations of the resolvent are given in terms of a weighted series over the contributions coming from loops of all possible lengths, from which we obtain the isolated eigenvalue as well as an analytical expression for the $O(1/N)$ correction to the continuous part of the spectrum. The comparison between this analytical formula and direct diagonalization results exhibits an excellent agreement, confirming the correctness of our expression.

preprint2014arXiv

Fractal free energy landscapes in structural glasses

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

preprint2014arXiv

Highly optimized simulations on single- and multi-GPU systems of 3D Ising spin glass

We present a highly optimized implementation of a Monte Carlo (MC) simulator for the three-dimensional Ising spin-glass model with bimodal disorder, i.e., the 3D Edwards-Anderson model running on CUDA enabled GPUs. Multi-GPU systems exchange data by means of the Message Passing Interface (MPI). The chosen MC dynamics is the classic Metropolis one, which is purely dissipative, since the aim was the study of the critical off-equilibrium relaxation of the system. We focused on the following issues: i) the implementation of efficient access patterns for nearest neighbours in a cubic stencil and for lagged-Fibonacci-like pseudo-Random Numbers Generators (PRNGs); ii) a novel implementation of the asynchronous multispin-coding Metropolis MC step allowing to store one spin per bit and iii) a multi-GPU version based on a combination of MPI and CUDA streams. We highlight how cubic stencils and PRNGs are two subjects of very general interest because of their widespread use in many simulation codes. Our code best performances ~3 and ~5 psFlip on a GTX Titan with our implementations of the MINSTD and MT19937 respectively.

preprint2014arXiv

Hopping and the Stokes-Einstein relation breakdown in simple glass formers

One of the most actively debated issues in the study of the glass transition is whether a mean-field description is a reasonable starting point for understanding experimental glass formers. Although the mean-field theory of the glass transition -- like that of other statistical systems -- is exact when the spatial dimension $d\rightarrow\infty$, the evolution of systems properties with $d$ may not be smooth. Finite-dimensional effects could dramatically change what happens in physical dimensions, $d=2,3$. For standard phase transitions finite-dimensional effects are typically captured by renormalization group methods, but for glasses the corrections are much more subtle and only partially understood. Here, we investigate hopping between localized cages formed by neighboring particles in a model that allows to cleanly isolate that effect. By bringing together results from replica theory, cavity reconstruction, void percolation, and molecular dynamics, we obtain insights into how hopping induces a breakdown of the Stokes--Einstein relation and modifies the mean-field scenario in experimental systems. Although hopping is found to supersede the dynamical glass transition, it nonetheless leaves a sizable part of the critical regime untouched. By providing a constructive framework for identifying and quantifying the role of hopping, we thus take an important step towards describing dynamic facilitation in the framework of the mean-field theory of glasses.

preprint2014arXiv

Large Deviations of Correlation Functions in Random Magnets

We present a large deviations theory of the spin-spin correlation functions in the Random Field Ising Model on the Bethe lattice, both at finite and zero temperature. Rare events of atypically correlated variables are particularly important at the critical point: the phase transition is driven by few pairs of strongly correlated spins, while the majority remains basically uncorrelated. At the zero temperature critical point the number of spin pairs correlated over a distance L is shown to be no longer exponential, but only linear in the spins separation.

preprint2014arXiv

Relations between Short Range and Long Range Ising models

We perform a numerical study of the long range (LR) ferromagnetic Ising model with power law decaying interactions ($J \propto r^{-d-σ}$) both on a one-dimensional chain ($d=1$) and on a square lattice ($d=2$). We use advanced cluster algorithms to avoid the critical slowing down. We first check the validity of the relation connecting the critical behavior of the LR model with parameters $(d,σ)$ to that of a short range (SR) model in an equivalent dimension $D$. We then study the critical behavior of the $d=2$ LR model close to the lower critical $σ$, uncovering that the spatial correlation function decays with two different power laws: the effect of the subdominant power law is much stronger than finite size effects and actually makes the estimate of critical exponents very subtle. By including this subdominant power law, the numerical data are consistent with the standard renormalization group (RG) prediction by Sak, thus making not necessary (and unlikely, according to Occam's razor) the recent proposal by Picco of having a new set of RG fixed points, in addition to the mean-field one and the SR one.

preprint2014arXiv

Replica Theory and Spin Glasses

These are notes from the lectures of Giorgio Parisi given at the autumn school "Statistical Physics, Optimization, Inference, and Message-Passing Algorithm", that took place in Les Houches, France from Monday September 30th, 2013, till Friday October 11th, 2013. The school was organized by Florent Krzakala from UPMC and ENS Paris, Federico Ricci-Tersenghi from "La Sapienza" Roma, Lenka Zdeborová from CEA Saclay and CNRS, and Riccardo Zecchina from Politecnico Torino. The first lecture contains an introduction to the replica method, along with a concrete application to the computation of the eigenvalue distribution of random matrices in the GOE. In the second lecture, the solution of the SK model is derived, along with the phenomenon of replica symmetry breaking (RSB). In the third part, the physical meaning of the RSB is explained. The ultrametricity of the space of pure states emerges as a consequence of the hierarchical RSB scheme. Moreover, it is shown how some low temperature properties of physical observables can be derived by invoking the stochastic stability principle. Lecture four contains some rigorous results on the SK model: the existence of the thermodynamic limit, and the proof of the exactness of the hierarchical RSB solution.

preprint2014arXiv

Scaling hypothesis for the Euclidean bipartite matching problem

We propose a simple yet very predictive form, based on a Poisson's equation, for the functional dependence of the cost from the density of points in the Euclidean bipartite matching problem. This leads, for quadratic costs, to the analytic prediction of the large $N$ limit of the average cost in dimension $d=1,2$ and of the subleading correction in higher dimension. A non-trivial scaling exponent, $γ_d=\frac{d-2}{d}$, which differs from the monopartite's one, is found for the subleading correction. We argue that the same scaling holds true for a generic cost exponent in dimension $d>2$.

preprint2014arXiv

Soft modes in jammed hard spheres (I): Mean field theory of the isostatic transition

In this paper we consider different models for soft modes in jammed hard spheres. We show how one can construct mean field models that can be solved analytically. The analytic solution of these models displays an excess of low energy soft modes that become more and more localized by decreasing the energy. A simple solution of these models is found in the infinite dimensional limit.

preprint2014arXiv

The crossover region between long-range and short-range interactions for the critical exponents

It is well know that systems with an interaction decaying as a power of the distance may have critical exponents that are different from those of short-range systems. The boundary between long-range and short-range is known, however the behavior in the crossover region is not well understood. In this paper we propose a general form for the crossover function and we compute it in a particular limit. We compare our predictions with the results of numerical simulations for two-dimensional long-range percolation.

preprint2013arXiv

Critical Exponents of the Random Field Hierarchical Model

We have studied the one dimensional Dyson hierarchical model in presence of a random field. This is a long range model where the interactions scale with the distance with a power law-like form J(r) ~ r^{-ρ} and we can explore mean field and non-mean field behavior by changing ρ. Thus, it can be used to approach the phase transitions in finite-dimensional disordered models. We studied the model at T=0 and we numerically computed its critical exponents in the non-mean field region for Gaussian disorder. We then computed an analytic expression for the critical exponent δ, that holds in the non-mean field region, and we noted an interesting relation between the critical exponents of the disordered model and the ones of the pure model, that seems to break down in the non-mean field region. We finally compare our results for the critical exponents with the expected ones in D-dimensional short range models and with the ones of the straightforward one dimensional long range model.

preprint2013arXiv

Critical Off-Equilibrium Dynamics in Glassy Systems

We consider off-equilibrium dynamics at the critical temperature in a class of glassy system. The off-equilibrium correlation and response functions obey a precise scaling form in the aging regime. The structure of the {\it equilibrium} replicated Gibbs free energy fixes the corresponding {\it off-equilibrium} scaling functions implicitly through two functional equations. The details of the model enter these equations only through the ratio $w_2/w_1$ of the cubic coefficients (proper vertexes) of the replicated Gibbs free energy. Therefore the off-equilibrium dynamical exponents are controlled by the very same parameter exponent $λ=w_2/w_1$ that determines equilibrium dynamics. We find approximate solutions to the equations and validate the theory by means of analytical computations and numerical simulations.

preprint2013arXiv

Dimensional dependence of the Stokes--Einstein relation and its violation

We generalize to higher spatial dimensions the Stokes--Einstein relation (SER) and the leading correction to diffusivity in periodic systems, and validate them using numerical simulations. Using these results, we investigate the evolution of the SER violation with dimension in simple hard sphere glass formers. The analysis suggests that the SER violation disappears around dimension d=8, above which SER is not violated. The critical exponent associated to the violation appears to evolve linearly in 8-d below d=8, as predicted by Biroli and Bouchaud [J. Phys.: Cond. Mat. 19, 205101 (2007)], but the linear coefficient is not consistent with their prediction. The SER violation evolution with d establishes a new benchmark for theory, and a complete description remains an open problem.

preprint2013arXiv

Ensemble renormalization group for disordered systems

We propose and study a renormalization group transformation that can be used also for models with strong quenched disorder, like spin glasses. The method is based on a mapping between disorder distributions, chosen such as to keep some physical properties (e.g., the ratio of correlations averaged over the ensemble) invariant under the transformation. We validate this ensemble renormalization group by applying it to the hierarchical model (both the diluted ferromagnetic version and the spin glass version), finding results in agreement with Monte Carlo simulations.

preprint2013arXiv

Entropy-Driven Phase Transitions of Entanglement

We study the behavior of bipartite entanglement at fixed von Neumann entropy. We look at the distribution of the entanglement spectrum, that is the eigenvalues of the reduced density matrix of a quantum system in a pure state. We report the presence of two continuous phase transitions, characterized by different entanglement spectra, which are deformations of classical eigenvalue distributions.

preprint2013arXiv

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

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

preprint2013arXiv

Finite size corrections to disordered systems on Erdös-Rényi random graphs

We study the finite size corrections to the free energy density in disorder spin systems on sparse random graphs, using both replica theory and cavity method. We derive an analytical expressions for the $O(1/N)$ corrections in the replica symmetric phase as a linear combination of the free energies of open and closed chains. We perform a numerical check of the formulae on the Random Field Ising Model at zero temperature, by computing finite size corrections to the ground state energy density.

preprint2013arXiv

Imry-Ma criterion for long-range random field Ising model: short-/long-range equivalence in a field

The Ising model in a random field and with power-law decaying ferromagnetic bonds is studied at zero temperature. Comparing the scaling of the energy contributions of the ferromagnetic domain wall flip and of the random field a la Imry-Ma we obtain a threshold value for the power $ρ$ of the long-range interaction, beyond which no critical behavior occurs. The critical threshold value is $ρ_c=3/2$, at a difference with the zero field model in which $ρ_c=2$. This prediction is confirmed by numerical computation of the ground states below, at, and above this threshold value. Some possible implications for the critical behavior of spin-glasses in a field are conjectured.

preprint2013arXiv

Jamming transition of randomly pinned systems

We consider a system of hard spheres close to jamming, where translation invariance is broken by pinning a randomly chosen set of particles. Using two different protocols, we generate two kinds of packings at the jamming point, isostatic and hyperstatic packings. In the case where the packings are isostatic, jamming transition is not affected by random pinning: the jamming density is only slightly reduced, a generalized isostaticity condition holds, the system is marginally rigid over the entire glass phase, and the typical structural signatures of jamming are unchanged. Besides, random pinning does not have effect on the vibrational modes of the amorphous system, at least on the time and length scale that we are able to probe. For packings that are hyperstatic at jamming, the microscopic properties of the system are strongly affected by the random pinning: the frozen degrees of freedom introduce an excess of constraints that drives the system far from the mechanical marginality conditions and stabilize the low frequency modes of the system, shifting them to higher frequencies. The distance from the isostatic point seems to be the only relevant parameter in the system in both cases. These two cases are in contrast with the behavior of the plane waves of a crystal: the whole spectrum in this case is shifted to higher frequencies as soon as some particles are pinned.

preprint2013arXiv

Multi-surface coding simulations of the restricted solid-on-solid model in four dimensions

We study the Restricted Solid on Solid (RSOS) model for surface growth in spatial dimension d=4 by means of a multi-surface coding technique that allows to analyze samples to analyze samples of size up to $256^4$ in the steady state regime. For such large systems we are able to achieve a controlled asymptotic regime where the typical scale of the fluctuations are larger than the lattice spacing used in the simulations. A careful finite-size scaling analysis of the critical exponents clearly indicate that d=4 is not the upper critical dimension of the model.

preprint2013arXiv

Quantitative field theory of the glass transition

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

preprint2013arXiv

Spatial correlation functions and dynamical exponents in very large samples of 4D spin glasses

The study of the low temperature phase of spin glass models by means of Monte Carlo simulations is a challenging task, because of the very slow dynamics and the severe finite size effects they show. By exploiting at the best the capabilities of standard modern CPUs (especially the SSE instructions), we have been able to simulate the four-dimensional (4D) Edwards-Anderson model with Gaussian couplings up to sizes $L=70$ and for times long enough to accurately measure the asymptotic behavior. By quenching systems of different sizes to the the critical temperature and to temperatures in the whole low temperature phase, we have been able to identify the regime where finite size effects are negligible: $ξ(t) \lesssim L/7$. Our estimates for the dynamical exponent ($z \simeq 1/T$) and for the replicon exponent ($α\simeq 1.0$ and $T$-independent), that controls the decay of the spatial correlation in the zero-overlap sector, are consistent with the RSB theory, but the latter differs from the theoretically conjectured value.

preprint2013arXiv

Static replica approach to critical correlations in glassy systems

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

preprint2013arXiv

Temperature chaos and quenched heterogeneities

We present a treatable generalization of the Sherrington Kirkpatrick model which introduces correlations in the elements of the coupling matrix through multiplicative disorder on the single variables and investigate the consequences on the phase diagram. We define a generalized $q_{EA}$ parameter and test the structural stability of the SK results in this correlated case evaluating the AT-line of the model. As main result we demonstrate the increase of temperature chaos effects due to heterogeneities.

preprint2013arXiv

The overlap in glassy systems

In this paper I will consider many of the various definitions of the overlap and of its probability distribution that have been introduced in the literature starting from the original papers of Edwards and Anderson; I will present also some of the most recent results on the probability distribution of the local overlap in spin glasses. These quantities are related to the fluctuation dissipation relations both in their local and in their global versions.

preprint2013arXiv

Universality Classes of Critical Points in Constrained Glasses

We analyze critical points that can be induced in glassy systems by the presence of constraints. These critical points are predicted by the Mean Field Thermodynamic approach and they are precursors of the standard glass transition in absence of constraints. Through a deep analysis of the soft modes appearing in the replica field theory we can establish the universality class of these points. In the case of the "annealed potential" of a symmetric coupling between two copies of the system, the critical point is in the Ising universality class. More interestingly, is the case of the "quenched potential" where the a single copy is coupled with an equilibrium reference configuration, or the "pinned particle" case where a fraction of particles is frozen in fixed positions. In these cases we find the Random Field Ising Model (RFIM) universality class. The effective random field is a "self-generated" disorder that reflects the random choice of the reference configuration. The RFIM representation of the critical theory predicts non-trivial relations governing the leading singular behavior of relevant correlation functions, that can be tested in numerical simulations.

preprint2012arXiv

Critical Dynamics in Glassy Systems

Critical dynamics in various glass models including those described by mode coupling theory is described by scale-invariant dynamical equations with a single non-universal quantity, i.e. the so-called parameter exponent that determines all the dynamical critical exponents. We show that these equations follow from the structure of the static replicated Gibbs free energy near the critical point. In particular the exponent parameter is given by the ratio between two cubic proper vertexes that can be expressed as six-point cumulants measured in a purely static framework.

preprint2012arXiv

Critical slowing down exponents in quenched disordered spin models for structural glasses: Random Orthogonal and related models

An important prediction of Mode-Coupling-Theory (MCT) is the relationship between the power- law decay exponents in the β regime. In the original structural glass context this relationship follows from the MCT equations that are obtained making rather uncontrolled approximations and λ has to be treated like a tunable parameter. It is known that a certain class of mean-field spin-glass models is exactly described by MCT equations. In this context, the physical meaning of the so called parameter exponent λ has recently been unveiled, giving a method to compute it exactly in a static framework. In this paper we exploit this new technique to compute the critical slowing down exponents in a class of mean-field Ising spin-glass models including, as special cases, the Sherrington-Kirkpatrick model, the p-spin model and the Random Orthogonal model.

preprint2012arXiv

Dimensional Study of the Caging Order Parameter at the Glass Transition

The glass problem is notoriously hard and controversial. Even at the mean-field level, little is agreed about how a fluid turns sluggish while exhibiting but unremarkable structural changes. It is clear, however, that the process involves self-caging, which provides an order parameter for the transition. It is also broadly assumed that this cage should have a Gaussian shape in the mean-field limit. Here we show that this ansatz does not hold. By performing simulations as a function of spatial dimension, we find the cage to keep a non-trivial form. Quantitative mean-field descriptions of the glass transition, such as mode-coupling theory, density functional theory, and replica theory, all miss this crucial element. Although the mean-field random first-order transition scenario of the glass transition is here qualitatively supported and non-mean-field corrections are found to remain small on decreasing dimension, reconsideration of its implementation is needed for it to result in a coherent description of experimental observations.

preprint2012arXiv

Dynamical critical exponents for the mean-field Potts glass

In this paper we study the critical behaviour of the fully-connected p-colours Potts model at the dynamical transition. In the framework of Mode Coupling Theory (MCT), the time autocorrelation function displays a two step relaxation, with two exponents governing the approach to the plateau and the exit from it. Exploiting a relation between statics and equilibrium dynamics which has been recently introduced, we are able to compute the critical slowing down exponents at the dynamical transition with arbitrary precision and for any value of the number of colours p. When available, we compare our exact results with numerical simulations. In addition, we present a detailed study of the dynamical transition in the large p limit, showing that the system is not equivalent to a random energy model.

preprint2012arXiv

Exact theory of dense amorphous hard spheres in high dimension. I. The free energy

We consider the theory of the glass transition and jamming of hard spheres in the large space dimension limit. Previous investigations were based on the assumption that the probability distribution within a "cage" is Gaussian, which is not fully consistent with numerical results. Here we perform a replica calculation without making any assumption on the cage shape. We show that thermodynamic functions turn out to be exact within the Gaussian ansatz -- provided one allows for arbitrary replica symmetry breaking --- and indeed agree well with numerical results. The actual structure function (the so-called non-ergodic parameter) is not Gaussian, an apparent paradox which we discuss. In this paper we focus on the free energy, future papers will present the results for the structure functions and a detailed comparison with numerical results.

preprint2012arXiv

Generalized Fluctuation-Dissipation Relation and Effective Temperature upon Heating a Deeply Supercooled Liquid

We show that a generalized fluctuation-dissipation relation applies upon instantaneously increasing the temperature of a deeply supercooled liquid. This has the same two-step shape of the relation found upon cooling the liquid, but with opposite violation, indicating an effective temperature that is lower than bath temperature. We show that the effective temperature exhibits some sensible time-dependence and that it retains its connection with the partitioned phase space visited in ageing. We underline the potential relevance of our numerical results for experimental studies of the fluctuation-dissipation relation in glassy systems.

preprint2012arXiv

Glassy Critical Points and Random Field Ising Model

We consider the critical properties of points of continuous glass transition as one can find in liquids in presence of constraints or in liquids in porous media. Through a one loop analysis we show that the critical Replica Field Theory describing these points can be mapped in the $ϕ^4$-Random Field Ising Model. We confirm our analysis studying the finite size scaling of the $p$-spin model defined on sparse random graph, where a fraction of variables is frozen such that the phase transition is of a continuous kind.

preprint2012arXiv

On the two-steps relaxation of mean-field glasses: p-spin model

Critical slowing down dynamics of supercooled glass-forming liquids is usually understood at the mean-field level in the framework of Mode Coupling Theory, providing a two-time relaxation scenario and power-law behaviors of the time correlation function at dynamic criticality. In this work we derive critical slowing down exponents of spin-glass models undergoing discontinuous transitions by computing their Gibbs free energy and connecting the dynamic behavior to static "in-state" properties. Both the spherical and Ising versions are considered and, in the simpler spherical case, a generalization to arbitrary schematic Mode Coupling kernels is presented. Comparison with dynamic results available in literature is performed. Analytical predictions for the Ising case are provided for any $p$.

preprint2012arXiv

Quasi-equilibrium in glassy dynamics: an algebraic view

We study a chain of identical glassy systems in a constrained equilibrium where each bond of the chain is forced to remain at a preassigned distance to the previous one. We apply this description to Mean Field Glassy systems in the limit of long chain where each bond is close to the previous one. We show that in specific conditions this pseudo-dynamic process can formally describe real relaxational dynamics the long time. In particular, in mean field spin glass models we can recover in this way the equations of Langevin dynamics in the long time limit at the dynamical transition temperature and below. We interpret the formal identity as an evidence that in these situations the configuration space is explored in a quasi-equilibrium fashion. Our general formalism, that relates dynamics to equilibrium puts slow dynamics in a new perspective and opens the way to the computation of new dynamical quantities in glassy systems.

preprint2012arXiv

Random Pinning Glass Model

Glass transition where viscosity of liquids increases dramatically upon decrease of temperature without any major change in structural properties, remains one of the most challenging problems in condensed matter physics (Cavagna, 2009; Berthier and Biroli, 2011) in spite of tremendous research efforts in last decades. On the other hand disordered freezing of spins in a magnetic materials with decreasing temperature, the so-called spin glass transition, is relatively better understood (Mezard, Parisi and Virasoro, 1987; Castellani and Cavagna, 2005). Previously found similarity between some spin glass models with the structural glasses (Kirkpatrick and Thirumalai, 1987; Kirkpatrick and Wolynes, 1987; Kirkpatrick and Wolynes, 1987; Franz and Parisi, 1999; Moore and Drossel, 2002) inspired development of theories of structural glasses (Kirkpatrick, Thirumalai and Wolynes, 1989; Barrat, Franz and Parisi, 1997; Mézard and Parisi, 1999; Lubchenko and Wolynes, 2007; Biroli and Bouchaud, 2012) based on the scenario of spin glass transition. This scenario though looks very appealing is still far from being well established. One of the main differences between standard spin systems to molecular systems is the absence of quenched disorder and the presence of translational invariance: it often assumed that this difference is not relevant, but this conjecture is still far from being established. The quantities, which are well defined and characterized for spin models, are not easily calculable for molecular glasses due to the lack of quenched disorder which breaks the translational invariance in the system and the characterization of the similarity between the spin and the structural glass transition remained an elusive subject still now. In this study we introduced a model structural glass with built in quenched disorder which alleviates this main difference between the spin and molecular glasses thereby helping us to compare these two systems: the possibility of producing a good thermalization at rather low temperatures is one of the advantages of this model.

preprint2012arXiv

Universal microstructure and mechanical stability of jammed packings

Jammed packings' mechanical properties depend sensitively on their detailed local structure. Here we provide a complete characterization of the pair correlation close to contact and of the force distribution of jammed frictionless spheres. In particular we discover a set of new scaling relations that connect the behavior of particles bearing small forces and those bearing no force but that are almost in contact. By performing systematic investigations for spatial dimensions d=3-10, in a wide density range and using different preparation protocols, we show that these scalings are indeed universal. We therefore establish clear milestones for the emergence of a complete microscopic theory of jamming. This description is also crucial for high-precision force experiments in granular systems.

preprint2011arXiv

A numerical study of the overlap probability distribution and its sample-to-sample fluctuations in a mean-field model

In this paper we study the fluctuations of the probability distributions of the overlap in mean field spin glasses in the presence of a magnetic field on the De Almeida-Thouless line. We find that there is a large tail in the left part of the distribution that is dominated by the contributions of rare samples. Different techniques are used to examine the data and to stress on different aspects of the contribution of rare samples.

preprint2011arXiv

Field Theory of Fluctuations in Glasses

We develop a field-theoretical description of dynamical heterogeneities and fluctuations in supercooled liquids close to the (avoided) MCT singularity. Using quasi-equilibrium arguments we eliminate time from the description and we completely characterize fluctuations in the beta regime. We identify different sources of fluctuations and show that the most relevant ones are associated to variations of "self-induced disorder" in the initial condition of the dynamics. It follows that heterogeneites can be describes through a cubic field theory with an effective random field term. The phenomenon of perturbative dimensional reduction ensues, well known in random field problems, which implies an upper critical dimension of the theory equal to 8. We apply our theory to finite size scaling for mean-field systems and we test its prediction against numerical simulations.

preprint2011arXiv

Glass transition and random close packing above three dimensions

Motivated by a recently identified severe discrepancy between a static and a dynamic theory of glasses, we numerically investigate the behavior of dense hard spheres in spatial dimensions 3 to 12. Our results are consistent with the static replica theory, but disagree with the dynamic mode-coupling theory, indicating that key ingredients of high-dimensional physics are missing from the latter. We also obtain numerical estimates of the random close packing density, which provides new insights into the mathematical problem of packing spheres in large dimension.

preprint2011arXiv

On the solution of a `solvable' model of an ideal glass of hard spheres displaying a jamming transition

We discuss the analytical solution through the cavity method of a mean field model that displays at the same time an ideal glass transition and a set of jamming points. We establish the equations describing this system, and we discuss some approximate analytical solutions and a numerical strategy to solve them exactly. We compare these methods and we get insight into the reliability of the theory for the description of finite dimensional hard spheres.

preprint2011arXiv

Replica symmetry breaking in and around six dimensions

Two, replica symmetry breaking specific, quantities of the Ising spin glass --- the breakpoint x1 of the order parameter function and the Almeida-Thouless line --- are calculated in six dimensions (the upper critical dimension of the replicated field theory used), and also below and above it. The results comfirm that replica symmetry breaking does exist below d=6, and also the tendency of its escalation for decreasing dimension continues. As a new feature, x1 has a nonzero and universal value for d<6 at criticality. Near six dimensions we have x1=3(6-d)+O[(6-d)^2]. A method to expand a generic theory with replica equivalence around the replica symmetric one is also demonstrated.

preprint2010arXiv

A renormalization group computation of the critical exponents of hierarchical spin glasses

The infrared behaviour of a non-mean field spin-glass system is analysed, and the critical exponent related to the divergence of the correlation length is computed at two loops within the epsilon-expansion technique with two independent methods. Both methods yield the same result confirming that the infrared behaviour of the theory if well-defined and the underlying ideas of the Renormalization Group hold also in such non-mean field disordered model. By pushing such calculation to high orders in epsilon, a consistent and predictive non-mean field theory for such disordered system could be established.

preprint2010arXiv

Critical behaviour of large scale dynamical heterogeneities in glasses: a complete theory

In this talk I will present a complete theory for the behaviour of large-scale dynamical heterogeneities in glasses. Following the work arXiv:1001.1746 I will show that we can write a (physically motivated) simple stochastic differential equation that is potentially able to explain the behaviour of large scale dynamical heterogeneities in glasses. It turns out that this behaviour is in the same universality class of the dynamics near the endpoint of a metastable phase in a disordered system, as far as reparametrization invariant quantities are concerned. Therefore Large scale dynamical heterogeneities in glasses have many points in contact with the Barkhausen noise. Numerical verifications of this theory have not yet done, but they are quite possible.

preprint2010arXiv

Interface Energy in the Edwards-Anderson model

We numerically investigate the spin glass energy interface problem in three dimensions. We analyze the energy cost of changing the overlap from -1 to +1 at one boundary of two coupled systems (in the other boundary the overlap is kept fixed to +1). We implement a parallel tempering algorithm that simulate finite temperature systems and work with both cubic lattices and parallelepiped with fixed aspect ratio. We find results consistent with a lower critical dimension $D_c=2.5$. The results show a good agreement with the mean field theory predictions.

preprint2010arXiv

Mean field theory of hard sphere glasses and jamming

Hard spheres are ubiquitous in condensed matter: they have been used as models for liquids, crystals, colloidal systems, granular systems, and powders. Packings of hard spheres are of even wider interest, as they are related to important problems in information theory, such as digitalization of signals, error correcting codes, and optimization problems. In three dimensions the densest packing of identical hard spheres has been proven to be the FCC lattice, and it is conjectured that the closest packing is ordered (a regular lattice, e.g, a crystal) in low enough dimension. Still, amorphous packings have attracted a lot of interest, because for polydisperse colloids and granular materials the crystalline state is not obtained in experiments for kinetic reasons. We review here a theory of amorphous packings, and more generally glassy states, of hard spheres that is based on the replica method: this theory gives predictions on the structure and thermodynamics of these states. In dimensions between two and six these predictions can be successfully compared with numerical simulations. We will also discuss the limit of large dimension where an exact solution is possible. Some of the results we present here have been already published, but others are original: in particular we improved the discussion of the large dimension limit and we obtained new results on the correlation function and the contact force distribution in three dimensions. We also try here to clarify the main assumptions that are beyond our theory and in particular the relation between our static computation and the dynamical procedures used to construct amorphous packings.

preprint2010arXiv

Properties of the perturbative expansion around the mode-coupling dynamical transition in glasses

In this letter we show how to perform a systematic perturbative approach for the mode-coupling theory. The results coincide with those obtained via the replica approach. The upper critical dimension turns out to be always 8 and the correlations have a double pole in momentum space in perturbations theory. Non-perturbative effects are found to be very important. We suggest a possible framework to compute these effects.

preprint2010arXiv

The Hierarchical Random Energy Model

We introduce a Random Energy Model on a hierarchical lattice where the interaction strength between variables is a decreasing function of their mutual hierarchical distance, making it a non-mean field model. Through small coupling series expansion and a direct numerical solution of the model, we provide evidence for a spin glass condensation transition similar to the one occuring in the usual mean field Random Energy Model. At variance with mean field, the high temperature branch of the free-energy is non-analytic at the transition point.

preprint2009arXiv

A theory of amorphous packings of binary mixtures of hard spheres

We extend our theory of amorphous packings of hard spheres to binary mixtures and more generally to multicomponent systems. The theory is based on the assumption that amorphous packings produced by typical experimental or numerical protocols can be identified with the infinite pressure limit of long lived metastable glassy states. We test this assumption against numerical and experimental data and show that the theory correctly reproduces the variation with mixture composition of structural observables, such as the total packing fraction and the partial coordination numbers.

preprint2009arXiv

Large Deviations of the Free-Energy in Diluted Mean-Field Spin-Glass

Sample-to-sample free energy fluctuations in spin-glasses display a markedly different behaviour in finite-dimensional and fully-connected models, namely Gaussian vs. non-Gaussian. Spin-glass models defined on various types of random graphs are in an intermediate situation between these two classes of models and we investigate whether the nature of their free-energy fluctuations is Gaussian or not. It has been argued that Gaussian behaviour is present whenever the interactions are locally non-homogeneous, i.e. in most cases with the notable exception of models with fixed connectivity and random couplings $J_{ij}=\pm \tilde{J}$. We confirm these expectation by means of various analytical results. In particular we unveil the connection between the spatial fluctuations of the populations of populations of fields defined at different sites of the lattice and the Gaussian nature of the free-energy fluctuations. On the contrary on locally homogeneous lattices the populations do not fluctuate over the sites and as a consequence the small-deviations of the free energy are non-Gaussian and scales as in the Sherrington-Kirkpatrick model.

preprint2009arXiv

On the structure of correlations in the three dimensional spin glasses

We investigate the low temperature phase of three-dimensional Edwards-Anderson model with Bernoulli random couplings. We show that at a fixed value $Q$ of the overlap the model fulfills the clustering property: the connected correlation functions between two local overlaps decay as a power whose exponent is independent of $Q$ for all $0\le |Q| < q_{EA}$. Our findings are in agreement with the RSB theory and show that the overlap is a good order parameter.

preprint2009arXiv

Phase diagram and large deviations in the free-energy of mean-field spin-glasses

We consider the probability distribution of large deviations in the spin-glass free energy for the Sherrington-Kirkpatrick mean field model, i.e. the exponentially small probability of finding a system with intensive free energy smaller than the most likely one. This result is obtained by computing $Φ(n,T)=T \bar{Z^n}/ n$, i.e. the average value of the partition function to the power $n$ as a function of $n$. We study in full details the phase diagram of $Φ(n,T)$ in the $(n,T)$ plane computing in particular the stability of the replica-symmetric solution. At low temperatures we compute $Φ(n,T)$ in series of $n$ and $τ=T_c-T$ at high orders using the standard hierarchical ansatz and confirm earlier findings on the $O(n^5)$ scaling. We prove that the $O(n^5)$ scaling is valid at all orders and obtain an exact expression for the coefficient in term of the function $q(x)$. Resumming the series we obtain the large deviations probability at all temperatures. At zero temperature the analytical prediction displays a remarkable quantitative agreement with the numerical data. A similar computation for the simpler spherical model is also performed and the connection between large and small deviations is discussed.

preprint2009arXiv

Replica method and finite volume corrections

In this note we introduce a method to calculate the finite volume corrections to the mean field results for the free energy when replica symmetry is broken at one-step. We find that the naive results are modified by the presence of additional corrections: these corrections can be interpreted as arising from fluctuations in the size of the blocks in the replica approach. The computation suggests a new approach for deriving the replica broken results in a rigorous way.

preprint2009arXiv

Scale-free correlations in bird flocks

From bird flocks to fish schools, animal groups often seem to react to environmental perturbations as if of one mind. Most studies in collective animal behaviour have aimed to understand how a globally ordered state may emerge from simple behavioural rules. Less effort has been devoted to understanding the origin of collective response, namely the way the group as a whole reacts to its environment. Yet collective response is the adaptive key to survivor, especially when strong predatory pressure is present. Here we argue that collective response in animal groups is achieved through scale-free behavioural correlations. By reconstructing the three-dimensional position and velocity of individual birds in large flocks of starlings, we measured to what extent the velocity fluctuations of different birds are correlated to each other. We found that the range of such spatial correlation does not have a constant value, but it scales with the linear size of the flock. This result indicates that behavioural correlations are scale-free: the change in the behavioural state of one animal affects and is affected by that of all other animals in the group, no matter how large the group is. Scale-free correlations extend maximally the effective perception range of the individuals, thus compensating for the short-range nature of the direct inter-individual interaction and enhancing global response to perturbations. Our results suggest that flocks behave as critical systems, poised to respond maximally to environmental perturbations.

preprint2009arXiv

Universality and Deviations in Disordered Systems

We compute the probability of positive large deviations of the free energy per spin in mean-field Spin-Glass models. The probability vanishes in the thermodynamic limit as $P(Δf) \propto \exp[-N^2 L_2(Δf)]$. For the Sherrington-Kirkpatrick model we find $L_2(Δf)=O(Δf)^{12/5}$ in good agreement with numerical data and with the assumption that typical small deviations of the free energy scale as $N^{1/6}$. For the spherical model we find $L_2(Δf)=O(Δf)^{3}$ in agreement with recent findings on the fluctuations of the largest eigenvalue of random Gaussian matrices. The computation is based on a loop expansion in replica space and the non-gaussian behaviour follows in both cases from the fact that the expansion is divergent at all orders. The factors of the leading order terms are obtained resumming appropriately the loop expansion and display universality, pointing to the existence of a single universal distribution describing the small deviations of any model in the full-Replica-Symmetry-Breaking class.

preprint2008arXiv

The STARFLAG handbook on collective animal behaviour: Part I, empirical methods

The most startling examples of collective animal behaviour are provided by very large and cohesive groups moving in three dimensions. Paradigmatic examples are bird flocks, fish schools and insect swarms. However, because of the sheer technical difficulty of obtaining 3D data, empirical studies conducted to date have only considered loose groups of a few tens of animals. Moreover, these studies were very seldom conducted in the field. Recently the STARFLAG project achieved the 3D reconstruction of thousands of birds under field conditions, thus opening the way to a new generation of quantitative studies of collective animal behaviour. Here, we review the main technical problems in 3D data collection of large animal groups and we outline some of the methodological solutions adopted by the STARFLAG project. In particular, we explain how to solve the stereoscopic correspondence - or matching - problem, which was the major bottleneck of all 3D studies in the past.

preprint2008arXiv

The STARFLAG handbook on collective animal behaviour: Part II, three-dimensional analysis

The study of collective animal behaviour must progress through a comparison between the theoretical predictions of numerical models and data coming from empirical observations. To this aim it is important to develop methods of three-dimensional (3D) analysis that are at the same time informative about the structure of the group and suitable to empirical data. In fact, empirical data are considerably noisier than numerical data, and they are subject to several constraints. We review here the tools of analysis used by the STARFLAG project to characterise the 3D structure of large flocks of starlings in the field. We show how to avoid the most common pitfalls i the quantitative analysis of 3D animal groups, with particular attention to the problem of the bias introduced by the border of the group. By means of practical examples, we demonstrate that neglecting border effects gives rise to artefacts when studying the 3D structure of a group. Moreover, we show that mathematical rigour is essential to distinguish important biological properties from trivial geometric features of animal groups.

preprint2007arXiv

Answer to Comment on "Ultrametricity in the Edwards-Anderson Model" arXiv:0709.0894

This reply shows that the argument presented in the comment by Jorg and Krzakala (cond mat 0709.0894) cannot be used to weaken the results presented in our paper on ultrametricity evidence in the 3d Edwards Anderson model (PRL 99, 057206, 2007; cond-mat/0607376). Our work in fact was mainly based on identifying the scaling law that governs the large volume approach to ultrametricity while NO asymptotic analysis has been done in (cond mat 0709.0894). We show here that the same method we used in our paper, when properly applied to the 2d case, reveals the expected lack of RSB picture at positive temperature, despite the fact that for a fixed finite volume some ultrametric features might still be seen in the joint overlap probability distribution. Those features disappear for increasing volume or when the system is away from the critical curve in the (T,d) plane.

preprint2007arXiv

Large Deviations in the Free-Energy of Mean-Field Spin-Glasses

We compute analytically the probability distribution of large deviations in the spin-glass free energy for the Sherrington-Kirkpatrick mean field model, i.e. we compute the exponentially small probability of finding a system with intensive free energy smaller than the most likely one. This result is obtained by computing the average value of the partition function to the power $n$ as a function of $n$. At zero temperature this absolute prediction displays a remarkable quantitative agreement with the numerical data.

preprint2006arXiv

On $k$-Core Percolation in Four Dimensions

The $k$-core percolation on the Bethe lattice has been proposed as a simple model of the jamming transition because of its hybrid first-order/second-order nature. We investigate numerically $k$-core percolation on the four-dimensional regular lattice. For $k=4$ the presence of a discontinuous transition is clearly established but its nature is strictly first order. In particular, the $k$-core density displays no singular behavior before the jump and its correlation length remains finite. For $k=3$ the transition is continuous.

preprint2004arXiv

Zero-Temperature Limit of the SUSY-breaking Complexity in Diluted Spin-Glass Models

We study the SUSY-breaking complexity of the Bethe Lattice Spin-Glass in the zero temperature limit. We consider both the Gaussian and the bimodal distribution of the coupling constants. For $J_{ij}=\pm 1$ the SUSY breaking theory yields fields distributions that concentrate on integer values at low temperatures, at variance with the unbroken SUSY theory. This concentration takes place both in the quenched as well as in the simpler annealed formulation.

preprint2000arXiv

Non trivial overlap distributions at zero temperature

We explore the consequences of Replica Symmetry Breaking at zero temperature. We introduce a repulsive coupling between a system and its unperturbed ground state. In the Replica Symmetry Breaking scenario a finite coupling induces a non trivial overlap probability distribution among the unperturbed ground state and the one in presence of the coupling. We find a closed formula for this probability for arbitrary ultrametric trees, in terms of the parameters defining the tree. The same probability is computed in numerical simulations of a simple model with many ground states, but no ultrametricity: polymers in random media in 1+1 dimension. This gives us an idea of what violation of our formula can be expected in cases when ultrametricity does not hold.

preprint1999arXiv

A first principle computation of the thermodynamics of glasses

We propose a first principle computation of the equilibrium thermodynamics of simple fragile glasses starting from the two body interatomic potential. A replica formulation translates this problem into that of a gas of interacting molecules, each molecule being built of m atoms, and having a gyration radius (related to the cage size) which vanishes at zero temperature. We use a small cage expansion, valid at low temperatures, which allows to compute the cage size, the specific heat (which follows the Dulong and Petit law), and the configurational entropy.

preprint1999arXiv

Physics of glassy systems

In this talk I present some of the recent theoretical results that have been obtained on glassy systems like spin glasses or structural glasses. The physical principles at the basis of the theory are explained in a simple language (without using replicas) and the results are compared with large scale numerical simulations. Finally we introduce the generalized fluctuation dissipation relation that can be directly tested in experiments with the present day technology.

preprint1998arXiv

Testing replica predictions in experiments

We review the predictions of the replica approach both for the statics and for the off-equilibrium dynamics. We stress the importance of the Cugliandolo-Kurchan off-equilibrium fluctuation-dissipation relation in providing a bridge between the statics and the dynamics. We present numerical evidence for the correctness of these relations. This approach allows an experimental determination of the basic parameters of the replica theory.

preprint1998arXiv

Thermodynamics of glasses: a first principle computation

We propose a first principle computation of the thermodynamics of simple fragile glasses starting from the two body interatomic potential. A replica formulation translates this problem into that of a gas of interacting molecules, each molecule being built of $m$ atoms, and having a gyration radius (related to the cage size) which vanishes at zero temperature. We use a small cage expansion, valid at low temperatures, which allows to compute the cage size, the specific heat (which follows the Dulong and Petit law), and the configurational entropy.

preprint1997arXiv

Effective potential in glassy systems: theory and simulations

We study the phase diagram of glassy systems in presence of an attractive coupling among real replicas. We find competition among a localized and a delocalized phase, that are separated by a coexistence line as in ordinary first order phase transitions. The coexistence line terminates in a critical point. We present numerical simulations for binary glasses in which show that this scenario is realized.

preprint1997arXiv

Numerical Simulations of Spin Glass Systems

We discuss the status of Monte Carlo simulations of (mainly finite dimensional) spin glass systems. After a short historical note and a brief theoretical introduction we start by discussing the (crucial) 3D case: the warm phase, the critical point and the cold phase, the ultrametric structure and the out of equilibrium dynamics. With the same style we discuss the cases of 4D and 2D. In a few appendices we give some details about the definition of states and about the tempering Monte Carlo approach.

preprint1997arXiv

On the mean field approach to glassy systems

In these lectures I will study some properties that are shared by many glassy systems. I will shown how some of these properties can be understood in the framework of the mean field approach based on the replica method and I will discuss which are the difficulties which we have to face when we apply the replica method to realistic short range models. We finally present a partially successfully application of the replica method to soft sphere glasses.

preprint1995arXiv

Replica Theory and Large D Josephson Junction Hypercubic Models

We study the statistical mechanics of a $D$-dimensional array of Josephson junctions in presence of a magnetic field on a lattice of side $2$. In the high temperature region the thermodynamical properties can be computed in the limit $D \to \infty$. A conjectural form of the thermodynamic properties in the low temperature phase is obtained by assuming that they are the same of an appropriate spin glass system, based on quenched disordered couplings. Numerical simulations show that this conjecture is very accurate in one regime of the magnetic field, while it is probably slightly inaccurate in a second regime.

preprint1992arXiv

Simulated Tempering: A New Monte Carlo Scheme

We propose a new global optimization method ({\em Simulated Tempering}) for simulating effectively a system with a rough free energy landscape (i.e. many coexisting states) at finite non-zero temperature. This method is related to simulated annealing, but here the temperature becomes a dynamic variable, and the system is always kept at equilibrium. We analyze the method on the Random Field Ising Model, and we find a dramatic improvement over conventional Metropolis and cluster methods. We analyze and discuss the conditions under which the method has optimal performances.