Emergence of metachronal waves in a chain of symmetrically beating filaments
Narina Jung
Korea Institute for Advanced Study, Seoul 02455, Korea
Won Kyu Kim
Korea Institute for Advanced Study, Seoul 02455, Korea
Changbong Hyeon
Korea Institute for Advanced Study, Seoul 02455, Korea
(June 27, 2024)
Abstract
Recent experiments have shown that metachronal waves (MCWs) can emerge from a chain of symmetrically beating nematodes aligned at the edge of sessile droplets.
Our study, employing a coupled elastohydrodynamic model of active filaments, elucidates that a misalignment caused by a tilt against the bounding wall disrupts the synchronization and generates a constant time lag between adjacent filaments, leading to MCWs.
The MCWs, enhancing the fluid circulation, achieve their maximum thermodynamic efficiency over the same range of tilt angles observed in the nematode experiments.
Introduction.
Collective dynamics observed at the scale of animal population are an outcome of local dynamics and coupling between individual agents
[1, 2, 3, 4, 5, 6].
Among them, long wavelength modulations of collective beating behavior called metachronal waves (MCWs), characterized by a constant phase or time lag between the dynamics of adjacent agents,
enhances the fluid flow and transport of nutrients essential for a host of microorganisms.
As a result, MCW has long been a topic of abiding interest [7, 8, 9, 10, 11, 12, 13].
Since the Taylor’s seminal works [7, 8],
long-range hydrodynamic coupling has widely been appreciated as the key interaction for the coordination and generation of MCWs among filamentous micro-swimmers in low Reynolds number environment [14, 15, 16, 17, 9, 10, 18].
Significance of short-range steric interactions was also highlighted in the gait synchronization and generation of collective waves with constant phase lags among C. elegans [19, 20].
Applying transverse external flow to symmetrically undulating filaments,
Guirao and Joanny have suggested that breaking the left-right (LR) symmetry in the filament dynamics as well as including the hydrodynamic interactions (HIs) is necessary to generate MCWs [21].
Peshkov et al. have recently demonstrated formation of MCWs in a chain of Turbatrix aceti (T. aceti) aligned along the contact line of a water droplet [22, 23].
In more broadly discussed MCWs in cilia, individual cilium exhibits two-phase asymmetric beating consisting of power and recovery strokes.
In contrast, despite its functional importance, the MCW emerging from symmetrically beating nematodes is physically not obvious and its origin remains unclear.
To elucidate how MCWs arise from the 1D array of nematodes,
we extend a phenomenological elastohydrodynamic model of a single active filament [24] to many filaments, and study the physics underlying the formation of MCWs and its effect on the flow of surrounding fluid.
Theoretical model. We consider that a 1D array of slender filaments with radius and length are aligned along the -axis, separated by (), with a tilt angle against -axis (see Fig. 1).
The force-balance equation for th filament coupled with its neighbors and the bounding wall is written in a non-dimensionalized form (see Supplemental Materials for details):
(1)
with the boundary conditions, and
.
Here is the transversal displacement (parallel to the -axis) at position along the centerline of th filament and at time .
In ,
the three terms in the parenthesis model the effect of an active force originating from a mechanical feedback of molecular motors [25, 26].
The parameter in the first term, denoting a dimensionless speed of wave along the centerline, reflects the filament activity.
To make the hydrodynamic forces, , acting on the th filament more tractable and interpretable, we consider the resistive force theory [27, 24], and simplify to
(2)
where denotes the position along the th filament.
The coefficient for the strength of HI is given by
(3)
with the first and the second terms from the Stokeslet and doublet, respectively [28, 29].
In our system with and , the coefficient is .
The terms and
model the repulsion between the filaments and the bounding wall effect on the filaments.
Inclusion of SIs in Eq. 1 is necessary in stabilizing
the in-phase synchronized state against otherwise the equally stable anti-phase synchronized state under hydrodynamic coupling, especially for an array of filaments satisfying the condition where (see SM Text and Fig. S1 for details).
We study the collective dynamics of by numerically solving Eq.1 with varying and tilt angle [30].
Figure 1: Coordinate system in two-dimension used to study the beating dynamics of filaments aligned in parallel to the -axis with a finite tilt angle towards the boundary line ().
The coordinate of the th filament in the reference frame is associated with the local coordinates with as
where is the size of mismatch.
The head position of the first filament () is set to the origin .
In what follows, we first characterize the single filament dynamics and next explore how the inter-filament coupling modulates the symmetrically beating dynamics of each filament
and contributes to the formation of MCWs.
Single filament dynamics.
For a single filament () without tilt (), the wall interaction (WI) is minimal, which reduces
Eq. 1 to
, and yields
a periodic solution of for , where the subscript denotes the waveform solution for a single filament () at which we use as a reference waveform throughout the paper (Fig. 2(b), see also SM Movie 1).
One of the key features of is that it preserves LR symmetry
such that both and are the solution of Eq.1,
while displaying a growing amplitude from the head to tail (Fig. 2(a)).
The -dependent beat period, ,
can either be read off from the kymograph (Fig. 2(b)) that satisfies ,
or be obtained as
from the power spectrum
.
The spectrum (Fig. 2(c))
shows a stripe at , effectively independent of .
Similarly, the wavelength at time is either estimated through or is
obtained as
from the power spectrum
(Fig. 2(d)).
The mean wavelength of a filament for is
with .
We confirm that the wave speed along the filament body calculated by using and is identical to the input wave speed (), i.e., .
Figure 2: Beat dynamics of a single filament at (see SM Movie 1).
(a) Filament configurations generated at in increasing time with the interval of (pale to dark).
(b) Kymograph () depicting the spatiotemporal evolution of filament dynamics.
(c), (d) Power spectra of in time and space:
(c) , and (d)
.
Figure 3: Dynamics of two interacting filaments with increasing mismatch.
(a) Configurations of the filaments (#1 (blue) and #2 (red)) and their power spectra, and , calculated for
, , and , and .
(b) Mean undulation period of each filament, , as a function of . See also SM Movie 2.
Dynamics of interacting filaments.
When multiple filaments interact with each other, additional force terms, , , and , come into play.
Based on the numerical solutions of Eq.1 at steady states, we posit that the presence of SI and WI terms renormalizes the coefficient for HIs is modified from to .
When the array of filaments are tilted at a finite angle
towards the bounding wall (Fig. 1), the filaments do not fully overlap, engendering a mismatch, whose size
is related to as (Fig. 1).
Using the resistive force theory that enables one to replace the hydrodynamic coupling with the velocity of transversal displacement (Eq. 2), the force-balance equation for each filament in the array made of filaments can be written as
(4)
The segment of th filament located at
interacts with the segment of th filament at .
When the filaments beat collectively with ,
the transversal displacements of the interacting segments satisfy the inequalities,
.
Thus, we assume that the transversal displacements of the neighboring filaments satisfy the relation , where
are the factors modulating the transversal displacement at .
In full synchrony without tilt (),
the filaments beat collectively, satisfying , and do not alter the waveform of adjacent filament (). Furthermore,
HIs become the most dominant among the interactions,
.
By summing the equations () in Eq. 4 and dividing it by we obtain an evolution equation for ,
.
Thus, one can consider either that is mapped to by rescaling the time as
, or that reduction of the drag coefficient occurs from to
.
This argument can be generalized to the case with a finite tilt (), enabling us to obtain the mean period of undulation for an array of filaments,
(5)
with .
For a pair of coupled filaments, we study the effect of misalignment between the adjacent filaments on the synchronization.
As the extent of mismatch () increases, three different dynamical behaviors emerge: (i) full synchronization; (ii) partial synchronization; (iii) decoupled dynamics.
Fig. 3(a) displays typical configurations of the filaments and their power spectra for three cases of .
(i) Sync regime ():
The beat dynamics of filaments are fully synchronized in the region of two filaments being overlapped, and there is a minor distinction in the pattern of power spectra between and . The corresponding mean period of undulation increases from to (see Fig. 3(b)).
This slight increase in the period is due to the reduced size of overlap between the two filaments.
In fact, for a filament characterized with a fixed wave speed,
both and tend to increase with decreasing filament size () [24].
(ii) Mixed regime ():
The tail part of the filament #1, whose transversal fluctuations is greater than that of head part, interacts with the head part of the filament #2, which engenders a mismatch in the velocity profiles, and .
Since the strength of HIs is proportional to the amplitude of transversal displacement for a fixed oscillatory period, i.e.,
(Eq. 2),
the mismatch leads to the left-right symmetry breaking of hydrodynamic forces, .
In this regime,
the synchronization and decoupling of filament dynamics coexist, and the power spectra and , exhibiting two main peaks at and , clearly differ from each other, suggestive of the broken LR symmetry in the beat dynamics of the two.
Figure 4:
MCWs generated from the 1D array of filaments at varying with periodic boundary (see SM Movie 3).
(a) Configurations of many filaments producing MCWs with increasing and the corresponding power spectra, and .
(b) Wavelength, period, and wave speed for the individual filaments in the array with increasing .
(c) The mean transversal fluctuations of filaments in the array with increasing : , and , where is the mean squared positive (or negative) fluctuations.
(left) The large amplitudes of at and is highlighted by the red arrows. (right) The ratio of left/right amplitudes, .
(d) Time evolution of the mass densities of the filaments in arrays projected on the -axis.
The speed of MCWs can be estimated from the slope of stripes: , and .
(iii) Decoupled regime (): The peak of is uniquely formed
at .
Although 40 of the filament body is subject to HIs for , the dynamics of two filaments are effectively decoupled, and the beat dynamics of each filament is restored to .
An increasing mismatch () weakens the strength of HIs (or reduces ), and it consequently tends to increase the mean period of undulation from the case with zero mismatch
to the case with full decoupling, i.e., .
As increases from the regime (i) to (iii),
the relative weight of the two main peaks in Fourier space ()
shifts from to and undergoes a first order type discrete transition (Fig. 3(a)).
Figure 5:
(a) Work production and dissipation, and (b) thermodynamic efficiency with varying . The inset shows a snapshot of nematode array forming from the Peshkov et al.’s experiment [22].
(c) The mean velocity of fluid flow calculated for an array of filaments () at .
We now analyze behaviors of 1D arrays of many filaments () at a fixed separation () by imposing a periodic boundary condition to eliminate the boundary effect (see SM Text and Fig. S2 for the case with varying without the periodic boundary).
First, for ,
The filaments beat in perfect in-phase synchrony at steady states, yielding the
waveform,
, whose wavelength and transversal fluctuation are identical to those of a single filament in isolation ( and ), yet the time is rescaled as
.
the oscillatory period is reduced to , which substantially accelerates the wave speed to .
This is in excellent agreement with a theoretical estimate, and , where we have used Eq. 5 with , (Eq. 3) and .
Individual filaments beat without phase or time lag between adjacent filaments and retain the LR symmetry.
Next, for , the waveforms for all the filaments are modified from to , while the individual filaments collectively undulate with a constant time lag against the adjacent filaments, generating MCWs that travel across the array of filaments (-axis), which differ from the body waves traveling along the centerline (-axis) of each filament at speed .
, calculated by analyzing ,
displays a clear increase with ( in Figs. 4(a) and 4(b)).
Although the wavelength () changes non-monotonically, maximized at ,
the speed of body wave ()
decreases with (Fig. 4(b)).
In light of the two-filament dynamics (Fig. 3(b)),
the cases with
and (Fig. 4(a)) belong to the mixed regime, and
the deviation of the wave characteristics from those of a single filament is evident (Fig. 4(b)).
Due to the increasing mismatch between the two filaments with , the average amplitude as well as changes non-monotonically from to , showing a tendency of recovering at
(Fig. 4(c)).
Importantly, for , the broken LR symmetry is apparent, such that the transversal displacement of filaments is significantly skewed towards the left side of the centerline, i.e., (see the red arrows in Fig. 4(c)).
Accordingly, the hydrodynamic coupling becomes stronger on the left side than on the right, breaking the LR symmetry of HIs.
The collective waveform (MCW), composed of the individual body waves, propagates along the array with a constant time lag between adjacent filaments.
The is
the time required for a position of a filament to travel from to and mimic the movement of the adjacent filament, whose head positions are separated by . Note that is related to as .
For , the speed of the MCW, , relates to as follows.
(6)
For different tilt angles, the speeds of MCWs are:
,
,
and
.
Notably, for all cases of .
Alternatively, the speed of the MCW can be determined by analyzing the spatiotemporal evolution of mass density projected on the axis,
.
As clearly shown in Fig. 4(d), is periodic in space and time as and , which enables to estimate the speed of MCW by using .
The estimated speeds,
, and ,
agree well with the values of obtained using Eq. 6.
Viscous dissipation, work production, and thermodynamic efficiency.
We investigate the thermodynamics associated with the generation of MCWs and the fluid flow.
First, given that the swimming velocity along the - and -directions resulting from the propulsion force of a filament are
and [31, 32],
the velocities along the - and -directions can be expressed as
and .
Then, the mean work production rate per filament averaged over a period, which
enhances the fluid flow along the -direction, is
.
Meanwhile, the mean viscous dissipation per filament is evaluated as .
Note that both and are maximized at (Fig. 5(a)).
On the other hand, the extent of the mechanical work transmitted to the surrounding fluid relative to the total energy consumption, quantified by calculating , is maximized in the range of (Fig. 5(b)).
Remarkably, we find that the array of T. aceti in the experiment [22] adopt effectively the same range of tilt angle observed in the experiment (, see the inset of Fig. 5(b)).
This suggests that the nematode population immersed in the fluid media is spontaneously adjusted to a configuration that maximizes the work production for a given energy resource, namely, the thermodynamic efficiency ().
The fluid flow generated by the collective motion of filaments is calculated by using the linear response between the propulsion force and the fluid velocity associated through a hydrodynamic tensor (Eq. S5, see SM Text for the details).
Fig. 5(c) shows the velocity fields of medium along -axis, generated by an array of beating filaments at , as a function of distance calculated by using a viscosity value corresponding to 1 cP.
It is noteworthy that the velocity field decreases with the distance from the -axis as for ,
reflecting a -decay characteristic of the fundamental solution of the hydrodynamic tensor (Eq. S5).
Concluding Remarks.
One of the key features of MCWs observed in the experiment on a chain of T. aceti is the modulation of the waveform of individual filaments [22].
To explore MCWs, previous studies
adopted a phase reduction approach, representing a filament dynamics by means of a frequency-locked phase oscillator
[33, 34, 35, 23, 36].
Such a strategy may only prove effective for perfectly aligned filaments whose waveform remains unaltered in synchrony (see Fig. S2) [37, 38, 39]; however, it overlooks
the physical reality that the wave characteristics of each filament is subject to change upon interacting with its neighbors and generating MCWs.
Our study makes it clear that the full evolution of waveforms is essential to evaluate the thermodynamics of a system accurately as well as to decipher the physics that gives rise to MCWs.
Unlike the asymmetrically beating dynamics of cilium, the dynamics of a single nematode is symmetric with respect to its centerline. Despite its biological importance of facilitating the fluid circulation, generation of MCWs from the nematode population appears unlikely at first sight.
Our study clarifies that a misalignment between adjacent filaments disrupts the synchronization and generates a constant time lag, which is the essential mechanism for the formation of MCWs.
Importantly, a large tilt () introduced to the array of filaments generates the energetically most efficient MCWs enhancing the circulation of surrounding fluid.
Acknowledgements.
This study was supported by the KIAS individual Grants, CG085701 (NJ), CG076002 (WKK), and CG035003 (CH) at Korea Institute for Advanced Study.
We thank the Center for Advanced Computation in KIAS for providing the computing resources.
Elgeti et al. [2015]J. Elgeti, R. G. Winkler, and G. Gompper, Physics of
microswimmers—single particle motion and collective behavior: a review, Rep. Prog. Phys. 78, 056601 (2015).
Cavagna and Giardina [2014]A. Cavagna and I. Giardina, Bird flocks as condensed
matter, Annu.
Rev. Condens. Matter Phys. 5, 183 (2014).
Moussaïd et al. [2011]M. Moussaïd, D. Helbing, and G. Theraulaz, How simple rules
determine pedestrian behavior and crowd disasters, Proc. Natl. Acad. Sci. U. S. A. 108, 6884 (2011).
Lee et al. [2018]S. Lee, C. Hyeon, and J. Jo, Thermodynamic uncertainty relation of interacting
oscillators in synchrony, Phys. Rev. E 98, 032119 (2018).
Hong et al. [2020]H. Hong, J. Jo, C. Hyeon, and H. Park, Thermodynamic cost of synchronizing a population of beating cilia, J. Stat. Mech. , 074001 (2020).
Taylor [1951]G. I. Taylor, Analysis of the swimming
of microscopic organisms, Proc. R. Soc. A: Math. Phys. Eng. Sci. 209, 447 (1951).
Taylor [1952]G. I. Taylor, Analysis of the swimming
of long and narrow animals, Proc. R. Soc. A: Math. Phys. Eng. Sci. 214, 158 (1952).
Gueron et al. [1997]S. Gueron, K. Levit-Gurevich, N. Liron, and J. J. Blum, Cilia internal mechanism and
metachronal coordination as the result of hydrodynamical coupling, Proc. Natl. Acad.
Sci. 94, 6001 (1997).
Gueron and Levit-Gurevich [1999]S. Gueron and K. Levit-Gurevich, Energetic
considerations of ciliary beating and the advantage of metachronal
coordination, Proc. Natl. Acad. Sci. U. S. A. 96, 12240 (1999).
Niedermayer et al. [2008]T. Niedermayer, B. Eckhardt, and P. Lenz, Synchronization, phase
locking, and metachronal wave formation in ciliary chains, Chaos 18 (2008).
Gu et al. [2020]H. Gu, Q. Boehler,
H. Cui, E. Secchi, G. Savorana, C. De Marco, S. Gervasoni, Q. Peyron, T.-Y. Huang, S. Pane, et al., Magnetic
cilia carpets with programmable metachronal waves, Nature Commun. 11, 2637 (2020).
Wang et al. [2024]T. Wang, T. ul Islam,
E. Steur, T. Homan, I. Aggarwal, P. R. Onck, J. M. den Toonder, and Y. Wang, Programmable metachronal motion of closely packed magnetic
artificial cilia, Lab on a Chip 24, 1573
(2024).
Hancock [1953]G. Hancock, The self-propulsion of
microscopic organisms through liquids, Proc. R. Soc. A: Math. Phys. Eng. Sci. 217, 96 (1953).
Yang et al. [2008]Y. Yang, J. Elgeti, and G. Gompper, Cooperation of sperm in two dimensions:
synchronization, attraction, and aggregation through hydrodynamic
interactions, Phys. Rev. E. 78, 061903 (2008).
Golestanian et al. [2011]R. Golestanian, J. M. Yeomans, and N. Uchida, Hydrodynamic
synchronization at low Reynolds number, Soft Matter 7, 3074 (2011).
Kim and Netz [2006]Y. W. Kim and R. R. Netz, Pumping fluids with periodically
beating grafted elastic filaments, Phys. Rev. Lett. 96, 158101 (2006).
Zheng et al. [2023]E. Zheng, M. Brandenbourger, L. Robinet, P. Schall,
E. Lerner, and C. Coulais, Self-oscillation and synchronization transitions in
elastoactive structures, Phys. Rev. Lett. 130, 178202 (2023).
Yuan et al. [2014]J. Yuan, D. M. Raizen, and H. H. Bau, Gait synchronization in
Caenorhabditis elegans, Proc. Natl. Acad. Sci. U. S. A. 111, 6865 (2014).
Chelakkot et al. [2021]R. Chelakkot, M. F. Hagan, and A. Gopinath, Synchronized
oscillations, traveling waves, and jammed clusters induced by steric
interactions in active filament arrays, Soft Matter 17, 1091 (2021).
Guirao and Joanny [2007]B. Guirao and J.-F. Joanny, Spontaneous creation of
macroscopic flow and metachronal waves in an array of cilia, Biophys. J. 92, 1900 (2007).
Peshkov et al. [2022]A. Peshkov, S. McGaffigan, and A. C. Quillen, Synchronized
oscillations in swarms of nematode Turbatrix aceti, Soft Matter 18, 1174 (2022).
Quillen et al. [2021]A. Quillen, A. Peshkov,
E. Wright, and S. McGaffigan, Metachronal waves in concentrations of swimming
Turbatrix aceti nematodes and an oscillator chain model for their
coordinated motions, Phys. Rev. E. 104, 014412 (2021).
Goldstein et al. [2016]R. E. Goldstein, E. Lauga,
A. I. Pesci, and M. R. Proctor, Elastohydrodynamic synchronization of
adjacent beating flagella, Phys. Rev. Fluids 1, 073201 (2016).
Heuser et al. [2009]T. Heuser, M. Raytchev,
J. Krell, M. E. Porter, and D. Nicastro, The dynein regulatory complex is the nexin link and a
major regulatory node in cilia and flagella, J. Cell Biol. 187, 921 (2009).
Jülicher and Prost [1997]F. Jülicher and J. Prost, Spontaneous oscillations of
collective molecular motors, Phys. Rev. Lett. 78, 4510 (1997).
Gray and Hancock [1955]J. Gray and G. Hancock, The propulsion of sea-urchin
spermatozoa, J.
Exp. Biol. 32, 802
(1955).
Rotne and Prager [1969]J. Rotne and S. Prager, Variational treatment of hydrodynamic
interaction in polymers, J. Chem. Phys. 50, 4831 (1969).
Yamakawa [1970]H. Yamakawa, Transport properties of
polymer chains in dilute solution: hydrodynamic interaction, J. Chem. Phys. 53, 436 (1970).
Tornberg and Shelley [2004]A.-K. Tornberg and M. J. Shelley, Simulating the dynamics
and interactions of flexible fibers in stokes flows, J. Comp. Phys. 196, 8 (2004).
Lauga and Powers [2009]E. Lauga and T. R. Powers, The hydrodynamics of
swimming microorganisms, Rep. Prog. Phys. 72, 096601 (2009).
Wiggins and Goldstein [1998]C. H. Wiggins and R. E. Goldstein, Flexive and propulsive
dynamics of elastica at low Reynolds number, Phys. Rev. Lett. 80, 3879 (1998).
Nakao [2016]H. Nakao, Phase reduction approach to
synchronisation of nonlinear oscillators, Contemporary Physics 57, 188 (2016).
Kawamura and Tsubaki [2018]Y. Kawamura and R. Tsubaki, Phase reduction approach
to elastohydrodynamic synchronization of beating flagella, Phys. Rev. E. 97, 022212 (2018).
Man and Kanso [2020]Y. Man and E. Kanso, Multisynchrony in active
microfilaments, Phys. Rev. Lett. 125, 148101 (2020).
Quillen [2023]A. Quillen, Robust formation of
metachronal waves in directional chains of phase oscillators, Phys. Rev. E 107, 034401 (2023).
Guo et al. [2018]H. Guo, L. Fauci, M. Shelley, and E. Kanso, Bistability in the synchronization of actuated microfilaments, J. Fluid Mech. 836, 304 (2018).
Young [2009]Y.-N. Young, Hydrodynamic interactions
between two semiflexible inextensible filaments in stokes flow, Phys. Rev. E. 79, 046317 (2009).
Vilfan and Jülicher [2006]A. Vilfan and F. Jülicher, Hydrodynamic flow
patterns and synchronization of beating cilia, Phys. Rev. Lett. 96, 058102 (2006).
Lighthill [1975]S. J. Lighthill, Mathematical biofluid
dynamics (SIAM, 1975).
Chakrabarti and Saintillan [2019]B. Chakrabarti and D. Saintillan, Spontaneous
oscillations, beating patterns, and hydrodynamics of active microfilaments, Phys. Rev.
Fluids 4, 043102
(2019).
Fang-Yen et al. [2010]C. Fang-Yen, M. Wyart,
J. Xie, R. Kawai, T. Kodger, S. Chen, Q. Wen, and A. D. Samuel, Biomechanical analysis of
gait adaptation in the nematode Caenorhabditis elegans, Proc. Natl. Acad. Sci. U. S. A. 107, 20323 (2010).
Sznitman et al. [2010]J. Sznitman, P. K. Purohit, P. Krajacic,
T. Lamitina, and P. E. Arratia, Material properties of Caenorhabditis
elegans swimming at low Reynolds number, Biophys. J. 98, 617 (2010).
Xu et al. [2016]G. Xu, K. S. Wilson,
R. J. Okamoto, J.-Y. Shao, S. K. Dutcher, and P. V. Bayly, Flexural rigidity and shear stiffness of flagella estimated from
induced bends and counterbends, Biophys. J. 110, 2759 (2016).
Gilpin et al. [2015]W. Gilpin, S. Uppaluri, and C. P. Brangwynne, Worms under pressure: bulk mechanical
properties of C. elegans are independent of the cuticle, Biophys. J. 108, 1887 (2015).
Elfring and Lauga [2009]G. J. Elfring and E. Lauga, Hydrodynamic phase locking
of swimming microorganisms, Phys. Rev. Lett. 103, 088101 (2009).
Leptos et al. [2013]K. C. Leptos, K. Y. Wan,
M. Polin, I. Tuval, A. I. Pesci, and R. E. Goldstein, Antiphase synchronization in a flagellar-dominance mutant
of chlamydomonas, Phys. Rev. Lett. 111, 158101 (2013).
I Supplemental Materials
II Supplemental Movies
•
SM Movie 1: Beat dynamics of a single filament at .
•
SM Movie 2: Dynamic of two interacting filaments at with varying sizes of mismatch .
•
SM Movie 3: Dynamics of filaments at
, , , and .
•
SM Movie 4: Demonstration of the effect of SI on two-filament dynamics starting from a V-shape configuration.
•
SM Movie 5: Many filament dynamics at with varying , without periodic boundary condition.
•
SM Movie 6: Many filament dynamics at with varying , without periodic boundary condition.
III Elastohydrodynamic model
We consider filaments, each with a length of and a cross-sectional radius of , tilted with while maintaining a body-to-body separation distance within a two-dimensional plane (Fig. 1).
The head of each filament is hinged along the -axis. The head-to-head distance depends on the tilt angle and is separated by for the given separation .
An aspect ratio of the filament is and the scale of used in the model is in the range of .
Under an assumption of small undulation compared to the body length, the curvilinear coordinate along the filament () is approximated to .
At low Reynolds number, the force balance on the th filament should be established among viscous drag (), bending (), active forces (), hydrodynamic couplings () and steric interactions () from the neighboring filaments, such that the sum of all the force components should vanish, .
Each force term can be related to the waveform of filaments .
The force associated with viscous drag is given as
(S1)
where is the drag coefficient for a slender body [40, 31].
The force arising from the bending penalty can be obtained from the functional derivative of the bending energy ,
(S2)
where is the coefficient of bending stiffness.
For the active force we adopt the phenomenological model by Goldstein et. al [24]
(S3)
where the first term () generates propagating waves along the centerline of filament (-axis in Fig. 1).
The second and the third term break head-tail symmetry of the amplitude of transversal undulations, which is magnified by the second term () and suppressed by the third term (), saturates to a finite value.
The hydrodynamic force from th to th filaments is defined via the fundamental singular solutions of Stokes equations:
(S4)
where , is the sum of Stokeslet and doublet
(S5)
with r denoting the relative position vector from to and .
The second term, corresponding to the doublet, corrects the boundary condition when the velocity is computed on the surface of the slender body near the singularity distribution. In Eq.S4, the term associated with viscous force of -th filament, , is also obtained from the resistive force theory,
.
It should be emphasized that for a non-zero tilt (),
the nearest point of adjacent chain to is for .
Taken together,
Eq.S4 is simplified to the following form [24]:
(S6)
For the steric repulsion between the adjacent filaments, we use the Weeks-Chandler-Anderson (WCA) potential
(S7)
where is a length scale at which and is a scale of interaction energy.
The repulsive force due to the steric interaction is obtained from
(S8)
For a single active elastic filament demonstrating a beating dynamics,
two dimensionless parameters can characterize its dynamics: (i) activity number, the active force relative to the bending force; (ii) sperm number, the beat frequency relative to the rate of the bending wave relaxation in a viscous fluid [32, 41].
A waveform generated from a collective dynamics of filaments in a dense suspension is also characterized by the two parameters of a single filament [15].
In our system, the activity and sperm numbers are given as
, .
Here, we assume that the characteristic time scale of active motion, , is the same as that of bending relaxation time scale, , and use it as a characteristic time scale of the system, .
The following expression that considers the WCA potential is numerically computed:
(S9)
for the steric repulsion between the neighboring filaments,
and
(S10)
for the wall repulsion,
where , , , and
.
The boundary conditions were imposed to create a hinged head () and a free tail ( and ).
As initial conditions, we used straight filaments in a parallel alignment, and studied the dynamics for various .
The governing equation (Eq.S11) was integrated using the finite difference and the semi-implicit scheme [30].
Using the relations, , ,
, and , which determine
the characteristic length and time scales of the system as
and ,
we non-dimensionalize the force balance equation for each filament as
(S11)
where , , , and
is a dimensionless speed of wave.
For the parameter values of Pasec, N [22], and Nmm2 (corresponding to the Young’s modulus of kPa) [42, 43, 44, 45], we obtain mm and
ms.
When , N/mm. We also set the values of and .
In Eq. 1, the hat symbol for non-dimensional length ( and ) and time () is dropped from the expression for simplicity.
The stability analysis on the linearized elastohydrodynamic equation,
with yields a dispersion relation,
(S12)
Thus, for the dynamics of the filament to be oscillatory,
the condition of imposes
the wavelength () of a single filament always greater than ,
(S13)
IV Effect of steric interactions on the synchronization of two hydrodynamically coupled filaments
The effect of SI on the two hydrodynamically coupled filaments are investigated for a system without tilt (), separated by .
To quantify the process of synchronization between two active filaments ( and ), we calculate a correlation function,
(S14)
Here, for perfect synchrony and when their dynamics are decoupled.
Figure S1: Dynamics of a pair of hydrodynamically-coupled filaments from V-shape initial configuration with and without SI (see SM Movie 4).
(a) Evolution of correlation functions () for the two filaments coupled either only through HIs (black) or through HIs and SIs (red).
Depicted in the insets are the steady-state configurations in (top) IP and (bottom) AP synchrony.
(b) The kymographs of two filaments ( and ) coupled through HIs and SIs.
(c) Phase diagrams of in-phase (IP) and anti-phase (AP) synchronization between two active filaments from the V-shaped initial condition for various values of and .
The steady state at given and is determined by using (Eq. S14), such that for the IP state, and
for the AP state.
The solid line in magenta plots the -dependent fluctuations of a single filament (, Eq. (S15)).
In the absence of SIs, depending on the initial condition,
the two hydrodynamically coupled filaments demonstrate either in-phase (IP) or anti-phase (AP) synchronization after a transient time [35].
Adding the short-range SI little affects the hydrodynamically coupled filaments already in IP synchrony; however, the same SI alters the two-filament configuration in AP synchrony into the one in IP synchrony, particularly, when the amplitude of undulation is greater than the inter-filament separation ().
Two-filament dynamics starting from a V-shape configuration demonstrates the effect of SI dramatically (see SM Movie 4).
The filaments, coupled only through HI,
converge to a configuration in AP synchrony (the black line in Fig. S1(a)).
However, an introduction of SI destabilizes the AP configuration, nudging the system towards the one in IP synchrony (the red graph in Fig. S1A) [46, 47, 35].
The kymographs of and depict the transition from AP () to IP () at , accompanied with the decrease of beat period from to (Fig. S1(b)).
The phase diagram of the steady state configurations of two filaments
(Fig. S1(c)) shows the stability of AP and IP modes with respect to and
[37, 24], along with their phase boundary which can be described using the -dependent amplitude of transversal displacement, , obtained
by assuming that a single traveling wave for the tail part is given as
.
Since should satisfy Eq. 1, we obtain the amplitude of transversal displacement [24]:
Figure S2: Dynamics among filaments at and with varying .
(a) Configuration of filaments in in-phase synchrony.
(b) , , calculated for filament.
For the demonstration, the filament in the middle is centered.
(c) Mean period of filament dynamics in synchrony, as a function of .
The red line is the fit using Eq. 5.
V Many filaments with non-periodic boundary conditions
We consider a case where
a group of active filaments with non-periodic boundary conditions to study the impact of the number of filaments on the collective dynamics at two tilt angles, and .
For a set of -filaments with (see Fig. S2(a) for a configuration of filaments), the filaments beat in perfect synchrony, so that the mean periods () and wavelengths () are identical for different filaments except for those at the ends in the case of and , whereas the amplitude of transversal displacement () is the largest for the filaments in the middle (Fig. S2(b)).
Notably, the wavelength of the waveform is effectively identical as , irrespective of ; however,
the mean period () decreases with (Fig. S2(c)). Thus, the wave propagation speed () along the centerline () increases with (see SM Movie 4).
Figure S3: Dynamics among filaments at and with varying .
(a) Configuration of N = 20 filaments in synchrony.
(b) Same as Fig. S2(b).
In fact, Eq. 5 with the fit parameters, , , and , well describes the periods of collective motion of filaments () for varying .
Note, however, that in the collective wave generated from the filaments, there is neither the phase lag between the adjacent filaments nor the significant change in the waveform.
The waveform of each filament bears the LR symmetry with no wave packet on average being transmitted to left or right along the -axis (see SM Movie 5).
In contrast, for , it is clearly seen that MCWs travel from left to right along the -axis, although the beat dynamics of filaments are only partially synchronized (see SM Movie 6).
For , the wave characteristics of the filaments for the two different tilt angles averaged over all the filaments excluding those contributing at the boundary are given as follows:
(i) and for ;
(ii) and for .
Thus, the wave propagation speeds along the centerline
are and .
Note that is greater than by 10 folds.
The transversal fluctuations, , are also greater as well for .
It is of note that more significant reduction is made in with increasing at than at , which is primarily due to the enhancement of HI in multiple filaments in full synchrony at in comparison with those in partial synchrony at .
VI Propulsion Force
With the resistive force theory [31], we compute the anisotropic drag forces on the th filament in the perpendicular and parallel directions with respect to the centerline of the filament body as follows.
(S16)
where is the angle formed between the filament orientation at and the centerline (-axis) and .
Then, the drag forces along the - and -axis are calculated as
(S17)
where we have used the condition of small undulation in comparison with
the body length (, , and ) to obtain the approximate expressions for the drag forces.
Thus, the total drag forces generated by the th filament along the and -axes are
(S18)
From the force balance, the propulsion force exerted to the fluid due to the beat motion of a filament
are and .
Thus, the total propulsion forces are obtained by summing over the contributions from filaments, and , and the tilt angle () dependent total propulsion forces along the - and -axes are
(S19)
The mean propulsion forces over the period of beat motion in the - and -directions are
Next, the velocity of the fluid flow at position at time generated by the undulatory motion of the th filament is
(S21)
where defined in Eq. (S17), is in Eq. (S5), , and .
The flow velocity along the and -axes are
(S22)
Therefore, the fluid flow generated by the collective motion of filaments is and .
Fig. 5 shows the mean values of at location from the -axis.
VII Viscous dissipation and work production
First, the mean viscous dissipation per filament (Fig. 5A) is evaluated as:
(S23)
Next, since the swimming velocity along the - and -directions [31, 32], resulting from the propulsion force of a filament are
and ,
the velocities along the - and -directions can be expressed as
and .
Consequently, the rate of work generated and transmitted to the fluid,
which enhances the fluid flow along the -direction, is: