Local micromechanics in a mean-field model of glasses reveal key properties of its non-equilibrium RSB phase

Technical Summary: Local Micromechanics in a Mean-Field Model of Glasses

Problem and Motivation
Glasses are out-of-equilibrium solids characterized by disorder and a lack of long-range order. While continuum-scale properties (e.g., elasticity) are well-understood, the microscopic origins of non-equilibrium phenomena—such as non-Debye vibrational spectra, yielding transitions, and mechanical memory—remain challenging. Recent advances in finite-dimensional computer glasses have identified "quasi-localized" low-frequency modes and a universal $\sim\omega^4$ tail in the vibrational density of states (VDoS), $D_G(\omega)$.

A recently formulated disordered Hamiltonian mean-field model [47, 48] has shown promise in unifying these glassy properties. This model features an equilibrium, zero-temperature Replica-Symmetry-Breaking (RSB) transition in certain parameter ranges. In the Replica-Symmetric (RS) phase approaching the RSB transition, the solution is described by an effective, self-consistent random potential for uncoupled degrees of freedom. However, the behavior of this model deep within the non-equilibrium RSB phase (following an instantaneous quench) and its connection to the effective potential description required further investigation. Specifically, there was a need to define "micromechanical" response functions in a fully connected model lacking spatial structure and to establish exact relations between local observables and global susceptibilities.

Methodology
The authors analyze a mean-field glass model defined by the Hamiltonian:
$$H = \frac{1}{2}\sum_{i=1}^N k_i x_i^2 + \frac{1}{4!}\sum_{i=1}^N x_i^4 + \sum_{i<j} J_{ij}x_i x_j - f(t)\sum_{i=1}^N x_i$$
where $k_i$ are initial stiffnesses drawn from a distribution (specifically $p(k) \propto k^{-1}$ on $[0,1]$ for the studied case), and $J_{ij}$ are Gaussian random couplings. The system is studied in the non-equilibrium RSB phase by minimizing $H$ (instantaneous quench) without thermalization.

The core methodological innovation is the definition of a force monopole response in this mean-field context. Analogous to force dipoles in finite-dimensional glasses (which screen long-wavelength phonons), the authors apply a local force monopole $f_i$ (a unit force on the $i$-th oscillator) and measure the linear response $z_i = M^{-1} \cdot f_i$, where $M$ is the Hessian matrix. They define a local monopole stiffness $\kappa_i$:
$$\kappa_i = \frac{f_i \cdot M^{-1} \cdot f_i}{f_i \cdot M^{-2} \cdot f_i} = \frac{(M^{-1})_{ii}}{(M^{-2})_{ii}}$$
Using Random Matrix Theory (RMT) techniques, the authors derive an exact mean-field identity relating these local stiffnesses to global linear ($\chi$) and spin-glass ($\chi_{SG}$) susceptibilities. They then utilize numerical simulations to validate these relations, reconstruct the full VDoS from monopole statistics, and identify a "stabilized sector" of oscillators based on a marginal stability criterion.

Key Contributions and Results

1. Exact Mean-Field Identity: The authors derive a rigorous relation connecting the local monopole stiffness $\kappa_i$ to global susceptibilities:
   $$\kappa_i = \frac{(k_i - J^2\chi) + \frac{1}{2}x_i^2}{1 + J^2\chi_{SG}}$$
   This identity demonstrates that entirely local properties are determined by global quantities, a manifestation of the model's fully-connected nature. Crucially, for $\kappa_i$ to be finite, the spin-glass susceptibility $\chi_{SG} = \int \omega^{-4} D_G(\omega) d\omega$ must be finite. This necessitates a pseudogapped VDoS with a tail exponent $>3$, consistent with the model's observed $\sim\omega^4$ tail.

2. Micromechanical Localization and Pseudogap: Numerical analysis reveals that the monopole response $z_i$ is strongly localized (participation number $\langle N_e \rangle \approx 4$, independent of system size $N$), mirroring quasi-localized modes in finite-dimensional glasses. Furthermore, the distribution of monopole stiffnesses $p(\kappa)$ is pseudogapped (vanishing as $\kappa \to 0$), indicating marginal stability.

3. Reconstruction of Vibrational Spectra: The authors demonstrate that the full nonphononic VDoS, $D_G(\omega)$, can be reconstructed solely from the statistics of the monopole stiffnesses $\{\kappa_i\}$. By solving a self-consistency equation for the resolvent using the derived stiffness statistics, they recover $D_G(\omega)$ with excellent quantitative agreement against direct Hessian diagonalization.

4. Characteristic Scales and the Boson Peak: The study identifies a single characteristic frequency scale controlling the nonphononic VDoS.

   • The authors define a "stabilized sector" of oscillators where the initial stiffness $k_i < J^2\chi$. These oscillators rely on the anharmonic displacement term ($\frac{1}{2}x_i^2$) to maintain positive stiffness.
   • The mean stiffness of this stabilized sector, $\langle \kappa \rangle_{st}$, scales as the square of the intrinsic boson-peak frequency: $\langle \kappa \rangle_{st} \sim \omega_p^2$.
   • Additionally, the characteristic frequency scale of the effective random potential, derived from force fluctuations $\Delta f$, scales as $\omega_p \sim (\Delta f)^{1/3}$.

5. Mesoscopic vs. Global Elasticity: The ratio of the mesoscopic stiffness (characterized by $\langle \kappa \rangle_{st}$) to the global shear modulus $G$ varies significantly with the disorder level $J$. Specifically, $\langle \kappa \rangle_{st}$ exhibits a broader range of variation than $G$, a behavior consistent with observations in finite-dimensional computer glasses.

Significance and Claims
The paper claims that the statistical micromechanics of the mean-field model, specifically the force monopole response, provide a concrete realization of the collective degrees of freedom described by the effective random potential in the equilibrium RS phase. The findings suggest that:

• The mathematical structure of the effective potential (derived for the equilibrium RS phase) remains valid deep within the non-equilibrium RSB phase.
• The "stabilized sector" of oscillators (those with $k_i < J^2\chi$) constitutes the physical manifestation of the soft glassy modes responsible for the $\omega^4$ tail and the boson peak.
• The model offers a unified framework bridging mean-field theory and finite-dimensional glass physics, where local micromechanical observables (monopole stiffness) directly encode global spectral properties and susceptibility limits.

The authors conclude that these results reinforce the theoretical suggestions that the model's energy landscape is inherently glassy and that the effective potential picture provides a unifying description for its non-equilibrium properties. They note that future work should explore sparse connectivity, anharmonic properties, and strongly driven dynamics.