Taming active turbulence with patterned soft interfaces

Nematic liquid crystal flow driven by time-varying active surface anchoring

We demonstrate the generation of diverse material flow regimes in nematic liquid cells as driven by time-variable active surface anchoring, including no-net flow, oscillatory flow, steady flow, and pulsating flow. Specifically, we numerically simulate a passive nematic fluid inside a cell bounded with two flat solid boundaries at which the time-dependent anchoring is applied with the dynamically variable surface anchoring easy axis. We show that different flow regimes emerge as the result of different anchoring driving directions (i.e. co-rotating or counter-rotating) and relative phase of anchoring driving. The flow magnitude is tunable by cell thickness and anchoring driving frequency. More generally, this work aims towards possible applications of responsive time-variable surfaces, including photonics or synthetic active matter.

Swarm Autonomy: From Agent Functionalization to Machine Intelligence

Swarm behaviors are common in nature, where individual organisms collaborate via perception, communication, and adaptation. Emulating these dynamics, large groups of active agents can self-organize through localized interactions, giving rise to complex swarm behaviors, which exhibit potential for applications across various domains. This review presents a comprehensive summary and perspective of synthetic swarms, to bridge the gap between the microscale individual agents and potential applications of synthetic swarms. It is begun by examining active agents, the fundamental units of synthetic swarms, to understand the origins of their motility and functionality in the presence of external stimuli. Then inter-agent communications and agent-environment communications that contribute to the swarm generation are summarized. Furthermore, the swarm behaviors reported to date and the emergence of machine intelligence within these behaviors are reviewed. Eventually, the applications enabled by distinct synthetic swarms are summarized. By discussing the emergent machine intelligence in swarm behaviors, insights are offered into the design and deployment of autonomous synthetic swarms for real-world applications.

Active nematic coherence probed under spatial patterns of distributed activity

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

Mechanochemical topological defects in an active nematic

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

Coexistence of Defect Morphologies in Three-Dimensional Active Nematics

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

Topological defects in multi-layered swarming bacteria

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

Oscillating edge current in polar active fluid

Dense bacterial suspensions exhibit turbulent behavior called bacterial turbulence. The behavior of the bulk unconstrained bacterial turbulence is described well by the Toner-Tu-Swift-Hohenberg (TTSH) equation for the velocity field. However, it remains unclear how we should treat boundary conditions on bacterial turbulence in contact with some boundaries (e.g., solid walls). To be more specific, although the importance of the edge current, the flow along the boundary, has been demonstrated in several experimental studies on confined bacterial suspensions, previous numerical studies based on the TTSH equation employ nonslip boundary conditions and do not seem to properly describe the behavior of bacteria near the boundaries. In this paper, we impose a slip boundary condition on the TTSH equation to describe the bacterial motion at boundaries. We develop a method to implement the slip boundary condition. Using this method, we have successfully produced edge current and discovered that the direction of the edge current temporally oscillates. The oscillation can be attributable to the advection term in the TTSH equation. Our paper demonstrates that boundary conditions could play an important role in the collective dynamics of active systems.

Probing active nematics with in situ microfabricated elastic inclusions

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

Maximally mixing active nematics

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

Vorticity phase separation and defect lattices in the isotropic phase of active liquid crystals

We use numerical simulations and linear stability analysis to study the dynamics of an active liquid crystal film on a substrate in the regime where the passive system would be isotropic. Extensile activity builds up local orientational order and destabilizes the quiescent isotropic state above a critical activity, eventually resulting in spatiotemporal chaotic dynamics akin to the one observed ubiquitously in the nematic state. Here we show that tuning substrate friction yields a variety of emergent structures at intermediate activity, including lattices of flow vortices with associated regular arrangements of topological defects and a new state where flow vortices trap pairs of +1/2 defect that chase each other's tail. These chiral units spontaneously pick the sense of rotation and organize in a hexagonal lattice, surrounded by a diffuse flow of opposite rotation to maintain zero net vorticity. The length scale of these emergent structures is set by the screening length of the flow, controlled by the shear viscosity η and the substrate friction Γ, and can be captured by simple mode selection of the vortical flows. We demonstrate that the emergence of coherent structures can be interpreted as a phase separation of vorticity, where friction plays a role akin to that of birth/death processes in breaking conservation of the phase separating species and selecting a characteristic scale for the patterns. Our work shows that friction provides an experimentally accessible tuning parameter for designing controlled active flows.

Measuring response functions of active materials from data

From flocks of birds to biomolecular assemblies, systems in which many individual components independently consume energy to perform mechanical work exhibit a wide array of striking behaviors. Methods to quantify the dynamics of these so-called active systems generally aim to extract important length or time scales from experimental fields. Because such methods focus on extracting scalar values, they do not wring maximal information from experimental data. We introduce a method to overcome these limitations. We extend the framework of correlation functions by taking into account the internal headings of displacement fields. The functions we construct represent the material response to specific types of active perturbation within the system. Utilizing these response functions we query the material response of disparate active systems composed of actin filaments and myosin motors, from model fluids to living cells. We show we can extract critical length scales from the turbulent flows of an active nematic, anticipate contractility in an active gel, distinguish viscous from viscoelastic dissipation, and even differentiate modes of contractility in living cells. These examples underscore the vast utility of this method which measures response functions from experimental observations of complex active systems.

Spontaneous flows and dynamics of full-integer topological defects in polar active matter

Polar active matter of self-propelled particles sustain spontaneous flows through the full-integer topological defects. We study theoretically the incompressible flow profiles around ±1 defects induced by polar and dipolar active forces. We show that dipolar forces induce vortical flows around the +1 defect, while the flow around the -1 defect has an 8-fold rotational symmetry. The vortical flow changes its chirality near the +1 defect core in the absence of the friction with a substrate. We show analytically that the flow induced by polar active forces is vortical near the +1 defect and is 4-fold symmetric near the -1 defect, while it becomes uniform in the far-field. For a pair of oppositely charged defects, this polar flow contributes to a mutual interaction force that depends only on the orientation of the defect pair relative to the background polarization, and that enhances defect pair annihilation. This is in contradiction with the effect of dipolar active forces which decay inversely proportional with the defect separation distance. As such, our analyses reveals a long-ranged mechanism for the pairwise interaction between topological defects in polar active matter.

Flow patterns and defect dynamics of active nematic liquid crystals under an electric field

The effects of an electric field on the flow patterns and defect dynamics of two-dimensional active nematic liquid crystals are numerically investigated. We found that field-induced director reorientation causes anisotropic active turbulence characterized by enhanced flow perpendicular to the electric field. The average flow speed and its anisotropy are maximized at an intermediate field strength. Topological defects in the anisotropic active turbulence are localized and show characteristic dynamics with simultaneous creation of two pairs of defects. A laning state characterized by stripe domains with alternating flow directions is found at a larger field strength near the transition to the uniformly aligned state. We obtained periodic oscillations between the laning state and active turbulence, which resembles an experimental observation of active nematics subject to anisotropic friction.

Complex-tensor theory of simple smectics

Matter self-assembling into layers generates unique properties, including structures of stacked surfaces, directed transport, and compact area maximization that can be highly functionalized in biology and technology. Smectics represent the paradigm of such lamellar materials - they are a state between fluids and solids, characterized by both orientational and partial positional ordering in one layering direction, making them notoriously difficult to model, particularly in confining geometries. We propose a complex tensor order parameter to describe the local degree of lamellar ordering, layer displacement and orientation of the layers for simple, lamellar smectics. The theory accounts for both dislocations and disclinations, by regularizing singularities within defect cores and so remaining continuous everywhere. The ability to describe disclinations and dislocation allows this theory to simulate arrested configurations and inclusion-induced local ordering. This tensorial theory for simple smectics considerably simplifies numerics, facilitating studies on the mesoscopic structure of topologically complex systems.

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

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

Active Bending of Disordered Microtubule Bundles by Kinesin Motors

Active networks of biopolymers and motor proteins in vitro self-organize and exhibit dynamic structures on length scales much larger than the interacting individual components of which they consist. How the dynamics is related across the range of length scales is still an open question. Here, we experimentally characterize and quantify the dynamic behavior of isolated microtubule bundles that bend due to the activity of motor proteins. At the motor level, we track and describe the motion features of kinesin-1 clusters stepping within the bending bundles. We find that there is a separation of length scales by at least 1 order of magnitude. At a run length of <1 μm, kinesin-1 activity leads to a bundle curvature in the range of tens of micrometers. We propose that the distribution of microtubule polarity plays a crucial role in the bending dynamics that we observe at both the bundle and motor levels. Our results contribute to the understanding of fundamental principles of vital intracellular processes by disentangling the multiscale dynamics in out-of-equilibrium active networks composed of cytoskeletal elements.

Active boundary layers in confined active nematics

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

Tubulation and dispersion of oil by bacterial growth on droplets

Bacteria on surfaces exhibit collective behaviors, such as active turbulence and active stresses, which result from their motility, growth, and interactions with their local surroundings. However, interfacial deformations on soft surfaces and liquid interfaces caused by active growth, particularly over long time scales, are not well understood. Here, we describe experimental observations on the emergence of tubular structures arising from the growth of rod-shaped bacteria at the interface of oil droplets in water. Using microfluidics and timelapse microscopy, the dimensions and extension rates of individual tubular structures as well as bulk bio-aggregate formation are quantified for hundreds of droplets over 72 hours. Tubular structures are comparable in length to the initial droplet radius and are composed of an outer shell of bacteria that stabilize an inner filament of oil. The oil filament breaks up into smaller microdroplets dispersed within the bacterial shell. This work provides insight into active stresses at deformable interfaces and improves our understanding of microbial oil biodegradation and its potential influence on the transport of droplets in the ocean water column.

Nematic order condensation and topological defects in inertial active nematics

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

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

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

Integer topological defects organize stresses driving tissue morphogenesis

Tissues acquire function and shape via differentiation and morphogenesis. Both processes are driven by coordinating cellular forces and shapes at the tissue scale, but general principles governing this interplay remain to be discovered. Here we report that self-organization of myoblasts around integer topological defects, namely spirals and asters, suffices to establish complex multicellular architectures. In particular, these arrangements can trigger localized cell differentiation or, alternatively, when differentiation is inhibited, they can drive the growth of swirling protrusions. Both localized differentiation and growth of cellular vortices require specific stress patterns. By analysing the experimental velocity and orientational fields through active gel theory, we show that integer topological defects can generate force gradients that concentrate compressive stresses. We reveal these gradients by assessing spatial changes in nuclear volume and deformations of elastic pillars. We propose integer topological defects as mechanical organizing centres controlling differentiation and morphogenesis.

Logic operations with active topological defects

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

Nano/Micromotors in Active Matter

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

Flow around topological defects in active nematic films

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

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.

Two-dimensional numerical model of Marangoni surfers: From single swimmer to crystallization

We numerically study the dynamics of an ensemble of Marangoni surfers in a two-dimensional and unconfined space. The swimmers are modeled as Gaussian sources of surfactant generating surface tension gradients and are shown to follow the Marangoni flow filtered at their spatial scale in the lubrication regime, an unstable situation leading to spontaneous motion as soon as the Marangoni effect is intense enough. As the system is fully unconstrained, it is possible to study the various dynamical regimes from single swimmer, two-body interaction, to the many-particles case characterized by an efficient particle dispersion. We show that, although the present model is very simple, it reproduces the experimentally observed transition between a regime of dispersion by random agitation when the number of swimmers is moderate to the regime of crystallization with imperfect hexagonal lattice at high density.

Orientational Correlations in Active and Passive Nematic Defects

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

Physics of liquid crystals in cell biology

The past decade has witnessed a rapid growth in understanding of the pivotal roles of mechanical stresses and physical forces in cell biology. As a result, an integrated view of cell biology is evolving, where genetic and molecular features are scrutinised hand in hand with physical and mechanical characteristics of cells. Physics of liquid crystals has emerged as a burgeoning new frontier in cell biology over the past few years, fuelled by an increasing identification of orientational order and topological defects in cell biology, spanning scales from subcellular filaments to individual cells and multicellular tissues. Here, we provide an account of the most recent findings and developments, together with future promises and challenges in this rapidly evolving interdisciplinary research direction.

Shear flow of active matter in thin channels

We study the shear flow of active filaments confined in a thin channel for extensile and contractile fibers. We apply the Ericksen-Leslie equations of liquid crystal flow with an activity source term. The dimensionless form of this system includes the Ericksen, activity, and Reynolds numbers, together with the aspect ratio of the channel, as the main driving parameters. We perform a normal mode stability analysis of the base shear flow. For both types of fibers, we arrive at a comprehensive description of the stability properties and their dependence on the parameters of the system. The transition to unstable frequencies in extensile fibers occurs at a positive threshold value of the activity parameter, whereas for contractile ones a complex behavior is found at low absolute value of the activity number. The latter might be an indication of the biologically relevant plasticity and phase transition issues. In contrast with extensile fibers, flows of contractile ones are also found to be highly sensitive to the Reynolds number. The work on extensile fibers is guided by experiments on active filaments in confined channels and aims at quantifying their findings in the prechaotic regime.

Multi-particle collision dynamics with a non-ideal equation of state. II. Collective dynamics of elongated squirmer rods

Simulations of flow fields around microscopic objects typically require methods that both solve the Navier-Stokes equations and also include thermal fluctuations. One such method popular in the field of soft-matter physics is the particle-based simulation method of multi-particle collision dynamics (MPCD). However, in contrast to the typically incompressible real fluid, the fluid of the traditional MPCD methods obeys the ideal-gas equation of state. This can be problematic because most fluid properties strongly depend on the fluid density. In a recent article, we proposed an extended MPCD algorithm and derived its non-ideal equation of state and an expression for the viscosity. In the present work, we demonstrate its accuracy and efficiency for the simulations of the flow fields of single squirmers and of the collective dynamics of squirmer rods. We use two exemplary squirmer-rod systems for which we compare the outcome of the extended MPCD method to the well-established MPCD version with an Andersen thermostat. First, we explicitly demonstrate the reduced compressibility of the MPCD fluid in a cluster of squirmer rods. Second, for shorter rods, we show the interesting result that in simulations with the extended MPCD method, dynamic swarms are more pronounced and have a higher polar order. Finally, we present a thorough study of the state diagram of squirmer rods moving in the center plane of a Hele-Shaw geometry. From a small to large aspect ratio and density, we observe a disordered state, dynamic swarms, a single swarm, and a jammed cluster, which we characterize accordingly.

Multiscale Microtubule Dynamics in Active Nematics

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

Submersed micropatterned structures control active nematic flow, topology, and concentration

Coupling between flows and material properties imbues rheological matter with its wide-ranging applicability, hence the excitement for harnessing the rheology of active fluids for which internal structure and continuous energy injection lead to spontaneous flows and complex, out-of-equilibrium dynamics. We propose and demonstrate a convenient, highly tunable method for controlling flow, topology, and composition within active films. Our approach establishes rheological coupling via the indirect presence of fully submersed micropatterned structures within a thin, underlying oil layer. Simulations reveal that micropatterned structures produce effective virtual boundaries within the superjacent active nematic film due to differences in viscous dissipation as a function of depth. This accessible method of applying position-dependent, effective dissipation to the active films presents a nonintrusive pathway for engineering active microfluidic systems.

Properties of twisted topological defects in 2D nematic liquid crystals

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

Flow coupling between active and passive fluids across water-oil interfaces

Active fluid droplets surrounded by oil can spontaneously develop circulatory flows. However, the dynamics of the surrounding oil and their influence on the active fluid remain poorly understood. To investigate interactions between the active fluid and the passive oil across their interface, kinesin-driven microtubule-based active fluid droplets were immersed in oil and compressed into a cylinder-like shape. The droplet geometry supported intradroplet circulatory flows, but the circulation was suppressed when the thickness of the oil layer surrounding the droplet decreased. Experiments with tracers and network structure analyses and continuum models based on the dynamics of self-elongating rods demonstrated that the flow transition resulted from flow coupling across the interface between active fluid and oil, with a millimeter-scale coupling length. In addition, two novel millifluidic devices were developed that could trigger and suppress intradroplet circulatory flows in real time: one by changing the thickness of the surrounding oil layer and the other by locally deforming the droplet. This work highlights the role of interfacial dynamics in the active fluid droplet system and shows that circulatory flows within droplets can be affected by millimeter-scale flow coupling across the interface between the active fluid and the oil.

Tuning the Properties of Active Microtubule Networks by Depletion Forces

Suspensions of microtubules and nonadsorbing particles form thick and long bundles due to depletion forces. Such interactions act at the nanometer scale and define the structural and dynamical properties of the resulting networks. In this study, we analyze the depletion forces exerted by two types of nonadsorbing particles, namely, the polymer, poly(ethylene glycol) (PEG), and the block copolymer, Pluronic. We characterize their effects both in passive and active networks by adding motor proteins to the suspensions. By exploiting its bundling effect via entropic forces, we observed that PEG generates a network with thick structures showing a nematic order and larger mesh size. On the other hand, Pluronic builds up a much denser gel-like network without a recognizable mesh structure. This difference is also reflected in the network activity. PEG networks show moderate contraction in lateral directions while Pluronic networks exhibit faster and isotropic contraction. Interestingly, by mixing the two nonadsorbing polymers in different ratios, we observed that the system showed a behavior that exhibited properties of both agents, leading to a robust and fast responsive structure compared to the single-depletant networks. In conclusion, we show how passive osmotic compression modifies the distribution of biopolymers. Its combination with active motors results in a new active material with potential for nanotechnological applications.

Fractal generation in a two-dimensional active-nematic fluid

Active fluids, composed of individual self-propelled agents, can generate complex large-scale coherent flows. A particularly important laboratory realization of such an active fluid is a system composed of microtubules, aligned in a quasi-two-dimensional (2D) nematic phase and driven by adenosine-triphosphate-fueled kinesin motor proteins. This system exhibits robust chaotic advection and gives rise to a pronounced fractal structure in the nematic contours. We characterize such experimentally derived fractals using the power spectrum and discover that the power spectrum decays as k^-β for large wavenumbers k. The parameter β is measured for several experimental realizations. Though β is effectively constant in time, it does vary with experimental parameters, indicating differences in the scale-free behavior of the microtubule-based active nematic. Though the fractal patterns generated in this active system are reminiscent of passively advected dye in 2D chaotic flows, the underlying mechanism for fractal generation is more subtle. We provide a simple, physically inspired mathematical model of fractal generation in this system that relies on the material being locally compressible, though the total area of the material is conserved globally. The model also requires that large-scale density variations are injected into the material periodically. The model reproduces the power-spectrum decay k^-β seen in experiments. Linearizing the model of fractal generation about the equilibrium density, we derive an analytic relationship between β and a single dimensionless quantity r, which characterizes the compressibility.

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

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

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.

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.

Machine learning forecasting of active nematics

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

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.

Stability analysis of flow of active extensile fibers in confined domains

In this article, we study shear flow of active extensile filaments confined in a narrow channel. They behave as nematic liquid crystals that we assumed are governed by the Ericksen-Leslie equations of balance of linear and angular momentum. The addition of an activity source term in the Leslie stress captures the role of the biofuel prompting the dynamics. The dimensionless form of the governing system includes the Ericksen, activity, and Reynolds numbers together with the aspect ratio of the channel as the main driving parameters affecting the stability of the system. The active system that guides our analysis is composed of microtubules concentrated in bundles, hundreds of microns long, placed in a narrow channel domain, of aspect ratios in the range between 10^-2 and 10^-3 dimensionless units, which are able to align due to the combination of adenosine triphosphate-supplied energy and confinement effects. Specifically, this work aims at studying the role of confinement on the behavior of active matter. It is experimentally observed that, at an appropriately low activity and channel width, the active flow is laminar, with the linear velocity profile and the angle of alignment analogous to those in passive shear, developing defects and becoming chaotic, at a large activity and a channel aspect ratio. The present work addresses the laminar regime, where defect formation does not play a role. We perform a normal mode stability analysis of the base shear flow. A comprehensive description of the stability properties is obtained in terms of the driving parameters of the system. Our main finding, in addition to the geometry and magnitude of the flow profiles, and also consistent with the experimental observations, is that the transition to instability of the uniformly aligned shear flow occurs at a threshold value of the activity parameter, with the transition also being affected by the channel aspect ratio. The role of the parameters on the vorticity and angular profiles of the perturbing flow is also analyzed and found to agree with the experimentally observed transition to turbulent regimes. A spectral method based on Chebyshev polynomials is used to solve the generalized eigenvalue problems arising in the stability analysis.

Optimal Control of Active Nematics

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

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.

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.

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.