Table of Contents
Fetching ...

Dispersion Relations for Active Undulators in Overdamped Environments

Christopher J. Pierce, Daniel Irvine, Lucinda Peng, Xuefei Lu, Hang Lu, Daniel I. Goldman

TL;DR

This work identifies a unifying dispersion relation between undulation frequency $\omega$ and wavenumber $k$ for overdamped locomotors, observed across nematodes and diverse environments. An active viscoelastic beam model, with internal and external dissipation and a phase-shifted neuromuscular input, yields $\operatorname{Re}\{\omega(k)\}$ and $\operatorname{Im}\{\omega(k)\}$ and reveals two scaling regimes: $\omega \propto k^{-2}$ in internal-dissipation-dominated conditions and $\omega \propto k^{2}$ when external dissipation dominates; stability imposes $k<1$. The model introduces regime constants $\alpha=(2\pi t_c s_c^2)^{-1}$ and $\beta=2\pi s_\eta^2/t_c$, linking frequency scales to wavelength $\lambda$ via $\omega \approx \alpha \lambda^2$ or $\omega \approx \beta/\lambda^2$, respectively. Experimental data on $\omega(k)$ from C. elegans, together with cross-species evidence from spermatozoa and fish larvae, support this unified framework, suggesting that environmental rheology and gait control co-determine how organisms adjust locomotor parameters in overdamped regimes.

Abstract

Organisms that locomote by propagating waves of body bending can maintain performance across heterogeneous environments by modifying their gait frequency $ω$ or wavenumber $k$. We identify a unifying relationship between these parameters for overdamped undulatory swimmers (including nematodes, spermatozoa, and mm-scale fish) moving in diverse environmental rheologies, in the form of an active `dispersion relation' $ω\propto k^{\pm2}$. A model treating the organisms as actively driven viscoelastic beams reproduces the experimentally observed scaling. The relative strength of rate-dependent dissipation in the body and the environment determines whether $k^2$ or $k^{-2}$ scaling is observed. The existence of these scaling regimes reflects the $k$ and $ω$ dependence of the various underlying force terms and how their relative importance changes with the external environment and the neuronally commanded gait.

Dispersion Relations for Active Undulators in Overdamped Environments

TL;DR

This work identifies a unifying dispersion relation between undulation frequency and wavenumber for overdamped locomotors, observed across nematodes and diverse environments. An active viscoelastic beam model, with internal and external dissipation and a phase-shifted neuromuscular input, yields and and reveals two scaling regimes: in internal-dissipation-dominated conditions and when external dissipation dominates; stability imposes . The model introduces regime constants and , linking frequency scales to wavelength via or , respectively. Experimental data on from C. elegans, together with cross-species evidence from spermatozoa and fish larvae, support this unified framework, suggesting that environmental rheology and gait control co-determine how organisms adjust locomotor parameters in overdamped regimes.

Abstract

Organisms that locomote by propagating waves of body bending can maintain performance across heterogeneous environments by modifying their gait frequency or wavenumber . We identify a unifying relationship between these parameters for overdamped undulatory swimmers (including nematodes, spermatozoa, and mm-scale fish) moving in diverse environmental rheologies, in the form of an active `dispersion relation' . A model treating the organisms as actively driven viscoelastic beams reproduces the experimentally observed scaling. The relative strength of rate-dependent dissipation in the body and the environment determines whether or scaling is observed. The existence of these scaling regimes reflects the and dependence of the various underlying force terms and how their relative importance changes with the external environment and the neuronally commanded gait.
Paper Structure (8 sections, 12 equations, 4 figures)

This paper contains 8 sections, 12 equations, 4 figures.

Figures (4)

  • Figure 1: Dispersion relation showing the approximate inverse quadratic scaling of $\omega \propto f$ with $k\propto 1/\lambda$ for C. elegans under changing environmental rheology. (a) Brightfield images of worm body shapes traveling through fluids of different viscosities, final positions are after 4.5s elapsed time. Violet curves illustrate full wave periods to show the changing wavelength $\lambda$. Inset, heatmap showing curvature $\kappa$ (in units of $mm^{-1}$) along the body coordinate $s$ (the arc length from head to tail in units of body length) and through time showing alternating bands of positive and negative curvature oscillating at a frequency $f$. (b) Postures of nematodes in buffer (b, i) and agar (b, ii) (scale bar, 1 mm). Linear (b, iii) and log plot (b, iii, inset) of the experimental dispersion relation for nematodes in diverse environments shown with the inverse quadratic (solid black line) and power law (dashed black line) fits, along with experimental data $\mathbf{\boldsymbol{*}}$ methylcellulose (0-3%), $\square$ Agar, $\boldsymbol{\times}$ PEG (1-5%), $\bullet$ dextran from Butler et al butler2015consistent, $\bigcirc$ buffer and $\boldsymbol{\square}$ agar from Fang-yen et al fang2010biomechanical. Power law exponent is $-2.28\pm0.30$.
  • Figure 2: An active damped beam model and resulting dispersion relations. (a) Model variables and internal and external forces for a generic undulator in a fluid. (b) The real part of the full dispersion relation (dashed, black curve) and the real part of the dispersion curve with the constraint of wave stability imposed (solid, gray curve). (b, inset) The same curves on a log-log scale. Dashed vertical lines indicate the value of $k$ separating quadratic and inverse scaling (purple) and the wave stability cutoff (green), for $s_c/s_\eta = 0.25$. Yellow shaded region shows the domain in $k$ which applies to the locomotion of nematodes.
  • Figure 3: Dispersion relation (\ref{['re']}) in natural units for different values of fluid drag coefficient $C$ taken from methylcellulose data. All the C. elegans data lie in the $k>s_c/s_{\eta}$ regime and are insensitive to the precise location of the peak at $k = s_c/s_{\eta}$.
  • Figure 4: Fish larvae and spermatozoa display dispersion relations where $f$ is an increasing function of $1/\lambda$, indicating that external fluid dissipation dominates internal dissipation. (a) Data from fish swimming meta-analysis with selected low-inertia swimmers separated by species (high-inertia swimmers plotted in black). Images of a mosquitofish G. amnis ($\sim$cm, Wikipedia) (b, i), zebrafish larva, D. rerio ($\sim$mm from Maes2012) (b, ii), an axolotl larva A. mexicanum ($\sim$cm photo credit John P. Clare) (b,iii), and a herring larva C. harengus ($\sim$mm from Fischbach2023) (b, vi). (c) Spermatozoa data from different species taken from Velho_Rodrigues2021 along with quadratic fit, $\mathbf{\boldsymbol{+}}$P. maxima (turbot), $\mathbf{\boldsymbol{\circ}}$A. curtula (beetle), $\mathbf{\boldsymbol{*}}$Lygaeus (milkweed bug), $\mathbf{\boldsymbol{\times}}$T. thynnus (tuna), $\mathbf{\boldsymbol{\square}}$M. merluccius (hake), $\mathbf{\boldsymbol{\diamond}}$C. capitata (fly), $\mathbf{\boldsymbol{\triangledown}}$ Cricket, $\mathbf{\boldsymbol{\triangle}}$B. marinus (toad), $\mathbf{\boldsymbol{\triangleright}}$Colobocentrotus (sea urchin), $\mathbf{\boldsymbol{\star}}$G. morhua (cod), $\mathbf{\boldsymbol{\triangleleft}}$Chaetopterus (annelid), $\mathbf{\boldsymbol{\davidsstar}}$Psammechinus (sea urchin), $\mathbf{\boldsymbol{+}}$G. domesticus (domestic fowl), $\mathbf{\boldsymbol{\circ}}$Ostrea (oyster), $\mathbf{\boldsymbol{*}}$Lytechinus (sea urchin), $\mathbf{\boldsymbol{\times}}$Ciona (tunicate), $\mathbf{\boldsymbol{\square}}$S. purpuratus (sea urchin), $\mathbf{\boldsymbol{\diamond}}$Myzostomus (worm), $\mathbf{\boldsymbol{\triangledown}}$C. mellus (midge), $\mathbf{\boldsymbol{\triangle}}$Homo (human), $\mathbf{\boldsymbol{\triangleright}}$Bos, bull, $\mathbf{\boldsymbol{\star}}$Ovis (ram), $\mathbf{\boldsymbol{\triangleleft}}$M. sclaris (fly), $\mathbf{\boldsymbol{\davidsstar}}$ Rabbit, $\mathbf{\boldsymbol{+}}$Mus (mouse), (d) Cricket sperm dispersion relation for multiple individuals with multiple measurements along each flagellum from Rikmenspoel1978 (different gray values represent different individuals), along with (e) an image of multiple waves on a single cricket flagellum ($\sim0.1-1$mm)Rikmenspoel1978.