Self-organization in a bimotility mixture of model microswimmers

The narrow escape problem of a chiral active particle (CAP): an optimal scheme

We report a simulation study on the narrow escape kinetics of a chiral active particle (CAP) confined to a circular domain with a narrow escape opening. The study's main objective is to optimize the CAP's escape chances as a function of the relevant parameters, such as translational and rotational speeds of the CAP, domain size, etc. We identified three regimes in the escape kinetics, namely the noise-dominated regime, the optimal regime, and the chiral activity-dominated regime. In particular, the optimal regime is characterized by an escape scheme that involves a direct passage to the domain boundary at first and then a unidirectional drift along the boundary towards the exit. Furthermore, we propose a non-dimensionalization approach to optimize the escape performance across microorganisms with varying motile characteristics. Additionally, we explore the influence of the translational and rotational noise on the CAP's escape kinetics.

Collective self-optimization of binary mixed heterogeneous populations

To maximize the survival chances of society members, collective self-organization must balance individual interests with promoting the collective welfare. Although situations where group members have equal optimal values are clear, how varying optimal values impacts group dynamics remains unclear. To address this gap, we conducted a self-optimization study of a binary system incorporating communication-enabled active particles with distinct optimal values. We demonstrate that similar particles will spontaneously aggregate and separate from each other to maximize their individual benefits during the process of self-optimization. Our research shows that both types of particles can produce the optimal field values at low density. However, only one type of particle can achieve the optimal field values at medium density. At high densities, neither type of particle is effective in reaching the optimal field values. Interestingly, we observed that during the self-optimization process, the mixture demixed spontaneously under certain circumstances of mixed particles. Particles with higher optimal values developed into larger clusters, while particles with lower optimal values migrated outside of these clusters, resulting in the separation of the mixture. To achieve this separation, suitable noise intensity, particle density, and the significant difference in optimal values were necessary. Our results provide a more profound comprehension of the self-optimization of synthetic or biological agents' communication and provide valuable insight into separating binary species and mixtures.

Large-scale demixing in a binary mixture of cells with rigidity disparity in biological tissues

Physical demixing on large scales of embryonic cell populations is fundamental to metazoan development, but whether a rigidity disparity alone is sufficient to driving large-scale demixing in a binary mixture of cell tissues is still an open question. To answer this question, we study mixing and demixing in a binary mixture of rigidity disparity cell tissues without heterotypic interactions using the Voronoi-based cellular model. Under suitable system parameters, the solid-like cells in the mixture can aggregate into a large cluster and the large-scale demixing occurs, which addresses that a rigidity disparity alone is sufficient to drive large-scale demixing. Remarkably, there exists an optimal temperature or rigidity disparity at which the binary mixture can be separated to the maximum extent. The necessary condition for the separation of mixtures is that the two types of cells are solid-like and liquid-like, respectively. The observation of robust demixing on large scales suggests that the sorting of progenitor cells may occur very early in the development process before robust heterotypic interfacial tensions are established. Our findings are relevant to understanding the mechanisms that drive cell sorting in confluent tissues.

Can playing Spirograph lead to an ordered structure in self-propelled particles?

Understanding the local dynamics of microorganisms infecting a cell could help us develop efficient strategies to counter their aggregation. In the present study we have introduced a simple model of self-propelled particles (SPPs) with constant linear velocity, both in 2 and 3 dimensions, which captures the essential features of a microorganism's aggregation as well the dynamics around an attractive point (AP). The static behavior shows the presence of an icosahedral structure for a finite number of SPPs, and a hexagonal closed packed structure for an infinite number of SPPs, which was confirmed using Steinhardt bond order parameters for a 3-dimensional model. For a single SPP the dynamic behaviour involves the formation of orbits around the AP, which can be categorised into three dynamical regions based on the strength of coupling between the AP and SPP. For weak coupling we observe a rosette-like trajectory reminiscent of the pattern formed by the Spirograph toy. For intermediate coupling, circular trajectories were observed, and for very strong coupling the SPP was static and was always aligned with the AP. The radial distance from the AP to SPP was determined by the angular velocities of the SPP for the rosette-like region whereas for the circular and static regions, it was determined by the coupling constant. Even for a finite number of SPPs we observed the same behavior as long as the SPPs could rotate around the AP without colliding with each other.

Binary mixtures of active and passive particles on a sphere

We study the cooperation and segregation dynamics of binary mixtures of active and passive particles on a sphere. According to the competition between rotational diffusion and polar alignment, we find three distinct phases: a mixed phase and two different demixed phases. When rotational diffusion dominates the dynamics, the demixing is due to the aggregation of passive particles, where active and passive particles respectively occupy two hemispheres. When polar alignment is dominated, the demixing is caused by the aggregation of active particles, where active particles occupy the equator of the sphere and passive particles occupy the two poles of the sphere. In this case, there exist a circulating band cluster and two cambered surface clusters, which is a purely curvature-driven effect with no equivalent in the planar model. When rotational diffusion and polar alignment are comparable, particles are completely mixed. Our findings are relevant to the experimental pursuit of segregation dynamics of binary mixtures on curved surfaces.

Dynamics near planar walls for various model self-phoretic particles

For chemically active particles suspended in a liquid solution and moving by self-phoresis, the dynamics near chemically inert, planar walls is studied theoretically by employing various choices for the activity function, i.e., the spatial distribution of the sites where various chemical reactions take place. We focus on the case of solutions composed of electrically neutral species. This analysis extends previous studies of the case that the chemical activity can be modeled effectively as the release of a "product" molecular species from parts of the surface of the particle by accounting for annihilation of the product molecules by chemical reactions, either on the rest of the surface of the particle or in the volume of the surrounding solution. We show that, for the models considered here, the emergence of "sliding" and "hovering" wall-bound states is a generic, robust feature. However, the details of these states, such as the range of parameters within which they occur, depend on the specific model for the activity function. Additionally, in certain cases there is a reversal of the direction of the motion compared to the one observed if the particle is far away from the wall. We have also studied the changes of the dynamics induced by a direct interaction between the particle and the wall by including a short-ranged repulsive component to the interaction in addition to the steric one (a procedure often employed in numerical simulations of active colloids). Upon increasing the strength of this additional component, while keeping its range fixed, significant qualitative changes occur in the phase portraits of the dynamics near the wall: for sufficiently strong short-ranged repulsion, the sliding steady states of the dynamics are transformed into hovering states. Furthermore, our studies provide evidence for an additional "oscillatory" wall-bound steady state of motion for chemically active particles due to a strong, short-ranged, and direct repulsion. This kind of particle translates along the wall at a distance from it which oscillates between a minimum and a maximum.

Cage dynamics leads to double relaxation of the intermediate scattering function in a binary colloidal system

A system of binary colloids where one fraction of particles is aggregating by forming irreversible bonds and the other fraction of particles only interacts as hard spheres, is simulated using Brownian cluster dynamics. These aggregating species always formed percolating clusters for the case of diffusing hard spheres while for the static case, formation of percolating clusters depended on the fraction of static hard spheres in the system. The dynamics of the hard spheres inside the percolating clusters was studied by restarting the simulation after the kinetics of aggregation was arrested. Two cases were studied, one where the percolated particles moved within the bonds or cage dynamics was allowed and another where the movement within the bonds was not allowed or the cages were static. The hard spheres showed anomalous diffusion in both cases. The mean square displacement showed that for the case of dynamic cages we always had diffusive hard spheres irrespective of the fraction of hard spheres for volume fractions below 0.49. Static cages, depending on the fraction of hard spheres, showed either diffusive or arrested behavior of hard spheres. The intermediate scattering function of only the hard sphere particles showed double relaxation similar to the colloidal glass system for low volume fraction, where the fraction of hard sphere particles was small. For higher fractions we observed only a single stretched exponential. We could differentiate between slow and fast particles for both static and dynamic cages. For the case of static cages the hard spheres were permanently stuck inside the cages while for the case of dynamic cages almost all the hard spheres were moving in and out of the cages.

Simulating squirmers with multiparticle collision dynamics

Multiparticle collision dynamics is a modern coarse-grained simulation technique to treat the hydrodynamics of Newtonian fluids by solving the Navier-Stokes equations. Naturally, it also includes thermal noise. Initially it has been applied extensively to spherical colloids or bead-spring polymers immersed in a fluid. Here, we review and discuss the use of multiparticle collision dynamics for studying the motion of spherical model microswimmers called squirmers moving in viscous fluids.