BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Matthieu Wyart

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

ResearcherUnclaimed source record
cond-mat.softcond-mat.stat-mechcond-mat.dis-nncond-mat.mtrl-sciMachine LearningBiological Physicsmath-phmath.MPnlin.AOphysics.flu-dynphysics.geo-phQuantitative Methods

Catalog footprint

What is connected

45works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

45 published item(s)

preprint2022arXiv

Avalanches and deformation in glasses and disordered systems

In this chapter, we discuss avalanches in glasses and disordered systems, and the macroscopic dynamical behavior that they mediate. We briefly review three classes of systems where avalanches are observed: depinning transition of disordered interfaces, yielding of amorphous materials, and the jamming transition. Without extensive formalism, we discuss results gleaned from theoretical approaches -- mean-field theory, scaling and exponent relations, the renormalization group, and a few results from replica theory. We focus both on the remarkably sophisticated physics of avalanches and on relatively new approaches to the macroscopic flow behavior exhibited past the depinning/yielding transition.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Failure and success of the spectral bias prediction for Kernel Ridge Regression: the case of low-dimensional data

Recently, several theories including the replica method made predictions for the generalization error of Kernel Ridge Regression. In some regimes, they predict that the method has a `spectral bias': decomposing the true function $f^*$ on the eigenbasis of the kernel, it fits well the coefficients associated with the O(P) largest eigenvalues, where $P$ is the size of the training set. This prediction works very well on benchmark data sets such as images, yet the assumptions these approaches make on the data are never satisfied in practice. To clarify when the spectral bias prediction holds, we first focus on a one-dimensional model where rigorous results are obtained and then use scaling arguments to generalize and test our findings in higher dimensions. Our predictions include the classification case $f(x)=$sign$(x_1)$ with a data distribution that vanishes at the decision boundary $p(x)\sim x_1^χ$. For $χ>0$ and a Laplace kernel, we find that (i) there exists a cross-over ridge $λ^*_{d,χ}(P)\sim P^{-\frac{1}{d+χ}}$ such that for $λ\gg λ^*_{d,χ}(P)$, the replica method applies, but not for $λ\llλ^*_{d,χ}(P)$, (ii) in the ridge-less case, spectral bias predicts the correct training curve exponent only in the limit $d\rightarrow\infty$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Mean-field description for the architecture of low-energy excitations in glasses

In amorphous materials, groups of particles can rearrange locally into a new stable configuration. Such elementary excitations are key as they determine the response to external stresses, as well as to thermal and quantum fluctuations. Yet, understanding what controls their geometry remains a challenge. Here we build a scaling description of the geometry and energy of low-energy excitations in terms of the distance to an instability, as predicted for instance at the dynamical transition in mean field approaches of supercooled liquids. We successfully test our predictions in ultrastable computer glasses, with a gapped and ungapped (regular) spectrum. Overall, our approach explains why excitations become less extended, with a higher energy and displacement scale upon cooling.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Asymptotic learning curves of kernel methods: empirical data v.s. Teacher-Student paradigm

How many training data are needed to learn a supervised task? It is often observed that the generalization error decreases as $n^{-β}$ where $n$ is the number of training examples and $β$ an exponent that depends on both data and algorithm. In this work we measure $β$ when applying kernel methods to real datasets. For MNIST we find $β\approx 0.4$ and for CIFAR10 $β\approx 0.1$, for both regression and classification tasks, and for Gaussian or Laplace kernels. To rationalize the existence of non-trivial exponents that can be independent of the specific kernel used, we study the Teacher-Student framework for kernels. In this scheme, a Teacher generates data according to a Gaussian random field, and a Student learns them via kernel regression. With a simplifying assumption -- namely that the data are sampled from a regular lattice -- we derive analytically $β$ for translation invariant kernels, using previous results from the kriging literature. Provided that the Student is not too sensitive to high frequencies, $β$ depends only on the smoothness and dimension of the training data. We confirm numerically that these predictions hold when the training points are sampled at random on a hypersphere. Overall, the test error is found to be controlled by the magnitude of the projection of the true function on the kernel eigenvectors whose rank is larger than $n$. Using this idea we predict relate the exponent $β$ to an exponent $a$ describing how the coefficients of the true function in the eigenbasis of the kernel decay with rank. We extract $a$ from real data by performing kernel PCA, leading to $β\approx0.36$ for MNIST and $β\approx0.07$ for CIFAR10, in good agreement with observations. We argue that these rather large exponents are possible due to the small effective dimension of the data.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Direct Coupling Analysis of Epistasis in Allosteric Materials

In allosteric proteins, the binding of a ligand modifies function at a distant active site. Such allosteric pathways can be used as target for drug design, generating considerable interest in inferring them from sequence alignment data. Currently, different methods lead to conflicting results, in particular on the existence of long-range evolutionary couplings between distant amino-acids mediating allostery. Here we propose a resolution of this conundrum, by studying epistasis and its inference in models where an allosteric material is evolved in silico to perform a mechanical task. We find in our model the four types of epistasis (Synergistic, Sign, Antagonistic, Saturation), which can be both short or long-range and have a simple mechanical interpretation. We perform a Direct Coupling Analysis (DCA) and find that DCA predicts well the cost of point mutations but is a rather poor generative model. Strikingly, it can predict short-range epistasis but fails to capture long-range epistasis, in consistence with empirical findings. We propose that such failure is generic when function requires subparts to work in concert. We illustrate this idea with a simple model, which suggests that other methods may be better suited to capture long-range effects.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Inferring the flow properties of epithelial tissues from their geometry

Amorphous materials exhibit complex material proprteties with strongly nonlinear behaviors. Below a yield stress they behave as plastic solids, while they start to yield above a critical stress $Σ_c$. A key quantity controlling plasticity which is, however, hard to measure is the density $P(x)$ of weak spots, where $x$ is the additional stress required for local plastic failure. In the thermodynamic limit $P(x)\sim x^θ$ is singular at $x= 0$ in the solid phase below the yield stress $Σ_c$. This singularity is related to the presence of system spannig avalanches of plastic events. Here we address the question if the density of weak spots and the flow properties of a material can be determined from the geometry of an amporphous structure alone. We show that a vertex model for cell packings in tissues exhibits the phenomenology of plastic amorphous systems. As the yield stress is approached from above, the strain rate vanishes and the avalanches size $S$ and their duration $τ$ diverge. We then show that in general, in materials where the energy functional depend on topology, the value $x$ is proportional to the length $L$ of a bond that vanishes in a plastic event. For this class of models $P(x)$ is therefore readily measurable from geometry alone. Applying this approach to a quantification of the cell packing geometry in the developing wing epithelium of the fruit fly, we find that in this tissue $P(L)$ exhibits a power law with exponents similar to those found numerically for a vertex model in its solid phase. This suggests that this tissue exhibits plasticity and non-linear material properties that emerge from collective cell behaviors and that these material properties govern developmental processes. Our approach based on the relation between topology and energetics suggests a new route to outstanding questions associated with the yielding transition.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Jamming with tunable roughness

We introduce a new model to study the effect of surface roughness on the jamming transition. By performing numerical simulations, we show that for a smooth surface, the jamming transition density and the contact number at the transition point both increase upon increasing asphericity, as for ellipsoids and spherocylinders. Conversely, for a rough surface, both quantities decrease, in quantitative agreement with the behavior of frictional particles. Furthermore, in the limit corresponding to the Coulomb friction law, the model satisfies a generalized isostaticity criterion proposed in previous studies. We introduce a counting argument that justifies this criterion and interprets it geometrically. Finally, we propose a simple theory to predict the contact number at finite friction from the knowledge of the force distribution in the infinite friction limit.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Perspective: A Phase Diagram for Deep Learning unifying Jamming, Feature Learning and Lazy Training

Deep learning algorithms are responsible for a technological revolution in a variety of tasks including image recognition or Go playing. Yet, why they work is not understood. Ultimately, they manage to classify data lying in high dimension -- a feat generically impossible due to the geometry of high dimensional space and the associated curse of dimensionality. Understanding what kind of structure, symmetry or invariance makes data such as images learnable is a fundamental challenge. Other puzzles include that (i) learning corresponds to minimizing a loss in high dimension, which is in general not convex and could well get stuck bad minima. (ii) Deep learning predicting power increases with the number of fitting parameters, even in a regime where data are perfectly fitted. In this manuscript, we review recent results elucidating (i,ii) and the perspective they offer on the (still unexplained) curse of dimensionality paradox. We base our theoretical discussion on the $(h,α)$ plane where $h$ is the network width and $α$ the scale of the output of the network at initialization, and provide new systematic measures of performance in that plane for MNIST and CIFAR 10. We argue that different learning regimes can be organized into a phase diagram. A line of critical points sharply delimits an under-parametrised phase from an over-parametrized one. In over-parametrized nets, learning can operate in two regimes separated by a smooth cross-over. At large initialization, it corresponds to a kernel method, whereas for small initializations features can be learnt, together with invariants in the data. We review the properties of these different phases, of the transition separating them and some open questions. Our treatment emphasizes analogies with physical systems, scaling arguments and the development of numerical observables to quantitatively test these results empirically.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Thermal origin of quasi-localised excitations in glasses

Key aspects of glasses are controlled by the presence of excitations in which a group of particles can rearrange. Surprisingly, recent observations indicate that their density is dramatically reduced and their size decreases as the temperature of the supercooled liquid is lowered. Some theories predict these excitations to cause a gap in the spectrum of quasi-localised modes of the Hessian that grows upon cooling, while others predict a pseudo-gap ${D_L(ω)} \sim ω^α$. To unify these views and observations, we generate glassy configurations of controlled gap magnitude $ω_c$ at temperature ${T=0}$, using so-called `breathing' particles, and study how such gapped states respond to thermal fluctuations. We find that \textit{(i)}~the gap always fills up at finite $T$ with ${D_L(ω) \approx A_4(T) \, ω^4}$ and ${A_4 \sim \exp(-E_a / T)}$ at low $T$, \textit{(ii)}~$E_a$ rapidly grows with $ω_c$, in reasonable agreement with a simple scaling prediction ${E_a\sim ω_c^4}$ and \textit{(iii)}~at larger $ω_c$ excitations involve fewer particles, as we rationalise, and eventually become string-like. We propose an interpretation of mean-field theories of the glass transition, in which the modes beyond the gap act as an excitation reservoir, from which a pseudo-gap distribution is populated with its magnitude rapidly decreasing at lower $T$. We discuss how this picture unifies the rarefaction as well as the decreasing size of excitations upon cooling, together with a string-like relaxation occurring near the glass transition.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Infinitesimal asphericity changes the universality of the jamming transition

The jamming transition of non-spherical particles is fundamentally different from the spherical case. Non-spherical particles are hypostatic at their jamming points, while isostaticity is ensured in the case of the jamming of spherical particles. This structural difference implies that the presence of asphericity affects the critical exponents related to the contact number and the vibrational density of states. Moreover, while the force and gap distributions of isostatic jamming present power-law behaviors, even an infinitesimal asphericity is enough to smooth out these singularities. In a recent work [PNAS 115(46), 11736], we have used a combination of marginal stability arguments and the replica method to explain these observations. We argued that systems with internal degrees of freedom, like the rotations in ellipsoids, or the variation of the radii in the case of the \textit{breathing} particles fall in the same universality class. In this paper, we review comprehensively the results about the jamming with internal degrees of freedom in addition to the translational degrees of freedom. We use a variational argument to derive the critical exponents of the contact number, shear modulus, and the characteristic frequencies of the density of states. Moreover, we present additional numerical data supporting the theoretical results, which were not shown in the previous work.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Effect of particle collisions in dense suspension flows

We study non-local effects associated with particle collisions in dense suspension flows, in the context of the affine solvent model known to capture various aspects of the jamming transition. We show that an individual collision changes significantly the velocity field on a characteristic volume $Ω_c\sim 1/δz$ that diverges as jamming is approached, where $δz$ is the deficit in coordination number required to jam the system. Such an event also affects the contact forces between particles on that same volume $Ω_c$, but this change is modest in relative terms, of order $f_{coll}\sim \bar{f}^{0.8}$, where $\bar{f}$ is the typical contact force scale. We then show that the requirement that coordination is stationary (such that a collision has a finite probability to open one contact elsewhere in the system) yields the scaling of the viscosity (or equivalently the viscous number) with coordination deficit $δz$. The same scaling result was derived in [E.~DeGiuli, G.~Düring, E.~Lerner, and M.~Wyart, Phys.~Rev.~E {\bf 91}, 062206 (2015)] via different arguments making an additional assumption. The present approach gives a mechanistic justification as to why the correct finite size scaling volume behaves as $1/δz$, and can be used to recover a marginality condition known to characterize the distributions of contact forces and gaps in jammed packings.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Evidence for marginal stability in emulsions

We report the first measurements of the effect of pressure on vibrational modes in emulsions, which serve as a model for soft frictionless spheres at zero temperature. As a function of the applied pressure, we find that the density of states D(omega) exhibits a low-frequency cutoff omega*, which scales linearly with the number of extra contacts per particle dz. Moreover, for omega<omega*, D(omega)~ omega^2/omega*^2; a quadratic behavior whose prefactor is larger than what is expected from Debye theory. This surprising result agrees with recent theoretical findings. Finally, the degree of localization of the softest low frequency modes increases with compression, as shown by the participation ratio as well as their spatial configurations. Overall, our observations show that emulsions are marginally stable and display non-plane-wave modes up to vanishing frequencies.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On Variational Arguments for Vibrational Modes near Jamming

Amorphous solids tend to present an abundance of soft elastic modes, which diminish their transport properties, generate heterogeneities in their elastic response, and affect non-linear processes like thermal activation of plasticity. This is especially true in packings of particles near their jamming transition, for which effective medium theory and variational arguments can both predict the density of vibrational modes. However, recent numerics support that one hypothesis of the variational argument does not hold. We provide a novel variational argument which overcomes this problem, and correctly predicts the scaling properties of soft modes near the jamming transition. Soft modes are shown to be related to the response to a local strain in more connected networks, and to be characterized by a volume $1/δz$, where $δz$ is the excess coordination above the Maxwell threshold. These predictions are verified numerically.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Scale-free channeling patterns near the onset of erosion of sheared granular beds

Erosion shapes our landscape and occurs when a sufficient shear stress is exerted by a fluid on a sedimented layer. What controls erosion at a microscopic level remains debated, especially near the threshold forcing where it stops. Here we study experimentally the collective dynamics of the moving particles, using a set-up where the system spontaneously evolves toward the erosion onset. We find that the spatial organization of the erosion flux is heterogeneous in space, and occurs along channels of local flux $σ$ whose distribution displays scaling near threshold and follows $P(σ)\sim J/σ$, where $J$ is the mean erosion flux. Channels are strongly correlated in the direction of forcing but not in the transverse direction. We show that these results quantitatively agree with a model where the dynamics is governed by the competition of disorder (which channels mobile particles) and particle interactions (which reduces channeling). These observations support that for laminar flows, erosion is a dynamical phase transition which shares similarity with the plastic depinning transition occurring in dirty superconductors. The methodology we introduce here could be applied to probe these systems as well.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Scaling Description of Non-Local Rheology

Non-locality is crucial to understand the plastic flow of an amorphous material, and has been successfully described by the fluidity, along with a cooperativity length scale ξ. We demonstrate, by applying the scaling hypothesis to the yielding transition, that non-local effects in non-uniform stress configurations can be explained within the framework of critical phenomena. From the scaling description, scaling relations between different exponents are derived, and collapses of strain rate profiles are made both in shear driven and pressure driven flow. We find that the cooperative length in non-local flow is governed by the same correlation length in finite dimensional homogeneous flow, excluding the mean field exponents. We also show that non-locality also affects the finite size scaling of the yield stress, especially the large finite size effects observed in pressure driven flow. Our theoretical results are nicely verified by the elasto-plastic model, and experimental data.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Unsteady flow and particle migration in dense, non-Brownian suspensions

We present experimental results on dense corn-starch suspensions as examples of non-Brownian, nearly-hard particles that undergo continuous and discontinuous shear thickening (CST and DST) at intermediate and high densities respectively. Our results offer strong support for recent theories involving a stress-dependent effective contact friction among particles. We show however that in the DST regime, where theory might lead one to expect steady-state shear bands oriented layerwise along the vorticity axis, the real flow is unsteady. To explain this, we argue that steady-state banding is generically ruled out by the requirement that, for hard non-Brownian particles, the solvent pressure and the normal-normal component of the particle stress must balance separately across the interface between bands. (Otherwise there is an unbalanced migration flux.) However, long-lived transient shear-bands remain possible.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A model for the erosion onset of a granular bed sheared by a viscous fluid

We study theoretically the erosion threshold of a granular bed forced by a viscous fluid. We first introduce a novel model of interacting particles driven on a rough substrate. It predicts a continuous transition at some threshold forcing $θ_c$, beyond which the particle current grows linearly $J\sim θ-θ_c$, in agreement with experiments. The stationary state is reached after a transient time $t_{\rm conv}$ which diverges near the transition as $t_{\rm conv}\sim |θ-θ_c|^{-z}$ with $z\approx 2.5$. The model also makes quantitative testable predictions for the drainage pattern: the distribution $P(σ)$ of local current is found to be extremely broad with $P(σ)\sim J/σ$, spatial correlations for the current are negligible in the direction transverse to forcing, but long-range parallel to it. We explain some of these features using a scaling argument and a mean-field approximation that builds an analogy with $q$-models. We discuss the relationship between our erosion model and models for the depinning transition of vortex lattices in dirty superconductors, where our results may also apply.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Adaptive Elastic Networks as models of supercooled liquids

The thermodynamics and dynamics of supercooled liquids correlate with their elasticity. In particular for covalent networks, the jump of specific heat is small and the liquid is {\it strong} near the threshold valence where the network acquires rigidity. By contrast, the jump of specific heat and the fragility are large away from this threshold valence. In a previous work [Proc. Natl. Acad. Sci. U.S.A., 110, 6307 (2013)], we could explain these behaviors by introducing a model of supercooled liquids in which local rearrangements interact via elasticity. However, in that model the disorder characterizing elasticity was frozen, whereas it is itself a dynamic variable in supercooled liquids. Here we study numerically and theoretically adaptive elastic network models where polydisperse springs can move on a lattice, thus allowing for the geometry of the elastic network to fluctuate and evolve with temperature. We show numerically that our previous results on the relationship between structure and thermodynamics hold in these models. We introduce an approximation where redundant constraints (highly coordinated regions where the frustration is large) are treated as an ideal gas, leading to analytical predictions that are accurate in the range of parameters relevant for real materials. Overall, these results lead to a description of supercooled liquids, in which the distance to the rigidity transition controls the number of directions in phase space that cost energy and the specific heat.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Criticality in the approach to failure in amorphous solids

Failure of amorphous solids is fundamental to various phenomena, including landslides and earthquakes. Recent experiments indicate that highly plastic regions form elongated structures that are especially apparent near the maximal shear stress $Σ_{\max}$ where failure occurs. This observation suggested that $Σ_{\max}$ acts as a critical point where the length scale of those structures diverges, possibly causing macroscopic transient shear bands. Here we argue instead that the entire solid phase ($Σ<Σ_{\max}$) is critical, that plasticity always involves system-spanning events, and that their magnitude diverges at $Σ_{\max}$ independently of the presence of shear bands. We relate the statistics and fractal properties of these rearrangements to an exponent $θ$ that captures the stability of the material, which is observed to vary continuously with stress, and we confirm our predictions in elastoplastic models.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dynamics and Correlations among Soft Excitations in Marginally Stable Glasses

Marginal stability is the notion that stability is achieved, but only barely so. This property constrains the ensemble of configurations explored at low temperature in a variety of systems, including spin, electron and structural glasses. A key feature of marginal states is a (saturated) pseudo-gap in the distribution of soft excitations. We study how such a pseudo-gap appears dynamically in the case of the Sherrington-Kirkpatrick (SK) spin glass. After revisiting and correcting the multi-spin-flip criterion for local stability, we show that stationarity along the hysteresis loop requires that soft spins are frustrated among each other, with a correlation that diverges as $C(λ)\sim 1/λ$, where $λ$ is the larger of two considered local fields. We explain how this arises spontaneously in a marginal system and develop an analogy between the spin dynamics in the SK model and random walks in two dimensions. We discuss the applicability of these findings to hard sphere packings.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Effects of coordination and pressure on sound attenuation, boson peak and elasticity in amorphous solids

Connectedness and applied stress strongly affect elasticity in solids. In various amorphous solids, mechanical stability can be lost either by reducing connectedness or by increasing pressure. We present an effective medium theory of elasticity that extends previous approaches by incorporating the effect of compression, of amplitude $e$, allowing one to describe quantitative features of sound propagation, transport, the boson peak, and elastic moduli near the elastic instability occurring at a compression $e_c$. The theory disentangles several frequencies characterizing the vibrational spectrum: the onset frequency $ω_0\sim \sqrt{e_c-e}$ where strongly-scattered modes appear in the vibrational spectrum, the pressure-independent frequency $ω_*$ where the density of states displays a plateau, the boson peak frequency $ω_{BP}$, and the Ioffe-Regel frequency $ω_{IR}$ where scattering length and wavelength become equal. We predict that sound attenuation crosses over from $ω^4$ to $ω^2$ behaviour at $ω_0$. We predict that a frequency-dependent length scale $l_s(ω)$ and speed of sound $ν(ω)$ characterize vibrational modes, and could be extracted from scattering data. One key result is the prediction of a flat diffusivity above $ω_0$, in agreement with previously unexplained observations. We find that the shear modulus does not vanish at the elastic instability, but drops by a factor of 2. We check our predictions in packings of soft particles and study the case of covalent networks and silica. Overall, our approach unifies sound attenuation, transport and length scales entering elasticity in a single framework where disorder is not the main parameter controlling the boson peak, in agreement with observations. This framework leads to a phase diagram where various glasses can be placed, connecting microscopic structure to vibrational properties.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Evolution of covalent networks under cooling: contrasting the rigidity window and jamming scenarios

We study the evolution of structural disorder under cooling in supercooled liquids, focusing on covalent networks. We introduce a model for the energy of networks that incorporates weak non-covalent interactions. We show that at low-temperature, these interactions considerably affect the network topology near the rigidity transition that occurs as the coordination increases. As a result, this transition becomes mean-field and does not present a line of critical points previously argued for, the "rigidity window". Vibrational modes are then not fractons, but instead are similar to the anomalous modes observed in packings of particles near jamming. These results suggest an alternative interpretation for the intermediate phase observed in chalcogenides.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Length scales and self-organization in dense suspension flows

Dense non-Brownian suspension flows of hard particles display mystifying properties: as the jamming threshold is approached, the viscosity diverges, as well as a length scale that can be identified from velocity correlations. To unravel the microscopic mechanism governing dissipation and its connection to the observed long-range correlations, we develop an analogy between suspension flows and the rigidity transition occurring when floppy networks are pulled -- a transition believed to be associated to the stress-stiffening of certain gels. After deriving the critical properties near the rigidity transition, we show numerically that suspensions flows lie close to it. We find that this proximity causes a decoupling between viscosity and the correlation length of velocities ξ, which scales as the length l_c characterizing the response of the velocity in flow to a local perturbation, previously predicted to follow l_c\sim 1/\sqrt{z_c-z}\sim p^{0.18} where p is the dimensionless particle pressure, z the coordination of the contact network made by the particles and z_c is twice the spatial dimension. We confirm these predictions numerically, predict the existence of a larger length scale l_r\sim 1/\sqrt{p} with mild effects on velocity correlation and the existence of a vanishing strain δγ\sim 1/p that characterizes de-correlation in flow.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Marginal Stability in Structural, Spin and Electron Glasses

We revisit the concept of marginal stability in glasses, and determine its range of applicability in the context of avalanche-type response to slow external driving. We argue that there is an intimate connection between a pseudo-gap in the distribution of local fields and crackling in systems with long-range interactions. We show how the principle of marginal stability offers a unifying perspective on the phenomenology of systems as diverse as spin and electron glasses, hard spheres, pinned elastic interfaces and the plasticity of soft amorphous solids.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Phonon gap and localization lengths in floppy materials

Gels of semi-flexible polymers, network glasses made of low valence elements, softly compressed ellipses and dense suspensions under flow are examples of floppy materials. These systems present collective motions with almost no restoring force. We study theoretically and numerically the frequency-dependence of the response of these materials, and the length scales that characterize their elasticity. We show that isotropic floppy elastic networks present a phonon gap for frequencies smaller than a frequency $ω^*$ governed by coordination, and that the elastic response is localized on a length scale $l_c\sim 1/\sqrt{ω^*}$ that diverges as the phonon gap vanishes (with a logarithmic correction in the two dimensional case). $l_c$ also characterizes velocity correlations under shear, whereas another length scale $l^*\sim 1/ω^*$ characterizes the effect of pinning boundaries on elasticity. We discuss the implications of our findings for suspensions flows, and the correspondence between floppy materials and amorphous solids near unjamming, where $l_c$ and $l^*$ have also been identified but where their roles are not fully understood.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Scaling description of the yielding transition in soft amorphous solids at zero temperature

Yield stress materials flow if a sufficiently large shear stress is ap- plied. Although such materials are ubiquitous and relevant for indus- try, there is no accepted microscopic description of how they yield, even in the simplest situations where temperature is negligible and where flow inhomogeneities such as shear bands or fractures are ab- sent. Here we propose a scaling description of the yielding transition in amorphous solids made of soft particles at zero temperature. Our description makes a connection between the Herschel-Bulkley expo- nent characterizing the singularity of the flow curve near the yield stress Σc, the extension and duration of the avalanches of plasticity observed at threshold, and the density P(x) of soft spots, or shear transformation zones, as a function of the stress increment x be- yond which they yield. We argue that the critical exponents of the yielding transition can be expressed in terms of three independent exponents θ, df and z, characterizing respectively the density of soft spots, the fractal dimension of the avalanches, and their duration. Our description shares some similarity with the depinning transition that occurs when an elastic manifold is driven through a random potential, but also presents some striking differences. We test our arguments in an elasto-plastic model, an automaton model similar to those used in depinning, but with a different interaction kernel, and find satisfying agreement with our predictions both in two and three dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Breakdown of continuum elasticity in amorphous solids

We show numerically that the response of simple amorphous solids (elastic networks and particle packings) to a local force dipole is characterized by a lengthscale $\ell_c$ that diverges as unjamming is approached as $\ell_c \sim (z - 2d)^{-1/2}$, where $z \ge 2d$ is the mean coordination, and $d$ is the spatial dimension, at odds with previous numerical claims. We also show how the magnitude of the lengthscale $\ell_c$ is amplified by the presence of internal stresses in the disordered solid. Our data suggests a divergence of $\ell_c\sim (p_c-p)^{-1/4}$ with proximity to a critical internal stress $p_c$ at which soft elastic modes become unstable.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Discontinuous shear thickening without inertia in dense non-Brownian suspensions

A consensus is emerging that discontinuous shear thickening (DST) in dense suspensions marks a transition from a flow state where particles remain well separated by lubrication layers, to one dominated by frictional contacts. We show here that reasonable assumptions about contact proliferation predict two distinct types of DST in the absence of inertia. The first occurs at densities above the jamming point of frictional particles; here the thickened state is completely jammed and (unless particles deform) cannot flow without inhomogeneity or fracture. The second regime shows strain- rate hysteresis and arises at somewhat lower densities where the thickened phase flows smoothly. DST is predicted to arise when finite-range repulsions defer contact formation until a characteristic stress level is exceeded.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Low-energy non-linear excitations in sphere packings

We study theoretically and numerically how hard frictionless particles in random packings can rearrange. We demonstrate the existence of two distinct unstable non-linear modes of rearrangement, both associated with the opening and the closing of contacts. Mode one, whose density is characterized by some exponent θ', corresponds to motions of particles extending throughout the entire system. Mode two, whose density is characterized by an exponent θ != θ', corresponds to the local buckling of a few particles. Mode one is shown to yield at a much higher rate than mode two when a stress is applied. We show that the distribution of contact forces follows P(f) f^{min(θ',θ)}, and that imposing that the packing cannot be densified further leads to the bounds γ >= 1/(2+θ') and γ >= (1-θ)/2, where γ characterizes the singularity of the pair distribution function g(r) at contact. These results extend the theoretical analysis of [M. Wyart, Phys. Rev. Lett 109, 125502 (2012)] where the existence of mode two was not considered. We perform numerics that support that these bounds are saturated with γ \approx 0.38, θ \approx 0.17 and θ' \approx 0.44. We measure systematically the stability of all such modes in packings, and confirm their marginal stability. The principle of marginal stability thus allows to make clearcut predictions on the ensemble of configurations visited in these out-of-equilibrium systems, and on the contact forces and pair distribution functions. It also reveals the excitations that need to be included in a description of plasticity or flow near jamming, and suggests a new path to study two-level systems and soft spots in simple amorphous solids of repulsive particles.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the density of shear transformation zones in amorphous solids

We study the stability of amorphous solids, focusing on the distribution P(x) of the local stress increase x that would lead to an instability. We argue that this distribution is singular P(x)x^θ, where the exponent θ is non-zero if the elastic interaction between rearranging regions is non-monotonic, and increases with the interaction range. For a class of finite dimensional models we show that stability implies a lower bound on θ, which is found to lie near saturation. For quadrupolar interactions these models yield θ ~ 0.6 for d=2 and θ~ 0.4 in d=3 where d is the spatial dimension, accurately capturing previously unresolved observations in atomistic models, both in quasi-static flow and after a fast quench.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Simulations of driven overdamped frictionless hard spheres

We introduce an event-driven simulation scheme for overdamped dynamics of frictionless hard spheres subjected to external forces, neglecting hydrodynamic interactions. Our event-driven approach is based on an exact equation of motion which relates the driving force to the resulting velocities through the geometric information characterizing the underlying network of contacts between the hard spheres. Our method allows for a robust extraction of the instantaneous coordination of the particles as well as contact force statistics and dynamics, under any chosen driving force, in addition to shear flow and compression. It can also be used for generating high-precision jammed packings under shear, compression, or both. We present a number of additional applications of our method.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Why glass elasticity affects the thermodynamics and fragility of super-cooled liquids

Super-cooled liquids are characterized by their fragility: the slowing down of the dynamics under cooling is more sudden and the jump of specific heat at the glass transition is generally larger in fragile liquids than in strong ones. Despite the importance of this quantity in classifying liquids, explaining what aspects of the microscopic structure controls fragility remains a challenge. Surprisingly, experiments indicate that the linear elasticity of the glass -- a purely local property of the free energy landscape -- is a good predictor of fragility. In particular, materials presenting a large excess of soft elastic modes, the so-called boson peak, are strong. This is also the case for network liquids near the rigidity percolation, known to affect elasticity. Here we introduce a model of the glass transition based on the assumption that particles can organize locally into distinct configurations, which are coupled spatially via elasticity. The model captures the mentioned observations connecting elasticity and fragility. We find that materials presenting an abundance of soft elastic modes have little elastic frustration: energy is insensitive to most directions in phase space, leading to a small jump of specific heat. In this framework strong liquids turn out to lie the closest to a critical point associated with a rigidity or jamming transition, and their thermodynamic properties are related to the problem of number partitioning and to Hopfield nets in the limit of small memory.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A unified framework for non-Brownian suspension flows and soft amorphous solids

While the rheology of non-Brownian suspensions in the dilute regime is well-understood, their behavior in the dense limit remains mystifying. As the packing fraction of particles increases, particle motion becomes more collective, leading to a growing length scale and scaling properties in the rheology as the material approaches the jamming transition. There is no accepted microscopic description of this phenomenon. However, in recent years it has been understood that the elasticity of simple amorphous solids is governed by a critical point, the unjamming transition where the pressure vanishes, and where elastic properties display scaling and a diverging length scale. The correspondence between these two transitions is at present unclear. Here we show that for a simple model of dense flow, which we argue captures the essential physics near the jamming threshold, a formal analogy can be made between the rheology of the flow and the elasticity of simple networks. This analogy leads to a new conceptual framework to relate microscopic structure to rheology. It enables us to define and compute numerically normal modes and a density of states. We find striking similarities between the density of states in flow, and that of amorphous solids near unjamming: both display a plateau above some frequency scale ω* ~ |z_c-z|, where z is the coordination of the network of particles in contact, z_c = 2D where D is the spatial dimension. However, a spectacular difference appears: the density of states in flow displays a single mode at another frequency scale ω_{min} << ω* governing the divergence of the viscosity.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Stability at Random Close Packing

The requirement that packings of hard particles, arguably the simplest structural glass, cannot be compressed by rearranging their network of contacts is shown to yield a new constraint on their microscopic structure. This constraint takes the form a bound between the distribution of contact forces P(f) and the pair distribution function g(r): if P(f) \sim f^θ and g(r) \sim (r-σ)^(-γ), where σ is the particle diameter, one finds that γ \geq 1/(2+θ). This bound plays a role similar to those found in some glassy materials with long-range interactions, such as the Coulomb gap in Anderson insulators or the distribution of local fields in mean-field spin glasses. There is ground to believe that this bound is saturated, offering an explanation for the presence of avalanches of rearrangements with power-law statistics observed in packings.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Toward a microscopic description of flow near the jamming threshold

We study the relationship between microscopic structure and viscosity in non-Brownian suspensions. We argue that the formation and opening of contacts between particles in flow effectively leads to a negative selection of the contacts carrying weak forces. We show that an analytically tractable model capturing this negative selection correctly reproduces scaling properties of flows near the jamming transition. In particular, we predict that (i) the viscosity η diverges with the coordination z as η ~ (z_c-z)^{-(3+θ)/(1+θ)}, (ii) the operator that governs flow displays a low-frequency mode that controls the divergence of viscosity, at a frequency ω_min\sim(z_c-z)^{(3+θ)/(2+2θ)}, and (iii) the distribution of forces displays a scale f* that vanishes near jamming as f*/<f>\sim(z_c-z)^{1/(1+θ)} where θ characterizes the distribution of contact forces P(f)\simf^θ at jamming, and where z_c is the Maxwell threshold for rigidity.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The jamming scenario - an introduction and outlook

The jamming scenario of disordered media, formulated about 10 years ago, has in recent years been advanced by analyzing model systems of granular media. This has led to various new concepts that are increasingly being explored in in a variety of systems. This chapter contains an introductory review of these recent developments and provides an outlook on their applicability to different physical systems and on future directions. The first part of the paper is devoted to an overview of the findings for model systems of frictionless spheres, focussing on the excess of low-frequency modes as the jamming point is approached. Particular attention is paid to a discussion of the cross-over frequency and length scales that govern this approach. We then discuss the effects of particle asphericity and static friction, the applicability to bubble models for wet foams in which the friction is dynamic, the dynamical arrest in colloids, and the implications for molecular glasses.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Correlations between vibrational entropy and dynamics in super-cooled liquids

A relation between vibrational entropy and particles mean square displacement is derived in super-cooled liquids, assuming that the main effect of temperature changes is to rescale the vibrational spectrum. Deviations from this relation, in particular due to the presence of a Boson Peak whose shape and frequency changes with temperature, are estimated. Using observations of the short-time dynamics in liquids of various fragility, it is argued that (i) if the crystal entropy is significantly smaller than the liquid entropy at $T_g$, the extrapolation of the vibrational entropy leads to the correlation $T_K\approx T_0$, where $T_K$ is the Kauzmann temperature and $T_0$ is the temperature extracted from the Vogel-Fulcher fit of the viscosity. (ii) The jump in specific heat associated with vibrational entropy is very small for strong liquids, and increases with fragility. The analysis suggests that these correlations stem from the stiffening of the Boson Peak under cooling, underlying the importance of this phenomenon on the dynamical arrest.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Scaling of phononic transport with connectivity in amorphous solids

The effect of coordination on transport is investigated theoretically using random networks of springs as model systems. An effective medium approximation is made to compute the density of states of the vibrational modes, their energy diffusivity (a spectral measure of transport) and their spatial correlations as the network coordination $z$ is varied. Critical behaviors are obtained as $z\to z_c$ where these networks lose rigidity. A sharp cross-over from a regime where modes are plane-wave-like toward a regime of extended but strongly-scattered modes occurs at some frequency $ω^*\sim z-z_c$, which does not correspond to the Ioffe-Regel criterion. Above $ω^*$ both the density of states and the diffusivity are nearly constant. These results agree remarkably with recent numerical observations of repulsive particles near the jamming threshold \cite{ning}. The analysis further predicts that the length scale characterizing the correlation of displacements of the scattered modes decays as $1/\sqrtω$ with frequency, whereas for $ω<<ω^*$ Rayleigh scattering is found with a scattering length $l_s\sim (z-z_c)^3/ω^4$. It is argued that this description applies to silica glass where it compares well with thermal conductivity data, and to transverse ultrasound propagation in granular matter.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Elasticity of soft particles and colloids near Random Close Packing

Assemblies of purely repulsive and frictionless particles, such as emulsions or hard spheres, display very curious properties near their jamming transition, which occurs at the random close packing for mono-disperse spheres. Although such systems do not contain the long and cross-linked polymeric chains characterizing a rubber, they behave macroscopically in a similar way: the shear modulus $G$ can become infinitely smaller than the bulk modulus $B$. After reviewing recent theoretical results on the structure of such packing (in particular their coordination) I will propose an explanation for the observed scaling of the elastic moduli, and explain why the arguments both apply to soft and hard particles.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Heat transport in model jammed solids

We calculate numerically the normal modes of vibrations in 3D jammed packings of soft spheres as a function of the packing fraction and obtain the energy diffusivity, a spectral measure of transport that controls sound propagation and thermal conductivity. The crossover frequency between weak and strong phonon scattering is controlled by the coordination and shifts to zero as the system is decompressed towards the critical packing fraction at which rigidity is lost. Below the crossover, the diffusivity displays a power-law divergence with inverse frequency, which suggests that the vibrational modes are primarily transverse waves, weakly scattered by disorder. Above it, a large number of modes appear whose diffusivity plateaus at a nearly constant value independent of the inter-particle potential, before dropping to zero above the Anderson localization frequency. The thermal conductivity of a marginally jammed solid just above the rigidity threshold is calculated and related to the one measured experimentally at room temperature for most glasses.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

On the dependence of the avalanche angle on the granular layer thickness

A layer of sand of thickness h flows down a rough surface if the inclination is larger than some threshold value theta which decreases with h. A tentative microscopic model for the dependence of theta with h is proposed for rigid frictional grains, based on the following hypothesis: (i) a horizontal layer of sand has some coordination z larger than a critical value z_c where mechanical stability is lost (ii) as the tilt angle is increased, the configurations visited present a growing proportion $_s of sliding contacts. Instability with respect to flow occurs when z-z_s=z_c. This criterion leads to a prediction for theta(h) in good agreement with empirical observations.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Rigidity-based approach to the boson peak in amorphous solids: from sphere packing to amorphous silica

Glasses have a large excess of low-frequency vibrational modes in comparison with continuous elastic body, the so-called Boson Peak, which appears to correlate with several crucial properties of glasses, such as transport or fragility. I review recent results showing that the Boson Peak is a necessary consequence of the weak connectivity of the solid. I explain why in assemblies repulsive spheres the boson peak shifts up to zero frequency as the pressure is lowered toward the jamming threshold, and derive the corresponding exponent. I show how these ideas capture the main low-frequency features of the vibrational spectrum of amorphous silica. These results extend arguments of Phillips on the presence of floppy modes in under-constrained covalent networks to glasses where the covalent network is rigid, or when interactions are purely radial.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access