Biphasic, lyotropic, active nematics

Interface morphodynamics in living tissues

Interfaces between distinct tissues or between tissues and environments are common in multicellular organisms. The evolution and stability of these interfaces are essential for tissue development, and their dysfunction can lead to diseases such as cancer. Mounting efforts, either theoretical or experimental, have been devoted to uncovering the morphodynamics of tissue interfaces. Here, we review the recent progress of studies on interface morphodynamics. The regulatory mechanisms governing interface evolution are dissected, with a focus on adhesion, cortical tension, cell activity, extracellular matrix, and microenvironment. We examine the methodologies used to study morphodynamics, emphasizing the characteristics of experimental techniques and theoretical models. Finally, we explore the broader implications of interface morphodynamics in tissue morphogenesis and diseases, offering a comprehensive perspective on this rapidly developing field.

Active phase separation: new phenomenology from non-equilibrium physics

In active systems, whose constituents have non-equilibrium dynamics at local level, fluid-fluid phase separation is widely observed. Examples include the formation of membraneless organelles within cells; the clustering of self-propelled colloidal particles in the absence of attractive forces, and some types of ecological segregation. A schematic understanding of such active phase separation was initially borrowed from what is known for the equilibrium case, in which detailed balance holds at microscopic level. However it has recently become clear that in active systems the absence of detailed balance, although it leave phase separation qualitatively unchanged in some regimes (for example domain growth driven by interfacial tension via Ostwald ripening), can in other regimes radically alter its phenomenology at mechanistic level. For example, microphase separation can be caused by reverse Ostwald ripening, a process that is hard to imagine from an equilibrium perspective. This and other new phenomena arise because, instead of having a single, positive interfacial tension like their equilibrium counterparts, the fluid-fluid interfaces created by active phase separation can have several distinct interfacial tensions governing different properties, some of which can be negative. These phenomena can be broadly understood by studying continuum field theories for a single conserved scalar order parameter (the fluid density), supplemented with a velocity field in cases where momentum conservation is also present. More complex regimes arise in systems described by multiple scalar order parameters (especially with nonreciprocal interactions between these); or when an order parameter undergoes both conserved and non-conserved dynamics (such that the combination breaks detailed balance); or in systems that support orientational long-range order in one or more of the coexisting phases. In this Review, we survey recent progress in understanding the specific role of activity in phase separation, drawing attention to many open questions. We focus primarily on continuum theories, especially those with a single scalar order parameter, reviewing both analytical and numerical work. We compare their predictions with particle-based models, which have mostly been studied numerically although a few have been explicitly coarse-grained to continuum level. We also compare, where possible, with experimental results. In the latter case, qualitative comparisons are broadly encouraging whereas quantitative ones are hindered by the dynamical complexity of most experimental systems relative that of simplified (particle-level or continuum) models of active matter.

Active Fluids Form System-Spanning Filamentary Networks

Recent experimental realizations of liquid-liquid phase separation of active liquid crystals have offered an insight into the interaction between phase separation, ubiquitous in soft matter and biology, and chaotic active flows. In this Letter, we use continuum theory to examine phase separation of an active liquid crystal and a passive fluid and report two new results. First, we provide an analytical derivation of the activity-induced suppression of the phase boundary of the coexistence region-a result first reported in simulations and experiments. We show that the shift in the critical point is a result of the balance between self-stirring active flows and phase-separating diffusive fluxes. Second, we show that this same balance is responsible for dramatically changing the morphology of the phase separated state, resulting in the emergence of a new mixed active phase consisting of a dynamical filamentous active network that invades the entire system area, trapping droplets of passive material. This structure exists even for very low volume fractions of active material. Our work provides an important step towards the goal of understanding how to use activity as a new handle for sculpting interfaces.

Turbulence-to-order transitions in activity-patterned active nematics

We numerically study two-dimensional active nematics with periodic activity patterning. For stripes of activity, we observe a transition from two-dimensional to one-dimensional active turbulence as the maximum active force and distance between activity stripes increases, followed by a transition to stable vortices ordered antiferromagnetically along the stripes and ferromagnetically transverse to the stripes. By comparing to a triangular lattice of activity circles, we find that transitions to two-dimensional active turbulence emerge from interplays between the active length scale and activity density, independent of the patterning geometry. The vortex ordering, however, is highly sensitive to patterning geometry, which we show by comparing the activity stripes to columns of activity circles, where the vortex ordering is lost. Our results provide a mechanism for inducing nonequilibrium phase transitions in active nematics using activity inhomogeneity, which can be further exploited to create activity patterned ordered phases.

Beyond Dipolar Activity: Quadrupolar Stress Drives Collapse of Nematic Order on Frictional Substrates

The field of active nematics has traditionally employed descriptions based on dipolar activity. However, it is theoretically predicted that interactions with a substrate, prevalent in most biological systems, lead to novel forms of activity, such as quadrupolar activity, that are governed by hydrodynamic screening. Here, combining experiments and numerical simulations, we show that upon light-induced solidification of the underlying medium, microtubule-kinesin mixtures undergo a transformation that leads to a biphasic active suspension. Using an active lyotropic model, we prove that the transition is governed by screening effects that alter the dominant form of active stress. Specifically, the combined effect of friction and quadrupolar activity leads to a hierarchical folding that follows the intrinsic bend instability of the active nematic layer. Our results demonstrate the dynamics of the collapse of orientational order in active nematics and present a new route for controlling active matter by modifying the activity through changing the surrounding environment.

Analytical model for the motion and interaction of two-dimensional active nematic defects

We develop an approximate, analytical model for the velocity of defects in active nematics by combining recent results for the velocity of topological defects in nematic liquid crystals with the flow field generated from individual defects in active nematics. Importantly, our model takes into account the long-range interactions between defects that result from the flows they produce as well as the orientational coupling between defects inherent in nematics. Our work complements previous studies of active nematic defect motion by introducing a linear approximation that allows us to treat defect interactions as two-body interactions and incorporates the hydrodynamic screening length as a tuning parameter. We show that the model can analytically predict bound states between two +1/2 winding number defects, effective attraction between two -1/2 defects, and the scaling of a critical unbinding length between ±1/2 defects with activity. The model also gives predictions for the trajectories of defects, such as the scattering of +1/2 defects by -1/2 defects at a critical impact parameter that depends on activity. In the presence of circular confinement, the model predicts a braiding motion for three +1/2 defects that was recently seen in experiments, as well as stable and ergodic trajectories for four or more defects.

Active nematic coherence probed under spatial patterns of distributed activity

A photoresponsive variant of the paradigmatic active nematic fluid made of microtubules and powered by kinesin motors is studied in a conventional two-dimensional interface under blue-light illumination. This advantageously permits the system's performance to be assessed under conditions of spatially distributed activity. Both turbulent and flow aligning conditions are separately analyzed. Under uniform illuminating conditions, active flows get enhanced, in accordance with previous observations. In contrast, patterning the activity appears to disturb the effective activity measured in terms of the vorticity of the elicited flows. We interpret this result as alternative evidence of the important role played by the active length scale in setting not only the textural and flow characteristics of the active nematic but also, most importantly, the range of material integrity. Our research continues to explore perspectives that should pave the way for an effective control of such an admirable material.

Mechanochemical topological defects in an active nematic

We propose a reaction-diffusion system that converts topological information of an active nematic into chemical signals. We show that a curvature-activated reaction dipole is sufficient for creating a system that dynamically senses topology by producing a concentration field possessing local extrema coinciding with ±1/2 defects. The enabling term is analogous to polarization charge density seen in dielectric materials. We demonstrate the ability of this system to identify defects in both passive and active nematics. Our results illustrate that a relatively simple feedback scheme, expressed as a system of partial differential equations, is capable of producing chemical signals in response to inherently nonlocal structures in anisotropic media. We posit that such coarse-grained systems can help generate testable hypotheses for regulated processes in biological systems, such as morphogenesis, and motivate the creation of bio-inspired materials that utilize dynamic coupling between nematic structure and biochemistry.

Traveling waves at the surface of active liquid crystals

Active liquid crystals exert nonequilibrium stresses on their surroundings through constant consumption of energy, giving rise to dynamical steady states not present in equilibrium. The paradigmatic example of an active liquid crystal is a suspension of microtubule bundles powered by kinesin motor proteins, which exhibits self-sustained spatiotemporal chaotic flows. This system has been modelled using continuum theories that couple the microtubule orientation to active flows. Recently the focus has shifted to the interfacial properties of mixtures of active liquid crystals and passive fluids. Active/passive interfaces have been shown to support propagating capillary waves in the absence of inertia and offer a promising route for relating experimental parameters to those of the continuum theory. In this paper we report the derivation of a minimal model that captures the linear dynamics of the interface between an active liquid crystal and a passive fluid. We show that the dynamics of the interface, although powered by active flows throughout the bulk, is qualitatively captured by equations that couple non-reciprocally interface height and nematic director at the interface. This minimal model reproduces the dynamical structure factor evaluated from numerical simulations and the qualitative form of the wave dispersion relation seen in experiments.

Phases and correlations in active nematic granular systems

We investigate the statistical behavior of a system comprising fore-aft symmetric rods confined between two vertically vibrating plates using numerical simulations closely resembling the experimental setup studied by V. Narayan et al., Science, 2007, 317, 105-108. Our focus lies in studying the phase transition of the system from an isotropic phase to a nematic phase as we increase either the rod length or the rod concentration. Our simulations confirm the presence of long-range ordering in the ordered phase. Furthermore, we identify a phase characterized by a periodic ordering profile that disrupts translation symmetry. We also provide a detailed analysis of the translational and rotational diffusive dynamics of the rods. Interestingly, the translational diffusivity of the rods is found to increase with the rod concentration.

Phase ordering in binary mixtures of active nematic fluids

We use a continuum, two-fluid approach to study a mixture of two active nematic fluids. Even in the absence of thermodynamically driven ordering, for mixtures of different activities we observe turbulent microphase separation, where domains form and disintegrate chaotically in an active turbulent background. This is a weak effect if there is no elastic nematic alignment between the two fluid components, but is greatly enhanced in the presence of an elastic alignment or substrate friction. We interpret the results in terms of relative flows between the two species which result from active anchoring at concentration gradients. Our results may have relevance in interpreting epithelial cell sorting and the dynamics of multispecies bacterial colonies.

Simulating dynamics of ellipsoidal particles using lattice Boltzmann method

Anisotropic particles are often encountered in different fields of soft matter and complex fluids. In this work, we present an implementation of the coupled hydrodynamics of solid ellipsoidal particles and the surrounding fluid using the lattice Boltzmann method. A standard link-based mechanism is used to implement the solid-fluid boundary conditions. We develop an implicit method to update the position and orientation of the ellipsoid. This exploits the relations between the quaternion which describes the ellipsoid's orientation and the ellipsoid's angular velocity to obtain a stable and robust dynamic update. The proposed algorithm is validated by looking at four scenarios: (i) the steady translational velocity of a spheroid subject to an external force in different orientations, (ii) the drift of an inclined spheroid subject to an imposed force, (iii) three-dimensional rotational motions in a simple shear flow (Jeffrey's orbits), and (iv) developed fluid flows and self-propulsion exhibited by a spheroidal microswimmer. In all cases the comparison of numerical results shows good agreement with known analytical solutions, irrespective of the choice of the fluid properties, geometrical parameters, and lattice Boltzmann model, thus demonstrating the robustness of the proposed algorithm.

Active nematic ratchet in asymmetric obstacle arrays

We numerically investigate the effect of an asymmetric periodic obstacle array in a two-dimensional active nematic. We find that activity in conjunction with the asymmetry leads to a ratchet effect or unidirectional flow of the fluid along the asymmetry direction. The directional flow is still present even in the active turbulent phase when the gap between obstacles is sufficiently small. We demonstrate that the dynamics of the topological defects transition from flow mirroring to smectic-like as the gap between obstacles is made smaller, and explain this transition in terms of the pinning of negative winding number defects between obstacles. This also leads to a nonmonotonic ratchet effect magnitude as a function of obstacle size, so that there is an optimal obstacle size for ratcheting at fixed activity.

Design rules for controlling active topological defects

Topological defects play a central role in the physics of many materials, including magnets, superconductors, and liquid crystals. In active fluids, defects become autonomous particles that spontaneously propel from internal active stresses and drive chaotic flows stirring the fluid. The intimate connection between defect textures and active flow suggests that properties of active materials can be engineered by controlling defects, but design principles for their spatiotemporal control remain elusive. Here, we propose a symmetry-based additive strategy for using elementary activity patterns, as active topological tweezers, to create, move, and braid such defects. By combining theory and simulations, we demonstrate how, at the collective level, spatial activity gradients act like electric fields which, when strong enough, induce an inverted topological polarization of defects, akin to a negative susceptibility dielectric. We harness this feature in a dynamic setting to collectively pattern and transport interacting active defects. Our work establishes an additive framework to sculpt flows and manipulate active defects in both space and time, paving the way to design programmable active and living materials for transport, memory, and logic.

[The Dynamic Model of the Active-Inactive Cell Interface]

Objective

To explore the morphodynamics of the active-inactive cell monolayer interfaces by using the active liquid crystal model.

Methods

A continuum mechanical model was established based on the active liquid crystal theory and the active-inactive cell monolayer interfaces were established by setting the activity difference of cell monolayers. The theoretical equations were solved numerically by the finite difference and the lattice Boltzmann method.

Results

The active-inactive cell interfaces displayed three typical morphologies, namely, flat interface, wavy interface, and finger-like interface. On the flat interfaces, the cells were oriented perpendicular to the interface, the -1/2 topological defects were clustered in the interfaces, and the interfaces were negatively charged. On the wavy interfaces, cells showed no obvious preference for orientation at the interfaces and the interfaces were neutrally charged. On the finger-like interfaces, cells were tangentially oriented at the interfaces, the +1/2 topological defects were collected at the interfaces, driving the growth of the finger-like structures, and the interfaces were positively charged.

Conclusion

The orientation of the cell alignment at the interface can significantly affect the morphologies of the active-inactive cell monolayer interfaces, which is closely associated with the dynamics of topological defects at the interfaces.

Vortex Lattices in Active Nematics with Periodic Obstacle Arrays

We numerically model a two-dimensional active nematic confined by a periodic array of fixed obstacles. Even in the passive nematic, the appearance of topological defects is unavoidable due to planar anchoring by the obstacle surfaces. We show that a vortex lattice state emerges as activity is increased, and that this lattice may be tuned from "ferromagnetic" to "antiferromagnetic" by varying the gap size between obstacles. We map the rich variety of states exhibited by the system as a function of distance between obstacles and activity, including a pinned defect state, motile defects, the vortex lattice, and active turbulence. We demonstrate that the flows in the active turbulent phase can be tuned by the presence of obstacles, and explore the effects of a frustrated lattice geometry on the vortex lattice phase.

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.

Controlling liquid-liquid phase behaviour with an active fluid

Demixing binary liquids is a ubiquitous transition explained using a well-established thermodynamic formalism that requires the equality of intensive thermodynamics parameters across phase boundaries. Demixing transitions also occur when binary fluid mixtures are driven away from equilibrium, but predicting and designing such out-of-equilibrium transitions remains a challenge. Here we study the liquid-liquid phase separation of attractive DNA nanostars driven away from equilibrium using a microtubule-based active fluid. We find that activity lowers the critical temperature and narrows the range of coexistence concentrations, but only in the presence of mechanical bonds between the liquid droplets and reconfiguring active fluid. Similar behaviours are observed in numerical simulations, suggesting that the activity suppression of the critical point is a generic feature of active liquid-liquid phase separation. Our work describes a versatile platform for building soft active materials with feedback control and providing an insight into self-organization in cell biology.

Active nematics with deformable particles

The hydrodynamic theory of active nematics has been often used to describe the spatio-temporal dynamics of cell flows and motile topological defects within soft confluent tissues. Those theories, however, often rely on the assumption that tissues consist of cells with a fixed, anisotropic shape and do not resolve dynamical cell shape changes due to flow gradients. In this paper we extend the continuum theory of active nematics to include cell shape deformability. We find that circular cells in tissues must generate sufficient active stress to overcome an elastic barrier to deforming their shape in order to drive tissue-scale flows. Above this threshold the systems enter a dynamical steady-state with regions of elongated cells and strong flows coexisting with quiescent regions of isotropic cells.

Geometrical control of interface patterning underlies active matter invasion

Interaction between active materials and the boundaries of geometrical confinement is key to many emergent phenomena in active systems. For living active matter consisting of animal cells or motile bacteria, the confinement boundary is often a deformable interface, and it has been unclear how activity-induced interface dynamics might lead to morphogenesis and pattern formation. Here, we studied the evolution of bacterial active matter confined by a deformable boundary. We found that an ordered morphological pattern emerged at the interface characterized by periodically spaced interfacial protrusions; behind the interfacial protrusions, bacterial swimmers self-organized into multicellular clusters displaying +1/2 nematic defects. Subsequently, a hierarchical sequence of transitions from interfacial protrusions to creeping branches allowed the bacterial active drop to rapidly invade surrounding space with a striking self-similar branch pattern. We found that this interface patterning is geometrically controlled by the local curvature of the interface, a phenomenon we denote as collective curvature sensing. Using a continuum active model, we revealed that the collective curvature sensing arises from enhanced active stresses near high-curvature regions, with the active length scale setting the characteristic distance between the interfacial protrusions. Our findings reveal a protrusion-to-branch transition as a unique mode of active matter invasion and suggest a strategy to engineer pattern formation of active materials.

Phase Separation Driven by Active Flows

We extend the continuum theories of active nematohydrodynamics to model a two-fluid mixture with separate velocity fields for each fluid component, coupled through a viscous drag. The model is used to study an active nematic fluid mixed with an isotropic fluid. We find microphase separation, and argue that this results from an interplay between active anchoring and active flows driven by concentration gradients. The results may be relevant to cell sorting and the formation of lipid rafts in cell membranes.

The crucial role of adhesion in the transmigration of active droplets through interstitial orifices

Active fluid droplets are a class of soft materials exhibiting autonomous motion sustained by an energy supply. Such systems have been shown to capture motility regimes typical of biological cells and are ideal candidates as building-block for the fabrication of soft biomimetic materials of interest in pharmacology, tissue engineering and lab on chip devices. While their behavior is well established in unconstrained environments, much less is known about their dynamics under strong confinement. Here, we numerically study the physics of a droplet of active polar fluid migrating within a microchannel hosting a constriction with adhesive properties, and report evidence of a striking variety of dynamic regimes and morphological features, whose properties crucially depend upon droplet speed and elasticity, degree of confinement within the constriction and adhesiveness to the pore. Our results suggest that non-uniform adhesion forces are instrumental in enabling the crossing through narrow orifices, in contrast to larger gaps where a careful balance between speed and elasticity is sufficient to guarantee the transition. These observations may be useful for improving the design of artificial micro-swimmers, of interest in material science and pharmaceutics, and potentially for cell sorting in microfluidic devices.

Activity-driven tissue alignment in proliferating spheroids

We extend the continuum theory of active nematic fluids to study cell flows and tissue dynamics inside multicellular spheroids, spherical, self-assembled aggregates of cells that are widely used as model systems to study tumour dynamics. Cells near the surface of spheroids have better access to nutrients and therefore proliferate more rapidly than those in the resource-depleted core. Using both analytical arguments and three-dimensional simulations, we find that the proliferation gradients result in flows and in gradients of activity both of which can align the orientation axis of cells inside the aggregates. Depending on environmental conditions and the intrinsic tissue properties, we identify three distinct alignment regimes: spheroids in which all the cells align either radially or tangentially to the surface throughout the aggregate and spheroids with angular cell orientation close to the surface and radial alignment in the core. The continuum description of tissue dynamics inside spheroids not only allows us to infer dynamic cell parameters from experimentally measured cell alignment profiles, but more generally motivates novel mechanisms for controlling the alignment of cells within aggregates which has been shown to influence the mechanical properties and invasive capabilities of tumors.

Collective rotational motion of freely expanding T84 epithelial cell colonies

Coordinated rotational motion is an intriguing, yet still elusive mode of collective cell migration, which is relevant in pathological and morphogenetic processes. Most of the studies on this topic have been carried out on epithelial cells plated on micropatterned substrates, where cell motion is confined in regions of well-defined shapes coated with extracellular matrix adhesive proteins. The driver of collective rotation in such conditions has not been clearly elucidated, although it has been speculated that spatial confinement can play an essential role in triggering cell rotation. Here, we study the growth of epithelial cell colonies freely expanding (i.e. with no physical constraints) on the surface of cell culture plates and focus on collective cell rotation in such conditions, a case which has received scarce attention in the literature. One of the main findings of our work is that coordinated cell rotation spontaneously occurs in cell clusters in the free growth regime, thus implying that cell confinement is not necessary to elicit collective rotation as previously suggested. The extent of collective rotation was size and shape dependent: a highly coordinated disc-like rotation was found in small cell clusters with a round shape, while collective rotation was suppressed in large irregular cell clusters generated by merging of different clusters in the course of their growth. The angular motion was persistent in the same direction, although clockwise and anticlockwise rotations were equally likely to occur among different cell clusters. Radial cell velocity was quite low as compared to the angular velocity, in agreement with the free expansion regime where cluster growth is essentially governed by cell proliferation. A clear difference in morphology was observed between cells at the periphery and the ones in the core of the clusters, the former being more elongated and spread out as compared to the latter. Overall, our results, to our knowledge, provide the first quantitative and systematic evidence that coordinated cell rotation does not require a spatial confinement and occurs spontaneously in freely expanding epithelial cell colonies, possibly as a mechanism for the system.

Active turbulence and spontaneous phase separation in inhomogeneous extensile active gels

We report numerical results for the hydrodynamics of inhomogeneous lyotropic and extensile active nematic gels. By simulating the coupled Cahn-Hilliard, Navier-Stokes, and Beris-Edwards equation for the evolution of the composition, flow and orientational order of an active nematic, we ask whether composition variations are important to determine its emergent physics. As in active gels of uniform composition, we find that increasing either activity or nematic tendency (e.g., overall active matter concentration) triggers a transition between an isotropic passive phase and an active nematic one. We show that composition inhomogeneities are important in the latter phase, where we find three types of possible dynamical regimes. First, we observe regular patterns with defects and vortices: these exist close to the passive-active transition. Second, for larger activity, or deeper in the nematic phase, we find active turbulence, as in active gels of uniform composition, but with exceedingly large composition variation. In the third regime, which is uniquely associated with inhomogeneity and occurs for large nematic tendency and low activity, we observe spontaneous microphase separation into active and passive domains. The microphase separated regime is notable in view of the absence of an explicit demixing term in the underlying free energy which we use, and we provide a theoretical analysis based on the common tangent construction which explains its existence. We hope this regime can be probed experimentally in the future.

Activity-Suppressed Phase Separation

We use a continuum model to examine the effect of activity on a phase-separating mixture of an extensile active nematic and a passive fluid. We highlight the distinct role of (i) previously considered interfacial active stresses and (ii) bulk active stresses that couple to liquid crystalline degrees of freedom. Interfacial active stresses can arrest phase separation, as previously demonstrated. Bulk extensile active stresses can additionally strongly suppress phase separation by sustained self-stirring of the fluid, substantially reducing the size of the coexistence region in the temperature-concentration plane relative to that of the passive system. The phase-separated state is a dynamical emulsion of continuously splitting and merging droplets, as suggested by recent experiments. Using scaling analysis and simulations, we identify various regimes for the dependence of droplet size on activity. These results can provide a criterion for identifying the mechanisms responsible for arresting phase separation in experiments.

Active boundary layers in confined active nematics

The role of boundary layers in conventional liquid crystals is commonly related to the mesogen anchoring on confining walls. In the classical view, anchoring enslaves the orientational field of the passive material under equilibrium conditions. In this work, we show that an active nematic can develop active boundary layers that topologically polarize the confining walls. We find that negatively-charged defects accumulate in the boundary layer, regardless of the wall curvature, and they influence the overall dynamics of the system to the point of fully controlling the behavior of the active nematic in situations of strong confinement. Further, we show that wall defects exhibit behaviors that are essentially different from those of their bulk counterparts, such as high motility or the ability to recombine with another defect of like-sign topological charge. These exotic behaviors result from a change of symmetry induced by the wall in the director field around the defect. Finally, we suggest that the collective dynamics of wall defects might be described in terms of a model equation for one-dimensional spatio-temporal chaos.

Self-mixing in microtubule-kinesin active fluid from nonuniform to uniform distribution of activity

Active fluids have applications in micromixing, but little is known about the mixing kinematics of systems with spatiotemporally-varying activity. To investigate, UV-activated caged ATP is used to activate controlled regions of microtubule-kinesin active fluid and the mixing process is observed with fluorescent tracers and molecular dyes. At low Péclet numbers (diffusive transport), the active-inactive interface progresses toward the inactive area in a diffusion-like manner that is described by a simple model combining diffusion with Michaelis-Menten kinetics. At high Péclet numbers (convective transport), the active-inactive interface progresses in a superdiffusion-like manner that is qualitatively captured by an active-fluid hydrodynamic model coupled to ATP transport. Results show that active fluid mixing involves complex coupling between distribution of active stress and active transport of ATP and reduces mixing time for suspended components with decreased impact of initial component distribution. This work will inform application of active fluids to promote micromixing in microfluidic devices.

Dispersion of activity at an active-passive nematic interface

Efficient nutrient mixing is crucial for the survival of bacterial colonies and other living systems known as active nematics. However, the dynamics of this mixing is non-trivial as there is a coupling between nutrients concentration and velocity field. To address this question, we solve the hydrodynamic equation for active nematics to model the bacterial swarms coupled to an advection-diffusion equation for the activity field, which is proportional to the concentration of nutrients. At the interface between active and passive nematics the activity field is transported by the interfacial flows and in turn it modifies them through the generation of active stresses. We find that the dispersion of this conserved activity field is subdiffusive due to the emergence of a barrier of negative defects at the active-passive interface, which hinders the propagation of the motile positive defects.

Activity gradients in two- and three-dimensional active nematics

We numerically investigate how spatial variations of extensile or contractile active stress affect bulk active nematic systems in two and three dimensions. In the absence of defects, activity gradients drive flows which re-orient the nematic director field and thus act as an effective anchoring force. At high activity, defects are created and the system transitions into active turbulence, a chaotic flow state characterized by strong vorticity. We find that in two-dimensional (2D) systems active torques robustly align +1/2 defects parallel to activity gradients, with defect heads pointing towards contractile regions. In three-dimensional (3D) active nematics disclination lines preferentially lie in the plane perpendicular to activity gradients due to active torques acting on line segments. The average orientation of the defect structures in the plane perpendicular to the line tangent depends on the defect type, where wedge-like +1/2 defects align parallel to activity gradients, while twist defects are aligned anti-parallel. Understanding the response of active nematic fluids to activity gradients is an important step towards applying physical theories to biology, where spatial variations of active stress impact morphogenetic processes in developing embryos and affect flows and deformations in growing cell aggregates, such as tumours.

Bacterial active matter

Bacteria are among the oldest and most abundant species on Earth. Bacteria successfully colonize diverse habitats and play a significant role in the oxygen, carbon, and nitrogen cycles. They also form human and animal microbiota and may become sources of pathogens and a cause of many infectious diseases. Suspensions of motile bacteria constitute one of the most studied examples of active matter: a broad class of non-equilibrium systems converting energy from the environment (e.g., chemical energy of the nutrient) into mechanical motion. Concentrated bacterial suspensions, often termed active fluids, exhibit complex collective behavior, such as large-scale turbulent-like motion (so-called bacterial turbulence) and swarming. The activity of bacteria also affects the effective viscosity and diffusivity of the suspension. This work reports on the progress in bacterial active matter from the physics viewpoint. It covers the key experimental results, provides a critical assessment of major theoretical approaches, and addresses the effects of visco-elasticity, liquid crystallinity, and external confinement on collective behavior in bacterial suspensions.

Quadrupolar active stress induces exotic patterns of defect motion in compressible active nematics

A wide range of living and artificial active matter exists in close contact with substrates and under strong confinement, where in addition to dipolar active stresses, quadrupolar active stresses can become important. Here, we numerically investigate the impact of quadrupolar non-equilibrium stresses on the emergent patterns of self-organisation in non-momentum conserving active nematics. Our results reveal that beyond having stabilising effects, the quadrupolar active forces can induce various modes of topological defect motion in active nematics. In particular, we find the emergence of both polar and nematic ordering of the defects, as well as new patterns of self-organisation that comprise topological defect chains and transient topological defect asters. The results contribute to further understanding of emergent patterns of collective motion and non-equilibrium self-organisation in active matter.

Self-rotation regulates interface evolution in biphasic active matter through taming defect dynamics

Chirality can endow nonequilibrium active matter with unique features and functions. Here, we explore the chiral dynamics in biphasic active nematics composed of self-rotating units that continuously inject energy and angular momentum at the microscale. We show that the self-rotation of units can regularize the boundaries between two phases, rendering sinusoidal-like interfaces, which allow lateral wave propagation and are characterized by chains of ordered antiferromagnetic cross-interface flow vortices. Through the spontaneous coordination of counter-rotating units across the interfaces, topological defects excited by activity are sorted spatiotemporally, where positive defects are locally trapped at the interfaces but, unexpectedly, are transported laterally in a unidirectional rather than wavy mode, whereas inertial negative defects remain spinning in the bulks. Our findings reveal that individual chirality could be harnessed to modulate interfacial morphodynamics in active systems and suggest a potential approach toward controlling topological defects for programmable microfluidics and logic operations.

Phase separation dynamics in deformable droplets

Phase separation can drive spatial organization of multicomponent mixtures. For instance in developing animal embryos, effective phase separation descriptions have been used to account for the spatial organization of different tissue types. Similarly, separation of different tissue types is also observed in stem cell aggregates, where the emergence of a polar organization can mimic early embryonic axis formation. Here, we describe such aggregates as deformable two-phase fluid droplets, which are suspended in a fluid environment (third phase). Using hybrid finite-volume Lattice-Boltzmann simulations, we numerically explore the out-of-equilibrium routes that can lead to the polar equilibrium state of such a droplet. We focus on the interplay between spinodal decomposition and advection with hydrodynamic flows driven by interface tensions, which we characterize by a Peclet number Pe. Consistent with previous work, for large Pe the coarsening process is generally accelerated. However, for intermediate Pe we observe long-lived, strongly elongated droplets, where both phases form an alternating stripe pattern. We show that these "croissant" states are close to mechanical equilibrium and coarsen only slowly through diffusive fluxes in an Ostwald-ripening-like process. Finally, we show that a surface tension asymmetry between both droplet phases leads to transient, rotationally symmetric states whose resolution leads to flows reminiscent of Marangoni flows. Our work highlights the importance of advection for the phase separation process in finite, deformable systems.

Singularity identification for the characterization of topology, geometry, and motion of nematic disclination lines

We introduce a characterization of disclination lines in three dimensional nematic liquid crystals as a tensor quantity related to the so called rotation vector around the line. This quantity is expressed in terms of the nematic tensor order parameter Q, and shown to decompose as a dyad involving the tangent vector to the disclination line and the rotation vector. Further, we derive a kinematic law for the velocity of disclination lines by connecting this tensor to a topological charge density as in the Halperin-Mazenko description of defects in vector models. Using this framework, analytical predictions for the velocity of interacting line disclinations and of self-annihilating disclination loops are given and confirmed through numerical computation.

Coupling Turing stripes to active flows

We numerically solve the active nematohydrodynamic equations of motion, coupled to a Turing reaction-diffusion model, to study the effect of active nematic flow on the stripe patterns resulting from a Turing instability. If the activity is uniform across the system, the Turing patterns dissociate when the flux from active advection balances that from the reaction-diffusion process. If the activity is coupled to the concentration of Turing morphogens, and neighbouring stripes have equal and opposite activity, the system self organises into a pattern of shearing flows, with stripes tending to fracture and slip sideways to join their neighbours. We discuss the role of active instabilities in controlling the crossover between these limits. Our results are of relevance to mechanochemical coupling in biological systems.

Activity-induced instabilities of brain organoids

We present an analytical and numerical investigation of the activity-induced hydrodynamic instabilities in model brain organoids. While several mechanisms have been introduced to explain the experimental observation of surface instabilities in brain organoids, the role of activity has been largely overlooked. Our results show that the active stress generated by the cells can be a, previously overlooked, contributor to the emergence of surface deformations in brain organoids.

Director alignment at the nematic-isotropic interface: elastic anisotropy and active anchoring

Activity in nematics drives interfacial flows that lead to preferential alignment that is tangential or planar for extensile systems (pushers) and perpendicular or homeotropic for contractile ones (pullers). This alignment is known as active anchoring and has been reported for a number of systems and described using active nematic hydrodynamic theories. The latter are based on the one-elastic constant approximation, i.e. they assume elastic isotropy of the underlying passive nematic. Real nematics, however, have different elastic constants, which lead to interfacial anchoring. In this paper, we consider elastic anisotropy in multiphase and multicomponent hydrodynamic models of active nematics and investigate the competition between the interfacial alignment driven by the elastic anisotropy of the passive nematic and the active anchoring. We start by considering systems with translational invariance to analyse the alignment at flat interfaces and, then, consider two-dimensional systems and active nematic droplets. We investigate the competition of the two types of anchoring over a wide range of the other parameters that characterize the system. The results of the simulations reveal that the active anchoring dominates except at very low activities, when the interfaces are static. In addition, we found that the elastic anisotropy does not affect the dynamics but changes the active length that becomes anisotropic. This article is part of the theme issue 'Progress in mesoscale methods for fluid dynamics simulation'.

The many behaviors of deformable active droplets

Active fluids consume fuel at the microscopic scale, converting this energy into forces that can drive macroscopic motions over scales far larger than their microscopic constituents. In some cases, the mechanisms that give rise to this phenomenon have been well characterized, and can explain experimentally observed behaviors in both bulk fluids and those confined in simple stationary geometries. More recently, active fluids have been encapsulated in viscous drops or elastic shells so as to interact with an outer environment or a deformable boundary. Such systems are not as well understood. In this work, we examine the behavior of droplets of an active nematic fluid. We study their linear stability about the isotropic equilibrium over a wide range of parameters, identifying regions in which different modes of instability dominate. Simulations of their full dynamics are used to identify their nonlinear behavior within each region. When a single mode dominates, the droplets behave simply: as rotors, swimmers, or extensors. When parameters are tuned so that multiple modes have nearly the same growth rate, a pantheon of modes appears, including zigzaggers, washing machines, wanderers, and pulsators.

Robustness of topological defects in discrete domains

Topological defects are singular points in vector fields, important in applications ranging from fingerprint detection to liquid crystals to biomedical imaging. In discretized vector fields, topological defects and their topological charge are identified by finite differences or finite-step paths around the tentative defect. As the topological charge is (half) integer, it cannot depend continuously on each input vector in a discrete domain. Instead, it switches discontinuously when vectors change beyond a certain amount, making the analysis of topological defects error prone in noisy data. We improve existing methods for the identification of topological defects by proposing a robustness measure for (i) the location of a defect, (ii) the existence of topological defects and the total topological charge within a given area, (iii) the annihilation of a defect pair, and (iv) the formation of a defect pair. Based on the proposed robustness measure, we show that topological defects in discrete domains can be identified with optimal trade-off between localization precision and robustness. The proposed robustness measure enables uncertainty quantification for topological defects in noisy discretized nematic fields (orientation fields) and polar fields (vector fields).

Optimal Control of Active Nematics

In this work we present the first systematic framework to sculpt active nematic systems, using optimal control theory and a hydrodynamic model of active nematics. We demonstrate the use of two different control fields, (i) applied vorticity and (ii) activity strength, to shape the dynamics of an extensile active nematic that is confined to a disk. In the absence of control inputs, the system exhibits two attractors, clockwise and counterclockwise circulating states characterized by two co-rotating topological +1/2 defects. We specifically seek spatiotemporal inputs that switch the system from one attractor to the other; we also examine phase-shifting perturbations. We identify control inputs by optimizing a penalty functional with three contributions: total control effort, spatial gradients in the control, and deviations from the desired trajectory. This work demonstrates that optimal control theory can be used to calculate nontrivial inputs capable of restructuring active nematics in a manner that is economical, smooth, and rapid, and therefore will serve as a guide to experimental efforts to control active matter.

Pattern formation in active model C with anchoring: bands, aster networks, and foams

We study the dynamics of pattern formation in a minimal model for active mixtures made of microtubules and molecular motors. We monitor the evolution of the (conserved) microtubule density and of the (non-conserved) nematic order parameter, focusing on the effects of an "anchoring" term that provides a direct coupling between the preferred microtubule direction and their density gradient. The key control parameter is the ratio between activity and elasticity. When elasticity dominates, the interplay between activity and anchoring leads to formation of banded structures that can undergo additional bending, rotational or splaying instabilities. When activity dominates, the nature of anchoring instead gives rise to a range of active cellular solids, including aster-like networks, disordered foams and spindle-like patterns. We speculate that the introduced "active model C" with anchoring is a minimal model to describe pattern formation in a biomimetic analogue of the microtubule cytoskeleton.

Collective chemotaxis of active nematic droplets

Collective chemotaxis plays a key role in the navigation of cell clusters in, e.g., embryogenesis and cancer metastasis. Using the active nematic continuum equations, coupled to a chemical field that regulates activity, we demonstrate and explain a physical mechanism that results in collective chemotaxis. The activity naturally leads to cell polarization at the cluster interface which induces outward flows. The chemical gradient then breaks the symmetry of the flow field, leading to a net motion. The velocity is independent of the cluster size, in agreement with experiment.

Pattern Formation and Defect Ordering in Active Chiral Nematics

Many biological systems display intriguing chiral patterns and dynamics. Here, we present an active nematic theory accounting for individual spin to explore the collective handedness in chiral rod-shaped aggregations. We show that coordinated individual spin and motility can engender a vortex-array pattern with chirality and drive ordering of topological defects. During this chiral process, the stationary trefoil-like defects self-organize into a periodic, hexagon-dominated polygonal network, which segregates persistently rotating cometlike defects in pairs within each polygon, leading to a translation symmetry at the global scale while a broken reflection symmetry at the local scale. Such defect ordering agrees exactly with the Voronoi tiling of two-dimensional space and the emergence of the hexagonal symmetry is deciphered in analogy with topological charge neutralization. We calculate energy barriers to the topological transition of the defect ordering and explain the existing metastable states with nonhexagonal polygons. Our findings shed light on the chiral morphodynamics in life processes and also suggest a potential route towards tuning self-organization in active materials.

Propagation of active nematic-isotropic interfaces on substrates

Motivated by results for the propagation of active-passive interfaces of bacterial Serratia marcescens swarms [Nat. Commun., 2018, 9, 5373], we used a hydrodynamic multiphase model to investigate the propagation of interfaces of active nematics on substrates. We characterized the active nematic phase of the model through the calculation of the spatial and temporal auto correlation functions and the energy spectrum and discussed its description of the statistical dynamics of the swarms reported in the experiment. We then studied the propagation of circular and flat active-passive interfaces. We found that the closing time of the circular passive domain decays quadratically with the activity and that the structure factor of the flat interface is similar to that reported for the swarms, with an activity dependent exponent. Finally, the effect of the substrate friction was investigated. We found an activity dependent threshold, above which the turbulent active nematic forms isolated islands that shrink until the system becomes isotropic and below which the active nematic expands, with a well defined propagating interface. We also found that the interface becomes static in the presence of a friction gradient.

Active nematics with anisotropic friction: the decisive role of the flow aligning parameter

We use continuum simulations to study the impact of anisotropic hydrodynamic friction on the emergent flows of active nematics. We show that, depending on whether the active particles align with or tumble in their collectively self-induced flows, anisotropic friction can result in markedly different patterns of motion. In a flow-aligning regime and at high anisotropic friction, the otherwise chaotic flows are streamlined into flow lanes with alternating directions, reproducing the experimental laning state that has been obtained by interfacing microtubule-motor protein mixtures with smectic liquid crystals. Within a flow-tumbling regime, however, we find that no such laning state is possible. Instead, the synergistic effects of friction anisotropy and flow tumbling can lead to the emergence of bound pairs of topological defects that align at an angle to the easy flow direction and navigate together throughout the domain. In addition to confirming the mechanism behind the laning states observed in experiments, our findings emphasise the role of the flow aligning parameter in the dynamics of active nematics.

Chiral stresses in nematic cell monolayers

Recent experiments on monolayers of spindle-like cells plated on adhesive stripe-shaped domains have provided a convincing demonstration that certain types of collective phenomena in epithelia are well described by active nematic hydrodynamics. While recovering some of the hallmark predictions of this framework, however, these experiments have also revealed a number of unexpected features that could be ascribed to the existence of chirality over length scales larger than the typical size of a cell. In this article we elaborate on the microscopic origin of chiral stresses in nematic cell monolayers and investigate how chirality affects the motion of topological defects, as well as the collective motion in stripe-shaped domains. We find that chirality introduces a characteristic asymmetry in the collective cellular flow, from which the ratio between chiral and non-chiral active stresses can be inferred by particle-image-velocimetry measurements. Furthermore, we find that chirality changes the nature of the spontaneous flow transition under confinement and that, for specific anchoring conditions, the latter has the structure of an imperfect pitchfork bifurcation.

Topology and Morphology of Self-Deforming Active Shells

We present a generic framework for modeling three-dimensional deformable shells of active matter that captures the orientational dynamics of the active particles and hydrodynamic interactions on the shell and with the surrounding environment. We find that the cross talk between the self-induced flows of active particles and dynamic reshaping of the shell can result in conformations that are tunable by varying the form and magnitude of active stresses. We further demonstrate and explain how self-induced topological defects in the active layer can direct the morphodynamics of the shell. These findings are relevant to understanding morphological changes during organ development and the design of bioinspired materials that are capable of self-organization.

Active matter invasion

Biologically active materials such as bacterial biofilms and eukaryotic cells thrive in confined micro-spaces. Here, we show through numerical simulations that confinement can serve as a mechanical guidance to achieve distinct modes of collective invasion when combined with growth dynamics and the intrinsic activity of biological materials. We assess the dynamics of the growing interface and classify these collective modes of invasion based on the activity of the constituent particles of the growing matter. While at small and moderate activities the active material grows as a coherent unit, we find that blobs of active material collectively detach from the cohort above a well-defined activity threshold. We further characterise the mechanical mechanisms underlying the crossovers between different modes of invasion and quantify their impact on the overall invasion speed.