Theory of the jamming transition at finite temperature

Unified understanding of the impact of semiflexibility, concentration, and molecular weight on macromolecular-scale ring diffusion

Conformationally fluctuating, globally compact macromolecules such as polymeric rings, single-chain nanoparticles, microgels, and many-arm stars display complex dynamic behaviors due to their rich topological structure and intermolecular organization. Synthetic rings are hybrid objects with conformations that display both ideal random walk and compact globular features, which can serve as models of genomic DNA. To date, emphasis has been placed on the effect of ring molecular weight on their unusual behaviors. Here, we combine simulations and a microscopic force-level theory to build a unified understanding for how key aspects of ring dynamics depend on different tunable molecular properties including backbone rigidity, monomer concentration, degree of traditional entanglement, and molecular weight. Our large-scale molecular dynamics simulations of ring melts with very different backbone stiffnesses reveal unanticipated behaviors which agree well with our generalized theory. This includes a universal master curve for center-of-mass diffusion constants as a function of molecular weight scaled by a chemistry and thermodynamic state-dependent critical molecular weight that generalizes the concept of an entanglement cross-over for linear chains. The key physics is how backbone rigidity and monomer concentration induced changes of the entanglement length, interring packing, degree of interpenetration, and liquid compressibility slow down space-time dynamic-force correlations on macromolecular scales. A power law decay of the center-of-mass diffusion constant with inverse molecular weight squared is the first consequence, followed by an ultraslow activated hopping transport regime. Our results set the stage to address slow dynamics and kinetic arrest in different families of compact synthetic and biological polymeric systems.

Local vs. cooperative: Unraveling glass transition mechanisms with SEER

Which phenomenon slows down the dynamics in supercooled liquids and turns them into glasses is a long-standing question of condensed matter. Most popular theories posit that as the temperature decreases, many events must occur in a coordinated fashion on a growing length scale for relaxation to occur. Instead, other approaches consider that local barriers associated with the elementary rearrangement of a few particles or "excitations" govern the dynamics. To resolve this conundrum, our central result is to introduce an algorithm, Systematic Excitation ExtRaction, which can systematically extract hundreds of excitations and their energy from any given configuration. We also provide a measurement of the activation energy, characterizing the liquid dynamics, based on fast quenching and reheating. We use these two methods in a popular liquid model of polydisperse particles. Such polydisperse models are known to capture the hallmarks of the glass transition and can be equilibrated efficiently up to millisecond time scales. The analysis reveals that cooperative effects do not control the fragility of such liquids: the change of energy of local barriers determines the change of activation energy. More generally, these methods can now be used to measure the degree of cooperativity of any liquid model.

Effective elastic moduli of composites with a strongly disordered host material

The local elastic properties of strongly disordered material are investigated using the theory of correlated random matrices. A significant increase in stiffness is shown in the interfacial region, the thickness of which depends on the strength of disorder. It is shown that this effect plays a crucial role in nanocomposites, in which interfacial regions are formed around each nanoparticle. The studied interfacial effect can significantly increase the influence of nanoparticles on the macroscopic stiffness of nanocomposites. The obtained thickness of the interfacial region is determined by the heterogeneity lengthscale and is of the same order as the lengthscale of the boson peak.

Hard-sphere jamming through the lens of linear optimization

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

Long-Range Anomalous Decay of the Correlation in Jammed Packings

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

Finite-size effects in the microscopic critical properties of jammed configurations: A comprehensive study of the effects of different types of disorder

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

Emergent solidity of amorphous materials as a consequence of mechanical self-organisation

Amorphous solids have peculiar properties distinct from crystals. One of the most fundamental mysteries is the emergence of solidity in such nonequilibrium, disordered state without the protection by long-range translational order. A jammed system at zero temperature, although marginally stable, has solidity stemming from the space-spanning force network, which gives rise to the long-range stress correlation. Here, we show that such nonlocal correlation already appears at the nonequilibrium glass transition upon cooling. This is surprising since we also find that the system suffers from giant anharmonic fluctuations originated from the fractal-like potential energy landscape. We reveal that it is the percolation of the force-bearing network that allows long-range stress transmission even under such circumstance. Thus, the emergent solidity of amorphous materials is a consequence of nontrivial self-organisation of the disordered mechanical architecture. Our findings point to the significance of understanding amorphous solids and nonequilibrium glass transition from a mechanical perspective.

Vibrational Properties of Hard and Soft Spheres Are Unified at Jamming

The unconventional thermal properties of jammed amorphous solids are directly related to their density of vibrational states. While the vibrational spectrum of jammed soft sphere solids has been fully described, the vibrational spectrum of hard spheres, a model glass former often related to physical colloidal glasses, is still unknown due to the difficulty of treating the nonanalytic interaction potential. We bypass this difficulty using the recently described effective interaction potential for the free energy of thermal hard spheres. By minimizing this effective free energy, we mimic the rapid compression of hard spheres and produce typical configurations of the thermal system. We measure the resulting vibrational spectrum and characterize its evolution toward the jamming point where configurations of hard and soft spheres are trivially unified. For densities approaching jamming from below, we observe low-frequency modes which agree with those found in numerical simulations of jammed soft spheres. Our measurements of the vibrational structure demonstrate that the jamming universality extends away from jamming: hard sphere thermal systems below jamming exhibit the same vibrational spectra as thermal and athermal soft sphere systems above the transition.

Soft-grain compression: Beyond the jamming point

We present the experimental studies of highly strained soft bidisperse granular systems made of hyperelastic and plastic particles. We explore the behavior of granular matter deep in the jammed state from local field measurement from the grain scale to the global scale. By means of a dedicated digital image correlation code and an accurate image recording method, we measure for each compression step the evolution of the particle geometries and their right Cauchy-Green strain tensor fields. We analyze the evolution of the usual macroscopic observables (stress, packing fraction, coordination, fraction of nonrattlers, etc.) along the compression process through the jamming point and far beyond. Analyzing the evolution of the local strain statistics, we evidence a crossover in the material behavior deep in the jammed state for both sorts of particles. We show that this crossover is due to a competition between material compression, dilation, and shear, so its position depends on the particle material. We argue that the strain field is a reliable observable to describe the evolution of a granular system through the jamming transition and deep in the dense packing state whatever the material behavior.

Gardner physics in amorphous solids and beyond

One of the most remarkable predictions to emerge out of the exact infinite-dimensional solution of the glass problem is the Gardner transition. Although this transition was first theoretically proposed a generation ago for certain mean-field spin glass models, its materials relevance was only realized when a systematic effort to relate glass formation and jamming was undertaken. A number of nontrivial physical signatures associated with the Gardner transition have since been considered in various areas, from models of structural glasses to constraint satisfaction problems. This perspective surveys these recent advances and discusses the novel research opportunities that arise from them.

Quasilocalized states of self stress in packing-derived networks

States of self stress (SSS) are assignments of forces on the edges of a network that satisfy mechanical equilibrium in the absence of external forces. In this work we show that a particular class of quasilocalized SSS in packing-derived networks, first introduced by D.M. Sussman, C.P. Goodrich, A.J. Liu (Soft Matter 12, 3982 (2016)), are characterized by a decay length that diverges as [Formula: see text] , where [Formula: see text] is the mean connectivity of the network, and [Formula: see text] is the Maxwell threshold in two dimensions, at odds with previous claims. Our results verify the previously proposed analogy between quasilocalized SSS and the mechanical response to a local dipolar force in random networks of relaxed Hookean springs. We show that the normalization factor that distinguishes between quasilocalized SSS and the response to a local dipole constitutes a measure of the mechanical coupling of the forced spring to the elastic network in which it is embedded. We further demonstrate that the lengthscale that characterizes quasilocalized SSS does not depend on its associated degree of mechanical coupling, but instead only on the network connectivity.

Higher-order corrections to the effective potential close to the jamming transition in the perceptron model

In view of the results achieved in a previously related work [A. Altieri, S. Franz, and G. Parisi, J. Stat. Mech. (2016) 093301]10.1088/1742-5468/2016/09/093301, regarding a Plefka-like expansion of the free energy up to the second order in the perceptron model, we improve the computation here focusing on the role of third-order corrections. The perceptron model is a simple example of constraint satisfaction problem, falling in the same universality class as hard spheres near jamming and hence allowing us to get exact results in high dimensions for more complex settings. Our method enables to define an effective potential (or Thouless-Anderson-Palmer free energy), namely a coarse-grained functional, which depends on the generalized forces and the effective gaps between particles. The analysis of the third-order corrections to the effective potential reveals that, albeit irrelevant in a mean-field framework in the thermodynamic limit, they might instead play a fundamental role in considering finite-size effects. We also study the typical behavior of generalized forces and we show that two kinds of corrections can occur. The first contribution arises since the system is analyzed at a finite distance from jamming, while the second one is due to finite-size corrections. We nevertheless show that third-order corrections in the perturbative expansion vanish in the jamming limit both for the potential and the generalized forces, in agreement with the isostaticity argument proposed by Wyart and coworkers. Finally, we analyze the relevant scaling solutions emerging close to the jamming line, which define a crossover regime connecting the control parameters of the model to an effective temperature.

Unjamming in models with analytic pairwise potentials

Canonical models for studying the unjamming scenario in systems of soft repulsive particles assume pairwise potentials with a sharp cutoff in the interaction range. The sharp cutoff renders the potential nonanalytic but makes it possible to describe many properties of the solid in terms of the coordination number z, which has an unambiguous definition in these cases. Pairwise potentials without a sharp cutoff in the interaction range have not been studied in this context, but should in fact be considered to understand the relevance of the unjamming phenomenology in systems where such a cutoff is not present. In this work we explore two systems with such interactions: an inverse power law and an exponentially decaying pairwise potential, with the control parameters being the exponent (of the inverse power law) for the former and the number density for the latter. Both systems are shown to exhibit the characteristic features of the unjamming transition, among which are the vanishing of the shear-to-bulk modulus ratio and the emergence of an excess of low-frequency vibrational modes. We establish a relation between the pressure-to-bulk modulus ratio and the distance to unjamming in each of our model systems. This allows us to predict the dependence of other key observables on the distance to unjamming. Our results provide the means for a quantitative estimation of the proximity of generic glass-forming models to the unjamming transition in the absence of a clear-cut definition of the coordination number and highlight the general irrelevance of nonaffine contributions to the bulk modulus.

Nontrivial Critical Fixed Point for Replica-Symmetry-Breaking Transitions

The transformation of the free-energy landscape from smooth to hierarchical is one of the richest features of mean-field disordered systems. A well-studied example is the de Almeida-Thouless transition for spin glasses in a magnetic field, and a similar phenomenon-the Gardner transition-has recently been predicted for structural glasses. The existence of these replica-symmetry-breaking phase transitions has, however, long been questioned below their upper critical dimension, d_{u}=6. Here, we obtain evidence for the existence of these transitions in d<d_{u} using a two-loop calculation. Because the critical fixed point is found in the strong-coupling regime, we corroborate the result by resumming the perturbative series with inputs from a three-loop calculation and an analysis of its large-order behavior. Our study offers a resolution of the long-lasting controversy surrounding phase transitions in finite-dimensional disordered systems.

Casimir effect between pinned particles in two-dimensional jammed systems

The Casimir effect arises when long-ranged fluctuations are geometrically confined between two surfaces, leading to a macroscopic force. Traditionally, these forces have been observed in quantum systems and near critical points in classical systems. Here we show the existence of Casimir-like forces between two pinned particles immersed in two-dimensional systems near the jamming transition. We observe two components to the total force: a short-ranged, depletion force and a long-ranged, repulsive Casimir-like force. The Casimir-like force dominates as the jamming transition is approached, and when the pinned particles are much larger than the ambient jammed particles. We show that this repulsive force arises due to a clustering of particles with strong contact forces around the perimeter of the pinned particles. As the separation between the pinned particles decreases, a region of high-pressure develops between them, leading to a net repulsive force.

Emergent interparticle interactions in thermal amorphous solids

Amorphous media at finite temperatures, be them liquids, colloids, or glasses, are made of interacting particles that move chaotically due to thermal energy, continuously colliding and scattering off each other. When the average configuration in these systems relaxes only at long times, one can introduce effective interactions that keep the mean positions in mechanical equilibrium. We introduce a framework to determine the effective force laws that define an effective Hessian that can be employed to discuss stability properties and the density of states of the amorphous system. We exemplify the approach with a thermal glass of hard spheres; these experience zero forces when not in contact and infinite forces when they touch. Close to jamming we recapture the effective interactions that at temperature T depend on the gap h between spheres as T/h [C. Brito and M. Wyart, Europhys. Lett. 76, 149 (2006)EULEEJ0295-507510.1209/epl/i2006-10238-x]. For hard spheres at lower densities or for systems whose binary bare interactions are longer ranged (at any density), the emergent force laws include ternary, quaternary, and generally higher-order many-body terms, leading to a temperature-dependent effective Hessian.

Understanding soft glassy materials using an energy landscape approach

Many seemingly different soft materials-such as soap foams, mayonnaise, toothpaste and living cells-display strikingly similar viscoelastic behaviour. A fundamental physical understanding of such soft glassy rheology and how it can manifest in such diverse materials, however, remains unknown. Here, by using a model soap foam consisting of compressible spherical bubbles, whose sizes slowly evolve and whose collective motion is simply dictated by energy minimization, we study the foam's dynamics as it corresponds to downhill motion on an energy landscape function spanning a high-dimensional configuration space. We find that these downhill paths, when viewed in this configuration space, are, surprisingly, fractal. The complex behaviour of our model, including power-law rheology and non-diffusive bubble motion and avalanches, stems directly from the fractal dimension and energy function of these paths. Our results suggest that ubiquitous soft glassy rheology may be a consequence of emergent fractal geometry in the energy landscapes of many complex fluids.

Finite-temperature mechanical instability in disordered lattices

Mechanical instability takes different forms in various ordered and disordered systems and little is known about how thermal fluctuations affect different classes of mechanical instabilities. We develop an analytic theory involving renormalization of rigidity and coherent potential approximation that can be used to understand finite-temperature mechanical stabilities in various disordered systems. We use this theory to study two disordered lattices: a randomly diluted triangular lattice and a randomly braced square lattice. These two lattices belong to two different universality classes as they approach mechanical instability at T=0. We show that thermal fluctuations stabilize both lattices. In particular, the triangular lattice displays a critical regime in which the shear modulus scales as G∼T(1/2), whereas the square lattice shows G∼T(2/3). We discuss generic scaling laws for finite-T mechanical instabilities and relate them to experimental systems.

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

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

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

Thermal fluctuations, mechanical response, and hyperuniformity in jammed solids

Jamming is a geometric phase transition occurring in dense particle systems in the absence of temperature. We use computer simulations to analyze the effect of thermal fluctuations on several signatures of the transition. We show that scaling laws for bulk and shear moduli only become relevant when thermal fluctuations are extremely small, and propose their relative ratio as a quantitative signature of jamming criticality. Despite the nonequilibrium nature of the transition, we find that thermally induced fluctuations and mechanical responses obey equilibrium fluctuation-dissipation relations near jamming, provided the appropriate fluctuating component of the particle displacements is analyzed. This shows that mechanical moduli can be directly measured from particle positions in mechanically unperturbed packings, and suggests that the definition of a "nonequilibrium index" is unnecessary for amorphous materials. We find that fluctuations of particle displacements are spatially correlated, and define a transverse and a longitudinal correlation length scale which both diverge as the jamming transition is approached. We analyze the frozen component of density fluctuations and find that it displays signatures of nearly hyperuniform behavior at large length scales. This demonstrates that hyperuniformity in jammed packings is unrelated to a vanishing compressibility and explains why it appears remarkably robust against temperature and density variations. Differently from jamming criticality, obstacles preventing the observation of hyperuniformity in colloidal systems do not originate from thermal fluctuations.