Gravity Mode (Astrophysics)

In stellar astrophysics and asteroseismology, a gravity mode (abbreviated g-mode) is a standing non-radial oscillation wave inside a star or planet for which buoyancy (the gravitational restoring force acting on density perturbations across a stably stratified medium, quantified by the Brunt–Väisälä frequency N) serves as the primary restoring force. First classified alongside acoustic pressure modes (p-modes) and surface fundamental modes (f-modes) by Thomas G. Cowling in 1941, astrophysical g-modes are internal fluid buoyancy waves governed by classical Newtonian gravity and are physically distinct from relativistic gravitational waves (ripples in the spacetime metric governed by Einstein's general relativity). Because buoyancy acts as a restoring force only where the square of the Brunt–Väisälä frequency is positive (N² > 0, satisfying the Ledoux criterion), g-modes propagate freely inside radiative zones and are evanescent (undergoing exponential decay) inside convective zones (N² < 0). Within astronomy, space-based asteroseismology of g-modes and mixed p–g modes with CoRoT, Kepler, and TESS provides the primary observational probe of convective overshooting, core helium burning, chemical composition gradients, and core-to-envelope differential rotation across the Hertzsprung–Russell diagram.

Spherical Harmonic Geometry (n, ℓ, m)

In a spherically symmetric star, the displacement vector ξ(r, θ, φ, t) of a linear normal oscillation mode is separable into a radial eigenfunction and a spherical harmonic Y_ℓ^m(θ, φ) oscillating with angular frequency ω = 2π ν:

ξ(r, θ, φ, t) = [ ξ_r(r) e_r + ξ_h(r) ( e_θ ∂/∂θ + [e_φ / sin θ] ∂/∂φ ) ] Y_ℓ^m(θ, φ) exp(−i ω t)

Each normal mode is indexed by three integer quantum numbers:

1. Radial order (n): The number of concentric radial nodes in the radial displacement eigenfunction ξ_r(r) (conventionally assigned negative integers n = −1, −2, −3, … for g-modes and positive integers n = +1, +2, … for p-modes).
2. Angular degree (ℓ ≥ 1): The total number of surface nodal lines on the sphere (note that radial modes with ℓ = 0 cannot exist as g-modes because buoyancy requires horizontal fluid displacement along an equipotential surface to create a lateral density contrast; hence all g-modes are non-radial with ℓ ≥ 1).
3. Azimuthal order (m ∈ [−ℓ, +ℓ]): The number of surface nodal lines passing through the stellar rotation axis.

The Characteristic Frequencies: Brunt–Väisälä (N) and Lamb (S_ℓ)

In the Cowling approximation (neglecting perturbations to the gravitational potential Φ') and the local WKB short-wavelength limit, the radial wavenumber k_r(r) of a stellar oscillation of angular frequency ω satisfies the dispersion relation:

k_r²(r) = ( ω² / c_s² ) ( 1 − S_ℓ² / ω² ) ( 1 − N² / ω² )

where c_s = (Γ₁ P / ρ)^1/2 is the adiabatic sound speed and the two characteristic cutoff frequencies of the stellar interior are:

• The Lamb (Acoustic) Frequency (S_ℓ): Defined as S_ℓ²(r) ≡ ℓ(ℓ + 1) c_s²(r) / r² = k_h² c_s², which represents the inverse time required for an acoustic sound wave to travel horizontally over one angular wavelength.
• The Brunt–Väisälä (Buoyancy) Frequency (N): Defined via the Ledoux criterion as N²(r) = (g / H_P) [ δ (∇_ad − ∇) + φ ∇_μ ], which represents the natural vertical oscillation frequency of a displaced fluid parcel under buoyancy.

For radial wave propagation to occur, k_r²(r) must be positive (k_r² > 0), which occurs in two distinct frequency domains:

Because g-modes require ω² < N², they cannot propagate inside convective zones (where N² < 0). In solar-type stars like the Sun, g-modes are trapped inside the radiative core (r < 0.713 R_☉) and must tunnel through the thick convective envelope to reach the photosphere, attenuating their surface amplitudes to sub-millimeter-per-second levels. Conversely, in stars with convective cores and radiative envelopes (such as γ Doradus and Slowly Pulsating B stars) or radiative sub-photospheric layers (such as pulsating white dwarfs and subdwarf B stars), g-modes reach the surface with readily observable photometric amplitudes.

Uniform Period Spacing (ΔΠ_ℓ)

Whereas high-order acoustic p-modes are asymptotically equally spaced in frequency (Δν = ν_n+1,ℓ − ν_n,ℓ ≈ [2 ∫ dr / c_s]⁻¹), Monique Tassoul proved in 1980 that high-order gravity modes (|n| ≫ ℓ) in a chemically homogeneous radiative cavity bounded between inner radius r₁ and outer radius r₂ are asymptotically equally spaced in period Π_n,ℓ = 2π / ω_n,ℓ:

Π_n,ℓ ≈ Π₀ ( |n| + ℓ/2 + α_g ) / √[ ℓ(ℓ + 1) ]

Consequently, the asymptotic period spacing ΔΠ_ℓ ≡ Π_n+1,ℓ − Π_n,ℓ between consecutive radial orders is constant for a given angular degree ℓ:

ΔΠ_ℓ = Π₀ / √[ ℓ(ℓ + 1) ], where Π₀ = 2π² / ∫_r1^r₂ ( N(r) / r ) dr

Measuring Π₀ directly measures the buoyancy integral ∫ (N / r) dr at the edge of the convective core, calibrating the mass of the convective core and the extent of convective overshooting (f_ov).

Chemical Gradient Mode Trapping and Rotational Tilt

1. Buoyancy Glitches and Mode Trapping: As a star burns hydrogen on the main sequence, its receding convective core leaves behind a steep mean molecular weight gradient ∇_μ > 0, creating a sharp spike in N_μ² at radius r_μ. Modes whose radial nodes align with r_μ are partially trapped in the μ-gradient zone, producing periodic oscillatory dips around the mean ΔΠ_ℓ baseline whose period and amplitude measure the Evolutionary age (X_c) and diffusive mixing profile (D_ov, rotational mixing) of the star.
2. Rotational Splitting and Core Spin Rates: In a rotating star with internal angular velocity Ω(r), the Coriolis force lifts the azimuthal (m) degeneracy and modifies g-mode periods according to the Traditional Approximation of Rotation (TAR). In a plot of period spacing ΔΠ_ℓ,m versus period Π, prograde modes (m > 0) exhibit a downward linear slope while retrograde modes (m < 0) exhibit an upward slope directly proportional to the near-core rotation frequency Ω_core.

• Slowly Pulsating B (SPB) Stars (M ≈ 3–9 M_☉) and γ Doradus Stars (M ≈ 1.3–1.8 M_☉): Main-sequence stars with convective cores and radiative envelopes whose high-order dipole (ℓ = 1) and quadrupole (ℓ = 2) g-modes (periods 0.4 to 3 days), excited by the iron-peak κ-mechanism or convective flux blocking, have yielded precise measurements of near-core rotation and convective overshooting from Kepler and TESS.
• Mixed p–g Modes in Red Giants and Red Clump** Stars:** In evolved red giants, the compact radiative helium core has such a high Brunt–Väisälä frequency N that core g-mode frequencies overlap the envelope p-mode frequencies. Quantum-like wave tunneling across the narrow evanescent zone creates mixed p–g modes that simultaneously reveal the core rotation rate (showing that red-giant cores rotate only 5–15 times faster than their envelopes, requiring strong angular momentum transport beyond pure hydrodynamic rotational mixing) and cleanly distinguish hydrogen-shell-burning red-giant-branch stars (ΔΠ₁ ≈ 50–100 s) from core-helium-burning red clump stars (ΔΠ₁ ≈ 250–350 s).
• Pulsating White Dwarfs** (ZZ Ceti / DAV, DBV, GW Vir) and Subdwarf B (sdB) Stars:** Compact remnants oscillating in non-radial g-modes with periods of 100 to 1,500 seconds, where mode trapping at the C/O core, He mantle, and H atmosphere boundaries measures layered shell masses and white-dwarf crystallization rates.

What is the difference between astrophysical gravity modes (g-modes) and relativistic gravitational waves?

Gravity modes (g-modes) are fluid-mechanical buoyancy waves propagating inside a star's plasma (analogous to internal ocean waves), where classical Newtonian buoyancy acts as the restoring force and the waves cannot travel through empty space. Gravitational waves are relativistic transverse ripples in the fabric of spacetime itself predicted by Einstein's general relativity and detected by LIGO/Virgo from merging black holes and neutron stars.

How do mixed modes distinguish a red-giant-branch star from a red-clump star?

A red-giant-branch star has a radiative, electron-degenerate helium core where g-modes propagate all the way to r → 0, making the integral ∫ (N / r) dr very large and the dipole period spacing small (ΔΠ₁ ≈ 50–100 s). Once helium ignites in the helium flash to form a red clump star, a central convective core (N² < 0) develops and excludes g-modes from the inner core, reducing ∫ (N / r) dr and jumping ΔΠ₁ up to 250–350 s.

1. Cowling, T. G. (1941). "The non-radial oscillations of polytropic stars." Monthly Notices of the Royal Astronomical Society, 101(8), 367–375. https://doi.org/10.1093/mnras/101.8.367
2. Tassoul, M. (1980). "Asymptotic approximations for stellar nonradial pulsations." The Astrophysical Journal Supplement Series, 43, 469–490. https://doi.org/10.1086/190678
3. Bedding, T. R., Mosser, B., Huber, D., Montalbán, J., Beck, P., Christensen-Dalsgaard, J., ... & White, T. R. (2011). "Gravity modes as a way to distinguish between hydrogen- and helium-burning red giant stars." Nature, 471(7340), 608–611. https://doi.org/10.1038/nature09935
4. Aerts, C. (2021). "Probing the interior physics of stars through asteroseismology." Reviews of Modern Physics, 93(1), 015001. https://doi.org/10.1103/RevModPhys.93.015001