---
title: Stellar Pulsation and the Instability Strip
draft: false
module: Stellar Evolution
moduleNumber: 7
lessonNumber: 4
order: 704
summary: >
  Radial pulsation is a standing sound wave whose period scales inversely with the
  square root of the mean density. The kappa mechanism, an opacity valve seated in
  the helium partial-ionization zone, turns a star into a heat engine that pumps the
  oscillation. Stars in the instability strip pulsate as Cepheids, RR Lyrae, and
  Mira variables, and the Cepheid period–luminosity relation calibrates the distance
  ladder.
topics: [Stellar Evolution]
sources:
  - book: Carroll & Ostlie
    ref: "Ch. 14 — Stellar Pulsation; §14.1 Observations of Pulsating Stars, §14.2 The Physics of Stellar Pulsation"
  - book: Maoz
    ref: "Ch. 2 — Stars: Basic Observations"
---

Some stars vary in brightness because they pulsate: the whole star expands and
contracts, brightening as it heats and dims as it cools. The oscillation is a
standing sound wave, its period fixed by the mean density, and it is sustained
against dissipation by an opacity valve in a partial-ionization layer of the
envelope. Only stars occupying a narrow, nearly vertical strip of the
Hertzsprung–Russell diagram pulsate this way, and their period–luminosity relation
makes them the measuring rods of extragalactic distance.[^co-puls]

## The period–mean-density relation

A radial pulsation is a standing acoustic wave in the star. Its period is the time a
sound wave takes to cross the star and return, so it depends on the sound speed and
the radius. Model the star crudely as a sphere of uniform density $\bar\rho$
oscillating adiabatically. A displaced shell feels a restoring pressure whose
adiabatic response is set by the first adiabatic exponent $\Gamma_1$, and the
linearized equation of motion has the fundamental period

$$
\Pi = \frac{2\pi}{\sqrt{\tfrac{4}{3}\pi G\bar\rho\,(3\Gamma_1 - 4)}}.
$$

For $\Gamma_1 = 5/3$ the factor $3\Gamma_1 - 4 = 1$, and this reduces to
$\Pi = \sqrt{3\pi/(G\bar\rho)}$. The essential content is the scaling

$$
\Pi\sqrt{\bar\rho} = Q = \text{const},
\qquad
\Pi\sqrt{\frac{\bar\rho}{\bar\rho_\odot}} = Q,
$$

the **period–mean-density relation**. The **pulsation constant** $Q$ is nearly the
same for all stars pulsating in the fundamental mode, about $0.04\ \text{days}$ for
classical Cepheids. Denser stars pulsate faster; the tenuous, luminous supergiants
have the longest periods. Note the requirement $\Gamma_1 > 4/3$ for a real period:
below it the restoring force vanishes and the star is dynamically unstable rather
than oscillatory.

The fundamental mode has the whole envelope moving in phase, with the center fixed
and the surface at maximum displacement. Overtones have one or more internal nodes
where the gas stays still while layers on either side move oppositely, and they have
shorter periods. Most Cepheids pulsate in the fundamental; many RR Lyrae stars
pulsate in an overtone, giving the double-peaked period distribution seen in globular
clusters.

$$
% caption: Radial pulsation modes shown as displacement against radius; the
% fundamental mode has all layers moving in phase with a single antinode at the
% surface, while the overtone has an interior node separating oppositely moving
% layers.
\begin{tikzpicture}[scale=1.0, font=\footnotesize]
\definecolor{acc}{HTML}{4A6FA5}
\draw[->, black] (0,2.3) -- (8.8,2.3) node[right, black!70] {radius};
\draw[->, black] (0,0.4) -- (0,4.4) node[above, black!70] {displacement};
\node[black, anchor=north east] at (0,2.3) {center};
\node[black, anchor=north] at (8.2,2.3) {surface};
% fundamental: monotone rise, all in phase
\draw[acc, very thick] (0,2.3) .. controls (3.5,2.5) and (6.0,3.4) .. (8.2,4.0);
\node[acc, anchor=south east] at (7.6,3.9) {fundamental};
% overtone: rises, node, dips below
\draw[black, very thick, densely dashed] (0,2.3) .. controls (2.0,3.2) and (3.0,3.0) .. (4.2,2.3)
  .. controls (5.4,1.6) and (6.6,0.9) .. (8.2,1.0);
\node[black, anchor=north] at (2.5,3.25) {overtone};
% node marker
\fill[black] (4.2,2.3) circle (1.6pt);
\node[black, anchor=south] at (4.2,2.35) {node};
\end{tikzpicture}
$$

## The kappa mechanism

A pulsation left to itself would damp as sound waves dissipate. To persist, the star
must do net work on the oscillation over each cycle, acting as a heat engine that
absorbs heat when compressed and releases it when expanded. Ordinary stellar material
does the opposite. Compression raises the temperature, and in the Kramers regime the
opacity $\kappa \propto \rho\,T^{-7/2}$ falls with temperature, so a compressed layer
becomes more transparent and lets heat leak out exactly when the engine would need to
retain it. Across most of a star, pulsations are damped.

The exception is a **partial-ionization zone**. Where a species is half-ionized,
compression does not raise the temperature much, because the energy goes into further
ionization rather than into heating the gas. With the temperature nearly fixed, the
density increase makes the opacity **rise** on compression. The layer dams the
outward flux while compressed, storing heat, and releases it on expansion. This is the
**kappa mechanism** (aided by the associated **gamma mechanism**, the same energy
diversion into ionization). The dominant valve in classical pulsators is the second
helium ionization zone, $\text{He}^+ \to \text{He}^{2+}$, near $4\times 10^4\ \text{K}$.

$$
% caption: The kappa-mechanism engine cycle; in the helium ionization zone opacity
% rises on compression, damming the outward flux and storing heat, which pushes the
% layer back out, whereupon opacity falls, heat escapes, and the layer falls back to
% compress again.
\begin{tikzpicture}[scale=1.0, font=\footnotesize]
\definecolor{acc}{HTML}{4A6FA5}
% four phases around a loop
\node[draw, very thick, align=center, minimum width=2.4cm, minimum height=1.0cm] (a) at (1.7,4.0) {compress};
\node[draw, very thick, align=center, minimum width=2.4cm, minimum height=1.0cm] (b) at (6.5,4.0) {opacity rises,\\heat stored};
\node[draw, very thick, align=center, minimum width=2.4cm, minimum height=1.0cm] (c) at (6.5,1.2) {push out,\\expand};
\node[draw, very thick, align=center, minimum width=2.4cm, minimum height=1.0cm] (d) at (1.7,1.2) {opacity falls,\\heat escapes};
\draw[black, very thick, ->] (a) -- (b);
\draw[black, very thick, ->] (b) -- (c);
\draw[black, very thick, ->] (c) -- (d);
\draw[black, very thick, ->] (d) -- (a);
\end{tikzpicture}
$$

## The instability strip

The kappa mechanism drives a star only when the helium ionization zone sits at the
right depth. Too hot, and the zone lies too near the surface, where it contains too
little mass to store significant heat, so the driving is feeble. Too cool, and the
outer envelope becomes convective; convection carries the energy and short-circuits
the valve, quenching the pulsation. Between these limits the zone lies at a depth
where its heat capacity is large and radiation still controls the flux, and the star
pulsates. The favorable range is nearly independent of luminosity, so it defines an
almost vertical band in the H-R diagram, the **instability strip**, tilted slightly
because the cool edge depends weakly on luminosity.

$$
% caption: The instability strip as a narrow, nearly vertical band crossing the H-R
% diagram; classical Cepheids occupy its luminous upper end, RR Lyrae stars its
% intersection with the horizontal branch, and Mira variables lie beyond on the
% asymptotic giant branch.
\begin{tikzpicture}[scale=1.0, font=\footnotesize]
\definecolor{acc}{HTML}{4A6FA5}
\draw[->, black] (9.0,0) -- (0,0) node[left, black!70] {hotter};
\draw[->, black] (9.0,0) -- (9.0,5.4) node[above, black!70] {luminosity};
% main sequence for reference
\draw[black, thick] (1.4,5.0) .. controls (4.0,3.0) and (6.5,1.4) .. (8.4,0.7);
\node[black, anchor=west, rotate=-27] at (5.4,1.9) {main sequence};
% instability strip band (tilted)
\fill[acc!10] (4.0,0.4) -- (4.9,0.4) -- (4.1,5.2) -- (3.2,5.2) -- cycle;
\draw[acc!60, densely dashed] (4.0,0.4) -- (3.2,5.2);
\draw[acc!60, densely dashed] (4.9,0.4) -- (4.1,5.2);
\node[acc, anchor=south, rotate=-80] at (3.55,2.6) {instability strip};
% Cepheids at top
\fill[black] (3.7,4.3) circle (2pt);
\node[black!70, anchor=west] at (3.85,4.4) {Cepheids};
% RR Lyrae mid
\fill[black] (4.3,2.2) circle (2pt);
\node[black!70, anchor=west] at (4.45,2.3) {RR Lyrae};
% Mira lower right (cool, luminous AGB - place upper right of strip)
\fill[black] (2.6,4.9) circle (2pt);
\node[black!70, anchor=west] at (2.7,4.9) {Mira};
\end{tikzpicture}
$$

The strip is populated by different stars at different luminosities:

- **Classical Cepheids** at the luminous top, massive stars ($3$–$20\,M_\odot$)
  crossing the strip on blue loops during core helium burning, with periods of one to
  a hundred days.
- **RR Lyrae stars** where the strip meets the horizontal branch, old low-mass stars
  of about $0.6\,M_\odot$ burning helium in the core, with periods near half a day.
- **Mira variables**, cool luminous stars on the asymptotic giant branch pulsating in
  the fundamental radial mode with periods of hundreds of days and large amplitudes.

## Light and velocity curves

A pulsating star is observed through its light curve and its radial-velocity curve.
The luminosity varies mainly with the surface temperature, which peaks slightly after
the star has passed through minimum radius and is expanding. The radial velocity,
measured from the Doppler shift of spectral lines, tracks the velocity of the surface
and is nearly the time derivative of the radius. Maximum brightness therefore lags the
epoch of maximum contraction, and the light and velocity curves are offset by roughly
a quarter of a period. For Cepheids the light curve is characteristically asymmetric,
rising fast and declining slowly.

$$
% caption: Light and radial-velocity curves of a Cepheid over one period; the
% brightness rises steeply and falls slowly, and the velocity curve, tracking the
% surface motion, is offset from the light curve by about a quarter period.
\begin{tikzpicture}[scale=1.0, font=\footnotesize]
\definecolor{acc}{HTML}{4A6FA5}
\draw[->, black] (0,2.6) -- (9.2,2.6) node[right, black!70] {phase};
\draw[->, black] (0,0.3) -- (0,5.0);
\node[acc, anchor=east] at (-0.05,4.2) {light};
\node[black, anchor=east] at (-0.05,1.3) {velocity};
% light curve: fast rise, slow decline (asymmetric), upper band
\draw[acc, very thick] (0.4,3.4) .. controls (0.8,4.6) and (1.1,4.6) .. (1.5,4.5)
  .. controls (3.0,4.1) and (4.3,3.3) .. (5.0,3.35)
  .. controls (5.6,3.4) and (6.0,4.55) .. (6.5,4.5)
  .. controls (8.0,4.1) and (8.6,3.4) .. (9.0,3.4);
% velocity curve: offset ~quarter period, lower band, dashed
\draw[black, very thick, densely dashed] (0.4,1.9) .. controls (1.4,1.0) and (2.2,1.0) .. (3.0,1.7)
  .. controls (3.8,2.3) and (4.4,2.3) .. (5.0,1.9)
  .. controls (5.9,1.0) and (6.7,1.0) .. (7.5,1.7)
  .. controls (8.2,2.2) and (8.7,2.2) .. (9.0,2.0);
\end{tikzpicture}
$$

## The period–luminosity relation

More luminous Cepheids are larger and less dense, so by the period–mean-density
relation they pulsate more slowly. This links period directly to luminosity: the
**period–luminosity relation**, or Leavitt law,

$$
M_V = a\,\log_{10}\!\left(\frac{\Pi}{\text{day}}\right) + b,
$$

with $a \approx -2.8$ and a small scatter set by the finite width of the instability
strip. Because the period is measured from the light curve alone, independent of
distance, the relation delivers the absolute magnitude, and comparison with the
apparent magnitude gives the distance modulus. Cepheids are luminous enough to be
resolved in galaxies tens of megaparsecs away, making them the rung that calibrates
the extragalactic scale in [the cosmic distance
ladder](/astrophysics-cosmology/observational-foundations/the-cosmic-distance-ladder).

$$
% caption: The Cepheid period–luminosity relation; the absolute magnitude brightens
% linearly with the logarithm of the pulsation period, with a narrow scatter set by
% the width of the instability strip.
\begin{tikzpicture}[scale=1.0, font=\footnotesize]
\definecolor{acc}{HTML}{4A6FA5}
\draw[->, black] (0,0) -- (9.0,0) node[right, black!70] {log period};
\draw[->, black] (0,0) -- (0,5.0) node[above, black!70] {brighter};
% P-L line (brighter = up)
\draw[acc, very thick] (0.7,0.7) -- (8.2,4.6);
\node[acc, anchor=south east] at (6.6,3.4) {slope near minus 2.8};
% scatter points about it
\fill[black] (1.4,1.2) circle (1.4pt);
\fill[black] (2.3,1.5) circle (1.4pt);
\fill[black] (3.1,2.1) circle (1.4pt);
\fill[black] (4.0,2.3) circle (1.4pt);
\fill[black] (5.0,3.1) circle (1.4pt);
\fill[black] (5.9,3.2) circle (1.4pt);
\fill[black] (6.8,3.9) circle (1.4pt);
\fill[black] (7.5,4.1) circle (1.4pt);
\end{tikzpicture}
$$

Pulsation converts a star's structure into an observable clock. The period fixes the
mean density and, through the instability-strip physics, the luminosity, and the
brightness variation is visible across intergalactic distances. The same instability
strip is crossed by massive stars on blue loops, treated in [the evolution of massive
stars](/astrophysics-cosmology/stellar-evolution/the-evolution-of-massive-stars), and
by horizontal-branch stars after the helium flash, treated in [post-main-sequence
evolution of low-mass
stars](/astrophysics-cosmology/stellar-evolution/post-main-sequence-low-mass-evolution).

[^co-puls]: Carroll & Ostlie, §14.1–14.2 — Stellar Pulsation: the period–mean-density relation and pulsation constant, radial modes, the kappa and gamma mechanisms seated in the helium ionization zone, the instability strip, and the Cepheid period–luminosity relation.