拥挤通道中活性铰链的对称性破缺运动：从颗粒介质约束到可调控超结构设计

《Soft Matter》：Symmetry-breaking motility of an active hinge in a crowded channel

本刊推荐：针对对称约束下活性物质如何自发打破对称实现持续运动的机制问题，研究者提出“活性铰链（active hinge）”模型，结合朗之万动力学与粗粒化状态分析，揭示堆积分数?与最大铰链角θmax对运动态寿命的非单调调控及被动颗粒构型突变效应，为活性-被动混合系统的动态模式切换提供理论框架。

　　

在自然界与人工系统中，从细胞集体的定向迁移到微型机器人的集群协同，微观尺度的“自主运动”始终是非平衡物理学与软物质研究的核心命题。这类由内部能量驱动、通过环境反馈形成有序结构的系统被称为活性物质（active matter）。然而，一个基础难题长期困扰学界：即便外部环境完全对称——例如均匀的几何边界、无偏置的势场——活性单元仍可能自发打破左右对称，产生稳定的单向运动（如顺时针或逆时针旋转的颗粒团簇）。这种“对称性破缺”（symmetry breaking）背后的正反馈机制，及其在复杂介质（如致密颗粒浴）中的稳定性规律，是理解活性系统自组织原理的关键，也为设计新一代可调控微纳器件提供了灵感。

近期，Son等人发表于《Soft Matter》的实验发现令人瞩目：在填充振摇颗粒介质的圆形边界内，少量椭圆状自驱动粒子（self-propelled particles, SPP）会自发聚集于壁面，形成沿边界持久单向移动的活团簇。这一现象源于团簇不对称性与周围颗粒响应的正反馈循环——前端卡住的被动圆盘维持团簇偏斜，而偏斜的推力又将圆盘压在壁面，直至被动颗粒偶然脱离导致对称恢复、运动停止。但实验体系的复杂性限制了参数空间的系统性探索：团簇有效倾角受限于边界曲率与粒子数量，难以厘清各物理变量对运动态稳定性的独立贡献。

为此，L. G. Rigon与Y. Baek在《Soft Matter》发表研究，提出“活性铰链（active hinge）”最小模型，将实验观测的边界团簇动力学简化为两段刚性杆连接的可滑动铰链，浸没于二维通道的被动布朗圆盘浴中，以此剥离冗余因素，聚焦对称性破缺的核心驱动力。他们的工作不仅重现了实验中的持续运动态，更揭示了运动寿命随颗粒堆积分数?的非单调变化，以及临界铰链角θ_max≈53°引发的被动颗粒构型突变——这些发现为“少量活性单元调控被动介质整体动力学”的超结构设计提供了定量蓝图。

研究者采用多尺度模拟策略：首先建立活性铰链的微观珠串模型——每根杆由M个重叠小圆盘构成，中央盘共享为铰点，通过朗之万方程描述其平动（x_act）与转动（θ_L, θ_R）；引入Weeks-Chandler-Andersen（WCA）排斥势刻画颗粒间及壁面相互作用，设定杆与壁最大夹角θ_max模拟实验中的曲率约束。其次，推导过阻尼极限下的有效阻力矩阵Γ，耦合平移与旋转自由度，并通过约束力算法确保角度在[0, θ_max]范围内。仿真以σ（被动盘直径）、k_BT（热力学能量）、τ=σ²/D（扩散时间）为单位，步长Δt=10^{?2τ，统计200次系综平均。为量化运动态寿命，定义首通时间（FPT）为单杆吸收（θ→0）耗时，利用生存概率指数衰减拟合估计均值（MFPT）；同时粗粒化“杆下被动盘数N”，解析不同N态的驻留时间T_N，揭示稳定构型的涌现机制。}

3. 对称性破缺的运动态

研究首先验证了活性铰链能否复现实验中的对称性破缺。从双杆对称初始态（θ_L=θ_R=θ_max）出发，模拟显示第一根杆被吸收的MFPT显著短于第二根——前者随堆积分数?增加缓慢增长，反映碰撞增多减缓角度归零；后者则呈现非单调依赖，暗示多重竞争机制。均方位移（MSD）进一步佐证：第一杆吸收前为正常扩散（MSD∝t），之后进入持久弹道运动（MSD∝t²），对应单杆驱动的稳态迁移，完美类比实验中“不对称团簇→单向运动”的正反馈锁模。

4. 运动态持续性的调控机制

4.1 堆积分数?的竞争效应

当活性力佩克莱数Pe、杆长l、最大角θ_max及杆盘迁移率μ均较大时，MFPT在中等?区间出现明显峰值，揭示“稳定运动态窗口”。粗粒化分析阐明本质：低?时，被动盘碰撞加速杆的转动扩散，缩短MFPT（机制ⅰ）；中低?时，单被动盘（N=1）易卡入杆-壁间隙，强活性推力使其在共动系中受指向铰点的等效力，形成“捕获增强不对称→更快运动→更强捕获”的正反馈（机制ⅱ），对应实验中团簇的自持机制；中高?时，上方颗粒压力挤压杆，减小杆-壁角，削弱推力与捕获力，一旦N=1态失稳，杆迅速被压至吸收态（机制ⅲ）；接近jamming（??0.7）时，全局动力学减速拉高MFPT（机制ⅳ）。若无足够Pe、l、θ_max、μ，正反馈无法激活，MFPT峰值消失，印证实验“团簇过大→对称化→运动猝灭”的阈值行为。

4.2 最大角θ_max的构型跃变

在?=0.7（近实验条件）下，MFPT随θ_max上升，但θ_max?53°且Pe=16、l=2时出现阶跃式增长。状态分解显示：θ_max?53°时，窄扇形区仅容单盘（N=1主导）；θ_max?53°后，三被动盘可形成稳定三角形簇（N=3态驻留极长），即使单盘逸出也易补充，而三角簇投影波动贡献N=2态时间。此构型突变精确映射实验现象：θ_max模拟团簇尺寸效应，其临界跃迁对应“增活性粒子数→有效倾角减→对称化→运动猝灭”的突发损失。

活性铰链模型的系统性仿真证明，对称环境中活性-被动耦合可自发涌现方向性运动，其持续性由四重竞争机制权衡：转动扩散（促吸收）、正反馈捕获（延寿命）、上方压力诱导失稳（缩寿命）与jamming减速（再延长）。临界参数（Pe、l、θ_max、μ）共同决定是否激活强对称破缺态，而θ_max跨越53°触发被动颗粒构型从单盘到三角簇的跃变，使运动寿命量级提升。这一成果不仅为Son等人的实验提供了可计算的最小机制模型，更指明“活性掺杂被动介质”超结构的设计路径：通过调节区域面积（控?）或活性力（控Pe），可切换直线与手性运动模式，为柔性微马达、靶向递送系统及自适应材料的动力学编程奠定理论基础。