Physical properties of collective motion in suspensions of bacteria

Vortex reversal is a precursor of confined bacterial turbulence

Active turbulence, or chaotic self-organized collective motion, is often observed in concentrated suspensions of motile bacteria and other systems of self-propelled interacting agents. To date, there is no fundamental understanding of how geometrical confinement orchestrates active turbulence and alters its physical properties. Here, by combining large-scale experiments, computer modeling, and analytical theory, we have identified a generic sequence of transitions occurring in bacterial suspensions confined in cylindrical wells of varying radii. With increasing the well's radius, we observed that persistent vortex motion gives way to periodic vortex reversals, four-vortex pulsations, and then well-developed active turbulence. Using computational modeling and analytical theory, we have shown that vortex reversal results from the nonlinear interaction of the first three azimuthal modes that become unstable with the radius increase. The analytical results account for our key experimental findings. To further validate our approach, we reconstructed equations of motion from experimental data. Our findings shed light on the universal properties of confined bacterial active matter and can be applied to various biological and synthetic active systems.

Collective motion of energy depot active disks

In the present work we have studied collectives of active disks with an energy depot, moving in the two-dimensional plane and interacting via an excluded volume. The energy depot accounts for the extraction of energy taking place at the level of each particle in order to perform self-propulsion, included in an underdamped Langevin dynamics. We show that this model undergoes a flocking transition, exhibiting some of the key features of the Vicsek model, namely, band formation and giant number fluctuations. These bands, either single or multiple, are dense and very strongly polarised propagating structures. Large density bands disappear as the activity is further increased, eventually reaching a homogeneous polar state. We unravel an effective alignment interaction at the level of two-particle collisions that can be controlled by activity and gives rise to flocking at large scales.

Dynamically induced spatial segregation in multispecies bacterial bioconvection

Active matter, from motile bacteria to animals, can exhibit striking collective and coherent behavior. Despite significant advances in understanding the behavior of homogeneous systems, little is known about the self-organization and dynamics of heterogeneous active matter, such as complex and diverse bacterial communities. Under oxygen gradients, many bacterial species swim towards air-liquid interfaces in auto-organized, directional bioconvective flows, whose spatial scales exceed the cell size by orders of magnitude. Here we show that multispecies bacterial suspensions undergoing oxytactic-driven bioconvection exhibit dynamically driven spatial segregation, despite the enhanced mixing of bioconvective flows, and the fact that these species coexist in their natural habitat. Segregation is observed as patterns of spatially interlocked domains, with local dominance of one of the constituent species in the suspension. Our findings suggest that segregation mechanisms are driven by species-specific motile behaviors under conditions of hydrodynamic flow, rather than biochemical repulsion. Thus, species with different motile characteristics in the same ecological context can enhance their access to limiting resources. This work provides novel insights on the role of heterogeneity in active matter, as well as on the dynamics of complex microbial communities, their spatial organization and their collective behavior.

Deviation of a system of nonreciprocally coupled harmonic oscillators from a conservative system

Discrete systems of coupled linear mechanical oscillators with nonreciprocal interaction are a model for a variety of physical systems. In general, the presence of nonreciprocal interactions renders their dynamics nonconservative, but under certain conditions it remains conservative. In this paper we show which thermodynamic properties induced by nonreciprocity can be observed in conservative systems and which are specific to nonconservative systems. To this end, we formulate a criterion for identifying conservative systems and construct a measure to quantify the deviation from conservativity.

Long-lived unidirectional flow of active particles within long narrow channels

Revealing the mechanism of directed transport of active matter is critical for advancing our fundamental understanding of non-equilibrium physics. Asymmetric microstructures are commonly used to rectify random movement of active particles. However, it remains unclear as to how to achieve unidirectional movement of active particles in long narrow channels. Here, we study the dynamics of active particles in a device which is divided into two chambers by V-shaped barriers and connected by a narrow channel. We find three distinct movement modes of active particles within this symmetric channel, including stochastic movement, self-sustained oscillation, and long-lived unidirectional flows. We demonstrate that the three movement modes are determined by the competition between the ratchet effect induced by the V-shaped barriers and the particle transport mediated by the long-narrow channel. Finally, we show that the unidirectional particle flow can serve as an "energy battery" to continuously supply energy for the directed transport of other objects. Our findings offer valuable insights into a unique approach for realizing unidirectional movement of active matter and open new avenues for application in microfluidics and material transport.

Strong Casimir-like Forces in Flocking Active Matter

Confining in space the equilibrium fluctuations of statistical systems with long-range correlations is known to result into effective forces on the boundaries. Here we demonstrate the occurrence of Casimir-like forces in the nonequilibrium context provided by flocking active matter. In particular, we consider a system of aligning self-propelled particles in two spatial dimensions that are transversally confined by reflecting or partially reflecting walls. We show that in the ordered flocking phase this confined active vectorial fluid is characterized by extensive boundary layers, as opposed to the finite ones usually observed in confined scalar active matter. Moreover, a finite-size, fluctuation-induced contribution to the pressure on the wall emerges, which decays slowly and algebraically upon increasing the distance between the walls. We explain our findings-which display a certain degree of universality-within a hydrodynamic description of the density and velocity fields.

Scaling Transition of Active Turbulence from Two to Three Dimensions

Turbulent flows are observed in low-Reynolds active fluids, which display similar phenomenology to the classical inertial turbulence but are of a different nature. Understanding the dependence of this new type of turbulence on dimensionality is a fundamental challenge in non-equilibrium physics. Real-space structures and kinetic energy spectra of bacterial turbulence are experimentally measured from two to three dimensions. The turbulence shows three regimes separated by two critical confinement heights, resulting from the competition of bacterial length, vortex size and confinement height. Meanwhile, the kinetic energy spectra display distinct universal scaling laws in quasi-2D and 3D regimes, independent of bacterial activity, length, and confinement height, whereas scaling exponents transition in two steps around the critical heights. The scaling behaviors are well captured by the hydrodynamic model we develop, which employs image systems to represent the effects of confining boundaries. The study suggests a framework for investigating the effect of dimensionality on non-equilibrium self-organized systems.

Exploring sperm cell rheotaxis in microfluidic channel: the role of flow and viscosity

Rheotaxis is a fundamental mechanism of sperm cells that guides them in navigating towards the oocyte. The present study investigates the phenomenon of sperm rheotaxis in Newtonian and non-Newtonian fluid media, which for the first time explores a viscosity range equivalent to that of the oviductal fluid of the female reproductive tract in rectilinear microfluidic channels. Three parameters, the progressive velocity while performing rheotaxis, the radius of rotation during rheotaxis, and the percentage of rheotactic sperm cells in the bulk and near-wall regions of the microfluidic channel were measured. Numerical simulations of the flow were conducted to estimate the shear rate, flow velocity, and the drag force acting on the sperm head at specific locations where the sperms undergo rheotaxis. Increasing the flow velocity resulted in a change in the position of rheotactic sperm from the bulk center to the near wall region, an increase and subsequent decrease in the sperm's upstream progressive velocity, and a decrease in the radius of rotation. We observed that with an increase in viscosity, rheotactic sperms migrate to the near wall regions at lower flow rates, the upstream progressive velocity of the sperm decreases for Newtonian and increases for non-Newtonian media, and the radius of rotation increases for Newtonian and decreases for non-Newtonian media. These results quantify the effects of fluid properties such as viscosity and flow rate on sperm rheotaxis and navigation, thereby paving the way for manipulating sperm behavior in microfluidic devices, potentially leading to advancements in assisted reproduction techniques.

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.

Nonequilibrium hyperuniform states in active turbulence

We demonstrate that the complex spatiotemporal structure in active fluids can feature characteristics of hyperuniformity. Using a hydrodynamic model, we show that the transition from hyperuniformity to nonhyperuniformity and antihyperuniformity depends on the strength of active forcing and can be related to features of active turbulence without and with scaling characteristics of inertial turbulence. Combined with identified signatures of Levy walks and nonuniversal diffusion in these systems, this allows for a biological interpretation and the speculation of nonequilibrium hyperuniform states in active fluids as optimal states with respect to robustness and strategies of evasion and foraging.

Active Darcy's Law

While bacterial swarms can exhibit active turbulence in vacant spaces, they naturally inhabit crowded environments. We numerically show that driving disorderly active fluids through porous media enhances Darcy's law. While purely active flows average to zero flux, hybrid active/driven flows display greater drift than purely pressure-driven flows. This enhancement is nonmonotonic with activity, leading to an optimal activity to maximize flow rate. We incorporate the active contribution into an active Darcy's law, which may serve to help understand anomalous transport of swarming in porous media.

Anomalous dynamics of a passive droplet in active turbulence

Motion of a passive deformable object in an active environment serves as a representative of both in-vivo systems such as intracellular particle motion in Acanthamoeba castellanii, or in-vitro systems such as suspension of beads inside dense swarms of Escherichia coli. Theoretical modeling of such systems is challenging due to the requirement of well resolved hydrodynamics which can explore the spatiotemporal correlations around the suspended passive object in the active fluid. We address this critical lack of understanding using coupled hydrodynamic equations for nematic liquid crystals with finite active stress to model the active bath, and a suspended nematic droplet with zero activity. The droplet undergoes deformation fluctuations and its movement shows periods of "runs" and "stays". At relatively low interfacial tension, the droplet begins to break and mix with the outer active bath. We establish that the motion of the droplet is influenced by the interplay of spatial correlations of the flow and the size of the droplet. The mean square displacement shows a transition from ballistic to normal diffusion which depends on the droplet size. We discuss this transition in relation to spatiotemporal scales associated with velocity correlations of the active bath and the droplet.

United under stress: High-speed transport network emerging at bacterial living edge

Individuals tend to move freely when there is enough room but would act collectively for their survival under external stress. In the case of living cells, for instance, when a drop of low-density flagellated bacterial solution is transferred onto the agar surface, the initially disordered movement of individual bacteria would be replaced with coordinated cell swarming after a lag phase of a few hours. Here, we study how such cooperation is established while overcoming the disorder at the onset of the lag phase with single nanoparticle tracking. Upon the spreading of the droplet, the bacteria in the solution cluster and align near the almost immobilized contact line confining the drop, forming a narrow ring of cells. As individual cells move in and out of the ring continuously, certain flow patterns emerge in the inter-bacterial fluid. We reveal high-speed long-distance unidirectional flows with definite chirality along the outside of the ring, along the inside of the ring and across the ring. We speculate that these flows enable the fast and efficient transport, facilitating the communication and unification of the bacterial community.

Pair Interactions of Self-Propelled SiO2-Pt Janus Colloids: Chemically Mediated Encounters

Driven by the necessity to achieve a thorough comprehension of the bottom-up fabrication process of functional materials, this experimental study investigates the pairwise interactions or collisions between chemically active SiO2-Pt Janus colloids. These collisions are categorized based on the Janus colloids' orientations before and after they make physical contact. In addition to the hydrodynamic interactions, the Janus colloids are also known to affect each other's chemical field, resulting in chemophoretic interactions, which depend on the degree of surface anisotropy in reactivity of Janus colloid and the solute-surface interaction at play. Our study reveals that these interactions lead to a noticeable decrease in particle speed and changes in orientation that correlate with the contact duration and yield different collision types. Distinct configurations of contact during collisions were found, whose mechanisms and likelihood are found to be dependent primarily on the chemical interactions. Such estimates of collision and their characterization in dilute suspensions shall have a key impact in determining the arrangement and time scales of dynamical structures and assemblies of denser suspensions and potentially the functional materials of the future.

Vortex Pattern Stabilization in Thin Films Resulting from Shear Thickening of Active Suspensions

The need for structuring on micrometer scales is abundant, for example, in view of phononic applications. We here outline a novel approach based on the phenomenon of active turbulence on the mesoscale. As we demonstrate, a shear-thickening carrier fluid of active microswimmers intrinsically stabilizes regular vortex patterns of otherwise turbulent active suspensions. The fluid self-organizes into a periodically structured nonequilibrium state. Introducing additional passive particles of intermediate size leads to regular spatial organization of these objects. Our approach opens a new path toward functionalization through patterning of thin films and membranes.

Fluid handling by foam wound dressings: From engineering theory to advanced laboratory performance evaluations

This article describes the contemporary bioengineering theory and practice of evaluating the fluid handling performance of foam-based dressings, with focus on the important and clinically relevant engineering structure-function relationships and on advanced laboratory testing methods for pre-clinical quantitative assessments of this common type of wound dressings. The effects of key wound dressing material-related and treatment-related physical factors on the absorbency and overall fluid handling of foam-based dressings are thoroughly and quantitively analysed. Discussions include exudate viscosity and temperature, action of mechanical forces and the dressing microstructure and associated interactions. Based on this comprehensive review, we propose a newly developed testing method, experimental metrics and clinical benchmarks that are clinically relevant and can set the standard for robust fluid handling performance evaluations. The purpose of this evaluative framework is to translate the physical characteristics and performance determinants of a foam dressing into achievable best clinical outcomes. These guiding principles are key to distinguishing desirable properties of a dressing that contribute to optimal performance in clinical settings.

Characterizing different motility-induced regimes in active matter with machine learning and noise

We examine motility-induced phase separation (MIPS) in two-dimensional run-and-tumble disk systems using both machine learning and noise fluctuation analysis. Our measures suggest that within the MIPS state there are several distinct regimes as a function of density and run time, so that systems with MIPS transitions exhibit an active fluid, an active crystal, and a critical regime. The different regimes can be detected by combining an order parameter extracted from principal component analysis with a cluster stability measurement. The principal component-derived order parameter is maximized in the critical regime, remains low in the active fluid, and has an intermediate value in the active crystal regime. We demonstrate that machine learning can better capture dynamical properties of the MIPS regimes compared to more standard structural measures such as the maximum cluster size. The different regimes can also be characterized via changes in the noise power of the fluctuations in the average speed. In the critical regime, the noise power passes through a maximum and has a broad spectrum with a 1/f^{1.6} signature, similar to the noise observed near depinning transitions or for solids undergoing plastic deformation.

Spontaneous shock waves in pulse-stimulated flocks of Quincke rollers

Active matter demonstrates complex spatiotemporal self-organization not accessible at equilibrium and the emergence of collective behavior. Fluids comprised of microscopic Quincke rollers represent a popular realization of synthetic active matter. Temporal activity modulations, realized by modulated external electric fields, represent an effective tool to expand the variety of accessible dynamic states in active ensembles. Here, we report on the emergence of shockwave patterns composed of coherently moving particles energized by a pulsed electric field. The shockwaves emerge spontaneously and move faster than the average particle speed. Combining experiments, theory, and simulations, we demonstrate that the shockwaves originate from intermittent spontaneous vortex cores due to a vortex meandering instability. They occur when the rollers' translational and rotational decoherence times, regulated by the electric pulse durations, become comparable. The phenomenon does not rely on the presence of confinement, and multiple shock waves continuously arise and vanish in the system.

Viscoelasticity enhances collective motion of bacteria

Bacteria form human and animal microbiota. They are the leading causes of many infections and constitute an important class of active matter. Concentrated bacterial suspensions exhibit large-scale turbulent-like locomotion and swarming. While the collective behavior of bacteria in Newtonian fluids is relatively well understood, many fundamental questions remain open for complex fluids. Here, we report on the collective bacterial motion in a representative biological non-Newtonian viscoelastic environment exemplified by mucus. Experiments are performed with synthetic porcine gastric mucus, natural cow cervical mucus, and a Newtonian-like polymer solution. We have found that an increase in mucin concentration and, correspondingly, an increase in the suspension's elasticity monotonously increases the length scale of collective bacterial locomotion. On the contrary, this length remains practically unchanged in Newtonian polymer solution in a wide range of concentrations. The experimental observations are supported by computational modeling. Our results provide insight into how viscoelasticity affects the spatiotemporal organization of bacterial active matter. They also expand our understanding of bacterial colonization of mucosal surfaces and the onset of antibiotic resistance due to swarming.

Scale-free correlations and potential criticality in weakly ordered populations of brain cancer cells

Collective behavior spans several orders of magnitude of biological organization, from cell colonies to flocks of birds. We used time-resolved tracking of individual glioblastoma cells to investigate collective motion in an ex vivo model of glioblastoma. At the population level, glioblastoma cells display weakly polarized motion in the (directional) velocities of single cells. Unexpectedly, fluctuations in velocities are correlated over distances many times the size of a cell. Correlation lengths scale linearly with the maximum end-to-end length of the population, indicating that they are scale-free and lack a characteristic decay scale other than the size of the system. Last, a data-driven maximum entropy model captures statistical features of the experimental data with only two free parameters: the effective length scale (nc) and strength (J) of local pairwise interactions between tumor cells. These results show that glioblastoma assemblies exhibit scale-free correlations in the absence of polarization, suggesting that they may be poised near a critical point.

Spontaneous symmetry-breaking of the active cluster drives the directed movement and self-sustained oscillation of symmetric rod-like passive particles

Active particles without detailed balance can rectify their random motions to drive the directed movement or rotation of asymmetric passive obstacles. However, whether they can drive the directed movement of symmetric passive obstacles is still unclear. Here, we show that a rod-like passive particle which is fixed to move along the x-axis in an active bath can keep long-lived directed movement at nearly constant speed due to the spontaneous symmetry breaking of the neighboring active particle cluster. If the passive particle is further confined by a harmonic potential, it may undergo self-sustained periodic oscillation for an appropriate length of the passive particle and self-propelled velocity of active particles. The restoring force from the harmonic potential will trigger the velocity jump-off and thus lead to self-sustained periodic oscillation. Remarkably, the relationship between the velocity of the passive particle and the external force shows that the effective viscosity of the active bath may become negative in some regime. Finally, we develop a minimum 1D theoretical model to further probe the mechanism underlying the directed movement and self-sustained oscillation of the passive particle. Our findings reveal the effect of the moving boundary on the active bath and demonstrate a novel method to extract practical mechanical work from the active bath to propel microdevices.

Phase behavior of active and passive dumbbells

We report phase separation in a mixture of "hot" and "cold" three-dimensional dumbbells which interact by Lennard-Jones potential. We also have studied the effect of asymmetry of dumbbells and the variation of ratio of "hot" and "cold" dumbbells on their phase separation. The ratio of the temperature difference between hot and cold dumbbells to the temperature of cold dumbbells is a measure of the activity χ of the system. From constant density simulations of symmetric dumbbells, we observe that the "hot" and "cold" dumbbells phase separate at higher activity ratio (χ>5.80) compared to that of a mixture of hot and cold Lennard-Jones monomers (χ>3.44). We find that, in the phase-separated system, the hot dumbbells have high effective volume and hence high entropy which is calculated by two-phase thermodynamic method. The high kinetic pressure of hot dumbbells forces the cold dumbbells to form dense clusters such that at the interface the high kinetic pressure of hot dumbbells is balanced by the virial pressure of cold dumbbells. We find that phase separation pushes the cluster of cold dumbbells to have solidlike ordering. Bond orientation order parameters reveal that the cold dumbbells form solidlike ordering consisting of predominantly face-centered cubic and hexagonal-close packing packing, but the individual dumbbells have random orientations. The simulation of the nonequilibrium system of symmetric dumbbells at different ratios of number of hot dumbbells to cold dumbbells reveals that the critical activity of phase separation decreases with increase in fraction of hot dumbbells. The simulation of equal mixture of hot and cold asymmetric dumbbells revealed that the critical activity of phase separation was independent of the asymmetry of dumbbells. We also observed that the clusters of cold asymmetric dumbbells showed both crystalline and noncrystalline order depending on the asymmetry of dumbbells.

Alignment rule and geometric confinement lead to stability of a vortex in active flow

Vortices are hallmarks of a wide range of nonequilibrium phenomena in fluids at multiple length scales. In this work, we numerically study the whirling motion of self-propelled soft point particles confined in circular domain, and aim at addressing the stability issue of the coherent vortex structure. By the combination of dynamical and statistical analysis at the individual particle level, we reveal the persistence of the whirling motion resulting from the subtle competition of activity and geometric confinement. In the stable whirling motion, the scenario of the coexistence of the irregular microscopic motions of individual particles and the regular global whirling motion is fundamentally different from the motion of a vortex in passive fluid. Possible orientational order coexisting with the whirling are further explored. This work shows the stability mechanism of vortical dynamics in active media under the alignment rule in confined space and may have implications in creating and harnessing macroscale coherent dynamical states by tuning the confining geometry.

Hydrodynamic Synchronization of Chiral Microswimmers

We study synchronization in bulk suspensions of spherical microswimmers with chiral trajectories using large scale numerics. The model is generic. It corresponds to the lowest order solution of a general model for self-propulsion at low Reynolds numbers, consisting of a nonaxisymmetric rotating source dipole. We show that both purely circular and helical swimmers can spontaneously synchronize their rotation. The synchronized state corresponds to velocity alignment with high orientational order in both the polar and azimuthal directions. Finally, we consider a racemic mixture of helical swimmers where intraspecies synchronization is observed while the system remains as a spatially uniform fluid. Our results demonstrate hydrodynamic synchronization as a natural collective phenomenon for microswimmers with chiral trajectories.

Flocking in one dimension: Asters and reversals

We study the one-dimensional active Ising model in which aligning particles undergo diffusion biased by the signs of their spins. The phase diagram obtained varying the density of particles, their hopping rate, and the temperature controlling the alignment shows a homogeneous disordered phase but no homogeneous ordered one, as well as two phases with localized dense structures. In the flocking phase, large ordered aggregates move ballistically and stochastically reverse their direction of motion. In what we termed the "aster" phase, dense immobile aggregates of opposite magnetization face each other, exchanging particles, without any net motion of the aggregates. Using a combination of numerical simulations and mean-field theory, we study the evolution of the shapes of the flocks, the statistics of their reversal times, and their coarsening dynamics. Solving exactly for the zero-temperature dynamics of an aster allows us to understand their coarsening, which shows extremal dynamics, while mean-field equations account for their shape.

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.

Asymptotics and Summation of the Effective Properties of Suspensions, Simple Liquids and Composites

We review the problem of summation for a very short truncation of a power series by means of special resummation techniques inspired by the field-theoretical renormalization group. Effective viscosity (EV) of active and passive suspensions is studied by means of a special algebraic renormalization approach applied to the first and second-order expansions in volume fractions of particles. EV of the 2D and 3D passive suspensions is analysed by means of various self-similar approximants such as iterated roots, exponential approximants, super-exponential approximants and root approximants. General formulae for all concentrations are derived. A brief introduction to the rheology of micro-swimmers is given. Microscopic expressions for the intrinsic viscosity of the active system of puller-like microswimmers are obtained. Special attention is given to the problem of the calculation of the critical indices and amplitudes of the EV and to the sedimentation rate in the vicinity of known critical points. Critical indices are calculated from the short truncation by means of minimal difference and minimal derivative conditions on the fixed points imposed directly on the critical properties. Accurate expressions are presented for the non-local diffusion coefficient of a simple liquid in the vicinity of a critical point. Extensions and corrections to the celebrated Kawasaki formula are discussed. We also discuss the effective conductivity for the classical analog of graphene and calculate the effective critical index for superconductivity dependent on the concentration of vacancies. Finally, we discuss the effective conductivity of a random 3D composite and calculate the superconductivity critical index of a random 3D composite.

Self-sustained non-equilibrium co-existence of fluid and solid states in a strongly coupled complex plasma system

A complex (dusty) plasma system is well known as a paradigmatic model for studying the kinetics of solid-liquid phase transitions in inactive condensed matter. At the same time, under certain conditions a complex plasma system can also display characteristics of an active medium with the micron-sized particles converting energy of the ambient environment into motility and thereby becoming active. We present a detailed analysis of the experimental complex plasmas system that shows evidence of a non-equilibrium stationary coexistence between a cold crystalline and a hot fluid state in the structure due to the conversion of plasma energy into the motion energy of microparticles in the central region of the system. The plasma mediated non-reciprocal interaction between the dust particles is the underlying mechanism for the enormous heating of the central subsystem, and it acts as a micro-scale energy source that keeps the central subsystem in the molten state. Accurate multiscale simulations of the system based on combined molecular dynamics and particle-in-cell approaches show that strong structural nonuniformity of the system under the action of electostatic trap makes development of instabilities a local process. We present both experimental tests conducted with a complex plasmas system in a DC glow discharge plasma and a detailed theoretical analysis.

Bistability in orbital trajectories of a chiral self-propelled particle interacting with an external field

In this paper, the dynamics of a self-propelled stochastic particle under the influence of an axisymmetric light field is experimentally studied. The particle under consideration has the main characteristic of carrying a light sensor in an eccentric location. For the chosen experimental conditions, the emerging trajectories are orbital, and, more interestingly, they suggest the existence of bistability. A mathematical model incorporating the key experimental components is introduced. By means of numerical simulations and theoretical analysis, it is found that, in addition to the orbiting behavior, the sensor location could produce trapped or diffusive behaviors. Furthermore, the study reveals that stochastic perturbation and the eccentric location of the sensor are responsible for inducing bistability in the orbital trajectories, supporting experimental observations.

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.

Hyperuniform Active Chiral Fluids with Tunable Internal Structure

Large density fluctuations observed in active systems and hyperuniformity are two seemingly incompatible phenomena. However, the formation of hyperuniform states has been recently predicted in nonequilibrium fluids formed by chiral particles performing circular motion with the same handedness. Here we report evidence of hyperuniformity realized in a chiral active fluid comprised of pear-shaped Quincke rollers of arbitrary handedness. We show that hyperuniformity and large density fluctuations, triggered by dynamic clustering, coexist in this system at different length scales. The system loses its hyperuniformity as the curvature of particles' motion increases, transforming them into localized spinners. Our results experimentally demonstrate a novel hyperuniform active fluid and provide new insights into an interplay between chirality, activity, and hyperuniformity.

Influence of confinement on the spreading of bacterial populations

The spreading of bacterial populations is central to processes in agriculture, the environment, and medicine. However, existing models of spreading typically focus on cells in unconfined settings—despite the fact that many bacteria inhabit complex and crowded environments, such as soils, sediments, and biological tissues/gels, in which solid obstacles confine the cells and thereby strongly regulate population spreading. Here, we develop an extended version of the classic Keller-Segel model of bacterial spreading via motility that also incorporates cellular growth and division, and explicitly considers the influence of confinement in promoting both cell-solid and cell-cell collisions. Numerical simulations of this extended model demonstrate how confinement fundamentally alters the dynamics and morphology of spreading bacterial populations, in good agreement with recent experimental results. In particular, with increasing confinement, we find that cell-cell collisions increasingly hinder the initial formation and the long-time propagation speed of chemotactic pulses. Moreover, also with increasing confinement, we find that cellular growth and division plays an increasingly dominant role in driving population spreading—eventually leading to a transition from chemotactic spreading to growth-driven spreading via a slower, jammed front. This work thus provides a theoretical foundation for further investigations of the influence of confinement on bacterial spreading. More broadly, these results help to provide a framework to predict and control the dynamics of bacterial populations in complex and crowded environments.

The spreading of bacteria through their environments critically impacts our everyday lives; it can be harmful, underlying the progression of infections and spoilage of foods, or can be beneficial, enabling the delivery of therapeutics, sustaining plant growth, and remediating polluted terrain. In all these cases, bacteria typically inhabit crowded environments, such as soils, sediments, and biological tissues/gels, in which solid obstacles confine the cells and regulate their spreading. However, existing models of spreading typically focus on cells in unconfined settings, and thus are frequently not applicable to cells in more complex environments. Here, we address this gap in knowledge by extending the classic Keller-Segel model of bacterial spreading via motility to also incorporate cellular growth and division, and explicitly consider the influence of confinement. Through numerical simulations of this extended model, we show how confinement fundamentally alters the dynamics and morphology of spreading bacterial populations—in particular, driving a transition from chemotactic spreading of motile cells to growth-driven spreading via a slower, jammed front. These results provide a foundation for further investigations of the influence of confinement on bacterial spreading, both by yielding testable predictions for future experiments, and by providing guidelines to predict and control the dynamics of bacterial populations in complex and crowded environments.

Activity-Induced Enhancement of Superdiffusive Transport in Bacterial Turbulence

Superdiffusion processes significantly promote the transport of tiny passive particles within biological fluids. Activity, one of the essential measures for living matter, however, is less examined in terms of how and to what extent it can improve the diffusivity of the moving particles. Here, bacterial suspensions are confined within the microfluidic channel at the state of bacterial turbulence, and are tuned to different activity levels by oxygen consumption in control. Systematic measurements are conducted to determine the superdiffusion exponent, which characterizes the diffusivity strength of tracer particles, depending on the continuously injecting energy converted to motile activity from swimming individuals. Higher activity is quantified to drastically enhance the superdiffusion process of passive tracers in the short-time regime. Moreover, the number density of the swimming bacteria is controlled to contribute to the field activity, and then to strengthen the super-diffusivity of tracers, distinguished by regimes with and without collective motion of interacting bacteria. Finally, the non-slip surfaces of the microfluidic channel lower the superdiffusion of immersed tracers due to the resistance, with the small diffusivity differing from the counterpart in the bulk. The findings here suggest ways of controlled diffusion and transport of substances within the living system with different levels of nutrition and resources and boundary walls, leading to efficient mixing, drug delivery and intracellular communications.

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.

Guiding Self-Assembly of Active Colloids by Temporal Modulation of Activity

Self-organization phenomena in ensembles of self-propelled particles open pathways to the synthesis of new dynamic states not accessible by traditional equilibrium processes. The challenge is to develop a set of principles that facilitate the control and manipulation of emergent active states. Here, we report that dielectric rolling colloids energized by a pulsating electric field self-organize into alternating square lattices with a lattice constant controlled by the parameters of the field. We combine experiments and simulations to examine spatiotemporal properties of the emergent collective patterns and investigate the underlying dynamics of the self-organization.We reveal the resistance of the dynamic lattices to compression and expansion stresses leading to a hysteretic behavior of the lattice constant. The general mechanism of pattern synthesis and control in active ensembles via temporal modulation of activity can be applied to other active colloidal systems.

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.

Flow transitions and length scales of a channel-confined active nematic

We perform lattice Boltzmann simulations of an active nematic fluid confined in a two-dimensional channel to study the range of flow states that are stabilised by the confinement: unidirectional flow, oscillatory flow, the dancing state, localised active turbulence and fully-developed active turbulence. We analyse the flows in Fourier space, and measure a range of different length scales which describe the flows. We argue that the different states occur as a result of flow instabilities inherent to the system. As a consequence the characteristic length scale for oscillatory flow, the dancing state and localised active turbulence is set by the channel width. Fully-developed active turbulence occurs only when the channel width is larger than the intrinsic, active length scale of the bulk fluid. The results clarify why the activity number is a control parameter for the flow transitions.

Surveying a Swarm: Experimental Techniques to Establish and Examine Bacterial Collective Motion

The survival and successful spread of many bacterial species hinges on their mode of motility. One of the most distinct of these is swarming, a collective form of motility where a dense consortium of bacteria employ flagella to propel themselves across a solid surface. Surface environments pose unique challenges, derived from higher surface friction/tension and insufficient hydration. Bacteria have adapted by deploying an array of mechanisms to overcome these challenges. Beyond allowing bacteria to colonize new terrain in the absence of bulk liquid, swarming also bestows faster speeds and enhanced antibiotic resistance to the collective. These crucial attributes contribute to the dissemination, and in some cases pathogenicity, of an array of bacteria. This mini-review highlights; 1) aspects of swarming motility that differentiates it from other methods of bacterial locomotion. 2) Facilitatory mechanisms deployed by diverse bacteria to overcome different surface challenges. 3) The (often difficult) approaches required to cultivate genuine swarmers. 4) The methods available to observe and assess the various facets of this collective motion, as well as the features exhibited by the population as a whole.

Emergence and dynamics of unconfined self-organised vortices in active magnetic roller liquids

Actively driven colloids demonstrate complex out-of-equilibrium dynamics often rivaling self-organized patterns and collective behavior observed in living systems. Recent studies revealed the emergence of steady macroscopic states with multiple interacting vortices in an unconfined environment that emerge from the coupling between microscale particle rotation and translation. Yet, insights into the microscopic behavior during the vortex emergence, growth, and formation of a multi-vortical state remain lacking. Here, we investigate in experiments and simulations how the microscale magnetic roller behavior leads to the emergence of seed vortices, their aggregation or annihilation, and the formation of stable large-scale vortical structures. We reveal that the coupling of roller-induced hydrodynamic flows guides the local self-densifications and self-organization of the micro-rollers into seed vortices. The resulting multi-vortical state is sensitive to the external magnetic field amplitude and allows tuning the rollers' number density in a vortex and its characteristic size.

Emergent probability fluxes in confined microbial navigation

When the motion of a motile cell is observed closely, it appears erratic, and yet the combination of nonequilibrium forces and surfaces can produce striking examples of organization in microbial systems. While most of our current understanding is based on bulk systems or idealized geometries, it remains elusive how and at which length scale self-organization emerges in complex geometries. Here, using experiments and analytical and numerical calculations, we study the motion of motile cells under controlled microfluidic conditions and demonstrate that probability flux loops organize active motion, even at the level of a single cell exploring an isolated compartment of nontrivial geometry. By accounting for the interplay of activity and interfacial forces, we find that the boundary's curvature determines the nonequilibrium probability fluxes of the motion. We theoretically predict a universal relation between fluxes and global geometric properties that is directly confirmed by experiments. Our findings open the possibility to decipher the most probable trajectories of motile cells and may enable the design of geometries guiding their time-averaged motion.

Active bath-induced localization and collapse of passive semiflexible polymers

The conformational and dynamical properties of a passive polymer embedded in a bath of active Brownian particles (ABPs) are studied by Langevin dynamics simulations. Various activities and ABP concentrations below and above the critical values for motility-induced phase separation (MIPS) are considered. In a homogeneous ABP fluid, the embedded polymer swells with increasing bath activity, with stronger swelling for larger densities. The polymer dynamics is enhanced, with the diffusion coefficient increasing by a power-law with increasing activity, where the exponent depends on the ABP concentration. For ABP concentrations in the MIPS regime, we observe a localization of the polymer in the low-density ABP phase associated with polymer collapse for moderate activities and a reswelling for high activities accompanied by a preferred localization in the high-density ABP phase. Localization and reswelling are independent of the polymer stiffness, with stiff polymers behaving similarly to flexible polymers. The polymer collapse is associated with a slowdown of its dynamics and a significantly smaller center-of-mass diffusion coefficient. In general, the polymer dynamics can only partially be described by an effective (bath) temperature. Moreover, the properties of a polymer embedded in a homogeneous active bath deviate quantitatively from those of a polymer composed of active monomers, i.e., linear chains of ABPs; however, such a polymer exhibits qualitatively similar activity-dependent features.

Rheology of bacterial superfluids in viscous environments

Viscous environments are ubiquitous in nature and in engineering applications, from mucus in lungs to oil recovery strategies in the earth's subsurface - and in all these environments, bacteria also thrive. The behavior of bacteria in viscous environments has been investigated for a single bacterium, but not for active suspensions. Dense populations of pusher-type bacteria are known to create superfluidic regimes where the effective viscosity of the entire suspension is reduced through collective motion, and the main purpose of this study is to investigate how a viscous environment will affect this behavior. Using a Couette rheometer, we measure shear stress as a function of the applied shear rate to define the effective viscosity of suspensions of Escherichia coli (E. coli), while varying both the bacterial density within the suspension and the viscosity of the suspending fluid. We document the remarkable observation that E. coli decreases the effective suspension viscosity to near-zero (superfluidic regime) for all solvent viscosities tested (1-17 mPa s). Specifically, we observe that the bacterial density needed to trigger this superfluidic regime and the maximum shear rate under which this regime can be sustained both decrease with increasing solvent viscosity. We find that the resulting rheograms can be well approximated by the Carreau-Yasuda law. Using this, we propose a constitutive model as a function of the solvent viscosity and the bacterial concentration only. This model captures the onset of the superfluidic regime and offers promising avenues for the modelling of flow of bacterial suspensions in viscous environments.

Odd Viscosity in Active Matter: Microscopic Origin and 3D Effects

In common fluids, viscosity is associated with dissipation. However, when time-reversal symmetry is broken a new type of nondissipative "viscosity" emerges. Recent theories and experiments on classical 2D systems with active spinning particles have heightened interest in "odd viscosity," but a microscopic theory for it in active materials is still absent. Here, we present such first-principles microscopic Hamiltonian theory, valid for both 2D and 3D, showing that odd viscosity is present in any system, even at zero temperature, with globally or locally aligned spinning components. Our work substantially extends the applicability of odd viscosity into 3D fluids, and specifically to internally driven active materials, such as living matter (e.g., actomyosin gels). We find intriguing 3D effects of odd viscosity such as propagation of anisotropic bulk shear waves and breakdown of Bernoulli's principle.

Hydrodynamics and direction change of tumbling bacteria

The bacterium Escherichia coli (E. coli) swims in viscous fluids by rotating several helical flagellar filaments, which are gathered in a bundle behind the cell during ‘runs’ wherein the cell moves steadily forward. In between runs, the cell undergoes quick ‘tumble’ events, during which at least one flagellum reverses its rotation direction and separates from the bundle, resulting in erratic motion in place and a random reorientation of the cell. Alternating between runs and tumbles allows cells to sample space by stochastically changing their propulsion direction after each tumble. The change of direction during a tumble is not uniformly distributed but is skewed towards smaller angles with an average of about 62°–68°, as first measured by Berg and Brown (1972). Here we develop a theoretical approach to model the angular distribution of swimming E. coli cells during tumbles. We first use past experimental imaging results to construct a kinematic description of the dynamics of the flagellar filaments during a tumble. We then employ low-Reynolds number hydrodynamics to compute the consequences of the kinematic model on the force and torque balance of the cell and to deduce the overall change in orientation. The results of our model are in good agreement with experimental observations. We find that the main change of direction occurs during the ‘bundling’ part of the process wherein, at the end of a tumble, the dispersed flagellar filaments are brought back together in the helical bundle, which we confirm using a simplified forced-sphere model.

Enhanced diffusion of a tracer particle in a lattice model of a crowded active system

Living systems at the subcellular, cellular, and multicellular levels are often crowded systems that contain active particles. The active motion of these particles can also propel passive particles, which typically results in enhanced effective diffusion of the passive particles. Here we study the diffusion of a passive tracer particle in such a dense system of active crowders using a minimal lattice model incorporating particles pushing each other. We show that the model exhibits several regimes of motility and quantify the enhanced diffusion as a function of density and activity of the active crowders. Moreover, we demonstrate an interplay of tracer diffusion and clustering of active particles, which suppresses the enhanced diffusion. Simulations of mixtures of passive and active crowders show that a rather small fraction of active particles is sufficient for the observation of enhanced diffusion.

Unraveling hidden interactions in complex systems with deep learning

Rich phenomena from complex systems have long intrigued researchers, and yet modeling system micro-dynamics and inferring the forms of interaction remain challenging for conventional data-driven approaches, being generally established by scientists with human ingenuity. In this study, we propose AgentNet, a model-free data-driven framework consisting of deep neural networks to reveal and analyze the hidden interactions in complex systems from observed data alone. AgentNet utilizes a graph attention network with novel variable-wise attention to model the interaction between individual agents, and employs various encoders and decoders that can be selectively applied to any desired system. Our model successfully captured a wide variety of simulated complex systems, namely cellular automata (discrete), the Vicsek model (continuous), and active Ornstein-Uhlenbeck particles (non-Markovian) in which, notably, AgentNet's visualized attention values coincided with the true variable-wise interaction strengths and exhibited collective behavior that was absent in the training data. A demonstration with empirical data from a flock of birds showed that AgentNet could identify hidden interaction ranges exhibited by real birds, which cannot be detected by conventional velocity correlation analysis. We expect our framework to open a novel path to investigating complex systems and to provide insight into general process-driven modeling.

Design of nematic liquid crystals to control microscale dynamics

The dynamics of small particles, both living such as swimming bacteria and inanimate, such as colloidal spheres, has fascinated scientists for centuries. If one could learn how to control and streamline their chaotic motion, that would open technological opportunities in the transformation of stored or environmental energy into systematic motion, with applications in micro-robotics, transport of matter, guided morphogenesis. This review presents an approach to command microscale dynamics by replacing an isotropic medium with a liquid crystal. Orientational order and associated properties, such as elasticity, surface anchoring, and bulk anisotropy, enable new dynamic effects, ranging from the appearance and propagation of particle-like solitary waves to self-locomotion of an active droplet. By using photoalignment, the liquid crystal can be patterned into predesigned structures. In the presence of the electric field, these patterns enable the transport of solid and fluid particles through nonlinear electrokinetics rooted in anisotropy of conductivity and permittivity. Director patterns command the dynamics of swimming bacteria, guiding their trajectories, polarity of swimming, and distribution in space. This guidance is of a higher level of complexity than a simple following of the director by rod-like microorganisms. Namely, the director gradients mediate hydrodynamic interactions of bacteria to produce an active force and collective polar modes of swimming. The patterned director could also be engraved in a liquid crystal elastomer. When an elastomer coating is activated by heat or light, these patterns produce a deterministic surface topography. The director gradients define an activation force that shapes the elastomer in a manner similar to the active stresses triggering flows in active nematics. The patterned elastomer substrates could be used to define the orientation of cells in living tissues. The liquid-crystal guidance holds a major promise in achieving the goal of commanding microscale active flows.

Chemotactic migration of bacteria in porous media

Chemotactic migration of bacteria-their ability to direct multicellular motion along chemical gradients-is central to processes in agriculture, the environment, and medicine. However, current understanding of migration is based on studies performed in bulk liquid, despite the fact that many bacteria inhabit tight porous media such as soils, sediments, and biological gels. Here, we directly visualize the chemotactic migration of Escherichia coli populations in well-defined 3D porous media in the absence of any other imposed external forcing (e.g., flow). We find that pore-scale confinement is a strong regulator of migration. Strikingly, cells use a different primary mechanism to direct their motion in confinement than in bulk liquid. Furthermore, confinement markedly alters the dynamics and morphology of the migrating population-features that can be described by a continuum model, but only when standard motility parameters are substantially altered from their bulk liquid values to reflect the influence of pore-scale confinement. Our work thus provides a framework to predict and control the migration of bacteria, and active matter in general, in complex environments.

Smoothed profile method for direct numerical simulations of hydrodynamically interacting particles

A general method is presented for computing the motions of hydrodynamically interacting particles in various kinds of host fluids for arbitrary Reynolds numbers. The method follows the standard procedure for performing direct numerical simulations (DNS) of particulate systems, where the Navier-Stokes equation must be solved consistently with the motion of the rigid particles, which defines the temporal boundary conditions to be satisfied by the Navier-Stokes equation. The smoothed profile (SP) method provides an efficient numerical scheme for coupling the continuum fluid mechanics with the dispersed moving particles, which are allowed to have arbitrary shapes. In this method, the sharp boundaries between solid particles and the host fluid are replaced with a smeared out thin shell (interfacial) region, which can be accurately resolved on a fixed Cartesian grid utilizing a SP function with a finite thickness. The accuracy of the SP method is illustrated by comparison with known exact results. In the present paper, the high degree of versatility of the SP method is demonstrated by considering several types of active and passive particle suspensions.

Spontaneous Motion and Rotation of Acid Droplets on the Surface of a Liquid Metal

Self-propulsion of droplets is of great significance in many fields. The spontaneous horizontal motion and self-jumping of droplets have been well realized in various ways. However, there is still a lack of an effective method to enable a droplet to rotate spontaneously and steadily. In this paper, by employing an acid droplet and a liquid metal, the spontaneous and steady rotation of droplets is achieved. For an acid droplet, it may spontaneously move when it is deposited on the surface of the liquid metal. By adjusting experimental parameters to the proper range, the self-rotation of droplet happens. This phenomenon originates from the fluctuation of the droplet boundary and the collective movement of bubbles that are generated by the chemical reactions between the acid droplet and liquid metal. This rotation has a simpler implementation method and more steady rotation state. Its angular velocity is much higher than that driven by other mechanisms. Moreover, the movements of acid droplets on the liquid metal are classified according to experimental conditions. The internal flow fields, the movements and distribution of bubbles, and the fluctuation of the droplet boundary are also explored and discussed. The theoretical model describing the rotational droplet is given. Our work may deepen the understanding of the physical system transition affected by chemical reactions and provide a new way for the design of potential applications, e.g., micro- and nanodevices.