Self-regulation in self-propelled nematic fluids

Hierarchical defect-induced condensation in active nematics

Topological defects play a central role in the formation and organization of various biological systems. Historically, such nonequilibrium defects have been mainly studied in the context of homogeneous active nematics. Phase-separated systems, in turn, are known to form dense and dynamic nematic bands, but typically lack topological defects. In this paper, we use agent-based simulations of weakly aligning, self-propelled polymers and demonstrate that contrary to the existing paradigm phase-separated active nematics form -1/2 defects. Moreover, these defects, emerging due to interactions among dense nematic bands, constitute a novel second-order collective state. We investigate the morphology of defects in detail and find that their cores correspond to a strong increase in density, associated with a condensation of nematic fluxes. Unlike their analogs in homogeneous systems, such condensed defects form and decay in a different way and do not involve positively charged partners. We additionally observe and characterize lateral arc-like structures that separate from a band's bulk and move in transverse direction. We show that the key control parameters defining the route from stable bands to the coexistence of dynamic lanes and defects are the total density of particles and their path persistence length. We introduce a hydrodynamic theory that qualitatively recapitulates all the main features of the agent-based model, and use it to show that the emergence of both defects and arcs can be attributed to the same anisotropic active fluxes. Finally, we present a way to artificially engineer and position defects, and speculate about experimental verification of the provided model.

Spontaneous flows and dynamics of full-integer topological defects in polar active matter

Polar active matter of self-propelled particles sustain spontaneous flows through the full-integer topological defects. We study theoretically the incompressible flow profiles around ±1 defects induced by polar and dipolar active forces. We show that dipolar forces induce vortical flows around the +1 defect, while the flow around the -1 defect has an 8-fold rotational symmetry. The vortical flow changes its chirality near the +1 defect core in the absence of the friction with a substrate. We show analytically that the flow induced by polar active forces is vortical near the +1 defect and is 4-fold symmetric near the -1 defect, while it becomes uniform in the far-field. For a pair of oppositely charged defects, this polar flow contributes to a mutual interaction force that depends only on the orientation of the defect pair relative to the background polarization, and that enhances defect pair annihilation. This is in contradiction with the effect of dipolar active forces which decay inversely proportional with the defect separation distance. As such, our analyses reveals a long-ranged mechanism for the pairwise interaction between topological defects in polar active matter.

Mesoscopic simulations of active nematics

Coarse-grained, mesoscale simulations are invaluable for studying soft condensed matter because of their ability to model systems in which a background solvent plays a substantial role but is not the primary interest. Such methods generally model passive solvents; however, far-from-equilibrium systems may also be composed of complex solutes suspended in an active fluid. Yet, few coarse-grained simulation methods exist to model an active medium. We introduce an algorithm to simulate active nematics, which builds on multiparticle collision dynamics (MPCD) for passive fluctuating nematohydrodynamics by introducing dipolar activity in the local collision operator. Active nematic MPCD (AN-MPCD) simulations not only exhibit the key characteristics of active nematic turbulence but, as a particle-based algorithm, also reproduce crucial attributes of active particle models. Thus, mesoscopic AN-MPCD is an approach that bridges microscopic and continuum descriptions, allowing simulations of composite active-passive systems.

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.

Rich complex behaviour of self-assembled nanoparticles far from equilibrium

A profoundly fundamental question at the interface between physics and biology remains open: what are the minimum requirements for emergence of complex behaviour from nonliving systems? Here, we address this question and report complex behaviour of tens to thousands of colloidal nanoparticles in a system designed to be as plain as possible: the system is driven far from equilibrium by ultrafast laser pulses that create spatiotemporal temperature gradients, inducing Marangoni flow that drags particles towards aggregation; strong Brownian motion, used as source of fluctuations, opposes aggregation. Nonlinear feedback mechanisms naturally arise between flow, aggregate and Brownian motion, allowing fast external control with minimal intervention. Consequently, complex behaviour, analogous to those seen in living organisms, emerges, whereby aggregates can self-sustain, self-regulate, self-replicate, self-heal and can be transferred from one location to another, all within seconds. Aggregates can comprise only one pattern or bifurcated patterns can coexist, compete, endure or perish.

Instabilities, defects, and defect ordering in an overdamped active nematic

We consider a phenomenological continuum theory for an extensile, overdamped active nematic liquid crystal, applicable in the dense regime. Constructed from general principles, the theory is universal, with parameters independent of any particular microscopic realization. We show that it exhibits two distinct instabilities, one of which arises due to shear forces, and the other due to active torques. Both lead to the proliferation of defects. We focus on the active torque bend instability and find three distinct nonequilibrium steady states including a defect-ordered nematic in which +½ disclinations develop polar ordering. We characterize the phenomenology of these phases and identify the relationship of this theoretical description to experimental realizations and other theoretical models of active nematics.

Role of the active viscosity and self-propelling speed in channel flows of active polar liquid crystals

We study channel flows of active polar liquid crystals (APLCs) focusing on the role played by the active viscosity (β) and the self-propelling speed (ω) on the formation and long time evolution of spontaneous flows using a continuum model. First, we study the onset of spontaneous flows by carrying out a linear stability analysis on two special steady states subject to various physical boundary conditions. We identify a single parameter b1, proportional to a linear combination of the active viscosity and the self-propelling speed, and inversely proportional to a Frank elastic constant, the solvent viscosity, and the liquid crystal relaxation time. We show that the active viscosity and the self-propelling speed influence the onset of spontaneous flows through b1 in that for any fixed value of the bulk activity parameter ζ, large enough |b1| can suppress the spontaneous flow. We then follow spontaneous flows in long time to further investigate the role of β and ω on spatial-temporal structures in the nonlinear regime numerically. The numerical study demonstrates a strong correlation between the most unstable eigenfunction obtained from the linear analysis and the terminal steady state or the persistent, traveling wave structure, revealing the genesis of flow and orientational structures in the active matter system. In the nonlinear regime, a nonzero b1 facilitates the formation of traveling waves in the case of boundary anchoring (the Dirichlet boundary condition) so long as the linear stability analysis predicts an onset of spontaneous flows; in the case of the free boundary condition (the Neumann boundary condition), a stable, spatially homogeneous tilted state always emerges in the presence of two active effects. Finally, we note that various fully out-of-plane spatio-temporal structures can emerge in long time dynamics depending on the boundary condition as well as the initial state of the polarity vector field.

Aspects of the density field in an active nematic

Active nematics are conceptually the simplest orientationally ordered phase of self-driven particles, but have proved to be a perennial source of surprises. We show here through numerical solution of coarse-grained equations for the order parameter and density that the growth of the active nematic phase from the isotropic phase is necessarily accompanied by a clumping of the density. The growth kinetics of the density domains is shown to be faster than the [Formula: see text] law expected for variables governed by a conservation law. Other results presented include the suppression of density fluctuations in the stationary ordered nematic by the imposition of an orienting field. We close by posing some open questions.

Phase separation and emergent structures in an active nematic fluid

We consider a phenomenological continuum theory for an active nematic fluid and show that there exists a universal, model-independent instability which renders the homogeneous nematic state unstable to order fluctuations. Using numerical and analytic tools we show that, in the vicinity of a critical point, this instability leads to a phase-separated state in which the ordered regions form bands in which the direction of nematic order is perpendicular to the direction of the density gradient. We argue that the underlying mechanism that leads to this phase separation is a universal feature of active fluids of different symmetries.

Capillary instability of axisymmetric, active liquid crystal jets

We study linear stability of an infinitely long, axisymmetric, cylindrical active liquid crystal (ALC) jet in a passive isotropic fluid matrix using a polar active liquid crystal (ALC) model. We identify three possible unstable modes (or mechanisms) as the result of the interaction between the flow and the active (or self-propelled) molecular motion. The first unstable mode is related to the polarity vector instability when coupled to the flow field in the presence of the molecular activity. It can be traced back to the inherent polarity vector instability in a bulk active liquid crystal flow. However, it can be grossly amplified in the ALC jet to encompass up to infinitely many unstable growth rates when the long range distortional elastic interaction is weak in certain parameter regimes; it can also be suppressed in other parameter regimes completely. The second unstable mode is related to the classical capillary or Rayleigh instability, which exists in a finite wave interval [0, k(cutoff)]. The new feature for this instability lies in the dependence of the cutoff wave number (k(cutoff)) on the activity of the active matter system. For ALC jets with sufficiently strong contractile activity, the instability can be completely suppressed though. The third unstable mode is due to the active viscous stress. This unstable mode can emerge in the intermediate wave number regime at a sufficiently strong active viscosity and even expand all the way to the zero wave number limit when the Rayleigh unstable mode is absent. It can also be suppressed in the regime of weak active viscous stress. At any given values of the model parameters, the three types of instabilities can show up either individually or in a certain combination, or be completely suppressed altogether. In this paper, we discuss the positive growth rates associated with the instabilities, windows of instability and their dependence on model parameters through extensive numerical computations aided by asymptotic analyses.

Topological structure dynamics revealing collective evolution in active nematics

Topological defects frequently emerge in active matter like bacterial colonies, cytoskeleton extracts on substrates, self-propelled granular or colloidal layers and so on, but their dynamical properties and the relations to large-scale organization and fluctuations in these active systems are seldom touched. Here we reveal, through a simple model for active nematics using self-driven hard elliptic rods, that the excitation, annihilation and transportation of topological defects differ markedly from those in non-active media. These dynamical processes exhibit strong irreversibility in active nematics in the absence of detailed balance. Moreover, topological defects are the key factors in organizing large-scale dynamic structures and collective flows, resulting in multi-spatial temporal effects. These findings allow us to control the self-organization of active matter through topological structures.

Topological defects are observed in a range of active systems, but their dynamical properties are largely unknown. Here, the authors use a simulation of self-propelled hard-rods to generate topological defects in active nematics, finding that their anomalous dynamics may lead to large-scale collective motions.

Casimir effect in swimmer suspensions

We show that the Casimir effect can emerge in microswimmer suspensions. In principle, two effects conspire against the development of Casimir effects in swimmer suspensions. First, at low Reynolds number, the force on any closed volume vanishes, but here the relevant effect is the drag by the flow produced by the swimmers, which can be finite. Second, the fluid velocity and the pressure are linear on the swimmer force dipoles, and averaging over the swimmer orientations would lead to a vanishing effect. However, being that the suspension is a discrete system, the noise terms of the coarse-grained equations depend on the density, which itself fluctuates, resulting in effective nonlinear dynamics. Applying the tools developed for other nonequilibrium systems to general coarse-grained equations for swimmer suspensions, the Casimir drag is computed on immersed objects, and it is found to depend on the correlation function between the rescaled density and dipolar density fields. By introducing a model correlation function with medium-range order, explicit expressions are obtained for the Casimir drag on a body. When the correlation length is much larger than the microscopic cutoff, the average drag is independent of the correlation length, with a range that depends only on the size of the immersed bodies.

Invasion-wave-induced first-order phase transition in systems of active particles

An instability near the transition to collective motion of self-propelled particles is studied numerically by Enskog-like kinetic theory. While hydrodynamics breaks down, the kinetic approach leads to steep solitonlike waves. These supersonic waves show hysteresis and lead to an abrupt jump of the global order parameter if the noise level is changed. Thus they provide a mean-field mechanism to change the second-order character of the phase transition to first order. The shape of the wave is shown to follow a scaling law and to quantitatively agree with agent-based simulations.

Nonlinear field equations for aligning self-propelled rods

We derive a set of minimal and well-behaved nonlinear field equations describing the collective properties of self-propelled rods from a simple microscopic starting point, the Vicsek model with nematic alignment. Analysis of their linear and nonlinear dynamics shows good agreement with the original microscopic model. In particular, we derive an explicit expression for density-segregated, banded solutions, allowing us to develop a more complete analytic picture of the problem at the nonlinear level.