Energy dissipation at the atomic scale explains how fracture energy depends on crack velocity in silica glass

Mediante simulaciones de dinámica molecular con un potencial aprendido por máquina, este estudio revela que la energía de fractura del vidrio de sílice aumenta hasta un 33% por debajo del umbral de ramificación debido a una combinación del aumento de la densidad de energía superficial intrínseca y la rugosidad a escala nanométrica, demostrando que la fractura dinámica crea una estructura superficial fundamentalmente diferente en lugar de simplemente aumentar el área superficial aparente.

Lo esencial

Imagina que intentas romper un trozo de vidrio grueso. Podrías pensar que la energía necesaria para romperlo es simplemente la energía requerida para cortar los diminutos enlaces atómicos que lo mantienen unido, como cortar un solo hilo de espagueti. Sin embargo, los científicos han sabido durante mucho tiempo que romper vidrio en realidad requiere mucho más energía de lo que sugiere ese cálculo simple. Es como si el vidrio se defendiera, exigiendo un esfuerzo extra para romperse.

Durante años, los investigadores creyeron que este "costo extra" se debía principalmente a que la grieta se volvía inestable y dentada a medida que aceleraba, creando un área de superficie más rugosa (como rasgar un trozo de papel en una tira dentada en lugar de una línea recta). Pero un nuevo estudio que utiliza simulaciones por computadora avanzadas ha revelado una historia más compleja.

Esto es lo que descubrió el artículo, explicado de forma sencilla:

1. La punta de la grieta "sobrecalentada"

Cuando una grieta se mueve muy rápido a través del vidrio, la punta de esa grieta se vuelve increíblemente caliente. El estudio encontró que a altas velocidades, los átomos justo en la punta de la grieta alcanzan temperaturas de 8.000 Kelvin (más calientes que la superficie del sol).

Piensa en la punta de la grieta no solo como un punto de ruptura, sino como un pequeño soplete microscópico. Este calor intenso no solo derrite el vidrio; cambia fundamentalmente la naturaleza de la superficie que se está creando.

2. Dos razones por las que cuesta más romper el vidrio

Los investigadores utilizaron un modelo informático superpreciso (como un microscopio digital que ve átomos individuales) para averiguar a dónde va toda esa energía extra. Descubrieron que el "costo extra" se divide aproximadamente 50/50 entre dos cosas:

• El factor de "Rugosidad" (Cantidad): A medida que la grieta acelera, la superficie que deja atrás no es perfectamente lisa. Se vuelve rugosa a escala nanoscópica, como una cordillera vista desde el espacio. Esto significa que la grieta en realidad crea más área de superficie de la que parece desde el exterior.
  • Analogía: Imagina rasgar un trozo de tela. Si lo rasgas lentamente, el borde es recto. Si lo rasgas rápido, el borde se deshilacha y se vuelve dentado. Has usado más tela para crear ese borde dentado.
• El factor de "Calidad" (Densidad de energía): Este es el nuevo descubrimiento. Incluso si alisaras esa superficie dentada, aún costaría más energía crearla que una superficie tranquila y lenta. El calor extremo en la punta de la grieta cambia la estructura atómica de la nueva superficie, haciéndola de "mayor energía" o más inestable.
  • Analogía: Imagina hornear un pastel. Un pastel cocinado lentamente tiene una textura estándar. Pero si lo bombardeas con un soplete, el exterior se quema y cambia químicamente. La superficie "quemada" es fundamentalmente diferente y requiere más energía para crearse que la versión lisa y cocinada lentamente.

3. La rugosidad "oculta"

Uno de los puntos más interesantes es que la "rugosidad" que encontró la computadora es tan pequeña (a escala atómica) que las herramientas estándar utilizadas por los ingenieros para medir el vidrio roto la pasarían por completo.

Si miraras un trozo de vidrio roto con un microscopio normal, verías una superficie lisa. Asumirías que toda la energía extra se fue a hacer la superficie "más caliente" o más energética. Pero este estudio muestra que una parte significativa de esa energía en realidad se fue a hacer la superficie físicamente más grande y más rugosa, solo que a una escala demasiado pequeña para que nuestros ojos la vean.

4. Corrigiendo las matemáticas sobre la velocidad de las grietas

El artículo también corrigió una fórmula de larga data utilizada para predecir la velocidad a la que se mueve una grieta basándose en la fuerza aplicada. La fórmula antigua (el "modelo de Freund") era como un mapa que se volvía un poco borroso a altas velocidades. El nuevo estudio encontró una fórmula mejor (una "relación de raíz cuadrada") que se ajusta perfectamente a los datos.

Esta corrección es importante porque ayuda a explicar por qué los experimentos anteriores que medían el calor de la rotura del vidrio (usando la luz emitida por la grieta, llamada fractoluminiscencia) no coincidían del todo con las predicciones de velocidad. Al usar la nueva fórmula, las velocidades y temperaturas predichas finalmente coinciden con lo que mostraron las simulaciones por computadora.

La conclusión

Romper vidrio no se trata solo de cortar enlaces. Cuando la grieta se mueve rápido, actúa como un pequeño láser sobrecalentado que:

1. Hace que la superficie sea físicamente más rugosa (creando más área).
2. Altera químicamente la superficie para hacerla más energética.

El estudio demuestra que la energía necesaria para romper vidrio no es un número fijo; cambia dependiendo de la velocidad a la que lo rompes, y está impulsada tanto por la forma de la ruptura como por el calor extremo en la punta.

Resumen técnico

Enunciado del Problema
Se observa que la energía de fractura de materiales frágiles, como el vidrio de sílice, aumenta con la velocidad de la grieta, un fenómeno que excede el costo termodinámico de la ruptura de enlaces. Si bien esta dependencia de la velocidad se atribuye a menudo al alisado superficial causado por inestabilidades dinámicas (como la microramificación), el comportamiento por debajo del umbral de ramificación sigue siendo poco comprendido. Estudios previos sobre polímeros (por ejemplo, PMMA) notaron un aumento de aproximadamente el 30% en la energía de fractura antes del inicio de la ramificación, atribuyéndolo a procesos de disipación no resueltos. Persiste una pregunta central: ¿este aumento de energía pre-ramificación surge de un aumento en la cantidad de área superficial (mediante rugosidad a nanoescala invisible para las mediciones estándar) o de un aumento en la calidad (densidad de energía intrínseca) de la superficie creada debido a condiciones extremas en la punta de la grieta? Además, los potenciales interatómicos empíricos utilizados en estudios previos de dinámica molecular (DM) de sílice han fallado en capturar con precisión eventos de ruptura de enlaces de alta energía, lo que ha llevado a velocidades de ramificación subestimadas y descripciones inexactas de la zona de proceso.

Metodología
Para abordar estas limitaciones, los autores emplearon extensas simulaciones de dinámica molecular de primeros principios de vidrio de sílice fundida.

• Potencial Interatómico: En lugar de potenciales empíricos, el estudio utilizó un Potencial Interatómico Aprendido por Máquina (MLIP) basado en la arquitectura Allegro. Este modelo fue entrenado con datos de Teoría del Funcional de la Densidad (DFT) utilizando el funcional r2SCAN, abarcando una vasta gama de estados que incluyen sílice fundida, estructuras vítreas y deformación mecánica hasta 15.000 K.
• Configuración de la Simulación: El sistema consistió en una tira de vidrio de sílice de 12,3 millones de átomos (300 × 8,5 × 75 nm³) con una entalla preexistente. Las simulaciones se realizaron a 300 K y 1 bar, seguidas de un desplazamiento de extensión rápido para inducir la propagación de la grieta.
• Parámetros: Se ejecutaron catorce simulaciones con tasas de desplazamiento variables, resultando en tensiones de tracción iniciales entre 3,9 GPa y 7 GPa.
• Análisis: Los autores distinguieron entre "energía de fractura estructural" (energía almacenada en la estructura dañada después del enfriamiento, normalizada por el área proyectada) y "densidad de energía superficial intrínseca" (energía normalizada por el área superficial real y rugosa medida mediante teselación de Delaunay). También analizaron las temperaturas en la punta de la grieta promediando la energía cinética en intervalos de 0,5 ps.

Resultados Clave

1. Energía de Fractura Dependiente de la Velocidad: Incluso por debajo del umbral de ramificación (definido en 0,72 de la velocidad de la onda de Rayleigh, $v_R$), la energía de fractura estructural aumenta hasta un 33% a medida que aumenta la velocidad de la grieta.
2. Descomposición del Aumento de Energía: Este aumento se divide aproximadamente por igual entre dos mecanismos:
   • Rugosidad a Nanoescala: Un aumento de aproximadamente el 16% en el área de la superficie de fractura real ($A_{real}$) en comparación con el área proyectada. Esta rugosidad crea una zona de daño más gruesa alrededor de la grieta, pero ocurre a escalas de longitud invisibles para el análisis óptico post-mortem estándar.
   • Densidad de Energía Superficial Intrínseca: Un aumento de aproximadamente el 15% en la densidad de energía superficial intrínseca ($\gamma_s$). Esto se atribuye a las condiciones térmicas extremas en la punta de la grieta, donde las temperaturas superan los 3.000 K dentro de una zona de proceso de solo 3–6 nm.
3. Velocidad de la Grieta e Inestabilidad: Las simulaciones revelan que la velocidad de la grieta sigue una relación de raíz cuadrada con el factor de intensidad de tensión ($v_c \propto \sqrt{1 - (K_{Ic}/K_I)^2}$), diferenciándose del modelo estándar de Freund. La ramificación ocurre a aproximadamente 2,4 km/s (0,72 $v_R$), consistente con las observaciones experimentales en sílice fundida.
4. Tamaño de la Zona de Proceso: El radio de la zona de proceso aumenta con la intensidad de la tensión, oscilando entre 3,1 nm para grietas planas y ~6,8 nm para grietas con ramificación. Este tamaño se alinea con las predicciones de Dugdale–Barenblatt y es menor que los valores derivados de simulaciones previas con potenciales empíricos (10 nm).
5. Resolución de la Discrepancia Térmica: El estudio resuelve una discrepancia entre las temperaturas de simulación y las predicciones experimentales basadas en fractoluminiscencia. Al aplicar la nueva relación velocidad-tensión de raíz cuadrada identificada, las temperaturas predichas en la punta de la grieta para una velocidad dada son más bajas, acercando los modelos analíticos a los datos de simulación.

Significado y Afirmaciones
El artículo afirma que la fractura dinámica en el vidrio de sílice aumenta la energía de fractura no simplemente creando más área superficial aparente, sino alterando fundamentalmente la superficie a nanoescala. Los autores demuestran que:

• La energía de fractura no es una propiedad material constante; depende de la velocidad incluso en ausencia de ramificación macroscópica.
• El análisis post-mortem estándar confundiría las contribuciones de la rugosidad a nanoescala y el aumento de la densidad de energía intrínseca, interpretando todo el aumento como un incremento en la densidad de energía superficial.
• Las simulaciones de DM de primeros principios utilizando MLIPs pueden ahora lograr precisión cuantitativa en el modelado de la fractura dinámica sin ajuste empírico, reproduciendo con éxito las propiedades elásticas, la tenacidad a la fractura y los umbrales de ramificación que coinciden con los datos experimentales.
• La relación velocidad-tensión revisada proporciona un marco más preciso para interpretar las temperaturas en la punta de la grieta inferidas de experimentos de fractoluminiscencia.

El estudio concluye que el aumento pre-ramificación en la energía de fractura, observado previamente pero sin explicación en materiales como el PMMA, surge de una combinación de rugosidad a nanoescala y una densidad de energía superficial intrínseca elevada impulsada por temperaturas locales extremas.