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

利用机器学习势函数的分子动力学模拟，本研究表明，在分叉阈值以下，二氧化硅玻璃的断裂能因本征表面能密度升高与纳米级粗糙化共同作用而最高增加 33%，这证明动态断裂形成了截然不同的表面结构，而不仅仅是增加了表观表面积。

要点

想象你正在试图折断一块厚玻璃。你可能会认为，将其折断所需的能量仅仅是切断将其维系在一起的微小原子键所需的能量，就像剪断一根意大利面一样。然而，科学家们早就知道，打破玻璃实际上所需的能量比这种简单计算所暗示的要多得多。这就像玻璃在反抗，需要额外的努力才能将其打破。

多年来，研究人员一直认为这种“额外成本”主要是由于裂纹在加速过程中变得不稳定且参差不齐，从而产生了更粗糙的表面积（就像将一张纸撕成锯齿状的条带，而不是直线）。但一项利用先进计算机模拟的新研究揭示了一个更为复杂的故事。

以下是该论文发现的简要说明：

1. “超热”的裂纹尖端

当裂纹在玻璃中高速移动时，其尖端会变得极其炽热。研究发现，在高速下，裂纹尖端的原子温度可达8,000 开尔文（比太阳表面还热！）。

请将裂纹尖端不仅仅视为一个断裂点，而将其视为一个微小的、微观的喷灯。这种 intense 的热量不仅熔化了玻璃，而且从根本上改变了所产生表面的性质。

2. 打破玻璃成本更高的两个原因

研究人员使用了一个超精确的计算机模型（就像能看见单个原子的数字显微镜）来弄清楚所有这些额外能量都去了哪里。他们发现，“额外成本”大致50/50地分配在以下两个方面：

• “粗糙度”因素（数量）： 随着裂纹加速，其留下的表面并非完全平滑。它在纳米尺度上变得粗糙，就像从太空看到的山脉。这意味着裂纹实际上创造了比从外部看起来更多的表面积。
  • 类比： 想象撕开一块布料。如果你慢慢撕，边缘是直的。如果你快速撕，边缘会变得毛糙且参差不齐。为了制造这种锯齿状边缘，你实际上使用了更多的布料。
• “质量”因素（能量密度）： 这是新发现。即使你抚平了那个锯齿状的表面，创造它所需的能量仍然比创造平静、缓慢的表面要多。裂纹尖端的极端热量改变了新表面的原子结构，使其具有“更高能量”或更不稳定。
  • 类比： 想象烤蛋糕。慢火烤制的蛋糕具有标准的质地。但如果你用喷灯猛烤，外部会变得焦黑并在化学上发生变化。“焦黑”的表面在本质上不同，并且比光滑、慢火烤制的版本需要更多的能量来创造。

3. “隐藏”的粗糙度

最有趣的观点之一是，计算机发现的这种“粗糙度”极其微小（处于原子尺度），以至于工程师用来测量破碎玻璃的标准工具会完全忽略它。

如果你用普通显微镜观察一块破碎的玻璃，你会看到一个光滑的表面。你会认为所有额外的能量都用于使表面“更热”或更具能量。但这项研究表明，其中相当一部分能量实际上用于使表面在物理上变得更大、更粗糙，只是这个尺度太小，我们的眼睛无法看见。

4. 修正裂纹移动速度的数学公式

该论文还修正了一个长期使用的公式，该公式用于根据施加的力来预测裂纹的移动速度。旧公式（“弗劳德模型”）就像一张在高速下变得有点模糊的地图。这项新研究发现了一个更好的公式（“平方根关系”），它能完美地拟合数据。

这一修正很重要，因为它有助于解释为什么之前测量破碎玻璃热量（利用裂纹发出的光，称为断裂发光）的实验与速度预测不太吻合。通过使用新公式，预测的速度和温度终于与计算机模拟显示的结果一致了。

核心结论

打破玻璃不仅仅是关于断开原子键。当裂纹快速移动时，它就像一个微小的、超热的激光，会：

1. 使表面在物理上变得更粗糙（创造更大的面积）。
2. 在化学上改变表面，使其更具能量。

这项研究证明，打破玻璃所需的能量不是一个固定数值；它会随着你打破它的速度而变化，并且是由断裂的形状以及尖端的极端热量共同驱动的。

技术摘要

问题陈述
脆性材料（如二氧化硅玻璃）的断裂能随裂纹速度增加而增大，这一现象超出了化学键断裂的热力学成本。虽然这种速度依赖性通常归因于动态不稳定性（如微分支）引起的表面粗糙化，但分支阈值以下的行为仍知之甚少。以往对聚合物（如 PMMA）的研究指出，在分支发生前断裂能增加了约 30%，并将其归因于未解析的耗散过程。一个核心问题依然存在：这种分支前的能量增加是源于表面积数量的增加（通过标准测量无法观测到的纳米级粗糙化），还是源于极端裂纹尖端条件下所创生表面的质量（本征能量密度）的提升？此外，以往二氧化硅分子动力学（MD）研究中使用的经验原子间势未能准确捕捉高能断键事件，导致分支速度被低估且过程区描述不准确。

方法论
为了解决这些局限性，作者对熔融二氧化硅玻璃进行了广泛的从头算分子动力学模拟。

• 原子间势：该研究未采用经验势，而是利用基于 Allegro 架构的机器学习原子间势（MLIP）。该模型使用 r2SCAN 泛函在密度泛函理论（DFT）数据上进行训练，涵盖了包括熔融二氧化硅、玻璃态结构以及高达 15,000 K 的机械变形在内的广泛状态。
• 模拟设置：系统由一个含有预制缺口的 1230 万个原子的二氧化硅玻璃条（300 × 8.5 × 75 nm³）组成。模拟在 300 K 和 1 bar 条件下进行，随后施加快速拉伸位移以诱发裂纹扩展。
• 参数：进行了十四次模拟，位移速率各不相同，产生的初始拉伸应力介于 3.9 GPa 至 7 GPa 之间。
• 分析：作者区分了“结构断裂能”（冷却后受损结构中储存的能量，按投影面积归一化）和“本征表面能密度”（按通过 Delaunay 剖分测量的真实粗糙表面积归一化的能量）。此外，他们还通过对 0.5 ps 时间间隔内的动能取平均值来分析裂纹尖端温度。

关键结果

1. 速度依赖的断裂能：即使在分支阈值以下（定义为瑞利波速 $v_R$ 的 0.72），随着裂纹速度的增加，结构断裂能也上升了高达 33%。
2. 能量增加的分解：这种上升大致均等地分配给两种机制：
   • 纳米级粗糙化：真实断裂表面积（$A_{real}$）比投影面积增加了约 16%。这种粗糙化在裂纹周围形成了一个更厚的损伤区，但其发生的尺度在标准的死后光学分析中是不可见的。
   • 本征表面能密度：本征表面能密度（$\gamma_s$）增加了约 15%。这归因于裂纹尖端的极端热条件，在仅 3–6 nm 的过程区内，温度超过了 3,000 K。
3. 裂纹速度与不稳定性：模拟显示，裂纹速度与应力强度因子遵循平方根关系（$v_c \propto \sqrt{1 - (K_{Ic}/K_I)^2}$），这与标准的 Freund 模型不同。分支发生在约 2.4 km/s（0.72 $v_R$），与熔融二氧化硅的实验观测一致。
4. 过程区尺寸：过程区半径随应力强度增加而增大，范围从平面裂纹的约 3.1 nm 到分支裂纹的约 6.8 nm。该尺寸与 Dugdale–Barenblatt 预测一致，且小于以往经验势模拟得出的数值（约 10 nm）。
5. 热差异的解决：该研究解决了模拟温度与基于断裂发光实验预测之间的差异。通过应用新发现的平方根速度 - 应力关系，给定速度下预测的裂纹尖端温度更低，从而使分析模型与模拟数据更加吻合。

意义与主张
该论文主张，二氧化硅玻璃的动态断裂不仅通过产生更多的表观表面积来增加断裂能，而且通过在纳米尺度上根本性地改变表面来实现。作者证明：

• 断裂能并非恒定的材料属性；即使在缺乏宏观分支的情况下，它也具有速度依赖性。
• 标准的死后分析会将纳米级粗糙化的贡献与增加的本征能量密度混为一谈，将整个增加量解释为表面能密度的上升。
• 使用 MLIP 的从头算 MD 模拟现在可以在无需经验调整的情况下实现动态断裂建模的定量精度，成功复现了与实验数据相匹配的弹性性能、断裂韧性和分支阈值。
• 修正后的速度 - 应力关系为解释从断裂发光实验中推断出的裂纹尖端温度提供了更准确的框架。

该研究得出结论，此前在 PMMA 等材料中观察到但无法解释的分支前断裂能增加，源于纳米级粗糙化与由极端局部温度驱动的本征表面能密度升高的共同作用。