Stellar pulsations and asteroseismology
Connect small stellar oscillations to sound travel, buoyancy, and structure while separating stable frequencies from excitation and damping.
On this page
The question for today
What you’ll understand
- Distinguish acoustic and buoyancy restoring forces in stellar normal modes.
- Compute an asymptotic acoustic frequency separation from a supplied travel time.
- State the homology and linearity assumptions behind density-frequency scaling.
A pulsating star offers information about its interior because waves sample pressure, density, sound speed, and stratification below the photosphere. The observed spectrum is a boundary-value problem, not merely a single oscillator attached to the surface. We begin with small perturbations of a spherical, slowly evolving hydrostatic model, neglecting rapid rotation and strong magnetic effects. Linear adiabatic theory supplies useful frequencies, while excitation, damping and nonlinear amplitude selection require additional physics. Our examples are hypothetical asymptotic spectra with supplied travel times and phase offsets, rather than fits to measured stars.
Identify the restoring force and perturbation domain
Pressure modes are predominantly acoustic, with pressure gradients providing the main restoring force. Gravity modes are predominantly restored by buoyancy in stably stratified regions; the word gravity here does not mean a gravitational wave of spacetime. Some evolved stars support mixed modes that sample both an acoustic envelope cavity and a buoyancy-dominated core cavity. Mode labels specify radial order and angular structure, but their precise sign and counting conventions vary, so a numerical comparison must declare how labels are assigned.
For adiabatic perturbations the sound speed squared is the adiabatic pressure derivative at fixed composition. It is not generally the isothermal derivative. Radial dynamical stability also depends on the pressure response during compression. A simplified homologous-compression argument gives a four-thirds threshold for a constant adiabatic exponent in a Newtonian model. A small outer region below that value does not alone prove global collapse; actual stability depends on the entire structure and mode equations. Relativistic corrections and nonadiabatic effects add further conditions.
Acoustic travel time organizes high-order modes
For high radial order and low angular degree in an appropriate acoustic cavity, the radial-order spacing is approximately the inverse of the round-trip sound travel time between center and surface. The asymptotic frequency formula includes a phase offset from turning regions and the surface. This pattern makes the travel-time integral diagnostically useful, but its assumptions can fail for low-order modes, mixed modes, sharp structural features, or rapidly varying surface phase. A constant phase offset is an explicit toy simplification.
The large separation scales roughly with the square root of mean density when comparing homologous structures with the same dimensionless sound-speed profile. It is not an exact relation across all stellar masses and stages. Deviations from a regular spacing, rotational splitting, and gravity-mode period patterns can provide additional information, but interpreting them requires mode identification and an appropriate structural model. In particular, an acoustic spacing alone does not uniquely reveal chemical composition or distinguish all possible core-burning histories.
Frequency, growth, and observed amplitude are different outputs
Writing perturbations with a minus-i-omega time convention makes a positive imaginary part of omega correspond to growth. In a damped mode the imaginary part is negative. Adiabatic stable-mode calculations give real frequencies but do not predict observed amplitudes. Nonadiabatic work from opacity variations, nuclear heating, radiative leakage, or convective heat exchange can drive or damp oscillations. Turbulent stochastic forcing can maintain a spectrum of otherwise damped modes. Strong amplitudes eventually require nonlinear dynamics and can invalidate linear eigenfunctions.
Observation adds geometric visibility and finite measurement duration: different angular patterns can cancel when the unresolved stellar disk is integrated. A broadened spectral peak is therefore interpreted through a model of excitation, damping, observation window, and noise rather than by identifying its width with a structural frequency spacing. The equal-frequency wave instrument can test coherent addition and cancellation, but it has no eigenfrequency search, spherical harmonics, internal cavities, or mode excitation. Its cancellation is an exact algebraic analogy, not evidence that a real star contains no oscillation energy.
Work it through
Read an ideal acoustic comb
A hypothetical acoustic model has one-way travel time τac=1000 s. Use the asymptotic formula with ε=1/2 and ℓ=0 to estimate Δν and the n=12 and n=13 frequencies. Treat the phase offset as fixed and use 1 μHz=10⁻⁶ Hz.
The round-trip time is 2000 s, so . The factor two is essential.
For n=12 the dimensionless factor is 12.5, giving .
For n=13 the factor is 13.5, giving 6750 μHz. Their difference is 500 μHz, consistently reproducing the supplied asymptotic spacing.
The toy large separation is 500 μHz, with consecutive estimated radial frequencies 6.25 and 6.75 mHz; no claim of exact low-order stellar eigenfrequencies is implied.
The living notebook
Explore the model
Interactive lab · waves
Coherent visibility can cancel a signal
Reset and set BOTH Phase difference to 0° and Amplitude ratio to 1; predict a sum twice either component. Then Reset and set BOTH Phase difference to 180° and Amplitude ratio to 1; predict a zero sum at every position and time. Play or Step to verify both. These are equal-frequency co-propagating sine waves, illustrating coherent cancellation only; the instrument has no stellar mode spectrum, standing-wave cavity, angular visibility integral, or damping calculation.
- Sum at x = 0
- 0.866
- Resultant amplitude
- 1.732
- Phase difference
- 60°
Two waves share the same frequency and wave number. Playback advances t in sin(2x − t) + a sin(2x − t + phase). Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.
Look forCancellation of a measured combination does not by itself imply the underlying components carry no oscillation.
Pause & reason
Check your understanding
Your turn
Try it yourself
Compare two hypothetical homologous acoustic models with the same dimensionless structure. Model 2 has twice the mass and twice the radius of model 1. Find the large-separation ratio and one-way acoustic-travel-time ratio. Explain the role of the homology assumption.
Take a moment to make a prediction or work it out on paper. Open a hint when you need one.
Give me a hint
Apply the square-root mean-density scaling before taking the inverse relation for travel time.
Compare with the solution
The mean-density ratio is , so the separation ratio is 1/2. Since Δν is approximately the inverse of twice travel time, the travel-time ratio is two. This follows because the dimensionless sound-speed profile is held fixed. Different core composition, stratification or surface phases can spoil the simple scaling, so the result is not an exact radius estimator for an arbitrary pair of stars.
Keep the idea
Take it with you
- Sound and buoyancy provide different restoring forces and can couple in mixed modes.
- The acoustic large separation measures a travel-time integral under asymptotic assumptions.
- Adiabatic frequencies, damping or growth rates, and observable amplitudes require distinct calculations.
Read further
Original sources and deeper explanations to follow at your own pace.
Optional progress, saved in this browser when available. Revisit any lesson at any time.
Lesson not marked complete.