Unveiling Structural Bottlenecks of Dynamic Disorder in a Density-Tunable Glass Former: From Strong to Fragile Regimes

Resumen Técnico: Revelando los Cuellos de Botella Estructurales del Desorden Dinámico en un Formador de Vidrio Sintonizable por Densidad

Planteamiento del Problema
Los mecanismos microscópicos que gobiernan la ralentización dramática de los líquidos formadores de vidrio al superenfriarse siguen siendo elusivos, particularmente en lo que respecta a la distinción entre comportamientos "fuertes" y "frágiles". Mientras que los líquidos fuertes exhiben una dependencia de temperatura casi Arrhenius y los líquidos frágiles muestran un crecimiento super-Arrhenius, no está claro si estos regímenes surgen del mismo mecanismo de relajación microscópica subyacente o si la naturaleza de la relajación misma cambia fundamentalmente. Estudios previos a menudo han atribuido las diferencias en la fragilidad a variaciones en la cooperatividad o la heterogeneidad dinámica, en lugar de a mecanismos distintos, pero la evidencia microscópica directa que vincule la fragilidad con las variables estructurales específicas que controlan las tasas de relajación es limitada. Un desafío importante para resolver esto es desentrañar los mecanismos universales de las características específicas del sistema (por ejemplo, la direccionalidad de los enlaces, la composición química) inherentes a la mayoría de los formadores de vidrio.

Metodología
Para abordar estas cuestiones, los autores emplearon una mezcla binaria de repulsión suave sintonizable por densidad (partículas A-B de 50:50 que interactúan mediante potenciales de esfera blanda armónica). Este modelo permite ajustar continuamente la fragilidad desde regímenes fuertes hasta frágiles únicamente mediante la variación de la densidad, evitando cambios en la topología de interacción o la composición química.

El estudio utilizó simulaciones de dinámica molecular mediante el Simulador de Masa Paralela Atómica/Molecular de Gran Escala (LAMLAMPS) a través de densidades ($\rho$) que oscilan entre 0.65 (fuerte) y 0.82 (frágil). El análisis se centró en eventos de salto resueltos por partícula, definidos por una función de salto $h(t)$, que distingue las vibraciones de la jaula localizadas de los saltos de ruptura de la jaula. Los marcos analíticos clave incluyeron:

1. Análisis de Desorden Dinámico: Caracterización de las distribuciones de tiempo de vida de la jaula ($\psi_{cage}$) y las probabilidades de supervivencia ($C_S$) para cuantificar las desviaciones de la estadística de Poisson utilizando el parámetro de aleatoriedad $R$ y el exponente de estiramiento $\beta$.
2. Identificación de Variables Estructurales: Uso de la divergencia de Kullback–Leibler (KL) y análisis de correlación de Pearson para identificar rangos de vecinos ($k$) específicos cuyas fluctuaciones de distancia están más fuertemente acopladas a los eventos de salto.
3. Parámetro de Lentitud Estructural (SSP): Construcción de coordenadas estructurales reducidas mediante el pesaje de las fluctuaciones de los rangos de vecinos identificados. Estos SSP se utilizaron para definir subestados para el cálculo de las probabilidades de supervivencia de "fluctuación lenta", probando si descriptores estructurales específicos podrían explicar la dinámica de salto no markoviana.
4. Descriptores Adicionales: Incorporación del volumen libre de Voronoi y parámetros de orden orientacional de enlace ($q_6$) para probar su necesidad en la descripción de la dinámica lenta en el régimen frágil.
5. Correlaciones Estáticas: Cálculo de longitudes de correlación punto a conjunto (PTS, $\xi_{PTS}$) para comparar la extensión espacial del orden amorfo estático con el rango espacial de las variables lentas identificadas.

Resultados Clave

1. Evolución del Desorden Dinámico: A medida que la fragilidad aumenta (de $\rho=0.65$ a $0.82$), el sistema exhibe distribuciones de tiempo de vida de la jaula progresivamente más amplias y probabilidades de supervivencia cada vez más no exponenciales. El parámetro de aleatoriedad $R$ y la desviación de $C_S(t)$ del límite de fluctuación rápida crecen más abruptamente a temperaturas más bajas en sistemas frágiles en comparación con los sistemas fuertes.
2. Extensión Espacial de los Cuellos de Botella Estructurales:
   • Régimen Fuerte ($\rho=0.65$): La dinámica de salto está gobernada por restricciones estructurales localizadas. Las variables lentas relevantes están confinadas a un conjunto compacto de rangos de vecinos (primera capa externa y el inicio de la segunda capa). Un único SSP basado en estas distancias captura con éxito la probabilidad de supervivencia de la fluctuación lenta.
   • Régimen Intermedio ($\rho=0.72$): Al enfriar, las fluctuaciones estructurales relevantes se extienden para incluir rangos de vecinos más distantes (más profundamente en la segunda capa). La extensión espacial de las variables lentas aumenta, requiriendo un SSP extendido para describir la dinámica.
   • Régimen Frágil ($\rho=0.82$): El cuello de botella estructural involucra reordenamientos de vecinos aún más extendidos (abarcando la segunda capa). Crucialmente, los descriptores basados únicamente en la distancia se vuelven insuficientes a las temperaturas más bajas. La probabilidad de supervivencia de la fluctuación lenta solo se reproduce con precisión cuando el SSP basado en la distancia extendida se combina con el volumen libre de Voronoi, indicando que las fluctuaciones de empaquetamiento local se convierten en un componente crítico de la variable lenta que controla la tasa.
3. Comparación con el Orden Estático: La longitud de correlación punto a conjunto ($\xi_{PTS}$) crece con la fragilidad, consistente con el crecimiento del orden amorfo estático. Sin embargo, la extensión espacial de las variables lentas que controlan la tasa (inferida a partir de los rangos de vecinos involucrados en el SSP) excede consistentemente a $\xi_{PTS}$. Esto sugiere que la escala de longitud basada en el solapamiento estático no captura plenamente el rango espacial del desorden dinámico que gobierna las tasas de salto en este sistema de repulsión suave.

Significancia y Reivindicaciones
El artículo afirma establecer un marco microscópico unificado que vincula la fragilidad con la naturaleza de las variables estructurales que regulan las fluctuaciones de la tasa microscópica. La contribución primaria es la demostración de que el "cuello de botella estructural" para el desorden dinámico evoluciona sistemáticamente a través del espectro fuerte-frágil:

• En líquidos fuertes, la relajación está controlada por fluctuaciones de distancia de vecinos localizadas.
• En líquidos frágiles, la relajación está controlada por reordenamientos de vecinos extendidos y fluctuaciones sensibles al empaquetamiento (capturadas por el volumen libre).

Los autores argumentan que la fragilidad no es simplemente una amplificación de la cooperatividad local, sino que refleja un cambio fundamental en la dimensionalidad y la naturaleza de los grados de libertad lentos. Específicamente, la transición al régimen frágil implica un aumento en la dimensionalidad efectiva del espacio de variables lentas, requiriendo descriptores tanto de distancia como de volumen. Además, el estudio destaca que la correspondencia entre el orden amorfo estático (PTS) y el desorden dinámico no es universal; en esta mezcla de repulsión suave, las fluctuaciones dinámicas se extienden más allá de la escala estática, lo que sugiere que el potencial de interacción moldea significativamente la relación entre la estructura estática y la heterogeneidad dinámica.