Sheared active fluids: thickening, thinning, and vanishing viscosity

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.

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.

Flow states of two dimensional active gels driven by external shear

Using a minimal hydrodynamic model, we theoretically and computationally study the Couette flow of active gels in straight and annular two-dimensional channels subject to an externally imposed shear. The gels are isotropic in the absence of externally- or activity-driven shear, but have nematic order that increases with shear rate. Using the finite element method, we determine the possible flow states for a range of activities and shear rates. Linear stability analysis of an unconfined gel in a straight channel shows that an externally imposed shear flow can stabilize an extensile fluid that would be unstable to spontaneous flow in the absence of the shear flow, and destabilize a contractile fluid that would be stable against spontaneous flow in the absence of shear flow. These results are in rough agreement with the stability boundaries between the base shear flow state and the nonlinear flow states that we find numerically for a confined active gel. For extensile fluids, we find three kinds of nonlinear flow states in the range of parameters we study: unidirectional flows, oscillatory flows, and dancing flows. To highlight the activity-driven spontaneous component of the nonlinear flows, we characterize these states by the average volumetric flow rate and the wall stress. For contractile fluids, we only find the linear shear flow and a nonlinear unidirectional flow in the range of parameters that we studied. For large magnitudes of the activity, the unidirectional contractile flow develops a boundary layer. Our analysis of annular channels shows how curvature of the streamlines in the base flow affects the transitions among flow states.

Motility-induced shear thickening in dense colloidal suspensions

Phase transitions and collective dynamics of active colloidal suspensions are fascinating topics in soft matter physics, particularly for out-of-equilibrium systems, which can lead to rich rheological behaviours in the presence of steady shear flow. Here the role of self-propulsion in the rheological response of a dense colloidal suspension is investigated by using particle-resolved Brownian dynamics simulations. First, the combined effect of activity and shear in the solid on the disordering transition of the suspension is analyzed. While both self-propulsion and shear destroy order and melt the system if critical values are exceeded, self-propulsion largely lowers the stress barrier needed to be overcome during the transition. We further explore the rheological response of the active sheared system once a steady state is reached. While passive suspensions show a solid-like behaviour, turning on particle motility fluidises the system. At low self-propulsion, the active suspension behaves in the steady state as a shear-thinning fluid. Increasing the self-propulsion changes the behaviour of the liquid from shear-thinning to shear-thickening. We attribute this to clustering in the sheared suspensions induced by motility. This new phenomenon of motility-induced shear thickening (MIST) can be used to tailor the rheological response of colloidal suspensions.

Coexistence of Active and Hydrodynamic Turbulence in Two-Dimensional Active Nematics

In active nematic liquid crystals, activity is able to drive chaotic spatiotemporal flows referred to as active turbulence. Active turbulence has been characterized through theoretical and experimental work as a low Reynolds number phenomenon. We show that, in two dimensions, the active forcing alone is able to trigger hydrodynamic turbulence leading to the coexistence of active and inertial turbulence. This type of flow develops for sufficiently active and extensile flow-aligning nematics. We observe that the combined effect of an extensile nematic and large values of the flow-aligning parameter leads to a broadening of the elastic energy spectrum that promotes a growth of kinetic energy able to trigger an inverse energy cascade.

Defect dynamics in active polar fluids vs. active nematics

Topological defects play a key role in two-dimensional active nematics, and a transient role in two-dimensional active polar fluids. Using a variational method, we study both the transient and long-time behavior of defects in two-dimensional active polar fluids in the limit of strong order and overdamped, compressible flow, and compare the defect dynamics with the corresponding active nematics model studied recently. One result is non-central interactions between defect pairs for active polar fluids, and by extending our analysis to allow orientation dynamics of defects, we find that the orientation of +1 defects, unlike that of ±1/2 defects in active nematics, is not locked to defect positions and relaxes to asters. Moreover, using a scaling argument, we explain the transient feature of active polar defects and show that in the steady state, active polar fluids are either devoid of defects or consist of a single aster. We argue that for contractile (extensile) active nematic systems, +1 vortices (asters) should emerge as bound states of a pair of +1/2 defects, which has been recently observed. Moreover, unlike the polar case, we show that for active nematics, a linear chain of equally spaced bound states of pairs of +1/2 defects can screen the activity term. A common feature in both models is the appearance of +1 defects (elementary in polar and composite in nematic) in the steady state.

Active nematic defects in compressible and incompressible flows

We study the dynamics of active nematic films on a substrate driven by active flows with or without the incompressible constraint. Through simulations and theoretical analysis, we show that arch patterns are stable in the compressible case, while they become unstable under the incompressibility constraint. For compressible flows at high enough activity, stable arches organize themselves into a smecticlike pattern, which induce an associated global polar ordering of +1/2 nematic defects. By contrast, divergence-free flows give rise to a local nematic order of the +1/2 defects, consisting of antialigned pairs of neighboring defects, as established in previous studies.

Active nematic gel with quenched disorder

With quenched disorder, we introduce two-dimensional active nematics suspended in an incompressible fluid. We write the coarse-grained hydrodynamic equations of motion for slow variables, viz. density, orientation, and flow fields. The quenched disorder is introduced such that it interacts with the local orientation at every point with some strength. Disorder strength is tuned from zero to large values. We numerically study the defect dynamics and system kinetics and find that the finite disorder slows the ordering. The presence of fluid induces large fluctuation in the orientation field, further disturbing the ordering. The large fluctuation in the orientation field due to the fluid is so dominant that it reduces the effect of the quenched disorder. We have also found that the disorder effect is almost the same for both the contractile and extensile nature of active stresses in the system. This study can help to understand the impact of quenched disorder on the ordering kinetics of active gels with nematic interaction among the constituent objects.

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.

Shear-induced orientational ordering in an active glass former

Dense assemblies of self-propelled particles that can form solid-like states also known as active or living glasses are abundant around us, covering a broad range of length scales and timescales: from the cytoplasm to tissues, from bacterial biofilms to vehicular traffic jams, and from Janus colloids to animal herds. Being structurally disordered as well as strongly out of equilibrium, these systems show fascinating dynamical and mechanical properties. Using extensive molecular dynamics simulation and a number of distinct dynamical and mechanical order parameters, we differentiate three dynamical steady states in a sheared model active glassy system: 1) a disordered state, 2) a propulsion-induced ordered state, and 3) a shear-induced ordered state. We supplement these observations with an analytical theory based on an effective single-particle Fokker-Planck description to rationalize the existence of the shear-induced orientational ordering behavior in an active glassy system without explicit aligning interactions of, for example, Vicsek type. This ordering phenomenon occurs in the large persistence time limit and is made possible only by the applied steady shear. Using a Fokker-Planck description with parameters that can be measured independently, we make testable predictions for the joint distribution of single-particle position and orientation. These predictions match well with the joint distribution measured from direct numerical simulation. Our results are of relevance for experiments exploring the rheological response of dense active colloids and jammed active granular matter systems.

Effect of polydispersity on the dynamics of active Brownian particles

We numerically study the dynamics and the phases of self-propelled disk-shaped particles of different sizes with soft repulsive potential in two dimensions. Size diversity is introduced by the polydispersity index (PDI) ε, which is the width of the uniform distribution of the particle's radius. The self-propulsion speed of the particles controls the activity v. We observe enhanced dynamics for large size diversity among the particles. We calculate the effective diffusion coefficient D_{eff} in the steady state. The system exhibits four distinct phases, jammed phase with small D_{eff} for small activity and liquid phase with enhanced D_{eff} for large activity. The number fluctuation is larger and smaller than the equilibrium limit in the liquid and jammed phases, respectively. Further, the jammed phase is of two types: solid jammed and liquid jammed for small and large PDI. Whereas the liquid phase is called motility induced phase separation (MIPS) liquid for small PDI and for large PDI, we find enhanced diffusivity and call it the pure liquid phase. The system is studied for three packing densities ϕ, and the response of the system for polydispersity is the same for all ϕ's. Our study can help understand the behavior of cells of various sizes in a tissue, artificial self-driven granular particles, or living organisms of different sizes in a dense environment.

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.

Fluctuations can induce local nematic order and extensile stress in monolayers of motile cells

Recent experiments in various cell types have shown that two-dimensional tissues often display local nematic order, with evidence of extensile stresses manifest in the dynamics of topological defects. Using a mesoscopic model where tissue flow is generated by fluctuating traction forces coupled to the nematic order parameter, we show that the resulting tissue dynamics can spontaneously produce local nematic order and an extensile internal stress. A key element of the model is the assumption that in the presence of local nematic alignment, cells preferentially crawl along the nematic axis, resulting in anisotropy of fluctuations. Our work shows that activity can drive either extensile or contractile stresses in tissue, depending on the relative strength of the contractility of the cortical cytoskeleton and tractions by cells on the extracellular matrix.

Cell Junction Mechanics beyond the Bounds of Adhesion and Tension

Cell-cell junctions, in particular adherens junctions, are major determinants of tissue mechanics during morphogenesis and homeostasis. In attempts to link junctional mechanics to tissue mechanics, many have utilized explicitly or implicitly equilibrium approaches based on adhesion energy, surface energy, and contractility to determine the mechanical equilibrium at junctions. However, it is increasingly clear that they have significant limitations, such as that it remains challenging to link the dynamics of the molecular components to the resulting physical properties of the junction, to its remodeling ability, and to its adhesion strength. In this perspective, we discuss recent attempts to consider the aspect of energy dissipation at junctions to draw contact points with soft matter physics where energy loss plays a critical role in adhesion theories. We set the grounds for a theoretical framework of the junction mechanics that bridges the dynamics at the molecular scale to the mechanics at the tissue scale.

Bacterial biomechanics-From individual behaviors to biofilm and the gut flora

Bacteria inhabit a variety of locations and play important roles in the environment and health. Our understanding of bacterial biomechanics has improved markedly in the last decade and has revealed that biomechanics play a significant role in microbial biology. The obtained knowledge has enabled investigation of complex phenomena, such as biofilm formation and the dynamics of the gut flora. A bottom-up strategy, i.e., from the cellular to the macroscale, facilitates understanding of macroscopic bacterial phenomena. In this Review, we first cover the biomechanics of individual bacteria in the bulk liquid and on surfaces as the base of complex phenomena. The collective behaviors of bacteria in simple environments are next introduced. We then introduce recent advances in biofilm biomechanics, in which adhesion force and the flow environment play crucial roles. We also review transport phenomena in the intestine and the dynamics of the gut flora, focusing on that in zebrafish. Finally, we provide an overview of the future prospects for the field.

Pressure-driven changes to spontaneous flow in active nematic liquid crystals

We consider the effects of a pressure gradient on the spontaneous flow of an active nematic liquid crystal in a channel, subject to planar anchoring and no-slip conditions on the boundaries of the channel. We employ a model based on the Ericksen-Leslie theory of nematics, with an additional active stress accounting for the activity of the fluid. By directly solving the flow equation, we consider an asymptotic solution for the director angle equation for large activity parameter values and predict the possible values of the director angle in the bulk of the channel. Through a numerical solution of the full nonlinear equations, we examine the effects of pressure on the branches of stable and unstable equilibria, some of which are disconnected from the no-flow state. In the absence of a pressure gradient, solutions are either symmetric or antisymmetric about the channel midpoint; these symmetries are changed by the pressure gradient. Considering the activity-pressure state space allows us to predict qualitatively the extent of each solution type and to show, for large enough pressure gradients, that a branch of non-trivial director angle solutions exists for all activity values.

The inhibition of concentrated active baths

Passive tracers in the active bath express fascinating behaviors. However, most studies are restricted to dilute active baths. Here, we use 2D simulation of suspensions consisting of active Brownian particles and a passive disk-shaped tracer to investigate tracers' diffusive behaviors in a wide range of volume fractions. Due to the competition between the thermal noise and collisions with active particles, tracers express a first transition from the normal diffusion to the superdiffusion at a short time scale and recur to normal diffusion at a long time scale. At a low volume fraction, infrequent active collisions retard the first transition of smaller tracers. At a high volume fraction, active particles with high activity aggregating around tracers induce a bimodal probability distribution function of tracer displacements during superdiffusion. Considering the enhancement of diffusion, the non-dimensional enhanced diffusivity increases asymptotically with the Peclet number. The asymptotic line gives an upper limit of non-dimensional enhanced diffusivity of tracers. Cases with lower enhanced diffusion have a high volume fraction and a low active velocity that indicates the inhibition of concentrated active baths. With the high negentropic work of these cases, the inhibition is explained as the change of the configuration of active baths for introducing tracers.

Different-shaped micro-objects driven by active particle aggregations

The dynamics of passive micro-objects in an active bath has been receiving much attention. However, the influence of the shapes of micro-objects remains unclear. Here, we use 2D simulation to investigate the interaction between active Brownian particles and different-shaped passive micro-objects. We show that active particles accumulate around micro-objects and self-assemble into living aggregations at a high active velocity and high volume fraction. The shapes of micro-objects affect the distributions of the aggregations. In turn, the different distribution of aggregations influences the motion of micro-objects and induces abnormal diffusive behaviors. We further demonstrate that polar distributed aggregations at a high active velocity and the inhibition of the active bath at a low active velocity induce the counterintuitive anisotropic enhanced diffusion of rods, and the steric interaction between active particles induces the reverse translation-rotation coupled diffusion of chevrons.

Darcy's Law without Friction in Active Nematic Rheology

We study the dynamics of a contractile active nematic fluid subjected to a Poiseuille flow. In a quasi-1D geometry, we find that the linear rheology of this material is reminiscent of Darcy's law in complex fluids, with a pluglike flow decaying to zero over a well-defined "permeation" length. As a result, the viscosity increases with size, but never diverges, thereby evading the yield stress predicted by previous theories. We find strong shear thinning controlled by an active Ericksen number quantifying the ratio between external pressure difference and internal active stresses. In 2D, the increase of linear regime viscosity with size only persists up to a critical length beyond which we observe active turbulent patterns, with very low apparent viscosity. The ratio between the critical and permeation length determining the stability of the Darcy regime can be made indefinitely large by varying the flow aligning parameter or magnitude of nematic order.

Rheology of Entangled Active Polymer-Like T. Tubifex Worms

We experimentally study the rheology of long, slender, and entangled living worms (Tubifex Tubifex). Their level of activity can be controlled by changing the temperature or by adding small amounts of alcohol to make the worms temporarily inactive. Performing classical rheology experiments on this entangled polymer-like system, we find that the rheology is qualitatively similar to that of usual polymers, but, quantitatively, (i) shear thinning is reduced by activity, (ii) the characteristic shear rate for the onset of shear-thinning is given by the time scale of the activity, and (iii) the low shear viscosity as a function of concentration shows a very different scaling from that of regular polymers. Our study paves the way towards a new experimental research field of active "polymer-like worms."

Negative viscosity of liquid crystals in the presence of turbulence: Conductivity dependence, phase diagram, and self-oscillation

Recently, we reported the discovery of enormous negative viscosity of a nematic liquid crystal in the presence of turbulence induced by ac electric fields, which enabled us to observe unique phenomena related to the negative viscosity, such as spontaneous shear flow, hysteresis in flow curves, and self-oscillation [Orihara et al., Phys. Rev. E 99, 012701 (2019)10.1103/PhysRevE.99.012701]. In the present paper, we report the rheological properties of another nematic liquid crystal, which is a homologue of the previous one. The properties of the present liquid crystal are strongly dependent on electrical conductivity. Three samples with different conductivities were prepared by changing the amount of an ionic dopant. It was found that the lowest-conductivity sample without dopant shows no negative viscosity whereas the other ion-doped samples exhibit negative viscosity with strong dependence on the frequency of the ac electric field, consistent with microscopic observations. Phase diagrams of the negative- and positive-viscosity states in the amplitude and frequency plane are constructed to show the conductivity effect. Furthermore, we propose a model to reproduce another type of self-oscillation found in the present study.

Swimmer Suspensions on Substrates: Anomalous Stability and Long-Range Order

We present a comprehensive theory of the dynamics and fluctuations of a two-dimensional suspension of polar active particles in an incompressible fluid confined to a substrate. We show that, depending on the sign of a single parameter, a state with polar orientational order is anomalously stable (or anomalously unstable), with a nonzero relaxation (or growth) rate for angular fluctuations, not parallel to the ordering direction, at zero wave number. This screening of the broken-symmetry mode in the stable state does lead to conventional rather than giant number fluctuations as argued by Bricard et al., Nature 503, 95 (2013), but their bend instability in a splay-stable flock does not exist and the polar phase has long-range order in two dimensions. Our theory also describes confined three-dimensional thin-film suspensions of active polar particles as well as dense compressible active polar rods, and predicts a flocking transition without a banding instability.

Rheology of active polar emulsions: from linear to unidirectional and inviscid flow, and intermittent viscosity

The rheological behaviour of an emulsion made of an active polar component and an isotropic passive fluid is studied by lattice Boltzmann methods. Different flow regimes are found by varying the values of the shear rate and extensile activity (occurring, e.g., in microtubule-motor suspensions). By increasing the activity, a first transition occurs from the linear flow regime to spontaneous persistent unidirectional macro-scale flow, followed by another transition either to a (low shear) intermittent flow regime with the coexistence of states with positive, negative, and vanishing apparent viscosity, or to a (high shear) symmetric shear thinning regime. The different behaviours can be explained in terms of the dynamics of the polarization field close to the walls. A maximum entropy production principle selects the most likely states in the intermittent regime.

Dry active turbulence in a model for microtubule-motor mixtures

We study the dynamics and phase behaviour of a dry suspension of microtubules and molecular motors. We obtain a set of continuum equations by rigorously coarse graining a microscopic model where motor-induced interactions lead to parallel or antiparallel ordering. Through numerical simulations, we show that this model generically creates either stable stripes, or a never-settling pattern where stripes periodically form, rotate and then split up. We derive a minimal model which displays the same instability as the full model, and clarifies the underlying physical mechanism. The necessary ingredients are an extensile flux arising from microtubule sliding and an interfacial torque favouring ordering along density gradients. We argue that our minimal model unifies various previous observations of chaotic behaviour in dry active matter into a general universality class.

Lattice Boltzmann methods and active fluids

We review the state of the art of active fluids with particular attention to hydrodynamic continuous models and to the use of Lattice Boltzmann Methods (LBM) in this field. We present the thermodynamics of active fluids, in terms of liquid crystals modelling adapted to describe large-scale organization of active systems, as well as other effective phenomenological models. We discuss how LBM can be implemented to solve the hydrodynamics of active matter, starting from the case of a simple fluid, for which we explicitly recover the continuous equations by means of Chapman-Enskog expansion. Going beyond this simple case, we summarize how LBM can be used to treat complex and active fluids. We then review recent developments concerning some relevant topics in active matter that have been studied by means of LBM: spontaneous flow, self-propelled droplets, active emulsions, rheology, active turbulence, and active colloids.

Phase Transition to Large Scale Coherent Structures in Two-Dimensional Active Matter Turbulence

The collective motion of microswimmers in suspensions induce patterns of vortices on scales that are much larger than the characteristic size of a microswimmer, attaining a state called bacterial turbulence. Hydrodynamic turbulence acts on even larger scales and is dominated by inertial transport of energy. Using an established modification of the Navier-Stokes equation that accounts for the small-scale forcing of hydrodynamic flow by microswimmers, we study the properties of a dense suspension of microswimmers in two dimensions, where the conservation of enstrophy can drive an inverse cascade through which energy is accumulated on the largest scales. We find that the dynamical and statistical properties of the flow show a sharp transition to the formation of vortices at the largest length scale. The results show that 2D bacterial and hydrodynamic turbulence are separated by a subcritical phase transition.

Exact results for sheared polar active suspensions with variable liquid crystalline order

We consider a confined sheared active polar liquid crystal with a uniform orientation and study the effect of variations in the magnitude of polarization. Restricting our analysis to one-dimensional geometries, we demonstrate that with asymmetric boundary conditions, this system is characterized, macroscopically, by a linear shear stress vs. shear strain relationship that does not pass through the origin: At a zero strain rate, the fluid sustains a non-zero stress. Analytic solutions for the polarization, density, and velocity fields are derived for asymptotically large or small systems and are shown by comparison with precise numerical solutions to be good approximations for finite-size systems.

Shear-Induced First-Order Transition in Polar Liquid Crystals

The hydrodynamic theory of polar liquid crystals is widely used to describe biological active fluids as well as passive molecular materials. Depending on the "shear-alignment parameter", in passive or weakly active polar fluids under external shear, the polar order parameter p is either inclined to the flow at a fixed (Leslie) angle, or rotates continuously. Here, we study the role of an additional "shear-elongation parameter" that has been neglected in the recent literature and causes |p| to change under flow. We show that this effect can give rise to a shear-induced first-order phase transition from isotropic to polar, and significantly change the rheological properties of both active and passive polar fluids.

Negative viscosity of a liquid crystal in the presence of turbulence

We report on the discovery of enormous negative viscosity in a nematic liquid crystal in the presence of turbulence induced by electric fields. As the negative viscosity in this system is so large, we are able to observe several phenomena originating from it. For example, we observe a spontaneous shear flow that rotates the upper disk of a rheometer, as well as the reversal of the rotational direction upon applying an external torque in the opposite direction. Hysteresis loops are also observed in the shear-stress-shear-rate curves, which is reminiscent of those seen for ferromagnetic and ferroelectric materials. The similarities between the phenomena observed for our system and ferroic materials are comprehensively demonstrated, although the two systems are fundamentally different in that the former is out of equilibrium. We elucidate the origin of the negative viscosity and propose a simple model that reproduces the phenomena observed in this active fluid.

Enhancement of Microorganism Swimming Speed in Active Matter

We study a swimming undulating sheet in the isotropic phase of an active nematic liquid crystal. Activity changes the effective shear viscosity, reducing it to zero at a critical value of activity. Expanding in the sheet amplitude, we find that the correction to the swimming speed due to activity is inversely proportional to the effective shear viscosity. Our perturbative calculation becomes invalid near the critical value of activity; using numerical methods to probe this regime, we find that activity enhances the swimming speed by an order of magnitude compared to the passive case.

Effective temperature of active fluids and sheared soft glassy materials

The dynamics within active fluids, driven by internal activity of the self-propelled particles, is a subject of intense study in non-equilibrium physics. These systems have been explored using simulations, where the motion of a passive tracer particle is followed. Similar studies have been carried out for a soft glassy material that is driven by shearing its boundaries. In both types of systems the non-equilibrium motion have been quantified by defining a set of "effective temperatures", using both the tracer particle kinetic energy and the fluctuation-dissipation relation. We demonstrate that these effective temperatures extracted from the many-body simulations fit analytical expressions that are obtained for a single active particle inside a visco-elastic fluid. This result provides testable predictions and suggests a unified description for the dynamics inside active systems.

Kinetic theory of pattern formation in mixtures of microtubules and molecular motors

In this study we formulate a theoretical approach, based on a Boltzmann-like kinetic equation, to describe pattern formation in two-dimensional mixtures of microtubular filaments and molecular motors. Following the previous work by Aranson and Tsimring [Phys. Rev. E 74, 031915 (2006)PLEEE81539-375510.1103/PhysRevE.74.031915] we model the motor-induced reorientation of microtubules as collision rules, and devise a semianalytical method to calculate the corresponding interaction integrals. This procedure yields an infinite hierarchy of kinetic equations that we terminate by employing a well-established closure strategy, developed in the pattern-formation community and based on a power-counting argument. We thus arrive at a closed set of coupled equations for slowly varying local density and orientation of the microtubules, and study its behavior by performing a linear stability analysis and direct numerical simulations. By comparing our method with the work of Aranson and Tsimring, we assess the validity of the assumptions required to derive their and our theories. We demonstrate that our approximation-free evaluation of the interaction integrals and our choice of a systematic closure strategy result in a rather different dynamical behavior than was previously reported. Based on our theory, we discuss the ensuing phase diagram and the patterns observed.

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.

Fluidization and Active Thinning by Molecular Kinetics in Active Gels

We derive the constitutive equations of an active polar gel from a model for the dynamics of elastic molecules that link polar elements. Molecular binding kinetics induces the fluidization of the material, giving rise to Maxwell viscoelasticity and, provided that detailed balance is broken, to the generation of active stresses. We give explicit expressions for the transport coefficients of active gels in terms of molecular properties, including nonlinear contributions on the departure from equilibrium. In particular, when activity favors linker unbinding, we predict a decrease of viscosity with activity-active thinning-of kinetic origin, which could explain some experimental results on the cell cortex. By bridging the molecular and hydrodynamic scales, our results could help understand the interplay between molecular perturbations and the mechanics of cells and tissues.

Probing the shear viscosity of an active nematic film

In vitro reconstituted active systems, such as the adenosine triphosphate (ATP)-driven microtubule bundle suspension developed by the Dogic group [T. Sanchez, D. T. Chen, S. J. DeCamp, M. Heymann, and Z. Dogic, Nature (London) 491, 431 (2012)10.1038/nature11591], provide a fertile testing ground for elucidating the phenomenology of active liquid crystalline states. Controlling such novel phases of matter crucially depends on our knowledge of their material and physical properties. In this Rapid Communication, we show that the shear viscosity of an active nematic film can be probed by varying its hydrodynamic coupling to a bounding oil layer. Using the motion of disclinations as intrinsic tracers of the flow field and a hydrodynamic model, we obtain an estimate for the shear viscosity of the nematic film. Knowing this now provides us with an additional handle for robust and precision tunable control of the emergent dynamics of active fluids.

Dynamics of ellipsoidal tracers in swimming algal suspensions

Enhanced diffusion of passive tracers immersed in active fluids is a universal feature of active fluids and has been extensively studied in recent years. Similar to microrheology for equilibrium complex fluids, the unusual enhanced particle dynamics reveal intrinsic properties of active fluids. Nevertheless, previous studies have shown that the translational dynamics of spherical tracers are qualitatively similar, independent of whether active particles are pushers or pullers-the two fundamental classes of active fluids. Is it possible to distinguish pushers from pullers by simply imaging the dynamics of passive tracers? Here, we investigated the diffusion of isolated ellipsoids in algal C. reinhardtii suspensions-a model for puller-type active fluids. In combination with our previous results on pusher-type E. coli suspensions [Peng et al., Phys. Rev. Lett. 116, 068303 (2016)PRLTAO0031-900710.1103/PhysRevLett.116.068303], we showed that the dynamics of asymmetric tracers show a profound difference in pushers and pullers due to their rotational degree of freedom. Although the laboratory-frame translation and rotation of ellipsoids are enhanced in both pushers and pullers, similar to spherical tracers, the anisotropic diffusion in the body frame of ellipsoids shows opposite trends in the two classes of active fluids. An ellipsoid diffuses fastest along its major axis when immersed in pullers, whereas it diffuses slowest along the major axis in pushers. This striking difference can be qualitatively explained using a simple hydrodynamic model. In addition, our study on algal suspensions reveals that the influence of the near-field advection of algal swimming flows on the translation and rotation of ellipsoids shows different ranges and strengths. Our work provides not only new insights into universal organizing principles of active fluids, but also a convenient tool for detecting the class of active particles.

Negative stiffness and modulated states in active nematics

We examine the dynamics of an active nematic liquid crystal on a frictional substrate. When frictional damping dominates over viscous dissipation, we eliminate flow in favor of active stresses to obtain a minimal dynamical model for the nematic order parameter, with elastic constants renormalized by activity. The renormalized elastic constants can become negative at large activity, leading to the selection of spatially inhomogeneous patterns via a mechanism analogous to that responsible for modulated phases arising at an equilibrium Lifshitz point. Tuning activity and the degree of nematic order in the passive system, we obtain a linear stability phase diagram that exhibits a nonequilibrium tricritical point where ordered, modulated and disordered phases meet. Numerical solution of the nonlinear equations yields a succession of spatial structures of increasing complexity with increasing activity, including kink walls and active turbulence, as observed in experiments on microtubule bundles confined at an oil-water interface. Our work provides a minimal model for an overdamped active nematic that reproduces all the nonequilibrium structures seen in simulations of the full active nematic hydrodynamics and provides a framework for understanding some of the mechanisms for selection of the nonequilibrium patterns in the language of equilibrium critical phenomena.

Emergent self-organization in active materials

Biological systems exhibit large-scale self-organized dynamics and structures which enable organisms to perform the functions of life. The field of active matter strives to develop and understand microscopically driven nonequilibrium materials, with emergent properties comparable to those of living systems. This review will describe two recently developed classes of active matter systems, in which simple building blocks-self-propelled colloidal particles or extensile rod-like particles-self-organize to form macroscopic structures with features not possible in equilibrium systems. We summarize the recent experimental and theoretical progress on each of these systems, and we present simple descriptions of the physics underlying their emergent behaviors.

Turning Bacteria Suspensions into Superfluids

The rheological response under simple shear of an active suspension of Escherichia coli is determined in a large range of shear rates and concentrations. The effective viscosity and the time scales characterizing the bacterial organization under shear are obtained. In the dilute regime, we bring evidence for a low-shear Newtonian plateau characterized by a shear viscosity decreasing with concentration. In the semidilute regime, for particularly active bacteria, the suspension displays a "superfluidlike" transition where the viscous resistance to shear vanishes, thus showing that, macroscopically, the activity of pusher swimmers organized by shear is able to fully overcome the dissipative effects due to viscous loss.

Physics of microswimmers--single particle motion and collective behavior: a review

Locomotion and transport of microorganisms in fluids is an essential aspect of life. Search for food, orientation toward light, spreading of off-spring, and the formation of colonies are only possible due to locomotion. Swimming at the microscale occurs at low Reynolds numbers, where fluid friction and viscosity dominates over inertia. Here, evolution achieved propulsion mechanisms, which overcome and even exploit drag. Prominent propulsion mechanisms are rotating helical flagella, exploited by many bacteria, and snake-like or whip-like motion of eukaryotic flagella, utilized by sperm and algae. For artificial microswimmers, alternative concepts to convert chemical energy or heat into directed motion can be employed, which are potentially more efficient. The dynamics of microswimmers comprises many facets, which are all required to achieve locomotion. In this article, we review the physics of locomotion of biological and synthetic microswimmers, and the collective behavior of their assemblies. Starting from individual microswimmers, we describe the various propulsion mechanism of biological and synthetic systems and address the hydrodynamic aspects of swimming. This comprises synchronization and the concerted beating of flagella and cilia. In addition, the swimming behavior next to surfaces is examined. Finally, collective and cooperate phenomena of various types of isotropic and anisotropic swimmers with and without hydrodynamic interactions are discussed.

The role of filament length, finite-extensibility and motor force dispersity in stress relaxation and buckling mechanisms in non-sarcomeric active gels

After relaxing some assumptions we apply a single-chain mean-field mathematical model recently introduced [RSC Adv. (2014)] to describe the role of molecular motors in the mechanical properties of active gels. The model allows physics that are not available in models postulated on coarser levels of description. Moreover it proposes a level of description that allows the prediction of observables at time scales too difficult to achieve in multi-chain simulations for realistic filament lengths and densities. We model the semiflexible filaments that compose the active gel as bead-spring chains; molecular motors are accounted for by using a mean-field approach, in which filaments undergo transitions of one motor attachment state depending on the state of the probe filament. The level of description includes the end-to-end distance and attachment state of the filaments, and the motor-generated forces, as stochastic state variables which evolve according to a proposed differential Chapman-Kolmogorov equation. The motor-generated forces are drawn from a stationary distribution of motor stall forces. We consider bead-spring chains with multiple beads, explore the effect of finite-extensibility of the strands and incorporate into the model motor force distributions that have been measured experimentally. The model can no longer be solved analytically but is amenable to numerical simulation. This version of the model allows a more quantitative description of buckling dynamics [Lenz et. al. PRL, 2012, 108, 238107] and the dynamic modulus of active gels. The effect of finite extensibility of the filament strands on the dynamic modulus was also found to be in agreement with the microrheology experiments of Mizuno et. al., [Science, 2007, 315, 370-373].

Defect dynamics in active nematics

Topological defects are distinctive signatures of liquid crystals. They profoundly affect the viscoelastic behaviour of the fluid by constraining the orientational structure in a way that inevitably requires global changes not achievable with any set of local deformations. In active nematic liquid crystals, topological defects not only dictate the global structure of the director, but also act as local sources of motion, behaving as self-propelled particles. In this article, we present a detailed analytical and numerical study of the mechanics of topological defects in active nematic liquid crystals.

Spontaneous motility of passive emulsion droplets in polar active gels

We study by computer simulations the dynamics of a droplet of passive, isotropic fluid, embedded in a polar active gel. The latter represents a fluid of active force dipoles, which exert either contractile or extensile stresses on their surroundings, modelling for instance a suspension of cytoskeletal filaments and molecular motors. When the polarisation of the active gel is anchored normal to the droplet at its surface, the nematic elasticity of the active gel drives the formation of a hedgehog defect; this defect then drives an active flow which propels the droplet forward. In an extensile gel, motility can occur even with tangential anchoring, which is compatible with a defect-free polarisation pattern. In this case, upon increasing activity the droplet first rotates uniformly, and then undergoes a discontinuous nonequilibrium transition into a translationally motile state, powered by bending deformations in the surrounding active medium.

Capillary instability of axisymmetric, active liquid crystal jets

We study linear stability of an infinitely long, axisymmetric, cylindrical active liquid crystal (ALC) jet in a passive isotropic fluid matrix using a polar active liquid crystal (ALC) model. We identify three possible unstable modes (or mechanisms) as the result of the interaction between the flow and the active (or self-propelled) molecular motion. The first unstable mode is related to the polarity vector instability when coupled to the flow field in the presence of the molecular activity. It can be traced back to the inherent polarity vector instability in a bulk active liquid crystal flow. However, it can be grossly amplified in the ALC jet to encompass up to infinitely many unstable growth rates when the long range distortional elastic interaction is weak in certain parameter regimes; it can also be suppressed in other parameter regimes completely. The second unstable mode is related to the classical capillary or Rayleigh instability, which exists in a finite wave interval [0, k(cutoff)]. The new feature for this instability lies in the dependence of the cutoff wave number (k(cutoff)) on the activity of the active matter system. For ALC jets with sufficiently strong contractile activity, the instability can be completely suppressed though. The third unstable mode is due to the active viscous stress. This unstable mode can emerge in the intermediate wave number regime at a sufficiently strong active viscosity and even expand all the way to the zero wave number limit when the Rayleigh unstable mode is absent. It can also be suppressed in the regime of weak active viscous stress. At any given values of the model parameters, the three types of instabilities can show up either individually or in a certain combination, or be completely suppressed altogether. In this paper, we discuss the positive growth rates associated with the instabilities, windows of instability and their dependence on model parameters through extensive numerical computations aided by asymptotic analyses.

Nonequilibrium structure and dynamics in a microscopic model of thin-film active gels

In the presence of adenosine triphosphate, molecular motors generate active force dipoles that drive suspensions of protein filaments far from thermodynamic equilibrium, leading to exotic dynamics and pattern formation. Microscopic modeling can help to quantify the relationship between individual motors plus filaments to organization and dynamics on molecular and supramolecular length scales. Here, we present results of extensive numerical simulations of active gels where the motors and filaments are confined between two infinite parallel plates. Thermal fluctuations and excluded-volume interactions between filaments are included. A systematic variation of rates for motor motion, attachment, and detachment, including a differential detachment rate from filament ends, reveals a range of nonequilibrium behavior. Strong motor binding produces structured filament aggregates that we refer to as asters, bundles, or layers, whose stability depends on motor speed and differential end detachment. The gross features of the dependence of the observed structures on the motor rate and the filament concentration can be captured by a simple one-filament model. Loosely bound aggregates exhibit superdiffusive mass transport, where filament translocation scales with lag time with nonunique exponents that depend on motor kinetics. An empirical data collapse of filament speed as a function of motor speed and end detachment is found, suggesting a dimensional reduction of the relevant parameter space. We conclude by discussing the perspectives of microscopic modeling in the field of active gels.