接触几何与面外弹性柔度调控石墨有限接触中载荷无关的结构超润滑

《CARBON》：Load-independent structural lubricity of graphite controlled by contact geometry and elastic compliance

　　

• 结构超润滑（structural lubricity）在压力跨越四个数量级范围内产生载荷无关摩擦。
• 有限接触边界通过边界贡献显著增强摩擦。
• 在有限接触中，低至中等压力下载荷无关性得以保持。
• 边界局域化石墨变形伴随载荷依赖性的出现。
• 增加石墨厚度可抑制面外变形。

该论文发表于《CARBON》，作者Hongyu Gao围绕石墨中载荷无关结构超润滑（structural lubricity）的接触几何与弹性柔度调控开展研究。背景上，传统Amontons-Coulomb摩擦近似认为摩擦力随法向载荷近似线性增加，宏观上源于微凸体接触集合响应，原子尺度上可由Prandtl-Tomlinson（PT）框架中法向载荷增强横向褶皱与滑移势垒解释。结构超润滑在有序、非公度晶体界面中因晶格失配导致横向力相消而可实现极低摩擦；有限接触中摩擦随接触面积亚线性增长，边界贡献相对增强。已有石墨/石墨烯异质结实验显示非常规载荷依赖甚至负摩擦系数，并被归因于moiré脊或边缘原子面外运动抑制及可能边缘脱附，建模也强调载荷诱导moiré褶皱展平；但结构超润滑能否使摩擦对法向载荷不敏感、有限接触中载荷无关性在何种条件下持续仍不清楚。有限石墨接触还引入接触线破坏横向平移对称和面外柔度，弹性可通过横向尺寸效应限制结构超润滑；因此研究人员在固定横向接触尺寸、改变法向压力的互补情形下，隔离边界介导面外柔度的作用。研究采用非公度Au(111)/石墨弱黏附界面作为受控参考，比较横向周期面积填充（AF，area-filling）与含显式接触线的有限接触（CL，contact-line）几何，系统改变法向载荷、滑动速度及几何和弹性参数。结论表明：AF界面因非公度性抑制相干横向力叠加，呈主要黏性响应且在四个数量级压力范围内基本载荷无关；CL引入边界贡献显著增强摩擦，低中压下载荷无关性仍保持，高压下伴随接触终止附近强非均匀面外变形而出现载荷依赖；增加石墨厚度抑制变形和高压摩擦。该研究的重要意义在于区分理想扩展界面与有限接触中的结构超润滑，提出边界几何与弹性柔度耦合的机械框架，说明载荷无关结构超润滑并非仅由非公度性决定，其鲁棒性取决于界面力相消与接触几何、弹性柔度的耦合。技术方法（≤250字）：采用分子动力学（MD，molecular dynamics）模拟，体系为刚性顶层Au(111)滑块在恒定速度下与石墨基底接触，法向压力转换为法向力并均匀施加；参考体系为三层石墨，比较横向周期AF和有限CL两种几何，压力0.1 MPa–1 GPa、速度0.1–100 m/s、时间步1 fs；Au用嵌入原子法（EAM，embedded-atom method）势、石墨用AIREBO（adaptive intermolecular reactive empirical bond order）势、Au-石墨用Morse势，并用半各向异性界面势（SAIP，semi-anisotropic interfacial potential）作稳健性检验；附加CL测试五层石墨、滑动方向30°变化、接触宽度和石墨跨度；滑动距离>200 nm，LAMMPS（Large-scale Atomic/Molecular Massively Parallel Simulator）完成。无样本队列来源。结果：3.1 面积填充与接触线几何：通过对比AF与CL的横向力空间功率谱密度（PSD，power spectral density）和力迹，AF无显著峰、滑动平滑，CL在q≈4.7 nm^?1处有峰，对应λ≈0.21 nm，接近石墨烯armchair周期（√3a/2，a=0.246 nm），力迹近周期振荡，说明CL为相干晶格周期力调制而非传统黏滑；CL摩擦显著高于AF，Morse与SAIP给出相同定性趋势。3.2 速度依赖：AF剪切应力随速度线性增加，符合Stokesian黏性响应，η_vis=0.365 kPa·s/m；CL呈次线性，τ(v)=η_visv+η_defv^c，η_vis=0.090 kPa·s/m、η_def=5.03 kPa·(s/m)^c、c=0.61，虽类似PT框架的弱于线性速度依赖，但无亚稳捕获或离散滑移证据，故为边界介导率依赖的经验描述；CL的η_vis较小是有效拟合参数，不表示黏性耦合更弱。3.3 载荷依赖和弹性柔度：AF在0.1 MPa–1 GPa剪切应力基本恒定；CL低中压近载荷无关，高压显著增加，并伴随接触终止附近石墨面外变形，三层CL的Δ_z由约0.20 ?增至1.57 ?，五层在1 GPa时Δ_z约1.12 ?、τ约149 kPa（三层约190 kPa）；层间距压缩约1.4%和1.3%较温和；经验式τ(Δ_z)=τ₀+K_zδ_z^βΘ(δ_z)拟合β≈2.2；改变滑动方向可在相似Δ_z下改变高压摩擦，增加接触宽度降低剪切应力，增加石墨跨度改变高压变形和摩擦；未观察到压力诱导横向晶格失配显著降低、塑性变形或结构转变；干净晶体界面为受控参考，实验接触还可能受污染、缺陷、边缘化学、粗糙度和多凸体影响。讨论部分总结：高压摩擦增强与接触终止附近强非均匀面外变形相关，而非均匀层间压缩或横向晶格失配大幅降低；石墨增厚同时抑制变形和高压摩擦；Δ_z是有用的弹性状态描述符但非普适摩擦控制参数；摩擦响应反映边界几何与层状基底空间非均匀弹性响应的共同作用。研究结论翻译：本工作确定接触几何与面外弹性柔度是控制有限晶体接触中载荷无关结构超润滑持续性的关键因素。在横向周期Au(111)/石墨界面中，非公度性抑制相干横向力叠加，并产生在四个数量级法向压力范围内基本不敏感的主要为黏性的响应。引入有限接触终止从根本上改变该行为：摩擦因边界贡献而增强，而载荷无关性仅在边界介导的弹性响应相对较弱时得以保持。高压下载荷依赖摩擦的出现与接触终止附近局域化的强非均匀面外变形相关，而非与均匀层间压缩或压力诱导的横向晶格失配显著降低相关。增加石墨厚度同时抑制变形和高压摩擦，而滑动方向、接触宽度和可用基底跨度的变化会改变其定量关系。因此，全局变形幅度Δ_z为弹性状态提供了有用描述符，但不是摩擦的普适控制参数。摩擦响应反映边界几何与层状基底空间非均匀弹性响应的共同影响。这些发现提示理想扩展界面与有限接触中的结构超润滑存在更广泛区别。非公度性可在界面机械响应受到足够约束时维持极弱载荷敏感性，但有限边界引入额外弹性自由度，使法向载荷可在不要求塑性变形、结构无序、静态钉扎或晶格失配显著损失的情况下改变滑动动力学。因此，载荷无关结构超润滑不应仅被视为非公度性的结果，而应被视为其鲁棒性取决于界面力相消如何与接触几何和弹性柔度耦合的状态。这为理解结构超润滑为何可在宽广加载范围内持续、但在有限石墨接触中最终变为载荷依赖提供了机械框架。