Coarse-graining amorphous plasticity: impact of rejuvenation and disorder

Resumen Técnico: Granularidad Gruesa de la Plasticidad Amorfa: Impacto del Rejuvenecimiento y el Desorden

Planteamiento del Problema
El comportamiento mecánico de los sólidos amorfos, caracterizado por fenómenos como el envejecimiento, el rejuvenecimiento, la formación de bandas de cizalla y el efecto Bauschinger, surge de reajustes estructurales localizados (Transformaciones de Cizalla) que interactúan a través de un medio elástico. Si bien las simulaciones atomísticas capturan estas características, son computacionalmente costosas. Por el contrario, los Modelos Elasto-Plásticos (EPM) mesoscópicos, que son modelos de coarse-graining, ofrecen un marco computacionalmente eficiente, pero a menudo dependen de suposiciones simplificadas que pueden fallar al intentar reproducir cuantitativamente los resultados atomísticos. El problema central abordado es el desafío del "coarse-graining" (granularidad gruesa): determinar cómo transferir parámetros de escalas atomísticas a modelos de red mesoscópicos para reproducir la física de la plasticidad amorfa, centrándose específicamente en el papel del paisaje de energía potencial (PEL) desordenado subyacente, la evolución del desorden bajo carga y el rejuvencimiento del estado vítreo.

Metodología
El estudio emplea un enfoque comparativo entre simulaciones atomísticas de Cuasi-Estática Atérmica (AQS) y un modelo EPM escalar mesoscópico simplificado.

1. Simulaciones Atomísticas:

   • Se prepararon tres muestras distintas de vidrio de Lennard-Jones binario 2D con historias térmicas contrastantes:
     • MCT: Enfriada instantáneamente desde una temperatura cercana a la Transición de Acoplamiento de Modo ($T_p \approx 0.38$).
     • ESL: Enfriada instantáneamente desde una temperatura más alta ($T_p = 0.351$), representando un vidrio mal recocido (blando).
     • GQ: Enfriada gradualmente desde $T_p = 0.351$, representando un vidrio bien recocido (duro).
   • Las muestras fueron sometidas a carga de cizalla AQS hasta una deformación global de $\gamma = 5$.
   • Las propiedades locales se caracterizaron mediante el "método de la matriz congelada", que aisla regiones circulares para medir los esfuerzos de fluencia locales ($\sigma_c$), las fluctuaciones de esfuerzo interno ($\delta\sigma$) y las caídas de esfuerzo ($\Delta\sigma$) a través de 18 direcciones de cizallamiento.

2. Modelo Mesoscópico:

   • Se implementó un modelo de red escalar en una rejilla cuadrada con condiciones de contorno bi-periódicas.
   • Dinámica: El modelo presenta una dinámica de umbral local donde una celda cede si el esfuerzo local excede un umbral $\sigma_c$. Los eventos plásticos activan una redistribución de esfuerzo elástico de largo alcance (núcleo de Eshelby) y una renovación del umbral local.
   • Parámetros Clave: El estudio se centra en dos parámetros adimensionales:
     • Fuerza del Desorden ($s$): La relación entre las fluctuaciones del esfuerzo de fluencia y la caída de esfuerzo inducida por un evento plástico.
     • Parámetro de Edad/Rejuvenecimiento ($r$): Definido como $r = (\sigma_{c}^{inicial} - \sigma_{c}^{renovado}) / \sigma_{c}^{renovado}$. Cuantifica la diferencia entre el estado vítreo inicial y el estado del material después de un reajuste plástico.
   • Estrategia de Up-scaling: Los parámetros para el modelo mesoscópico (distribuciones de esfuerzos de fluencia, incrementos de deslizamiento y distribuciones de renovación) se extrajeron directamente de los datos atomísticos. El estudio probó dos hipótesis específicas:
     • Escenario Desordenado: La distribución de renovación de los esfuerzos de fluencia es idéntica a la distribución inicial ($r=0$).
     • Escenario de Rejuvenecimiento: La distribución de renovación corresponde al vidrio MCT (un estado más "fresco"), lo que implica que $r > 0$.
     • Hipótesis de Correlación: Se introdujo una correlación directa entre el esfuerzo de fluencia local y la magnitud del incremento de deslizamiento subsiguiente ($\varepsilon_0 \propto \sigma_c^2$), reflejando la estructura jerárquica del PEL.

Contribuciones Clave y Resultados

1. Comportamiento Esfuerzo-Deformación y Transitorios:

   • Los modelos simples sin correlación entre el esfuerzo de fluencia y los incrementos de deslizamiento fallaron en reproducir el pico de esfuerzo transitorio observado en las simulaciones atomísticas, produciendo únicamente curvas monotónicas.
     implantando una correlación entre el esfuerzo de fluencia local y la amplitud del incremento de deslizamiento, el modelo mesoscópico pudo reproducir la curva completa de esfuerzo-deformación, incluyendo el pico transitorio y el esfuerzo de flujo en estado estacionario, tanto para vidrios blandos (ESL) como duros (GQ).
   • Se encontró que el esfuerzo de flujo en estado estacionario está controlado principalmente por el parámetro de fuerza del desorden $s$ (vinculado a la amplitud media del deslizamiento), más que por los detalles específicos del escenario de rejuvenecimiento.

2. Rejuvenecimiento y Evolución del Esfuerzo de Fluencia Local:

   • Vidrio Blando (ESL): Las simulaciones atomísticas mostraron un ligero endurecimiento (aumento del esfuerzo de fluencia local medio) seguido de una meseta. El modelo mesoscópico con rejuvenecimiento ($r \approx 0.14$) logró un acuerdo cuasi-cuantitativo con esta evolución. Sin el rejuvenecimiento, el modelo sobreestimaba el esfuerzo de fluencia.
   • Vidrio Duro (GQ): Las simulaciones atomísticas mostraron un endurecimiento inicial seguido de un régimen de ablandamiento asociado con la formación de bandas de cizalla. El modelo mesoscópico sin rejuvenecimiento no pudo reproducir la fase de ablandamiento, mostrando solo un endurecimiento continuo. La inclusión del escenario de rejuvenecimiento ($r \approx 0.65$) fue crucial para recuperar cualitativamente la disminución del esfuerzo de fluencia medio, aunque el acuerdo cuantitativo siguió siendo imperfecto.

3. Localización y Bandas de Cizalla:

   • Para el vidrio duro (GQ), las simulaciones atomísticas exhibieron bandas de cizalla claras.
   • El modelo mesoscópico con rejuvencimiento y correlación reprodujo con éxito la evolución del índice de localización: un aumento rápido durante la nucleación de la banda seguido de una disminución gradual a medida que la banda se ensancha.
   • Los modelos sin rejuvencimiento mostraron una localización significativa, pero subestimaron la intensidad y persistencia de la banda de cizalla.

4. Fluctuaciones de Deformación y Esfuerzo:

   • Si bien el modelo mesoscópico reprodujo la tendencia de la evolución de la varianza de la deformación, sobreestimó sistemáticamente la magnitud de las fluctuaciones por un factor de ~1.5 en comparación con los resultados atomísticos.
   • Las fluctuaciones del campo de esfuerzo interno fueron sobreestimadas aún más significativamente (por un factor de 2–3). Los autores atribuyen esta discrepancia a la naturaleza escalar del modelo (ignorando el desorden orientacional de las direcciones de deslizamiento) y a la negligencia de la heterogeneidad elástica y los efectos no afines en la formulación mesoscópica.

Significancia y Reivindicaciones
El artículo sostiene que un modelo escalar de granularidad gruesa puede reproducir las características salientes de la plasticidad amorfa, siempre que se incorpore información estructural específica de la escala atomística. La principal significancia radica en la identificación de dos ingredientes críticos para el acuerdo cuantitativo:

1. Correlación: La necesidad de correlacionar las caídas de esfuerzo locales con los esfuerzos de fluencia locales, lo que refleja la naturaleza jerárquica del Paisaje de Energía Potencial (pozos más profundos corresponden a deslizamientos mayores).
2. Rejuvenecimiento: La validación del escenario de rejuvenecimiento, donde los reajustes plásticos transforman la estructura vítrea local en un estado que se asemeja a un vidrio enfriado cerca de la Transición de Acoplamiento de Modo. Este mecanismo es esencial para reproducir el comportamiento de ablandamiento y las bandas de cizalla observados en vidrios bien recocidos.

Los autores reconocen modestamente que el modelo es "demasiado simple y caricaturesco" para un acuerdo cuantitativo completo, particularmente respecto a las amplitudes de las fluctuaciones. Concluyen que, si bien el modelo escalar captura la física esencial del flujo y la localización, el trabajo futuro debe abordar la heterogeneidad elástica, el desorden tensorial y los efectos de memoria para resolver completamente las discrepancias en la estadística de las fluctuaciones.