Fluid dynamics of bacterial turbulence

Hydrodynamic Clustering of Human Sperm in Viscoelastic Fluids

We have numerically investigated sperm clustering behaviours, modelling cells as superpositions of regularised flow singularities, coarse-grained from experimentally obtained digital microscopy of human sperm, both in watery medium and a highly viscous-weakly elastic, methylcellulose medium. We find that the cell yaw and cell pulling dynamics inhibit clustering in low viscosity media. In contrast clustering is readily visible in simulations modelling sperm within a methylcellulose medium, in line with previous observations that bovine sperm clustering is much more prominent in a rheological polyacrylamide medium. Furthermore, the fine-scale details of sperm flagellar movement substantially impact large-scale collective behaviours, further motivating the need for the digital microscopy and characterization of sperm to understand their dynamics.

Active nematics

Active matter extracts energy from its surroundings at the single particle level and transforms it into mechanical work. Examples include cytoskeleton biopolymers and bacterial suspensions. Here, we review experimental, theoretical and numerical studies of active nematics - a type of active system that is characterised by self-driven units with elongated shape. We focus primarily on microtubule-kinesin mixtures and the hydrodynamic theories that describe their properties. An important theme is active turbulence and the associated motile topological defects. We discuss ways in which active turbulence may be controlled, a pre-requisite to harvesting energy from active materials, and we consider the appearance, and possible implications, of active nematics and topological defects to cellular systems and biological processes.

Collective behavior of self-propelled rods with quorum sensing

Active agents-like phoretic particles, bacteria, sperm, and cytoskeletal filaments in motility assays-show a large variety of motility-induced collective behaviors, such as aggregation, clustering, and phase separation. The behavior of dense suspensions of engineered phoretic particles and of bacteria during biofilm formation is determined by two qualitatively different physical mechanisms: (i) volume exclusion (short-range steric repulsion) and (ii) quorum sensing (longer-range reduced propulsion due to alteration of the local chemical environment). To systematically characterize such systems, we study semi-penetrable self-propelled rods in two dimensions, with a propulsion force that decreases with increasing local rod density, by employing Brownian dynamics simulations. Volume exclusion and quorum sensing both lead to phase separation; however, the structure of the systems and the rod dynamics vastly differ. Quorum sensing enhances the polarity of the clusters, induces perpendicularity of rods at the cluster borders, and enhances cluster formation. For systems where the rods essentially become passive at high densities, formation of asters and stripes is observed. Systems of rods with larger aspect ratios show more ordered structures compared to those with smaller aspect ratios, due to their stronger alignment, with almost circular asters for strongly density-dependent propulsion force. With increasing range of the quorum-sensing interaction, the local density decreases, asters become less stable, and polar hedgehog clusters and clusters with domains appear.

Colloidal transport through trap arrays controlled by active microswimmers

We investigate the dynamics of a binary mixture consisting of active and passive colloidal particles diffusing in a 2D array of truncated harmonic wells, or traps. We explore the possibility of using a small fraction of active particles to manipulate a much larger fraction of passive particles, for instance, to confine them in or extract them from the traps. The results of our study have potential application in biology and medical sciences, for example, to remove dead cells or undesired contaminants from biological systems by means of self-propelled nano-robots.

State diagram of a three-sphere microswimmer in a channel

Geometric confinements are frequently encountered in soft matter systems and in particular significantly alter the dynamics of swimming microorganisms in viscous media. Surface-related effects on the motility of microswimmers can lead to important consequences in a large number of biological systems, such as biofilm formation, bacterial adhesion and microbial activity. On the basis of low-Reynolds-number hydrodynamics, we explore the state diagram of a three-sphere microswimmer under channel confinement in a slit geometry and fully characterize the swimming behavior and trajectories for neutral swimmers, puller- and pusher-type swimmers. While pushers always end up trapped at the channel walls, neutral swimmers and pullers may further perform a gliding motion and maintain a stable navigation along the channel. We find that the resulting dynamical system exhibits a supercritical pitchfork bifurcation in which swimming in the mid-plane becomes unstable beyond a transition channel height while two new stable limit cycles or fixed points that are symmetrically disposed with respect to the channel mid-height emerge. Additionally, we show that an accurate description of the averaged swimming velocity and rotation rate in a channel can be captured analytically using the method of hydrodynamic images, provided that the swimmer size is much smaller than the channel height.

Symmetric shear banding and swarming vortices in bacterial superfluids

Bacterial suspensions-a premier example of active fluids-show an unusual response to shear stresses. Instead of increasing the viscosity of the suspending fluid, the emergent collective motions of swimming bacteria can turn a suspension into a superfluid with zero apparent viscosity. Although the existence of active superfluids has been demonstrated in bulk rheological measurements, the microscopic origin and dynamics of such an exotic phase have not been experimentally probed. Here, using high-speed confocal rheometry, we study the dynamics of concentrated bacterial suspensions under simple planar shear. We find that bacterial superfluids under shear exhibit unusual symmetric shear bands, defying the conventional wisdom on shear banding of complex fluids, where the formation of steady shear bands necessarily breaks the symmetry of unsheared samples. We propose a simple hydrodynamic model based on the local stress balance and the ergodic sampling of nonequilibrium shear configurations, which quantitatively describes the observed symmetric shear-banding structure. The model also successfully predicts various interesting features of swarming vortices in stationary bacterial suspensions. Our study provides insights into the physical properties of collective swarming in active fluids and illustrates their profound influences on transport processes.

Turbulent Dynamics of Epithelial Cell Cultures

We investigate the large length and long time scales collective flows and structural rearrangements within in vitro human bronchial epithelial cell (HBEC) cultures. Activity-driven collective flows result in ensembles of vortices randomly positioned in space. By analyzing a large population of vortices, we show that their area follows an exponential law with a constant mean value and their rotational frequency is size independent, both being characteristic features of the chaotic dynamics of active nematic suspensions. Indeed, we find that HBECs self-organize in nematic domains of several cell lengths. Nematic defects are found at the interface between domains with a total number that remains constant due to the dynamical balance of nucleation and annihilation events. The mean velocity fields in the vicinity of defects are well described by a hydrodynamic theory of extensile active nematics.

Anomalous Chained Turbulence in Actively Driven Flows on Spheres

Recent experiments demonstrate the importance of substrate curvature for actively forced fluid dynamics. Yet, the covariant formulation and analysis of continuum models for nonequilibrium flows on curved surfaces still poses theoretical challenges. Here, we introduce and study a generalized covariant Navier-Stokes model for fluid flows driven by active stresses in nonplanar geometries. The analytical tractability of the theory is demonstrated through exact stationary solutions for the case of a spherical bubble geometry. Direct numerical simulations reveal a curvature-induced transition from a burst phase to an anomalous turbulent phase that differs distinctly from externally forced classical 2D Kolmogorov turbulence. This new type of active turbulence is characterized by the self-assembly of finite-size vortices into linked chains of antiferromagnetic order, which percolate through the entire fluid domain, forming an active dynamic network. The coherent motion of the vortex chain network provides an efficient mechanism for upward energy transfer from smaller to larger scales, presenting an alternative to the conventional energy cascade in classical 2D turbulence.

Hydrodynamic Choreographies of Microswimmers

We unveil orbital topologies of two nearby swimming microorganisms using an artificial microswimmer, called Quadroar. Depending on the initial conditions of the microswimmers, we find diverse families of attractors including dynamical equilibria, bound orbits, braids, and pursuit–evasion games. We also observe a hydrodynamic slingshot effect: a system of two hydrodynamically interacting swimmers moving along braids can advance in space faster than non-interacting swimmers that have the same actuation parameters and initial conditions as the interacting ones. Our findings suggest the existence of complex collective behaviors of microswimmers, from equilibrium to rapidly streaming states.

Activity-induced instability of phonons in 1D microfluidic crystals

One-dimensional crystals of passively-driven particles in microfluidic channels exhibit collective vibrational modes reminiscent of acoustic 'phonons'. These phonons are induced by the long-range hydrodynamic interactions among the particles and are neutrally stable at the linear level. Here, we analyze the effect of particle activity - self-propulsion - on the emergence and stability of these phonons. We show that the direction of wave propagation in active crystals is sensitive to the intensity of the background flow. We also show that activity couples, at the linear level, transverse waves to the particles' rotational motion, inducing a new mode of instability that persists in the limit of large background flow, or, equivalently, vanishingly small activity. We then report a new phenomenon of phonons switching back and forth between two adjacent crystals in both passively-driven and active systems, similar in nature to the wave switching observed in quantum mechanics, optical communication, and density stratified fluids. These findings could have implications for the design of commercial microfluidic systems and the self-assembly of passive and active micro-particles into one-dimensional structures.

Vortex dynamics and Lagrangian statistics in a model for active turbulence

Cellular suspensions such as dense bacterial flows exhibit a turbulence-like phase under certain conditions. We study this phenomenon of "active turbulence" statistically by using numerical tools. Following Wensink et al. (Proc. Natl. Acad. Sci. U.S.A. 109, 14308 (2012)), we model active turbulence by means of a generalized Navier-Stokes equation. Two-point velocity statistics of active turbulence, both in the Eulerian and the Lagrangian frame, is explored. We characterize the scale-dependent features of two-point statistics in this system. Furthermore, we extend this statistical study with measurements of vortex dynamics in this system. Our observations suggest that the large-scale statistics of active turbulence is close to Gaussian with sub-Gaussian tails.

Dynamically generated patterns in dense suspensions of active filaments

We use Langevin dynamics simulations to study dynamical behavior of a dense planar layer of active semiflexible filaments. Using the strength of active force and the thermal persistence length as parameters, we map a detailed phase diagram and identify several nonequilibrium phases in this system. In addition to a slowly flowing melt phase, we observe that, for sufficiently high activity, collective flow accompanied by signatures of local polar and nematic order appears in the system. This state is also characterized by strong density fluctuations. Furthermore, we identify an activity-driven crossover from this state of coherently flowing bundles of filaments to a phase with no global flow, formed by individual filaments coiled into rotating spirals. This suggests a mechanism where the system responds to activity by changing the shape of active agents, an effect with no analog in systems of active particles without internal degrees of freedom.

Derivation of a hydrodynamic theory for mesoscale dynamics in microswimmer suspensions

In this paper, we systematically derive a fourth-order continuum theory capable of reproducing mesoscale turbulence in a three-dimensional suspension of microswimmers. We start from overdamped Langevin equations for a generic microscopic model (pushers or pullers), which include hydrodynamic interactions on both small length scales (polar alignment of neighboring swimmers) and large length scales, where the solvent flow interacts with the order parameter field. The flow field is determined via the Stokes equation supplemented by an ansatz for the stress tensor. In addition to hydrodynamic interactions, we allow for nematic pair interactions stemming from excluded-volume effects. The results here substantially extend and generalize earlier findings [S. Heidenreich et al., Phys. Rev. E 94, 020601 (2016)2470-004510.1103/PhysRevE.94.020601], in which we derived a two-dimensional hydrodynamic theory. From the corresponding mean-field Fokker-Planck equation combined with a self-consistent closure scheme, we derive nonlinear field equations for the polar and the nematic order parameter, involving gradient terms of up to fourth order. We find that the effective microswimmer dynamics depends on the coupling between solvent flow and orientational order. For very weak coupling corresponding to a high viscosity of the suspension, the dynamics of mesoscale turbulence can be described by a simplified model containing only an effective microswimmer velocity.

Control of active turbulence through addressable soft interfaces

We present an experimental study of a kinesin/tubulin active nematic formed at different oil interfaces. By tuning the interfacial rheology of the contacting oil, we have been able to condition and control the seemingly chaotic motion that characterizes the self-sustained active flows in our preparations. The active nematic is inherently unstable and spontaneously develops defects from an initial homogeneous state. We show that the steady state and, in particular, the density and dynamics of the defects strongly depends on the rheology of the contacting oil. Using a smectic-A thermotropic liquid crystal as the oil phase, we pattern the interface thanks to the anisotropy of the shear viscosity in this material. The geometry of the active nematic adapts to the boundary conditions at the interface by changing from the so-called active turbulent regime to laminar flows along the easy flow directions. The latter can be either a lattice of self-assembled circular paths or reconfigurable homogeneous orientations that can be addressed by means of an external magnetic field. We show that, under all confinement conditions, the spatiotemporal modes exhibited by the active liquid are consistent with a single intrinsic length scale, which can be tuned by the material parameters, and obey basic topological requirements imposed on the defects that drive the active flows. Future control strategies, including a tunable depleting agent, are discussed.

Active turbulence in a gas of self-assembled spinners

Colloidal particles subject to an external periodic forcing exhibit complex collective behavior and self-assembled patterns. A dispersion of magnetic microparticles confined at the air-liquid interface and energized by a uniform uniaxial alternating magnetic field exhibits dynamic arrays of self-assembled spinners rotating in either direction. Here, we report on experimental and simulation studies of active turbulence and transport in a gas of self-assembled spinners. We show that the spinners, emerging as a result of spontaneous symmetry breaking of clock/counterclockwise rotation of self-assembled particle chains, generate vigorous vortical flows at the interface. An ensemble of spinners exhibits chaotic dynamics due to self-generated advection flows. The same-chirality spinners (clockwise or counterclockwise) show a tendency to aggregate and form dynamic clusters. Emergent self-induced interface currents promote active diffusion that could be tuned by the parameters of the external excitation field. Furthermore, the erratic motion of spinners at the interface generates chaotic fluid flow reminiscent of 2D turbulence. Our work provides insight into fundamental aspects of collective transport in active spinner materials and yields rules for particle manipulation at the microscale.

Flow properties and hydrodynamic interactions of rigid spherical microswimmers

We analyze a minimal model for a rigid spherical microswimmer and explore the consequences of its extended surface on the interplay between its self-propulsion and flow properties. The model is the first order representation of microswimmers, such as bacteria and algae, with rigid bodies and flexible propelling appendages. The flow field of such a microswimmer at finite distances significantly differs from that of a point-force (Stokeslet) dipole. For a suspension of microswimmers, we derive the grand mobility matrix that connects the motion of an individual swimmer to the active and passive forces and torques acting on all the swimmers. Our investigation of the mobility tensors reveals that hydrodynamic interactions among rigid-bodied microswimmers differ considerably from those among the corresponding point-force dipoles. Our results are relevant for the study of collective behavior of hydrodynamically interacting microswimmers by means of Stokesian dynamics simulations at moderate concentrations.

Intermittent turbulence in flowing bacterial suspensions

Dense suspensions of motile bacteria, possibly including the human gut microbiome, exhibit collective dynamics akin to those observed in classic, high Reynolds number turbulence with important implications for chemical and biological transport, yet this analogy has remained primarily qualitative. Here, we present experiments in which a dense suspension of Bacillus subtilis bacteria was flowed through microchannels and the velocity statistics of the flowing suspension were quantified using a recently developed velocimetry technique coupled with vortex identification methods. Observations revealed a robust intermittency phenomenon, whereby the average velocity profile of the suspension fluctuated between a plug-like flow and a parabolic flow profile. This intermittency is a hallmark of the onset of classic turbulence and Lagrangian tracking revealed that it here originates from the presence of transient vortices in the active, collective motion of the bacteria locally reinforcing the externally imposed flow. These results link together two entirely different manifestations of turbulence and show the potential of the microfluidic approach to mimic the environment characteristic of certain niches of the human microbiome.

Motility of active nematic films driven by "active anchoring"

We provide a minimal model for an active nematic film in contact with both a solid substrate and a passive isotropic fluid, and explore its dynamics in one and two dimensions using a combination of hybrid Lattice Boltzmann simulations and analytical calculations. By imposing nematic anchoring at the substrate while active flows induce a preferred alignment at the interface ("active anchoring"), we demonstrate that directed fluid flow spontaneously emerges in cases where the two anchoring types are opposing. In one dimension, our model reduces to an analogue of a loaded elastic column. Here, the transition from a stationary to a motile state is akin to the buckling bifurcation, but offers the possibility to reverse the flow direction for a given set of parameters and boundary conditions solely by changing initial conditions. The two-dimensional variant of our model allows for additional tangential instabilities, and it is found that undulations form in the interface above a threshold activity. Our results might be relevant for the design of active microfluidic geometries or curvature-guided self-assembly.

Taming active turbulence with patterned soft interfaces

Active matter embraces systems that self-organize at different length and time scales, often exhibiting turbulent flows apparently deprived of spatiotemporal coherence. Here, we use a layer of a tubulin-based active gel to demonstrate that the geometry of active flows is determined by a single length scale, which we reveal in the exponential distribution of vortex sizes of active turbulence. Our experiments demonstrate that the same length scale reemerges as a cutoff for a scale-free power law distribution of swirling laminar flows when the material evolves in contact with a lattice of circular domains. The observed prevalence of this active length scale can be understood by considering the role of the topological defects that form during the spontaneous folding of microtubule bundles. These results demonstrate an unexpected strategy for active systems to adapt to external stimuli, and provide with a handle to probe the existence of intrinsic length and time scales.

Active nematics consist of self-driven components that develop orientational order and turbulent flow. Here Guillamat et al. investigate an active nematic constrained in a quasi-2D geometrical setup and show that there exists an intrinsic length scale that determines the geometry in all forcing regimes.

Role of tumbling in bacterial swarming

Typical wild-type bacteria swimming in sparse suspensions exhibit a movement pattern called "run and tumble," characterized by straight trajectories (runs) interspersed by shorter, random reorientation (tumbles). This is achieved by rotating their flagella counterclockwise, or clockwise, respectively. The chemotaxis signaling network operates in controlling the frequency of tumbles, enabling navigation toward or away from desired regions in the medium. In contrast, while in dense populations, flagellated bacteria exhibit collective motion and form large dynamic clusters, whirls, and jets, with intricate dynamics that is fundamentally different than trajectories of sparsely swimming cells. Although collectively swarming cells do change direction at the level of the individual cell, often exhibiting reversals, it has been suggested that chemotaxis does not play a role in multicellular colony expansion, but the change in direction stems from clockwise flagellar rotation. In this paper, the effects of cell rotor switching (i.e., the ability to tumble) and chemotaxis on the collective statistics of swarming bacteria are studied experimentally in wild-type Bacillus subtilis and two mutants-one that does not tumble and one that tumbles independently of the chemotaxis system. We show that while several of the parameters examined are similar between the strains, other collective and individual characteristics are significantly different. The results demonstrate that tumbling and/or flagellar directional rotor switching has an important role on the dynamics of swarming, and imply that swarming models of self-propelled rods that do not take tumbling and/or rotor switching into account may be oversimplified.

Role of Correlations in the Collective Behavior of Microswimmer Suspensions

In this Letter, we study the collective behavior of a large number of self-propelled microswimmers immersed in a fluid. Using unprecedentedly 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 nonperturbative in the swimmer density. To provide an experimentally accessible measure of correlations, we calculate the diffusivity of passive tracers and reveal its nontrivial 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.

Intrinsic flow structure and multifractality in two-dimensional bacterial turbulence

The active interaction between the bacteria and fluid generates turbulent structures even at zero Reynolds number. The velocity of such a flow obtained experimentally has been quantitatively investigated based on streamline segment analysis. There is a clear transition at about 16 times the organism body length separating two different scale regimes, which may be attributed to the different influence of the viscous effect. Surprisingly the scaling extracted from the streamline segment indicates the existence of scale similarity even at the zero Reynolds number limit. Moreover, the multifractal feature can be quantitatively described via a lognormal formula with the Hurst number H=0.76 and the intermittency parameter μ=0.20, which is coincidentally in agreement with the three-dimensional hydrodynamic turbulence result. The direction of cascade is measured via the filter-space technique. An inverse energy cascade is confirmed. For the enstrophy, a forward cascade is observed when r/R≤3, and an inverse one is observed when r/R>3, where r and R are the separation distance and the bacteria body size, respectively. Additionally, the lognormal statistics is verified for the coarse-grained energy dissipation and enstrophy, which supports the lognormal formula to fit the measured scaling exponent.

Anomalous Fluctuations in the Orientation and Velocity of Swarming Bacteria

Simultaneous acquisition of phase-contrast light microscopy and fluorescently labeled bacteria, moving within a dense swarm, reveals the intricate interactions between cells and the collective flow around them. By comparing wild-type and immotile cells embedded in a dense wild-type swarm, the effect of the active thrust generated by the flagella can be singled out. It is shown that while the distribution of angles among cell velocity, cell orientation, and the local flow around it is Gaussian-like for immotile bacteria, wild-type cells exhibit anomalous non-Gaussian deviations and are able to move in trajectories perpendicular to the collective flow. Thus, cells can maneuver or switch between local streams and jets. A minimal model describing bacteria as hydrodynamic force dipoles shows that steric effects, hydrodynamics interactions, and local alignments all have to be taken into account to explain the observed dynamics. These findings shed light on the physical mechanisms underlying bacterial swarming and the balance between individual and collective dynamics.

Dancing disclinations in confined active nematics

The spontaneous emergence of collective flows is a generic property of active fluids and often leads to chaotic flow patterns characterised by swirls, jets, and topological disclinations in their orientation field. However, the ability to achieve structured flows and ordered disclinations is of particular importance in the design and control of active systems. By confining an active nematic fluid within a channel, we find a regular motion of disclinations, in conjunction with a well defined and dynamic vortex lattice. As pairs of moving disclinations travel through the channel, they continually exchange partners producing a dynamic ordered state, reminiscent of Ceilidh dancing. We anticipate that this biomimetic ability to self-assemble organised topological disclinations and dynamically structured flow fields in engineered geometries will pave the road towards establishing new active topological microfluidic devices.

Onset of meso-scale turbulence in active nematics

Meso-scale turbulence is an innate phenomenon, distinct from inertial turbulence, that spontaneously occurs at low Reynolds number in fluidized biological systems. This spatiotemporal disordered flow radically changes nutrient and molecular transport in living fluids and can strongly affect the collective behaviour in prominent biological processes, including biofilm formation, morphogenesis and cancer invasion. Despite its crucial role in such physiological processes, understanding meso-scale turbulence and any relation to classical inertial turbulence remains obscure. Here we show how the motion of active matter along a micro-channel transitions to meso-scale turbulence through the evolution of locally disordered patches (active puffs) from an ordered vortex-lattice flow state. We demonstrate that the stationary critical exponents of this transition to meso-scale turbulence in a channel coincide with the directed percolation universality class. This finding bridges our understanding of the onset of low-Reynolds-number meso-scale turbulence and traditional scale-invariant turbulence in confinement.

Elastic turbulence in entangled semi-dilute DNA solutions measured with optical coherence tomography velocimetry

The flow instabilities of solutions of high molecular weight DNA in the entangled semi-dilute concentration regime were investigated using optical coherence tomography velocimetry, a technique that provides high spatial (probe volumes of 3.4 pL) and temporal resolution (sub μs) information on the flow behaviour of complex fluids in a rheometer. The velocity profiles of the opaque DNA solutions (high and low salt) were measured as a function of the distance across the gap of a parallel plate rheometer, and their evolution over time was measured. At lower DNA concentrations and low shear rates, the velocity fluctuations were well described by Gaussian functions and the velocity gradient was uniform across the rheometer gap, which is expected for Newtonian flows. As the DNA concentration and shear rate were increased there was a stable wall slip regime followed by an evolving wall slip regime, which is finally followed by the onset of elastic turbulence. Strain localization (shear banding) is observed on the boundaries of the flows at intermediate shear rates, but decreases in the high shear elastic turbulence regime, where bulk strain localization occurs. A dynamic phase diagram for non-linear flow was created to describe the different behaviours.

Active matter logic for autonomous microfluidics

Chemically or optically powered active matter plays an increasingly important role in materials design, but its computational potential has yet to be explored systematically. The competition between energy consumption and dissipation imposes stringent physical constraints on the information transport in active flow networks, facilitating global optimization strategies that are not well understood. Here, we combine insights from recent microbial experiments with concepts from lattice-field theory and non-equilibrium statistical mechanics to introduce a generic theoretical framework for active matter logic. Highlighting conceptual differences with classical and quantum computation, we demonstrate how the inherent non-locality of incompressible active flow networks can be utilized to construct universal logical operations, Fredkin gates and memory storage in set-reset latches through the synchronized self-organization of many individual network components. Our work lays the conceptual foundation for developing autonomous microfluidic transport devices driven by bacterial fluids, active liquid crystals or chemically engineered motile colloids.

Influence of Physical Effects on the Swarming Motility of Pseudomonas aeruginosa

Many species of bacteria can spread over a moist surface via a particular form of collective motion known as "surface swarming". This form of motility is typically studied by inoculating bacteria on a gel formed by 0.4-1.5% agar, which contains essential nutrients for their growth and proliferation. Using Pseudomonas aeruginosa and its pili-less mutant, ΔPilA, we investigate physical factors that either facilitate or restrict the swarming motility, measured by the rate of increase in area covered by a spreading bacterial colony, i.e., a swarm. The wild-type colony spreads over the agar surface in highly branched structures. The pili-less mutant fills up the area more fully as it spreads, but it also produces numerous and fragmented branches, or tendrils, at the swarm front. Whereas additional surfactants enhance swarming, increasing the agar percentage, adding extra salt or sugar or incorporating viscous agents in the agar matrix all decrease swarming, supporting the conclusion that swarming motility is restricted by the surface tension at the swarm front and swarm growth is limited by the rate of water supply from within the agar gel. The physical basis elaborated through this study provides a useful framework for understanding the swarming behavior of numerous species 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.

Non-Gaussian limit fluctuations in active swimmer suspensions

We investigate the hydrodynamic fluctuations in suspensions of swimming microorganisms (Chlamydomonas) by observing the probe particles dispersed in the media. Short-term fluctuations of probe particles were superdiffusive and displayed heavily tailed non-Gaussian distributions. The analytical theory that explains the observed distribution was derived by summing the power-law-decaying hydrodynamic interactions from spatially distributed field sources (here, swimming microorganisms). The summing procedure, which we refer to as the physical limit operation, is applicable to a variety of physical fluctuations to which the classical central limiting theory does not apply. Extending the analytical formula to compare to experiments in active swimmer suspensions, we show that the non-Gaussian shape of the observed distribution obeys the analytic theory concomitantly with independently determined parameters such as the strength of force generations and the concentration of Chlamydomonas. Time evolution of the distributions collapsed to a single master curve, except for their extreme tails, for which our theory presents a qualitative explanation. Investigations thereof and the complete agreement with theoretical predictions revealed broad applicability of the formula to dispersions of active sources of fluctuations.

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.

Hydrodynamic interaction between particles near elastic interfaces

We present an analytical calculation of the hydrodynamic interaction between two spherical particles near an elastic interface such as a cell membrane. The theory predicts the frequency dependent self- and pair-mobilities accounting for the finite particle size up to the 5th order in the ratio between particle diameter and wall distance as well as between diameter and interparticle distance. We find that particle motion towards a membrane with pure bending resistance always leads to mutual repulsion similar as in the well-known case of a hard-wall. In the vicinity of a membrane with shearing resistance, however, we observe an attractive interaction in a certain parameter range which is in contrast to the behavior near a hard wall. This attraction might facilitate surface chemical reactions. Furthermore, we show that there exists a frequency range in which the pair-mobility for perpendicular motion exceeds its bulk value, leading to short-lived superdiffusive behavior. Using the analytical particle mobilities we compute collective and relative diffusion coefficients. The appropriateness of the approximations in our analytical results is demonstrated by corresponding boundary integral simulations which are in excellent agreement with the theoretical predictions.

Lattice-Boltzmann simulations of microswimmer-tracer interactions

Hydrodynamic interactions in systems composed of self-propelled particles, such as swimming microorganisms and passive tracers, have a significant impact on the tracer dynamics compared to the equivalent "dry" sample. However, such interactions are often difficult to take into account in simulations due to their computational cost. Here, we perform a systematic investigation of swimmer-tracer interaction using an efficient force-counterforce-based lattice-Boltzmann (LB) algorithm [De Graaf et al., J. Chem. Phys. 144, 134106 (2016)JCPSA60021-960610.1063/1.4944962] in order to validate its ability to capture the relevant low-Reynolds-number physics. We show that the LB algorithm reproduces far-field theoretical results well, both in a system with periodic boundary conditions and in a spherical cavity with no-slip walls, for which we derive expressions here. The force-lattice coupling of the LB algorithm leads to a "smearing out" of the flow field, which strongly perturbs the tracer trajectories at close swimmer-tracer separations, and we analyze how this effect can be accurately captured using a simple renormalized hydrodynamic theory. Finally, we show that care must be taken when using LB algorithms to simulate systems of self-propelled particles, since its finite momentum transport time can lead to significant deviations from theoretical predictions based on Stokes flow. These insights should prove relevant to the future study of large-scale microswimmer suspensions using these methods.

Geometric control of active collective motion

Recent experimental studies have shown that confinement can profoundly affect self-organization in semi-dilute active suspensions, leading to striking features such as the formation of steady and spontaneous vortices in circular domains and the emergence of unidirectional pumping motions in periodic racetrack geometries. Motivated by these findings, we analyze the two-dimensional dynamics in confined suspensions of active self-propelled swimmers using a mean-field kinetic theory where conservation equations for the particle configurations are coupled to the forced Navier-Stokes equations for the self-generated fluid flow. In circular domains, a systematic exploration of the parameter space casts light on three distinct states: equilibrium with no flow, stable vortex, and chaotic motion, and the transitions between these are explained and predicted quantitatively using a linearized theory. In periodic racetracks, similar transitions from equilibrium to net pumping to traveling waves to chaos are observed in agreement with experimental observations and are also explained theoretically. Our results underscore the subtle effects of geometry on the morphology and dynamics of emerging patterns in active suspensions and pave the way for the control of active collective motion in microfluidic devices.

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.

Classical Nucleation Theory Description of Active Colloid Assembly

Nonaligning self-propelled particles with purely repulsive excluded volume interactions undergo athermal motility-induced phase separation into a dilute gas and a dense cluster phase. Here, we use enhanced sampling computational methods and analytic theory to examine the kinetics of formation of the dense phase. Despite the intrinsically nonequilibrium nature of the phase transition, we show that the kinetics can be described using an approach analogous to equilibrium classical nucleation theory, governed by an effective free energy of cluster formation with identifiable bulk and surface terms. The theory captures the location of the binodal, nucleation rates as a function of supersaturation, and the cluster size distributions below the binodal, while discrepancies in the metastable region reveal additional physics about the early stages of active crystal formation. The success of the theory shows that a framework similar to equilibrium thermodynamics can be obtained directly from the microdynamics of an active system, and can be used to describe the kinetics of evolution toward nonequilibrium steady states.

Correlation lengths in hydrodynamic models of active nematics

We examine the scaling with activity of the emergent length scales that control the nonequilibrium dynamics of an active nematic liquid crystal, using two popular hydrodynamic models that have been employed in previous studies. In both models we find that the chaotic spatio-temporal dynamics in the regime of fully developed active turbulence is controlled by a single active scale determined by the balance of active and elastic stresses, regardless of whether the active stress is extensile or contractile in nature. The observed scaling of the kinetic energy and enstrophy with activity is consistent with our single-length scale argument and simple dimensional analysis. Our results provide a unified understanding of apparent discrepancies in the previous literature and demonstrate that the essential physics is robust to the choice of model.

Collective Motion of Microorganisms in a Viscoelastic Fluid

We study the collective motion of a suspension of rodlike microswimmers in a two-dimensional film of viscoelastic fluids. We find that the fluid elasticity has a small effect on a suspension of pullers, while it significantly affects the pushers. The attraction and orientational ordering of the pushers are enhanced in viscoelastic fluids. The induced polymer stresses break down the large-scale flow structures and suppress velocity fluctuations. In addition, the energy spectra and induced mixing in the suspension of pushers are greatly modified by fluid elasticity.

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.

Active Brownian motion of emulsion droplets: Coarsening dynamics at the interface and rotational diffusion

A micron-sized droplet of bromine water immersed in a surfactant-laden oil phase can swim (S. Thutupalli, R. Seemann, S. Herminghaus, New J. Phys. 13 073021 (2011). The bromine reacts with the surfactant at the droplet interface and generates a surfactant mixture. It can spontaneously phase-separate due to solutocapillary Marangoni flow, which propels the droplet. We model the system by a diffusion-advection-reaction equation for the mixture order parameter at the interface including thermal noise and couple it to fluid flow. Going beyond previous work, we illustrate the coarsening dynamics of the surfactant mixture towards phase separation in the axisymmetric swimming state. Coarsening proceeds in two steps: an initially slow growth of domain size followed by a nearly ballistic regime. On larger time scales thermal fluctuations in the local surfactant composition initiates random changes in the swimming direction and the droplet performs a persistent random walk, as observed in experiments. Numerical solutions show that the rotational correlation time scales with the square of the inverse noise strength. We confirm this scaling by a perturbation theory for the fluctuations in the mixture order parameter and thereby identify the active emulsion droplet as an active Brownian particle.

Stability and dynamics of anisotropically tumbling chemotactic swimmers

Microswimmers such as bacteria perform random walks known as run-and-tumbles to move up chemoattractant gradients and as a result aggregate with others. It is also known that such micro-swimmers can self-organize into macroscopic patterns due to interactions with neighboring cells through the fluidic environment they live in. While the pattern formation resulting from chemotactic and hydrodynamic interactions separately and together have been previously investigated, the effect of the anisotropy in the tumbles of microswimmers has been unexplored. Here we show through linear analysis and full nonlinear simulations that the slight anisotropy in the individual swimmer tumbles can alter the collective pattern formation in nontrivial ways. We show that tumbling anisotropy diminishes the magnitude of the chemotactic aggregates but may result in more such aggregation peaks.

Intermittency measurement in two-dimensional bacterial turbulence

In this paper, an experimental velocity database of a bacterial collective motion, e.g., Bacillus subtilis, in turbulent phase with volume filling fraction 84% provided by Professor Goldstein at Cambridge University (UK), was analyzed to emphasize the scaling behavior of this active turbulence system. This was accomplished by performing a Hilbert-based methodology analysis to retrieve the scaling property without the β-limitation. A dual-power-law behavior separated by the viscosity scale ℓ_{ν} was observed for the qth-order Hilbert moment L_{q}(k). This dual-power-law belongs to an inverse-cascade since the scaling range is above the injection scale R, e.g., the bacterial body length. The measured scaling exponents ζ(q) of both the small-scale (k>k_{ν}) and large-scale (k<k_{ν}) motions are convex, showing the multifractality. A log-normal formula was put forward to characterize the multifractal intensity. The measured intermittency parameters are μ_{S}=0.26 and μ_{L}=0.17, respectively, for the small- and large-scale motions. It implies that the former cascade is more intermittent than the latter one, which is also confirmed by the corresponding singularity spectrum f(α) versus α. Comparison with the conventional two-dimensional Ekman-Navier-Stokes equation, a continuum model indicates that the origin of the multifractality could be a result of some additional nonlinear interaction terms, which deservers a more careful investigation.

Dynamics of a bacterial flagellum under reverse rotation

To initiate tumbling of an E. coli, one of the helical flagella reverses its sense of rotation. It then transforms from its normal form first to the transient semicoiled state and subsequently to the curly-I state. The dynamics of polymorphism is effectively modeled by describing flagellar elasticity through an extended Kirchhoff free energy. However, the complete landscape of the free energy remains undetermined because the ground state energies of the polymorphic forms are not known. We investigate how variations in these ground state energies affect the dynamics of a reversely rotated flagellum of a swimming bacterium. We find that the flagellum exhibits a number of distinct dynamical states and comprehensively summarize them in a state diagram. As a result, we conclude that tuning the landscape of the extended Kirchhoff free energy alone cannot generate the intermediate full-length semicoiled state. However, our model suggests an ad hoc method to realize the sequence of polymorphic states as observed for a real bacterium. Since the elastic properties of bacterial flagella are similar, our findings can easily be extended to other peritrichous bacteria.

Dynamics of inert spheres in active suspensions of micro-rotors

Inert particles suspended in active fluids of self-propelled particles are known to often exhibit enhanced diffusion and novel coherent structures. Here we numerically investigate the dynamical behavior and self-organization in a system consisting of passive and actively rotating spheres of the same size. The particles interact through direct collisions and the fluid flows generated as they move. In the absence of passive particles, three states emerge in a binary mixture of spinning spheres depending on particle fraction: a dilute gas-like state where the rotors move chaotically, a phase-separated state where like-rotors move in lanes or vortices, and a jammed state where crystals continuously assemble, melt and move (K. Yeo, E. Lushi, and P. M. Vlahovska, Phys. Rev. Lett., 2015, 114, 188301). Passive particles added to the rotor suspension modify the system dynamics and pattern formation: while states identified in the pure active suspension still emerge, they occur at different densities and mixture proportions. The dynamical behavior of the inert particles is also non-trivially dependent on the system composition.

Stochastic cycle selection in active flow networks

Active biological flow networks pervade nature and span a wide range of scales, from arterial blood vessels and bronchial mucus transport in humans to bacterial flow through porous media or plasmodial shuttle streaming in slime molds. Despite their ubiquity, little is known about the self-organization principles that govern flow statistics in such nonequilibrium networks. Here we connect concepts from lattice field theory, graph theory, and transition rate theory to understand how topology controls dynamics in a generic model for actively driven flow on a network. Our combined theoretical and numerical analysis identifies symmetry-based rules that make it possible to classify and predict the selection statistics of complex flow cycles from the network topology. The conceptual framework developed here is applicable to a broad class of biological and nonbiological far-from-equilibrium networks, including actively controlled information flows, and establishes a correspondence between active flow networks and generalized ice-type models.

Active micromachines: Microfluidics powered by mesoscale turbulence

Dense active matter, from bacterial suspensions and microtubule bundles driven by motor proteins to cellular monolayers and synthetic Janus particles, is characterized by mesoscale turbulence, which is the emergence of chaotic flow structures. By immersing an ordered array of symmetric rotors in an active fluid, we introduce a microfluidic system that exploits spontaneous symmetry breaking in mesoscale turbulence to generate work. The lattice of rotors self-organizes into a spin state where neighboring discs continuously rotate in permanent alternating directions due to combined hydrodynamic and elastic effects. Our virtual prototype demonstrates a new research direction for the design of micromachines powered by the nematohydrodynamic properties of active turbulence.

Lagrangian model of copepod dynamics: Clustering by escape jumps in turbulence

Planktonic copepods are small crustaceans that have the ability to swim by quick powerful jumps. Such an aptness is used to escape from high shear regions, which may be caused either by flow perturbations, produced by a large predator (i.e., fish larvae), or by the inherent highly turbulent dynamics of the ocean. Through a combined experimental and numerical study, we investigate the impact of jumping behavior on the small-scale patchiness of copepods in a turbulent environment. Recorded velocity tracks of copepods displaying escape response jumps in still water are here used to define and tune a Lagrangian copepod (LC) model. The model is further employed to simulate the behavior of thousands of copepods in a fully developed hydrodynamic turbulent flow obtained by direct numerical simulation of the Navier-Stokes equations. First, we show that the LC velocity statistics is in qualitative agreement with available experimental observations of copepods in turbulence. Second, we quantify the clustering of LC, via the fractal dimension D_{2}. We show that D_{2} can be as low as ∼2.3 and that it critically depends on the shear-rate sensitivity of the proposed LC model, in particular it exhibits a minimum in a narrow range of shear-rate values. We further investigate the effect of jump intensity, jump orientation, and geometrical aspect ratio of the copepods on the small-scale spatial distribution. At last, possible ecological implications of the observed clustering on encounter rates and mating success are discussed.