How glass breaks -- Damage explains the difference between surface and fracture energies in amorphous silica

Resumen Técnico: Cómo se rompe el vidrio – El daño explica la diferencia entre las energías de superficie y de fractura en el sílice amorfo

Planteamiento del problema
En la mecánica de la fractura, existe una discrepancia fundamental entre la energía de superficie libre ($\gamma$), que cuantifica la energía necesaria para crear nuevas superficies, y la energía de fractura ($G_c$), que representa la energía total disipada durante la propagación de una grieta. En materiales frágiles ideales, $G_c$ debería ser igual a $2\gamma$. Sin embargo, las mediciones experimentales en vidrios de silicato muestran consistentemente que $G_c$ excede a $2\gamma$ por un factor de aproximadamente 5. Esta energía excedente se atribuye a menudo a la disipación plástica en la punta de la grieta, pero esta explicación sigue siendo controvertida para el sílice amorfo. La evidencia de plasticidad es difícil de obtener experimentalmente a escala atómica, y los modelos de plasticidad estándar no se pueden aplicar a los silicatos debido a su falta de endurecimiento. En consecuencia, el origen físico de esta discrepancia de energía en los silicatos frágiles amorfos ha permanecido sin resolver.

Metodología
Los autores emplearon un enfoque de simulación multiescala combinando dinámica molecular (MD) a escala atómica con un modelado de campo de fase continuo para investigar la fractura en el sílice amorfo.

• Simulaciones Atómicas: Se realizaron simulaciones de MD utilizando el potencial BKS para modelar las interacciones atómicas y la truncación de Wolf para las fuerzas coulombianas. Las muestras fueron generadas mediante la colocación aleatoria de átomos, equilibradas y enfriadas a $10^{-5}$ K. La propagación de la fractura se simuló bajo condiciones cuasiestáticas y atermales utilizando una condición de contorno de desplazamiento de campo K para impulsar el crecimiento de la grieta.
• Granularidad (Coarse-Graining): Para tender un puente entre las escalas atómica y continua, se aplicó una técnica de convolución gaussiana basada físicamente. Este método generó campos homogéneos de desplazamiento, deformación, tensión y energía potencial, conservando la masa y la energía.
• Cuantificación de Daño y Energía:
  • Daño: Definido mecánicamente como la pérdida de rigidez, calculado comparando la energía elástica real con la energía de deformación no dañada.
  • Energía de Superficie Libre (FSE): Aislada mediante el cálculo de la diferencia entre la energía potencial total y la energía de deformación elástica ($\psi_{FSE} = \psi_{pot} - \psi_{\epsilon}$), separando efectivamente la energía de formación de la superficie de la deformación elástica.
  • Calibración de Campo de Fase: Se utilizó un esquema de Actualización de Elementos Finitos (FEMU) para ajustar el campo de daño obtenido de las simulaciones de MD a un modelo de campo de fase continuo. Esto permitió la extracción de parámetros de material, específicamente la tenacidad a la fractura ($g_c$) y la escala de longitud interna ($l_c$), directamente de instantáneas atómicas.
• Análisis Estructural: Los cambios estructurales locales se cuantificaron rastreando la evolución de las estructuras de anillos (usando el algoritmo de Le Roux & Jund) y el número de coordinación de los átomos de silicio.

Resultos Clave

1. Daño Difuso vs. Plasticidad: El estudio demuestra que la disipación de energía alrededor de la punta de la grieta se debe principalmente al daño difuso en lugar de a la plasticidad. La zona de daño se extiende aproximadamente de 16 a 23 Å desde la trayectoria de la grieta. El análisis de la respuesta tensión-deformación y del balance de energía confirma que la disipación plástica y la elasticidad no lineal son insignificantes; el material pierde rigidez (daño) en lugar de retenerla (plasticidad).
2. Separación de las Contribuciones de Energía: Los autores separaron con éxito la energía de fractura total en dos componentes distintos:
   • Energía de Superficie Libre: Localizada estrictamente en la superficie, asociada con la ruptura de enlaces y la relajación atómica.
   • Energía de Daño: Extendiéndose más profundamente en el material (~20 Å), asociada con reordenamientos estructurales.
     La relación entre la energía de superficie libre y la energía de fractura total ($2\gamma/g_c$) se calculó en aproximadamente un 23%, lo cual es consistente con las observaciones experimentales donde la energía superficial representa solo cerca del 20% de la energía total.
3. Firmas Estructurales:
   • Daño Difuso: Correlaciona fuertemente con los cambios en la estructura de anillos de la red vítrea, la cual posee una extensión espacial más amplia.
   • Energía de Superficie Libre: Correlaciona con los cambios en el número de coordinación de los átomos de silicio, que están más localizados cerca de la superficie.
4. Validación: Los valores de tenacidad a la fractura derivados del ajuste de campo de fase, el balance de energía global y el método de la integral J fueron consistentes. La escala de longitud interna calculada ($l_c \approx 12\text{--}20$ Å) se alinea con el ancho observado de la zona de daño.

Significado y Reivindicaciones
El artículo afirma proporcionar un marco totalmente consistente para comprender la fractura en el sílice amorfo al resolver la discrepancia de larga data entre las energías de superficie y de fractura. La contribución principal es la identificación del daño difuso —distinto de la plasticidad— como el mecanismo dominante para la disipación de energía en los vidrios de silicato frágiles. Al vincular este daño con reordenamientos estructurales específicos (ruptura de la estructura de anillos) y cuantificar su extensión espacial, el estudio ofrece una explicación estructural para la naturaleza "frágil" pero disipativa de la fractura del sílice.

Los autores aseveran que su enfoque valida el concepto de daño difuso utilizado en las técnicas de fractura de campo de fase, proporcionando una demostración directa de su base física a la escala atómica. Además, los resultados sugieren que las visiones tradicionales de la fractura frágil pueden pasar por alto la sutil disipación mecánica inherente en los materiales amorfos. El artículo concluye modestamente, señalando que, si bien estos hallazgos aclaran el mecanismo en los silicatos, se requiere mayor investigación para determinar si estos conceptos se aplican a otros sistemas amorfos, como los vidrios metálicos, o a escenarios de carga dinámica y de fatiga.