Source author record

L. Mahadevan

L. Mahadevan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

83works

40topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

Stigmergic optimal transport

Efficient navigation in swarms often relies on the emergence of decentralized approaches that minimize traversal time or energy. Stigmergy, where agents modify a shared environment that then modifies their behavior, is a classic mechanism that can encode this strategy. We develop a theoretical framework for stigmergic transport by casting it as a stochastic optimal control problem: agents (collectively) lay and (individually) follow trails while minimizing expected traversal time. Simulations and analysis reveal two emergent behaviors: path straightening in homogeneous environments and path refraction at material interfaces, both consistent with experimental observations of insect trails. While reminiscent of Fermat's principle, our results show how local, noisy agent+field interactions can give rise to geodesic trajectories in heterogeneous environments, without centralized coordination or global knowledge, relying instead on an embodied slow fast dynamical mechanism.

preprint2022arXiv

How to grow a flat leaf

Growing a flat lamina such as a leaf is almost impossible without some feedback to stabilize long wavelength modes that are easy to trigger since they are energetically cheap. Here we combine the physics of thin elastic plates with feedback control theory to explore how a leaf can remain flat while growing. We investigate both in-plane (metric) and out-of-plane (curvature) growth variation and account for both local and nonlocal feedback laws. We show that a linearized feedback theory that accounts for both spatially nonlocal and temporally delayed effects suffices to suppress long wavelength fluctuations effectively and explains recently observed statistical features of growth in tobacco leaves. Our work provides a framework for understanding the regulation of the shape of leaves and other laminar objects.

preprint2022arXiv

Modular representation and control of floppy networks

Geometric graph models of systems as diverse as proteins, robots, and mechanical structures from DNA assemblies to architected materials point towards a unified way to represent and control them in space and time. While much work has been done in the context of characterizing the behavior of these networks close to critical points associated with bond and rigidity percolation, isostaticity, etc., much less is known about floppy, under-constrained networks that are far more common in nature and technology. Here we combine geometric rigidity and algebraic sparsity to provide a framework for identifying the zero-energy floppy modes via a representation that illuminates the underlying hierarchy and modularity of the network, and thence the control of its nestedness and locality. Our framework allows us to demonstrate a range of applications of this approach that include robotic reaching tasks with motion primitives, and predicting the linear and nonlinear response of elastic networks based solely on infinitesimal rigidity and sparsity, which we test using physical experiments. Our approach is thus likely to be of use broadly in dissecting the geometrical properties of floppy networks using algebraic sparsity to optimize their function and performance.

preprint2022arXiv

Quasicrystal kirigami

Kirigami, the art of introducing cuts in thin sheets to enable articulation and deployment, has become an inspiration for a novel class of mechanical metamaterials with unusual properties. Here we complement the use of periodic tiling patterns for kirigami designs by showing that quasicrystals can also serve as the basis for designing deployable kirigami structures, and analyze the geometrical, topological and mechanical properties of these aperiodic kirigami structures.

preprint2021arXiv

Active nematic defects and epithelial morphogenesis

Inspired by recent experiments that highlight the role of nematic defects in the morphogenesis of epithelial tissues, we develop a minimal framework to study the dynamics of an active curved surface driven by its nematic texture. Allowing the surface to evolve via relaxational dynamics leads to a theory linking nematic defect dynamics, cellular division rates and Gaussian curvature. Regions of large positive (negative) curvature and positive (negative) growth are colocalized with the presence of positive (negative) defects. Applying this framework to the dynamics of cultured murine neural progenitor cells (NPCs) in an ex-vivo setting, we find that cells accumulate at positive defects and are depleted at negative defects. In contrast, applying this to the dynamics of a basal marine invertebrate \emph{Hydra} in an in-vivo setting, we show that activity stabilizes a bound $+1$ defect state by creating an incipient tentacle, while a bound $+1$ defect state surrounded by two $-1/2$ defects can create a stationary ring configuration of tentacles, consistent with observations.

preprint2021arXiv

Combing a double helix

Combing hair involves brushing away the topological tangles in a collective curl. Using a combination of experiment and computation, we study this problem that naturally links topology, geometry and mechanics. Observations show that the dominant interactions in hair are those of a two-body nature, corresponding to a braided homochiral double helix. Using this minimal model, we study the detangling of an elastic double helix via a single stiff tine that moves along it, leaving two untangled filaments in its wake. Our results quantify how the forces of detangling are correlated with the magnitude and spatial extent of the link density, a topological quantity, that propagates ahead of the tine. This in turn provides a measure of the maximum characteristic length of a single combing stroke, and thus the trade-offs between comfort, efficiency and speed of combing in the many-body problem on a head of hair.

preprint2021arXiv

Geometric control of topological dynamics in a singing saw

The common handsaw can be converted into a bowed musical instrument capable of producing exquisitely sustained notes when its blade is appropriately bent. Acoustic modes localized at an inflection point are known to underlie the saw's sonorous quality, yet the origin of localization has remained mysterious. Here we uncover a topological basis for the existence of localized modes, that relies on and is protected by spatial curvature. By combining experimental demonstrations, theory and computation, we show how spatial variations in blade curvature control the localization of these trapped states, allowing the saw to function as a geometrically tunable high quality oscillator. Our work establishes an unexpected connection between the dynamics of thin shells and topological insulators, and offers a robust principle to design high quality resonators across scales, from macroscopic instruments to nanoscale devices, simply through geometry.

preprint2021arXiv

Geometrical dynamics of edge-driven surface growth

Accretion of mineralized thin wall-like structures via localized growth along their edges is observed in a range of physical and biological systems ranging from molluscan and brachiopod shells to carbonate-silica composite precipitates. To understand the shape of these mineralized structures, we develop a mathematical framework that treats the thin-walled shells as a smooth surface left in the wake of the growth front that can be described as an evolving space curve. Our theory then takes an explicit geometric form for the prescription of the velocity of the growth front curve, along with some compatibility relations and a closure equation related to the nature of surface curling. The result is a set of equations for the geometrical dynamics of a curve that leaves behind a compatible surface. Solutions of these equations capture a range of geometric precipitate patterns seen in abiotic and biotic forms across scales. In addition to providing a framework for the growth and form of these thin-walled morphologies, our theory suggests a new class of dynamical systems involving moving space curves that are compatible with non-Euclidean embeddings of surfaces.

preprint2021arXiv

Instability-induced patterning of a jelling jet

When a thin stream of aqueous sodium alginate is extruded into a reacting calcium chloride bath, it polymerizes into a soft elastic tube that spontaneously forms helical coils due to the ambient fluid drag. We quantify the onset of this drag-induced instability and its nonlinear evolution using experiments, and explain the results using a combination of scaling, theory and simulations. By co-extruding a second (internal) liquid within the aqueous sodium alginate jet and varying the rates of co-extrusion of the two liquids, as well as the diameter of the jet, we show that we can tune the local composition of the composite filament and the nature of the ensuing instabilities to create soft filaments of variable relative buoyancy, shape and mechanical properties. All together, by harnessing the fundamental varicose (jetting) and sinuous (buckling) instabilities associated with the extrusion of a jelling filament, we show that it is possible to print complex three-dimensional filamentous structures in the ambient fluid.

preprint2021arXiv

Prestrain-induced contraction in 1D random elastic chains

Prestrained elastic networks arise in a number of biological and technological systems ranging from the cytoskeleton of cells to tensegrity structures. To understand the response of such a network as a function of the prestrain, we consider a minimal model in one dimension. We do this by considering a chain (1D network) of elastic springs upon which a random, zero mean, finite variance prestrain is imposed. Numerical simulations and analytical predictions quantify the magnitude of the contraction as a function of the variance of the prestrain, and show that the chain always shrinks. To test these predictions, we vary the topology of the chain and consider more complex connectivity and show that our results are relatively robust to these changes.

preprint2021arXiv

Static adhesion hysteresis in elastic structures

Adhesive interactions between elastic structures such as graphene sheets, carbon nanotubes, and microtubules have been shown to exhibit hysteresis due to irrecoverable energy loss associated with bond breakage, even in static (rate-independent) experiments. To understand this phenomenon, we start with a minimal theory for the peeling of a thin sheet from a substrate, coupling the local event of bond breaking to the nonlocal elastic relaxation of the sheet and show that this can drive static adhesion hysteresis over a bonding/debonding cycle. Using this model we quantify hysteresis in terms of the adhesion and elasticity parameters of the system. This allows us to derive a scaling relation that preserves hysteresis at different levels of granularity while resolving a seeming paradox of lattice trapping in the continuum limit of a discrete fracture process. Finally, to verify our theory, we use new experiments to demonstrate and measure adhesion hysteresis in bundled microtubules.

preprint2020arXiv

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.

preprint2020arXiv

Models of benthic bipedalism

Walking is a common bipedal and quadrupedal gait and is often associated with terrestrial and aquatic organisms. Inspired by recent evidence of the neural underpinnings of primitive aquatic walking in the little skate Leucoraja erinacea, we introduce a theoretical model of aquatic walking that reveals robust and efficient gaits with modest requirements for body morphology and control. The model predicts undulatory behavior of the system body with a regular foot placement pattern which is also observed in the animal, and additionally predicts the existence of gait bistability between two states, one with a large energetic cost for locomotion and another associated with almost no energetic cost. We show that these can be discovered using a simple reinforcement learning scheme. To test these theoretical frameworks, we built a bipedal robot and show that its behaviors are similar to those of our minimal model: its gait is also periodic and exhibits bistability, with a low efficiency gait separated from a high efficiency gait by a "jump" transition. Overall, our study highlights the physical constraints on the evolution of walking and provides a guide for the design of efficient biomimetic robots.

preprint2020arXiv

Optimal policies for mitigating pandemic costs

Several non-pharmaceutical interventions have been proposed to control the spread of the COVID-19 pandemic. On the large scale, these empirical solutions, often associated with extended and complete lockdowns, attempt to minimize the costs associated with mortality, economic losses and social factors, while being subject to constraints such as finite hospital capacity. Here we pose the question of how to mitigate pandemic costs subject to constraints by adopting the language of optimal control theory. This allows us to determine top-down policies for the nature and dynamics of social contact rates given an age-structured model for the dynamics of the disease. Depending on the relative weights allocated to life and socioeconomic losses, we see that the optimal strategies range from long-term social-distancing only for the most vulnerable, to partial lockdown to ensure not over-running hospitals, to alternating-shifts with significant reduction in life and/or socioeconomic losses. Crucially, commonly used strategies that involve long periods of broad lockdown are almost never optimal, as they are highly unstable to reopening and entail high socioeconomic costs. Using parameter estimates from data available for Germany and the USA, we quantify these policies and use sensitivity analysis in the relevant model parameters and initial conditions to determine the range of robustness of our policies. Finally we also discuss how bottom-up behavioral changes can also change the dynamics of the pandemic and show how this in tandem with top-down control policies can mitigate pandemic costs even more effectively.

preprint2020arXiv

Self-excited motions of volatile drops on swellable sheets

When a volatile droplet is deposited on a floating swellable sheet, it becomes asymmetric, lobed and mobile. We describe and quantify this phenomena that involves nonequilibrium swelling, evaporation and motion, working together to realize a self-excitable spatially extended oscillator. Solvent penetration causes the film to swell locally and eventually buckle, changing its shape and the drop responds by moving. Simultaneously, solvent evaporation from the swollen film causes it to regain its shape once the droplet has moved away. The process repeats and leads to complex pulsatile spinning and/or sliding movements. We use a one-dimensional experiment to highlight the slow swelling of and evaporation from the film and the fast motion of the drop, a characteristic of excitable systems. Finally, we provide a phase diagram for droplet excitability as a function of drop size and film thickness and scaling laws for the motion of the droplet.

preprint2019arXiv

Random sequential adsorption of spheres on a cylinder

Inspired by observations of beads packed on a thin string in such systems as sea-grapes and dental plaque, we study the random sequential adsorption of spheres on a cylinder. We determine the asymptotic fractional coverage of the cylinder as a function of the sole parameter in the problem, the ratio of the sphere radius to the cylinder radius (for a very long cylinder) using a combination of analysis and numerical simulations. Examining the asymptotic structures, we find weak chiral ordering on sufficiently small spatial scales. Experiments involving colloidal microspheres that can attach irreversibly to a silica wire via electrostatic forces or DNA hybridization allow us to verify our predictions for the asymptotic coverage.

preprint2018arXiv

Programming shape using kirigami tessellations

Kirigami tessellations, regular planar patterns formed by cutting flat, thin sheets, have attracted recent scientific interest for their rich geometries, surprising material properties and promise for technologies. Here we pose and solve the inverse problem of designing the number, size, and orientation of cuts that allows us to convert a closed, compact regular kirigami tessellation of the plane into a deployment that conforms approximately to any prescribed target shape in two and three dimensions. We do this by first identifying the constraints on the lengths and angles of generalized kirigami tessellations which guarantee that their reconfigured face geometries can be contracted from a non-trivial deployed shape to a novel planar cut pattern. We encode these conditions in a flexible constrained optimization framework which allows us to deform the geometry of periodic kirigami tesselations with three, four, and sixfold symmetry, among others, into generalized kirigami patterns that deploy to a wide variety of prescribed boundary target shapes. Physically fabricated models verify our inverse design approach and allow us to determine the tunable material response of the resulting structures. We then extend our framework to create generalized kirigami patterns that deploy to approximate curved surfaces in $\mathbb{R}^3$. Altogether, this work illustrates a novel framework for designing complex, shape-changing sheets from simple cuts showing the power of kirigami tessellations as flexible mechanical metamaterials.

preprint2017arXiv

Active elastohydrodynamics of vesicles in narrow, blind constrictions

Fluid-resistance limited transport of vesicles through narrow constrictions is a recurring theme in many biological and engineering applications. Inspired by the motor-driven movement of soft membrane-bound vesicles into closed neuronal dendritic spines, here we study this problem using a combination of passive three-dimensional simulations and a simplified semi-analytical theory for active transport of vesicles that are forced through such constrictions by molecular motors. We show that the motion of these objects is characterized by two dimensionless quantities related to the geometry and the strength of forcing relative to the vesicle elasticity. We use numerical simulations to characterize the transit time for a vesicle forced by fluid pressure through a constriction in a channel, and find that relative to an open channel, transport into a blind end leads to the formation of an effective lubrication layer that strongly impedes motion. When the fluid pressure forcing is complemented by forces due to molecular motors that are responsible for vesicle trafficking into dendritic spines, we find that the competition between motor forcing and fluid drag results in multistable dynamics reminiscent of the real system. Our study highlights the role of non-local hydrodynamic effects in determining the kinetics of vesicular transport in constricted geometries.

preprint2016arXiv

Self-sustained lift and low friction via soft lubrication

Relative motion between soft wet solids arises in a number of applications in natural and artificial settings, and invariably couples elastic deformation and fluid flow. We explore this in a minimal setting by considering a fluid-immersed negatively-buoyant cylinder moving along a soft inclined wall. Our experiments show that there is an emergent robust steady-state sliding regime of the cylinder with an effective friction that is significantly reduced relative to that of rigid fluid-lubricated contacts. A simple scaling approach that couples the cylinder-induced flow to substrate deformation allows us to explain the emergence of an elastohydrodynamic lift that underlies the self-sustained lubricated motion of the cylinder, consistent with recent theoretical predictions. Our results suggest an explanation for a range of effects such as reduced wear in animal joints and long-runout landslides, and can be couched as a design principle for low-friction interfaces.

preprint2015arXiv

Elastic Cheerios effect: self-assembly of cylinders on a soft solid

A rigid cylinder placed on a soft gel deforms its surface. When multiple cylinders are placed on the surface, they interact with each other via the topography of the deformed gel which serves as an energy landscape; as they move, the landscape changes which in turn changes their interaction. We use a combination of experiments, simple scaling estimates and numerical simulations to study the self-assembly of cylinders in this elastic analog of the Cheerios effect for capillary interactions on a fluid interface. Our results show that the effective two body interaction can be well described by an exponential attraction potential as a result of which the dynamics also show an exponential behavior with respect to the separation distance. When many cylinders are placed on the gel, the cylinders cluster together if they are not too far apart; otherwise their motion gets elastically arrested.

preprint2015arXiv

Elastohydrodynamics and kinetics of protein patterning in the immunological synapse

The cellular basis for the adaptive immune response during antigen recognition relies on a specialized protein interface known as the immunological synapse (IS). Understanding the biophysical basis for protein patterning by deciphering the quantitative rules for their formation and motion is an important aspect of characterizing immune cell recognition and thence the rules for immune system activation. We propose a minimal mathematical model for the physical basis of membrane protein patterning in the IS, which encompass membrane mechanics, protein binding kinetics and motion, and fluid flow in the synaptic cleft. Our theory leads to simple predictions for the spatial and temporal scales of protein cluster formation, growth and arrest as a function of membrane stiffness, rigidity and kinetics of the adhesive proteins, and the fluid in the synaptic cleft. Numerical simulations complement these scaling laws by quantifying the nucleation, growth and stabilization of proteins domains on the size of the cell. Direct comparison with experiment shows that passive elastohydrodynamics and kinetics of protein binding in the synaptic cleft can describe the short-time formation and organization of protein clusters, without evoking any active processes in the cytoskeleton. Despite the apparent complexity of the process, our analysis highlights the role of just two dimensionless parameters that characterize the spatial and temporal evolution of the protein pattern: a ratio of membrane elasticity to protein elasticity, and the ratio of a hydrodynamic time scale for fluid flow relative to the protein binding rate, and we present a simple phase diagram that encompasses the variety of patterns that can arise.

preprint2015arXiv

Elastohydrodynamics of a sliding, spinning and sedimenting cylinder near a soft wall

We consider the motion of a fluid-immersed negatively buoyant particle in the vicinity of a thin compressible elastic wall, a situation that arises in a variety of technological and natural settings. We use scaling arguments to establish different regimes of sliding, and complement these estimates using thin-film lubrication dynamics to determine an asymptotic theory for the sedimentation, sliding, and spinning motions of a cylinder. The resulting theory takes the form of three coupled nonlinear singular-differential equations. Numerical integration of the resulting equations confirms our scaling relations and further yields a range of unexpected behaviours. Despite the low-Reynolds feature of the flow, we demonstrate that the particle can spontaneously oscillate when sliding, can generate lift via a Magnus-like effect, can undergo a spin-induced reversal effect, and also shows an unusual sedimentation singularity. Our description also allows us to address a sedimentation-sliding transition that can lead to the particle coasting over very long distances, similar to certain geophysical phenomena. Finally, we show that a small modification of our theory allows to generalize the results to account for additional effects such as wall poroelasticity.

preprint2015arXiv

Elastohydrodynamics of contact in adherent sheets

Adhesive contact between a thin elastic sheet and a substrate in a liquid environment arises in a range of biological, physical and technological applications. By considering the dynamics of this process that naturally couples fluid flow, long wavelength elastic deformations and microscopic adhesion, and solving the resulting partial differential equation numerically, we uncover the short-time dynamics of the onset of adhesion and the long-time dynamics of a steady propagating adhesion front. Simple scaling laws corroborate our results for characteristic waiting-time for adhesive contact, as well as the speed of the adhesion front. A similarity analysis of the governing partial differential equation further allows us to determine the shape of a fluid filled bump ahead of the adhesion zone. Finally, our analysis yields the boundary conditions for the apparent elastohydrodynamic contact line, generalizing the well known conditions for static elastic contact while highlighting how microscale physics regularizes the dynamics of contact.

preprint2015arXiv

Generalized Erdos Numbers for network analysis

In this paper we consider the concept of `closeness' between nodes in a weighted network that can be defined topologically even in the absence of a metric. The Generalized Erdős Numbers (GENs) satisfy a number of desirable properties as a measure of topological closeness when nodes share a finite resource between nodes as they are real-valued and non-local, and can be used to create an asymmetric matrix of connectivities. 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 a new global measure of centrality and is highly correlated with Page Rank. The relative asymmetry of the GENs (due to their non-metric definition) is linked also to the asymmetry in the mean first passage time between nodes in a random walk, and we use a linearized form of the GENs to develop a continuum model for `closeness' in spatial networks. As an example of their practicality, we deploy them to characterize the structure of static networks and show how it relates to dynamics on networks in such situations as the spread of an epidemic.

preprint2015arXiv

Gyrification from constrained cortical expansion

The exterior of the mammalian brain - the cerebral cortex - has a conserved layered structure whose thickness varies little across species. However, selection pressures over evolutionary time scales have led to cortices that have a large surface area to volume ratio in some organisms, with the result that the brain is strongly convoluted into sulci and gyri. Here we show that the gyrification can arise as a nonlinear consequence of a simple mechanical instability driven by tangential expansion of the gray matter constrained by the white matter. A physical mimic of the process using a layered swelling gel captures the essence of the mechanism, and numerical simulations of the brain treated as a soft solid lead to the formation of cusped sulci and smooth gyri similar to those in the brain. The resulting gyrification patterns are a function of relative cortical expansion and relative thickness (compared with brain size), and are consistent with observations of a wide range of brains, ranging from smooth to highly convoluted. Furthermore, this dependence on two simple geometric parameters that characterize the brain also allows us to qualitatively explain how variations in these parameters lead to anatomical anomalies in such situations as polymicrogyria, pachygyria, and lissencephalia.

preprint2015arXiv

Optimal switching strategies for stochastic geocentric/egocentric navigation

Animals use a combination of egocentric navigation driven by the internal integration of environmental cues, interspersed with geocentric course correction and reorientation, often with uncertainty in sensory acquisition of information, planning and execution. Inspired directly by observations of dung beetle navigational strategies that show switching between geocentric and egocentric strategies, we consider the question of optimal strategies for the navigation of an agent along a preferred direction in the presence of multiple sources of noise. We address this using a model that takes the form of a correlated random walk at short time scales that is interspersed with reorientation events that yields a biased random walks at long time scales. We identify optimal alternation schemes and characterize their robustness in the context of noisy sensory acquisition, and performance errors linked with variations in environmental conditions and agent-environment interactions.

preprint2015arXiv

Persistence of activity in noisy motor-filament assemblies

Long, elastic filaments cross-linked and deformed by active molecular motors occur in various natural settings. The overall macroscopic mechanical response of such a composite network depends on the coupling between the active and the passive properties of the underlying constituents and nonlocal interactions between different parts of the composite. In a simple one dimensional system, using a mean field model, it has been shown that the combination of motor activity and finite filament extensibility yields a persistence length scale over which strain decays. Here we study a similar system, in the complementary limit of strong noise and moderate extensibility, using Brownian multi-particle collision dynamics-based numerical simulations that includes the coupling between motor kinetics and local filament extensibility. While the numerical model shows deviations from the mean field predictions due to the presence of strong active noise caused by the variations in individual motor activity, several qualitative features are still retained. Specifically, for fixed motor attachment and detachment rates, the decay is length is set by the ratio of the passive elasticity to the active shear resistance generated by attached motors. Our study generalizes the notion of persistence in passive thermal systems to actively driven systems with testable predictions.

preprint2015arXiv

Similarity and singularity in adhesive elastohydrodynamic touchdown

We consider the touchdown of an elastic sheet as it adheres to a wall, which has a dynamics that is limited by the viscous resistance provided by the squeeze flow of the intervening liquid trapped between the two solid surfaces. The dynamics of the sheet is described mathematically by elastohydrodynamic lubrication theory, coupling the elastic deformation of the sheet, the microscopic van der Waals adhesion and the viscous thin film flow. We use a combination of numerical simulations of the governing partial differential equation and a scaling analysis to describe the self-similar solution of the touchdown of the sheet as it approaches the wall. An analysis of the equation satisfied by the similarity variables in the vicinity of the touchdown event shows that an entire sequence of solutions are allowed. However, a comparison of these shows that only the fundamental similarity solution is observed in the time-dependent numerical simulations, consistent with the fact that it alone is stable. Our analysis generalizes similar approaches for rupture in capillary thin film hydrodynamics and suggests experimentally verifiable predictions for a new class of singular flows linking elasticity, hydrodynamics and adhesion.

preprint2015arXiv

Solid friction between soft filaments

Any macroscopic deformation of a filamentous bundle is necessarily accompanied by local sliding and/or stretching of the constituent filaments. Yet the nature of the sliding friction between two aligned filaments interacting through multiple contacts remains largely unexplored. Here, by directly measuring the sliding forces between two bundled F-actin filaments, we show that these frictional forces are unexpectedly large, scale logarithmically with sliding velocity as in solid-like friction, and exhibit complex dependence on the filaments' overlap length. We also show that a reduction of the frictional force by orders of magnitude, associated with a transition from solid-like friction to Stokes' drag, can be induced by coating F-actin with polymeric brushes. Furthermore, we observe similar transitions in filamentous microtubules and bacterial flagella. Our findings demonstrate how altering a filament's elasticity, structure and interactions can be used to engineer interfilament friction and thus tune the properties of fibrous composite materials.

preprint2015arXiv

The Organization and Control of an Evolving Interdependent Population

Starting with Darwin, biologists have asked how populations evolve from a low fitness state that is evolutionarily stable to a high fitness state that is not. Specifically of interest is the emergence of cooperation and multicellularity where the fitness of individuals often appears in conflict with that of the population. Theories of social evolution and evolutionary game theory have produced a number of fruitful results employing two-state two-body frameworks. In this study we depart from this tradition and instead consider a multi-player, multi-state evolutionary game, in which the fitness of an agent is determined by its relationship to an arbitrary number of other agents. We show that populations organize themselves in one of four distinct phases of interdependence depending on one parameter, selection strength. Some of these phases involve the formation of specialized large-scale structures. We then describe how the evolution of independence can be manipulated through various external perturbations.

preprint2014arXiv

2D and Axisymmetric Incompressible Elastic Green's Functions

We compile a list of 2-D and axisymmetric Green's functions for isotropic full and half spaces, to complement our letter "Linear elasticity of incompressible solids". We also extend the isotropic exactly incompressible linear theory from our letter to include isotropic neo-Hookean solids subject to a large pre-strain, and present Green's functions in these cases. The Green's function for a pre-strained half-space reproduces the Biot instability.

preprint2014arXiv

Dynamics of evaporative colloidal patterning

Drying suspensions often leave behind complex patterns of particulates, as might be seen in the coffee stains on a table. Here we consider the dynamics of periodic band or uniform solid film formation on a vertical plate suspended partially in a drying colloidal solution. Direct observations allow us to visualize the dynamics of the band and film deposition, and the transition in between when the colloidal concentration is varied. A minimal theory of the liquid meniscus motion along the plate reveals the dynamics of the banding and its transition to the filming as a function of the ratio of deposition and evaporation rates. We also provide a complementary multiphase model of colloids dissolved in the liquid, which couples the inhomogeneous evaporation at the evolving meniscus to the fluid and particulate flows and the transition from a dilute suspension to a porous plug. This allows us to determine the concentration dependence of the bandwidth and the deposition rate. Together, our findings allow for the control of drying-induced patterning as a function of the colloidal concentration and evaporation rate.

preprint2014arXiv

Elastocapillary coalescence of plates and pillars

When a fluid-immersed array of lamellae or filaments that is attached to a substrate is dried, evaporation leads to the formation of menisci on the tips of the plates or pillars that bring them together. Similarly, when hair dries it clumps together due to capillary forces induced by the liquid menisci between the flexible hairs. Building on prior experimental observations, we use a combination of theory and computation to understand the nature of this instability and its evolution in both the two-dimensional and three-dimensional setting of the problem. For the case of lamellae, we explicitly derive the interaction torques based on the relevant physical parameters. A Bloch-wave analysis for our periodic mechanical system captures the critical volume of the liquid and the 2-plate-collapse eigenmode at the onset of instability. We study the evolution of clusters and their arrest using numerical simulations to explain the hierarchical cluster formation and characterize the sensitive dependence of the final structures on the initial perturbations. We then generalize our analysis to treat the problem of pillar collapse in 3D, where the fluid domain is completely connected and the interface is a surface with the uniform mean curvature. Our theory and simulations capture the salient features of both previous experimental observations and our own in terms of the key parameters that can be used to control the kinetics of the process.

preprint2014arXiv

Evaporation-driven ring and film deposition from colloidal droplets

Evaporating suspensions of colloidal particles lead to the formation of a variety of patterns, ranging from a left-over ring of a dried coffee drop to uniformly distributed solid pigments left behind wet paint. To characterize the transition between single rings, multiple concentric rings, broad bands, and uniform deposits, we investigate the dynamics of a drying droplet via a multiphase model of colloidal particles in a solvent. Our theory couples the inhomogeneous evaporation at the evolving droplet interface to the dynamics inside the drop, i.e. the liquid flow, local variations of the particle concentration, and the propagation of the deposition front where the solute forms an incompressible porous medium at high concentrations. A dimensionless parameter combining the capillary number and the droplet aspect ratio captures the formation conditions of different pattern types.

preprint2014arXiv

Exactly isochoric deformations of soft solids

Many materials of contemporary interest, such as gels, biological tissues and elastomers, are easily deformed but essentially incompressible. Traditional linear theory of elasticity implements incompressibility only to first order and thus permits some volume changes, which become problematically large even at very small strains. Using a mixed coordinate transformation originally due to Gauss, we enforce the constraint of isochoric deformations exactly to develop a linear theory with perfect volume conservation that remains valid until strains become geometrically large. We demonstrate the utility of this approach by calculating the response of an infinite soft isochoric solid to a point force that leads to a nonlinear generalization of the Kelvin solution. Our approach naturally generalizes to a range of problems involving deformations of soft solids and interfaces in 2 dimensional and axisymmetric geometries, which we exemplify by determining the solution to a distributed load that mimics muscular contraction within the bulk of a soft solid.

preprint2014arXiv

Feedback induced phase transitions in active porous media

Flow through passive porous media is typically described in terms of a linear theory relating current fluxes and driving forces, in the presence of a prescribed heterogeneous permeability. However, many porous systems such as glacial drainage networks, erosional river bed networks, vascular networks, social insect swarms and animal architectures such as termite mounds are continuously remodeled by the flow and thence modify the flow, i.e. they are active. Here we consider a minimal model for an active porous medium where flow and resistance are coupled to each other. Using numerical simulations, we show that this results in both channelization and wall-building transitions depending on the form of the feedback. A continuum model allows us to understand the qualitative features of the resulting phase diagram, and suggests ways to realize complex architectures using simple rules in engineered systems.

preprint2014arXiv

Fluid driven fingering instability of a confined elastic meniscus

When a fluid is pumped into a cavity in a confined elastic layer, at a critical pressure, destabilizing fingers of fluid invade the elastic solid along its meniscus (Saintyves, Dauchot, and Bouchaud, 2013). These fingers occur without fracture or loss of adhesion and are reversible, disappearing when the pressure is decreased. We develop an asymptotic theory of pressurized highly elastic layers trapped between rigid bodies to explain these observations, with predictions for the critical fluid pressure for fingering, and the finger wavelength. We also show that the theory links this fluid-driven fingering with a similar transition driven instead by transverse stretching of the elastic layer. We further verify these predictions by using finite-element simulations on the two systems which show that, in both cases, the fingering transition is first-order (sudden) and hence has a region of bistability. Our predictions are in good agreement with recent observations of this elastic analog of the classical Saffman-Taylor interfacial instability in hydrodynamics.

preprint2014arXiv

Forced tearing of ductile and brittle thin sheets

Tearing a thin sheet by forcing a rigid object through it leads to complex crack morphologies; a single oscillatory crack arises when a tool is driven laterally through a brittle sheet, while two diverging cracks and a series of concertina-like folds forms when a tool is forced laterally through a ductile sheet. On the other hand, forcing an object perpendicularly through the sheet leads to radial petal-like tears in both ductile and brittle materials. To understand these different regimes we use a combination of experiments, simulations and simple theories. In particular, we describe the transition from brittle oscillatory tearing via a single crack to ductile concertina tearing with two tears by deriving laws that describe the crack paths and wavelength of the concertina folds and provide a simple phase diagram for the morphologies in terms of the material properties of the sheet and the relative size of the tool.

preprint2014arXiv

How a Blister Heals

We use experiments to study the dynamics of the healing of a blister, a localized bump in a thin elastic layer that is adhered to a soft substrate everywhere except at the bump. We create a blister by gently placing a glass cover slip on a PDMS substrate. The pressure jump across the elastic layer drives fluid flow through micro-channels that form at the interface between the layer and the substrate; these channels coalesce at discrete locations as the blister heals and eventually disappear at a lower critical radius. The spacing of the channel follows a simple scaling law that can be theoretically justified, and the kinetics of healing is rate limited by fluid flow, but with a non-trivial dependence on the substrate thickness that likely arises due to channelization. Our study is relevant to a variety of soft adhesion scenarios.

preprint2014arXiv

How the velvet worm squirts slime

The rapid squirt of a proteinaceous slime jet endows the ancient velvet worms (Onychophora) with a unique mechanism for defense from predators and for capturing prey by entangling them in a disordered web that immobilizes their target. However, to date neither qualitative nor quantitative descriptions have been provided for this unique adaptation. Here we investigate the fast oscillatory motion of the oral papillae and the exiting liquid jet that oscillates with frequencies $f\sim 30-60$ Hz. Using anatomical images, high speed videography, theoretical analysis and a physical simulacrum we show that this fast oscillatory motion is the result of an elastohydrodynamic instability driven by the interplay between the elasticity of oral papillae and the fast unsteady flow during squirting. Our results demonstrate how passive strategies can be cleverly harnessed by organisms, while suggesting future oscillating micro-fluidic devices as well as novel ways for micro and nano fiber production using bioinspired strategies.

preprint2014arXiv

Models for elastic shells with incompatible strains

The three-dimensional shapes of thin lamina such as leaves, flowers, feathers, wings etc, are driven by the differential strain induced by the relative growth. The growth takes place through variations in the Riemannian metric, given on the thin sheet as a function of location in the central plane and also across its thickness. The shape is then a consequence of elastic energy minimization on the frustrated geometrical object. Here we provide a rigorous derivation of the asymptotic theories for shapes of residually strained thin lamina with nontrivial curvatures, i.e. growing elastic shells in both the weakly and strongly curved regimes, generalizing earlier results for the growth of nominally flat plates. The different theories are distinguished by the scaling of the mid-surface curvature relative to the inverse thickness and growth strain, and also allow us to generalize the classical Föppl-von Kármán energy to theories of prestrained shallow shells.

preprint2014arXiv

Neuromimetic Circuits with Synaptic Devices based on Strongly Correlated Electron Systems

Strongly correlated electron systems such as the rare-earth nickelates (RNiO3, R = rare-earth element) can exhibit synapse-like continuous long term potentiation and depression when gated with ionic liquids; exploiting the extreme sensitivity of coupled charge, spin, orbital, and lattice degrees of freedom to stoichiometry. We present experimental real-time, device-level classical conditioning and unlearning using nickelate-based synaptic devices in an electronic circuit compatible with both excitatory and inhibitory neurons. We establish a physical model for the device behavior based on electric-field driven coupled ionic-electronic diffusion that can be utilized for design of more complex systems. We use the model to simulate a variety of associate and non-associative learning mechanisms, as well as a feedforward recurrent network for storing memory. Our circuit intuitively parallels biological neural architectures, and it can be readily generalized to other forms of cellular learning and extinction. The simulation of neural function with electronic device analogues may provide insight into biological processes such as decision making, learning and adaptation, while facilitating advanced parallel information processing in hardware.

preprint2014arXiv

Optimal control of plate shape with incompatible strain fields

A flat plate can bend into a curved surface if it experiences an inhomogeneous growth field. In this article a method is described that numerically determines the optimal growth field giving rise to an arbitrary target shape, optimizing for closeness to the target shape and for growth field smoothness. Numerical solutions are presented, for the full non-symmetric case as well as for simplified one-dimensional and axisymmetric geometries. This system can also be solved semi-analytically by positing an ansatz for the deformation and growth fields in a circular disk with given thickness profile. Paraboloidal, cylindrical and saddle-shaped target shapes are presented as examples, of which the last two exemplify a soft mode arising from a non-axisymmetric deformation of a structure with axisymmetric material properties.

preprint2014arXiv

Surface sulci in squeezed soft solids

The squeezing of soft solids, the constrained growth of biological tissues, and the swelling of soft elastic solids such as gels can generate large compressive stresses at their surfaces. This causes the otherwise smooth surface of such a solid to becomes unstable when its stress exceeds a critical value. Previous analyses of the surface instability have assumed two-dimensional plane-strain conditions, but in experiments isotropic stresses often lead to complex three-dimensional sulcification patterns. Here we show how such diverse morphologies arise by numerically modeling the lateral compression of a rigidly clamped elastic layer. For incompressible solids, close to the instability threshold, sulci appear as I-shaped lines aligned orthogonally with their neighbors; at higher compressions they are Y-shaped and prefer a hexagonal arrangement. In contrast, highly compressible solids when squeezed show only one sulcified phase characterized by a hexagonal sulcus network.

preprint2014arXiv

Synchronization in a stochastic Hebbian network of phase oscillators

We investigate a generalized Kuramoto phase-oscillator model with Hebb-like couplings that evolve according to a stochastic differential equation on various topologies. Numerical simulations show that even with identical oscillators, there is a regime in the nearest-neighbor coupling topologies and a complex network topology where oscillators move between an in phase and anti-phase state. Phase diagrams show the transition probabilities as a function of the noise strength and rate of evolution of network coupling. A minimal theoretical model allows us to understand these transitions.

preprint2014arXiv

The dynamics of sperm cooperation in a competitive environment

Sperm cooperation has evolved in a variety of taxa and is often considered a response to sperm competition, yet the benefit of this form of collective movement remains unclear. Here we use fine-scale imaging and a minimal mathematical model to study sperm aggregation in the rodent genus $Peromyscus$. We demonstrate that as the number of sperm cells in an aggregate increase, the group moves with more persistent linearity but without increasing speed; this benefit, however, is offset in larger aggregates as the geometry of the group forces sperm to swim against one another. The result is a non-monotonic relationship between aggregate size and average velocity with both a theoretically predicted and empirically observed optimum of 6-7 sperm/aggregate. To understand the role of sexual selection in driving these sperm group dynamics, we compared two sister-species with divergent mating systems and find that sperm of $P.\,maniculatus$ (highly promiscuous), which have evolved under intense competition, form optimal-sized aggregates more often than sperm of $P.\,polionotus$ (strictly monogamous), which lack competition. Our combined mathematical and experimental study of coordinated sperm movement reveals the importance of geometry, motion and group size on sperm velocity and suggests how these physical variables interact with evolutionary selective pressures to regulate cooperation in competitive environments.

preprint2013arXiv

Collective thermoregulation in bee clusters

Swarming is an essential part of honeybee behaviour, wherein thousands of bees cling onto each other to form a dense cluster that may be exposed to the environment for several days. This cluster has the ability to maintain its core temperature actively without a central controller, and raises the question of how this is achieved. We suggest that the swarm cluster is akin to an active porous structure whose functional requirement is to adjust to outside conditions by varying its porosity to control its core temperature. Using a continuum model that takes the form of a set of advection-diffusion equations for heat transfer in a mobile porous medium, we show that the equalization of an effective "behavioural pressure", which propagates information about the ambient temperature through variations in density, leads to effective thermoregulation. Our model extends and generalizes previous models by focusing the question of mechanism on the form and role of the behavioural pressure, and allows us to explain the vertical asymmetry of the cluster (as a consequence of buoyancy driven flows), the ability of the cluster to overpack at low ambient temperatures without breaking up at high ambient temperatures, and the relative insensitivity to large variations in the ambient temperature. Finally, our theory makes testable hypotheses for how the cluster bee density should respond to externally imposed temperature inhomogeneities, and suggests strategies for biomimetic thermoregulation.

preprint2013arXiv

Digital instability of a confined elastic meniscus

Thin soft elastic layers serving as joints between relatively rigid bodies may function as sealants, thermal, electrical, or mechanical insulators, bearings, or adhesives. When such a joint is stressed, even though perfect adhesion is maintained, the exposed free meniscus in the thin elastic layer becomes unstable, leading to the formation of spatially periodic digits of air that invade the elastic layer, reminiscent of viscous fingering in a thin fluid layer. How- ever, the elastic instability is reversible and rate-independent, dis- appearing when the joint is unstressed. We use theory, experiments, and numerical simulations to show that the transition to the digital state is sudden (first-order), the wavelength and amplitude of the fingers are proportional to the thickness of the elastic layer, and the required separation to trigger the instability is inversely proportional to the in-plane dimension of the layer. Our study reveals the energetic origin of this instability and has implications for the strength of polymeric adhesives; it also suggests a method for patterning thin films reversibly with any arrangement of localized fingers in a digital elastic memory, which we confirm experimentally.

preprint2013arXiv

Elastic Platonic Shells

On microscopic scales, the crystallinity of flexible tethered or cross linked membranes determines their mechanical response. We show that by controlling the type, number and distribution of defects on a spherical elastic shell, it is possible to direct the morphology of these structures. Our numerical simulations show that by deflating a crystalline shell with defects, we can create elastic shell analogs of the classical Platonic solids. These morphologies arise via a sharp buckling transition from the sphere which is strongly hysteretic in loading-unloading. We construct a minimal Landau theory for the transition using quadratic and cubic invariants of the spherical harmonic modes. Our approach suggests methods to engineer shape into soft spherical shells using a frozen defect topology.

preprint2013arXiv

Increased Network Interdependency Leads to Aging

Although species longevity is subject to a diverse range of selective forces, the mortality curves of a wide variety of organisms are rather similar. We argue that aging and its universal characteristics may have evolved by means of a gradual increase in the systemic interdependence between a large collection of biochemical or mechanical components. Modeling the organism as a dependency network which we create using a constructive evolutionary process, we age it by allowing nodes to be broken or repaired according to a probabilistic algorithm that accounts for random failures/repairs and dependencies. Our simulations show that the network slowly accumulates damage and then catastrophically collapses. We use our simulations to fit experimental data for the time dependent mortality rates of a variety of multicellular organisms and even complex machines such as automobiles. Our study suggests that aging is an emergent finite-size effect in networks with dynamical dependencies and that the qualitative and quantitative features of aging are not sensitively dependent on the details of system structure.

preprint2013arXiv

Statistical Mechanics of Multistability in Microscopic Shells

Unlike macroscopic multistable mechanical systems such as snap bracelets or elastic shells that must be physically manipulated into various conformations, microscopic systems can undergo spontaneous conformation switching between multistable states due to thermal fluctuations. Here we investigate the statistical mechanics of shape transitions in small elastic elliptical plates and shells driven by noise. By assuming that the effects of edges are small, which we justify exactly for plates and shells with a lenticular section, we decompose the shapes into a few geometric modes whose dynamics are easy to follow. We use Monte Carlo simulations to characterize the shape transitions between conformation minima as a function of noise strength, and corroborate our results using a Fokker-Planck formalism to study the stationary distribution and the mean first passage time problem. Our results are applicable to objects such as such as graphene flakes or protein β-sheets, where fluctuations, geometry and finite size effects are important.

preprint2013arXiv

Swarming, swirling and stasis in sequestered bristle-bots

The collective ability of organisms to move coherently in space and time is ubiquitous in any group of autonomous agents that can move and sense each other and the environment. Here we investigate the origin of collective motion and its loss using macroscopic self-propelled Bristle-Bots, simple automata made from a toothbrush and powered by an onboard cell phone vibrator-motor, that can sense each other through shape-dependent local interactions, and can also sense the environment non-locally via the effects of confinement and substrate topography. We show that when Bristle-Bots are confined to a limited arena with a soft boundary, increasing the density drives a transition from a disordered and uncoordinated motion to organized collective motion either as a swirling cluster or a collective dynamical stasis. This transition is regulated by a single parameter, the relative magnitude of spinning and walking in a single automaton. We explain this using quantitative experiments and simulations that emphasize the role of the agent shape, environment, and confinement via boundaries. Our study shows how the behavioral repertoire of these physically interacting automatons controlled by one parameter translates into the mechanical intelligence of swarms.

preprint2012arXiv

A simple model for nanofiber formation by rotary jet-spinning

Nanofibers are microstructured materials that span a broad range of applications from tissue engineering scaffolds to polymer transistors. An efficient method of nanofiber production is Rotary Jet-Spinning (RJS), consisting of a perforated reservoir rotating at high speeds along its axis of symmetry, which propels a liquid, polymeric jet out of the reservoir orifice. These jets undergo stretching and solidify forming nanoscale fibers. We report a minimal scaling framework complemented by a semi-analytic and numerical approach to characterize the regimes of nanofiber production using RJS. Our theoretical model is validated for the fiber radius as a function of experimentally tunable parameters. We summarize our findings in a phase diagram for the design space of continuous nanofibers as a function of process parameters, in good agreement with experiments and with natural implications for the production rates as well as in the morphological quality of fibers.

preprint2012arXiv

Flagellar dynamics of a connected chain of active, Brownian particles

Eukaryotic flagella are active structures with a complex architecture of microtubules, motor proteins and elastic links. They are capable of whiplike motions driven by motors sliding along filaments that are themselves constrained at an end. Here, we show that active, self-propelled particles that are connected together to form a single chain that is anchored at one end can produce the graceful beating motions of flagella. We use a combination of numerical simulations, scaling analysis and mean field continuum elastic theory to demarcate the phase diagram for this type of oscillatory motion as a function of the filament length, passive elasticity, propulsion force and longitudinal persistence of propulsion directions. Depending on the nature of the anchoring, we show that filament either undergoes flagella-like beating or assumes a steadily rotating coiled conformation. Our system is simpler than its biological inspiration, and thus could be experimentally realized using a variety of self-propelled particles.

preprint2012arXiv

Geometric Mechanics of Curved Crease Origami

Folding a sheet of paper along a curve can lead to structures seen in decorative art and utilitarian packing boxes. Here we present a theory for the simplest such structure: an annular circular strip that is folded along a central circular curve to form a three-dimensional buckled structure driven by geometrical frustration. We quantify this shape in terms of the radius of the circle, the dihedral angle of the fold and the mechanical properties of the sheet of paper and the fold itself. When the sheet is isometrically deformed everywhere except along the fold itself, stiff folds result in creases with constant curvature and oscillatory torsion. However, relatively softer folds inherit the broken symmetry of the buckled shape with oscillatory curvature and torsion. Our asymptotic analysis of the isometrically deformed state is corroborated by numerical simulations which allow us to generalize our analysis to study multiply folded structures.

preprint2012arXiv

Geometric Mechanics of Periodic Pleated Origami

Origami is the archetype of a structural material with unusual mechanical properties that arise almost exclusively from the geometry of its constituent folds and forms the basis for mechanical metamaterials with an extreme deformation response. Here we consider a simple periodically folded structure Miura-ori, which is composed of identical unit cells of mountain and valley folds with four-coordinated ridges, defined completely by 2 angles and 2 lengths. We use the geometrical properties of a Miura-ori plate to characterize its elastic response to planar and non-planar piece- wise isometric deformations and calculate the two-dimensional stretching and bending response of a Miura-ori sheet, and show that the in-plane and out-of-plane Poisson's ratios are equal in magnitude, but opposite in sign. Our geometric approach also allows us to solve the inverse design problem of determining the geometric parameters that achieve the optimal geometric and mechanical response of such structures.

preprint2012arXiv

Macroscopic magnetic frustration

Although geometrical frustration transcends scale, it has primarily been evoked in the micro and mesoscopic realm to characterize such phases as spin-ice liquids and glasses and to explain the behavior of such materials as multiferroics, high temperature superconductors, colloids and copolymers. Here we introduce a system of macroscopic ferromagnetic rotors arranged in a planar lattice capable of out-of-plane movement that exhibit the characteristic honeycomb spin ice rules studied and seen so far only in its mesoscopic manifestation. We find that a polarized initial state of this system settles into the honeycomb spin ice phase with relaxation on multiple time scales. We explain this relaxation process using a minimal classical mechanical model which includes Coulombic interactions between magnetic charges located at the ends of the magnets and viscous dissipation at the hinges. Our study shows how macroscopic frustration arises in a purely classical setting that is amenable to experiment, easy manipulation, theory and computation, and shows phenomena that are not visible in their microscopic counterparts.

preprint2012arXiv

Minimal surfaces bounded by elastic lines

In mathematics, the classical Plateau problem consists of finding the surface of least area that spans a given rigid boundary curve. A physical realization of the problem is obtained by dipping a stiff wire frame of some given shape in soapy water and then removing it; the shape of the spanning soap film is a solution to the Plateau problem. But what happens if a soap film spans a loop of inextensible but flexible wire? We consider this simple query that couples Plateau's problem to Euler's Elastica: a special class of twist-free curves of given length that minimize their total squared curvature energy. The natural marriage of two of the oldest geometrical problems linking physics and mathematics leads to a quest for the shape of a minimal surface bounded by an elastic line: the Euler-Plateau problem. We use a combination of simple physical experiments with soap films that span soft filaments, scaling concepts, exact and asymptotic analysis combined with numerical simulations to explore some of the richness of the shapes that result. Our study raises questions of intrinsic interest in geometry and its natural links to a range of disciplines including materials science, polymer physics, architecture and even art.

preprint2012arXiv

Scale and Nature of Sulcification Patterns

Sulci are surface folds commonly seen in strained soft elastomers and form via a strongly subsubcriticalcritical, yet scale-free instability. Treating the threshold for nonlinear instability as a nonlinear critical point, we explain the nature of sulcus patterns in terms of the scale and translation symmetries which are broken by the formation of an isolated, small sulcus. Our perturbative theory and simulations show that sulcus formation in a thick, compressed slab can arise either as a supercritical or as a weakly subcritical bifurcation relative to this nonlinear critical point, depending on the boundary conditions. An infinite number of competing, equilibrium patterns simultaneously emerge at this critical point, but the one selected has the lowest energy. We give a simple, physical explanation for the formation of these sulcification patterns using an analogy to a solid-solid phase transition with a finite energy of transformation.

preprint2012arXiv

The size, shape, and dynamics of cellular blebs

A cellular bleb grows when a portion of the cell membrane detaches from the underlying cortex under the influence of a cytoplasmic pressure. We develop a quantitative model for the growth and dynamics of these objects in a simple two-dimensional setting. In particular, we first find the minimum cytoplasmic pressure and minimum unsupported membrane length for a stationary bleb to nucleate and grow as a function of the membrane-cortex adhesion. We next show how a bleb may travel around the periphery of the cell when the cytoplasmic pressure varies in space and time in a prescribed way and find that the traveling speed is governed by the speed of the pressure change induced by local cortical contraction while the shape of the traveling bleb is governed by the speed of cortical healing. Finally, we relax the assumption that the pressure change is prescribed and couple it hydrodynamically to the cortical contraction and membrane deformation. By quantifying the phase space of bleb formation and dynamics, our framework serves to delineate the range and scope of bleb-associated cell motility and synthesizes a variety of experimental observations.

preprint2012arXiv

Twisting Graphene Nanoribbons into Carbon Nanotubes

Although carbon nanotubes consist of honeycomb carbon, they have never been fabricated from graphene directly. Here, it is shown by quantum molecular-dynamics simulations and classical continuum-elasticity modeling, that graphene nanoribbons can, indeed, be transformed into carbon nanotubes by means of twisting. The chiralities of the tubes thus fabricated can be not only predicted but also externally controlled. This twisting route is an opportunity for nanofabrication, and is easily generalizable to ribbons made of other planar nanomaterials.

preprint2011arXiv

Banding, Excitability and Chaos in Active Nematic Suspensions

Motivated by the observation of highly unstable flowing states in suspensions of microtubules and kinesin, we analyze a model of mutually-propelled filaments suspended in a solvent. The system undergoes a mean-field isotropic-nematic transition for large enough filament concentrations when the nematic order parameter is allowed to vary in space and time. We analyze the model in two contexts: a quasi-one-dimensional channel with no-slip walls and a two-dimensional box with periodic boundaries. Using stability analysis and numerical calculations we show that the interplay between non-uniform nematic order, activity, and flow results in a variety of complex scenarios that include spontaneous banded laminar flow, relaxation oscillations, and chaos.

preprint2011arXiv

Buckling instability of a thin-layer rectilinear Couette flow

We analyse the buckling stability of a thin, viscous sheet when subject to simple shear, providing conditions for the onset of the dominant out-of-plane modes using two models: (i) an asymptotic theory for the dynamics of a viscous plate and (ii) the full Stokes equations. In either case, the plate is stabilised by a combination of viscous resistance, surface tension and buoyancy relative to an underlying denser fluid. In the limit of vanishing thickness, plates buckle at a shear rate $γ/(μd)$ independent of buoyancy, where 2d is the plate thickness, $γ$ is the average surface tension between the upper and lower surfaces and $μ$ is the fluid viscosity. For thicker plates stabilised by an equal surface tension at the upper and lower surfaces, at and above onset, the most unstable mode has moderate wavelength, is stationary in the frame of the centre-line, spans the width of the plate with crests and troughs aligned at approximately $45^\circ$ to the walls and closely resembles elastic shear modes. The thickest plates that can buckle have an aspect ratio (thickness/width) approximately $\EC$ and are stabilised only by internal viscous resistance. We show that the viscous plate model can only accurately describe the onset of buckling for vanishingly thin plates but provides an excellent description of the most unstable mode above onset. Finally, we show that by modifying the plate model to incorporate advection and make the model material frame-invariant, it is possible to extend its predictive power to describe relatively short, travelling waves.

preprint2011arXiv

Dynamic instability of a growing adsorbed polymorphic filament

The intermittent transition between slow growth and rapid shrinkage in polymeric assemblies is termed dynamic instability, a feature observed in a variety of biochemically distinct assemblies including microtubules, actin and their bacterial analogs. The existence of this labile phase of a polymer has many functional consequences in cytoskeletal dynamics, and its repeated appearance suggests that it is relatively easy to evolve. Here, we consider the minimal ingredients for the existence of dynamic instability by considering a single polymorphic filament that grows by binding to a substrate, undergoes a conformation change, and may unbind as a consequence of the residual strains induced by this change. We identify two parameters that control the phase space of possibilities for the filament: a structural mechanical parameter that characterizes the ratio of the bond strengths along the filament to those with the substrate (or equivalently the ratio of longitudinal to lateral interactions in an assembly), and a kinetic parameter that characterizes the ratio of time scales for growth and conformation change. In the deterministic limit, these parameters serve to demarcate a region of uninterrupted growth from that of collapse. However, in the presence of disorder in either the structural or the kinetic parameter the growth and collapse phases can coexist where the filament can grow slowly, shrink rapidly, and transition between these phases, thus exhibiting dynamic instability. We exhibit the window for the existence of dynamic instability in a phase diagram that allows us to quantify the evolvability of this labile phase.

preprint2011arXiv

Multistability of free spontaneously-curved anisotropic strips

Multistable structures are objects with more than one stable conformation, exemplified by the simple switch. Continuum versions are often elastic composite plates or shells, such as the common measuring tape or the slap bracelet, both of which exhibit two stable configurations: rolled and unrolled. Here we consider the energy landscape of a general class of multistable anisotropic strips with spontaneous Gaussian curvature. We show that while strips with non-zero Gaussian curvature can be bistable, strips with positive spontaneous curvature are always bistable, independent of the elastic moduli, strips of spontaneous negative curvature are bistable only in the presence of spontaneous twist and when certain conditions on the relative stiffness of the strip in tension and shear are satisfied. Furthermore, anisotropic strips can become tristable when their bending rigidity is small. Our study complements and extends the theory of multistability in anisotropic shells and suggests new design criteria for these structures.

preprint2011arXiv

Probability and dynamics in the toss of a non-bouncing thick coin

When a thick cylindrical coin is tossed in the air and lands without bouncing on an inelastic substrate, it ends up on its face or its side. We account for the rigid body dynamics of spin and precession and calculate the probability distribution of heads, tails, and sides for a thick coin as a function of its dimensions and the distribution of its initial conditions. Our theory yields a simple expression for the aspect ratio of homogeneous coins with a prescribed frequency of heads/tails compared to sides, which we validate by tossing experiments using coins of different aspect ratios.

preprint2011arXiv

Reply to Comment on "Statistical mechanics of developable ribbons"

Starostin and van der Heijden [arXiv:1111.2029] question the validity of the results reported in our Letter [Phys. Rev. Lett. 104, 238104 (2010)] regarding the existence of spontaneous helical structure in finite-temperature developable ribbons. According to the authors of the Comment, our setting is affected by a systematic error in the calculation of the torsion τ that results in a bias of the Monte Carlo sampling toward helical configurations, and originates from the fact that the Frenet frame, used by us to describe the centerline of the ribbons, rotates by 180^{\circ} about the direction of the tangent vector every time an inflection point (i.e. a point where the curvature κ changes in sign) is encountered. While the rotation of the frame at infection points is a well known property of Frenet frame, the conclusions of Starostin and van der Heijden on how this phenomenon would affect the Monte Carlo sampling are incorrect, as are many of the other statements that they make.

preprint2011arXiv

Shock driven jamming and periodic fracture of particulate rafts

A tenuous monolayer of hydrophobic particles at the air-water interface often forms a scum or raft. When such a monolayer is disturbed by the localized introduction of a surfactant droplet, a radially divergent surfactant shock front emanates from the surfactant origin and packs the particles into a jammed, compact, annular band with a packing fraction that saturates at a peak packing fraction $ϕ^*$. As the resulting two-dimensional, disordered elastic band grows with time and is driven radially outwards by the surfactant, it fractures to form periodic triangular cracks with robust geometrical features. We find the number of cracks $N$ and the compaction band radius $R^*$ at fracture onset vary monotonically with the initial packing fraction ($ϕ_{init}$). However, its width $W^*$ is constant for all $ϕ_{init}$. A simple geometric theory that treats the compaction band as an elastic annulus, and accounts for mass conservation allows us to deduce that $N \simeq 2πR^*/W^* \simeq 4πϕ_{RCP}/ϕ_{init}$, a result we verify both experimentally and numerically. We show the essential ingredients for this phenomenon are an initially low enough particulate packing fraction that allows surfactant driven advection to cause passive jamming and eventual fracture of the hydrophobic particulate interface.

preprint2011arXiv

Skating on a Film of Air: Drops Impacting on a Surface

Drops impacting on a surface are ubiquitous in our everyday experience. This impact is understood within a commonly accepted hydrodynamic picture: it is initiated by a rapid shock and a subsequent ejection of a sheet leading to beautiful splashing patterns. However, this picture ignores the essential role of the air that is trapped between the impacting drop and the surface. Here we describe a new imaging modality that is sensitive to the behavior right at the surface. We show that a very thin film of air, only a few tens of nanometers thick, remains trapped between the falling drop and the surface as the drop spreads. The thin film of air serves to lubricate the drop enabling the fluid to skate on the air film laterally outward at surprisingly high velocities, consistent with theoretical predictions. Eventually this thin film of air must break down as the fluid wets the surface. We suggest that this occurs in a spinodal-like fashion, and causes a very rapid spreading of a wetting front outwards; simultaneously the wetting fluid spreads inward much more slowly, trapping a bubble of air within the drop. Our results show that the dynamics of impacting drops are much more complex than previously thought and exhibit a rich array of unexpected phenomena that require rethinking classical paradigms.

preprint2011arXiv

Unfolding the Sulcus

Sulci are localized furrows on the surface of soft materials that form by a compression-induced instability. We unfold this instability by breaking its natural scale and translation invariance, and compute a limiting bifurcation diagram for sulcfication showing that it is a scale-free, sub-critical {\em nonlinear} instability. In contrast with classical nucleation, sulcification is {\em continuous}, occurs in purely elastic continua and is structurally stable in the limit of vanishing surface energy. During loading, a sulcus nucleates at a point with an upper critical strain and an essential singularity in the linearized spectrum. On unloading, it quasi-statically shrinks to a point with a lower critical strain, explained by breaking of scale symmetry. At intermediate strains the system is linearly stable but nonlinearly unstable with {\em no} energy barrier. Simple experiments confirm the existence of these two critical strains.

preprint2010arXiv

Excitable Patterns in Active Nematics

We analyze a model of mutually-propelled filaments suspended in a two-dimensional solvent. The system undergoes a mean-field isotropic-nematic transition for large enough filament concentrations and the nematic order parameter is allowed to vary in space and time. We show that the interplay between non-uniform nematic order, activity and flow results in spatially modulated relaxation oscillations, similar to those seen in excitable media. In this regime the dynamics consists of nearly stationary periods separated by "bursts" of activity in which the system is elastically distorted and solvent is pumped throughout. At even higher activity the dynamics becomes chaotic.

preprint2010arXiv

Flow-induced channelization in a porous medium

We propose a theory for erosional channelization induced by fluid flow in a saturated granular porous medium. When the local fluid flow-induced stress is larger than a critical threshold, grains are dislodged and carried away so that the porosity of the medium is altered by erosion. This in turn affects the local hydraulic conductivity and pressure in the medium and results in the growth and development of channels that preferentially conduct the flow. Our multiphase model involves a dynamical porosity field that evolves along with the volume fraction of the mobile and immobile grains in response to fluid flow that couples the spatiotemporal dynamics of the three phases. Numerical solutions of the resulting initial boundary value problem show how channels form in porous media and highlights how heterogeneity in the erosion threshold dictates the form of the patterns and thus the ability to control them.

preprint2010arXiv

Generalized Erdos Numbers

We propose a simple real-valued generalization of the well known integer-valued Erdos number as a topological, non-metric measure of the `closeness' felt between two nodes in an undirected, weighted graph. These real-valued Erdos numbers are asymmetric and are able to distinguish between network topologies that standard distance metrics view as identical. We use this measure to study some simple analytically tractable networks, and show the utility of our measure to devise a ratings scheme based on the generalized Erdos number that we deploy on the data from the NetFlix prize, and find a significant improvement in our ratings prediction over a baseline.

preprint2010arXiv

Geometry, mechanics and electronics of singular structures and wrinkles in graphene

As the thinnest atomic membrane, graphene presents an opportunity to combine geometry, elasticity and electronics at the limits of their validity. The availability of reliable atomistic potentials for graphene allows unprecedented precise simulations of the mechanical response of atomic membranes. Here we describe the transport and electronic structure in the neighbourhood of conical singularities, the elementary excitations of the ubiquitous wrinkled and crumpled graphene that occur in non-epitaxial suspended samples where shear stresses are unavoidable. We use a combination of atomistic mechanical simulations, analytical geometry and transport calculations in curved graphene, and exact diagonalization of the electronic spectrum to calculate the effects of geometry on electronic structure, transport and mobility in suspended samples. We also point out how the geometry-generated pseudo-magnetic/electric fields might disrupt Landau quantization under a magnetic field.

preprint2010arXiv

Statistical Mechanics of Developable Ribbons

We investigate the statistical mechanics of long developable ribbons of finite width and very small thickness. The constraint of isometric deformations in these ribbon-like structures that follows from the geometric separation of scales introduces a coupling between bending and torsional degrees of freedom. Using analytical techniques and Monte Carlo simulations, we find that the tangent-tangent correlation functions always exhibits an oscillatory decay at any finite temperature implying the existence of an underlying helical structure even in absence of a preferential zero-temperature twist. In addition the persistence length is found to be over three times larger than that of a wormlike chain having the same bending rigidity. Our results are applicable to many ribbon-like objects in polymer physics and nanoscience that are not described by the classical worm-like chain model.

preprint2010arXiv

The von Karman equations for plates with residual strain

We provide a derivation of the Foppl-von Karman equations for the shape of and stresses in an elastic plate with residual strains. These might arise from a range of causes: inhomogeneous growth, plastic deformation, swelling or shrinkage driven by solvent absorption. Our analysis gives rigorous bounds on the convergence of the three dimensional equations of elasticity to the low-dimensional description embodied in the plate-like description of laminae and thus justifies a recent formulation of the problem to the shape of growing leaves. It also formalizes a procedure that can be used to derive other low-dimensional descriptions of active materials.

preprint2009arXiv

The Shape and Motion of a Ruck in a Rug

The motion of a ruck in a rug is used as an analogy to explain the role of dislocations in the deformation of crystalline solids. We take the analogy literally and study the shape and motion of a bump, wrinkle or ruck in a thin sheet in partial contact with a rough substrate in a gravitational field. Using a combination of experiments, scaling analysis and numerical solutions of the governing equations, we first quantify the static shape of a ruck on a horizontal plane. When the plane is inclined, the ruck becomes asymmetric and moves by rolling only when the the inclination of the plane reaches a critical angle. We find that the angle at which this first occurs is larger than the angle at which the ruck stops, i.e. static rolling friction is larger than dynamic rolling friction. Once the ruck is in motion, it travels at a constant speed proportional to the sine of the angle of inclination, a result that we rationalize in terms of a simple power balance. We conclude with a simple implication of our study for the onset of rolling motion at soft interfaces.

preprint2008arXiv

A coordination-based approach to elasticity of floppy and stiff random networks

We study the role of connectivity on the linear and nonlinear elastic behavior of amorphous systems using a two-dimensional random network of harmonic springs as a model system. A natural characterization of these systems arises in terms of the network coordination relative to that of an isostatic network $δz$; a floppy network has $δz<0$, while a stiff network has $δz>0$. Under the influence of an externally applied load we observe that the response of both floppy and rigid network are controlled by the same critical point, corresponding to the onset of rigidity. We use numerical simulations to compute the exponents which characterize the shear modulus, the amplitude of non-affine displacements, and the network stiffening as a function of $δz$, derive these theoretically and make predictions for the mechanical response of glasses and fibrous networks.

preprint2006arXiv

Curvature condensation and bifurcation in an elastic shell

We study the formation and evolution of localized geometrical defects in an indented cylindrical elastic shell using a combination of experiment and numerical simulation. We find that as a symmetric localized indentation on a semi-cylindrical shell increases, there is a transition from a global mode of deformation to a localized one which leads to the condensation of curvature along a symmetric parabolic crease. This process introduces a soft mode in the system, converting a load-bearing structure into a hinged, kinematic mechanism. Further indentation leads to twinning wherein the parabolic crease bifurcates into two creases that move apart on either side of the line of symmetry. A qualitative theory captures the main features of the phenomena and leads to sharper questions about the nucleation of these defects.

preprint2006arXiv

Mechanics of Interfacial Composite Materials

Recent experiments and simulations have demonstrated that particle-covered interfaces can exist in stable non-spherical shapes as a result of the steric jamming of the interfacially trapped particles, which confers the interface with solid-like properties. We provide an experimental and theoretical characterization of the mechanical properties of these armored objects, with attention given to the two-dimensional granular state of the interface. Small inhomogeneous stresses produce a plastic response while homogeneous stresses produce a weak elastic response. Shear-driven particle-scale rearrangements explain the basic threshold needed to obtain the near-perfect plastic deformation that is observed. Furthermore, the inhomogeneous stress state of the interface is exhibited experimentally by using surfactants to destabilize the particles on the surface. Since the interfacially trapped particles retain their individual characteristics, armored interfaces can be recognized as a kind of composite material with distinct chemical, structural and mechanical properties.

preprint2006arXiv

Persistence of a pinch in a pipe

The response of low-dimensional solid objects combines geometry and physics in unusual ways, exemplified in structures of great utility such as a thin-walled tube that is ubiquitous in nature and technology. Here we provide a particularly surprising consequence of this confluence of geometry and physics in tubular structures: the anomalously large persistence of a localized pinch in an elastic pipe whose effect decays very slowly as an oscillatory exponential with a persistence length that diverges as the thickness of the tube vanishes, which we confirm experimentally. The result is more a consequence of geometry than material properties, and is thus equally applicable to carbon nanotubes as it is to oil pipelines.

preprint2004arXiv

Dynamics of poroelastic filaments

We investigate the stability and geometrically non-linear dynamics of slender rods made of a linear isotropic poroelastic material. Dimensional reduction leads to the evolution equation for the shape of the poroelastica where, in addition to the usual terms for the bending of an elastic rod, we find a term that arises from fluid-solid interaction. Using the poroelastica equation as a starting point, we consider the load controlled and displacement controlled planar buckling of a slender rod, as well as the closely related instabilities of a rod subject to twisting moments and compression when embedded in an elastic medium. This work has applications to the active and passive mechanics of thin filaments and sheets made from gels, plant organs such as stems, roots and leaves, sponges, cartilage layers and bones.

preprint2000arXiv

Rippling Instability of a Collapsing Bubble

When a bubble of air rises to the top of a highly viscous liquid, it forms a dome-shaped protuberance on the free surface. Unlike a soap bubble, it bursts so slowly as to collapse under its own weight simultaneously, and folds into a striking wavy structure. This rippling effect occurs in fact for both elastic and viscous sheets, and a theory for its onset is formulated. The growth of the corrugation is governed by the competition between gravitational and bending (shearing) forces and is exhibited for a range of densities, stiffnesses (viscosities), and sizes -- a result which arises less from dynamics than from geometry, suggesting a wide validity. A quantitative expression for the number of ripples is presented, together with experimental results which are in agreement with the theoretical predictions.

L. Mahadevan

What is connected

Connect this record

See the researcher in context

Building this map preview

83 published item(s)

Stigmergic optimal transport

How to grow a flat leaf

Modular representation and control of floppy networks

Quasicrystal kirigami

Active nematic defects and epithelial morphogenesis

Combing a double helix

Geometric control of topological dynamics in a singing saw

Geometrical dynamics of edge-driven surface growth

Instability-induced patterning of a jelling jet

Prestrain-induced contraction in 1D random elastic chains

Static adhesion hysteresis in elastic structures

Deterministic and stochastic control of kirigami topology

Models of benthic bipedalism

Optimal policies for mitigating pandemic costs

Self-excited motions of volatile drops on swellable sheets

Random sequential adsorption of spheres on a cylinder

Programming shape using kirigami tessellations

Active elastohydrodynamics of vesicles in narrow, blind constrictions

Self-sustained lift and low friction via soft lubrication

Elastic Cheerios effect: self-assembly of cylinders on a soft solid

Elastohydrodynamics and kinetics of protein patterning in the immunological synapse

Elastohydrodynamics of a sliding, spinning and sedimenting cylinder near a soft wall

Elastohydrodynamics of contact in adherent sheets

Generalized Erdos Numbers for network analysis

Gyrification from constrained cortical expansion

Optimal switching strategies for stochastic geocentric/egocentric navigation

Persistence of activity in noisy motor-filament assemblies

Similarity and singularity in adhesive elastohydrodynamic touchdown

Solid friction between soft filaments

The Organization and Control of an Evolving Interdependent Population

2D and Axisymmetric Incompressible Elastic Green's Functions

Dynamics of evaporative colloidal patterning

Elastocapillary coalescence of plates and pillars

Evaporation-driven ring and film deposition from colloidal droplets

Exactly isochoric deformations of soft solids

Feedback induced phase transitions in active porous media

Fluid driven fingering instability of a confined elastic meniscus

Forced tearing of ductile and brittle thin sheets

How a Blister Heals

How the velvet worm squirts slime

Models for elastic shells with incompatible strains

Neuromimetic Circuits with Synaptic Devices based on Strongly Correlated Electron Systems

Optimal control of plate shape with incompatible strain fields

Surface sulci in squeezed soft solids

Synchronization in a stochastic Hebbian network of phase oscillators

The dynamics of sperm cooperation in a competitive environment

Collective thermoregulation in bee clusters

Digital instability of a confined elastic meniscus

Elastic Platonic Shells

Increased Network Interdependency Leads to Aging

Statistical Mechanics of Multistability in Microscopic Shells

Swarming, swirling and stasis in sequestered bristle-bots

A simple model for nanofiber formation by rotary jet-spinning

Flagellar dynamics of a connected chain of active, Brownian particles

Geometric Mechanics of Curved Crease Origami

Geometric Mechanics of Periodic Pleated Origami

Macroscopic magnetic frustration

Minimal surfaces bounded by elastic lines

Scale and Nature of Sulcification Patterns

The size, shape, and dynamics of cellular blebs

Twisting Graphene Nanoribbons into Carbon Nanotubes

Banding, Excitability and Chaos in Active Nematic Suspensions

Buckling instability of a thin-layer rectilinear Couette flow

Dynamic instability of a growing adsorbed polymorphic filament

Multistability of free spontaneously-curved anisotropic strips

Probability and dynamics in the toss of a non-bouncing thick coin

Reply to Comment on "Statistical mechanics of developable ribbons"

Shock driven jamming and periodic fracture of particulate rafts

Skating on a Film of Air: Drops Impacting on a Surface

Unfolding the Sulcus

Excitable Patterns in Active Nematics

Flow-induced channelization in a porous medium

Generalized Erdos Numbers

Geometry, mechanics and electronics of singular structures and wrinkles in graphene