Hydrodynamic Synchronization of Spontaneously Beating Filaments

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.

Spontaneous oscillation of an active filament under viscosity gradients

We investigate the effects of uniform viscosity gradients on the spontaneous oscillations of an elastic, active filament in viscous fluids. Combining numerical simulations and linear stability analysis, we demonstrate that a viscosity gradient increasing from the filament's base to tip destabilises the system, facilitating its self-oscillation. This effect is elucidated through a reduced-order model, highlighting the delicate balance between destabilising active forces and stabilising viscous forces. Additionally, we reveal that while a perpendicular viscosity gradient to the filament's orientation minimally affects instability, it induces asymmetric ciliary beating, thus generating a net flow along the gradient. Our findings offer new insights into the complex behaviours of biological and artificial filaments in complex fluid environments, contributing to the broader understanding of filament dynamics in heterogeneous viscous media.

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.

The reaction-diffusion basis of animated patterns in eukaryotic flagella

The flagellar beat of bull spermatozoa and C. Reinhardtii are modelled by a minimal, geometrically exact, reaction-diffusion system. Spatio-temporal animated patterns describe flagellar waves, analogous to chemical-patterns from classical reaction-diffusion systems, with sliding-controlled molecular motor reaction-kinetics. The reaction-diffusion system is derived from first principles as a consequence of the high-internal dissipation by the flagellum relative to the external hydrodynamic dissipation. Quantitative comparison with nonlinear, large-amplitude simulations shows that animated reaction-diffusion patterns account for the experimental beating of both bull sperm and C. Reinhardtii. Our results suggest that a unified mechanism may exist for motors controlled by sliding, without requiring curvature-sensing, and uninfluenced by hydrodynamics. High-internal dissipation instigates autonomous travelling waves independently of the external fluid, enabling progressive swimming, otherwise not possible, in low viscosity environments, potentially critical for external fertilizers and aquatic microorganisms. The reaction-diffusion system may prove a powerful tool for studying pattern formation of movement on animated structures.

Mobility-induced order in active XY spins on a substrate

We elucidate that the nearly phase-ordered active XY spins in contact with a conserved, diffusing species on a substrate can be stable. For wide-ranging model parameters, it has stable uniform phases robust against noises. These are distinguished by generalized quasi-long-range (QLRO) orientational order logarithmically stronger or weaker than the well-known QLRO in equilibrium, together with miniscule (i.e., hyperuniform) or giant number fluctuations, respectively. This illustrates a direct correspondence between the two. The scaling of both phase and density fluctuations in the stable phase-ordered states is nonuniversal: they depend on the nonlinear dynamical couplings. For other parameters, it has no stable uniformly ordered phase. Our model, a theory for active spinners, provides a minimal framework for wide-ranging systems, e.g., active superfluids on substrates, synchronization of oscillators, active carpets of cilia and bacterial flagella, and active membranes.

Active XY model on a substrate: Density fluctuations and phase ordering

We explore the generic long-wavelength properties of an active XY model on a substrate, consisting of a collection of nearly phase-ordered active XY spins in contact with a diffusing, conserved species, as a representative system of active spinners with a conservation law. The spins rotate actively in response to the local density fluctuations and local phase differences, on a solid substrate. We investigate this system by Monte Carlo simulations of an agent-based model, which we set up, complemented by the hydrodynamic theory for the system. We demonstrate that this system can phase-synchronize without any hydrodynamic interactions. Our combined numerical and analytical studies show that this model, when stable, displays hitherto unstudied scaling behavior: As a consequence of the interplay between the mobility, active rotation, and number conservation, such a system can be stable over a wide range of the model parameters characterized by a novel correspondence between the phase and density fluctuations. In different regions of the phase space where the phase-ordered system is stable, it displays generalized quasi-long-range order (QLRO): It shows phase ordering which is generically either logarithmically stronger than the conventional QLRO found in its equilibrium limit, together with "miniscule number fluctuations," or logarithmically weaker than QLRO along with "giant number fluctuations," showing a novel one-to-one correspondence between phase ordering and density fluctuations in the ordered states. Intriguingly, these scaling exponents are found to depend explicitly on the model parameters. We further show that in other parameter regimes there are no stable, ordered phases. Instead, two distinct types of disordered states with short-range phase order are found, characterized by the presence or absence of stable clusters of finite sizes. In a surprising connection, the hydrodynamic theory for this model also describes the fluctuations in a Kardar-Parisi-Zhang (KPZ) surface with a conserved species on it, or an active fluid membrane with a finite tension, without momentum conservation and a conserved species living on it. This implies the existence of stable fluctuating surfaces that are only logarithmically smoother or rougher than the Edward-Wilkinson surface at two dimensions (2D) can exist, in contrast to the 2D pure KPZ-like "rough" surfaces.

Hydrodynamic interactions hinder transport of flow-driven colloidal particles

The flow-driven transport of interacting micron-sized particles occurs in many soft matter systems spanning from the translocation of proteins to moving emulsions in microfluidic devices. Here we combine experiments and theory to investigate the collective transport properties of colloidal particles along a rotating ring of optical traps. In the corotating reference frame, the particles are driven by a vortex flow of the surrounding fluid. When increasing the depth of the optical potential, we observe a jamming behavior that manifests itself in a strong reduction of the current with increasing particle density. We show that this jamming is caused by hydrodynamic interactions that enhance the energetic barriers between the optical traps. This leads to a transition from an over- to an under-critical tilting of the potential in the corotating frame. Based on analytical considerations, the enhancement effect is estimated to increase with increasing particle size or decreasing radius of the ring of traps. Measurements for different ring radii and Stokesian dynamics simulations for corresponding particle sizes confirm this. The enhancement of potential barriers in the flow-driven system is contrasted to the reduction of barriers in a force-driven one. This diverse behavior demonstrates that hydrodynamic interactions can have a very different impact on the collective dynamics of many-body systems. Applications to soft matter and biological systems require careful consideration of the driving mechanism and of the role of hydrodynamic interactions.

Self-sustaining oscillation of two axonemal microtubules based on a stochastic bonding model between microtubules and dynein

The motility of cilia and flagella plays important physiological roles, and there has been a great deal of research on the mechanisms underlying the motility of molecular motors. Although recent molecular structural analyses have revealed the components of the ciliary axoneme, the mechanisms involved in the regulation of dynein activity are still unknown, and how multiple dyneins coordinate their movements remains unclear. In particular, the mode of binding for axonemal dynein has not been elucidated. In this study, we constructed a thermodynamic stochastic model of microtubule-dynein coupling and reproduced the experiments of Aoyama and Kamiya on the minimal component of axonemal microtubule-dynein. We then identified the binding mode of axonemal dynein and clarified the relationship between dynein activity distribution and axonemal movement. Based on our numerical results, the slip-bond mechanism agrees quantitatively with the experimental results in terms of amplitude, frequency, and propagation velocity, implying that axial microtubule-dynein coupling may follow a slip-bond mechanism. Moreover, the frequency and propagation velocity decayed in proportion to the fourth power of microtubule length, and the critical load of the trigger for the oscillation agreed well with Euler's critical load.

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.

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.

System size synchronization

In this paper we bring out the existence of a kind of synchronization associated with the size of a complex system. A dichotomic random jump process associated with the dynamics of an externally driven stochastic system with N coupled units is constructed. We define an output frequency and phase diffusion coefficient. System size synchronization occurs when the average output frequency is locked to the external one and the average phase diffusion coefficient shows a very deep minimum for a range of system sizes. Analytical and numerical procedures are introduced to study the phenomenon, and the results describe successfully the existence of system size synchronization.

Active droploids

Active matter comprises self-driven units, such as bacteria and synthetic microswimmers, that can spontaneously form complex patterns and assemble into functional microdevices. These processes are possible thanks to the out-of-equilibrium nature of active-matter systems, fueled by a one-way free-energy flow from the environment into the system. Here, we take the next step in the evolution of active matter by realizing a two-way coupling between active particles and their environment, where active particles act back on the environment giving rise to the formation of superstructures. In experiments and simulations we observe that, under light-illumination, colloidal particles and their near-critical environment create mutually-coupled co-evolving structures. These structures unify in the form of active superstructures featuring a droplet shape and a colloidal engine inducing self-propulsion. We call them active droploids-a portmanteau of droplet and colloids. Our results provide a pathway to create active superstructures through environmental feedback.

Energetics of synchronization for model flagella and cilia

Synchronization is often observed in the swimming of flagellated cells, either for multiple appendages on the same organism or between the flagella of nearby cells. Beating cilia are also seen to synchronize their dynamics. In 1951, Taylor showed that the observed in-phase beating of the flagella of coswimming spermatozoa was consistent with minimization of the energy dissipated in the surrounding fluid. Here we revisit Taylor's hypothesis for three models of flagella and cilia: (1) Taylor's waving sheets with both longitudinal and transverse modes, as relevant for flexible flagella, (2) spheres orbiting above a no-slip surface to model interacting flexible cilia, and (3) whirling rods above a no-slip surface to address the interaction of nodal cilia. By calculating the flow fields explicitly, we show that the rate of working of the model flagella or cilia is minimized in our three models for (1) a phase difference depending on the separation of the sheets and precise waving kinematics, (2) in-phase or opposite-phase motion depending on the relative position and orientation of the spheres, and (3) in-phase whirling of the rods. These results will be useful in future models probing the dynamics of synchronization in these setups.

Hydrodynamic synchronization in strong confinement

Cellular appendages conferring motility, such as flagella and cilia, are known to synchronise their periodic beats. The origin of synchronization is a combination of long-range hydrodynamic interactions with physical mechanisms allowing the phases of these biological oscillators to evolve. Two of such mechanisms have been identified by previous work, the elastic compliance of the periodic orbit or oscillations driven by phase-dependent biological forcing, both of which can lead generically to stable phase locking. In order to help uncover the physical mechanism for hydrodynamic synchronization most essential overall in biology, we theoretically investigate in this paper the effect of strong confinement on the effectiveness of hydrodynamic synchronization. Following past work, we use minimal models of cilia where appendages are modeled as rigid spheres forced to move along circular trajectories near a rigid surface. Strong confinement is modeled by adding a second nearby surface, parallel to the first one, where the distance between the surfaces is much smaller than the typical distance between the cilia, which results in a qualitative change in the nature of hydrodynamic interactions. We calculate separately the impact of hydrodynamic confinement on the synchronization dynamics of the elastic compliance and the force modulation mechanisms and compare our results to the usual case with a single surface. Applying our results to the biologically relevant situation of nodal cilia, we show that force modulation is a mechanism that leads to phase-locked states under strong confinement that are very similar to those without confinement as a difference with the elastic compliance mechanism. Our results point therefore to the robustness of force modulation for synchronization, an important feature for biological dynamics that therefore suggests it could be the most essential physical mechanism overall in arrays of nodal cilia. We further examine the distinct biologically relevant situation of primary cilia and show in that case that the difference in robustness of the mechanisms is not as pronounced but still favors the force modulation.

Minimal model of the hydrodynamical coupling of flagella on a spherical body with application to Volvox

Flagella are hairlike appendages attached to microorganisms that allow the organisms to traverse their fluid environment. The algae Volvox are spherical swimmers with thousands of individual flagella on their surface, and their coordination is not fully understood. In this work, a previously developed minimal model of flagella synchronization is extended to the outer surface of a sphere submerged in a fluid. Each beating flagellum tip is modeled as a small sphere, elastically bound to a circular orbit just above the spherical surface and a regularized image system for Stokes flow outside of a sphere is used to enforce the no-slip condition. Biologically relevant distributions of rotors results in a rapidly developing and robust symplectic metachronal wave traveling from the anterior to the posterior of the spherical Volvox body.

Multisynchrony in Active Microfilaments

Biological microfilaments exhibit a variety of synchronization modes. Recent experiments observed that a pair of isolated eukaryotic flagella, coupled solely via the fluid medium, display synchrony at nontrivial phase lags in addition to in-phase and antiphase synchrony. Using an elastohydrodynamic filament model in conjunction with numerical simulations and a Floquet-type theoretical analysis, we demonstrate that it is possible to reach multiple synchronization states by varying the intrinsic activity of the filament and the strength of hydrodynamic coupling between the two filaments. Then, we derive an evolution equation for the phase difference between the two filaments at weak coupling, and use a Kuramoto-style phase sensitivity analysis to reveal the nature of the bifurcations underlying the transitions between these different synchronized states.

Internal friction controls active ciliary oscillations near the instability threshold

Ciliary oscillations driven by molecular motors cause fluid motion at micron scale. Stable oscillations require a substantial source of dissipation to balance the energy input of motors. Conventionally, it stems from external fluid. We show, in contrast, that external fluid friction is negligible compared to internal elastic stress through a simultaneous measurement of motion and flow field of an isolated and active Chlamydomonas cilium beating near the instability threshold. Consequently, internal friction emerges as the sole source of dissipation for ciliary oscillations. We combine these experimental insights with theoretical modeling of active filaments to show that an instability to oscillations takes place when active stresses are strain softening and shear thinning. Together, our results reveal a counterintuitive mechanism of ciliary beating and provide a general experimental and theoretical methodology to analyze other active filaments, both biological and synthetic ones.

Human sperm uses asymmetric and anisotropic flagellar controls to regulate swimming symmetry and cell steering

Flagellar beating drives sperm through the female reproductive tract and is vital for reproduction. Flagellar waves are generated by thousands of asymmetric molecular components; yet, paradoxically, forward swimming arises via symmetric side-to-side flagellar movement. This led to the preponderance of symmetric flagellar control hypotheses. However, molecular asymmetries must still dictate the flagellum and be manifested in the beat. Here, we reconcile molecular and microscopic observations, reconnecting structure to function, by showing that human sperm uses asymmetric and anisotropic controls to swim. High-speed three-dimensional (3D) microscopy revealed two coactive transversal controls: An asymmetric traveling wave creates a one-sided stroke, and a pulsating standing wave rotates the sperm to move equally on all sides. Symmetry is thus achieved through asymmetry, creating the optical illusion of bilateral symmetry in 2D microscopy. This shows that the sperm flagellum is asymmetrically controlled and anisotropically regularized by fast-signal transduction. This enables the sperm to swim forward.