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

En utilisant des simulations de dynamique moléculaire avec un potentiel appris par machine, cette étude révèle que l'énergie de fracture du verre de silice augmente jusqu'à 33 % en dessous du seuil de ramification en raison d'une combinaison de l'augmentation de la densité d'énergie de surface intrinsèque et du rugosissement à l'échelle nanométrique, démontrant ainsi que la fracture dynamique crée une structure de surface fondamentalement différente plutôt que d'augmenter simplement la surface apparente.

L'essentiel

Imaginez que vous essayez de casser un morceau de verre épais. Vous pourriez penser que l'énergie nécessaire pour le briser est simplement celle requise pour rompre les minuscules liaisons atomiques qui le maintiennent ensemble, comme couper un seul brin de spaghetti. Cependant, les scientifiques savent depuis longtemps que briser du verre nécessite en réalité beaucoup plus d'énergie que ce que ce calcul simple suggère. C'est comme si le verre se défendait, exigeant un effort supplémentaire pour se briser.

Pendant des années, les chercheurs ont cru que ce « coût supplémentaire » était principalement dû au fait que la fissure devenait vacillante et irrégulière en accélérant, créant une surface plus rugueuse (comme déchirer un morceau de papier en une bande irrégulière plutôt qu'en une ligne droite). Mais une nouvelle étude utilisant des simulations informatiques avancées a révélé une histoire plus complexe.

Voici ce que l'article a découvert, expliqué simplement :

1. L'extrémité de la fissure « surchauffée »

Lorsqu'une fissure se déplace très rapidement à travers le verre, l'extrémité de cette fissure devient incroyablement chaude. L'étude a révélé qu'à haute vitesse, les atomes situés juste à l'extrémité de la fissure atteignent des températures de 8 000 Kelvin (plus chaudes que la surface du soleil !).

Imaginez l'extrémité de la fissure non pas seulement comme un point de rupture, mais comme un minuscule chalumeau microscopique. Cette chaleur intense ne fait pas simplement fondre le verre ; elle modifie fondamentalement la nature de la surface en cours de création.

2. Deux raisons pour lesquelles le verre coûte plus cher à briser

Les chercheurs ont utilisé un modèle informatique ultra-précis (comme un microscope numérique capable de voir les atomes individuels) pour déterminer où toute cette énergie supplémentaire est dépensée. Ils ont découvert que le « coût supplémentaire » est réparti approximativement 50/50 entre deux facteurs :

• Le facteur « Rugosité » (Quantité) : Lorsque la fissure accélère, la surface qu'elle laisse derrière elle n'est pas parfaitement lisse. Elle devient nanoscopiquement rugueuse, comme une chaîne de montagnes vue depuis l'espace. Cela signifie que la fissure crée en réalité plus de surface que ce qui apparaît de l'extérieur.
  • Analogie : Imaginez déchirer un morceau de tissu. Si vous le déchirez lentement, le bord est droit. Si vous le déchirez rapidement, le bord devient effiloché et irrégulier. Vous avez utilisé plus de tissu pour créer ce bord irrégulier.
• Le facteur « Qualité » (Densité d'énergie) : C'est la nouvelle découverte. Même si vous lissiez cette surface irrégulière, il en coûterait toujours plus d'énergie pour la créer que pour une surface calme et lente. La chaleur extrême à l'extrémité de la fissure modifie la structure atomique de la nouvelle surface, la rendant « à plus haute énergie » ou plus instable.
  • Analogie : Imaginez cuire un gâteau. Un gâteau cuit lentement a une texture standard. Mais si vous le bombardez avec un chalumeau, l'extérieur devient carbonisé et chimiquement différent. La surface « carbonisée » est fondamentalement différente et nécessite plus d'énergie pour être créée que la version lisse et cuite lentement.

3. La « rugosité cachée »

L'un des points les plus intéressants est que la « rugosité » découverte par l'ordinateur est si minuscule (à l'échelle des atomes) que les outils standards utilisés par les ingénieurs pour mesurer le verre brisé la manqueraient complètement.

Si vous regardiez un morceau de verre brisé avec un microscope normal, vous verriez une surface lisse. Vous supposeriez que toute l'énergie supplémentaire a servi à rendre la surface « plus chaude » ou plus énergétique. Mais cette étude montre qu'une part importante de cette énergie a en réalité servi à rendre la surface physiquement plus grande et plus rugueuse, simplement à une échelle trop petite pour que nos yeux puissent la voir.

4. Corriger les mathématiques sur la vitesse de propagation des fissures

L'article a également corrigé une formule de longue date utilisée pour prédire la vitesse de propagation d'une fissure en fonction de la force appliquée. L'ancienne formule (le « modèle Freund ») était comme une carte qui devenait un peu floue à haute vitesse. La nouvelle étude a trouvé une meilleure formule (une « relation en racine carrée ») qui correspond parfaitement aux données.

Cette correction est importante car elle aide à expliquer pourquoi les expériences précédentes mesurant la chaleur dégagée par la rupture du verre (en utilisant la lumière émise par la fissure, appelée fractoluminescence) ne correspondaient pas tout à fait aux prédictions de vitesse. En utilisant la nouvelle formule, les vitesses et températures prédites correspondent enfin à ce que les simulations informatiques ont montré.

L'essentiel

Briser du verre ne consiste pas seulement à rompre des liaisons. Lorsque la fissure se déplace rapidement, elle agit comme un minuscule laser surchauffé qui :

1. Rend la surface physiquement plus rugueuse (créant plus de surface).
2. Modifie chimiquement la surface pour la rendre plus énergétique.

L'étude prouve que l'énergie nécessaire pour briser le verre n'est pas un nombre fixe ; elle change en fonction de la vitesse à laquelle vous le brisez, et elle est pilotée à la fois par la forme de la rupture et la chaleur extrême à l'extrémité.

Résumé technique

Énoncé du problème
L'énergie de fracture des matériaux fragiles, tels que le verre de silice, augmente avec la vitesse de fissure, un phénomène qui dépasse le coût thermodynamique de la rupture des liaisons. Bien que cette dépendance à la vitesse soit souvent attribuée au rugosité de surface causée par des instabilités dynamiques (telles que la micro-ramification), le comportement en dessous du seuil de ramification reste mal compris. Des études antérieures sur les polymères (par exemple, le PMMA) ont noté une augmentation d'environ 30 % de l'énergie de fracture avant le début de la ramification, l'attribuant à des processus de dissipation non résolus. Une question centrale persiste : cette augmentation d'énergie pré-ramification provient-elle d'une augmentation de la quantité de surface (via un rugosité nanométrique invisible aux mesures standard) ou d'une augmentation de la qualité (densité d'énergie intrinsèque) de la surface créée due aux conditions extrêmes à la pointe de la fissure ? De plus, les potentiels interatomiques empiriques utilisés dans les études antérieures de dynamique moléculaire (DM) sur la silice ont échoué à capturer avec précision les événements de rupture de liaisons à haute énergie, conduisant à une sous-estimation des vitesses de ramification et à des descriptions imprécises de la zone de processus.

Méthodologie
Pour remédier à ces limitations, les auteurs ont employé des simulations étendues de dynamique moléculaire à partir des premiers principes du verre de silice fondue.

• Potentiel interatomique : Au lieu de potentiels empiriques, l'étude a utilisé un Potentiel Interatomique Appris par Machine (MLIP) basé sur l'architecture Allegro. Ce modèle a été entraîné sur des données de la théorie de la fonctionnelle de la densité (DFT) utilisant le fonctionnel r2SCAN, couvrant un vaste éventail d'états incluant la silice fondue, les structures vitreuses et la déformation mécanique jusqu'à 15 000 K.
• Configuration de la simulation : Le système consistait en une bande de verre de silice de 12,3 millions d'atomes (300 × 8,5 × 75 nm³) avec une entaille préexistante. Les simulations ont été menées à 300 K et 1 bar, suivies d'un déplacement de traction rapide pour induire la propagation de la fissure.
• Paramètres : Quatorze simulations ont été exécutées avec des taux de déplacement variables, entraînant des contraintes de traction initiales comprises entre 3,9 GPa et 7 GPa.
• Analyse : Les auteurs ont distingué l'"énergie de fracture structurelle" (énergie stockée dans la structure endommagée après refroidissement, normalisée par la surface projetée) et la "densité d'énergie de surface intrinsèque" (énergie normalisée par la surface réelle, rugueuse, mesurée via la tessellation de Delaunay). Ils ont également analysé les températures à la pointe de la fissure en moyennant l'énergie cinétique sur des intervalles de 0,5 ps.

Résultats clés

1. Énergie de fracture dépendante de la vitesse : Même en dessous du seuil de ramification (défini à 0,72 de la vitesse de l'onde de Rayleigh, $v_R$), l'énergie de fracture structurelle augmente jusqu'à 33 % à mesure que la vitesse de fissure augmente.
2. Décomposition de l'augmentation d'énergie : Cette hausse est répartie à peu près également entre deux mécanismes :
   • Rugosité nanométrique : Une augmentation d'environ 16 % de la surface réelle de fracture ($A_{real}$) par rapport à la surface projetée. Ce rugosité crée une zone d'endommagement plus épaisse autour de la fissure mais se produit à des échelles de longueur invisibles pour l'analyse optique post-mortem standard.
   • Densité d'énergie de surface intrinsèque : Une augmentation d'environ 15 % de la densité d'énergie de surface intrinsèque ($\gamma_s$). Cela est attribué aux conditions thermiques extrêmes à la pointe de la fissure, où les températures dépassent 3 000 K dans une zone de processus de seulement 3 à 6 nm.
3. Vitesse de fissure et instabilité : Les simulations révèlent que la vitesse de fissure suit une relation racine carrée avec le facteur d'intensité de contrainte ($v_c \propto \sqrt{1 - (K_{Ic}/K_I)^2}$), différant du modèle Freund standard. La ramification se produit à environ 2,4 km/s (0,72 $v_R$), ce qui est cohérent avec les observations expérimentales sur la silice fondue.
4. Taille de la zone de processus : Le rayon de la zone de processus augmente avec l'intensité de la contrainte, variant d'environ 3,1 nm pour les fissures planes à environ 6,8 nm pour les fissures ramifiées. Cette taille s'aligne sur les prédictions de Dugdale–Barenblatt et est inférieure aux valeurs dérivées des simulations antérieures utilisant des potentiels empiriques (~10 nm).
5. Résolution de la divergence thermique : L'étude résout une divergence entre les températures de simulation et les prédictions expérimentales basées sur la fractoluminescence. En appliquant la nouvelle relation vitesse–contrainte à racine carrée identifiée, les températures prédites à la pointe de la fissure pour une vitesse donnée sont plus basses, rapprochant les modèles analytiques des données de simulation.

Signification et revendications
L'article revendique que la fracture dynamique dans le verre de silice augmente l'énergie de fracture non seulement en créant plus de surface apparente, mais en altérant fondamentalement la surface à l'échelle nanométrique. Les auteurs démontrent que :

• L'énergie de fracture n'est pas une propriété matérielle constante ; elle dépend de la vitesse même en l'absence de ramification macroscopique.
• L'analyse post-mortem standard confondrait les contributions du rugosité nanométrique et de l'augmentation de la densité d'énergie intrinsèque, interprétant l'ensemble de l'augmentation comme une hausse de la densité d'énergie de surface.
• Les simulations DM à partir des premiers principes utilisant des MLIP peuvent désormais atteindre une précision quantitative dans la modélisation de la fracture dynamique sans réglage empirique, reproduisant avec succès les propriétés élastiques, la ténacité à la rupture et les seuils de ramification qui correspondent aux données expérimentales.
• La relation vitesse–contrainte révisée fournit un cadre plus précis pour interpréter les températures à la pointe de la fissure déduites des expériences de fractoluminescence.

L'étude conclut que l'augmentation pré-ramification de l'énergie de fracture, précédemment observée mais inexpliquée dans des matériaux comme le PMMA, provient d'une combinaison de rugosité nanométrique et d'une densité d'énergie de surface intrinsèque élevée, pilotée par des températures locales extrêmes.