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

Resumen Técnico: Micromecánica Local en un Modelo de Campo Medio de Vidrios

Problema y Motivación
Los vidrios son sólidos fuera del equilibrio caracterizados por el desorden y la falta de orden de largo alcance. Mientras que las propiedades a escala de continuo (por ejemplo, la elasticidad) están bien comprendidas, los orígenes microscópicos de los fenómenos fuera del equilibrio —como los espectros vibracionales no Debye, las transiciones de fluencia (yielding) y la memoria mecánica— siguen siendo un desafío. Avances recientes en vidrios computacionales de dimensión finita han identificado modos de baja frecuencia "cuasi-localizados" y una cola universal de $\sim\omega^4$ en la densidad de estados vibracionales (VDoS), $D_G(\omega)$.

Un modelo de campo medio de Hamiltoniano desordenado recientemente formulado [47, 48] ha mostrado ser prometedor para unificar estas propiedades vítreas. Este modelo presenta una transición de Ruptura de Simetría de Réplica (RSB) de equilibrio a temperatura cero en ciertos rangos de parámetros. En la fase de Simetría de Réplica (RS) que se aproxima a la transición RSB, la solución se describe mediante un potencial aleatorio efectivo y auto-consistente para grados de libertad desacoplados. Sin embargo, el comportamiento de este modelo profundamente dentro de la fase RSB fuera del equilibrio (tras un enfriamiento instantáneo o quench) y su conexión con la descripción del potencial efectivo requerían mayor investigación. Específicamente, existía la necesidad de definir funciones de respuesta "micromecánicas" en un modelo de conectividad total que carece de estructura espacial, y de establecer relaciones exactas entre observables locales y susceptibilidades globales.

Metodología
Los autores analizan un modelo de vidrio de campo medio definido por el Hamiltoniano:
$$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$$
donde $k_i$ son rigideces iniciales extraídas de una distribución (específicamente $p(k) \propto k^{-1}$ en $[0,1]$ para el caso estudiado), y $J_{ij}$ son acoplamientos aleatorios gaussianos. El sistema se estudia en la fase RSB fuera del equilibrio minimizando $H$ (un quench instantáneo) sin termalización.

La innovación metodológica central es la definición de una respuesta de monopolo de fuerza en este contexto de campo medio. Análogamente a los dipolos de fuerza en vidrios de dimensión finita (que blindan los fonones de longitud de onda larga), los autores aplican un monopolo de fuerza local $f_i$ (una fuerza unitaria sobre el oscilador $i$) y miden la respuesta lineal $z_i = M^{-1} \cdot f_i$, donde $M$ es la matriz Hessiana. Definen una rigidez de monopolo local $\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}}$$
Utilizando técnicas de Teoría de Matrices Aleatorias (RMT), los autores derivan una identidad exacta de campo medio que relaciona estas rigideces locales con las susceptibilidades lineales ($\chi$) y de vidrio de espín ($\chi_{SG}$) globales. Posteriormente, utilizan simulaciones numéricas para validar estas relaciones, reconstruir la VDoS completa a partir de las estadísticas de los monopolos e identificar un "sector estabilizado" de osciladores basado en un criterio de estabilidad marginal.

Contribuciones Clave y Resultados

1. Identidad Exacta de Campo Medio: Los autores derivan una relación rigurosa que conecta la rigidez de monopolo local $\kappa_i$ con las susceptibilidades globales:
   $$\kappa_i = \frac{(k_i - J^2\chi) + \frac{1}{2}x_i^2}{1 + J^2\chi_{SG}}$$
   Esta identidad demuestra que las propiedades enteramente locales están determinadas por cantidades globales, una manifestación de la naturaleza de conectividad total del modelo. Crucialmente, para que $\kappa_i$ sea finito, la susceptibilidad de vidrio de espín $\chi_{SG} = \int \omega^{-4} D_G(\omega) d\omega$ debe ser finita. Esto requiere una VDoS con un pseudogap con un exponente de cola $>3$, consistente con la cola $\sim\omega^4$ observada en el modelo.

2. Localización Micromecánica y Pseudogap: El análisis numérico revela que la respuesta de monopolo $z_i$ está fuertemente localizada (número de participación $\langle N_e \rangle \approx 4$, independiente del tamaño del sistema $N$), reflejando los modos cuasi-localizados en vidrios de dimensión finita. Además, la distribución de las rigideces de monopolo $p(\kappa)$ presenta un pseudogap (desvaneciéndose cuando $\kappa \to 0$), lo que indica estabilidad marginal.

3. Reconstrucción de Espectros Vibracionales: Los autores demuestran que la VDoS completa no fonónica, $D_G(\omega)$, puede reconstruirse únicamente a partir de las estadísticas de las rigideces de monopolo $\{\kappa_i\}$. Al resolver una ecuación de auto-consistencia para el resolvente utilizando las estadísticas de rigidez derivadas, recuperan $D_G(\omega)$ con una excelente concordancia cuantitativa frente a la diagonalización directa de la Hessiana.

4. Escalas Características y el Pico de Bosón: El estudio identifica una única escala de frecuencia característica que controla la VDoS no fonónica.

   • Los autores definen un "sector estabilizado" de osciladores donde la rigidez inicial $k_i < J^2\chi$. Estos osciladores dependen del término de desplazamiento anharmónico ($\frac{1}{2}x_i^2$) para mantener una rigidez positiva.
   • La rigidez media de este sector estabilizado, $\langle \kappa \rangle_{st}$, escala como el cuadrado de la frecuencia intrínseca del pico de bosón: $\langle \kappa \rangle_{st} \sim \omega_p^2$.
   • Adicionalmente, la escala de frecuencia característica del potencial aleatorio efectivo, derivada de las fluctuaciones de fuerza $\Delta f$, escala como $\omega_p \sim (\Delta f)^{1/3}$.

5. Elasticidad Mesoscópica vs. Global: La relación entre la rigidez mesoscópica (caracterizada por $\langle \kappa \rangle_{st}$) y el módulo de corte global $G$ varía significativamente con el nivel de desorden $J$. Específicamente, $\langle \kappa \rangle_{st}$ exhibe una variación más amplia que $G$, un comportamiento consistente con las observaciones en vidrios computacionales de dimensión finita.

Significancia y Reivindicaciones
El artículo sostiene que la micromecánica estadística del modelo de campo medio, específicamente la respuesta de monopolo de fuerza, proporciona una realización concreta de los grados de libertad colectivos descritos por el potencial aleatorio efectivo en la fase RS de equilibrio. Los hallazgos sugieren que:

• La estructura matemática del potencial efectivo (derivada para la fase RS de equilibrio) permanece válida profundamente dentro de la fase RSB fuera del equilibrio.
• El "sector estabilizado" de osciladores (aquellos con $k_i < J^2\chi$) constituye la manifestación física de los modos vítreos blandos responsables de la cola $\omega^4$ y del pico de bosón.
• El modelo ofrece un marco unificado que une la teoría de campo medio y la física de vidrios de dimensión finita, donde los observables micromecánicos locales (rigidez de monopolo) codifican directamente las propiedades espectrales globales y los límites de susceptibilidad.

Los autores concluyen que estos resultados refuerzan las sugerencias teóricas de que el paisaje de energía del modelo es inherentemente vítreo y que la descripción del potencial efectivo proporciona una descripción unificada para sus propiedades fuera del equilibrio. Señalan que trabajos futuros deberían explorar la conectividad dispersa (sparse), las propiedades anharmónicas y la dinámica fuertemente impulsada.