## New Research In

### Physical Sciences

### Social Sciences

#### Featured Portals

#### Articles by Topic

### Biological Sciences

#### Featured Portals

#### Articles by Topic

- Agricultural Sciences
- Anthropology
- Applied Biological Sciences
- Biochemistry
- Biophysics and Computational Biology
- Cell Biology
- Developmental Biology
- Ecology
- Environmental Sciences
- Evolution
- Genetics
- Immunology and Inflammation
- Medical Sciences
- Microbiology
- Neuroscience
- Pharmacology
- Physiology
- Plant Biology
- Population Biology
- Psychological and Cognitive Sciences
- Sustainability Science
- Systems Biology

# Why glass elasticity affects the thermodynamics and fragility of supercooled liquids

Edited by Giorgio Parisi, University of Rome, Rome, Italy, and approved March 1, 2013 (received for review January 9, 2013)

## Abstract

Supercooled 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 that 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.

When a liquid is cooled rapidly to avoid crystallization, its viscosity increases up to the glass transition where the material becomes solid. Although this phenomenon was already used in ancient times to mold artifacts, the nature of the glass transition and the microscopic cause for the slowing down of the dynamics remain controversial. Glass-forming liquids are characterized by their fragility (1, 2); the least fragile liquids are called strong, and their characteristic time scale τ follows approximatively an Arrhenius law , where the activation energy is independent of temperature. Instead, in fragile liquids the activation energy grows as the temperature decreases, leading to a sudden slowing down of the dynamics. The fragility of liquids strongly correlates with their thermodynamic properties (3, 4): The jump in the specific heat that characterizes the glass transition is large in fragile liquids and moderate in strong ones. Various theoretical works (5⇓⇓–8), starting with Adam and Gibbs, have proposed explanations for such correlations. By contrast, a few propositions (for example, refs. 9⇓–11), have been made to understand which aspects of the microscopic structure of a liquid determine its fragility and the amplitude of the jump in the specific heat at the transition.

Observations indicate that the linear elasticity of the glass is a key factor determining fragility— a fact a priori surprising because linear elasticity is a local property of the energy landscape, whereas fragility is a nonlocal property characterizing transition between metastable states. In particular, (*i*) glasses are known to present an excess of soft elastic modes with respect to Debye vibrations at low frequencies, the so-called boson peak that appears in scattering measurements (12). The amplitude of the boson peak is strongly anticorrelated with fragility, both in network and molecular liquids: Structures presenting an abundance of soft elastic modes tend to be strong (13, 14). (*ii*) In network glasses, where particles interact via covalent bonds and via the much weaker Van der Waals interactions, the microscopic structure and the elasticity can be monitored by continuously changing the composition of compounds (15⇓–17). As the average valence *r* is increased toward some threshold , the covalent networks display a rigidity transition (18, 19) where the number of covalent bonds is just sufficient to guarantee mechanical stability. Rigidity percolation has striking effects on the thermal properties of supercooled liquids: In its vicinity, liquids are strong (20) and the jump of specific heat is small (15), whereas they become fragile with a large jump in specific heat both when the valence is increased and decreased (15, 20). It was argued (9) that fragility should decrease with valence, at least when the valence is small. There is no explanation, however, why increasing the valence affects the glass transition properties in a nonmonotonic way, and why such properties are extremal when the covalent network acquires rigidity (21).

Recently it has been shown that the presence of soft modes in various amorphous materials, including granular media (22⇓⇓–25), Lennard–Jones glasses (23, 26), colloidal suspensions (27⇓–29), and silica glass (23, 30) was controlled by the proximity of a jamming transition (31), a sort of rigidity transition that occurs, for example, when purely repulsive particles are decompressed toward vanishing pressure (24). Near the jamming transition spatial fluctuations play a limited role and simple theoretical arguments (22, 23) capture the connection between elasticity and structure. They imply that soft modes must be abundant near the transition, suggesting a link between observations (*i*) and (*ii*). However, these results apply to linear elasticity and cannot explain intrinsically nonlinear phenomena such as those governing fragility or the jump of specific heat. In this article we propose to bridge that gap by introducing a model for the structural relaxation in supercooled liquids. Our starting assumption is that particles can organize locally into distinct configurations, which are coupled at different points in space via elasticity. We study what is perhaps the simplest model realizing this idea and show numerically that it captures qualitatively the relationships between elasticity, rigidity, thermodynamics, and fragility. The thermodynamic properties of this model can be treated theoretically within a good accuracy in the temperature range we explore. Our key result is the following physical picture: When there is an abundance of soft elastic modes, elastic frustration vanishes, in the sense that a limited number of directions in phase space cost energy. Only those directions contribute to the specific heat, which is thus small. Away from the critical point, elastic frustration increases: More degrees of freedom contribute to the jump of specific heat, which increases while the boson peak is reduced.

## Model

Our main assumption is that in a supercooled liquid nearby particles can organize themselves into a few distinct configurations. Consider, for example, covalent networks sketched in Fig. 1, where we use the label to indicate the existence of a covalent bond between particles *i* and *j*. If two covalent bonds and are adjacent, there exists locally another configuration for which these bonds are broken, and where the bonds and are formed instead. These two configurations do not have the same energy in general. Moreover, going from one configuration to the other generates a local strain, which creates an elastic stress that propagates in space. In turn, this stress changes the energy difference between local configurations elsewhere in the system. This process leads to an effective interaction between local configurations at different locations.

Our contention is that even a simple description of the local configurations—in our case we will consider two-level systems, and we will make the approximation that the elastic properties do not depend on the levels—can capture several unexplained aspects of supercooled liquids, as long as the salient features of the elasticity of amorphous materials are taken into account. To incorporate in particular the presence of soft modes in the vibrational spectrum we consider random elastic networks. The elasticity of three types of networks has been studied extensively: networks of springs randomly deposited on a lattice (32), on-lattice self-organized networks (33), and off-lattice random networks with small spatial fluctuations of coordination (34⇓–36). We shall consider the last class of networks, which are known to capture correctly the scaling properties of elasticity near jamming and can be treated analytically (22, 35, 37). It is straightforward to extend our model to self-organized networks. If the network is assumed to present a rigidity window for which a subpart of the system is critical, we expect that within this window thermodynamics and fragility would behave as the critical point in the present model. In our model two kinds of springs connect the *N* nodes of the network: strong ones, of stiffness *k* and coordination *z*, and weak ones, of stiffness and coordination . These networks undergo a rigidity transition as *z* crosses , where *d* is the spatial dimension. For elastic stability is guaranteed by the presence of the weak springs. As indicated in Fig. 1, this situation is similar to covalent networks, where the weak Van der Waals interactions are required to ensure stability when the valence *r* is smaller than its critical value .

Initially when our network is built, every spring is at rest; the rest length follows , where is the initial position of the node *i*. To allow for local changes of configurations we shall consider that any strong spring can switch between two rest lengths: , where is a spin variable. There are thus two types of variables: the spin variables , which we shall denote using the ket notation , and the coordinates of the particles denoted by . The elastic energy is a function of both types of variables. The inherent structure energy associated with any configuration is defined as

where we have introduced the dimensionless Hamiltonian . We shall consider the limit of small ϵ, where the vibrational energy is simply that of harmonic oscillators. In this limit all of the relevant information is contained in the inherent structures energy, because including the vibrational energy would increase the specific heat by a constant, which does not contribute to the jump of that quantity at the glass transition. In this limit, linear elasticity implies the form

where *γ* and *β* label strong springs and is the Green function describing how a dipole of force applied on the contact *γ* changes the amplitude of the force in the contact *β*. Note that models in which some kind of defects interact elastically, leading to Hamiltonians similar in spirit to that of Eq. **2**, have been proposed to investigate the low-temperature properties of glasses (38) and supercooled liquids (10). These models, however, assume continuum elasticity, unlike our model, which can incorporate the effects of a rigidity transition and the presence of a boson peak.

is computed in *Supporting Information, A. Stiffness and Coupling Matrices* and reads

where is the identity matrix and and the dimensionless stiffness matrix are standard linear operators connecting forces and displacements in elastic networks (39). They can be formally written as

where , indicate a summation over the strong springs () or the weak springs (). is a matrix that projects any displacement field onto the contact space of strong or weak springs. The components of this linear operator are uniquely determined by the unit vectors directed along the contacts and point toward the node *i*.

Finally, note that the topology of the elastic network is frozen in our model. This addition of frozen disorder is obviously an approximation, because the topology itself should evolve as local configurations change. Our model thus misses the evolution of elasticity with temperature that presumably affects the slowing down of the dynamics (8) and gives a vibrational contribution to the jump of specific heat (8, 40). Building models that incorporate this possibility, although still tractable numerically and theoretically, remains a challenge.

## Simulation

### Network Structure.

Random networks with weak spatial fluctuations of coordination can be generated from random packings of compressed soft particles (34⇓–36). We consider packings with periodic boundary conditions. The centers of the particles correspond to the nodes of the network, of unit mass , and unstretched springs of stiffness are put between particles in contact. Then springs are removed, preferably where the local coordination is high, so as to achieve the desired coordination *z*. In a second phase, weak springs are added between the closest unconnected pairs of nodes. The relative effect of those weak springs is best characterized by , which we modulate by fixing and changing . Note that an order of magnitude estimate of *α* in covalent glasses can be obtained by comparing the behavior of the shear modulus *G* in the elastic networks (34) and in network glasses near the rigidity transition. As shown in Fig. S1, this comparison yields the estimate that discussed in *Supporting Information*, *B. Shear Modulus of the Random Elastic Networks*.

### Thermodynamics.

We introduce the rescaled temperature , where is the temperature. To equilibrate the system, we perform a one-spin-flip Monte Carlo algorithm. The energy of configurations is computed using Eq. **2**. We use five networks of nodes in two dimensions and in three dimensions, each run with 10 different initial configurations. Thus, our results are averaged on these 50 realizations. We perform Monte Carlo steps at each *T*. The time-average inherent structure energy is calculated as a function of temperature, together with the specific heat . The intensive quantity is represented in Fig. 2 for various excess coordination and . We observe that the specific heat increases under cooling, until the glass transition temperature , where rapidly vanishes, indicating that the system falls out of equilibrium.

The amplitude of just above thus corresponds to the jump of specific heat and is shown in Fig. 3. In our approach the absolute value of the specific heat will depend on the minimal number of particles needed to generate distinct local configurations, and on the number of configurations such a group can generate. If those two numbers are of order one, our model predicts that is of order one per particle or “beads,” as observed near the glass transition (6).

Our key finding is that as the coordination increases varies nonmonotonically and is minimal in the vicinity of the rigidity transition for all values of *α* investigated, as observed experimentally (15). This behavior seems to result from a sharp asymmetric transition at . For we observe that . The jump in specific heat thus vanishes as , where the system can be called “perfectly strong.” For , is very rapidly of order one. When *α* increases, this sharp transition becomes a cross-over, marked by a minimum of at some coordination larger but close to .

### Dynamics.

To characterize the dynamics we compute the correlation function , which decays to zero at long time in the liquid phase. We define the relaxation time *τ* as and the glass transition temperature as . Finite size effects on *τ* seem to be weak, as shown in *Supporting Information, C. Finite Size Effects on Fragility*. The Angell plot representing the dependence of *τ* with inverse rescaled temperature is shown in Fig. 4 *Inset*. It is found that the dynamics follow an Arrhenius behavior for and . Away from the rigidity transition, the slowing down of the dynamics is faster than Arrhenius. To quantify this effect we compute the fragility , whose variation with coordination is presented in Fig. 4. Our key finding is that for all weak interaction amplitudes *α* studied, the fragility depends nonmonotonically on coordination and is minimal near the rigidity transition, again as observed empirically (20) in covalent liquids. As was the case for the thermodynamic properties, the fragility seems to be controlled by a critical point present at and where the liquid is strong, and the dynamic is simply Arrhenius. As the coordination changes and increases, the liquid becomes more fragile. The rapid change of fragility near the rigidity transition is smoothed over when the amplitude of the weak interaction *α* is increased.

### Correlating Boson Peak and Fragility.

The presence of soft elastic modes in glasses can be analyzed by considering the maximum of (12), where is the vibrational density of states and the Debye model for this quantity. The maximum of , defines the boson peak frequency (12) and its normalized amplitude (13). The inverse of was shown to strongly correlate with fragility (13, 14) both in molecular liquids and covalent networks.

To test whether our model can capture this behavior we compute the density of states via a direct diagonalization of the stiffness matrix (Eq. **4**). To compute the Debye density of states is estimated as , where *G* is the shear modulus (bulk and shear moduli scale identically in this model; e.g., ref. 37). Then we extract the maximum . The dependence of is represented in Fig. 5 and shows a minimum near the rigidity transition, and even a cusp in the limit . This behavior can be explained in terms of previous theoretical results on the density of states near the rigidity transition that support that when . When and , and (22), whereas (23), leading to . For , (34). However, the boson peak is governed by the fraction of floppy modes, which gain a finite frequency (35), and thus we expect and .

Fig. 5 shows that and the liquid fragility *m* are correlated in our model, thus capturing observations in molecular liquids. The model also predicts that and the jump of specific heat are correlated. Note that the correlation between fragility and is not perfect and that two branches, for glasses with low and with high coordinations, are clearly distinguishable. In general we expect physical properties to depend on the full structure of the density of states, as will be made clear for the thermodynamics of our model below. The variable , which is a single number, cannot fully capture this relationship. In our framework it is a useful quantity, however, because it characterizes well the proximity of the jamming transition.

## Theory

### Thermodynamics in the Absence of Weak Interactions ().

In the absence of weak springs the thermodynamics is nontrivial if ; otherwise, the inherent structure energies are all zero. Then Eq. **4** implies , and Eq. **3** leads to . Inspection of this expression indicates that is a projector on the kernel of , which is generically of dimension . This kernel corresponds to all of the sets of contact forces that balance forces on each node (23). We denote by an orthonormal basis of this space. We may then rewrite and Eq. **2** as

Eq. **5** is a key result, because it implies that near the rigidity transition the number of directions of phase space that cost energy vanishes. Only those directions can contribute to the specific heat, which must thus vanish linearly in as the rigidity transition is approached from above.

Eq. **5** also makes a connection between strong liquids in our framework and well-known problems in statistical mechanics. In particular, Eq. **5** is similar to that describing Hopfield nets (41) used to store memories consisting of the spin states . The key difference is the sign: In Hopfield nets memories correspond to metastables states, whereas in our model the vectors correspond to maxima of the energy. A particularly interesting case is , the closest point to the jamming transition that is nontrivial. In this situation the sum in Eq. **5** contains only one term: . This Hamiltonian corresponds to the NP-complete partitioning problem (42), where given a list of numbers (the ) one must partition this list into two groups whose sums are as identical as possible. Thermodynamically this problem is known (43) to map into the random energy model (44), where energy levels are randomly distributed.

It is in general very difficult to compute the thermodynamic functions of the problem defined by Eq. **5** because the vectors present spatial correlations, as must be the case because the amplitude of the interaction kernel must decay with distance. However, this effect is expected to be mild near the rigidity transition. Indeed, there exists a diverging length scale at the transition (ref. 35 includes a recent discussion), below which is dominated by fluctuations and decays mildly with distance. Beyond this length scale presents a dipolar structure, as in a standard continuous elastic medium. We shall thus assume that are random unitary vectors, an approximation of mean-field character expected to be good near the rigidity transition.

Within this approximation, the thermodynamic properties can be derived for any spectrum of (45). If the orthogonality of the vectors is preserved, the Hamiltonian of Eq. **5** corresponds to the random orthogonal model (ROM) whose thermodynamic properties have been derived (45) as well as some aspects of the dynamics (46). A comparison of the specific heat of our model and the ROM predictions of ref. 45 is shown in Fig. 3, and they are found to be very similar. For sake of simplicity, in what follows we shall also relax the orthogonal condition on the vectors . This approximation allows for a straightforward analytical treatment in the general case , and is also very accurate near the rigidity transition because the number of vectors is significantly smaller than the dimension of the space they live in, making random vectors effectively orthogonal. Under these assumptions we recover the random Hopfield model with negative temperature.

In the parameter range of interest, the Hopfield free energy [here represents the disorder average on the ] is approximated very precisely by the annealed free energy (this is obviously true for the number partitioning problem that maps to the random energy model), which can be easily calculated. Indeed, in our approximations the quantities are independent Gaussian random variables of variance one, and

Performing the Gaussian integrals we find

The Kautzman temperature defined as is found to follow . Eq. **7** evaluated at is tested against the numerics in Fig. 3 and performs remarkably well for the range of coordination probed.

### General Case ().

To solve our model analytically in the presence of weak interactions, we make the additional approximation that the associated coordination , while keeping constant. In this limit weak springs lead to an effective interaction between each node and the center of mass of the system, which is motionless. Thus, the restoring force stemming from weak interactions follows , leading to a simple expression in the stiffness matrix Eq. **4** for the weak spring contribution . It is useful to perform the eigenvalue decompositionwhere is the vibrational mode of frequency *ω* in the elastic network without weak interactions. We introduce the orthonormal eigenvectors in contact space defined for . For these vectors form a complete basis of that space, of dimension . When , however, this set is of dimension , and it must be completed by the kernel of , that is, the set of the previously introduced. Using this decomposition in Eqs. **3** and **4** we findwhere the first term exists only for . Using the mean field approximation that the set of and are random Gaussian vectors, the annealed free energy is readily computed, as shown in *Supporting Information, E. Canonical Ensemble Approach with Weak Springs*. We find in particular for the specific heatwhere is the unitary step function. To compare this prediction with our numerics without fitting parameters, we compute numerically the vibrational frequencies for each value of the coordination. Our results are again in excellent agreement with our observations, as seen in Figs. 2 and 3.

To obtain the asymptotic behavior near jamming, we replace the summation over frequencies in Eq. **10** by an integral . The associated density of vibrational modes in such networks has been computed theoretically (22, 35, 37). These results allow us to compute the scaling behavior of thermodynamic properties near the rigidity transition (*Supporting Information, F. Analytical Results of Continuous Density of States Limit*). We find that the specific heat increases monotonically with decreasing temperature. Its value at the Kautzman temperature thus yields an upper bound on the jump of specific heat. In the limit , we find that a sudden discontinuity of the jump of specific heat occurs at the rigidity transition:

Eq. **11** states that adding weak interactions is not a singular perturbation for , and we recover Eq. **7**. However, for the energy of inherent structures is zero in the absence of weak springs, which thus have a singular effect. The relevant scale of temperature is then a function of *α*. In particular, we find that the Kautzman temperature is sufficiently low that all of the terms in the second sum of Eq. **10** contribute significantly to the specific heat, which is therefore large as Eq. **12** implies. Thus, as the coordination decreases below the rigidity transition, one goes discontinuously from a regime where at the relevant temperature scale the energy landscape consists of a vanishing number of costly directions in phase space, whose cost is governed by the strong interaction *k*, to a regime where the weak interaction *α* is the relevant one, and where at the relevant temperature scale all directions in phase space contribute to the specific heat.

Note that although the sharp change of thermodynamic behavior that occurs at the rigidity transition is important conceptually, empirically a smooth cross-over will always be observed. This is the case because (*i*) *α* is small but finite. As *α* increases this sharp discontinuity is replaced by a cross-over at a coordination (*Supporting Information, F. Analytical Results of Continuous Density of States Limit*) where is minimal, as indicated in Fig. 3 *Inset*. (*ii*) The Kautzman temperature range is not accessible dynamically, that is, near the rigidity transition. Comparing Fig. 3 with its inset, our theory predicts that the minimum of is closer to and more pronounced than at .

## Discussion

Previous work (31) has shown that well-coordinated glasses must have a small boson peak, which increases as the coordination (or valence for network glasses) is decreased toward the jamming (or rigidity) transition. Here we have argued that as this process occurs elastic frustration vanishes: Thanks to the abundance of soft modes, any configuration (conceived here as a set of local arrangements of the particles) can relax more and more of its energy as jamming is approached from above. As a result, the effective number of degrees of freedom that cost energy and contribute to the jump of specific heat at the glass transition vanishes. As the coordination is decreased further below the rigidity transition, the scale of energy becomes governed by the weak interactions (such as Van der Waals) responsible for the finite elasticity of the glass. At that scale, all direction in phase space has a significant cost and the specific heat increases. This view potentially explains why linear elasticity strongly correlates to key aspects of the energy landscape in network and molecular glasses (13⇓⇓–16). This connection we propose between structure and dynamics can also be tested numerically. Our results are consistent with the observation that soft particles become more fragile when compressed away from the random close packing (47). A quantitative treatment of the vicinity of random close packing would require including the dependence of elasticity on temperature neglected here. Another interesting parameter to manipulate is the amplitude of weak interactions, which can be increased by adding long-range forces to the interaction potential (23, 26). According to our analysis, doing so should increase fragility, in agreement with existing observations (48).

The model of the glass transition we introduced turns out to be a spin glass model, with the specificity that (*i*) the interaction is dipolar in the far field and that (*ii*) the sign of the interaction is approximatively random below some length scale that diverges near jamming, where the coupling matrix has a vanishingly small rank. Applying spin glass models to structural glasses have a long history. In particular the random first order transition theory (6) is based on mean-field spin glass models that display a thermodynamic transition at some where the entropy vanishes. A phenomenological description of relaxation in liquids near based on the nucleation of random configurations leads to a diverging time scale and length scale ξ at (6, 7). One limitation of this approach is that no finite dimensional spin models with two-body interactions have been shown to follow this scenario so far (49), and it would thus be important to know whether our model does display a critical point at finite temperature. Our model will also allow one to investigate the generally neglected role of the action at a distance allowed by elasticity, characterized by a scale . In supercooled liquids heterogeneities of elasticity (that correlate to irreversible rearrangements) can be rather extended (50), suggesting that is large. This length scale may thus play an important role in a description of relaxation in liquids, and in deciphering the relationship between elastic and dynamical heterogeneities.

## Acknowledgments

We thank A. Grosberg, P. Hohenberg, E. Lerner, D. Levine, D. Pine, E. Vanden-Eijnden, M. Vucelja., and A. Lefèvre for discussions and E. Lerner for discussions leading to Eq. **S1**. This work was supported primarily by National Science Foundation (NSF) Chemical, Bioengineering, Environmental and Transport Systems Division Grant 1236378 and partially by the Sloan Fellowship, NSF Division of Materials Research Grant 1105387, and Petroleum Research Fund Grant 52031-DNI9.

## Footnotes

- ↵
^{1}To whom correspondence should be addressed. E-mail: mw135{at}nyu.edu.

Author contributions: L.Y., G.D., and M.W. designed research, performed research, analyzed data, and wrote the paper.

The authors declare no conflict of interest.

This article is a PNAS Direct Submission.

This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1300534110/-/DCSupplemental.

## References

- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- Dyre JC

- ↵
- ↵
- ↵
- ↵
- Phillips WA

- Anderson AC

- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- Phillips JC

- ↵
- ↵
- Böhmer R,
- Angell CA

- ↵
- ↵
- ↵
- Wyart M

- ↵
- ↵
- ↵
- ↵
- ↵
- Ghosh A,
- Chikkadi VK,
- Schall P,
- Kurchan J,
- Bonn D

- ↵
- ↵
- ↵
- Berthier L,
- Biroli G,
- Bouchaud JP,
- van Saarloos W

- Liu AJ,
- Nagel SR,
- van Saarloos W,
- Wyart M

- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- Hopfield JJ

- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- Berthier L,
- Tarjus G

- ↵Cammarota C, Biroli G, Tarzia M, Tarjus G (2012) On the fragility of the mean-field scenario of structural glasses for finite-dimensional disordered spin models. arXiv:1210.2941.
- ↵

## Citation Manager Formats

## Sign up for Article Alerts

## Article Classifications

- Physical Sciences
- Physics