Velocity correlations in an active nematic

Evidence of universal conformal invariance in living biological matter

The emergent dynamics of collective cellular movement are typically thought to depend on how cells interact with one another and the mechanisms used to drive motility, both of which exhibit remarkable diversity across different biological systems. Here we report experimental evidence of a universal feature in the patterns of flow that spontaneously emerge in groups of collectively moving cells. Specifically, we demonstrate that the flows generated by collectively moving dog kidney cells, human breast cancer cells and two different strains of pathogenic bacteria exhibit robust conformal invariance. We also show that the precise form of invariance in all four systems is described by the Schramm-Loewner evolution-a family of planar curves defined by a single parameter-and belongs to the percolation universality class. The presence of universal conformal invariance reveals that the macroscopic features of living biological matter exhibit universal translational, rotational and scale symmetries that are independent of the microscopic properties of its constituents. Our results show that flow patterns generated by different systems are highly conserved and that biological systems can be used to experimentally test predictions from the theories for conformally invariant structures.

Inertia Drives Concentration-Wave Turbulence in Swimmer Suspensions

We discover an instability mechanism in suspensions of self-propelled particles that does not involve active stress. Instead, it is driven by a subtle interplay of inertia, swimmer motility, and concentration fluctuations, through a crucial time lag between the velocity and the concentration field. The resulting time-persistent state seen in our high-resolution numerical simulations consists of self-sustained waves of concentration and orientation, transiting from regular oscillations to wave turbulence. We analyze the statistical features of this active turbulence, including an intriguing connection to the Batchelor spectrum of passive scalars.

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

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

Surface Topography Induces and Orients Nematic Swarms of Active Filaments: Considerations for Lab-On-A-Chip Devices

Surface-bound molecular motors can drive the collective motion of cytoskeletal filaments in the form of nematic bands and polar flocks in reconstituted gliding assays. Although these "swarming transitions" are an emergent property of active filament collisions, they can be controlled and guided by tuning the surface chemistry or topography of the substrate. To date, the impact of surface topography on collective motion in active nematics is only partially understood, with most experimental studies focusing on the escape of a single filament from etched channels. Since the late 1990s, significant progress has been made to utilize the nonequilibrium properties of active filaments and create a range of functional nanodevices relevant to biosensing and parallel computation; however, the complexity of these swarming transitions presents a challenge when attempting to increase filament surface concentrations. In this work, we etch shallow, linear trenches into glass substrates to induce the formation of swarming nematic bands and investigate the mechanisms by which surface topography regulates the two-dimensional (2D) collective motion of driven filamentous actin (F-actin). We demonstrate that nematic swarms only appear at intermediate trench spacings and vanish if the trenches are made too narrow, wide, or tortuous. To rationalize these results, we simulate the F-actin as self-propelled, semiflexible chains subject to a soft, spatially modulated potential that encodes the energetic cost of bending a filament along the edge of a trench. In our model, we hypothesize that an individual filament experiences a penalty when its projected end-to-end distance is smaller than the trench spacing ("bending and turning"). However, chains that span the channel width glide above the trenches in a force- and torque-free manner ("crowd-surfing"). Our simulations demonstrate that collections of filaments form nematic bands only at intermediate trench spacings, consistent with our experimental findings.

Multifractal fluctuations in zebrafish (Danio rerio) polarization time series

In this work, we study the polarization time series obtained from experimental observation of a group of zebrafish (Danio rerio) confined in a circular tank. The complex dynamics of the individual trajectory evolution lead to the appearance of multiple characteristic scales. Employing the Multifractal Detrended Fluctuation Analysis (MF-DFA), we found distinct behaviors according to the parameters used. The polarization time series are multifractal at low fish densities and their average scales with ρ - 1 / 4 . On the other hand, they tend to be monofractal, and their average scales with ρ - 1 / 2 for high fish densities. These two regimes overlap at critical density ρ c , suggesting the existence of a phase transition separating them. We also observed that for low densities, the polarization velocity shows a non-Gaussian behavior with heavy tails associated with long-range correlation and becomes Gaussian for high densities, presenting an uncorrelated regime.

Active Darcy's Law

While bacterial swarms can exhibit active turbulence in vacant spaces, they naturally inhabit crowded environments. We numerically show that driving disorderly active fluids through porous media enhances Darcy's law. While purely active flows average to zero flux, hybrid active/driven flows display greater drift than purely pressure-driven flows. This enhancement is nonmonotonic with activity, leading to an optimal activity to maximize flow rate. We incorporate the active contribution into an active Darcy's law, which may serve to help understand anomalous transport of swarming in porous media.

Novel turbulence and coarsening arrest in active-scalar fluids

We uncover a new type of turbulence - activity-induced homogeneous and isotropic turbulence - in a model that has been employed to investigate motility-induced phase separation (MIPS) in a system of microswimmers. The active Cahn-Hilliard-Navier-Stokes (CHNS) equations, also called active model H, provide a natural theoretical framework for our study. In this CHNS model, a single scalar order parameter ϕ, positive (negative) in regions of high (low) microswimmer density, is coupled with the velocity field u. The activity of the microswimmers is governed by an activity parameter ζ that is positive for extensile swimmers and negative for contractile swimmers. With extensile swimmers, this system undergoes complete phase separation, which is similar to that in binary-fluid mixtures. By carrying out pseudospectral direct numerical simulations (DNSs), we show, for the first time, that (a) this model develops an emergent nonequilibrium, but statistically steady, state (NESS) of active turbulence, for the case of contractile swimmers, if ζ is sufficiently large and negative, and (b) this turbulence arrests the phase separation. We quantify this suppression by showing how the coarsening-arrest length scale does not grow indefinitely, with time t, but saturates at a finite value at large times. We characterise the statistical properties of this active-scalar turbulence by employing energy spectra and fluxes and the spectrum of ϕ. For sufficiently high Reynolds numbers, the energy spectrum (k) displays an inertial range, with a power-law dependence on the wavenumber k. We demonstrate that, in this range, the flux Π(k) assumes a nearly constant, negative value, which indicates that the system shows an inverse cascade of energy, even though energy injection occurs over a wide range of wavenumbers in our active-CHNS model.

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

Objective

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

Methods

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

Results

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

Conclusion

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

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.

Maximally mixing active nematics

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

Active nematic liquid crystals simulated by particle-based mesoscopic methods

Two Multi-particle collision dynamics algorithms that simulate nematic liquid crystals are generalised to reproduce active behaviour. One of the algorithms is due to Shendruk and Yeomans and is based on particles that carry an orientation vector ordered by a mean-field energy [T. N. Shendruk and J. M. Yeomans, Soft Matter, 2015, 11, 5101]. In the other algorithm, due to Mandal and Mazza, particles possess an order parameter tensor which evolves according to the Qian-Sheng model of nematohydrodynamics [S. Mandal and M. G. Mazza, Phys. Rev. E, 2019, 99, 063319]. For both methods activity is incorporated through a force proportional to the divergence of the local average order parameter tensor. Both implementations produce disclination curves in the nematic fluid that undergo nucleation and self-annihilation dynamics. Topological defects are found to be consistent with those observed in recent experiments of three-dimensional active nematics. Results permit to compare the length-scales over which the different nematic Multi-particle collision dynamics methods operate. The structure and dynamics of the orientation and flow fields agree with those obtained recently in numerical studies of continuum three-dimensional active nematics. Overall, our results open the opportunity to use mesoscopic particle-based approaches to study active liquid crystals in situations such as nonequilibrium states driven by flow or colloidal particles in active anisotropic solvents.

Controlling liquid-liquid phase behaviour with an active fluid

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

Physically informed data-driven modeling of active nematics

A continuum description is essential for understanding a variety of collective phenomena in active matter. However, building quantitative continuum models of active matter from first principles can be extremely challenging due to both the gaps in our knowledge and the complicated structure of nonlinear interactions. Here, we use a physically informed data-driven approach to construct a complete mathematical model of an active nematic from experimental data describing kinesin-driven microtubule bundles confined to an oil-water interface. We find that the structure of the model is similar to the Leslie-Ericksen and Beris-Edwards models, but there are appreciable and important differences. Rather unexpectedly, elastic effects are found to play no role in the experiments considered, with the dynamics controlled entirely by the balance between active stresses and friction stresses.

Microdomains and stress distributions in bacterial monolayers on curved interfaces

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

Active turbulence and spontaneous phase separation in inhomogeneous extensile active gels

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

Data-Driven Discovery of Active Nematic Hydrodynamics

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

Defect dynamics in active polar fluids vs. active nematics

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

Bacterial active matter

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

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

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

A topological fluctuation theorem

Fluctuation theorems specify the non-zero probability to observe negative entropy production, contrary to a naive expectation from the second law of thermodynamics. For closed particle trajectories in a fluid, Stokes theorem can be used to give a geometric characterization of the entropy production. Building on this picture, we formulate a topological fluctuation theorem that depends only by the winding number around each vortex core and is insensitive to other aspects of the force. The probability is robust to local deformations of the particle trajectory, reminiscent of topologically protected modes in various classical and quantum systems. We demonstrate that entropy production is quantized in these strongly fluctuating systems, and it is controlled by a topological invariant. We demonstrate that the theorem holds even when the probability distributions are non-Gaussian functions of the generated heat.

Activity-Induced Enhancement of Superdiffusive Transport in Bacterial Turbulence

Superdiffusion processes significantly promote the transport of tiny passive particles within biological fluids. Activity, one of the essential measures for living matter, however, is less examined in terms of how and to what extent it can improve the diffusivity of the moving particles. Here, bacterial suspensions are confined within the microfluidic channel at the state of bacterial turbulence, and are tuned to different activity levels by oxygen consumption in control. Systematic measurements are conducted to determine the superdiffusion exponent, which characterizes the diffusivity strength of tracer particles, depending on the continuously injecting energy converted to motile activity from swimming individuals. Higher activity is quantified to drastically enhance the superdiffusion process of passive tracers in the short-time regime. Moreover, the number density of the swimming bacteria is controlled to contribute to the field activity, and then to strengthen the super-diffusivity of tracers, distinguished by regimes with and without collective motion of interacting bacteria. Finally, the non-slip surfaces of the microfluidic channel lower the superdiffusion of immersed tracers due to the resistance, with the small diffusivity differing from the counterpart in the bulk. The findings here suggest ways of controlled diffusion and transport of substances within the living system with different levels of nutrition and resources and boundary walls, leading to efficient mixing, drug delivery and intracellular communications.

Spontaneous flow created by active topological defects

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

Logic operations with active topological defects

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

Role of Advective Inertia in Active Nematic Turbulence

Suspensions of active agents with nematic interactions exhibit complex spatiotemporal dynamics such as mesoscale turbulence. Since the Reynolds number of microscopic flows is very small on the scale of individual agents, inertial effects are typically excluded in continuum theories of active nematic turbulence. Whether active stresses can collectively excite inertial flows is currently unclear. To address this question, we investigate a two-dimensional continuum theory for active nematic turbulence. In particular, we compare mesoscale turbulence with and without the effects of advective inertia. We find that inertial effects can influence the flow already close to the onset of the turbulent state and, moreover, give rise to large-scale fluid motion for strong active driving. A detailed analysis of the kinetic energy budget reveals an energy transfer to large scales mediated by inertial advection. While this transfer is small in comparison to energy injection and dissipation, its effects accumulate over time. The inclusion of friction, which is typically present in experiments, can compensate for this effect. The findings suggest that the inclusion of inertia and friction may be necessary for dynamically consistent theories of active nematic turbulence.

Flow transitions and length scales of a channel-confined active nematic

We perform lattice Boltzmann simulations of an active nematic fluid confined in a two-dimensional channel to study the range of flow states that are stabilised by the confinement: unidirectional flow, oscillatory flow, the dancing state, localised active turbulence and fully-developed active turbulence. We analyse the flows in Fourier space, and measure a range of different length scales which describe the flows. We argue that the different states occur as a result of flow instabilities inherent to the system. As a consequence the characteristic length scale for oscillatory flow, the dancing state and localised active turbulence is set by the channel width. Fully-developed active turbulence occurs only when the channel width is larger than the intrinsic, active length scale of the bulk fluid. The results clarify why the activity number is a control parameter for the flow transitions.

Behavior of chiral active nematics confined to nanoscopic circular region

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

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

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

Collective excitations in active fluids: Microflows and breakdown in spectral equipartition of kinetic energy

The effect of particle activity on collective excitations in active fluids of microflyers is studied. With an in silico study, we observed an oscillating breakdown of equipartition (uniform spectral distribution) of kinetic energy in reciprocal space. The phenomenon is related to short-range velocity-velocity correlations that were realized without forming of long-lived mesoscale vortices in the system. This stands in contrast to well-known mesoscale turbulence operating in active nematic systems (bacterial or artificial) and reveals the features of collective dynamics in active fluids, which should be important for structural transitions and glassy dynamics in active matter.

Spatiotemporal control of liquid crystal structure and dynamics through activity patterning

Active materials are capable of converting free energy into mechanical work to produce autonomous motion, and exhibit striking collective dynamics that biology relies on for essential functions. Controlling those dynamics and transport in synthetic systems has been particularly challenging. Here, we introduce the concept of spatially structured activity as a means of controlling and manipulating transport in active nematic liquid crystals consisting of actin filaments and light-sensitive myosin motors. Simulations and experiments are used to demonstrate that topological defects can be generated at will and then constrained to move along specified trajectories by inducing local stresses in an otherwise passive material. These results provide a foundation for the design of autonomous and reconfigurable microfluidic systems where transport is controlled by modulating activity with light.

Design of nematic liquid crystals to control microscale dynamics

The dynamics of small particles, both living such as swimming bacteria and inanimate, such as colloidal spheres, has fascinated scientists for centuries. If one could learn how to control and streamline their chaotic motion, that would open technological opportunities in the transformation of stored or environmental energy into systematic motion, with applications in micro-robotics, transport of matter, guided morphogenesis. This review presents an approach to command microscale dynamics by replacing an isotropic medium with a liquid crystal. Orientational order and associated properties, such as elasticity, surface anchoring, and bulk anisotropy, enable new dynamic effects, ranging from the appearance and propagation of particle-like solitary waves to self-locomotion of an active droplet. By using photoalignment, the liquid crystal can be patterned into predesigned structures. In the presence of the electric field, these patterns enable the transport of solid and fluid particles through nonlinear electrokinetics rooted in anisotropy of conductivity and permittivity. Director patterns command the dynamics of swimming bacteria, guiding their trajectories, polarity of swimming, and distribution in space. This guidance is of a higher level of complexity than a simple following of the director by rod-like microorganisms. Namely, the director gradients mediate hydrodynamic interactions of bacteria to produce an active force and collective polar modes of swimming. The patterned director could also be engraved in a liquid crystal elastomer. When an elastomer coating is activated by heat or light, these patterns produce a deterministic surface topography. The director gradients define an activation force that shapes the elastomer in a manner similar to the active stresses triggering flows in active nematics. The patterned elastomer substrates could be used to define the orientation of cells in living tissues. The liquid-crystal guidance holds a major promise in achieving the goal of commanding microscale active flows.

Fluctuations can induce local nematic order and extensile stress in monolayers of motile cells

Recent experiments in various cell types have shown that two-dimensional tissues often display local nematic order, with evidence of extensile stresses manifest in the dynamics of topological defects. Using a mesoscopic model where tissue flow is generated by fluctuating traction forces coupled to the nematic order parameter, we show that the resulting tissue dynamics can spontaneously produce local nematic order and an extensile internal stress. A key element of the model is the assumption that in the presence of local nematic alignment, cells preferentially crawl along the nematic axis, resulting in anisotropy of fluctuations. Our work shows that activity can drive either extensile or contractile stresses in tissue, depending on the relative strength of the contractility of the cortical cytoskeleton and tractions by cells on the extracellular matrix.

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.

Confinement-induced self-organization in growing bacterial colonies

We investigate the emergence of global alignment in colonies of dividing rod-shaped cells under confinement. Using molecular dynamics simulations and continuous modeling, we demonstrate that geometrical anisotropies in the confining environment give rise to an imbalance in the normal stresses, which, in turn, drives a collective rearrangement of the cells. This behavior crucially relies on the colony's solid-like mechanical response at short time scales and can be recovered within the framework of active hydrodynamics upon modeling bacterial colonies as growing viscoelastic gels characterized by Maxwell-like stress relaxation.

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.

Driven topological transitions in active nematic films

The topological properties of many materials are central to their behavior. In intrinsically out-of-equilibrium active materials, the dynamics of topological defects can be particularly important. In this paper, local manipulation of the order, dynamics, and topological properties of microtubule-based active nematic films is demonstrated in a joint experimental and simulation study. Hydrodynamic stresses created by magnetically actuated rotation of disk-shaped colloids in proximity to the films compete with internal stresses in the active nematic, influencing the local motion of +1/2 charge topological defects that are intrinsic to the nematic order in the spontaneously turbulent active films. Sufficiently large applied stresses drive the formation of +1 charge topological vortices through the merger of two +1/2 defects. The directed motion of the defects is accompanied by ordering of the vorticity and velocity of the active flows within the film that is qualitatively unlike the response of passive viscous films. Many features of the film's response to the stress are captured by lattice Boltzmann simulations, providing insight into the anomalous viscoelastic nature of the active nematic. The topological vortex formation is accompanied by a rheological instability in the film that leads to significant increase in the flow velocities. Comparison of the velocity profile in vicinity of the vortex with fluid-dynamics calculations provides an estimate of the film viscosity.

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.

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.

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.

Darcy's Law without Friction in Active Nematic Rheology

We study the dynamics of a contractile active nematic fluid subjected to a Poiseuille flow. In a quasi-1D geometry, we find that the linear rheology of this material is reminiscent of Darcy's law in complex fluids, with a pluglike flow decaying to zero over a well-defined "permeation" length. As a result, the viscosity increases with size, but never diverges, thereby evading the yield stress predicted by previous theories. We find strong shear thinning controlled by an active Ericksen number quantifying the ratio between external pressure difference and internal active stresses. In 2D, the increase of linear regime viscosity with size only persists up to a critical length beyond which we observe active turbulent patterns, with very low apparent viscosity. The ratio between the critical and permeation length determining the stability of the Darcy regime can be made indefinitely large by varying the flow aligning parameter or magnitude of nematic order.

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.

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.

Self-straining of actively crosslinked microtubule networks

Cytoskeletal networks are foundational examples of active matter and central to self-organized structures in the cell. In vivo, these networks are active and densely crosslinked. Relating their large-scale dynamics to the properties of their constituents remains an unsolved problem. Here, we study an in vitro active gel made from aligned microtubules and XCTK2 kinesin motors. Using photobleaching, we demonstrate that the gel's aligned microtubules, driven by motors, continually slide past each other at a speed independent of the local microtubule polarity and motor concentration. This phenomenon is also observed, and remains unexplained, in spindles. We derive a general framework for coarse graining microtubule gels crosslinked by molecular motors from microscopic considerations. Using microtubule-microtubule coupling through a force-velocity relationship for kinesin, this theory naturally explains the experimental results: motors generate an active strain rate in regions of changing polarity, which allows microtubules of opposite polarities to slide past each other without stressing the material.

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.

Particle-resolved lattice Boltzmann simulations of 3-dimensional active turbulence

Collective behaviour in suspensions of microswimmers is often dominated by the impact of long-ranged hydrodynamic interactions. These phenomena include active turbulence, where suspensions of pusher bacteria at sufficient densities exhibit large-scale, chaotic flows. To study this collective phenomenon, we use large-scale (up to N = 3 × 106) particle-resolved lattice Boltzmann simulations of model microswimmers described by extended stresslets. Such system sizes enable us to obtain quantitative information about both the transition to active turbulence and characteristic features of the turbulent state itself. In the dilute limit, we test analytical predictions for a number of static and dynamic properties against our simulation results. For higher swimmer densities, where swimmer-swimmer interactions become significant, we numerically show that the length- and timescales of the turbulent flows increase steeply near the predicted finite-system transition density.

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.