物理学快报(Physical Review Letters)发表研究揭示活性物质界面普适类的新机制
主要作者及发表信息
本研究由巴黎-萨克雷大学(Université Paris-Saclay)、法国国家科学研究中心(CNRS)、固体物理实验室(Laboratoire de Physique des Solides)的Raphaël Maire(通讯作者)、Andrea Plati、Frank Smallenburg及Giuseppe Foffi(通讯作者)共同完成。该论文于2026年4月6日在线发表于《物理评论快报》(Physical Review Letters)第136卷,并被编辑推荐为亮点文章。
研究背景与动机
在热力学平衡态下,不同体系中相分离界面的统计行为表现出极强的普适性。无论是胶体悬浮液、分子流体还是二维伊辛模型(Ising model),其界面在大尺度上均受毛细波哈密顿量(capillary-wave Hamiltonian)支配,静态特性与动力学特性存在解耦。然而,在远离平衡态的活性物质系统中,这种普适性被认为会失效。理论预测指出,根据主体相动力学中守恒量的不同,非生长型界面应存在三种不同的非平衡普适类:JQKPZ普适类(jqkpz universality class)、湿JQKPZ普适类(wet-jqkpz universality class)以及超均匀界面普适类(hyperuniform interface universality class)。这三类在一维界面中的静态粗糙度指数(roughness exponent)χ和动力学指数(dynamical exponent)z均不相同,打破了平衡态的普适规律。
然而,尽管理论预言丰富,实验和数值模拟中观察到的活性物质界面波动却长期表现出惊人的类平衡特性,其静态高度关联函数仍近似满足S(k) ∼ |k|⁻²。JQKPZ行为仅在极少数强受限系统中被报道过,而湿JQKPZ行为此前从未被观测到。这一理论与实验的鸿沟构成了该领域的核心难题。本研究旨在构建一个能够清晰展现全部三种理论预言的普适类的微观模型,并探索此前被忽视的、由主体相慢动力学引发的新界面行为。
研究流程与模型设计
为实现上述目标,研究团队设计了一个受振动颗粒系统启发的“活性碰撞驱动硬碟模型”(hard-disk model driven by active collisions)。模型包含N个质量为m、直径为σ的二维硬碟,置于周期性边界条件的长方形盒子中。每个粒子的运动遵循朗之万方程(Langevin equation),其中阻尼项γ和温度项T决定了系统是否具备动量守恒。该模型的核心创新在于引入了非平衡碰撞机制:当两个粒子接触时,它们的动能并非简单耗散或弹性反弹,而是根据一套内置的“充能”规则进行改变。
具体而言,碰撞能量变化δE_ij由两部分组成:一部分是与碰撞相对速度相关的耗散项,由恢复系数α(coefficient of restitution)控制;另一部分是能量注入项δE⁺_ij,其大小取决于每个粒子自上次碰撞以来所经过的时间τ_i和τ_j。注入能量在最小值2δe0和最大值2(δe0 + δe)之间,通过一个特征充电时间τr和陡度β控制的插值函数进行调节。这种设计巧妙地将准二维振动颗粒单层的物理过程——垂直振动板的能量间歇性地通过碰撞转化为水平运动——映射为了粒子内部状态的更新。在这种粗粒化描述下,活性并非源于持续的自我推进,而是通过碰撞介导的能量转移涌现出来,这与自推进粒子模型(如活性布朗粒子)有本质区别。研究采用分子动力学模拟,结合时间步进和事件驱动算法,在细长的盒子中初始化一个阶梯密度分布,并追踪界面高度h(y, t)随时间的演化。
主要结果与数据支撑
通过精确调控模型参数,研究团队首次在同一个模型框架内,清晰地观测到了全部三种理论预言的界面普适类。
(一)JQKPZ普适类:当同时存在有限温度(T ≠ 0)和有限阻尼(γ ≠ 0)时,朗之万热浴耗散了动量,使得密度成为唯一的守恒场。模拟结果显示,稠密相温度显著低于稀疏相,这是强烈的非平衡特征。稳态高度关联函数呈现S_h(k) ∼ k⁻¹⋅⁵的标度,对应粗糙度指数χ ≃ 0.25,显著低于平衡态的χ_eq = 1/2。在粗化阶段,界面宽度平方w²(t) ∼ t⁰⋅¹⁸增长,与JQKPZ类理论预测的慢动力学一致。结合χ值,导出动力学指数z ≃ 2.78,通过有限尺寸标度分析得到的稳态宽度w²_ss和稳态建立时间τ_ss与系统尺寸L的关系也进一步证实了这些指数。
(二)湿JQKPZ普适类:当阻尼γ = 0时,由于缺少朗之万热浴,动量成为局域守恒量。此条件下,界面显著更粗糙。S_h(k) ∼ k⁻³和w²_ss ∼ L²均给出χ ≃ 1,与理论预期χ_wet = 1完美吻合。然而,其动力学表现出异常粗化行为,界面宽度未呈现清晰的幂律增长,这与直接模拟湿JQKPZ场方程时报告的异常粗化现象一致。
(三)超均匀界面普适类:当温度T = 0而阻尼γ ≠ 0,即动量被阻尼耗散但碰撞过程仍守恒质心系总动量时,界面变得极其平滑。S_h(k) ∼ k⁻¹,证实了理论预言的χ_hu = 0。由于界面宽度仅在粗化阶段呈对数增长,无法从早期动力学提取z,但通过测量弛豫时间,得到z ≃ 3,与理论值完全一致。
在完成液-气界面的研究后,团队进一步探索了主体相慢动力学对界面标度的影响。首先,对于液-固界面,当稠密相堆积分数达到约0.7以上时,系统出现六角序,可通过键取向序参量ψ₆定量表征。无论是动量守恒还是非动量守恒系统,一旦稠密相发生结构化,静态高度关联的标度指数χ均趋向于0,表明界面变得超均匀光滑,偏离了其液-气状态下的JQKPZ或湿JQKPZ行为。其次,对于液-玻璃界面,通过使用大小粒子二元混合物抑制结晶,同样发现当稠密相发生动力学停滞(表现为均方位移出现平台)时,χ从湿JQKPZ的1单调递减至约0.25。这两项发现共同揭示了一个全新的、此前被忽视的机制:无论具体来源是结构序还是动力学停滞,主体相中慢动力学的出现会普遍性地压平界面,导致其静态起伏大幅降低。
结论与价值
这项研究的贡献具有双重性。其一,它通过一个精巧设计的、直接映射到可实现的粒状实验系统的微观模型,首次在同一框架内清晰地实现了活性物质非生长型界面全部三种理论预言的普适类,解决了该领域长期存在的理论与实验脱节问题。其二,它揭示了由主体相固相有序化或玻璃化转变引发的慢动力学可以深刻重塑界面静态特性,并由此发现了全新的普适类。
本研究的科学价值在于,它为跨越非平衡界面物理的理论与实验鸿沟提供了关键的缺失环节。该模型不同于许多需要极长持续长度或界面极化才能观测到普适标度的典型活性物质系统,它可能为实验学家提供一条更容易通往渐近标度行为的途径,从而推动对非平衡普适类的定量验证。在应用层面,活性物质中的液-固和液-玻璃界面广泛存在于生物学体系中,例如细胞内生物分子凝聚体与细胞骨架的相互作用、组织力学等领域,因此,理解慢动力学如何影响界面起伏对理解这些生命过程具有潜在意义。此外,该模型框架和发现为未来的理论研究(如发展能关联主体相弛豫和界面指数的解析理论)和实验探索(如在准二维振动颗粒系统中寻找这些普适类)开辟了新的方向。