Particle-resolved lattice Boltzmann simulations of 3-dimensional active turbulence

Collective behavior of squirmers in thin films

Bacteria in biofilms form complex structures and can collectively migrate within mobile aggregates, which is referred to as swarming. This behavior is influenced by a combination of various factors, including morphological characteristics and propulsive forces of swimmers, their volume fraction within a confined environment, and hydrodynamic and steric interactions between them. In our study, we employ the squirmer model for microswimmers and the dissipative particle dynamics method for fluid modeling to investigate the collective motion of swimmers in thin films. The film thickness permits a free orientation of non-spherical squirmers, but constraints them to form a two-layered structure at maximum. Structural and dynamic properties of squirmer suspensions confined within the slit are analyzed for different volume fractions of swimmers, motility types (e.g., pusher, neutral squirmer, puller), and the presence of a rotlet dipolar flow field, which mimics the counter-rotating flow generated by flagellated bacteria. Different states are characterized, including a gas-like phase, swarming, and motility-induced phase separation, as a function of increasing volume fraction. Our study highlights the importance of an anisotropic swimmer shape, hydrodynamic interactions between squirmers, and their interaction with the walls for the emergence of different collective behaviors. Interestingly, the formation of collective structures may not be symmetric with respect to the two walls. Furthermore, the presence of a rotlet dipole significantly mitigates differences in the collective behavior between various swimmer types. These results contribute to a better understanding of the formation of bacterial biofilms and the emergence of collective states in confined active matter.

Harnessing density to control the duration of intermittent Lévy walks in bacterial turbulence

Dense bacterial suspensions display collective motion exhibiting coherent flow structures reminiscent of turbulent flows. However, in contrast to inertial turbulence, the microscopic dynamics underlying bacterial turbulence is only beginning to be understood. Here, we report experiments revealing correlations between microscopic dynamics and the emergence of collective motion in bacterial suspensions. Our results demonstrate the existence of three microscopic dynamical regimes: initial ballistic dynamics followed by an intermittent Lévy walk before the intriguing decay to random Gaussian fluctuations. Our experiments capture that the fluid correlation time earmarks the transition from Lévy to Gaussian fluctuations demonstrating the microscopic reason underlying the observation. By harnessing the flow activity via bacterial concentration, we reveal systematic control over the flow correlation timescales, which, in turn, allows controlling the duration of the Lévy walk.

Active turbulence and spontaneous phase separation in inhomogeneous extensile active gels

We report numerical results for the hydrodynamics of inhomogeneous lyotropic and extensile active nematic gels. By simulating the coupled Cahn-Hilliard, Navier-Stokes, and Beris-Edwards equation for the evolution of the composition, flow and orientational order of an active nematic, we ask whether composition variations are important to determine its emergent physics. As in active gels of uniform composition, we find that increasing either activity or nematic tendency (e.g., overall active matter concentration) triggers a transition between an isotropic passive phase and an active nematic one. We show that composition inhomogeneities are important in the latter phase, where we find three types of possible dynamical regimes. First, we observe regular patterns with defects and vortices: these exist close to the passive-active transition. Second, for larger activity, or deeper in the nematic phase, we find active turbulence, as in active gels of uniform composition, but with exceedingly large composition variation. In the third regime, which is uniquely associated with inhomogeneity and occurs for large nematic tendency and low activity, we observe spontaneous microphase separation into active and passive domains. The microphase separated regime is notable in view of the absence of an explicit demixing term in the underlying free energy which we use, and we provide a theoretical analysis based on the common tangent construction which explains its existence. We hope this regime can be probed experimentally in the future.

Alignment and propulsion of squirmer pusher-puller dumbbells

The properties of microswimmer dumbbells composed of pusher-puller pairs are investigated by mesoscale hydrodynamic simulations employing the multiparticle collision dynamics approach for the fluid. An individual microswimmer is represented by a squirmer, and various active-stress combinations in a dumbbell are considered. The squirmers are connected by a bond, which does not impose any geometrical restriction on the individual rotational motion. Our simulations reveal a strong influence of the squirmers' flow fields on the orientation of their propulsion directions, their fluctuations, and the swimming behavior of a dumbbell. The properties of pusher-puller pairs with equal magnitude of the active stresses dependent only weakly on the stress magnitude. This is similar to dumbbells of microswimmers without hydrodynamic interactions. However, for non-equal stress magnitudes, the active stress implies strong orientational correlations of the swimmers' propulsion directions with respect to each other as well as the bond vector. The orientational coupling is most pronounced for pairs with large differences of the active stress magnitude. The alignment of the squirmer propulsion directions with respect to each other is preferentially orthogonal in dumbbells with a strong pusher and weak puller, and antiparallel in the opposite case when the puller dominates. These strong correlations affect the active motion of dumbbells which is faster for strong pushers and slower for strong pullers.

Density fluctuations and energy spectra of 3D bacterial suspensions

Giant number fluctuations are often considered as a hallmark of the emergent nonequilibrium dynamics of active fluids. However, these anomalous density fluctuations have only been reported experimentally in two-dimensional dry active systems heretofore. Here, we investigate density fluctuations of bulk Escherichia coli suspensions, a paradigm of three-dimensional (3D) wet active fluids. Our experiments demonstrate the existence and quantify the scaling relation of giant number fluctuations in 3D bacterial suspensions. Surprisingly, the anomalous scaling persists at small scales in low-concentration suspensions well before the transition to active turbulence, reflecting the long-range nature of hydrodynamic interactions of 3D wet active fluids. To illustrate the origin of the density fluctuations, we measure the energy spectra of suspension flows and explore the density-energy coupling in both the steady and transient states of active turbulence. A scale-invariant density-independent correlation between density fluctuations and energy spectra is uncovered across a wide range of length scales. In addition, our experiments show that the energy spectra of bacterial turbulence exhibit the scaling of 3D active nematic fluids, challenging the common view of dense bacterial suspensions as active polar fluids.

Imaging the emergence of bacterial turbulence: Phase diagram and transition kinetics

We experimentally study the emergence of collective bacterial swimming, a phenomenon often referred to as bacterial turbulence. A phase diagram of the flow of 3D Escherichia coli suspensions spanned by bacterial concentration, the swimming speed of bacteria, and the number fraction of active swimmers is systematically mapped, which shows quantitative agreement with kinetic theories and demonstrates the dominant role of hydrodynamic interactions in bacterial collective swimming. We trigger bacterial turbulence by suddenly increasing the swimming speed of light-powered bacteria and image the transition to the turbulence in real time. Our experiments identify two unusual kinetic pathways, i.e., the one-step transition with long incubation periods near the phase boundary and the two-step transition driven by long-wavelength instabilities deep inside the turbulent phase. Our study provides not only a quantitative verification of existing theories but also insights into interparticle interactions and transition kinetics of bacterial turbulence.

Symmetric Mixtures of Pusher and Puller Microswimmers Behave as Noninteracting Suspensions

Suspensions of rear- and front-actuated microswimmers immersed in a fluid, known respectively as "pushers" and "pullers," display qualitatively different collective behaviors: beyond a characteristic density, pusher suspensions exhibit a hydrodynamic instability leading to collective motion known as active turbulence, a phenomenon which is absent for pullers. In this Letter, we describe the collective dynamics of a binary pusher-puller mixture using kinetic theory and large-scale particle-resolved simulations. We derive and verify an instability criterion, showing that the critical density for active turbulence moves to higher values as the fraction χ of pullers is increased and disappears for χ≥0.5. We then show analytically and numerically that the two-point hydrodynamic correlations of the 1∶1 mixture are equal to those of a suspension of noninteracting swimmers. Strikingly, our numerical analysis furthermore shows that the full probability distribution of the fluid velocity fluctuations collapses onto the one of a noninteracting system at the same density, where swimmer-swimmer correlations are strictly absent. Our results thus indicate that the fluid velocity fluctuations in 1∶1 pusher-puller mixtures are exactly equal to those of the corresponding noninteracting suspension at any density, a surprising cancellation with no counterpart in equilibrium long-range interacting systems.