Diffusive transport without detailed balance in motile bacteria: does microbiology need statistical physics?

Density and polarization of active Brownian particles in curved activity landscapes

Suspensions of motile active particles with space-dependent activity form characteristic polarization and density patterns. Recent single-particle studies for planar activity landscapes identified several quantities associated with emergent density-polarization patterns that are solely determined by bulk variables. Naive thermodynamic intuition suggests that these results might hold for arbitrary activity landscapes mediating bulk regions, and thus could be used as benchmarks for simulations and theories. However, the considered system operates in a nonequilibrium steady state and we prove by construction that the quantities in question lose their simple form for curved activity landscapes. Specifically, we provide a detailed analytical study of polarization and density profiles induced by radially symmetric activity steps, and of the total polarization for the case of a general radially symmetric activity landscape. While the qualitative picture is similar to the planar case, all the investigated variables depend not only on bulk variables but also comprise geometry-induced contributions. We verified that all our analytical results agree with exact numerical calculations.

Condensation transition in the late-time position of a run-and-tumble particle

We study the position distribution P(R[over ⃗],N) of a run-and-tumble particle (RTP) in arbitrary dimension d, after N runs. We assume that the constant speed v>0 of the particle during each running phase is independently drawn from a probability distribution W(v) and that the direction of the particle is chosen isotropically after each tumbling. The position distribution is clearly isotropic, P(R[over ⃗],N)→P(R,N) where R=|R[over ⃗]|. We show that, under certain conditions on d and W(v) and for large N, a condensation transition occurs at some critical value of R=R_{c}∼O(N) located in the large-deviation regime of P(R,N). For R<R_{c} (subcritical fluid phase), all runs are roughly of the same size in a typical trajectory. In contrast, an RTP trajectory with R>R_{c} is typically dominated by a "condensate," i.e., a large single run that subsumes a finite fraction of the total displacement (supercritical condensed phase). Focusing on the family of speed distributions W(v)=α(1-v/v_{0})^{α-1}/v_{0}, parametrized by α>0, we show that, for large N, P(R,N)∼exp[-Nψ_{d,α}(R/N)], and we compute exactly the rate function ψ_{d,α}(z) for any d and α. We show that the transition manifests itself as a singularity of this rate function at R=R_{c} and that its order depends continuously on d and α. We also compute the distribution of the condensate size for R>R_{c}. Finally, we study the model when the total duration T of the RTP, instead of the total number of runs, is fixed. Our analytical predictions are confirmed by numerical simulations, performed using a constrained Markov chain Monte Carlo technique, with precision ∼10^{-100}.

Phase behavior and surface tension of soft active Brownian particles

We study quasi two-dimensional, monodisperse systems of active Brownian particles (ABPs) for a range of activities, stiffnesses, and densities. We develop a microscopic, analytical method for predicting the dense phase structure formed after motility-induced phase separation (MIPS) has occurred, including the dense cluster's area fraction, interparticle pressure, and radius. Our predictions are in good agreement with our Brownian dynamics simulations. We, then, derive a continuum model to investigate the relationship between the predicted interparticle pressure, the swim pressure, and the macroscopic pressure in the momentum equation. We find that formulating the point-wise macroscopic pressure as the interparticle pressure and modeling the particle activity through a spatially variant body force - as opposed to a volume-averaged swim pressure - results in consistent predictions of pressure in both the continuum model and the microscopic theory. This formulation of pressure also results in nearly zero surface tension for the phase separated domains, irrespective of activity, stiffness, and area fraction. Furthermore, using Brownian dynamics simulations and our continuum model, we showed that both the interface width and surface tension, are intrinsic characteristics of the system. On the other hand, if we were to exclude the body force induced by activity, we find that the resulting surface tension values are linearly dependent on the size of the simulation, in contrast to the statistical mechanical definition of surface tension.

Thermotaxis of Janus particles

The interactions of autonomous microswimmers play an important role for the formation of collective states of motile active matter. We study them in detail for the common microswimmer-design of two-faced Janus spheres with hemispheres made from different materials. Their chemical and physical surface properties may be tailored to fine-tune their mutual attractive, repulsive or aligning behavior. To investigate these effects systematically, we monitor the dynamics of a single gold-capped Janus particle in the external temperature field created by an optically heated metal nanoparticle. We quantify the orientation-dependent repulsion and alignment of the Janus particle and explain it in terms of a simple theoretical model for the induced thermoosmotic surface fluxes. The model reveals that the particle's angular velocity is solely determined by the temperature profile on the equator between the Janus particle's hemispheres and their phoretic mobility contrast. The distortion of the external temperature field by their heterogeneous heat conductivity is moreover shown to break the apparent symmetry of the problem.

How Activity Landscapes Polarize Microswimmers without Alignment Forces

Active-particle suspensions exhibit distinct polarization-density patterns in activity landscapes, even without anisotropic particle interactions. Such polarization without alignment forces is at work in motility-induced phase separation and betrays intrinsic microscopic activity to mesoscale observers. Using stable long-term confinement of a single thermophoretic microswimmer in a dedicated force-free particle trap, we examine the polarized interfacial layer at a motility step and confirm that it does not exert pressure onto the bulk. Our observations are quantitatively explained by an analytical theory that can also guide the analysis of more complex geometries and many-body effects.

Polarization-density patterns of active particles in motility gradients

The colocalization of density modulations and particle polarization is a characteristic emergent feature of motile active matter in activity gradients. We employ the active-Brownian-particle model to derive precise analytical expressions for the density and polarization profiles of a single Janus-type swimmer in the vicinity of an abrupt activity step. Our analysis allows for an optional (but not necessary) orientation-dependent propulsion speed, as often employed in force-free particle steering. The results agree well with measurement data for a thermophoretic microswimmer presented in the companion paper [Söker et al., Phys. Rev. Lett. 126, 228001 (2021)10.1103/PhysRevLett.126.228001], and they can serve as a template for more complex applications, e.g., to motility-induced phase separation or studies of physical boundaries. The essential physics behind our formal results is robustly captured and elucidated by a schematic two-species "run-and-tumble" model.

Chemotaxis of Cargo-Carrying Self-Propelled Particles

Active particles with their characteristic feature of self-propulsion are regarded as the simplest models for motility in living systems. The accumulation of active particles in low activity regions has led to the general belief that chemotaxis requires additional features and at least a minimal ability to process information and to control motion. We show that self-propelled particles display chemotaxis and move into regions of higher activity if the particles perform work on passive objects, or cargo, to which they are bound. The origin of this cooperative chemotaxis is the exploration of the activity gradient by the active particle when bound to a load, resulting in an average excess force on the load in the direction of higher activity. Using a new theoretical model, we capture the most relevant features of these active-passive dimers, and in particular we predict the crossover between antichemotactic and chemotactic behavior. Moreover, we show that merely connecting active particles to chains is sufficient to obtain the crossover from antichemotaxis to chemotaxis with increasing chain length. Such an active complex is capable of moving up a gradient of activity such as provided by a gradient of fuel and to accumulate where the fuel concentration is at its maximum. The observed transition is of significance to protoforms of life, enabling them to locate a source of nutrients even in the absence of any supporting sensomotoric apparatus.

Exploring dissipative sources of non-Markovian biochemical reaction systems

Many biological processes including important intracellular processes are governed by biochemical reaction networks. Usually, these reaction systems operate far from thermodynamic equilibrium, implying free-energy dissipation. On the other hand, single reaction events happen often in a memory manner, leading to non-Markovian kinetics. A question then arises: how do we calculate free-energy dissipation (defined as the entropy production rate) in this physically real case? We derive an analytical formula for calculating the energy consumption of a general reaction system with molecular memory characterized by nonexponential waiting-time distributions. It shows that this dissipation is composed of two parts: one from broken detailed balance of an equivalent Markovian system with the same topology and substrates, and the other from the direction-time dependence of waiting-time distributions. But, if the system is in a detailed balance and the waiting-time distribution is direction-time independent, there is no energy dissipation even in the non-Markovian case. These general results provide insights into the physical mechanisms underlying nonequilibrium processes. A continuous-time random-walk model and a generalized model of stochastic gene expression are chosen to clearly show dissipative sources and the relationship between energy dissipation and molecular memory.

Chain stiffness boosts active nanoparticle transport in polymer networks

Recent advances in technologies such as nanomanufacturing and nanorobotics have opened new pathways for the design of active nanoparticles (NPs) capable of penetrating biolayers for biomedical applications, e.g., for drug delivery. The coupling and feedback between active NP motility (with large stochastic increments relative to passive NPs) and the induced nonequilibrium deformation and relaxation responses of the polymer network, spanning scales from the NP to the local structure of the network, remain to be clarified. Using molecular dynamics simulations, combined with a Rouse mode analysis of network chains and position and velocity autocorrelation functions of the NPs, we demonstrate that the mobility of active NPs within cross-linked, concentrated polymer networks is a monotonically increasing function of chain stiffness, contrary to passive NPs, for which chain stiffness suppresses mobility. In flexible networks, active NPs exhibit a behavior similar to passive NPs, with a boost in mobility proportional to the self-propulsion force. These results are suggestive of design strategies for active NP penetration of stiff biopolymer matrices.

Phase separation of an active colloidal suspension via quorum-sensing

We present the Brownian dynamics simulation of an active colloidal suspension in two dimensions, where the self-propulsion speed of a colloid is regulated according to the local density sensed by it. The role of concentration-dependent motility in the phase-separation of colloids and their dynamics is investigated in detail. Interestingly, the system phase separates at a very low packing fraction (Φ≈ 0.125) at higher self-propulsion speeds (Pe), into a dense phase coexisting with a homogeneous phase and attains a long-range crystalline order beyond the transition point. The transition point is quantified here from the local density profiles and local and global-bond order parameters. We have shown that the characteristics of the phase diagram are qualitatively akin to the active Brownian particle (ABP) model. Moreover, our investigation reveals that the density-dependent motility amplifies the slow-down of the directed speed, which facilitates phase-separation even at low packing fractions. The effective diffusivity shows a crossover from quadratic rise to a power-law behavior of exponent 3/2 with Pe in the phase-separated regime. Furthermore, we have shown that the effective diffusion decreases exponentially with packing fraction in the phase-separated regime, while it shows a linear decrease in the single phase regime.

Swarming bacteria undergo localized dynamic phase transition to form stress-induced biofilms

Self-organized multicellular behaviors enable cells to adapt and tolerate stressors to a greater degree than isolated cells. However, whether and how cellular communities alter their collective behaviors adaptively upon exposure to stress is largely unclear. Here, we investigate this question using Bacillus subtilis, a model system for bacterial multicellularity. We discover that, upon exposure to a spatial gradient of kanamycin, swarming bacteria activate matrix genes and transit to biofilms. The initial stage of this transition is underpinned by a stress-induced multilayer formation, emerging from a biophysical mechanism reminiscent of motility-induced phase separation (MIPS). The physical nature of the process suggests that stressors which suppress the expansion of swarms would induce biofilm formation. Indeed, a simple physical barrier also induces a swarm-to-biofilm transition. Based on the gained insight, we propose a strategy of antibiotic treatment to inhibit the transition from swarms to biofilms by targeting the localized phase transition.

Migration and accumulation of bacteria with chemotaxis and chemokinesis

Bacteria can chemotactically migrate up attractant gradients by controlling run-and-tumble motility patterns. In addition to this well-known chemotactic behaviour, several soil and marine bacterial species perform chemokinesis; they adjust their swimming speed according to the local concentration of chemoeffector, with higher speed at higher concentration. A field of attractant then induces a spatially varying swimming speed, which results in a drift towards lower attractant concentrations-contrary to the drift created by chemotaxis. Here, to explore the biological benefits of chemokinesis and investigate its impact on the chemotactic response, we extend a Keller-Segel-type model to include chemokinesis. We apply the model to predict the dynamics of bacterial populations capable of chemokinesis and chemotaxis in chemoeffector fields inspired by microfluidic and agar plate migration assays. We find that chemokinesis combined with chemotaxis not only may enhance the population response with respect to pure chemotaxis, but also modifies it qualitatively. We conclude presenting predictions for bacteria around dynamic finite-size nutrient sources, simulating, e.g. a marine particle or a root. We show that chemokinesis can reduce the measuring bias that is created by a decaying attractant gradient.

Emergent Field-Driven Robot Swarm States

We present an ecology-inspired form of active matter consisting of a robot swarm. Each robot moves over a planar dynamic resource environment represented by a large light-emitting diode array in search of maximum light intensity; the robots deplete (dim) locally by their presence the local light intensity and seek maximum light intensity. Their movement is directed along the steepest local light intensity gradient; we call this emergent symmetry breaking motion "field drive." We show there emerge dynamic and spatial transitions similar to gas, crystalline, liquid, glass, and jammed states as a function of robot density, resource consumption rates, and resource recovery rates. Paradoxically the nongas states emerge from smooth, flat resource landscapes, not rough ones, and each state can directly move to a glassy state if the resource recovery rate is slow enough, at any robot density.

Statistical mechanics of active Ornstein-Uhlenbeck particles

We study the statistical properties of active Ornstein-Uhlenbeck particles (AOUPs). In this simplest of models, the Gaussian white noise of overdamped Brownian colloids is replaced by a Gaussian colored noise. This suffices to grant this system the hallmark properties of active matter, while still allowing for analytical progress. We study in detail the steady-state distribution of AOUPs in the small persistence time limit and for spatially varying activity. At the collective level, we show AOUPs to experience motility-induced phase separation both in the presence of pairwise forces or due to quorum-sensing interactions. We characterize both the instability mechanism leading to phase separation and the resulting phase coexistence. We probe how, in the stationary state, AOUPs depart from their thermal equilibrium limit by investigating the emergence of ratchet currents and entropy production. In the small persistence time limit, we show how fluctuation-dissipation relations are recovered. Finally, we discuss how the emerging properties of AOUPs can be characterized from the dynamics of their collective modes.

Stochastic resetting of active Brownian particles with Lorentz force

Equilibrium properties of a system of passive diffusing particles in an external magnetic field are unaffected by Lorentz force. In contrast, active Brownian particles exhibit steady-state phenomena that depend on both the strength and the polarity of the applied magnetic field. The intriguing effects of the Lorentz force, however, can only be observed when out-of-equilibrium density gradients are maintained in the system. To this end, we use the method of stochastic resetting on active Brownian particles in two dimensions by resetting them to the line x = 0 at a constant rate and periodicity in the y direction. Under stochastic resetting, an active system settles into a nontrivial stationary state which is characterized by an inhomogeneous density distribution, polarization and bulk fluxes perpendicular to the density gradients. We show that whereas for a uniform magnetic field the properties of the stationary state of the active system can be obtained from its passive counterpart, novel features emerge in the case of an inhomogeneous magnetic field which have no counterpart in passive systems. In particular, there exists an activity-dependent threshold rate such that for smaller resetting rates, the density distribution of active particles becomes non-monotonic. We also study the mean first-passage time to the x axis and find a surprising result: it takes an active particle more time to reach the target from any given point for the case when the magnetic field increases away from the axis. The theoretical predictions are validated using Brownian dynamics simulations.

Underdamped active Brownian heat engine

Active Brownian engines rectify energy from reservoirs composed of self-propelling nonequilibrium molecules into work. We consider a class of such engines based on an underdamped Brownian particle trapped in a power-law potential. The energy they transform has thermodynamic properties of heat only if the nonequilibrium reservoir can be assigned a suitable effective temperature consistent with the second law and thus yielding an upper bound on the engine efficiency. The effective temperature exists if the total force exerted on the particle by the bath is not correlated with the particle position. In general, this occurs if the noise autocorrelation function and the friction kernel are proportional as in the fluctuation-dissipation theorem. But even if the proportionality is broken, the effective temperature can be defined in restricted, fine-tuned, parameter regimes, as we demonstrate on a specific example with harmonic potential.

Hard core run and tumble particles on a one-dimensional lattice

We study the large scale behavior of a collection of hard core run and tumble particles on a one-dimensional lattice with periodic boundary conditions. Each particle has persistent motion in one direction decided by an associated spin variable until the direction of spin is reversed. We map the run and tumble model to a mass transfer model with fluctuating directed bonds. We calculate the steady-state single-site mass distribution in the mass model within a mean field approximation for larger spin-flip rates and by analyzing an appropriate coalescence-fragmentation model for small spin-flip rates. We also calculate the hydrodynamic coefficients of diffusivity and conductivity for both large and small spin-flip rates and show that the Einstein relation is violated in both regimes. We also show how the nongradient nature of the process can be taken into account in a systematic manner to calculate the hydrodynamic coefficients.

Periodic oscillations in a string of camphor infused disks

The rhythmic beating motion of autonomously motile filaments has many practical applications. Here, we present an experimental study on a filament made of camphor infused paper disks, stitched together adjacent to each other using nylon thread. The filament displays spontaneous translatory motion when it is placed on the surface of water due to the surface tension gradients created by camphor molecules on the water surface. When this filament is clamped on one end, we obtain regular oscillatory motion instead of translation. The filament shows qualitatively different dynamics at different activity levels, which is controlled by the amount of camphor infused into the paper disks. For a better physical understanding of the filament dynamics, we develop a minimal numerical model involving a semi-flexible filament made of active polar disks, where the polarity is coupled to the instantaneous velocity of the particle. This model qualitatively reproduces different oscillatory modes of the filament. Moreover, our model reveals a rich dynamical state diagram of the system, as a function of filament activity and the coupling strength.

Lamellar to Micellar Phases and Beyond: When Tactic Active Systems Admit Free Energy Functionals

We consider microscopic models of active particles whose velocities, rotational diffusivities, and tumbling rates depend on the gradient of a local field that is either externally imposed or depends on all particle positions. Despite the fundamental differences between active and passive dynamics at the microscopic scale, we show that a large class of such tactic active systems admit fluctuating hydrodynamics equivalent to those of interacting Brownian colloids in equilibrium. We exploit this mapping to show how taxis may lead to the lamellar and micellar phases observed for soft repulsive colloids. In the context of chemotaxis, we show how the competition between chemoattractant and chemorepellent may lead to a bona fide equilibrium liquid-gas phase separation in which a loss of thermodynamic stability of the fluid signals the onset of a chemotactic collapse.

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.

Isomorph theory beyond thermal equilibrium

This paper generalizes isomorph theory to systems that are not in thermal equilibrium. The systems are assumed to be R-simple, i.e., to have a potential energy that as a function of all particle coordinates R obeys the hidden-scale-invariance condition U(R_a) < U(R_b) ⇒ U(λR_a) < U(λR_b). "Systemic isomorphs" are introduced as lines of constant excess entropy in the phase diagram defined by density and systemic temperature, which is the temperature of the equilibrium state point with the average potential energy equal to U(R). The dynamics is invariant along a systemic isomorph if there is a constant ratio between the systemic and the bath temperature. In thermal equilibrium, the systemic temperature is equal to the bath temperature and the original isomorph formalism is recovered. The new approach rationalizes within a consistent framework previously published observations of isomorph invariance in simulations involving nonlinear steady-state shear flows, zero-temperature plastic flows, and glass-state isomorphs. This paper relates briefly to granular media, physical aging, and active matter. Finally, we discuss the possibility that the energy unit defining the reduced quantities should be based on the systemic rather than the bath temperature.

Universal properties of a run-and-tumble particle in arbitrary dimension

We consider an active run-and-tumble particle (RTP) in d dimensions, starting from the origin and evolving over a time interval [0,t]. We examine three different models for the dynamics of the RTP: the standard RTP model with instantaneous tumblings, a variant with instantaneous runs and a general model in which both the tumblings and the runs are noninstantaneous. For each of these models, we use the Sparre Andersen theorem for discrete-time random walks to compute exactly the probability that the x component does not change sign up to time t, showing that it does not depend on d. As a consequence of this result, we compute exactly other x-component properties, namely, the distribution of the time of the maximum and the record statistics, showing that they are universal, i.e., they do not depend on d. Moreover, we show that these universal results hold also if the speed v of the particle after each tumbling is random, drawn from a generic probability distribution. Our findings are confirmed by numerical simulations. Some of these results have been announced in a recent Letter [Phys. Rev. Lett. 124, 090603 (2020)10.1103/PhysRevLett.124.090603].

Single active particle engine utilizing a nonreciprocal coupling between particle position and self-propulsion

We recently argued that a self-propelled particle is formally equivalent to a system consisting of two subsystems coupled by a nonreciprocal interaction [Phys. Rev. E 100, 050603(R) (2019)2470-004510.1103/PhysRevE.100.050603]. Here, we show that this nonreciprocal coupling allows us to extract useful work from a single self-propelled particle maintained at constant temperature, by using an aligning interaction to control correlations between the particle's position and self-propulsion.

Haloarchaea swim slowly for optimal chemotactic efficiency in low nutrient environments

Archaea have evolved to survive in some of the most extreme environments on earth. Life in extreme, nutrient-poor conditions gives the opportunity to probe fundamental energy limitations on movement and response to stimuli, two essential markers of living systems. Here we use three-dimensional holographic microscopy and computer simulations to reveal that halophilic archaea achieve chemotaxis with power requirements one hundred-fold lower than common eubacterial model systems. Their swimming direction is stabilised by their flagella (archaella), enhancing directional persistence in a manner similar to that displayed by eubacteria, albeit with a different motility apparatus. Our experiments and simulations reveal that the cells are capable of slow but deterministic chemotaxis up a chemical gradient, in a biased random walk at the thermodynamic limit.

Universality in incompressible active fluid: Effect of nonlocal shear stress

Phase transitions in active fluids attracted significant attention within the last decades. Recent results show [L. Chen et al., New J. Phys. 17, 042002 (2015)10.1088/1367-2630/17/4/042002] that an order-disorder phase transition in incompressible active fluids belongs to a new universality class. In this work, we further investigate this type of phase transition and focus on the effect of long-range interactions. This is achieved by introducing a nonlocal shear stress into the hydrodynamic description, which leads to superdiffusion of the velocity field, and can be viewed as a result of the active particles performing Lévy walks. The universal properties in the critical region are derived by performing a perturbative renormalization group analysis of the corresponding response functional within the one-loop approximation. We show that the effect of nonlocal shear stress decreases the upper critical dimension of the model, and can lead to the irrelevance of the active fluid self-advection with the resulting model belonging to an unusual long-range Model A universality class not reported before, to our knowledge. Moreover, when the degree of nonlocality is sufficiently high all nonlinearities become irrelevant and the mean-field description is valid in any spatial dimension.

Characterization of active matter in dense suspensions with heterodyne laser Doppler velocimetry

We present a novel approach for characterizing the properties and performance of active matter in dilute suspension as well as in crowded environments. We use Super-Heterodyne Laser-Doppler-Velocimetry (SH-LDV) to study large ensembles of catalytically active Janus particles moving under UV illumination. SH-LDV facilitates a model-free determination of the swimming speed and direction, with excellent ensemble averaging. In addition, we obtain information on the distribution of the catalytic activity. Moreover, SH-LDV operates away from walls and permits a facile correction for multiple scattering contributions. It thus allows for studies of concentrated suspensions of swimmers or of systems where swimmers propel actively in an environment crowded by passive particles. We demonstrate the versatility and the scope of the method with a few selected examples. We anticipate that SH-LDV complements established methods and paves the way for systematic measurements at previously inaccessible boundary conditions.

Testing a thermodynamic approach to collective animal behavior in laboratory fish schools

Collective behaviors displayed by groups of social animals are observed frequently in nature. Understanding and predicting the behavior of complex biological systems is dependent on developing effective descriptions and models. While collective animal systems are characteristically nonequilibrium, we can employ concepts from equilibrium statistical mechanics to motivate the measurement of material-like properties in laboratory animal aggregates. Here, we present results from a new set of experiments that utilize high speed footage of two-dimensional schooling events, particle tracking, and projected static and dynamic light fields to observe and control the behavior of negatively phototaxic fish schools (Hemigrammus bleheri). First, we use static light fields consisting of dark circular regions to produce visual stimuli that confine the schools to a range of areas. We find that schools have a maximum density which is independent of group size, and that a swim pressurelike quantity, Π increases linearly with number density, suggesting that unperturbed schools exist on an isotherm. Next, we use dynamic light fields where the radius of the dark region shrinks linearly with time to compress the schools. We find that an effective temperature parameter depends on the compression time and our results are thus consistent with the school having a constant heat flux. These findings further evidence the utility of effective thermodynamic descriptions of nonequilibrium systems in collective animal behavior.

An attraction-repulsion transition of force on two asymmetric wedges induced by active particles

Effective interaction between two asymmetric wedges immersed in a two-dimensional active bath is investigated by computer simulations. The attraction–repulsion transition of effective force between two asymmetric wedges is subjected to the relative position of two wedges, the wedge-to-wedge distance, the active particle density, as well as the apex angle of two wedges. By exchanging the position of the two asymmetric wedges in an active bath, firstly a simple attraction–repulsion transition of effective force occurs, completely different from passive Brownian particles. Secondly the transition of effective force is symmetric for the long-range distance between two asymmetric wedges, while it is asymmetric for the short-range case. Our investigations may provide new possibilities to govern the motion and assembly of microscopic objects by taking advantage of the self-driven behaviour of active particles.

Information and motility exchange in collectives of active particles

We examine the interplay of motility and information exchange in a model of run-and-tumble active particles where the particle's motility is encoded as a bit of information that can be exchanged upon contact according to the rules of AND and OR logic gates in a circuit. Motile AND particles become non-motile upon contact with a non-motile particle. Conversely, motile OR particles remain motile upon collision with their non-motile counterparts. AND particles that have become non-motile additionally "reawaken", i.e., recover their motility, at a fixed rate μ, as in the SIS (susceptible, infected, susceptible) model of epidemic spreading, where an infected agent can become healthy again, but keeps no memory of the recent infection, hence it is susceptible to a renewed infection. For μ = 0, both AND and OR particles relax irreversibly to absorbing states of all non-motile or all motile particles, respectively. The relaxation kinetics is, however, faster for OR particles that remain active throughout the process. At finite μ, the AND dynamics is controlled by the interplay between reawakening and collision rates. The system evolves to a state of all motile particles (an absorbing state in the language of absorbing phase transitions) for μ > μc and to a mixed state with coexisting motile and non-motile particles (an active state in the language of absorbing phase transitions) for μ < μc. The final state exhibits a rich structure controlled by motility-induced aggregation. Our work can be relevant to biochemical signaling in motile bacteria, the spreading of epidemics and of social consensus, as well as light-controlled organization of active colloids.

Symmetric Mixtures of Pusher and Puller Microswimmers Behave as Noninteracting Suspensions

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

Run-and-tumble particles in two dimensions: Marginal position distributions

We study a set of run-and-tumble particle (RTP) dynamics in two spatial dimensions. In the first case of the orientation θ of the particle can assume a set of n possible discrete values, while in the second case θ is a continuous variable. We calculate exactly the marginal position distributions for n=3,4 and the continuous case and show that in all cases the RTP shows a crossover from a ballistic to diffusive regime. The ballistic regime is a typical signature of the active nature of the systems and is characterized by nontrivial position distributions which depend on the specific model. We also show that the signature of activity at long times can be found in the atypical fluctuations, which we also characterize by computing the large deviation functions explicitly.

Binary mixtures of active and passive particles on a sphere

We study the cooperation and segregation dynamics of binary mixtures of active and passive particles on a sphere. According to the competition between rotational diffusion and polar alignment, we find three distinct phases: a mixed phase and two different demixed phases. When rotational diffusion dominates the dynamics, the demixing is due to the aggregation of passive particles, where active and passive particles respectively occupy two hemispheres. When polar alignment is dominated, the demixing is caused by the aggregation of active particles, where active particles occupy the equator of the sphere and passive particles occupy the two poles of the sphere. In this case, there exist a circulating band cluster and two cambered surface clusters, which is a purely curvature-driven effect with no equivalent in the planar model. When rotational diffusion and polar alignment are comparable, particles are completely mixed. Our findings are relevant to the experimental pursuit of segregation dynamics of binary mixtures on curved surfaces.

Formation of stable and responsive collective states in suspensions of active colloids

Many animal species organise into disordered swarms, polarised flocks or swirls to protect from predators or optimise foraging. Previous studies suggest that such collective states are related to a critical point, which could explain their balance between robustness to noise and high responsiveness regarding external perturbations. Here we provide experimental evidence for this idea by investigating the stability of swirls formed by light-responsive active colloids which adjust their individual motion to positions and orientations of neighbours. Because their behaviour can be precisely tuned, controlled changes between different collective states can be achieved. During the transition between stable swirls and swarms we observe a maximum of the group's susceptibility indicating the vicinity of a critical point. Our results support the idea of system-independent organisation principles of collective states and provide useful strategies for the realisation of responsive yet stable ensembles in microrobotic systems.

Steady-state distributions and nonsteady dynamics in nonequilibrium systems

We search for steady states in a class of fluctuating and driven physical systems that exhibit sustained currents. We find that the physical concept of a steady state, well known for systems at equilibrium, must be generalized to describe such systems. In these, the generalization of a steady state is associated with a stationary probability density of microstates and a deterministic dynamical system whose trajectories the system follows on average. These trajectories are a manifestation of nonstationary macroscopic currents observed in these systems. We determine precise conditions for the steady state to exist as well as the requirements for it to be stable. We illustrate this with some examples.

Oscillating collective motion of active rotors in confinement

Due to its inherent out-of-equilibrium nature, active matter in confinement may exhibit collective behavior absent in unconfined systems. Extensive studies have indicated that hydrodynamic or steric interactions between active particles and boundary play an important role in the emergence of collective behavior. However, besides introducing external couplings at the single-particle level, the confinement also induces an inhomogeneous density distribution due to particle-position correlations, whose effect on collective behavior remains unclear. Here, we investigate this effect in a minimal chiral active matter composed of self-spinning rotors through simulation, experiment, and theory. We find that the density inhomogeneity leads to a position-dependent frictional stress that results from interrotor friction and couples the spin to the translation of the particles, which can then drive a striking spatially oscillating collective motion of the chiral active matter along the confinement boundary. Moreover, depending on the oscillation properties, the collective behavior has three different modes as the packing fraction varies. The structural origins of the transitions between the different modes are well identified by the percolation of solid-like regions or the occurrence of defect-induced particle rearrangement. Our results thus show that the confinement-induced inhomogeneity, dynamic structure, and compressibility have significant influences on collective behavior of active matter and should be properly taken into account.

The 2020 motile active matter roadmap

Activity and autonomous motion are fundamental in living and engineering systems. This has stimulated the new field of 'active matter' in recent years, which focuses on the physical aspects of propulsion mechanisms, and on motility-induced emergent collective behavior of a larger number of identical agents. The scale of agents ranges from nanomotors and microswimmers, to cells, fish, birds, and people. Inspired by biological microswimmers, various designs of autonomous synthetic nano- and micromachines have been proposed. Such machines provide the basis for multifunctional, highly responsive, intelligent (artificial) active materials, which exhibit emergent behavior and the ability to perform tasks in response to external stimuli. A major challenge for understanding and designing active matter is their inherent nonequilibrium nature due to persistent energy consumption, which invalidates equilibrium concepts such as free energy, detailed balance, and time-reversal symmetry. Unraveling, predicting, and controlling the behavior of active matter is a truly interdisciplinary endeavor at the interface of biology, chemistry, ecology, engineering, mathematics, and physics. The vast complexity of phenomena and mechanisms involved in the self-organization and dynamics of motile active matter comprises a major challenge. Hence, to advance, and eventually reach a comprehensive understanding, this important research area requires a concerted, synergetic approach of the various disciplines. The 2020 motile active matter roadmap of Journal of Physics: Condensed Matter addresses the current state of the art of the field and provides guidance for both students as well as established scientists in their efforts to advance this fascinating area.

Pseudo-chemotaxis of active Brownian particles competing for food

Active Brownian particles (ABPs) are physical models for motility in simple life forms and easily studied in simulations. An open question is to what extent an increase of activity by a gradient of fuel, or food in living systems, results in an evolutionary advantage of actively moving systems such as ABPs over non-motile systems, which rely on thermal diffusion only. It is an established fact that within confined systems in a stationary state, the activity of ABPs generates density profiles that are enhanced in regions of low activity, which is thus referred to as ‘anti-chemotaxis’. This would suggest that a rather complex sensoric subsystem and information processing is a precondition to recognize and navigate towards a food source. We demonstrate in this work that in non-stationary setups, for instance as a result of short bursts of fuel/food, ABPs do in fact exhibit chemotactic behavior. In direct competition with inactive, but otherwise identical Brownian particles (BPs), the ABPs are shown to fetch a larger amount of food. We discuss this result based on simple physical arguments. From the biological perspective, the ability of primitive entities to move in direct response to the available amount of external energy would, even in absence of any sensoric devices, encompass an evolutionary advantage.

Inferring Entropy Production from Short Experiments

We provide a strategy for the exact inference of the average as well as the fluctuations of the entropy production in nonequilibrium systems in the steady state, from the measurements of arbitrary current fluctuations. Our results are built upon the finite-time generalization of the thermodynamic uncertainty relation, and require only very short time series data from experiments. We illustrate our results with exact and numerical solutions for two colloidal heat engines.

Universal Survival Probability for a d-Dimensional Run-and-Tumble Particle

We consider an active run-and-tumble particle (RTP) in d dimensions and compute exactly the probability S(t) that the x component of the position of the RTP does not change sign up to time t. When the tumblings occur at a constant rate, we show that S(t) is independent of d for any finite time t (and not just for large t), as a consequence of the celebrated Sparre Andersen theorem for discrete-time random walks in one dimension. Moreover, we show that this universal result holds for a much wider class of RTP models in which the speed v of the particle after each tumbling is random, drawn from an arbitrary probability distribution. We further demonstrate, as a consequence, the universality of the record statistics in the RTP problem.

Is Graphene Oxide a Chemoattractant?

The use of nanomaterials to regulate cell surface receptors is considered a novel strategy to manipulate cell behaviors. However, recognition is important to drive nanoparticle-cell complex formation. Here, we report a novel approach that uses graphene oxide (GO) as a chemoattractant to lure bacteria to nanosurface, facilitating complex formation. The amount of Escherichia coli (E. coli) cells attracted into capillaries containing 20 mg/L GO was more than 8.6-fold higher than that attracted into capillaries containing 20 mg/L glucose. The inherent mechanism involved interference with transmembrane chemoreceptors and activation of the chemotactic system via GO attachment and a subsequent increase in cell aggregation and migration via self-secreted quorum sensing molecules. The key feature of this strategy is the potential to improve the efficiency of the nanoparticle-cell recognition pattern and to expedite the development of surface-contact-related nanotechnology.

A combined rheometry and imaging study of viscosity reduction in bacterial suspensions

Suspending self-propelled "pushers" in a liquid lowers its viscosity. We study how this phenomenon depends on system size in bacterial suspensions using bulk rheometry and particle-tracking rheoimaging. Above the critical bacterial volume fraction needed to decrease the viscosity to zero, [Formula: see text], large-scale collective motion emerges in the quiescent state, and the flow becomes nonlinear. We confirm a theoretical prediction that such instability should be suppressed by confinement. Our results also show that a recent application of active liquid-crystal theory to such systems is untenable.

Nonquenched rotators ease flocking and memorize it

We introduce a minimal model for a two-dimensional polar flock with nonquenched rotators and show that the rotators make the usual macroscopic long-range order of the flock more robust than the clean system. The rotators memorize the flock-information which helps in establishing the robustness. Moreover, the memory of the rotators assists in probing the moving flock. We also formulate a hydrodynamic framework for the microscopic model that makes our study comprehensive. Using linearized hydrodynamics, it is shown that the presence of such nonquenched heterogeneities increases the sound speeds of the flock. The enhanced sound speeds lead to faster convection of information and consequently the robust ordering in the system. We argue that similar nonquenched heterogeneities may be useful in monitoring and controlling large crowds.

Run-and-tumble motion in one dimension with space-dependent speed

We consider a particle performing run-and-tumble dynamics with space-dependent speed. The model has biological relevance as it describes motile bacteria or cells in heterogeneous environments. We give exact expression for the probability density function in the case of free motion in unbounded space. We then analyze the case of a particle moving in a confined interval in the presence of partially absorbing boundaries, reporting the probability density in the Laplace (time) domain and the mean time to absorption. We also discuss the relaxation to the steady state in the case of confinement with reflecting boundaries and drift effects due to direction-dependent tumbling rates, modeling taxis phenomena of cells. The case of diffusive particles with spatially variable diffusivity is obtained as a limiting case.

Rotating oil droplets driven by motile bacteria at interfaces

We show that oil droplets suspended near a liquid-solid interface can be driven to rotate by motile bacteria adhered to the droplet surface. Droplets rotate clockwise when viewed from the liquid side, due to symmetry-breaking hydrodynamic interactions of bacteria with the interface. The angular speed of rotation for droplets decreases as their size is increased. Differences in the speed of rotation driven by Escherichia coli, Shewanella haliotis, and Halomonas titanicae bacteria reflects differences in the number of bacteria adhered at the droplet surface and their interfacial affinity. Adding surfactant reduces the number of adherent bacteria and hence lowers the speed of rotation. Together, these results demonstrate that bacterial activity can be used to manipulate suspended droplets.

Linear Response Theory and Green-Kubo Relations for Active Matter

We address the question of how interacting active systems in a nonequilibrium steady state respond to an external perturbation. We establish an extended fluctuation-dissipation theorem for active Brownian particles (ABP), which highlights the role played by the local violation of detailed balance due to activity. By making use of a Markovian approximation we derive closed Green-Kubo expressions for the diffusivity and mobility of ABP and quantify the deviations from the Stokes-Einstein relation. We compute the linear response function to an external force using unperturbed simulations of ABP and compare the results with the analytical predictions of the transport coefficients. Our results show the importance of the interplay between activity and interactions in the departure from equilibrium linear response.

Long-time position distribution of an active Brownian particle in two dimensions

We study the late-time dynamics of a single active Brownian particle in two dimensions with speed v_{0} and rotation diffusion constant D_{R}. We show that at late times t≫D_{R}^{-1}, while the position probability distribution P(x,y,t) in the x-y plane approaches a Gaussian form near its peak describing the typical diffusive fluctuations, it has non-Gaussian tails describing atypical rare fluctuations when sqrt[x^{2}+y^{2}]∼v_{0}t. In this regime, the distribution admits a large deviation form, P(x,y,t)∼exp{-tD_{R}Φ[sqrt[x^{2}+y^{2}]/(v_{0}t)]}, where we compute the rate function Φ(z) analytically and also numerically using an importance sampling method. We show that the rate function Φ(z), encoding the rare fluctuations, still carries the trace of activity even at late times. Another way of detecting activity at late times is to subject the active particle to an external harmonic potential. In this case we show that the stationary distribution P_{stat}(x,y) depends explicitly on the activity parameter D_{R}^{-1} and undergoes a crossover, as D_{R} increases, from a ring shape in the strongly active limit (D_{R}→0) to a Gaussian shape in the strongly passive limit (D_{R}→∞).

Stochastic thermodynamics for self-propelled particles

We propose a generalization of stochastic thermodynamics to systems of active particles, which move under the combined influence of stochastic internal self-propulsions (activity) and a heat bath. The main idea is to consider joint trajectories of particles' positions and self-propulsions. It is then possible to exploit formal similarity of an active system and a system consisting of two subsystems interacting with different heat reservoirs and coupled by a nonsymmetric interaction. The resulting thermodynamic description closely follows the standard stochastic thermodynamics. In particular, total entropy production, Δs_{tot}, can be decomposed into housekeeping, Δs_{hk}, and excess, Δs_{ex}, parts. Both Δs_{tot} and Δs_{hk} satisfy fluctuation theorems. The average rate of the steady-state housekeeping entropy production can be related to the violation of the fluctuation-dissipation theorem via a Harada-Sasa relation. The excess entropy production enters into a Hatano-Sasa-like relation, which leads to a generalized Clausius inequality involving the change of the system's entropy and the excess entropy production. Interestingly, although the evolution of particles' self-propulsions is free and uncoupled from that of their positions, nontrivial steady-state correlations between these variables lead to the nonzero excess dissipation in the reservoir coupled to the self-propulsions.

Phase coexistence of active Brownian particles

We investigate motility-induced phase separation of active Brownian particles, which are modeled as purely repulsive spheres that move due to a constant swim force with freely diffusing orientation. We develop on the basis of power functional concepts an analytical theory for nonequilibrium phase coexistence and interfacial structure. Theoretical predictions are validated against Brownian dynamics computer simulations. We show that the internal one-body force field has four nonequilibrium contributions: (i) isotropic drag and (ii) interfacial drag forces against the forward motion, (iii) a superadiabatic spherical pressure gradient, and (iv) the quiet life gradient force. The intrinsic spherical pressure is balanced by the swim pressure, which arises from the polarization of the free interface. The quiet life force opposes the adiabatic force, which is due to the inhomogeneous density distribution. The balance of quiet life and adiabatic forces determines bulk coexistence via equality of two bulk state functions, which are independent of interfacial contributions. The internal force fields are kinematic functionals which depend on density and current but are independent of external and swim forces, consistent with power functional theory. The phase transition originates from nonequilibrium repulsion, with the agile gas being more repulsive than the quiet liquid.