Dynamic clustering and chemotactic collapse of self-phoretic active particles

Artificial chemotaxis under electrodiffusiophoresis

HYPOTHESIS

Through a large parameter space, electric fields can tune colloidal interactions and forces leading to diverse static and dynamical structures. So far, however, field-driven interactions have been limited to dipole-dipole and hydrodynamic contributions. Nonetheless, in this work, we propose that under the right conditions, electric fields can also induce interactions based on local chemical fields and diffusiophoretic flows.

EXPERIMENTS

Herein, we present a strategy to generate and measure 3D chemical gradients under electric fields. In this approach, faradaic reactions at electrodes induce global pH gradients that drive long-range transport through electrodiffusiophoresis. Simultaneously, the electric field induces local pH gradients by driving the particle's double layer far from equilibrium.

FINDINGS

As a result, while global pH gradients lead to 2D focusing away from electrodes, local pH gradients induce aggregation in the third dimension. Evidence points to a mechanism of interaction based on diffusiophoresis. Interparticle interactions display a strong dependence on surface chemistry, zeta potential and diameter of particles. Furthermore, pH gradients can be readily tuned by adjusting the voltage and frequency of the electric field. For large Péclet numbers, we observed a collective chemotactic-like collapse of particles. Remarkably, such collapse occurs without reactions at a particle's surface. By mixing particles with different sizes, we also demonstrate, through experiments and Brownian dynamics simulations, the emergence of non-reciprocal interactions, where small particles are more drawn towards large ones.

Chemotactic behavior for a self-phoretic Janus particle near a patch source of fuel

Many biological microswimmers are capable of chemotaxis, i.e., they can sense an ambient chemical gradient and adjust their mechanism of motility to move towards or away from the source of the gradient. Synthetic active colloids endowed with chemotactic behavior hold considerable promise for targeted drug delivery and the realization of programmable and reconfigurable materials. Here, we study the chemotactic behavior of a Janus particle, which converts "fuel" molecules, released at an axisymmetric chemical patch located on a planar wall, into "product" molecules at its catalytic cap and moves by self-phoresis induced by the product. The chemotactic behavior is characterized as a function of the interplay between the rates of release (at the patch) and the consumption (at the particle) of fuel, as well as of details of the phoretic response of the particle (i.e., its phoretic mobility). Among other results, we find that, under certain conditions, the particle is attracted to a stable "hovering state" in which it aligns its axis normal to the wall and rests (positions itself) at an activity-dependent distance above the center of the patch.

Active string fluids and gels formed by dipolar active Brownian particles in 3D

Self-propelled particles possessing permanent magnetic dipole moments occur naturally in magnetotactic bacteria and can be built into man-made systems such as active colloids or micro-robots. Yet, the interplay between self-propulsion and anisotropic dipole-dipole interactions on dynamic self-assembly in three dimensions (3D) remains poorly understood. We conduct Brownian dynamics simulations of active dipolar particles in 3D, focusing on the low-density regime, where dipolar hard spheres tend to form chain-like aggregates and percolated networks with increasing dipolar coupling strength. We find that strong active forces override dipolar attractions, effectively inhibiting chain-like aggregation and network formation. Conversely, activating particles with low to moderate forces results in a fluid composed of active chains and rings. At strong dipolar coupling strengths, this active fluid transitions into an active gel, consisting of a percolated network of active chains. Although the overall structure of the active gel remains interconnected, the network experiences more frequent configurational rearrangements due to the reduced bond lifetime of active dipolar particles. Consequently, particles exhibit enhanced translational and rotational diffusion within the active fluid of strings and active gels compared to their passive counterparts. We quantify the influence of activity on aggregate topology as they transition from branched structures to unconnected chains and rings. Our findings are summarized in a state diagram, delineating the impact of dipolar coupling strength and active force magnitude on the system.

Spatiotemporal control of structure and dynamics in a polar active fluid

We apply optimal control theory to a model of a polar active fluid (the Toner-Tu model), with the objective of driving the system into particular emergent dynamical behaviors or programming switching between states on demand. We use the effective self-propulsion speed as the control parameter (i.e. the means of external actuation). We identify control protocols that achieve outcomes such as relocating asters to targeted positions, forcing propagating solitary waves to reorient to a particular direction, and switching between stationary asters and propagating fronts. We analyze the solutions to identify generic principles for controlling polar active fluids. Our findings have implications for achieving spatiotemporal control of active polar systems in experiments, particularly in vitro cytoskeletal systems. Additionally, this research paves the way for leveraging optimal control methods to engineer the structure and dynamics of active fluids more broadly.

Anomalous Fluctuations in a Droplet of Chemically Active Colloids or Enzymes

Chemically active colloids or enzymes cluster into dense droplets driven by their phoretic response to collectively generated chemical gradients. Employing Brownian dynamics simulation techniques, our study of the dynamics of such a chemically active droplet uncovers a rich variety of structures and dynamical properties, including the full range of fluidlike to solidlike behavior, and non-Gaussian positional fluctuations. Our work sheds light on the complex dynamics of the active constituents of metabolic clusters, which are the main drivers of nonequilibrium activity in living systems.

Clustering induces switching between phoretic and osmotic propulsion in active colloidal rafts

Active particles driven by chemical reactions are the subject of intense research to date due to their rich physics, being intrinsically far from equilibrium, and their multiple technological applications. Recent attention in this field is now shifting towards exploring the fascinating dynamics of active and passive mixtures. Here we realize active colloidal rafts, composed of a single catalytic particle encircled by several shells of passive microspheres, and assembled via light-activated chemophoresis. We show that the cluster propulsion mechanism transits from diffusiophoretic to diffusioosmotic as the number of colloidal shells increases. Using the Lorentz reciprocal theorem, we demonstrate that in large clusters self-propulsion emerges by considering the hydrodynamic flow via the diffusioosmotic response of the substrate. The dynamics in our active colloidal rafts are governed by the interplay between phoretic and osmotic effects. Thus, our work highlights their importance in understanding the rich physics of active catalytic systems.

Chemotactic particles as strong electrolytes: Debye-Hückel approximation and effective mobility law

We consider a binary mixture of chemically active particles that produce or consume solute molecules and that interact with each other through the long-range concentration fields they generate. We analytically calculate the effective phoretic mobility of these particles when the mixture is submitted to a constant, external concentration gradient, at leading order in the overall concentration. Relying on an analogy with the modeling of strong electrolytes, we show that the effective phoretic mobility decays with the square root of the concentration: our result is, therefore, a nonequilibrium counterpart to the celebrated Kohlrausch and Debye-Hückel-Onsager conductivity laws for electrolytes, which are extended here to particles with long-range nonreciprocal interactions. The effective mobility law we derive reveals the existence of a regime of maximal mobility and could find applications in the description of nanoscale transport phenomena in living cells.

Pair Interactions of Self-Propelled SiO2-Pt Janus Colloids: Chemically Mediated Encounters

Driven by the necessity to achieve a thorough comprehension of the bottom-up fabrication process of functional materials, this experimental study investigates the pairwise interactions or collisions between chemically active SiO2-Pt Janus colloids. These collisions are categorized based on the Janus colloids' orientations before and after they make physical contact. In addition to the hydrodynamic interactions, the Janus colloids are also known to affect each other's chemical field, resulting in chemophoretic interactions, which depend on the degree of surface anisotropy in reactivity of Janus colloid and the solute-surface interaction at play. Our study reveals that these interactions lead to a noticeable decrease in particle speed and changes in orientation that correlate with the contact duration and yield different collision types. Distinct configurations of contact during collisions were found, whose mechanisms and likelihood are found to be dependent primarily on the chemical interactions. Such estimates of collision and their characterization in dilute suspensions shall have a key impact in determining the arrangement and time scales of dynamical structures and assemblies of denser suspensions and potentially the functional materials of the future.

Controlling active turbulence by activity patterns

By patterning activity in space, one can control active turbulence. To show this, we use Doi's hydrodynamic equations of a semidilute solution of active rods. A linear stability analysis reveals the resting isotropic fluid to be unstable above an absolute pusher activity. The emergent activity-induced paranematic state displays active turbulence, which we characterize by different quantities including the energy spectrum, which shows the typical power-law decay with exponent -4. Then, we control the active turbulence by a square lattice of circular spots where activity is switched off. In the parameter space lattice constant versus surface-to-surface distance of the spots, we identify different flow states. Most interestingly, for lattice constants below the vorticity correlation length and for spot distances smaller than the nematic coherence length, we observe a multi-lane flow state, where flow lanes with alternating flow directions are separated by a street of vortices. The flow pattern displays pronounced multistability and also appears transiently at the transition to the isotropic active-turbulence state. At larger lattice constants a trapped vortex state is identified with a non-Gaussian vorticity distribution due to the low flow vorticity at the spots. It transitions to conventional active turbulence for increasing spot distance.

Self-organization of active colloids mediated by chemical interactions

Self-propelled colloidal particles exhibit rich non-equilibrium phenomena and have promising applications in fields such as drug delivery and self-assembled active materials. Previous experimental and theoretical studies have shown that chemically active colloids that consume or produce a chemical can self-organize into clusters with diverse characteristics depending on the effective phoretic interactions. In this paper, we investigate self-organization in systems with multiple chemical species that undergo a network of reactions and multiple colloidal species that participate in different reactions. Active colloids propelled by complex chemical reactions with potentially nonlinear kinetics can be realized using enzymatic reactions that occur on the surface of enzyme-coated particles. To demonstrate how the self-organizing behavior depends on the chemical reactions active colloids catalyze and their chemical environment, we consider first a single type of colloid undergoing a simple catalytic reaction, and compare this often-studied case with self-organization in binary mixtures of colloids with sequential reactions, and binary mixtures with nonlinear autocatalytic reactions. Our results show that in general active colloids at low particle densities can form localized clusters in the presence of bulk chemical reactions and phoretic attractions. The characteristics of the clusters, however, depend on the reaction kinetics in the bulk and on the particles and phoretic coefficients. With one or two chemical species that only undergo surface reactions, the space for possible self-organizations are limited. By considering the additional system parameters that enter the chemical reaction network involving reactions on the colloids and in the fluid, the design space of colloidal self-organization can be enlarged, leading to a variety of non-equilibrium structures.

Network Effects Lead to Self-Organization in Metabolic Cycles of Self-Repelling Catalysts

Mixtures of particles that interact through phoretic effects are known to aggregate if they belong to species that exhibit attractive self-interactions. We study self-organization in a model metabolic cycle composed of three species of catalytically active particles that are chemotactic toward the chemicals that define their connectivity network. We find that the self-organization can be controlled by the network properties, as exemplified by a case where a collapse instability is achieved by design for self-repelling species. Our findings highlight a possibility for controlling the intricate functions of metabolic networks by taking advantage of the physics of phoretic active matter.

Chemotactic Motility-Induced Phase Separation

Collectives of actively moving particles can spontaneously separate into dilute and dense phases-a fascinating phenomenon known as motility-induced phase separation (MIPS). MIPS is well-studied for randomly moving particles with no directional bias. However, many forms of active matter exhibit collective chemotaxis, directed motion along a chemical gradient that the constituent particles can generate themselves. Here, using theory and simulations, we demonstrate that collective chemotaxis strongly competes with MIPS-in some cases, arresting or completely suppressing phase separation, or in other cases, generating fundamentally new dynamic instabilities. We establish principles describing this competition, thereby helping to reveal and clarify the rich physics underlying active matter systems that perform chemotaxis, ranging from cells to robots.

Alignment interaction and band formation in assemblies of autochemorepulsive walkers

Chemotaxis refers to the motion of an organism induced by chemical stimuli and is a motility mode shared by many living species that has been developed by evolution to optimize certain biological processes such as foraging or immune response. In particular, autochemotaxis refers to chemotaxis mediated by a cue produced by the chemotactic particle itself. Here, we investigate the collective behavior of autochemotactic particles that are repelled by the cue and therefore migrate preferentially towards low-concentration regions. To this end, we introduce a lattice model inspired by the true self-avoiding walk which reduces to the Keller-Segel model in the continuous limit, for which we describe the rich phase behavior. We first rationalize the chemically mediated alignment interaction between walkers in the limit of stationary concentration fields, and then describe the various large-scale structures that can spontaneously form and the conditions for them to emerge, among which we find stable bands traveling at constant speed in the direction transverse to the band.

Dynamics of Motility-Induced Clusters: Coarsening beyond Ostwald Ripening

We study the dynamics of clusters of active Brownian disks generated by motility-induced phase separation, by applying an algorithm that we devised to track cluster trajectories. We identify an aggregation mechanism that goes beyond Ostwald ripening but also yields a dynamic exponent characterizing the cluster growth z=3, in the timescales explored numerically. Clusters of mass M self-propel with enhanced diffusivity D∼Pe^{2}/sqrt[M]. Their fast motion drives aggregation into large fractal structures, which are patchworks of diverse hexatic orders, and coexist with regular, orientationally uniform, smaller ones. To bring out the impact of activity, we perform a comparative study of a passive system that evidences major differences with the active case.

Self-organization of primitive metabolic cycles due to non-reciprocal interactions

One of the greatest mysteries concerning the origin of life is how it has emerged so quickly after the formation of the earth. In particular, it is not understood how metabolic cycles, which power the non-equilibrium activity of cells, have come into existence in the first instances. While it is generally expected that non-equilibrium conditions would have been necessary for the formation of primitive metabolic structures, the focus has so far been on externally imposed non-equilibrium conditions, such as temperature or proton gradients. Here, we propose an alternative paradigm in which naturally occurring non-reciprocal interactions between catalysts that can partner together in a cyclic reaction lead to their recruitment into self-organized functional structures. We uncover different classes of self-organized cycles that form through exponentially rapid coarsening processes, depending on the parity of the cycle and the nature of the interaction motifs, which are all generic but have readily tuneable features.

Hydrodynamic interactions between squirmers near walls: far-field dynamics and near-field cluster stability

Confinement increases contacts between microswimmers in dilute suspensions and affects their interactions. In particular, boundaries have been shown experimentally to lead to the formation of clusters that would not occur in bulk fluids. To what extent does hydrodynamics govern these boundary-driven encounters between microswimmers? We consider theoretically the symmetric boundary-mediated encounters of model microswimmers under gravity through far-field interaction of a pair of weak squirmers, as well as the lubrication interactions occurring after contact between two or more squirmers. In the far field, the orientation of microswimmers is controlled by the wall and the squirming parameter. The presence of a second swimmer influences the orientation of the original squirmer, but for weak squirmers, most of the interaction occurs after contact. We thus analyse next the near-field reorientation of circular groups of squirmers. We show that a large number of swimmers and the presence of gravity can stabilize clusters of pullers, while the opposite is true for pushers; to be stable, clusters of pushers thus need to be governed by other interactions (e.g. phoretic). This simplified approach to the phenomenon of active clustering enables us to highlight the hydrodynamic contribution, which can be hard to isolate in experimental realizations.

The interplay between chemo-phoretic interactions and crowding in active colloids

Many motile microorganisms communicate with each other and their environments via chemical signaling which leads to long-range interactions mediated by self-generated chemical gradients. However, consequences of the interplay between crowding and chemotactic interactions on their collective behavior remain poorly understood. In this work, we use Brownian dynamics simulations to investigate the effect of packing fraction on the formation of non-equilibrium structures in a monolayer of diffusiophoretic self-propelled colloids as a model for chemically active particles. Focusing on the case when a chemical field induces attractive positional and repulsive orientational interactions, we explore dynamical steady-states of active colloids of varying packing fractions and degrees of motility. In addition to collapsed, active gas, and dynamical clustering steady-states reported earlier for low packing fractions, a new phase-separated state emerges. The phase separation results from a competition between long-range diffusiophoretic interactions and motility and is observed at moderate activities and a wide range of packing fractions. Our analysis suggests that the fraction of particles in the largest cluster is a suitable order parameter for capturing the transition from an active gas and dynamical clustering states to a phase-separated state.

Kinetic modeling of the chemotactic process in run-and-tumble bacteria

The chemotactic process of run-and-tumble bacteria results from modulating the tumbling rate in response to changes in chemoattractant gradients felt by the bacteria. The response has a characteristic memory time and is subject to important fluctuations. These ingredients are considered in a kinetic description of chemotaxis, allowing the computation of the stationary mobility and the relaxation times needed to reach the steady state. For large memory times, these relaxation times become large, implying that finite-time measurements give rise to nonmonotonic currents as a function of the imposed chemoattractant gradient, contrary to the stationary regime where the response is monotonic. The case of an inhomogeneous signal is analyzed. Contrary to the usual Keller-Segel model, the response is nonlocal, and the bacterial profile is smoothed with a characteristic length that grows with the memory time. Finally, the case of traveling signals is considered, where appreciable differences appear compared to memoryless chemotactic descriptions.

Collective motion of active particles exhibiting non-reciprocal orientational interactions

We present a Brownian dynamics study of a 2D bath of active particles interacting among each other through usual steric interactions and, additionally, via non-reciprocal avoidant orientational interactions. We motivate them by the fact that the two flagella of the alga Chlamydomonas interact sterically with nearby surfaces such that a torque acts on the alga. As expected, in most cases such interactions disrupt the motility-induced particle clustering in active baths. Surprisingly, however, we find that the active particles can self-organize into collectively moving flocks if the range of non-reciprocal interactions is close to that of steric interactions. We observe that the flocking motion can manifest itself through a variety of structural forms, spanning from single dense bands to multiple moderately-dense stripes, which are highly dynamic. The flocking order parameter is found to be only weakly dependent on the underlying flock structure. Together with the variance of the local-density distribution, one can clearly group the flocking motion into the two separate band and dynamic-stripes states.

Clustering of self-thermophilic asymmetric dimers: the relevance of hydrodynamics

Self-thermophilic dimers are characterized by a net phoretic attraction which, in combination with hydrodynamic interactions, results in the formation of crystalline-like aggregates. To distinguish the effect of the different contributions is frequently an important challenge. We present a simulation investigation done in parallel in the presence and the absence of hydrodynamic interactions for the case of asymmetric self-thermophoretic dimers. In the absence of hydrodynamics, the clusters have the standard heads-in configurations. In contrast, in the presence of hydrodynamics, clusters with heads-in conformation are being formed, in which dimers with their propulsion velocity pointing out of the cluster are assembled and stabilized by strong hydrodynamic osmotic flows. Significant variation in the material properties is to be expected from such differences in the collective behavior, whose understanding and control is of great relevance for the development of new synthetic active materials.

Analytic Solution of an Active Brownian Particle in a Harmonic Well

We provide an analytical solution for the time-dependent Fokker-Planck equation for a two-dimensional active Brownian particle trapped in an isotropic harmonic potential. Using the passive Brownian particle as basis states we show that the Fokker-Planck operator becomes lower diagonal, implying that the eigenvalues are unaffected by the activity. The propagator is then expressed as a combination of the equilibrium eigenstates with weights obeying exact iterative relations. We show that for the low-order correlation functions, such as the positional autocorrelation function, the recursion terminates at finite order in the Péclet number, allowing us to generate exact compact expressions and derive the velocity autocorrelation function and the time-dependent diffusion coefficient. The nonmonotonic behavior of latter quantities serves as a fingerprint of the nonequilibrium dynamics.

Nonequilibrium thermodynamics of Janus particle self-assembly

We compute the energetic cost of formation of Janus particle structures. Using an approach that couples particle dynamics to the evolution of fuel concentration in the medium, which we consider to be initially inhomogeneous, we show the different types of emerging structures. The energy dissipated in the formation of such structures is obtained from the entropy production rate, which is a non-monotonic function of the fraction of assembled particles and, thus, different in each self-assembly regime. An analysis of the free energy of these particles allows us to establish a thermodynamic criterion of structure formation based on the behavior of chemical potential as a function of the fraction of assembled particles.

Chemotactic self-caging in active emulsions

The out-of-equilibrium dynamics of chemotactic active matter—be it animate or inanimate—is closely coupled to the environment, a chemical landscape shaped by secretions from the motile agents, fuel uptake, or autochemotactic signaling. This gives rise to complex collective effects, which can be exploited by the agents for colony migration strategies or pattern formation. We study such effects using an idealized experimental system: self-propelled microdroplets that communicate via chemorepulsive trails. We present a comprehensive experimental analysis that involves direct probing of the diffusing chemical trails and the trail–droplet interactions and use it to construct a generic theoretical model. We connect these repulsive autochemotactic interactions to the collective dynamics in emulsions, demonstrating a state of dynamical arrest: chemotactic self-caging.

A common feature of biological self-organization is how active agents communicate with each other or their environment via chemical signaling. Such communications, mediated by self-generated chemical gradients, have consequences for both individual motility strategies and collective migration patterns. Here, in a purely physicochemical system, we use self-propelling droplets as a model for chemically active particles that modify their environment by leaving chemical footprints, which act as chemorepulsive signals to other droplets. We analyze this communication mechanism quantitatively both on the scale of individual agent–trail collisions as well as on the collective scale where droplets actively remodel their environment while adapting their dynamics to that evolving chemical landscape. We show in experiment and simulation how these interactions cause a transient dynamical arrest in active emulsions where swimmers are caged between each other’s trails of secreted chemicals. Our findings provide insight into the collective dynamics of chemically active particles and yield principles for predicting how negative autochemotaxis shapes their navigation strategy.

Self-phoretic Brownian dynamics simulations

A realistic and effective model to simulate phoretic Brownian dynamics swimmers based on the general form of the thermophoretic force is here presented. The collective behavior of self-phoretic dimers is investigated with this model and compared with two simpler versions, allowing the understanding of the subtle interplay of steric interactions, propulsion, and phoretic effects. The phoretic Brownian dynamics method has control parameters which can be tuned to closely map the properties of experiments or simulations with explicit solvent, in particular those performed with multiparticle collision dynamics. The combination of the phoretic Brownian method and multiparticle collision dynamics is a powerful tool to precisely identify the importance of hydrodynamic interactions in systems of self-phoretic swimmers.

Generic Rules for Distinguishing Autophoretic Colloidal Motors

Distinguishing the operating mechanisms of nano- and micromotors powered by chemical gradients, i.e. "autophoresis", holds the key for fundamental and applied reasons. In this article, we propose and experimentally confirm that the speeds of a self-diffusiophoretic colloidal motor scale inversely to its population density but not for self-electrophoretic motors, because the former is an ion source and thus increases the solution ionic strength over time while the latter does not. They also form clusters in visually distinguishable and quantifiable ways. This pair of rules is simple, powerful, and insensitive to the specific material composition, shape or size of a colloidal motor, and does not require any measurement beyond typical microscopy. These rules are not only useful in clarifying the operating mechanisms of typical autophoretic micromotors, but also in predicting the dynamics of unconventional ones that are yet to be experimentally realized, even those involving enzymes.

Active particles in explicit solvent: Dynamics of clustering for alignment interaction

We study the dynamics of clustering in systems containing active particles that are immersed in an explicit solvent. For this, we have adopted a hybrid simulation method, consisting of molecular dynamics and multiparticle collision dynamics. In our model, the overlap-avoiding passive interaction of an active particle with another active particle or a solvent particle has been taken care of via variants of the Lennard-Jones potential. Dynamic interactions among the active particles have been incorporated via a Vicsek-like alignment rule in self-propulsion that facilitates clustering. We quantify the effects of activity and importance of hydrodynamics on the dynamics of clustering via variations of relevant system parameters. Results are obtained for low overall density of active particles, for which the state point is close to the vapor branch of the coexistence curve, and thus the morphology consists of disconnected clusters. In such a situation, the mechanism of growth switches among particle diffusion, diffusive coalescence, and ballistic aggregation, depending upon the presence or absence of active and hydrodynamic interactions providing different kinds of mobilities to the clusters. Corresponding growth laws have been quantified and discussed in the context of appropriate theoretical pictures. Our results suggest that multiparticle collision dynamics is an effective method for the investigation of hydrodynamic phenomena in phase-separating active matter systems.

Hunting active Brownian particles: Learning optimal behavior

We numerically study active Brownian particles that can respond to environmental cues through a small set of actions (switching their motility and turning left or right with respect to some direction) which are motivated by recent experiments with colloidal self-propelled Janus particles. We employ reinforcement learning to find optimal mappings between the state of particles and these actions. Specifically, we first consider a predator-prey situation in which prey particles try to avoid a predator. Using as reward the squared distance from the predator, we discuss the merits of three state-action sets and show that turning away from the predator is the most successful strategy. We then remove the predator and employ as collective reward the local concentration of signaling molecules exuded by all particles and show that aligning with the concentration gradient leads to chemotactic collapse into a single cluster. Our results illustrate a promising route to obtain local interaction rules and design collective states in active matter.

Interactions in active colloids

The past two decades have seen a remarkable progress in the development of synthetic colloidal agents which are capable of creating directed motion in an unbiased environment at the microscale. These self-propelling particles are often praised for their enormous potential to self-organize into dynamic nonequilibrium structures such as living clusters, synchronized super-rotor structures or self-propelling molecules featuring a complexity which is rarely found outside of the living world. However, the precise mechanisms underlying the formation and dynamics of many of these structures are still barely understood, which is likely to hinge on the gaps in our understanding of how active colloids interact. In particular, besides showing comparatively short-ranged interactions which are well known from passive colloids (Van der Waals, electrostatic etc), active colloids show novel hydrodynamic interactions as well as phoretic and substrate-mediated 'osmotic' cross-interactions which hinge on the action of the phoretic field gradients which are induced by the colloids on other colloids in the system. The present article discusses the complexity and the intriguing properties of these interactions which in general are long-ranged, non-instantaneous, non-pairwise and non-reciprocal and which may serve as key ingredients for the design of future nonequilibrium colloidal materials. Besides providing a brief overview on the state of the art of our understanding of these interactions a key aim of this review is to emphasize open key questions and corresponding open challenges.

Multi-particle collision dynamics with a non-ideal equation of state. II. Collective dynamics of elongated squirmer rods

Simulations of flow fields around microscopic objects typically require methods that both solve the Navier-Stokes equations and also include thermal fluctuations. One such method popular in the field of soft-matter physics is the particle-based simulation method of multi-particle collision dynamics (MPCD). However, in contrast to the typically incompressible real fluid, the fluid of the traditional MPCD methods obeys the ideal-gas equation of state. This can be problematic because most fluid properties strongly depend on the fluid density. In a recent article, we proposed an extended MPCD algorithm and derived its non-ideal equation of state and an expression for the viscosity. In the present work, we demonstrate its accuracy and efficiency for the simulations of the flow fields of single squirmers and of the collective dynamics of squirmer rods. We use two exemplary squirmer-rod systems for which we compare the outcome of the extended MPCD method to the well-established MPCD version with an Andersen thermostat. First, we explicitly demonstrate the reduced compressibility of the MPCD fluid in a cluster of squirmer rods. Second, for shorter rods, we show the interesting result that in simulations with the extended MPCD method, dynamic swarms are more pronounced and have a higher polar order. Finally, we present a thorough study of the state diagram of squirmer rods moving in the center plane of a Hele-Shaw geometry. From a small to large aspect ratio and density, we observe a disordered state, dynamic swarms, a single swarm, and a jammed cluster, which we characterize accordingly.

Crystallisation and polymorph selection in active Brownian particles

We explore crystallisation and polymorph selection in active Brownian particles with numerical simulation. In agreement with previous work (Wysocki et al. in Europhys Lett 105:48004, 2014), we find that crystallisation is suppressed by activity and occurs at higher densities with increasing Péclet number ([Formula: see text]). While the nucleation rate decreases with increasing activity, the crystal growth rate increases due to the accelerated dynamics in the melt. As a result of this competition, we observe the transition from a nucleation and growth regime at high [Formula: see text] to "spinodal nucleation" at low [Formula: see text]. Unlike the case of passive hard spheres, where preference for FCC over HCP polymorphs is weak, activity causes the annealing of HCP stacking faults, thus strongly favouring the FCC symmetry at high [Formula: see text]. When freezing occurs more slowly, in the nucleation and growth regime, this tendency is much reduced and we see a trend towards the passive case of little preference for either polymorph.

Optimal Non-Markovian Search Strategies with n-Step Memory

Stochastic search processes are ubiquitous in nature and are expected to become more efficient when equipped with a memory, where the searcher has been before. A natural realization of a search process with long-lasting memory is a migrating cell that is repelled from the diffusive chemotactic signal that it secretes on its way, denoted as an autochemotactic searcher. To analyze the efficiency of this class of non-Markovian search processes, we present a general formalism that allows one to compute the mean first-passage time (MFPT) for a given set of conditional transition probabilities for non-Markovian random walks on a lattice. We show that the optimal choice of the n-step transition probabilities decreases the MFPT systematically and substantially with an increasing number of steps. It turns out that the optimal search strategies can be reduced to simple cycles defined by a small parameter set and that mirror-asymmetric walks are more efficient. For the autochemotactic searcher, we show that an optimal coupling between the searcher and the chemical reduces the MFPT to 1/3 of the one for a Markovian random walk.

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.

Mesoscopic clustering behavior of thermophoretic-type active particle suspension under quasi-one-dimensional microconfinement

We experimentally investigate the mesoscopic clustering behavior of thermophoretic-type active particle suspension under quasi-one-dimensional spatial confinement (high aspect ratio microchannel). The microchannel enhances the viscous dissipation to operate the system in subpropulsion regime. We find that, in the subpropulsion regime, the steady-state configuration of active particle suspension exhibits a transition from homogeneous state to sausagelike clustering bundle located at the channel center, quasiperiodic isolated clusters at the channel center, aperiodic isolated cluster deviated from channel center, and finally to the typical propulsion-induced accumulation around the channel boundary as increasing the excitation laser intensity. The formation of those patterns is under the interplay of outward-pointing mesoscopic scaled thermophoretic force and the applied spatial confinement. The finding of those special patterns can provide some further possibilities of particle manipulation at mesoscopic scale.

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.

Collective behavior of thermophoretic dimeric active colloids in three-dimensional bulk

Colloids driven by phoresis constitute one of the main avenues for the design of synthetic microswimmers. For these swimmers, the specific form of the phoretic and hydrodynamic interactions dramatically influences their dynamics. Explicit solvent simulations allow the investigation of the different behaviors of dimeric Janus active colloids. The phoretic character is modified from thermophilic to thermophobic, and this, together with the relative size of the beads, strongly influences the resulting solvent velocity fields. Hydrodynamic flows can change from puller-type to pusher-type, although the actual flows significantly differ from these standard flows. Such hydrodynamic interactions combined with phoretic interactions between dimers result in several interesting phenomena in three-dimensional bulk conditions. Thermophilic dimeric swimmers are attracted to each other and form large and stable aggregates. Repulsive phoretic interactions among thermophobic dimeric swimmers hinder such clustering and lead, together with long- and short-ranged attractive hydrodynamic interactions, to short-lived, aligned swarming structures.

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.

Disorder-Induced Long-Ranged Correlations in Scalar Active Matter

We study the impact of quenched random potentials and torques on scalar active matter. Microscopic simulations reveal that motility-induced phase separation is replaced in two dimensions by an asymptotically homogeneous phase with anomalous long-ranged correlations and nonvanishing steady-state currents. Using a combination of phenomenological models and a field-theoretical treatment, we show the existence of a lower-critical dimension d_{c}=4, below which phase separation is only observed for systems smaller than an Imry-Ma length scale. We identify a weak-disorder regime in which the structure factor scales as S(q)∼1/q^{2}, which accounts for our numerics. In d=2, we predict that, at larger scales, the behavior should cross over to a strong-disorder regime. In d>2, these two regimes exist separately, depending on the strength of the potential.

Target Search of Active Agents Crossing High Energy Barriers

Target search by active agents in rugged energy landscapes has remained a challenge because standard enhanced sampling methods do not apply to irreversible dynamics. We overcome this nonequilibrium rare-event problem by developing an algorithm generalizing transition-path sampling to active Brownian dynamics. This method is exemplified and benchmarked for a paradigmatic two-dimensional potential with a high barrier. We find that even in such a simple landscape the structure and kinetics of the ensemble of transition paths changes drastically in the presence of activity. Indeed, active Brownian particles reach the target more frequently than passive Brownian particles, following longer and counterintuitive search patterns.

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.

Active Assembly of Spheroidal Photocatalytic BiVO4 Microswimmers

We create single-component photocatalytic bismuth vanadate (BiVO₄) microswimmers with a spheroidal shape that move individually upon irradiation without any asymmetrization step. These particles form active assemblies which we investigate combining an experimental approach with numerical simulations and analytical calculations. We systematically explore the speed and assembly of the swimmers into clusters of up to four particles and find excellent agreement between experiment and theory, which helps us to understand motion patterns and speed trends. Moreover, different batches of particles can be functionalized separately, making them ideal candidates to fulfill a multitude of tasks, such as sensing or environmental remediation. To exemplify this, we coat our swimmers with silica (SiO₂) and selectively couple some of their modules to fluorophores in a way which does not inhibit self-propulsion. The present work establishes spheroidal BiVO₄ microswimmers as a versatile platform to design multifunctional microswimmers.

Unravelling the role of phoretic and hydrodynamic interactions in active colloidal suspensions

Active fluids comprise a variety of systems composed of elements immersed in a fluid environment which can convert some form of energy into directed motion; as such they are intrinsically out-of-equilibrium in the absence of any external force. A fundamental problem in the physics of active matter concerns the understanding of how the characteristics of autonomous propulsion and agent-agent interactions determine the collective dynamics of the system. We study numerically the suspensions of self-propelled diffusiophoretic colloids, in (quasi)-2d configurations, accounting for both dynamically resolved solute-mediated phoretic interactions and solvent-mediated hydrodynamic interactions. Our results show that the system displays different scenarios at changing the colloid-solute affinity and it develops a cluster phase in the chemoattractive case. We study the statistics of cluster sizes and cluster morphologies for different magnitudes of colloidal activity. Finally, we provide evidences that hydrodynamics plays a relevant role in the aggregation kinetics and cluster morphology, significantly hindering cluster growth.

Competing active and passive interactions drive amoebalike crystallites and ordered bands in active colloids

Swimmers and self-propelled particles are physical models for the collective behavior and motility of a wide variety of living systems, such as bacteria colonies, bird flocks, and fish schools. Such artificial active materials are amenable to physical models which reveal the microscopic mechanisms underlying the collective behavior. Here we study colloids in a dc electric field. Our quasi-two-dimensional system of electrically driven particles exhibits a rich and exotic phase behavior exhibiting passive crystallites, motile crystallites, an active gas, and banding. Amongst these are two mesophases, reminiscent of systems with competing interactions. At low field strengths activity suppresses demixing, leading to motile crystallites. Meanwhile, at high field strengths, activity drives partial demixing to traveling bands. We parametrize a particulate simulation model which reproduces the experimentally observed phases.

Active Brownian Motion with Orientation-Dependent Motility: Theory and Experiments

Combining experiments on active colloids, whose propulsion velocity can be controlled via a feedback loop, and the theory of active Brownian motion, we explore the dynamics of an overdamped active particle with a motility that depends explicitly on the particle orientation. In this case, the active particle moves faster when oriented along one direction and slower when oriented along another, leading to anisotropic translational dynamics which is coupled to the particle's rotational diffusion. We propose a basic model of active Brownian motion for orientation-dependent motility. On the basis of this model, we obtain analytical results for the mean trajectories, averaged over the Brownian noise for various initial configurations, and for the mean-square displacements including their non-Gaussian behavior. The theoretical results are found to be in good agreement with the experimental data. Orientation-dependent motility is found to induce significant anisotropy in the particle displacement, mean-square displacement, and non-Gaussian parameter even in the long-time limit. Our findings establish a methodology for engineering complex anisotropic motilities of active Brownian particles, with a potential impact in the study of the swimming behavior of microorganisms subjected to anisotropic driving fields.

Shape- and Material-Dependent Self-Propulsion of Photocatalytic Active Colloids, Interfacial Effects, and Dynamic Interparticle Interactions

Active colloids powered by self-generated, local chemical concentration gradients exhibit dynamics that are a function of the particles' morphology and material properties. These characteristics also govern how the active colloids interact with surfaces, including other particles and nearby walls. Thus, by targeted design, the dynamic behavior, on average, can be engineered, despite a lack of "external" control such as an applied magnetic field. This allows for the development of new applications and the investigation of novel effects that arise when self-propelled active colloids have complex shapes and material composition. Here, we explore some of our recent work on this topic including the dynamics and interactions of photoactivated, self-propelled colloids with such multifaceted properties. We also delve into some special cases, such as a new variety of active particle-particle interaction that we recently developed, in which direct contact between the active colloids is forbidden, and the direction of propulsion for pairs of particles is correlated. The unifying theme of the work highlighted herein is the relationship between the physical, chemical, and material properties of active colloids and their motive behavior, the understanding of which opens up a wide range of new possibilities as we move toward the ultimate goal of realizing functional, man-made micro- and nanomachinery.

Chemically Active Particles: From One to Few on the Way to Many

Chemically active particles suspended in a liquid solution can achieve self-motility by locally changing the chemical composition of the solution via catalytic reactions at their surfaces. They operate intrinsically out of equilibrium, continuously extracting free energy from the environment to power the dissipative self-motility. The effective interactions involving active particles are, in general, nonreciprocal and anisotropic, even if the particles have simple shapes (e.g., Janus spheres). Accordingly, for chemically active particles a very rich behavior of collective motion and self-assembly may be expected to emerge, including phenomena such as microphase separation in the form of kinetically stable, finite-sized aggregates. Here, I succinctly review a number of recent experimental studies that demonstrate the self-assembly of structures, involving chemically active Janus particles, which exhibit various patterns of motion. These examples illustrate concepts such as "motors made out of motors" (as suggestively named by Fischer [Fischer, P. Nat. Phys. 2018, 14, 1072]). The dynamics of assembly and structure formation observed in these systems can provide benchmark, in-depth testing of the current understanding of motion and effective interactions produced by chemical activity. Finally, one notes that these significant achievements are likely just the beginning of the field. Recently reported particles endowed with time-dependent chemical activity or switchable reaction mechanisms open the way for exciting developments, such as periodic reshaping of self-assembled structures based on man-made internal clocks.

Emergent collective dynamics of bottom-heavy squirmers under gravity

We present the results of hydrodynamic simulations using the method of multi-particle collision dynamics for a system of squirmer microswimmers moving under the influence of gravity at low Reynolds numbers. In addition, the squirmers are bottom-heavy so that they experience a torque which aligns them along the vertical. The squirmers interact hydrodynamically by the flow fields of a stokeslet and rotlet, which are initiated by the acting gravitational force and torque, respectively, and by their own flow fields. By varying the ratio of swimming to bulk sedimentation velocity and the torque, we determine state diagrams for the emergent collective dynamics of neutral squirmers as well as strong pushers and pullers. For low swimming velocity and torque we observe conventional sedimentation, while the sedimentation profile becomes inverted when their values are increased. For neutral squirmers we discover convective rolls of circulating squirmers between both sedimentation states, which sit at the bottom of the system and are fed by plumes made of collectively sinking squirmers. At larger torques porous clusters occur that spawn single squirmers. The two latter states can also occur transiently starting from a uniform squirmer distribution and then disappear in the long-time limit. For strong pushers and pullers only weak plume formation is observed.

Exact Phoretic Interaction of Two Chemically Active Particles

We study the nonequilibrium interaction of two isotropic chemically active particles taking into account the exact near-field chemical interactions as well as hydrodynamic interactions. We identify regions in the parameter space wherein the dynamical system describing the two particles can have a fixed point-a phenomenon that cannot be captured under the far-field approximation. We find that, due to near-field effects, the particles may reach a stable equilibrium at a nonzero gap size or make a complex that can dissociate in the presence of sufficiently strong noise. We explicitly show that the near-field effects originate from a self-generated neighbor-reflected chemical gradient, similar to interactions of a self-propelling phoretic particle and a flat substrate.

Swarm Hunting and Cluster Ejections in Chemically Communicating Active Mixtures

A large variety of microorganisms produce molecules to communicate via complex signaling mechanisms such as quorum sensing and chemotaxis. The biological diversity is enormous, but synthetic inanimate colloidal microswimmers mimic microbiological communication (synthetic chemotaxis) and may be used to explore collective behaviour beyond the one-species limit in simpler setups. In this work we combine particle based and continuum simulations as well as linear stability analyses, and study a physical minimal model of two chemotactic species. We observed a rich phase diagram comprising a “hunting swarm phase”, where both species self-segregate and form swarms, pursuing, or hunting each other, and a “core-shell-cluster phase”, where one species forms a dense cluster, which is surrounded by a (fluctuating) corona of particles from the other species. Once formed, these clusters can dynamically eject their core such that the clusters almost turn inside out. These results exemplify a physical route to collective behaviours in microorganisms and active colloids, which are so-far known to occur only for comparatively large and complex animals like insects or crustaceans.

The mesoscopic collective motion of self-propelling active particle suspension confined in two-dimensional micro-channel

The mesoscopic collective motion of self-propelling active particle suspension confined in high aspect ratio two-dimensional micro-channel is numerically studied through coupled [Formula: see text] equation by considering background thermal fluctuation, inter-particle interaction, self-propulsion and micro-channel confinement. Both the self-propulsion and micro-channel confinement are the factors driving the system away from equilibrium and sustaining heterogeneous motion. In such system, the propulsion induced particle accumulation around the channel walls is a universal phenomenon with spatial heterogeneity, where large fraction of particles are caged inside the accumulated cluster with local oscillation coexisting with few fast propelling particles in the center region. Although the formation mechanism of the induced accumulation is well studied, post the cluster formation, how the cluster evolves and its dynamical properties is rarely discussed. Based on the merits of [Formula: see text] equation, the dynamical evolution of induced accumulation is revealed by particle trajectories. It is found that the induced accumulation can be dissociated through the slow re-orientation process of few jammed particles. By using the idea of force chain network, how the transverse confinement couples the transverse displacement with the longitudinal displacement is evidenced. It is further verified by the statistical measurement of correlation probability between transverse and longitudinal displacements. The suppressed displacements in both directions is the origin leading to the slow dynamics of cluster evolution. Temporally, within the orientational relaxation time, this system exhibits non-trivial anomalous diffusion under the competition between the counter effects of self-propulsion (enhanced diffusion) and micro-channel confinement (suppressed diffusion). Additionally, by considering the orientational coupling, the deep hysteresis of accumulation has been found even for very weak orientational coupling strength.