Chemical Physics of Active Matter

Interplay of chemotactic force, Péclet number, and dimensionality dictates the dynamics of auto-chemotactic chiral active droplets

In living and synthetic active matter systems, the constituents can self-propel and interact with each other and with the environment through various physicochemical mechanisms. Among these mechanisms, chemotactic and auto-chemotactic effects are widely observed. The impact of (auto-)chemotactic effects on achiral active matter has been a recent research focus. However, the influence of these effects on chiral active matter remains elusive. Here, we develop a Brownian dynamics model coupled with a diffusion equation to examine the dynamics of auto-chemotactic chiral active droplets in both quasi-two-dimensional (2D) and three-dimensional (3D) systems. By quantifying the droplet trajectory as a function of the dimensionless Péclet number and chemotactic strength, our simulations well reproduce the curling and helical trajectories of nematic droplets in a surfactant-rich solution reported by Krüger et al. [Phys. Rev. Lett. 117, 048003 (2016)]. The modeled curling trajectory in 2D exhibits an emergent chirality, also consistent with the experiment. We further show that the geometry of the chiral droplet trajectories, characterized by the pitch and diameter, can be used to infer the velocities of the droplet. Interestingly, we find that, unlike the achiral case, the velocities of chiral active droplets show dimensionality dependence: its mean instantaneous velocity is higher in 3D than in 2D, whereas its mean migration velocity is lower in 3D than in 2D. Taken together, our particle-based simulations provide new insights into the dynamics of auto-chemotactic chiral active droplets, reveal the effects of dimensionality, and pave the way toward their applications, such as drug delivery, sensors, and micro-reactors.

Synchronized rotations of active particles on chemical substrates

Many microorganisms use chemical 'signaling' - a quintessential self-organizing strategy in non-equilibrium - that can induce spontaneous aggregation and coordinated motion. Using synthetic signaling as a design principle, we construct a minimal model of active Brownian particles (ABPs) having soft repulsive interactions on a chemically quenched patterned substrate. The interplay between chemo-phoretic interactions and activity is numerically investigated for a proposed variant of the Keller-Segel model for chemotaxis. Such competition not only results in a chemo-motility-induced phase-separated state, but also results in a new cohesive clustering phase with synchronized rotations. Our results suggest that rotational order can emerge in systems by virtue of activity and repulsive interactions alone without an explicit alignment interaction. These rotations can also be exploited by designing mechanical devices that can generate reorienting torques using active particles.

Dynamic phase transition induced by active molecules in a supercooled liquid

The purpose of this work is to use active particles to investigate the effect of facilitation on supercooled liquids. To this end we examine the behavior of a model supercooled liquid that is doped with a mixture of active particles and slowed particles. To simulate the facilitation mechanism, the activated particles are subjected to a force that follows the mobility of their most mobile neighboring molecule, while the slowed particles experience a friction force. Upon activation, we observe a fluidization of the entire medium along with a significant increase in dynamic heterogeneity. This effect is reminiscent of the fluidization observed experimentally when introducing molecular motors into soft materials. Interestingly, when the characteristic time τ_{μ}, used to define the mobility in the facilitation mechanism, matches the physical time t^{*} that characterizes the spontaneous cooperativity of the material, we observe a phase transition accompanied by structural aggregation of the active molecules. This transition is characterized by a sharp increase in fluidization and dynamic heterogeneity.

Bridging scales in a multiscale pattern-forming system

Self-organized pattern formation is vital for many biological processes. Reaction-diffusion models have advanced our understanding of how biological systems develop spatial structures, starting from homogeneity. However, biological processes inherently involve multiple spatial and temporal scales and transition from one pattern to another over time, rather than progressing from homogeneity to a pattern. To deal with such multiscale systems, coarse-graining methods are needed that allow the dynamics to be reduced to the relevant degrees of freedom at large scales, but without losing information about the patterns at small scales. Here, we present a semiphenomenological approach which exploits mass conservation in pattern formation, and enables reconstruction of information about patterns from the large-scale dynamics. The basic idea is to partition the domain into distinct regions (coarse grain) and determine instantaneous dispersion relations in each region, which ultimately inform about local pattern-forming instabilities. We illustrate our approach by studying the Min system, a paradigmatic model for protein pattern formation. By performing simulations, we first show that the Min system produces multiscale patterns in a spatially heterogeneous geometry. This prediction is confirmed experimentally by in vitro reconstitution of the Min system. Using a recently developed theoretical framework for mass-conserving reaction-diffusion systems, we show that the spatiotemporal evolution of the total protein densities on large scales reliably predicts the pattern-forming dynamics. Our approach provides an alternative and versatile theoretical framework for complex systems where analytical coarse-graining methods are not applicable, and can, in principle, be applied to a wide range of systems with an underlying conservation law.

Quantitative kinetic theory of flocking with three-particle closure

We consider aligning self-propelled particles in two dimensions. Their motion is given by Langevin equations and includes nonadditive N-particle interactions. The qualitative behavior is as for the famous Vicsek model. We develop a kinetic theory of flocking beyond mean field. In particular, we self-consistently take into account the full pair correlation function. We find excellent quantitative agreement of the pair correlations with direct agent-based simulations within the disordered regime. Furthermore we use a closure relation to incorporate spatial correlations of three particles. In that way we achieve good quantitative agreement of the onset of flocking with direct simulations. Compared to mean-field theory, the flocking transition is shifted significantly toward lower noise because directional correlations favor disorder. We compare our theory with a recently developed Landau-kinetic theory.

Characterization of MIPS in a suspension of repulsive active Brownian particles through dynamical features

We study a two-dimensional system composed by Active Brownian Particles (ABPs), focusing on the onset of Motility Induced Phase Separation (MIPS), by means of molecular dynamics simulations. For a pure hard-disk system with no translational diffusion, the phase diagram would be completely determined by their density and Péclet number. In our model, two additional effects are present: translational noise and the overlap of particles; we study the effects of both in the phase space. As we show, the second effect can be mitigated if we use, instead of the standard Weeks-Chandler-Andersen potential, a stiffer potential: the pseudo-hard sphere potential. Moreover, in determining the boundary of our phase space, we explore different approaches to detect MIPS and conclude that observing dynamical features, via the non-Gaussian parameter, is more efficient than observing structural ones, such as through the local density distribution function. We also demonstrate that the Vogel-Fulcher equation successfully reproduces the decay of the diffusion as a function of density, with the exception of very high densities. Thus, in this regard, the ABP system behaves similar to a fragile glass.

Orientation of motion of a flat folding nano-swimmer in soft matter

It is well established that anisotropic molecules do have a preferential direction of motion at short time scales that is washed out at larger times by Brownian noise. Anisotropic molecular motors are able to move at lower temperatures when Brownian noise is smaller suggesting the possibility of oriented motion for larger time scales. We use molecular dynamics simulations to investigate that possibility, calculating the displacements of a simple flat folding molecular nano-swimmer embedded in soft matter. We find actually that the motor displacement is oriented in the direction of its length. We note that the observed orientation of the displacement explains the experimental polarization effect in surface relief gratings formation in agreement with the caterpillar model for azobenzene SRG formation mechanism. That result also suggests a simple route for the creation of molecular motors with oriented displacements.

Pattern-induced local symmetry breaking in active-matter systems

The emergence of macroscopic order and patterns is a central paradigm in systems of (self-)propelled agents and a key component in the structuring of many biological systems. The relationships between the ordering process and the underlying microscopic interactions have been extensively explored both experimentally and theoretically. While emerging patterns often show one specific symmetry (e.g., nematic lane patterns or polarized traveling flocks), depending on the symmetry of the alignment interactions patterns with different symmetries can apparently coexist. Indeed, recent experiments with an actomysin motility assay suggest that polar and nematic patterns of actin filaments can interact and dynamically transform into each other. However, theoretical understanding of the mechanism responsible remains elusive. Here, we present a kinetic approach complemented by a hydrodynamic theory for agents with mixed alignment symmetries, which captures the experimentally observed phenomenology and provides a theoretical explanation for the coexistence and interaction of patterns with different symmetries. We show that local, pattern-induced symmetry breaking can account for dynamically coexisting patterns with different symmetries. Specifically, in a regime with moderate densities and a weak polar bias in the alignment interaction, nematic bands show a local symmetry-breaking instability within their high-density core region, which induces the formation of polar waves along the bands. These instabilities eventually result in a self-organized system of nematic bands and polar waves that dynamically transform into each other. Our study reveals a mutual feedback mechanism between pattern formation and local symmetry breaking in active matter that has interesting consequences for structure formation in biological systems.

Universal reshaping of arrested colloidal gels via active doping

Colloids that interact via a short-range attraction serve as the primary building blocks for a broad range of self-assembled materials. However, one of the well-known drawbacks to this strategy is that these building blocks rapidly and readily condense into a metastable colloidal gel. Using computer simulations, we illustrate how the addition of a small fraction of purely repulsive self-propelled colloids, a technique referred to as active doping, can prevent the formation of this metastable gel state and drive the system toward its thermodynamically favored crystalline target structure. The simplicity and robust nature of this strategy offers a systematic and generic pathway to improving the self-assembly of a large number of complex colloidal structures. We discuss in detail the process by which this feat is accomplished and provide quantitative metrics for exploiting it to modulate the self-assembly. We provide evidence for the generic nature of this approach by demonstrating that it remains robust under a number of different anisotropic short-ranged pair interactions in both two and three dimensions. In addition, we report on a novel microphase in mixtures of passive and active colloids. For a broad range of self-propelling velocities, it is possible to stabilize a suspension of fairly monodisperse finite-size crystallites. Surprisingly, this microphase is also insensitive to the underlying pair interaction between building blocks. The active stabilization of these moderately sized monodisperse clusters is quite remarkable and should be of great utility in the design of hierarchical self-assembly strategies. This work further bolsters the notion that active forces can play a pivotal role in directing colloidal self-assembly.