Gyrochronology

Gyrochronology (from the Ancient Greek gyros, "rotation," chronos, "time," and -logia, "study") is an empirical astrophysical age-dating method formulated by Sydney A. Barnes in 2003 and 2007 that derives the ages of low-mass main-sequence stars (spectral types late-F, G, K, and M, with masses M ≲ 1.3 M_☉) from two directly measurable observables: their axial rotation period P_rot and their intrinsic color index (such as Johnson B − V or Gaia G_BP − G_RP, which serves as a proxy for stellar mass and convective turnover time). Whereas traditional isochrone fitting on the Hertzsprung–Russell diagram yields large age uncertainties (±50–100%) for main-sequence field G, K, and M dwarfs because their luminosities and effective temperatures change by less than a few percent per gigayear during core hydrogen burning, their rotation periods lengthen steadily via magnetic braking according to the Skumanich law (P_rot ∝ t^1/2). Within astronomy, gyrochronology calibrated against coeval open clusters by the Kepler, K2, and TESS space missions provides ages with 10% to 15% precision for exoplanet host stars and Galactic disk populations between ~100 million and ~4 billion years old.

Why Isochrones Fail and Rotation Succeeds on the Lower Main Sequence

For stars cooler than the Kraft break (T_eff ≲ 6,250 K, M ≲ 1.3 M_☉), the nuclear main-sequence lifetime spans 4 Gyr (for an F6V star) to over 100 Gyr (for an M dwarf). Because nuclear evolution is so slow on the lower main sequence, a 1-gigayear-old K dwarf and a 6-gigayear-old K dwarf occupy virtually identical positions on a color–magnitude diagram.

By contrast, every star below the Kraft break possesses an outer convective envelope that powers a turbulent magnetic dynamo and launches an ionized coronal Parker wind. As coronal magnetic field lines force the wind to corotate out to the Alfvén radius R_A, magnetic braking extracts angular momentum at a rate proportional to the cube of the rotation rate (dJ / dt ∝ −Ω³). Because this cubic negative feedback decelerates initially fast rotators far more strongly than initially slow rotators, stars of a given mass erase their birth spin dispersion within 100 to 600 Myr and converge onto a single deterministic surface in three-dimensional (Period – Color/Mass – Age) space.

Barnes's C-Sequence, Gap, and Asymptotic I-Sequence

In his 2003 classification of open-cluster color–period diagrams (plotting rotation period P_rot against intrinsic color B − V), Sydney Barnes identified two distinct morphological branches separated by a transitional rotational gap:

1. The C-Sequence (Convective / Fast Sequence): Populated by young stars (t ≲ 50–150 Myr for G dwarfs; up to ~500 Myr for late-K and M dwarfs) with short rotation periods (P_rot ≈ 0.2 to 2 days) whose magnetic dynamos operate in the saturated regime (Rossby number Ro ≡ P_rot / τ_c < 0.13). In this regime, angular momentum loss scales only linearly with rotation (dJ / dt ∝ −Ω), and the outer convective envelope is rotationally decoupled from the radiative core.
2. The I-Sequence (Interface / Slow Sequence): Populated by stars that have spun down into the unsaturated dynamo regime (Ro ≳ 0.13), where the radiative core and convective envelope are tightly coupled across the tachocline interface and angular momentum loss obeys dJ / dt ∝ −Ω³. Once a star crosses the transitional gap from the C-sequence onto the I-sequence, its subsequent rotation period is a tight, single-valued function of its age t and color/mass.

Separable Empirical Factorizations (Barnes 2007; Mamajek & Hillenbrand 2008)

For stars that have converged onto the I-sequence, Barnes (2007) expressed the rotation period P_rot (in days) as a separable product of an age function g(t) (matching the Skumanich law tⁿ with n ≈ 0.5) and a color/mass dependence f(B − V):

P_rot(t, B − V) = g(t) × f(B − V) = tⁿ × a ( B − V − c )^b

Inverting this equation yields the gyrochronological age t_gyro (in millions of years):

t_gyro = [ P_rot / ( a [ B − V − c ]^b ) ]^{1 / n}

where c ≈ 0.40–0.495 marks the blueward color boundary of the Kraft break (where convective envelopes vanish and f → 0). Standard empirical calibrations by Barnes (2007), Eric E. Mamajek and Lynne A. Hillenbrand (2008), and Sören Meibom et al. (2009–2015) yield:

• Age index: n = 0.519 ± 0.007 (or n = 0.566 in Mamajek & Hillenbrand 2008)
• Color amplitude: a = 0.773 ± 0.011 (or a = 0.407 for c = 0.495)
• Color power-law slope: b = 0.601 ± 0.024 (or b = 0.325 for c = 0.495)

Non-Separable Convective Turnover Time Models (Barnes 2010)

To unify both the young C-sequence and the mature I-sequence without assuming strict separability, Barnes (2010) and Barnes & Kim (2010) formulated a differential spin-down equation expressed directly in terms of the convective turnover timescale τ_c(M) (which increases from ~4 days in late-F stars to ~15 days in the Sun and >70 days in M dwarfs) and the dimensionless Rossby number Ro = P_rot / τ_c:

t = ( τ_c / k_C ) ln( P_rot / P₀ ) + ( k_I / [ 2 τ_c ] ) ( P_rot² − P₀² )

where k_C ≈ 0.646 day Myr⁻¹ and k_I ≈ 452 Myr day⁻¹ are universal dimensionless constants. For young fast rotators (P_rot small), the first logarithmic term dominates (C-sequence); for older stars (P_rot ≫ P₀), the quadratic second term dominates, recovering P_rot ∝ (τ_c t)^1/2 on the I-sequence.

Space photometry from Kepler, K2, and TESS combined with Gaia parallaxes has uncovered two physical departures from simple power-law gyrochronology that are now incorporated into modern gyrochronology software packages (such as gyro-interp by Luke Bouma et al., 2023):

1. K-Dwarf Temporary Spin-Down Stalling (~0.6 to 1.5 Gyr): Comparisons between Praesepe (670 Myr) and NGC 6811 (1.0 Gyr) by Jason Curtis, Marcel Agüeros, and Collaborators (2019–2020) showed that while G dwarfs spin down uninterruptedly as t^1/2, lower-mass K dwarfs (M ≈ 0.6–0.85 M_☉) experience a temporary ~700 Myr epoch of stalled surface spin-down. During this phase, angular momentum stored in the rapidly rotating radiative core is transported outward to the convective envelope at a rate matching wind losses at the surface; once core–envelope recoupling completes after ~1.5 Gyr, K dwarfs resume Skumanich spin-down.
2. Weakened Magnetic Braking at High Rossby Number (Ro > Ro_crit ≈ 2.0): By combining Kepler asteroseismic ages with rotation periods for old field stars, Jennifer van Saders et al. (2016) and Oliver Hall et al. (2021) discovered that once a star spins down past a threshold Rossby number Ro_crit ≈ 2.0 (reached at t ≈ 2–3 Gyr in F dwarfs, ~4.5 Gyr in G2V stars like the Sun, and >7 Gyr in K dwarfs), the coronal magnetic field shifts from a large-scale dipole to higher-order multipoles, reducing the Alfvén lever arm R_A and causing rotation periods to stall near a constant ceiling for the remainder of the main sequence.

How are stellar rotation periods measured for gyrochronology?

As a magnetically active star rotates, dark photospheric starspots and bright faculae cross the visible hemisphere, producing periodic quasi-sinusoidal brightness modulations of 0.01% to 5% in time-series light curves from Kepler, K2, TESS, or ground-based surveys, whose periodogram peak directly yields the surface rotation period P_rot.

Why can't standard gyrochronology be applied to close binary stars?

In close binaries with orbital periods shorter than ~10–15 days, gravitational tidal torques synchronize the stars' spins with the orbit, overriding natural wind-driven spin-down and making the rotation period reflect binary orbital dynamics rather than stellar age.

1. Barnes, S. A. (2003). "On the Rotational Evolution of Solar- and Late-Type Stars, Its Magnetic Origins, and the Possibility of Stellar Gyrochronology." The Astrophysical Journal, 586(1), 464–479. https://doi.org/10.1086/367639
2. Barnes, S. A. (2007). "Ages for Illustrative Field Stars Using Gyrochronology: Viability, Limitations, and Errors." The Astrophysical Journal, 669(2), 1167–1189. https://doi.org/10.1086/519295
3. Mamajek, E. E., & Hillenbrand, L. A. (2008). "Improved Age Estimation for Solar-Type Dwarfs Using Activity-Rotation Diagnostics." The Astrophysical Journal, 687(2), 1264–1293. https://doi.org/10.1086/591785
4. Curtis, J. L., Agüeros, M. A., Matt, S. P., Covey, K. R., Douglas, S. T., Angus, R., ... & Wright, J. T. (2020). "When Do Stalled Stars Resume Spinning Down? Advancing Gyrochronology with Ruprecht 147." The Astrophysical Journal, 904(2), 140. https://doi.org/10.3847/1538-4357/abbf58
5. Bouma, L. G., Palumbo, E. K., & Hillenbrand, L. A. (2023). "The Empirical Limits of Gyrochronology." The Astrophysical Journal Letters, 947(1), L3. https://doi.org/10.3847/2041-8213/acc589