<?xml version="1.0" encoding="utf-8" standalone="yes"?>
<rss version="2.0" xmlns:atom="http://www.w3.org/2005/Atom" xmlns:content="http://purl.org/rss/1.0/modules/content/">
  <channel>
    <title>Cats on Sebastian Spicker</title>
    <link>https://sebastianspicker.github.io/tags/cats/</link>
    <description>Recent content in Cats on Sebastian Spicker</description>
    <image>
      <title>Sebastian Spicker</title>
      <url>https://sebastianspicker.github.io/og-image.png</url>
      <link>https://sebastianspicker.github.io/og-image.png</link>
    </image>
    <generator>Hugo -- 0.160.0</generator>
    <language>en</language>
    <lastBuildDate>Mon, 23 Jun 2025 00:00:00 +0000</lastBuildDate>
    <atom:link href="https://sebastianspicker.github.io/tags/cats/index.xml" rel="self" type="application/rss+xml" />
    <item>
      <title>The Cat&#39;s Eye: Slit Pupils, Thin-Film Mirrors, and 135-Fold Dynamic Range</title>
      <link>https://sebastianspicker.github.io/posts/cat-eyes-slit-pupils-tapetum/</link>
      <pubDate>Mon, 23 Jun 2025 00:00:00 +0000</pubDate>
      <guid>https://sebastianspicker.github.io/posts/cat-eyes-slit-pupils-tapetum/</guid>
      <description>A cat&amp;rsquo;s eye contains two distinct optical technologies that human engineers have copied — one consciously, one not. The slit pupil achieves a dynamic range of 135:1 in light transmission, nearly ten times that of the human circular pupil. The tapetum lucidum is a multilayer thin-film reflector of crystalline rodlets, producing constructive interference at the peak of scotopic sensitivity and sending light through the retina twice. Banks et al. (Science Advances, 2015) showed why the slit geometry specifically evolved in ambush predators; Percy Shaw&amp;rsquo;s 1934 Catseye road reflector borrowed the principle directly.</description>
      <content:encoded><![CDATA[<p><em>Flash photography of cats produces glowing eyes. This is familiar enough that
most people do not find it strange. But the physics that produces it — a
biological multilayer interference reflector built from crystalline rodlets of
riboflavin and zinc, tuned to the peak of night-vision sensitivity, sending returning photons through
the retina for a second pass — is not familiar at all. I started thinking about
this after photographing our cats at dusk — through the doorway; they are indoor
cats now, for health reasons — and finding their eyes lit up a colour
that depends on the angle: greenish from straight ahead, golden from the side.
The angle-dependence is a direct consequence of the thin-film interference
condition, and the different colours correspond to different constructive
interference wavelengths at different angles of incidence.</em></p>
<p><em>The eye contains two optical solutions — pupil geometry and tapetum — that
address different aspects of the same problem: how to function across a very
large range of light levels, from bright midday sun to the dim luminance of a
starlit field.</em></p>
<hr>
<h2 id="the-dynamic-range-problem">The Dynamic Range Problem</h2>
<p>A crepuscular predator — active around dawn and dusk — must function visually
across a light-level range of roughly $10^8$:$1$. The sun on a bright day
produces retinal illuminance of around $10^5\,\mathrm{photons}/(\mu\mathrm{m}^2\cdot\mathrm{s})$;
a moonless night produces roughly $10^{-3}$ in the same units. The ratio is
approximately $10^8$.</p>
<p>The pupil is the variable aperture that controls how much light reaches the
retina. The larger the pupil area, the more light admitted; the smaller the
area, the less. For the human eye, the pupil diameter ranges from approximately
$2\,\mathrm{mm}$ (bright light) to $8\,\mathrm{mm}$ (darkness), giving a
maximum area ratio of:</p>
$$\frac{A_\mathrm{max}}{A_\mathrm{min}} = \left(\frac{8}{2}\right)^2 = 16.$$<p>This is a dynamic range of 16:1 from the pupil alone. The remaining
$10^8 / 16 \approx 6 \times 10^6$ factor in adaptation comes from neural
and photochemical mechanisms in the retina itself (photopigment bleaching,
dark adaptation of rods vs. cones, lateral inhibition).</p>
<p>For a domestic cat, the same measurement gives something different.</p>
<hr>
<h2 id="the-slit-pupil-1351-dynamic-range">The Slit Pupil: 135:1 Dynamic Range</h2>
<p>Banks, Sprague, Schmoll, Parnell, and Love published &ldquo;Why do animal eyes have
pupils of different shapes?&rdquo; in <em>Science Advances</em> in 2015 (1:7, e1500391).
They analysed pupil shape and size data from 214 terrestrial species and
correlated pupil geometry with ecological niche.</p>
<p>Their principal finding for slit pupils: the domestic cat pupil, a vertical
slit, achieves an area ratio of approximately <strong>135:1</strong> between maximum dilation
and maximum constriction. Numerically:</p>
$$\frac{A_\mathrm{max}}{A_\mathrm{min}} \approx 135.$$<p>The mechanism that makes this possible is geometrical. A circular pupil&rsquo;s
minimum area is limited by diffraction: constricting a circular aperture below
about $2\,\mathrm{mm}$ diameter produces diffraction rings that degrade image
quality. A slit, by contrast, can be made arbitrarily narrow in one direction
while retaining a larger dimension in the other, limiting diffraction in only
one axis. The vertical slit in a cat pupil can constrict to a width of
$\sim 0.3\,\mathrm{mm}$ while retaining a height of $\sim 15\,\mathrm{mm}$,
giving an area of roughly $0.3 \times 15 / (3.14 \times (8/2)^2) \times
A_\mathrm{max}$ — approximately 135 times smaller.</p>
<p>The 135:1 ratio is nearly <strong>nine times</strong> the dynamic range achievable by the
human circular pupil (16:1). This allows the cat&rsquo;s pupil to do substantially
more of the work of light adaptation, reducing the load on the slower neural
and photochemical mechanisms.</p>
<hr>
<h2 id="why-vertical-the-ecological-correlation">Why Vertical? The Ecological Correlation</h2>
<p>Banks et al. found a striking correlation between pupil geometry and predator
ecology:</p>
<ul>
<li><strong>Vertical slit pupils</strong> correlate with <em>ambush predators whose eyes are
close to the ground</em> — animals with shoulder height below approximately
$42\,\mathrm{cm}$.</li>
<li><strong>Horizontal slit pupils</strong> correlate with <em>prey animals and grazing
herbivores</em> (horses, goats, sheep, deer). The horizontal slit, when the
animal lowers its head to graze, rotates to remain approximately horizontal
(the eye counterrotates in the orbit), providing a wide panoramic field
of view for detecting approaching predators.</li>
<li><strong>Circular pupils</strong> correlate with <em>pursuit predators</em> (humans, dogs, large
raptors) that hunt at larger distances where the precise vertical depth
cues provided by the slit geometry are less critical.</li>
</ul>
<p>The functional advantage of a <strong>vertical slit for a low-to-the-ground ambush
predator</strong> is depth estimation by <em>blur circles</em>. The slit geometry produces
strong defocus blur in the horizontal direction but sharp focus in the vertical
direction. An ambush predator lying in grass needs to estimate the horizontal
distance to prey accurately; the defocus differential between horizontal and
vertical blur provides a stereoscopic-like depth cue even with one eye. This
is a form of <strong>astigmatic blur ranging</strong>: the degree of horizontal blur for a
given focal setting encodes the object&rsquo;s distance.</p>
<p>The correlation across 214 species is not perfect, but it is statistically
robust: slit pupils in ground-level ambush predators is not coincidence, it is
selection pressure.</p>
<hr>
<h2 id="the-tapetum-lucidum-a-biological-dielectric-mirror">The Tapetum Lucidum: A Biological Dielectric Mirror</h2>
<p>Behind the retina, most nocturnal and crepuscular mammals possess a reflective
layer called the <em>tapetum lucidum</em> (literally: &ldquo;bright carpet&rdquo;). Light that
passes through the retina without being absorbed by a photoreceptor strikes
the tapetum and is reflected back through the retina for a second absorption
opportunity. This roughly doubles the effective optical path length through
the photoreceptor layer, substantially increasing the probability of photon
capture at low light levels.</p>
<p>The cat tapetum is a <strong>tapetum cellulosum</strong>: a layer of specialised cells
whose cytoplasm contains dense arrays of rod-shaped crystalline inclusions
composed primarily of riboflavin (vitamin B$_2$) and zinc. (This is distinct
from the guanine-crystal tapeta found in fish and some reptiles.) The
crystalline rodlets have a refractive index of approximately $n_1 \approx 1.8$;
they alternate with layers of cytoplasm with refractive index $n_2 \approx
1.33$ (close to water). The rodlet arrays form a multilayer thin-film
reflector.</p>
<hr>
<h2 id="thin-film-interference-the-physics-of-the-reflection">Thin-Film Interference: The Physics of the Reflection</h2>
<p>The physics of the tapetum is identical to the physics of anti-reflection
coatings on camera lenses and dielectric mirrors in laser cavities.</p>
<p>Consider a single thin film of thickness $d$ and refractive index $n_1$
embedded between media of index $n_2 < n_1$. Light of wavelength $\lambda$
(in vacuum) incident at angle $\theta$ to the normal undergoes partial
reflection at both interfaces. The two reflected beams interfere
constructively when their optical path difference is a multiple of the
wavelength:</p>
$$\Delta = 2 n_1 d \cos\theta = m\lambda, \quad m = 1, 2, 3, \ldots$$<p>For the tapetum, typical rodlet diameter is $d \approx 100$–$120\,\mathrm{nm}$.
With $n_1 \approx 1.8$ and $\theta \approx 0°$ (normal incidence), the first
constructive interference maximum for a single layer occurs at:</p>
$$\lambda_\mathrm{peak} = 2 n_1 d = 2 \times 1.8 \times 100\,\mathrm{nm}
\approx 360\,\mathrm{nm}.$$<p>Wait — that is in the ultraviolet. The tapetum must have multiple layers.</p>
<p>For a stack of $N$ rodlet layers, the reflectance is strongly enhanced
(approaching unity for large $N$) and the peak wavelength of the fundamental
reflection maximum shifts. The relevant periodicity is the combined optical
thickness of one rodlet layer plus one cytoplasm layer:</p>
$$d_\mathrm{eff} = n_1 d_1 + n_2 d_2,$$<p>where $d_1 \approx 100\,\mathrm{nm}$ is the rodlet diameter and
$d_2 \approx 50$–$100\,\mathrm{nm}$ is the cytoplasm spacing. Taking
$d_2 \approx 60\,\mathrm{nm}$:</p>
$$d_\mathrm{eff} = 1.8 \times 100 + 1.33 \times 60 \approx 180 + 80
= 260\,\mathrm{nm}.$$<p>Constructive interference (quarter-wave condition for a multilayer stack) at
$m = 1$:</p>
$$\lambda_\mathrm{peak} = 2 d_\mathrm{eff} \approx 520\,\mathrm{nm}.$$<p>This is green — close to the peak of the scotopic (rod) sensitivity
curve at $\lambda_\mathrm{max,rod} = 498\,\mathrm{nm}$. The tapetum is tuned
to reflect the wavelengths that the night-vision photoreceptors are most
sensitive to. (The exact peak depends on rodlet spacing, which varies across
the tapetum; this produces the observed variation from green to yellow.)</p>
<p>The angle-dependence of the peak wavelength follows from the interference
condition: at angle $\theta$ to the normal, $\lambda_\mathrm{peak}(\theta)
= 2 d_\mathrm{eff} \cos\theta$. At $\theta = 30°$, $\cos 30° \approx 0.87$,
giving $\lambda_\mathrm{peak} \approx 450\,\mathrm{nm}$ — blue. At
$\theta = 60°$, $\cos 60° = 0.5$, giving $\lambda \approx 260\,\mathrm{nm}$ —
ultraviolet, invisible. The colour of eyeshine in a flash photograph therefore
depends on the angle between the camera and the eye, exactly as observed.</p>
<hr>
<h2 id="reflectance-of-a-multilayer-stack">Reflectance of a Multilayer Stack</h2>
<p>For $N$ identical bilayers (each of optical thickness $n_1 d_1 + n_2 d_2$),
the reflectance at the design wavelength is given by the transfer matrix
method. For the cat tapetum with $N \approx 10$–$15$ bilayers:</p>
$$R = \left(\frac{1 - (n_2/n_1)^{2N}}{1 + (n_2/n_1)^{2N}}\right)^2
\approx 1 - 4\left(\frac{n_2}{n_1}\right)^{2N}.$$<p>With $n_2/n_1 = 1.33/1.8 \approx 0.739$ and $N = 15$:</p>
$$(0.739)^{30} \approx 1.1 \times 10^{-4}.$$<p>The reflectance is approximately $1 - 4 \times 1.1 \times 10^{-4} \approx
0.9996$ — essentially $100\%$ at the design wavelength for a sufficiently thick
stack. The tapetum is a near-perfect reflector in a narrow wavelength band,
a biological dielectric mirror.</p>
<hr>
<h2 id="photon-statistics-at-low-light">Photon Statistics at Low Light</h2>
<p>The tapetum&rsquo;s function becomes clearest when framed in terms of photon
statistics. A single rod photoreceptor has an absorption probability of
approximately $\eta_\mathrm{single} \approx 25\%$ for a photon passing through
it once at $\lambda = 500\,\mathrm{nm}$.</p>
<p>With the tapetum reflecting the photon back for a second pass, the total
absorption probability becomes:</p>
$$\eta_\mathrm{total} = \eta + (1 - \eta)\, R\, \eta,$$<p>where $R \approx 1$ is the tapetum reflectance. For $\eta = 0.25$ and $R =
0.98$:</p>
$$\eta_\mathrm{total} = 0.25 + (0.75)(0.98)(0.25) = 0.25 + 0.184 \approx 0.43.$$<p>The double pass increases the photon detection efficiency from $25\%$ to
approximately $43\%$ — a factor of $1.7\times$.</p>
<p>At extremely low light levels, photon detection becomes a counting problem
governed by Poisson statistics. If a mean of $\bar{n}$ photons reaches a
single photoreceptor per integration time, the probability of detecting at
least one photon (and hence registering the presence of light) is:</p>
$$P(\text{detection}) = 1 - e^{-\bar{n}\,\eta_\mathrm{total}}.$$<p>For very dim stimuli where $\bar{n} \approx 1$–$3$ photons per rod per
integration time (close to the absolute threshold of cat vision at around
$7 \times 10^{-7}\,\mathrm{lux}$), increasing $\eta$ by a factor of $\sim
1.7$ has a significant effect on detection probability. The tapetum is not a
luxury
at low light levels; it is a biophysical necessity for sub-threshold light
detection.</p>
<hr>
<h2 id="percy-shaw-and-the-road-catseye">Percy Shaw and the Road Catseye</h2>
<p>In 1934, Percy Shaw, a road-mender from Halifax, applied for a British patent
for a retroreflective road stud that he called the &ldquo;Catseye.&rdquo; Shaw&rsquo;s stated
inspiration was the reflection of his car headlights from a cat&rsquo;s eyes while
driving on an unlit road at night. Whether this story is entirely accurate is
unclear, but the name and the inspiration are both documented in period sources.</p>
<p>Shaw&rsquo;s device uses a different retroreflection mechanism from the tapetum. The
tapetum produces specular (mirror-like) reflection in the back-focal plane of
the eye&rsquo;s lens — light returning along its incident path because the lens
refocuses it. Shaw&rsquo;s Catseye uses glass hemisphere retroreflectors (or, in
later versions, corner-cube retroreflectors) that return light toward its
source by total internal reflection rather than thin-film interference.</p>
<p>The corner-cube geometry guarantees retroreflection: any ray entering a trihedral
corner (three mutually perpendicular surfaces) reflects from all three surfaces
and exits parallel to the incident direction, regardless of the angle of
incidence. The mathematical proof is that the product of three reflections in
mutually perpendicular planes is the identity transformation on vectors up to
a sign change — the direction vector $\hat{v}$ exits as $-\hat{v}$, which is
exactly retroreflection.</p>
$$\hat{v}_\mathrm{out} = -\hat{v}_\mathrm{in}.$$<p>Shaw&rsquo;s road Catseye became standard equipment on British roads during the Second World War,
credited with a significant reduction in road fatalities during blackouts and
foggy conditions. The biological original was a multilayer interference mirror;
the engineering copy is a corner-cube retroreflector. Different physics, same
function, same name.</p>
<hr>
<h2 id="two-optical-solutions-to-one-problem">Two Optical Solutions to One Problem</h2>
<p>The cat&rsquo;s eye contains two distinct optical technologies:</p>
<ol>
<li>
<p><strong>The slit pupil</strong> — a variable aperture with 135:1 dynamic range, optimised
for depth estimation by astigmatic blur in a low-to-the-ground ambush predator.</p>
</li>
<li>
<p><strong>The tapetum lucidum</strong> — a multilayer thin-film reflector of riboflavin
crystalline rodlets, tuned to the scotopic sensitivity peak, achieving
near-100% reflectance at design wavelength and increasing photon detection
efficiency by a factor of approximately $1.7\times$.</p>
</li>
</ol>
<p>Both solutions were arrived at by natural selection over millions of years of
low-light hunting. Both have been copied — one consciously (Shaw&rsquo;s road
reflectors), one as a model for engineered multilayer reflectors in telescopes,
laser cavities, and narrowband optical filters.</p>
<p>When I photograph our cats at dusk and their eyes glow green, I am seeing
the thin-film interference of a biological photonic crystal — riboflavin
rodlets in cytoplasm — wavelength-selected to send green photons back through
rod cells for a second chance at absorption.
The green is not cosmetic. It is functional, and it is physics.</p>
<hr>
<h2 id="references">References</h2>
<ul>
<li>
<p>Banks, M.S., Sprague, W.W., Schmoll, J., Parnell, J.A.Q., &amp; Love, G.D.
(2015). Why do animal eyes have pupils of different shapes? <em>Science Advances</em>,
1(7), e1500391. <a href="https://doi.org/10.1126/sciadv.1500391">https://doi.org/10.1126/sciadv.1500391</a></p>
</li>
<li>
<p>Ollivier, F.J., Samuelson, D.A., Brooks, D.E., Lewis, P.A., Kallberg, M.E.,
&amp; Komaromy, A.M. (2004). Comparative morphology of the tapetum lucidum
(among selected species). <em>Veterinary Ophthalmology</em>, 7(1), 11–22.
<a href="https://doi.org/10.1111/j.1463-5224.2004.00318.x">https://doi.org/10.1111/j.1463-5224.2004.00318.x</a></p>
</li>
<li>
<p>Born, M., &amp; Wolf, E. (1999). <em>Principles of Optics</em> (7th ed.). Cambridge
University Press. (Chapters 1, 7 on thin-film interference and multilayer
coatings.)</p>
</li>
<li>
<p>Shaw, P. (1934). <em>Improvements in Studs for Roads and like Surfaces.</em> British
Patent 436,290. Applied 3 April 1934.</p>
</li>
<li>
<p>Warrant, E.J. (1999). Seeing better at night: Life style, eye design and the
optimum strategy of spatial and temporal summation. <em>Vision Research</em>, 39(9),
1611–1630. <a href="https://doi.org/10.1016/S0042-6989(98)00262-4">https://doi.org/10.1016/S0042-6989(98)00262-4</a></p>
</li>
</ul>
<hr>
<h2 id="changelog">Changelog</h2>
<ul>
<li><strong>2025-12-15</strong>: Corrected the adoption date of Percy Shaw&rsquo;s road Catseyes from &ldquo;from 1945 onward&rdquo; to &ldquo;during the Second World War&rdquo; (widespread adoption began under wartime blackout conditions, not after the war ended). Removed the Machan, Gu, &amp; Bharthuar (2020) reference, which could not be confirmed in available databases.</li>
</ul>
]]></content:encoded>
    </item>
    <item>
      <title>Why Cats Purr at 25 Hz: Vocal Fold Pads and the Physics of Self-Sustained Oscillation</title>
      <link>https://sebastianspicker.github.io/posts/purring-physics-vocal-fold-pads/</link>
      <pubDate>Mon, 09 Sep 2024 00:00:00 +0000</pubDate>
      <guid>https://sebastianspicker.github.io/posts/purring-physics-vocal-fold-pads/</guid>
      <description>For decades, the mechanism of purring was disputed. A 2023 paper in Current Biology showed that cat larynges purr without any neural input: airflow alone drives a self-sustained oscillation. The secret is connective tissue pads embedded in the vocal folds that increase effective mass and lower the resonant frequency to 25–30 Hz — the same range used clinically for bone- density stimulation and fracture healing under Wolff&amp;rsquo;s law.</description>
      <content:encoded><![CDATA[<p><em>The first thing either of our cats did when I sat still long enough was purr.
Not after food, not during play — the purr arrived when I sat down and held
still and they settled against me, and it arrived as a physical fact, a vibration
felt through the sternum and the ribs, not merely heard. The frequency was low:
around 25–30 cycles per second, which you can feel as a buzz rather than hear
as a tone. This is, I later confirmed, not far from the frequency at which
clinical devices stimulate bone growth. They are indoor cats now, on our vet&rsquo;s
recommendation — they find this unreasonable, but sitting still and being purred
on has become a regular feature of working from home.</em></p>
<p><em>The physics of how the larynx produces that frequency is, as of 2023, finally
resolved — and the mechanism is more elegant than anyone suspected.</em></p>
<hr>
<h2 id="the-frequency-and-its-peculiarity">The Frequency and Its Peculiarity</h2>
<p>Domestic cats purr at approximately $25$–$30\,\mathrm{Hz}$. This is
remarkably low for an animal of cat size. A human vocal fold — roughly
comparable in size — vibrates at $85$–$255\,\mathrm{Hz}$ for normal speech.
A cat&rsquo;s larynx is smaller than a human&rsquo;s, not larger, which makes the low
frequency surprising: in a simple spring-mass oscillator model, smaller and
lighter vocal folds should vibrate <em>faster</em>, not slower.</p>
<p>The frequency range $25$–$50\,\mathrm{Hz}$ has clinical significance in a
different field. Therapeutic vibration platforms used in sports medicine and
osteoporosis treatment operate in exactly this range, exploiting Wolff&rsquo;s law
(bone remodelling under mechanical stress) to increase bone density and
accelerate fracture repair. The coincidence is suggestive. It was first
noted quantitatively by von Muggenthaler (2001, <em>Journal of the Acoustical
Society of America</em> 110, 2666), who recorded purrs from 44 felids and
found that all produced dominant frequencies between $25$ and $150\,\mathrm{Hz}$.</p>
<p>Whether cats deliberately exploit this frequency for self-healing is a separate
biological question. The physics question is simpler: how does the larynx
produce it?</p>
<hr>
<h2 id="flow-induced-vocal-fold-oscillation">Flow-Induced Vocal Fold Oscillation</h2>
<p>Vocal fold oscillation in mammals is a flow-induced, self-sustained mechanical
phenomenon. The Bernoulli effect and elastic restoring forces create a
feedback loop that keeps the folds oscillating as long as subglottal air
pressure is maintained.</p>
<p>The mechanism is as follows. The lungs supply a steady subglottal pressure
$p_\mathrm{sub}$. This drives airflow through the glottis (the gap between the
vocal folds). As the folds are pushed apart by the pressure, the airflow
velocity in the narrowed glottis increases; by Bernoulli&rsquo;s principle,</p>
$$p + \tfrac{1}{2}\rho v^2 = \mathrm{const},$$<p>the pressure drops, drawing the folds back together. The folds&rsquo; elastic
restoring force adds to this: they spring back when displaced. The result is
an oscillation — the folds open and close periodically, chopping the airflow
into pressure pulses that we perceive as sound (or vibration, for low
frequencies).</p>
<p>The fundamental frequency is approximately:</p>
$$f_0 \approx \frac{1}{2L}\sqrt{\frac{T}{\rho_s}},$$<p>where $L$ is the vibrating length of the vocal fold, $T$ is the longitudinal
tension, and $\rho_s$ is the surface density (mass per unit area). This is
the same formula as for a vibrating string — and the physics is closely
related.</p>
<p>For a cat-sized larynx with $L \approx 1\,\mathrm{cm}$, realistic tissue
tension, and tissue density $\rho_s \sim 1\,\mathrm{kg/m}^2$, this formula
gives $f_0$ in the hundreds of hertz — far above the observed purring
frequency of $25$–$30\,\mathrm{Hz}$.</p>
<p>Something is missing from the model.</p>
<hr>
<h2 id="the-long-standing-controversy">The Long-Standing Controversy</h2>
<p>Until 2023, the dominant explanation for the low purring frequency was the
<strong>Active Muscular Contraction (AMC) hypothesis</strong>: the laryngeal muscles
contract rhythmically at the purring frequency, mechanically driving the
vocal folds rather than relying on passive aeroelastic oscillation. On this
view, purring is more like a drumming than a singing — the neural drive at
$25$–$30\,\mathrm{Hz}$ sets the frequency, overriding the natural aeroelastic
frequency.</p>
<p>The AMC hypothesis was difficult to test directly because the larynx is
inaccessible in a live, purring cat without interfering with the purr.
Electromyographic recordings from laryngeal muscles of purring cats showed
rhythmic activity consistent with the AMC hypothesis, but causality was unclear:
were the muscles driving the oscillation, or responding to it?</p>
<p>The alternative hypothesis — that purring is passive, driven purely by
aeroelastic forces — faced the problem noted above: the aeroelastic frequency
of a cat-sized larynx should be far too high to explain $25$–$30\,\mathrm{Hz}$.
Unless something was being added to the vocal folds to lower their effective
resonant frequency.</p>
<hr>
<h2 id="herbst-et-al-2023-the-mass-loading-mechanism">Herbst et al. 2023: The Mass-Loading Mechanism</h2>
<p>In October 2023, Christian Herbst and colleagues at the University of Vienna
published &ldquo;Domestic cat larynges can produce purring frequencies without neural
input&rdquo; (<em>Current Biology</em> 33, 4727–4732). The experiment was decisive.</p>
<p>The team excised larynges from domestic cats (post-mortem, within a short time
window to preserve tissue properties) and mounted them in a flow bench: a
controlled airflow was supplied to the subglottal side, and the larynges were
held at physiologically realistic tension and hydration.</p>
<p><strong>The result</strong>: all eight excised larynges produced self-sustained oscillations
at $25$–$30\,\mathrm{Hz}$ — the normal purring frequency — without any neural
input whatsoever. No muscular contraction was present (no motor neurons, no
calcium signalling, no ATP). The oscillation was purely passive, driven by the
airflow and maintained by the tissue mechanics.</p>
<p>This ruled out the AMC hypothesis. The neural drive is not needed to sustain
the oscillation; it may modulate it, start or stop it, but the fundamental
frequency is set by the tissue mechanics, not the neural firing rate.</p>
<p>The follow-up finding was the key to the physics: histological analysis of the
vocal fold tissue revealed <strong>connective tissue pads</strong> embedded in the vocal
fold mucosa, up to $4\,\mathrm{mm}$ thick. These pads are not present in the
vocal folds of humans or other mammals that do not purr. They increase the
effective mass of the oscillating tissue significantly, without adding
corresponding stiffness.</p>
<hr>
<h2 id="the-mass-loading-physics">The Mass-Loading Physics</h2>
<p>The fundamental frequency of a harmonic oscillator is:</p>
$$f_0 = \frac{1}{2\pi}\sqrt{\frac{k}{m}},$$<p>where $k$ is the effective stiffness and $m$ is the effective mass. Adding mass
(at constant stiffness) lowers the frequency as $f_0 \propto m^{-1/2}$.</p>
<p>For the vocal folds, the spring constant $k$ is set by tissue tension and
elasticity — properties that the tissue pads do not significantly alter. But
the pads add a substantial mass $\Delta m$ to the oscillating system. The
purring frequency becomes:</p>
$$f_\mathrm{purr} = \frac{1}{2\pi}\sqrt{\frac{k}{m_0 + \Delta m}},$$<p>where $m_0$ is the baseline vocal fold mass and $\Delta m$ is the added mass
from the pads.</p>
<p>As a rough estimate: if the unloaded aeroelastic frequency were in the
range $f_\mathrm{normal} \approx 200$–$400\,\mathrm{Hz}$ (the range of
cat meow fundamental frequencies), lowering it to $f_\mathrm{purr} \approx
25\,\mathrm{Hz}$ would require a mass increase by a factor of</p>
$$\frac{m_0 + \Delta m}{m_0} = \left(\frac{f_\mathrm{normal}}{f_\mathrm{purr}}\right)^2
\approx 64\text{–}256.$$<p>This is a large factor, but not implausible for pads up to 4 mm thick
embedded in a mucosal membrane that is itself very thin. The simple
harmonic oscillator model is an idealisation — the actual frequency reduction
also involves changes in vibration mode shape, tissue coupling, and
aerodynamic loading — but the mass-loading effect is the dominant mechanism.
The tissue pads are, in effect, frequency dividers: they convert a
high-frequency aeroelastic oscillator into a low-frequency vibration
generator.</p>
<p>This is the same principle used in engineering to lower the natural frequency
of mechanical structures: add mass without changing stiffness. Tuned mass
dampers in skyscrapers work on the same principle. So do the heavy flywheel
weights added to engines to suppress rotational vibration.</p>
<p>The cat&rsquo;s larynx evolved this solution independently, and with a mass ratio
that would impress a structural engineer.</p>
<hr>
<h2 id="the-self-sustained-oscillation-criterion">The Self-Sustained Oscillation Criterion</h2>
<p>Not every mass-loaded oscillator will self-sustain under airflow. The
Bernoulli-elastic feedback loop must overcome the viscous damping of the
tissue. A dimensional scaling estimate for the critical subglottal pressure is:</p>
$$p^* \sim \eta_\mathrm{tissue} \cdot \frac{v}{L} \sim \eta_\mathrm{tissue} \cdot f_0,$$<p>where $\eta_\mathrm{tissue}$ is the tissue viscosity, $v \sim f_0 L$ is the
characteristic mucosal wave velocity, and $L$ is the fold length. (The full
phonation threshold pressure, as derived by Titze (2006), depends on
additional geometric and aerodynamic parameters.) For typical laryngeal tissue properties and the observed purring
frequency, this critical pressure is of order $100$–$200\,\mathrm{Pa}$ —
low enough to be sustained by the respiratory system without extraordinary
effort.</p>
<p>This is consistent with the observation that cats can purr both during
inhalation and exhalation, maintaining a continuous acoustic output throughout
the breathing cycle. The oscillation threshold is low enough that normal
respiration can maintain it.</p>
<hr>
<h2 id="wolffs-law-and-the-25-hz-coincidence">Wolff&rsquo;s Law and the 25 Hz Coincidence</h2>
<p>Julius Wolff (1892) proposed that bone remodels in response to mechanical
loading: osteoblasts (bone-building cells) are stimulated by cyclic compressive
stress, while osteoclasts (bone-resorbing cells) dominate in the absence of
loading. This principle — now called Wolff&rsquo;s law — underpins the use of
therapeutic vibration in orthopaedics.</p>
<p>The optimal frequency for osteoblast stimulation, determined empirically in
clinical studies, is $20$–$50\,\mathrm{Hz}$. Vibration at these frequencies,
applied at amplitudes of $0.2$–$1.0\,g$ (where $g$ is gravitational
acceleration), produces measurable increases in bone mineral density, accelerates
fracture healing, and reduces bone loss in microgravity. The frequency range
is not a narrow resonance; it reflects the natural frequencies of cellular
mechanotransduction pathways involving focal adhesion kinase (FAK) and
integrin signalling.</p>
<p>Cat purring produces vibration in the frequency range $25$–$50\,\mathrm{Hz}$
at the body surface. Whether this is sufficient to produce meaningful bone
stimulation — and whether cats evolved purring partly as a bone-maintenance
mechanism — is not yet resolved by controlled experiments. The hypothesis is
physiologically plausible: cats conserve metabolic energy by resting for up
to 16 hours per day, and during this rest period, bone would normally be
unstressed and subject to resorption. A continuous low-frequency vibration
during rest could counteract this.</p>
<p>This is speculative at the level of evolutionary causation. What is not
speculative is that the purring frequency overlaps precisely with the
therapeutic vibration range, and that this overlap is not obviously accidental.</p>
<hr>
<h2 id="across-felid-species">Across Felid Species</h2>
<p>Von Muggenthaler&rsquo;s 2001 survey of 44 felids found that most domestic
cats purr in the range $25$–$30\,\mathrm{Hz}$, with harmonics at $50$,
$75\,\mathrm{Hz}$, and so on. Cheetahs purr at $20$–$25\,\mathrm{Hz}$;
pumas (mountain lions) at $20$–$30\,\mathrm{Hz}$; servals and ocelots at
$22$–$28\,\mathrm{Hz}$.</p>
<p>The large roaring cats — lions, tigers, leopards, jaguars — do not purr in
the continuous sense that domestic cats do. Their enlarged hyoid apparatus
allows roaring by a different mechanism (a modified laryngeal pad that
allows very low-frequency, high-intensity sound production). Some large cats
produce purr-like sounds during exhalation but not the continuous through-
inhalation-and-exhalation purring of smaller felids.</p>
<p>The vocal fold pad mechanism appears to be specific to the non-roaring felids,
though detailed histological comparisons across species are still sparse.</p>
<hr>
<h2 id="what-i-hear">What I Hear</h2>
<p>When one of our cats purrs while settled against me, what I am feeling is the
mechanical resonance of a mass-loaded aeroelastic oscillator at approximately
$25\,\mathrm{Hz}$, the frequency having been lowered by connective tissue pads
from a natural aeroelastic frequency several hundred hertz higher. The pads
evolved, we think, to produce exactly this frequency — sustained under normal
respiratory airflow pressure with no additional muscular energy. The acoustic
output is a byproduct of a vibration.</p>
<p>Whether the vibration serves a direct physiological function in the cat&rsquo;s own
bones is, as of this writing, still an open question. What seems clear is that
the 2023 paper settled the mechanism question conclusively: the frequency is
set by mass loading, not neural drive. The larynx purrs by itself when you
blow air through it.</p>
<p>I find this reassuring. The physics is in the cat, not in its nervous system.
The cat purrs the way a tuning fork rings — not because it decides to, but
because that is what it does when the conditions are right.</p>
<hr>
<h2 id="references">References</h2>
<ul>
<li>
<p>Herbst, C.T., Prigge, T., Garcia, M., Hampala, V., Hofer, R., Weissengruber,
G.E., Svec, J.G., &amp; Fitch, W.T. (2023). Domestic cat larynges can produce
purring frequencies without neural input. <em>Current Biology</em>, 33(22),
4727–4732.e4. <a href="https://doi.org/10.1016/j.cub.2023.09.014">https://doi.org/10.1016/j.cub.2023.09.014</a></p>
</li>
<li>
<p>von Muggenthaler, E. (2001). The felid purr: A healing mechanism?
<em>Journal of the Acoustical Society of America</em>, 110(5), 2666.
<a href="https://doi.org/10.1121/1.4777098">https://doi.org/10.1121/1.4777098</a></p>
</li>
<li>
<p>Titze, I.R. (2006). <em>The Myoelastic Aerodynamic Theory of Phonation.</em>
National Center for Voice and Speech.</p>
</li>
<li>
<p>Wolff, J. (1892). <em>Das Gesetz der Transformation der Knochen.</em> A. Hirschwald.
(English translation: Maquet, P., &amp; Furlong, R., 1986. <em>The Law of Bone
Remodelling.</em> Springer.)</p>
</li>
<li>
<p>Rubin, C.T., &amp; Lanyon, L.E. (1984). Regulation of bone formation by applied
dynamic loads. <em>Journal of Bone and Joint Surgery</em>, 66(3), 397–402.
<a href="https://doi.org/10.2106/00004623-198466030-00012">https://doi.org/10.2106/00004623-198466030-00012</a></p>
</li>
<li>
<p>Christiansen, P. (2008). Evolution of skull and mandible shape in cats
(Carnivora: Felidae). <em>PLOS ONE</em>, 3(7), e2807.
<a href="https://doi.org/10.1371/journal.pone.0002807">https://doi.org/10.1371/journal.pone.0002807</a></p>
</li>
</ul>
]]></content:encoded>
    </item>
    <item>
      <title>How Cats Drink: Inertia, Gravity, and the Froude Number at the Tip of a Tongue</title>
      <link>https://sebastianspicker.github.io/posts/how-cats-drink-froude-number/</link>
      <pubDate>Mon, 22 Jul 2024 00:00:00 +0000</pubDate>
      <guid>https://sebastianspicker.github.io/posts/how-cats-drink-froude-number/</guid>
      <description>Cats do not scoop water with their tongues — they exploit a delicate balance between inertia and gravity at the air-water interface. The tip of the tongue just touches the surface; rapid withdrawal pulls a fluid column upward; the jaw closes at exactly the moment the column peaks. Reis, Jung, Aristoff, and Stocker (Science, 2010) showed that the lapping frequency of all felids — from domestic cats to lions — is tuned so that the Froude number at the tongue tip is approximately unity.</description>
      <content:encoded><![CDATA[<p><em>I have spent a non-trivial amount of time watching our cats drink — they are
indoor-only cats, on our vet&rsquo;s advice, which gives them few distractions and
gives me ample opportunity to observe. This is not entirely voluntary. Once you have noticed that something is happening at the
water bowl that does not look right — the tongue moves too fast, the water
column is pulled upward rather than scooped, the jaw closes before the tongue
returns — you find yourself crouching beside the bowl with your phone propped
against a chair, filming at 240 frames per second and feeling that you have
perhaps chosen an unusual way to spend a Tuesday morning.</em></p>
<p><em>Pedro Reis, Sunghwan Jung, Jeffrey Aristoff, and Roman Stocker had the same
impulse, with better equipment. Their 2010 paper in Science, &ldquo;How Cats Lap:
Water Uptake by Felis catus,&rdquo; is one of the more elegant pieces of dimensional
analysis in recent biology.</em></p>
<hr>
<h2 id="how-cats-do-not-drink">How Cats Do Not Drink</h2>
<p>The simplest hypothesis — that cats curl the tongue into a spoon and scoop
water into the mouth — is false. High-speed photography shows that the cat&rsquo;s
tongue does not form a cup shape. Instead, the cat extends the tongue tip
downward toward the water surface and then rapidly retracts it. The motion is
fast — too fast for normal video — and the tongue barely contacts the surface.</p>
<p>The contrast with dogs is instructive. Dogs <em>do</em> scoop: the tongue curls
backward (not forward), forming a ladle shape that scoops water upward and
backwards into the mouth. The mechanism is vigorous and inefficient — a
significant fraction of the water misses the mouth entirely, which is why
drinking dogs produce splashing and dogs often have wet chins. The mechanism
works but is inelegant.</p>
<p>Cats produce almost no splash. The mechanism is different in kind.</p>
<hr>
<h2 id="the-physical-mechanism">The Physical Mechanism</h2>
<p>Reis et al. (2010) used high-speed photography (1000 frames per second) to
resolve the cat&rsquo;s lapping motion. Their observations:</p>
<ol>
<li>
<p>The cat extends the tongue tip downward until the <em>dorsal surface</em> (the top
side) just touches the water surface. The ventral surface (the smooth
underside) does not contact the water.</p>
</li>
<li>
<p>The cat then rapidly retracts the tongue upward. The tongue tip is moving
at roughly $v \approx 0.7\,\mathrm{m/s}$ during this retraction.</p>
</li>
<li>
<p>As the tongue tip pulls away from the surface, a column of liquid is pulled
upward by the adhesion between the liquid and the retreating tongue. The
column rises against gravity.</p>
</li>
<li>
<p>The column eventually stalls — inertia is overcome by gravity — and begins
to fall back. The cat closes its jaw at exactly the moment of maximum column
height, capturing the peak volume of water.</p>
</li>
<li>
<p>The cat then extends the tongue for the next lap.</p>
</li>
</ol>
<p>The cat closes its jaw before the tongue fully retracts. This is important:
the jaw closure captures the water column, not the water adhering to the tongue.
The tongue is the mechanism that <em>creates</em> the column; the jaw captures it.</p>
<hr>
<h2 id="dimensional-analysis-the-froude-number">Dimensional Analysis: The Froude Number</h2>
<p>The relevant competition is between <strong>inertia</strong> (which drives the column
upward) and <strong>gravity</strong> (which pulls it back down). Surface tension plays a
role in stabilising the column but is not the primary factor governing the
column height.</p>
<p>The balance between inertia and gravity for a fluid column moving at speed
$v$ and of characteristic length scale $L$ (here, the diameter of the tongue
tip, $L \approx 5\,\mathrm{mm}$ for a domestic cat) is captured by the
<strong>Froude number</strong>:</p>
$$\mathrm{Fr} = \frac{v}{\sqrt{gL}},$$<p>where $g = 9.81\,\mathrm{m/s}^2$ is gravitational acceleration.</p>
<p>When $\mathrm{Fr} \ll 1$: gravity dominates, inertia is insufficient to pull a
significant column of water upward. Very slow tongue motion would lift almost
no water.</p>
<p>When $\mathrm{Fr} \gg 1$: inertia dominates, the column rises far above the
surface but the jaw must be closed quickly before the large amount of water
falls back. Very fast tongue motion wastes water and requires rapid jaw closure.</p>
<p>The optimal lapping frequency — maximising captured volume per lap — occurs
near $\mathrm{Fr} \approx 1$, where inertial and gravitational forces are
comparable and the column height is matched to the jaw closure dynamics.</p>
<h3 id="checking-the-numbers-for-a-domestic-cat">Checking the Numbers for a Domestic Cat</h3>
<p>For a domestic cat:</p>
<ul>
<li>Tongue tip diameter: $L \approx 5\,\mathrm{mm} = 5 \times 10^{-3}\,\mathrm{m}$</li>
<li>Characteristic tongue tip speed: $v \approx 0.7\,\mathrm{m/s}$</li>
</ul>
$$\mathrm{Fr} = \frac{0.7}{\sqrt{9.81 \times 5 \times 10^{-3}}}
= \frac{0.7}{\sqrt{0.049}} = \frac{0.7}{0.22} \approx 3.2.$$<p>Reis et al. found Fr of order unity — inertial and gravitational forces
comparable — confirming that the lapping speed is tuned to the inertia-gravity
balance. (The exact numerical value depends on the choice of characteristic
length scale; using the tongue tip diameter as above gives Fr in the range
1–3, squarely in the regime where neither force dominates.)</p>
<hr>
<h2 id="scaling-across-felids">Scaling Across Felids</h2>
<p>The Froude number prediction yields a scaling law for lapping frequency across
felid species of different sizes. If all felids lap at $\mathrm{Fr} \approx 1$,
then the characteristic speed scales as $v \sim \sqrt{gL}$, and the lapping
frequency scales as:</p>
$$f = \frac{v}{d} \sim \frac{\sqrt{gL}}{d},$$<p>where $d$ is the distance the tongue travels per lap (roughly proportional to
tongue length, which scales with body size). Since $L \sim d$ scales with body
size, we get:</p>
$$f \sim \frac{\sqrt{g \cdot d}}{d} = \sqrt{\frac{g}{d}} \propto d^{-1/2}.$$<p>Larger cats have longer tongues and lap more slowly. The prediction is that
lapping frequency scales as the square root of inverse tongue length — or,
equivalently, as the inverse square root of body mass (since linear dimensions
scale as mass$^{1/3}$):</p>
$$f \propto m^{-1/6}.$$<p>Reis et al. tested this against high-speed footage of large felids. A domestic
cat laps at approximately $4\,\mathrm{Hz}$; a lion laps at approximately
$1.2\,\mathrm{Hz}$; a tiger at roughly $1\,\mathrm{Hz}$. The scaling is
consistent with $f \propto m^{-1/6}$ across three orders of magnitude in
body mass.</p>
<p>The table below shows the predicted versus observed scaling:</p>
<table>
  <thead>
      <tr>
          <th>Species</th>
          <th>Body mass (kg)</th>
          <th>Predicted $f$ relative to cat</th>
          <th>Predicted $f$ (Hz)</th>
          <th>Observed $f$ (Hz)</th>
      </tr>
  </thead>
  <tbody>
      <tr>
          <td>Domestic cat</td>
          <td>4</td>
          <td>1.0</td>
          <td>4.0</td>
          <td>~4.0</td>
      </tr>
      <tr>
          <td>Jaguar</td>
          <td>80</td>
          <td>$\left(\frac{4}{80}\right)^{1/6} \approx 0.61$</td>
          <td>2.4</td>
          <td>~2.0</td>
      </tr>
      <tr>
          <td>Lion</td>
          <td>200</td>
          <td>$\left(\frac{4}{200}\right)^{1/6} \approx 0.52$</td>
          <td>2.1</td>
          <td>~1.5</td>
      </tr>
      <tr>
          <td>Tiger</td>
          <td>220</td>
          <td>$\left(\frac{4}{220}\right)^{1/6} \approx 0.51$</td>
          <td>2.1</td>
          <td>~1.0</td>
      </tr>
  </tbody>
</table>
<p>The $m^{-1/6}$ scaling captures the correct trend — larger cats lap more
slowly — though the predicted frequencies for the largest cats somewhat
overestimate the observed values. The discrepancy may reflect the limitations
of the simple allometric assumption (that all linear dimensions scale as
$m^{1/3}$) and the fact that tongue geometry does not scale isometrically
across the full range of felid body sizes.</p>
<hr>
<h2 id="why-not-just-lick">Why Not Just Lick?</h2>
<p>A natural question: why not simply allow the tongue to fully submerge and
absorb water through the papillae, as the tongue already contacts water when
lapping? Several answers:</p>
<ol>
<li>
<p><strong>Papillae are not sponges.</strong> Feline papillae are hollow and scoop-shaped
(filiform papillae with hollow tips), optimised for grooming and food
manipulation, not passive absorption. Active wicking is limited.</p>
</li>
<li>
<p><strong>The cat cannot breathe with its mouth submerged.</strong> A lapping mechanism
that keeps the mouth mostly closed except for the brief jaw-closure moment
allows continuous breathing through the nose during drinking.</p>
</li>
<li>
<p><strong>Speed and efficiency.</strong> The inertial column mechanism delivers significantly
more water per jaw movement than surface tension adhesion alone. At 4 laps
per second, a domestic cat takes in roughly $0.14\,\mathrm{mL}$ per lap,
for a total of roughly $34\,\mathrm{mL/min}$ — comparable to sipping rates
in animals that use more direct intake mechanisms.</p>
</li>
</ol>
<p>The cat has converged on a hydrodynamically optimal strategy under the
constraint of keeping the oral cavity mostly sealed during the intake cycle.</p>
<hr>
<h2 id="the-robotic-tongue">The Robotic Tongue</h2>
<p>Reis et al. constructed a robotic cat tongue to verify the mechanism: a smooth
glass disc lowered to the water surface and retracted at controlled speeds.
The column height as a function of speed followed the predicted inertia-gravity
balance, confirming that the mechanism does not depend on any specifically
biological property of the tongue — it is a fluid dynamics result that applies
to any surface moving away from a water interface at the right speed.</p>
<p>The robot lapped at the same Froude number as the cat.</p>
<hr>
<h2 id="dogs-horses-and-the-comparison">Dogs, Horses, and the Comparison</h2>
<p>Dogs cup the tongue <em>caudally</em> (backwards) rather than ventrally, forming a
ladle. The mechanism is faster and delivers more water per stroke but is
messy — the ladle is formed outside the mouth, and water sloshes freely. Dogs
lap at roughly $3\,\mathrm{Hz}$ with a tongue tip speed significantly higher
than cats, producing Fr well above unity. The excess inertia is why dog
drinking generates splashing.</p>
<p>Horses, by contrast, create a near-seal with their lips and use suction —
a fundamentally different mechanism that requires no tongue projection at all.
The lapping mechanism of felids is phylogenetically specific and appears to
have evolved under selection pressure for both efficiency and noise suppression,
consistent with the ambush-predator lifestyle. A cat that splashed while
drinking would alert prey at a water source. A cat that laps near-silently
does not.</p>
<hr>
<h2 id="a-note-on-the-measurement">A Note on the Measurement</h2>
<p>Getting reliable high-speed footage of a cat drinking is harder than it sounds.
Our cats drink at different times of day, in different moods, and the presence
of a camera tripod next to the water bowl is regarded as grounds for drinking
elsewhere. Pedro Reis et al. solved this by filming their laboratory cat, Cutta
Cutta, in a controlled setting. Their footage is available online and is
genuinely beautiful: a slow-motion waterfall in miniature, rising improbably
from the tongue tip and held there by the balance between upward momentum and
downward gravity, until the jaw swings shut.</p>
<p>The physics is in the timing.</p>
<hr>
<h2 id="references">References</h2>
<ul>
<li>
<p>Reis, P.M., Jung, S., Aristoff, J.M., &amp; Stocker, R. (2010). How cats lap:
Water uptake by <em>Felis catus</em>. <em>Science</em>, 330(6008), 1231–1234.
<a href="https://doi.org/10.1126/science.1195421">https://doi.org/10.1126/science.1195421</a></p>
</li>
<li>
<p>Aristoff, J.M., Stocker, R., Jung, S., &amp; Reis, P.M. (2011). On the water
lapping of felines and the water running of lizards. <em>Communicative &amp;
Integrative Biology</em>, 4(2), 213–215.</p>
</li>
<li>
<p>Vogel, S. (1994). <em>Life in Moving Fluids: The Physical Biology of Flow</em>
(2nd ed.). Princeton University Press.</p>
</li>
</ul>
<hr>
<h2 id="changelog">Changelog</h2>
<ul>
<li><strong>2025-12-15</strong>: Updated water intake per lap from 0.04 mL to 0.14 mL (Reis et al. report ~0.14 +/- 0.04 mL per lap; the previous value was the standard deviation), and updated the intake rate accordingly (~34 mL/min). Updated the papillae location from ventral to dorsal surface. Updated the Aristoff et al. reference to the correct 2011 <em>Communicative &amp; Integrative Biology</em> article. Removed the Jung &amp; Kim (2012) PRL reference (article number 034501 resolves to a different paper).</li>
</ul>
]]></content:encoded>
    </item>
    <item>
      <title>Are Cats Liquid? The Deborah Number and the Rheology of Cats</title>
      <link>https://sebastianspicker.github.io/posts/liquid-cats-deborah-number/</link>
      <pubDate>Wed, 03 Apr 2024 00:00:00 +0000</pubDate>
      <guid>https://sebastianspicker.github.io/posts/liquid-cats-deborah-number/</guid>
      <description>Marc-Antoine Fardin won the 2017 Ig Nobel Prize in Physics for proving, rigorously, that cats are liquid. The argument rests on the Deborah number De = τ/T: if the material&amp;rsquo;s relaxation time τ is shorter than the observation time T, the material behaves as a fluid. A cat filling a sink (De ≈ 0.008) is a liquid. A cat bouncing off a table (De ≫ 1) is a solid. The classification is not a joke — it is standard rheology, applied to an unusual substrate.</description>
      <content:encoded><![CDATA[<p><em>One of our strays discovered, sometime in her first winter indoors — they are
strictly indoor cats now, on our vet&rsquo;s recommendation — that she could fit into
a salad bowl. Not sit beside it, not rest her head on its rim: fit into it,
curled into a precise sphere with her tail tucked under her chin and her ears
folded flat, filling the bowl as liquid fills a container. The bowl has a
diameter of 22 centimetres. I did not find this as surprising as perhaps I
should have: there is a quantity in materials science that determines, rigorously,
whether a given material in a given situation should be classified as a solid or
a liquid. For a cat in a bowl, this quantity is comfortably below one.</em></p>
<p><em>The material is a liquid. The material is also a cat.</em></p>
<hr>
<h2 id="the-definition-of-a-fluid">The Definition of a Fluid</h2>
<p>The intuitive distinction between solids and liquids is that solids hold their
shape and liquids conform to their container. But this distinction is one of
timescale, not of material identity.</p>
<p>A classic demonstration: place a ball of silly putty on a table. Over the
course of an hour, it flows slowly outward, taking the shape of the table
surface — clearly a liquid. Strike it sharply with a hammer and it shatters —
clearly a solid. The material has not changed. The timescale of the
interaction has.</p>
<p>The same principle applies to glass (contrary to popular myth, medieval window
glass is not thicker at the bottom because it has flowed — the variation is
from the manufacturing process, and the relaxation time of soda-lime glass at
room temperature is of order $10^{23}$ years — but at elevated temperatures
near the glass transition, silicate glass flows readily). It applies
to mantle rock, which is solid on the scale of earthquake waves and liquid on
the scale of continental drift. It applies to pitch, to ice sheets, to asphalt
on a hot day.</p>
<p>The formal tool for capturing this is the <strong>Deborah number</strong>.</p>
<hr>
<h2 id="the-deborah-number">The Deborah Number</h2>
<p>The Deborah number was introduced by Marcus Reiner in 1964, in a short note
in <em>Physics Today</em> (Reiner 1964). It is defined as:</p>
$$\mathrm{De} = \frac{\tau}{T},$$<p>where $\tau$ is the <strong>relaxation time</strong> of the material — roughly, the
characteristic time over which it can rearrange its internal structure and
relieve stress — and $T$ is the <strong>observation time</strong> or the timescale of the
imposed deformation.</p>
<ul>
<li>$\mathrm{De} \ll 1$: The material relaxes quickly relative to the timescale
of observation. Internal stresses are continuously relieved. The material
behaves as a <strong>fluid</strong>.</li>
<li>$\mathrm{De} \gg 1$: The material relaxes slowly relative to the observation
timescale. Internal stresses persist. The material behaves as a <strong>solid</strong>.</li>
<li>$\mathrm{De} \sim 1$: The material is in a viscoelastic regime — partly
fluid, partly solid, exhibiting time-dependent behaviour that is neither.</li>
</ul>
<p>The name comes from the prophetess Deborah, who sang in Judges 5:5: <em>&ldquo;The
mountains flowed before the Lord.&rdquo;</em> At the timescale of a divine perspective,
mountains are liquid. At the timescale of a human lifetime, they are not.
Reiner&rsquo;s point was that the solid-liquid distinction is not a property of
the material but of the relationship between the material&rsquo;s internal
dynamics and the observer&rsquo;s timescale.</p>
<p>For Newtonian fluids (water, air at ordinary conditions), $\tau \to 0$ and
$\mathrm{De} \to 0$ for any finite observation time — they are always liquid.
For a perfectly elastic solid (an ideal spring), $\tau \to \infty$ and
$\mathrm{De} \to \infty$ for any finite observation time — always solid. Real
materials lie between these extremes.</p>
<hr>
<h2 id="the-maxwell-viscoelastic-model">The Maxwell Viscoelastic Model</h2>
<p>The simplest model of a material with a finite relaxation time is the Maxwell
element: a spring (elastic, spring constant $G$) in series with a dashpot
(viscous, viscosity $\eta$). Under a step stress $\sigma_0$ applied at time
$t = 0$, the strain evolves as:</p>
$$\epsilon(t) = \frac{\sigma_0}{G} + \frac{\sigma_0}{\eta}\,t,$$<p>where $\tau = \eta / G$ is the Maxwell relaxation time. The first term is the
instantaneous elastic deformation of the spring; the second is the linear
viscous creep of the dashpot. For $t \ll \tau$, the elastic strain dominates
and the material behaves as a solid; for $t \gg \tau$, the viscous flow
dominates and the material behaves as a liquid. The material &ldquo;decides&rdquo; whether
to be solid or liquid depending on the ratio of $\tau$ to the duration of the
applied stress — which is precisely the Deborah number.</p>
<p>The <strong>creep compliance</strong> $J(t) = \epsilon(t)/\sigma_0 = t/\eta + 1/G$ grows
linearly with time for $t \gg \tau$, confirming liquid behaviour on long
timescales. The <strong>relaxation modulus</strong> $G(t) = \sigma(t)/\epsilon_0 = G
e^{-t/\tau}$ decays exponentially to zero, confirming that the material
cannot sustain a permanent stress — again, liquid behaviour on long timescales.</p>
<hr>
<h2 id="on-the-rheology-of-cats">On the Rheology of Cats</h2>
<p>In 2014, Marc-Antoine Fardin, a physicist at the ENS Lyon,
published &ldquo;On the Rheology of Cats&rdquo; in the <em>Rheology Bulletin</em> 83(2), 16–17.
The paper asked whether cats satisfy the defining rheological criterion for
liquids, using the Deborah number as the test. Fardin was awarded the 2017
Ig Nobel Prize in Physics — which is awarded for research that &ldquo;makes you
laugh, then makes you think&rdquo; — for this work.</p>
<p>The paper is not a joke. It is standard rheology applied to an unusual material,
with appropriately hedged conclusions and correct citations to the primary
literature on viscoelastic flow. The humour is in the application; the physics
is serious.</p>
<h3 id="estimating-the-cats-relaxation-time">Estimating the Cat&rsquo;s Relaxation Time</h3>
<p>The relaxation time $\tau$ of a cat is the time scale over which the cat&rsquo;s
body deforms to fill a container. This is observable. A cat placed near a
suitable container — a salad bowl, a cardboard box, a bathroom sink —
adopts a conformed shape on a timescale of roughly 5–30 seconds. The initial
posture (stiff, alert) gives way to a relaxed conformation as the cat
assesses the container and adjusts. Fardin estimated $\tau \approx 1$–$30$
seconds, with the exact value depending on the container&rsquo;s attractiveness
to the specific cat.</p>
<p>This is the material&rsquo;s characteristic relaxation time. The fact that it is
finite — that the cat does eventually conform to the container — is the
essential observation.</p>
<h3 id="computing-the-deborah-number-for-various-situations">Computing the Deborah Number for Various Situations</h3>
<p><strong>Scenario 1: Cat in a sink.</strong>
A cat taking ten minutes to settle into a bathroom sink. Observation time
$T = 600\,\mathrm{s}$, relaxation time $\tau \approx 5\,\mathrm{s}$.</p>
$$\mathrm{De}_\mathrm{sink} = \frac{5}{600} \approx 0.008 \ll 1.$$<p>The cat is unambiguously a <strong>liquid</strong>.</p>
<p><strong>Scenario 2: Cat in a cardboard box.</strong>
Conformation over approximately 30 minutes, $\tau \approx 20\,\mathrm{s}$.</p>
$$\mathrm{De}_\mathrm{box} = \frac{20}{1800} \approx 0.011 \ll 1.$$<p><strong>Liquid.</strong></p>
<p><strong>Scenario 3: Cat dropping from a bookshelf.</strong>
Contact time during a jump approximately $T \approx 0.05\,\mathrm{s}$,
relaxation time still $\tau \approx 5\,\mathrm{s}$.</p>
$$\mathrm{De}_\mathrm{jump} = \frac{5}{0.05} = 100 \gg 1.$$<p><strong>Solid.</strong> The cat does not deform into the shape of the bookshelf during the
jump; it rebounds elastically.</p>
<p><strong>Scenario 4: Cat startled by a loud noise.</strong>
Reaction time $T \approx 0.3\,\mathrm{s}$, $\tau \approx 5\,\mathrm{s}$.</p>
$$\mathrm{De}_\mathrm{startle} = \frac{5}{0.3} \approx 17 \gg 1.$$<p><strong>Solid.</strong> On short timescales, cats behave as elastic materials — they spring,
they bounce, they do not flow.</p>
<p>The cat is neither permanently solid nor permanently liquid. It is a
<strong>viscoelastic material</strong> whose phase classification depends on the timescale
of the interaction. This is not a loose analogy; it is the definition of
viscoelasticity.</p>
<hr>
<h2 id="non-newtonian-behaviour-and-flow-instabilities">Non-Newtonian Behaviour and Flow Instabilities</h2>
<p>Fardin noted an additional complication: cat flow is not Newtonian. A Newtonian
fluid has a viscosity $\eta$ that is independent of the applied shear rate
$\dot\gamma$. Many real materials are <strong>shear-thinning</strong> (viscosity decreases
with increasing shear rate — ketchup, blood, many polymer solutions) or
<strong>shear-thickening</strong> (viscosity increases with increasing shear rate —
cornstarch suspension, some dense suspensions). Cats, Fardin observed, appear
to be shear-thinning: the more rapidly you attempt to move a relaxed cat from
its current position, the more &ldquo;liquid&rdquo; (accommodating, compliant) it becomes,
up to a point at which the cat transitions to solid behaviour (claws, teeth).</p>
<p>This is, formally, the behaviour of a <strong>yield-stress fluid</strong>: a material that
behaves as a solid below a critical stress $\sigma_y$ and flows above it. The
Herschel–Bulkley model describes such fluids:</p>
$$\sigma = \sigma_y + k \dot\gamma^n, \quad \sigma > \sigma_y,$$<p>where $k$ is the flow consistency index and $n < 1$ for shear-thinning. The
challenge of fitting $k$, $n$, and $\sigma_y$ for a specific cat is
experimental, and Fardin acknowledged this was left to future work.</p>
<p>The <strong>Deborah number</strong> and the <strong>yield stress</strong> together provide a two-parameter
phase diagram for cat rheology:</p>
<ul>
<li>Low stress, short timescale: solid (De ≫ 1 or σ &lt; σ_y)</li>
<li>Low stress, long timescale: liquid (De ≪ 1)</li>
<li>High stress: yield, followed by flow</li>
</ul>
<hr>
<h2 id="flow-instabilities-the-rayleigh-plateau-connection">Flow Instabilities: The Rayleigh-Plateau Connection</h2>
<p>Fardin also noted that cats confined to containers thinner than their body
diameter can exhibit flow instabilities. A cat attempting to fit into a glass
too narrow for its body will sometimes adopt a helical or coiled configuration —
an instability reminiscent of the <strong>Rayleigh–Plateau instability</strong> of a liquid
jet.</p>
<p>The Rayleigh–Plateau instability occurs when a cylindrical fluid jet of radius
$r_0$ is subject to perturbations of wavelength $\lambda > 2\pi r_0$. Modes
with wavelength longer than the cylinder&rsquo;s circumference are unstable and grow,
breaking the jet into droplets. The dispersion relation for growth rate $\sigma$
as a function of wavenumber $k = 2\pi/\lambda$ (for an inviscid jet) is:</p>
$$\sigma^2 = \frac{\gamma}{\rho r_0^3}\, k r_0 \bigl(1 - k^2 r_0^2\bigr)
I_1(kr_0)/I_0(kr_0),$$<p>where $\gamma$ is surface tension and $I_0, I_1$ are modified Bessel functions.
The analogy with a cat is inexact — surface tension is not the dominant
restoring force — but the qualitative instability mechanism (a long cylinder of
material is unstable to perturbations whose wavelength exceeds the cylinder&rsquo;s
circumference) appears to apply, suggesting that very elongated cats in very
narrow containers should be unstable to coiling. This is, again, left to future
experimental work.</p>
<hr>
<h2 id="why-the-deborah-number-matters-outside-of-cat-physics">Why the Deborah Number Matters (Outside of Cat Physics)</h2>
<p>The Deborah number is not a curiosity; it is a central dimensionless number
in engineering and materials science.</p>
<p><strong>Polymer processing</strong>: The flow of polymer melts through injection-moulding
channels involves De in the range $10^{-2}$–$10^2$. Too high a De leads to
elastic instabilities, melt fracture, and surface defects in the finished part.</p>
<p><strong>Blood rheology</strong>: Blood is a non-Newtonian viscoelastic fluid. In the large
arteries (low shear rate), red blood cells aggregate into <em>rouleaux</em> and
blood behaves as a shear-thinning fluid. In the capillaries (high shear rate),
rouleaux break up and individual cells deform to fit through vessels smaller
than their resting diameter — liquid behaviour on short length scales.</p>
<p><strong>Geophysics</strong>: The mantle is an elastic solid for seismic waves ($T \sim$
seconds, De ≫ 1) and a viscous fluid for convection ($T \sim 10^8$–$10^9$
years, De ≪ 1). The same material. Different Deborah numbers.</p>
<p><strong>Glaciology</strong>: Ice is an elastic solid for rapid fracture (calving of icebergs)
and a viscous fluid for glacier flow. The transition occurs at timescales of
years to decades, depending on temperature and stress.</p>
<p>The cat is in good company.</p>
<hr>
<h2 id="references">References</h2>
<ul>
<li>
<p>Fardin, M.-A. (2014). On the rheology of cats. <em>Rheology Bulletin</em>, 83(2),
16–17.</p>
</li>
<li>
<p>Reiner, M. (1964). The Deborah number. <em>Physics Today</em>, 17(1), 62.
<a href="https://doi.org/10.1063/1.3051374">https://doi.org/10.1063/1.3051374</a></p>
</li>
<li>
<p>Barnes, H.A., Hutton, J.F., &amp; Walters, K. (1989). <em>An Introduction to
Rheology.</em> Elsevier (Rheology Series, Vol. 3).</p>
</li>
<li>
<p>Bird, R.B., Armstrong, R.C., &amp; Hassager, O. (1987). <em>Dynamics of Polymeric
Liquids, Vol. 1: Fluid Mechanics</em> (2nd ed.). Wiley-Interscience.</p>
</li>
<li>
<p>Eggers, J. (1997). Nonlinear dynamics and breakup of free-surface flows.
<em>Reviews of Modern Physics</em>, 69(3), 865–930.
<a href="https://doi.org/10.1103/RevModPhys.69.865">https://doi.org/10.1103/RevModPhys.69.865</a></p>
</li>
</ul>
<hr>
<h2 id="changelog">Changelog</h2>
<ul>
<li><strong>2025-12-15</strong>: Fixed Deborah number in summary from 0.08 to 0.008 (matching the body calculation: 5/600 = 0.00833).</li>
<li><strong>2025-12-15</strong>: Corrected Fardin&rsquo;s institutional affiliation from &ldquo;Paris Diderot University&rdquo; to &ldquo;ENS Lyon&rdquo; — his affiliation on the 2014 <em>Rheology Bulletin</em> paper is Université de Lyon / ENS Lyon (CNRS UMR 5672). He moved to Paris Diderot later in 2014, after the paper was published.</li>
</ul>
]]></content:encoded>
    </item>
    <item>
      <title>Zero Angular Momentum: The Falling Cat and the Geometry of Shape Space</title>
      <link>https://sebastianspicker.github.io/posts/falling-cat-geometric-phase/</link>
      <pubDate>Tue, 03 Oct 2023 00:00:00 +0000</pubDate>
      <guid>https://sebastianspicker.github.io/posts/falling-cat-geometric-phase/</guid>
      <description>A cat dropped upside-down rotates 180° and lands on its feet, despite having zero angular momentum throughout. This is not a trick and not a violation of physics. The explanation took physicists from 1894 to 1993 to fully work out, and the answer — a geometric phase arising from the holonomy of a fiber bundle — is the same mathematics that governs the Berry phase in quantum mechanics and the Aharonov-Bohm effect in electrodynamics. We adopted two strays this year. They fall beautifully.</description>
      <content:encoded><![CDATA[<p><em>We adopted two stray cats in 2023. They had been living under a garden shed and
had strong opinions about most things, including the correct height from which to
leap onto a bookshelf and whether landing was optional. They are indoor cats now,
for health reasons — a vet&rsquo;s recommendation they find unconvincing but have largely
accepted. Watching one of them drop
from a windowsill — always feet-first, always orientated correctly, from heights
that would leave me reconsidering my life choices — I found myself thinking about
a problem I had first encountered in a mechanics course and had never fully
resolved to my satisfaction.</em></p>
<p><em>How does a cat rotate with zero angular momentum?</em></p>
<hr>
<h2 id="the-problem">The Problem</h2>
<p>When a cat is dropped from an inverted position — upside-down, held by a
practised experimenter, then released — it rotates approximately 180° and
lands on its feet. The drop takes around 0.3 seconds. The cat begins with
negligible angular momentum (the experimenter can release it with almost no
spin), and there are no external torques during free fall. By conservation of
angular momentum, the total angular momentum of the cat must remain constant
throughout the fall.</p>
<p>The total angular momentum is therefore approximately zero throughout the
fall.</p>
<p>And yet the cat rotates 180°.</p>
<p>This is the falling cat problem. It was first documented quantitatively by
Étienne-Jules Marey in 1894 using chronophotography — among the first
high-speed photography of any biological motion — and it has
occupied physicists, mathematicians, neuroscientists, and roboticists ever
since.</p>
<p>The problem is not exotic. Every cat owner has seen it. What requires
explanation is why our intuitions about angular momentum fail here, and what
replaces them.</p>
<hr>
<h2 id="why-the-obvious-answers-do-not-work">Why the Obvious Answers Do Not Work</h2>
<p>There are two naive explanations for the cat&rsquo;s righting reflex, both wrong.</p>
<p><strong>Explanation 1: The cat uses initial angular momentum.</strong> The experimenter
gives the cat a small spin before releasing it; the cat amplifies this to
achieve the full 180°. This fails because controlled experiments (and Marey&rsquo;s
original photographs) confirm that cats can right themselves even when
released with zero initial spin. Careful experimenters have verified this
explicitly.</p>
<p><strong>Explanation 2: The cat pushes against the air.</strong> A falling cat could, in
principle, use aerodynamic forces to push against the air and generate a
reaction. This fails because the angular impulse from air drag over 0.3
seconds is far too small to account for the observed 180° rotation. Marey&rsquo;s
chronophotographs already showed that the motion begins immediately on
release, before air resistance could contribute meaningfully.</p>
<p>Both explanations appeal to external torques. The correct explanation requires
none.</p>
<hr>
<h2 id="marey-and-the-photographic-evidence">Marey and the Photographic Evidence</h2>
<p>Étienne-Jules Marey published his chronophotographic sequence of a falling
cat in <em>La Nature</em> on 10 November 1894. The images, taken at 60 frames per
second, show the following clearly:</p>
<ol>
<li>The front and rear halves of the cat move <em>asymmetrically</em>. The front half
rotates in one direction; the rear half rotates by a smaller angle in the
opposite direction.</li>
<li>The cat pulls its front legs in close to its body (reducing the moment of
inertia of the front half) while extending its rear legs (increasing the
moment of inertia of the rear half).</li>
<li>The front half then rotates rapidly (large angle, small moment of inertia);
the rear half rotates slowly in the opposite direction (small angle, large
moment of inertia).</li>
<li>The cat then extends its front legs and pulls in its rear legs, and reverses
the process.</li>
</ol>
<p>The net effect: the cat&rsquo;s body orientation rotates by 180° even though the
<em>total</em> angular momentum — computed as the sum of both halves — remains
constant. The key word is <em>sum</em>. Individual parts can exchange angular momentum
through internal torques; the sum is conserved.</p>
<p>This mechanism — internal redistribution of angular momentum without changing
its total — is correct but not complete. It explains <em>that</em> rotation is
possible, not <em>how much</em> rotation is achieved per cycle of shape change. For
that, we need the mathematics.</p>
<hr>
<h2 id="kane-and-scher-the-two-cylinder-model">Kane and Scher: The Two-Cylinder Model</h2>
<p>The first rigorous mechanical model was published by T.R. Kane and M.P. Scher
in 1969 (<em>International Journal of Solids and Structures</em> 5, 663–670).</p>
<p>They modelled the cat as two rigid axisymmetric cylinders — a front half and
a rear half — connected at a joint that allows relative bending and twisting.
The joint constraint imposes that the relative twist between the two halves is
zero (a &ldquo;no-twist&rdquo; condition: the cylinders cannot spin relative to each other
at their connection). The total angular momentum of the system is held fixed
at zero.</p>
<p>Let the two cylinders have moments of inertia $I_1$ and $I_2$ about their
symmetry axes, and let $\phi$ be the bend angle between them and $\psi$ the
twist angle. The zero-angular-momentum constraint, combined with the no-twist
condition, gives a system of equations that can be integrated numerically to
find the net body rotation as a function of the shape-change trajectory
$(\phi(t), \psi(t))$.</p>
<p>Kane and Scher showed that a specific sequence of shape changes — one complete
cycle in the $(\phi, \psi)$ plane — produces a net rotation of approximately
90–100°. A second cycle gives the rest. The calculation was the first to
confirm, from mechanics alone, that the righting manoeuvre requires no external
torques and is entirely consistent with conservation of angular momentum.</p>
<p>What the Kane–Scher model does not explain is <em>why</em> the net rotation per cycle
depends on the area enclosed by the trajectory in shape space — or why the
same mathematical structure appears in quantum mechanics. For that, we need
Montgomery&rsquo;s formulation.</p>
<hr>
<h2 id="montgomery-fiber-bundles-and-geometric-holonomy">Montgomery: Fiber Bundles and Geometric Holonomy</h2>
<p>In 1993, Richard Montgomery published a reformulation of the falling cat problem
using gauge theory (<em>Dynamics and Control of Mechanical Systems</em>, Fields
Institute Communications, AMS, pp. 193–218). The reformulation is the
definitive mathematical treatment, and it connects the cat to one of the deepest
structures in modern physics.</p>
<h3 id="the-configuration-space">The Configuration Space</h3>
<p>The full configuration space of the cat — the space of all possible positions
and orientations — is</p>
$$Q = SO(3) \times \mathcal{S},$$<p>where $SO(3)$ is the rotation group (describing the cat&rsquo;s overall orientation
in space) and $\mathcal{S}$ is the <em>shape space</em> (describing the internal
geometry: the bend angle, the twist, the position of each limb relative to the
body).</p>
<p>The angular momentum constraint $\mathbf{L} = 0$ defines a <em>horizontal
distribution</em> on $Q$ — a preferred subspace of tangent vectors at each point
that correspond to shape changes at zero angular momentum. This distribution is
not integrable (it does not come from a foliation), which is the mathematical
signature that holonomy is possible.</p>
<h3 id="the-fiber-bundle">The Fiber Bundle</h3>
<p>The projection</p>
$$\pi \colon Q \to \mathcal{S}, \qquad (R, s) \mapsto s,$$<p>makes $Q$ into a principal fiber bundle over $\mathcal{S}$ with structure group
$SO(3)$. The fiber above each shape $s \in \mathcal{S}$ is the set of all
orientations the cat can have with that shape.</p>
<p>A <em>connection</em> on this bundle is a rule for &ldquo;lifting&rdquo; paths in the base
$\mathcal{S}$ to horizontal paths in the total space $Q$ — that is, paths
along which the angular momentum constraint is satisfied. This connection
$\mathcal{A}$ is a one-form on $\mathcal{S}$ taking values in the Lie algebra
$\mathfrak{so}(3)$.</p>
<h3 id="holonomy-the-geometric-phase">Holonomy: The Geometric Phase</h3>
<p>When the cat executes a closed loop $\gamma$ in shape space — a sequence of
shape changes that returns it to its initial shape — the <em>holonomy</em> of the
connection $\mathcal{A}$ around $\gamma$ gives the net rotation:</p>
$$R_\gamma = \mathrm{Hol}_\mathcal{A}(\gamma) \in SO(3).$$<p>For the full non-Abelian case ($SO(3)$), the holonomy is a path-ordered
exponential along $\gamma$ and its relationship to the curvature involves
non-Abelian corrections. But the essential geometric intuition is captured
by the Abelian case — rotation about a single axis — where Stokes&rsquo;s theorem
gives the net rotation directly:</p>
$$\theta_\gamma = \iint_{\Sigma} F,$$<p>where $\Sigma$ is a surface bounded by $\gamma$ and $F = d\mathcal{A}$
is the curvature 2-form. The cat&rsquo;s net rotation per cycle is the integral
of the curvature over the area enclosed by its shape-change loop in
$\mathcal{S}$. For small loops, the curvature $F_\mathcal{A} = d\mathcal{A}</p>
<ul>
<li>\mathcal{A} \wedge \mathcal{A}$ determines the holonomy to leading order
in both the Abelian and non-Abelian cases.</li>
</ul>
<p>The rotation is <em>geometric</em>: it depends on the shape of the loop, not on the
speed at which the loop is traversed. A cat executing the same shape-change
sequence twice as fast achieves the same rotation in half the time.</p>
<hr>
<h2 id="the-connection-to-berry-phase">The Connection to Berry Phase</h2>
<p>The gauge structure of the falling cat problem is not an isolated curiosity.
It is the same mathematical structure that governs several central phenomena
in modern physics.</p>
<p><strong>The Berry phase</strong> (Berry 1984, <em>Proceedings of the Royal Society A</em>) arises
when a quantum system is transported adiabatically around a closed loop $C$ in
parameter space. The state acquires a phase</p>
$$\gamma_B = \oint_C \mathbf{A} \cdot d\mathbf{R},$$<p>where $\mathbf{A} = i\langle n(\mathbf{R}) | \nabla_\mathbf{R} | n(\mathbf{R}) \rangle$
is the Berry connection — a gauge field on parameter space. The Berry phase is
the holonomy of this connection, which is to say: the cat righting itself and
a quantum state accumulating a geometric phase are instances of the <em>same
mathematical theorem</em>.</p>
<p>Shapere and Wilczek (1989) made this connection explicit for deformable bodies,
noting that the net rotation of a swimming microorganism or a falling cat is
the holonomy of a gauge connection on shape space — exactly the Berry phase,
expressed in the language of classical mechanics.</p>
<p><strong>The Foucault pendulum</strong> precesses at a rate of $2\pi\sin\phi$ per sidereal
day, where $\phi$ is the latitude. The holonomy of the Levi-Civita connection
on $S^2$ for parallel transport around the circle of latitude is the solid
angle of the enclosed polar cap, $\Omega = 2\pi(1 - \sin\phi)$. The
lab-frame precession $2\pi\sin\phi = 2\pi - \Omega$ is the complementary
angle — the two sum to a full rotation because the local frame itself
completes one circuit per sidereal day. It is another geometric phase.</p>
<p><strong>The Aharonov-Bohm effect</strong> (1959) produces a phase shift for electrons
circling a solenoid, even when the electrons travel only through field-free
regions. The phase is the holonomy of the electromagnetic vector potential
$\mathbf{A}$ around the loop — a Berry phase for the electromagnetic field.</p>
<p>All four phenomena — the falling cat, the Berry phase, the Foucault pendulum,
the Aharonov-Bohm effect — are manifestations of the same structure: a
connection on a fiber bundle, and holonomy as the geometric consequence of
traversing a closed loop.</p>
<p>Batterman (2003, <em>Studies in History and Philosophy of Modern Physics</em> 34,
527–557) gives a particularly clear account of this unification, drawing out
the common mathematical skeleton and its physical implications.</p>
<hr>
<h2 id="high-rise-syndrome-terminal-velocity-and-the-parachute-cat">High-Rise Syndrome: Terminal Velocity and the Parachute Cat</h2>
<p>There is a grounding empirical footnote to the elegant geometry above. Whitney
and Mehlhaff (1987, <em>Journal of the American Veterinary Medical Association</em>
191, 1399–1403) analysed 132 cats brought to a Manhattan veterinary clinic after
falling from buildings of two to thirty-two stories. Their finding was
counterintuitive:</p>
<p>Cats falling from above seven stories had a <em>lower</em> injury rate than cats
falling two to six stories. Overall, 90% of the cats in the study survived,
with injuries paradoxically less severe at greater heights.</p>
<p>The explanation involves two phases. Below seven stories, the cat is still
accelerating: it is tense, its legs are extended to brace for impact, and it
absorbs the force of landing poorly. Above seven stories, the cat reaches
terminal velocity — approximately $100\,\mathrm{km/h}$ for a falling cat — and
then, apparently, <em>relaxes</em>. The vestibular system, having identified that the
fall is not ending imminently, switches from the righting reflex to a
parachute posture: legs spread horizontally, body flattened, increasing the
cross-sectional area and hence air resistance.</p>
<p>Terminal velocity is reached when the drag force equals the gravitational force:</p>
$$mg = \frac{1}{2} C_D \rho A v_t^2, \qquad
v_t = \sqrt{\frac{2mg}{C_D \rho A}}.$$<p>For a spread-eagle cat ($m \approx 4\,\mathrm{kg}$, $A \approx 0.06\,\mathrm{m}^2$,
$C_D \approx 1.0$, $\rho_\mathrm{air} \approx 1.2\,\mathrm{kg/m}^3$):</p>
$$v_t \approx \sqrt{\frac{2 \times 4 \times 9.8}{1.0 \times 1.2 \times 0.06}}
\approx 33\,\mathrm{m/s} \approx 120\,\mathrm{km/h}.$$<p>(The exact value depends on posture and fur drag; empirical estimates for
cats in the parachute posture are lower, roughly $25$–$30\,\mathrm{m/s}$,
because the effective area increases when the limbs are spread.)</p>
<p>A human in free-fall has terminal velocity around $55\,\mathrm{m/s}$
($200\,\mathrm{km/h}$) — faster, because the mass-to-area ratio is higher.
The cat, with its low mass and high drag relative to body weight, hits a
gentler terminal velocity and distributes the impact more effectively.</p>
<p>The study is sometimes cited as evidence that cats are invincible. A significant
caveat is <strong>survivorship bias</strong>: cats that died on impact were likely not brought
to the veterinary clinic, so the dataset underrepresents fatal outcomes,
especially for higher falls. The apparent decrease in injury rate above seven
stories may partly reflect the fact that the most severely injured cats from
those heights never entered the study. The aerodynamic posture explanation is
plausible, but the data do not cleanly separate it from the sampling bias.</p>
<hr>
<h2 id="robotics-and-spacecraft">Robotics and Spacecraft</h2>
<p>The falling cat problem has practical applications beyond veterinary statistics.</p>
<p><strong>Spacecraft attitude control</strong>: Astronauts in free fall can change their
body orientation without thrusters, using the same gauge-theoretic mechanism
as the cat. NASA and ESA have studied cat-inspired reorientation manoeuvres
for astronauts and satellites.</p>
<p><strong>Robotics</strong>: The two-cylinder model inspired early robot designs capable of
reorienting in free fall — useful for robots deployed from aircraft or
spacecraft. Subsequent work (including a 2022 review in <em>IEEE Transactions on
Robotics</em>) has produced legged robots that can right themselves after being
knocked over using shape-change sequences derived from the Montgomery connection.</p>
<p><strong>Gymnastics and diving</strong>: Human athletes performing somersaults and twists
exploit the same gauge structure, though without articulating the mathematics.
A tuck increases rotation rate (smaller $I$, constant $L$ → larger $\omega$);
a layout decreases it. Changing the tuck–layout timing mid-rotation produces
a net twist — holonomy in the shape space of a human body.</p>
<hr>
<h2 id="the-view-from-a-windowsill">The View from a Windowsill</h2>
<p>My cats have no opinion about fiber bundles. When one of them drops from the
top of the bookcase, she is not solving the variational problem</p>
$$\min_{\gamma \in \Omega} \int_\gamma |\dot{s}|^2 \, dt,
\quad \text{subject to } \mathrm{Hol}_\mathcal{A}(\gamma) = R_{180°},$$<p>she is executing a motor program refined over millions of years of feline
evolution. The vestibular system provides continuous feedback on body
orientation; the cerebellum coordinates the shape-change sequence; the whole
manoeuvre is over in a third of a second.</p>
<p>What physics tells us is that the manoeuvre is <em>possible</em> — that no law of
nature forbids a body with zero angular momentum from reorienting — and gives
the precise geometric reason: the curvature of a connection on shape space is
non-zero, which means the holonomy of closed loops is non-trivial.</p>
<p>The same curvature that allows a cat to right itself allows a quantum state to
accumulate a geometric phase, allows the Foucault pendulum to precess, and
allows the Aharonov-Bohm effect to shift an interference fringe without a local
field. These are not analogies. They are the same theorem, applied to different
physical systems in different mathematical languages.</p>
<p>I find this more remarkable than the cat.</p>
<hr>
<h2 id="references">References</h2>
<ul>
<li>
<p>Batterman, R.W. (2003). Falling cats, parallel parking, and polarized light.
<em>Studies in History and Philosophy of Modern Physics</em>, 34(4), 527–557.
<a href="https://doi.org/10.1016/S1355-2198(03)00062-5">https://doi.org/10.1016/S1355-2198(03)00062-5</a></p>
</li>
<li>
<p>Berry, M.V. (1984). Quantal phase factors accompanying adiabatic changes.
<em>Proceedings of the Royal Society A</em>, 392, 45–57.
<a href="https://doi.org/10.1098/rspa.1984.0023">https://doi.org/10.1098/rspa.1984.0023</a></p>
</li>
<li>
<p>Gbur, G.J. (2019). <em>Falling Felines and Fundamental Physics.</em> Yale University
Press.</p>
</li>
<li>
<p>Kane, T.R., &amp; Scher, M.P. (1969). A dynamical explanation of the falling cat
phenomenon. <em>International Journal of Solids and Structures</em>, 5(7), 663–670.
<a href="https://doi.org/10.1016/0020-7683(69)90086-9">https://doi.org/10.1016/0020-7683(69)90086-9</a></p>
</li>
<li>
<p>Marey, É.-J. (1894). Des mouvements que certains animaux exécutent pour
retomber sur leurs pieds lorsqu&rsquo;ils sont précipités d&rsquo;un lieu élevé. <em>La
Nature</em>, 10 November 1894.</p>
</li>
<li>
<p>Montgomery, R. (1993). Gauge theory of the falling cat. In M. Enos (Ed.),
<em>Dynamics and Control of Mechanical Systems</em> (Fields Institute Communications,
Vol. 1, pp. 193–218). American Mathematical Society.</p>
</li>
<li>
<p>Shapere, A., &amp; Wilczek, F. (Eds.). (1989). <em>Geometric Phases in Physics.</em>
World Scientific.</p>
</li>
<li>
<p>Whitney, W.O., &amp; Mehlhaff, C.J. (1987). High-rise syndrome in cats. <em>Journal
of the American Veterinary Medical Association</em>, 191(11), 1399–1403.</p>
</li>
</ul>
<hr>
<h2 id="changelog">Changelog</h2>
<ul>
<li><strong>2025-12-15</strong>: Corrected the Marey publication date from 22 November 1894 to 10 November 1894 (in text and in reference). Updated the Whitney &amp; Mehlhaff (1987) statistics to reflect that the 90% survival rate applies to all cats in the study, as reported in the paper, rather than specifically to those falling from above seven stories.</li>
</ul>
]]></content:encoded>
    </item>
  </channel>
</rss>
