Superfluid Behavior of Active Suspensions from Diffusive Stretching

Microrheology of active suspensions

We study the microrheology of active suspensions through direct hydrodynamic simulations using model pusher-like microswimmers. We demonstrate that the friction coefficient of a probe particle is notably reduced by hydrodynamic interactions (HIs) among a moving probe and the swimmers. When a swimmer approaches a probe from the rear (front) side, the repulsive HIs between them are weakened (intensified), which results in a slight front-rear asymmetry in swimmer orientation distribution around the probe, creating a significant additional net driving force acting on the probe from the rear side. The present drag-reduction mechanism qualitatively differs from that of the viscosity-reduction observed in sheared bulk systems and depends on probing details. This study provides insights into our fundamental knowledge of hydrodynamic effects in active suspensions and serves as a practical example illuminating distinctions between micro- and macrorheology measurements.

Thinning by cluster breaking: Active matter and shear flows share thinning mechanisms

Among many unexpected phenomena of active matter is the recently observed superfluid-like thinning (viscosity drop) behavior of bacteria suspensions. Understanding this peculiar self-propelled thinning by active matter is of theoretical and practical importance. Here, we find that, although distinct in driving mechanisms, active matter and shear flows exhibit similar thinning behaviors upon the increase of self-propulsion and shear forces, respectively. Our structural characterizations reveal that they actually share the same cluster-breaking mechanism of thinning. How fast and how shattered the cluster is broken determines the (dis)continuity of the thinning. This explains why adding active particles to Newtonian fluids can cause thinning, in which rotation of active particles play a key role in breaking clusters. Our work proposes a mechanism of self-propelled thinning and further establishes the underlying connections between active matter and shear flows.

Colloidal Self-Assembly: From Passive to Active Systems

Self-assembly fundamentally implies the organization of small sub-units into large structures or patterns without the intervention of specific local interactions. This process is commonly observed in nature, occurring at various scales ranging from atomic/molecular assembly to the formation of complex biological structures. Colloidal particles may serve as micrometer-scale surrogates for studying assembly, particularly for the poorly understood kinetic and dynamic processes at the atomic scale. Recent advances in colloidal self-assembly have enabled the programmable creation of novel materials with tailored properties. We here provide an overview and comparison of both passive and active colloidal self-assembly, with a discussion on the energy landscape and interactions governing both types. In the realm of passive colloidal assembly, many impressive and important structures have been realized, including colloidal molecules, one-dimensional chains, two-dimensional lattices, and three-dimensional crystals. In contrast, active colloidal self-assembly, driven by optical, electric, chemical, or other fields, involves more intricate dynamic processes, offering more flexibility and potential new applications. A comparative analysis underscores the critical distinctions between passive and active colloidal assemblies, highlighting the unique collective behaviors emerging in active systems. These behaviors encompass collective motion, motility-induced phase segregation, and exotic properties arising from out-of-equilibrium thermodynamics. Through this comparison, we aim to identify the future opportunities in active assembly research, which may suggest new application domains.

Active Brownian particles in external force fields: Field-theoretical models, generalized barometric law, and programmable density patterns

We investigate the influence of external forces on the collective dynamics of interacting active Brownian particles in two as well as three spatial dimensions. Via explicit coarse graining, we derive predictive models, i.e., models that give a direct relation between the models' coefficients and the bare parameters of the system, that are applicable for space- and time-dependent external force fields. We study these models for the cases of gravity and harmonic traps. In particular, we derive a generalized barometric formula for interacting active Brownian particles under gravity that is valid for low to high concentrations and activities of the particles. Furthermore, we show that one can use an external harmonic trap to induce motility-induced phase separation in systems that, without external fields, remain in a homogeneous state. This finding makes it possible to realize programmable density patterns in systems of active Brownian particles. Our analytic predictions are found to be in very good agreement with Brownian dynamics simulations.

Active XY model on a substrate: Density fluctuations and phase ordering

We explore the generic long-wavelength properties of an active XY model on a substrate, consisting of a collection of nearly phase-ordered active XY spins in contact with a diffusing, conserved species, as a representative system of active spinners with a conservation law. The spins rotate actively in response to the local density fluctuations and local phase differences, on a solid substrate. We investigate this system by Monte Carlo simulations of an agent-based model, which we set up, complemented by the hydrodynamic theory for the system. We demonstrate that this system can phase-synchronize without any hydrodynamic interactions. Our combined numerical and analytical studies show that this model, when stable, displays hitherto unstudied scaling behavior: As a consequence of the interplay between the mobility, active rotation, and number conservation, such a system can be stable over a wide range of the model parameters characterized by a novel correspondence between the phase and density fluctuations. In different regions of the phase space where the phase-ordered system is stable, it displays generalized quasi-long-range order (QLRO): It shows phase ordering which is generically either logarithmically stronger than the conventional QLRO found in its equilibrium limit, together with "miniscule number fluctuations," or logarithmically weaker than QLRO along with "giant number fluctuations," showing a novel one-to-one correspondence between phase ordering and density fluctuations in the ordered states. Intriguingly, these scaling exponents are found to depend explicitly on the model parameters. We further show that in other parameter regimes there are no stable, ordered phases. Instead, two distinct types of disordered states with short-range phase order are found, characterized by the presence or absence of stable clusters of finite sizes. In a surprising connection, the hydrodynamic theory for this model also describes the fluctuations in a Kardar-Parisi-Zhang (KPZ) surface with a conserved species on it, or an active fluid membrane with a finite tension, without momentum conservation and a conserved species living on it. This implies the existence of stable fluctuating surfaces that are only logarithmically smoother or rougher than the Edward-Wilkinson surface at two dimensions (2D) can exist, in contrast to the 2D pure KPZ-like "rough" surfaces.

Boundary design regulates the diffusion of active matter in heterogeneous environments

Physical boundaries play a key role in governing the overall transport properties of nearby self-propelled particles. In this work, we develop dispersion theories and conduct Brownian dynamics simulations to predict the coupling between surface accumulation and effective diffusivity of active particles in boundary-rich media. We focus on three models that are well-understood for passive systems: particle transport in (i) an array of fixed volume-excluding obstacles; (ii) a pore with spatially heterogeneous width; and (iii) a tortuous path with kinks and corners. While the impact of these entropic barriers on passive particle transport is well established, we find that these classical models of porous media flows break down due to the unique interplay between activity and the microstructure of the internal geometry. We study the activity-induced slowdown of effective diffusivity by formulating a Smoluchowski description of long-time self diffusivity which contains contributions from the density and fluctuation fields of the active particles. Particle-based and finite element simulations corroborate this perspective and reveal important nonequilibrium considerations of active transport.

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.

Hydrodynamics and surface properties influence biofilm proliferation

A biofilm is an interface-associated colloidal dispersion of bacterial cells and excreted polymers in which microorganisms find protection from their environment. Successful colonization of a surface by a bacterial community is typically a detriment to human health and property. Insight into the biofilm life-cycle provides clues on how their proliferation can be suppressed. In this review, we follow a cell through the cycle of attachment, growth, and departure from a colony. Among the abundance of factors that guide the three phases, we focus on hydrodynamics and stratum properties due to the synergistic effect such properties have on bacteria rejection and removal. Cell motion, whether facilitated by the environment via medium flow or self-actuated by use of an appendage, drastically improves the survivability of a bacterium. Once in the vicinity of a stratum, a single cell is exposed to near-surface interactions, such as van der Waals, electrostatic and specific interactions, similarly to any other colloidal particle. The success of the attachment and the potential for detachment is heavily influenced by surface properties such as material type and topography. The growth of the colony is similarly guided by mainstream flow and the convective transport throughout the biofilm. Beyond the growth phase, hydrodynamic traction forces on a biofilm can elicit strongly non-linear viscoelastic responses from the biofilm soft matter. As the colony exhausts the means of survival at a particular location, a set of trigger signals activates mechanisms of bacterial release, a life-cycle phase also facilitated by fluid flow. A review of biofilm-relevant hydrodynamics and startum properties provides insight into future research avenues.

Predictive local field theory for interacting active Brownian spheres in two spatial dimensions

We present a predictive local field theory for the nonequilibrium dynamics of interacting active Brownian particles with a spherical shape in two spatial dimensions. The theory is derived by a rigorous coarse-graining starting from the Langevin equations that describe the trajectories of the individual particles. For high accuracy and generality of the theory, it includes configurational order parameters and derivatives up to infinite order. In addition, we discuss possible approximations of the theory and present reduced models that are easier to apply. We show that our theory contains popular models such as Active Model B+  as special cases and that it provides explicit expressions for the coefficients occurring in these and other, often phenomenological, models. As a further outcome, the theory yields an analytical expression for the density-dependent mean swimming speed of the particles. To demonstrate an application of the new theory, we analyze a simple reduced model of the lowest nontrivial order in derivatives, which is able to predict the onset of motility-induced phase separation of the particles. By a linear stability analysis, an analytical expression for the spinodal corresponding to motility-induced phase separation is obtained. This expression is evaluated for the case of particles interacting repulsively by a Weeks-Chandler-Andersen potential. The analytical predictions for the spinodal associated with these particles are found to be in very good agreement with the results of Brownian dynamics simulations that are based on the same Langevin equations as our theory. Furthermore, the critical point predicted by our analytical results agrees excellently with recent computational results from the literature.

Time-reversal symmetry breaking in two-dimensional nonequilibrium viscous fluids

We study the rheological signatures of departure from equilibrium in two-dimensional viscous fluids with and without internal spin. Under the assumption of isotropy, we provide the most general linear constitutive relations for stress and couple stress in terms of the velocity and spin fields. Invoking Onsager's regression hypothesis for fluctuations about steady states, we derive the Green-Kubo formulas relating the transport coefficients to time-correlation functions of the fluctuating stress. In doing so, we show that one of the nonequilibrium transport coefficients, the odd viscosity, requires time-reversal symmetry breaking in the case of systems without internal spin. However, the Green-Kubo relations for systems with internal spin also show that there is a possibility for nonvanishing odd viscosity even when time-reversal symmetry is preserved. Furthermore, we find that breakdown of equipartition in nonequilibrium steady states results in the decoupling of the two rotational viscosities relating the vorticity and the internal spin.

Microscopic origins of the swim pressure and the anomalous surface tension of active matter

The unique pressure exerted by active particles-the "swim" pressure-has proven to be a useful quantity in explaining many of the seemingly confounding behaviors of active particles. However, its use has also resulted in some puzzling findings including an extremely negative surface tension between phase separated active particles. Here, we demonstrate that this contradiction stems from the fact that the swim pressure is not a true pressure. At a boundary or interface, the reduction in particle swimming generates a net active force density-an entirely self-generated body force. The pressure at the boundary, which was previously identified as the swim pressure, is in fact an elevated (relative to the bulk) value of the traditional particle pressure that is generated by this interfacial force density. Recognizing this unique mechanism for stress generation allows us to define a much more physically plausible surface tension. We clarify the utility of the swim pressure as an "equivalent pressure" (analogous to those defined from electrostatic and gravitational body forces) and the conditions in which this concept can be appropriately applied.

Hydrodynamics control shear-induced pattern formation in attractive suspensions

Dilute suspensions of repulsive particles exhibit a Newtonian response to flow that can be accurately predicted by the particle volume fraction and the viscosity of the suspending fluid. However, such a description fails when the particles are weakly attractive. In a simple shear flow, suspensions of attractive particles exhibit complex, anisotropic microstructures and flow instabilities that are poorly understood and plague industrial processes. One such phenomenon, the formation of log-rolling flocs, which is ubiquitously observed in suspensions of attractive particles that are sheared while confined between parallel plates, is an exemplar of this phenomenology. Combining experiments and discrete element simulations, we demonstrate that this shear-induced structuring is driven by hydrodynamic coupling between the flocs and the confining boundaries. Clusters of particles trigger the formation of viscous eddies that are spaced periodically and whose centers act as stable regions where particles aggregate to form flocs spanning the vorticity direction. Simulation results for the wavelength of the periodic pattern of stripes formed by the logs and for the log diameter are in quantitative agreement with experimental observations on both colloidal and noncolloidal suspensions. Numerical and experimental results are successfully combined by means of rescaling in terms of a Mason number that describes the strength of the shear flow relative to the rupture force between contacting particles in the flocs. The introduction of this dimensionless group leads to a universal stability diagram for the log-rolling structures and allows for application of shear-induced structuring as a tool for assembling and patterning suspensions of attractive particles.

Active Suspensions have Nonmonotonic Flow Curves and Multiple Mechanical Equilibria

We point out unconventional mechanical properties of confined active fluids, such as bacterial suspensions, under shear. Using a minimal model of an active liquid crystal with no free parameters, we predict the existence of a window of bacteria concentration for which a suspension of E. Coli effectively behaves, at steady-state, as a negative viscosity fluid and reach a quantitative agreement with experimental measurements. Our theoretical analysis further shows that a negative apparent viscosity is due to a nonmonotonic local velocity profile, and it is associated with a nonmonotonic stress versus strain rate flow curve. This implies that fixed stress and fixed strain rate ensembles are not equivalent for active fluids.

Symmetric shear banding and swarming vortices in bacterial superfluids

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

Do hydrodynamic interactions affect the swim pressure?

We study the motion of a spherical active Brownian particle (ABP) of size a, moving with a fixed speed U0, and reorienting on a time scale τR in the presence of a confining boundary. Because momentum is conserved in the embedding fluid, we show that the average force per unit area on the boundary equals the bulk mechanical pressure P∞ = p∞f + Π∞, where p∞f is the fluid pressure and Π∞ is the particle pressure; this is true for active and passive particles alike regardless of how the particles interact with the boundary. As an example, we investigate how hydrodynamic interactions (HI) change the particle-phase pressure at the wall, and find that Πwall = n∞(kBT + ζ(Δ)U0l(Δ)/6), where ζ is the (Stokes) drag on the swimmer, l = U0τR is the run length, and Δ is the minimum gap size between the particle and the wall; as Δ → ∞ this is the familiar swim pressure [Takatori et al., Phys. Rev. Lett., 2014, 113, 1-5].

Effective viscosity of a suspension of flagellar-beating microswimmers: Three-dimensional modeling

Micro-organisms usually can swim in their liquid environment by flagellar or ciliary beating. In this numerical work, we analyze the influence of flagellar beating on the orbits of a swimming cell in a shear flow. We also calculate the effect of the flagellar beating on the rheology of a dilute suspension of microswimmers. A three-dimensional model is proposed for Chlamydomonas Reinhardtii swimming with a breaststroke-like beating of two anterior flagella modeled by two counter-rotating fore beads. The active swimmer model reveals unusual angular orbits in a linear shear flow. Namely, the swimmer sustains orientations transiently across the flow. Such behavior is a result of the interplay between shear flow and the swimmer's periodic beating motion of flagella, which exert internal torques on the cell body. This peculiar behavior has some significant consequences on the rheological properties of the suspension. We calculate Einstein's viscosity of the suspension composed of such isolated modeled microswimmers (dilute case) in a shear flow. We use numerical simulations based on a Rotne-Prager-like approximation for hydrodynamic interaction between simplified flagella and the cell body. The results show an increased intrinsic viscosity for active swimmer suspensions in comparison to nonactive ones as well as a shear thinning behavior in accordance with previous experimental measurements [Phys. Rev. Lett. 104, 098102 (2010)10.1103/PhysRevLett.104.098102].

Effective viscosity of a two-dimensional suspension of interacting active particles

Suspensions of hydrodynamical active particles exhibit interesting rheological properties. For a dilute suspension of microswimmers, it has been shown that the effective viscosity of the suspension depends on the volume fraction of swimmers, and it behaves differently for pushers and pullers. Here we develop a theoretical framework to study the rheological properties of an interacting suspension. Taking into account the hydrodynamic interaction between swimmers and considering the small Péclet number condition, we calculate the effective viscosity of a two-dimensional suspension. For a dilute suspension, a perturbative result is obtained up to the second order of the surface fraction of swimmers. Our results show that the effective viscosity for the suspension can be very different for pushers and pullers.