机器学习势分子动力学模拟揭示柔性限域通道中的固液界面摩擦新机制

该研究利用机器学习势分子动力学揭示，柔性二维通道的热涨落通过改变界面静态力波动与局部剪切梯度，使固液摩擦系数呈现随通道尺度、应变和弯曲刚度可调的有限尺寸标度。

被忽略的自由度：二维通道壁并非刚性平面

二维材料广泛用于海水淡化、能量收集、生物传感和离子分离。传统纳米流体模型通常把通道壁视为固定边界，但根据 Mermin–Wagner 定理，原子级薄膜在有限温度下必然产生面外热涨落。实验片层的微米或亚微米尺度远大于典型热相干长度，长波 ZA 弯曲模因而不可忽略。

膜统计力学给出热相干长度与均方根起伏幅值的尺度估计：

其中 \(D\) 为弯曲刚度，\(E_{2D}\) 为二维杨氏模量，\(L\) 为通道横向尺寸。\(\bar h\propto L/\sqrt D\) 表明：通道越宽、壁越薄软，热起伏越强；其高度可与水分子直径 2.78 Å 相当，足以改变局部流动边界。

用 Green–Kubo 关系测量界面动量耗散

研究训练石墨烯、六方氮化硼、二硫化钼和蛭石四类通道的 NEP 势，并借助 GPUMD 在第一性原理精度下开展大尺寸、长时间分子动力学模拟。平衡态下的固液摩擦系数由线性响应理论给出：

这里 \(A=L_xL_y\) 为界面面积，\(F_{\alpha}\) 是通道壁受到的瞬时切向总力。该式把摩擦理解为界面力涨落的时间积分，即液体动量向固体耗散的统计度量。

当 \(L\approx5\) nm 时，石墨烯、h-BN、MoS₂ 和蛭石通道的 \(\lambda\) 跨越约两个数量级。更出人意料的是，柔性通道的 \(\lambda\) 随 \(L\) 增大而升高，并与 \(1/L\) 呈线性关系；刚性通道中这一规律消失。

机械应变调控 ZA 模：把柔性变成控制旋钮

面内应变可通过声子硬化或软化改变低频 ZA 模。研究对石墨烯施加 \(-0.7\%\) 至 \(1.4\%\) 的双轴预应变，发现拉伸应变抑制热起伏并降低摩擦，压缩应变则增强起伏并提高摩擦。

负斜率意味着通道尺寸增大时 \(\lambda\) 逐渐趋近更高的无限尺寸极限；随拉伸应变增加，斜率绝对值减小，说明机械加载削弱了长波涨落对摩擦的贡献。

摩擦为何改变：静态力涨落主导，而非相关时间

为分离界面摩擦的静态与动态来源，Green–Kubo 关系可重写为：

模拟显示，石墨烯的力相关时间 \(\tau_F\) 在不同尺寸和应变下基本保持在 30–50 fs；真正随通道尺度和应变显著变化的是单位面积静态力方差 \(\langle F^2\rangle/A\)。这说明柔性壁的动态曲率主要通过改变局部水—壁作用力和剪切梯度，打开新的动量耗散路径。

最终标度：弯曲刚度决定尺寸效应强弱

论文将总摩擦分解为刚性晶格粗糙度贡献 \(\lambda_0\) 与热起伏贡献 \(\lambda_{\mathrm{rip}}\)。结合连续介质流体力学与膜统计力学，可将有限尺寸效应写成相对无限尺寸极限的摩擦亏损，并得到其斜率：

图 5 的归一化结果验证了 \(|k|\propto\eta/\sqrt D\)。弯曲刚度越小，面外涨落越强，尺寸效应越突出；石墨烯的无限尺寸热起伏贡献可达刚性摩擦的约 350%，h-BN 约为 50%，而 MoS₂ 和蛭石中几乎可以忽略。

从摩擦到滑移：对纳米流动设计的意义

摩擦系数直接决定滑移长度 \(l_s\)，二者满足：

因此，若将通道简单冻结为刚性平面，可能高估滑移并低估真实能量耗散。通道尺寸、层数、弯曲刚度和预应变不再只是几何或材料参数，而成为调节纳米流量的力学自由度；未来还可通过外加应变或光—声子耦合选择性驱动壁面模态，实现主动流动控制。

该经典分子动力学框架揭示的是由几何涨落引起的界面力学贡献；涉及固体电子激发与水介电谱耦合的量子摩擦不在本文模拟范围内。