Asymmetry-Induced Chiral Dynamics in Coupled Self-Propelled Robots: Spinning and Circular Motion

技术摘要：耦合自驱动机器人中不对称诱导的手性动力学

问题陈述
微观生物游泳者，如大肠杆菌（E. coli）、莱氏衣藻（Chlamydomonas reinhardtii）和精子，表现出多样化的运动策略，包括趋跑与翻转（run-and-tumble, RnT）运动、圆形轨迹和旋转。这些行为通常由结构不对称性（例如，偏心的鞭毛附着或形态不对称）驱动，从而将推进力与净力矩耦合。虽然之前的极简机器人模型（如 Paramanick 等人）通过两个由刚性杆连接的自驱动机器人成功复制了 RnT 运动，但一个关键问题仍然存在：通过控制机械不对称性，在具有弹性耦合的系统中，是否可以产生圆形运动和其他独特的动力学机制？此外，弹性顺应性（而非刚性约束）在切换运动行为中的功能作用尚未得到充分探索。本文研究了两对由弹簧耦合、具有不对称推进方向的自驱动布朗机器人的动力学，以解决这些空白。

方法论
作者开发了一个理论模型，通过将刚性杆替换为劲度系数为 $k$ 的胡克弹簧，扩展了 Paramanick 等人的刚性杆框架。该系统由两个主动布朗机器人组成，每个机器人具有以下特征：

• 自推进： 一个恒定的主动力 $F_a$ 沿着取向向量 $n_i$ 作用。
• 几何不对称性： 枢轴偏移量 $\delta$（机器人中心与弹簧附着点之间的距离）和推进角 $\alpha_i$（推进方向与连接枢轴与中心连线的夹角）。
• 耦合： 机器人通过一个静止长度为 $l$、劲度系数为 $k$ 的弹簧连接。
• 噪声： 引入了旋转扩散 $D_r$ 以模拟随机重定向。

运动方程是在过阻尼极限下推导的，忽略了惯性项。该模型考虑了各向异性的平移摩擦（平行于机器人取向的摩擦系数不同）。研究利用数值模拟分析了系统在推进角（$\alpha_1, \alpha_2$）、弹簧劲度（$K$）和旋转噪声（$D_r$）参数空间中的行为。关键观测指标包括平均自旋角速度（$\langle \omega_s \rangle$）、轨道角速度（$\omega_c$）、平均质心速度（$\langle v_c \rangle$）以及取向自相关函数（$C_0(\tau)$ 和 $C(\tau)$）。

主要贡献与结果

1. 涌现出三种动力学机制：
   系统表现出三种截然不同的机制，这主要取决于几何配置（$\alpha_1, \alpha_2$）：

• 趋跑与翻转 (RnT)： 发生在对称配置（$\alpha_1 = \alpha_2$）中。系统在直线趋跑和重定向事件之间交替，类似于刚性杆模型。
• 旋转运动 (Spinning Motion)： 当 $\alpha_1 + \alpha_2 = 360^\circ$（推进方向反向对齐）时发生。在这种机制下，净主动力消失，几何力矩诱导快速的原地旋转，且平移几乎可以忽略不计。手性（顺时针或逆时针）由 $\alpha_1 - \alpha_2$ 的符号决定。
• 圆形运动 (Circular Motion)： 发生在不对称配置（$\alpha_alpha_1 + \alpha_2 \neq 360^\circ$ 且 $\alpha_1 \neq \alpha_2$）中。在此机制下，非零的净主动力与净力矩并存，驱动系统的质心沿持续的曲线轨迹运动。曲率半径和手性受推进角控制。

2. 几何不对称性与不动点分析：
   作者分析了 $(\theta_+, \theta_-)$ 平面（其中 $\theta_\pm$ 是机器人取向的组合）中的确定性流。他们识别出了对应于旋转状态和圆形运动状态的稳定不动点。

• 对于 $\alpha_1 + \alpha_2 = 360^\circ$，稳定不动点对应于完全反向对齐的机器人，导致纯旋转。
• 对于 $\alpha_1 + \alpha_2 \neq 360^\circ$，稳定不动点涉及非完全反向对齐的机器人，产生净力使质心漂移，从而导致圆形运动。
• 研究发现，由于存在替代稳定不动点，在某些小参数区域内，运动的手性相对于 $\alpha_1 - \alpha_2$ 的符号会发生反转。

3. 旋转噪声 ($D_r$) 的影响：
   增加旋转噪声会逐渐使推进方向随机化，从而削弱几何力矩。

• 旋转运动： 随着 $D_r$ 增加，$\langle \omega_s \rangle$ 减小。高噪声会将系统从稳定的旋转不动点推向趋跑状态。具有较大 $\alpha_1$（接近 $90^\circ$）的配置比较小 $\alpha_1$ 的配置对噪声更具鲁棒性。
• 圆形运动： 圆形运动仅在噪声水平低于阈值（$D_r \approx 0.15 \, \text{rad}^2 \text{s}^{-1}$）时持续存在。超过此阈值后，轨道角速度 $\omega_c$ 降至零，运动变为扩散运动。

4. 弹簧劲度 ($K$) 的影响：
   弹簧劲度作为机械顺应性的调节参数。

• $\langle \omega_s \rangle$ 和 $\omega_c$ 均随 $K$ 单调增加，并在 $K \to \infty$ 时趋于饱和，恢复刚性杆极限。
• 在低 $K$（高顺应性）时，弹簧无法维持必要的约束力来有效传递力矩，从而抑制了动力学响应。
• 这一发现表明，弹性顺应性能连续地在被抑制的动力学与充分发展的集体运动之间进行插值。

5. 均方位移 (MSD)：

• 圆形机制： MSD 在短时间内表现出弹道式增长，随后在中间时间段出现显著振荡（轨道限制的特征），以及受抑制的长时间扩散。
• 旋转机制： MSD 表现出弹道式增长，随后平滑过渡到扩散，没有振荡，反映了缺乏轨道限制。

意义与主张
本文声称，几何不对称性和弹性顺应性足以产生生物微游泳者所观察到的全谱系集体运动行为——从直线趋跑与翻转到手性圆形轨迹和旋转——而无需需要流体动力相互作用、化学梯度或主动感知。

作者将其工作定位为人工活性系统与生物微游泳者之间的桥梁。通过展示一种简单的机械设计（两个由弹簧耦合、具有不对称推进方向的机器人）可以复制复杂的行为，如莱氏衣藻中通过调节鞭毛摆动来实现手性切换的行为，本研究确立了一个极简的机械设计原则。该原则强调，结构不对称性结合机械约束，可以通过纯粹的机械相互作用产生多样化的运动行为。研究结果表明，生物系统中的弹性耦合（例如莱氏衣藻中的远端纹状纤维或细菌鞭毛钩）可能在调节运动机制方面发挥功能性作用，这一假设得到了模型的支持，因为该模型能够通过调节劲度来实现这些机制的转换。