Defect annihilation and proliferation in active nematics

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.

Emergence of dynamic vortex glasses in disordered polar active fluids

In equilibrium, disorder conspires with topological defects to redefine the ordered states of matter in systems as diverse as crystals, superconductors, and liquid crystals. Far from equilibrium, however, the consequences of quenched disorder on active condensed matter remain virtually uncharted. Here, we reveal a state of strongly disordered active matter with no counterparts in equilibrium: a dynamical vortex glass. Combining microfluidic experiments and theory, we show how colloidal flocks collectively cruise through disordered environments without relaxing the topological singularities of their flows. The resulting state is highly dynamical but the flow patterns, shaped by a finite density of frozen vortices, are stationary and exponentially degenerated. Quenched isotropic disorder acts as a random gauge field turning active liquids into dynamical vortex glasses. We argue that this robust mechanism should shape the collective dynamics of a broad class of disordered active matter, from synthetic active nematics to collections of living cells exploring heterogeneous media.

Machine learning active-nematic hydrodynamics

Hydrodynamic theories effectively describe many-body systems out of equilibrium in terms of a few macroscopic parameters. However, such parameters are difficult to determine from microscopic information. Seldom is this challenge more apparent than in active matter, where the hydrodynamic parameters are in fact fields that encode the distribution of energy-injecting microscopic components. Here, we use active nematics to demonstrate that neural networks can map out the spatiotemporal variation of multiple hydrodynamic parameters and forecast the chaotic dynamics of these systems. We analyze biofilament/molecular-motor experiments with microtubule/kinesin and actin/myosin complexes as computer vision problems. Our algorithms can determine how activity and elastic moduli change as a function of space and time, as well as adenosine triphosphate (ATP) or motor concentration. The only input needed is the orientation of the biofilaments and not the coupled velocity field which is harder to access in experiments. We can also forecast the evolution of these chaotic many-body systems solely from image sequences of their past using a combination of autoencoders and recurrent neural networks with residual architecture. In realistic experimental setups for which the initial conditions are not perfectly known, our physics-inspired machine-learning algorithms can surpass deterministic simulations. Our study paves the way for artificial-intelligence characterization and control of coupled chaotic fields in diverse physical and biological systems, even in the absence of knowledge of the underlying dynamics.

Geometry and mechanics of disclination lines in 3D nematic liquid crystals

In 3D nematic liquid crystals, disclination lines have a range of geometric structures. Locally, they may resemble +1/2 or -1/2 defects in 2D nematic phases, or they may have 3D twist. Here, we analyze the structure in terms of the director deformation modes around the disclination, as well as the nematic order tensor inside the disclination core. Based on this analysis, we construct a vector to represent the orientation of the disclination, as well as tensors to represent higher-order structure. We apply this method to simulations of a 3D disclination arch, and determine how the structure changes along the contour length. We then use this geometric analysis to investigate three types of forces acting on a disclination: Peach-Koehler forces due to external stress, interaction forces between disclination lines, and active forces. These results apply to the motion of disclination lines in both conventional and active liquid crystals.

Integer topological defects of cell monolayers: Mechanics and flows

Monolayers of anisotropic cells exhibit long-ranged orientational order and topological defects. During the development of organisms, orientational order often influences morphogenetic events. However, the linkage between the mechanics of cell monolayers and topological defects remains largely unexplored. This holds specifically at the timescales relevant for tissue morphogenesis. Here, we build on the physics of liquid crystals to determine material parameters of cell monolayers. In particular, we use a hydrodynamical description of an active polar fluid to study the steady-state mechanical patterns at integer topological defects. Our description includes three distinct sources of activity: traction forces accounting for cell-substrate interactions as well as anisotropic and isotropic active nematic stresses accounting for cell-cell interactions. We apply our approach to C2C12 cell monolayers in small circular confinements, which form isolated aster or spiral topological defects. By analyzing the velocity and orientational order fields in spirals as well as the forces and cell number density fields in asters, we determine mechanical parameters of C2C12 cell monolayers. Our work shows how topological defects can be used to fully characterize the mechanical properties of biological active matter.

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.

Role of Friction in Multidefect Ordering

We use continuum simulations to study the impact of friction on the ordering of defects in an active nematic. Even in a frictionless system, +1/2 defects tend to align side by side and orient antiparallel reflecting their propensity to form, and circulate with, flow vortices. Increasing friction enhances the effectiveness of the defect-defect interactions, and defects form dynamically evolving, large-scale, positionally, and orientationally ordered structures, which can be explained as a competition between hexagonal packing, preferred by the -1/2 defects, and rectangular packing, preferred by the +1/2 defects.

Active microfluidic transport in two-dimensional handlebodies

Unlike traditional nematic liquid crystals, which adopt ordered equilibrium configurations compatible with the topological constraints imposed by the boundaries, active nematics are intrinsically disordered because of their self-sustained internal flows. Controlling the flow patterns of active nematics remains a limiting step towards their use as functional materials. Here we show that confining a tubulin-kinesin active nematic to a network of connected annular microfluidic channels enables controlled directional flows and autonomous transport. In single annular channels, for narrow widths, the typically chaotic streams transform into well-defined circulating flows, whose direction or handedness can be controlled by introducing asymmetric corrugations on the channel walls. The dynamics is altered when two or three annular channels are interconnected. These more complex topologies lead to scenarios of synchronization, anti-correlation, and frustration of the active flows, and to the stabilisation of high topological singularities in both the flow field and the orientational field of the material. Controlling textures and flows in these microfluidic platforms opens unexplored perspectives towards their application in biotechnology and materials science.

Flow States and Transitions of an Active Nematic in a Three-Dimensional Channel

We use active nematohydrodynamics to study the flow of an active fluid in a 3D microchannel, finding a transition between active turbulence and regimes where there is a net flow along the channel. We show that the net flow is only possible if the active nematic is flow aligning and that, in agreement with experiments, the appearance of the net flow depends on the aspect ratio of the channel cross section. We explain our results in terms of when the hydrodynamic screening due to the channel walls allows the emergence of vortex rolls across the channel.

Wrinkling Instability in 3D Active Nematics

In nature, interactions between biopolymers and motor proteins give rise to biologically essential emergent behaviors. Besides cytoskeleton mechanics, active nematics arise from such interactions. Here we present a study on 3D active nematics made of microtubules, kinesin motors, and depleting agent. It shows a rich behavior evolving from a nematically ordered space-filling distribution of microtubule bundles toward a flattened and contracted 2D ribbon that undergoes a wrinkling instability and subsequently transitions into a 3D active turbulent state. The wrinkle wavelength is independent of the ATP concentration and our theoretical model describes its relation with the appearance time. We compare the experimental results with a numerical simulation that confirms the key role of kinesin motors in cross-linking and sliding the microtubules. Our results on the active contraction of the network and the independence of wrinkle wavelength on ATP concentration are important steps forward for the understanding of these 3D systems.

Pattern Formation and Defect Ordering in Active Chiral Nematics

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

Dynamics of active nematic defects on the surface of a sphere

A nematic liquid crystal confined to the surface of a sphere exhibits topological defects of total charge +2 due to the topological constraint. In equilibrium, the nematic field forms four +1/2 defects, located at the corners of a regular tetrahedron inscribed within the sphere, since this minimizes the Frank elastic energy. If additionally the individual nematogens exhibit self-driven directional motion, the resulting active system creates large-scale flow that drives it out of equilibrium. In particular, the defects now follow complex dynamic trajectories which, depending on the strength of the active forcing, can be periodic (for weak forcing) or chaotic (for strong forcing). In this paper we derive an effective particle theory for this system, in which the topological defects are the degrees of freedom, whose exact equations of motion we subsequently determine. Numerical solutions of these equations confirm previously observed characteristics of their dynamics and clarify the role played by the time dependence of their global rotation. We also show that Onsager's variational principle offers an exceptionally transparent way to derive these dynamical equations, and we explain the defect mobility at the hydrodynamics level.

Phases and excitations of active rod-bead mixtures: simulations and experiments

We present a large-scale numerical study, supplemented by experimental observations, on a quasi-two-dimensional active system of polar rods and spherical beads confined between two horizontal plates and energised by vertical vibration. For a low rod concentration Φr, our observations are consistent with a direct phase transition, as bead concentration Φb is increased, from the isotropic phase to a homogeneous flock. For Φr above a threshold value, an ordered band dense in both rods and beads occurs between the disordered phase and the homogeneous flock, in both experiments and simulations. Within the size ranges accessible, we observe only a single band, whose width increases with Φr. Deep in the ordered state, we observe broken-symmetry "sound" modes and giant number fluctuations. The direction-dependent sound speeds and the scaling of fluctuations are consistent with the predictions of field theories of flocking; sound damping rates show departures from such theories, but the range of wavenumbers explored is modest. At very high densities, we see phase separation into rod-rich and bead-rich regions, both of which move coherently.

Pressure-driven changes to spontaneous flow in active nematic liquid crystals

We consider the effects of a pressure gradient on the spontaneous flow of an active nematic liquid crystal in a channel, subject to planar anchoring and no-slip conditions on the boundaries of the channel. We employ a model based on the Ericksen-Leslie theory of nematics, with an additional active stress accounting for the activity of the fluid. By directly solving the flow equation, we consider an asymptotic solution for the director angle equation for large activity parameter values and predict the possible values of the director angle in the bulk of the channel. Through a numerical solution of the full nonlinear equations, we examine the effects of pressure on the branches of stable and unstable equilibria, some of which are disconnected from the no-flow state. In the absence of a pressure gradient, solutions are either symmetric or antisymmetric about the channel midpoint; these symmetries are changed by the pressure gradient. Considering the activity-pressure state space allows us to predict qualitatively the extent of each solution type and to show, for large enough pressure gradients, that a branch of non-trivial director angle solutions exists for all activity values.

Steady-state distributions and nonsteady dynamics in nonequilibrium systems

We search for steady states in a class of fluctuating and driven physical systems that exhibit sustained currents. We find that the physical concept of a steady state, well known for systems at equilibrium, must be generalized to describe such systems. In these, the generalization of a steady state is associated with a stationary probability density of microstates and a deterministic dynamical system whose trajectories the system follows on average. These trajectories are a manifestation of nonstationary macroscopic currents observed in these systems. We determine precise conditions for the steady state to exist as well as the requirements for it to be stable. We illustrate this with some examples.

Topology control of human fibroblast cells monolayer by liquid crystal elastomer

Eukaryotic cells in living tissues form dynamic patterns with spatially varying orientational order that affects important physiological processes such as apoptosis and cell migration. The challenge is how to impart a predesigned map of orientational order onto a growing tissue. Here, we demonstrate an approach to produce cell monolayers of human dermal fibroblasts with predesigned orientational patterns and topological defects using a photoaligned liquid crystal elastomer (LCE) that swells anisotropically in an aqueous medium. The patterns inscribed into the LCE are replicated by the tissue monolayer and cause a strong spatial variation of cells phenotype, their surface density, and number density fluctuations. Unbinding dynamics of defect pairs intrinsic to active matter is suppressed by anisotropic surface anchoring allowing the estimation of the elastic characteristics of the tissues. The demonstrated patterned LCE approach has potential to control the collective behavior of cells in living tissues, cell differentiation, and tissue morphogenesis.

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

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

Chiral stresses in nematic cell monolayers

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

Microfluidic control over topological states in channel-confined nematic flows

Compared to isotropic liquids, orientational order of nematic liquid crystals makes their rheological properties more involved, and thus requires fine control of the flow parameters to govern the orientational patterns. In microfluidic channels with perpendicular surface alignment, nematics discontinuously transition from perpendicular structure at low flow rates to flow-aligned structure at high flow rates. Here we show how precise tuning of the driving pressure can be used to stabilize and manipulate a previously unresearched topologically protected chiral intermediate state which arises before the homeotropic to flow-aligned transition. We characterize the mechanisms underlying the transition and construct a phenomenological model to describe the critical behaviour and the phase diagram of the observed chiral flow state, and evaluate the effect of a forced symmetry breaking by introduction of a chiral dopant. Finally, we induce transitions on demand through channel geometry, application of laser tweezers, and careful control of the flow rate.

It is interesting phenomenon that chiral order can emerge in intrinsically achiral liquid crystals. Here Čopar et al. demonstrate achiral-to-chiral transition of the nematic liquid crystals flow in microfluidic channels and their behaviour, stability, and dependence on geometric and material parameters.

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.

Active matter invasion

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

Dynamic clustering and re-dispersion in concentrated colloid-active gel composites

We study the dynamics of quasi-two-dimensional concentrated suspensions of colloidal particles in active gels by computer simulations. Remarkably, we find that activity induces a dynamic clustering of colloids even in the absence of any preferential anchoring of the active nematic director at the particle surface. When such an anchoring is present, active stresses instead compete with elastic forces and re-disperse the aggregates observed in passive colloid-liquid crystal composites. Our quasi-two-dimensional "inverse" dispersions of passive particles in active fluids (as opposed to the more common "direct" suspensions of active particles in passive fluids) provide a promising route towards the self-assembly of new soft materials.

Nonsphere Drop Impact Assembly of Graphene Oxide Liquid Crystals

Creating long-lived topological textured liquid crystals (LCs) in confined nonspherical space is of significance in both generations of structures and fundamental studies of topological physics. However, it remains a great challenge due to the fluid character of LCs and the unstable tensional state of transient nonspheres. Here, we realize a rich series of topological textures confined in nonspherical geometries by drop impact assembly (DIA) of graphene oxide (GO) aqueous LCs. Various highly curved nonspherical morphologies of LCs were captured by gelator bath, generating distinct out-of-equilibrium yet long-lived macroscopic topological textures in 3D confinement. Our hydrodynamic investigations on DIA processes reveal that the shear-thinning fluid behavior of LCs and the arrested GO alignments mainly contribute to the topological richness in DIA. Utilizing the shaping behavior of GO LCs compared to other conventional linear polymers such as alginate, we further extend the DIA methodology to design more complex yet highly controllable functional composites and hybrids. This work thus reveals the potential to scale production of uniform yet anisotropic materials with rich topologic textures and tailored composition.

Lattice Boltzmann methods and active fluids

We review the state of the art of active fluids with particular attention to hydrodynamic continuous models and to the use of Lattice Boltzmann Methods (LBM) in this field. We present the thermodynamics of active fluids, in terms of liquid crystals modelling adapted to describe large-scale organization of active systems, as well as other effective phenomenological models. We discuss how LBM can be implemented to solve the hydrodynamics of active matter, starting from the case of a simple fluid, for which we explicitly recover the continuous equations by means of Chapman-Enskog expansion. Going beyond this simple case, we summarize how LBM can be used to treat complex and active fluids. We then review recent developments concerning some relevant topics in active matter that have been studied by means of LBM: spontaneous flow, self-propelled droplets, active emulsions, rheology, active turbulence, and active colloids.

Activity Driven Orientational Order in Active Nematic Liquid Crystals on an Anisotropic Substrate

We investigate the effect of an anisotropic substrate on the turbulent dynamics of a simulated two-dimensional active nematic. This is introduced as an anisotropic friction and an effective anisotropic viscosity, with the orientation of the anisotropy being defined by the substrate. In this system, we observe the emergence of global nematic order of topological defects that is controlled by the degree of anisotropy in the viscosity and the magnitude of the active stress. No global defect alignment is seen in passive liquid crystals with anisotropic viscosity or friction confirming that ordering is driven by the active stress. We then closely examine the active flow generated by a single defect to show that the net kinetic energy of the flow is dependent on the orientation of the defect relative to the substrate, resulting in a torque on the defect to align it with the anisotropy in the substrate.

Self-propulsion of an active polar drop

We investigate the self-propulsive motion of a drop containing an active polar field. The drop demonstrates spontaneous symmetry breaking from a uniform orientational order into a splay or bend instability depending on the types of active stress, namely, contractile or extensile, respectively. We develop an analytical theory of the mechanism of this instability, which has been observed only in numerical simulations. We show that both contractile and extensile active stresses result in the instability and self-propulsive motion. We also discuss asymmetry between contractile and extensile stresses and show that extensile active stress generates chaotic motion even under a simple model of the polarity field coupled with motion and deformation of the drop.

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.

Geometrical Control of Active Turbulence in Curved Topographies

We investigate the turbulent dynamics of a two-dimensional active nematic liquid crystal constrained to a curved surface. Using a combination of hydrodynamic and particle-based simulations, we demonstrate that the fundamental structural features of the fluid, such as the topological charge density, the defect number density, the nematic order parameter, and defect creation and annihilation rates, are approximately linear functions of the substrate Gaussian curvature, which then acts as a control parameter for the chaotic flow. Our theoretical predictions are then compared with experiments on microtubule-kinesin suspensions confined on toroidal droplets, finding excellent qualitative agreement.

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.

Self-organized dynamics and the transition to turbulence of confined active nematics

We study how confinement transforms the chaotic dynamics of bulk microtubule-based active nematics into regular spatiotemporal patterns. For weak confinements in disks, multiple continuously nucleating and annihilating topological defects self-organize into persistent circular flows of either handedness. Increasing confinement strength leads to the emergence of distinct dynamics, in which the slow periodic nucleation of topological defects at the boundary is superimposed onto a fast procession of a pair of defects. A defect pair migrates toward the confinement core over multiple rotation cycles, while the associated nematic director field evolves from a distinct double spiral toward a nearly circularly symmetric configuration. The collapse of the defect orbits is punctuated by another boundary-localized nucleation event, that sets up long-term doubly periodic dynamics. Comparing experimental data to a theoretical model of an active nematic reveals that theory captures the fast procession of a pair of [Formula: see text] defects, but not the slow spiral transformation nor the periodic nucleation of defect pairs. Theory also fails to predict the emergence of circular flows in the weak confinement regime. The developed confinement methods are generalized to more complex geometries, providing a robust microfluidic platform for rationally engineering 2D autonomous flows.

Active matter invasion of a viscous fluid: Unstable sheets and a no-flow theorem

We investigate the dynamics of a dilute suspension of hydrodynamically interacting motile or immotile stress-generating swimmers or particles as they invade a surrounding viscous fluid. Colonies of aligned pusher particles are shown to elongate in the direction of particle orientation and undergo a cascade of transverse concentration instabilities, governed at small times by an equation that also describes the Saffman-Taylor instability in a Hele-Shaw cell, or the Rayleigh-Taylor instability in a two-dimensional flow through a porous medium. Thin sheets of aligned pusher particles are always unstable, while sheets of aligned puller particles can either be stable (immotile particles), or unstable (motile particles) with a growth rate that is nonmonotonic in the force dipole strength. We also prove a surprising "no-flow theorem": a distribution initially isotropic in orientation loses isotropy immediately but in such a way that results in no fluid flow everywhere and for all time.

Active transport in a channel: stabilisation by flow or thermodynamics

Recent experiments on active materials, such as dense bacterial suspensions and microtubule-kinesin motor mixtures, show a promising potential for achieving self-sustained flows. However, to develop active microfluidics it is necessary to understand the behaviour of active systems confined to channels. Therefore here we use continuum simulations to investigate the behaviour of active fluids in a two-dimensional channel. Motivated by the fact that most experimental systems show no ordering in the absence of activity, we concentrate on temperatures where there is no nematic order in the passive system, so that any nematic order is induced by the active flow. We systematically analyze the results, identify several different stable flow states, provide a phase diagram and show that the key parameters controlling the flow are the ratio of channel width to the length scale of active flow vortices, and whether the system is flow aligning or flow tumbling.

Optimizing active work: Dynamical phase transitions, collective motion, and jamming

Active work measures how far the local self-forcing of active particles translates into real motion. Using population Monte Carlo methods, we investigate large deviations in the active work for repulsive active Brownian disks. Minimizing the active work generically results in dynamical arrest; in contrast, despite the lack of aligning interactions, trajectories of high active work correspond to a collectively moving, aligned state. We use heuristic and analytic arguments to explain the origin of dynamical phase transitions separating the arrested, typical, and aligned regimes.

Emergence of Active Nematic Behavior in Monolayers of Isotropic Cells

There is now growing evidence of the emergence and biological functionality of liquid crystal features, including nematic order and topological defects, in cellular tissues. However, how such features that intrinsically rely on particle elongation emerge in monolayers of cells with isotropic shapes is an outstanding question. In this Letter, we present a minimal model of cellular monolayers based on cell deformation and force transmission at the cell-cell interface that explains the formation of topological defects and captures the flow-field and stress patterns around them. By including mechanical properties at the individual cell level, we further show that the instability that drives the formation of topological defects, and leads to active turbulence, emerges from a feedback between shape deformation and active driving. The model allows us to suggest new explanations for experimental observations in tissue mechanics, and to propose designs for future experiments.

Theory of defect motion in 2D passive and active nematic liquid crystals

The motion of topological defects is an important feature of the dynamics of all liquid crystals, and is especially conspicuous in active liquid crystals. Understanding defect motion is a challenging theoretical problem, because the dynamics of orientational order is coupled with backflow of the fluid, and because a liquid crystal has several distinct viscosity coefficients. Here, we suggest a coarse-grained, variational approach, which describes the motion of defects as effective "particles". For passive liquid crystals, the theory shows how the drag depends on defect orientation, and shows the coupling between translational and rotational motion. For active liquid crystals, the theory provides an alternative way to describe motion induced by the activity coefficient.

Dynamical subclasses of dry active nematics

We show that the dominant mode of alignment plays an important role in dry active nematics, leading to two dynamical subclasses defined by the nature of the instability of the nematic bands that characterize, in these systems, the coexistence phase separating the isotropic and fluctuating nematic states. In addition to the well-known instability inducing long undulations along the band, another stronger instability leading to the breakup of the band in many transversal segments may arise. We elucidate the origin of this strong instability for a realistic model of self-propelled rods and determine the high-order nonlinear terms responsible for it at the hydrodynamic level.

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.

Harnessing Medium Anisotropy To Control Active Matter

Active matter is a wide class of nonequilibrium systems consisting of interacting self-propelled agents transducing the energy stored in the environment into mechanical motion. Numerous examples range from microscopic cytoskeletal filaments and swimming organisms (bacteria and unicellular algae), synthetic catalytic nanomotors, colloidal self-propelled Janus particles, to macroscopic bird flocks, fish schools, and even human crowds. Active matter demonstrates a remarkable tendency toward self-organization and development of collective states with the long-range spatial order. Furthermore, active materials exhibit properties that are not present in traditional materials like plastics or ceramics: self-repair, shape change, and adaptation. A suspension of microscopic swimmers, such as motile bacteria or self-propelled colloids (active suspensions), is possibly the simplest and the most explored realization of active matter. Recent studies of active suspensions revealed a wealth of unexpected behaviors, from a dramatic reduction of the effective viscosity, enhanced mixing and self-diffusion, rectification of chaotic motion, to artificial rheotaxis (drift against the imposed flow) and cross-stream migration. To date, most of the studies of active matter are performed in isotropic suspending medium, like water with the addition of some "fuel", e.g., nutrient for bacteria or H₂O₂ for catalytic bimetallic AuPt nanorods. A highly structured anisotropic suspending medium represented by lyotropic liquid crystal (water-soluble) opens enormous opportunities to control and manipulate active matter. Liquid crystals exhibit properties intermediate between solid and liquids; they may flow like a liquid but respond to deformations as a solid due to a crystal-like orientation of molecules. Liquid crystals doped by a small amount of active component represent a new class of composite materials (living liquid crystals or LLCs) with unusual mechanical and optical properties. LLCs demonstrate a variety of highly organized dynamic collective states, spontaneous formation of dynamic textures of topological defects (singularities of local molecular orientation), controlled and reconfigurable transport of cargo particles, manipulation of individual trajectories of microswimmers, and many others. Besides insights into fundamental mechanisms governing active materials, living liquid crystals may have intriguing applications, such as the design of new classes of soft adaptive bioinspired materials capable to respond to physical and chemical stimuli, such as light, magnetic, and electric fields, mechanical shear, airborne pollutants, and bacterial toxins. This Account details the most recent developments in the field of LLCs and discusses how the anisotropy of liquid crystals can be harnessed to control and manipulate active materials.

A growing bacterial colony in two dimensions as an active nematic

How a single bacterium becomes a colony of many thousand cells is important in biomedicine and food safety. Much is known about the molecular and genetic bases of this process, but less about the underlying physical mechanisms. Here we study the growth of single-layer micro-colonies of rod-shaped Escherichiacoli bacteria confined to just under the surface of soft agarose by a glass slide. Analysing this system as a liquid crystal, we find that growth-induced activity fragments the colony into microdomains of well-defined size, whilst the associated flow orients it tangentially at the boundary. Topological defect pairs with charges \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pm {\textstyle{1 \over 2}}$$\end{document}±12 are produced at a constant rate, with the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$+ {\textstyle{1 \over 2}}$$\end{document}+12 defects being propelled to the periphery. Theoretical modelling suggests that these phenomena have different physical origins from similar observations in other extensile active nematics, and a growing bacterial colony belongs to a new universality class, with features reminiscent of the expanding universe.

Rod-shaped bacteria are an example of active matter. Here the authors find that a growing bacterial colony harbours internal cellular flows affecting orientational ordering in its interior and at the boundary. Results suggest this system may belong to a new active matter universality class.

Tunable structure and dynamics of active liquid crystals

Active materials are capable of converting free energy into directional motion, giving rise to notable dynamical phenomena. Developing a general understanding of their structure in relation to the underlying nonequilibrium physics would provide a route toward control of their dynamic behavior and pave the way for potential applications. The active system considered here consists of a quasi-two-dimensional sheet of short (≈1 μm) actin filaments driven by myosin II motors. By adopting a concerted theoretical and experimental strategy, new insights are gained into the nonequilibrium properties of active nematics over a wide range of internal activity levels. In particular, it is shown that topological defect interactions can be led to transition from attractive to repulsive as a function of initial defect separation and relative orientation. Furthermore, by examining the +1/2 defect morphology as a function of activity, we found that the apparent elastic properties of the system (the ratio of bend-to-splay elastic moduli) are altered considerably by increased activity, leading to an effectively lower bend elasticity. At high levels of activity, the topological defects that decorate the material exhibit a liquid-like structure and adopt preferred orientations depending on their topological charge. Together, these results suggest that it should be possible to tune internal stresses in active nematic systems with the goal of designing out-of-equilibrium structures with engineered dynamic responses.

Edges impose planar alignment in nematic monolayers by directing cell elongation and enhancing migration

Boundaries play an important role in the emergence of nematic order in classical liquid crystal systems; we explore their importance in adhesive cells that form active nematics. In particular, we study how cells are affected by an edge, which in our experiments is a boundary between adhesive and non-adhesive domains on a planar surface. We find that such edges induce elongation and direct the migration of isolated fibroblasts. In confluent monolayers, these elongated cells co-align and migrate to form an active, two-dimensional nematic structure in which edges enforce planar alignment and provide local slip to streams of cells that move along them. On an adhesive square island of dimensions 1 mm × 1 mm, cells near the edges in confluent nematic monolayers have enhanced alignment and velocity. The corners of the adhesive island seed defects with signs that depend on the direction of the motion of the streams of cells that meet there. Distortions emerge with rotations of -π/2 to form a -1/4 defect for streams that move clockwise or counterclockwise, and +π/2 to form a +1/4 defect for converging streams. We explore how cells transmit alignment information to each other in the absence of an edge by studying cell pairs and find that while such pairs do co-align, this alignment is only transient and short lived. These results shed light on the importance of edges in imposing nematic order in confluent monolayers and how edges can be used as tools to pattern the long-range organization of cells for tissue engineering applications.

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.

Frustration-induced phases in migrating cell clusters

Certain malignant cancer cells form clusters in a chemoattractant gradient, which can spontaneously show three different phases of motion: translational, rotational, and random. Guided by our experiments on the motion of two-dimensional clusters in vitro, we developed an agent-based model in which the cells form a cohesive cluster due to attractive and alignment interactions. We find that when cells at the cluster rim are more motile, all three phases of motion coexist, in agreement with our observations. Using the model, we show that the transitions between different phases are driven by competition between an ordered rim and a disordered core accompanied by the creation and annihilation of topological defects in the velocity field. The model makes specific predictions, which we verify with our experimental data. Our results suggest that heterogeneous behavior of individuals, based on local environment, can lead to novel, experimentally observed phases of collective motion.

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.

Twist-induced crossover from two-dimensional to three-dimensional turbulence in active nematics

While studies of active nematics in two dimensions have shed light on various aspects of the flow regimes and topology of active matter, three-dimensional properties of topological defects and chaotic flows remain unexplored. By confining a film of active nematics between two parallel plates, we use continuum simulations and analytical arguments to demonstrate that the crossover from quasi-two-dimensional (quasi-2D) to three-dimensional (3D) chaotic flows is controlled by the morphology of the disclination lines. For small plate separations, the active nematic behaves as a quasi-2D material, with straight topological disclination lines spanning the height of the channel and exhibiting effectively 2D active turbulence. Upon increasing channel height, we find a crossover to 3D chaotic flows due to the contortion of disclinations above a critical activity. Above this critical activity highly contorted disclination lines and disclination loops are formed. We further show that these contortions are engendered by twist perturbations producing a sharp change in the curvature of disclinations.

A nonequilibrium force can stabilize 2D active nematics

Suspensions of actively driven anisotropic objects exhibit distinctively nonequilibrium behaviors, and current theories predict that they are incapable of sustaining orientational order at high activity. By contrast, here we show that nematic suspensions on a substrate can display order at arbitrarily high activity due to a previously unreported, potentially stabilizing active force. This force moreover emerges inevitably in theories of active orientable fluids under geometric confinement. The resulting nonequilibrium ordered phase displays robust giant number fluctuations that cannot be suppressed even by an incompressible solvent. Our results apply to virtually all experimental assays used to investigate the active nematic ordering of self-propelled colloids, bacterial suspensions, and the cytoskeleton and have testable implications in interpreting their nonequilibrium behaviors.

Spontaneously Flowing Crystal of Self-Propelled Particles

We experimentally and numerically study the structure and dynamics of a monodisperse packing of spontaneously aligning self-propelled hard disks. The packings are such that their equilibrium counterparts form perfectly ordered hexagonal structures. Experimentally, we first form a perfect crystal in a hexagonal arena which respects the same crystalline symmetry. Frustration of the hexagonal order, obtained by removing a few particles, leads to the formation of a rapidly diffusing "droplet." Removing more particles, the whole system spontaneously forms a macroscopic sheared flow, while conserving an overall crystalline structure. This flowing crystalline structure, which we call a "rheocrystal," is made possible by the condensation of shear along localized stacking faults. Numerical simulations very well reproduce the experimental observations and allow us to explore the parameter space. They demonstrate that the rheocrystal is induced neither by frustration nor by noise. They further show that larger systems flow faster while still remaining ordered.

Turbulent Dynamics of Epithelial Cell Cultures

We investigate the large length and long time scales collective flows and structural rearrangements within in vitro human bronchial epithelial cell (HBEC) cultures. Activity-driven collective flows result in ensembles of vortices randomly positioned in space. By analyzing a large population of vortices, we show that their area follows an exponential law with a constant mean value and their rotational frequency is size independent, both being characteristic features of the chaotic dynamics of active nematic suspensions. Indeed, we find that HBECs self-organize in nematic domains of several cell lengths. Nematic defects are found at the interface between domains with a total number that remains constant due to the dynamical balance of nucleation and annihilation events. The mean velocity fields in the vicinity of defects are well described by a hydrodynamic theory of extensile active nematics.