Collective behavior of squirmers in thin films

本研究采用 squirmer 模型和耗散粒子动力学，研究了游泳者形状、体积分数、流体动力相互作用以及转子偶极子如何影响受限薄膜中细菌的集体行为——涵盖从类气态相到集群现象及运动诱导相分离的过程——并揭示了不对称的结构形成以及转子偶极子对不同类型游泳者之间差异的缓解作用。

要点

想象一个拥挤的舞池，但舞者不是人，而是微小的细菌。这些可不是普通的舞者，它们是“自驱动”的，这意味着它们拥有自己的小引擎，可以自主游动。科学家们称之为“蠕动者”（squirmers）。

这篇论文就像是对这些游泳细菌在被挤进一个非常薄、非常扁的水平水层（像一层三明治一样薄）中会发生什么的某种高科技观察。这个水层薄到它们只能在两个层级中移动，就像双层巴士一样，但仍有足够的空间进行旋转。

以下是研究人员发现的内容，用简单的语言进行了解释：

1. 三种类型的舞者

研究人员研究了三种不同的“性格”的游泳细菌，其区别在于它们如何通过推动水来移动：

• 推手型（如大肠杆菌）： 它们通过将水从尾部推开来向前移动。想象一下，它们就像是在人群中推搡着通过门口的人。
• 拉手型： 它们通过将水向头部拉动来向前移动。想象一下，它们就像是使用绳子把自己拉过人群的人。
• 中性型： 它们只是滑行，不会产生明显的推力或拉力。

2. “薄膜”舞池

科学家们将这些游泳者放在两堵墙之间的一个狭窄间隙中。

• 结果： 细菌自然地形成了两层，一层靠近顶壁，一层靠近底壁。它们并没有堆叠成高塔，而是保持在这些扁平的薄片中。
• 取向： 大多数情况下，细菌平行于墙壁游泳，就像在高速公路上行驶的汽车。然而，“拉手型”稍微有点叛逆；在人群稀疏时，它们喜欢几乎垂直于墙壁站立（生长），就像从花盆里长出的花朵。

3. 人群如何改变舞蹈

研究人员改变了盒子中细菌的数量（体积占比），以观察人群的行为：

• 低密度人群（气相）： 当细菌很少时，它们只是随机游动，就像在大型空旷公园里闲逛的人。
• 中等密度人群（集群）： 随着加入更多细菌，它们开始形成协同移动的群体。这被称为“集群”（swarming）。这就像鱼群或鸟群同步移动一样。
  • 转折点： “推手型”和具有特殊旋转特征（称为“转子偶极子”，rotlet dipole）的细菌非常擅长集群。而没有这种旋转特征的“拉手型”不太擅长集群；它们更喜欢聚集成紧密的、静止的团块。
• 高密度人群（交通堵塞）： 当盒子被塞得满满的时，细菌会陷入巨大的、无法移动的簇群中。它们无法再移动了。这被称为“运动诱导相分离”（MIPS）。这就像一场交通堵塞，每个人都卡在一个巨大的、无法移动的堆积物中。

4. “陀螺效应”（转子偶极子）

涉及到一个特定的流场——“转子偶极子”——的研究非常有趣。想象一个细菌不仅向前游动，还在移动时像陀螺一样旋转身体。

• 神奇之处： 当研究人员加入这种旋转运动时，它就像是一个“通用平衡器”。无论细菌是推手型、拉手型还是中性型，它们表现得都一样。
• 结果： 旋转让它们变得活跃得多。它们不再形成那些紧密的、被困住的团块，而是保持活跃移动。它们也更频繁地在顶层和底层之间跳跃，就像在高速公路上变换车道以避开交通拥堵一样。

5. 为什么这对于生物膜很重要

生物膜是你在牙齿（牙菌斑）或溪流中的岩石上看到的粘性细菌层。它们从表面的一层细菌开始生长。

• 核心问题： 它们是如何长成厚厚的、多层结构的？
• 答案： 研究表明，如果细菌是“推手型”（如大肠杆菌）或者具有那种“陀螺旋转”的运动，它们非常擅长从底层跳跃到顶层。这使得它们能够快速构建出多层结构。
• 例外情况： “拉手型”倾向于停留在它们特定的模式中，不像其他类型那样容易切换层级，这可能会减慢它们生物膜增厚的速度。

总结

简而言之，这篇论文展示了细菌的形状、它们推动水的方式以及它们在游泳时是否旋转，共同决定了它们会形成混乱的人群、协调的集群，还是被困住的交通堵塞。“旋转”这一特征特别强大，因为它能让细菌保持移动并防止被困，从而帮助它们构建出更厚、更复杂的结构。

技术摘要

问题陈述
微型游泳者（如生物膜中的细菌）的集体运动是一个复杂的现象，由形态特征、推进力、体积分数以及流体动力学和位阻相互作用共同驱动。虽然在准二维（quasi-2D）单层和三维（3D）周期系统中，运动诱导相分离（MIPS）、集群行为（swarming）和主动湍流已被广泛研究，但对于允许多层结构的薄膜中游泳者的行为研究仍不够深入。具体而言，需要填补高度受限的准二维系统与本体三维系统之间的空白，以理解细菌如何从表面单层过渡到多层生物膜结构。关键问题包括：薄膜中的游泳者更倾向于平行还是垂直于壁面的取向？在什么条件下它们会自发离开壁面形成多层结构？以及它们的集体行为与准二维系统有何不同？

方法论
作者采用 squirmer 模型来表示微型游泳者，并使用耗散粒子动力学（DPD）方法进行流体建模。系统由位于两个平行壁面之间的球形 squirmer 组成，其薄膜厚度（$L_y$）足以允许完全的旋转自由度，但又受到限制以确保最多只能形成两层结构。

• Squirmer 模型： 游泳者被建模为长宽比为 2（$b_z/b_x = 2$）的球体，模拟大肠杆菌（E. coli）。其表面滑移速度包含自推进项（$B_1$）、用于区分推力型（$\beta < 0$）、中性型（$\beta = 0$）和拉力型（$\beta > 0$）游泳者的主动应力项（$\beta$），以及一个用于模拟鞭毛细菌反向旋转流的转子偶极子项（$\lambda$）。
• 模拟设置： 模拟在 $x$ 和 $z$ 方向具有周期性边界条件、并在 $y$ 方向受壁面限制的区域内进行。研究通过改变体积分数（$\phi \in \{0.18, 0.35, 0.44, 0.56\}$）、主动应力（$\beta \in \{-5, 0, 5\}$）和转子偶极子强度（$\tilde{\lambda} \in \{0, 133.5\}$）来进行。
• 分析： 通过簇大小分布、径向分布函数和取向分布来分析结构特性。动态特性通过有效旋转扩散、均方位移（MSD）和平均游泳速度来表征。

主要贡献与结果
研究确定了随着体积分数增加和游泳者类型的变化，系统呈现出不同的集体状态——气态相、集群行为和 MIPS。

1. 相行为与体积分数：
   • 在低体积分数（$\phi \approx 0.18$）下，系统通常表现为小簇的气态相，但没有转子偶极子的拉力型游泳者由于吸引性的流体动力学相互作用会形成大簇（MIPS）。
   • 在中等体积分数（$\phi \approx 0.35$）下，推力型游泳者和带有转子偶极子的游泳者表现出集群行为。没有转子偶极子的拉力型则趋向于 MIPS。
   • 在高体积分数（$\phi \ge 0.44$）下，所有系统都向以大型、近乎静止的簇为特征的 MIPS 相转变。
2. 游泳者类型与壁面相互作用的作用：
   • 推力型和中性型游泳者： 倾向于平行于壁面排列。推力型经常在两层之间切换，促进了多层结构的形成。
   • 拉力型游泳者： 在没有转子偶极子的情况下，拉力型倾向于近乎垂直于壁面的取向，导致在低 $\phi$ 时形成“花状”结构（一个中心拉力型垂直于壁面，周围环绕着“花瓣”状拉力型）。这种取向促进了早期簇的成核和 MIPS。拉力型很少自发切换层级。
3. 转子偶极子的影响：
   • 转子偶极子的存在（$\tilde{\lambda} = 133.5$）显著减轻了推力型、中性型和拉力型游泳者之间集体行为的差异。
   • 存在转子偶极子时，所有类型的游泳者表现出相似的行为：它们平行于壁面排列，在中等 $\phi$ 时表现出集群行为，并在高 $\phi$ 时向 MIPS 转变。
   • 转子偶极子诱导了壁面附近的圆形轨迹和频繁的层间切换，从而增强了即使在中等体积分数下的有效旋转扩散和迁移率。
4. 结构不对称性： 集体结构的形成并不一定相对于两个壁面对称，尤其是在缺乏转子偶极子的情况下。

意义
本文声称提供了连接极度受限的准二维系统与较少受限的三维集体行为的第一步。研究结果强调了游泳者的各向异性、流体动力学相互作用以及壁面相互作用对于决定集体状态的关键作用。具体而言，研究表明：

• 游泳者的长宽比及其与壁面的流体动力学相互作用对于不同集体行为的出现至关重要。
• 转子偶极子可以使不同类型的游泳者行为趋于一致，减少推力型、中性型和拉力型游泳者之间的区别。
• 类推力型游泳者（类似于大肠杆菌）由于频繁的层间切换，能够自发地从单层配置转变为多层配置，而拉力型则倾向于保持特定的取向以促进簇的成核。

作者指出，虽然 squirmer 模型能很好地捕捉远场流体动力学，但对于真实细菌的近场细节（如鞭毛几何形状和翻滚行为）进行了简化。因此，直接与实验进行定量比较需要进行校准，但这些模拟为受限活性物质提供了稳健的定性表征。