联手王中林院士！北京理工大学，最新Nature Reviews Chemistry！

DFT与MD把接触起电拆成三条路径：电子、离子与材料转移

接触起电长期被视为有害现象，可引发静电放电、火灾和器件故障，如今又被用于微纳能源收集、传感和催化，因此理解其内在机制对静电防护和高效能量转换都至关重要。已有研究指出电子、离子和带电材料碎片都可能贡献接触起电，但传统表面表征主要提供接触后的表面电位和化学态信息，难以直接捕获实时电荷转移、瞬态断键和亚纳米界面变形；实验验证和深入研究因此受制于高空间与时间分辨界面信息的获取困难。

结合实验发现与微观模拟，北京理工大学秦勇教授、中国科学院北京纳米能源与系统研究所王中林院士合作，提出接触起电可能由两种表面之间电子转移、离子转移和材料转移在不同条件下的竞争与耦合所主导。文章强调，DFT和MD为解决定量影响不明的问题提供了可能：电子转移由功函数变化或前沿轨道匹配驱动，并受应变诱导电子态重构和挠曲电效应调控；离子转移由毛细水桥形成促进，并受界面水层氢键网络破坏调控；材料转移源于机械应力。相关论文以“Advances in understanding contact electrification via DFT and molecular dynamics simulations”为题，发表在Nature Reviews Chemistry上。

综述首先用一张总览示意图把接触起电拆成离子转移、电子转移和材料转移三条路径，分别对应水桥中的离子迁移、表面间电子转移和材料碎片转移，从而确立全文以“多机制竞争与耦合”组织证据的框架（图1）。

图1 | 接触起电现象的密度泛函理论和分子动力学研究。

随后，实验机制图把三种路径并置：电子转移模型强调接触时电子转移和表面态对电荷的稳定（图2Aa、2Ab）；离子转移既可由离子聚合物接触时离子作为电荷载流子（图2Ba、2Bb），也可由OH⁻和H₃O⁺沿水桥在局域电场与浓度梯度驱动下迁移（图2Bc）；材料转移则用碎片转移（图2Ca）、马赛克电荷图案（图2Cb）和表面电荷密度与杨氏模量的反相关（图2Cc）说明软聚合物表面电荷密度更高，正文进一步将其归因于软聚合物更易机械诱导断裂并转移带电碎片。由此得出的结论是，摩擦电序列只能作为指南，三种机制均已在实验中出现，但定量贡献仍难分离。

图2 | 实验观察到的三种电荷转移机制

实验机制图之后，综述用双电层两步模型解释固-液界面电子转移与离子转移的共存：液体接触原始固体表面后，水分子和离子在热运动与液体压力驱动下随机碰撞表面，此过程中可发生电子转移；固体表面通过电子转移和/或表面电离反应带电，因此表面结合离子位点应与液体中可移动离子区分；随后液体中抗衡离子受静电吸引向带电表面聚集，形成电双层，紧邻带电表面的紧凑抗衡离子层为Stern层（图3a–d）。正文进一步指出，该模型由超快电子转移（<10⁻¹² s）建立初始界面电位，随后抗衡离子组织成结构化双电层（约10⁻⁶–10⁻³ s），从而稳定表面电荷，并将短时间尺度电子动力学与实验可及的离子响应联系起来。

图3 | 电双层两步物理模型。

在电子转移路径上，功函数与前沿轨道图用能级图（图4A）、金属功函数与界面能级随分子相互作用的演化（图4Ba）、MXene表面偶极与功函数关系（图4Bb）、火山型关系（图4Bc）以及PTFE不同终止结构的LUMO分布（图4Db），说明DFT可以把经典功函数差和聚合物前沿轨道能级匹配统一到界面电子转移框架中；文章据此指出，无机晶体可用功函数预测电荷转移倾向，聚合物则需辅以前沿轨道能级与空间分布（图4C、4Db），且界面轨道隙约低于3 eV时纯电子图像最合理，超过该阈值则机械化学断键主导的材料转移更有利；PVA中引入MOF-303可提高HOMO、窄化供体—受体隙并增强电子转移（图4Da）。

图4 | 功函数与前沿轨道。

在应变调控方面，应变多尺度图把一维PTFE链（图5Aa、5Ab）、二维波纹2H-MoTe₂接触（图5Ba、5Bb）和三维PTFE-Al界面（图5Ca、5Cb）并置，展示应变如何调制前沿轨道空间分布、电荷密度差和隧穿势垒；其结论是，应变不仅改变带隙、功函数和载流子有效质量，还可通过挠曲电极化和非线性势垒驱动不可逆电子隧穿，从而把机械变形直接写入电子转移路径。

图5 | 应变对原子和电子结构影响的多尺度分析。

在非电子路径上，离子与材料转移机制图用热电子发射区分电子与“粘性”离子（图6Aa）、用空气—水界面离子分布（图6Ab）和MD水桥模型计算H₃O⁺/OH⁻在两聚合物表面的分配自由能（图6Ac），并用能量判据（图6Ba）、内聚能密度矩阵（图6Bb）和FEM-RMD（图6Bc）展示断键与碎片转移；文章据此认为，界面氢键网络破坏是离子载流子的关键来源，而共价断键驱动的材料转移是聚合物接触起电的普遍机制，相同聚合物接触时表面形貌造成的拉伸/压缩不对称可解释光滑面倾向负电、粗糙面带正电。

图6 | 离子转移和材料转移机制。

总结图进一步以原子间相互作用势串起接触—分离循环：长程吸引区由水桥介导离子迁移（图7Ba），平衡接触区由功函数差和前沿轨道匹配促成电子转移（图7Bb），短程压缩区由电子云重叠、隧穿和挠曲电局域场增强电荷输运（图7Bc），分离区则留下局域电子、断键碎片和再吸附水分子造成的持续电荷不平衡（图7Bd–f）；接触起电是高度耦合的多尺度非平衡界面现象，电子、离子和材料转移在不同接触—分离区间各自占主导（图7A、7B）。

图7 | 固-固接触-分离循环中三种主要电荷转移路径及其相关物理化学机制的时空演化和可能排序。

结论与展望：接触起电的三种机制在不同接触—分离区间各自占主导，但常产生难以区分的宏观信号；当前挑战包括势相互作用剖面具有材料特异性、表面原子尺度粗糙导致位点差异、瞬态效应量化困难，以及CE本质上的多尺度非平衡特征。未来需要把原位表征与多尺度模拟紧密结合，并谨慎评估DFT泛函误差、经典MD对断键和极化的缺失、AIMD尺度限制以及机器学习力场对训练数据覆盖度和标注精度的依赖，以建立能连接原子尺度电荷转移与宏观电荷弛豫的可靠框架。