No Arabic abstract
Cilia are elastic hairlike protuberances of the cell membrane found in various unicellular organisms and in several tissues of most living organisms. In some tissues such as the airway tissues of the lung, the coordinated beating of cilia induce a fluid flow of crucial importance as it allows the continuous cleaning of our bronchia, known as mucociliary clearance. While most of the models addressing the question of collective dynamics and metachronal wave consider homogeneous carpets of cilia, experimental observations rather show that cilia clusters are heterogeneously distributed over the tissue surface. The purpose of this paper is to investigate the role of spatial heterogeneity on the coherent beating of cilia using a very simple one dimensional model for cilia known as the rower model. We systematically study systems consisting of a few rowers to hundreds of rowers and we investigate the conditions for the emergence of collective beating. When considering a small number of rowers, a phase drift occurs, hence a bifurcation in beating frequency is observed as the distance between rowers clusters is changed. In the case of many rowers, a distribution of frequencies is observed. We found in particular the pattern of the patchy structure that shows the best robustness in collective beating behavior, as the density of cilia is varied over a wide range.
In this Letter, we study the collective behaviour of a large number of self-propelled microswimmers immersed in a fluid. Using unprecedently large-scale lattice Boltzmann simulations, we reproduce the transition to bacterial turbulence. We show that, even well below the transition, swimmers move in a correlated fashion that cannot be described by a mean-field approach. We develop a novel kinetic theory that captures these correlations and is non-perturbative in the swimmer density. To provide an experimentally accessible measure of correlations, we calculate the diffusivity of passive tracers and reveal its non-trivial density dependence. The theory is in quantitative agreement with the lattice Boltzmann simulations and captures the asymmetry between pusher and puller swimmers below the transition to turbulence.
Synchronization among arrays of beating cilia is one of the emergent phenomena in biological processes at meso-scopic scales. Strong inter-ciliary couplings modify the natural beating frequencies, $omega$, of individual cilia to produce a collective motion that moves around a group frequency $omega_m$. Here we study the thermodynamic cost of synchronizing cilia arrays by mapping their dynamics onto a generic phase oscillator model. The model suggests that upon synchronization the mean heat dissipation rate is decomposed into two contributions, dissipation from each ciliums own natural driving force and dissipation arising from the interaction with other cilia, the latter of which can be interpreted as the one produced by a potential with a time-dependent protocol in the framework of our model. The spontaneous phase-synchronization of beating dynamics of cilia induced by strong inter-ciliary coupling is always accompanied with a significant reduction of dissipation for the cilia population, suggesting that organisms as a whole expend less energy by attaining a temporal order. At the level of individual cilia, however, a population of cilia with $|omega|< omega_m$ expend more amount of energy upon synchronization.
We calculate the hydrodynamic flow field generated far from a cilium which is attached to a surface and beats periodically. In the case of two beating cilia, hydrodynamic interactions can lead to synchronization of the cilia, which are nonlinear oscillators. We present a state diagram where synchronized states occur as a function of distance of cilia and the relative orientation of their beat. Synchronized states occur with different relative phases. In addition, asynchronous solutions exist. Our work could be relevant for the synchronized motion of cilia generating hydrodynamic flows on the surface of cells.
Ventricular tachycardia (VT) and ventricular fibrillation (VF) are lethal rhythm disorders, which are associated with the occurrence of abnormal electrical scroll waves in the heart. Given the technical limitations of imaging and probing, the in situ visualization of these waves inside cardiac tissue remains a challenge. Therefore, we must, perforce, rely on in-silico simulations of scroll waves in mathematical models for cardiac tissue to develop an understanding of the dynamics of these waves in mammalian hearts. We use direct numerical simulations of the Hund-Rudy-Dynamic (HRD) model, for canine ventricular tissue, to examine the interplay between electrical scroll-waves and conduction and ionic inhomogeneities, in anatomically realistic canine ventricular geometries with muscle-fiber architecture. We find that millimeter-sized, distributed, conduction inhomogeneities cause a substantial decrease in the scroll wavelength, thereby increasing the probability for wave breaks; by contrast, single, localized, medium-sized ($simeq $ cm) conduction inhomogeneities, exhibit the potential to suppress wave breaks or enable the self-organization of wave fragments into stable, intact scrolls. We show that ionic inhomogeneities, both distributed or localised, suppress scroll-wave break up. The dynamics of a stable rotating wave is not affected significantly by such inhomogeneities, except at high concentrations of distributed inhomogeneities, which can cause a partial break up of scroll waves. Our results indicate that inhomogeneities in the canine ventricular tissue are less arrhythmogenic than inhomogeneities in porcine ventricular tissue, for which an earlier in silico study has shown that the inhomogeneity-induced suppression of scroll waves is a rare occurrence.
We investigate the heterogeneity of dynamics, the breakdown of the Stokes-Einstein relation and fragility in a model glass forming liquid, a binary mixture of soft spheres with a harmonic interaction potential, for spatial dimensions from 3 to 8. Dynamical heterogeneity is quantified through the dynamical susceptibility $chi_4$, and the non-Gaussian parameter $alpha_2$. We find that the fragility, the degree of breakdown of the Stokes-Einstein relation, as well as heterogeneity of dynamics, decrease with increasing spatial dimensionality. We briefly describe the dependence of fragility on density, and use it to resolve an apparent inconsistency with previous results.