Activity driven buckling and pattern formation in shells of oriented solids

Technical Summary: Activity Driven Buckling and Pattern Formation in Shells of Oriented Solids

Problem Statement
Mechanical instabilities are a fundamental mechanism for structure and function generation in materials, ranging from passive wrinkling membranes to shape-morphing metamaterials. While the interplay between geometry and elasticity in passive systems is well-understood, the mechanics of active solids—materials where constituents continuously generate internal stresses within an elastic network—remain largely unexplored, particularly regarding deformable shells. Unlike active fluids, active solids couple orientational order (nematic alignment) directly to mechanical deformation and shape changes. This paper addresses the gap in theoretical understanding of how activity, orientational order, and shell geometry interact to drive instabilities and pattern formation in deformable cylindrical shells of active nematic solids.

Methodology
The authors develop a continuum mechanical model for a thin cylindrical shell populated by oriented active particles. The system is governed by a free energy functional $F$ that includes:

1. Shell elasticity: Bending energy (flexural rigidity $D$) and stretching energy (Young's modulus $E$, Poisson's ratio $\nu$).
2. Nematic ordering: Landau-de Gennes free energy terms ($A$ for ordering strength, $L$ for Frank elastic constant).
3. Coupling: A term ($\lambda$) coupling the nematic order parameter $Q_{ij}$ to the strain tensor $\epsilon_{ij}$.
4. Active stress: A contribution $-\tilde{\zeta}Q_{ij}$ representing internally generated stresses, where the sign of $\tilde{\zeta}$ distinguishes extensile ($\zeta > 0$) from contractile ($\zeta < 0$) activity.

The dynamics are modeled in the overdamped limit, describing in-plane and radial displacements ($u_j, w$) and the evolution of the nematic tensor $Q_{ij}$ driven by a molecular field $H_{ij}$.

The study employs two primary analytical and numerical approaches:

• Linear Stability Analysis: The authors perturb a pre-buckled steady state (defined by uniform nematic orientation $\psi$ and applied loads) to derive a growth rate $\sigma$ for buckling modes characterized by axial ($\alpha$) and azimuthal ($n$) wavenumbers. This yields an analytical expression for the critical activity required to trigger instability.
• Nonlinear Simulations: Full numerical simulations are performed to validate linear predictions and explore the long-time behavior of the system beyond the onset of instability, capturing self-organized patterns and dynamic states.

Key Contributions and Results

1. New Class of Activity-Driven Buckling Instabilities:
   The linear stability analysis reveals that active stresses drive buckling modes absent in passive shells. The unstable mode is selected by the interplay between the nematic orientation ($\psi$) and the sign of activity ($\zeta$).

   • Mode Selection: The system exhibits three distinct deformation modes:
     • Axial modes: Occur when the director is perpendicular to the cylinder axis (contractile activity).
     • Circumferential modes: Occur when the director is parallel to the axis (contractile activity).
     • Helical modes: Occur at intermediate angles.
   • Arbitrary Small Activity Instability: A remarkable finding is that circumferential modes can become unstable at arbitrarily small activity levels. This occurs because these modes involve no axial variation, thereby incurring no stretching (membrane) energy cost; the destabilizing active term (scaling as $k^2$) dominates the stabilizing bending term (scaling as $k^4$) at long wavelengths.

2. Critical Activity and Phase Diagrams:
   The authors derive a critical activity threshold $\zeta_c$ that depends on the director angle $\psi$ and shell geometry. The stability boundary forms a "rhodonea" or rose-curve structure in the $(\zeta, \psi)$ plane. Numerical simulations confirm that the analytically predicted stability boundaries align closely with the transition between stable and unstable states observed in simulations.

3. Emergent Nonlinear Patterns and Dynamics:
   Beyond the linear regime, nonlinear simulations uncover rich spatiotemporal behaviors driven by the feedback between strain and nematic alignment:

   • Diamond Patterns: For positive coupling ($\lambda \ge 0$), the system settles into stationary, self-organized diamond-shaped patterns.
   • Dynamic States: For negative coupling ($\lambda < 0$), the feedback creates frustration between the preferred direction of active stress and the nematic alignment. This destabilizes steady states, leading to:
     • Bulk Oscillations: Periodic switching of surface deformation and nematic orientation.
     • Traveling Waves: Propagating domain walls and waves along the cylinder axis, driven by a persistent phase lag between nematic reorientation and mechanical deformation.

Significance
The paper establishes fundamental buckling modes and emergent patterns in shells of active oriented solid materials. The authors claim these results have potential relevance to two main domains:

• Biological Systems: The findings offer a mechanical explanation for shape formation and dynamics in biological tissues that behave as active solids, such as epithelial tissues, endothelial tubes, and cytoskeletal structures. Specifically, the mechanism of activity-orientation-geometry interplay may contribute to phenomena like peristalsis in tubular organs without invoking complex biochemical pathways.
• Engineered Materials: The work provides a systematic map (encoded in the stability diagram) for designing responsive shells and soft actuators. It demonstrates that axial, circumferential, or helical deformations can be programmatically selected by tuning the nematic orientation and activity level, offering a design principle for materials where oriented active elements are embedded in deformable elastic surfaces.

The study concludes that extending this framework to non-cylindrical geometries (e.g., heart, bladder, developing brain) could uncover further non-equilibrium morphodynamic phenomena at the intersection of active matter physics and tissue mechanics.