Aspect-ratio-dependent void formation in active rhomboidal and elliptical particle systems

技术摘要：活性菱形与椭圆形颗粒系统中依赖于长宽比的空隙形成

问题陈述
活性向列相系统（由将化学能转化为动能的自驱动颗粒组成）表现出复杂的时空有序性和密度不均匀性。虽然先前的研究已经确定了拓扑缺陷驱动这些动力学过程，并且在流体动力学模型中颗粒的长宽比会影响集体行为，但关于颗粒几何形状（超出简单长宽比之外）如何决定“干”活性向列相系统（其中相互作用由直接碰撞和排斥体积效应主导）中“空隙区域”（即局部颗粒密度接近于零的区域）的形成机制，目前仍不完全明确。具体而言，在中等密度系统中，空隙区域的出现——即局部密度几乎为零的区域——其尺寸和稳定性是如何受颗粒形状（不仅是长宽比）以及直接接触力学的调节作用影响的，目前尚未得到充分表征。

研究方法
作者对一个具有周期性边界条件的二维正方形系统中的活性向列相模型进行了数值模拟。该系统由 $N=4140$ 个仅通过排斥体积效应进行相互作用的自驱动颗粒组成。

• 颗粒形状： 对比了两种不同的几何形状：菱形颗粒和椭圆形颗粒。
• 参数： 保持颗粒面积（$s = \pi/23$）不变，同时将长宽比（$a$）从 1.5 变化到 6.0。这要求根据每种形状相应地调整长轴长度。
• 动力学： 颗粒运动通过包含自驱动、随机旋转和随机位移的尝试步进行建模。重叠被拒绝，以强制执行排斥体积约束。一个记忆变量（$\Lambda$）用于追踪重叠历史，以影响随后的运动。
• 分析： 研究分析了局部数密度、局部向列相序、拓扑缺陷缠绕数以及空隙区域特征。空隙区域被定义为局部密度低于系统平均值 30% 的网格点连通簇。关键指标包括平均空隙面积、空隙尺寸分布，以及全局和针对空隙周围颗粒计算的取向相关长度。

主要结果

1. 空隙出现： 只有当长宽比超过临界阈值时（菱形约为 $a \approx 3.0$，椭圆形约为 $a \approx 4.0$），两个系统才会出现空隙区域。局部密度方差和空隙区域的总面积随长宽比增加而增加，并最终趋于饱和。
2. 形状依赖的空隙特征：
   • 菱形系统： 表现出特征性的空隙尺寸，该尺寸随长宽比单调增加。空隙尺寸分布较窄，且总空隙面积在 $a=5.0$ 左右趋于饱和。
   • 椭圆形系统： 显示出没有明确特征尺度的宽广空隙尺寸分布。椭圆形系统的空隙数密度低于菱形系统，且总空隙面积相对于长宽比并未表现出与菱形系统相同的明显饱和行为。
3. 拓扑缺陷： 两个系统都由 $\pm 1/2$ 缠绕数的缺陷主导。然而，菱形系统包含更高密度的缺陷。空隙区域在两个系统中通常与 $+1$ 缠绕数相关联，并且周围的颗粒通常沿空隙边缘切向排列。
4. 取向相关性： 取向相关长度（$\lambda$）在两个系统中都随长宽比增加。至关重要的是，专门位于空隙周围颗粒的取向相关长度（$\lambda^{(void)}$）在菱形系统中显著高于椭圆形系统。颗粒沿空隙边缘的排列程度（由 $\langle \bar{\beta} \rangle$ 量化）在菱形系统中更强。

意义与主张
本文声称，在干活性向列相系统中，颗粒的具体几何形状（菱形对比椭圆形）在决定宏观集体模式方面起着至关重要的作用，这与单纯的长宽比效应不同。

• 空隙定标机制： 作者提出，在菱形系统中，空隙区域的特征尺寸是由空隙周围颗粒场的取向相关长度决定的。这得到了以下发现的支持：对于菱形颗粒，观测到的平均空隙面积与基于取向相关长度估算的面积之比保持近乎常数。
• 椭圆形系统的对比： 与之相反，椭圆形颗粒的空隙尺寸分布很宽，这表明空隙尺寸可能由系统尺寸而非局部取向相关长度决定。
• 模型启示： 研究结果表明，仅依赖长宽比的活性向列相模型（在流体动力学模型中很常见）对于受直接碰撞主导的系统可能是不充分的。颗粒几何形状的细节（它影响了颗粒在碰撞过程中如何堆积和排列）是空隙形成等时空模式出现的关键因素。

研究结论认为，虽然长宽比和颗粒形状都会影响集体运动，但由于取向相关性和堆积约束的差异，决定空隙尺寸的具体机制在菱形颗粒和椭圆形颗粒之间存在本质区别。