Defect structures and transitions in active nematic membranes

We investigate the dynamics of active nematic liquid crystals on deformable membranes, focusing on the interplay between active stress and anisotropic curvature coupling. Using a minimal model, we simulate the coupled evolution of the nematic order parameter and membrane height. We demonstrate a continuous transition from a curvature-dominated regime, where topological defects are trapped by local deformation, to an activity-dominated regime exhibiting active turbulence. A scaling analysis reveals that the critical activity threshold $ζ_c$ scales as $α^2/κ$, where $α$ and $κ$ are the coupling constant and bending stiffness, respectively; this relationship is confirmed by our numerical results. Furthermore, we find that significant correlations between the orientational pattern and membrane geometry persist even in the turbulent regime. Specifically, we identify that “walls” in the director field induce characteristic wave-like curvature profiles, providing a mechanism for dynamic coupling between order and shape. These results offer a physical framework for understanding defect-mediated deformation in nonequilibrium biological membranes.

💡 Research Summary

In this work the authors investigate how active nematic liquid crystals interact with deformable membranes, focusing on the coupling between active stresses and an anisotropic curvature–director interaction. They formulate a minimal continuum model in the Monge gauge, describing the nematic order by a symmetric traceless tensor Q and the membrane shape by a height field h(x,y). The total free energy consists of a Landau–de Gennes term, a bending term (κ∇⁴h), and a curvature–nematic coupling α Q : ∇∇h. The dynamics are governed by low‑Reynolds‑number Stokes flow, the Q‑tensor transport equation with flow‑alignment and rotational viscosity, an active stress Σ^(a)=−ζ Q, and a purely relaxational equation for h that includes the fourth‑order bending operator and the coupling source α∇·∇·Q.

A scaling analysis of the static height equation shows that the membrane curvature induced by nematic textures scales as ∇²h ∼ αS₀/κ, where S₀ is the typical magnitude of the order parameter. Balancing the active driving force on defects (∼ζ S₀) against the curvature‑induced trapping force (∼α² S₀/κ) yields a critical activity ζ_c ∝ α²/κ.

Numerical simulations are performed on a 128 × 128 lattice using a fourth‑order Runge‑Kutta scheme for the Q‑tensor and height fields, while the incompressible Stokes equation is solved in Fourier space at each time step. Fixed parameters are A = −0.5, B = 2.0, M = 1.0, λ = 0.1, Γ_Q = 0.1, Γ_h = 0.01, κ = 100, while α is varied from 0 to 10 and ζ from 0 to 1. The initial condition is a flat membrane with random Q components.

The simulations reveal two distinct regimes. At low activity (ζ ≲ 0.3) the curvature coupling dominates: topological ±½ defects become trapped in local membrane deformations, the bending energy ⟨f_bend⟩ remains high, and the height correlation length ξ_H grows much more slowly than the nematic correlation length ξ_Q. In this regime the ratio ξ_Q/ξ_H approaches ≈0.9 in the late stage, indicating that membrane relaxation limits defect coarsening.

When ζ exceeds a threshold ζ_c that follows the predicted α²/κ scaling, the active stress overwhelms the curvature trapping. The bending energy drops abruptly, the membrane flattens, and the mean squared fluid velocity ⟨v²⟩ rises linearly with ζ, signalling the onset of active turbulence. Defects unbind and move freely, and the system reaches a statistically steady turbulent state.

A particularly striking observation is the formation of “walls” in the director field—extended regions where the orientation rotates sharply and +½ defects accumulate. These walls generate characteristic wave‑like curvature profiles in h, a direct consequence of the α Q : ∇∇h term. Even in the turbulent regime, correlations between the nematic pattern and membrane shape persist, demonstrating that the curvature–director coupling remains dynamically relevant.

The authors compile a phase diagram in the (ζ, α) plane, using the time‑averaged bending energy as an order parameter. The boundary between the curvature‑dominated “trapped” phase and the activity‑dominated turbulent phase aligns well with the theoretical ζ_c ∝ α²/κ line, confirming the scaling analysis.

Overall, the paper provides a clear physical picture of how anisotropic curvature coupling controls defect dynamics and membrane deformation in active nematic membranes. By establishing the simple scaling law ζ_c ∝ α²/κ and demonstrating its validity through extensive simulations, the work offers a valuable framework for interpreting defect‑mediated morphogenesis in biological membranes and for designing synthetic active materials where shape and order are tightly coupled.