Growing, Buckling, and Swirling: motility from polymerization

技术摘要：生长、失稳与旋涡：由聚合作用驱动的运动性

问题陈述
本文研究了由锚定在表面的细胞骨架丝件生长所驱动的一类独特的微型游泳者（microswimmer）推进机制。传统的微型游泳模型（如鞭毛旋转或纤毛摆动）通常规定边界条件或表面滑移速度，而本研究则探讨了主动力是如何通过聚合动力学、丝件定向动力学以及流体动力学相互作用自洽地产生的。该研究受到生物系统（如产生纤维素纤维并因失稳而产生运动性的 Acetobacter xylinum）以及由肌动蛋白聚合驱动的合成胶体颗粒的启发。核心挑战在于建立一个理论框架，能够同时解决低雷诺数（斯托克斯流）环境下，丝件动力学、取向动力学与流固耦合之间的关系，其中运动是由丝件地毯（filament carpet）的自组织动力学而非预定义的边界条件决定的。

方法论
作者开发了一个将微观聚合动力学与宏观流体动力学耦合的自洽连续介质框架：

1. 平均场动力学：系统被建模为由密集的丝件组成的“地毯”，这些丝件经历着成核、恒定速度生长 ($V_g$) 以及随机解聚（灾变）。分布函数 $\Psi(\ell, t)$ 描述了丝件长度，而广义分布 $\psi(y, \ell, p, t)$ 则描述了空间位置、长度和取向 $p$。
2. 时间尺度分离：模型假设聚合与解聚发生的时间尺度比丝件取向重定向的时间尺度更快。这允许对长度动力学进行绝热消除，从而将问题简化为取向分布 $\tilde{g}(y, p, t)$ 和极化场 $n(y, t)$ 的演化。
3. 丝件动力学：单个丝件被建模为通过扭转弹簧（刚度为 $k_\alpha$）锚定在表面的刚性棒。其重定向受杰弗里方程（Jeffery's equation，流体剪切）和恢复力矩控制，并与聚合作用产生的对抗粘性流体的有效力相平衡。
4. 流体动力学耦合：主动力被粗粒化为作用在偏置表面 $\partial S$（位于刚性本体上方平均丝件长度 $\bar{\ell}$ 处）上的有效牵引跳跃。该牵引力驱动斯托克斯流。刚性本体被视为无力且无力矩的，其平移速度 ($U$) 和角速度 ($\Omega$) 是通过自洽方式确定的。
5. 数值实现：
   • 闭合方案：使用极化 Bingham 闭合方案来近似高阶取向矩（$\langle pp \rangle, \langle ppp \rangle$），将其表示为极化向量 $n$ 的函数，从而避免了求解完整福克-普朗克方程（Fokker-Planck equation）的高昂计算成本。
   • 求解器：使用边界积分法（BIM）结合五阶求积方案以及扩展求积法（QBX）处理近奇异相互作用来求解斯托克斯方程。系统通过迭代使用 GMRES 进行求解，以处理表面牵引力与速度梯度之间的隐式耦合。

关键结果

• 线性稳定性分析：在半空间几何结构中，作者证明了垂直生长的均匀丝件地毯对长波扰动是不稳定的。聚合诱导的压缩力驱动了失稳（buckling instability），从而打破对称性，导致丝件自发倾斜并产生大规模流动。该失稳受控于破坏性的主动应力与稳定性的弹性恢复力矩及解聚速率之间的竞争。
• 固定球体动力学：对于具有均匀成核机制的刚性球体，失稳会导致一种旋涡模式，其特征是存在两个对跖的 $+1$ 型拓扑缺陷。这种构型在球体外部产生相干旋转流，类似于果蝇（Drosophila）卵母细胞中的胞质流。
• 自由球体动力学（自发旋转）：当球体可以自由移动时，对称性破缺产生的旋涡流会产生净流体动力矩，导致球体绕着自发选择的轴线进行持续的刚体旋转。
• 图案化成核与运动性：通过打破成核位点的空间对称性，作者展示了对力和力矩平衡的控制：
  • 赤道带：产生净力矩但零净力，导致纯旋转。
  • 极区斑块：产生净力和净力矩，导致圆周运动（平移与旋转结合）。
  • 四重对称性：四斑块排列抵消了净力矩同时叠加了合力，产生定向的、持久的平移。
  • 手性图案：手性排列的成核位点将旋转与平移耦合，产生同时旋转并直线前进的轨迹，模拟了自然界中的手性推进机制（如银杏精子 Ginkgo sperm）。
• 几何形状的作用：通过将分析扩展到长球体，作者发现物体形状会显著改变运动轨迹。由于阻力各向异性，带有手性图案的长球体会追踪圆周螺旋线而非直线，因为净牵引力与主轴之间的偏差会导致推进方向发生连续的进动。

意义与主张
本文为由聚合丝件集体动力学驱动的运动性建立了通用的理论基础。它证明了聚合本身就是一种通用的自推进途径，仅通过丝件取向与生长的自组织，即可产生多样化的行为（旋转、定向运动、手性游泳）。

作者声称，其框架为合成微米级游泳者提供了一种新的设计原则，即“表面代码”（成核位点的空间图案化）决定了涌现出的运动策略。该工作弥合了微观聚合动力学与宏观流体动力学之间的鸿沟，提供了一个适用于聚合驱动和马达驱动主动地毯的统一连续介质描述。此外，结果表明，通过这些机制可以实现 $0.1\text{--}1 \, \mu\text{m/s}$ 量级的游泳速度。论文最后指出，虽然目前的模型假设了“薄地毯”近似，但未来的工作需要解决丝件长度超过本体尺寸的机制，这需要混合离散-连续的方法。