Geometric phase predicts locomotion performance in undulating living systems across scales.

Self-propelling organisms locomote via generation of patterns of self-deformation. Despite the diversity of body plans, internal actuation schemes and environments in limbless vertebrates and invertebrates, such organisms often use similar traveling waves of axial body bending for movement. Delineating how self-deformation parameters lead to locomotor performance (e.g. speed, energy, turning capabilities) remains challenging. We show that a geometric framework, replacing laborious calculation with a diagrammatic scheme, is well-suited to discovery and comparison of effective patterns of wave dynamics in diverse living systems. We focus on a regime of undulatory locomotion, that of highly damped environments, which is applicable not only to small organisms in viscous fluids, but also larger animals in frictional fluids (sand) and on frictional ground. We find that the traveling wave dynamics used by mm-scale nematode worms and cm-scale desert dwelling snakes and lizards can be described by time series of weights associated with two principal modes. The approximately circular closed path trajectories of mode weights in a self-deformation space enclose near-maximal surface integral (geometric phase) for organisms spanning two decades in body length. We hypothesize that such trajectories are targets of control (which we refer to as "serpenoid templates"). Further, the geometric approach reveals how seemingly complex behaviors such as turning in worms and sidewinding snakes can be described as modulations of templates. Thus, the use of differential geometry in the locomotion of living systems generates a common description of locomotion across taxa and provides hypotheses for neuromechanical control schemes at lower levels of organization.

Field-mediated locomotor dynamics on highly deformable surfaces.

Studies of active matter-systems consisting of individuals or ensembles of internally driven and damped locomotors-are of interest to physicists studying nonequilibrium dynamics, biologists interested in individuals and swarm locomotion, and engineers designing robot controllers. While principles governing active systems on hard ground or within fluids are well studied, another class of systems exists at deformable interfaces. Such environments can display mixes of fluid-like and elastic features, leading to locomotor dynamics that are strongly influenced by the geometry of the surface, which, in itself, can be a dynamical entity. To gain insight into principles by which locomotors are influenced via a deformation field alone (and can influence other locomotors), we study robot locomotion on an elastic membrane, which we propose as a model of active systems on highly deformable interfaces. As our active agent, we use a differential driven wheeled robotic vehicle which drives straight on flat homogeneous surfaces, but reorients in response to environmental curvature. We monitor the curvature field-mediated dynamics of a single vehicle interacting with a fixed deformation as well as multiple vehicles interacting with each other via local deformations. Single vehicles display precessing orbits in centrally deformed environments, while multiple vehicles influence each other by local deformation fields. The active nature of the system facilitates a differential geometry-inspired mathematical mapping from the vehicle dynamics to those of test particles in a fictitious "spacetime," allowing further understanding of the dynamics and how to control agent interactions to facilitate or avoid multivehicle membrane-induced cohesion.

Functional consequences of convergently evolved microscopic skin features on snake locomotion.

The small structures that decorate biological surfaces can significantly affect behavior, yet the diversity of animal-environment interactions essential for survival makes ascribing functions to structures challenging. Microscopic skin textures may be particularly important for snakes and other limbless locomotors, where substrate interactions are mediated solely through body contact. While previous studies have characterized ventral surface features of some snake species, the functional consequences of these textures are not fully understood. Here, we perform a comparative study, combining atomic force microscopy measurements with mathematical modeling to generate predictions that link microscopic textures to locomotor performance. We discover an evolutionary convergence in the ventral skin structures of a few sidewinding specialist vipers that inhabit sandy deserts-an isotropic texture that is distinct from the head-to-tail-oriented, micrometer-sized spikes observed on a phylogenetically broad sampling of nonsidewinding vipers and other snakes from diverse habitats and wide geographic range. A mathematical model that relates structural directionality to frictional anisotropy reveals that isotropy enhances movement during sidewinding, whereas anisotropy improves movement during slithering via lateral undulation of the body. Our results highlight how an integrated approach can provide quantitative predictions for structure-function relationships and insights into behavioral and evolutionary adaptations in biological systems.

Mechanical diffraction reveals the role of passive dynamics in a slithering snake.

Limbless animals like snakes inhabit most terrestrial environments, generating thrust to overcome drag on the elongate body via contacts with heterogeneities. The complex body postures of some snakes and the unknown physics of most terrestrial materials frustrates understanding of strategies for effective locomotion. As a result, little is known about how limbless animals contend with unplanned obstacle contacts. We studied a desert snake, Chionactis occipitalis, which uses a stereotyped head-to-tail traveling wave to move quickly on homogeneous sand. In laboratory experiments, we challenged snakes to move across a uniform substrate and through a regular array of force-sensitive posts. The snakes were reoriented by the array in a manner reminiscent of the matter-wave diffraction of subatomic particles. Force patterns indicated the animals did not change their self-deformation pattern to avoid or grab the posts. A model using open-loop control incorporating previously described snake muscle activation patterns and body-buckling dynamics reproduced the observed patterns, suggesting a similar control strategy may be used by the animals. Our results reveal how passive dynamics can benefit limbless locomotors by allowing robust transit in heterogeneous environments with minimal sensing.

Surprising simplicities and syntheses in limbless self-propulsion in sand.

Animals moving on and in fluids and solids move their bodies in diverse ways to generate propulsion and lift forces. In fluids, animals can wiggle, stroke, paddle or slap, whereas on hard frictional terrain, animals largely engage their appendages with the substrate to avoid slip. Granular substrates, such as desert sand, can display complex responses to animal interactions. This complexity has led to locomotor strategies that make use of fluid-like or solid-like features of this substrate, or combinations of the two. Here, we use examples from our work to demonstrate the diverse array of methods used and insights gained in the study of both surface and subsurface limbless locomotion in these habitats. Counterintuitively, these seemingly complex granular environments offer certain experimental, theoretical, robotic and computational advantages for studying terrestrial movement, with the potential for providing broad insights into morphology and locomotor control in fluids and solids, including neuromechanical control templates and morphological and behavioral evolution. In particular, granular media provide an excellent testbed for a locomotion framework called geometric mechanics, which was introduced by particle physicists and control engineers in the last century, and which allows quantitative analysis of alternative locomotor patterns and morphology to test for control templates, optimality and evolutionary alternatives. Thus, we posit that insights gained from movement in granular environments can be translated into principles that have broader applications across taxa, habitats and movement patterns, including those at microscopic scales.

Divergence of Voronoi Cell Anisotropy Vector: A Threshold-Free Characterization of Local Structure in Amorphous Materials.

Characterizing structural inhomogeneity is an essential step in understanding the mechanical response of amorphous materials. We introduce a threshold-free measure based on the field of vectors pointing from the center of each particle to the centroid of the Voronoi cell in which the particle resides. These vectors tend to point in toward regions of high free volume and away from regions of low free volume, reminiscent of sinks and sources in a vector field. We compute the local divergence of these vectors, where positive values correspond to overpacked regions and negative values identify underpacked regions within the material. Distributions of this divergence are nearly Gaussian with zero mean, allowing for structural characterization using only the moments of the distribution. We explore how the standard deviation and skewness vary with the packing fraction for simulations of bidisperse systems and find a kink in these moments that coincides with the jamming transition.

Tunable capillary-induced attraction between vertical cylinders.

Deformation of a fluid interface caused by the presence of objects at the interface can lead to large lateral forces between objects. We explore these fluid-mediated attractive force between partially submerged vertical cylinders. Forces are experimentally measured by slowly separating cylinder pairs and cylinder triplets after capillary rise is initially established for cylinders in contact. For cylinder pairs, numerical computations and a theoretical model are found to be in good agreement with measurements. The model provides insight into the relative importance of the contributions to the total force. For small separations, the lateral force is dominated by the fluid pressure acting over the wetted cylinder surfaces. At large separations, the surface tension acting along the contact line dominates the lateral force. A crossover between the two regimes occurs at a separation of around half of a capillary length. The experimentally measured forces between cylinder triplets are also in good agreement with numerical computations, and we show that pairwise contributions account for nearly all of the attractive force between triplets. For cylinders with an equilibrium capillary rise height greater than the height of the cylinder, we find that the attractive force depends on the height of the cylinders above the submersion level, which provides a means to create precisely controlled tunable cohesive forces between objects deforming a fluid interface.

Rheology of soft colloids across the onset of rigidity: scaling behavior, thermal, and non-thermal responses.

We study the rheological behavior of colloidal suspensions composed of soft sub-micron-size hydrogel particles across the liquid-solid transition. The measured stress and strain-rate data, when normalized by thermal stress and time scales, suggest our systems reside in a regime wherein thermal effects are important. In a different vein, critical point scaling predictions for the jamming transition, typical in athermal systems, are tested. Near dynamic arrest, the suspensions exhibit scaling exponents similar to those reported in Nordstrom et al., Phys. Rev. Lett., 2010, 105, 175701. The observation suggests that our system exhibits a glass transition near the onset of rigidity, but it also exhibits a jamming-like scaling further from the transition point. These observations are thought-provoking in light of recent theoretical and simulation findings, which show that suspension rheology across the full range of microgel particle experiments can exhibit both thermal and athermal mechanisms.

A general locomotion control framework for multi-legged locomotors.

Serially connected robots are promising candidates for performing tasks in confined spaces such as search and rescue in large-scale disasters. Such robots are typically limbless, and we hypothesize that the addition of limbs could improve mobility. However, a challenge in designing and controlling such devices lies in the coordination of high-dimensional redundant modules in a way that improves mobility. Here we develop a general framework to discover templates to control serially connected multi-legged robots. Specifically, we combine two approaches to build a general shape control scheme which can provide baseline patterns of self-deformation ('gaits') for effective locomotion in diverse robot morphologies. First, we take inspiration from a dimensionality reduction and a biological gait classification scheme to generate cyclic patterns of body deformation and foot lifting/lowering, which facilitate the generation of arbitrary substrate contact patterns. Second, we extend geometric mechanics, which was originally introduced to study swimming at low Reynolds numbers, to frictional environments, allowing the identification of optimal body-leg coordination in this common terradynamic regime. Our scheme allows the development of effective gaits on flat terrain with diverse numbers of limbs (4, 6, 16, and even 0 limbs) and backbone actuation. By properly coordinating the body undulation and leg placement, our framework combines the advantages of both limbless robots (modularity and narrow profile) and legged robots (mobility). Our framework can provide general control schemes for the rapid deployment of general multi-legged robots, paving the way toward machines that can traverse complex environments. In addition, we show that our framework can also offer insights into body-leg coordination in living systems, such as salamanders and centipedes, from a biomechanical perspective.

Side-impact collision: mechanics of obstacle negotiation in sidewinding snakes.

Snakes excel at moving through cluttered environments, and heterogeneities can be used as propulsive contacts for snakes performing lateral undulation. However, sidewinding, which is often associated with sandy deserts, cuts a broad path through its environment that may increase its vulnerability to obstacles. Our prior work demonstrated that sidewinding can be represented as a pair of orthogonal body waves (vertical and horizontal) that can be independently modulated to achieve high maneuverability and incline ascent, suggesting that sidewinders may also use template modulations to negotiate obstacles. To test this hypothesis, we recorded overhead video of four sidewinder rattlesnakes (Crotalus cerastes) crossing a line of vertical pegs placed in the substrate. Snakes used three methods to traverse the obstacles: a Propagate Through behavior in which the lifted moving portion of the snake was deformed around the peg and dragged through as the snake continued sidewinding (115/160 runs), Reversal turns that reorient the snake entirely (35/160), or switching to Concertina locomotion (10/160). The Propagate Through response was only used if the anterior-most region of static contact would propagate along a path anterior to the peg, or if a new region of static contact could be formed near the head to satisfy this condition; otherwise, snakes could only use Reversal turns or switch to Concertina locomotion. Reversal turns allowed the snake to re-orient and either escape without further peg contact or re-orient into a posture amenable to using the Propagate Through response. We developed an algorithm to reproduce the Propagate Through behavior in a robophysical model using a modulation of the two-wave template. This range of behavioral strategies provides sidewinders with a versatile range of options for effectively negotiating obstacles in their natural habitat, as well as provide insights into the design and control of robotic systems dealing with heterogeneous habitats.

Mitigating memory effects during undulatory locomotion on hysteretic materials.

While terrestrial locomotors often contend with permanently deformable substrates like sand, soil, and mud, principles of motion on such materials are lacking. We study the desert-specialist shovel-nosed snake traversing a model sand and find body inertia is negligible despite rapid transit and speed dependent granular reaction forces. New surface resistive force theory (RFT) calculation reveals how wave shape in these snakes minimizes material memory effects and optimizes escape performance given physiological power limitations. RFT explains the morphology and waveform-dependent performance of a diversity of non-sand-specialist snakes but overestimates the capability of those snakes which suffer high lateral slipping of the body. Robophysical experiments recapitulate aspects of these failure-prone snakes and elucidate how re-encountering previously deformed material hinders performance. This study reveals how memory effects stymied the locomotion of a diversity of snakes in our previous studies (Marvi et al., 2014) and indicates avenues to improve all-terrain robots.

Dynamics of scattering in undulatory active collisions.

Natural and artificial self-propelled systems must manage environmental interactions during movement. In complex environments, these interactions include active collisions, in which propulsive forces create persistent contacts with heterogeneities. Due to the driven and dissipative nature of these systems, such collisions are fundamentally different from those typically studied in classical physics. Here we experimentally and numerically study the effects of active collisions on a laterally undulating sensory-deprived robophysical model, whose dynamics are relevant to self-propelled systems across length scales and environments. Interactions with a single rigid post scatter the robot, and this deflection is dominated by head-post contact. These results motivate a model which reduces the snake to a circular particle with two key features: The collision dynamics are set by internal driving subject to the geometric constraints of the post, and the particle has an effective length equal to the wavelength of the snake. Interactions with a single row of evenly spaced posts (with interpost spacing d) produce distributions reminiscent of far-field diffraction patterns: As d decreases, distinct secondary peaks emerge as large deflections become more likely. Surprisingly, we find that the presence of multiple posts does not change the nature of individual collisions; instead, multimodal scattering patterns arise from multiple posts altering the likelihood of individual collisions to occur. As d decreases, collisions near the leading edges of the posts become more probable, and we find that these interactions are associated with larger deflections. Our results, which highlight the surprising dynamics that can occur during active collisions of self-propelled systems, can inform control principles for locomotors in complex terrain and facilitate design of task-capable active matter.

Tail use improves performance on soft substrates in models of early vertebrate land locomotors.

In the evolutionary transition from an aquatic to a terrestrial environment, early tetrapods faced the challenges of terrestrial locomotion on flowable substrates, such as sand and mud of variable stiffness and incline. The morphology and range of motion of appendages can be revealed in fossils; however, biological and robophysical studies of modern taxa have shown that movement on such substrates can be sensitive to small changes in appendage use. Using a biological model (the mudskipper), a physical robot model, granular drag measurements, and theoretical tools from geometric mechanics, we demonstrate how tail use can improve robustness to variable limb use and substrate conditions. We hypothesize that properly coordinated tail movements could have provided a substantial benefit for the earliest vertebrates to move on land.

Deformation-driven diffusion and plastic flow in amorphous granular pillars.

We report a combined experimental and simulation study of deformation-induced diffusion in compacted quasi-two-dimensional amorphous granular pillars, in which thermal fluctuations play a negligible role. The pillars, consisting of bidisperse cylindrical acetal plastic particles standing upright on a substrate, are deformed uniaxially and quasistatically by a rigid bar moving at a constant speed. The plastic flow and particle rearrangements in the pillars are characterized by computing the best-fit affine transformation strain and nonaffine displacement associated with each particle between two stages of deformation. The nonaffine displacement exhibits exponential crossover from ballistic to diffusive behavior with respect to the cumulative deviatoric strain, indicating that in athermal granular packings, the cumulative deviatoric strain plays the role of time in thermal systems and drives effective particle diffusion. We further study the size-dependent deformation of the granular pillars by simulation, and find that different-sized pillars follow self-similar shape evolution during deformation. In addition, the yield stress of the pillars increases linearly with pillar size. Formation of transient shear lines in the pillars during deformation becomes more evident as pillar size increases. The width of these elementary shear bands is about twice the diameter of a particle, and does not vary with pillar size.