NAVIS：一种用于高效计算纳米通道速度滑移与热界面滑移的LAMMPS-Python框架

《Computer Physics Communications》：NAVIS: A LAMMPS-Python framework for efficient computation of nanochannel velocity and thermal interfacial slip

　　

研究人员提出NAVIS（NAnochannel Velocity and thermal Interfacial Slip，纳米通道速度与热界面滑移），一种基于LAMMPS-Python脚本的工具包，用于计算任意固-液界面的Navier（流体力学）摩擦系数与Kapitza（热）阻。NAVIS基于平衡分子动力学（EMD）方法计算界面的线性响应摩擦与热阻，以及相应的速度与温度滑移。该方法依托前期研究（Hansen等人，Phys. Rev. E 84, 016313 (2011)；Varghese等人，J. Chem. Phys. 154, 184707 (2021)；Alosious等人，J. Chem. Phys. 151, 194502 (2019)；Alosious等人，Langmuir 37, 2355–2361 (2021)），并在本工作中提供教学式框架，在两类体系上实现：水-石墨烯体系（流体力学滑移）与水-碳纳米管（CNT）体系（热滑移）。研究人员给出利用LAMMPS包执行EMD模拟及利用Python模块处理输出以获取目标量的详细流程。该工具包有望服务于研究界面摩擦与热输运的计算学者，这两者是纳米流体系统高效实用化的关键因素。

研究背景方面，固-液界面的摩擦与传热近年成为研究焦点，其关乎减摩润滑与换热效率提升，对缓解气候变化下的能耗压力具有重要意义。当前计算界面输运性质主要分为平衡分子动力学（EMD）与非平衡分子动力学（NEMD）两类；NEMD在实验相关外场下难以获得统计可靠的流场，需外推至实验可达场强，且对高滑移体系偏差明显，因此基于线性响应理论的EMD更受青睐。但经典Green-Kubo类方法在受限（非零波矢）体系中因关联函数长时发散导致积分消失，给出偏低的摩擦系数与偏高的滑移长度；Bocquet-Barrat方法只能得到两壁混合的有效摩擦，且把流体本体黏性阻力混入界面本征摩擦。针对这些争议，Hansen等人提出不依赖体积积分的EMD替代方案，Varghese等人其后改进统计精度；热滑移侧类似地，Barrat-Chiaruttini的Green-Kubo类积分也长时发散，Alosious等人提出非积分时间关联法更稳定。在此背景下，研究人员开展NAVIS框架研究，旨在将前述理论以可复现的教学式代码落地，在LAMMPS与Python间搭建流水线，得出可分别计算Navier摩擦系数与Kapitza阻的实用结论，论文发表于《Computer Physics Communications》，对纳米流体界面输运计算具有方法学推广意义。

关键技术方法上，研究人员采用三类样本体系：水-石墨烯狭缝通道（测流体力学滑移）、水-碳纳米管（CNT）体系（测热滑移），均取x、y方向周期、z向受限构型。计算内核沿用Hansen-Varghese改进的非体积积分EMD方案提取界面本征Navier摩擦系数ξ₀，即由近壁流体薄层速度u_x(t)与壁-流切向力F_x(t)构造关联C_uxFx(t)=?u_x(0)F_x(t)?与C_uxux(t)=?u_x(0)u_x(t)?，经Laplace域关系ξ₀=-k_BT/lim_s→0 C?_uxFx(s)/C?_uxux(s)求取；热侧沿用Alosious等人的非积分时间关联法算Kapitza阻R_K，避免长时发散。模拟端用LAMMPS跑EMD产出轨迹与力/速度时间序列，后端用Python模块做分块平均、关联估计与Laplace变换数值反演，整套流程在GitHub以GPL-v3发布。

研究结果部分，按论文主体小节归纳如下。1. Introduction一节指出23%全球能耗源于摩擦接触，其中五分之一用于克服摩擦，凸显固-液界面知识缺口；并列明EMD相对NEMD在线性区的优势及Green-Kubo发散根源。2. Theory小节给出ξ₀与R_K的局部界面定义：二者皆由紧邻壁面的流体层与壁原子相互作用主导，通过本构式定义的相关量时间关联函数在EMD中抽取。2.1 Navier Friction coefficient给出u_x(t)、F_x(t)及两个关联函数定义，并引入Laplace变换L{f(t)}=∫₀^∞f(t)e^-stdt=??(s)，导出零频摩擦系数公式。2.2 Kapitza resistance类比给出热流-温度关联的非积分稳定表达式。3. Simulation procedure描述水-石墨烯与水-CNT的盒型、力场（SPC/E水、AIREBO或TIP3P配CHARMM类壁）、控温控压与EMD采样长度设置。4. Post-processing with Python说明分窗求速度剖面、力自相关、块平均降噪与Laplace域比值外推s→0的具体步骤。5. Results and discussion中，水-石墨烯算出ξ₀对应滑移长度与文献NEMD外推值吻合；水-CNT得R_K随管径与手性变化趋势，验证Alosious法稳定性。结论部分翻译为：研究人员证实NAVIS能以纯EMD流程、免长时积分截断争议，分别给出双壁各异的本征Navier摩擦与Kapitza阻；框架对水-石墨烯与水-CNT的示范表明，其输出与既有外推/NEMD结果一致，可作为纳米流体界面速度滑移与热滑移的标准化计算入口。

讨论总结方面，原文强调NAVIS的价值不在提出新物理，而在把Hansen-Varghese与Alosious等分散方法学收束为带示例数据的可复现LAMMPS-Python教具，回避了Bocquet-Barrat体积积分混入本体黏性、Green-Kubo长时发散须人为截断两处痛点；同时单轮EMD即可分辩两不同壁的本征系数，优于只能出混合有效值的旧法。局限上原文未声称覆盖非线性和大驱动场，也未替代TTCF-NEMD在强非平衡场景的角色，定位仍是线性响应下的界面本征输运抽取器。