Dynamics and density distribution of strongly confined noninteracting nonaligning self-propelled particles in a nonconvex boundary

Kinetically arrested periodic clusters in active filament arrays

We study the dynamics and pattern formation of ordered arrays of active semi-flexible filaments, each of which is pinned at one end and free at the other. The filaments are modeled as connected chains of polar active particles with activity incorporated through local follower forces acting along their local tangent. Using Brownian dynamics simulations in two dimensions, we show that for a range of activity and filament separation, the filament array self-assembles into regularly spaced, kinetically arrested compact clusters. Activity, array geometry, filament elasticity, and grafting density are each seen to crucially influence the show size, shape, and spacing of these emergent clusters. Furthermore, cluster shapes for different grafting densities can be rescaled into self-similar forms with activity-dependent scaling exponents. We derive theoretical expressions that relate the number of filaments in a cluster and the spacing between adjacent clusters to filament activity, filament elasticity, and grafting density. Our results provide insight into the physical mechanisms involved in the initiation of clustering and suggest that steric contact forces and friction balance active forces and filament elasticity to shape and stabilize emergence clusters.

Achieving designed texture and flows in bulk active nematics using optimal control theory

Being intrinsically nonequilibrium, active materials can potentially perform functions that would be thermodynamically forbidden in passive materials. However, active systems have diverse local attractors that correspond to distinct dynamical states, many of which exhibit chaotic turbulent-like dynamics and thus cannot perform work or useful functions. Designing such a system to choose a specific dynamical state is a formidable challenge. Motivated by recent advances enabling optogenetic control of experimental active materials, we describe an optimal control theory framework that identifies a spatiotemporal sequence of light-generated activity that drives an active nematic system toward a prescribed dynamical steady state. Active nematics are unstable to spontaneous defect proliferation and chaotic streaming dynamics in the absence of control. We demonstrate that optimal control theory can compute activity fields that redirect the dynamics into a variety of alternative dynamical programs and functions. This includes dynamically reconfiguring between states, selecting and stabilizing emergent behaviors that do not correspond to attractors, and are hence unstable in the uncontrolled system. Our results provide a roadmap to leverage optical control methods to rationally design structure, dynamics, and function in a wide variety of active materials.

Strong Casimir-like Forces in Flocking Active Matter

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

Motility-Induced Clustering of Active Particles under Soft Confinement

We investigate the structural and dynamic properties of active Brownian particles (APs) confined within a soft annulus-shaped channel. Depending on the strength of the confinement and the Péclet number, we observe a novel reentrant behavior that is not present in unconfined systems. Our findings are substantiated by numerical simulations and analytical considerations, revealing that this behavior arises from the strong coupling between the Péclet number and the effective confining dimensionality of the APs. Our work highlights the peculiarities of soft boundaries for APs and how clogging can be avoided under such conditions.

Motion of microswimmers in cylindrical microchannels

Biological and artificial microswimmers often have to propel through a variety of environments, ranging from heterogeneous suspending media to strong geometrical confinement. Under confinement, local flow fields generated by microswimmers, and steric and hydrodynamic interactions with their environment determine the locomotion. We propose a squirmer-like model to describe the motion of microswimmers in cylindrical microchannels, where propulsion is generated by a fixed surface slip velocity. The model is studied using an approximate analytical solution for cylindrical swimmer shapes, and by numerical hydrodynamics simulations for spherical and spheroidal shapes. For the numerical simulations, we employ the dissipative particle dynamics method for modelling fluid flow. Both the analytical model and simulations show that the propulsion force increases with increasing confinement. However, the swimming velocity under confinement remains lower than the swimmer speed without confinement for all investigated conditions. In simulations, different swimming modes (i.e. pusher, neutral, puller) are investigated, and found to play a significant role in the generation of propulsion force when a swimmer approaches a dead end of a capillary tube. Propulsion generation in confined systems is local, such that the generated flow field generally vanishes beyond the characteristic size of the swimmer. These results contribute to a better understanding of microswimmer force generation and propulsion under strong confinement, including the motion in porous media and in narrow channels.

Inertia suppresses signatures of activity of active Brownian particles in a harmonic potential

A harmonically trapped active Brownian particle exhibits two types of positional distributions-one has a single peak and the other has a single well-that signify steady-state dynamics with low and high activity, respectively. Adding inertia to the translational motion preserves this strict classification of either single-peak or single-well densities but shifts the dividing boundary between the states in the parameter space. We characterize this shift for the dynamics in one spatial dimension using the static Fokker-Planck equation for the full joint distribution of the state space. We derive local results analytically with a perturbation method for a small rotational velocity and then extend them globally with a numerical approach.

Chiral active particles are sensitive reporters to environmental geometry

Chiral active particles (CAPs) are self-propelling particles that break time-reversal symmetry by orbiting or spinning, leading to intriguing behaviors. Here, we examined the dynamics of CAPs moving in 2D lattices of disk obstacles through active Brownian dynamics simulations and granular experiments with grass seeds. We find that the effective diffusivity of the CAPs is sensitive to the structure of the obstacle lattice, a feature absent in achiral active particles. We further studied the transport of CAPs in obstacle arrays under an external field and found a reentrant directional locking effect, which can be used to sort CAPs with different activities. Finally, we demonstrated that parallelogram lattices of obstacles without mirror symmetry can separate clockwise and counter-clockwise CAPs. The mechanisms of the above three novel phenomena are qualitatively explained. As such, our work provides a basis for designing chirality-based tools for single-cell diagnosis and separation, and active particle-based environmental sensors.

Passive objects in confined active fluids: A localization transition

We study how walls confining active fluids interact with asymmetric passive objects placed in their bulk. We show that the objects experience nonconservative long-ranged forces mediated by the active bath. To leading order, these forces can be computed using a generalized image theorem. The walls repel asymmetric objects, irrespective of their microscopic properties or their orientations. For circular cavities, we demonstrate how this may lead to the localization of asymmetric objects in the center of the cavity, something impossible for symmetric ones.

Computational exploration of treadmilling and protrusion growth observed in fire ant rafts

Collective living systems regularly achieve cooperative emergent functions that individual organisms could not accomplish alone. The rafts of fire ants (Solenopsis invicta) are often studied in this context for their ability to create aggregated structures comprised entirely of their own bodies, including tether-like protrusions that facilitate exploration of and escape from flooded environments. While similar protrusions are observed in cytoskeletons and cellular aggregates, they are generally dependent on morphogens or external gradients leaving the isolated role of local interactions poorly understood. Here we demonstrate through an ant-inspired, agent-based numerical model how protrusions in ant rafts may emerge spontaneously due to local interactions. The model is comprised of a condensed structural network of agents that represents the monolayer of interconnected worker ants, which floats on the water and gives ant rafts their form. Experimentally, this layer perpetually contracts, which we capture through the pairwise contraction of all neighboring structural agents at a strain rate of [Formula: see text]. On top of the structural layer, we model a dispersed, on-lattice layer of motile agents that represents free ants, which walk on top of the floating network. Experimentally, these self-propelled free ants walk with some mean persistence length and speed that we capture through an ant-inspired phenomenological model. Local interactions occur between neighboring free ants within some radius of detection, R, and the persistence length of freely active agents is tuned through a noise parameter, η as introduced by the Vicsek model. Both R and η where fixed to match the experimental trajectories of free ants. Treadmilling of the raft occurs as agents transition between the structural and free layers in accordance with experimental observations. Ultimately, we demonstrate how phases of exploratory protrusion growth may be induced by increased ant activity as characterized by a dimensionless parameter, [Formula: see text]. These results provide an example in which functional morphogenesis of a living system may emerge purely from local interactions at the constituent length scale, thereby providing a source of inspiration for the development of decentralized, autonomous active matter and swarm robotics.

Absorption of self-propelled particles into a dense porous medium

We study the absorption of self-propelled particles into a finite-size dense porous medium, which is mimicked by an obstacle array. We find that, depending on the competition of the propelling strength versus the repulsive barrier formed by obstacles and the contrast between the characteristic time scales of permeation and propelling persistence, the absorption process exhibits three distinct types of behavior. In Type I and II behavior, the propelling strength is not large enough to surmount the barrier, and hence particles transport in the medium by barrier-hopping dynamics. The initial permeation of particles toward the medium center is phenomenologically similar to a normal slow diffusion process. But, surprisingly, after the initial permeation process, a concentrated nucleus of particle aggregates forms and grows at the medium center in Type I, due to the long propelling persistence. Such an abnormal "nucleation" phenomenon does not appear in Type II, in which the propelling persistence is low. When the propelling strength is very high (Type III), particles transport smoothly in the medium, hence the initial slow diffusion process disappears and small particle clusters form and merge randomly in the medium. Our results provide a foundation for applications of active objects in a complex environment and also suggest the possible usage of a porous medium, for example, in the selection or sorting of active matter.

Rotation reversal of a ratchet gear powered by active particles

Rotation of a gear powered by active particles is numerically investigated in a circular chamber. Due to the nonequilibrium properties of active particles, net gear rotation is achieved in a bath composed of self-propelling particles. Our setup can convert the random motion of active particles into the directional rotation of the ratchet gear. The direction of rotation is determined by the asymmetry of the gear and the persistence length (the ratio of the self-propulsion speed to the rotation diffusion coefficient) of active particles. Remarkably, the direction of rotation for large persistence length is opposite to the direction of rotation for small persistence length. Therefore, for a given asymmetric gear, we can observe the rotation reversal when tuning the system parameters (e.g., the self-propulsion speed, the rotation diffusion coefficient, and the packing fraction of active particles). Our findings are relevant to the experimental pursuit of rectifying random motion to directional motion in active matter.

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.

Treadmilling and dynamic protrusions in fire ant rafts

Fire ants (Solenopsis invicta) are exemplary for their formation of cohered, buoyant and dynamic structures composed entirely of their own bodies when exposed to flooded environments. Here, we observe tether-like protrusions that emerge from aggregated fire ant rafts when docked to stationary, vertical rods. Ant rafts comprise a floating, structural network of interconnected ants on which a layer of freely active ants walk. We show here that sustained shape evolution is permitted by the competing mechanisms of perpetual raft contraction aided by the transition of bulk structural ants to the free active layer and outward raft expansion owing to the deposition of free ants into the structural network at the edges, culminating in global treadmilling. Furthermore, we see that protrusions emerge as a result of asymmetries in the edge deposition rate of free ants. Employing both experimental characterization and a model for self-propelled particles in strong confinement, we interpret that these asymmetries are likely to occur stochastically owing to wall accumulation effects and directional motion of active ants when strongly confined by the protrusions' relatively narrow boundaries. Together, these effects may realize the cooperative, yet spontaneous formation of protrusions that fire ants sometimes use for functional exploration and to escape flooded environments.

Active particles with polar alignment in ring-shaped confinement

We study steady-state properties of active, nonchiral and chiral Brownian particles with polar alignment and steric interactions confined within a ring-shaped confinement (annulus) in two dimensions. Exploring possible interplays between polar interparticle alignment, geometric confinement and the surface curvature, being incorporated here on minimal levels, we report a surface-population reversal effect, whereby active particles migrate from the outer concave boundary of the annulus to accumulate on its inner convex boundary. This contrasts the conventional picture, implying stronger accumulation of active particles on concave boundaries relative to the convex ones. The population reversal is caused by both particle alignment and surface curvature, disappearing when either of these factors is absent. We explore the ensuing consequences for the chirality-induced current and swim pressure of active particles and analyze possible roles of system parameters, such as the mean number density of particles and particle self-propulsion, chirality, and alignment strengths.

Confinement-induced alternating interactions between inclusions in an active fluid

In a system of colloidal inclusions suspended in an equilibrium bath of smaller particles, the particulate bath engenders effective, short-ranged, primarily attractive interactions between the inclusions, known as depletion interactions, that originate from the steric depletion of bath particles from the immediate vicinity of the inclusions. In a bath of active (self-propelled) particles, the nature of such bath-mediated interactions can qualitatively change from attraction to repulsion, and they become stronger in magnitude and range of action as compared with typical equilibrium depletion interactions, especially as the bath activity (particle self-propulsion) is increased. We study effective interactions mediated by a bath of active Brownian particles between two fixed, impenetrable, and disk-shaped inclusions in a planar (channel) confinement in two dimensions. Confinement is found to strongly influence the effective interaction between the inclusions, specifically by producing alternating interaction profiles with possible attractive and repulsive regions in sufficiently narrow channels. We study the dependence of the ensuing interactions on different system parameters and the orientational (parallel versus perpendicular) configuration of the inclusion pair relative to the channel walls. The confinement effects are largely regulated by the layering of active particles next to the surface boundaries, both of the inclusions and the channel walls that counteract one another in accumulating the active particles in their own proximities. In narrow channels, the combined effects due to the channel walls and the inclusions lead to peculiar structuring of active particles (reminiscent of wavelike interference patterns) within the channel.

Transport of self-propelled particles across a porous medium: trapping, clogging, and the Matthew effect

We study the transport of self-propelled particles from one free chamber to another across two stripe-like areas of dense porous medium. The medium is mimicked by arrays of obstacles. We find that active motion could greatly speed up the transport of particles. However, more and more particles become trapped in the obstacle arrays with the enhancement of activity. At high persistence (low rotational diffusion rate) and moderate particle concentration, we observe the Matthew effect in the aggregation of particles in the two obstacle arrays. This effect is weakened by introduction of randomness or deformability into the obstacle arrays. Moreover, the dependence on deformability shows the characteristics of first-order phase transition. In rare situations, the system could be stuck in a dynamic unstable state, e.g. the particles alternatively gather more in one of the two obstacle arrays, exhibiting oscillation of particle number between the arrays. Our results reveal new features in the transport of active objects in a complex medium and have implications for manipulating their collective assembly.

Vortex formation of spherical self-propelled particles around a circular obstacle

A vortex is a common ratchet phenomenon in active systems. The spatial symmetry is usually broken by introducing asymmetric shapes or spontaneously by collective motion in the presence of hydrodynamic interactions or other alignment effects. Unexpectedly, we observe, by simulations, the formation of a vortex in the simplest model of a circular obstacle immersed in a bath of spherical self-propelled particles. No symmetry-breaking factors mentioned above are included in this model. The vortex forms only when the particle activity is high, i.e. large persistence. The obstacle size is also a key factor and the vortex only forms in a limited range of obstacle sizes. The sustainment of the vortex originates from the bias of the rotating particle cluster around the obstacle in accepting the incoming particles based on their propelling directions. Our results provide new understanding of and insights into the spontaneous symmetry-breaking and ratchet phenomena in active matter.

Rectification of twitching bacteria through narrow channels: A numerical simulations study

Bacteria living on surfaces use different types of motility mechanisms to move on the surface in search of food or to form microcolonies. Twitching is one such form of motility employed by bacteria such as Neisseria gonorrhoeae, in which the polymeric extensions known as type IV pili mediate its movement. Pili extending from the cell body adhere to the surface and pull the bacteria by retraction. The bacterial movement is decided by the two-dimensional tug-of-war among the pili attached to the surface. Natural surfaces on which these microcrawlers dwell are generally spatially inhomogeneous and have varying surface properties. Their motility is known to be affected by the topography of the surfaces. Therefore, it is possible to control bacterial movement by designing structured surfaces which can be potentially utilized for controlling biofilm architecture. In this paper, we numerically investigate the twitching motility in a two-dimensional corrugated channel. The bacterial movement is simulated by two different models: (a) a detailed tug-of-war model which extensively describe the twitching motility of bacteria assisted by pili and (b) a coarse-grained run-and-tumble model which depicts the motion of wide-ranging self-propelled particles. The simulation of bacterial motion through asymmetric corrugated channels using the above models show rectification. The bacterial transport depends on the geometric parameters of the channel and inherent system parameters such as persistence length and self-propelled velocity. In particular, the variation of the particle current with the geometric parameters of the microchannels shows that one can optimize the particle current for specific values of these parameters.

Topotaxis of active Brownian particles

Recent experimental studies have demonstrated that cellular motion can be directed by topographical gradients, such as those resulting from spatial variations in the features of a micropatterned substrate. This phenomenon, known as topotaxis, has been extensively studied for topographical gradients at the subcellular scale, but can also occur in the presence of a spatially varying density of cell-sized features. Such a large-scale topotaxis has recently been observed in highly motile cells that persistently crawl within an array of obstacles with smoothly varying lattice spacing. We introduce a toy model of large-scale topotaxis, based on active Brownian particles. Using numerical simulations and analytical arguments, we demonstrate that topographical gradients introduce a spatial modulation of the particles' persistence, leading to directed motion toward regions of higher persistence. Our results demonstrate that persistent motion alone is sufficient to drive large-scale topotaxis and could serve as a starting point for more detailed studies on self-propelled particles and cells.

Response of active Brownian particles to boundary driving

We computationally study the behavior of underdamped active Brownian particles in a sheared channel geometry. Due to their underdamped dynamics, the particles carry momentum a characteristic distance away from the boundary before it is dissipated into the substrate. We correlate this distance with the persistence of particle trajectories, determined jointly by their friction and self-propulsion. Within this characteristic length, we observe counterintuitive phenomena stemming from the interplay of activity, interparticle interactions, and the boundary driving. Depending on the values of friction and self-propulsion, interparticle interactions can either aid or hinder momentum transport. More dramatically, in certain cases we observe a flow reversal near the wall, which we correlate with an induced polarization of the particle self-propulsion directions. We rationalize these results in terms of a simple kinetic picture of particle trajectories.

Self-propelled particle in a nonconvex external potential: Persistent limit in one dimension

Equilibrium mapping techniques for nonaligning self-propelled particles have made it possible to predict the density profile of an active ideal gas in a wide variety of external potentials. However, they fail when the self-propulsion is very persistent and the potential is nonconvex, which is precisely when the most uniquely active phenomena occur. Here, we show how to predict the density profile of a 1D active Ornstein-Uhlenbeck particle in an arbitrary external potential in the persistent limit and discuss the consequences of the potential's nonconvexity on the structure of the solution, including the central role of the potential's inflection points and the nonlocal dependence of the density profile on the potential.

Geometric control of bacterial surface accumulation

Controlling and suppressing bacterial accumulation at solid surfaces is essential for preventing biofilm formation and biofouling. Whereas various chemical surface treatments are known to reduce cell accumulation and attachment, the role of complex surface geometries remains less well understood. Here, we report experiments and simulations that explore the effects of locally varying boundary curvature on the scattering and accumulation dynamics of swimming Escherichia coli bacteria in quasi-two-dimensional microfluidic channels. Our experimental and numerical results show that a concave periodic boundary geometry can decrease the average cell concentration at the boundary by more than 50% relative to a flat surface.

Microrotor of a chain-grafted colloidal disk immersed in the active bath: The impact of particle concentration, grafting density, and chain rigidity

In an earlier work, we discussed the possibility to realize a microrotor by immersing a chain-grafted colloidal disk in a thin film of active-particle suspension. Under certain conditions, the colloidal disk rotates unidirectionally driven by the bath active particles. Here we systematically study the role of active-particle concentration, grafting density, and chain rigidity in the phenomenon of the spontaneous symmetry breaking of the chain configurations and the unidirectional rotation of the disk. We find that high chain rigidity can help stabilize both the collective asymmetric chain configurations and the rotation of the disk, while it has a weak impact on the rotational speed/efficiency. Increasing the number of grafted chains can also stabilize the rotation but has a negative impact on the rotational speed/efficiency. Active particles power the rotation of the colloidal disk, yet their contribution saturates beyond a certain concentration. Our work provides new insights into the active systems with chain-structured objects and the design of soft/deformable micromachines.

Current reversals of active particles in time-oscillating potentials

Rectification of interacting active particles is numerically investigated in a two-dimensional time-oscillating potential. It is found that the oscillation of the potential and the self-propulsion of active particles are two different types of nonequilibrium driving, which can induce net currents with opposite directions. For a given asymmetry of the potential, the direction of the transport is determined by the competition of the self-propulsion and the oscillation of the potential. There exists an optimal oscillating angular frequency (or self-propulsion speed) at which the average velocity takes its maximal positive or negative value. Remarkably, when the oscillation of the potential competes with the self-propulsion, the average velocity can change direction several times due to the change in the oscillating frequency. Especially, particles with different self-propulsion velocities will move in opposite directions and can be separated. Our results provide a novel and convenient method for controlling and manipulating the transport (or separation) of active particles.

Active fluids at circular boundaries: swim pressure and anomalous droplet ripening

We investigate the swim pressure exerted by non-chiral and chiral active particles on convex or concave circular boundaries. Active particles are modeled as non-interacting and non-aligning self-propelled Brownian particles. The convex and concave circular boundaries are used to model a fixed inclusion immersed in an active bath and a cavity (or container) enclosing the active particles, respectively. We first present a detailed analysis of the role of convex versus concave boundary curvature and of the chirality of active particles in their spatial distribution, chirality-induced currents, and the swim pressure they exert on the bounding surfaces. The results will then be used to predict the mechanical equilibria of suspended fluid enclosures (generically referred to as 'droplets') in a bulk with active particles being present either inside the bulk fluid or within the suspended droplets. We show that, while droplets containing active particles behave in accordance with standard capillary paradigms when suspended in a normal bulk, those containing a normal fluid exhibit anomalous behaviors when suspended in an active bulk. In the latter case, the excess swim pressure results in non-monotonic dependence of the inside droplet pressure on the droplet radius; hence, revealing an anomalous regime of behavior beyond a threshold radius, in which the inside droplet pressure increases upon increasing the droplet size. Furthermore, for two interconnected droplets, mechanical equilibrium can occur also when the droplets have different sizes. We thus identify a regime of anomalous droplet ripening, where two unequal-sized droplets can reach a final state of equal size upon interconnection, in stark contrast with the standard Ostwald ripening phenomenon, implying shrinkage of the smaller droplet in favor of the larger one.

Active motion on curved surfaces

A theoretical analysis of active motion on curved surfaces is presented in terms of a generalization of the telegrapher equation. Such a generalized equation is explicitly derived as the polar approximation of the hierarchy of equations obtained from the corresponding Fokker-Planck equation of active particles diffusing on curved surfaces. The general solution to the generalized telegrapher equation is given for a pulse with vanishing current as initial data. Expressions for the probability density and the mean squared geodesic displacement are given in the limit of weak curvature. As an explicit example of the formulated theory, the case of active motion on the sphere is presented, where oscillations observed in the mean squared geodesic displacement are explained.

Clogging and depinning of ballistic active matter systems in disordered media

We numerically examine ballistic active disks driven through a random obstacle array. Formation of a pinned or clogged state occurs at much lower obstacle densities for the active disks than for passive disks. As a function of obstacle density, we identify several distinct phases including a depinned fluctuating cluster state, a pinned single-cluster or jammed state, a pinned multicluster state, a pinned gel state, and a pinned disordered state. At lower active disk densities, a drifting uniform liquid forms in the absence of obstacles, but when even a small number of obstacles are introduced, the disks organize into a pinned phase-separated cluster state in which clusters nucleate around the obstacles, similar to a wetting phenomenon. We examine how the depinning threshold changes as a function of disk or obstacle density and find a crossover from a collectively pinned cluster state to a disordered plastic depinning transition as a function of increasing obstacle density. We compare this to the behavior of nonballistic active particles and show that as we vary the activity from completely passive to completely ballistic, a clogged phase-separated state appears in both the active and passive limits, while for intermediate activity, a readily flowing liquid state appears and there is an optimal activity level that maximizes the flux through the sample.

Anomalous Chained Turbulence in Actively Driven Flows on Spheres

Recent experiments demonstrate the importance of substrate curvature for actively forced fluid dynamics. Yet, the covariant formulation and analysis of continuum models for nonequilibrium flows on curved surfaces still poses theoretical challenges. Here, we introduce and study a generalized covariant Navier-Stokes model for fluid flows driven by active stresses in nonplanar geometries. The analytical tractability of the theory is demonstrated through exact stationary solutions for the case of a spherical bubble geometry. Direct numerical simulations reveal a curvature-induced transition from a burst phase to an anomalous turbulent phase that differs distinctly from externally forced classical 2D Kolmogorov turbulence. This new type of active turbulence is characterized by the self-assembly of finite-size vortices into linked chains of antiferromagnetic order, which percolate through the entire fluid domain, forming an active dynamic network. The coherent motion of the vortex chain network provides an efficient mechanism for upward energy transfer from smaller to larger scales, presenting an alternative to the conventional energy cascade in classical 2D turbulence.

Active Brownian particles near straight or curved walls: Pressure and boundary layers

Unlike equilibrium systems, active matter is not governed by the conventional laws of thermodynamics. Through a series of analytic calculations and Langevin dynamics simulations, we explore how systems cross over from equilibrium to active behavior as the activity is increased. In particular, we calculate the profiles of density and orientational order near straight or circular walls and show the characteristic width of the boundary layers. We find a simple relationship between the enhancements of density and pressure near a wall. Based on these results, we determine how the pressure depends on wall curvature and hence make approximate analytic predictions for the motion of curved tracers, as well as the rectification of active particles around small openings in confined geometries.

Hierarchical collective motion of a mixture of active dipolar Janus particles and passive charged colloids in two dimensions

We use computer simulations to study the behavior of a mixture of large passive charged colloids in a suspension of smaller active dipolar Janus particles. We find that when a single charged colloid is present in solution, it acquires a rotational or translational motion depending on how the active dipoles self-assemble on its surface to form active complexes. The collective behavior of these complexes is quite remarkable, and includes swarming behavior and coherent macroscopic motion. We detail how the variety of different phenomenologies emerging in this system can ultimately be controlled by the strength of the active forces and the relative concentration of the two species.

Memory effects in active particles with exponentially correlated propulsion

The Ornstein-Uhlenbeck particle (OUP) model imagines a microscopic swimmer propelled by an active force which is correlated with itself on a finite time scale. Here we investigate the influence of external potentials on an ideal suspension of OUPs, in both one and two spatial dimensions, with particular attention paid to the pressure exerted on "confining walls." We employ a mathematical connection between the local density of OUPs and the statistics of their propulsion force to demonstrate the existence of an equation of state in one dimension. In higher dimensions we show that active particles generate a nonconservative force field in the surrounding medium. A simplified far-from-equilibrium model is proposed to account for OUP behavior in the vicinity of potentials. Building on this, we interpret simulations of OUPs in more complicated situations involving asymmetrical and spatially curved potentials, and characterize the resulting inhomogeneous stresses in terms of competing active length scales.

Van't Hoff's law for active suspensions: the role of the solvent chemical potential

We extend Van't Hoff's law for the osmotic pressure to a suspension of active Brownian particles. The propelled particles exert a net reaction force on the solvent, and thereby either drive a measurable solvent flow from the connecting solvent reservoir through the semipermeable membrane, or increase the osmotic pressure and cause the suspension to rise to heights as large as micrometers for experimentally realized microswimmers described in the literature. The increase in osmotic pressure is caused by the background solvent being, in contrast to passive suspensions, no longer at the chemical potential of the solvent reservoir. The difference in solvent chemical potentials depends on the colloid-membrane interaction potential, which implies that the osmotic pressure is a state function of a state that itself is influenced by the membrane potential.

Enhanced Dynamics of Confined Cytoskeletal Filaments Driven by Asymmetric Motors

Cytoskeletal filaments and molecular motors facilitate the micron-scale force generation necessary for the distribution of organelles and the restructuring of the cytoskeleton within eukaryotic cells. Although the mesoscopic structure and the dynamics of such filaments have been studied in vitro and in vivo, their connection with filament polarity-dependent motor-mediated force generation is not well understood. Using 2D Brownian dynamics simulations, we study a dense, confined mixture of rigid microtubules (MTs) and active springs that have arms that cross-link neighboring MT pairs and move unidirectionally on the attached MT. We simulate depletion interactions between MTs using an attractive potential. We show that dimeric motors, with a motile arm on only one of the two MTs, produce large polarity-sorted MT clusters, whereas tetrameric motors, with motile arms on both microtubules, produce bundles. Furthermore, dimeric motors induce, on average, higher velocities between antialigned MTs than tetrameric motors. Our results, where MTs move faster near the confining wall, are consistent with experimental observations in Drosophila oocytes where enhanced microtubule activity is found close to the confining plasma membrane.

Forces in inhomogeneous open active-particle systems

We study the force that noninteracting pointlike active particles apply to a symmetric inert object in the presence of a gradient of activity and particle sources and sinks. We consider two simple patterns of sources and sinks that are common in biological systems. We analytically solve a one-dimensional model designed to emulate higher-dimensional systems, and study a two-dimensional model by numerical simulations. We specify when the particle flux due to the creation and annihilation of particles can act to smooth the density profile that is induced by a gradient in the velocity of the active particles, and find the net resultant force due to both the gradient in activity and the particle flux. These results are compared qualitatively to observations of nuclear motion inside the oocyte, that is driven by a gradient in activity of actin-coated vesicles.

Spontaneous symmetry breaking induced unidirectional rotation of a chain-grafted colloidal particle in the active bath

Exploiting the energy of randomly moving active agents such as bacteria is a fascinating way to power a microdevice. Here we show, by simulations, that a chain-grafted disk-like colloidal particle can rotate unidirectionally and hence output work when immersed in a thin film of active particle suspension. The collective spontaneous symmetry breaking of chain configurations is the origin of the unidirectional rotation. Long persistence time, large propelling force and/or small rotating friction are keys to sustaining the collective broken symmetry and realizing the rotation. In the rotating state, we find very simple linear relations, e.g. between the mean angular speed and the propelling force. The time-evolving asymmetry of chain configurations reveals that there are two types of non-rotating state. The basic phenomena are also observed in the macroscopic granular experiments, implying the generic nature of these phenomena. Our findings provide new insights into the collective spontaneous symmetry breaking in active systems with flexible objects and also open the way to conceive new soft/deformable microdevices.

Generalized Archimedes' principle in active fluids

We show how a gradient in the motility properties of noninteracting pointlike active particles can cause a pressure gradient that pushes a large inert object. We calculate the force on an object inside a system of active particles with position-dependent motion parameters, in one and two dimensions, and show that a modified Archimedes' principle is satisfied. We characterize the system, both in terms of the model parameters and in terms of experimentally measurable quantities: the spatial profiles of the density, velocity and pressure. This theoretical analysis is motivated by recent experiments, which showed that the nucleus of a mouse oocyte (immature egg cell) moves from the cortex to the center due to a gradient of activity of vesicles propelled by molecular motors; it more generally applies to artificial systems of controlled localized activity.

Equilibrium mappings in polar-isotropic confined active particles

Despite their fundamentally nonequilibrium nature, the individual and collective behavior of active systems with polar propulsion and isotropic interactions (polar-isotropic active systems) are remarkably well captured by equilibrium mapping techniques. Here we examine two signatures of equilibrium systems --the existence of a local free energy function and the independence of the coarse-grained behavior on the details of the microscopic dynamics-- in polar-isotropic active particles confined by hard walls of arbitrary geometry at the one-particle level. We find that boundaries that possess concave regions make the density profile strongly dynamics-dependent and give it a nonlocal dependence on the geometry of the confining box. This in turn constrains the scope of equilibrium mapping techniques in polar-isotropic active systems.

Transport of underdamped active particles in ratchet potentials

We study the rectified transport of underdamped active noninteracting particles in an asymmetric periodic potential. It is found that the ratchet effect of active noninteracting particles occurs in a single direction (along the easy direction of the substrate asymmetry) in the overdamped limit. However, when the inertia is considered, it is possible to observe reversals of the ratchet effect, where the motion is along the hard direction of the substrate asymmetry. By changing the friction coefficient or the self-propulsion force, the average velocity can change its direction several times. Therefore, by suitably tailoring the parameters, underdamped active particles with different self-propulsion forces can move in different directions and can be separated.

Effects of hydrodynamic interactions on rectified transport of self-propelled particles

Directed transport of self-propelled particles is numerically investigated in a three-dimensional asymmetric potential. Beside the steric repulsive forces, hydrodynamic interactions between particles have been taken into account in an approximate way. From numerical simulations, we find that hydrodynamic interactions can strongly affect the rectified transport of self-propelled particles. Hydrodynamic interactions enhance the performance of the rectified transport when particles can easily pass across the barrier of the potential, and reduce the rectified transport when particles are mainly trapped in the potential well.

Gaussian theory for spatially distributed self-propelled particles

Obtaining a reduced description with particle and momentum flux densities outgoing from the microscopic equations of motion of the particles requires approximations. The usual method, we refer to as truncation method, is to zero Fourier modes of the orientation distribution starting from a given number. Here we propose another method to derive continuum equations for interacting self-propelled particles. The derivation is based on a Gaussian approximation (GA) of the distribution of the direction of particles. First, by means of simulation of the microscopic model, we justify that the distribution of individual directions fits well to a wrapped Gaussian distribution. Second, we numerically integrate the continuum equations derived in the GA in order to compare with results of simulations. We obtain that the global polarization in the GA exhibits a hysteresis in dependence on the noise intensity. It shows qualitatively the same behavior as we find in particles simulations. Moreover, both global polarizations agree perfectly for low noise intensities. The spatiotemporal structures of the GA are also in agreement with simulations. We conclude that the GA shows qualitative agreement for a wide range of noise intensities. In particular, for low noise intensities the agreement with simulations is better as other approximations, making the GA to an acceptable candidates of describing spatially distributed self-propelled particles.

Kinetics of motility-induced phase separation and swim pressure

Active Brownian particles (ABPs) represent a minimal model of active matter consisting of self-propelled spheres with purely repulsive interactions and rotational noise. Here we examine the pressure of ABPs in two dimensions in both closed boxes and systems with periodic boundary conditions and show that its nonmonotonic behavior with density is a general property of ABPs and is not the result of finite-size effects. We correlate the time evolution of the mean pressure towards its steady-state value with the kinetics of motility-induced phase separation. For parameter values corresponding to phase-separated steady states, we identify two dynamical regimes. The pressure grows monotonically in time during the initial regime of rapid cluster formation, overshooting its steady-state value and then quickly relaxing to it, and remains constant during the subsequent slower period of cluster coalescence and coarsening. The overshoot is a distinctive feature of active systems.

Active Brownian particles and run-and-tumble particles separate inside a maze

A diverse range of natural and artificial self-propelled particles are known and are used nowadays. Among them, active Brownian particles (ABPs) and run-and-tumble particles (RTPs) are two important classes. We numerically study non-interacting ABPs and RTPs strongly confined to different maze geometries in two dimensions. We demonstrate that by means of geometrical confinement alone, ABPs are separable from RTPs. By investigating Matryoshka-like mazes with nested shells, we show that a circular maze has the best filtration efficiency. Results on the mean first-passage time reveal that ABPs escape faster from the center of the maze, while RTPs reach the center from the rim more easily. According to our simulations and a rate theory, which we developed, ABPs in steady state accumulate in the outermost region of the Matryoshka-like mazes, while RTPs occupy all locations within the maze with nearly equal probability. These results suggest a novel technique for separating different types of self-propelled particles by designing appropriate confining geometries without using chemical or biological agents.

Shape and Displacement Fluctuations in Soft Vesicles Filled by Active Particles

We investigate numerically the dynamics of shape and displacement fluctuations of two-dimensional flexible vesicles filled with active particles. At low concentration most of the active particles accumulate at the boundary of the vesicle where positive particle number fluctuations are amplified by trapping, leading to the formation of pinched spots of high density, curvature and pressure. At high concentration the active particles cover the vesicle boundary almost uniformly, resulting in fairly homogeneous pressure and curvature, and nearly circular vesicle shape. The change between polarized and spherical shapes is driven by the number of active particles. The center-of-mass of the vesicle performs a persistent random walk with a long time diffusivity that is strongly enhanced for elongated active particles due to orientational correlations in their direction of propulsive motion. In our model shape-shifting induces directional sensing and the cell spontaneously migrate along the polarization direction.

Active Particles with Soft and Curved Walls: Equation of State, Ratchets, and Instabilities

We study, from first principles, the pressure exerted by an active fluid of spherical particles on general boundaries in two dimensions. We show that, despite the nonuniform pressure along curved walls, an equation of state is recovered upon a proper spatial averaging. This holds even in the presence of pairwise interactions between particles or when asymmetric walls induce ratchet currents, which are accompanied by spontaneous shear stresses on the walls. For flexible obstacles, the pressure inhomogeneities lead to a modulational instability as well as to the spontaneous motion of short semiflexible filaments. Finally, we relate the force exerted on objects immersed in active baths to the particle flux they generate around them.

Ideal bulk pressure of active Brownian particles

The extent to which active matter might be described by effective equilibrium concepts like temperature and pressure is currently being discussed intensely. Here, we study the simplest model, an ideal gas of noninteracting active Brownian particles. While the mechanical pressure exerted onto confining walls has been linked to correlations between particles' positions and their orientations, we show that these correlations are entirely controlled by boundary effects. We also consider a definition of local pressure, which describes interparticle forces in terms of momentum exchange between different regions of the system. We present three pieces of analytical evidence which indicate that such a local pressure exists, and we show that its bulk value differs from the mechanical pressure exerted on the walls of the system. We attribute this difference to the fact that the local pressure in the bulk does not depend on boundary effects, contrary to the mechanical pressure. We carefully examine these boundary effects using a channel geometry, and we show a virial formula for the pressure correctly predicts the mechanical pressure even in finite channels. However, this result no longer holds in more complex geometries, as exemplified for a channel that includes circular obstacles.

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.

Swim pressure on walls with curves and corners

The concept of swim pressure quantifies the average force exerted by microswimmers on confining walls in nonequilibrium. Here we explore how the swim pressure depends on the wall curvature and on the presence of sharp corners in the wall. For active Brownian particles at high dilution, we present a coherent framework which describes the force and torque on passive particles of arbitrary shape, in the limit of large particles compared to the persistence length of the swimmer trajectories. The resulting forces can be used to derive, for example, the activity-induced depletion interaction between two disks, as well as to optimize the shape of a tracer particle for high swimming velocity. Our predictions are verifiable in experiments on passive obstacles exposed to a bath of bacteria or artificial microswimmers.