Deterministic and stochastic control of kirigami topology

Kirigami, the creative art of paper cutting, is a promising paradigm for mechanical metamaterials. However, to make kirigami-inspired structures a reality requires controlling the topology of kirigami to achieve connectivity and rigidity. We address this question by deriving the maximum number of cuts (minimum number of links) that still allow us to preserve global rigidity and connectivity of the kirigami. A deterministic hierarchical construction method yields an efficient topological way to control both the number of connected pieces and the total degrees of freedom. A statistical approach to the control of rigidity and connectivity in kirigami with random cuts complements the deterministic pathway, and shows that both the number of connected pieces and the degrees of freedom show percolation transitions as a function of the density of cuts (links). Together, this provides a general framework for the control of rigidity and connectivity in planar kirigami.

Collective ventilation in honeybee nests

European honey bees ( Apis mellifera) live in large congested nest cavities with a single opening that limits passive ventilation. When the local air temperature exceeds a threshold, the nests are actively ventilated by bees fanning their wings at the nest entrance. Here, we show that colonies with relatively large nest entrances use an emergent ventilation strategy where fanning bees self-organize to form groups, separating regions of continuous inflow and outflow. The observed spatio-temporal patterns correlate the air velocity and air temperature along the entrances to the distribution of fanning bees. A mathematical model that couples these variables to known fanning behaviour of individuals recapitulates their collective dynamics. Additionally, the model makes predictions about the temporal stability of the fanning group as a function of the temperature difference between the environment and the nest. Consistent with these predictions, we observe that the fanning groups drift, cling to the entrance boundaries, break-up and reform as the ambient temperature varies over a period of days. Overall, our study shows how honeybees use flow-mediated communication to self-organize into a steady state in fluctuating environments.

Topology, Geometry, and Mechanics of Strongly Stretched and Twisted Filaments: Solenoids, Plectonemes, and Artificial Muscle Fibers

Soft elastic filaments that can be stretched, bent, and twisted exhibit a range of topologically and geometrically complex morphologies. Recently, a number of experiments have shown how to use these building blocks to create filament-based artificial muscles that use the conversion of writhe to extension or contraction, exposing the connection between topology, geometry, and mechanics. Here, we combine numerical simulations of soft elastic filaments that account for geometric nonlinearities and self-contact to map out the basic structures underlying artificial muscle fibers in a phase diagram that is a function of the extension and twist density. We then use ideas from computational topology to track the interconversion of link, twist, and writhe in these geometrically complex physical structures to explain the physical principles underlying artificial muscle fibers and provide guidelines for their design.

Computational analysis of size, shape and structure of insect wings

The size, shape and structure of insect wings are intimately linked to their ability to fly. However, there are few systematic studies of the variability of the natural patterns in wing morphology across insects. We have assembled a dataset of 789 insect wings with representatives from 25 families and performed a comprehensive computational analysis of their morphology using topological and geometric notions in terms of (i) wing size and contour shape, (ii) vein topology, and (iii) shape and distribution of wing membrane domains. These morphospaces are complementary to existing methods for quantitatively characterizing wing morphology and are likely to be useful for investigating wing function and evolution. This Methods and Techniques paper is accompanied by a set of computational tools for open use.

This article has an associated First Person interview with the first author of the paper.

Summary: We provide a set of simple quantitative measures to compare morphological variation in size, shape, and structure of insect wings across species, families and orders.

Geometric localization in supported elastic struts

Localized deformation patterns are a common motif in morphogenesis and are increasingly finding applications in materials science and engineering, in such instances as mechanical memories. Here, we describe the emergence of spatially localized deformations in a minimal mechanical system by exploring the impact of growth and shear on the conformation of a semi-flexible filament connected to a pliable shearable substrate. We combine numerical simulations of a discrete rod model with theoretical analysis of the differential equations recovered in the continuum limit to quantify (in the form of scaling laws) how geometry, mechanics and growth act together to give rise to such localized structures in this system. We find that spatially localized deformations along the filament emerge for intermediate shear modulus and increasing growth. Finally, we use experiments on a 3D-printed multi-material model system to demonstrate that external control of the amount of shear and growth may be used to regulate the spatial extent of the localized strain texture.

Size control of the inner ear via hydraulic feedback

Animals make organs of precise size, shape, and symmetry but how developing embryos do this is largely unknown. Here, we combine quantitative imaging, physical theory, and physiological measurement of hydrostatic pressure and fluid transport in zebrafish to study size control of the developing inner ear. We find that fluid accumulation creates hydrostatic pressure in the lumen leading to stress in the epithelium and expansion of the otic vesicle. Pressure, in turn, inhibits fluid transport into the lumen. This negative feedback loop between pressure and transport allows the otic vesicle to change growth rate to control natural or experimentally-induced size variation. Spatiotemporal patterning of contractility modulates pressure-driven strain for regional tissue thinning. Our work connects molecular-driven mechanisms, such as osmotic pressure driven strain and actomyosin tension, to the regulation of tissue morphogenesis via hydraulic feedback to ensure robust control of organ size.

Editorial note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

Programming shape using kirigami tessellations

Kirigami tessellations, regular planar patterns formed by partially cutting flat, thin sheets, allow compact shapes to morph into open structures with rich geometries and unusual material properties. However, geometric and topological constraints make the design of such structures challenging. Here we pose and solve the inverse problem of determining the number, size and orientation of cuts that enables the deployment of a closed, compact regular kirigami tessellation to conform approximately to any prescribed target shape in two or three dimensions. We first identify the constraints on the lengths and angles of generalized kirigami tessellations that guarantee that their reconfigured face geometries can be contracted from a non-trivial deployed shape to a compact, non-overlapping planar cut pattern. We then encode these conditions into a flexible constrained optimization framework to obtain generalized kirigami patterns derived from various periodic tesselations of the plane that can be deployed into a wide variety of prescribed shapes. A simple mechanical analysis of the resulting structure allows us to determine and control the stability of the deployed state and control the deployment path. Finally, we fabricate physical models that deploy in two and three dimensions to validate this inverse design approach. Altogether, our approach, combining geometry, topology and optimization, highlights the potential for generalized kirigami tessellations as building blocks for shape-morphing mechanical metamaterials.

The effect of step size on straight-line orientation

Moving along a straight path is a surprisingly difficult task. This is because, with each ensuing step, noise is generated in the motor and sensory systems, causing the animal to deviate from its intended route. When relying solely on internal sensory information to correct for this noise, the directional error generated with each stride accumulates, ultimately leading to a curved path. In contrast, external compass cues effectively allow the animal to correct for errors in its bearing. Here, we studied straight-line orientation in two different sized dung beetles. This allowed us to characterize and model the size of the directional error generated with each step, in the absence of external visual compass cues (motor error) as well as in the presence of these cues (compass and motor errors). In addition, we model how dung beetles balance the influence of internal and external orientation cues as they orient along straight paths under the open sky. We conclude that the directional error that unavoidably accumulates as the beetle travels is inversely proportional to the step size of the insect, and that both beetle species weigh the two sources of directional information in a similar fashion.

Dynamics of Growth and Form in Prebiotic Vesicles

The growth, form, and division of prebiotic vesicles, membraneous bags of fluid of varying components and shapes is hypothesized to have served as the substrate for the origin of life. The dynamics of these out-of-equilibrium structures is controlled by physicochemical processes that include the intercalation of amphiphiles into the membrane, fluid flow across the membrane, and elastic deformations of the membrane. To understand prebiotic vesicular forms and their dynamics, we construct a minimal model that couples membrane growth, deformation, and fluid permeation, ultimately couched in terms of two dimensionless parameters that characterize the relative rate of membrane growth and the membrane permeability. Numerical simulations show that our model captures the morphological diversity seen in extant precursor mimics of cellular life, and might provide simple guidelines for the synthesis of these complex shapes from simple ingredients.

Spatial control of irreversible protein aggregation

Liquid cellular compartments form in the cyto- or nucleoplasm and can regulate aberrant protein aggregation. Yet, the mechanisms by which these compartments affect protein aggregation remain unknown. Here, we combine kinetic theory of protein aggregation and liquid-liquid phase separation to study the spatial control of irreversible protein aggregation in the presence of liquid compartments. We find that even for weak interactions aggregates strongly partition into the liquid compartment. Aggregate partitioning is caused by a positive feedback mechanism of aggregate nucleation and growth driven by a flux maintaining the phase equilibrium between the compartment and its surrounding. Our model establishes a link between specific aggregating systems and the physical conditions maximizing aggregate partitioning into the compartment. The underlying mechanism of aggregate partitioning could be used to confine cytotoxic protein aggregates inside droplet-like compartments but may also represent a common mechanism to spatially control irreversible chemical reactions in general.

A multiphase theory for spreading microbial swarms and films

Bacterial swarming and biofilm formation are collective multicellular phenomena through which diverse microbial species colonize and spread over water-permeable tissue. During both modes of surface translocation, fluid uptake and transport play a key role in shaping the overall morphology and spreading dynamics. Here we develop a generalized two-phase thin-film model that couples bacterial growth, extracellular matrix swelling, fluid flow, and nutrient transport to describe the expansion of both highly motile bacterial swarms, and sessile bacterial biofilms. We show that swarm expansion corresponds to steady-state solutions in a nutrient-rich, capillarity dominated regime. In contrast, biofilm colony growth is described by transient solutions associated with a nutrient-limited, extracellular polymer stress driven limit. We apply our unified framework to explain a range of recent experimental observations of steady and unsteady expansion of microbial swarms and biofilms. Our results demonstrate how the physics of flow and transport in slender geometries serve to constrain biological organization in microbial communities.

Bacteria can grow and thrive in many different environments. Although we usually think of bacteria as single-celled organisms, they are not always solitary; they can also form groups containing large numbers of individuals. These aggregates work together as one super-colony, allowing the bacteria to feed and protect themselves more efficiently than they could as isolated cells.

These colonies move and grow in characteristic patterns as they respond to their environment. They can form swarms, like insects, or biofilms, which are thin, flat structures containing both cells and a film-like substance that the cells secrete. Availability of food and water influences the way colonies spread; however, since movement and growth are accompanied by mechanical forces, physical constraints are also important. These include the ability of the bacteria to change the water balance and their local mechanical environment, and the forces they create as they grow and move.

Previous research has used a variety of experimental and theoretical approaches to explain the dynamics of bacterial swarms and biofilms as separate phenomena. However, while they do differ biologically, they also share many physical characteristics.

Srinivasan et al. wanted to exploit these similarities, and use them to predict the growth and shape of biofilms and bacterial swarms under different conditions. To do this, a unified mathematical model for the growth of both swarms and biofilms was created. The model accounted for various factors, such as the transport of nutrients into the colony, the movement of water between the colony and the surface on which it grew, and mechanical changes in the environment (e.g. swelling/softening). The theoretical results were then compared with results from experimental measurements of different bacterial aggregates grown on a soft, hydrated gel. For both swarms and biofilms, the model correctly predicted how fast the colony expanded overall, as well as the shape and location of actively growing regions.

Biofilms and other bacterial aggregates can cause diseases and increase inflammation in tissues, and also hinder industrial processes by damage to submerged surfaces, such as ships and waterpipes. The results described here may open up new approaches to restrict the spreading of bacterial aggregates by focusing on their physical constraints.

Morphogenesis of termite mounds

Several species of millimetric-sized termites across Africa, Asia, Australia, and South America collectively construct large, meter-sized, porous mound structures that serve to regulate mound temperature, humidity, and gas concentrations. These mounds display varied yet distinctive morphologies that range widely in size and shape. To explain this morphological diversity, we introduce a mathematical model that couples environmental physics to insect behavior: The advection and diffusion of heat and pheromones through a porous medium are modified by the mound geometry and, in turn, modify that geometry through a minimal characterization of termite behavior. Our model captures the range of naturally observed mound shapes in terms of a minimal set of dimensionless parameters and makes testable hypotheses for the response of mound morphology to external temperature oscillations and internal odors. Our approach also suggests mechanisms by which evolutionary changes in odor production rate and construction behavior coupled to simple physical laws can alter the characteristic mound morphology of termites.

Mechanics of biomimetic 4D printed structures

Recent progress in additive manufacturing and materials engineering has led to a surge of interest in shape-changing plate and shell-like structures. Such structures are typically printed in a planar configuration and, when exposed to an ambient stimulus such as heat or humidity, swell into a desired three-dimensional geometry. Viewed through the lens of differential geometry and elasticity, the application of the physical stimulus can be understood as a local change in the metric of a two dimensional surface embedded in three dimensions. To relieve the resulting elastic frustration, the structure will generally bend and buckle out-of-plane. Here, we propose a numerical approach to convert the discrete geometry of filament bilayers, associated with print paths of inks with given material properties, into continuous plates with inhomogeneous growth patterns and thicknesses. When subject to prescribed growth anisotropies, we can then follow the evolution of the shapes into their final form. We show that our results provide a good correspondence between experiments and simulations, and lead to a framework for the prediction and design of shape-changing structures.

Meniscus instabilities in thin elastic layers

We consider meniscus instabilities in thin elastic layers perfectly adhered to, and confined between, much stiffer bodies. When the free boundary associated with the meniscus of the elastic layer recedes into the layer, for example by pulling the stiffer bodies apart or injecting air between them, then the meniscus will eventually undergo a purely elastic instability in which fingers of air invade the layer. Here we show that the form of this instability is identical in a range of different loading conditions, provided only that the thickness of the meniscus, a, is small compared to the in-plane dimensions and to two emergent in-plane length scales that arise if the substrate is soft or if the layer is compressible. In all such situations, we predict that the instability will occur when the meniscus has receded by approximately 1.27a, and that the instability will have wavelength λ ≈ 2.75a. We illustrate this by also calculating the threshold for fingering in a thin wedge of elastic material bonded to two rigid plates that are pried apart, and the threshold for fingering when a flexible plate is peeled from an elastic layer that glues the plate to a rigid substrate.

Planar morphometrics using Teichmüller maps

Inspired by the question of quantifying wing shape, we propose a computational approach for analysing planar shapes. We first establish a correspondence between the boundaries of two planar shapes with boundary landmarks using geometric functional data analysis and then compute a landmark-matching curvature-guided Teichmüller mapping with uniform quasi-conformal distortion in the bulk. This allows us to analyse the pair-wise difference between the planar shapes and construct a similarity matrix on which we deploy methods from network analysis to cluster shapes. We deploy our method to study a variety of Drosophila wings across species to highlight the phenotypic variation between them, and Lepidoptera wings over time to study the developmental progression of wings. Our approach of combining complex analysis, computation and statistics to quantify, compare and classify planar shapes may be usefully deployed in other biological and physical systems.

The statistical shape of geometric reasoning

Geometric reasoning has an inherent dissonance: its abstract axioms and propositions refer to perfect, idealized entities, whereas its use in the physical world relies on dynamic perception of objects. How do abstract Euclidean concepts, dynamics, and statistics come together to support our intuitive geometric reasoning? Here, we address this question using a simple geometric task – planar triangle completion. An analysis of the distribution of participants’ errors in localizing a fragmented triangle’s missing corner reveals scale-dependent deviations from a deterministic Euclidean representation of planar triangles. By considering the statistical physics of the process characterized via a correlated random walk with a natural length scale, we explain these results and further predict participants’ estimates of the missing angle, measured in a second task. Our model also predicts the results of a categorical reasoning task about changes in the triangle size and shape even when such completion strategies need not be invoked. Taken together, our findings suggest a critical role for noisy physical processes in our reasoning about elementary Euclidean geometry.

Generalized Erdős numbers for network analysis

The identification of relationships in complex networks is critical in a variety of scientific contexts. This includes the identification of globally central nodes and analysing the importance of pairwise relationships between nodes. In this paper, we consider the concept of topological proximity (or 'closeness') between nodes in a weighted network using the generalized Erdős numbers (GENs). This measure satisfies a number of desirable properties for networks with nodes that share a finite resource. These include: (i) real-valuedness, (ii) non-locality and (iii) asymmetry. We show that they can be used to define a personalized measure of the importance of nodes in a network with a natural interpretation that leads to new methods to measure centrality. We show that the square of the leading eigenvector of an importance matrix defined using the GENs is strongly correlated with well-known measures such as PageRank, and define a personalized measure of centrality that is also well correlated with other existing measures. The utility of this measure of topological proximity is demonstrated by showing the asymmetries in both the dynamics of random walks and the mean infection time in epidemic spreading are better predicted by the topological definition of closeness provided by the GENs than they are by other measures.

Multifunctional ferrofluid-infused surfaces with reconfigurable multiscale topography

Developing adaptive materials with geometries that change in response to external stimuli provides fundamental insights into the links between the physical forces involved and the resultant morphologies and creates a foundation for technologically relevant dynamic systems^1,2. In particular, reconfigurable surface topography as a means to control interfacial properties³ has recently been explored using responsive gels⁴, shape-memory polymers⁵, liquid crystals^6-8 and hybrid composites^9-14, including magnetically active slippery surfaces^12-14. However, these designs exhibit a limited range of topographical changes and thus a restricted scope of function. Here we introduce a hierarchical magneto-responsive composite surface, made by infiltrating a ferrofluid into a microstructured matrix (termed ferrofluid-containing liquid-infused porous surfaces, or FLIPS). We demonstrate various topographical reconfigurations at multiple length scales and a broad range of associated emergent behaviours. An applied magnetic-field gradient induces the movement of magnetic nanoparticles suspended in the ferrofluid, which leads to microscale flow of the ferrofluid first above and then within the microstructured surface. This redistribution changes the initially smooth surface of the ferrofluid (which is immobilized by the porous matrix through capillary forces) into various multiscale hierarchical topographies shaped by the size, arrangement and orientation of the confining microstructures in the magnetic field. We analyse the spatial and temporal dynamics of these reconfigurations theoretically and experimentally as a function of the balance between capillary and magnetic pressures^15-19 and of the geometric anisotropy of the FLIPS system. Several interesting functions at three different length scales are demonstrated: self-assembly of colloidal particles at the micrometre scale; regulated flow of liquid droplets at the millimetre scale; and switchable adhesion and friction, liquid pumping and removal of biofilms at the centimetre scale. We envision that FLIPS could be used as part of integrated control systems for the manipulation and transport of matter, thermal management, microfluidics and fouling-release materials.

Differential Activity-Driven Instabilities in Biphasic Active Matter

Active stresses can cause instabilities in contractile gels and living tissues. Here we provide a generic hydrodynamic theory that treats these systems as a mixture of two phases of varying activity and different mechanical properties. We find that differential activity between the phases causes a uniform mixture to undergo a demixing instability. We follow the nonlinear evolution of the instability and characterize a phase diagram of the resulting patterns. Our study complements other instability mechanisms in mixtures driven by differential adhesion, differential diffusion, differential growth, and differential motion.

Microtubules soften due to cross-sectional flattening

We use optical trapping to continuously bend an isolated microtubule while simultaneously measuring the applied force and the resulting filament strain, thus allowing us to determine its elastic properties over a wide range of applied strains. We find that, while in the low-strain regime, microtubules may be quantitatively described in terms of the classical Euler-Bernoulli elastic filament, above a critical strain they deviate from this simple elastic model, showing a softening response with increasingdeformations. A three-dimensional thin-shell model, in which the increased mechanical compliance is caused by flattening and eventual buckling of the filament cross-section, captures this softening effect in the high strain regime and yields quantitative values of the effective mechanical properties of microtubules. Our results demonstrate that properties of microtubules are highly dependent on the magnitude of the applied strain and offer a new interpretation for the large variety in microtubule mechanical data measured by different methods.

Forward and inverse problems in the mechanics of soft filaments

Soft slender structures are ubiquitous in natural and artificial systems, in active and passive settings and across scales, from polymers and flagella, to snakes and space tethers. In this paper, we demonstrate the use of a simple and practical numerical implementation based on the Cosserat rod model to simulate the dynamics of filaments that can bend, twist, stretch and shear while interacting with complex environments via muscular activity, surface contact, friction and hydrodynamics. We validate our simulations by solving a number of forward problems involving the mechanics of passive filaments and comparing them with known analytical results, and extend them to study instabilities in stretched and twisted filaments that form solenoidal and plectonemic structures. We then study active filaments such as snakes and other slender organisms by solving inverse problems to identify optimal gaits for limbless locomotion on solid surfaces and in bulk liquids.

Solar-powered ventilation of African termite mounds

How termite mounds function to facilitate climate control is still only partially understood. Recent experimental evidence in the mounds of a single species, the south Asian termite Odontotermes obesus, suggests that the daily oscillations of radiant heating associated with diurnal insolation patterns drive convective flow within them. How general this mechanism is remains unknown. To probe this, we consider the mounds of the African termite Macrotermes michaelseni, which thrives in a very different environment. By directly measuring air velocities and temperatures within the mound, we see that the overall mechanisms and patterns involved are similar to that in the south Asian species. However, there are also some notable differences between the physiology of these mounds associated with the temporal variations in radiant heating patterns and CO₂ dynamics. Because of the difference between direct radiant heating driven by the position of the sun in African conditions, and the more shaded south Asian environments, we see changes in the convective flows in the two types of mounds. Furthermore, we also see that the south Asian mounds show a significant overturning of stratified gases, once a day, while the African mounds have a relatively uniform concentration of CO₂ Overall, our observations show that despite these differences, termite architectures can harness periodic solar heating to drive ventilation inside them in very different environments, functioning as an external lung, with clear implications for human engineering.

Competing failure modes in finite adhesive pads

Thin adhesive pads used to attach objects to each other often fail catastrophically. Here we consider the nature of failure of such a pad under loading parallel to the adhesive substrate. To determine the modes of failure of the pad and to understand what limits its load bearing capacity, we conduct experiments with finite pads composed of a soft adhesive layer with a stiff backing and load them parallel to the surface of adhesion. We find that two different peeling mechanisms emerge as a function of the slenderness of the adhesive pad: an interfacial peeling mechanism that starts close to the pulling end for very long pads, and an unstable curling mechanism that starts at the opposite end for relatively short pads. A minimal theoretical framework allows us to explain our observations and reveals the adhesive bond stiffness as a dominant parameter in defining the peeling mode. A phase diagram that delineates the different regimes of peeling modes brings our experiments and theory together. Our results suggest that unstable peeling by curling may be more common than previously thought, and could perhaps occur naturally in such examples as the gecko foot.