Designing highly efficient interlocking interactions in anisotropic active particles.

Cluster formation of microscopic swimmers is key to the formation of biofilms and colonies, efficient motion and nutrient uptake, but, in the absence of other interactions, requires high swimmer concentrations to occur. Here we experimentally and numerically show that cluster formation can be dramatically enhanced by an anisotropic swimmer shape. We analyze a class of model microswimmers with a shape that can be continuously tuned from spherical to bent and straight rods. In all cases, clustering can be described by Michaelis-Menten kinetics governed by a single scaling parameter that depends on particle density and shape only. We rationalize these shape-dependent dynamics from the interplay between interlocking probability and cluster stability. The bent rod shape promotes assembly in an interlocking fashion even at vanishingly low particle densities and we identify the most efficient shape to be a semicircle. Our work provides key insights into how shape can be used to rationally design out-of-equilibrium self-organization, key to creating active functional materials and processes that require two-component assembly with high fidelity.

Lipid membranes supported by polydimethylsiloxane substrates with designed geometry.

The membrane curvature of cells and intracellular compartments continuously adapts to enable cells to perform vital functions, from cell division to signal trafficking. Understanding how membrane geometry affects these processes in vivo is challenging because of the biochemical and geometrical complexity as well as the short time and small length scales involved in cellular processes. By contrast, in vitro model membranes with engineered curvature would provide a versatile platform for this investigation and applications to biosensing and biocomputing. Here, we present a strategy that allows fabrication of lipid membranes with designed shape by combining 3D micro-printing and replica-molding lithography with polydimethylsiloxane to create curved micrometer-sized scaffolds with virtually any geometry. The resulting supported lipid membranes are homogeneous and fluid. We demonstrate the versatility of the system by fabricating structures of interesting combinations of mean and Gaussian curvature. We study the lateral phase separation and how local curvature influences the effective diffusion coefficient. Overall, we offer a bio-compatible platform for understanding curvature-dependent cellular processes and developing programmable bio-interfaces for living cells and nanostructures.

Hydrodynamics and multiscale order in confluent epithelia.

We formulate a hydrodynamic theory of confluent epithelia: i.e. monolayers of epithelial cells adhering to each other without gaps. Taking advantage of recent progresses toward establishing a general hydrodynamic theory of p-atic liquid crystals, we demonstrate that collectively migrating epithelia feature both nematic (i.e. p = 2) and hexatic (i.e. p = 6) orders, with the former being dominant at large and the latter at small length scales. Such a remarkable multiscale liquid crystal order leaves a distinct signature in the system's structure factor, which exhibits two different power-law scaling regimes, reflecting both the hexagonal geometry of small cells clusters and the uniaxial structure of the global cellular flow. We support these analytical predictions with two different cell-resolved models of epithelia - i.e. the self-propelled Voronoi model and the multiphase field model - and highlight how momentum dissipation and noise influence the range of fluctuations at small length scales, thereby affecting the degree of cooperativity between cells. Our construction provides a theoretical framework to conceptualize the recent observation of multiscale order in layers of Madin-Darby canine kidney cells and pave the way for further theoretical developments.

Hexanematic crossover in epithelial monolayers depends on cell adhesion and cell density.

Changes in tissue geometry during developmental processes are associated with collective migration of cells. Recent experimental and numerical results suggest that these changes could leverage on the coexistence of nematic and hexatic orientational order at different length scales. How this multiscale organization is affected by the material properties of the cells and their substrate is presently unknown. In this study, we address these questions in monolayers of Madin-Darby canine kidney cells having various cell densities and molecular repertoires. At small length scales, confluent monolayers are characterized by a prominent hexatic order, independent of the presence of E-cadherin, monolayer density, and underlying substrate stiffness. However, all three properties affect the meso-scale tissue organization. The length scale at which hexatic order transits to nematic order, the "hexanematic" crossover scale, strongly depends on cell-cell adhesions and correlates with monolayer density. Our study demonstrates how epithelial organization is affected by mechanical properties, and provides a robust description of tissue organization during developmental processes.

Tuneable defect-curvature coupling and topological transitions in active shells.

Recent experimental observations have suggested that topological defects can facilitate the creation of sharp features in developing embryos. Whereas these observations echo established knowledge about the interplay between geometry and topology in two-dimensional passive liquid crystals, the role of activity has mostly remained unexplored. In this article we focus on deformable shells consisting of either polar or nematic active liquid crystals and demonstrate that activity renders the mechanical coupling between defects and curvature much more involved and versatile than previously thought. Using a combination of linear stability analysis and three-dimensional computational fluid dynamics, we demonstrate that such a coupling can in fact be tuned, depending on the type of liquid crystal order, the specific structure of the defect (i.e. asters or vortices) and the nature of the active forces. In polar systems, this can drive a spectacular transition from spherical to toroidal topology, in the presence of large extensile activity. Our analysis strengthens the idea that defects could serve as topological morphogens and provides a number of predictions that could be tested in in vitro studies, for instance in the context of organoids.

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.

Hydrodynamic theory of p-atic liquid crystals.

We formulate a comprehensive hydrodynamic theory of two-dimensional liquid crystals with generic p-fold rotational symmetry, also known as p-atics, of which nematics (p=2) and hexatics (p=6) are the two best known examples. Previous hydrodynamic theories of p-atics are characterized by continuous O(2) rotational symmetry, which is higher than the discrete rotational symmetry of p-atic phases. By contrast, here we demonstrate that the discrete rotational symmetry allows the inclusion of additional terms in the hydrodynamic equations, which, in turn, lead to novel phenomena, such as the possibility of flow alignment at high shear rates, even for p>2. Furthermore, we show that any finite imposed shear will induce long-ranged orientational order in any p-atic liquid crystal, in contrast to the quasi-long-ranged order that occurs in the absence of shear. The induced order parameter scales like a nonuniversal power of the applied shear rate at small shear rates.

Long-Ranged Order and Flow Alignment in Sheared p-atic Liquid Crystals.

We formulate a hydrodynamic theory of p-atic liquid crystals, namely, two-dimensional anisotropic fluids endowed with generic p-fold rotational symmetry. Our approach, based on an order parameter tensor that directly embodies the discrete rotational symmetry of p-atic phases, allows us to unveil several unknown aspects of flowing p-atics, that previous theories, characterized by O(2) rotational symmetry, could not account for. This includes the onset of long-ranged orientational order in the presence of a simple shear flow of arbitrary shear rate, as opposed to the standard quasi-long-ranged order of two-dimensional liquid crystals, and the possibility of flow alignment at large shear rates.

Statistics and topology of fluctuating ribbons.

Ribbons are a class of slender structures whose length, width, and thickness are widely separated from each other. This scale separation gives a ribbon unusual mechanical properties in athermal macroscopic settings, for example, it can bend without twisting, but cannot twist without bending. Given the ubiquity of ribbon-like biopolymers in biology and chemistry, here we study the statistical mechanics of microscopic inextensible, fluctuating ribbons loaded by forces and torques. We show that these ribbons exhibit a range of topologically and geometrically complex morphologies exemplified by three phases-a twist-dominated helical phase (HT), a writhe-dominated helical phase (HW), and an entangled phase-that arise as the applied torque and force are varied. Furthermore, the transition from HW to HT phases is characterized by the spontaneous breaking of parity symmetry and the disappearance of perversions (that correspond to chirality-reversing localized defects). This leads to a universal response curve of a topological quantity, the link, as a function of the applied torque that is similar to magnetization curves in second-order phase transitions.

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.

Self-Assembly Dynamics of Reconfigurable Colloidal Molecules.

Colloidal molecules are designed to mimic their molecular analogues through their anisotropic shape and interactions. However, current experimental realizations are missing the structural flexibility present in real molecules thereby restricting their use as model systems. We overcome this limitation by assembling reconfigurable colloidal molecules from silica particles functionalized with mobile DNA linkers in high yields. We achieve this by steering the self-assembly pathway toward the formation of finite-sized clusters by employing high number ratios of particles functionalized with complementary DNA strands. The size ratio of the two species of particles provides control over the overall cluster size, i.e., the number of bound particles N, as well as the degree of reconfigurability. The bond flexibility provided by the mobile linkers allows the successful assembly of colloidal clusters with the geometrically expected maximum number of bound particles and shape. We quantitatively examine the self-assembly dynamics of these flexible colloidal molecules by a combination of experiments, agent-based simulations, and an analytical model. Our "flexible colloidal molecules" are exciting building blocks for investigating and exploiting the self-assembly of complex hierarchical structures, photonic crystals, and colloidal metamaterials.

Confinement-induced self-organization in growing bacterial colonies.

We investigate the emergence of global alignment in colonies of dividing rod-shaped cells under confinement. Using molecular dynamics simulations and continuous modeling, we demonstrate that geometrical anisotropies in the confining environment give rise to an imbalance in the normal stresses, which, in turn, drives a collective rearrangement of the cells. This behavior crucially relies on the colony's solid-like mechanical response at short time scales and can be recovered within the framework of active hydrodynamics upon modeling bacterial colonies as growing viscoelastic gels characterized by Maxwell-like stress relaxation.

Faceting and Flattening of Emulsion Droplets: A Mechanical Model.

When cooled down, emulsion droplets stabilized by a frozen interface of alkane molecules and surfactants have been observed to undergo a spectacular sequence of morphological transformations: from spheres to faceted liquid icosahedra, down to flattened liquid platelets. While generally ascribed to the interplay between the elasticity of the frozen interface and surface tension, the physical mechanisms underpinning these transitions have remained elusive, despite different theoretical pictures having been proposed in recent years. In this Letter, we introduce a comprehensive mechanical model of morphing emulsion droplets, which quantitatively accounts for various experimental observations, including the size scaling behavior of the faceting transition. Our analysis highlights the role of gravity and the spontaneous curvature of the frozen interface in determining the specific transition pathway.

Measuring Gaussian Rigidity Using Curved Substrates.

The Gaussian (saddle splay) rigidity of fluid membranes controls their equilibrium topology but is notoriously difficult to measure. In lipid mixtures, typical of living cells, linear interfaces separate liquid ordered (LO) from liquid disordered (LD) bilayer phases at subcritical temperatures. Here, we consider such membranes supported by curved substrates that thereby control the membrane curvatures. We show how spectral analysis of the fluctuations of the LO-LD interface provides a novel way of measuring the difference in Gaussian rigidity between the two phases. We provide a number of conditions for such interface fluctuations to be both experimentally measurable and sufficiently sensitive to the value of the Gaussian rigidity, while remaining in the perturbative regime of our analysis.

Geometric pinning and antimixing in scaffolded lipid vesicles.

Previous studies on the phase behaviour of multicomponent lipid bilayers found an intricate interplay between membrane geometry and its composition, but a fundamental understanding of curvature-induced effects remains elusive. Thanks to a combination of experiments on lipid vesicles supported by colloidal scaffolds and theoretical work, we demonstrate that the local geometry and global chemical composition of the bilayer determine both the spatial arrangement and the amount of mixing of the lipids. In the mixed phase, a strong geometrical anisotropy can give rise to an antimixed state, where the lipids are mixed, but their relative concentration varies across the membrane. After phase separation, the bilayer organizes in multiple lipid domains, whose location is pinned in specific regions, depending on the substrate curvature and the bending rigidity of the lipid domains. Our results provide critical insights into the phase separation of cellular membranes and, more generally, two-dimensional fluids on curved substrates.

Dislocation screening in crystals with spherical topology.

Whereas disclination defects are energetically prohibitive in two-dimensional flat crystals, their existence is necessary in crystals with spherical topology, such as viral capsids, colloidosomes, or fullerenes. Such a geometrical frustration gives rise to large elastic stresses, which render the crystal unstable when its size is significantly larger than the typical lattice spacing. Depending on the compliance of the crystal with respect to stretching and bending deformations, these stresses are alleviated either by a local increase of the intrinsic curvature in proximity of the disclinations or by the proliferation of excess dislocations, often organized in the form of one-dimensional chains known as "scars." The associated strain field of the scars is such as to counterbalance the one resulting from the isolated disclinations. Here we develop a continuum theory of dislocation screening in two-dimensional closed crystals with genus one. Upon modeling the flux of scars emanating from a given disclination as an independent scalar field, we demonstrate that the elastic energy of closed two-dimensional crystals with various degrees of asphericity can be expressed as a simple quadratic function of the screened topological charge of the disclinations, at both zero and finite temperature. This allows us to predict the optimal density of the excess dislocations as well as the minimal stretching energy attained by the crystal.

Mechanical interplay between cell shape and actin cytoskeleton organization.

We investigate the mechanical interplay between the spatial organization of the actin cytoskeleton and the shape of animal cells adhering on micropillar arrays. Using a combination of analytical work, computer simulations and in vitro experiments, we demonstrate that the orientation of the stress fibers strongly influences the geometry of the cell edge. In the presence of a uniformly aligned cytoskeleton, the cell edge can be well approximated by elliptical arcs, whose eccentricity reflects the degree of anisotropy of the cell's internal stresses. Upon modeling the actin cytoskeleton as a nematic liquid crystal, we further show that the geometry of the cell edge feeds back on the organization of the stress fibers by altering the length scale at which these are confined. This feedback mechanism is controlled by a dimensionless number, the anchoring number, representing the relative weight of surface-anchoring and bulk-aligning torques. Our model allows to predict both cellular shape and the internal structure of the actin cytoskeleton and is in good quantitative agreement with experiments on fibroblastoid (GDß1, GDß3) and epithelioid (GEß1, GEß3) cells.

Lipid exchange enhances geometric pinning in multicomponent membranes on patterned substrates.

Experiments on supported lipid bilayers featuring liquid ordered/disordered domains have shown that the spatial arrangement of the lipid domains and their chemical composition are strongly affected by the curvature of the substrate. Furthermore, theoretical predictions suggest that both these effects are intimately related with the closed topology of the bilayer. In this work, we test this hypothesis by fabricating supported membranes consisting of colloidal particles of various shapes lying on a flat substrate. A single lipid bilayer coats both colloids and substrate, allowing local lipid exchange between them, thus rendering the system thermodynamically open, i.e. able to exchange heat and molecules with an external reservoir in the neighborhood of the colloid. By reconstructing the Gibbs phase diagram for this system, we demonstrate that the free-energy landscape is directly influenced by the geometry of the colloid. In addition, we find that local lipid exchange enhances the pinning of the liquid disordered phase in highly curved regions. This allows us to provide estimates of the bending moduli difference of the domains. Finally, by combining experimental and numerical data, we forecast the outcome of possible experiments on catenoidal and conical necks and show that these geometries could greatly improve the precision of the current estimates of the bending moduli.