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

Résumé Technique : Comment le verre se brise – L'endommagement explique la différence entre l'énergie de surface et l'énergie de fracture dans la silice amorphe

Énoncé du problème
En mécanique de la rupture, une divergence fondamentale existe entre l'énergie de surface libre ($\gamma$), qui quantifie l'énergie requise pour créer de nouvelles surfaces, et l'énergie de fracture ($G_c$), qui représente l'énergie totale dissipée lors de la propagation d'une fissure. Pour les matériaux fragiles idéaux, $G_c$ devrait être égal à $2\gamma$. Cependant, les mesures expérimentales sur les verres de silicate montrent systématiquement que $G_c$ dépasse $2\gamma$ d'un facteur d'environ 5. Cette énergie excédentaire est souvent attribuée à la dissipation plastique à la pointe de la fissure, mais cette explication reste controversée pour la silice amorphe. L'obtention de preuves de plasticité est difficile expérimentalement à l'échelle atomique, et les modèles de plasticité standards ne s'appliquent pas aux silicates en raison de leur absence d'écrouissage. Par conséquent, l'origine physique de cette divergence énergétique dans les silicates amorphes fragiles est restée non résolue.

Méthodologie
Les auteurs ont employé une approche de simulation multi-échelle combinant la dynamique moléculaire (MD) à l'échelle atomique avec une modélisation par champ de phase (phase-field) continue pour étudier la fracture de la silice amorphe.

• Simulations atomiques : Les simulations MD ont été réalisées en utilisant le potentiel BKS pour modéliser les interactions atomiques et la troncature de Wolf pour les forces coulombiennes. Les échantillons ont été générés par placement aléatoire d'atomes, équilibrés, puis refroidis à $10^{-5}$ K. La propagation de la fissure a été simulée sous des conditions quasi-statiques et athermiques en utilisant une condition limite de déplacement de type champ-K pour piloter la croissance de la fissure.
• Moyennage (Coarse-Graining) : Pour faire le pont entre les échelles atomique et continue, une technique de convolution gaussienne basée sur la physique a été appliquée. Cette méthode a généré des champs homogènes de déplacement, de déformation, de contrainte et d'énergie potentielle tout en conservant la masse et l'énergie.
• Quantification de l'endommagement et de l'énergie :
  • Endommagement : Défini mécaniquement comme la perte de rigidité, calculé en comparant l'énergie élastique réelle à l'énergie de déformation non endommagée.
  • Énergie de surface libre (FSE) : Isolée en calculant la différence entre l'énergie potentielle totale et l'énergie de déformation élastique ($\psi_{FSE} = \psi_{pot} - \psi_{\epsilon}$), séparant ainsi l'énergie de formation de surface de la déformation élastique.
  • Calibrage par champ de phase : Un schéma de mise à jour par éléments finis (FEMU) a été utilisé pour ajuster le champ d'endommagement obtenu par MD à un modèle de champ de phase continu. Cela a permis d'extraire les paramètres de matériau, spécifiquement la ténacité de rupture ($g_c$) et l'échelle de longueur interne ($l_c$), directement à partir de clichés atomiques.
• Analyse structurelle : Les changements structurels locaux ont été quantifiés en suivant l'évolution des structures de cycles (en utilisant l'algorithme de Le Roux & Jund) et le nombre de coordination des atomes de silicium.

Résultats clés

1. Endommagement diffus vs Plasticité : L'étude démontre que la dissipation d'énergie autour de la pointe de la fissure est principalement due à un endommagement diffus plutôt qu'à la plasticité. La zone d'endommagement s'étend sur environ 16 à 23 Å à partir du chemin de la fissure. L'analyse de la réponse contrainte-déformation et du bilan énergétique confirme que la dissipation plastique et l'élasticité non linéaire sont négligeables ; le matériau perd sa rigidité (endommagement) plutôt que de la conserver (plasticité).
2. Séparation des contributions énergétiques : Les auteurs ont réussi à séparer l'énergie de fracture totale en deux composantes distinctes :
   • Énergie de surface libre : Localisée strictement à la surface, associée à la rupture des liaisons et à la relaxation atomique.
   • Énergie d'endommagement : S'étendant plus profondément dans le matériau (~20 Å), associée aux réarrangements structurels.
     Le rapport entre l'énergie de surface libre et l'énergie de fracture totale ($2\gamma/g_c$) a été calculé à environ 23 %, ce qui est cohérent avec les observations expérimentales où l'énergie de surface ne représente qu'environ 20 % de l'énergie totale.
3. Signatures structurelles :
   • Endommagement diffus : Corrèle fortement avec les changements dans la structure de cycle du réseau de verre, qui possède une étendue spatiale plus large.
   • Énergie de surface libre : Corrèle avec les changements du nombre de coordination des atomes de silicium, qui sont plus localisés près de la surface.
4. Validation : Les valeurs de ténacité de rupture dérivées de l'ajustement par champ de phase, du bilan énergétique global et de la méthode de l'intégrale J étaient cohérentes. L'échelle de longueur interne calculée ($l_c \approx 12\text{--}20$ Å) concorde avec la largeur observée de la zone d'endommagement.

Signification et affirmations
L'article affirme fournir un cadre pleinement cohérent pour comprendre la fracture dans la silice amorphe en résolvant la divergence de longue date entre les énergies de surface et de fracture. La contribution primaire est l'identification de l'endommagement diffus — distinct de la plasticité — comme le mécanisme dominant de dissipation d'énergie dans les verres de silicate fragiles. En liant cet endommagement à des réarrangements structurels spécifiques (rupture de la structure de cycle) et en quantifiant son étendue spatiale, l'étude offre une explication structurelle de la nature à la fois « fragile » et dissipative de la fracture de la silice.

Les auteurs affirment que leur approche valide le concept d'endommagement diffus utilisé dans les techniques de fracture par champ de phase, fournissant une démonstration directe de sa base physique à l'échelle atomique. De plus, les résultats suggèrent que les vues traditionnelles de la fracture fragile pourraient négliger la subtile dissipation mécanique inhérente aux matériaux amorphes. L'article conclut avec modestie, notant que si ces conclusions clarifient le mécanisme dans les silicates, des recherches supplémentaires sont nécessaires pour déterminer si ces concepts s'appliquent à d'autres systèmes amorphes, tels que les verres métalliques, ou aux scénarios de chargement dynamique et de fatigue.