Magnetic Braking (Astronomy)
Magnetic braking is the magnetohydrodynamic process by which a rotating star with a convective envelope loses spin or binary orbital angular momentum through an ionized stellar wind forced to corotate along open coronal magnetic field lines out to the Alfvén radius.
Magnetic braking is the magnetohydrodynamic angular momentum loss mechanism by which rotating stars—and tidally synchronized close binary systems—shed angular momentum via an ionized stellar wind threaded by a large-scale coronal magnetic field. First proposed by the French astrophysicist Evry Schatzman in 1962 to explain why low-mass main-sequence stars cooler than spectral type F5 (M ≲ 1.3 M☉) rotate orders of magnitude more slowly than hotter O, B, and A stars, and formalized analytically in the equatorial magnetohydrodynamic solar wind solution of Edmund J. Weber and Leverett Davis Jr. in 1967, magnetic braking relies on the magnetic tension of open field lines to force escaping wind plasma ([Parker wind]) into quasi-rigid corotation out to a large Alfvén radius RA ≫ R*. Because specific angular momentum scales with the square of the lever arm (j = Ω RA2), even a minuscule mass-loss rate (such as the present solar wind, dM / dt ≈ 2 × 10−14 M☉ yr−1) extracts substantial angular momentum over gigayear timescales. Within astronomy, magnetic braking governs the Skumanich law (veq ∝ t−1/2) and stellar age dating via gyrochronology, shapes the pre-main-sequence spin evolution of T Tauri stars, and drives the secular orbital inspiral of cataclysmic variables and low-mass X-ray binaries above the 2–3 hour orbital period gap.
The Kraft Break and Convective Dynamo Threshold
Why Hot Stars Spin Fast and Cool Stars Spin Slowly
Spectroscopic surveys of projected rotational velocities (v sin i) along the main sequence reveal a sharp transition discovered by Robert P. Kraft in 1967—the Kraft break—near effective temperature Teff ≈ 6,250 K (spectral type F5V, stellar mass M ≈ 1.3 M☉):
- Above the Kraft Break (M ≳ 1.3 M☉, O, B, A, and early-F stars): Stars possess convective cores surrounded by stably stratified radiative envelopes. Lacking a deep sub-photospheric convective zone to power a solar-type αΩ magnetic dynamo or a hot, magnetically channeled corona, these stars retain most of their natal angular momentum from Kelvin-Helmholtz contraction and rotate rapidly (veq ≈ 100 to 300 km s−1).
- Below the Kraft Break (M ≲ 1.3 M☉, late-F, G, K, and M dwarfs): Hydrogen ionization opacity near the photosphere creates a sub-surface convective envelope whose thickness increases toward lower stellar mass (becoming fully convective below M ≈ 0.35 M☉). Cyclonic convective turbulence combined with differential rotation generates a self-sustaining magnetic dynamo (B ~ 1 to 103 G) and heats a ~106 K X-ray corona that launches a continuous ionized Parker wind, spinning the star down via magnetic braking to veq ≈ 1 to 15 km s−1 within a few hundred million years.
Weber–Davis Theory and the Alfvén Lever Arm
Corotation Out to the Alfvén Surface
In a non-magnetic star of radius R* rotating at angular velocity Ω, an escaping stellar wind detaches directly from the photosphere, carrying away a small specific angular momentum j0 = (2/3) Ω R*2. To reduce the Sun's angular momentum by a factor of 10 via an unmagnetized wind would require shedding more than 90% of the Sun's mass.
When the star possesses an open coronal magnetic field B, however, the conductive ionized wind plasma is frozen onto the magnetic field lines by Alfvén's theorem. Close to the star (r < RA), the magnetic energy density B2 / (8π) exceeds the kinetic energy density (1/2) ρ vr2 of the outflowing wind. Magnetic tension forces the outflowing plasma to corotate almost rigidly at angular velocity Ω with the stellar surface, while magnetic stresses transfer angular momentum outward through the sub-Alfvénic wind.
As the wind accelerates outward according to the Parker wind equation and the magnetic field strength B(r) decays with radius, the radial wind velocity vr(r) eventually reaches the local radial Alfvén speed vA(r):
vr(RA) = vA(RA) ≡ Br(RA) / √(4π ρ(RA))
The spherical (or axisymmetric) locus where vr = vA is the Alfvén radius RA. Beyond RA, the wind kinetic energy dominates over magnetic tension, the field lines are wound into an Archimedean Parker spiral, and the wind coasts outward conserving its specific angular momentum.
Weber–Davis Angular Momentum Loss Formula
As proved by Edmund J. Weber and Leverett Davis Jr. (1967) from the conservation of mass flux and magnetohydrodynamic azimuthal momentum (including both plasma azimuthal motion and Maxwell magnetic stress Br Bφ / 4π), the total torque dJ / dt exerted on the star equals the mass-loss rate dM / dt multiplied by the specific angular momentum of rigid corotation at the Alfvén radius:
dJ / dt = − (2/3) ( dM / dt ) Ω RA2 = ( dJ / dt )unmag × ( RA / R* )2
In the present-day Sun, the mean Alfvén radius measured directly by NASA's Parker Solar Probe when it crossed into the sub-Alfvénic solar corona in 2021 lies at RA ≈ 15 to 20 R☉ (0.07–0.09 au). Consequently, the magnetic lever arm amplifies the angular momentum removed per gram of solar wind by a factor of (RA / R☉)2 ≈ 225 to 400! In young, rapidly rotating T Tauri stars and zero-age main-sequence G/K dwarfs with kilogauss fields, RA reaches 30 to 50 R*, amplifying the braking torque by >103.
Kawaler's Wind Braking Law and Dynamo Saturation
Dependence on Field Geometry (n) and Dynamo Scaling (a)
In 1988, Steven D. Kawaler generalized the Mestel (1968) and Weber–Davis torque into a unified parametric formulation for stellar evolution codes. Assuming a radial magnetic field variation B(r) = B0 (R* / r)n—where n = 2 corresponds to a pure radial split-monopole field (B ∝ r−2) and n = 3 corresponds to a closed dipole field (B ∝ r−3)—and a dynamo law scaling surface field strength with rotation rate as B0 ∝ Ωa, Kawaler derived:
dJ / dt = − KW Ω1 + (4a n / 3) R*2 − n ( dM / dt )1 − (2n / 3) M*−n / 3
Setting a linear dynamo law (a = 1) and an intermediate coronal field geometry n = 3/2 (reflecting a mixture of closed dipolar loops near the equator and open radial streamers at higher latitudes, as confirmed by Zeeman–Doppler imaging and 3D MHD wind simulations by Sean Matt et al., 2012) yields dJ / dt ∝ −Ω3. Integrating I (dΩ / dt) = −K Ω3 at constant moment of inertia I recovers the classical Skumanich law:
Ω(t) ∝ t−1/2 (or rotation period Prot(t) ∝ t+1/2)
| Braking Regime | Rossby Number (Ro = Prot / τc) | Torque Scaling (dJ / dt) | Period Evolution (Prot(t)) | Physical Mechanism |
|---|---|---|---|---|
| Saturated Dynamo Regime (Ultrafast Rotators) | Ro < Rosat ≈ 0.13 | dJ / dt ∝ −Ω | Exponential spin-down (Prot ∝ exp(t / τ)) | Coronal magnetic flux and X-ray emission saturate (LX / Lbol ≈ 10−3); independent of Ω |
| Unsaturated Skumanich Regime | 0.13 ≲ Ro ≲ 2.0 | dJ / dt ∝ −Ω3 | Skumanich law: Prot ∝ t1/2 | Unsaturated αΩ dynamo (B ∝ Ro−1 ∝ Ω); standard gyrochronology sequence |
| Weakened Magnetic Braking (Old Stars) | Ro > Rocrit ≈ 2.0 | dJ / dt drops sharply | Prot stalls near constant value | Shift from dipole (n = 2–3) to higher-order multipole magnetic topology reduces RA (van Saders et al., 2016) |
Binary Orbital Evolution and the Cataclysmic Variable Period Gap
In close binary systems—such as cataclysmic variables (a white dwarf accreting from a Roche-lobe-filling K/M-dwarf donor), low-mass X-ray binaries (a neutron star or black hole primary), and W Ursae Majoris contact binaries—strong tidal forces lock the convective secondary star's spin period Prot to the binary orbital period Porb (Ωspin = Ωorb).
Because magnetic braking continuously removes spin angular momentum from the secondary star's convective envelope while tidal torques replenish that spin strictly at the expense of the binary orbit, magnetic braking extracts orbital angular momentum (dJorb / dt < 0), driving stable Roche-lobe mass transfer at dM / dt ~ 10−9–10−8 M☉ yr−1 for orbital periods Porb > 3 hours. In the standard disrupted magnetic braking model (Rappaport, Verbunt, & Joss, 1983), when the donor star is whittled down to M2 ≈ 0.25–0.30 M☉ at Porb ≈ 3.18 hours, it becomes fully convective, disrupting the tachocline dynamo and shutting down magnetic braking; the bloated donor shrinks inside its Roche lobe and detaches until gravitational-wave radiation shrinks the orbit across the 2–3 hour period gap to re-ignite mass transfer at Porb ≈ 2.15 hours.
Frequently Asked Questions
Why does magnetic braking extract orbital angular momentum from a binary if only one star has a magnetic wind?
Tidal friction in a close binary acts to keep the donor star's rotation synchronized with the orbital revolution (Ωspin = Ωorb). Whenever the donor's magnetized wind tries to slow the donor's spin, tidal torques promptly transfer angular momentum from the binary orbit into the donor's spin to restore synchronization, causing the two stars to spiral closer together and actually spin up to a shorter orbital period!
How did the Parker Solar Probe directly verify the solar Alfvén radius?
On April 28, 2021, NASA's Parker Solar Probe crossed below the Alfvén surface for the first time at a heliocentric distance of 18.8 solar radii (R☉ ≈ 0.087 au), measuring solar wind flow speeds slower than the local Alfvén wave speed (MA = vr / vA < 1) and directly confirming the sub-Alfvénic magnetic lever arm predicted by Weber and Davis (1967).
References
- Schatzman, E. (1962). "A theory of the role of magnetic activity during star formation." Annales d'Astrophysique, 25, 18–29. https://ui.adsabs.harvard.edu/abs/1962AnAp...25...18S
- Weber, E. J., & Davis, L., Jr. (1967). "The Angular Momentum of the Solar Wind." The Astrophysical Journal, 148, 217–227. https://doi.org/10.1086/149138
- Kawaler, S. D. (1988). "Angular momentum loss in low-mass stars." The Astrophysical Journal, 333, 236–247. https://doi.org/10.1086/166740
- Kasper, J. C., Klein, K. G., Lichko, E., Huang, J., Chen, C. H. K., Badman, S. T., ... & Zank, G. P. (2021). "Parker Solar Probe Enters the Magnetically Dominated Solar Corona." Physical Review Letters, 127(25), 255101. https://doi.org/10.1103/PhysRevLett.127.255101