Predictive local field theory for interacting active Brownian spheres in two spatial dimensions

Statistical properties of microphase and bubbly phase-separated active fluids

In phase-separated active fluids, the Ostwald process can go into reverse, leading to either microphase separation or bubbly phase separation. We show that the latter is formed of two macroscopic regions that are occupied by the homogeneous fluid and by the microphase separated one. Within the microphase-separated fluid, the relative rate of the Ostwald process, coalescence, and nucleation determines whether the size distribution of mesoscopic domains is narrowly peaked or displays a broad range of sizes before attaining a cutoff independent of system size. Our results are obtained via large-scale simulations of a minimal field theory for active phase separation and reproduced by an effective model in which the degrees of freedom are the locations and sizes of the microphase-separated domains.

Influence of anomalous agents on the dynamics of an active system

Swarming behavior in systems of self-propelled particles, whether biological or artificial, has received increased attention in recent years. Here, we show that even a small number of particles with anomalous behavior can change dramatically collective dynamics of the swarming system and can impose unusual behavior and transitions between dynamic states. Our results pave the way to practical approaches and concepts of multiagent dynamics in groups of flocking animals: birds, insects, and fish, i.e., active and living soft matter.

Active crystallization from power functional theory

We address the gas, liquid, and crystal phase behaviors of active Brownian particles in three dimensions. The nonequilibrium force balance at coexistence leads to equality of state functions for which we use power functional approximations. Motility-induced phase separation starts at a critical point and quickly becomes metastable against active freezing for Péclet numbers above a nonequilibrium triple point. The mean swim speed acts as a state variable, similar to the density of depletion agents in colloidal demixing. We obtain agreement with recent simulation results and correctly predict the strength of particle number fluctuations in active fluids.

Orientation-Dependent Propulsion of Active Brownian Spheres: From Self-Advection to Programmable Cluster Shapes

Applications of active particles require a method for controlling their dynamics. While this is typically achieved via direct interventions, indirect interventions based, e.g., on an orientation-dependent self-propulsion speed of the particles, become increasingly popular. In this Letter, we investigate systems of interacting active Brownian spheres in two spatial dimensions with orientation-dependent propulsion using analytical modeling and Brownian dynamics simulations. It is found that the orientation dependence leads to self-advection, circulating currents, and programmable cluster shapes.

Active Brownian particles in external force fields: Field-theoretical models, generalized barometric law, and programmable density patterns

We investigate the influence of external forces on the collective dynamics of interacting active Brownian particles in two as well as three spatial dimensions. Via explicit coarse graining, we derive predictive models, i.e., models that give a direct relation between the models' coefficients and the bare parameters of the system, that are applicable for space- and time-dependent external force fields. We study these models for the cases of gravity and harmonic traps. In particular, we derive a generalized barometric formula for interacting active Brownian particles under gravity that is valid for low to high concentrations and activities of the particles. Furthermore, we show that one can use an external harmonic trap to induce motility-induced phase separation in systems that, without external fields, remain in a homogeneous state. This finding makes it possible to realize programmable density patterns in systems of active Brownian particles. Our analytic predictions are found to be in very good agreement with Brownian dynamics simulations.

Dependence of the acoustic propulsion of nano- and microcones on their orientation and aspect ratio

Recent research revealed the orientation-dependent propulsion of a cone-shaped colloidal particle that is exposed to a planar traveling ultrasound wave. Here, we extend the previous research by considering nano- and microcones with different aspect ratios and studying how the propulsion of a particle depends on its orientation and aspect ratio. We also study how the orientation-averaged propulsion of a cone-shaped particle, which corresponds to an isotropic ultrasound field, depends on its aspect ratio and identify an aspect ratio of 1/2 where the orientation-averaged propulsion is particularly strong. To make our simulation results easier reusable for follow-up research, we provide a corresponding simple analytic representation.

Equation of state for active matter

We characterize the steady states of a suspension of two-dimensional active Brownian particles (ABPs). By approximating the first-order correction to the steady-state probability distribution to lowest order in Peclet number, we show that macroscopic quantities can be calculated in analogous way to equilibrium systems using this probability distribution. We then derive expressions for the macroscopic pressure and position-orientation correlation functions. We check our results by direct comparison with extensive numerical simulations. A key finding is the importance of many-body effective interactions even at very low densities.

Inertia changes evolution of motility-induced phase separation in active matter across particle activity

The effects of inertia in active matter and motility-induced phase separation (MIPS) have attracted growing interest but still remain poorly studied. We studied MIPS behavior in the Langevin dynamics across a broad range of particle activity and damping rate values with molecular dynamic simulations. Here we show that the MIPS stability region across particle activity values consists of several domains separated by discontinuous or sharp changes in susceptibility of mean kinetic energy. These domain boundaries have fingerprints in the system's kinetic energy fluctuations and characteristics of gas, liquid, and solid subphases, such as the number of particles, densities, or the power of energy release due to activity. The observed domain cascade is most stable at intermediate damping rates but loses its distinctness in the Brownian limit or vanishes along with phase separation at lower damping values.

From a microscopic inertial active matter model to the Schrödinger equation

Active field theories, such as the paradigmatic model known as 'active model B+', are simple yet very powerful tools for describing phenomena such as motility-induced phase separation. No comparable theory has been derived yet for the underdamped case. In this work, we introduce active model I+, an extension of active model B+ to particles with inertia. The governing equations of active model I+ are systematically derived from the microscopic Langevin equations. We show that, for underdamped active particles, thermodynamic and mechanical definitions of the velocity field no longer coincide and that the density-dependent swimming speed plays the role of an effective viscosity. Moreover, active model I+ contains an analog of the Schrödinger equation in Madelung form as a limiting case, allowing one to find analoga of the quantum-mechanical tunnel effect and of fuzzy dark matter in active fluids. We investigate the active tunnel effect analytically and via numerical continuation.

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.

Analytical approach to chiral active systems: suppressed phase separation of interacting Brownian circle swimmers

We consider chirality in active systems by exemplarily studying the phase behavior of planar systems of interacting Brownian circle swimmers with a spherical shape. For this purpose, we derive a predictive field theory that is able to describe the collective dynamics of circle swimmers. The theory yields a mapping between circle swimmers and noncircling active Brownian particles and predicts that the angular propulsion of the particles leads to a suppression of their motility-induced phase separation, being in line with recent simulation results. In addition, the theory provides analytical expressions for the spinodal corresponding to the onset of motility-induced phase separation and the associated critical point as well as for their dependence on the angular propulsion of the circle swimmers. We confirm our findings by Brownian dynamics simulations. The agreement between results from theory and simulations is found to be good.

Orientation-Dependent Propulsion of Triangular Nano- and Microparticles by a Traveling Ultrasound Wave

Previous studies on ultrasound-propelled nano- and microparticles have considered only systems in which the particle orientation is perpendicular to the direction of propagation of the ultrasound. However, in future applications of these particles, they will typically be able to attain other orientations. Therefore, using direct acoustofluidic simulations, here we study how the propulsion of triangular nano- and microparticles, which are known to have a particularly efficient acoustic propulsion and are therefore promising candidates for future applications, depends on their orientation relative to the propagation direction of a traveling ultrasound wave. Our results reveal that the propulsion of the particles depends strongly on their orientation relative to the direction of wave propagation and that the particles tend to orient perpendicularly to the wave direction. We also address the orientation-averaged translational and angular velocities of the particles, which correspond to the particles' effective propulsion for an isotropic exposure to ultrasound. Our results allow assessment of how free ultrasound-propelled colloidal particles move in three spatial dimensions and thus constitute an important step toward the realization of envisaged future applications of such particles.

Phase separation of self-propelled disks with ferromagnetic and nematic alignment

We present a comprehensive study of a model system of repulsive self-propelled disks in two dimensions with ferromagnetic and nematic velocity alignment interactions. We characterize the phase behavior of the system as a function of the alignment and self-propulsion strength, featuring orientational order for strong alignment and motility-induced phase separation (MIPS) at moderate alignment but high enough self-propulsion. We derive a microscopic theory for these systems yielding a closed set of hydrodynamic equations from which we perform a linear stability analysis of the homogenous disordered state. This analysis predicts MIPS in the presence of aligning torques. The nature of the continuum theory allows for an explicit quantitative comparison with particle-based simulations, which consistently shows that ferromagnetic alignment fosters phase separation, while nematic alignment does not alter either the nature or the location of the instability responsible for it. In the ferromagnetic case, such behavior is due to an increase of the imbalance of the number of particle collisions along different orientations, giving rise to the self-trapping of particles along their self-propulsion direction. On the contrary, the anisotropy of the pair correlation function, which encodes this self-trapping effect, is not significantly affected by nematic torques. Our work shows the predictive power of such microscopic theories to describe complex active matter systems with different interaction symmetries and sheds light on the impact of velocity-alignment interactions in motility-induced phase separation.

Acoustically propelled nano- and microcones: fast forward and backward motion

We focus on cone-shaped nano- and microparticles, which have recently been found to show particularly strong propulsion when they are exposed to a traveling ultrasound wave, and study based on direct acoustofluidic computer simulations how their propulsion depends on the cones' aspect ratio. The simulations reveal that the propulsion velocity and even its sign are very sensitive to the aspect ratio, where short particles move forward whereas elongated particles move backward. Furthermore, we identify a cone shape that allows for a particularly large propulsion speed. Our results contribute to the understanding of the propulsion of ultrasound-propelled colloidal particles, suggest a method for separation and sorting of nano- and microcones concerning their aspect ratio, and provide useful guidance for future experiments and applications.

Dynamic overlap concentration scale of active colloids

By introducing the notion of a dynamic overlap concentration scale, we identify additional universal features of the mechanical properties of active colloids. We codify these features by recognizing that the characteristic length scale of an active particle's trajectory, the run length, introduces a concentration scale ϕ^{*}. Large-scale simulations of repulsive active Brownian particles (ABPs) confirm that this run-length dependent concentration, the trajectory-space analog of the overlap concentration in polymer solutions, delineates distinct concentration regimes in which interparticle collisions alter particle trajectories. Using ϕ^{*} and concentration scales associated with colloidal jamming, the mechanical equation of state for ABPs collapses onto a set of principal curves that contain several overlooked features. The inclusion of these features qualitatively alters previous predictions of the behavior for active colloids, as we demonstrate by computing the spinodal for a suspension of purely repulsive ABPs. Our findings suggest that dynamic overlap concentration scales should help unravel the behavior of active and driven systems.

Ordering kinetics in the active model B

We undertake a detailed numerical study of the Active Model B proposed by Wittkowski et al., [Nature Commun. 5, 4351 (2014)]2041-172310.1038/ncomms5351. We find that the introduction of activity has a drastic effect on the ordering kinetics. First, the domain growth law shows a crossover from the usual Lifshitz-Slyozov growth law for phase separation (L∼t^{1/3}, where t is the time) to a novel growth law (L∼t^{1/4}) at late times. Second, the equal-time correlation function of the density field exhibits dynamical scaling for a given activity strength λ, but the scaling function depends on λ.

Statistical mechanics of active Ornstein-Uhlenbeck particles

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

Pair correlation of dilute active Brownian particles: From low-activity dipolar correction to high-activity algebraic depletion wings

We study the pair correlation of active Brownian particles at low density using numerical simulations and analytical calculations. We observe a winged pair correlation: While particles accumulate in front of an active particle as expected, the depletion wake consists of two depletion wings. In the limit of soft particles, we obtain a closed equation for the pair correlation, allowing us to characterize the depletion wings. In particular, we unveil two regimes at high activity, where the wings adopt a self-similar profile and decay algebraically. We also perform experiments of self-propelled Janus particles and indeed observe the depletion wings.

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.

Effects of social distancing and isolation on epidemic spreading modeled via dynamical density functional theory

For preventing the spread of epidemics such as the coronavirus disease COVID-19, social distancing and the isolation of infected persons are crucial. However, existing reaction-diffusion equations for epidemic spreading are incapable of describing these effects. In this work, we present an extended model for disease spread based on combining a susceptible-infected-recovered model with a dynamical density functional theory where social distancing and isolation of infected persons are explicitly taken into account. We show that the model exhibits interesting transient phase separation associated with a reduction of the number of infections, and allows for new insights into the control of pandemics.

On the shape-dependent propulsion of nano- and microparticles by traveling ultrasound waves

We address the propulsion mechanism of ultrasound-propelled nano- and microparticles that are exposed to a traveling ultrasound wave. Based on direct computational fluid dynamics simulations, we study the effect of two important aspects of the particle shape on the propulsion: rounded vs. pointed and filled vs. hollow shapes. We also study the flow field generated around such particles. Our results reveal that pointedness leads to an increase of the propulsion speed, whereas it is not significantly affected by hollowness. Furthermore, we show that the flow field near to ultrasound-propelled particles can look similar to the flow field generated by pusher squirmers.

Pair-distribution function of active Brownian spheres in two spatial dimensions: Simulation results and analytic representation

We investigate the full pair-distribution function of a homogeneous suspension of spherical active Brownian particles interacting by a Weeks-Chandler-Andersen potential in two spatial dimensions. The full pair-distribution function depends on three coordinates describing the relative positions and orientations of two particles, the Péclet number specifying the activity of the particles, and their mean packing density. This five-dimensional function is obtained from Brownian dynamics simulations. We discuss its structure taking into account all of its degrees of freedom. In addition, we present an approximate analytic expression for the product of the full pair-distribution function and the interparticle force. We find that the analytic expression, which is typically needed when deriving analytic models for the collective dynamics of active Brownian particles, is in good agreement with the simulation results. The results of this work can thus be expected to be helpful for the further theoretical investigation of active Brownian particles as well as nonequilibrium statistical physics in general.