Collective deformation of anisotropic particles with internal pulsation

本文研究了各向异性椭圆颗粒的内部脉动如何驱动致密组装体中同步变形波与多样化集体模式的涌现，并提出了一个流体动力学框架，该框架成功捕捉了受心脏组织行为启发的这些自组织动力学。

要点

想象一个拥挤的舞池，每个人都试图移动，但每个人不仅是在行走，还在不断改变自己的体型。有些人伸展得又高又瘦，而另一些人则蜷缩成又小又圆的球体。现在，想象舞池中的每个人都在像心跳一样有节奏地脉动，并试图与邻居同步他们的动作。

这就是本文研究人员正在探索的世界。他们研究的是“活性物质”——由微小单元组成的系统，这些单元利用自身的能量进行移动和改变形状，就像活体中的细胞一样。具体来说，他们观察了当这些单元是椭圆形（而非简单的圆形）并以两种不同方式脉动时会发生什么。

以下是他们发现的简要说明，使用了简单的类比：

两种类型的“舞者”

研究人员创建了这两种椭圆形粒子的模型，以观察它们在人群中的行为：

1. “挤压者”：想象一个椭圆形，当它脉动到最大尺寸时，变成一个完美的圆形。但当它收缩到最小尺寸时，它会变得非常细且拉长（像一根面条）。
2. “拉伸者”：想象一个椭圆形，当它收缩到最小尺寸时，变成一个完美的圆形。但当它生长到最大尺寸时，它会拉伸成一根又长又细的面条。

人群的三种主要“情绪”

当这些粒子在密集的人群中聚集在一起时，它们不会静止不动。根据它们的紧密程度以及它们倾听邻居的程度，整个群体会进入三种截然不同的模式之一：

• “停滞”状态（冻结的人群）：如果人群过于密集，粒子就会卡住。它们无法彼此穿过，因此它们的脉动节奏被锁定。每个人都无法有效地改变形状，整个系统冻结。
• “循环”状态（同步舞蹈）：如果空间稍多，且粒子善于彼此倾听，它们就会完美同步地脉动。它们一起膨胀和收缩，就像一个巨大的生物在呼吸。
• “波浪”状态（体育场波浪）：在中间地带，事情变得混乱但美丽。粒子并非完全同时脉动。相反，变形的波浪在人群中荡漾。想象一下“体育场波浪”，人们一个接一个地站起来和坐下，形成移动的涟漪。在这个模型中，“站起来”就是粒子拉伸或挤压。

惊喜：形状决定秩序

最有趣的发现发生在挤压者（那些在变小时成为细长面条的粒子）身上。

当挤压者的人群变得非常密集时，发生了一些特别的事情。因为它们的最小形状是一根又长又细的面条，它们自然地想要彼此排成一行，就像一盒未煮过的意大利面。这创造了一种向列序状态。

• 类比：想象一盒铅笔。如果你摇晃它们，它们可能会指向随机方向。但如果你将它们非常紧密地包装，它们会自然地并排对齐。
• 结果：挤压者在高密度下完美地彼此对齐。然而，拉伸者（那些在变小时成为圆形球体的粒子）并没有这样做。当它们变小时，它们只是圆球，因此没有理由排成一行。它们保持无序。

“流体动力学”地图

研究人员不仅观察了粒子，还构建了一个数学“地图”（流体动力学理论）来预测这些行为。将此地图想象成人群的天气预报。它成功预测了：

1. 你可以得到波浪、停滞或同步循环。
2. 只有“挤压者”在人群非常密集时会自然排成一行（形成向列序）。

为什么这很重要（根据论文）

该论文表明，这有助于我们理解活体组织（如心肌）如何工作。心肌细胞是椭圆形的，并且沿着它们的长轴收缩（挤压）。研究人员发现，这种特定类型的“挤压”形状变化很可能是帮助这些细胞组织起来并产生健康心跳所需的波浪的原因，即使它们没有物理地从一处移动到另一处。

简而言之：形状即命运。一个脉动粒子是“挤压者”还是“拉伸者”，不仅决定了它如何移动，还决定了它能否组织成协调的波浪状模式或完美对齐的直线。

技术摘要

技术摘要：具有内部脉动的各向异性粒子的集体变形

问题陈述
本文旨在解决捕捉收缩性活体组织中变形波涌现的挑战，特别关注那些细胞运动（自推进）可忽略不计、但内部形状变化驱动集体动力学的系统。虽然先前的脉动活性物质模型大多考虑了各向同性粒子，但生物细胞（例如心肌细胞）通常是各向异性的，并沿特定轴改变形状。作者研究了脉动粒子的各向异性变形——特别是向列序与同步内部脉动之间的相互作用——如何影响致密组装体的集体状态。

方法论
本研究采用微观数值模拟与宏观流体动力学理论相结合的双重方法。

1. 微观模型：

   • 系统由二维空间中的 $N$ 个过阻尼各向异性粒子（椭圆）组成。
   • 动力学： 粒子位置（$\mathbf{r}_i$）和取向（$\theta_i$）在排斥相互作用和热噪声的作用下演化。
   • 变形： 通过内部相位变量 $\phi_i$ 引入两种不同的变形模式：
     • 挤压： 短轴变化，长轴保持不变。粒子在最大尺寸时呈圆形，在最小尺寸时呈椭圆形。
     • 拉伸： 长轴变化，短轴保持不变。粒子在最大尺寸时呈椭圆形，在最小尺寸时呈圆形。
   • 内部驱动： 相位 $\phi_i$ 根据一个驱动方程演化，该方程包含自然频率 $\omega$、邻居间的同步耦合 $\epsilon$ 以及噪声。
   • 相互作用： 使用粒子间距离的近似解析表达式来计算排斥力，同时考虑椭圆的取向和瞬时形状。

2. 流体动力学理论：

   • 利用随机微积分对微观运动方程应用粗粒化程序。
   • 该理论推导出一组关于密度场 $\rho$、同步场 $f_\phi$（与相位相干性相关）和向列场 $f_\theta$（与取向序相关）的封闭流体动力学方程。
   • 推导过程涉及高阶模式的绝热消除，以及关于耦合参数 $\nu$ 的微扰展开，该参数将粒子形状各向异性与向列排列联系起来。

主要贡献与结果

• 相图与集体状态：
  微观模型的数值模拟揭示了取决于密度（$\rho_0$）和同步强度（$\epsilon$）的三种主要集体状态：

  1. 阻滞态： 在高密度下，排斥力抵消了驱动力，导致相位同步（$R_\phi \approx 1$），但电流消失（$\Upsilon \approx 0$）。
  2. 循环态： 在中等密度和高同步下，粒子表现出同步相位，且变形电流最大（$\Upsilon \approx 1$）。
  3. 变形波： 在阻滞与循环之间的区域（低同步、高电流），系统表现出不稳定性，导致传播的变形波。这些波可以形成包含拓扑缺陷（顺时针或逆时针）的螺旋图案，或不包含缺陷的圆形图案。

• 各向异性与向列序的作用：
  一个关键发现是，在极高密度下出现了向列有序态（$R_\theta \approx 1$），但这仅发生在挤压变形模型中。在此区域，粒子达到其最小尺寸，对应于椭圆形，从而允许空间位阻排列诱导向列序。相反，拉伸模型在高密度下导致粒子呈圆形（最小尺寸），因此不产生向列序。

• 流体动力学验证：
  推导出的流体动力学方程在定性上重现了模拟中观察到的相图。该理论成功捕捉了：

  • 无序态、阻滞态和循环态之间的转变。
  • 向列序在挤压情况下高密度时的特定涌现。
  • 在循环态与阻滞态之间同步场中传播波和缺陷的形成。
  • 微扰分析表明，相位同步与向列序之间的耦合是驱动挤压模型与拉伸模型之间不同行为的机制。

意义
本文声称，其模型提供了关键见解，阐明了活性各向异性变形如何产生自组织成各种动力学模式的波。通过将脉动活性物质的框架扩展至包含各向异性形状，作者证明了仅靠形状变化（无需自推进）即可驱动复杂的集体行为，包括向列排序和波传播。

作者指出，这些发现对于建模同时表现出面积和形状变化的生物组织（如心脏组织和上皮层）具有相关性。该工作阐明了内部脉动与形状各向异性如何相互作用以产生此类系统中观察到的集体动力学，为理解心脏心律失常等背景下细胞运动并非主要驱动力的变形波提供了理论基础。