行波表面波驱动下石墨烯上粒子的速度锁定效应

摘要

纳米尺度上已检测到热涨落和振动诱导的快速扩散现象。本文通过分子动力学（MD）模拟和理论模型，研究行波表面波作用下粒子在石墨烯层上的运动。证实粒子在特定前提条件下将以波速持续运动，即“速度锁定效应”。通过将粒子与波动表面的范德华（vdW）势表示为曲率函数，基于相对波坐标系中的势能“水坑”阐明机制。提出两个前提条件：粒子初始位置位于势能水坑内，且初始动能不足以使其跳出水坑。参数分析显示：波长减小、振幅增大、vdW 势增强均扩大速度锁定区域。该工作为基于势能水坑理论的层状材料上粒子相干运动提供了新视角，可解释纳米尺度快速扩散现象。

引言

近年在微纳尺度检测到一系列表面波/声子诱导的快速输运和扩散现象。热迁移现象在碳纳米管 (CNT) 和石墨烯带中被广泛研究：轴向温度梯度驱动 CNT 内连续水流；大尺度石墨烯衬底温度梯度可控制石墨烯纳米薄片运动；热驱动水滴在石墨烯和 h-BN 表面输运。这些现象与特定声子模式相关：呼吸模式、屈曲行波声子等。此外，振动也能驱动粒子和液滴在 CNT 内外输运：横向声波传递动量推动纳米液滴（速度 ~30 nm/ns）；振动悬臂离心力驱动水分子定向流动；振动电荷与不对称表面能 CNT 实现单向水泵；多米诺波可将内部分子射出至 1 km/s。尽管现象丰富，但波与粒子运动的耦合机制尚不统一。主流解释为表面动量传递给粒子，但难以量化振幅、频率、相互作用势与输运的关系。有研究指出石墨烯曲率符号偏好是吸附物快速扩散的关键，暗示波动表面形貌诱导的相互作用势至关重要。因此，探究波动表面与外部粒子的相互作用对于理解纳米尺度快速输运机制具有核心意义。

方法

MD 模拟用 LAMMPS 实施。石墨烯原子面密度 ρ = 3.85×10¹⁹ m⁻²，初始平铺 (z=0)，长 6344 Å (x 向)，6000 原子；y 向周期性边界 (周期 12.2 Å)。球形粒子质量 m = 0.83×10⁻²⁵ kg，初始位置 z=7 Å, x=200 Å，初速 z 向 -50 m/s，x 向 ~2000 m/s，通过控制 z 向初速启动时间调节落点。石墨烯原子间用 REBO 势；粒子-碳原子间用 LJ 势：u(R)=ε[(σ/R)¹²-(σ/R)⁶]，ε=5.92×10⁻²¹ J, σ=4×10⁻¹⁰ m。平衡高度 h=4.2×10⁻¹⁰ m（法向力为零条件及模拟确定）。行波函数：y = A sin(2πx/λ - ωt)，A=1×10⁻⁹ m, λ=21.75×10⁻⁹ m, ω=2π/10⁻¹² s⁻¹ (周期 10 ps)，波速 v_wave=2175 m/s。左端 10 Å 石墨烯按上述振幅频率扰动触发行波，右端 x>6010 Å 固定。平坦石墨烯模拟时解除钳制并施加弱弹簧锚定 (0.0938 eV/Ų)，A=0。初始温度 5 K（抑制热激发波纹），NVE 集成演化，温度近乎不变。

结果与讨论

图 1 对比波动/平坦石墨烯上粒子轨迹（时间间隔取波周期）。波动表面上粒子相对波峰/波谷位置不变，即速度锁定于波速；平坦表面上粒子因摩擦迅速减速。不同温度、波参数下的模拟及 Xe 原子、C₆₀ 分子验证了普适性（见附件）。

理论模型（图 2）：建立波动表面 S 与粒子 P 相互作用模型。波长 λ、振幅 A，最近距离 h，面密度 ρ_s。基于 LJ 对势，推导出基于平均曲率 H 和高斯曲率 K 的相互作用势（式 3），经 MD 验证吻合。考虑摩擦（平坦石墨烯模拟估得 f = -5.2×10⁻¹³ N），构建随波移动的相对波坐标系 PXY（红色），波在该系“冻结”。粒子相对石墨烯右移，摩擦恒向左，相对势能 P = U_n + f·x（式 4），形成周期性势能水坑（红区 λ₁）和势垒（蓝区 λ₂）。

锁定条件：(1) 粒子初始位置落在势能水坑 λ₁ 内；(2) 相对波动能 ½m(v₀ - v_wave)² ≤ ΔU（水坑深度）。满足则粒子被困在水坑内振荡，绝对坐标下整体速度锁定于波速。否则粒子落入蓝区或动能过大，将跨越水坑向左滑移，最终滞后于波直至力平衡（可能静止）。文献 [21] 指出速度差会导致波散射产生推进力，故终速不为零。

运动方程（式 5-6）：法向力 F_n=∂U/∂h，切向力 F_t=∂U/∂H ∇H + ∂U/∂K ∇K。法向力零点确定 h。x 向运动方程 mẍ = F_x - f。正弦波面高斯曲率 K=0，平均曲率 H=κ，数值求解得轨迹。图 3 对比锁定/解锁案例：理论轨迹与 MD 高度吻合；锁定区粒子速度振荡于波速附近，整体等于波速；解锁区粒子逐渐停止，但波动表面仍比平坦表面增强运动。

参数效应分析（图 4，假设摩擦不变）：(a) 波长 λ 增大，水坑深度减小，超临界 λ 无锁定区；频率 ω 降低（λ 增大）亦缩小锁定区。(b) 振幅 A 增大，水坑加深，锁定效应增强；A 过小则消失。(c) 固定 λ/A 比例增大尺度，水坑深度迅速下降，表明大尺度波动下锁定效应消失。(d) LJ 参数 ε 增大，锁定区拓宽；ε 过弱则消失。参数分析与文献现象一致：涟漪振幅增大增强扩散系数（水坑加深）；涟漪谷势能低于峰（对应水坑位置）；行波纳米流道流速随振幅频率增大而增加。

潜在应用：调节振幅频率可实现近全锁定/解锁，驱动粒子定向运动或停止，垂直波动转化为横向运动类似棘轮效应，可用于纳米机电系统；不同粒子-表面相互作用导致轨迹差异，可实现相分离。

结论

本文揭示了粒子与行波石墨烯表面间独特的速度锁定现象。MD 模拟证实粒子速度可锁定于波速（需满足条件）。理论模型基于势能水坑机制阐明了锁定条件：初始位置在水坑内、相对动能不足以逃逸。理论预测轨迹与 MD 结果高度一致。波长、振幅、LJ 势参数对锁定区影响规律明确。该工作为理解波动表面上快速扩散/输运提供了新视角，并对相分离等潜在应用具有启示意义。