Confinement Controls the Bend Instability of Three-Dimensional Active Liquid Crystals

Active nematics in corrugated channels

Active nematic fluids exhibit complex dynamics in both bulk and in simple confining geometries. However, complex confining geometries could have substantial impact on active spontaneous flows. Using multiparticle collision dynamics simulations adapted for active nematic particles, we study the dynamic behaviour of an active nematic fluid confined in a corrugated channel. The transition from a quiescent state to a spontaneous flow state occurs from a weak swirling flow to a strong coherent flow due to the presence of curved-wall induced active flows. We show that the active nematic fluid flows in corrugated channels can be understood in two different ways: (i) as the result of an early or delayed flow transition when compared with that in a flat-walled channel of appropriate width and (ii) boundary-induced active flows in the corrugations providing an effective slip velocity to the coherent flows in the bulk. Thus, our work illustrates the crucial role of corrugations of the confining boundary in dictating the flow transition and flow states of active fluids.

Scaling Transition of Active Turbulence from Two to Three Dimensions

Turbulent flows are observed in low-Reynolds active fluids, which display similar phenomenology to the classical inertial turbulence but are of a different nature. Understanding the dependence of this new type of turbulence on dimensionality is a fundamental challenge in non-equilibrium physics. Real-space structures and kinetic energy spectra of bacterial turbulence are experimentally measured from two to three dimensions. The turbulence shows three regimes separated by two critical confinement heights, resulting from the competition of bacterial length, vortex size and confinement height. Meanwhile, the kinetic energy spectra display distinct universal scaling laws in quasi-2D and 3D regimes, independent of bacterial activity, length, and confinement height, whereas scaling exponents transition in two steps around the critical heights. The scaling behaviors are well captured by the hydrodynamic model we develop, which employs image systems to represent the effects of confining boundaries. The study suggests a framework for investigating the effect of dimensionality on non-equilibrium self-organized systems.

Heterogeneous distribution of kinesin-streptavidin complexes revealed by mass photometry

Kinesin-streptavidin complexes are widely used in microtubule-based active-matter studies. The stoichiometry of the complexes is empirically tuned but experimentally challenging to determine. Here, mass photometry measurements reveal heterogenous distributions of kinesin-streptavidin complexes. Our binding model indicates that heterogeneity arises from both the kinesin-streptavidin mixing ratio and the kinesin-biotinylation efficiency.

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.

Flow states of two dimensional active gels driven by external shear

Using a minimal hydrodynamic model, we theoretically and computationally study the Couette flow of active gels in straight and annular two-dimensional channels subject to an externally imposed shear. The gels are isotropic in the absence of externally- or activity-driven shear, but have nematic order that increases with shear rate. Using the finite element method, we determine the possible flow states for a range of activities and shear rates. Linear stability analysis of an unconfined gel in a straight channel shows that an externally imposed shear flow can stabilize an extensile fluid that would be unstable to spontaneous flow in the absence of the shear flow, and destabilize a contractile fluid that would be stable against spontaneous flow in the absence of shear flow. These results are in rough agreement with the stability boundaries between the base shear flow state and the nonlinear flow states that we find numerically for a confined active gel. For extensile fluids, we find three kinds of nonlinear flow states in the range of parameters we study: unidirectional flows, oscillatory flows, and dancing flows. To highlight the activity-driven spontaneous component of the nonlinear flows, we characterize these states by the average volumetric flow rate and the wall stress. For contractile fluids, we only find the linear shear flow and a nonlinear unidirectional flow in the range of parameters that we studied. For large magnitudes of the activity, the unidirectional contractile flow develops a boundary layer. Our analysis of annular channels shows how curvature of the streamlines in the base flow affects the transitions among flow states.

Microscopic interactions control a structural transition in active mixtures of microtubules and molecular motors

Microtubules and molecular motors are essential components of the cellular cytoskeleton, driving fundamental processes in vivo, including chromosome segregation and cargo transport. When reconstituted in vitro, these cytoskeletal proteins serve as energy-consuming building blocks to study the self-organization of active matter. Cytoskeletal active gels display rich emergent dynamics, including extensile flows, locally contractile asters, and bulk contraction. However, it is unclear how the protein-protein interaction kinetics set their contractile or extensile nature. Here, we explore the origin of the transition from extensile bundles to contractile asters in a minimal reconstituted system composed of stabilized microtubules, depletant, adenosine 5'-triphosphate (ATP), and clusters of kinesin-1 motors. We show that the microtubule-binding and unbinding kinetics of highly processive motor clusters set their ability to end-accumulate, which can drive polarity sorting of the microtubules and aster formation. We further demonstrate that the microscopic time scale of end-accumulation sets the emergent time scale of aster formation. Finally, we show that biochemical regulation is insufficient to fully explain the transition as generic aligning interactions through depletion, cross-linking, or excluded volume interactions can drive bundle formation despite end-accumulating motors. The extensile-to-contractile transition is well captured by a simple self-assembly model where nematic and polar aligning interactions compete to form either bundles or asters. Starting from a five-dimensional organization phase space, we identify a single control parameter given by the ratio of the different component concentrations that dictates the material-scale organization. Overall, this work shows that the interplay of biochemical and mechanical tuning at the microscopic level controls the robust self-organization of active cytoskeletal materials.

Viscoelastic confinement induces periodic flow reversals in active nematics

We use linear stability analysis and hybrid lattice Boltzmann simulations to study the dynamical behavior of an active nematic confined in a channel made of viscoelastic material. We find that the quiescent, ordered active nematic is unstable above a critical activity. The transition is to a steady flow state for high elasticity of the channel surroundings. However, below a threshold elastic modulus, the system produces spontaneous oscillations with periodic flow reversals. We provide a phase diagram that highlights the region where time-periodic oscillations are observed and explain how they are produced by the interplay of activity and viscoelasticity. Our results suggest experiments to study the role of viscoelastic confinement in the spatiotemporal organization and control of active matter.

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.

Light-activated microtubule-based two-dimensional active nematic

We assess the ability of two light responsive kinesin motor clusters to drive dynamics of microtubule-based active nematics: opto-K401, a processive motor, and opto-K365, a non-processive motor. Measurements reveal an order of magnitude improvement in the contrast of nematic flow speeds between maximally- and minimally-illuminated states for opto-K365 motors when compared to opto-K401 construct. For opto-K365 nematics, we characterize both the steady-state flow and defect density as a function of applied light. We also examine the transient behavior as the system switches between steady-states upon changes in light intensities. Although nematic flows reach a steady state within tens of seconds, the defect density exhibits transient behavior for up to 10 minutes, showing a separation between small-scale active flows and system-scale structural states. Our work establishes an experimental platform that can exploit spatiotemporally-heterogeneous patterns of activity to generate targeted dynamical states.

Compressive forces stabilize microtubules in living cells

Microtubules are cytoskeleton components with unique mechanical and dynamic properties. They are rigid polymers that alternate phases of growth and shrinkage. Nonetheless, the cells can display a subset of stable microtubules, but it is unclear whether microtubule dynamics and mechanical properties are related. Recent in vitro studies suggest that microtubules have mechano-responsive properties, being able to stabilize their lattice by self-repair on physical damage. Here we study how microtubules respond to cycles of compressive forces in living cells and find that microtubules become distorted, less dynamic and more stable. This mechano-stabilization depends on CLASP2, which relocates from the end to the deformed shaft of microtubules. This process seems to be instrumental for cell migration in confined spaces. Overall, these results demonstrate that microtubules in living cells have mechano-responsive properties that allow them to resist and even counteract the forces to which they are subjected, being a central mediator of cellular mechano-responses.

Spatio-temporal patterning of extensile active stresses in microtubule-based active fluids

Microtubule-based active fluids exhibit turbulent-like autonomous flows, which are driven by the molecular motor powered motion of filamentous constituents. Controlling active stresses in space and time is an essential prerequisite for controlling the intrinsically chaotic dynamics of extensile active fluids. We design single-headed kinesin molecular motors that exhibit optically enhanced clustering and thus enable precise and repeatable spatial and temporal control of extensile active stresses. Such motors enable rapid, reversible switching between flowing and quiescent states. In turn, spatio-temporal patterning of the active stress controls the evolution of the ubiquitous bend instability of extensile active fluids and determines its critical length dependence. Combining optically controlled clusters with conventional kinesin motors enables one-time switching from contractile to extensile active stresses. These results open a path towards real-time control of the autonomous flows generated by active fluids.

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.

Competing instabilities reveal how to rationally design and control active crosslinked gels

How active stresses generated by molecular motors set the large-scale mechanics of the cell cytoskeleton remains poorly understood. Here, we combine experiments and theory to demonstrate how the emergent properties of a biomimetic active crosslinked gel depend on the properties of its microscopic constituents. We show that an extensile nematic elastomer exhibits two distinct activity-driven instabilities, spontaneously bending in-plane or buckling out-of-plane depending on its composition. Molecular motors play a dual antagonistic role, fluidizing or stiffening the gel depending on the ATP concentration. We demonstrate how active and elastic stresses are set by each component, providing estimates for the active gel theory parameters. Finally, activity and elasticity were manipulated in situ with light-activable motor proteins, controlling the direction of the instability optically. These results highlight how cytoskeletal stresses regulate the self-organization of living matter and set the foundations for the rational design and optogenetic control of active materials.

Catapulting of topological defects through elasticity bands in active nematics

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

Active Microphase Separation in Mixtures of Microtubules and Tip-Accumulating Molecular Motors

Mixtures of filaments and molecular motors form active materials with diverse dynamical behaviors that vary based on their constituents' molecular properties. To develop a multiscale of these materials, we map the nonequilibrium phase diagram of microtubules and tip-accumulating kinesin-4 molecular motors. We find that kinesin-4 can drive either global contractions or turbulentlike extensile dynamics, depending on the concentrations of both microtubules and a bundling agent. We also observe a range of spatially heterogeneous nonequilibrium phases, including finite-sized radial asters, 1D wormlike chains, extended 2D bilayers, and system-spanning 3D active foams. Finally, we describe intricate kinetic pathways that yield microphase separated structures and arise from the inherent frustration between the orientational order of filamentous microtubules and the positional order of tip-accumulating molecular motors. Our work reveals a range of novel active states. It also shows that the form of active stresses is not solely dictated by the properties of individual motors and filaments, but is also contingent on the constituent concentrations and spatial arrangement of motors on the filaments.

Helical flow states in active nematics

We show that confining extensile nematics in three-dimensional (3D) channels leads to the emergence of two self-organized flow states with nonzero helicity. The first is a pair of braided antiparallel streams-this double helix occurs when the activity is moderate, anchoring negligible, and reduced temperature high. The second consists of axially aligned counter-rotating vortices-this grinder train arises between spontaneous axial streaming and the vortex lattice. These two unanticipated helical flow states illustrate the potential of active fluids to break symmetries and form complex but organized spatiotemporal structures in 3D fluidic devices.

Defect loops in three-dimensional active nematics as active multipoles

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

Crosslinking and depletion determine spatial instabilities in cytoskeletal active matter

Active gels made of cytoskeletal proteins are valuable materials with attractive non-equilibrium properties such as spatial self-organization and self-propulsion. At least four typical routes to spatial patterning have been reported to date in different types of cytoskeletal active gels: bending and buckling instabilities in extensile systems, and global and local contraction instabilities in contractile gels. Here we report the observation of these four instabilities in a single type of active gel and we show that they are controlled by two parameters: the concentrations of ATP and depletion agent. We demonstrate that as the ATP concentration decreases, the concentration of passive motors increases until the gel undergoes a gelation transition. At this point, buckling is selected against bending, while global contraction is favored over local ones. Our observations are coherent with a hydrodynamic model of a viscoelastic active gel where the filaments are crosslinked with a characteristic time that diverges as the ATP concentration decreases. Our work thus provides a unified view of spatial instabilities in cytoskeletal active matter.

Active Extensile Stress Promotes 3D Director Orientations and Flows

We use numerical simulations and linear stability analysis to study an active nematic layer where the director is allowed to point out of the plane. Our results highlight the difference between extensile and contractile systems. Contractile stress suppresses the flows perpendicular to the layer and favors in-plane orientations of the director. By contrast extensile stress promotes instabilities that can turn the director out of the plane, leaving behind a population of distinct, in-plane regions that continually elongate and divide. This supports extensile forces as a mechanism for the initial stages of layer formation in living systems, and we show that a planar drop with extensile (contractile) activity grows into three dimensions (remains in two dimensions). The results also explain the propensity of disclination lines in three dimensional active nematics to be of twist type in extensile or wedge type in contractile materials.

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.

Effects of confinement on the dynamics and correlation scales in kinesin-microtubule active fluids

We study the influence of solid boundaries on dynamics and structure of kinesin-driven microtubule active fluids as the height of the container, H, increases from hundreds of micrometers to several millimeters. By three-dimensional tracking of passive tracers dispersed in the active fluid, we observe that the activity level, characterized by velocity fluctuations, increases as system size increases and retains a small-scale isotropy. Concomitantly, as the confinement level decreases, the velocity-velocity temporal correlation develops a strong positive correlation at longer times, suggesting the establishment of a "memory". We estimate the characteristic size of the flow structures from the spatial correlation function and find that, as the confinement becomes weaker, the correlation length, l_{c}, saturates at approximately 400 microns. This saturation suggests an intrinsic length scale which, along with the small-scale isotropy, demonstrates the multiscale nature of this kinesin-driven bundled microtubule active system.

Spontaneous circulation of active microtubules confined by optical traps

We propose an experiment to demonstrate spontaneous ordering and symmetry breaking of kinesin-driven microtubules confined to an optical trap. Calculations involving the feasibility of such an experiment are first performed which analyze the power needed to confine microtubules and address heating concerns. We then present the results of first-principles simulations of active microtubules confined in such a trap and analyze the types of motion observed by the microtubules as well as the velocity of the surrounding fluid, both near the trap and in the far-field. We find three distinct phases characterized by breaking of distinct symmetries and also analyze the power spectrum of the angular momenta of polymers to further quantify the differences between these phases. Under the correct conditions, microtubules were found to spontaneously align with one another and circle the trap in one direction.

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

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

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.

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.