advanced35 min studyLesson 899 in Physics

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

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.

cs2=(∂P∂ρ)s,Xi=Γ1Pρ,Γ1=(∂ln⁡P∂ln⁡ρ)s,Xic_s^2=\left(\frac{\partial P}{\partial\rho}\right)_{s,X_i}=\Gamma_1\frac{P}{\rho},\qquad \Gamma_1=\left(\frac{\partial\ln P}{\partial\ln\rho}\right)_{s,X_i}
cs is m s⁻¹ and Γ₁ is dimensionless. These derivatives concern a displaced material element. The background polytropic exponent need not equal Γ₁. A local sound speed and buoyancy frequency enter a global eigenvalue problem with center and surface conditions.↔ Scroll wide equations horizontally.

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.

τac=∫0Rdrcs,Δν≃12τac,νnℓ≃Δν(n+ℓ2+ε),Δν∝MR3 (homology)\tau_{\rm ac}=\int_0^R\frac{dr}{c_s},\qquad \Delta\nu\simeq\frac1{2\tau_{\rm ac}},\qquad \nu_{n\ell}\simeq\Delta\nu\left(n+\frac{\ell}{2}+\varepsilon\right),\qquad \Delta\nu\propto\sqrt{\frac{M}{R^3}}\ \text{(homology)}
τac is one-way acoustic travel time in s; Δν and ν are cyclic frequencies in Hz, not angular frequencies. n is a chosen high radial-order integer and ℓ angular degree. ε absorbs the phase convention. The mean-density scaling assumes homologous dimensionless structure and is approximate outside that family.↔ Scroll wide equations horizontally.

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.

ξ(r,t)=Re⁡ ⁣[ξ^(r)e−iωt],ω=ωr+iγ,∣ξ∣envelope∝eγt\xi(\mathbf r,t)=\operatorname{Re}\!\left[\hat\xi(\mathbf r)e^{-i\omega t}\right],\qquad \omega=\omega_r+i\gamma,\qquad |\xi|_{\rm envelope}\propto e^{\gamma t}
ωr and γ have units s⁻¹. This convention gives γ>0 growth and γ<0 damping; the real cyclic frequency is ωr/(2π). The envelope statement describes linear exponential evolution of a mode amplitude, not nonlinear saturation or an inferred observational probability.↔ Scroll wide equations horizontally.

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.

  1. The round-trip time is 2000 s, so Δν=1/2000 Hz=500 μHz\Delta\nu=1/2000\ {\rm Hz}=500\ \mu{\rm Hz}. The factor two is essential.

  2. For n=12 the dimensionless factor is 12.5, giving ν12,0=6250 μHz=6.25 mHz\nu_{12,0}=6250\ \mu{\rm Hz}=6.25\ {\rm mHz}.

  3. For n=13 the factor is 13.5, giving 6750 μHz. Their difference is 500 μHz, consistently reproducing the supplied asymptotic spacing.

The result

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.

Explore the model

Interactive lab

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.

Coherent visibility can cancel a signalTwo waves share the same frequency and wave number. Playback advances t in sin(2x − t) + a sin(2x − t + phase). Current readings: Sum at x = 0: 0.866; Resultant amplitude: 1.732; Phase difference: 60°. All parameters are adjustable with the labeled controls below.−3.14−1.5701.573.14−3.4−1.701.73.4x · radiansAmplitude
SumFirst componentSecond component
Sum at x = 0
0.866
Resultant amplitude
1.732
Phase difference
60°
t = 0 s

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.

Check your understanding

What sets the approximate large separation of high-order, low-degree acoustic modes?

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 2/23=1/42/2^3=1/4, 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.

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.