Hydrodynamic length-scale selection in microswimmer suspensions

Oscillating edge current in polar active fluid

Dense bacterial suspensions exhibit turbulent behavior called bacterial turbulence. The behavior of the bulk unconstrained bacterial turbulence is described well by the Toner-Tu-Swift-Hohenberg (TTSH) equation for the velocity field. However, it remains unclear how we should treat boundary conditions on bacterial turbulence in contact with some boundaries (e.g., solid walls). To be more specific, although the importance of the edge current, the flow along the boundary, has been demonstrated in several experimental studies on confined bacterial suspensions, previous numerical studies based on the TTSH equation employ nonslip boundary conditions and do not seem to properly describe the behavior of bacteria near the boundaries. In this paper, we impose a slip boundary condition on the TTSH equation to describe the bacterial motion at boundaries. We develop a method to implement the slip boundary condition. Using this method, we have successfully produced edge current and discovered that the direction of the edge current temporally oscillates. The oscillation can be attributable to the advection term in the TTSH equation. Our paper demonstrates that boundary conditions could play an important role in the collective dynamics of active systems.

Self-enhanced mobility enables vortex pattern formation in living matter

Ranging from subcellular organelle biogenesis to embryo development, the formation of self-organized structures is a hallmark of living systems. Whereas the emergence of ordered spatial patterns in biology is often driven by intricate chemical signalling that coordinates cellular behaviour and differentiation1-4, purely physical interactions can drive the formation of regular biological patterns such as crystalline vortex arrays in suspensions of spermatozoa5 and bacteria6. Here we discovered a new route to self-organized pattern formation driven by physical interactions, which creates large-scale regular spatial structures with multiscale ordering. Specifically we found that dense bacterial living matter spontaneously developed a lattice of mesoscale, fast-spinning vortices; these vortices each consisted of around 104-105 motile bacterial cells and were arranged in space at greater than centimetre scale and with apparent hexagonal order, whereas individual cells in the vortices moved in coordinated directions with strong polar and vortical order. Single-cell tracking and numerical simulations suggest that the phenomenon is enabled by self-enhanced mobility in the system-that is, the speed of individual cells increasing with cell-generated collective stresses at a given cell density. Stress-induced mobility enhancement and fluidization is prevalent in dense living matter at various scales of length7-9. Our findings demonstrate that self-enhanced mobility offers a simple physical mechanism for pattern formation in living systems and, more generally, in other active matter systems10 near the boundary of fluid- and solid-like behaviours11-17.

Stochastic motion under nonlinear friction representing shear thinning

We study stochastic motion under a nonlinear frictional force that levels off with increasing velocity. Specifically, our frictional force is of the so-called Coulomb-tanh type. At small speed, it increases approximately linearly with velocity, while at large speed, it approaches a constant magnitude, similarly to solid (dry, Coulomb) friction. In one spatial dimension, a formal analogy between the associated Fokker-Planck equation and the Schrödinger equation for a quantum mechanical oscillator in a nonharmonic Pöschl-Teller potential is revealed. Then, the stationary velocity statistics can be treated analytically. From such analytical considerations, we determine associated diffusion coefficients, which we confirm by agent-based simulations. Moreover, from such simulations and from numerically solving the associated Fokker-Planck equation, we find that the spatial distribution function, starting from an initial Gaussian shape, develops tails that appear exponential at intermediate timescales. At small magnitudes of stochastic driving, the velocity distribution resembles the case of linear friction, while at large magnitudes, it rather approaches the case of solid (dry, Coulomb) friction. The same is true for diffusion coefficients. In a certain sense thus interpolating between the two extreme cases of linear friction and solid (dry, Coulomb) friction, our approach should be useful to describe several cases of practical relevance. For instance, a reduced increase in friction with increasing relative speed is typical of shear-thinning behavior. Therefore, driven motion in shear-thinning environments is one specific example to which our description may be applied.

Flocking turbulence of microswimmers in confined domains

We extensively study the Toner-Tu-Swift-Hohenberg model of motile active matter by means of direct numerical simulations in a two-dimensional confined domain. By exploring the space of parameters of the model we investigate the emergence of a new state of active turbulence which occurs when the aligning interactions and the self-propulsion of the swimmers are strong. This regime of flocking turbulence is characterized by a population of few strong vortices, each surrounded by an island of coherent flocking motion. The energy spectrum of flocking turbulence displays a power-law scaling with an exponent which depends weakly on the model parameters. By increasing the confinement we observe that the system, after a long transient characterized by power-law-distributed transition times, switches to the ordered state of a single giant vortex.

Inertial and geometrical effects of self-propelled elliptical Brownian particles

Active particles that self-propel by transforming energy into mechanical motion represent a growing area of research in mathematics, physics, and chemistry. Here we investigate the dynamics of nonspherical inertial active particles moving in a harmonic potential, introducing geometric parameters which take into account the role of eccentricity for nonspherical particles. A comparison between the overdamped and underdamped models for elliptical particles is performed. The model of overdamped active Brownian motion has been used to describe most of the basic aspects of micrometer-sized particles moving in a liquid ("microswimmers"). We consider active particles by extending the active Brownian motion model to incorporate translation and rotation inertia and account for the role of eccentricity. We show how the overdamped and the underdamped models behave in the same way for small values of activity (Brownian case) if eccentricity is equal to zero, but increasing eccentricity leads the two dynamics to substantially depart from each other-in particular, the action of a torque induced by external forces, induced a marked difference close to the walls of the domain if eccentricity is high. Effects induced by inertia include an inertial delay time of the self-propulsion direction from the particle velocity, and the differences between the overdamped and underdamped systems are particularly evident in the first and second moments of the particle velocities. Comparison with the experimental results of vibrated granular particles shows good agreement and corroborates the notion that self-propelling massive particles moving in gaseous media are dominated by inertial effects.

Giant vortex dynamics in confined bacterial turbulence

We report the numerical evidence of a new state of bacterial turbulence in confined domains. By means of extensive numerical simulations of the Toner-Tu-Swift-Hohenberg model for dense bacterial suspensions in circular geometry, we discover the formation a stable, ordered state in which the angular momentum symmetry is broken. This is achieved by self-organization of a turbulent-like flow into a single, giant vortex of the size of the domain. The giant vortex is surrounded by an annular region close to the boundary, characterized by small-scale, radial vorticity streaks. The average radial velocity profile of the vortex is found to be in agreement with a simple analytical prediction. We also provide an estimate of the temporal and spatial scales of a suitable experimental setup comparable with our numerical findings.

Emergent collective dynamics of pusher and puller squirmer rods: swarming, clustering, and turbulence

We study the interplay of steric and hydrodynamic interactions in suspensions of elongated microswimmers by simulating the full hydrodynamics of squirmer rods in the quasi two-dimensional geometry of a Hele-Shaw cell. To create pusher or puller-type squirmer rods, we concentrate the surface slip-velocity field more to the back or to the front of the rod and thereby are able to tune the rod's force-dipole strength. We study a wide range of aspect ratios and area fractions and provide corresponding state diagrams. The flow field of pusher-type squirmer rods destabilizes ordered structures and favors the disordered state at small area fractions and aspect ratios. Only when steric interactions become relevant, we observe a turbulent and dynamic cluster state, while for large aspect ratios a single swarm and jammed cluster occurs. The power spectrum of the turbulent state shows two distinct energy cascades at small and large wave numbers with power-law scaling and non-universal exponents. Pullers show a strong tendency to form swarms instead of the disordered state found for neutral and pusher rods. At large area fractions a dynamic cluster is observed and at larger aspect ratio a single swarm or jammed cluster occurs.

Statistics for an object actively driven by spontaneous symmetry breaking into reversible directions

Propulsion of otherwise passive objects is achieved by mechanisms of active driving. We concentrate on cases in which the direction of active drive is subject to spontaneous symmetry breaking. In our case, this direction will be maintained until a large enough impulse by an additional stochastic force reverses it. Examples may be provided by self-propelled droplets, gliding bacteria stochastically reversing their propulsion direction, or nonpolar vibrated hoppers. The magnitude of active forcing is regarded as constant, and we include the effect of inertial contributions. Interestingly, this situation can formally be mapped to stochastic motion under (dry, solid) Coulomb friction, however, with a negative friction parameter. Diffusion coefficients are calculated by formal mapping to the situation of a quantum-mechanical harmonic oscillator exposed to an additional repulsive delta-potential. Results comprise a ditched or double-peaked velocity distribution and spatial statistics showing outward propagating maxima when starting from initially concentrated arrangements.

Bacterial active matter

Bacteria are among the oldest and most abundant species on Earth. Bacteria successfully colonize diverse habitats and play a significant role in the oxygen, carbon, and nitrogen cycles. They also form human and animal microbiota and may become sources of pathogens and a cause of many infectious diseases. Suspensions of motile bacteria constitute one of the most studied examples of active matter: a broad class of non-equilibrium systems converting energy from the environment (e.g., chemical energy of the nutrient) into mechanical motion. Concentrated bacterial suspensions, often termed active fluids, exhibit complex collective behavior, such as large-scale turbulent-like motion (so-called bacterial turbulence) and swarming. The activity of bacteria also affects the effective viscosity and diffusivity of the suspension. This work reports on the progress in bacterial active matter from the physics viewpoint. It covers the key experimental results, provides a critical assessment of major theoretical approaches, and addresses the effects of visco-elasticity, liquid crystallinity, and external confinement on collective behavior in bacterial suspensions.

Analytical approach to chiral active systems: suppressed phase separation of interacting Brownian circle swimmers

We consider chirality in active systems by exemplarily studying the phase behavior of planar systems of interacting Brownian circle swimmers with a spherical shape. For this purpose, we derive a predictive field theory that is able to describe the collective dynamics of circle swimmers. The theory yields a mapping between circle swimmers and noncircling active Brownian particles and predicts that the angular propulsion of the particles leads to a suppression of their motility-induced phase separation, being in line with recent simulation results. In addition, the theory provides analytical expressions for the spinodal corresponding to the onset of motility-induced phase separation and the associated critical point as well as for their dependence on the angular propulsion of the circle swimmers. We confirm our findings by Brownian dynamics simulations. The agreement between results from theory and simulations is found to be good.

Ising-like Critical Behavior of Vortex Lattices in an Active Fluid

Turbulent vortex structures emerging in bacterial active fluids can be organized into regular vortex lattices by weak geometrical constraints such as obstacles. Here we show, using a continuum-theoretical approach, that the formation and destruction of these patterns exhibit features of a continuous second-order equilibrium phase transition, including long-range correlations, divergent susceptibility, and critical slowing down. The emerging vorticity field can be mapped onto a two-dimensional (2D) Ising model with antiferromagnetic nearest-neighbor interactions by coarse graining. The resulting effective temperature is found to be proportional to the strength of the nonlinear advection in the continuum model.

Turbulence-induced clustering in compressible active fluids

We study a novel phase of active polar fluids, which is characterized by the continuous creation and destruction of dense clusters due to self-sustained turbulence. This state arises due to the interplay between self-advection of the aligned swimmers and their defect topology. The typical cluster size is determined by the characteristic vortex size. Our results are obtained by investigating a continuum model of compressible polar active fluids, which incorporates typical experimental observations in bacterial suspensions, in particular a non-monotone dependence of speed on density.

Multi-particle collision dynamics with a non-ideal equation of state. II. Collective dynamics of elongated squirmer rods

Simulations of flow fields around microscopic objects typically require methods that both solve the Navier-Stokes equations and also include thermal fluctuations. One such method popular in the field of soft-matter physics is the particle-based simulation method of multi-particle collision dynamics (MPCD). However, in contrast to the typically incompressible real fluid, the fluid of the traditional MPCD methods obeys the ideal-gas equation of state. This can be problematic because most fluid properties strongly depend on the fluid density. In a recent article, we proposed an extended MPCD algorithm and derived its non-ideal equation of state and an expression for the viscosity. In the present work, we demonstrate its accuracy and efficiency for the simulations of the flow fields of single squirmers and of the collective dynamics of squirmer rods. We use two exemplary squirmer-rod systems for which we compare the outcome of the extended MPCD method to the well-established MPCD version with an Andersen thermostat. First, we explicitly demonstrate the reduced compressibility of the MPCD fluid in a cluster of squirmer rods. Second, for shorter rods, we show the interesting result that in simulations with the extended MPCD method, dynamic swarms are more pronounced and have a higher polar order. Finally, we present a thorough study of the state diagram of squirmer rods moving in the center plane of a Hele-Shaw geometry. From a small to large aspect ratio and density, we observe a disordered state, dynamic swarms, a single swarm, and a jammed cluster, which we characterize accordingly.

Emergence and melting of active vortex crystals

Melting of two-dimensional (2D) equilibrium crystals is a complex phenomenon characterized by the sequential loss of positional and orientational order. In contrast to passive systems, active crystals can self-assemble and melt into an active fluid by virtue of their intrinsic motility and inherent non-equilibrium stresses. Currently, the non-equilibrium physics of active crystallization and melting processes is not well understood. Here, we establish the emergence and investigate the melting of self-organized vortex crystals in 2D active fluids using a generalized Toner-Tu theory. Performing extensive hydrodynamic simulations, we find rich transition scenarios. On small domains, we identify a hysteretic transition as well as a transition featuring temporal coexistence of active vortex lattices and active turbulence, both of which can be controlled by self-propulsion and active stresses. On large domains, an active vortex crystal with solid order forms within the parameter range corresponding to active vortex lattices. The melting of this crystal proceeds through an intermediate hexatic phase. Generally, these results highlight the differences and similarities between crystalline phases in active fluids and their equilibrium counterparts.

Anomalous Diffusion and Lévy Walks Distinguish Active from Inertial Turbulence

Bacterial swarms display intriguing dynamical states like active turbulence. Now, using a hydrodynamic model, we show that such dense active suspensions manifest superdiffusion, via Lévy walks, which masquerades as a crossover from ballistic to diffusive scaling in measurements of mean-squared displacements, and is tied to the emergence of hitherto undetected oscillatory streaks in the flow. Thus, while laying the theoretical framework of an emergent advantageous strategy in the collective behavior of microorganisms, our Letter underlines the essential differences between active and inertial turbulence.

Towards an analytical description of active microswimmers in clean and in surfactant-covered drops

Geometric confinements are frequently encountered in the biological world and strongly affect the stability, topology, and transport properties of active suspensions in viscous flow. Based on a far-field analytical model, the low-Reynolds-number locomotion of a self-propelled microswimmer moving inside a clean viscous drop or a drop covered with a homogeneously distributed surfactant, is theoretically examined. The interfacial viscous stresses induced by the surfactant are described by the well-established Boussinesq-Scriven constitutive rheological model. Moreover, the active agent is represented by a force dipole and the resulting fluid-mediated hydrodynamic couplings between the swimmer and the confining drop are investigated. We find that the presence of the surfactant significantly alters the dynamics of the encapsulated swimmer by enhancing its reorientation. Exact solutions for the velocity images for the Stokeslet and dipolar flow singularities inside the drop are introduced and expressed in terms of infinite series of harmonic components. Our results offer useful insights into guiding principles for the control of confined active matter systems and support the objective of utilizing synthetic microswimmers to drive drops for targeted drug delivery applications.

Squirmer rods as elongated microswimmers: flow fields and confinement

Microswimmers or active elements, such as bacteria and active filaments, have an elongated shape, which determines their individual and collective dynamics. There is still a need to identify what role long-range hydrodynamic interactions play in their fascinating dynamic structure formation. We construct rods of different aspect ratios using several spherical squirmer model swimmers. With the help of the mesoscale simulation method of multi-particle collision dynamics we analyze the flow fields of these squirmer rods both in a bulk fluid and in Hele-Shaw geometries of different slab widths. Based on the hydrodynamic multipole expansion either for bulk or confinement between two parallel plates, we categorize the different multipole contributions of neutral as well as pusher-type squirmer rods. We demonstrate how confinement alters the radial decay of the flow fields for a given force or source multipole moment compared to the bulk fluid.

Active phase field crystal systems with inertial delay and underdamped dynamics

Active matter systems often are well approximated as overdamped, meaning that any inertial momentum is immediately dissipated by the environment. On the other hand, especially for macroscopic systems but also for many mesoscopic ones particle mass can become relevant for the dynamics. For such systems we recently proposed an underdamped continuum model which captures translationally inertial dynamics via two contributions. First, convection and second a damping time scale of inertial motion. In this paper, we ask how both of these features influence the collective behavior compared to overdamped dynamics by studying the example of the active phase field crystal model. We first focus on the case of suppressed convection to study the role of the damping time. We quantify that the relaxation process to the steady collective motion state is considerably prolonged with damping time due to the increasing occurrence of transient groups of circularly moving density peaks. Finally, we illustrate the fully underdamped case with convection. Instead of collective motion of density peaks we then find a coexistence of constant high and low density phases reminiscent of motility-induced phase separation.

Mean field approach of dynamical pattern formation in underdamped active matter with short-ranged alignment and distant anti-alignment interactions

Many active matter systems, especially on the microscopic scale, are well approximated as overdamped, meaning that any inertial momentum is immediately dissipated by the environment. On the other hand, especially for macroscopic active systems but also for many mesoscopic systems the time scale of translational inertial motion can become large enough to be relevant for the dynamics. This raises the question how collective dynamics and the resulting states in active matter are influenced by inertia. Therefore, we propose a coarse-grained continuum model for underdamped active matter based on a mean field description for passive systems. Furthermore, as an example, we apply the model to a system with interactions that support an alignment on short distances and an anti-alignment on longer length scales as known in the context of pattern formation due to orientational interactions. Our numerical calculations of the under- and overdamped dynamics both predict a structured laning state. However, activity induced convective flows that are only present in the underdamped model destabilize this state when the anti-alignment is weakened, leading to a collective motion state which does not occur in the overdamped limit. A turbulent transition regime between the two states can be characterized by strong density fluctuations and the absence of global ordering.

Clustering and phase separation in mixtures of dipolar and active particles

The self-assembly of colloidal particles in dynamic environments has become an important field of study because of potential applications in fabricating out-of-equilibrium materials. We investigate the phase behavior of mixtures of passive dipolar colloids and active soft spheres using Brownian dynamics simulations in two dimensions. The phase behaviors exhibited include dipolar percolated network, dipolar string-fluid, isotropic fluid, and a phase-separated state. We find that the clustering of dipolar colloids is enhanced in the presence of slow-moving active particles compared to the clustering of dipolar particles mixed with passive particles. When the active particle motility is high, the chains of dipolar particles are either broken into short chains or pushed into dense clusters. Motility-induced phase separation into dense and dilute phases is also present. The area fraction of particles in the dilute phase increases as the fraction of active particles in the system decreases, while the area fraction of particles in the dense phase remains constant. Our findings are relevant to the development of reconfigurable self-assembled materials.

Self-straining of actively crosslinked microtubule networks

Cytoskeletal networks are foundational examples of active matter and central to self-organized structures in the cell. In vivo, these networks are active and densely crosslinked. Relating their large-scale dynamics to the properties of their constituents remains an unsolved problem. Here, we study an in vitro active gel made from aligned microtubules and XCTK2 kinesin motors. Using photobleaching, we demonstrate that the gel's aligned microtubules, driven by motors, continually slide past each other at a speed independent of the local microtubule polarity and motor concentration. This phenomenon is also observed, and remains unexplained, in spindles. We derive a general framework for coarse graining microtubule gels crosslinked by molecular motors from microscopic considerations. Using microtubule-microtubule coupling through a force-velocity relationship for kinesin, this theory naturally explains the experimental results: motors generate an active strain rate in regions of changing polarity, which allows microtubules of opposite polarities to slide past each other without stressing the material.

Frequency-dependent higher-order Stokes singularities near a planar elastic boundary: Implications for the hydrodynamics of an active microswimmer near an elastic interface

The emerging field of self-driven active particles in fluid environments has recently created significant interest in the biophysics and bioengineering communities owing to their promising future for biomedical and technological applications. These microswimmers move autonomously through aqueous media, where under realistic situations they encounter a plethora of external stimuli and confining surfaces with peculiar elastic properties. Based on a far-field hydrodynamic model, we present an analytical theory to describe the physical interaction and hydrodynamic couplings between a self-propelled active microswimmer and an elastic interface that features resistance toward shear and bending. We model the active agent as a superposition of higher-order Stokes singularities and elucidate the associated translational and rotational velocities induced by the nearby elastic boundary. Our results show that the velocities can be decomposed in shear and bending related contributions which approach the velocities of active agents close to a no-slip rigid wall in the steady limit. The transient dynamics predict that contributions to the velocities of the microswimmer due to bending resistance are generally more pronounced than those due to shear resistance. Bending can enhance (suppress) the velocities resulting from higher-order singularities whereas the shear related contribution decreases (increases) the velocities. Most prominently, we find that near an elastic interface of only energetic resistance toward shear deformation, such as that of an elastic capsule designed for drug delivery, a swimming bacterium undergoes rotation of the same sense as observed near a no-slip wall. In contrast to that, near an interface of only energetic resistance toward bending, such as that of a fluid vesicle or liposome, we find a reversed sense of rotation. Our results provide insight into the control and guidance of artificial and synthetic self-propelling active microswimmers near elastic confinements.

Which interactions dominate in active colloids?

Despite mounting evidence that the same gradients, which active colloids use for swimming, induce important cross-interactions (phoretic interactions), they are still ignored in most many-body descriptions, perhaps to avoid complexity and a zoo of unknown parameters. Here we derive a simple model, which reduces phoretic far-field interactions to a pair-interaction whose strength is mainly controlled by one genuine parameter (swimming speed). The model suggests that phoretic interactions are generically important for autophoretic colloids (unless effective screening of the phoretic fields is strong) and should dominate over hydrodynamic interactions for the typical case of half-coating and moderately nonuniform surface mobilities. Unlike standard minimal models, but in accordance with canonical experiments, our model generically predicts dynamic clustering in active colloids at a low density. This suggests that dynamic clustering can emerge from the interplay of screened phoretic attractions and active diffusion.

Memory-based mediated interactions between rigid particulate inclusions in viscoelastic environments

Many practically relevant materials combine properties of viscous fluids and elastic solids to viscoelastic behavior. Our focus is on the induced dynamic behavior of damped finite-sized particulate inclusions in such substances. We explicitly describe history-dependent interactions that emerge between the embedded particles. These interactions are mediated by the viscoelastic surroundings. They result from the flows and distortions of the viscoelastic medium when induced by the rigid inclusions. Both viscoelastic environments of terminal fluidlike flow and of completely reversible damped elastic behavior are covered. For illustration and to highlight the role of the formalism in potential applications, we briefly address the relevant examples of dragging a rigid sphere through a viscoelastic environment together with subsequent relaxation dynamics, the switching dynamics of magnetic fillers in elastic gel matrices, and the swimming behavior of active microswimmers in viscoelastic solutions. The approach provides a basis for more quantitative and extended investigations of these and related systems in the future.

Delayed feedback control of active particles: a controlled journey towards the destination

We explore theoretically the navigation of an active particle based on delayed feedback control. The delayed feedback enters in our expression for the particle orientation which, for an active particle, determines (up to noise) the direction of motion in the next time step. Here we estimate the orientation by comparing the delayed position of the particle with the actual one. This method does not require any real-time monitoring of the particle orientation and may thus be relevant also for controlling sub-micron sized particles, where the imaging process is not easily feasible. We apply the delayed feedback strategy to two experimentally relevant situations, namely, optical trapping and photon nudging. To investigate the performance of our strategy, we calculate the mean arrival time analytically (exploiting a small-delay approximation) and by simulations.

A statistical physics view of swarming bacteria

Bacterial swarming is a collective mode of motion in which cells migrate rapidly over surfaces, forming dynamic patterns of whirls and jets. This review presents a physical point of view of swarming bacteria, with an emphasis on the statistical properties of the swarm dynamics as observed in experiments. The basic physical principles underlying the swarm and their relation to contemporary theories of collective motion and active matter are reviewed and discussed in the context of the biological properties of swarming cells. We suggest a paradigm according to which bacteria have optimized some of their physical properties as a strategy for rapid surface translocation. In other words, cells take advantage of favorable physics, enabling efficient expansion that enhances survival under harsh conditions.

Clustering of microswimmers: interplay of shape and hydrodynamics

The spatiotemporal dynamics in systems of active self-propelled particles is controlled by the propulsion mechanism in combination with various direct interactions, such as steric repulsion and hydrodynamics. These direct interactions are typically anisotropic, and come in different "flavors", such as spherical and elongated particle shapes, pusher and puller flow fields, etc. The combination of the various aspects is expected to lead to new emergent behavior. However, it is a priori not evident whether shape and hydrodynamics act synergistically or antagonistically to generate motility-induced clustering (MIC) and phase separation (MIPS). We employ a model of prolate spheroidal microswimmers-called squirmers-in quasi-two-dimensional confinement to address this issue by mesoscale hydrodynamic simulations. For comparison, non-hydrodynamic active Brownian particles (ABPs) are considered to elucidate the contribution of hydrodynamic interactions. For spherical particles, the comparison between ABPs and hydrodynamic-squirmer ensembles reveals a suppression of MIPS due to hydrodynamic interactions. Yet, our analysis shows that dynamic clusters exist, with a broad size distribution. The fundamental difference between ABPs and squirmers is attributed to an increased reorientation of squirmers by hydrodynamic torques during their collisions. In contrast, for elongated squirmers, hydrodynamics interactions enhance MIPS. The transition to a phase-separated state strongly depends on the nature of the swimmer's flow field-with an increased tendency toward MIPS for pullers, and a reduced tendency for pushers. Thus, hydrodynamic interactions show opposing effects on MIPS for spherical and elongated microswimmers, and details of the propulsion mechanism of biological microswimmers may be very important to determine their collective behavior.

Clustering and phase separation of circle swimmers dispersed in a monolayer

We perform Brownian dynamics simulations in two dimensions to study the collective behavior of circle swimmers, which are driven by both, an (effective) translational and rotational self-propulsion, and interact via steric repulsion. We find that active rotation generally opposes motility-induced clustering and phase separation, as demonstrated by a narrowing of the coexistence region upon increase of the propulsion angular velocity. Moreover, although the particles are intrinsically assigned to rotate counterclockwise, a novel state of clockwise vortices emerges at an optimal value of the effective propulsion torque. We propose a simple gear-like model to capture the underlying mechanism of the clockwise vortices.

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.

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.

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.

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.

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.