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

Mittels Molekulardynamiksimulationen mit einem maschinell erlernten Potenzial zeigt diese Studie, dass die Bruchenergie von Silikatglas unterhalb des Verzweigungsschwellenwerts um bis zu 33 % ansteigt, bedingt durch eine Kombination aus steigender intrinsischer Oberflächenenergiedichte und nanoskaliger Rauheit, was belegt, dass dynamische Rissausbreitung eine grundlegend andere Oberflächenstruktur erzeugt, statt lediglich die scheinbare Oberfläche zu vergrößern.

Das Wesentliche

Stellen Sie sich vor, Sie versuchen, ein Stück dickes Glas zu zerbrechen. Man könnte denken, dass die Energie, die zum Brechen erforderlich ist, einfach die Energie entspricht, die benötigt wird, um die winzigen atomaren Bindungen zu durchtrennen, die es zusammenhalten – ähnlich wie das Durchschneiden eines einzelnen Spaghetti-Strangs. Wissenschaftler wissen jedoch schon lange, dass das Brechen von Glas tatsächlich viel mehr Energie erfordert, als diese einfache Rechnung nahelegt. Es ist, als würde das Glas sich wehren und zusätzliche Anstrengung verlangen, um gebrochen zu werden.

Jahrelang glaubten Forscher, dass diese „zusätzlichen Kosten" hauptsächlich darauf zurückzuführen waren, dass der Riss beim Beschleunigen wackelig und gezackt wurde und so eine rauere Oberfläche erzeugte (wie das Zerreißen eines Papierstücks in einen gezackten Streifen statt in eine gerade Linie). Doch eine neue Studie mit fortschrittlichen Computersimulationen hat eine komplexere Geschichte enthüllt.

Hier ist das, was die Studie entdeckt hat, einfach erklärt:

1. Die „überhitzte" Rissspitze

Wenn ein Riss sehr schnell durch Glas läuft, wird die Spitze dieses Risses unglaublich heiß. Die Studie ergab, dass bei hohen Geschwindigkeiten die Atome direkt an der Rissspitze Temperaturen von 8.000 Kelvin erreichen (heißer als die Oberfläche der Sonne!).

Stellen Sie sich die Rissspitze nicht nur als Bruchpunkt vor, sondern als winzigen, mikroskopischen Brenner. Diese intensive Hitze schmilzt das Glas nicht nur; sie verändert grundlegend die Natur der entstehenden Oberfläche.

2. Zwei Gründe, warum Glas mehr Energie zum Brechen benötigt

Die Forscher nutzten ein hochpräzises Computermodell (wie ein digitales Mikroskop, das einzelne Atome sieht), um herauszufinden, wohin all diese zusätzliche Energie fließt. Sie stellten fest, dass die „zusätzlichen Kosten" ungefähr 50/50 zwischen zwei Faktoren aufgeteilt sind:

• Der „Rauheits"-Faktor (Menge): Wenn der Riss schneller wird, ist die hinterlassene Oberfläche nicht perfekt glatt. Sie wird nanoskopisch rau, wie ein Gebirgszug aus dem Weltraum betrachtet. Das bedeutet, dass der Riss tatsächlich mehr Oberfläche erzeugt, als es von außen scheint.
  • Vergleich: Stellen Sie sich vor, Sie reißen ein Stück Stoff. Wenn Sie langsam reißen, ist die Kante gerade. Reißen Sie schnell, wird die Kante ausgefranst und gezackt. Sie haben mehr Stoff verbraucht, um diesen gezackten Rand zu erzeugen.
• Der „Qualitäts"-Faktor (Energiedichte): Dies ist die neue Entdeckung. Selbst wenn Sie diese gezackte Oberfläche glätten würden, würde es immer noch mehr Energie kosten, sie zu erzeugen, als eine ruhige, langsame Oberfläche. Die extreme Hitze an der Rissspitze verändert die atomare Struktur der neuen Oberfläche und macht sie „energiereicher" oder instabiler.
  • Vergleich: Stellen Sie sich vor, Sie backen einen Kuchen. Ein langsam gegarter Kuchen hat eine Standardtextur. Wenn Sie ihn jedoch mit einem Brenner beschossen, wird die Außenseite verkohlt und chemisch verändert. Die „verkohlte" Oberfläche ist grundlegend anders und erfordert mehr Energie zur Erzeugung als die glatte, langsam gegarte Version.

3. Die „versteckte" Rauheit

Einer der interessantesten Punkte ist, dass die vom Computer gefundene „Rauheit" so winzig ist (im Maßstab von Atomen), dass Standardwerkzeuge, die Ingenieure zur Messung von zerbrochenem Glas verwenden, sie völlig übersehen würden.

Wenn Sie ein zerbrochenes Glasstück mit einem normalen Mikroskop betrachten, würden Sie eine glatte Oberfläche sehen. Sie würden annehmen, dass die gesamte zusätzliche Energie darin investiert wurde, die Oberfläche „heißer" oder energiereicher zu machen. Doch diese Studie zeigt, dass ein erheblicher Teil dieser Energie tatsächlich darin investiert wurde, die Oberfläche physisch größer und rauer zu machen, nur in einem Maßstab, der für unser Auge zu klein ist.

4. Die Mathematik zur Geschwindigkeit von Rissen korrigieren

Die Studie korrigierte auch eine langjährige Formel, die verwendet wurde, um vorherzusagen, wie schnell ein Riss sich basierend auf der ausgeübten Kraft bewegt. Die alte Formel (das „Freund-Modell") war wie eine Karte, die bei hohen Geschwindigkeiten etwas verschwommen wurde. Die neue Studie fand eine bessere Formel (eine „Wurzel-Beziehung"), die perfekt zu den Daten passt.

Diese Korrektur ist wichtig, weil sie hilft zu erklären, warum frühere Experimente zur Messung der Hitze beim Brechen von Glas (unter Verwendung des vom Riss emittierten Lichts, genannt Fraktolumineszenz) nicht ganz mit den Geschwindigkeitsvorhersagen übereinstimmten. Durch die Verwendung der neuen Formel stimmen die vorhergesagten Geschwindigkeiten und Temperaturen endlich mit dem überein, was die Computersimulationen zeigten.

Das Fazit

Das Brechen von Glas geht nicht nur darum, Bindungen zu zerreißen. Wenn sich der Riss schnell bewegt, wirkt er wie ein winziger, überhitzter Laser, der:

1. Die Oberfläche physisch rauer macht (und mehr Fläche erzeugt).
2. Die Oberfläche chemisch verändert, um sie energiereicher zu machen.

Die Studie beweist, dass die zum Brechen von Glas erforderliche Energie keine feste Zahl ist; sie ändert sich je nachdem, wie schnell Sie brechen, und wird sowohl durch die Form des Bruchs als auch durch die extreme Hitze an der Spitze getrieben.

Technische Zusammenfassung

Problemstellung
Die Bruchenergie spröder Materialien, wie z. B. von Siliziumdioxidglas, nimmt mit der Rissgeschwindigkeit zu, ein Phänomen, das die thermodynamischen Kosten des Bindungsbruchs übersteigt. Während diese Geschwindigkeitsabhängigkeit häufig auf eine durch dynamische Instabilitäten (wie Mikroverzweigung) verursachte Oberflächenrauheit zurückgeführt wird, bleibt das Verhalten unterhalb der Verzweigungsschwelle weitgehend unverstanden. Frühere Studien an Polymeren (z. B. PMMA) verzeichneten einen Anstieg der Bruchenergie um etwa 30 % vor dem Einsetzen der Verzweigung und führten dies auf nicht aufgelöste Dissipationsprozesse zurück. Eine zentrale Frage bleibt bestehen: Entsteht dieser energieerhöhende Effekt vor der Verzweigung aus einer Zunahme der Menge der Oberfläche (durch nanoskalige Rauheit, die für Standardmessungen unsichtbar ist) oder aus einer Zunahme der Qualität (intrinsische Energiedichte) der erzeugten Oberfläche aufgrund extremer Bedingungen an der Rissspitze? Darüber hinaus haben empirische interatomare Potentiale, die in früheren Molekulardynamik-(MD-)Studien zu Siliziumdioxid verwendet wurden, hochenergetische Bindungsbrüche nicht korrekt erfasst, was zu unterschätzten Verzweigungsgeschwindigkeiten und ungenauen Beschreibungen der Prozesszone führte.

Methodik
Um diese Einschränkungen zu überwinden, wandten die Autoren umfangreiche Molekulardynamik-Simulationen aus ersten Prinzipien auf geschmolzenes Siliziumdioxidglas an.

• Interatomares Potential: Anstelle empirischer Potentiale nutzte die Studie ein maschinell erlerntes interatomares Potential (MLIP) auf Basis der Allegro-Architektur. Dieses Modell wurde auf Daten der Dichtefunktionaltheorie (DFT) unter Verwendung des r2SCAN-Funktionals trainiert und deckte einen weiten Bereich von Zuständen ab, einschließlich geschmolzenem Siliziumdioxid, glasartigen Strukturen und mechanischer Verformung bis zu 15.000 K.
• Simulationsaufbau: Das System bestand aus einem 12,3 Millionen Atome umfassenden Streifen aus Siliziumdioxidglas (300 × 8,5 × 75 nm³) mit einem vorhandenen Kerb. Die Simulationen wurden bei 300 K und 1 bar durchgeführt, gefolgt von einer schnellen Dehnungsverschiebung, um die Rissausbreitung auszulösen.
• Parameter: Vierzehn Simulationen wurden mit variierenden Verschiebungsraten durchgeführt, was zu anfänglichen Zugspannungen zwischen 3,9 GPa und 7 GPa führte.
• Analyse: Die Autoren unterschieden zwischen „struktureller Bruchenergie" (Energie, die nach dem Abkühlen in der beschädigten Struktur gespeichert ist, normiert auf die projizierte Fläche) und „intrinsischer Oberflächenenergiedichte" (Energie, normiert auf die reale, raue Oberfläche, gemessen mittels Delaunay-Tessellation). Zudem analysierten sie die Temperatur an der Rissspitze durch Mittelung der kinetischen Energie über Intervalle von 0,5 ps.

Hauptergebnisse

1. Geschwindigkeitsabhängige Bruchenergie: Selbst unterhalb der Verzweigungsschwelle (definiert bei 0,72 der Rayleigh-Wellengeschwindigkeit, $v_R$) steigt die strukturelle Bruchenergie mit zunehmender Rissgeschwindigkeit um bis zu 33 %.
2. Zerlegung des Energieanstiegs: Dieser Anstieg wird annähernd gleichmäßig auf zwei Mechanismen verteilt:
   • Nanoskalige Rauheit: Eine Zunahme der realen Bruchfläche ($A_{real}$) um etwa 16 % im Vergleich zur projizierten Fläche. Diese Rauheit erzeugt eine dickere Schadenszone um den Riss herum, tritt jedoch auf Längenskalen auf, die für eine standardmäßige optische Nachuntersuchung unsichtbar sind.
   • Intrinsische Oberflächenenergiedichte: Eine Zunahme der intrinsischen Oberflächenenergiedichte ($\gamma_s$) um etwa 15 %. Dies wird auf die extremen thermischen Bedingungen an der Rissspitze zurückgeführt, wo die Temperaturen innerhalb einer Prozesszone von nur 3–6 nm 3.000 K überschreiten.
3. Rissgeschwindigkeit und Instabilität: Die Simulationen zeigen, dass die Rissgeschwindigkeit eine Wurzelbeziehung zum Spannungsintensitätsfaktor folgt ($v_c \propto \sqrt{1 - (K_{Ic}/K_I)^2}$), was sich vom Standard-Freund-Modell unterscheidet. Die Verzweigung tritt bei etwa 2,4 km/s (0,72 $v_R$) auf, was mit experimentellen Beobachtungen an geschmolzenem Siliziumdioxid übereinstimmt.
4. Größe der Prozesszone: Der Radius der Prozesszone nimmt mit der Spannungsintensität zu und reicht von etwa 3,1 nm für ebene Risse bis zu etwa 6,8 nm für verzweigte Risse. Diese Größe stimmt mit Dugdale–Barenblatt-Vorhersagen überein und ist kleiner als Werte, die aus früheren Simulationen mit empirischen Potentialen abgeleitet wurden (~10 nm).
5. Auflösung thermischer Diskrepanzen: Die Studie löst eine Diskrepanz zwischen Simulationstemperaturen und auf Fraktolumineszenz basierenden experimentellen Vorhersagen auf. Durch Anwendung der neu identifizierten Wurzelbeziehung zwischen Geschwindigkeit und Spannung sind die vorhergesagten Temperaturen an der Rissspitze für eine gegebene Geschwindigkeit niedriger, was analytische Modelle näher an die Simulationsdaten heranführt.

Bedeutung und Behauptungen
Die Arbeit behauptet, dass die dynamische Bruchfestigkeit in Siliziumdioxidglas die Bruchenergie nicht nur durch die Schaffung einer größeren scheinbaren Oberfläche erhöht, sondern die Oberfläche auf nanoskaliger Ebene fundamental verändert. Die Autoren zeigen Folgendes:

• Die Bruchenergie ist keine konstante Materialeigenschaft; sie ist auch in Abwesenheit makroskopischer Verzweigung geschwindigkeitsabhängig.
• Eine standardmäßige Nachuntersuchung würde die Beiträge der nanoskaligen Rauheit und der erhöhten intrinsischen Oberflächenenergiedichte vermischen und den gesamten Anstieg als Zunahme der Oberflächenenergiedichte interpretieren.
• MD-Simulationen aus ersten Prinzipien unter Verwendung von MLIPs können nun eine quantitative Genauigkeit bei der Modellierung dynamischer Brüche ohne empirische Anpassung erreichen und reproduzieren erfolgreich elastische Eigenschaften, Bruchzähigkeit und Verzweigungsschwellen, die mit experimentellen Daten übereinstimmen.
• Die revidierte Geschwindigkeits-Spannungs-Beziehung bietet einen genaueren Rahmen für die Interpretation der aus Fraktolumineszenz-Experimenten abgeleiteten Temperaturen an der Rissspitze.

Die Studie kommt zu dem Schluss, dass der Anstieg der Bruchenergie vor der Verzweigung, der zuvor bei Materialien wie PMMA beobachtet, aber nicht erklärt wurde, aus einer Kombination von nanoskaliger Rauheit und erhöhter intrinsischer Oberflächenenergiedichte resultiert, die durch extreme lokale Temperaturen getrieben wird.