Defect Localization by Stress Anisotropy in Active Nematic Turbulence

Technical Summary: Defect Localization by Stress Anisotropy in Active Nematic Turbulence

Problem Statement
Active systems, ranging from cytoskeletal filaments to cellular monolayers, exhibit complex collective behaviors such as pattern formation and active turbulence. While liquid crystal hydrodynamics (specifically active nematics) has successfully modeled the orientational order and topological defects ($\pm 1/2$ charges) in these systems, a critical gap remains in understanding the relationship between these orientational defects and the underlying mechanical stress fields. Specifically, it is unclear how the components of the stress tensor relate to the nematic director field, whether stress fields exhibit their own topological singularities, and how these stress defects localize within the turbulent regime. Addressing this is crucial because in biological contexts (e.g., confluent cell layers), mechanical stress measurements via traction force microscopy are often more accessible than detailed cell shape characterization.

Methodology
The authors employ a generic continuum model of two-dimensional active nematics suspended in an incompressible fluid. The system is governed by coupled hydrodynamic equations describing the evolution of the nematic order parameter tensor ($Q_{ij}$) and the momentum of the fluid ($v$).

• Model Equations: The dynamics include a generalized advection term ($S_{ij}$), rotational diffusivity ($\Gamma$), and a molecular field ($H_{ij}$) derived from a Landau-de Gennes free energy. The total stress tensor ($\sigma_{ij}$) comprises passive, viscous, and active components ($\sigma^a_{ij} = -\zeta Q_{ij}$), where $\zeta$ represents the activity strength (positive for extensile, negative for contractile).
• Simulation: The coupled equations are solved numerically using a hybrid Lattice-Boltzmann algorithm. Simulations are conducted in a turbulent regime with periodic boundary conditions for various activity strengths ($\zeta$).
• Analysis: The study focuses on the maximal principal stress ($\sigma_{max}$), its orientation ($\theta_p$), and the second invariant of the deviatoric stress tensor ($J_2$). The authors analyze the alignment between the nematic director ($n$) and the principal stress director ($n_s$), compute cross-correlations, and map the spatial distribution of topological defects in both fields.

Key Results

1. Alignment of Stress and Director: The study establishes a systematic relationship between the maximal principal stress and the nematic director. In extensile systems ($\zeta > 0$), the maximal principal stress aligns perpendicular to the nematic director. Conversely, in contractile systems ($\zeta < 0$), the maximal principal stress aligns parallel to the director. This alignment is robust across different activity strengths.
2. Topological Defects in Stress Fields: The orientation field of the maximal principal stress forms topological defects with charges $+1/2$ and $-1/2$, identical to those found in the nematic director field. However, the internal structure differs: the region of tension (positive isotropic stress) is centered at the core of nematic $+1/2$ defects but is displaced away from the core of stress $+1/2$ defects.
3. Defect Localization via Stress Anisotropy: A primary finding is the identification of a distinct isoline defined by the condition $\sigma_{xx} - \sigma_{yy} = 0$ (where the anisotropic stress component vanishes). Along this isoline, the second invariant of the deviatoric stress ($J_2$) reaches a minimum. The authors demonstrate that all $\pm 1/2$ defects (both nematic and stress) are localized precisely at the points where $J_2$ vanishes or is minimized.
4. Universality: This localization feature is robust; it remains unchanged regardless of the magnitude of activity or whether the system is extensile or contractile. The authors attribute this to the dominant role of the active stress term, which shares a traceless symmetric structure with the nematic tensor, causing the zeros of their respective anisotropic components to coincide.

Significance and Claims
The paper claims to provide a new framework for probing the mechanical and rheological properties of confluent cell layers. By demonstrating that stress defects are not random but organized along specific geometric backbones (the $J_2$ minima isoline), the authors suggest that mechanical stress measurements can serve as a proxy for understanding topological dynamics in active matter.

The authors emphasize that this approach offers a complementary description to traditional nematic order analysis. Since principal stresses can be experimentally inferred using techniques like traction force microscopy and monolayer stress microscopy, this stress-based defect localization provides a direct link between theoretical active nematic models and measurable mechanical quantities in biological tissues. The work does not propose new experimental protocols but rather interprets existing mechanical observables through the lens of active nematic hydrodynamics to reveal underlying topological organization.