Water cavitation results from the kinetic competition of bulk, surface and surface-defect nucleation events

本文提出了一种经分子动力学模拟验证的动力学模型，该模型通过阐明体相、表面及表面缺陷成核路径之间的竞争（尤其是纳米尺度疏水缺陷的主导作用）如何决定水在负压下的亚稳态，从而解释了实验观测到的水空化压力存在广泛差异的原因。

要点

想象一下处于负压状态的水，就像一根被紧紧拉伸的橡皮筋。它渴望弹回松弛的状态，却拼命维持着紧绷。最终，它放弃抵抗，通过在水中内部形成微小的气泡而“崩断”。这种崩断事件被称为空化。

长期以来，科学家们一直困惑于为何这种崩断会在不同强度（压力）下发生，具体取决于情境。有时，水能坚持到承受巨大应力（极端的负压）时才崩断，而有时它几乎会立即崩断。

本文如同一篇侦探故事，揭示了水决定在何处以及如何崩断的谜团。作者构建了一个计算机模型来模拟箱体内的水，发现水实际上有三种不同的断裂方式，并且它们彼此之间不断竞争：

1. “房间中央”断裂（体相空化）

想象一个完美洁净、空旷的房间，墙壁光滑且湿润。如果你将橡皮筋（水）拉得足够紧，它最终会在房间正中央、远离任何墙壁的地方崩断。

• 结果：这需要极端的应力。水必须被拉伸至约**-100 MPa**（巨大的负压）时，才会在中央崩断。这是最“纯粹”的断裂形式，但极难实现，因为现实中的水很少是绝对纯净的。

2. “墙壁”断裂（表面空化）

现在，想象房间的墙壁并非完美湿润，而是有点“油腻”或具有排斥性（疏水性）。水不喜欢接触这些墙壁。

• 类比：想象水试图拥抱一面它不喜欢的墙壁。如果墙壁过于“排斥”（具体而言，如果接触角大于50°至 60°），水就会放弃依附墙壁，直接在表面形成气泡，而不是等待在中央崩断。
• 结果：这种情况发生起来容易得多。水在低得多的应力水平下就会崩断，大约在**-30 MPa**。墙壁的“粘性”决定了是否会发生这种情况。如果墙壁亲水性很强（易湿润），水就会保持原位；如果墙壁具有排斥性，气泡就会提前形成。

3. “隐藏陷阱”断裂（缺陷空化）

这是最戏剧性的场景。想象墙壁上有一道微小的划痕、一个凹坑，或是一点超级油腻的灰尘（一个“纳米级缺陷”）。

• 类比：将这个缺陷想象成一个预先制作好的活板门。即使房间其余部分是完美湿润的表面，这个微小的油腻凹坑也会像磁铁一样吸引气泡。它如此有效，以至于即使水只受到轻微应力，气泡也能在那里几乎瞬间形成。
• 结果：单个微小的缺陷（小至几纳米）就能主导整个过程。它显著提高了“断裂点”，意味着水会在比完美系统中高得多的压力下崩断（更接近零甚至为正压）。

全局视角：这为何重要？

本文解释了为何实验会显示出如此多样的结果。

• 如果你拥有超纯水并置于完美光滑且湿润的容器中，它将坚持直到达到极端的**-100 MPa**极限（体相空化）。
• 如果你拥有普通水且表面略带油腻，它会在早得多的时候崩断，大约在**-30 MPa**（表面空化）。
• 如果你拥有脏水或表面带有微小划痕/凹坑，它几乎会立即崩断（缺陷空化）。

核心结论：
作者创建了一个“规则手册”（动力学模型），将这三种情景结合起来。他们发现，这场竞争的“赢家”取决于两个主要因素：

1. 表面的排斥程度：如果表面过于“油腻”（接触角 > 60°），气泡就会在表面形成。
2. 微小陷阱的存在：即使只有一个微小的缺陷，也能劫持整个过程，使水在远早于完美水物理预测的时间点崩断。

简而言之，水并非随机崩断；它会在可用的“最薄弱环节”崩断，无论是液体中央、墙壁，还是墙壁上的微小划痕。这解释了为何自然界和工程系统在受压的水中会观察到如此不同的行为。

技术摘要

技术摘要：水空化结果源于体相、表面及表面缺陷成核事件的动力学竞争

问题陈述
处于负压状态的水存在于亚稳态，直到通过空化（即蒸汽气泡的成核）达到平衡。实验上，空化阈值压力（$p_{cav}$）表现出巨大的变异性，具体取决于水的纯度、表面接触角以及表面质量。虽然经典成核理论（CNT）预测洁净体相水能够承受远低于 $-100$ MPa 的负压，但实验观察通常显示空化发生在高得多的压力下（例如 $-30$ MPa，甚至在标准容器中为正压）。这种差异表明空化路径存在竞争：体相中的均匀成核与表面及表面缺陷处的非均匀成核。此前缺乏一个统一的动力学框架来描述这三种路径（体相、表面和缺陷）之间的竞争，阻碍了对合成系统和生物系统中空化现象的理解。

方法论
作者通过构建一个将经典成核理论（CNT）与原子分子动力学（MD）模拟相结合的动力学模型，填补了这一空白。

• 模拟设置：研究采用由两层自组装单分子层（SAMs）组成的平面表面模型，该单分子层由羟基化烷烃构成。通过缩放头基的部分电荷来调节表面极性，进而调节水接触角（$\theta$）。系统用水（SPC/E 模型）溶剂化并施加负压。
• 动力学协议：为了确定成核尝试频率（$\kappa$），作者采用压力斜坡模拟协议，其中负压随时间线性降低（$p(t) = \dot{p}t$）。当模拟盒长度增加 50% 时，判定为空化发生。
• 理论框架：总空化速率（$k_{tot}$）被建模为体相（$k_{3D}$）、光滑表面（$k_{2D}$）和表面缺陷（$k_{def}$）独立速率之和。该模型假设这些路径互不干扰。作者推导了各路径的自由能势垒（$G^*$）和尝试频率的表达式，并利用 MD 数据对难以通过解析方法获得的尝试频率（$\kappa_{3D}$ 和 $\kappa_{2D}$）进行参数化。

主要贡献与结果

1. 尝试频率的量化：作者通过 MD 模拟确定了体相尝试频率密度（$\kappa_{3D}$）和表面尝试频率密度（$\kappa_{2D}$）随接触角 $\theta$ 的变化关系。研究发现 $\kappa_{2D}$ 随 $\theta$ 变化，在约 $75^\circ$ 处出现最大值，但这种变化相对于自由能势垒的指数依赖性而言是次要的。
2. 动力学竞争与交叉：模型揭示了体相空化与表面空化之间存在急剧的动力学转变。
   • 无缺陷表面：对于接触角低于临界阈值 $\theta^* \approx 50^\circ - 60^\circ$ 的亲水表面，体相空化占主导地位，其 $p_{cav}$ 与 CNT 预测一致（约 $-100$ MPa）。
   • 疏水表面：当 $\theta$ 超过 $\theta^*$ 时，表面空化成为主导路径。对于强疏水表面，典型的空化压力显著升高至约 $-30$ MPa。
   • 系统尺寸依赖性：交叉角 $\theta^*$ 仅微弱依赖于系统尺寸和观测时间，变化范围狭窄。
3. 表面缺陷的作用：研究表明，纳米级疏水表面缺陷充当高效的空化成核位点。
   • 即使单个纳米级缺陷也能主导宏观系统中的空化动力学。
   • 承载预存蒸汽气泡的缺陷可将临界空化压力提升至远高于光滑表面的水平。
   • 该模型提供了一个统一的框架，用于根据缺陷尺寸（$R_{def}$）、缺陷接触角（$\theta_{def}$）和缺陷数量（$N_{def}$）预测 $p_{cav}$。
4. 验证：包含单个疏水凹坑（缺陷）系统的显式 MD 模拟证实了理论预测。与无缺陷表面（$-95$ MPa）相比，半径为 1 nm 和 2 nm 的缺陷系统分别在显著更高的压力下（分别为 $-66$ MPa 和 $-45$ MPa）发生空化，在定性上与模型相符。

意义
本文声称其统一的动力学模型解释了不同系统中实验观测到的空化压力的巨大变异性。通过量化体相、表面和缺陷成核之间的竞争，研究结果强调：

• “洁净体相”极限（$p_{cav} \approx -100$ MPa）仅能在无缺陷且高度亲水的环境中实现（例如微观石英空腔）。
• 在大多数实际场景中，空化由表面性质和缺陷主导，而非体相热力学。
• 该模型为理解从喷墨打印和超声波清洗到生物现象（如树木中树液的上升和鼓虾的力学机制）等应用中的空化现象提供了预测工具。
• 研究结果表明，表面吸附（例如两亲性分子）可以通过改变接触角或覆盖缺陷来抑制空化，这为控制空化稳定性提供了潜在的机制。