Orienting-Field Effects on Instability and Mode Selection in Active Nematics

Technical Summary: Orienting-Field Effects on Instability and Mode Selection in Active Nematics

Problem Statement
Active nematics are non-equilibrium systems composed of motile agents that exhibit orientational order, often modeled by incorporating active stresses into the hydrodynamics of passive nematics. While the spontaneous-flow transition—where a quiescent active nematic becomes unstable above a critical activity threshold—is well-established, the impact of external orienting fields on these confined systems remains theoretically underexplored. Biological and synthetic active matter frequently operate in confined environments subject to external cues (e.g., electric, magnetic, or optical fields, or anisotropic substrates). This work addresses the gap in understanding how such orienting fields interact with activity and confinement to modify instability onset, bifurcation structures, and mode selection in active nematics.

Methodology
The authors employ a minimal Ericksen–Leslie framework adapted for low Reynolds number nemato-hydrodynamics. The model describes the dynamics of the director field $\mathbf{n}$ and velocity field $\mathbf{v}$ in a layer of active nematic confined between two parallel plates ($z = \pm d/2$) with infinite planar anchoring.

• Governing Equations: The system includes active stresses ($\zeta$) and a field-induced torque term ($\chi F^2$) in the angular momentum balance. The orienting field is modeled to align the director either parallel to the substrate normal (z-direction) or parallel to the substrate anchoring (x-direction).
• Linear Stability Analysis: The authors linearize the governing equations about the trivial planar state ($\theta=0, v=0$). They derive exact and approximate analytic forms for the stability boundaries, identifying two distinct instability modes:
  • S-modes: Antisymmetric director distortions about the layer mid-point.
  • D-modes: Symmetric director distortions about the layer mid-point.
• Numerical Simulation: To verify the linear predictions and explore nonlinear evolution, the authors perform direct numerical simulations of the full nonlinear equations using COMSOL Multiphysics. These simulations test steady states at various points in the parameter space defined by activity ($\zeta$) and field strength ($F$).

Key Contributions and Results

1. Breaking of Mode Degeneracy: In the absence of an external field, the spontaneous-flow transition at the critical activity $\zeta_c=1$ involves a double pitchfork bifurcation where both S and D modes destabilize simultaneously, with the S-mode typically growing fastest. The introduction of an orienting field breaks this degeneracy, shifting the critical thresholds for S and D modes independently and restructuring the bifurcation landscape.

2. Cooperative Instability (Destabilizing Field): When the orienting field aligns the director perpendicular to the substrate anchoring (normal to the plates), it acts cooperatively with activity:

   • Lowered Thresholds: Instability can occur when both the activity and field strength are below their respective classical critical values ($\zeta < 1$ and $F < 1$). The field lowers the activity threshold, and activity lowers the field threshold for the Fréedericksz-type transition.
   • Mode Selection: The field strength determines the symmetry of the fastest-growing mode. At low field strengths (but above the cooperative instability threshold), the D-mode (symmetric) dominates. As the field strength increases, the system transitions to a regime where the S-mode (antisymmetric) becomes the fastest-growing mode. This allows for the selection of specific flow symmetries by tuning the field strength alone.

3. Stabilization (Aligning Field): When the orienting field aligns the director parallel to the substrate anchoring, it stabilizes the trivial planar state. In this regime, a sufficiently strong field can suppress instabilities entirely, even at high activity levels. Any remaining instability is dominated by the S-mode.

4. Validation of Linear Theory: Numerical solutions of the full nonlinear system confirm that the linear stability analysis correctly predicts the symmetries of the long-time states. The simulations demonstrate that the system evolves into the steady states (S-mode or D-mode) predicted by the relative growth rates derived from the linear analysis, even in the cooperative instability region where both parameters are sub-critical.

Significance
The paper claims that orienting fields do more than simply shift the instability threshold; they qualitatively restructure the bifurcation landscape of confined active nematics. By removing the mode degeneracy inherent in the zero-field case, orienting fields provide a mechanism to control the symmetry of the emerging flow states. This offers a physical framework for interpreting experiments on field-responsive active matter and suggests strategies for tuning instability onset and selecting modes in confined active nematic systems without altering material properties or geometry. The work highlights a distinct difference between activity-driven transitions and field-driven transitions (Fréedericksz), noting that the "active Fréedericksz" transition possesses a unique bifurcation structure where activity and field can cooperatively drive instability below classical thresholds.