Source author record

Giulio Biroli

Giulio Biroli 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

72works
21topics
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

72 published item(s)

preprint2026arXiv

High-Dimensional Analysis of Gradient Flow for Extensive-Width Quadratic Neural Networks

We study the high-dimensional training dynamics of a shallow neural network with quadratic activation in a teacher-student setup. We focus on the extensive-width regime, where the teacher and student network widths scale proportionally with the input dimension, and the sample size grows quadratically. This scaling aims to describe overparameterized neural networks in which feature learning still plays a central role. In the high-dimensional limit, we derive a dynamical characterization of the gradient flow, in the spirit of dynamical mean-field theory (DMFT). Under l2-regularization, we analyze these equations at long times and characterize the performance and spectral properties of the resulting estimator. This result provides a quantitative understanding of the effect of overparameterization on learning and generalization, and reveals a double descent phenomenon in the presence of label noise, where generalization improves beyond interpolation. In the small regularization limit, we obtain an exact expression for the perfect recovery threshold as a function of the network widths, providing a precise characterization of how overparameterization influences recovery.

preprint2026arXiv

The critical slowing down in diffusion models

Computational sampling has been central to the sciences since the mid-20th century. While machine-learning-based approaches have recently enabled major advances, their behavior remains poorly understood, with limited theoretical control over when and why they succeed. Here we provide such insight for diffusion models-a class of generative schemes highly effective in practice-by analyzing their application to the $O(n)$ model of statistical field theory in the Gaussian limit $n \to \infty$. In this analytically tractable setting, we show that training a score model with a one-layer network architecture matching the exact solution exhibits a form of critical slowing down in parameter learning. This slowing down also impacts the generation process, indicating that the well-known difficulties of sampling near criticality persist even for learned generative models. To overcome this bottleneck, we demonstrate the power of combining architectural depth with physical locality. We find that using a two-layer architecture drastically reduces the critical slowing down, with the training time scaling logarithmically rather than quadratically with system size. By introducing a local score approximation we show that this acceleration in training time can be achieved without increasing the number of neural network parameters. Taken together, these results demonstrate that diffusion models can overcome the critical slowing down through appropriate architectural design, and establish a controlled framework for understanding and improving learned sampling methods in statistical physics and beyond.

preprint2023arXiv

Landscape Complexity for the Empirical Risk of Generalized Linear Models

We present a method to obtain the average and the typical value of the number of critical points of the empirical risk landscape for generalized linear estimation problems and variants. This represents a substantial extension of previous applications of the Kac-Rice method since it allows to analyze the critical points of high dimensional non-Gaussian random functions. Under a technical hypothesis, we obtain a rigorous explicit variational formula for the annealed complexity, which is the logarithm of the average number of critical points at fixed value of the empirical risk. This result is simplified, and extended, using the non-rigorous Kac-Rice replicated method from theoretical physics. In this way we find an explicit variational formula for the quenched complexity, which is generally different from its annealed counterpart, and allows to obtain the number of critical points for typical instances up to exponential accuracy.

preprint2022arXiv

Equilibrium Fluctuations in Mean-field Disordered Models

Mean-field models of glasses that present a random first order transition exhibit highly non-trivial fluctuations. Building on previous studies that focused on the critical scaling regime, we here obtain a fully quantitative framework for all equilibrium conditions. By means of the replica method we evaluate Gaussian fluctuations of the overlaps around the thermodynamic limit, decomposing them in thermal fluctuations inside each state and heterogeneous fluctuations between different states. We first test and compare our analytical results with numerical simulation results for the p-spin spherical model and the random orthogonal model, and then analyze the random Lorentz gas. In all cases, a strong quantitative agreement is obtained. Our analysis thus provides a robust scheme for identifying the key finite-size (or finite-dimensional) corrections to the mean-field treatment of these paradigmatic glass models.

preprint2022arXiv

Local dynamical heterogeneity in glass formers

We study the local dynamical fluctuations in glass-forming models of particles embedded in $d$-dimensional space, in the mean-field limit of $d\to\infty$. Our analytical calculation reveals that single-particle observables, such as squared particle displacements, display divergent fluctuations around the dynamical (or mode-coupling) transition, due to the emergence of nontrivial correlations between displacements along different directions. This effect notably gives rise to a divergent non-Gaussian parameter, $α_2$. The $d\to\infty$ local dynamics therefore becomes quite rich upon approaching the glass transition. The finite-$d$ remnant of this phenomenon further provides a long sought-after, first-principle explanation for the growth of $α_2$ around the glass transition that is \emph{not based on multi-particle correlations}.

preprint2022arXiv

Optimal learning rate schedules in high-dimensional non-convex optimization problems

Learning rate schedules are ubiquitously used to speed up and improve optimisation. Many different policies have been introduced on an empirical basis, and theoretical analyses have been developed for convex settings. However, in many realistic problems the loss-landscape is high-dimensional and non convex -- a case for which results are scarce. In this paper we present a first analytical study of the role of learning rate scheduling in this setting, focusing on Langevin optimization with a learning rate decaying as $η(t)=t^{-β}$. We begin by considering models where the loss is a Gaussian random function on the $N$-dimensional sphere ($N\rightarrow \infty$), featuring an extensive number of critical points. We find that to speed up optimization without getting stuck in saddles, one must choose a decay rate $β<1$, contrary to convex setups where $β=1$ is generally optimal. We then add to the problem a signal to be recovered. In this setting, the dynamics decompose into two phases: an \emph{exploration} phase where the dynamics navigates through rough parts of the landscape, followed by a \emph{convergence} phase where the signal is detected and the dynamics enter a convex basin. In this case, it is optimal to keep a large learning rate during the exploration phase to escape the non-convex region as quickly as possible, then use the convex criterion $β=1$ to converge rapidly to the solution. Finally, we demonstrate that our conclusions hold in a common regression task involving neural networks.

preprint2022arXiv

Rare events and disorder control the brittle yielding of well-annealed amorphous solids

We use atomistic computer simulations to provide a microscopic description of the brittle failure of amorphous materials, and we assess the role of rare events and quenched disorder. We argue that brittle yielding originates at rare soft regions, similarly to Griffiths effects in disordered systems. We numerically demonstrate how localized plastic events in such soft regions trigger macroscopic failure via the propagation of a shear band. This physical picture, which no longer holds in poorly annealed ductile materials, allows us to discuss the role of finite size effects in brittle yielding and reinforces the similarities between yielding and other disorder-controlled nonequilibrium phase transitions.

preprint2022arXiv

The RFOT Theory of Glasses: Recent Progress and Open Issues

The Random First Order Transition (RFOT) theory started with the pioneering work of Kirkpatrick, Thirumalai and Wolynes. It leverages the methods and advances of the theory of disordered systems. It fares remarkably well at reproducing the salient experimental facts of super-cooled liquids. Yet, direct and indisputable experimental validations are missing. In this short survey, we will review recent investigations that broadly support all static aspects of RFOT, but also those for which the standard dynamical extension of the theory appears to be struggling, in particular in relation with facilitation effects. We discuss possible solutions and open issues.

preprint2021arXiv

Amorphous Order & Non-linear Susceptibilities in Glassy Materials

We review 15 years of theoretical and experimental work on the non-linear response of glassy systems. We argue that an anomalous growth of the peak value of non-linear susceptibilities is a signature of growing "amorphous order" in the system, with spin-glasses as a case in point. Experimental results on supercooled liquids are fully compatible with the RFOT prediction of compact "glassites" of increasing volume as temperature is decreased, or as the system ages. We clarify why such a behaviour is hard to explain within purely kinetic theories of glass formation, despite recent claims to the contrary.

preprint2021arXiv

Critical behavior of the Anderson model on the Bethe lattice via a large-deviation approach

We present a new large-deviation approach to investigate the critical properties of the Anderson model on the Bethe lattice close to the localization transition in the thermodynamic limit. Our method allows us to study accurately the distribution of the local density of states (LDoS) down to very small probability tails as small as $10^{-50}$ which are completely out of reach for standard numerical techniques. We perform a thorough analysis of the functional form and of the tails of the probability distribution of the LDoS which yields for the first time a direct, transparent, and precise estimation of the correlation volume close to the Anderson transition. Such correlation volume is found to diverge exponentially when the localization is approached from the delocalized regime, in a singular way that is in agreement with the analytic predictions of the supersymmetric treatment.

preprint2021arXiv

Effects of intraspecific cooperative interactions in large ecosystems

We analyze the role of the Allee effect, a positive correlation between population density and mean individual fitness, for ecological communities formed by a large number of species. Our study is performed using the generalized Lotka-Volterra model with random interactions between species. We obtain the phase diagram and analyze the nature of the multiple equilibria phase. Remarkable differences emerge with respect to the case of the logistic growth case, thus revealing the major role played by the functional response in determining aggregate behaviours of ecosystems.

preprint2020arXiv

An analytic theory of shallow networks dynamics for hinge loss classification

Neural networks have been shown to perform incredibly well in classification tasks over structured high-dimensional datasets. However, the learning dynamics of such networks is still poorly understood. In this paper we study in detail the training dynamics of a simple type of neural network: a single hidden layer trained to perform a classification task. We show that in a suitable mean-field limit this case maps to a single-node learning problem with a time-dependent dataset determined self-consistently from the average nodes population. We specialize our theory to the prototypical case of a linearly separable dataset and a linear hinge loss, for which the dynamics can be explicitly solved. This allow us to address in a simple setting several phenomena appearing in modern networks such as slowing down of training dynamics, crossover between rich and lazy learning, and overfitting. Finally, we asses the limitations of mean-field theory by studying the case of large but finite number of nodes and of training samples.

preprint2020arXiv

Anomalous dynamics in the ergodic side of the Many-Body Localization transition and the glassy phase of Directed Polymers in Random Media

Using the non-interacting Anderson tight-binding model on the Bethe lattice as a toy model for the many-body quantum dynamics, we propose a novel and transparent theoretical explanation of the anomalously slow dynamics that emerges in the bad metal phase preceding the Many-Body Localization transition. By mapping the time-decorrelation of many-body wave-functions onto Directed Polymers in Random Media, we show the existence of a glass transition within the extended regime separating a metallic-like phase at small disorder, where delocalization occurs on an exponential number of paths, from a bad metal-like phase at intermediate disorder, where resonances are formed on rare, specific, disorder dependent site orbitals on very distant generations. The physical interpretation of subdiffusion and non-exponential relaxation emerging from this picture is complementary to the Griffiths one, although both scenarios rely on the presence of heavy-tailed distribution of the escape times. We relate the dynamical evolution in the glassy phase to the depinning transition of Directed Polymers, which results in macroscopic and abrupt jumps of the preferred delocalizing paths when a parameter like the energy is varied, and produce a singular behavior of the overlap correlation function between eigenstates at different energies. By comparing the quantum dynamics on loop-less Cayley trees and Random Regular Graphs we discuss the effect of loops, showing that in the latter slow dynamics and apparent power-laws extend on a very large time-window but are eventually cut-off on a time-scale that diverges at the MBL transition.

preprint2020arXiv

Attractive versus truncated repulsive supercooled liquids: The dynamics is encoded in the pair correlation function

We compare glassy dynamics in two liquids that differ in the form of their interaction potentials. Both systems have the same repulsive interactions but one has also an attractive part in the potential. These two systems exhibit very different dynamics despite having nearly identical pair correlation functions. We demonstrate that a properly weighted integral of the pair correlation function, which amplifies the subtle differences between the two systems, correctly captures their dynamical differences. The weights are obtained from a standard machine learning algorithm.

preprint2020arXiv

Complex Dynamics in Simple Neural Networks: Understanding Gradient Flow in Phase Retrieval

Despite the widespread use of gradient-based algorithms for optimizing high-dimensional non-convex functions, understanding their ability of finding good minima instead of being trapped in spurious ones remains to a large extent an open problem. Here we focus on gradient flow dynamics for phase retrieval from random measurements. When the ratio of the number of measurements over the input dimension is small the dynamics remains trapped in spurious minima with large basins of attraction. We find analytically that above a critical ratio those critical points become unstable developing a negative direction toward the signal. By numerical experiments we show that in this regime the gradient flow algorithm is not trapped; it drifts away from the spurious critical points along the unstable direction and succeeds in finding the global minimum. Using tools from statistical physics we characterize this phenomenon, which is related to a BBP-type transition in the Hessian of the spurious minima.

preprint2020arXiv

Double Trouble in Double Descent : Bias and Variance(s) in the Lazy Regime

Deep neural networks can achieve remarkable generalization performances while interpolating the training data perfectly. Rather than the U-curve emblematic of the bias-variance trade-off, their test error often follows a "double descent" - a mark of the beneficial role of overparametrization. In this work, we develop a quantitative theory for this phenomenon in the so-called lazy learning regime of neural networks, by considering the problem of learning a high-dimensional function with random features regression. We obtain a precise asymptotic expression for the bias-variance decomposition of the test error, and show that the bias displays a phase transition at the interpolation threshold, beyond which it remains constant. We disentangle the variances stemming from the sampling of the dataset, from the additive noise corrupting the labels, and from the initialization of the weights. Following up on Geiger et al. 2019, we first show that the latter two contributions are the crux of the double descent: they lead to the overfitting peak at the interpolation threshold and to the decay of the test error upon overparametrization. We then quantify how they are suppressed by ensemble averaging the outputs of K independently initialized estimators. When K is sent to infinity, the test error remains constant beyond the interpolation threshold. We further compare the effects of overparametrizing, ensembling and regularizing. Finally, we present numerical experiments on classic deep learning setups to show that our results hold qualitatively in realistic lazy learning scenarios.

preprint2020arXiv

Dynamical Mean-Field Theory and Aging Dynamics

Dynamical Mean-Field Theory (DMFT) replaces the many-body dynamical problem with one for a single degree of freedom in a thermal bath whose features are determined self-consistently. By focusing on models with soft disordered $p$-spin interactions, we show how to incorporate the mean-field theory of aging within dynamical mean-field theory. We study cases with only one slow time-scale, corresponding statically to the one-step replica symmetry breaking (1RSB) phase, and cases with an infinite number of slow time-scales, corresponding statically to the full replica symmetry breaking (FRSB) phase. For the former, we show that the effective temperature of the slow degrees of freedom is fixed by requiring critical dynamical behavior on short time-scales, i.e. marginality. For the latter, we find that aging on an infinite number of slow time-scales is governed by a stochastic equation where the clock for dynamical evolution is fixed by the change of effective temperature, hence obtaining a dynamical derivation of the stochastic equation at the basis of the FRSB phase. Our results extend the realm of the mean-field theory of aging to all situations where DMFT holds.

preprint2020arXiv

Finding the Needle in the Haystack with Convolutions: on the benefits of architectural bias

Despite the phenomenal success of deep neural networks in a broad range of learning tasks, there is a lack of theory to understand the way they work. In particular, Convolutional Neural Networks (CNNs) are known to perform much better than Fully-Connected Networks (FCNs) on spatially structured data: the architectural structure of CNNs benefits from prior knowledge on the features of the data, for instance their translation invariance. The aim of this work is to understand this fact through the lens of dynamics in the loss landscape. We introduce a method that maps a CNN to its equivalent FCN (denoted as eFCN). Such an embedding enables the comparison of CNN and FCN training dynamics directly in the FCN space. We use this method to test a new training protocol, which consists in training a CNN, embedding it to FCN space at a certain ``relax time'', then resuming the training in FCN space. We observe that for all relax times, the deviation from the CNN subspace is small, and the final performance reached by the eFCN is higher than that reachable by a standard FCN of same architecture. More surprisingly, for some intermediate relax times, the eFCN outperforms the CNN it stemmed, by combining the prior information of the CNN and the expressivity of the FCN in a complementary way. The practical interest of our protocol is limited by the very large size of the highly sparse eFCN. However, it offers interesting insights into the persistence of architectural bias under stochastic gradient dynamics. It shows the existence of some rare basins in the FCN loss landscape associated with very good generalization. These can only be accessed thanks to the CNN prior, which helps navigate the landscape during the early stages of optimization.

preprint2020arXiv

How to iron out rough landscapes and get optimal performances: Averaged Gradient Descent and its application to tensor PCA

In many high-dimensional estimation problems the main task consists in minimizing a cost function, which is often strongly non-convex when scanned in the space of parameters to be estimated. A standard solution to flatten the corresponding rough landscape consists in summing the losses associated to different data points and obtain a smoother empirical risk. Here we propose a complementary method that works for a single data point. The main idea is that a large amount of the roughness is uncorrelated in different parts of the landscape. One can then substantially reduce the noise by evaluating an empirical average of the gradient obtained as a sum over many random independent positions in the space of parameters to be optimized. We present an algorithm, called Averaged Gradient Descent, based on this idea and we apply it to tensor PCA, which is a very hard estimation problem. We show that Averaged Gradient Descent over-performs physical algorithms such as gradient descent and approximate message passing and matches the best algorithmic thresholds known so far, obtained by tensor unfolding and methods based on sum-of-squares.

preprint2020arXiv

Marvels and Pitfalls of the Langevin Algorithm in Noisy High-dimensional Inference

Gradient-descent-based algorithms and their stochastic versions have widespread applications in machine learning and statistical inference. In this work we perform an analytic study of the performances of one of them, the Langevin algorithm, in the context of noisy high-dimensional inference. We employ the Langevin algorithm to sample the posterior probability measure for the spiked matrix-tensor model. The typical behaviour of this algorithm is described by a system of integro-differential equations that we call the Langevin state evolution, whose solution is compared with the one of the state evolution of approximate message passing (AMP). Our results show that, remarkably, the algorithmic threshold of the Langevin algorithm is sub-optimal with respect to the one given by AMP. We conjecture this phenomenon to be due to the residual glassiness present in that region of parameters. Finally we show how a landscape-annealing protocol, that uses the Langevin algorithm but violate the Bayes-optimality condition, can approach the performance of AMP.

preprint2020arXiv

Role of fluctuations in the yielding transition of two-dimensional glasses

We numerically study yielding in two-dimensional glasses which are generated with a very wide range of stabilities by swap Monte-Carlo simulations and then slowly deformed at zero temperature. We provide strong numerical evidence that stable glasses yield via a nonequilibrium discontinuous transition in the thermodynamic limit. A critical point separates this brittle yielding from the ductile one observed in less stable glasses. We find that two-dimensional glasses yield similarly to their three-dimensional counterparts but display larger sample-to-sample disorder-induced fluctuations, stronger finite-size effects, and rougher spatial wandering of the observed shear bands. These findings strongly constrain effective theories of yielding.

preprint2020arXiv

Searching for the Gardner transition in glassy glycerol

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

preprint2020arXiv

Triple descent and the two kinds of overfitting: Where & why do they appear?

A recent line of research has highlighted the existence of a "double descent" phenomenon in deep learning, whereby increasing the number of training examples $N$ causes the generalization error of neural networks to peak when $N$ is of the same order as the number of parameters $P$. In earlier works, a similar phenomenon was shown to exist in simpler models such as linear regression, where the peak instead occurs when $N$ is equal to the input dimension $D$. Since both peaks coincide with the interpolation threshold, they are often conflated in the litterature. In this paper, we show that despite their apparent similarity, these two scenarios are inherently different. In fact, both peaks can co-exist when neural networks are applied to noisy regression tasks. The relative size of the peaks is then governed by the degree of nonlinearity of the activation function. Building on recent developments in the analysis of random feature models, we provide a theoretical ground for this sample-wise triple descent. As shown previously, the nonlinear peak at $N\!=\!P$ is a true divergence caused by the extreme sensitivity of the output function to both the noise corrupting the labels and the initialization of the random features (or the weights in neural networks). This peak survives in the absence of noise, but can be suppressed by regularization. In contrast, the linear peak at $N\!=\!D$ is solely due to overfitting the noise in the labels, and forms earlier during training. We show that this peak is implicitly regularized by the nonlinearity, which is why it only becomes salient at high noise and is weakly affected by explicit regularization. Throughout the paper, we compare analytical results obtained in the random feature model with the outcomes of numerical experiments involving deep neural networks.

preprint2020arXiv

Who is Afraid of Big Bad Minima? Analysis of Gradient-Flow in a Spiked Matrix-Tensor Model

Gradient-based algorithms are effective for many machine learning tasks, but despite ample recent effort and some progress, it often remains unclear why they work in practice in optimising high-dimensional non-convex functions and why they find good minima instead of being trapped in spurious ones. Here we present a quantitative theory explaining this behaviour in a spiked matrix-tensor model. Our framework is based on the Kac-Rice analysis of stationary points and a closed-form analysis of gradient-flow originating from statistical physics. We show that there is a well defined region of parameters where the gradient-flow algorithm finds a good global minimum despite the presence of exponentially many spurious local minima. We show that this is achieved by surfing on saddles that have strong negative direction towards the global minima, a phenomenon that is connected to a BBP-type threshold in the Hessian describing the critical points of the landscapes.

preprint2019arXiv

A jamming transition from under- to over-parametrization affects loss landscape and generalization

We argue that in fully-connected networks a phase transition delimits the over- and under-parametrized regimes where fitting can or cannot be achieved. Under some general conditions, we show that this transition is sharp for the hinge loss. In the whole over-parametrized regime, poor minima of the loss are not encountered during training since the number of constraints to satisfy is too small to hamper minimization. Our findings support a link between this transition and the generalization properties of the network: as we increase the number of parameters of a given model, starting from an under-parametrized network, we observe that the generalization error displays three phases: (i) initial decay, (ii) increase until the transition point --- where it displays a cusp --- and (iii) slow decay toward a constant for the rest of the over-parametrized regime. Thereby we identify the region where the classical phenomenon of over-fitting takes place, and the region where the model keeps improving, in line with previous empirical observations for modern neural networks.

preprint2019arXiv

Numerical implementation of dynamical mean field theory for disordered systems: application to the Lotka-Volterra model of ecosystems

Dynamical mean field theory (DMFT) is a tool that allows to analyze the stochastic dynamics of $N$ interacting degrees of freedom in terms of a self-consistent $1$-body problem. In this work, focusing on models of ecosystems, we present the derivation of DMFT through the dynamical cavity method, and we develop a method for solving it numerically. Our numerical procedure can be applied to a large variety of systems for which DMFT holds. We implement and test it for the generalized random Lotka-Volterra model, and show that complex dynamical regimes characterized by chaos and aging can be captured and studied by this framework.

preprint2019arXiv

Scaling description of generalization with number of parameters in deep learning

Supervised deep learning involves the training of neural networks with a large number $N$ of parameters. For large enough $N$, in the so-called over-parametrized regime, one can essentially fit the training data points. Sparsity-based arguments would suggest that the generalization error increases as $N$ grows past a certain threshold $N^{*}$. Instead, empirical studies have shown that in the over-parametrized regime, generalization error keeps decreasing with $N$. We resolve this paradox through a new framework. We rely on the so-called Neural Tangent Kernel, which connects large neural nets to kernel methods, to show that the initialization causes finite-size random fluctuations $\|f_{N}-\bar{f}_{N}\|\sim N^{-1/4}$ of the neural net output function $f_{N}$ around its expectation $\bar{f}_{N}$. These affect the generalization error $ε_{N}$ for classification: under natural assumptions, it decays to a plateau value $ε_{\infty}$ in a power-law fashion $\sim N^{-1/2}$. This description breaks down at a so-called jamming transition $N=N^{*}$. At this threshold, we argue that $\|f_{N}\|$ diverges. This result leads to a plausible explanation for the cusp in test error known to occur at $N^{*}$. Our results are confirmed by extensive empirical observations on the MNIST and CIFAR image datasets. Our analysis finally suggests that, given a computational envelope, the smallest generalization error is obtained using several networks of intermediate sizes, just beyond $N^{*}$, and averaging their outputs.

preprint2016arXiv

Breakdown of Elasticity in Amorphous Solids

What characterises a solid is its way to respond to external stresses. Ordered solids, such crystals, display an elastic regime followed by a plastic one, both well understood microscopically in terms of lattice distortion and dislocations. For amorphous solids the situation is instead less clear, and the microscopic understanding of the response to deformation and stress is a very active research topic. Several studies have revealed that even in the elastic regime the response is very jerky at low temperature, resembling very much the one of disordered magnetic materials. Here we show that in a very large class of amorphous solids this behaviour emerges by decreasing the temperature as a phase transition where standard elastic behaviour breaks down. At the transition all non-linear elastic modulii diverge and standard elasticity theory does not hold anymore. Below the transition the response to deformation becomes history and time-dependent.

preprint2016arXiv

Field-induced superdiffusion and dynamical heterogeneity

By analyzing two Kinetically Constrained Models of supercooled liquids we show that the anomalous transport of a driven tracer observed in supercooled liquids is another facet of the phenomenon of dynamical heterogeneity. We focus on the Fredrickson-Andersen and the Bertin-Bouchaud-Lequeux models. By numerical simulations and analytical arguments we demonstrate that the violation of the Stokes-Einstein relation and the observed field-induced superdiffusion have the same physical origin: while a fraction of probes do not move, others jump repeatedly because they are close to local mobile regions. The anomalous fluctuations observed out of equilibrium in presence of a pulling force $ε$, $σ_x^2(t) = \langle x_ε^2(t) \rangle - \langle x_ε(t) \rangle^2 \sim t^{3/2}$, which are accompanied by the asymptotic decay $α_ε(t)\sim t^{-1/2}$ of the non-Gaussian parameter from non-trivial values to zero, are due to the splitting of the probes population in the two (mobile and immobile) groups and to dynamical correlations, a mechanism expected to happen generically in supercooled liquids.

preprint2016arXiv

From non-ergodic eigenvectors to local resolvent statistics and back: a random matrix perspective

We study the statistics of the local resolvent and non-ergodic properties of eigenvectors for a generalised Rosenzweig-Porter $N\times N$ random matrix model, undergoing two transitions separated by a delocalised non-ergodic phase. Interpreting the model as the combination of on-site random energies $\{a_i\}$ and a structurally disordered hopping, we found that each eigenstate is delocalised over $N^{2-γ}$ sites close in energy $|a_j-a_i|\leq N^{1-γ}$ in agreement with Kravtsov \emph{et al}, arXiv:1508.01714. Our other main result, obtained combining a recurrence relation for the resolvent matrix with insights from Dyson's Brownian motion, is to show that the properties of the non-ergodic delocalised phase can be probed studying the statistics of the local resolvent in a non-standard scaling limit.

preprint2016arXiv

Level statistics and localization transitions of Lévy matrices

This work provide a thorough study of Lévy or heavy-tailed random matrices (LM). By analysing the self-consistent equation on the probability distribution of the diagonal elements of the resolvent we establish the equation determining the localisation transition and obtain the phase diagram of LMs. Using arguments based on super-symmetric field theory and Dyson Brownian motion we show that the eigenvalue statistics is the same one of the Gaussian Orthogonal Ensemble in the whole delocalised phase and is Poisson in the localised phase. Our numerics confirms these findings, valid in the limit of infinitely large LMs, but also reveals that the characteristic scale governing finite size effects diverges much faster than a power law approaching the transition and is already very large far from it. This leads to a very wide cross-over region in which the system looks as if it were in a mixed phase. Our results, together with the ones obtained previously, provide now a complete theory of Lévy matrices.

preprint2016arXiv

Spinodals with Disorder: from Avalanches in Random Magnets to Glassy Dynamics

We revisit the phenomenon of spinodals in the presence of quenched disorder and develop a complete theory for it. We focus on the spinodal of an Ising model in a quenched random field (RFIM), which has applications in many areas from materials to social science. By working at zero temperature in the quasi-statically driven RFIM, thermal fluctuations are eliminated and one can give a rigorous content to the notion of spinodal. We show that the latter is due to the depinning and the subsequent expansion of rare droplets. We work out the associated critical behavior, which, in any finite dimension, is very different from the mean-field one: the characteristic length diverges exponentially and the thermodynamic quantities display very mild nonanalyticities much like in a Griffith phenomenon. From the recently established connection between the spinodal of the RFIM and glassy dynamics, our results also allow us to conclusively assess the physical content and the status of the dynamical transition predicted by the mean-field theory of glass-forming liquids.

preprint2015arXiv

Nonperturbative fluctuations and metastability in a simple model: from observables to microscopic theory and back

Slow dynamics in glassy systems is often interpreted as due to thermally activated events between "metastable" states. This emphasizes the role of nonperturbative fluctuations, which is especially dramatic when these fluctuations destroy a putative phase transition predicted at the mean-field level. To gain insight into such hard problems, we consider the implementation of a generic back-and-forth process, between microscopic theory and observable behavior via effective theories, in a toy model that is simple enough to allow for a thorough investigation: the one-dimensional $φ^4$ theory at low temperature. We consider two ways of restricting the extent of the fluctuations, which both lead to a nonconvex effective potential (or free energy) : either through a finite-size system or by means of a running infrared cutoff within the nonperturbative Renormalization Group formalism. We discuss the physical insight one can get and the ways to treat strongly nonperturbative fluctuations in this context.

preprint2015arXiv

Slow Relaxations and Non-Equilibrium Dynamics in Classical and Quantum Systems

The aim of these lectures, given at the Les Houches Summer School of Physics "Strongly Interacting Quantum Systems Out of Equilibrium", is providing an introduction to several important and interesting facets of out of equilibrium dynamics. In recent years, there has been a boost in the research on quantum systems out of equilibrium. If fifteen years ago hard condensed matter and classical statistical physics remained rather separate research fields, now the focus on several kinds of out of equilibrium dynamics is making them closer and closer. The aim of my lectures was to present to the students the richness of this topic, insisting on the common concepts and showing that there is much to gain in considering and learning out of equilibrium dynamics as a whole research field.

preprint2015arXiv

Spin Glass in a Field: a New Zero-Temperature Fixed Point in Finite Dimensions

By using real space renormalisation group (RG) methods we show that spin-glasses in a field display a new kind of transition in high dimensions. The corresponding critical properties and the spin-glass phase are governed by two non-perturbative zero temperature fixed points of the RG flow. We compute the critical exponents, discuss the RG flow and its relevance for three dimensional systems. The new spin-glass phase we discovered has unusual properties, which are intermediate between the ones conjectured by droplet and full replica symmetry breaking theories. These results provide a new perspective on the long-standing debate about the behaviour of spin-glasses in a field.

preprint2015arXiv

The Super-Potts glass: a new disordered model for glass-forming liquids

We introduce a new disordered system, the Super-Potts model, which is a more frustrated version of the Potts glass. Its elementary degrees of freedom are variables that can take M values and are coupled via pair-wise interactions. Its exact solution on a completely connected lattice demonstrates that for large enough M it belongs to the class of mean-field systems solved by a one step replica symmetry breaking Ansatz. Numerical simulations by the parallel tempering technique show that in three dimensions it displays a phenomenological behaviour similar to the one of glass-forming liquids. The Super-Potts glass is therefore the first long-sought disordered model allowing one to perform extensive and detailed studies of the Random First Order Transition in finite dimensions. We also discuss its behaviour for small values of M, which is similar to the one of spin-glasses in a field.

preprint2014arXiv

Critical dynamical heterogeneities close to continuous second-order glass transitions

We analyse, using Inhomogenous Mode-Coupling Theory, the critical scaling behaviour of the dynamical susceptibility at a distance epsilon from continuous second-order glass transitions. We find that the dynamical correlation length xi behaves generically as epsilon^{-1/3} and that the upper critical dimension is equal to six. More surprisingly, we find activated dynamic scaling, where xi grows with time as [ln(t)]^2 exactly at criticality. All these results suggest a deep analogy between the glassy behaviour of attractive colloids or randomly pinned supercooled liquids and that of the Random Field Ising Model.

preprint2014arXiv

The Gardner transition in finite dimensions

Recent works on hard spheres in the limit of infinite dimensions revealed that glass states, envisioned as meta-basins in configuration space, can break up in a multitude of separate basins at low enough temperature or high enough pressure, leading to the emergence of new kinds of soft-modes and unusual properties. In this paper we study by perturbative renormalisation group techniques the critical properties of this transition, which has been discovered in disordered mean-field models in the '80s. We find that the upper critical dimension $d_u$ above which mean-field results hold is strictly larger than six and apparently non-universal, i.e. system dependent. Below $d_u$, we do not find any perturbative attractive fixed point (except for a tiny region of the 1RSB breaking parameter), thus showing that the transition in three dimensions either is governed by a non-perturbative fixed point unrelated to the Gaussian mean-field one or becomes first order or does not exist. We also discuss possible relationships with the behavior of spin glasses in a field.

preprint2013arXiv

Comparison of Static Length-Scales Characterizing the Glass Transition

The dramatic dynamic slowing down associated with the glass transition is considered by many to be related to the existence of a static length scale that grows when temperature decreases. Defining, identifying and measuring such a length is a subtle and non-trivial problem. Recently, two proposals, based on very different insights regarding the relevant physics, were put forward. One approach is based on the point-to-set correlation technique and the other on the scale where the lowest eigenvalue of the Hessian matrix becomes sensitive to disorder. In this Letter we present numerical evidence that the two approaches might result in the same identical length scale. This provides further mutual support to their relevance and, at the same time, raise interesting theoretical questions, discussed in the conclusion, concerning the fundamental reason for their relationship.

preprint2013arXiv

Dynamic Length Scales in Glass-Forming Liquids: A Inhomogeneous Molecular Dynamics Simulation Approach

In this work we numerically investigate a new method for the characterization of growing length scales associated with spatially heterogeneous dynamics of glass-forming liquids. This approach, motivated by the formulation of the inhomogeneous mode-coupling theory (IMCT) [Biroli G. et al. Phys. Rev. Lett. 2006 97, 195701], utilizes inhomogeneous molecular dynamics simulations in which the system is perturbed by a spatially modulated external potential. We show that the response of the two-point correlation function to the external field allows one to probe dynamic correlations. We examine the critical properties shown by this function, in particular the associated dynamic correlation length, that is found to be comparable to the one extracted from standardly-employed four-point correlation functions. Our numerical results are in qualitative agreement with IMCT predictions but suggest that one has to take into account fluctuations not included in this mean-field approach in order to reach quantitative agreement. Advantages of our approach over the more conventional one based on four-point correlation functions are discussed.

preprint2013arXiv

Perspective: The Glass Transition

We provide here a brief perspective on the glass transition field. It is an assessment, written from the point of view of theory, of where the field is and where it seems to be heading. We first give an overview of the main phenomenological characteristics, or "stylised facts", of the glass transition problem, i.e. the central observations that a theory of the physics of glass formation should aim to explain in a unified manner. We describe recent developments, with a particular focus on real space properties, including dynamical heterogeneity and facilitation, the search for underlying spatial or structural correlations, and the relation between the thermal glass transition and athermal jamming. We then discuss briefly how competing theories of the glass transition have adapted and evolved to account for such real space issues. We consider in detail two conceptual and methodological approaches put forward recently, that aim to access the fundamental critical phenomenon underlying the glass transition, be it thermodynamic or dynamic in origin, by means of biasing of ensembles, of configurations in the thermodynamic case, or of trajectories in the dynamic case. We end with a short outlook.

preprint2013arXiv

Quantum quenches, dynamical transitions and off-equilibrium quantum criticality

Several mean-field computations have revealed the existence of an out of equilibrium dynamical transition induced by quantum quenching an isolated system starting from its symmetry broken phase. In this work we focus on the quantum phi^4 N-component field theory. By taking into account dynamical fluctuations at the Hartree-Fock level, corresponding to the leading order of the 1/N expansion, we derive the critical properties of the dynamical transition beyond mean-field theory (including at finite temperature). We find diverging time and length-scales, dynamic scaling and aging. Finally, we unveil a relationship with coarsening, an off-equilibrium dynamical regime that can be induced by quenching from the symmetric toward the symmetry broken phase.

preprint2013arXiv

Random-field-like criticality in glass-forming liquids

We introduce an approach to derive an effective scalar field theory for the glass transition; the fluctuating field is the overlap between equilibrium configurations. We apply it to the case of constrained liquids for which the introduction of a conjugate source to the overlap field was predicted to lead to an equilibrium critical point. We show that the long-distance physics in the vicinity of this critical point is in the same universality class as that of a paradigmatic disordered model: the random-field Ising model. The quenched disorder is provided here by a reference equilibrium liquid configuration. We discuss to what extent this field-theoretical description and the mapping to the random field Ising model hold in the whole supercooled liquid regime, in particular near the glass transition.

preprint2012arXiv

A tentative replica theory of glassy Helium 4

We develop a quantum replica method for interacting particle systems and use it to estimate the location of the glass transition line in Helium 4. Although we do not fully succeed in taking into account all quantum effects, we make a thorough semiclassical analysis. We confirm previous suggestions that quantum fluctuations promote the formation of the glass and give a quantitative estimate of this effect at high density. Finally, we discuss the difficulties that are met when one tries to extend the calculation to the region of low densities and low temperatures, where quantum effects are strong and the semiclassical expansion breaks down.

preprint2012arXiv

Dynamic criticality at the jamming transition

We characterize vibrational motion occurring at low temperatures in dense suspensions of soft repulsive spheres over a broad range of volume fractions encompassing the jamming transition at (T = 0, phi = phi_J). We find that characteristic time and length scales of thermal vibrations obey critical scaling in the vicinity of the jamming transition. We show in particular that the amplitude and the time scale of dynamic fluctuations diverge symmetrically on both sides of the transition, and directly reveal a diverging correlation length. The critical region near phi_J is divided in three different regimes separated by a characteristic temperature scale T*(phi) that vanishes quadratically with the distance to phi_J. While two of them, (T < T*(phi), phi > phi_J) and (T < T*(phi), phi < phi_J), are described by harmonic theories developed in the zero temperature limit, the third one for T > T*(phi) is inherently anharmonic and displays new critical properties. We find that the quadratic scaling of T*(phi) is due to nonperturbative anharmonic contributions, its amplitude being orders of magnitude smaller than the perturbative prediction based on the expansion to quartic order in the interactions. Our results show that thermal vibrations in colloidal assemblies directly reveal the critical nature of the jamming transition. The critical region, however, is very narrow and has not yet been attained experimentally, even in recent specifically-dedicated experiments.

preprint2012arXiv

Finite size effects in the dynamics of glass-forming liquids

We present a comprehensive theoretical study of finite size effects in the relaxation dynamics of glass-forming liquids. Our analysis is motivated by recent theoretical progress regarding the understanding of relevant correlation length scales in liquids approaching the glass transition. We obtain predictions both from general theoretical arguments and from a variety of specific perspectives: mode-coupling theory, kinetically constrained and defect models, and random first order transition theory. In the latter approach, we predict in particular a non-monotonic evolution of finite size effects across the mode-coupling crossover due to the competition between mode-coupling and activated relaxation. We study the role of competing relaxation mechanisms in giving rise to non-monotonic finite size effects by devising a kinetically constrained model where the proximity to the mode-coupling singularity can be continuously tuned by changing the lattice topology. We use our theoretical findings to interpret the results of extensive molecular dynamics studies of four model liquids with distinct structures and kinetic fragilities. While the less fragile model only displays modest finite size effects, we find a more significant size dependence evolving with temperature for more fragile models, such as Lennard-Jones particles and soft spheres. Finally, for a binary mixture of harmonic spheres we observe the predicted non-monotonic temperature evolution of finite size effects near the fitted mode-coupling singularity, suggesting that the crossover from mode-coupling to activated dynamics is more pronounced for this model. Finally, we discuss the close connection between our results and the recent report of a non-monotonic temperature evolution of a dynamic length scale near the mode-coupling crossover in harmonic spheres.

preprint2012arXiv

Patch-repetition correlation length in glassy systems

We obtain the patch-repetition entropy Sigma within the Random First Order Transition theory (RFOT) and for the square plaquette system, a model related to the dynamical facilitation theory of glassy dynamics. We find that in both cases the entropy of patches of linear size l, Sigma(l), scales as s_c l^d+A l^{d-1} down to length-scales of the order of one, where A is a positive constant, s_c is the configurational entropy density and d the spatial dimension. In consequence, the only meaningful length that can be defined from patch-repetition is the cross-over length xi=A/s_c. We relate xi to the typical length-scales already discussed in the literature and show that it is always of the order of the largest static length. Our results provide new insights, which are particularly relevant for RFOT theory, on the possible real space structure of super-cooled liquids. They suggest that this structure differs from a mosaic of different patches having roughly the same size.

preprint2012arXiv

Random Pinning Glass Transition: Hallmarks, Mean-Field Theory and Renormalization Group Analysis

We present a detailed analysis of glass transitions induced by pinning particles at random from an equilibrium configuration. We first develop a mean-field analysis based on the study of p-spin spherical disordered models and then obtain the three dimensional critical behavior by the Migdal-Kadanoff real space renormalization group method. We unveil the important physical differences with the case in which particles are pinned from a random (or very high temperature) configuration. We contrast the pinning particles approach to the ones based on biasing dynamical trajectories with respect to their activity and on coupling to equilibrium configurations. Finally, we discuss numerical and experimental tests.

preprint2012arXiv

Social interaction, noise and antibiotic-mediated switches in the intestinal microbiota

The intestinal microbiota plays important roles in digestion and resistance against entero-pathogens. As with other ecosystems, its species composition is resilient against small disturbances but strong perturbations such as antibiotics can affect the consortium dramatically. Antibiotic cessation does not necessarily restore pre-treatment conditions and disturbed microbiota are often susceptible to pathogen invasion. Here we propose a mathematical model to explain how antibiotic-mediated switches in the microbiota composition can result from simple social interactions between antibiotic-tolerant and antibiotic-sensitive bacterial groups. We build a two-species (e.g. two functional-groups) model and identify regions of domination by antibiotic-sensitive or antibiotic-tolerant bacteria, as well as a region of multistability where domination by either group is possible. Using a new framework that we derived from statistical physics, we calculate the duration of each microbiota composition state. This is shown to depend on the balance between random fluctuations in the bacterial densities and the strength of microbial interactions. The singular value decomposition of recent metagenomic data confirms our assumption of grouping microbes as antibiotic-tolerant or antibiotic-sensitive in response to a single antibiotic. Our methodology can be extended to multiple bacterial groups and thus it provides an ecological formalism to help interpret the present surge in microbiome data.

preprint2012arXiv

The Generalized Arrhenius Law in Out of Equilibrium Systems

In this work we provide a comprehensive analysis of the activation problem out of equilibrium. We generalize the Arrhenius law for systems driven by non conservative time independent forces, subjected to retarded friction and non-Markovian noise. The role of the energy function is now played by the out of equilibrium potential ϕ = -lim_{T{\to}0} T log P_s, with P_s being the steady state probability distribution and T the strength of the noise. We unveil the relationship between the generalized Arrhenius law and a time-reversal transformation discussed in the context of fluctuations theorems out of equilibrium. Moreover, we characterize the noise-activated trajectories by obtaining their explicit expressions and identifying their irreversible nature. Finally, we discuss a real biological application that illustrates our results.

preprint2011arXiv

Aging and relaxation near Random Pinning Glass Transitions

Pinning particles at random in supercooled liquids is a promising route to make substantial progress on the glass transition problem. Here we develop a mean-field theory by studying the equilibrium and non-equilibrium dynamics of the spherical p-spin model in presence of a fraction c of pinned spins. Our study shows the existence of two dynamic critical lines: one corresponding to usual Mode Coupling transitions and the other one to dynamic spinodal transitions. Quenches in the portion of the c - T phase diagram delimited by those two lines leads to aging. By extending our results to finite dimensional systems we predict non-interrupted aging only for quenches on the ideal glass transition line and two very different types of equilibrium relaxations for quenches below and above it.

preprint2011arXiv

Ideal Glass Transitions by Random Pinning

We study the effect of freezing the positions of a fraction $c$ of particles from an equilibrium configuration of a supercooled liquid at a temperature $T$. We show that within the Random First-Order Transition theory pinning particles leads to an ideal glass transition for a critical fraction $c=c_{K}(T)$ even for moderate super-cooling, e.g. close to the Mode-Coupling transition temperature. We first derive the phase diagram in the $T-c$ plane by mean field approximations. Then, by applying a real-space renormalization group method, we obtain the critical properties for $|c-c_{K}(T)|\rightarrow 0$, in particular the divergence of length and time scales. These are dominated by two zero-temperature fixed points. We also show that for $c=c_{K}(T)$ the typical distance between frozen particles is related to the static point-to-set lengthscale of the unconstrained liquid. We discuss what are the main differences when particles are frozen in other geometries and not from an equilibrium configuration. Finally, we explain why the glass transition induced by freezing particles provides a new and very promising avenue of research to probe the glassy state and ascertain, or disprove, the validity of the theories of the glass transition.

preprint2011arXiv

Leggett's bound for amorphous solids

We investigate the constraints on the superfluid fraction of an amorphous solid following from an upper bound derived by Leggett. In order to accomplish this, we use as input density profiles generated for amorphous solids in a variety of different manners including by investigating Gaussian fluctuations around classical results. These rough estimates suggest that, at least at the level of the upper bound, there is not much difference in terms of superfluidity between a glass and a crystal characterized by the same Lindemann ratio. Moreover, we perform Path Integral Monte Carlo simulations of distinguishable Helium 4 rapidly quenched from the liquid phase to very lower temperature, at the density of the freezing transition. We find that the system crystallizes very quickly, without any sign of intermediate glassiness. Overall our results suggest that the experimental observations of large superfluid fractions in Helium 4 after a rapid quench correspond to samples evolving far from equilibrium, instead of being in a stable glass phase. Other scenarios and comparisons to other results on the super-glass phase are also discussed.

preprint2011arXiv

Theoretical perspective on the glass transition and amorphous materials

We provide a theoretical perspective on the glass transition in molecular liquids at thermal equilibrium, on the spatially heterogeneous and aging dynamics of disordered materials, and on the rheology of soft glassy materials. We start with a broad introduction to the field and emphasize its connections with other subjects and its relevance. The important role played by computer simulations to study and understand the dynamics of systems close to the glass transition at the molecular level is spelled out. We review the recent progress on the subject of the spatially heterogeneous dynamics that characterizes structural relaxation in materials with slow dynamics. We then present the main theoretical approaches describing the glass transition in supercooled liquids, focusing on theories that have a microscopic, statistical mechanics basis. We describe both successes and failures, and critically assess the current status of each of these approaches. The physics of aging dynamics in disordered materials and the rheology of soft glassy materials are then discussed, and recent theoretical progress is described. For each section, we give an extensive overview of the most recent advances, but we also describe in some detail the important open problems that, we believe, will occupy a central place in this field in the coming years.

preprint2010arXiv

Anderson model on Bethe lattices: density of states, localization properties and isolated eigenvalue

We revisit the Anderson localization problem on Bethe lattices, putting in contact various aspects which have been previously only discussed separately. For the case of connectivity 3 we compute by the cavity method the density of states and the evolution of the mobility edge with disorder. Furthermore, we show that below a certain critical value of the disorder the smallest eigenvalue remains delocalized and separated by all the others (localized) ones by a gap. We also study the evolution of the mobility edge at the center of the band with the connectivity, and discuss the large connectivity limit.

preprint2010arXiv

Coarsening of Disordered Quantum Rotors under a Bias Voltage

We solve the dynamics of an ensemble of interacting rotors coupled to two leads at different chemical potential letting a current flow through the system and driving it out of equilibrium. We show that at low temperature the coarsening phase persists under the voltage drop up to a critical value of the applied potential that depends on the characteristics of the electron reservoirs. We discuss the properties of the critical surface in the temperature, voltage, strength of quantum fluctuations and coupling to the bath phase diagram. We analyze the coarsening regime finding, in particular, which features are essentially quantum mechanical and which are basically classical in nature. We demonstrate that the system evolves via the growth of a coherence length with the same time-dependence as in the classical limit, $R(t) \simeq t^{1/2}$ -- the scalar curvature driven universality class. We obtain the scaling function of the correlation function at late epochs in the coarsening regime and we prove that it coincides with the classical one once a prefactor that encodes the dependence on all the parameters is factorized. We derive a generic formula for the current flowing through the system and we show that, for this model, it rapidly approaches a constant that we compute.

preprint2010arXiv

Elementary Excitation Modes in a Granular Glass above Jamming

The dynamics of granular media in the jammed, glassy region is described in terms of "modes", by applying a Principal Component Analysis (PCA) to the covariance matrix of the position of individual grains. We first demonstrate that this description is justified and gives sensible results in a regime of time/densities such that a metastable state can be observed on long enough timescale to define the reference configuration. For small enough times/system sizes, or at high enough packing fractions, the spectral properties of the covariance matrix reveals large, collective fluctuation modes that cannot be explained by a Random Matrix benchmark where these correlations are discarded. We then present a first attempt to find a link between the softest modes of the covariance matrix during a certain "quiet" time interval and the spatial structure of the rearrangement event that ends this quiet period. The motion during these cracks is indeed well explained by the soft modes of the dynamics before the crack, but the number of cracks preceded by a "quiet" period strongly reduces when the system unjams, questioning the relevance of a description in terms of modes close to the jamming transition, at least for frictional grains.

preprint2010arXiv

Kibble-Zurek mechanism and infinitely slow annealing through critical points

We revisit the Kibble-Zurek mechanism by analyzing the dynamics of phase ordering systems during an infinitely slow annealing across a second order phase transition. We elucidate the time and cooling rate dependence of the typical growing length and we use it to predict the number of topological defects left over in the symmetry broken phase as a function of time, both close and far from the critical region. Our results extend the Kibble-Zurek mechanism and reveal its limitations.

preprint2010arXiv

Overview of different characterisations of dynamic heterogeneity

Dynamic heterogeneity is now recognised as a central aspect of structural relaxation in disordered materials with slow dynamics, and was the focus of intense research in the last decade. Here we describe how initial, indirect observations of dynamic heterogeneity have recently evolved into well-defined, quantitative, statistical characterisations, in particular through the use of high-order correlation and response functions. We highlight both recent progress and open questions about the characterisation of dynamic heterogeneity in glassy materials. We also discuss the limits of available tools and describe a few candidates for future research in order to gain deeper understanding of the origin and nature of glassiness in disordered systems.

preprint2010arXiv

Quantum quenches and off-equilibrium dynamical transition in the infinite-dimensional Bose-Hubbard model

We study the off-equilibrium dynamics of the infinite dimensional Bose Hubbard Model after a quantum quench. The dynamics can be analyzed exactly by mapping it to an effective Newtonian evolution. For integer filling, we find a dynamical transition separating regimes of small and large quantum quenches starting from the superfluid state. This transition is very similar to the one found for the fermionic Hubbard model by mean field approximations.

preprint2010arXiv

Renormalization group analysis of the random first order transition

We consider the approach describing glass formation in liquids as a progressive trapping in an exponentially large number of metastable states. To go beyond the mean-field setting, we provide a real-space renormalization group (RG) analysis of the associated replica free-energy functional. The present approximation yields in finite dimensions an ideal glass transition similar to that found in mean field. However, we find that along the RG flow the properties associated with metastable glassy states, such as the configurational entropy, are only defined up to a characteristic length scale that diverges as one approaches the ideal glass transition. The critical exponents characterizing the vicinity of the transition are the usual ones associated with a first-order discontinuity fixed point.

preprint2010arXiv

Statistical properties of the spectrum the extended Bose-Hubbard model

Motivated by the role that spectral properties play for the dynamical evolution of a quantum many-body system, we investigate the level spacing statistic of the extended Bose-Hubbard model. In particular, we focus on the distribution of the ratio of adjacent level spacings, useful at large interaction, to distinguish between chaotic and non-chaotic regimes. After revisiting the bare Bose-Hubbard model, we study the effect of two different perturbations: next-nearest neighbor hopping and nearest-neighbor interaction. The system size dependence is investigated together with the effect of the proximity to integrable points or lines. Lastly, we discuss the consequences of a cutoff in the number of onsite bosons onto the level statistics.

preprint2010arXiv

Symmetries of generating functionals of Langevin processes with colored multiplicative noise

We present a comprehensive study of the symmetries of the generating functionals of generic Langevin processes with multiplicative colored noise. We treat both Martin-Siggia-Rose-Janssen-deDominicis and supersymmetric formalisms. We summarize the relations between observables that they imply including fluctuation relations, fluctuation-dissipation theorems, and Schwinger-Dyson equations. Newtonian dynamics and their invariances follow in the vanishing friction limit.

preprint2009arXiv

Dynamical Heterogeneity in Lattice Glass Models

In this paper we consider in detail the properties of dynamical heterogeneity in lattice glass models (LGMs). LGMs are lattice models whose dynamical rules are based on thermodynamic, as opposed to purely kinetic, considerations. We devise a LGM that is not prone to crystallization and displays properties of a fragile glass-forming liquid. Particle motion in this model tends to be locally anisotropic on intermediate time scales even though the rules governing the model are isotropic. The model demonstrates violations of the Stokes-Einstein relation and the growth of various length scales associated with dynamical heterogeneity. We discuss future avenues of research comparing the predictions of lattice glass models and kinetically constrained models to atomistic systems.

preprint2006arXiv

Critical fluctuations and breakdown of Stokes-Einstein relation in the Mode-Coupling Theory of glasses

We argue that the critical dynamical fluctuations predicted by the mode-coupling theory (MCT) of glasses provide a natural mechanism to explain the breakdown of the Stokes-Einstein relation. This breakdown, observed numerically and experimentally in a region where MCT should hold, is one of the major difficulty of the theory, for which we propose a natural resolution based on the recent interpretation of the MCT transition as a bona fide critical point with a diverging length scale. We also show that the upper critical dimension of MCT is d_c=8.

preprint2006arXiv

Jamming percolation and glassy dynamics

We present a detailed physical analysis of the dynamical glass-jamming transition which occurs for the so called Knight models recently introduced and analyzed in a joint work with D.S.Fisher \cite{letterTBF}. Furthermore, we review some of our previous works on Kinetically Constrained Models. The Knights models correspond to a new class of kinetically constrained models which provide the first example of finite dimensional models with an ideal glass-jamming transition. This is due to the underlying percolation transition of particles which are mutually blocked by the constraints. This jamming percolation has unconventional features: it is discontinuous (i.e. the percolating cluster is compact at the transition) and the typical size of the clusters diverges faster than any power law when $ρ\nearrowρ_c$. These properties give rise for Knight models to an ergodicity breaking transition at $ρ_c$: at and above $ρ_{c}$ a finite fraction of the system is frozen. In turn, this finite jump in the density of frozen sites leads to a two step relaxation for dynamic correlations in the unjammed phase, analogous to that of glass forming liquids. Also, due to the faster than power law divergence of the dynamical correlation length, relaxation times diverge in a way similar to the Vogel-Fulcher law.

preprint2006arXiv

On the top eigenvalue of heavy-tailed random matrices

We study the statistics of the largest eigenvalue lambda_max of N x N random matrices with unit variance, but power-law distributed entries, P(M_{ij})~ |M_{ij}|^{-1-mu}. When mu > 4, lambda_max converges to 2 with Tracy-Widom fluctuations of order N^{-2/3}. When mu < 4, lambda_max is of order N^{2/mu-1/2} and is governed by Fréchet statistics. The marginal case mu=4 provides a new class of limiting distribution that we compute explicitely. We extend these results to sample covariance matrices, and show that extreme events may cause the largest eigenvalue to significantly exceed the Marcenko-Pastur edge. Connections with Directed Polymers are briefly discussed.

preprint2004arXiv

What do we learn from the shape of the dynamical susceptibility of glass-formers?

We compute analytically and numerically the four-point correlation function that characterizes non-trivial cooperative dynamics in glassy systems within several models of glasses: elasto-plastic deformations, mode-coupling theory (MCT), collectively rearranging regions (CRR), diffusing defects and kinetically constrained models (KCM). Some features of the four-point susceptibility chi_4(t) are expected to be universal. at short times we expect an elastic regime characterized by a t or sqrt{t} growth. We find both in the beta, and the early alpha regime that chi_4 sim t^mu, where mu is directly related to the mechanism responsible for relaxation. This regime ends when a maximum of chi_4 is reached at a time t=t^* of the order of the relaxation time of the system. This maximum is followed by a fast decay to zero at large times. The height of the maximum also follows a power-law, chi_4(t^*) sim t^{*lambda}. The value of the exponents mu and lambda allows one to distinguish between different mechanisms. For example, freely diffusing defects in d=3 lead to mu=2 and lambda=1, whereas the CRR scenario rather predicts either mu=1 or a logarithmic behaviour depending on the nature of the nucleation events, and a logarithmic behaviour of chi_4(t^*). MCT leads to mu=b and lambda =1/gamma, where b and gamma are the standard MCT exponents. We compare our theoretical results with numerical simulations on a Lennard-Jones and a soft-sphere system. Within the limited time-scales accessible to numerical simulations, we find that the exponent mu is rather small, mu < 1, with a value in reasonable agreement with the MCT predictions.

preprint2000arXiv

From inherent structures to pure states: some simple remarks and examples

The notions of pure states and inherent structures, i.e. stable configurations against 1-spin flip are discussed. We explain why these different concepts accidentally coincide in mean-field models with infinite connectivity and present an exactly solvable unidimensional model where they do not. At zero temperature pure states are to some extent related to k-spin flip stable configurations with k going to infinity after the thermodynamical limit has been taken. This relationship is supported by an explicit analysis of the TAP equations and calculation of the number of pure states and k-spin flips stable configurations in a mean-field model with finite couplings. Finally we discuss the relevance of the concepts of pure states and inherent structures in finite dimensional glassy systems.

preprint2000arXiv

Metastable states in glassy systems

Truly stable metastable states are an artifact of the mean-field approximation or the zero temperature limit. If such appealing concepts in glass theory as configurational entropy are to have a meaning beyond these approximations, one needs to cast them in a form involving states with finite lifetimes. Starting from elementary examples and using results of Gaveau and Schulman, we propose a simple expression for the configurational entropy and revisit the question of taking flat averages over metastable states. The construction is applicable to finite dimensional systems, and we explicitly show that for simple mean-field glass models it recovers, justifies and generalises the known results. The calculation emphasises the appearance of new dynamical order parameters.