Geometrical control of interface patterning underlies active matter invasion.

Interaction between active materials and the boundaries of geometrical confinement is key to many emergent phenomena in active systems. For living active matter consisting of animal cells or motile bacteria, the confinement boundary is often a deformable interface, and it has been unclear how activity-induced interface dynamics might lead to morphogenesis and pattern formation. Here, we studied the evolution of bacterial active matter confined by a deformable boundary. We found that an ordered morphological pattern emerged at the interface characterized by periodically spaced interfacial protrusions; behind the interfacial protrusions, bacterial swimmers self-organized into multicellular clusters displaying +1/2 nematic defects. Subsequently, a hierarchical sequence of transitions from interfacial protrusions to creeping branches allowed the bacterial active drop to rapidly invade surrounding space with a striking self-similar branch pattern. We found that this interface patterning is geometrically controlled by the local curvature of the interface, a phenomenon we denote as collective curvature sensing. Using a continuum active model, we revealed that the collective curvature sensing arises from enhanced active stresses near high-curvature regions, with the active length scale setting the characteristic distance between the interfacial protrusions. Our findings reveal a protrusion-to-branch transition as a unique mode of active matter invasion and suggest a strategy to engineer pattern formation of active materials.

Cell sorting by active forces in a phase-field model of cell monolayers.

Cell sorting, the segregation of cells with different properties into distinct domains, is a key phenomenon in biological processes such as embryogenesis. We use a phase-field model of a confluent cell layer to study the role of activity in cell sorting. We find that a mixture of cells with extensile or contractile dipolar activity, and which are identical apart from their activity, quickly sort into small, elongated patches which then grow slowly in time. We interpret the sorting as driven by the different diffusivity of the extensile and contractile cells, mirroring the ordering of Brownian particles connected to different hot and cold thermostats. We check that the free energy is not changed by either partial or complete sorting, thus confirming that activity can be responsible for the ordering even in the absence of thermodynamic mechanisms.

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

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

Active Forces in Confluent Cell Monolayers.

We use a computational phase-field model together with analytical analysis to study how intercellular active forces can mediate individual cell morphology and collective motion in a confluent cell monolayer. We explore the regime where intercellular forces dominate the tissue dynamics, and polar forces are negligible. Contractile intercellular interactions lead to cell elongation, nematic ordering, and active turbulence characterized by motile topological defects. Extensile interactions result in frustration, and perpendicular cell orientations become more prevalent. Furthermore, we show that contractile behavior can change to extensile behavior if anisotropic fluctuations in cell shape are considered.

Phase Separation Driven by Active Flows.

We extend the continuum theories of active nematohydrodynamics to model a two-fluid mixture with separate velocity fields for each fluid component, coupled through a viscous drag. The model is used to study an active nematic fluid mixed with an isotropic fluid. We find microphase separation, and argue that this results from an interplay between active anchoring and active flows driven by concentration gradients. The results may be relevant to cell sorting and the formation of lipid rafts in cell membranes.

Shape-Tension Coupling Produces Nematic Order in an Epithelium Vertex Model.

We study the vertex model for epithelial tissue mechanics extended to include coupling between the cell shapes and tensions in cell-cell junctions. This coupling represents an active force which drives the system out of equilibrium and leads to the formation of nematic order interspersed with prominent, long-lived +1 defects. The defects in the nematic ordering are coupled to the shape of the cell tiling, affecting cell areas and coordinations. This intricate interplay between cell shape, size, and coordination provides a possible mechanism by which tissues could spontaneously develop long-range polarity through local mechanical forces without resorting to long-range chemical patterning.

Activity-driven tissue alignment in proliferating spheroids.

We extend the continuum theory of active nematic fluids to study cell flows and tissue dynamics inside multicellular spheroids, spherical, self-assembled aggregates of cells that are widely used as model systems to study tumour dynamics. Cells near the surface of spheroids have better access to nutrients and therefore proliferate more rapidly than those in the resource-depleted core. Using both analytical arguments and three-dimensional simulations, we find that the proliferation gradients result in flows and in gradients of activity both of which can align the orientation axis of cells inside the aggregates. Depending on environmental conditions and the intrinsic tissue properties, we identify three distinct alignment regimes: spheroids in which all the cells align either radially or tangentially to the surface throughout the aggregate and spheroids with angular cell orientation close to the surface and radial alignment in the core. The continuum description of tissue dynamics inside spheroids not only allows us to infer dynamic cell parameters from experimentally measured cell alignment profiles, but more generally motivates novel mechanisms for controlling the alignment of cells within aggregates which has been shown to influence the mechanical properties and invasive capabilities of tumors.

Active nematics with deformable particles.

The hydrodynamic theory of active nematics has been often used to describe the spatio-temporal dynamics of cell flows and motile topological defects within soft confluent tissues. Those theories, however, often rely on the assumption that tissues consist of cells with a fixed, anisotropic shape and do not resolve dynamical cell shape changes due to flow gradients. In this paper we extend the continuum theory of active nematics to include cell shape deformability. We find that circular cells in tissues must generate sufficient active stress to overcome an elastic barrier to deforming their shape in order to drive tissue-scale flows. Above this threshold the systems enter a dynamical steady-state with regions of elongated cells and strong flows coexisting with quiescent regions of isotropic cells.

Steering self-organisation through confinement.

Self-organisation is the spontaneous emergence of spatio-temporal structures and patterns from the interaction of smaller individual units. Examples are found across many scales in very different systems and scientific disciplines, from physics, materials science and robotics to biology, geophysics and astronomy. Recent research has highlighted how self-organisation can be both mediated and controlled by confinement. Confinement is an action over a system that limits its units' translational and rotational degrees of freedom, thus also influencing the system's phase space probability density; it can function as either a catalyst or inhibitor of self-organisation. Confinement can then become a means to actively steer the emergence or suppression of collective phenomena in space and time. Here, to provide a common framework and perspective for future research, we examine the role of confinement in the self-organisation of soft-matter systems and identify overarching scientific challenges that need to be addressed to harness its full scientific and technological potential in soft matter and related fields. By drawing analogies with other disciplines, this framework will accelerate a common deeper understanding of self-organisation and trigger the development of innovative strategies to steer it using confinement, with impact on, e.g., the design of smarter materials, tissue engineering for biomedicine and in guiding active matter.

Topological defects in epithelia govern cell death and extrusion.

Epithelial tissues (epithelia) remove excess cells through extrusion, preventing the accumulation of unnecessary or pathological cells. The extrusion process can be triggered by apoptotic signalling, oncogenic transformation and overcrowding of cells. Despite the important linkage of cell extrusion to developmental, homeostatic and pathological processes such as cancer metastasis, its underlying mechanism and connections to the intrinsic mechanics of the epithelium are largely unexplored. We approach this problem by modelling the epithelium as an active nematic liquid crystal (that has a long range directional order), and comparing numerical simulations to strain rate and stress measurements within monolayers of MDCK (Madin Darby canine kidney) cells. Here we show that apoptotic cell extrusion is provoked by singularities in cell alignments in the form of comet-shaped topological defects. We find a universal correlation between extrusion sites and positions of nematic defects in the cell orientation field in different epithelium types. The results confirm the active nematic nature of epithelia, and demonstrate that defect-induced isotropic stresses are the primary precursors of mechanotransductive responses in cells, including YAP (Yes-associated protein) transcription factor activity, caspase-3-mediated cell death, and extrusions. Importantly, the defect-driven extrusion mechanism depends on intercellular junctions, because the weakening of cell-cell interactions in an α-catenin knockdown monolayer reduces the defect size and increases both the number of defects and extrusion rates, as is also predicted by our model. We further demonstrate the ability to control extrusion hotspots by geometrically inducing defects through microcontact printing of patterned monolayers. On the basis of these results, we propose a mechanism for apoptotic cell extrusion: spontaneously formed topological defects in epithelia govern cell fate. This will be important in predicting extrusion hotspots and dynamics in vivo, with potential applications to tissue regeneration and the suppression of metastasis. Moreover, we anticipate that the analogy between the epithelium and active nematic liquid crystals will trigger further investigations of the link between cellular processes and the material properties of epithelia.

Investigating the nature of active forces in tissues reveals how contractile cells can form extensile monolayers.

Actomyosin machinery endows cells with contractility at a single-cell level. However, within a monolayer, cells can be contractile or extensile based on the direction of pushing or pulling forces exerted by their neighbours or on the substrate. It has been shown that a monolayer of fibroblasts behaves as a contractile system while epithelial or neural progentior monolayers behave as an extensile system. Through a combination of cell culture experiments and in silico modelling, we reveal the mechanism behind this switch in extensile to contractile as the weakening of intercellular contacts. This switch promotes the build-up of tension at the cell-substrate interface through an increase in actin stress fibres and traction forces. This is accompanied by mechanotransductive changes in vinculin and YAP activation. We further show that contractile and extensile differences in cell activity sort cells in mixtures, uncovering a generic mechanism for pattern formation during cell competition, and morphogenesis.

Active Extensile Stress Promotes 3D Director Orientations and Flows.

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

Activity gradients in two- and three-dimensional active nematics.

We numerically investigate how spatial variations of extensile or contractile active stress affect bulk active nematic systems in two and three dimensions. In the absence of defects, activity gradients drive flows which re-orient the nematic director field and thus act as an effective anchoring force. At high activity, defects are created and the system transitions into active turbulence, a chaotic flow state characterized by strong vorticity. We find that in two-dimensional (2D) systems active torques robustly align +1/2 defects parallel to activity gradients, with defect heads pointing towards contractile regions. In three-dimensional (3D) active nematics disclination lines preferentially lie in the plane perpendicular to activity gradients due to active torques acting on line segments. The average orientation of the defect structures in the plane perpendicular to the line tangent depends on the defect type, where wedge-like +1/2 defects align parallel to activity gradients, while twist defects are aligned anti-parallel. Understanding the response of active nematic fluids to activity gradients is an important step towards applying physical theories to biology, where spatial variations of active stress impact morphogenetic processes in developing embryos and affect flows and deformations in growing cell aggregates, such as tumours.

Coupling Turing stripes to active flows.

We numerically solve the active nematohydrodynamic equations of motion, coupled to a Turing reaction-diffusion model, to study the effect of active nematic flow on the stripe patterns resulting from a Turing instability. If the activity is uniform across the system, the Turing patterns dissociate when the flux from active advection balances that from the reaction-diffusion process. If the activity is coupled to the concentration of Turing morphogens, and neighbouring stripes have equal and opposite activity, the system self organises into a pattern of shearing flows, with stripes tending to fracture and slip sideways to join their neighbours. We discuss the role of active instabilities in controlling the crossover between these limits. Our results are of relevance to mechanochemical coupling in biological systems.

Memory effects, arches and polar defect ordering at the cross-over from wet to dry active nematics.

We use analytic arguments and numerical solutions of the continuum, active nematohydrodynamic equations to study how friction alters the behaviour of active nematics. Concentrating on the case where there is nematic ordering in the passive limit, we show that, as the friction is increased, memory effects become more prominent and +1/2 topological defects leave increasingly persistent trails in the director field as they pass. The trails are preferential sites for defect formation and they tend to impose polar order on any new +1/2 defects. In the absence of noise and for high friction, it becomes very difficult to create defects, but trails formed by any defects present at the beginning of the simulations persist and organise into parallel arch-like patterns in the director field. We show aligned arches of equal width are approximate steady state solutions of the equations of motion which co-exist with the nematic state. We compare our results to other models in the literature, in particular dry systems with no hydrodynamics, where trails, arches and polar defect ordering have also been observed.

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.

Flow States and Transitions of an Active Nematic in a Three-Dimensional Channel.

We use active nematohydrodynamics to study the flow of an active fluid in a 3D microchannel, finding a transition between active turbulence and regimes where there is a net flow along the channel. We show that the net flow is only possible if the active nematic is flow aligning and that, in agreement with experiments, the appearance of the net flow depends on the aspect ratio of the channel cross section. We explain our results in terms of when the hydrodynamic screening due to the channel walls allows the emergence of vortex rolls across the channel.

Role of Friction in Multidefect Ordering.

We use continuum simulations to study the impact of friction on the ordering of defects in an active nematic. Even in a frictionless system, +1/2 defects tend to align side by side and orient antiparallel reflecting their propensity to form, and circulate with, flow vortices. Increasing friction enhances the effectiveness of the defect-defect interactions, and defects form dynamically evolving, large-scale, positionally, and orientationally ordered structures, which can be explained as a competition between hexagonal packing, preferred by the -1/2 defects, and rectangular packing, preferred by the +1/2 defects.

Degenerate states, emergent dynamics and fluid mixing by magnetic rotors.

We investigate the collective motion of magnetic rotors suspended in a viscous fluid under a uniform rotating magnetic field. The rotors are positioned on a square lattice, and low Reynolds hydrodynamics is assumed. For a 3 × 3 array of magnets, we observe three characteristic dynamical patterns as the external field strength is varied: a synchronized pattern, an oscillating pattern, and a chessboard pattern. The relative stability of these depends on the competition between the energy due to the external magnetic field and the energy of the magnetic dipole-dipole interactions among the rotors. We argue that the chessboard pattern can be understood as an alternation in the stability of two degenerate states, characterized by striped and spin-ice configurations, as the applied magnetic field rotates. For larger arrays, we observe propagation of slip waves that are similar to metachronal waves. The rotor arrays have potential as microfluidic devices that can mix fluids and create vortices of different sizes.

Active nematics with anisotropic friction: the decisive role of the flow aligning parameter.

We use continuum simulations to study the impact of anisotropic hydrodynamic friction on the emergent flows of active nematics. We show that, depending on whether the active particles align with or tumble in their collectively self-induced flows, anisotropic friction can result in markedly different patterns of motion. In a flow-aligning regime and at high anisotropic friction, the otherwise chaotic flows are streamlined into flow lanes with alternating directions, reproducing the experimental laning state that has been obtained by interfacing microtubule-motor protein mixtures with smectic liquid crystals. Within a flow-tumbling regime, however, we find that no such laning state is possible. Instead, the synergistic effects of friction anisotropy and flow tumbling can lead to the emergence of bound pairs of topological defects that align at an angle to the easy flow direction and navigate together throughout the domain. In addition to confirming the mechanism behind the laning states observed in experiments, our findings emphasise the role of the flow aligning parameter in the dynamics of active nematics.