No Arabic abstract
We study, from first principles, the pressure exerted by an active fluid of spherical particles on general boundaries in two dimensions. We show that, despite the non-uniform pressure along curved walls, an equation of state is recovered upon a proper spatial averaging. This holds even in the presence of pairwise interactions between particles or when asymmetric walls induce ratchet currents, which are accompanied by spontaneous shear stresses on the walls. For flexible obstacles, the pressure inhomogeneities lead to a modulational instability as well as to the spontaneous motion of short semi-flexible filaments. Finally, we relate the force exerted on objects immersed in active baths to the particle flux they generate around them.
Active Brownian particles (ABPs) and Run-and-Tumble particles (RTPs) both self-propel at fixed speed $v$ along a body-axis ${bf u}$ that reorients either through slow angular diffusion (ABPs) or sudden complete randomisation (RTPs). We compare the physics of these two model systems both at microscopic and macroscopic scales. Using exact results for their steady-state distribution in the presence of external potentials, we show that they both admit the same effective equilibrium regime perturbatively that breaks down for stronger external potentials, in a model-dependent way. In the presence of collisional repulsions such particles slow down at high density: their propulsive effort is unchanged, but their average speed along ${bf u}$ becomes $v(rho) < v$. A fruitful avenue is then to construct a mean-field description in which particles are ghost-like and have no collisions, but swim at a variable speed $v$ that is an explicit function or functional of the density $rho$. We give numerical evidence that the recently shown equivalence of the fluctuating hydrodynamics of ABPs and RTPs in this case, which we detail here, extends to microscopic models of ABPs and RTPs interacting with repulsive forces.
We develop a dynamic mean-field theory for polar active particles that interact through a self-generated field, in particular one generated through emitting a chemical signal. While being a form of chemotactic response, it is different from conventional chemotaxis in that particles discontinuously change their motility when the local concentration surpasses a threshold. The resulting coupled equations for density and polarization are linear and can be solved analytically for simple geometries, yielding inhomogeneous density profiles. Specifically, here we consider a planar and circular interface. Our theory thus explains the observed coexistence of dense aggregates with an active gas. There are, however, differences to the more conventional picture of liquid-gas coexistence based on a free energy, most notably the absence of a critical point. We corroborate our analytical predictions by numerical simulations of active particles under confinement and interacting through volume exclusion. Excellent quantitative agreement is reached through an effective translational diffusion coefficient. We finally show that an additional response to the chemical gradient direction is sufficient to induce vortex clusters. Our results pave the way to engineer motility responses in order to achieve aggregation and collective behavior even at unfavorable conditions.
The equilibrium properties of a system of passive diffusing particles in an external magnetic field are unaffected by the Lorentz force. In contrast, active Brownian particles exhibit steady-state phenomena that depend on both the strength and the polarity of the applied magnetic field. The intriguing effects of the Lorentz force, however, can only be observed when out-of-equilibrium density gradients are maintained in the system. To this end, we use the method of stochastic resetting on active Brownian particles in two dimensions by resetting them to the line $x=0$ at a constant rate and periodicity in the $y$ direction. Under stochastic resetting, an active system settles into a nontrivial stationary state which is characterized by an inhomogeneous density distribution, polarization and bulk fluxes perpendicular to the density gradients. We show that whereas for a uniform magnetic field the properties of the stationary state of the active system can be obtained from its passive counterpart, novel features emerge in the case of an inhomogeneous magnetic field which have no counterpart in passive systems. In particular, there exists an activity-dependent threshold rate such that for smaller resetting rates, the density distribution of active particles becomes non-monotonic. We also study the mean first-passage time to the $x$ axis and find a surprising result: it takes an active particle more time to reach the target from any given point for the case when the magnetic field increases away from the axis. The theoretical predictions are validated using Brownian dynamics simulations.
The dynamics of dissipative soft-sphere gases obeys Newtons equation of motion which are commonly solved numerically by (force-based) Molecular Dynamics schemes. With the assumption of instantaneous, pairwise collisions, the simulation can be accelerated considerably using event-driven Molecular Dynamics, where the coefficient of restitution is derived from the interaction force between particles. Recently it was shown, however, that this approach may fail dramatically, that is, the obtained trajectories deviate significantly from the ones predicted by Newtons equations. In this paper, we generalize the concept of the coefficient of restitution and derive a numerical scheme which, in the case of dilute systems and frictionless interaction, allows us to perform highly efficient event-driven Molecular Dynamics simulations even for non-instantaneous collisions. We show that the particle trajectories predicted by the new scheme agree perfectly with the corresponding (force-based) Molecular Dynamics, except for a short transient period whose duration corresponds to the duration of the contact. Thus, the new algorithm solves Newtons equations of motion like force-based MD while preserving the advantages of event-driven simulations.
We study two interacting identical run and tumble particles (RTPs) in one dimension. Each particle is driven by a telegraphic noise, and in some cases, also subjected to a thermal white noise with a corresponding diffusion constant $D$. We are interested in the stationary bound state formed by the two RTPs in the presence of a mutual attractive interaction. The distribution of the relative coordinate $y$ indeed reaches a steady state that we characterize in terms of the solution of a second-order differential equation. We obtain the explicit formula for the stationary probability $P(y)$ of $y$ for two examples of interaction potential $V(y)$. The first one corresponds to $V(y) sim |y|$. In this case, for $D=0$ we find that $P(y)$ contains a delta function part at $y=0$, signaling a strong clustering effect, together with a smooth exponential component. For $D>0$, the delta function part broadens, leading instead to weak clustering. The second example is the harmonic attraction $V(y) sim y^2$ in which case, for $D=0$, $P(y)$ is supported on a finite interval. We unveil an interesting relation between this two-RTP model with harmonic attraction and a three-state single RTP model in one dimension, as well as with a four-state single RTP model in two dimensions. We also provide a general discussion of the stationary bound state, including examples where it is not unique, e.g., when the particles cannot cross due to an additional short-range repulsion.