Topological structure dynamics revealing collective evolution in active nematics

A machine learning approach to robustly determine director fields and analyze defects in active nematics

Active nematics are dense systems of rodlike particles that consume energy to drive motion at the level of the individual particles. They exist in natural systems like biological tissues and artificial materials such as suspensions of self-propelled colloidal particles or synthetic microswimmers. Active nematics have attracted significant attention in recent years due to their spectacular nonequilibrium collective spatiotemporal dynamics, which may enable applications in fields such as robotics, drug delivery, and materials science. The director field, which measures the direction and degree of alignment of the local nematic orientation, is a crucial characteristic of active nematics and is essential for studying topological defects. However, determining the director field is a significant challenge in many experimental systems. Although director fields can be derived from images of active nematics using traditional imaging processing methods, the accuracy of such methods is highly sensitive to the settings of the algorithms. These settings must be tuned from image to image due to experimental noise, intrinsic noise of the imaging technology, and perturbations caused by changes in experimental conditions. This sensitivity currently limits automatic analysis of active nematics. To address this, we developed a machine learning model for extracting reliable director fields from raw experimental images, which enables accurate analysis of topological defects. Application of the algorithm to experimental data demonstrates that the approach is robust and highly generalizable to experimental settings that are different from those in the training data. It could be a promising tool for investigating active nematics and may be generalized to other active matter systems.

Symmetry Breaking of Self-Propelled Topological Defects in Thin-Film Active Chiral Nematics

Active nematics represent a range of dense active matter systems which can engender spontaneous flows and self-propelled topological defects. Two-dimensional (2D) active nematic theory and simulation have been successful in explaining many quasi-2D experiments in which self-propelled +1/2 defects are observed to move along their symmetry axis. However, many active liquid crystals are essentially chiral nematic, but their twist mode becomes irrelevant under the 2D assumption. Here, we use theory and simulation to examine a three-dimensional active chiral nematic confined to a thin film, thus forming a quasi-2D system. We predict that the self-propelled +1/2 disclination in a curved thin film can break its mirror symmetry by moving circularly. Our prediction is confirmed by hydrodynamic simulations of thin spherical-shell and thin cylindrical-shell systems. In the spherical-shell confinement, the four emerged +1/2 disclinations exhibit rich dynamics as a function of activity and chirality. As such, we have proposed a new symmetry-breaking scenario in which self-propelled defects in quasi-2D active nematics can acquire an active angular velocity, greatly enriching their dynamics for finer control and emerging applications.

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.

Mesoscopic simulations of active nematics

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

Motile Topological Defects Hinder Dynamical Arrest in Dense Liquids of Active Ellipsoids

Recent numerical studies have identified the persistence time of active motion as a critical parameter governing glassy dynamics in dense active matter. Here we studied dynamics in liquids of granular active ellipsoids with tunable persistence and velocity. We show that increasing the persistence time at moderate supercooling is equivalent to increasing the strength of attraction in equilibrium liquids and results in reentrant dynamics not just in the translational degrees of freedom, as anticipated, but also in the orientational ones. However, at high densities, motile topological defects, unique to active liquids of elongated particles, hindered dynamical arrest. Most remarkably, for the highest activity, we observed intermittent dynamics due to the jamming-unjamming of these defects for the first time.

Absolute negative mobility of active polymer chains in steady laminar flows

We investigate the transport of active polymer chains in steady laminar flows in the presence of thermal noise and an external constant force. In the model, the polymer chain is worm-like and is propelled by active forces along its tangent vectors. Compared with inertial Brownian particles, active polymer chains in steady laminar flows exhibit richer movement patterns due to their specific spatial structures. The simulation results show that the velocity-force relation is strongly dependent on the system parameters such as the chain length, bending rigidity, active force and so on. The polymer chain may move in some preferential movement directions and exhibits absolute negative mobility within appropriate parameter regimes, i.e., the polymer chain can move in a direction opposite to the external constant force. In particular, we can observe giant negative mobility in a broad range of parameter regimes.

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.

Topology of Orientational Defects in Confined Smectic Liquid Crystals

We propose a general formalism to characterize orientational frustration of smectic liquid crystals in confinement by interpreting the emerging networks of grain boundaries as objects with a topological charge. In a formal idealization, this charge is distributed in pointlike units of quarter-integer magnitude, which we identify with tetratic disclinations located at the end points and nodes. This coexisting nematic and tetratic order is analyzed with the help of extensive Monte Carlo simulations for a broad range of two-dimensional confining geometries as well as colloidal experiments, showing how the observed defect networks can be universally reconstructed from simple building blocks. We further find that the curvature of the confining wall determines the anchoring behavior of grain boundaries, such that the number of nodes in the emerging networks and the location of their end points can be tuned by changing the number and smoothness of corners, respectively.

Molecular Simulation Studies on the Interactions of Bilirubin at Different States with a Lipid Bilayer

The unconjugated bilirubin (BR) may penetrate through the cell membrane and cause a severe cytotoxicity. However, the molecular mechanism underlying the penetration of BR into the cell membrane is still largely unknown. In this work, we systematically investigate the interaction of BR and a lipid bilayer under different conditions by using all-atom molecular dynamics simulations. It is found that BR at the Z,Z conformation can easily enter into the interior of the lipid bilayer due to its hydrophobicity. However, when BR transforms from the Z,Z conformation to the E,E conformation (after the blue-light emission), its penetration ability is greatly reduced (especially at its ionized state). This study may offer useful physical insights into the effect of phototherapy on the penetration behavior and the cytotoxicity of the unconjugated BR.

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.

Absorption of self-propelled particles into a dense porous medium

We study the absorption of self-propelled particles into a finite-size dense porous medium, which is mimicked by an obstacle array. We find that, depending on the competition of the propelling strength versus the repulsive barrier formed by obstacles and the contrast between the characteristic time scales of permeation and propelling persistence, the absorption process exhibits three distinct types of behavior. In Type I and II behavior, the propelling strength is not large enough to surmount the barrier, and hence particles transport in the medium by barrier-hopping dynamics. The initial permeation of particles toward the medium center is phenomenologically similar to a normal slow diffusion process. But, surprisingly, after the initial permeation process, a concentrated nucleus of particle aggregates forms and grows at the medium center in Type I, due to the long propelling persistence. Such an abnormal "nucleation" phenomenon does not appear in Type II, in which the propelling persistence is low. When the propelling strength is very high (Type III), particles transport smoothly in the medium, hence the initial slow diffusion process disappears and small particle clusters form and merge randomly in the medium. Our results provide a foundation for applications of active objects in a complex environment and also suggest the possible usage of a porous medium, for example, in the selection or sorting of active matter.

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.

Barrier-controlled nonequilibrium criticality in reactive particle systems

Nonequilibrium critical phenomena generally exist in many dynamic systems, like chemical reactions and some driven-dissipative reactive particle systems. Here, by using computer simulation and theoretical analysis, we demonstrate the crucial role of the activation barrier on the criticality of dynamic phase transitions in a minimal reactive hard-sphere model. We find that at zero thermal noise, with increasing the activation barrier, the type of transition changes from a continuous conserved directed percolation into a discontinuous dynamic transition by crossing a tricritical point. A mean-field theory combined with field simulation is proposed to explain this phenomenon. The possibility of Ising-type criticality in the nonequilibrium system at finite thermal noise is also discussed.

Alignment of a topological defect by an activity gradient

As a method for controlling active materials, researchers have suggested designing patterns of activity on a substrate, which should guide the motion of topological defects. To investigate this concept, we model the behavior of a single defect of topological charge +1/2, moving in an activity gradient. This modeling uses three methods: (1) approximate analytic solution of hydrodynamic equations, (2) macroscopic, symmetry-based theory of the defect as an effective oriented particle, and (3) numerical simulation. All three methods show that an activity gradient aligns the defect orientation, and hence should be useful to control defect motion.

Memory effects, arches and polar defect ordering at the cross-over from wet to dry active nematics

We use analytic arguments and numerical solutions of the continuum, active nematohydrodynamic equations to study how friction alters the behaviour of active nematics. Concentrating on the case where there is nematic ordering in the passive limit, we show that, as the friction is increased, memory effects become more prominent and +1/2 topological defects leave increasingly persistent trails in the director field as they pass. The trails are preferential sites for defect formation and they tend to impose polar order on any new +1/2 defects. In the absence of noise and for high friction, it becomes very difficult to create defects, but trails formed by any defects present at the beginning of the simulations persist and organise into parallel arch-like patterns in the director field. We show aligned arches of equal width are approximate steady state solutions of the equations of motion which co-exist with the nematic state. We compare our results to other models in the literature, in particular dry systems with no hydrodynamics, where trails, arches and polar defect ordering have also been observed.

Nonequilibrium pattern formation in circularly confined two-dimensional systems with competing interactions

We numerically investigate the nonequilibrium behaviors of classic particles with competing interactions confined in a two-dimensional logarithmic trap. We reveal a quench-induced surprising dynamics exhibiting rich dynamic patterns depending upon confinement strength and trap size, which is attributed to the time-dependent competition between interparticle repulsions and attractions under a circular confinement. Moreover, in the collectively diffusive motions of the particles, we find that the emergence of dynamic structure transformation coincides with a diffusive mode transition from superdiffusion to subdiffusion. These findings are likely useful in understanding the pattern selection and evolution in various chemical and biological systems in addition to modulated systems, and add a new route to tailoring the morphology of pattern-forming systems.

Machine learning forecasting of active nematics

Active nematics are a class of far-from-equilibrium materials characterized by local orientational order of force-generating, anisotropic constitutes. Traditional methods for predicting the dynamics of active nematics rely on hydrodynamic models, which accurately describe idealized flows and many of the steady-state properties, but do not capture certain detailed dynamics of experimental active nematics. We have developed a deep learning approach that uses a Convolutional Long-Short-Term-Memory (ConvLSTM) algorithm to automatically learn and forecast the dynamics of active nematics. We demonstrate our purely data-driven approach on experiments of 2D unconfined active nematics of extensile microtubule bundles, as well as on data from numerical simulations of active nematics.

Preferred penetration of active nano-rods into narrow channels and their clustering

In a channel connected to a reservoir, passive particles prefer staying in the reservoir than the channel due to the entropic effect, as the size of the particles is comparable to that of the channel. Self-propelled rods can exhibit out-of-equilibrium phenomena, and their partition behavior may differ from that of passive rods due to their persistent swimming ability. In this work, the distribution of active nano-rods between the nanoscale channel and reservoir is explored using dissipative particle dynamics. The ratio of the nano-rod concentration in the slit to that in the reservoir, defined as the partition ratio Ψ, is a function of active force, channel width, and rod length. Although passive nano-rods prefer staying in bulk (Ψ < 1), active rods can overcome the entropic barrier and show favorable partition toward narrow channels (Ψ > 1). As the slit width decreases to about the rod's width, active rods entering the slit behave like a quasi-two-dimensional system dynamically. At sufficiently high concentrations and Peclet numbers, nano-rods tend to align and move together in the same direction for a certain time. The distribution (PM) of the cluster size (M) follows a power law, PM ∝ M-2, for small clusters.

Active nematics with quenched disorder

We introduce a two-dimensional active nematic with quenched disorder. We write the coarse-grained hydrodynamic equations of motion for slow variables, viz. density and orientation. Disorder strength is tuned from zero to large values. Results from the numerical solution of equations of motion as well as the calculation of two-point orientation correlation function using linear approximation shows that the ordered steady state follows a disorder dependent crossover from quasi-long-range order to short-range order. Such crossover is due to the pinning of ±1/2 topological defects in the presence of finite disorder, which breaks the system in uncorrelated domains. Finite disorder slows the dynamics of +1/2 defect, and it leads to slower growth dynamics. The two-point correlation functions for the density and orientation fields show good dynamic scaling but no static scaling for the different disorder strengths. Our findings can motivate experimentalists to verify the results and find applications in living and artificial apolar systems in the presence of a quenched disorder.

Controlling ion transport in a C2N-based nanochannel with tunable interlayer spacing

Selective ion transport through a nanochannel formed by stacked two-dimensional materials plays a key role in water desalination, nanofiltration, and ion separation. Although there have been many functional nanomaterials used in these applications, how to well control ion transport in a laminar structure so as to obtain the desired selectivity still remains a challenging problem. In the present work, the transport of ions through a C₂N-based nanochannel is investigated by using all-atom molecular dynamics simulation. It is found that C₂N-based nanochannels with different interlayer spacing posses diverse ion selectivity, which is mainly attributed to the distinct loading capability among ions and the different velocity of ions inside the nanochannel. Moreover, we also find that the ion selectivity is dependent on the electric field, but nearly independent of the salt concentration. The present study may provide some physical insights into the experimental design of C₂N-based nanodevices in nanofiltration.

Transport of self-propelled particles across a porous medium: trapping, clogging, and the Matthew effect

We study the transport of self-propelled particles from one free chamber to another across two stripe-like areas of dense porous medium. The medium is mimicked by arrays of obstacles. We find that active motion could greatly speed up the transport of particles. However, more and more particles become trapped in the obstacle arrays with the enhancement of activity. At high persistence (low rotational diffusion rate) and moderate particle concentration, we observe the Matthew effect in the aggregation of particles in the two obstacle arrays. This effect is weakened by introduction of randomness or deformability into the obstacle arrays. Moreover, the dependence on deformability shows the characteristics of first-order phase transition. In rare situations, the system could be stuck in a dynamic unstable state, e.g. the particles alternatively gather more in one of the two obstacle arrays, exhibiting oscillation of particle number between the arrays. Our results reveal new features in the transport of active objects in a complex medium and have implications for manipulating their collective assembly.

Collective Dynamics of Model Pili-Based Twitcher-Mode Bacilliforms

Pseudomonas aeruginosa, like many bacilliforms, are not limited only to swimming motility but rather possess many motility strategies. In particular, twitching-mode motility employs hair-like pili to transverse moist surfaces with a jittery irregular crawl. Twitching motility plays a critical role in redistributing cells on surfaces prior to and during colony formation. We combine molecular dynamics and rule-based simulations to study twitching-mode motility of model bacilliforms and show that there is a critical surface coverage fraction at which collective effects arise. Our simulations demonstrate dynamic clustering of twitcher-type bacteria with polydomains of local alignment that exhibit spontaneous correlated motions, similar to rafts in many bacterial communities.

Filamentous active matter: Band formation, bending, buckling, and defects

Motor proteins drive persistent motion and self-organization of cytoskeletal filaments. However, state-of-the-art microscopy techniques and continuum modeling approaches focus on large length and time scales. Here, we perform component-based computer simulations of polar filaments and molecular motors linking microscopic interactions and activity to self-organization and dynamics from the filament level up to the mesoscopic domain level. Dynamic filament cross-linking and sliding and excluded-volume interactions promote formation of bundles at small densities and of active polar nematics at high densities. A buckling-type instability sets the size of polar domains and the density of topological defects. We predict a universal scaling of the active diffusion coefficient and the domain size with activity, and its dependence on parameters like motor concentration and filament persistence length. Our results provide a microscopic understanding of cytoplasmic streaming in cells and help to develop design strategies for novel engineered active materials.

Vortex formation of spherical self-propelled particles around a circular obstacle

A vortex is a common ratchet phenomenon in active systems. The spatial symmetry is usually broken by introducing asymmetric shapes or spontaneously by collective motion in the presence of hydrodynamic interactions or other alignment effects. Unexpectedly, we observe, by simulations, the formation of a vortex in the simplest model of a circular obstacle immersed in a bath of spherical self-propelled particles. No symmetry-breaking factors mentioned above are included in this model. The vortex forms only when the particle activity is high, i.e. large persistence. The obstacle size is also a key factor and the vortex only forms in a limited range of obstacle sizes. The sustainment of the vortex originates from the bias of the rotating particle cluster around the obstacle in accepting the incoming particles based on their propelling directions. Our results provide new understanding of and insights into the spontaneous symmetry-breaking and ratchet phenomena in active matter.

Topological structure and dynamics of three-dimensional active nematics

Topological structures are effective descriptors of the nonequilibrium dynamics of diverse many-body systems. For example, motile, point-like topological defects capture the salient features of two-dimensional active liquid crystals composed of energy-consuming anisotropic units. We dispersed force-generating microtubule bundles in a passive colloidal liquid crystal to form a three-dimensional active nematic. Light-sheet microscopy revealed the temporal evolution of the millimeter-scale structure of these active nematics with single-bundle resolution. The primary topological excitations are extended, charge-neutral disclination loops that undergo complex dynamics and recombination events. Our work suggests a framework for analyzing the nonequilibrium dynamics of bulk anisotropic systems as diverse as driven complex fluids, active metamaterials, biological tissues, and collections of robots or organisms.

Understanding Dense Active Nematics from Microscopic Models

We study dry, dense active nematics at both particle and continuous levels. Specifically, extending the Boltzmann-Ginzburg-Landau approach, we derive well-behaved hydrodynamic equations from a Vicsek-style model with nematic alignment and pairwise repulsion. An extensive study of the phase diagram shows qualitative agreement between the two levels of description. We find in particular that the dynamics of topological defects strongly depends on parameters and can lead to "arch" solutions forming a globally polar, smecticlike arrangement of Néel walls. We show how these configurations are at the origin of the defect ordered states reported previously. This work offers a detailed understanding of the theoretical description of dense active nematics directly rooted in their microscopic dynamics.

Defect dynamics in clusters of self-propelled rods in circular confinement

Rod-shaped active micro/nano-particles, such as bacterial and bipolar metallic micro/nano-motors, demonstrate novel collective phenomena far from the equilibrium state compared to passive particles. We apply a simulation approach --dissipative particle dynamics (DPD)-- to explore the collectively ordered states of self-propelled rods (SPRs). The SPRs are confined in a finite circular zone and repel each other when two rods touch each other. It is found that for a long enough rods system, the global vortex patterns, dynamic pattern oscillation between hedgehog pattern and vortex pattern, and hedgehog patterns are observed successively with increasing active force F_a. For the vortex pattern, the total interaction energy between the rods U is linear with active force F_a, i.e., U ∼ F_a . While the relation U ∼ F_a² is obtained for the hedgehog structure. It is observed that a new hedgehog pattern with one defect core is created by two ejections of polar cluster in opposite directions from the original hedgehog pattern, and then merges into one through the diffusion of the two aggregates, i.e., the creation and annihilation of topological charges.

Advances in Biological Liquid Crystals

Biological liquid crystals, a rich set of soft materials with rod-like structures widely existing in nature, possess typical lyotropic liquid crystalline phase properties both in vitro (e.g., cellulose, peptides, and protein assemblies) and in vivo (e.g., cellular lipid membrane, packed DNA in bacteria, and aligned fibroblasts). Given the ability to undergo phase transition in response to various stimuli, numerous practices are exercised to spatially arrange biological liquid crystals. Here, a fundamental understanding of interactions between rod-shaped biological building blocks and their orientational ordering across multiple length scales is addressed. Discussions are made with regard to the dependence of physical properties of nonmotile objects on the first-order phase transition and the coexistence of multi-phases in passive liquid crystalline systems. This work also focuses on how the applied physical stimuli drives the reorganization of constituent passive particles for a new steady-state alignment. A number of recent progresses in the dynamics behaviors of active liquid crystals are presented, and particular attention is given to those self-propelled animate elements, like the formation of motile topological defects, active turbulence, correlation of orientational ordering, and cellular functions. Finally, future implications and potential applications of the biological liquid crystalline materials are discussed.

Unfolding of a diblock chain and its anomalous diffusion induced by active particles

We study the structural and dynamical behavior of an A-B diblock chain in the bath of active Brownian particles (ABPs) by Brownian dynamics simulations in two dimensions. We are interested in the situation that the effective interaction between the A segments is attractive, while that between the B segments is repulsive. Therefore, in thermal (nonactive) equilibrium, the A block "folds" into a compact globule, while the B block is in the expanded coil state. Interestingly, we find that the A block could "unfold" sequentially like unknitting a sweater, driven by the surrounding ABPs when the propelling strength on them is beyond a certain value. This threshold value decreases and then levels off as the length of the B block increases. We also find a simple power-law relation between the unfolding time of the A block and the self-propelling strength and an exponential relation between the unfolding time and the length of the B block. Finally, we probe the translational and rotational diffusion of the chain and find that both of them show "super-diffusivity" in a large time window, especially when the self-propelling strength is small and the A block is in the folded state. Such super-diffusivity is due to the strong asymmetric distribution of ABPs around the chain. Our work provides new insights into the behavior of a polymer chain in the environment of active objects.

Statistical properties of autonomous flows in 2D active nematics

We study the dynamics of a tunable 2D active nematic liquid crystal composed of microtubules and kinesin motors confined to an oil-water interface. Kinesin motors continuously inject mechanical energy into the system through ATP hydrolysis, powering the relative microscopic sliding of adjacent microtubules, which in turn generates macroscale autonomous flows and chaotic dynamics. We use particle image velocimetry to quantify two-dimensional flows of active nematics and extract their statistical properties. In agreement with the hydrodynamic theory, we find that the vortex areas comprising the chaotic flows are exponentially distributed, which allows us to extract the characteristic system length scale. We probe the dependence of this length scale on the ATP concentration, which is the experimental knob that tunes the magnitude of the active stress. Our data suggest a possible mapping between the ATP concentration and the active stress that is based on the Michaelis-Menten kinetics that governs the motion of individual kinesin motors.

Data-driven quantitative modeling of bacterial active nematics

Active matter comprises individual units that convert energy into mechanical motion. In many examples, such as bacterial systems and biofilament assays, constituent units are elongated and can give rise to local nematic orientational order. Such "active nematics" systems have attracted much attention from both theorists and experimentalists. However, despite intense research efforts, data-driven quantitative modeling has not been achieved, a situation mainly due to the lack of systematic experimental data and to the large number of parameters of current models. Here, we introduce an active nematics system made of swarming filamentous bacteria. We simultaneously measure orientation and velocity fields and show that the complex spatiotemporal dynamics of our system can be quantitatively reproduced by a type of microscopic model for active suspensions whose important parameters are all estimated from comprehensive experimental data. This provides unprecedented access to key effective parameters and mechanisms governing active nematics. Our approach is applicable to different types of dense suspensions and shows a path toward more quantitative active matter research.

The interplay between activity and filament flexibility determines the emergent properties of active nematics

Active nematics are microscopically driven liquid crystals that exhibit dynamical steady states characterized by the creation and annihilation of topological defects. Motivated by differences between previous simulations of active nematics based on rigid rods and experimental realizations based on semiflexible biopolymer filaments, we describe a large-scale simulation study of a particle-based computational model that explicitly incorporates filament semiflexibility. We find that energy injected into the system at the particle scale preferentially excites bend deformations, reducing the apparent filament bend modulus. The emergent characteristics of the active nematic depend on activity and flexibility only through this activity-renormalized bend 'modulus', demonstrating that apparent values of material parameters, such as the Frank 'constants', depend on activity. Thus, phenomenological parameters within continuum hydrodynamic descriptions of active nematics must account for this dependence. Further, we present a systematic way to estimate these parameters from observations of deformation fields and defect shapes in experimental or simulation data.

Defect Unbinding in Active Nematics

We formulate the statistical dynamics of topological defects in the active nematic phase, formed in two dimensions by a collection of self-driven particles on a substrate. An important consequence of the nonequilibrium drive is the spontaneous motility of strength +1/2 disclinations. Starting from the hydrodynamic equations of active nematics, we derive an interacting particle description of defects that includes active torques. We show that activity, within perturbation theory, lowers the defect-unbinding transition temperature, determining a critical line in the temperature-activity plane that separates the quasi-long-range ordered (nematic) and disordered (isotropic) phases. Below a critical activity, defects remain bound as rotational noise decorrelates the directed dynamics of +1/2 defects, stabilizing the quasi-long-range ordered nematic state. This activity threshold vanishes at low temperature, leading to a reentrant transition. At large enough activity, active forces always exceed thermal ones and the perturbative result fails, suggesting that in this regime activity will always disorder the system. Crucially, rotational diffusion being a two-dimensional phenomenon, defect unbinding cannot be described by a simplified one-dimensional model.

Active nematics

Active matter extracts energy from its surroundings at the single particle level and transforms it into mechanical work. Examples include cytoskeleton biopolymers and bacterial suspensions. Here, we review experimental, theoretical and numerical studies of active nematics - a type of active system that is characterised by self-driven units with elongated shape. We focus primarily on microtubule-kinesin mixtures and the hydrodynamic theories that describe their properties. An important theme is active turbulence and the associated motile topological defects. We discuss ways in which active turbulence may be controlled, a pre-requisite to harvesting energy from active materials, and we consider the appearance, and possible implications, of active nematics and topological defects to cellular systems and biological processes.

Pair creation, motion, and annihilation of topological defects in two-dimensional nematic liquid crystals

We present a framework for the study of disclinations in two-dimensional active nematic liquid crystals and topological defects in general. The order tensor formalism is used to calculate exact multiparticle solutions of the linearized static equations inside a planar uniformly aligned state so that the total charge has to vanish. Topological charge conservation then requires that there is always an equal number of q=1/2 and q=-1/2 charges. Starting from a set of hydrodynamic equations, we derive a low-dimensional dynamical system for the parameters of the static solutions, which describes the motion of a half-disclination pair or of several pairs. Within this formalism, we model defect production and annihilation, as observed in experiments. Our dynamics also provide an estimate for the critical density at which production and annihilation rates are balanced.

Dynamical patterns in nematic active matter on a sphere

Using simulations of self-propelled agents with short-range repulsion and nematic alignment, we explore the dynamical phases of a dense active nematic confined to the surface of a sphere. We map the nonequilibrium phase diagram as a function of curvature, alignment strength, and activity. Our model reproduces several phases seen in recent experiments on active microtubule bundles confined the surfaces of vesicles. At low driving, we recover the equilibrium nematic ground state with four +1/2 defects. As the driving is increased, geodesic forces drive the transition to a polar band wrapping around an equator, with large empty spherical caps corresponding to two +1 defects at the poles. Upon further increasing activity, the bands fold onto themselves, and the system eventually transitions to a turbulent state marked by the proliferation of pairs of topological defects. We highlight the key role of the nematic persistence length in controlling pattern formation in these confined systems with positive Gaussian curvature.

Dynamically generated patterns in dense suspensions of active filaments

We use Langevin dynamics simulations to study dynamical behavior of a dense planar layer of active semiflexible filaments. Using the strength of active force and the thermal persistence length as parameters, we map a detailed phase diagram and identify several nonequilibrium phases in this system. In addition to a slowly flowing melt phase, we observe that, for sufficiently high activity, collective flow accompanied by signatures of local polar and nematic order appears in the system. This state is also characterized by strong density fluctuations. Furthermore, we identify an activity-driven crossover from this state of coherently flowing bundles of filaments to a phase with no global flow, formed by individual filaments coiled into rotating spirals. This suggests a mechanism where the system responds to activity by changing the shape of active agents, an effect with no analog in systems of active particles without internal degrees of freedom.

Control of active turbulence through addressable soft interfaces

We present an experimental study of a kinesin/tubulin active nematic formed at different oil interfaces. By tuning the interfacial rheology of the contacting oil, we have been able to condition and control the seemingly chaotic motion that characterizes the self-sustained active flows in our preparations. The active nematic is inherently unstable and spontaneously develops defects from an initial homogeneous state. We show that the steady state and, in particular, the density and dynamics of the defects strongly depends on the rheology of the contacting oil. Using a smectic-A thermotropic liquid crystal as the oil phase, we pattern the interface thanks to the anisotropy of the shear viscosity in this material. The geometry of the active nematic adapts to the boundary conditions at the interface by changing from the so-called active turbulent regime to laminar flows along the easy flow directions. The latter can be either a lattice of self-assembled circular paths or reconfigurable homogeneous orientations that can be addressed by means of an external magnetic field. We show that, under all confinement conditions, the spatiotemporal modes exhibited by the active liquid are consistent with a single intrinsic length scale, which can be tuned by the material parameters, and obey basic topological requirements imposed on the defects that drive the active flows. Future control strategies, including a tunable depleting agent, are discussed.

Taming active turbulence with patterned soft interfaces

Active matter embraces systems that self-organize at different length and time scales, often exhibiting turbulent flows apparently deprived of spatiotemporal coherence. Here, we use a layer of a tubulin-based active gel to demonstrate that the geometry of active flows is determined by a single length scale, which we reveal in the exponential distribution of vortex sizes of active turbulence. Our experiments demonstrate that the same length scale reemerges as a cutoff for a scale-free power law distribution of swirling laminar flows when the material evolves in contact with a lattice of circular domains. The observed prevalence of this active length scale can be understood by considering the role of the topological defects that form during the spontaneous folding of microtubule bundles. These results demonstrate an unexpected strategy for active systems to adapt to external stimuli, and provide with a handle to probe the existence of intrinsic length and time scales.

Active nematics consist of self-driven components that develop orientational order and turbulent flow. Here Guillamat et al. investigate an active nematic constrained in a quasi-2D geometrical setup and show that there exists an intrinsic length scale that determines the geometry in all forcing regimes.

On the mechanism of long-range orientational order of fibroblasts

Long-range alignment ordering of fibroblasts have been observed in the vicinity of cancerous tumors and can be recapitulated with in vitro experiments. However, the mechanisms driving their ordering are not understood. Here, we show that local collision-driven nematic alignment interactions among fibroblasts are insufficient to explain observed long-range alignment. One possibility is that there exists another orientation field coevolving with the cells and reinforcing their alignment. We propose that this field reflects the mechanical cross-talk between the fibroblasts and the underlying fibrous material on which they move. We show that this long-range interaction can give rise to high nematic order and to the observed patterning of the cancer microenvironment.

Order-disorder transition in active nematic: A lattice model study

We introduce a lattice model for active nematic composed of self-propelled apolar particles, study its different ordering states in the density-temperature parameter space, and compare with the corresponding equilibrium model. The active particles interact with their neighbours within the framework of the Lebwohl-Lasher model, and move anisotropically along their orientation to an unoccupied nearest neighbour lattice site. An interplay of the activity, thermal fluctuations and density gives rise distinct states in the system. For a fixed temperature, the active nematic shows a disordered isotropic state, a locally ordered inhomogeneous mixed state, and bistability between the inhomogeneous mixed and a homogeneous globally ordered state in different density regime. In the low temperature regime, the isotropic to the inhomogeneous mixed state transition occurs with a jump in the order parameter at a density less than the corresponding equilibrium disorder-order transition density. Our analytical calculations justify the shift in the transition density and the jump in the order parameter. We construct the phase diagram of the active nematic in the density-temperature plane.

Topological defects control collective dynamics in neural progenitor cell cultures

Cultured stem cells have become a standard platform not only for regenerative medicine and developmental biology but also for biophysical studies. Yet, the characterization of cultured stem cells at the level of morphology and of the macroscopic patterns resulting from cell-to-cell interactions remains largely qualitative. Here we report on the collective dynamics of cultured murine neural progenitor cells (NPCs), which are multipotent stem cells that give rise to cells in the central nervous system. At low densities, NPCs moved randomly in an amoeba-like fashion. However, NPCs at high density elongated and aligned their shapes with one another, gliding at relatively high velocities. Although the direction of motion of individual cells reversed stochastically along the axes of alignment, the cells were capable of forming an aligned pattern up to length scales similar to that of the migratory stream observed in the adult brain. The two-dimensional order of alignment within the culture showed a liquid-crystalline pattern containing interspersed topological defects with winding numbers of +1/2 and -1/2 (half-integer due to the nematic feature that arises from the head-tail symmetry of cell-to-cell interaction). We identified rapid cell accumulation at +1/2 defects and the formation of three-dimensional mounds. Imaging at the single-cell level around the defects allowed us to quantify the velocity field and the evolving cell density; cells not only concentrate at +1/2 defects, but also escape from -1/2 defects. We propose a generic mechanism for the instability in cell density around the defects that arises from the interplay between the anisotropic friction and the active force field.

Probing the shear viscosity of an active nematic film

In vitro reconstituted active systems, such as the adenosine triphosphate (ATP)-driven microtubule bundle suspension developed by the Dogic group [T. Sanchez, D. T. Chen, S. J. DeCamp, M. Heymann, and Z. Dogic, Nature (London) 491, 431 (2012)10.1038/nature11591], provide a fertile testing ground for elucidating the phenomenology of active liquid crystalline states. Controlling such novel phases of matter crucially depends on our knowledge of their material and physical properties. In this Rapid Communication, we show that the shear viscosity of an active nematic film can be probed by varying its hydrodynamic coupling to a bounding oil layer. Using the motion of disclinations as intrinsic tracers of the flow field and a hydrodynamic model, we obtain an estimate for the shear viscosity of the nematic film. Knowing this now provides us with an additional handle for robust and precision tunable control of the emergent dynamics of active fluids.

Mesoscale pattern formation of self-propelled rods with velocity reversal

We study self-propelled particles with velocity reversal interacting by uniaxial (nematic) alignment within a coarse-grained hydrodynamic theory. Combining analytical and numerical continuation techniques, we show that the physics of this active system is essentially controlled by the reversal frequency. In particular, we find that elongated, high-density, ordered patterns, called bands, emerge via subcritical bifurcations from spatially homogeneous states. Our analysis reveals further that the interaction of bands is weakly attractive and, consequently, bands fuse upon collision in analogy with nonequilibrium nucleation processes. Moreover, we demonstrate that a renormalized positive line tension can be assigned to stable bands below a critical reversal rate, beyond which they are transversally unstable. In addition, we discuss the kinetic roughening of bands as well as their nonlinear dynamics close to the threshold of transversal instability. Altogether, the reduction of the multiparticle system onto the dynamics of bands provides a unified framework to understand the emergence and stability of nonequilibrium patterns in this self-propelled particle system. In this regard, our results constitute a proof of principle in favor of the hypothesis in microbiology that velocity reversal of gliding rod-shaped bacteria regulates the transitions between various self-organized patterns observed during the bacterial life cycle.

Generalized Swift-Hohenberg models for dense active suspensions

In describing the physics of living organisms, a mathematical theory that captures the generic ordering principles of intracellular and multicellular dynamics is essential for distinguishing between universal and system-specific features. Here, we compare two recently proposed nonlinear high-order continuum models for active polar and nematic suspensions, which aim to describe collective migration in dense cell assemblies and the ordering processes in ATP-driven microtubule-kinesin networks, respectively. We discuss the phase diagrams of the two models and relate their predictions to recent experiments. The satisfactory agreement with existing experimental data lends support to the hypothesis that non-equilibrium pattern formation phenomena in a wide range of active systems can be described within the same class of higher-order partial differential equations.

Dressed active particles in spherical crystals

We investigate the dynamics of an active particle in two-dimensional spherical crystals, which provide an ideal environment to illustrate the interplay between active particles and crystallographic defects. A moving active particle is observed to be surrounded by localized topological defects, becoming a dressed active particle. Such a physical picture characterizes both the lattice distortion around the moving particle and the healing of the distorted lattice in its trajectory. We find that the dynamical behaviors of an active particle in both random and ballistic motions uniformly conform to this featured scenario, whether the particle is initially a defect or not. We further observe that the defect pattern around a dressed ballistic active particle randomly oscillates between two well-defined wing-like defect motifs regardless of its speed. The established physical picture of dressed active particles in this work partially deciphers the complexity of the intriguing nonequilibrium behaviors in active crystals, and opens the promising possibility of introducing the activity to engineer defects, which has strong connections with the design of materials.

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

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

Control of active liquid crystals with a magnetic field

Living cells sense the mechanical features of their environment and adapt to it by actively remodeling their peripheral network of filamentary proteins, known as cortical cytoskeleton. By mimicking this principle, we demonstrate an effective control strategy for a microtubule-based active nematic in contact with a hydrophobic thermotropic liquid crystal. By using well-established protocols for the orientation of liquid crystals with a uniform magnetic field, and through the mediation of anisotropic shear stresses, the active nematic reversibly self-assembles with aligned flows and textures that feature orientational order at the millimeter scale. The turbulent flow, characteristic of active nematics, is in this way regularized into a laminar flow with periodic velocity oscillations. Once patterned, the microtubule assembly reveals its intrinsic length and time scales, which we correlate with the activity of motor proteins, as predicted by existing theories of active nematics. The demonstrated commanding strategy should be compatible with other viable active biomaterials at interfaces, and we envision its use to probe the mechanics of the intracellular matrix.

Mathematical model suitable for efficient simulation of thin semi-flexible polymers in complex environments

We present an alternative approach to simulations of semi-flexible polymers. In contrast with the usual bead-rod compromise between bead-spring and rigid rod models, we use deformable cylindrical segments as basic units of the polymer. The length of each segment is not preserved with end points diffusing under constraints keeping the polymer chain nature intact. The model allows the simulation of tension transport and elasticity properties. In particular we describe a new cooperative regime in the relaxation of the polymer from its fully elongated configuration.