Replica theory of the rigidity of structural glasses

Abstract: We present a first principle scheme to compute the rigidity, i. e. the shear-modulus of structural glasses at finite temperatures using the cloned liquid theory, which combines the replica theory and the liquid theory. With the aid of the replica method which enables disentanglement of thermal fluctuations in liquids into intra-state and inter-state fluctuations, we extract the rigidity of metastable amorphous solid states in the supercooled liquid and glass phases. The result can be understood intuitively without replicas. As a test case, we apply the scheme to the supercooled and glassy state of a binary mixture of soft-spheres. The result compares well with the shear-modulus obtained by a previous molecular dynamic simulation. The rigidity of metastable states is significantly reduced with respect to the instantaneous rigidity, namely the Born term, due to non-affine responses caused by displacements of particles inside cages at all temperatures down to T=0. It becomes nearly independent of temperature below the Kauzmann temperature T_K. At higher temperatures in the supercooled liquid state, the non-affine correction to the rigidity becomes stronger suggesting melting of the metastable solid state. Inter-state part of the static response implies jerky, intermittent stress-strain curves with static analogue of yielding at mesoscopic scales.

• The paper demonstrates that a cloned liquid replica approach successfully computes the shear modulus by separating intra-state and inter-state contributions.
• It uses a systematic cage expansion technique to account for non-affine corrections and captures temperature-dependent rigidity near the Kauzmann transition.
• The methodology links metastable state dynamics with measurable viscoelastic responses, validated against molecular dynamics benchmarks.

Replica Theory Approach to the Rigidity of Structural Glasses

Introduction and Motivation

The computation of mechanical rigidity—specifically, the shear modulus—of structural glasses from microscopic principles remains a central problem in condensed matter theory, with wide-ranging implications for both theoretical understanding and material design. The lack of long-range order in amorphous solids precludes direct application of standard elasticity theory as used in crystals. Moreover, the emergent rigidity in glasses arises from complex structural correlations and metastable states generated via supercooling, undergoing visco-elastic responses with pronounced $\alpha$ and $\beta$ regimes. The "Replica theory of the rigidity of structural glasses" (1203.2341) systematically addresses the computation of shear modulus in glasses using the cloned liquid formalism, combining replica and liquid theories to disentangle intra-state and inter-state thermal fluctuations. This formalism provides a path to statically define and compute the effective rigidity associated with metastable amorphous states at finite temperature.

Theoretical Framework and Methodology

The central theoretical advance is the use of the cloned liquid (replicated liquid) approach, grounded in the random first order transition (RFOT) framework. The method constructs $m$ non-interacting replicas ("clones") of the system, biased to fall into the same metastable free energy basin. This is implemented as an effective "molecular" liquid, with each molecule corresponding to the coordinates of a particle across all replicas, caged near its center of mass with variance (cage size) $A$. $A$ functions as the glass order parameter—analogous to the Edwards-Anderson parameter in spin glasses—signaling ergodicity breaking and the emergence of amorphous order.

A key theoretical component is the decomposition of static response functions into intra-state (metabasin) and inter-state (inter-metabasin) contributions via the replica structure:

• The intra-state part (effective rigidity $\hat{\mu}$) captures the response within a typical metastable state, corresponding to $\beta$-relaxation.
• The inter-state part embodies transitions between valleys (configurational states), which restore macroscopic translational symmetry in the $N\to\infty$ limit, leading to an eventual vanishing of bulk rigidity in static equilibrium.

The physical shear modulus is then extracted through a fluctuation formula based on the second derivative of the free energy with respect to infinitesimal shear, with explicit microscopic implementation in the presence of local cages. A systematic cage expansion (in powers of $A$) is used to calculate corrections to the Born (affine) term, culminating in robust expressions for the shear modulus in terms of cage size and (renormalized) liquid structure factors.

Analytical Formulation and Hierarchy of Rigidities

The static shear modulus in the glass phase is obtained as

$$\mu = \langle b \rangle - N\beta \left( \langle \sigma^2 \rangle - \langle \sigma \rangle^2 \right)$$

where $b$ is the "Born" term (affine response) and the second term is the non-affine correction due to thermal and configurational stress fluctuations. Under the cloned liquid formalism, the thermal average splits via replica symmetry breaking structure into intra-state and inter-state contributions. The intra-state modulus $\hat{\mu}$ is:

$$\hat{\mu} = b - N\beta \left(\langle \sigma^2 \rangle_\text{MS} - \langle \sigma \rangle_\text{MS}^2 \right)$$

(Here, $\langle \cdot \rangle_\text{MS}$ denotes averaging within a metabasin).

A hierarchy of rigidities emerges:

• The instantaneous Born term $b$,
• The rigidity at the level of inherent structures (energy minima, $\mu_\text{IS}$),
• The rigidity of metastable states (metabasins, $\hat{\mu}$), i.e., the effective $\beta$-regime rigidity,
• The vanishing $N\to\infty$ static modulus after including all inter-state fluctuations.

It is shown analytically and in application that $\hat{\mu} < \mu_\text{IS} < b$, with the difference arising from inelastic corrections due to cage fluctuations and transitions between ISs inside a metabasin—non-affine processes capturing the collective and localized plasticity unique to glassy states.

Key Results and Numerical Benchmarks

Application to a binary mixture of soft spheres allows for concrete evaluation:

• The full cage expansion yields explicit analytic expressions for $\hat{\mu}$ as a function of $A$, structural correlation functions, and interaction derivatives, with non-affine corrections evaluated via standard techniques (hypernetted-chain).
• Comparison to reference molecular dynamics results by Barrat et al. reveals quantitative agreement across the glass and supercooled regimes: the computed $\hat{\mu}$ remains finite and nearly temperature-independent below the computed Kauzmann temperature $T_K$, while it sharply drops at higher $T$ at and above $T_K$, in concordance with the simulated shear modulus. The position of the vanishing of $\hat{\mu}$ closely tracks the mode-coupling crossover.
• The Born term contrasts sharply: it increases with $T$ (a direct result of stronger particle collisions), demonstrating that the physical shear modulus is dominated at all $T$ by significant non-affine reductions arising from local cage dynamics and collective rearrangements.

Notably, the clear temperature dependence and eventual collapse of $\hat{\mu}$ above $T_K$ strongly supports the theoretical RFOT scenario: the rigidity is governed by a discontinuous transition of the order parameter (cage size), with mean-field analysis predicting a square-root singularity as $T$ approaches the dynamical transition from below. However, finite-dimensional corrections smooth the transition, and full vanishing is only expected via decorrelation of mesoscopic metastable regions.

Theoretical Consequences: Rigidity Crisis and Metastability

The work provides robust mathematical support for the view that rigidity in glasses is not set by a putative continuous yielding (as in crystals) but by a critical-like scenario governed by the order parameter, exhibiting a discontinuous drop at the (mean-field) glass transition. The decomposition via replicas clarifies the essential role of plastic (non-affine) reorganizations even deep in the glass phase, with an explicit link between shear modulus and the "cage" order parameter.

Shear-induced level crossings between metastable states induce intermittent stress-strain behavior at mesoscopic scales. This static analog of yielding arises even in the absence of macroscopic flow and is theoretically universal in the presence of configurational multiplicity, with predicted scaling for event spacings.

Integration over cell-level shear fluctuations in coarse-grained functional theories predicts further coupling and potential first-order-like melting of metastability as non-affinity grows with temperature, refining global elasticity beyond mean-field Ginzburg-Landau approaches.

Practical and Conceptual Implications

• The formalism delivers a first-principles framework for computing mechanical response in amorphous solids, with direct connections to observable rheological features (elastic lags, onset of yielding, plastic events).
• The approach is sufficiently general for extension to other types of disordered solids (e.g., gels, granular packings) with appropriate generalizations of the cloned liquid machinery, and provides a foundation for more accurate treatments of the elastic and non-elastic constants of real glasses.
• The identification of intra-state susceptibilities as robust quasi-equilibrium measures opens a principled route for connecting theoretical quantities to experimental observables in out-of-equilibrium protocols.

Conclusion

This work rigorously establishes the replica-based cloned liquid as a powerful microscopic backbone for the calculation of rigidity in structural glasses (1203.2341). The analysis clarifies the role of metastable state structure, the necessity of non-affine corrections rooted in local and collective plasticity, and delineates the theoretical predictions for viscoelastic response and melting. The ability to quantitatively match two-body theory predictions with simulation benchmarks underscores the soundness of the framework. The implications extend to the design of coarse-grained field theories for amorphous matter, inform on the universality of rigidity transitions, and set the stage for refined treatments of finite-dimensional effects and dynamic heterogeneity in glassy materials.