Search

A cornerstone of the directed motion of microscopic self-propelling particles is an asymmetric particle structure defining a polarity axis along which these tiny machines move. This structural asymmetry ties the orientational Brownian motion to the microswimmers directional motion, limiting their persistence and making the long time motion effectively diffusive. Here, we demonstrate a completely symmetric thermoplasmonic microswimmer, which is propelled by laser-induced self-thermophoresis. The propulsion direction is imprinted externally to the particle by the heating laser position. The orientational Brownian motion, thus, becomes irrelevant for the propulsion, allowing enhanced control over the particles dynamics with almost arbitrary steering capability. We characterize the particle motion in experiments and simulations and also theoretically. The analysis reveals additional noise appearing in these systems, which is conjectured to be relevant for biological systems. Our experimental results show that even very small particles can be precisely controlled, enabling more advanced applications of these micromachines.

From the organization of animal flocks to the emergence of swarming behaviors in bacterial suspension, populations of motile organisms at all scales display coherent collective motion. Recent studies showed that the anisotropic interaction between active particles plays a key role in the phase behaviors. Here we investigate the collective behaviors of based-active Janus particles that experience an anisotropic interaction of which the orientation is opposite to the direction of active force by using Langevin dynamics simulations in two dimensional space. Interestingly, the system shows emergence of collective swarming states upon increasing the total area fraction of particles, which is not observed in systems without anisotropic interaction or activity. The threshold for emergence of swarming states decreases as particle activity or interaction strength increases. We have also performed basic kinetic analysis to reproduce the essential features of the simulation results. Our results demonstrate that anisotropic interactions at the individual level are sufficient to set homogeneous active particles into stable directed motion.

In a numerical experiment, it is found that at least three localized stationary breather type modes can exist in a chain of Rayleigh oscillators with a nonlinear coupling between them: immobile, “slow,” and “fast” dissipative breathers. The dynamics of synchronous motions of chain elements depends on the relation of characteristic time scales of the system related to the frequency of the oscillators and rigidity of the coupling between them; it is multistable.

This paper is devoted to a mathematical proof of the continuous dependence on the initial data for the discrete thermostatted kinetic framework, for all T > 0. This is a versatile model for describing the time-evolution of a biological complex system which is composed by a large number of interacting entities, called active particles, and is subjected to an external force field due to the environment. A thermostat term is introduced in order to keep the 2nd-order moment of the system, corresponding to the physical global activation energy, constant in time. This model is expressed by a system of nonlinear ordinary differential equations with quadratic nonlinearity.

We examine the interplay of motility and information exchange in a model of run-and-tumble active particles where the particle's motility is encoded as a bit of information that can be exchanged upon contact according to the rules of AND and OR logic gates in a circuit. Motile AND particles become non-motile upon contact with a non-motile particle. Conversely, motile OR particles remain motile upon collision with their non-motile counterparts. AND particles that have become non-motile additionally "reawaken", i.e., recover their motility, at a fixed rate $\mu$, as in the SIS (Susceptible, Infected, Susceptible) model of epidemic spreading, where an infected agent can become healthy again, but keeps no memory of the recent infection, hence it is susceptible to a renewed infection. For $\mu=0$, both AND and OR particles relax irreversibly to absorbing states of all non-motile or all motile particles, respectively. The relaxation kinetics is, however, faster for OR particles that remain active throughout the process. At finite $\mu$, the AND dynamics is controlled by the interplay between reawakening and collision rates. The system evolves to a state of all motile particles (an absorbing state in the language of absorbing phase transitions) for $\mu>\mu_c$ and to a mixed state with coexisting motile and non-motile particles (an active state in the language of absorbing phase transitions) for $\mu

Dipolar active particles describe a class of self-propelled, biological or artificial particles equipped with an internal (typically magnetic) dipole moment. Because of the interplay between self-propulsion and dipole-dipole interactions, complex collective behavior is expected to emerge in systems of such particles. Here, we use Brownian dynamics simulations to explore this collective behavior. We focus on the structures that form in small systems in spatial confinement. We quantify the type of structures that emerge and how they depend on the self-propulsion speed and the dipolar (magnetic) strength of the particles. We observe that the dipolar active particles self-assemble into chains and rings. The dominant configuration is quantified with an order parameter for chain and ring formation and shown to depend on the self-propulsion speed and the dipolar magnetic strength of the particles. In addition, we show that the structural configurations are also affected by the confining walls. To that end, we compare different confining geometries and study the impact of a reorienting 'wall torque' upon collisions of a particle with a wall. Our results indicate that dipolar interactions can further enhance the already rich variety of collective behaviors of active particles.

Confinement is known to have profound effects on the collective dynamics of many active systems. Here, we investigate a modeled active system in circular confinement consisting of biased active particles, where the direction of active force deviates a biased angle from the principle orientation of the anisotropic interaction. We find that such particles can spontaneously form a spiral vortex with two concentric and counter-rotating regions near the boundary. The emerged vortex can be measured by the vortex order parameter which shows nonmonotonic dependencies on both the biased angle and the strength of the anisotropic interaction. Our work can provide an understanding of such dynamic behaviors and enable different strategies for designing ordered collective behaviors.

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 are originated from a self-generated neighbor-reflected chemical gradient, similar to interactions of a self-propelling phoretic particle and a flat substrate.

We present a modified conserved lattice gas (CLG) model in two dimensions with a thermal noise parameter, P = exp(−1/T), where T is the effective temperature. In particular, we focus on the interplay of thermal noise P and the total density of particles ρ in the CLG model. The universality class of the case of 0

For a large class of nonequilibrium systems, thermodynamic notions like work, heat, and, in particular, entropy production can be identified on the level of fluctuating dynamical trajectories. Within stochastic thermodynamics various fluctuation theorems relating these quantities have been proven. Their application to experimental systems requires that all relevant mesostates are accessible. Recent advances address the typical situation that only partial, or coarse-grained, information about a system is available. Thermodynamic inference as a general strategy uses consistency constraints derived from stochastic thermodynamics to infer otherwise hidden properties of nonequilibrium systems. An important class in this respect are active particles, for which we resolve the conflicting strategies that have been proposed to identify entropy production. As a paradigm for thermodynamic inference, the thermodynamic uncertainty relation provides a lower bound on the entropy production through measurements of the dispersion of any current in the system. Likewise, it quantifies the cost of precision for biomolecular processes. Generalizations and ramifications allow the inference of, inter alia, model-free upper bounds on the efficiency of molecular motors and of the minimal number of intermediate states in enzymatic networks.

The presence of active forces in various biological and artificial systems may change how those systems behave when forced. We present a minimal model of a suspension of passive or active swimmers driven on the boundaries by time-dependent forcing. In particular, we consider a time-periodic drive from which we determine the linear response functions of the suspension. The meaning of these response functions are interpreted in terms of the storage and dissipation of energy through the particles within the system. We find that while a slowly driven active system responds in a way similar to a passive system with a re-defined diffusion constant, a rapidly driven active system exhibits a novel behavior related to a change in the motoring activity of the particles due to the external drive.

Abstract Lattice models allow for a computationally efficient investigation of motility-induced phase separation (MIPS) compared to off-lattice systems. Simulations are less demanding, and thus, bigger systems can be accessed with higher accuracy and better statistics. In equilibrium, lattice and off-lattice models with comparable interactions belong to the same universality class. Whether concepts of universality also hold for active particles is still a controversial and open question. Here, we examine two recently proposed active lattice systems that undergo MIPS and investigate numerically their critical behavior. In particular, we examine the claim that these systems and MIPS in general belong to the Ising universality class. We also take a more detailed look on the influence and role of rotational diffusion and active velocity in these systems. Graphic Abstract

We study kinetic Monte-Carlo (KMC) descriptions of active particles. By relying on large discrete time steps, KMC algorithms accelerate the relaxational dynamics of active systems towards their steady-state. We show, however, that their continuous-time limit is ill-defined, leading to the vanishing of trademark behaviors of active matter such as the motility-induced phase separation, ratchet effects, as well as to a diverging mechanical pressure. We show how mixing passive steps with active ones regularizes this behavior, leading to a well-defined continuous-time limit. We propose new AKMC algorithms whose continuous-time limits lead to the active dynamics of Active-Ornstein Uhlenbeck, Active Brownian, and Run-and-Tumbles particles.

We present theory and experiments demonstrating the existence of invariant manifolds that impede the motion of microswimmers in two-dimensional fluid flows. One-way barriers are apparent in a hyperbolic fluid flow that block the swimming of both smooth-swimming and run-and-tumble \emph{Bacillus subtilis} bacteria. We identify key phase-space structures, called swimming invariant manifolds (SwIMs), that serve as separatrices between different regions of long-time swimmer behavior. When projected into $xy$-space, the edges of the SwIMs act as one-way barriers, consistent with the experiments., Comment: 11 pages, 8 figures

Suspensions of motile active particles with space-dependent activity form characteristic polarization and density patterns. Recent single-particle studies for planar activity landscapes identified several quantities associated with emergent density-polarization patterns that are solely determined by bulk variables. Naive thermodynamic intuition suggests that these results might hold for arbitrary activity landscapes mediating bulk regions, and thus could be used as benchmarks for simulations and theories. However, the considered system operates in a nonequilibrium steady state and we prove by construction that the quantities in question lose their simple form for curved activity landscapes. Specifically, we provide a detailed analytical study of polarization and density profiles induced by radially symmetric activity steps, and of the total polarization for the case of a general radially symmetric activity landscape. While the qualitative picture is similar to the planar case, all the investigated variables depend not only on bulk variables but also comprise geometry-induced contributions. We verified that all our analytical results agree with exact numerical calculations.

We investigate coarse-grained models of suspended self-thermophoretic microswimmers. Upon heating, the Janus spheres, with hemispheres made of different materials, induce a heterogeneous local solvent temperature that causes the self-phoretic particle propulsion. Starting from atomistic molecular dynamics simulations, we verify the coarse-grained description of the fluid in terms of a local molecular temperature field, and its role for the particle's thermophoretic self-propulsion and hot Brownian motion. The latter is governed by effective nonequilibrium temperatures, which are measured from simulations by confining the particle position and orientation. They are theoretically shown to remain relevant for any further spatial coarse-graining towards a hydrodynamic description of the entire suspension as a homogeneous complex fluid., Comment: 20 pages, 5 figures

We study the dynamics of one-dimensional active particles confined in a double-well potential, focusing on the escape properties of the system, such as the mean escape time from a well. We first consider a single-particle both in near and far-from-equilibrium regimes by varying the persistence time of the active force and the swim velocity. A non-monotonic behavior of the mean escape time is observed with the persistence time of the activity, revealing the existence of an optimal choice of the parameters favoring the escape process. For small persistence times, a Kramers-like formula with an effective potential obtained within the unified colored noise approximation is shown to hold. Instead, for large persistence times, we developed a simple theoretical argument based on the first passage theory, which explains the linear dependence of the escape time with the persistence of the active force. In the second part of the work, we consider the escape on two active particles mutually repelling. Interestingly, the subtle interplay of active and repulsive forces may lead to a correlation between particles, favoring the simultaneous jump across the barrier. This mechanism cannot be observed in the escape process of two passive particles. Finally, we find that in the small persistence regime, the repulsion favors the escape, such as in passive systems, in agreement with our theoretical predictions, while for large persistence times, the repulsive and active forces produce an effective attraction, which hinders the barrier crossing.

Hydrodynamic interactions (HIs) are important in biophysics research because they influence both the collective and the individual behaviour of microorganisms and self-propelled particles. For instance, HIs at the micro-swimmer level determine the attraction or repulsion between individuals, and hence their collective behaviour. Meanwhile, HIs between swimming appendages (e.g. cilia and flagella) influence the emergence of swimming gaits, synchronised bundles and metachronal waves. In this study, we address the issue of HIs between slender filaments separated by a distance larger than their contour length (d>L) by means of asymptotic calculations and numerical simulations. We first derive analytical expressions for the extended resistance matrix of two arbitrarily-shaped rigid filaments as a series expansion in inverse powers of d/L>1. The coefficients in our asymptotic series expansion are then evaluated using two well-established methods for slender filaments, resistive-force theory (RFT) and slender-body theory (SBT), and our asymptotic theory is verified using numerical simulations based on SBT for the case of two parallel helices. The theory captures the qualitative features of the interactions in the regime d/L>1, which opens the path to a deeper physical understanding of hydrodynamically governed phenomena such as inter-filament synchronisation and multiflagellar propulsion. To demonstrate the usefulness of our results, we next apply our theory to the case of two helices rotating side-by-side, where we quantify the dependence of all forces and torques on the distance and phase difference between them. Using our understanding of pairwise HIs, we then provide physical intuition for the case of a circular array of rotating helices. Our theoretical results will be useful for the study of HIs between bacterial flagella, nodal cilia, and slender microswimmers.

One of the major challenges in stochastic thermodynamics is to compute the distributions of stochastic observables for small-scale systems for which fluctuations play a significant role. Hitherto much theoretical and experimental research has focused on systems composed of passive Brownian particles. In this paper, we study the heat fluctuations in a system of interacting active particles. Specifically we consider a one-dimensional harmonic chain of $N$ active Ornstein-Uhlenbeck particles, with the chain ends connected to heat baths of different temperatures. We compute the moment-generating function for the heat flow in the steady state. We employ our general framework to explicitly compute the moment-generating function for two example single-particle systems. Further, we analytically obtain the scaled cumulants for the heat flow for the chain. Numerical Langevin simulations confirm the long-time analytical expressions for first and second cumulants for the heat flow for a two-particle chain., Comment: 23 pages, 7 figures, Major revision

Experimental data on the active Brownian motion of single particles in the RF (radio-frequency) discharge plasma under the influence of thermophoretic force, induced by laser radiation, depending on the material and type of surface of the particle, are presented. Unlike passive Brownian particles, active Brownian particles, also known as micro-swimmers, move directionally. It was shown that different dust particles in gas discharge plasma can convert the energy of a surrounding medium (laser radiation) into the kinetic energy of motion. The movement of the active particle is a superposition of chaotic motion and self-propulsion.

We study the escape problem for interacting, self-propelled particles confined to a disk, where particles can exit through one open slot on the circumference. Within a minimal two-dimensional Vicsek model, we numerically study the statistics of escape events when the self-propelled particles can be in a flocking state. We show that while an exponential survival probability is characteristic for noninteracting self-propelled particles at all times, the interacting particles have an initial exponential phase crossing over to a subexponential late-time behavior. We use a phenomenological model based on nonstationary Poisson processes which includes the Allee effect to explain this subexponential trend and perform numerical simulations for various noise intensities.

I study deterministic dynamics of chiral active particles in two dimensions. Particles are considered as discs interacting with elastic repulsive forces. An ensemble of particles, started from random initial conditions, demonstrates chaotic collisions resulting in their normal diffusion. This chaos is transient, as rather abruptly a synchronous collisionless state establishes. The life time of chaos grows exponentially with the number of particles. External forcing (periodic or chaotic) is shown to facilitate the synchronization transition.

Collective cell migration is crucial to many physiological and pathological processes such as embryo development, wound healing, and cancer invasion. Recent experimental studies have indicated that the active traction forces generated by migrating cells in a fibrous extracellular matrix (ECM) can mechanically remodel the ECM, giving rise to bundlelike mesostructures bridging individual cells. Such fiber bundles also enable long-range propagation of cellular forces, leading to correlated migration dynamics regulated by the mechanical communication among the cells. Motivated by these experimental discoveries, we develop an active-particle model with polarized effective attractions (APPA) to investigate emergent multicellular migration dynamics resulting from ECM-mediated mechanical communications. In particular, the APPA model generalizes the classic active-Brownian-particle (ABP) model by imposing a pairwise polarized attractive force between the particles, which depends on the instantaneous dynamic states of the particles and mimics the effective mutual pulling between the cells via the fiber bundle bridge. The APPA system exhibits enhanced aggregation behaviors compared to the classic ABP system, and the contrast is more apparent at lower particle densities and higher rotational diffusivities. Importantly, in contrast to the classic ABP system where the particle velocities are not correlated for all particle densities, the high-density phase of the APPA system exhibits strong dynamic correlations, which are characterized by the slowly decaying velocity correlation functions with a correlation length comparable to the linear size of the high-density phase domain (i.e., the cluster of particles). The strongly correlated multicellular dynamics predicted by the APPA model is subsequently verified in in vitro experiments using MCF-10A cells. Our studies indicate the importance of incorporating ECM-mediated mechanical coupling among the migrating cells for appropriately modeling emergent multicellular dynamics in complex microenvironments.

Researchers make systems of self-propelled particles produce the same large-scale dynamics as passive-particle systems.

The co-localization of density modulations and particle polarization is a characteristic emergent feature of motile active matter in activity gradients. It can therefore play the role of a smoking gun for the mesoscale detection of intrinsic microscopic activity. We employ the active-Brownian-particle (ABP) model to derive precise analytical expressions for the density and polarization profiles of a single Janus-type swimmer in the vicinity of an abrupt activity step, including situations with an orientation-dependent propulsion speed. Our results agree well with measurement data for a hot microswimmer presented in the companion paper [1] and can serve as a template for more complex applications, e.g., to motility-induced phase separation or studies of physical boundaries. The essential physics behind our formal results is robustly captured and elucidated by a schematic two-species ``run-and-tumble'' model., 21 pages, 11 figures

We report a novel state of active matter—a swirlonic state. It is comprised of swirlons, formed by groups of active particles orbiting their common center of mass. These quasi-particles demonstrate a surprising behavior: In response to an external load they move with a constant velocity proportional to the applied force, just as objects in viscous media. The swirlons attract each other and coalesce forming a larger, joint swirlon. The coalescence is extremely slow, decelerating process, resulting in a rarified state of immobile quasi-particles. In addition to the swirlonic state, we observe gaseous, liquid and solid states, depending on the inter-particle and self-driving forces. Interestingly, in contrast to molecular systems, liquid and gaseous states of active matter do not coexist. We explain this unusual phenomenon by the lack of fast particles in active matter. We perform extensive numerical simulations and theoretical analysis. The predictions of the theory agree qualitatively and quantitatively with the simulation results.

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 minimalistic theoretical model, we capture the most relevant features of these active-passive dimers and in particular we predict the crossover between anti-chemotactic and chemotactic behaviour. Moreover we show that merely connecting active particles to chains is sufficient to obtain the crossover from anti-chemotaxis 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 proto-forms of life enabling them to locate a source of nutrients even in the absence of any supporting sensomotoric apparatus.

This note provides a list of errata and their correction for Reference [1].

We report on the emergence of stable self-propelled bound defects in monolayers of active nematics, which form virtual full-integer topological defects in the form of vortices and asters. Through numerical simulations and analytical arguments, we identify the phase space of the bound defect formation in active nematic monolayers. It is shown that an intricate synergy between the nature of active stresses and the flow-aligning behavior of active particles can stabilize the motion of self-propelled positive half-integer defects into specific bound structures. Our findings therefore, uncover conditions for the formation of full integer topological defects in active nematics with the potential for triggering further experiments and theories.

Effective interaction between two asymmetric wedges immersed in a two-dimensional active bath is investigated by computer simulations. The attraction–repulsion transition of effective force between two asymmetric wedges is subjected to the relative position of two wedges, the wedge-to-wedge distance, the active particle density, as well as the apex angle of two wedges. By exchanging the position of the two asymmetric wedges in an active bath, firstly a simple attraction–repulsion transition of effective force occurs, completely different from passive Brownian particles. Secondly the transition of effective force is symmetric for the long-range distance between two asymmetric wedges, while it is asymmetric for the short-range case. Our investigations may provide new possibilities to govern the motion and assembly of microscopic objects by taking advantage of the self-driven behaviour of active particles.

It is now well established that microswimmers can be sorted or segregated fabricating suitable microfluidic devices or using external fields. A natural question is how these techniques can be employed for dividing swimmers of different motility. In this paper, using numerical simulations in the dilute limit, we investigate how motility parameters (time of persistence and velocity) impact the narrow-escape time of active particles from circular domains. We show that the escape time undergoes a crossover between two asymptotic regimes. The control parameters of the crossover is the ratio between the persistence length of the active motion and the typical length scale of the circular domain. We explore the possibility of taking advantage of this finding for sorting active particles by motility parameters.