Spontaneous motion in hierarchically assembled active matter

Biomimetic swarm of active particles with coupled passive-active interactions

We study the universal behavior of a class of active colloids whose design is inspired by the collective dynamics of natural systems such as schools of fish and flocks of birds. These colloids, with off-center repulsive interaction sites, self-organize into polar swarms exhibiting long-range order and directional motion without significant hydrodynamic interactions. Our simulations show that the system transitions from motile perfect crystals to solid-like, liquid-like, and gas-like states depending on noise levels, repulsive interaction strength, and particle density. By analyzing swarm polarity and hexatic bond order parameters, we demonstrate that effective volume fractions based on force-range and torque-range interactions explain the system's universal behavior. This work lays a groundwork for biomimetic applications utilizing the cooperative polar dynamics of active colloids.

Morphogenesis of confined biofilms: how mechanical interactions determine cellular patterning and global geometry

Biofilms are surface-attached bacterial communities encased within extracellular matrices that play significant roles in health and society and serve as prototypical examples of proliferating active nematics. Recent advances in fluorescence microscopy have facilitated an unprecedented view of biofilm development at the single-cell level, thus providing the opportunity to develop a mechanistic understanding of how biofilm development is influenced by mechanical interactions with the environment. Here, we review recent studies that examined the morphogenesis of Vibrio cholerae biofilms under confinement at both single-cell and continuum levels. We describe how biofilms under different confinement modes-embedded and interstitial-can acquire various global geometries and forms of cell orientational ordering different from those in unconfined biofilms, and we demonstrate how these properties arise from the mechanical interplay between the biofilm and its confining medium. We also discuss how this interplay is fundamentally governed by the extracellular matrix, which facilitates the transmission of mechanical stress from the medium into the biofilm via adhesion and friction, and serves as the key feature that distinguishes biofilms from classical bacterial colonies. These studies lay the groundwork for many potential future directions, all of which will contribute to the establishment of a new "developmental biology" of biofilms.

Topology controls flow patterns in active double emulsions

Active emulsions and liquid crystalline shells are intriguing and experimentally realisable types of topological matter. Here we numerically study the morphology and spatiotemporal dynamics of a double emulsion, where one or two passive small droplets are embedded in a larger active droplet. We find activity introduces a variety of rich and nontrivial nonequilibrium states in the system. First, a double emulsion with a single active droplet becomes self-motile, and there is a transition between translational and rotational motion: both of these regimes remain defect-free, hence topologically trivial. Second, a pair of particles nucleate one or more disclination loops, with conformational dynamics resembling a rotor or chaotic oscillator, accessed by tuning activity. In the first state a single, topologically charged, disclination loop powers the rotation. In the latter state, this disclination stretches and writhes in 3D, continuously undergoing recombination to yield an example of an active living polymer. These emulsions can be self-assembled in the lab, and provide a pathway to form flow and topological patterns in active matter in a controlled way, as opposed to bulk systems that typically yield active turbulence.

Direct detection of deformation modes on varying length scales in active biopolymer networks

Correlated flows and forces that emerge from active matter orchestrate complex processes such as shape regulation and deformations in biological cells and tissues. The active materials central to cellular mechanics are cytoskeletal networks, where molecular motor activity drives deformations and remodeling. Here, we investigate deformation modes in contractile actin networks driven by the molecular motor myosin II through quantitative fluorescence microscopy. We examine the deformation anisotropy at different length scales in networks of sparsely cross-linked and bundled actin. In sparsely cross-linked networks, we find myosin-dependent biaxial buckling modes across length scales. Interestingly, both long and short-wavelength buckling may contribute to network contractility. In cross-linked bundled networks, uniaxial contraction predominates on long length scales, while the uniaxial or biaxial nature of the deformation depends on bundle microstructure at shorter length scales. The anisotropy of deformations may provide insight to the mechanical origins of contractility in actin networks and regulation of collective behavior in a variety of active materials.

Noise-induced transitions from contractile to extensile active stress in isotropic fluids

Tissues of living cells are a prime example of active fluids. There is experimental evidence that tissues generate extensile active stress even though their constituting cells are contractile. Fluctuating forces that could result from cell-substrate interactions have been proposed to be able to induce a transition from contractile to extensile active stress. We define the notion of contractile and extensile active stress in isotropic and anisotropic active matter. Through analytic calculations and numerical computations, we then show that in isotropic active fluids, nonlinearities and coupling between fluctuating forces and fluid density are necessary for such a transition to occur. Here, both transitions from extensile to contractile and vice versa are possible.

Active droplet driven by collective chemotaxis

Surfactant-laden fluid interfaces of soft colloids, such as bubbles and droplets, are ubiquitously seen in various natural phenomena and industrial settings. In canonical systems where microparticles are driven in hydrodynamic flows, convection of the surfactant changes local surface tension. Subsequently, the interplay of Marangoni and hydrodynamic stresses leads to rich interfacial dynamics that directly impact the particle motions. Here we introduce a new mechanism for self-propelled droplets, driven by a thin layer of chemically active microparticles situated at the interface of a suspended droplet, which is a direct extension of the planar collective surfing model by Masoud and Shelley (H. Masoud and M. J. Shelley, Phys. Rev. Lett., 2014, 112, 128304). These particles can generate chemicals locally, leading to spontaneous Marangoni flows that drive the self-aggregation of microparticles. This process, in turn, creates a polarized surfactant distribution, which induces collective chemotaxis and dipolar bulk flows, ultimately breaking the symmetry. By assuming the local surfactant production to be either proportional to particle density or saturated at a high particle density, we observe that the system can be chemotactically diverging or approach a steady state with constant migration velocity. The system is studied analytically in the linear region for the initial transient dynamics, yielding critical numbers and familiar patterns, as well as numerically for larger amplitudes and over a long time using spectral methods.

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.

Force balance of opposing diffusive motors generates polarity-sorted microtubule patterns

The internal organization of cells is largely determined by the architecture and orientation of the microtubule network. Microtubules serve as polar tracks for the selective transport of specific molecular motors toward either their plus or minus ends. How both motors reciprocally move microtubules and organize the network's arrangement and polarity is unknown. Here, we combined experiments on reconstituted systems and theory to study the interaction of microtubules with both plus- and minus-end directed motors bound to a fluid membrane. Depending on motor concentrations, the system could lead either to the constant transport of microtubules or to their alignment, stacking, and immobilization in regular bands that separate motors into domains of opposite polarities. In bands, microtubules shared the same polarity and segregated the two opposing motors accordingly. These regular patterns resulted from the balance of forces produced by the two motors as they walked in opposite directions along microtubules. The system was maintained in a dynamic steady state in which the directional transport of microtubule-bound motors compensates for the random diffusion of lipid-bound motors. The size of motor domains depended on their respective concentrations. The constant flow of motors allowed the system to respond to variations in motor concentrations by moving microtubules to adapt to the new force balance. The polar sorting and linear arrangement of microtubules associated with the segregation of motors of opposite polarity are typical of cellular architectures, which these data may help to better understand.

Topological transition in filamentous cyanobacteria: from motion to structure

Many active systems are capable of forming intriguing patterns at scales significantly larger than the size of their individual constituents. Cyanobacteria are one of the most ancient and important phyla of organisms that has allowed the evolution of more complex life forms. Despite its importance, the role of motility on the pattern formation of their colonies is not understood. Here, we investigate the large-scale collective effects and rich dynamics of gliding filamentous cyanobacteria colonies, while still retaining information about the individual constituents' dynamics and their interactions. We investigate both the colony's transient and steady-state dynamics and find good agreement with experiments. We furthermore show that the Péclet number and aligning interaction strength govern the system's topological transition from an isotropic distribution to a state of large-scale reticulate patterns. Although the system is topologically non-trivial, the parallel and perpendicular pair correlation functions provide structural information about the colony, and thus can be used to extract information about the early stages of biofilm formation. Finally, we find that the effects of the filaments' length cannot be reduced to a system of interacting points. Our model proves to reproduce both cyanobacteria colonies and systems of biofilaments where curvature is transported by motility.

Spontaneous flows and quantum analogies in heterogeneous active nematic films

Incorporating the inherent heterogeneity of living systems into models of active nematics is essential to provide a more realistic description of biological processes such as bacterial growth, cell dynamics and tissue development. Spontaneous flow of a confined active nematic is a fundamental feature of these systems, in which the role of heterogeneity has not yet been considered. We therefore determine the form of spontaneous flow transition for an active nematic film with heterogeneous activity, identifying a correspondence between the unstable director modes and solutions to Schrödinger's equation. We consider both activity gradients and steps between regions of distinct activity, finding that such variations can change the signature properties of the flow. The threshold activity required for the transition can be raised or lowered, the fluid flux can be reduced or reversed and interfaces in activity induce shear flows. In a biological context fluid flux influences the spread of nutrients while shear flows affect the behaviour of rheotactic microswimmers and can cause the deformation of biofilms. All the effects we identify are found to be strongly dependent on not simply the types of activity present in the film but also on how they are distributed.

Multi-scaled regulation for cholesteric organization of cellulose nanocrystals based on internal and external factors

Cellulose nanocrystals (CNCs), as one of the most promising bio-sourced materials, have been drawing increasing attention as they offer an attractive route for the rational design and sustainable manufacturing of photonic materials owing to their cholesteric self-assembly from the suspension to solid state. Such an organization process can be readily regulated depending on either internal factors or external forces. In this review, recent advances in the control over the self-assembly process and photonic organizations of CNCs are summarized. Furthermore, an in-depth understanding of diverse factors affecting the nano-scaled periodicity and micro-scaled alignment of CNC cholesteric organization is obtained from perspectives of bulk building blocks, solution environment, extra additives, and external forces. Additionally, the roles of the multi-sized photonic architecture associated with photonic-photonic coupling and the macrogeometry related to the complex confined self-assembly are highlighted for sophisticated CNC optical materials. Finally, insights into the future challenges in the field of CNC photonic materials, regarding the precise mechanism of CNC self-assembly and translation of CNC photonic technology from academia to industry, are proposed.

Controlling wall-particle interactions with activity

We theoretically determine the effective forces on hard disks near walls embedded inside active nematic liquid crystals. When the disks are sufficiently close to the wall and the flows are sufficiently slow, we can obtain exact expressions for the effective forces. We find these forces and the dynamics of disks near the wall depend both on the properties of the active nematic and on the anchoring conditions on the disks and the wall. Our results show that the presence of active stresses attract planar anchored disks to walls if the activity is extensile, and repel them if contractile. For normal anchored disks the reverse is true; they are attracted in contractile systems, and repelled in extensile ones. By choosing the activity and anchoring, these effects may be helpful in controlling the self assembly of active nematic colloids.

Regulation of minimal spindle midzone organization by mitotic kinases

During cell division, the microtubule cytoskeleton undergoes dramatic cell cycle-driven reorganizations of its architecture. Coordinated by changes in the phosphorylation patterns of a multitude of microtubule associated proteins, the mitotic spindle first self-assembles to capture the chromosomes and then reorganizes in anaphase as the chromosomes are segregated. A key protein for this reorganization is PRC1 which is differentially phosphorylated by the mitotic kinases CDK1 and PLK1. How the phosphorylation state of PRC1 orchestrates spindle reorganization is not understood mechanistically. Here, we reconstitute in vitro the transition between metaphase and anaphase-like microtubule architectures triggered by the changes in PRC1 phosphorylation. We find that whereas PLK1 regulates its own recruitment by PRC1, CDK1 controls the affinity of PRC1 for antiparallel microtubule binding. Dephosphorylation of CDK1-phosphorylated PRC1 is required and sufficient to trigger the reorganization of a minimal anaphase midzone in the presence of the midzone length controlling kinesin KIF4A. These results demonstrate how phosphorylation-controlled affinity changes regulate the architecture of active microtubule networks, providing new insight into the mechanistic underpinnings of the cell cycle-driven reorganization of the central spindle during mitosis.

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.

Collective self-caging of active filaments in virtual confinement

Motility coupled to responsive behavior is essential for many microorganisms to seek and establish appropriate habitats. One of the simplest possible responses, reversing the direction of motion, is believed to enable filamentous cyanobacteria to form stable aggregates or accumulate in suitable light conditions. Here, we demonstrate that filamentous morphology in combination with responding to light gradients by reversals has consequences far beyond simple accumulation: Entangled aggregates form at the boundaries of illuminated regions, harnessing the boundary to establish local order. We explore how the light pattern, in particular its boundary curvature, impacts aggregation. A minimal mechanistic model of active flexible filaments resembles the experimental findings, thereby revealing the emergent and generic character of these structures. This phenomenon may enable elongated microorganisms to generate adaptive colony architectures in limited habitats or guide the assembly of biomimetic fibrous materials.

Traveling waves at the surface of active liquid crystals

Active liquid crystals exert nonequilibrium stresses on their surroundings through constant consumption of energy, giving rise to dynamical steady states not present in equilibrium. The paradigmatic example of an active liquid crystal is a suspension of microtubule bundles powered by kinesin motor proteins, which exhibits self-sustained spatiotemporal chaotic flows. This system has been modelled using continuum theories that couple the microtubule orientation to active flows. Recently the focus has shifted to the interfacial properties of mixtures of active liquid crystals and passive fluids. Active/passive interfaces have been shown to support propagating capillary waves in the absence of inertia and offer a promising route for relating experimental parameters to those of the continuum theory. In this paper we report the derivation of a minimal model that captures the linear dynamics of the interface between an active liquid crystal and a passive fluid. We show that the dynamics of the interface, although powered by active flows throughout the bulk, is qualitatively captured by equations that couple non-reciprocally interface height and nematic director at the interface. This minimal model reproduces the dynamical structure factor evaluated from numerical simulations and the qualitative form of the wave dispersion relation seen in experiments.

Anomalous diffusion of active Brownian particles in responsive elastic gels: Nonergodicity, non-Gaussianity, and distributions of trapping times

Understanding actual transport mechanisms of self-propelled particles (SPPs) in complex elastic gels-such as in the cell cytoplasm, in in vitro networks of chromatin or of F-actin fibers, or in mucus gels-has far-reaching consequences. Implications beyond biology/biophysics are in engineering and medicine, with a particular focus on microrheology and on targeted drug delivery. Here, we examine via extensive computer simulations the dynamics of SPPs in deformable gellike structures responsive to thermal fluctuations. We treat tracer particles comparable to and larger than the mesh size of the gel. We observe distinct trapping events of active tracers at relatively short times, leading to subdiffusion; it is followed by an escape from meshwork-induced traps due to the flexibility of the network, resulting in superdiffusion. We thus find crossovers between different transport regimes. We also find pronounced nonergodicity in the dynamics of SPPs and non-Gaussianity at intermediate times. The distributions of trapping times of the tracers escaping from "cages" in our quasiperiodic gel often reveal the existence of two distinct timescales in the dynamics. At high activity of the tracers these timescales become comparable. Furthermore, we find that the mean waiting time exhibits a power-law dependence on the activity of SPPs (in terms of their Péclet number). Our results additionally showcase both exponential and nonexponential trapping events at high activities. Extensions of this setup are possible, with the factors such as anisotropy of the particles, different topologies of the gel network, and various interactions between the particles (also of a nonlocal nature) to be considered.

Supramolecular fibrillation in coacervates and other confined systems towards biomimetic function

As in natural cytoskeletons, the cooperative assembly of fibrillar networks can be hosted inside compartments to engineer biomimetic functions, such as mechanical actuation, transport, and reaction templating. Coacervates impose an optimal liquid-liquid phase separation within the aqueous continuum, functioning as membrane-less compartments that can organise such self-assembling processes as well as the exchange of information with their environment. Furthermore, biological fibrillation can often be controlled or assisted by intracellular compartments. Thus, the reconstitution of analogues of natural filaments in simplified artificial compartments, such as coacervates, offer a suitable model to unravel, mimic, and potentially exploit cellular functions. This perspective summarises the latest developments towards assembling fibrillar networks under confinement inside coacervates and related compartments, including a selection of examples ranging from biological to fully synthetic monomers. Comparative analysis between coacervates, lipid vesicles, and droplet emulsions showcases the interplay between supramolecular fibres and the boundaries of the corresponding compartment. Combining inspiration from natural systems and the custom properties of tailored synthetic fibrillators, rational monomer and compartment design will contribute towards engineering increasingly complex and more realistic artificial protocells.

Deep-learning optical flow for measuring velocity fields from experimental data

Deep learning-based optical flow (DLOF) extracts features in adjacent video frames with deep convolutional neural networks. It uses those features to estimate the inter-frame motions of objects. We evaluate the ability of optical flow to quantify the spontaneous flows of microtubule (MT)-based active nematics under different labeling conditions, and compare its performance to particle image velocimetry (PIV). We obtain flow velocity ground truths either by performing semi-automated particle tracking on samples with sparsely labeled filaments, or from passive tracer beads. DLOF produces more accurate velocity fields than PIV for densely labeled samples. PIV cannot reliably distinguish contrast variations at high densities, particularly along the nematic director. DLOF overcomes this limitation. For sparsely labeled samples, DLOF and PIV produce comparable results, but DLOF gives higher-resolution fields. Our work establishes DLOF as a versatile tool for measuring fluid flows in a broad class of active, soft, and biophysical systems.

Phases and correlations in active nematic granular systems

We investigate the statistical behavior of a system comprising fore-aft symmetric rods confined between two vertically vibrating plates using numerical simulations closely resembling the experimental setup studied by V. Narayan et al., Science, 2007, 317, 105-108. Our focus lies in studying the phase transition of the system from an isotropic phase to a nematic phase as we increase either the rod length or the rod concentration. Our simulations confirm the presence of long-range ordering in the ordered phase. Furthermore, we identify a phase characterized by a periodic ordering profile that disrupts translation symmetry. We also provide a detailed analysis of the translational and rotational diffusive dynamics of the rods. Interestingly, the translational diffusivity of the rods is found to increase with the rod concentration.

Morphological control of bundled actin networks subject to fixed-mass depletion

Depletion interactions are thought to significantly contribute to the organization of intracellular structures in the crowded cytosol. The strength of depletion interactions depends on physical parameters such as the depletant number density and the depletant size ratio. Cells are known to dynamically regulate these two parameters by varying the copy number of proteins of a wide distribution of sizes. However, mammalian cells are also known to keep the total protein mass density remarkably constant, to within 0.5% throughout the cell cycle. We thus ask how the strength of depletion interactions varies when the total depletant mass is held fixed, a.k.a. fixed-mass depletion. We answer this question via scaling arguments, as well as by studying depletion effects on networks of reconstituted semiflexible actin in silico and in vitro. We examine the maximum strength of the depletion interaction potential U∗ as a function of q, the size ratio between the depletant and the matter being depleted. We uncover a scaling relation U∗ ∼ qζ for two cases: fixed volume fraction φ and fixed mass density ρ. For fixed volume fraction, we report ζ < 0. For the fixed mass density case, we report ζ > 0, which suggests that the depletion interaction strength increases as the depletant size ratio is increased. To test this prediction, we prepared our filament networks at fixed mass concentrations with varying sizes of the depletant molecule poly(ethylene glycol) (PEG). We characterize the depletion interaction strength in our simulations via the mesh size. In experiments, we observe two distinct actin network morphologies, which we call weakly bundled and strongly bundled. We identify a mass concentration where different PEG depletant sizes lead to weakly bundled or strongly bundled morphologies. For these conditions, we find that the mesh size and intra-bundle spacing between filaments across the different morphologies do not show significant differences, while the dynamic light scattering relaxation time and storage modulus between the two states do show significant differences. Our results demonstrate the ability to tune actin network morphology and mechanics by controlling depletant size and give insights into depletion interaction mechanisms under the fixed-depletant-mass constraint relevant to living cells.

Motor-driven microtubule diffusion in a photobleached dynamical coordinate system

Motor-driven cytoskeletal remodeling in cellular systems can often be accompanied by a diffusive-like effect at local scales, but distinguishing the contributions of the ordering process, such as active contraction of a network, from this active diffusion is difficult to achieve. Using light-dimerizable kinesin motors to spatially control the formation and contraction of a microtubule network, we deliberately photobleach a grid pattern onto the filament network serving as a transient and dynamic coordinate system to observe the deformation and translation of the remaining fluorescent squares of microtubules. We find that the network contracts at a rate set by motor speed but is accompanied by a diffusive-like spread throughout the bulk of the contracting network with effective diffusion constant two orders of magnitude lower than that for a freely-diffusing microtubule. We further find that on micron scales, the diffusive timescale is only a factor of ≈ 3 slower than that of advection regardless of conditions, showing that the global contraction and long-time relaxation from this diffusive behavior are both motor-driven but exhibit local competition within the network bulk.

The Casimir-like effect induced by active nematics

We consider an active nematic phase and use hydrodynamical equations of it to model the activity as an internal field. The interaction of this field with the nematic director in a confined geometry is included in the Hamiltonian of the system. Based on this model Hamiltonian and the standard field theoretical approach, the Casimir-like force induced between the boundaries of such a confined film is discussed. The force depends on the geometrical shape and the dynamical character of the constituents of our active phase, as well as the anchoring conditions. For homeotropically aligned rod-like particles which in principle tend to align along a planar flow field, extensile activity enhances the attraction present in a thin nematic film. As the film thickness increases the force reduces. Beyond a critical thickness, a planar flow field instantaneous to a bend distortion sets in. Near but below the threshold of this activity-induced instability, the force crosses zero and repulsively diverges right at the critical threshold of this so-called flow instability. For contractile rods, in the same geometry as above, the structure is stable and the Casimir-like force diminishes by an exponential factor as a function of the film thickness. On the other side for a planar director alignment, rod-like contractile particles can induce opposite shear flows at the boundaries creating a splay distortion for the director between the plates. In this configuration, we obtain the same universal pretransitional behavior for the force as above. Vice versa, for extensile rod-like particles in this geometry, the director fluctuations become massive and the Casimir-like force diminishes again by an exponential factor as the film thickness increases. The effect of the active field on thermal fluctuations of the director and the fluctuation-induced Casimir force per area is derived through a "semi"-dynamical approach as well. However, the results of the calculation due to a mathematical sum over the fluctuating modes do not lead to an approved closed form.

Aspect-ratio-dependent void formation in active rhomboidal and elliptical particle systems

We execute a numerical simulation of active nematics with particles interacting by an excluded-volume effect. Systems with rhomboidal particles and with elliptical particles are considered in order to investigate the effect of the direct contact of particles. In our simulation, the void regions, where the local number density is almost zero, appear in both systems when the aspect ratio of the particles is high. We focus on the relationship between the void regions and the particle orientation of the bulk. The particle number density, particle orientation, topological defects, and void regions are analyzed for different aspect ratios in both systems. The systems with rhomboidal particles have characteristic void sizes, which increase with an increase in the aspect ratio. In contrast, the distribution of the void-region size in the systems with elliptical particles is broad. The present results suggest that the void size in the systems with rhomboidal particles is determined by the correlation length of the particle orientational field around the void regions, while the void size might be determined by the system size in the systems with elliptical particles.

Computational tools for cellular scale biophysics

Mathematical models are indispensable for disentangling the interactions through which biological components work together to generate the forces and flows that position, mix, and distribute proteins, nutrients, and organelles within the cell. To illuminate the ever more specific questions studied at the edge of biological inquiry, such models inevitably become more complex. Solving, simulating, and learning from these more realistic models requires the development of new analytic techniques, numerical methods, and scalable software. In this review, we discuss some recent developments in tools for understanding how large numbers of cytoskeletal filaments, driven by molecular motors and interacting with the cytoplasm and other structures in their environment, generate fluid flows, instabilities, and material deformations which help drive crucial cellular processes.

Structure and dynamics of motor-driven microtubule bundles

Connecting the large-scale emergent behaviors of active cytoskeletal materials to the microscopic properties of their constituents is a challenge due to a lack of data on the multiscale dynamics and structure of such systems. We approach this problem by studying the impact of depletion attraction on bundles of microtubules and kinesin-14 molecular motors. For all depletant concentrations, kinesin-14 bundles generate comparable extensile dynamics. However, this invariable mesoscopic behavior masks the transition in the microscopic motion of microtubules. Specifically, with increasing attraction, we observe a transition from bi-directional sliding with extension to pure extension with no sliding. Small-angle X-ray scattering shows that the transition in microtubule dynamics is concurrent with a structural rearrangement of microtubules from an open hexagonal to a compressed rectangular lattice. These results demonstrate that bundles of microtubules and molecular motors can display the same mesoscopic extensile behaviors despite having different internal structures and microscopic dynamics. They provide essential information for developing multiscale models of active matter.

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.

Interplay of chemotactic force, Péclet number, and dimensionality dictates the dynamics of auto-chemotactic chiral active droplets

In living and synthetic active matter systems, the constituents can self-propel and interact with each other and with the environment through various physicochemical mechanisms. Among these mechanisms, chemotactic and auto-chemotactic effects are widely observed. The impact of (auto-)chemotactic effects on achiral active matter has been a recent research focus. However, the influence of these effects on chiral active matter remains elusive. Here, we develop a Brownian dynamics model coupled with a diffusion equation to examine the dynamics of auto-chemotactic chiral active droplets in both quasi-two-dimensional (2D) and three-dimensional (3D) systems. By quantifying the droplet trajectory as a function of the dimensionless Péclet number and chemotactic strength, our simulations well reproduce the curling and helical trajectories of nematic droplets in a surfactant-rich solution reported by Krüger et al. [Phys. Rev. Lett. 117, 048003 (2016)]. The modeled curling trajectory in 2D exhibits an emergent chirality, also consistent with the experiment. We further show that the geometry of the chiral droplet trajectories, characterized by the pitch and diameter, can be used to infer the velocities of the droplet. Interestingly, we find that, unlike the achiral case, the velocities of chiral active droplets show dimensionality dependence: its mean instantaneous velocity is higher in 3D than in 2D, whereas its mean migration velocity is lower in 3D than in 2D. Taken together, our particle-based simulations provide new insights into the dynamics of auto-chemotactic chiral active droplets, reveal the effects of dimensionality, and pave the way toward their applications, such as drug delivery, sensors, and micro-reactors.

Orientational order and topological defects in a dilute solutions of rodlike polymers at low Reynolds number

The relationship between the polymer orientation and the chaotic flow, in a dilute solution of rigid rodlike polymers at low Reynolds number, is investigated by means of direct numerical simulations. It is found that the rods tend to align with the velocity field in order to minimize the friction with the solvent fluid, while regions of rotational disorder are related to strong vorticity gradients, and therefore to the chaotic flow. The "turbulent-like" behavior of the system is therefore associated with the emergence and interaction of topological defects of the mean director field, similarly to active nematic turbulence. The analysis has been carried out in both two and three spatial dimensions.

Self-Timed and Spatially Targeted Delivery of Chemical Cargo by Motile Liquid Crystal

A key challenge underlying the design of miniature machines is encoding materials with time- and space-specific functional behaviors that require little human intervention. Dissipative processes that drive materials beyond equilibrium and evolve continuously with time and location represent one promising strategy to achieve such complex functions. This work reports how internal nonequilibrium states of liquid crystal (LC) emulsion droplets undergoing chemotaxis can be used to time the delivery of a chemical agent to a targeted location. During ballistic motion, hydrodynamic shear forces dominate LC elastic interactions, dispersing microdroplet inclusions (microcargo) within double emulsion droplets. Scale-dependent colloidal forces then hinder the escape of dispersed microcargo from the propelling droplet. Upon arrival at the targeted location, a circulatory flow of diminished strength allows the microcargo to cluster within the LC elastic environment such that hydrodynamic forces grow to exceed colloidal forces and thus trigger the escape of the microcargo. This work illustrates the utility of the approach by using microcargo that initiate polymerization upon release through the outer interface of the carrier droplet. These findings provide a platform that utilizes nonequilibrium strategies to design autonomous spatial and temporal functions into active materials.

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.

Percolation of nonequilibrium assemblies of colloidal particles in active chiral liquids

The growing interest in the non-equilibrium assembly of colloidal particles in active liquids is driven by the motivation to create novel structures endowed with tunable properties unattainable within the confines of equilibrium systems. Here, we present an experimental investigation of the structural features of colloidal assemblies in active liquids of chiral E. coli. The colloidal particles form dynamic clusters due to the effective interaction mediated by active media. The activity and chirality of the swimmers strongly influence the dynamics and local ordering of colloidal particles, resulting in clusters with persistent rotation, whose structure differs significantly from those in equilibrium systems with attractive interactions, such as colloid-polymer mixtures. Our colloid-bacteria mixture displays several hallmark features of a percolation transition at a critical density, where the clusters span the system size. A closer examination of the critical exponents associated with cluster size distribution, the average cluster size, and the correlation length in the vicinity of the critical density shows deviations from the prediction of the standard continuum percolation model. Therefore, our experiments reveal a richer phase behavior of colloidal assemblies in active liquids.

Microswimmers Knead Nematics into Cholesterics

The hydrodynamic stresses created by active particles can destabilize orientational order present in the system. This is manifested, for example, by the appearance of a bend instability in active nematics or in quasi-two-dimensional living liquid crystals consisting of swimming bacteria in thin nematic films. Using large-scale hydrodynamics simulations, we study a system consisting of spherical microswimmers within a three-dimensional nematic liquid crystal. We observe a spontaneous chiral symmetry breaking, where the uniform nematic state is kneaded into a continuously twisting state, corresponding to a helical director configuration akin to a cholesteric liquid crystal. The transition arises from the hydrodynamic coupling between the liquid crystalline elasticity and the swimmer flow fields, leading to a twist-bend instability of the nematic order. It is observed for both pusher (extensile) and puller (contractile) swimmers. Further, we show that the liquid crystal director and particle trajectories are connected: in the cholesteric state the particle trajectories become helicoidal.

Ordering kinetics in the active Ising model

We undertake a numerical study of the ordering kinetics in the two-dimensional (2D) active Ising model (AIM), a discrete flocking model with a conserved density field coupled to a nonconserved magnetization field. We find that for a quench into the liquid-gas coexistence region and in the ordered liquid region, the characteristic length scale of both the density and magnetization domains follows the Lifshitz-Cahn-Allen growth law, R(t)∼t^{1/2}, consistent with the growth law of passive systems with scalar order parameter and nonconserved dynamics. The system morphology is analyzed with the two-point correlation function and its Fourier transform, the structure factor, which conforms to the well-known Porod's law, a manifestation of the coarsening of compact domains with smooth boundaries. We also find the domain growth exponent unaffected by different noise strengths and self-propulsion velocities of the active particles. However, transverse diffusion is found to play the most significant role in the growth kinetics of the AIM. We extract the same growth exponent by solving the hydrodynamic equations of the AIM.

Active nematic ratchet in asymmetric obstacle arrays

We numerically investigate the effect of an asymmetric periodic obstacle array in a two-dimensional active nematic. We find that activity in conjunction with the asymmetry leads to a ratchet effect or unidirectional flow of the fluid along the asymmetry direction. The directional flow is still present even in the active turbulent phase when the gap between obstacles is sufficiently small. We demonstrate that the dynamics of the topological defects transition from flow mirroring to smectic-like as the gap between obstacles is made smaller, and explain this transition in terms of the pinning of negative winding number defects between obstacles. This also leads to a nonmonotonic ratchet effect magnitude as a function of obstacle size, so that there is an optimal obstacle size for ratcheting at fixed activity.

Periodic orbits, pair nucleation, and unbinding of active nematic defects on cones

Geometric confinement and topological constraints present promising means of controlling active materials. By combining analytical arguments derived from the Born-Oppenheimer approximation with numerical simulations, we investigate the simultaneous impact of confinement together with curvature singularity by characterizing the dynamics of an active nematic on a cone. Here, the Born-Oppenheimer approximation means that textures can follow defect positions rapidly on the timescales of interest. Upon imposing strong anchoring boundary conditions at the base of a cone, we find a rich phase diagram of multidefect dynamics, including exotic periodic orbits of one or two +1/2 flank defects, depending on activity and nonquantized geometric charge at the cone apex. By characterizing the transitions between these ordered dynamical states, we present detailed understanding of (i) defect unbinding, (ii) defect absorption, and (iii) defect pair nucleation at the apex. Numerical simulations confirm theoretical predictions of not only the nature of the circular orbits but also defect unbinding from the apex.

Reconstitution of cytolinker-mediated crosstalk between actin and vimentin

Cell shape and motility are determined by the cytoskeleton, an interpenetrating network of actin filaments, microtubules, and intermediate filaments. The biophysical properties of each filament type individually have been studied extensively by cell-free reconstitution. By contrast, the interactions between the three cytoskeletal networks are relatively unexplored. They are coupled via crosslinkers of the plakin family such as plectin. These are challenging proteins for reconstitution because of their giant size and multidomain structure. Here we engineer a recombinant actin-vimentin crosslinker protein called 'ACTIF' that provides a minimal model system for plectin, recapitulating its modular design with actin-binding and intermediate filament-binding domains separated by a coiled-coil linker for dimerisation. We show by fluorescence and electron microscopy that ACTIF has a high binding affinity for vimentin and actin and creates mixed actin-vimentin bundles. Rheology measurements show that ACTIF-mediated crosslinking strongly stiffens actin-vimentin composites. Finally, we demonstrate the modularity of this approach by creating an ACTIF variant with the intermediate filament binding domain of Adenomatous Polyposis Coli. Our protein engineering approach provides a new cell-free system for the biophysical characterization of intermediate filament-binding crosslinkers and for understanding the mechanical synergy between actin and vimentin in mesenchymal cells.

Phase separation kinetics and cluster dynamics in two-dimensional active dumbbell systems

Molecular dynamics simulations were employed to investigate the phase separation process of a two-dimensional active Brownian dumbbell model. We evaluated the time dependence of the typical size of the dense component using the scaling properties of the structure factor, along with the averaged number of clusters and their radii of gyration. The growth observed is faster than in active disk models, and this effect is further enhanced under stronger activity. Next, we focused on studying the hexatic order of the clusters. The length associated with the orientational order increases algebraically with time and faster than for spherical active Brownian particles. Under weak active forces, most clusters exhibit a uniform internal orientational order. However, under strong forces, large clusters consist of domains with different orientational orders. We demonstrated that the latter configurations are not stable, and given sufficient time to evolve, they eventually achieve homogeneous configurations as well. No gas bubbles are formed within the clusters, even when there are patches of different hexatic order. Finally, attention was directed towards the geometry and motion of the clusters themselves. By employing a tracking algorithm, we showed that clusters smaller than the typical size at the observation time exhibit regular shapes, while larger ones display fractal characteristics. In between collisions or break-ups, the clusters behave as solid bodies. Their centers of mass undergo circular motion, with radii increasing with the cluster size. The angular velocity of the center of mass equals that of the constituents with respect to their center of mass. These observations were rationalised with a simple mechanical model.

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.

Programming tunable active dynamics in a self-propelled robot

We present a scheme for producing tunable active dynamics in a self-propelled robotic device. The robot moves using the differential drive mechanism where two wheels can vary their instantaneous velocities independently. These velocities are calculated by equating robot's equations of motion in two dimensions with well-established active particle models and encoded into the robot's microcontroller. We demonstrate that the robot can depict active Brownian, run and tumble, and Brownian dynamics with a wide range of parameters. The resulting motion analyzed using particle tracking shows excellent agreement with the theoretically predicted trajectories. Later, we show that its motion can be switched between different dynamics using light intensity as an external parameter. Intriguingly, we demonstrate that the robot can efficiently navigate through many obstacles by performing stochastic reorientations driven by the gradient in light intensity towards a desired location, namely the target. This work opens an avenue for designing tunable active systems with the potential of revealing the physics of active matter and its application for bio- and nature-inspired robotics.

Design rules for controlling active topological defects

Topological defects play a central role in the physics of many materials, including magnets, superconductors, and liquid crystals. In active fluids, defects become autonomous particles that spontaneously propel from internal active stresses and drive chaotic flows stirring the fluid. The intimate connection between defect textures and active flow suggests that properties of active materials can be engineered by controlling defects, but design principles for their spatiotemporal control remain elusive. Here, we propose a symmetry-based additive strategy for using elementary activity patterns, as active topological tweezers, to create, move, and braid such defects. By combining theory and simulations, we demonstrate how, at the collective level, spatial activity gradients act like electric fields which, when strong enough, induce an inverted topological polarization of defects, akin to a negative susceptibility dielectric. We harness this feature in a dynamic setting to collectively pattern and transport interacting active defects. Our work establishes an additive framework to sculpt flows and manipulate active defects in both space and time, paving the way to design programmable active and living materials for transport, memory, and logic.

Stress and alignment response to curved obstacles in growing bacterial monolayers

Monolayers of growing bacteria, confined within channel geometries, exhibit self-organization into a highly aligned laminar state along the axis of the channel. Although this phenomenon has been observed in experiments and simulations under various boundary conditions, the underlying physical mechanism driving this alignment remains unclear. In this study, we conduct simulations of growing bacteria in two-dimensional channel geometries perturbed by fixed obstacles, either circular or arc shaped, placed at the channel's center. Our findings reveal that even sizable obstacles cause only short-ranged disruptions to the baseline laminar state. These disruptions arise from a competition between local planar anchoring and bulk laminar alignment. At smaller obstacle sizes, bulk alignment fully dominates, while at larger sizes planar anchoring induces increasing local disruptions. Furthermore, our analysis indicates that the resulting configurations of the bacterial system display a striking resemblance to the arrangement of hard-rod smectic liquid crystals around circular obstacles. This suggests that modeling hard-rod bacterial monolayers as smectic, rather than nematic, liquid crystals may yield successful outcomes. The insights gained from our study contribute to the expanding body of research on bacterial growth in channels. Our work provides perspectives on the stability of the laminar state and extends our understanding to encompass more intricate confinement schemes.

Weak Ion-Exchange Based Magnetic Swarm for Targeted Drug Delivery and Chemotherapy

Swimming microrobots that are actuated by multiple stimuli/fields display various intriguing collective behaviors, ranging from phase separation to clustering and giant number fluctuation; however, it is still chanllenging to achieve multiple responses and functionalities within one colloidal system to emulate high environmental adaptability and improved tasking capability of natural swarms. In this work, a weak ion-exchange based swarm is presented that can self-organize and reconfigure by chemical, light, and magnetic fields, showing living crystal, amorphous glass, liquid, chain, and wheel-like structures. By changing the frequency and strength of the rotating magnetic field, various well-controlled and fast transformations are obtained. Experiments show the high adaptability and functionality of the microrobot swarm in delivering drugs in confined spaces, such as narrow channels with turns or obstacles. The drug-carrying swarm exhibits excellent chemtherapy for Hela and CT26 cells due to the pH-enhanced drug release and locomotion. This reconfigurable microswarm provides a new platform for biomedical and environmental applications.

Non-equilibrium structural and dynamic behaviors of active polymers in complex and crowded environments

Active matter systems, which convert internal chemical energy or energy from the environment into directed motion, are ubiquitous in nature and exhibit a range of emerging non-equilibrium behaviors. However, most of the current works on active matter have been devoted to particles, and the study of active polymers has only recently come into the spotlight due to their prevalence within living organisms. The intricate interplay between activity and conformational degrees of freedom gives rise to novel structural and dynamical behaviors of active polymers. Research in active polymers remarkably broadens diverse concepts of polymer physics, such as molecular architecture, dynamics, scaling and so on, which is of significant importance for the development of new polymer materials with unique performance. Furthermore, active polymers are often found in strongly interacting and crowded systems and in complex environments, so that the understanding of this behavior is essential for future developments of novel polymer-based biomaterials. This review thereby focuses on the study of active polymers in complex and crowded environments, and aims to provide insights into the fundamental physics underlying the adaptive and collective behaviors far from equilibrium, as well as the open challenges that the field is currently facing.

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.

Machine learning topological defects in confluent tissues

Active nematics is an emerging paradigm for characterizing biological systems. One aspect of particularly intense focus is the role active nematic defects play in these systems, as they have been found to mediate a growing number of biological processes. Accurately detecting and classifying these defects in biological systems is, therefore, of vital importance to improving our understanding of such processes. While robust methods for defect detection exist for systems of elongated constituents, other systems, such as epithelial layers, are not well suited to such methods. Here, we address this problem by developing a convolutional neural network to detect and classify nematic defects in confluent cell layers. Crucially, our method is readily implementable on experimental images of cell layers and is specifically designed to be suitable for cells that are not rod shaped, which we demonstrate by detecting defects on experimental data using the trained model. We show that our machine learning model outperforms current defect detection techniques and that this manifests itself in our method as requiring less data to accurately capture defect properties. This could drastically improve the accuracy of experimental data interpretation while also reducing costs, advancing the study of nematic defects in biological systems.

Motor crosslinking augments elasticity in active nematics

In active materials, uncoordinated internal stresses lead to emergent long-range flows. An understanding of how the behavior of active materials depends on mesoscopic (hydrodynamic) parameters is developing, but there remains a gap in knowledge concerning how hydrodynamic parameters depend on the properties of microscopic elements. In this work, we combine experiments and multiscale modeling to relate the structure and dynamics of active nematics composed of biopolymer filaments and molecular motors to their microscopic properties, in particular motor processivity, speed, and valency. We show that crosslinking of filaments by both motors and passive crosslinkers not only augments the contributions to nematic elasticity from excluded volume effects but dominates them. By altering motor kinetics we show that a competition between motor speed and crosslinking results in a nonmonotonic dependence of nematic flow on motor speed. By modulating passive filament crosslinking we show that energy transfer into nematic flow is in large part dictated by crosslinking. Thus motor proteins both generate activity and contribute to nematic elasticity. Our results provide new insights for rationally engineering active materials.

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.

Coherent light scattering from cellular dynamics in living tissues

This review examines the biological physics of intracellular transport probed by the coherent optics of dynamic light scattering from optically thick living tissues. Cells and their constituents are in constant motion, composed of a broad range of speeds spanning many orders of magnitude that reflect the wide array of functions and mechanisms that maintain cellular health. From the organelle scale of tens of nanometers and upward in size, the motion inside living tissue is actively driven rather than thermal, propelled by the hydrolysis of bioenergetic molecules and the forces of molecular motors. Active transport can mimic the random walks of thermal Brownian motion, but mean-squared displacements are far from thermal equilibrium and can display anomalous diffusion through Lévy or fractional Brownian walks. Despite the average isotropic three-dimensional environment of cells and tissues, active cellular or intracellular transport of single light-scattering objects is often pseudo-one-dimensional, for instance as organelle displacement persists along cytoskeletal tracks or as membranes displace along the normal to cell surfaces, albeit isotropically oriented in three dimensions. Coherent light scattering is a natural tool to characterize such tissue dynamics because persistent directed transport induces Doppler shifts in the scattered light. The many frequency-shifted partial waves from the complex and dynamic media interfere to produce dynamic speckle that reveals tissue-scale processes through speckle contrast imaging and fluctuation spectroscopy. Low-coherence interferometry, dynamic optical coherence tomography, diffusing-wave spectroscopy, diffuse-correlation spectroscopy, differential dynamic microscopy and digital holography offer coherent detection methods that shed light on intracellular processes. In health-care applications, altered states of cellular health and disease display altered cellular motions that imprint on the statistical fluctuations of the scattered light. For instance, the efficacy of medical therapeutics can be monitored by measuring the changes they induce in the Doppler spectra of livingex vivocancer biopsies.

Controlling active turbulence by activity patterns

By patterning activity in space, one can control active turbulence. To show this, we use Doi's hydrodynamic equations of a semidilute solution of active rods. A linear stability analysis reveals the resting isotropic fluid to be unstable above an absolute pusher activity. The emergent activity-induced paranematic state displays active turbulence, which we characterize by different quantities including the energy spectrum, which shows the typical power-law decay with exponent -4. Then, we control the active turbulence by a square lattice of circular spots where activity is switched off. In the parameter space lattice constant versus surface-to-surface distance of the spots, we identify different flow states. Most interestingly, for lattice constants below the vorticity correlation length and for spot distances smaller than the nematic coherence length, we observe a multi-lane flow state, where flow lanes with alternating flow directions are separated by a street of vortices. The flow pattern displays pronounced multistability and also appears transiently at the transition to the isotropic active-turbulence state. At larger lattice constants a trapped vortex state is identified with a non-Gaussian vorticity distribution due to the low flow vorticity at the spots. It transitions to conventional active turbulence for increasing spot distance.