Defect Localization by Stress Anisotropy in Active Nematic Turbulence

技术摘要：活性向列相湍流中由应力各向异性驱动的缺陷定位

问题陈述
从细胞骨架丝到细胞单层，活性系统表现出复杂的集体行为，如模式形成和活性湍流。虽然液晶流体力学（特别是活性向列相）已成功模拟了这些系统中定向序和拓扑缺陷（$\pm 1/2$ 电荷）的演化，但仍存在一个关键空白，即如何理解这些定向缺陷与底层机械应力场之间的关系。具体而言，目前尚不清楚应力张量的分量如何与向列相指向矢场相关联，以及应力场是否具有自身的拓扑奇异性，以及这些应力缺陷如何在湍流状态下定位。解决这一问题对于生物学背景（例如，汇合细胞层）至关重要，因为在这些背景下，通过牵引力显微术进行的机械应力测量通常比详细的细胞形状表征更容易实现。

方法论
作者采用了一个通用的二维活性向列相连续介质模型，该模型悬浮在不可压缩流体中。该系统由描述向列相序参数（$Q_{ij}$）和流体动量（$v$）演化的耦合流体力学方程控制。

• 模型方程： 动力学包括广义平流项（$S_{ij}$）、旋转扩散率（$\Gamma$）和源自 Landau-de Gennes 自由能的分子场（$H_{ij}$）。总应力张量（$\sigma_{ij}$）包含被动、粘性和活性分量（$\sigma^a_{ij} = -\zeta Q_{ij}$），其中 $\zeta$ 代表活性强度（正值表示伸展型，负值表示收缩型）。
• 模拟： 使用混合格子玻尔兹曼算法对耦合方程进行数值求解。在不同活性强度（$\zeta$）下，在周期性边界条件下对处于湍流状态的系统进行模拟。
• 分析： 研究重点在于最大主应力（$\sigma_{max}$）、其方向（$\theta_p$）以及偏应力张量的第二不变性（$J_2$）。作者分析了向列相指向矢（$n$）与主应力指向矢（$n_s$）之间的对齐情况，计算了互相关，并绘制了两个场中拓扑缺陷的空间分布图。

主要结果

1. 应力与指向矢的对齐： 研究确立了最大主应力与向列相指向矢之间的系统性关系。在伸展型系统（$\zeta > 0$）中，最大主应力与向列相指向矢呈垂直关系。相反，在收缩型系统（$\zeta < 0$）中，最大主应力与指向矢呈平行关系。这种对齐在不同的活性强度下都是稳健的。
2. 应力场中的拓扑缺陷： 最大主应力的方向场形成了与向列相指向矢场相同的 $+1/2$ 和 $-1/2$ 电荷拓扑缺陷。然而，其内部结构有所不同：张力区域（正各向同性应力）中心位于向列相 $+1/2$ 缺陷的核心处，但相对于应力 $+1/2$ 缺陷的核心则发生位移。
3. 通过应力各向异性进行缺陷定位： 一个主要发现是确定了一个由条件 $\sigma_{xx} - \sigma_{yy} = 0$（即各向异性应力分量消失处）定义的独特等值线。沿着该等值线，偏应力的第二不变性（$J_2$）达到极小值。作者证明了所有 $\pm 1/2$ 缺陷（包括向列相和应力缺陷）都精确地定位在 $J_2$ 消失或最小化的点上。
4. 普适性： 这种定位特征是稳健的；无论活性大小或系统是伸展型还是收缩型，该特征均保持不变。作者将其归因于活性应力项的主导作用，该项与向列相张量具有相同的无迹对称结构，导致它们各自各向异性分量的零点相重合。

意义与主张
本文声称为探测汇合细胞层的机械和流变性质提供了一个新框架。通过证明应力缺陷并非随机分布，而是组织在特定的几何骨架（$J_2$ 极小值等值线）之上，作者表明机械应力测量可以作为理解活性物质中拓扑动力学的代理指标。

作者强调，这种方法为传统的向列相序分析提供了一种补充描述。由于主应力可以通过牵引力显微术和单层应力显微术等技术进行实验推断，因此这种基于应力的缺陷定位为理论活性向列相模型与可测量的生物组织机械量之间建立了直接联系。这项工作并非提出新的实验方案，而是通过活性向列相流体力学的视角来解释现有的机械观测量，从而揭示底层的拓扑组织。