Collective forces in scalar active matter

Partial and complete wetting of droplets of active Brownian particles

We study wetting droplets formed of active Brownian particles in contact with a repulsive potential barrier, in a wedge geometry. Our numerical results demonstrate a transition between partially wet and completely wet states, as a function of the barrier height, analogous to the corresponding surface phase transition in passive fluids. We analyse partially wet configurations characterised by a nonzero contact angle θ between the droplet surface and the barrier including the average density profile and its fluctuations. These findings are compared with two equilibrium systems: a Lennard-Jones fluid and a simple contour model for a liquid-vapour interface. We locate the wetting transition where cos(θ) = 1, and the neutral state where cos(θ) = 0. We discuss the implications of these results for possible definitions of surface tensions in active fluids.

Nonequilibrium structure formation in electrohydrodynamic emulsions

Application of an electric field across the interface of two fluids with low, but non-zero, conductivities gives rise to a sustained electrohydrodynamic (EHD) fluid flow. In the presence of neighboring drops, drops interact via the EHD flows of their neighbors, as well as through a dielectrophoretic (DEP) force, a consequence of drops encountering disturbance electric fields around their neighbors. We explore the collective dynamics of emulsions with drops undergoing EHD and DEP interactions. The interplay between EHD and DEP results in a rich set of emergent behaviors. We simulate the collective behavior of large numbers of drops; in two dimensions, where drops are confined to a plane; and three dimensions. In monodisperse emulsions, drops in two dimensions cluster or crystallize depending on the relative strengths of EHD and DEP, and form spaced clusters when EHD and DEP balance. In three dimensions, chain formation observed under DEP alone is suppressed by EHD, and lost entirely when EHD dominates. When a second population of drops are introduced, such that the electrical conductivity, permittivity, or viscosity are different from the first population of drops, the interaction between the drops becomes non-reciprocal, an apparent violation of Newton's Third Law. The breadth of consequences due to these non-reciprocal interactions are vast: we show selected cases in two dimensions, where drops cluster into active dimers, trimers, and larger clusters that continue to translate and rotate over long timescales; and three dimensions, where drops form stratified chains, or combine into a single dynamic sheet.

Force Generation in Confined Active Fluids: The Role of Microstructure

We experimentally determine the force exerted by a bath of active particles onto a passive probe as a function of its distance to a wall and compare it to the measured averaged density distribution of active particles around the probe. Within the framework of an active stress, we demonstrate that both quantities are-up to a factor-directly related to each other. Our results are in excellent agreement with a minimal numerical model and confirm a general and system-independent relationship between the microstructure of active particles and transmitted forces.

Anomalous Transport of Tracers in Active Baths

We derive the long-time dynamics of a tracer immersed in a one-dimensional active bath. In contrast to previous studies, we find that the damping and noise correlations possess long-time tails with exponents that depend on the tracer symmetry. For generic tracers, shape asymmetry induces ratchet effects that alter fluctuations and lead to superdiffusion and friction that grows with time when the tracer is dragged at a constant speed. In the singular limit of a completely symmetric tracer, we recover normal diffusion and finite friction. Furthermore, for small symmetric tracers, the active contribution to the friction becomes negative: active particles enhance motion rather than oppose it. These results show that, in low-dimensional systems, the motion of a passive tracer in an active bath cannot be modeled as a persistent random walker with a finite correlation time.

Active depletion torque between two passive rods

The active depletion torque experienced by two anisotropic objects in an active bath is a conceptional generalization of the equilibrium entropic torque. Using Brownian dynamics simulations, we compute the active depletion torque suffered by two passive rods immersed in an ensemble of active Brownian particles. Our results demonstrate that the active depletion torque is qualitatively different from its passive counterpart. Interestingly, we find that the active depletion torque can be greatly affected by the external constraint applied on the rotational degree of freedom of the rods, and even the direction may be changed with the orientational constraint, which is in contrast to the equilibrium depletion torque. The main reason for the remarkable features of the active depletion torque is that the active particles can significantly accumulate in the vicinity of the rods due to persistent self-propulsion, which is sensitively dependent on the constraint strength and the rod configurations. Our findings could be relevant for understanding the self-assembly and dynamics of anisotropic macromolecules in living environments.

Critical behavior of active Brownian particles: Connection to field theories

We explore the relation between active Brownian particles, a minimal particle-based model for active matter, and scalar field theories. Both show a liquid-gas-like phase transition toward stable coexistence of a dense liquid with a dilute active gas that terminates in a critical point. However, a comprehensive mapping between the particle-based model parameters and the effective coefficients governing the field theories has not been established yet. We discuss conflicting recent numerical results for the critical exponents of active Brownian particles in two dimensions. Starting from the intermediate effective hydrodynamic equations, we then present a construction for a scalar order parameter for active Brownian particles that yields the active model B+. We argue that a crucial ingredient is the coupling between density and polarization in the particle current. The renormalization flow close to two dimensions exhibits a pair of perturbative fixed points that limit the attractive basin of the Wilson-Fisher fixed point, with the perspective that the critical behavior of active Brownian particles in two dimensions is governed by a strong-coupling fixed point different from Wilson-Fisher and not necessarily corresponding to Ising universality.

Wetting Transition of Active Brownian Particles on a Thin Membrane

We study nonequilibrium analogues of surface phase transitions in a minimal model of active particles in contact with a purely repulsive potential barrier that mimics a thin porous membrane. Under conditions of bulk motility-induced phase separation, the interaction strength ϵ_{w} of the barrier controls the affinity of the dense phase for the barrier region. We uncover clear signatures of a wetting phase transition as ϵ_{w} is varied. In common with its equilibrium counterpart, the character of this transition depends on the system dimensionality: a continuous transition with large density fluctuations and gas bubbles is uncovered in 2D while 3D systems exhibit a sharp transition absent of large correlations.

Noncentral forces mediated between two inclusions in a bath of active Brownian rods

Using Brownian Dynamics simulations, we study effective interactions mediated between two identical and impermeable disks (inclusions) immersed in a bath of identical, active (self-propelled), Brownian rods in two spatial dimensions, by assuming that the self-propulsion axis of the rods may generally deviate from their longitudinal axis. When the self-propulsion is transverse (perpendicular to the rod axis), the accumulation of active rods around the inclusions is significantly enhanced, causing a more expansive steric layering (ring formation) of the rods around the inclusions, as compared with the reference case of longitudinally self-propelling rods. As a result, the transversally self-propelling rods also mediate a significantly longer ranged effective interaction between the inclusions. The bath-mediated interaction arises due to the overlaps between the active-rod rings formed around the inclusions, as they are brought into small separations. When the self-propulsion axis is tilted relative to the rod axis, we find an asymmetric imbalance of active-rod accumulation around the inclusion dimer. This leads to a noncentral interaction, featuring an anti-parallel pair of transverse force components and, hence, a bath-mediated torque on the dimer.

Active matter shepherding and clustering in inhomogeneous environments

We consider a mixture of active and passive run-and-tumble disks in an inhomogeneous environment where only half of the sample contains quenched disorder or pinning. The disks are initialized in a fully mixed state of uniform density. We identify several distinct dynamical phases as a function of motor force and pinning density. At high pinning densities and high motor forces, there is a two-step process initiated by a rapid accumulation of both active and passive disks in the pinned region, which produces a large density gradient in the system. This is followed by a slower species phase separation process where the inactive disks are shepherded by the active disks into the pin-free region, forming a nonclustered fluid and producing a more uniform density with species phase separation. For higher pinning densities and low motor forces, the dynamics becomes very slow and the system maintains a strong density gradient. For weaker pinning and large motor forces, a floating clustered state appears, and the time-averaged density of the system is uniform. We illustrate the appearance of these phases in a dynamic phase diagram.

Clogging, dynamics, and reentrant fluid for active matter on periodic substrates

We examine the collective states of run-and-tumble active matter disks driven over a periodic obstacle array. When the drive is applied along a symmetry direction of the array, we find a clog-free uniform liquid state for low activity, while at higher activity, the density becomes increasingly heterogeneous and an active clogged state emerges in which the mobility is strongly reduced. For driving along nonsymmetry or incommensurate directions, there are two different clogging behaviors consisting of a drive-dependent clogged state in the low activity thermal limit and a drive-independent clogged state at high activity. These regimes are separated by a uniform flowing liquid at intermediate activity. There is a critical activity level above which the thermal clogged state does not occur, as well as an optimal activity level that maximizes the disk mobility. Thermal clogged states are dependent on the driving direction while active clogged states are not. In the low activity regime, diluting the obstacles produces a monotonic increase in the mobility; however, for large activities, the mobility is more robust against obstacle dilution. We also examine the velocity-force curves for driving along nonsymmetry directions and find that they are linear when the activity is low or intermediate but become nonlinear at high activity and show behavior similar to that found for the plastic depinning of solids. At higher drives, the active clustering is lost. For low activity, we also find a reentrant fluid phase, where the system transitions from a high mobility fluid at low drives to a clogged state at higher drives and then back into another fluid phase at very high drives. We map the regions in which the thermally clogged, partially clogged, active uniform fluid, clustered fluid, active clogged, and directionally locked states occur as a function of disk density, drift force, and activity.

Motility-induced inter-particle correlations and dynamics: a microscopic approach for active Brownian particles

Amongst the theoretical approaches towards the dynamics and phase behaviour of suspensions of active Brownian particles (ABPs), no attempt has been made to specify the motility-induced inter-particle correlations as quantified by the pair-correlation function. Here, we derive expressions for the pair-correlation function for ABPs with very short-ranged direct interactions for small and large swimming velocities and low concentrations. The pair-correlation function is the solution of a differential equation that is obtained from the Fokker-Planck equation for the probability density function of the positions and orientations of the ABPs. For large swimming Peclet numbers, λ, the pair-correlation function is highly asymmetric. The pair-correlation function attains a large value, ∼λ, within a small region of spatial extent, ∼1/λ, near contact of the ABPs when the ABPs approach each other. The pair-correlation function is small within a large region of spatial extent, ∼λ1/3, when the ABPs move apart, with a contact value that is essentially zero. From the explicit expressions for the pair-correlation function, Fick's diffusion equation is generalized to include motility. It is shown that mass transport, in case of large swimming velocities, is dominated by a preferred swimming direction that is induced by concentration gradients. The expression for the pair-correlation function derived in this paper could serve as a starting point to obtain approximate results for high concentrations, which could then be employed in a first-principles analysis of the dynamics and phase behaviour of ABPs at higher concentrations.

Design of nematic liquid crystals to control microscale dynamics

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

Vorticity Determines the Force on Bodies Immersed in Active Fluids

When immersed into a fluid of active Brownian particles, passive bodies might start to undergo linear or angular directed motion depending on their shape. Here we exploit the divergence theorem to relate the forces responsible for this motion to the density and current induced by-but far away from-the body. In general, the force is composed of two contributions: due to the strength of the dipolar field component and due to particles leaving the boundary, generating a nonvanishing vorticity of the polarization. We derive and numerically corroborate results for periodic systems, which are fundamentally different from unbounded systems with forces that scale with the area of the system. We demonstrate that vorticity is localized close to the body and to points at which the local curvature changes, enabling the rational design of particle shapes with desired propulsion properties.

Shaping the gradients driving phoretic micro-swimmers: influence of swimming speed, budget of carbonic acid and environment

pH gradient-driven modular micro-swimmers are investigated as a model for a large variety of quasi-two-dimensional chemi-phoretic self-propelled entities. Using three-channel micro-photometry, we obtain a precise large field mapping of pH at a spatial resolution of a few microns and a pH resolution of [Formula: see text] units for swimmers of different velocities propelling on two differently charged substrates. We model our results in terms of solutions of the three-dimensional advection-diffusion equation for a 1:1 electrolyte, i.e. carbonic acid, which is produced by ion exchange and consumed by equilibration with dissolved [Formula: see text]. We demonstrate the dependence of gradient shape and steepness on swimmer speed, diffusivity of chemicals, as well as the fuel budget. Moreover, we experimentally observe a subtle, but significant feedback of the swimmer's immediate environment in terms of a substrate charge-mediated solvent convection. We discuss our findings in view of different recent results from other micro-fluidic or active matter investigations. We anticipate that they are relevant for quantitative modelling and targeted applications of diffusio-phoretic flows in general and artificial micro-swimmers in particular.

Wetting of a solid surface by active matter

A lattice model is used to study repulsive active particles at a planar surface. A rejection-free Kinetic Monte Carlo method is employed to characterize the wetting behaviour. The model predicts a motility-induced phase separation of active particles, and the bulk coexistence of dense liquid-like and dilute vapour-like steady states is determined. An "ensemble", with a varying number of particles, analogous to a grand canonical ensemble in equilibrium, is introduced. The formation and growth of the liquid film between the solid surface and the vapour phase is investigated. At constant activity, as the system is brought towards coexistence from the vapour side, the thickness of the adsorbed film exhibits a divergent behaviour regardless of the activity. This suggests a complete wetting scenario along the full coexistence curve.

Coexistence of active Brownian disks: van der Waals theory and analytical results

At thermal equilibrium, intensive quantities like temperature and pressure have to be uniform throughout the system, restricting inhomogeneous systems composed of different phases. The paradigmatic example is the coexistence of vapor and liquid, a state that can also be observed for active Brownian particles steadily driven away from equilibrium. Recently, a strategy has been proposed that allows to predict phase equilibria of active particles [Solon et al., Phys. Rev. E 97, 020602(R) (2018)2470-004510.1103/PhysRevE.97.020602]. Here we elaborate on this strategy and formulate it in the framework of a van der Waals theory for active disks. For a given equation of state, we derive the effective free energy analytically and show that it yields coexisting densities in very good agreement with numerical results. We discuss the interfacial tension and the relation to Cahn-Hilliard models.

Spatiotemporal dynamics of a self-propelled system with opposing alignment and repulsive forces

Effect of concurrent alignment and repulsion is studied in the purview of a confined active matter system using a modified force-based Vicsek model. On alteration of the alignment and the repulsive force parameters, a low alignment random phase, a midrange alignment milling phase, and a high alignment oscillatory phase are identified. Based on the particle aggregations, the milling phase is further classified into three subphases, two of which are spatial patterns: one consisting of compact ring-shaped mills and the other incorporating both rings and clusters. A correlation function based on the inner product of spatial velocity fluctuations of the particles shows a high correlation length for the ringed milling and the rings-clusters hybrid milling state. On analyzing temporal velocity fluctuations of particles through chaos detection techniques, low alignment and high alignment states are indicative of chaos, while the middle order alignment is symbolic of periodicity. The extent of synchronization of the particles' motion is analyzed through a Hilbert transform-based mean frequency approach, leading to the detection of a weak chimera state in the case of the spatial structures. The ringed milling state shows a unique category of weak chimera consisting of multiple oscillator groups showcasing different synchronization frequencies coexisting with desynchronized oscillators.

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.

Directional locking effects for active matter particles coupled to a periodic substrate

Directional locking occurs when a particle moving over a periodic substrate becomes constrained to travel along certain substrate symmetry directions. Such locking effects arise for colloids and superconducting vortices moving over ordered substrates when the direction of the external drive is varied. Here we study the directional locking of run-and-tumble active matter particles interacting with a periodic array of obstacles. In the absence of an external biasing force, we find that the active particle motion locks to various symmetry directions of the substrate when the run time between tumbles is large. The number of possible locking directions depends on the array density and on the relative sizes of the particles and the obstacles. For a square array of large obstacles, the active particle only locks to the x, y, and 45^{∘} directions, while for smaller obstacles, the number of locking angles increases. Each locking angle satisfies θ=arctan(p/q), where p and q are integers, and the angle of motion can be measured using the ratio of the velocities or the velocity distributions in the x and y directions. When a biasing driving force is applied, the directional locking behavior is affected by the ratio of the self-propulsion force to the biasing force. For large biasing, the behavior resembles that found for directional locking in passive systems. For large obstacles under biased driving, a trapping behavior occurs that is nonmonotonic as a function of increasing run length or increasing self-propulsion force, and the trapping diminishes when the run length is sufficiently large.

Effective interactions mediated between two permeable disks in an active fluid

We study steady-state properties of a bath of active Brownian particles (ABPs) in two dimensions in the presence of two fixed, permeable (hollow) disklike inclusions, whose interior and exterior regions can exhibit mismatching motility (self-propulsion) strengths for the ABPs. We show that such a discontinuous motility field strongly affects spatial distribution of ABPs and thus also the effective interaction mediated between the inclusions through the active bath. Such net interactions arise from soft interfacial repulsions between ABPs that sterically interact with and/or pass through permeable membranes assumed to enclose the inclusions. Both regimes of repulsion and attractive (albeit with different mechanisms) are reported and summarized in overall phase diagrams.

Universal reshaping of arrested colloidal gels via active doping

Colloids that interact via a short-range attraction serve as the primary building blocks for a broad range of self-assembled materials. However, one of the well-known drawbacks to this strategy is that these building blocks rapidly and readily condense into a metastable colloidal gel. Using computer simulations, we illustrate how the addition of a small fraction of purely repulsive self-propelled colloids, a technique referred to as active doping, can prevent the formation of this metastable gel state and drive the system toward its thermodynamically favored crystalline target structure. The simplicity and robust nature of this strategy offers a systematic and generic pathway to improving the self-assembly of a large number of complex colloidal structures. We discuss in detail the process by which this feat is accomplished and provide quantitative metrics for exploiting it to modulate the self-assembly. We provide evidence for the generic nature of this approach by demonstrating that it remains robust under a number of different anisotropic short-ranged pair interactions in both two and three dimensions. In addition, we report on a novel microphase in mixtures of passive and active colloids. For a broad range of self-propelling velocities, it is possible to stabilize a suspension of fairly monodisperse finite-size crystallites. Surprisingly, this microphase is also insensitive to the underlying pair interaction between building blocks. The active stabilization of these moderately sized monodisperse clusters is quite remarkable and should be of great utility in the design of hierarchical self-assembly strategies. This work further bolsters the notion that active forces can play a pivotal role in directing colloidal self-assembly.