Topology and dynamics of active nematic vesicles

Hydrodynamic Enhancement of p-atic Defect Dynamics

We investigate numerically and analytically the effects of hydrodynamics on the dynamics of topological defects in p-atic liquid crystals, i.e., two-dimensional liquid crystals with p-fold rotational symmetry. Importantly, we find that hydrodynamics fuels a generic passive self-propulsion mechanism for defects of winding number s=(p-1)/p and arbitrary p. Strikingly, we discover that hydrodynamics always accelerates the annihilation dynamics of pairs of ±1/p defects and that, contrary to expectations, this effect increases with p. Our Letter paves the way toward understanding cell intercalation and other remodeling events in epithelial layers.

Curvature in Biological Systems: Its Quantification, Emergence, and Implications across the Scales

Surface curvature both emerges from, and influences the behavior of, living objects at length scales ranging from cell membranes to single cells to tissues and organs. The relevance of surface curvature in biology is supported by numerous experimental and theoretical investigations in recent years. In this review, first, a brief introduction to the key ideas of surface curvature in the context of biological systems is given and the challenges that arise when measuring surface curvature are discussed. Giving an overview of the emergence of curvature in biological systems, its significance at different length scales becomes apparent. On the other hand, summarizing current findings also shows that both single cells and entire cell sheets, tissues or organisms respond to curvature by modulating their shape and their migration behavior. Finally, the interplay between the distribution of morphogens or micro-organisms and the emergence of curvature across length scales is addressed with examples demonstrating these key mechanistic principles of morphogenesis. Overall, this review highlights that curved interfaces are not merely a passive by-product of the chemical, biological, and mechanical processes but that curvature acts also as a signal that co-determines these processes.

The crucial role of adhesion in the transmigration of active droplets through interstitial orifices

Active fluid droplets are a class of soft materials exhibiting autonomous motion sustained by an energy supply. Such systems have been shown to capture motility regimes typical of biological cells and are ideal candidates as building-block for the fabrication of soft biomimetic materials of interest in pharmacology, tissue engineering and lab on chip devices. While their behavior is well established in unconstrained environments, much less is known about their dynamics under strong confinement. Here, we numerically study the physics of a droplet of active polar fluid migrating within a microchannel hosting a constriction with adhesive properties, and report evidence of a striking variety of dynamic regimes and morphological features, whose properties crucially depend upon droplet speed and elasticity, degree of confinement within the constriction and adhesiveness to the pore. Our results suggest that non-uniform adhesion forces are instrumental in enabling the crossing through narrow orifices, in contrast to larger gaps where a careful balance between speed and elasticity is sufficient to guarantee the transition. These observations may be useful for improving the design of artificial micro-swimmers, of interest in material science and pharmaceutics, and potentially for cell sorting in microfluidic devices.

Polarity and chirality control of an active fluid by passive nematic defects

Much like passive materials, active systems can be affected by the presence of imperfections in their microscopic order, called defects, that influence macroscopic properties. This suggests the possibility to steer collective patterns by introducing and controlling defects in an active system. Here we show that a self-assembled, passive nematic is ideally suited to control the pattern formation process of an active fluid. To this end, we force microtubules to glide inside a passive nematic material made from actin filaments. The actin nematic features self-assembled half-integer defects that steer the active microtubules and lead to the formation of macroscopic polar patterns. Moreover, by confining the nematic in circular geometries, chiral loops form. We find that the exact positioning of nematic defects in the passive material deterministically controls the formation and the polarity of the active flow, opening the possibility of efficiently shaping an active material using passive defects.

Regularization of dielectric tensor tomography

Dielectric tensor tomography reconstructs the three-dimensional dielectric tensors of microscopic objects and provides information about the crystalline structure orientations and principal refractive indices. Because dielectric tensor tomography is based on transmission measurement, it suffers from the missing cone problem, which causes poor axial resolution, underestimation of the refractive index, and halo artifacts. In this study, we study the application of total variation and positive semi-definiteness regularization to three-dimensional tensor distributions. In particular, we demonstrate the reduction of artifacts when applied to dielectric tensor tomography.

DisGUVery: A Versatile Open-Source Software for High-Throughput Image Analysis of Giant Unilamellar Vesicles

Giant unilamellar vesicles (GUVs) are cell-sized aqueous compartments enclosed by a phospholipid bilayer. Due to their cell-mimicking properties, GUVs have become a widespread experimental tool in synthetic biology to study membrane properties and cellular processes. In stark contrast to the experimental progress, quantitative analysis of GUV microscopy images has received much less attention. Currently, most analysis is performed either manually or with custom-made scripts, which makes analysis time-consuming and results difficult to compare across studies. To make quantitative GUV analysis accessible and fast, we present DisGUVery, an open-source, versatile software that encapsulates multiple algorithms for automated detection and analysis of GUVs in microscopy images. With a performance analysis, we demonstrate that DisGUVery's three vesicle detection modules successfully identify GUVs in images obtained with a wide range of imaging sources, in various typical GUV experiments. Multiple predefined analysis modules allow the user to extract properties such as membrane fluorescence, vesicle shape, and internal fluorescence from large populations. A new membrane segmentation algorithm facilitates spatial fluorescence analysis of nonspherical vesicles. Altogether, DisGUVery provides an accessible tool to enable high-throughput automated analysis of GUVs, and thereby to promote quantitative data analysis in synthetic cell research.

Active morphogenesis of patterned epithelial shells

Shape transformations of epithelial tissues in three dimensions, which are crucial for embryonic development or in vitro organoid growth, can result from active forces generated within the cytoskeleton of the epithelial cells. How the interplay of local differential tensions with tissue geometry and with external forces results in tissue-scale morphogenesis remains an open question. Here, we describe epithelial sheets as active viscoelastic surfaces and study their deformation under patterned internal tensions and bending moments. In addition to isotropic effects, we take into account nematic alignment in the plane of the tissue, which gives rise to shape-dependent, anisotropic active tensions and bending moments. We present phase diagrams of the mechanical equilibrium shapes of pre-patterned closed shells and explore their dynamical deformations. Our results show that a combination of nematic alignment and gradients in internal tensions and bending moments is sufficient to reproduce basic building blocks of epithelial morphogenesis, including fold formation, budding, neck formation, flattening, and tubulation.

Topological defect-mediated morphodynamics of active-active interfaces

Physical interfaces widely exist in nature and engineering. Although the formation of passive interfaces is well elucidated, the physical principles governing active interfaces remain largely unknown. Here, we combine simulation, theory, and cell-based experiment to investigate the evolution of an active-active interface. We adopt a biphasic framework of active nematic liquid crystals. We find that long-lived topological defects mechanically energized by activity display unanticipated dynamics nearby the interface, where defects perform "U-turns" to keep away from the interface, push the interface to develop local fingers, or penetrate the interface to enter the opposite phase, driving interfacial morphogenesis and cross-interface defect transport. We identify that the emergent interfacial morphodynamics stems from the instability of the interface and is further driven by the activity-dependent defect-interface interactions. Experiments of interacting multicellular monolayers with extensile and contractile differences in cell activity have confirmed our predictions. These findings reveal a crucial role of topological defects in active-active interfaces during, for example, boundary formation and tissue competition that underlie organogenesis and clinically relevant disorders.

Active Bending of Disordered Microtubule Bundles by Kinesin Motors

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

Motility of acoustically powered micro-swimmers in a liquid crystalline environment

Suspensions of microswimmers in liquid crystals demonstrate remarkably complex dynamics and serve as a model system for studying active nematics. So far, experimental realization of microswimmers suspended in liquid crystalline media has relied on biological microorganisms that impose strict limitations on the compatible media and makes it difficult to regulate activity. Here, we demonstrate that acoustically powered bubble microswimmers can efficiently self-propel in a lyotropic liquid crystal. The velocity of the swimmers is controlled by the amplitude of the acoustic field. Histograms of swimming directions with respect to the local nematic field reveal a bimodal distribution: the swimmers tend to either fully align with or swim perpendicular to the director field of the liquid crystal, occasionally switching between these two states. The bubble-induced streaming from a swimmer locally melts the liquid crystal and produces topological defects at the tail of the swimmer. We show that the defect proliferation rate increases with the angle between the swimmer's velocity and the local orientation of the director field.

Defect dynamics in active polar fluids vs. active nematics

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

Collective dynamics and shattering of disturbed two-dimensional Lennard-Jones crystals

Elucidating collective dynamics in crystalline systems is a common scientific question in multiple fields. In this work, by combination of high-precision numerical approach and analytical normal mode analysis, we systematically investigate the dynamical response of two-dimensional Lennard-Jones crystal as a purely classical mechanical system under random disturbance of varying strength, and reveal rich microscopic dynamics. Specifically, we observe highly symmetric velocity field composed of sharply divided coherent and disordered regions, and identify the order-disorder dynamical transition of the velocity field. Under stronger disturbance, we reveal the vacancy-driven shattering of the crystal. This featured disruption mode is fundamentally different from the dislocation-unbinding scenario in two-dimensional melting. We also examine the healing dynamics associated with vacancies of varying size. The results in this work advance our understanding about the formation of collective dynamics and crystal disruption, and may have implications in elucidating relevant non-equilibrium behaviors in a host of crystalline systems. Microscopic dynamics and underlying topological defects in the disruption of the Lennard-Jones lattice under random disturbance.

Active nematic defects in compressible and incompressible flows

We study the dynamics of active nematic films on a substrate driven by active flows with or without the incompressible constraint. Through simulations and theoretical analysis, we show that arch patterns are stable in the compressible case, while they become unstable under the incompressibility constraint. For compressible flows at high enough activity, stable arches organize themselves into a smecticlike pattern, which induce an associated global polar ordering of +1/2 nematic defects. By contrast, divergence-free flows give rise to a local nematic order of the +1/2 defects, consisting of antialigned pairs of neighboring defects, as established in previous studies.

Tubulation and dispersion of oil by bacterial growth on droplets

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

Non-equilibrium shapes and dynamics of active vesicles

Active vesicles, constructed through the confinement of self-propelled particles (SPPs) inside a lipid membrane shell, exhibit a large variety of non-equilibrium shapes, ranging from the formation of local tethers and dendritic conformations, to prolate and bola-like structures. To better understand the behavior of active vesicles, we perform simulations of membranes modelled as dynamically triangulated surfaces enclosing active Brownian particles. A systematic analysis of membrane deformations and SPP clustering, as a function of SPP activity and volume fraction inside the vesicle is carried out. Distributions of membrane local curvature, and the clustering and mobility of SPPs obtained from simulations of active vesicles are analysed. There exists a feedback mechanism between the enhancement of membrane curvature, the formation of clusters of active particles, and local or global changes in vesicle shape. The emergence of active tension due to the activity of SPPs can well be captured by the Young-Laplace equation. Furthermore, a simple numerical method for tether detection is presented and used to determine correlations between the number of tethers, their length, and local curvature. We also provide several geometrical arguments to explain different tether characteristics for various conditions. These results contribute to the future development of steerable active vesicles or soft micro-robots whose behaviour can be controlled and used for potential applications.

Monodomain Liquid Crystals of Two-Dimensional Sheets by Boundary-Free Sheargraphy

Eliminating topological defects to achieve monodomain liquid crystals is highly significant for the fundamental studies of soft matter and building long-range ordered materials. However, liquid crystals are metastable and sensitive to external stimuli, such as flow, confinement, and electromagnetic fields, which cause their intrinsic polycrystallinity and topological defects. Here, we achieve the monodomain liquid crystals of graphene oxide over 30 cm through boundary-free sheargraphy. The obtained monodomain liquid crystals exhibit large-area uniform alignment of sheets, which has the same optical polarized angle and intensity. The monodomain liquid crystals provide bidirectionally ordered skeletons, which can be applied as lightweight thermal management materials with bidirectionally high thermal and electrical conductivity. Furthermore, we extend the controllable topology of two-dimensional colloids by introducing singularities and disclinations in monodomain liquid crystals. Topological structures with defect strength from - 2 to + 2 were realized. This work provides a facile methodology to study the structural order of soft matter at a macroscopic level, facilitating the fabrication of metamaterials with tunable and highly anisotropic architectures.

Response of vesicle shapes to dense inner active matter

We consider experimentally the Takatori-Sahu model of vesicle shape fluctuations induced by enclosed active matter, a model till present tested only in the absence of collective motion because few enclosed bacteria were used to generate the desired active motion (S. C. Takatori and A. Sahu, Phys. Rev. Lett., 2020, 124, 158102). Using deformable giant unilamellar vesicles (GUVs) and phase contrast microscopy, we extract the mode-dependence of GUV shape fluctuations when hundreds of E. coli bacteria are contained within each GUV. In the microscope focal plane, patterns of collective bacteria flow include vortex flow, dipolar flow, and chaotic motion, all of which influence the GUV shapes. The Takatori-Sahu model generalizes well to this situation if one considers the moving element to be the experimentally-determined size of the collecively-moving flock.

Emergent collective dynamics of pusher and puller squirmer rods: swarming, clustering, and turbulence

We study the interplay of steric and hydrodynamic interactions in suspensions of elongated microswimmers by simulating the full hydrodynamics of squirmer rods in the quasi two-dimensional geometry of a Hele-Shaw cell. To create pusher or puller-type squirmer rods, we concentrate the surface slip-velocity field more to the back or to the front of the rod and thereby are able to tune the rod's force-dipole strength. We study a wide range of aspect ratios and area fractions and provide corresponding state diagrams. The flow field of pusher-type squirmer rods destabilizes ordered structures and favors the disordered state at small area fractions and aspect ratios. Only when steric interactions become relevant, we observe a turbulent and dynamic cluster state, while for large aspect ratios a single swarm and jammed cluster occurs. The power spectrum of the turbulent state shows two distinct energy cascades at small and large wave numbers with power-law scaling and non-universal exponents. Pullers show a strong tendency to form swarms instead of the disordered state found for neutral and pusher rods. At large area fractions a dynamic cluster is observed and at larger aspect ratio a single swarm or jammed cluster occurs.

Defect absorption and emission for p-atic liquid crystals on cones

We investigate the ground-state configurations of two-dimensional liquid crystals with p-fold rotational symmetry (p-atics) on fixed curved surfaces. We focus on the intrinsic geometry and show that isothermal coordinates are particularly convenient as they explicitly encode a geometric contribution to the elastic potential. In the special case of a cone with half-angle β, the apex develops an effective topological charge of -χ, where 2πχ=2π(1-sinβ) is the deficit angle of the cone, and a topological defect of charge σ behaves as if it had an effective topological charge Q_{eff}=(σ-σ^{2}/2) when interacting with the apex. The effective charge of the apex leads to defect absorption and emission at the cone apex as the deficit angle of the cone is varied. For total topological defect charge 1, e.g., imposed by tangential boundary conditions at the edge, we find that for a disk the ground-state configuration consists of p defects each of charge +1/p lying equally spaced on a concentric ring of radius d=(p-1/3p-1)^{1/2p}R, where R is the radius of the disk. In the case of a cone with tangential boundary conditions at the base, we find three types of ground-state configurations as a function of cone angle: (i) for sharp cones, all of the +1/p defects are absorbed by the apex; (ii) at intermediate cone angles, some of the +1/p defects are absorbed by the apex and the rest lie equally spaced along a concentric ring on the flank; and (iii) for nearly flat cones, all of the +1/p defects lie equally spaced along a concentric ring on the flank. Here the defect positions and the absorption transitions depend intricately on p and the deficit angle, which we analytically compute. We check these results with numerical simulations for a set of commensurate cone angles and find excellent agreement.

Biotropic liquid crystal phase transformations in cellulose-producing bacterial communities

Biological functionality is often enabled by a fascinating variety of physical phenomena that emerge from orientational order of building blocks, a defining property of nematic liquid crystals that is also pervasive in nature. Out-of-equilibrium, "living" analogs of these technological materials are found in biological embodiments ranging from myelin sheath of neurons to extracellular matrices of bacterial biofilms and cuticles of beetles. However, physical underpinnings behind manifestations of orientational order in biological systems often remain unexplored. For example, while nematiclike birefringent domains of biofilms are found in many bacterial systems, the physics behind their formation is rarely known. Here, using cellulose-synthesizing Acetobacter xylinum bacteria, we reveal how biological activity leads to orientational ordering in fluid and gel analogs of these soft matter systems, both in water and on solid agar, with a topological defect found between the domains. Furthermore, the nutrient feeding direction plays a role like that of rubbing of confining surfaces in conventional liquid crystals, turning polydomain organization within the biofilms into a birefringent monocrystal-like order of both the extracellular matrix and the rod-like bacteria within it. We probe evolution of scalar orientational order parameters of cellulose nanofibers and bacteria associated with fluid-gel and isotropic-nematic transformations, showing how highly ordered active nematic fluids and gels evolve with time during biological-activity-driven, disorder-order transformation. With fluid and soft-gel nematics observed in a certain range of biological activity, this mesophase-exhibiting system is dubbed "biotropic," analogously to thermotropic nematics that exhibit solely orientational order within a temperature range, promising technological and fundamental-science applications.

Bacterial active matter

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

Pattern Propagation Driven by Surface Curvature

Pattern dynamics on curved surfaces are found everywhere in nature. The geometry of surfaces has been shown to influence dynamics and play a functional role, yet a comprehensive understanding is still elusive. Here, we report for the first time that a static Turing pattern on a flat surface can propagate on a curved surface, as opposed to previous studies, where the pattern is presupposed to be static irrespective of the surface geometry. To understand such significant changes on curved surfaces, we investigate reaction-diffusion systems on axisymmetric curved surfaces. Numerical and theoretical analyses reveal that both the symmetries of the surface and pattern participate in the initiation of pattern propagation. This study provides a novel and generic mechanism of pattern propagation that is caused by surface curvature, as well as insights into the general role of surface geometry.

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

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

Topological fine structure of smectic grain boundaries and tetratic disclination lines within three-dimensional smectic liquid crystals

Observing and characterizing the complex ordering phenomena of liquid crystals subjected to external constraints constitutes an ongoing challenge for chemists and physicists alike. To elucidate the delicate balance appearing when the intrinsic positional order of smectic liquid crystals comes into play, we perform Monte-Carlo simulations of rod-like particles in a range of cavities with a cylindrical symmetry. Based on recent insights into the topology of smectic orientational grain boundaries in two dimensions, we analyze the emerging three-dimensional defect structures from the perspective of tetratic symmetry. Using an appropriate three-dimensional tetratic order parameter constructed from the Steinhardt order parameters, we show that those grain boundaries can be interpreted as a pair of tetratic disclination lines that are located on the edges of the nematic domain boundary. Thereby, we shed light on the fine structure of grain boundaries in three-dimensional confined smectics.

Activity-induced polar patterns of filaments gliding on a sphere

Active matter systems feature the ability to form collective patterns as observed in a plethora of living systems, from schools of fish to swimming bacteria. While many of these systems move in a wide, three-dimensional environment, several biological systems are confined by a curved topology. The role played by a non-Euclidean geometry on the self-organization of active systems is not yet fully understood, and few experimental systems are available to study it. Here, we introduce an experimental setup in which actin filaments glide on the inner surface of a spherical lipid vesicle, thus embedding them in a curved geometry. We show that filaments self-assemble into polar, elongated structures and that, when these match the size of the spherical geometry, both confinement and topological constraints become relevant for the emergent patterns, leading to the formation of polar vortices and jammed states. These results experimentally demonstrate that activity-induced complex patterns can be shaped by spherical confinement and topology.

Integer topological defects organize stresses driving tissue morphogenesis

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

Fractional defect charges in liquid crystals with p-fold rotational symmetry on cones

Conical surfaces, with a δ function of Gaussian curvature at the apex, are perhaps the simplest example of geometric frustration. We study two-dimensional liquid crystals with p-fold rotational symmetry (p-atics) on the surfaces of cones. For free boundary conditions at the base, we find both the ground state(s) and a discrete ladder of metastable states as a function of both the cone angle and the liquid crystal symmetry p. We find that these states are characterized by a set of fractional defect charges at the apex and that the ground states are in general frustrated due to effects of parallel transport along the azimuthal direction of the cone. We check our predictions for the ground-state energies numerically for a set of commensurate cone angles (corresponding to a set of commensurate Gaussian curvatures concentrated at the cone apex), whose surfaces can be polygonized as a perfect triangular or square mesh, and find excellent agreement with our theoretical predictions.

Assembling Microtubule-Based Active Matter

Studied for more than a century, equilibrium liquid crystals provided insight into the properties of ordered materials, and led to commonplace applications such as display technology. Active nematics are a new class of liquid crystal materials that are driven out of equilibrium by continuous motion of the constituent anisotropic units. A versatile experimental realization of active nematic liquid crystals is based on rod-like cytoskeletal filaments that are driven out of equilibrium by molecular motors. We describe protocols for assembling microtubule-kinesin based active nematic liquid crystals and associated isotropic fluids. We describe the purification of each protein and the assembly process of a two-dimensional active nematic on a water-oil interface. Finally, we show examples of nematic formation and describe methods for quantifying their non-equilibrium dynamics.

Theory of defect-mediated morphogenesis

Growing experimental evidence indicates that topological defects could serve as organizing centers in the morphogenesis of tissues. Here, we provide a quantitative explanation for this phenomenon, rooted in the buckling theory of deformable active polar liquid crystals. Using a combination of linear stability analysis and computational fluid dynamics, we demonstrate that active layers, such as confined cell monolayers, are unstable to the formation of protrusions in the presence of disclinations. The instability originates from an interplay between the focusing of the elastic forces, mediated by defects, and the renormalization of the system's surface tension by the active flow. The posttransitional regime is also characterized by several complex morphodynamical processes, such as oscillatory deformations, droplet nucleation, and active turbulence. Our findings offer an explanation of recent observations on tissue morphogenesis and shed light on the dynamics of active surfaces in general.

Active elastocapillarity in soft solids with negative surface tension

Active solids consume energy to allow for actuation, shape change, and wave propagation not possible in equilibrium. Whereas active interfaces have been realized across many experimental systems, control of three-dimensional (3D) bulk materials remains a challenge. Here, we develop continuum theory and microscopic simulations that describe a 3D soft solid whose boundary experiences active surface stresses. The competition between active boundary and elastic bulk yields a broad range of previously unexplored phenomena, which are demonstrations of so-called active elastocapillarity. In contrast to thin shells and vesicles, we discover that bulk 3D elasticity controls snap-through transitions between different anisotropic shapes. These transitions meet at a critical point, allowing a universal classification via Landau theory. In addition, the active surface modifies elastic wave propagation to allow zero, or even negative, group velocities. These phenomena offer robust principles for programming shape change and functionality into active solids, from robotic metamaterials down to shape-shifting nanoparticles.

Engineering stability, longevity, and miscibility of microtubule-based active fluids

Microtubule-based active matter provides insight into the self-organization of motile interacting constituents. We describe several formulations of microtubule-based 3D active isotropic fluids. Dynamics of these fluids is powered by three types of kinesin motors: a processive motor, a non-processive motor, and a motor which is permanently linked to a microtubule backbone. Another modification uses a specific microtubule crosslinker to induce bundle formation instead of a non-specific polymer depletant. In comparison to the already established system, each formulation exhibits distinct properties. These developments reveal the temporal stability of microtubule-based active fluids while extending their reach and the applicability.

Logic operations with active topological defects

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

Shape multistability in flexible tubular crystals through interactions of mobile dislocations

Crystalline sheets rolled up into cylinders occur in diverse biological and synthetic systems, including carbon nanotubes, biofilaments of the cellular cytoskeleton, and packings of colloidal particles. In this work, we show, computationally, that such tubular crystals can be programmed with reconfigurable shapes, due to motions of defects that interrupt the periodicity of the crystalline lattice. By identifying and exploiting stable patterns of these defects, we cause tubular crystals to relax into desired target geometries, a design principle that could guide the creation of versatile colloidal analogues to nanotubes. Our results suggest routes to tunable and switchable material properties in ordered, soft materials on deformable surfaces.

We study avenues to shape multistability and shape morphing in flexible crystalline membranes of cylindrical topology, enabled by glide mobility of dislocations. Using computational modeling, we obtain states of mechanical equilibrium presenting a wide variety of tubular crystal deformation geometries, due to an interplay of effective defect interactions with out-of-tangent-plane deformations that reorient the tube axis. Importantly, this interplay often stabilizes defect configurations quite distinct from those predicted for a two-dimensional crystal confined to the surface of a rigid cylinder. We find that relative and absolute stability of competing states depend strongly on control parameters such as bending rigidity, applied stress, and spontaneous curvature. Using stable dislocation pair arrangements as building blocks, we demonstrate that targeted macroscopic three-dimensional conformations of thin crystalline tubes can be programmed by imposing certain sparse patterns of defects. Our findings reveal a broad design space for controllable and reconfigurable colloidal tube geometries, with potential relevance also to architected carbon nanotubes and microtubules.

Exact Coherent Structures and Phase Space Geometry of Preturbulent 2D Active Nematic Channel Flow

Confined active nematics exhibit rich dynamical behavior, including spontaneous flows, periodic defect dynamics, and chaotic "active turbulence." Here, we study these phenomena using the framework of exact coherent structures, which has been successful in characterizing the routes to high Reynolds number turbulence of passive fluids. Exact coherent structures are stationary, periodic, quasiperiodic, or traveling wave solutions of the hydrodynamic equations that, together with their invariant manifolds, serve as an organizing template of the dynamics. We compute the dominant exact coherent structures and connecting orbits in a preturbulent active nematic channel flow, which enables a fully nonlinear but highly reduced-order description in terms of a directed graph. Using this reduced representation, we compute instantaneous perturbations that switch the system between disparate spatiotemporal states occupying distant regions of the infinite-dimensional phase space. Our results lay the groundwork for a systematic means of understanding and controlling active nematic flows in the moderate- to high-activity regime.

Self-Organized Liquid Crystal Droplets as Phototunable Softmasks

A single-step self-organized pathway is harnessed to generate large-area and high-density liquid-crystal (LC) microdroplets via rapid spreading of an LC-laden volatile liquid film on an aqueous surfactant bath. The surfactant loading on the water bath and LC loading in the solvent fluid help in tuning the size, periodicity, and ordering of LC microdroplets. Remarkably, the experiments reveal a transition from a spinodal to heterogeneous nucleation pathway of dewetting when the surfactant loading is modulated from below to beyond the critical micellar concentration in the aqueous phase. In the process, a host of unprecedented drop formation modes, such as dewetting and contact-line instability, random ejection, and "fire cracker" toroid splitting, have been uncovered. Subsequently, the LC microdroplets on the air-water interface are employed as photomasks suitable for soft-photolithography applications. Such masks help in the decoration of a host of mesoscale three-dimensional features on the films of photoresists when photons are guided through the LC droplets. In such a scenario, phase transition of LC droplets under solvent vapor annealing is employed to control the movement of photons through drops and subsequently modulate the light exposure on the photoresist surface. Such a simple soft-photolithography setup leads to an array of flattened droplets on a positive resist, while donut features are observed on the negative tone. Remarkably, the orientation of nematogens within 4-cyano-4'-pentylbiphenyl droplets and at the three-phase contact-line provides additional handles in controlling the transmission of photons, which facilitates such a unique pattern formation. A number of low-cost and simple strategies are also discussed to order such soft-photolithography patterns. Importantly, with a minor modification to the same experimental setup, we could also measure the variation in the order parameter of the LC droplet during its phase transitions from the nematic to isotropic state.

Vesicle shape transformations driven by confined active filaments

In active matter systems, deformable boundaries provide a mechanism to organize internal active stresses. To study a minimal model of such a system, we perform particle-based simulations of an elastic vesicle containing a collection of polar active filaments. The interplay between the active stress organization due to interparticle interactions and that due to the deformability of the confinement leads to a variety of filament spatiotemporal organizations that have not been observed in bulk systems or under rigid confinement, including highly-aligned rings and caps. In turn, these filament assemblies drive dramatic and tunable transformations of the vesicle shape and its dynamics. We present simple scaling models that reveal the mechanisms underlying these emergent behaviors and yield design principles for engineering active materials with targeted shape dynamics.

Pattern formation, localized and running pulsation on active spherical membranes

Active force generation by an actin-myosin cortex coupled to a cell membrane allows the cell to deform, respond to the environment, and mediate cell motility and division. Several membrane-bound activator proteins move along it and couple to the membrane curvature. Besides, they can act as nucleating sites for the growth of filamentous actin. Actin polymerization can generate a local outward push on the membrane. Inward pull from the contractile actomyosin cortex can propagate along the membrane via actin filaments. We use coupled evolution of fields to perform linear stability analysis and numerical calculations. As activity overcomes the stabilizing factors such as surface tension and bending rigidity, the spherical membrane shows instability towards pattern formation, localized pulsation, and running pulsation between poles. We present our results in terms of phase diagrams and evolutions of the coupled fields. They have relevance for living cells and can be verified in experiments on artificial cell-like constructs.

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.

Orientational Correlations in Active and Passive Nematic Defects

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

Metamachines of pluripotent colloids

Machines enabled the Industrial Revolution and are central to modern technological progress: A machine's parts transmit forces, motion, and energy to one another in a predetermined manner. Today's engineering frontier, building artificial micromachines that emulate the biological machinery of living organisms, requires faithful assembly and energy consumption at the microscale. Here, we demonstrate the programmable assembly of active particles into autonomous metamachines using optical templates. Metamachines, or machines made of machines, are stable, mobile and autonomous architectures, whose dynamics stems from the geometry. We use the interplay between anisotropic force generation of the active colloids with the control of their orientation by local geometry. This allows autonomous reprogramming of active particles of the metamachines to achieve multiple functions. It permits the modular assembly of metamachines by fusion, reconfiguration of metamachines and, we anticipate, a shift in focus of self-assembly towards active matter and reprogrammable materials.

Designing architectural materials: from granular form to functional granular material

Designed granular materials are a novel class of architectural material system. Following one of the key paradigms of designed matter, material form and material function are closely interrelated in these systems. In this context, the article aims to contribute a parametric particle design model as an interface for this interrelation. A granular material is understood as an aggregation of large numbers of individual particles between which only short-range repulsive contact forces are acting. Granular materials are highly pertinent material systems for architecture. Due to the fact that they can act both as a solid and a liquid, they can be recycled and reconfigured multiple times and are thus highly sustainable. Designed granular materials have the added potential that the function of the granular material can be calibrated through the definition of the particles' form. Research on the design of granular materials in architecture is nascent. In physics they have been explored mainly with respect to different particle shapes. However, no coherent parametric particle design model of designed particle shapes for granular material systems in architecture has yet been established which considers both fabrication constraints and simulation requirements. The parametric particle design model proposed in this article has been based on a design system which has been developed through feasibility tests and simulations conducted in research and teaching. Based on this design system the parametric particle design model is developed integrating both fabrication constraints for architecture-scale particle systems and the geometric requirements of established simulation methods for granular materials. Initially the design system and related feasibility tests are presented. The parametric particle design model resulting from that is then described in detail. Directions of further research are discussed especially with respect to the integration of the parametric particle design model in 'inverse' design methods.

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

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

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.