Defect annihilation and proliferation in active nematics

Control of Modular Tissue Flows Shaping the Embryo in Avian Gastrulation

Avian gastrulation requires coordinated flows of thousands of cells to form the body plan. We quantified these flows using their fundamental kinematic units: one attractor and two repellers constituting its Dynamic Morphoskeleton (DM). We have also elucidated the mechanistic origin of the attractor, marking the primitive streak (PS), and controlled its shape, inducing gastrulation flows in the chick embryo that are typical of other vertebrates. However, the origins of repellers and dynamic embryo shape remain unclear. Here, we address these questions using active matter physics and experiments. Repeller 1, separating the embryo proper (EP) from extraembryonic (EE) tissues, arises from the tug-of-war between EE epiboly and EP isotropic myosin-induced active stress. Repeller 2, bisecting the anterior and posterior PS and associated with embryo shape change, arises from anisotropic myosin-induced active intercalation in the mesendoderm. Combining mechanical confinement with inhibition of mesendoderm induction, we eliminated either one or both repellers, as predicted by our model. Our results reveal a remarkable modularity of avian gastrulation flows delineated by the DM, uncovering the mechanistic roles of EE epiboly, EP active constriction, mesendoderm intercalation and ingression. These findings offer a new perspective for deconstructing morphogenetic flows, uncovering their modular origin, and aiding synthetic morphogenesis.

Coexistence of Defect Morphologies in Three-Dimensional Active Nematics

We establish how active stress globally affects the morphology of disclination lines of a three-dimensional active nematic liquid crystal under chaotic flow. Thanks to a defect detection algorithm based on the local nematic orientation, we show that activity selects a crossover length scale in between the size of small defect loops and that of long and tangled defect lines of fractal dimension 2. This length scale crossover is consistent with the scaling of the average separation between defects as a function of activity. Moreover, on the basis of numerical simulation in a 3D periodic geometry, we show the presence of a network of regular defect loops, contractible onto the 3-torus, always coexisting with wrapping defect lines. While the length of regular defects scales linearly with the emerging active length scale, it verifies an inverse quadratic dependence for wrapping defects. The shorter the active length scale, the more the defect lines wrap around the periodic boundaries, resulting in extremely long and buckled structures.

Analysis of the number of topological defects in active nematic fluids under applied shear flow

The number of topological defects in the shear flow of active nematic fluids is numerically investigated in this study. The evolution of the flow state of extensile active nematic fluids is explored by increasing the activity of active nematic fluids. Evidently, medium-activity active nematic fluids exhibit a highly ordered vortex lattice fluid state. However, high-activity active nematic fluids exhibit a meso-scale turbulent flow accompanied by topological defects. The number of topological defects (Ndef) increases with increasing shear Reynolds number (Res). Fluid viscosity strongly influences Ndef, while the influence of fluid density is relatively weak. Ndef decreases with increasing activity length scale (lζ) value. A small Res value strongly influences Ndef, whereas a large lζ value only weakly influences Ndef. As the activity increases, Ndef in contractile active nematic fluids becomes larger than that of extensile active nematic fluids.

Ordering kinetics in the active Ising model

We undertake a numerical study of the ordering kinetics in the two-dimensional (2D) active Ising model (AIM), a discrete flocking model with a conserved density field coupled to a nonconserved magnetization field. We find that for a quench into the liquid-gas coexistence region and in the ordered liquid region, the characteristic length scale of both the density and magnetization domains follows the Lifshitz-Cahn-Allen growth law, R(t)∼t^{1/2}, consistent with the growth law of passive systems with scalar order parameter and nonconserved dynamics. The system morphology is analyzed with the two-point correlation function and its Fourier transform, the structure factor, which conforms to the well-known Porod's law, a manifestation of the coarsening of compact domains with smooth boundaries. We also find the domain growth exponent unaffected by different noise strengths and self-propulsion velocities of the active particles. However, transverse diffusion is found to play the most significant role in the growth kinetics of the AIM. We extract the same growth exponent by solving the hydrodynamic equations of the AIM.

Topological defects in multi-layered swarming bacteria

Topological defects, which are singular points in a director field, play a major role in shaping active systems. Here, we experimentally study topological defects and the flow patterns around them, that are formed during the highly rapid dynamics of swarming bacteria. The results are compared to the predictions of two-dimensional active nematics. We show that, even though some of the assumptions underlying the theory do not hold, the swarm dynamics is in agreement with two-dimensional nematic theory. In particular, we look into the multi-layered structure of the swarm, which is an important feature of real, natural colonies, and find a strong coupling between layers. Our results suggest that the defect-charge density is hyperuniform, i.e., that long range density-fluctuations are suppressed.

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.

Machine learning topological defects in confluent tissues

Active nematics is an emerging paradigm for characterizing biological systems. One aspect of particularly intense focus is the role active nematic defects play in these systems, as they have been found to mediate a growing number of biological processes. Accurately detecting and classifying these defects in biological systems is, therefore, of vital importance to improving our understanding of such processes. While robust methods for defect detection exist for systems of elongated constituents, other systems, such as epithelial layers, are not well suited to such methods. Here, we address this problem by developing a convolutional neural network to detect and classify nematic defects in confluent cell layers. Crucially, our method is readily implementable on experimental images of cell layers and is specifically designed to be suitable for cells that are not rod shaped, which we demonstrate by detecting defects on experimental data using the trained model. We show that our machine learning model outperforms current defect detection techniques and that this manifests itself in our method as requiring less data to accurately capture defect properties. This could drastically improve the accuracy of experimental data interpretation while also reducing costs, advancing the study of nematic defects in biological systems.

Probing active nematics with in situ microfabricated elastic inclusions

In this work, we report a direct measurement of the forces exerted by a tubulin/kinesin active nematic gel as well as its complete rheological characterization, including the quantification of its shear viscosity, η, and its activity parameter, α. For this, we develop a method that allows us to rapidly photo-polymerize compliant elastic inclusions in the continuously remodeling active system. Moreover, we quantitatively settle long-standing theoretical predictions, such as a postulated relationship encoding the intrinsic time scale of the active nematic in terms of η and α. In parallel, we infer a value for the nematic elasticity constant, K, by combining our measurements with the theorized scaling of the active length scale. On top of the microrheology capabilities, we demonstrate strategies for defect encapsulation, quantification of defect mechanics, and defect interactions, enabled by the versatility of the microfabrication strategy that allows to combine elastic motifs of different shapes and stiffnesses that are fabricated in situ.

Maximally mixing active nematics

Active nematics are an important new paradigm in soft condensed matter systems. They consist of rodlike components with an internal driving force pushing them out of equilibrium. The resulting fluid motion exhibits chaotic advection, in which a small patch of fluid is stretched exponentially in length. Using simulation, this paper shows that this system can exhibit stable periodic motion when confined to a sufficiently small square with periodic boundary conditions. Moreover, employing tools from braid theory, we show that this motion is maximally mixing, in that it optimizes the (dimensionless) "topological entropy"-the exponential stretching rate of a material line advected by the fluid. That is, this periodic motion of the defects, counterintuitively, produces more chaotic mixing than chaotic motion of the defects. We also explore the stability of the periodic state. Importantly, we show how to stabilize this orbit into a larger periodic tiling, a critical necessity for it to be seen in future experiments.

Mechanochemical Active Feedback Generates Convergence Extension in Epithelial Tissue

Convergence extension, the simultaneous elongation of tissue along one axis while narrowing along a perpendicular axis, occurs during embryonic development. A fundamental process that contributes to shaping the organism, it happens in many different species and tissue types. Here, we present a minimal continuum model, that can be directly linked to the controlling microscopic biochemistry, which shows spontaneous convergence extension. It is comprised of a 2D viscoelastic active material with a mechanochemical active feedback mechanism coupled to a substrate via friction. Robust convergent extension behavior emerges beyond a critical value of the activity parameter and is controlled by the boundary conditions and the coupling to the substrate. Oscillations and spatial patterns emerge in this model when internal dissipation dominates over friction, as well as in the active elastic limit.

Non-spherical assemblies of chitin nanocrystals by drop impact assembly

Preparing complex non-spherical assemblies of elongated nanoparticles and exploring their topological conformations is a challenge due to liquid crystals' mobility and elastic distortion. Here, we fabricated a variety of non-spherical liquid crystal assemblies of chitin nanocrystals (ChNCs) in a coagulation bath containing sodium triphosphate (STP) by drop impact assembly method, and the forming mechanism and internal topology were systematically investigated. The collection height, ChNCs concentration, and STP concentration have significant influence on the shape and size of the assembled structures. Long-range ordered structures and long-lived topological textures of the ChNCs liquid crystal can be obtained since a molecular interaction of hydrogen bonding and electrostatic attractions between ChNCs and STP occur during the impact assembly. Rheological and kinetic analysis suggested the shear thinning behavior of the ChNCs liquid crystals and the rapid gelation phenomenon of ChNCs induced by STP. Morphology results showed that the rod-like ChNCs in the non-spherical assemblies were orderly and closely arranged with periodic repetition and layered structure. The non-spherical assemblies of ChNCs liquid crystals can be used as carriers of carbon nanotubes, magnetic Fe3O4 nanoparticles, synthesized polymers, and anticancer drugs for functional composite applications. The drop impact assembly method of ChNCs liquid crystal structure is highly controllable on the composition, morphology, and function, which shows promising applications in energy, environmental-friendly, and bioactive materials.

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.

Molecular-scale substrate anisotropy, crowding and division drive collective behaviours in cell monolayers

The ability of cells to reorganize in response to external stimuli is important in areas ranging from morphogenesis to tissue engineering. While nematic order is common in biological tissues, it typically only extends to small regions of cells interacting via steric repulsion. On isotropic substrates, elongated cells can co-align due to steric effects, forming ordered but randomly oriented finite-size domains. However, we have discovered that flat substrates with nematic order can induce global nematic alignment of dense, spindle-like cells, thereby influencing cell organization and collective motion and driving alignment on the scale of the entire tissue. Remarkably, single cells are not sensitive to the substrate's anisotropy. Rather, the emergence of global nematic order is a collective phenomenon that requires both steric effects and molecular-scale anisotropy of the substrate. To quantify the rich set of behaviours afforded by this system, we analyse velocity, positional and orientational correlations for several thousand cells over days. The establishment of global order is facilitated by enhanced cell division along the substrate's nematic axis, and associated extensile stresses that restructure the cells' actomyosin networks. Our work provides a new understanding of the dynamics of cellular remodelling and organization among weakly interacting cells.

Rectified Rotational Dynamics of Mobile Inclusions in Two-Dimensional Active Nematics

We investigate the dynamics of mobile inclusions embedded in 2D active nematics. The interplay between the inclusion shape, boundary-induced nematic order, and autonomous flows powers the inclusion motion. Disks and achiral gears exhibit unbiased rotational motion, but with distinct dynamics. In comparison, chiral gear-shaped inclusions exhibit long-term rectified rotation, which is correlated with dynamics and polarization of nearby +1/2 topological defects. The chirality of defect polarities and the active nematic texture around the inclusion correlate with the inclusion's instantaneous rotation rate. Inclusions provide a promising tool for probing the rheological properties of active nematics and extracting ordered motion from their inherently chaotic motion.

Microdomains and stress distributions in bacterial monolayers on curved interfaces

Monolayers of growing non-motile rod-shaped bacteria act as active nematic materials composed of hard particles rather than the flexible components of other commonly studied active nematics. The organization of these granular monolayers has been studied on flat surfaces but not on curved surfaces, which are known to change the behavior of other active nematics. We use molecular dynamics simulations to track alignment and stress in growing monolayers fixed to curved surfaces, and investigate how these vary with changing surface curvature and cell aspect ratio. We find that the length scale of alignment (measured by average microdomain size) increases with cell aspect ratio and decreases with curvature. Additionally, we find that alignment controls the distribution of extensile stresses in the monolayer by concentrating stress in negative-order regions. These results connect active nematic physics to bacterial monolayers and can be applied to model bacteria growing on droplets, such as oil-degrading marine bacteria.

Hydrodynamic Enhancement of p-atic Defect Dynamics

We investigate numerically and analytically the effects of hydrodynamics on the dynamics of topological defects in p-atic liquid crystals, i.e., two-dimensional liquid crystals with p-fold rotational symmetry. Importantly, we find that hydrodynamics fuels a generic passive self-propulsion mechanism for defects of winding number s=(p-1)/p and arbitrary p. Strikingly, we discover that hydrodynamics always accelerates the annihilation dynamics of pairs of ±1/p defects and that, contrary to expectations, this effect increases with p. Our Letter paves the way toward understanding cell intercalation and other remodeling events in epithelial layers.

Defect self-propulsion in active nematic films with spatially varying activity

We study the dynamics of topological defects in active nematic films with spatially varying activity and consider two set-ups: (i) a constant activity gradient and (ii) a sharp jump in activity. A constant gradient of extensile (contractile) activity endows the comet-like +1/2 defect with a finite vorticity that drives the defect to align its nose in the direction of decreasing (increasing) gradient. A constant gradient does not, however, affect the known self-propulsion of the +1/2 defect and has no effect on the -1/2 that remains a non-motile particle. A sharp jump in activity acts like a wall that traps the defects, affecting the translational and rotational motion of both charges. The +1/2 defect slows down as it approaches the interface and the net vorticity tends to reorient the defect polarization so that it becomes perpendicular to the interface. The -1/2 defect acquires a self-propulsion towards the activity interface, while the vorticity-induced active torque tends to align the defect to a preferred orientation. This effective attraction of the negative defects to the wall is consistent with the observation of an accumulation of negative topological charge at both active/passive interfaces and physical boundaries.

Emergence of macroscopic directional motion of deformable active cells in confined structures

There is now growing evidence of collective turbulentlike motion of cells in dense tissues. However, how to control and harness this collective motion is an open question. We investigate the transport of deformable active cells in a periodically asymmetric channel by using a phase-field model. We demonstrate that collective turbulent-like motion of cells can power and steer the macroscopic directional motion through the ratchet channel. The active intercellular forces proportional to the deformation of cells can break thermodynamical equilibrium and induce the directional motion. This directional motion is caused by the ratchet effect rather than the spontaneous symmetry breaking. The motion direction is determined by the asymmetry of the channel. Remarkably, there exits an optimal nonequilibrium driving (depending on the active strength, the elasticity, and the packing fraction) at which the average velocity reaches the maximum. In addition, the optimized packing fraction and the optimized minimum width of the channel can facilitate the directional motion of cells. Our findings are relevant to understanding how macroscopic directional motion relates to the local force transmission mediated by cell-cell contacts in cellular monolayers.

Data-Driven Discovery of Active Nematic Hydrodynamics

Active nematics can be modeled using phenomenological continuum theories that account for the dynamics of the nematic director and fluid velocity through partial differential equations (PDEs). While these models provide a statistical description of the experiments, the relevant terms in the PDEs and their parameters are usually identified indirectly. We adapt a recently developed method to automatically identify optimal continuum models for active nematics directly from spatiotemporal data, via sparse regression of the coarse-grained fields onto generic low order PDEs. After extensive benchmarking, we apply the method to experiments with microtubule-based active nematics, finding a surprisingly minimal description of the system. Our approach can be generalized to gain insights into active gels, microswimmers, and diverse other experimental active matter systems.

Continuous generation of topological defects in a passively driven nematic liquid crystal

Synthetic active matter is emerging as the prime route for the realisation of biological mechanisms such as locomotion, active mixing, and self-organisation in soft materials. In particular, passive nematic complex fluids are known to form out-of-equilibrium states with topological defects, but their locomotion, activation and experimental realization has been developed and understood to only a limited extent. Here, we report that the concentration-driven flow of small molecules triggers turbulent flow in the thin film of a nematic liquid crystal that continuously generates pairs of topological defects with an integer topological charge. The diffusion results in the formation of counter-rotating vortex rolls in the liquid crystal, which above a velocity threshold transform into a turbulent flow with continuous generation and annihilation of the defect pairs. The pairs of defects are created by the self-amplifying splay instability between the vortices, until a pair of oppositely charged defects is formed.

It has been known that spontaneous defect formation and annihilation can be triggered by turbulent flows in active nematic liquid crystals. Here, Mur et al. show a complementary mechanism induced by the flow of foreign organic molecules into the liquid crystal following the concentration gradient.

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.

Defect dynamics in active polar fluids vs. active nematics

Topological defects play a key role in two-dimensional active nematics, and a transient role in two-dimensional active polar fluids. Using a variational method, we study both the transient and long-time behavior of defects in two-dimensional active polar fluids in the limit of strong order and overdamped, compressible flow, and compare the defect dynamics with the corresponding active nematics model studied recently. One result is non-central interactions between defect pairs for active polar fluids, and by extending our analysis to allow orientation dynamics of defects, we find that the orientation of +1 defects, unlike that of ±1/2 defects in active nematics, is not locked to defect positions and relaxes to asters. Moreover, using a scaling argument, we explain the transient feature of active polar defects and show that in the steady state, active polar fluids are either devoid of defects or consist of a single aster. We argue that for contractile (extensile) active nematic systems, +1 vortices (asters) should emerge as bound states of a pair of +1/2 defects, which has been recently observed. Moreover, unlike the polar case, we show that for active nematics, a linear chain of equally spaced bound states of pairs of +1/2 defects can screen the activity term. A common feature in both models is the appearance of +1 defects (elementary in polar and composite in nematic) in the steady state.

Active nematic defects in compressible and incompressible flows

We study the dynamics of active nematic films on a substrate driven by active flows with or without the incompressible constraint. Through simulations and theoretical analysis, we show that arch patterns are stable in the compressible case, while they become unstable under the incompressibility constraint. For compressible flows at high enough activity, stable arches organize themselves into a smecticlike pattern, which induce an associated global polar ordering of +1/2 nematic defects. By contrast, divergence-free flows give rise to a local nematic order of the +1/2 defects, consisting of antialigned pairs of neighboring defects, as established in previous studies.

Fluctuation-induced dynamics of nematic topological defects

Topological defects are increasingly being identified in various biological systems, where their characteristic flow fields and stress patterns are associated with continuous active stress generation by biological entities. Here, using numerical simulations of continuum fluctuating nematohydrodynamics, we show that even in the absence of any specific form of active stresses associated with self-propulsion, mesoscopic fluctuations in either orientational alignment or hydrodynamics can independently result in flow patterns around topological defects that resemble the ones observed in active systems. Our simulations further show the possibility of extensile- and contractile-like motion of fluctuation-induced positive half-integer topological defects. Remarkably, isotropic stress fields also reproduce the experimentally measured stress patterns around topological defects in epithelia. Our findings further reveal that extensile- or contractile-like flow and stress patterns around fluctuation-induced defects are governed by passive elastic stresses and flow-aligning behavior of the nematics.

Active nematic gel with quenched disorder

With quenched disorder, we introduce two-dimensional active nematics suspended in an incompressible fluid. We write the coarse-grained hydrodynamic equations of motion for slow variables, viz. density, orientation, and flow fields. The quenched disorder is introduced such that it interacts with the local orientation at every point with some strength. Disorder strength is tuned from zero to large values. We numerically study the defect dynamics and system kinetics and find that the finite disorder slows the ordering. The presence of fluid induces large fluctuation in the orientation field, further disturbing the ordering. The large fluctuation in the orientation field due to the fluid is so dominant that it reduces the effect of the quenched disorder. We have also found that the disorder effect is almost the same for both the contractile and extensile nature of active stresses in the system. This study can help to understand the impact of quenched disorder on the ordering kinetics of active gels with nematic interaction among the constituent objects.

Active Nematic Defects and Epithelial Morphogenesis

Inspired by recent experiments that highlight the role of nematic defects in the morphogenesis of epithelial tissues, we develop a minimal framework to study the dynamics of an active curved surface driven by its nematic texture. Allowing the surface to evolve via relaxational dynamics leads to a theory linking nematic defect dynamics, cellular division rates, and Gaussian curvature. Regions of large positive (negative) curvature and positive (negative) growth are colocalized with the presence of positive (negative) defects. In an ex-vivo setting of cultured murine neural progenitor cells, we show that our framework is consistent with the observed cell accumulation at positive defects and depletion at negative defects. In an in-vivo setting, we show that the defect configuration consisting of a bound +1 defect state, which is stabilized by activity, surrounded by two -1/2 defects can create a stationary ring configuration of tentacles, consistent with observations of a basal marine invertebrate Hydra.

Catapulting of topological defects through elasticity bands in active nematics

Active materials are those in which individual, uncoordinated local stresses drive the material out of equilibrium on a global scale. Examples of such assemblies can be seen across scales from schools of fish to the cellular cytoskeleton and underpin many important biological processes. Synthetic experiments that recapitulate the essential features of such active systems have been the object of study for decades as their simple rules allow us to elucidate the physical underpinnings of collective motion. One system of particular interest has been active nematic liquid crystals (LCs). Because of their well understood passive physics, LCs provide a rich platform to interrogate the effects of active stress. The flows and steady state structures that emerge in an active LCs have been understood to result from a competition between nematic elasticity and the local activity. However most investigations of such phenomena consider only the magnitude of the elastic resistance and not its peculiarities. Here we investigate a nematic liquid crystal and selectively change the ratio of the material's splay and bend elasticities. We show that increases in the nematic's bend elasticity specifically drives the material into an exotic steady state where elongated regions of acute bend distortion or "elasticity bands" dominate the structure and dynamics. We show that these bands strongly influence defect dynamics, including the rapid motion or "catapulting" along the disintegration of one of these bands thus converting bend distortion into defect transport. Thus, we report a novel dynamical state resultant from the competition between nematic elasticity and active stress.

Nematic order condensation and topological defects in inertial active nematics

Living materials at different length scales manifest active nematic features such as orientational order, nematic topological defects, and active nematic turbulence. Using numerical simulations we investigate the impact of fluid inertia on the collective pattern formation in active nematics. We show that an incremental increase in inertial effects due to reduced viscosity results in gradual melting of nematic order with an increase in topological defect density before a discontinuous transition to a vortex-condensate state. The emergent vortex-condensate state at low enough viscosities coincides with nematic order condensation within the giant vortices and the drop in the density of topological defects. We further show flow field around topological defects is substantially affected by inertial effects. Moreover, we demonstrate the strong dependence of the kinetic energy spectrum on the inertial effects, recover the Kolmogorov scaling within the vortex-condensate phase, but find no evidence of universal scaling at higher viscosities. The findings reveal complexities in active nematic turbulence and emphasize the important cross-talk between active and inertial effects in setting flow and orientational organization of active particles.

Bridging microscopic cell dynamics to nematohydrodynamics of cell monolayers

It is increasingly being realized that liquid-crystalline features can play an important role in the properties and dynamics of cell monolayers. Here, we present a cell-based model of cell layers, based on the phase-field formulation, that connects cell-cell interactions specified at the single cell level to large-scale nematic and hydrodynamic properties of the tissue. In particular, we present a minimal formulation that reproduces the well-known bend and splay hydrodynamic instabilities of the continuum nemato-hydrodynamic formulation of active matter, together with an analytical description of the instability threshold in terms of activity and elasticity of the cells. Furthermore, we provide a quantitative characterisation and comparison of flows and topological defects for extensile and contractile stress generation mechanisms, and demonstrate activity-induced heterogeneity and spontaneous formation of gaps within a confluent monolayer. Together, these results contribute to bridging the gap between cell-scale dynamics and tissue-scale collective cellular organisation.

Physics-informed graph neural networks enhance scalability of variational nonequilibrium optimal control

When a physical system is driven away from equilibrium, the statistical distribution of its dynamical trajectories informs many of its physical properties. Characterizing the nature of the distribution of dynamical observables, such as a current or entropy production rate, has become a central problem in nonequilibrium statistical mechanics. Asymptotically, for a broad class of observables, the distribution of a given observable satisfies a large deviation principle when the dynamics is Markovian, meaning that fluctuations can be characterized in the long-time limit by computing a scaled cumulant generating function. Calculating this function is not tractable analytically (nor often numerically) for complex, interacting systems, so the development of robust numerical techniques to carry out this computation is needed to probe the properties of nonequilibrium materials. Here, we describe an algorithm that recasts this task as an optimal control problem that can be solved variationally. We solve for optimal control forces using neural network ansätze that are tailored to the physical systems to which the forces are applied. We demonstrate that this approach leads to transferable and accurate solutions in two systems featuring large numbers of interacting particles.

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.

Defect loops in three-dimensional active nematics as active multipoles

We develop a description of defect loops in three-dimensional active nematics based on a multipole expansion of the far-field director and show how this leads to a self-dynamics dependent on the loop's geometric type. The dipole term leads to active stresses that generate a global self-propulsion for splay and bend loops. The quadrupole moment is nonzero only for nonplanar loops and generates a net "active torque," such that defect loops are both self-motile and self-orienting. Our analysis identifies right- and left-handed twist loops as the only force- and torque-free geometries, suggesting a mechanism for generating an excess of twist loops. Finally, we determine the Stokesian flows created by defect loops and describe qualitatively their hydrodynamics.

Assembling Microtubule-Based Active Matter

Studied for more than a century, equilibrium liquid crystals provided insight into the properties of ordered materials, and led to commonplace applications such as display technology. Active nematics are a new class of liquid crystal materials that are driven out of equilibrium by continuous motion of the constituent anisotropic units. A versatile experimental realization of active nematic liquid crystals is based on rod-like cytoskeletal filaments that are driven out of equilibrium by molecular motors. We describe protocols for assembling microtubule-kinesin based active nematic liquid crystals and associated isotropic fluids. We describe the purification of each protein and the assembly process of a two-dimensional active nematic on a water-oil interface. Finally, we show examples of nematic formation and describe methods for quantifying their non-equilibrium dynamics.

Spontaneous flow created by active topological defects

Topological defects are at the root of the large-scale organization of liquid crystals. In two-dimensional active nematics, two classes of topological defects of charges [Formula: see text] are known to play a major role due to active stresses. Despite this importance, few analytical results have been obtained on the flow-field and active-stress patterns around active topological defects. Using the generic hydrodynamic theory of active systems, we investigate the flow and stress patterns around these topological defects in unbounded, two-dimensional active nematics. Under generic assumptions, we derive analytically the spontaneous velocity and stall force of self-advected defects in the presence of both shear and rotational viscosities. Applying our formalism to the dynamics of monolayers of elongated cells at confluence, we show that the non-conservation of cell number generically increases the self-advection velocity and could provide an explanation for their observed role in cellular extrusion and multilayering. We finally investigate numerically the influence of the Ericksen stress. Our work paves the way to a generic study of the role of topological defects in active nematics, and in particular in monolayers of elongated cells.

Large-scale orientational order in bacterial colonies during inward growth

During colony growth, complex interactions regulate the bacterial orientation, leading to the formation of large-scale ordered structures, including topological defects, microdomains, and branches. These structures may benefit bacterial strains, providing invasive advantages during colonization. Active matter dynamics of growing colonies drives the emergence of these ordered structures. However, additional biomechanical factors also play a significant role during this process. Here, we show that the velocity profile of growing colonies creates strong radial orientation during inward growth when crowded populations invade a closed area. During this process, growth geometry sets virtual confinement and dictates the velocity profile. Herein, flow-induced alignment and torque balance on the rod-shaped bacteria result in a new stable orientational equilibrium in the radial direction. Our analysis revealed that the dynamics of these radially oriented structures, also known as aster defects, depend on bacterial length and can promote the survival of the longest bacteria around localized nutritional hotspots. The present results indicate a new mechanism underlying structural order and provide mechanistic insights into the dynamics of bacterial growth on complex surfaces.

Logic operations with active topological defects

Logic operations performed by semiconductor-based transistors are the basis of modern computing. There is considerable interest in creating autonomous materials systems endowed with the capability to make decisions. In this work, we introduce the concept of using topological defects in active matter to perform logic operations. When an extensile active stress in a nematic liquid crystal is turned on, +1/2 defects can self-propel, in analogy to electron transport under a voltage gradient. By relying on hydrodynamic simulations of active nematics, we demonstrate that patterns of activity, when combined with surfaces imparting certain orientations, can be used to control the formation and transport of +1/2 defects. We further show that asymmetric high- and low-activity patterns can be used to create effective defect gates, tunnels, and amplifiers. The proposed active systems offer the potential to perform computations and transmit information in active soft materials, including actin-, tubulin-, and cell-based systems.

Nano/Micromotors in Active Matter

Nano/micromotors (NMMs) are tiny objects capable of converting energy into mechanical motion. Recently, a wealth of active matter including synthetic colloids, cytoskeletons, bacteria, and cells have been used to construct NMMs. The self-sustained motion of active matter drives NMMs out of equilibrium, giving rise to rich dynamics and patterns. Alongside the spontaneous dynamics, external stimuli such as geometric confinements, light, magnetic field, and chemical potential are also harnessed to control the movements of NMMs, yielding new application paradigms of active matter. Here, we review the recent advances, both experimental and theoretical, in exploring biological NMMs. The unique dynamical features of collective NMMs are focused on, along with some possible applications of these intriguing systems.

Understanding topological defects in fluidized dry active nematics

Dense assemblies of self-propelling rods (SPRs) may exhibit fascinating collective behaviors and anomalous physical properties that are far away from equilibrium. Using large-scale Brownian dynamics simulations, we investigate the dynamics of disclination defects in 2D fluidized swarming motions of dense dry SPRs (i.e., without hydrodynamic effects) that form notable local positional topological structures that are reminiscent of smectic order. We find the deformations of smectic-like rod layers can create unique polar structures that lead to slow translations and rotations of ±1/2-order defects, which are fundamentally different from the fast streaming defect motions observed in wet active matter. We measure and characterize the statistical properties of topological defects and reveal their connections with the coherent structures. Furthermore, we construct a bottom-up active-liquid-crystal model to analyze the instability of polar lanes, which effectively leads to defect formation between interlocked polar lanes and serves as the origin of the large-scale swarming motions.

Unified analysis of topological defects in 2D systems of active and passive disks

We provide a comprehensive quantitative analysis of localized and extended topological defects in the steady state of 2D passive and active repulsive Brownian disk systems. We show that, both in and out-of-equilibrium, the passage from the solid to the hexatic is driven by the unbinding of dislocations, in quantitative agreement with the KTHNY singularity. Instead, extended clusters of defects largely dominate below the solid-hexatic critical line. The latter percolate in the liquid phase very close to the hexatic-liquid transition, both for continuous and discontinuous transitions, in the homogeneous liquid regime. At critical percolation the clusters of defects are fractal with statistical and geometric properties that are independent of the activity and compatible with the universality class of uncorrelated critical percolation. We also characterize the spatial organization of point-like defects and we show that the disclinations are not free, but rather always very near more complex defect structures. At high activity, the bulk of the dense phase generated by Motility-Induced Phase Separation is characterized by a density of point-like defects, and statistics and morphology of defect clusters, set by the amount of activity and not the packing fraction. Hexatic domains within the dense phase are separated by grain-boundaries along which a finite network of topological defects resides, interrupted by gas bubbles in cavitation. This structure is dynamic in the sense that the defect network allows for an unzipping mechanism that leaves free space for gas bubbles to appear, close, and even be released into the dilute phase.

Flow around topological defects in active nematic films

We study the active flow around isolated defects and the self-propulsion velocity of +1/2 defects in an active nematic film with both viscous dissipation (with viscosity η) and frictional damping Γ with a substrate. The interplay between these two dissipation mechanisms is controlled by the hydrodynamic dissipation length ℓd=η/Γ that screens the flows. For an isolated defect, in the absence of screening from other defects, the size of the shear vorticity around the defect is controlled by the system size R. In the presence of friction that leads to a finite value of ℓd, the vorticity field decays to zero on the lengthscales larger than ℓd. We show that the self-propulsion velocity of +1/2 defects grows with R in small systems where R<ℓd, while in the infinite system limit or when R≫ℓd, it approaches a constant value determined by ℓd.

Orientational Correlations in Active and Passive Nematic Defects

We investigate the emergence of orientational order among +1/2 disclinations in active nematic liquid crystals. Using a combination of theoretical and experimental methods, we show that +1/2 disclinations have short-range antiferromagnetic alignment, as a consequence of the elastic torques originating from their polar structure. The presence of intermediate -1/2 disclinations, however, turns this interaction from antialigning to aligning at scales that are smaller than the typical distance between like-sign defects. No long-range orientational order is observed. Strikingly, these effects are insensitive to material properties and qualitatively similar to what is found for defects in passive nematic liquid crystals.

Multiscale Microtubule Dynamics in Active Nematics

In microtubule-based active nematics, motor-driven extensile motion of microtubule bundles powers chaotic large-scale dynamics. We quantify the interfilament sliding motion both in isolated bundles and in a dense active nematic. The extension speed of an isolated microtubule pair is comparable to the molecular motor stepping speed. In contrast, the net extension in dense 2D active nematics is significantly slower; the interfilament sliding speeds are widely distributed about the average and the filaments exhibit both contractile and extensile relative motion. These measurements highlight the challenge of connecting the extension rate of isolated bundles to the multimotor and multifilament interactions present in a dense 2D active nematic. They also provide quantitative data that is essential for building multiscale models.

Behavior of chiral active nematics confined to nanoscopic circular region

We performed molecular dynamic simulations of a model active nematic confined to a two-dimensional nanoscopic circular region under both tangential and radial anchoring boundary conditions. This active material is assumed to be composed of elongated chiral particles which interact with each other by means of isotropic Lennard-Jones and anisotropic Maier-Saupe-like potentials. These particles have the lateral appendage emitting a jet of some substance generated by a certain internal chemical reaction. As a result, such elongated particles are exposed to both the reactive self-propelled force and the torque that provide an additional translational movement of particles and a self-rotation with respect to their geometric centers. It has been found that the chiral active nematics under consideration form time-dependent vortex-like structures with two +1/2 topological defects which are similar to experimentally observed structures in active materials.

Topological braiding and virtual particles on the cell membrane

Braiding of topological structures in complex matter fields provides a robust framework for encoding and processing information, and it has been extensively studied in the context of topological quantum computation. In living systems, topological defects are crucial for the localization and organization of biochemical signaling waves, but their braiding dynamics remain unexplored. Here, we show that the spiral wave cores, which organize the Rho-GTP protein signaling dynamics and force generation on the membrane of starfish egg cells, undergo spontaneous braiding dynamics. Experimentally measured world line braiding exponents and topological entropy correlate with cellular activity and agree with predictions from a generic field theory. Our analysis further reveals the creation and annihilation of virtual quasi-particle excitations during defect scattering events, suggesting phenomenological parallels between quantum and living matter.

Defect Superdiffusion and Unbinding in a 2D XY Model of Self-Driven Rotors

We consider a nonequilibrium extension of the 2D XY model, equivalent to the noisy Kuramoto model of synchronization with short-range coupling, where rotors sitting on a square lattice are self-driven by random intrinsic frequencies. We study the static and dynamic properties of topological defects (vortices) and establish how self-spinning affects the Berezenskii-Kosterlitz-Thouless phase transition scenario. The nonequilibrium drive breaks the quasi-long-range ordered phase of the 2D XY model into a mosaic of ordered domains of controllable size and results in self-propelled vortices that generically unbind at any temperature, featuring superdiffusion ⟨r^{2}(t)⟩∼t^{3/2} with a Gaussian distribution of displacements. Our work provides a simple framework to investigate topological defects in nonequilibrium matter and sheds new light on the problem of synchronization of locally coupled oscillators.

Properties of twisted topological defects in 2D nematic liquid crystals

Topological defects are one of the most conspicuous features of liquid crystals. In two dimensional nematics, they have been shown to behave effectively as particles with both charge and orientation, which dictate their interactions. Here, we study "twisted" defects that have a radially dependent orientation. We find that twist can be partially relaxed through the creation and annihilation of defect pairs. By solving the equations for defect motion and calculating the forces on defects, we identify four distinct elements that govern the relative relaxational motion of interacting topological defects, namely attraction, repulsion, co-rotation and co-translation. The interaction of these effects can lead to intricate defect trajectories, which can be controlled by setting relevant timescales.

Active liquid crystals powered by force-sensing DNA-motor clusters

Cytoskeletal active nematics exhibit striking nonequilibrium dynamics that are powered by energy-consuming molecular motors. To gain insight into the structure and mechanics of these materials, we design programmable clusters in which kinesin motors are linked by a double-stranded DNA linker. The efficiency by which DNA-based clusters power active nematics depends on both the stepping dynamics of the kinesin motors and the chemical structure of the polymeric linker. Fluorescence anisotropy measurements reveal that the motor clusters, like filamentous microtubules, exhibit local nematic order. The properties of the DNA linker enable the design of force-sensing clusters. When the load across the linker exceeds a critical threshold, the clusters fall apart, ceasing to generate active stresses and slowing the system dynamics. Fluorescence readout reveals the fraction of bound clusters that generate interfilament sliding. In turn, this yields the average load experienced by the kinesin motors as they step along the microtubules. DNA-motor clusters provide a foundation for understanding the molecular mechanism by which nanoscale molecular motors collectively generate mesoscopic active stresses, which in turn power macroscale nonequilibrium dynamics of active nematics.

Defect Spirograph: Dynamical Behavior of Defects in Spatially Patterned Active Nematics

Topological defects in active liquid crystals can be confined by introducing gradients of activity. Here, we examine the dynamical behavior of two defects confined by a sharp gradient of activity that separates an active circular region and a surrounding passive nematic material. Continuum simulations are used to explain how the interplay among energy injection into the system, hydrodynamic interactions, and frictional forces governs the dynamics of topologically required self-propelling +1/2 defects. Our findings are rationalized in terms of a phase diagram for the dynamical response of defects in terms of activity and frictional damping strength. Different regions of the underlying phase diagram correspond to distinct dynamical modes, namely immobile defects, steady rotation of defects, bouncing defects, bouncing-cruising defects, dancing defects, and multiple defects with irregular dynamics. These dynamic states raise the prospect of generating synchronized defect arrays for microfluidic applications.