Synchronized oscillations, traveling waves, and jammed clusters induced by steric interactions in active filament arrays

Multiflagellate Swimming Controlled by Hydrodynamic Interactions

Many eukaryotic microorganisms propelled by multiple flagella can swim very rapidly with distinct gaits. Here, we model a three-dimensional mutiflagellate swimmer, resembling the microalgae. When the flagella are actuated synchronously, the swimming efficiency can be enhanced or reduced by interflagella hydrodynamic interactions (HIs), determined by the intrinsic tilting angle of the flagella. The asynchronous gait with a phase difference between neighboring flagella can reduce oscillatory motion via the basal mechanical coupling. In the presence of a spherical body, simulations taking into account the flagella-body interactions reveal the advantage of anterior configuration compared with posterior configuration, where in the latter case an optimal flagella number arises. Apart from understanding the role of HIs in the multiflagellate microorganisms, this work could also guide laboratory fabrications of novel microswimmers.

Mechanisms and functions of multiciliary coordination

Ciliated organisms are present in virtually every branch of the eukaryotic tree of life. In diverse systems, cilia operate in a coordinated manner to drive fluid flows, or even propel entire organisms. How do groups of motile cilia coordinate their activity within a cell or across a tissue to fulfil essential functions of life? In this review, we highlight the latest developments in our understanding of the mechanisms and functions of multiciliary coordination in diverse systems. We explore new and emerging trends in bioimaging, analytical, and computational methods, which together with their application in new model systems, have conspired to deliver important insights into one of the most fundamental questions in cellular dynamics.

Active Spaghetti: Collective Organization in Cyanobacteria

Filamentous cyanobacteria can show fascinating examples of nonequilibrium self-organization, which, however, are not well understood from a physical perspective. We investigate the motility and collective organization of colonies of these simple multicellular lifeforms. As their area density increases, linear chains of cells gliding on a substrate show a transition from an isotropic distribution to bundles of filaments arranged in a reticulate pattern. Based on our experimental observations of individual behavior and pairwise interactions, we introduce a nonreciprocal model accounting for the filaments' large aspect ratio, fluctuations in curvature, motility, and nematic interactions. This minimal model of active filaments recapitulates the observations, and rationalizes the appearance of a characteristic length scale in the system, based on the Péclet number of the cyanobacteria filaments.

Nonreciprocal interactions give rise to fast cilium synchronization in finite systems

Motile cilia beat in an asymmetric fashion in order to propel the surrounding fluid. When many cilia are located on a surface, their beating can synchronize such that their phases form metachronal waves. Here, we computationally study a model where each cilium is represented as a spherical particle, moving along a tilted trajectory with a position-dependent active driving force and a position-dependent internal drag coefficient. The model thus takes into account all the essential broken symmetries of the ciliary beat. We show that taking into account the near-field hydrodynamic interactions, the effective coupling between cilia even over an entire beating cycle can become nonreciprocal: The phase of a cilium is more strongly affected by an adjacent cilium on one side than by a cilium at the same distance in the opposite direction. As a result, synchronization starts from a seed at the edge of a group of cilia and propagates rapidly across the system, leading to a synchronization time that scales proportionally to the linear dimension of the system. We show that a ciliary carpet is characterized by three different velocities: the velocity of fluid transport, the phase velocity of metachronal waves, and the group velocity of order propagation. Unlike in systems with reciprocal coupling, boundary effects are not detrimental for synchronization, but rather enable the formation of the initial seed.

Fluid circulation driven by collectively organized metachronal waves in swimming T. aceti nematodes

Recent experiments have shown that the nematode T. aceti can assemble into collectively undulating groups at the edge of fluid drops. This coordinated state consists of metachronal waves and drives fluid circulation inside the drop. We find that the circulation velocity is about 2 mm/s and nearly half the speed of the metachronal wave. We develop a quasi-two-dimensional hydrodynamics model using the Stokes flow approximation. The periodic motion of the nematodes constitute our moving boundary condition that drives the flow. Our model suggests that large-amplitude excursions of the nematode tails produce the fluid circulation. We discuss the constraints on containers that would enhance fluid motion, which could be used in the future design of on demand flow generating systems.

Elastohydrodynamic Synchronization of Rotating Bacterial Flagella

To rotate continuously without jamming, the flagellar filaments of bacteria need to be locked in phase. While several models have been proposed for eukaryotic flagella, the synchronization of bacterial flagella is less well understood. Starting from a reduced model of flexible and hydrodynamically coupled bacterial flagella, we rigorously coarse grain the equations of motion using the method of multiple scales, and hence show that bacterial flagella generically synchronize to zero phase difference via an elastohydrodynamic mechanism. Remarkably, the far-field rate of synchronization is maximized at an intermediate value of elastic compliance, with surprising implications for bacteria.

Lateral contact yields longitudinal cohesion in active undulatory systems

Many animals and robots move using undulatory motion of their bodies. When the bodies are in close proximity undulatory motion can lead to novel collective behavior such as gait synchronization, spatial reconfiguration, and clustering. Here we study the role of contact interactions between model undulatory swimmers: three-link robots in experiment and multilink swimmers in simulation. The undulatory gait of each swimmer is generated through a time-dependent sinusoidal-like waveform which has a fixed phase offset, ϕ. By varying the phase relationship between neighboring swimmers we seek to study how contact forces and planar configurations are governed by the phase difference between neighboring swimmers. We find that undulatory actuation in close proximity drives neighboring swimmers into planar equilibrium configurations that depend on the actuation phase difference. We propose a model for stable planar configurations of nearest-neighbor undulatory swimmers which we call the gait compatibility condition, which is the set of planar and phase configurations in which no collisions occur. Robotic experiments with two, three, and four swimmers exhibit good agreement with the compatibility model. To study the contact forces and the time-averaged equilibrium between undulatory systems we perform simulations. To probe the interaction potential between undulatory swimmers we apply a small force to each swimmer longitudinally to separate them from the compatible configuration and we measure their steady-state displacement. These studies reveal that undulatory swimmers in close proximity exhibit attractive longitudinal interaction forces that drive the swimmers from incompatible to compatible configurations. This system of undulatory swimmers provides new insight into active-matter systems which move through body undulation. In addition to the importance of velocity and orientation coherence in active-matter swarms, we demonstrate that undulatory phase coherence is also important for generating stable, cohesive group configurations.

A multiscale biophysical model gives quantized metachronal waves in a lattice of beating cilia

Motile cilia are slender, hair-like cellular appendages that spontaneously oscillate under the action of internal molecular motors and are typically found in dense arrays. These active filaments coordinate their beating to generate metachronal waves that drive long-range fluid transport and locomotion. Until now, our understanding of their collective behavior largely comes from the study of minimal models that coarse grain the relevant biophysics and the hydrodynamics of slender structures. Here we build on a detailed biophysical model to elucidate the emergence of metachronal waves on millimeter scales from nanometer-scale motor activity inside individual cilia. Our study of a one-dimensional lattice of cilia in the presence of hydrodynamic and steric interactions reveals how metachronal waves are formed and maintained. We find that, in homogeneous beds of cilia, these interactions lead to multiple attracting states, all of which are characterized by an integer charge that is conserved. This even allows us to design initial conditions that lead to predictable emergent states. Finally, and very importantly, we show that, in nonuniform ciliary tissues, boundaries and inhomogeneities provide a robust route to metachronal waves.

Metachronal waves in concentrations of swimming Turbatrix aceti nematodes and an oscillator chain model for their coordinated motions

At high concentration, free swimming nematodes known as vinegar eels (Turbatrix aceti), collectively exhibit metachronal waves near a boundary. We find that the frequency of the collective traveling wave is lower than that of the freely swimming organisms. We explore models based on a chain of oscillators with nearest-neighbor interactions that inhibit oscillator phase velocity. The phase of each oscillator represents the phase of the motion of the eel's head back and forth about its mean position. A strongly interacting directed chain model mimicking steric repulsion between organisms robustly gives traveling wave states and can approximately match the observed wavelength and oscillation frequency of the observed traveling wave. We predict body shapes assuming that waves propagate down the eel body at a constant speed. The phase oscillator model that impedes eel head overlaps also reduces close interactions throughout the eel bodies.

Three-dimensional nonlinear dynamics of prestressed active filaments: Flapping, swirling, and flipping

Initially straight slender elastic filaments or rods with constrained ends buckle and form stable two-dimensional shapes when prestressed by bringing the ends together. Beyond a critical value of this prestress, rods can also deform off plane and form twisted three-dimensional equilibrium shapes. Here, we analyze the three-dimensional instabilities and dynamics of such deformed filaments subject to nonconservative active follower forces and fluid drag. We find that softly constrained filaments that are clamped at one end and pinned at the other exhibit stable two-dimensional planar flapping oscillations when active forces are directed toward the clamped end. Reversing the directionality of the forces quenches the instability. For strongly constrained filaments with both ends clamped, computations reveal an instability arising from the twist-bend-activity coupling. Planar oscillations are destabilized by off-planar perturbations resulting in twisted three-dimensional swirling patterns interspersed with periodic flipping or reversal of the swirling direction. These striking swirl-flip transitions are characterized by two distinct timescales: the time period for a swirl (rotation) and the time between flipping events. We interpret these reversals as relaxation oscillation events driven by accumulation of torsional energy. Each cycle is initiated by a fast jump in torsional deformation with a subsequent slow decrease in net torsion until the next cycle. Our work reveals the rich tapestry of spatiotemporal patterns when weakly inertial strongly damped rods are deformed by nonconservative active forces. Taken together, our results suggest avenues by which prestress, elasticity, and activity may be used to design synthetic macroscale pumps or mixers.