Active Filaments on Curved Surfaces: From Single Filaments to Dilute Suspensions

1. Problem Statement

The paper addresses the gap in understanding how surface curvature influences the dynamics of active semiflexible filaments (e.g., cytoskeletal filaments like actin or microtubules). While curvature is known to organize active matter in biological systems (e.g., tissue morphogenesis) and synthetic systems (e.g., microtubules on curved interfaces), previous studies were largely limited to:

• Point-particle models: Which ignore the extended nature of filaments.
• Constant curvature surfaces: Primarily spheres, which do not capture the effects of spatially heterogeneous or sign-changing curvature (positive vs. negative Gaussian curvature).
• Passive systems: Lacking the interplay between self-propulsion (activity) and elastic bending rigidity.

The authors aim to determine how spatially varying curvature interacts with filament flexibility and activity to produce new dynamical regimes, specifically focusing on the transition from single-filament behavior to collective dynamics in dilute suspensions.

2. Methodology

The study employs large-scale Langevin dynamics simulations to model active semiflexible chains confined to smooth, closed curved surfaces in the "dry limit" (ignoring long-range hydrodynamic interactions).

• Simulation Model:
  • Filaments: Modeled as chains of $N_b = 80$ beads.
  • Forces: Includes active propulsion ($f_{act}$), bonded interactions (harmonic angle potential for bending, tether potential for stretching), non-bonded steric repulsion (WCA potential), and thermal noise.
  • Confinement: Beads are projected onto the surface, with velocities and forces restricted to the local tangent plane.
• Geometries Investigated:
  • Gaussian Bump: An axisymmetric surface with a peak of positive curvature and a base of negative curvature, decaying to a flat plane.
  • Cassini Ovals: Surfaces of revolution generated by rotating Cassini ovals. These allow for a continuous deformation from a sphere ($a=0$) to a prolate spheroid, and finally to a "peanut" shape with two lobes connected by a narrow neck of negative Gaussian curvature ($a \approx 0.99b$).
• Key Parameters:
  • Activity ($f_p$): Magnitude of self-propulsion.
  • Bending Rigidity ($\kappa$): Resistance to bending.
  • Péclet Number ($Pe$): A geometry-dependent dimensionless number defined as $Pe = f_p l_s / \kappa$, where $l_s = 1/\sqrt{|K|_{max}}$ is the characteristic length scale of curvature. This quantifies the competition between active propulsion timescales and elastic deformation timescales.

3. Key Contributions

• Theoretical Framework: Extends the analysis of filament bending energy on curved surfaces using the Darboux frame, highlighting that geodesics (shortest paths) are not always energy-minimizing configurations unless specific symmetry conditions (like vanishing geodesic torsion) are met.
• New Dimensionless Parameter: Introduces a curvature-dependent Péclet number ($Pe$) to systematically map the regimes where activity dominates over elasticity versus where elasticity dominates.
• Discovery of Trapping Mechanisms: Identifies that regions of negative Gaussian curvature act as traps for active filaments in collective systems, a phenomenon not observed on constant-curvature spheres.

4. Key Results

A. Single Filament Dynamics

• Geodesic Alignment vs. Elastic Deviation:
  • High Activity / Low Rigidity (High $Pe$): Filaments closely follow geodesic paths. The active force drives the filament along the surface geometry faster than elastic forces can deform it.
  • Low Activity / High Rigidity (Low $Pe$): Filaments deviate from geodesics. The bending energy penalty prevents the filament from following the curvature of the surface, leading to trajectories that minimize elastic energy rather than following geodesics.
• Gaussian Bump: On a Gaussian bump, filaments with high rigidity deviate significantly from the initial geodesic when crossing the bump, whereas floppy filaments ($\kappa=0$) follow the geodesic perfectly.
• Cassini Ovals: Even on highly deformed "peanut" shapes, single filaments (in the absence of noise and crowding) explore the entire surface and do not get trapped in specific lobes, provided they are initialized on a geodesic that minimizes bending energy.

B. Collective Dynamics (Dilute Suspensions)

• Rotating Bands: On spherical and mildly deformed (prolate/peanut) surfaces, filaments self-organize into rotating polar bands near the equator. This mirrors experimental observations of actin on spherical vesicles.
• Curvature-Induced Trapping:
  • On highly deformed Cassini ovals ($a=0.99b$) with a narrow neck of negative curvature, a trapping transition occurs.
  • At high densities, filaments become confined to a single lobe and fail to cross the neck. This is a crowding effect: the neck is too narrow to accommodate the collective flow, and the negative curvature creates an effective barrier.
  • At low densities, filaments can collectively cross the neck, exploring the full surface.
• Activity Dependence: Increasing the Péclet number ($Pe$) at fixed high density allows filaments to overcome the trapping barrier, enabling full surface exploration. However, at high $Pe$, the motion becomes disordered, and the coherent rotating bands break down.

5. Significance and Conclusion

• Geometric Control: The study demonstrates that surface geometry is a powerful, tunable mechanism to control the organization and transport of active matter. By designing surfaces with specific curvature profiles (e.g., introducing negative curvature necks), one can induce trapping or guide collective flow.
• Biological Relevance: The findings offer insights into biological processes where curvature and activity compete, such as organelle transport, cell migration on curved tissues, and the formation of cytoskeletal structures in non-spherical cells.
• Design Principles: The results suggest that for microfluidic devices or bio-inspired materials, one can manipulate the "flocking" behavior of active agents by engineering the substrate's curvature and topology, rather than just relying on chemical gradients or external fields.

In summary, the paper establishes that while activity drives filaments toward geodesics, the interplay with bending rigidity and surface topology (specifically negative curvature regions) can lead to complex emergent behaviors, including curvature-induced trapping and the breakdown of collective order.