Rotational Mixing (Astrophysics)

Rotational mixing is the collective set of fluid-dynamical instabilities and meridional circulation currents that transport chemical species and angular momentum across stably stratified radiative zones in rotating stars. In a non-rotating, non-magnetic star, radiative envelopes satisfy the Ledoux criterion and remain chemically unmixed outside zones of convective overshooting or double-diffusive instability. When a star rotates, however, centrifugal distortion destroys spherical symmetry—producing radiative flux divergence on oblate equipotential surfaces as demonstrated by the von Zeipel theorem—and drives large-scale Eddington-Sweet circulation alongside vertical differential rotation dΩ / dr that triggers dynamical and secular shear instabilities. Within astronomy, rotational mixing accounts for the observed surface enrichment of nitrogen (¹⁴N) and helium accompanied by carbon and boron depletion in main-sequence O and B stars, lowers the mass threshold for Wolf-Rayet star formation at low metallicity, and drives chemically homogeneous evolution (CHE) in rapidly spinning high-mass stars that produce long gamma-ray bursts and merging binary black holes.

Anisotropic Horizontal vs. Vertical Turbulence

Because the restoring Brunt–Väisälä buoyancy force (N² > 0) acts parallel to the effective gravity vector g_eff (perpendicular to isobaric surfaces P = constant) while no buoyant restoring force opposes fluid motions along horizontal isobaric surfaces, turbulence in a rotating radiative zone is strongly anisotropic. As formalized by Jean-Paul Zahn (1992) and developed by André Maeder and Georges Meynet (1998–2012), the horizontal turbulent eddy diffusivity D_h is orders of magnitude larger than the vertical turbulent diffusivity D_v (D_h ≫ D_v).

This vigorous horizontal turbulence rapidly erases latitudinal variations in angular velocity Ω(r, θ) and chemical composition along each isobaric shell, enforcing shellular rotation in which angular velocity depends almost exclusively on the mean radius r of the isobar:

Ω(r, θ) ≈ Ω(r)

Coupled Advection–Diffusion Transport of Angular Momentum and Chemicals

Under shellular rotation, the radial transport of specific angular momentum satisfies an advection–diffusion partial differential equation (as distinct from pure diffusion), because Eddington-Sweet circulation with radial velocity profile U₂(r) P₂(cos θ) advects angular momentum across shells:

ρ (d / dt)[ r² Ω ]_Mr = (1 / 5 r²) (∂ / ∂r)[ ρ r⁴ Ω U₂(r) ] + (1 / r²) (∂ / ∂r)[ ρ D_shear r⁴ (∂Ω / ∂r) ]

By contrast, as shown by Chaboyer and Zahn (1992), the interaction between meridional advection U₂(r) and strong horizontal turbulence D_h reduces the net radial transport of chemical species (mass fraction X_i) to an effective diffusion coefficient D_eff added to the vertical shear diffusivity D_shear:

D_rot = D_eff + D_shear, where D_eff = | r U₂(r) |² / (30 D_h)

Because D_h appears in the denominator of D_eff, stronger horizontal turbulence actually homogenizes chemical composition along horizontal shells before meridional currents can carry it far vertically, whereas vertical shear turbulence D_shear directly mixes composition across isobars.

The Richardson Criterion and Secular Shear Mixing

Differential rotation (S ≡ r |dΩ / dr| ≠ 0) stores free rotational kinetic energy in adjacent spherical shells moving at different velocities. Whether Kelvin–Helmholtz shear eddies can overturn the stably stratified radiative zone depends on the dimensionless Richardson number Ri, which compares the potential energy required to lift a parcel against the Brunt–Väisälä buoyancy frequency N² = N_T² + N_μ² to the kinetic energy available in the vertical velocity shear:

Ri ≡ N² / (r dΩ / dr)² < Ri_crit ≈ 1/4

1. Dynamical Shear Instability: If the shear is so steep that the unfiltered Richardson criterion Ri < 1/4 is satisfied directly, Kelvin–Helmholtz billows grow on the dynamical timescale.
2. Secular (Thermal) Shear Instability: In most stellar radiative zones, the thermal buoyancy N_T² is large enough that Ri > 1/4 for adiabatic displacements. However, just as in semiconvection and stellar thermohaline mixing, radiative thermal diffusion K_T weakens the thermal restoring force N_T² on small eddy scales without erasing the composition barrier N_μ². Accounting for radiative heat leakage (Zahn, 1992; Maeder, 1997; Talon & Zahn, 1997) yields the secular shear diffusivity:

D_shear ≈ f_energ [ H_P / (g δ) ] [ K_T r² (dΩ / dr)² ] / [ ∇_ad − ∇ + (K_T / D_h) (φ / δ) ∇_μ ]

where f_energ ≈ Ri_crit = 1/4. Additional baroclinic instabilities—including the Goldreich–Schubert–Fricke (GSF) instability (which develops when ∂Ω / ∂z ≠ 0 parallel to the rotation axis) and the ABCD instability—plus magnetic torques from the Tayler–Spruit dynamo further couple the core and envelope.

Surface Nitrogen, Carbon, and Boron Abundances in OB Stars

In main-sequence stars with masses M ≳ 2 M_☉, core hydrogen fusion proceeds via the CNO cycle, which converts most core carbon (¹²C) and oxygen (¹⁶O) into nitrogen-14 (¹⁴N) because the proton-capture reaction ¹⁴N(p,γ)¹⁵O has the smallest nuclear cross section at the Gamow peak. Simultaneously, fragile light elements such as boron (¹⁰B, ¹¹B) are destroyed by proton capture at shallow depths where T ≳ 5 × 10⁶ K.

High-resolution spectroscopy from the VLT-FLAMES Survey of Massive Stars across the Milky Way and Magellanic Clouds confirms that rapidly rotating main-sequence B and O stars (v sin i ≳ 150–300 km s⁻¹) exhibit surface ¹⁴N/H enrichments of 0.3 to 0.8 dex correlated with strong surface boron depletion and modest ¹²C depletion long before the main-sequence turnoff, validating rotational diffusion across the radiative envelope.

Low-Metallicity Evolution and Chemically Homogeneous Evolution (CHE)

Rotational mixing becomes markedly more efficient at low metallicity (Z ≪ Z_☉) for two reasons: (1) metal-poor stars have lower line-driven radiative opacities and are more compact (smaller radius R and steeper internal gradients for a given angular momentum), and (2) weaker line-driven stellar winds preserve high initial spin rates throughout the main sequence:

1. Higher Luminosity and Wolf-Rayet Populations: Rotational mixing feeds fresh hydrogen into the convective core and helium into the radiative envelope, increasing main-sequence lifetimes by 20–30% and lowering the minimum initial mass required to form stripped-envelope Wolf-Rayet stars from ~35–40 M_☉ down to ~22–25 M_☉.
2. Chemically Homogeneous Evolution (CHE): When a high-mass star (M ≳ 20–30 M_☉) rotates above ~40–50% of its breakup velocity (either from birth at low metallicity or via tidal synchronization in a tight binary), rotational mixing becomes faster than the nuclear fusion timescale (τ_mix < τ_nuc). No composition barrier ∇_μ ever forms at the core boundary; the entire star converts homogeneously into a compact helium Wolf-Rayet star that evolves blueward on the Hertzsprung–Russell diagram without ever expanding into a red supergiant. This CHE channel avoids common envelope evolution and forms heavy binary black holes (30–50 M_☉) that merge within a Hubble time, as detected by LIGO/Virgo/KAGRA, as well as collapsar progenitors of long-duration gamma-ray bursts.

Why is boron depletion used as a clean test of rotational mixing in B-type stars?

Unlike nitrogen—which is synthesized deep inside the CNO-burning core at T > 2 × 10⁷ K and must be transported across the entire radiative zone—boron (¹⁰B and ¹¹B) is destroyed by proton capture in the outer radiative layers at T ≈ 5–6 × 10⁶ K. Consequently, stars experiencing shallow rotational mixing show surface boron depletion early on the main sequence before core-synthesized ¹⁴N reaches the photosphere, mapping the radial depth of the mixing profile.

How do internal magnetic fields modify rotational mixing?

When a magnetic dynamo or fossil field exerts strong Maxwell stresses across the radiative zone, it suppresses radial differential rotation (dΩ / dr → 0). While this reduces Kelvin–Helmholtz shear mixing (D_shear), solid-body rotation strengthens the centrifugal perturbation in the outer envelope and accelerates Eddington-Sweet circulation, producing distinct asteroseismic gravity-mode splitting patterns.

1. Zahn, J.-P. (1992). "Circulation and turbulence in rotating stars." Astronomy and Astrophysics, 265, 115–132. https://ui.adsabs.harvard.edu/abs/1992A%26A...265..115Z
2. Maeder, A., & Meynet, G. (2000). "The Evolution of Rotating Stars." Annual Review of Astronomy and Astrophysics, 38, 143–190. https://doi.org/10.1146/annurev.astro.38.1.143
3. Hunter, I., Brott, I., Langer, N., Lennon, D. J., Dufton, P. L., Howarth, I. D., ... & Evans, C. J. (2009). "The VLT-FLAMES Survey of Massive Stars: constraints on stellar evolution from the chemical composition of rapidly rotating Galactic and Magellanic Cloud B-type stars." Astronomy & Astrophysics, 496(3), 841–853. https://doi.org/10.1051/0004-6361/200810024
4. Mandel, I., & de Mink, S. E. (2016). "Merging binary black holes formed through chemically homogeneous evolution in short-period stellar binaries." Monthly Notices of the Royal Astronomical Society, 458(3), 2634–2647. https://doi.org/10.1093/mnras/stw379