Physical properties of collective motion in suspensions of bacteria

Effect of Cell Aspect Ratio on Swarming Bacteria

Swarming bacteria collectively migrate on surfaces using flagella, forming dynamic whirls and jets that consist of millions of individuals. Because some swarming bacteria elongate prior to actual motion, cell aspect ratio may play a significant role in the collective dynamics. Extensive research on self-propelled rodlike particles confirms that elongation promotes alignment, strongly affecting the dynamics. Here, we study experimentally the collective dynamics of variants of swarming Bacillus subtilis that differ in length. We show that the swarming statistics depends on the aspect ratio in a critical, fundamental fashion not predicted by theory. The fastest motion was obtained for the wild-type and variants that are similar in length. However, shorter and longer cells exhibit anomalous, non-Gaussian statistics and nonexponential decay of the autocorrelation function, indicating lower collective motility. These results suggest that the robust mechanisms to maintain aspect ratios may be important for efficient swarming motility. Wild-type cells are optimal in this sense.

Emergence of Lévy walks in systems of interacting individuals

We propose a model of superdiffusive Lévy walk as an emergent nonlinear phenomenon in systems of interacting individuals. The aim is to provide a qualitative explanation of recent experiments [G. Ariel et al., Nat. Commun. 6, 8396 (2015)2041-172310.1038/ncomms9396] revealing an intriguing behavior: swarming bacteria fundamentally change their collective motion from simple diffusion into a superdiffusive Lévy walk dynamics. We introduce microscopic mean-field kinetic equations in which we combine two key ingredients: (1) alignment interactions between individuals and (2) non-Markovian effects. Our interacting run-and-tumble model leads to the superdiffusive growth of the mean-squared displacement and the power-law distribution of run length with infinite variance. The main result is that the superdiffusive behavior emerges as a cooperative effect without using the standard assumption of the power-law distribution of run distances from the inception. At the same time, we find that the collision and repulsion interactions lead to the density-dependent exponential tempering of power-law distributions. This qualitatively explains the experimentally observed transition from superdiffusion to the diffusion of mussels as their density increases [M. de Jager et al., Proc. R. Soc. B 281, 20132605 (2014)PRLBA40962-845210.1098/rspb.2013.2605].

Spatial confinement of active microtubule networks induces large-scale rotational cytoplasmic flow

Collective behaviors of motile units through hydrodynamic interactions induce directed fluid flow on a larger length scale than individual units. In cells, active cytoskeletal systems composed of polar filaments and molecular motors drive fluid flow, a process known as cytoplasmic streaming. The motor-driven elongation of microtubule bundles generates turbulent-like flow in purified systems; however, it remains unclear whether and how microtubule bundles induce large-scale directed flow like the cytoplasmic streaming observed in cells. Here, we adopted Xenopus egg extracts as a model system of the cytoplasm and found that microtubule bundle elongation induces directed flow for which the length scale and timescale depend on the existence of geometrical constraints. At the lower activity of dynein, kinesins bundle and slide microtubules, organizing extensile microtubule bundles. In bulk extracts, the extensile bundles connected with each other and formed a random network, and vortex flows with a length scale comparable to the bundle length continually emerged and persisted for 1 min at multiple places. When the extracts were encapsulated in droplets, the extensile bundles pushed the droplet boundary. This pushing force initiated symmetry breaking of the randomly oriented bundle network, leading to bundles aligning into a rotating vortex structure. This vortex induced rotational cytoplasmic flows on the length scale and timescale that were 10- to 100-fold longer than the vortex flows emerging in bulk extracts. Our results suggest that microtubule systems use not only hydrodynamic interactions but also mechanical interactions to induce large-scale temporally stable cytoplasmic flow.

Spontaneous mirror-symmetry breaking induces inverse energy cascade in 3D active fluids

Classical turbulence theory assumes that energy transport in a 3D turbulent flow proceeds through a Richardson cascade whereby larger vortices successively decay into smaller ones. By contrast, an additional inverse cascade characterized by vortex growth exists in 2D fluids and gases, with profound implications for meteorological flows and fluid mixing. The possibility of a helicity-driven inverse cascade in 3D fluids had been rejected in the 1970s based on equilibrium-thermodynamic arguments. Recently, however, it was proposed that certain symmetry-breaking processes could potentially trigger a 3D inverse cascade, but no physical system exhibiting this phenomenon has been identified to date. Here, we present analytical and numerical evidence for the existence of an inverse energy cascade in an experimentally validated 3D active fluid model, describing microbial suspension flows that spontaneously break mirror symmetry. We show analytically that self-organized scale selection, a generic feature of many biological and engineered nonequilibrium fluids, can generate parity-violating Beltrami flows. Our simulations further demonstrate how active scale selection controls mirror-symmetry breaking and the emergence of a 3D inverse cascade.

A criterion of collective behavior of bacteria

It was established in the previous works that hydrodynamic interactions between the swimmers can lead to collective motion. Its implicit evidences were confirmed by reduction in the effective viscosity. We propose a new quantitative criterion to detect such a collective behavior. Our criterion is based on a new computationally effective RVE (representative volume element) theory based on the basic statistic moments (e-sums or generalized Eisenstein-Rayleigh sums). The criterion can be applied to various two-phase dispersed media (biological systems, composites etc). The locations of bacteria are modeled by short segments having a small width randomly embedded in medium without overlapping. We compute the e-sums of the simulated disordered sets and of the observed experimental locations of Bacillus subtilis. The obtained results show a difference between these two sets that demonstrates the collective motion of bacteria.

Long-range nematic order and anomalous fluctuations in suspensions of swimming filamentous bacteria

We study the collective dynamics of elongated swimmers in a very thin fluid layer by devising long filamentous nontumbling bacteria. The strong confinement induces weak nematic alignment upon collision, which, for large enough density of cells, gives rise to global nematic order. This homogeneous but fluctuating phase, observed on the largest experimentally accessible scale of millimeters, exhibits the properties predicted by standard models for flocking, such as the Vicsek-style model of polar particles with nematic alignment: true long-range nematic order and nontrivial giant number fluctuations.

Weak synchronization and large-scale collective oscillation in dense bacterial suspensions

Collective oscillatory behaviour is ubiquitous in nature, having a vital role in many biological processes from embryogenesis and organ development to pace-making in neuron networks. Elucidating the mechanisms that give rise to synchronization is essential to the understanding of biological self-organization. Collective oscillations in biological multicellular systems often arise from long-range coupling mediated by diffusive chemicals, by electrochemical mechanisms, or by biomechanical interaction between cells and their physical environment. In these examples, the phase of some oscillatory intracellular degree of freedom is synchronized. Here, in contrast, we report the discovery of a weak synchronization mechanism that does not require long-range coupling or inherent oscillation of individual cells. We find that millions of motile cells in dense bacterial suspensions can self-organize into highly robust collective oscillatory motion, while individual cells move in an erratic manner, without obvious periodic motion but with frequent, abrupt and random directional changes. So erratic are individual trajectories that uncovering the collective oscillations of our micrometre-sized cells requires individual velocities to be averaged over tens or hundreds of micrometres. On such large scales, the oscillations appear to be in phase and the mean position of cells typically describes a regular elliptic trajectory. We found that the phase of the oscillations is organized into a centimetre-scale travelling wave. We present a model of noisy self-propelled particles with strictly local interactions that accounts faithfully for our observations, suggesting that self-organized collective oscillatory motion results from spontaneous chiral and rotational symmetry breaking. These findings reveal a previously unseen type of long-range order in active matter systems (those in which energy is spent locally to produce non-random motion). This mechanism of collective oscillation may inspire new strategies to control the self-organization of active matter and swarming robots.

Dynamics of ellipsoidal tracers in swimming algal suspensions

Enhanced diffusion of passive tracers immersed in active fluids is a universal feature of active fluids and has been extensively studied in recent years. Similar to microrheology for equilibrium complex fluids, the unusual enhanced particle dynamics reveal intrinsic properties of active fluids. Nevertheless, previous studies have shown that the translational dynamics of spherical tracers are qualitatively similar, independent of whether active particles are pushers or pullers-the two fundamental classes of active fluids. Is it possible to distinguish pushers from pullers by simply imaging the dynamics of passive tracers? Here, we investigated the diffusion of isolated ellipsoids in algal C. reinhardtii suspensions-a model for puller-type active fluids. In combination with our previous results on pusher-type E. coli suspensions [Peng et al., Phys. Rev. Lett. 116, 068303 (2016)PRLTAO0031-900710.1103/PhysRevLett.116.068303], we showed that the dynamics of asymmetric tracers show a profound difference in pushers and pullers due to their rotational degree of freedom. Although the laboratory-frame translation and rotation of ellipsoids are enhanced in both pushers and pullers, similar to spherical tracers, the anisotropic diffusion in the body frame of ellipsoids shows opposite trends in the two classes of active fluids. An ellipsoid diffuses fastest along its major axis when immersed in pullers, whereas it diffuses slowest along the major axis in pushers. This striking difference can be qualitatively explained using a simple hydrodynamic model. In addition, our study on algal suspensions reveals that the influence of the near-field advection of algal swimming flows on the translation and rotation of ellipsoids shows different ranges and strengths. Our work provides not only new insights into universal organizing principles of active fluids, but also a convenient tool for detecting the class of active particles.

Generalized Swift-Hohenberg models for dense active suspensions

In describing the physics of living organisms, a mathematical theory that captures the generic ordering principles of intracellular and multicellular dynamics is essential for distinguishing between universal and system-specific features. Here, we compare two recently proposed nonlinear high-order continuum models for active polar and nematic suspensions, which aim to describe collective migration in dense cell assemblies and the ordering processes in ATP-driven microtubule-kinesin networks, respectively. We discuss the phase diagrams of the two models and relate their predictions to recent experiments. The satisfactory agreement with existing experimental data lends support to the hypothesis that non-equilibrium pattern formation phenomena in a wide range of active systems can be described within the same class of higher-order partial differential equations.

Hydrodynamic length-scale selection in microswimmer suspensions

A universal characteristic of mesoscale turbulence in active suspensions is the emergence of a typical vortex length scale, distinctly different from the scale invariance of turbulent high-Reynolds number flows. Collective length-scale selection has been observed in bacterial fluids, endothelial tissue, and active colloids, yet the physical origins of this phenomenon remain elusive. Here, we systematically derive an effective fourth-order field theory from a generic microscopic model that allows us to predict the typical vortex size in microswimmer suspensions. Building on a self-consistent closure condition, the derivation shows that the vortex length scale is determined by the competition between local alignment forces, rotational diffusion, and intermediate-range hydrodynamic interactions. Vortex structures found in simulations of the theory agree with recent measurements in Bacillus subtilis suspensions. Moreover, our approach yields an effective viscosity enhancement (reduction), as reported experimentally for puller (pusher) microorganisms.

Mode instabilities and dynamic patterns in a colony of self-propelled surfactant particles covering a thin liquid layer

We consider a colony of point-like self-propelled surfactant particles (swimmers) without direct interactions that cover a thin liquid layer on a solid support. The particles predominantly swim normal to the free film surface with only a small component parallel to the film surface. The coupled dynamics of the swimmer density and film height profile is captured in a long-wave model allowing for diffusive and convective transport of the swimmers (including rotational diffusion). The dynamics of the film height profile is determined by i) the upward pushing force of the swimmers onto the liquid-gas interface, ii) the solutal Marangoni force due to gradients in the swimmer concentration, and iii) the rotational diffusion of the swimmers together with the in-plane active motion. After reviewing and extending the analysis of the linear stability of the uniform state, we analyse the fully nonlinear dynamic equations and show that point-like swimmers, which only interact via long-wave deformations of the liquid film, self-organise in highly regular (standing, travelling, and modulated waves) and various irregular patterns.

Fabricating large two-dimensional single colloidal crystals by doping with active particles

Using simulations we explore the behaviour of two-dimensional colloidal (poly)crystals doped with active particles. We show that these active dopants can provide an elegant new route to removing grain boundaries in polycrystals. Specifically, we show that active dopants both generate and are attracted to defects, such as vacancies and interstitials, which leads to clustering of dopants at grain boundaries. The active particles both broaden and enhance the mobility of the grain boundaries, causing rapid coarsening of the crystal domains. The remaining defects recrystallize upon turning off the activity of the dopants, resulting in a large-scale single-domain crystal.

Ferromagnetic and antiferromagnetic order in bacterial vortex lattices

Despite their inherent non-equilibrium nature1, living systems can self-organize in highly ordered collective states2,3 that share striking similarities with the thermodynamic equilibrium phases4,5 of conventional condensed matter and fluid systems. Examples range from the liquid-crystal-like arrangements of bacterial colonies6,7, microbial suspensions8,9 and tissues10 to the coherent macro-scale dynamics in schools of fish11 and flocks of birds12. Yet, the generic mathematical principles that govern the emergence of structure in such artificial13 and biological6-9,14 systems are elusive. It is not clear when, or even whether, well-established theoretical concepts describing universal thermostatistics of equilibrium systems can capture and classify ordered states of living matter. Here, we connect these two previously disparate regimes: Through microfluidic experiments and mathematical modelling, we demonstrate that lattices of hydrodynamically coupled bacterial vortices can spontaneously organize into distinct phases of ferro- and antiferromagnetic order. The preferred phase can be controlled by tuning the vortex coupling through changes of the inter-cavity gap widths. The emergence of opposing order regimes is tightly linked to the existence of geometry-induced edge currents15,16, reminiscent of those in quantum systems17-19. Our experimental observations can be rationalized in terms of a generic lattice field theory, suggesting that bacterial spin networks belong to the same universality class as a wide range of equilibrium systems.

Effective Rheological Properties in Semi-dilute Bacterial Suspensions

Interactions between swimming bacteria have led to remarkable experimentally observable macroscopic properties such as the reduction in the effective viscosity, enhanced mixing, and diffusion. In this work, we study an individual-based model for a suspension of interacting point dipoles representing bacteria in order to gain greater insight into the physical mechanisms responsible for the drastic reduction in the effective viscosity. In particular, asymptotic analysis is carried out on the corresponding kinetic equation governing the distribution of bacteria orientations. This allows one to derive an explicit asymptotic formula for the effective viscosity of the bacterial suspension in the limit of bacterium non-sphericity. The results show good qualitative agreement with numerical simulations and previous experimental observations. Finally, we justify our approach by proving existence, uniqueness, and regularity properties for this kinetic PDE model.

Rapid expulsion of microswimmers by a vortical flow

Interactions of microswimmers with their fluid environment are exceptionally complex. Macroscopic shear flow alters swimming trajectories in a highly nontrivial way and results in dramatic reduction of viscosity and heterogeneous bacterial distributions. Here we report on experimental and theoretical studies of rapid expulsion of microswimmers, such as motile bacteria, by a vortical flow created by a rotating microparticle. We observe a formation of a macroscopic depletion area in a high-shear region, in the vicinity of a microparticle. The rapid migration of bacteria from the shear-rich area is caused by a vortical structure of the flow rather than intrinsic random fluctuations of bacteria orientations, in stark contrast to planar shear flow. Our mathematical model reveals that expulsion is a combined effect of motility and alignment by a vortical flow. Our findings offer a novel approach for manipulation of motile microorganisms and shed light on bacteria–flow interactions.

The control of microswimmers such as bacteria is important for emerging applications of active bioinspired materials. Here, the authors demonstrate the use of vortical shear to expel suspended motile bacteria from the vicinity of a rotating microparticle.

Diffusion of Ellipsoids in Bacterial Suspensions

Active fluids such as swarming bacteria and motile colloids exhibit exotic properties different from conventional equilibrium materials. As a peculiar example, a spherical tracer immersed inside active fluids shows an enhanced translational diffusion, orders of magnitude stronger than its intrinsic Brownian motion. Here, rather than spherical tracers, we investigate the diffusion of isolated ellipsoids in a quasi-two-dimensional bacterial bath. Our study shows a nonlinear enhancement of both translational and rotational diffusions of ellipsoids. More importantly, we uncover an anomalous coupling between particles' translation and rotation that is strictly prohibited in Brownian diffusion. The coupling reveals a counterintuitive anisotropic particle diffusion, where an ellipsoid diffuses fastest along its minor axis in its body frame. Combining experiments with theoretical modeling, we show that such an anomalous diffusive behavior arises from the generic straining flow of swimming bacteria. Our work illustrates an unexpected feature of active fluids and deepens our understanding of transport processes in microbiological systems.

Shape matters: Near-field fluid mechanics dominate the collective motions of ellipsoidal squirmers

Microswimmers show a variety of collective motions. Despite extensive study, questions remain regarding the role of near-field fluid mechanics in collective motion. In this paper, we describe precisely the Stokes flow around hydrodynamically interacting ellipsoidal squirmers in a monolayer suspension. The results showed that various collective motions, such as ordering, aggregation, and whirls, are dominated by the swimming mode and the aspect ratio. The collective motions are mainly induced by near-field fluid mechanics, despite Stokes flow propagation over a long range. These results emphasize the importance of particle shape in collective motion.

New class of turbulence in active fluids

Turbulence is a fundamental and ubiquitous phenomenon in nature, occurring from astrophysical to biophysical scales. At the same time, it is widely recognized as one of the key unsolved problems in modern physics, representing a paradigmatic example of nonlinear dynamics far from thermodynamic equilibrium. Whereas in the past, most theoretical work in this area has been devoted to Navier-Stokes flows, there is now a growing awareness of the need to extend the research focus to systems with more general patterns of energy injection and dissipation. These include various types of complex fluids and plasmas, as well as active systems consisting of self-propelled particles, like dense bacterial suspensions. Recently, a continuum model has been proposed for such "living fluids" that is based on the Navier-Stokes equations, but extends them to include some of the most general terms admitted by the symmetry of the problem [Wensink HH, et al. (2012) Proc Natl Acad Sci USA 109:14308-14313]. This introduces a cubic nonlinearity, related to the Toner-Tu theory of flocking, which can interact with the quadratic Navier-Stokes nonlinearity. We show that as a result of the subtle interaction between these two terms, the energy spectra at large spatial scales exhibit power laws that are not universal, but depend on both finite-size effects and physical parameters. Our combined numerical and analytical analysis reveals the origin of this effect and even provides a way to understand it quantitatively. Turbulence in active fluids, characterized by this kind of nonlinear self-organization, defines a new class of turbulent flows.

Comparison between Smoluchowski and Boltzmann approaches for self-propelled rods

Considering systems of self-propelled polar particles with nematic interactions ("rods"), we compare the continuum equations describing the evolution of polar and nematic order parameters, derived either from Smoluchowski or Boltzmann equations. Our main goal is to understand the discrepancies between the continuum equations obtained so far in both frameworks. We first show that, in the simple case of point-like particles with only alignment interactions, the continuum equations obtained have the same structure in both cases. We further study, in the Smoluchowski framework, the case where an interaction force is added on top of the aligning torque. This clarifies the origin of the additional terms obtained in previous works. Our observations lead us to emphasize the need for a more involved closure scheme than the standard normal form of the distribution when dealing with active systems.

Swarming bacteria migrate by Lévy Walk

Individual swimming bacteria are known to bias their random trajectories in search of food and to optimize survival. The motion of bacteria within a swarm, wherein they migrate as a collective group over a solid surface, is fundamentally different as typical bacterial swarms show large-scale swirling and streaming motions involving millions to billions of cells. Here by tracking trajectories of fluorescently labelled individuals within such dense swarms, we find that the bacteria are performing super-diffusion, consistent with Lévy walks. Lévy walks are characterized by trajectories that have straight stretches for extended lengths whose variance is infinite. The evidence of super-diffusion consistent with Lévy walks in bacteria suggests that this strategy may have evolved considerably earlier than previously thought.

Turbulence of swarming sperm

Collective motion of self-sustained swarming flows has recently provided examples of small-scale turbulence arising where viscous effects are dominant. We report the first observation of universal enstrophy cascade in concentrated swarming sperm consistent with a body of evidence built from various independent measurements. We found a well-defined k^{-3} power-law decay of a velocity field power spectrum and relative dispersion of small beads consistent with theoretical predictions in 2D turbulence. Concentrated living sperm displays long-range, correlated whirlpool structures of a size that provides an integral scale of turbulence. We propose a consistent explanation for this quasi-2D turbulence based on self-structured laminated flow forced by steric interactions and alignment, a state of active matter that we call "swarming liquid crystal." We develop scaling arguments consistent with this interpretation.

Flexibility of bacterial flagella in external shear results in complex swimming trajectories

Many bacteria use rotating helical flagella in swimming motility. In the search for food or migration towards a new habitat, bacteria occasionally unbundle their flagellar filaments and tumble, leading to an abrupt change in direction. Flexible flagella can also be easily deformed by external shear flow, leading to complex bacterial trajectories. Here, we examine the effects of flagella flexibility on the navigation of bacteria in two fundamental shear flows: planar shear and Poiseuille flow realized in long channels. On the basis of slender body elastodynamics and numerical analysis, we discovered a variety of non-trivial effects stemming from the interplay of self-propulsion, elasticity and shear-induced flagellar bending. We show that in planar shear flow the bacteria execute periodic motion, whereas in Poiseuille flow, they migrate towards the centre of the channel or converge towards a limit cycle. We also find that even a small amount of random reorientation can induce a strong response of bacteria, leading to overall non-periodic trajectories. Our findings exemplify the sensitive role of flagellar flexibility and shed new light on the navigation of bacteria in complex shear flows.

Turning Bacteria Suspensions into Superfluids

The rheological response under simple shear of an active suspension of Escherichia coli is determined in a large range of shear rates and concentrations. The effective viscosity and the time scales characterizing the bacterial organization under shear are obtained. In the dilute regime, we bring evidence for a low-shear Newtonian plateau characterized by a shear viscosity decreasing with concentration. In the semidilute regime, for particularly active bacteria, the suspension displays a "superfluidlike" transition where the viscous resistance to shear vanishes, thus showing that, macroscopically, the activity of pusher swimmers organized by shear is able to fully overcome the dissipative effects due to viscous loss.

Suppressing turbulence of self-propelling rods by strongly coupled passive particles

The strong turbulence suppression, mainly for large-scale modes, of two-dimensional self-propelling rods, by increasing the long-range coupling strength Γ of low-concentration passive particles, is numerically demonstrated. It is found that large-scale collective rod motion in forms of swirls or jets is mainly contributed from well-aligned dense patches, which can push small poorly aligned rod patches and uncoupled passive particles. The more efficient momentum transfer and dissipation through increasing passive particle coupling leads to the formation of a more ordered and slowed down network of passive particles, which competes with coherent dense active rod clusters. The frustration of active rod alignment ordering and coherent motion by the passive particle network, which interrupt the inverse cascading of forming large-scale swirls, is the key for suppressing collective rod motion with scales beyond the interpassive distance, even in the liquid phase of passive particles. The loosely packed active rods are weakly affected by increasing passive particle coupling due to the weak rod-particle interaction. They mainly contribute to the small-scale modes and high-speed motion.

Antibiotic-induced anomalous statistics of collective bacterial swarming

Under sublethal antibiotics concentrations, the statistics of collectively swarming Bacillus subtilis transitions from normal to anomalous, with a heavy-tailed speed distribution and a two-step temporal correlation of velocities. The transition is due to changes in the properties of the bacterial motion and the formation of a motility-defective subpopulation that self-segregates into regions. As a result, both the colonial expansion and the growth rate are not affected by antibiotics. This phenomenon suggests a new strategy bacteria employ to fight antibiotic stress.

Individual behavior and pairwise interactions between microswimmers in anisotropic liquid

A motile bacterium swims by generating flow in its surrounding liquid. Anisotropy of the suspending liquid significantly modifies the swimming dynamics and corresponding flow signatures of an individual bacterium and impacts collective behavior. We study the interactions between swimming bacteria in an anisotropic environment exemplified by lyotropic chromonic liquid crystal. Our analysis reveals a significant localization of the bacteria-induced flow along a line coaxial with the bacterial body, which is due to strong viscosity anisotropy of the liquid crystal. Despite the fact that the average viscosity of the liquid crystal is two to three orders of magnitude higher than the viscosity of pure water, the speed of bacteria in the liquid crystal is of the same order of magnitude as in water. We show that bacteria can transport a cargo (a fluorescent particle) along a predetermined trajectory defined by the direction of molecular orientation of the liquid crystal. We demonstrate that while the hydrodynamic interaction between flagella of two close-by bacteria is negligible, the observed convergence of the swimming speeds as well as flagella waves' phase velocities may occur due to viscoelastic interaction between the bacterial bodies.

Vortex arrays and mesoscale turbulence of self-propelled particles

Inspired by the Turing mechanism for pattern formation, we propose a simple self-propelled particle model with short-range alignment and antialignment at larger distances. It is able to produce orientationally ordered states, periodic vortex patterns, and mesoscale turbulence, which resembles observations in dense suspensions of swimming bacteria. The model allows a systematic derivation and analysis of a kinetic theory as well as hydrodynamic equations for density and momentum fields. A phase diagram with regions of pattern formation as well as orientational order is obtained from a linear stability analysis of these continuum equations. Microscopic Langevin simulations of self-propelled particles are in agreement with these findings.

Active nematic materials with substrate friction

Active turbulence in dense active systems is characterized by high vorticity on a length scale that is large compared to that of individual entities. We describe the properties of active turbulence as momentum propagation is screened by frictional damping. As friction is increased, the spacing between the walls in the nematic director field decreases as a consequence of the more rapid velocity decays. This leads to, first, a regime with more walls and an increased number of topological defects, and then to a jammed state in which the walls deliminate bands of opposing flow, analogous to the shear bands observed in passive complex fluids.

Vorticity, defects and correlations in active turbulence

We describe a numerical investigation of a continuum model of an active nematic, concentrating on the regime of active turbulence. Results are presented for the effect of three parameters, activity, elastic constant and rotational diffusion constant, on the order parameter and flow fields. Defects and distortions in the director field act as sources of vorticity, and thus vorticity is strongly correlated to the director field. In particular, the characteristic length of decay of vorticity and order parameter correlations is controlled by the defect density. By contrast, the decay of velocity correlations is determined by a balance between activity and dissipation. We highlight the role of microscopic flow generation mechanisms in determining the flow patterns and characteristic scales of active turbulence and contrast the behaviour of extensile and contractile active nematics.

Subdiffusion-reaction processes with A→B reactions versus subdiffusion-reaction processes with A+B→B reactions

We consider the subdiffusion-reaction process with reactions of a type A+B→B (in which particles A are assumed to be mobile, whereas B are assumed to be static) in comparison to the subdiffusion-reaction process with A→B reactions which was studied by Sokolov, Schmidt, and Sagués [Phys. Rev. E 73, 031102 (2006)]. In both processes a rule that reactions can only occur between particles which continue to exist is taken into account. Although in both processes a probability of the vanishing of particle A due to a reaction is independent of both time and space variables (assuming that in the system with the A+B→B reactions, particles B are distributed homogeneously), we show that subdiffusion-reaction equations describing these processes as well as their Green's functions are qualitatively different. The reason for this difference is as follows. In the case of the former reaction, particles A and B have to meet with some probability before the reaction occurs in contradiction with the case of the latter reaction. For the subdiffusion process with the A+B→B reactions we consider three models which differ in some details concerning a description of the reactions. We base the method considered in this paper on a random walk model in a system with both discrete time and discrete space variables. Then the system with discrete variables is transformed into a system with both continuous time and continuous space variables. Such a method seems to be convenient in analyzing subdiffusion-reaction processes with partially absorbing or partially reflecting walls. The reason is that within this method we can determine Green's functions without a necessity of solving a fractional differential subdiffusion-reaction equation with boundary conditions at the walls. As an example, we use the model to find the Green's functions for a subdiffusive reaction system (with the reactions mentioned above), which is bounded by a partially absorbing wall. This example shows how the model can be used to analyze the subdiffusion-reaction process in a system with partially absorbing or reflecting thin membranes. Employing a simple phenomenological model, we also derive equations related to the reaction parameters used in the considered models.

Large-scale chaos and fluctuations in active nematics

We show that dry active nematics, e.g., collections of shaken elongated granular particles, exhibit large-scale spatiotemporal chaos made of interacting dense, ordered, bandlike structures in a parameter region including the linear onset of nematic order. These results are obtained from the study of both the well-known (deterministic) hydrodynamic equations describing these systems and of the self-propelled particle model they were derived from. We prove, in particular, that the chaos stems from the generic instability of the band solution of the hydrodynamic equations. Revisiting the status of the strong fluctuations and long-range correlations in the particle model, we show that the giant number fluctuations observed in the chaotic phase are a trivial consequence of density segregation. However anomalous, curvature-driven number fluctuations are present in the homogeneous quasiordered nematic phase and characterized by a nontrivial scaling exponent.

Fluid flows created by swimming bacteria drive self-organization in confined suspensions

Concentrated suspensions of swimming microorganisms and other forms of active matter are known to display complex, self-organized spatiotemporal patterns on scales that are large compared with those of the individual motile units. Despite intensive experimental and theoretical study, it has remained unclear the extent to which the hydrodynamic flows generated by swimming cells, rather than purely steric interactions between them, drive the self-organization. Here we use the recent discovery of a spiral-vortex state in confined suspensions of Bacillus subtilis to study this issue in detail. Those experiments showed that if the radius of confinement in a thin cylindrical chamber is below a critical value, the suspension will spontaneously form a steady single-vortex state encircled by a counter-rotating cell boundary layer, with spiral cell orientation within the vortex. Left unclear, however, was the flagellar orientation, and hence the cell swimming direction, within the spiral vortex. Here, using a fast simulation method that captures oriented cell-cell and cell-fluid interactions in a minimal model of discrete particle systems, we predict the striking, counterintuitive result that in the presence of collectively generated fluid motion, the cells within the spiral vortex actually swim upstream against those flows. This prediction is then confirmed by the experiments reported here, which include measurements of flagella bundle orientation and cell tracking in the self-organized state. These results highlight the complex interplay between cell orientation and hydrodynamic flows in concentrated suspensions of microorganisms.

Nonequilibrium glassy dynamics of self-propelled hard disks

We analyze the collective dynamics of self-propelled particles in the large-density regime where passive particles undergo a kinetic arrest to an amorphous glassy state. We capture the competition between self-propulsion and crowding effects using a two-dimensional model of self-propelled hard disks, which we study using Monte Carlo simulations. Although the activity drives the system far from equilibrium, self-propelled particles undergo a kinetic arrest, which we characterize in detail and compare with its equilibrium counterpart. In particular, the critical density for dynamic arrest continuously shifts to larger densities with increasing activity, and the relaxation time is surprisingly well described by an algebraic divergence resulting from the emergence of highly collective dynamics. These results show that dense assemblies of active particles undergo a nonequilibrium glass transition that is profoundly affected by self-propulsion mechanisms.

Gait synchronization in Caenorhabditis elegans

Collective motion is observed in swarms of swimmers of various sizes, ranging from self-propelled nanoparticles to fish. The mechanisms that govern interactions among individuals are debated, and vary from one species to another. Although the interactions among relatively large animals, such as fish, are controlled by their nervous systems, the interactions among microorganisms, which lack nervous systems, are controlled through physical and chemical pathways. Little is known, however, regarding the mechanism of collective movements in microscopic organisms with nervous systems. To attempt to remedy this, we studied collective swimming behavior in the nematode Caenorhabditis elegans, a microorganism with a compact nervous system. We evaluated the contributions of hydrodynamic forces, contact forces, and mechanosensory input to the interactions among individuals. We devised an experiment to examine pair interactions as a function of the distance between the animals and observed that gait synchronization occurred only when the animals were in close proximity, independent of genes required for mechanosensation. Our measurements and simulations indicate that steric hindrance is the dominant factor responsible for motion synchronization in C. elegans, and that hydrodynamic interactions and genotype do not play a significant role. We infer that a similar mechanism may apply to other microscopic swimming organisms and self-propelled particles.

Transport powered by bacterial turbulence

We demonstrate that collective turbulentlike motion in a bacterial bath can power and steer the directed transport of mesoscopic carriers through the suspension. In our experiments and simulations, a microwedgelike "bulldozer" draws energy from a bacterial bath of varied density. We obtain that an optimal transport speed is achieved in the turbulent state of the bacterial suspension. This apparent rectification of random motion of bacteria is caused by polar ordered bacteria inside the cusp region of the carrier, which is shielded from the outside turbulent fluctuations.

Living liquid crystals

Collective motion of self-propelled organisms or synthetic particles, often termed "active fluid," has attracted enormous attention in the broad scientific community because of its fundamentally nonequilibrium nature. Energy input and interactions among the moving units and the medium lead to complex dynamics. Here, we introduce a class of active matter--living liquid crystals (LLCs)--that combines living swimming bacteria with a lyotropic liquid crystal. The physical properties of LLCs can be controlled by the amount of oxygen available to bacteria, by concentration of ingredients, or by temperature. Our studies reveal a wealth of intriguing dynamic phenomena, caused by the coupling between the activity-triggered flow and long-range orientational order of the medium. Among these are (i) nonlinear trajectories of bacterial motion guided by nonuniform director, (ii) local melting of the liquid crystal caused by the bacteria-produced shear flows, (iii) activity-triggered transition from a nonflowing uniform state into a flowing one-dimensional periodic pattern and its evolution into a turbulent array of topological defects, and (iv) birefringence-enabled visualization of microflow generated by the nanometers-thick bacterial flagella. Unlike their isotropic counterpart, the LLCs show collective dynamic effects at very low volume fraction of bacteria, on the order of 0.2%. Our work suggests an unorthodox design concept to control and manipulate the dynamic behavior of soft active matter and opens the door for potential biosensing and biomedical applications.

Controlling active self-assembly through broken particle-shape symmetry

Many structural properties of conventional passive materials are known to arise from the symmetries of their microscopic constituents. By contrast, it is largely unclear how the interplay between particle shape and self-propulsion controls the meso- and macroscale behavior of active matter. Here we use large-scale simulations of homo- and heterogeneous self-propelled particle systems to identify generic effects of broken particle-shape symmetry on collective motion. We find that even small violations of fore-aft symmetry lead to fundamentally different collective behaviors, which may facilitate demixing of differently shaped species as well as the spontaneous formation of stable microrotors. These results suggest that variation of particle shape yields robust physical mechanisms to control self-assembly of active matter, with possibly profound implications for biology and materials design.

Collective motion of spherical bacteria

A large variety of motile bacterial species exhibit collective motions while inhabiting liquids or colonizing surfaces. These collective motions are often characterized by coherent dynamic clusters, where hundreds of cells move in correlated whirls and jets. Previously, all species that were known to form such motion had a rod-shaped structure, which enhances the order through steric and hydrodynamic interactions. Here we show that the spherical motile bacteria Serratia marcescens exhibit robust collective dynamics and correlated coherent motion while grown in suspensions. As cells migrate to the upper surface of a drop, they form a monolayer, and move collectively in whirls and jets. At all concentrations, the distribution of the bacterial speed was approximately Rayleigh with an average that depends on concentration in a non-monotonic way. Other dynamical parameters such as vorticity and correlation functions are also analyzed and compared to rod-shaped bacteria from the same strain. Our results demonstrate that self-propelled spherical objects do form complex ordered collective motion. This opens a door for a new perspective on the role of cell aspect ratio and alignment of cells with regards to collective motion in nature.

Active matter ratchets with an external drift

When active matter particles such as swimming bacteria are placed in an asymmetric array of funnels, it has been shown that a ratchet effect can occur even in the absence of an external drive. Here we examine active ratchets for two-dimensional arrays of funnels or L shapes where there is also an externally applied dc drive or drift. We show that for certain conditions the ratchet effect can be strongly enhanced and it is possible to have conditions under which run-and-tumble particles with one run length move in the opposite direction from particles with a different run length. For the arrays of L shapes, we find that the application of a drift force can enhance a transverse rectification in the direction perpendicular to the drift. When particle-particle steric interactions are included, we find that the ratchet effects can be either enhanced or suppressed depending on barrier geometry, particle run length, and particle density.

Correlation properties of collective motion in bacterial suspensions

The study of collective motion in bacterial suspensions has been of significant recent interest. To better understand the non-trivial spatio-temporal correlations emerging in the course of collective swimming in suspensions of motile bacteria, a simple model is employed: a bacterium is represented as a force dipole with size, through the use of a short-range repelling potential, and shape. The model emphasizes two fundamental mechanisms: dipolar hydrodynamic interactions and short-range bacterial collisions. Using direct particle simulations validated by a dedicated experiment, we show that changing the swimming speed or concentration alters the time scale of sustained collective motion, consistent with experiment. Also, the correlation length in the collective state is almost constant as concentration and swimming speed change even though increasing each greatly increases the input of energy to the system. We demonstrate that the particle shape is critical for the onset of collective effects. In addition, new experimental results are presented illustrating the onset of collective motion with an ultrasound technique. This work exemplifies the delicate balance between various physical mechanisms governing collective motion in bacterial suspensions and provides important insights into its mesoscopic nature.

Bacterial turbulence reduction by passive magnetic particle chains

We report the experimental observation of the bacterial turbulence reduction in dense E. coli suspensions by increasing the coupling of passive particle additives (paramagnetic particles). Applying an external magnetic field induces magnetic dipoles for particles and causes the formation of vertical chain bundles, which are hard for bacterial flows to tilt and break. The larger effective drag coefficient of chains causes slow horizontal motion of chains, which in turn form obstacles to suppress bacterial flows through the strong correlation in coherent bacterial clusters and intercluster interaction. The interruption of the upward energy flow from individual self-propelling bacteria to the larger scale in the bacterial turbulence with multiscaled coherent flow by the chain bundle leads to more severe suppression in the low frequency (wave number) regimes of the power spectra.

Confinement stabilizes a bacterial suspension into a spiral vortex

Confining surfaces play crucial roles in dynamics, transport, and order in many physical systems, but their effects on active matter, a broad class of dynamically self-organizing systems, are poorly understood. We investigate here the influence of global confinement and surface curvature on collective motion by studying the flow and orientational order within small droplets of a dense bacterial suspension. The competition between radial confinement, self-propulsion, steric interactions, and hydrodynamics robustly induces an intriguing steady single-vortex state, in which cells align in inward spiraling patterns accompanied by a thin counterrotating boundary layer. A minimal continuum model is shown to be in good agreement with these observations.