--- title: "热润滑与速率理论的边界" description: "扫得足够慢，针尖会在两个谷之间来回跳，摩擦正比于速度并趋于零。估算这在什么速度出现，再检查速率理论在哪里不再可靠。" canonical_url: "https://sandboxphysics.com/zh/stories/nanofriction/thermolubricity-and-the-limits-of-rate-theory/" language: "zh-CN" page_type: "story" alternate_languages: en: "https://sandboxphysics.com/stories/nanofriction/thermolubricity-and-the-limits-of-rate-theory/" last_modified: "2026-10-30" --- [纳米摩擦](https://sandboxphysics.com/zh/stories/nanofriction/) L3 · 第 5 章 # 热润滑与速率理论的边界 Sandbox Physics · 2026-10-30 上一章的对数律假设针尖只往前跳。如果扫描器慢到针尖在一个周期里能来回跳很多次，它就始终处在热平衡附近，摩擦变成与速度成正比，速度趋于零时消失。Vanossi 等人引用 Krylov 等人的工作，把这称为热润滑。这一章估计它从多慢的速度开始，用 Lab 里一个弱起伏针尖把它演示出来，再看前两章的速率理论在哪些条件下成立。 ## 开始前需要什么 第四章的 Kramers 速率和对数律；指数函数的数量级估算；知道 Kolmogorov 距离用来比较两个分布会有帮助，但不是必需的。 - 说出 Vanossi 等人用来区分热润滑和黏滑的无量纲比，并估算 Sang 等人的针尖在室温下需要多慢才会热润滑。 - 在弱起伏的针尖上认出热润滑：来回跳跃、平均摩擦接近零并随速度增长。 - 列出强阻尼假设、闭式近似和时间步三项检验的结果，并说明模型无法解释的实验现象。 ## 两个区域和分开它们的数 Vanossi 等人引入一个无量纲比 $\delta$：拖动的速率 $v/a$ 除以热激活跳过势垒的速率。$\delta\ll1$ 时，针尖有足够的时间在势垒两侧来回跳，摩擦按式 (4) 正比于速度；$\delta\gg1$ 时，热涨落只偶尔帮针尖提前越过势垒，就是上一章的黏滑区域。 他们写道，热润滑一般出现在极低速度或高温下，室温下典型地要求 $v<1\,\mathrm{nm\,s^{-1}}$。Sang 等人用他们的参数给出更严格的界：线性蠕变只在低于 0.1 nm/s 时出现。 判断针尖能否来回跳，要看两个谷能量相等时的势垒。这发生在扫描器位于两个格点中间，$X=a/2$。课程模型给出这个势垒 $\Delta E_{1/2}$，再用 Kramers 速率估算来回跳跃的速度尺度。 $$ \delta=\frac{v/a}{\omega_0\,e^{-\Delta E/k_BT}},\qquad F_k=\alpha(T)\,v\quad(\delta\ll1) $$ Vanossi 等人的式 (4) 及其适用条件。$\omega_0$ 是尝试频率，$\alpha(T)$ 是针尖感受到的平衡阻尼，与模型里人为选的 $\gamma$ 无关。 $$ v_{\mathrm{hop}}\approx a\,\frac{\Omega^{2}}{2\pi\gamma}\,e^{-\Delta E_{1/2}/k_BT} $$ 针尖每跳一次前进一个周期；扫描速度低于这个量级时，针尖来得及来回跳。 来源：A. Vanossi · Y. Sang ## 在 Lab 里看到热润滑 Lab 的“热润滑”预设把起伏降到 0.06 eV，以 10 nm/s 拖动。$\eta=2.55$，零温下仍会黏滑，准静态摩擦 0.046 nN，折点力 0.139 nN。但两个谷能量相等时，势垒只有 0.0385 eV，即 1.52 $k_BT$。 课程的模拟在 80 个周期里每 50 个时间步记一次针尖所在的谷：向前跨越 3288 次，向后 3208 次，平均每个周期往返约 40 次。平均摩擦是 $(0.2\pm3.5)\times10^{-4}\,\mathrm{nN}$，在统计误差内为零，比零温时小两个数量级以上。 速度提高，往返来不及完成，摩擦随之出现：100、300、1000 nm/s 时分别是 0.0044、0.0098、0.027 nN，误差约 0.001 nN。1000 nm/s 时它已经接近零温值的六成，针尖正在回到黏滑。 这里不能用 Kramers 公式求准确的跳跃速率：它要求势垒远大于 $k_BT$，而 1.52 $k_BT$ 的势垒下针尖在谷里还没达到平衡就跳出去了。公式给出约 $8\times10^{3}\,\mathrm{s^{-1}}$，只能当数量级看。 来源：A. Vanossi · Y. Sang ## 速率理论在哪里成立 阻尼。Sang 等人的式 (5) 是强阻尼形式，他们也指出它对应 $\gamma\gg\Omega$ 的极限。但他们的阻尼 $8.9\times10^{5}\,\mathrm{s^{-1}}$ 小于谷底的振动角频率 $1.15\times10^{6}\,\mathrm{s^{-1}}$。课程把阻尼加大十倍、速度降为十分之一，$v^{*}$ 保持不变：模拟摩擦 0.3437 nN，速率积分 0.3448 nN，相差 0.3%。在 Sang 的阻尼下，同一个 $v^{*}$ 的模拟是 0.362 nN，比速率积分高 5%，多出的部分来自惯性。 闭式结果。两种情形的闭式结果都是 0.358 nN，因为它只依赖 $v^{*}$。在强阻尼下它比准确的速率积分高 4%，这是只保留对数领头项的代价。上一章 1.1% 的符合，有一部分来自两项误差相互抵消。 滑移力的分布。用 Kolmogorov 距离比较：Sang 的阻尼下，30 nm/s 的 400 次滑移与式 (8) 相距 0.29，与速率积分相距 0.21，都远超 1% 显著性的临界值 0.08。阻尼加大十倍后，300 次滑移与速率积分相距 0.039，低于临界值 0.094；与式 (8) 相距 0.102，略高于临界值。Sang 等人说他们图 3 的分布符合得很好；我们没有数字化那张图，无法判断它对应的条件。 时间步。1000 nm/s 时，默认时间步给出 0.4125 nN，时间步减半给出 0.4115 nN，差别小于两者的统计误差。 适用速度的上限。Sang 等人的脚注 18 说，斜坡蠕变只在速度低于 $F_c/(M\gamma)$ 的量级时成立；用他们的参数约为每秒 70 微米。 来源：Y. Sang · A. Vanossi ## 完整例题：Sang 等人的针尖要多慢才会热润滑 题目。$U_0=0.27\,\mathrm{eV}$，$T=293\,\mathrm{K}$，其余参数同前。估算针尖开始来回跳的扫描速度，并与 Sang 等人和 Vanossi 等人给出的界比较。 第一步，对称点的势垒。扫描器位于 $X=a/2$ 时，课程模型给出 $\Delta E_{1/2}=0.433\,\mathrm{eV}=17.2\,k_BT$。 第二步，跳跃速率。指数前因子 $2.15\times10^{5}\,\mathrm{s^{-1}}$，$e^{-17.16}=3.5\times10^{-8}$，速率 $7.6\times10^{-3}\,\mathrm{s^{-1}}$。17 个 $k_BT$ 的势垒足够高，这里 Kramers 公式是可靠的。 第三步，速度尺度。$v_{\mathrm{hop}}\approx0.4\,\mathrm{nm}\times7.6\times10^{-3}\,\mathrm{s^{-1}}=3\times10^{-3}\,\mathrm{nm\,s^{-1}}$。 对照。这比 Sang 等人的界 0.1 nm/s 低约 30 倍，也远低于 Vanossi 等人说的典型值 1 nm/s；三者并不矛盾，因为后两者是上限或典型量级，不是这组参数的估算。Sang 等人模拟的速度从 5 nm/s 起，都远在这个尺度之上，所以他们的针尖处在黏滑区域。 解释。势垒每降低一个 $k_BT$，这个速度尺度就增大约 2.7 倍。把起伏降到 0.06 eV，势垒只剩 1.5 $k_BT$，同样的估算给出微米每秒的量级，所以 Lab 的 10 nm/s 已经深入热润滑区。 来源：Y. Sang · A. Vanossi ## 到实验里试一试 ### [来回跳的针尖](https://sandboxphysics.com/zh/labs/afm-stick-slip-friction/) 1. 先预测：起伏 0.06 eV、10 nm/s、室温，平均摩擦会接近零温的 0.046 nN，还是接近零？ 2. 选“热润滑”预设，看能量图里针尖在两个谷之间的往返，读出平均摩擦力。 3. 把温度拖到 0 K 再看，解释为什么同一个起伏在室温下几乎没有摩擦，而在零温下有锯齿。 ### [加大阻尼，看速率理论何时准确](https://sandboxphysics.com/zh/labs/thermal-friction-velocity-scan/) 1. 先预测：阻尼加大十倍，模拟摩擦会向速率理论曲线靠近还是远离？ 2. 打开“高级：阻尼”，把阻尼率从 $8.9\times10^{5}$ 调到约 $8.9\times10^{6}\,\mathrm{s^{-1}}$，并把扫描速度降为原来的十分之一，读出模拟摩擦和速率理论摩擦。 3. 与默认设置比较两者之差，并用惯性解释为什么阻尼大时速率理论更准。 ## 练习与解答 1. 把温度升到 373 K，Sang 等人的针尖开始来回跳的速度尺度变成多少？（$\Delta E_{1/2}=0.433\,\mathrm{eV}$，指数前因子不变。） 展开解答 $k_BT=0.0321\,\mathrm{eV}$，势垒为 $13.5\,k_BT$。速率 $2.15\times10^{5}\times1.4\times10^{-6}=0.30\,\mathrm{s^{-1}}$，速度尺度 $0.4\,\mathrm{nm}\times0.30\,\mathrm{s^{-1}}=0.12\,\mathrm{nm\,s^{-1}}$，比 293 K 时大 40 倍。 2. 为什么在起伏 0.06 eV 的针尖上不能用 Kramers 公式求准确的跳跃速率？ 展开解答 Kramers 公式假设针尖在谷里先达到热平衡，再偶尔越过高势垒，这要求势垒远大于 $k_BT$。0.0385 eV 只有 1.52 $k_BT$，针尖几乎是在两个浅谷之间自由扩散，公式只能给出数量级。 3. 闭式结果在“阻尼 $\gamma$、300 nm/s”和“阻尼 $10\gamma$、30 nm/s”两种情形下都给出 0.358 nN，而模拟分别给出 0.362 和 0.344 nN。为什么闭式结果相同？模拟的差别说明什么？ 展开解答 闭式结果只通过 $v^{*}\propto\gamma v$ 依赖阻尼和速度，两种情形的 $v^{*}$ 相同。模拟的差别来自惯性：阻尼小时针尖在谷里带着惯性运动，速率理论的强阻尼前因子随之失准。 [上一章：热激活为什么让摩擦随速度的对数增长](https://sandboxphysics.com/zh/stories/nanofriction/thermal-activation-and-the-speed-of-friction/) [下一章：一片晶体的摩擦怎样由各个原子相加而来](https://sandboxphysics.com/zh/stories/nanofriction/from-one-tip-to-a-flake/) ## 来源与复算范围 1. [A. Vanossi, N. Manini, M. Urbakh, S. Zapperi and E. Tosatti, Colloquium: Modeling friction: from nanoscale to mesoscale, arXiv:1112.3234v4; Reviews of Modern Physics 85, 529 (2013)](https://arxiv.org/abs/1112.3234v4) 2013 年 4 月 3 日 v4。第 I 节：Coulomb–Amontons 摩擦定律。第 II.A 节：Prandtl–Tomlinson 模型（式 1）、黏滑参数、滑移开始的折点、在折点附近按二分之三次方缩小的势垒、该参数大于 4.604 时的多重滑移，以及低速摩擦等于能量落差的规则。第 II.C 节：带噪声的 Langevin 方程（式 2、3）、拉动速率与跳跃速率之比、热润滑（式 4）、对数规律（式 5），以及与低温实验的分歧。第 II.D 节：Frenkel–Kontorova 链（式 6）与 Aubry 转变。第 II.E 节：超润滑、对这个名称的批评、文中转述的石墨薄片实验（多数取向低于 50 pN，接近相称取向时约 250 pN，以及薄片转回相称取向），图 4 与图 5，以及“摩擦峰的角宽度约为薄片直径的倒数”这条规则。第 III.A 节：热浴，以及没有公认取法的阻尼率。第 III.B 节：模拟能达到的时间与尺寸，以及实验与模拟中黏滑消失速度之间最多六个数量级的差距；文中印出的“每个模拟日约 1 微秒”比它自己的算式小约 9 倍。第 VI 节：结论。我们用自己的模型复算了黏滑参数、折点力和角宽度，没有数字化任何图。 2. [Y. Sang, M. Dubé and M. Grant, Thermal effects on atomic friction, arXiv:cond-mat/0109168v1; Physical Review Letters 87, 174301 (2001)](https://arxiv.org/abs/cond-mat/0109168v1) 2001 年 9 月 10 日 v1。式 (1)：斜坡蠕变：摩擦的下降量正比于温度的三分之二次方乘“速度与温度之比的对数”的三分之二次方。式 (2)、(3)：Langevin 方程，以及极小值落在晶格格点上的余弦针尖势。式 (4)：折点位置。拟合 Gnecco 等人 NaCl 数据得到的参数：晶格周期 0.4 nm，等效刚度 0.86 N/m，起伏 0.27 eV，阻尼率 8.9 乘 10 的 5 次方每秒，针尖质量 8.7 乘 10 的负 12 次方 kg；脚注 21 给出弹簧常数 0.93 N/m。正文把弹簧共振频率印成 52 MHz；用这里的质量和弹簧常数算是 52.0 kHz，所以我们把 MHz 看作笔误。式 (5)：Kramers 速率。式 (6)、(7)：折点附近的曲率与势垒。式 (8)：滑移位置分布。式 (9)：标度速度。式 (10) 至 (12)：平均摩擦的闭合式及其常数，文中给出理论值 0.55 nN，模拟与实验拟合为 0.54 nN，不确定度 0.01 nN。我们用自己的模型复算了标度速度、这些常数和平均摩擦，没有数字化图 1 至图 3。