Breaking gyrochronology through the collapse of coronal winds

Basic Information

• Paper ID: 2510.14164
• Title: Breaking gyrochronology through the collapse of coronal winds
• Authors: Michaël Lévesque, Paul Charbonneau (Université de Montréal)
• Classification: astro-ph.SR (Solar and Stellar Astrophysics), physics.comp-ph
• Publication Date: October 17, 2025 (Draft version)
• Paper Link: https://arxiv.org/abs/2510.14164

Abstract

Gyrochronology is a method for determining the age of field stars based on their rotation rates, but recent observations reveal that this method fails for stars older than the Sun. The conventional explanation attributes this failure to the shutdown of the stellar dynamo or changes in its mode, resulting in a sharp decrease in angular momentum loss from magnetized coronal winds. This paper explores an alternative possibility: the collapse of the wind itself due to reduced coronal heating. The study demonstrates that in the low coronal temperature limit, even at solar-like low rotation rates and coronal magnetic field strengths, magnetic centrifugal effects become important. When temperatures drop below approximately 1.5 MK, mass and angular momentum loss rates cannot be expressed as power-law functions of temperature or rotation rate.

Research Background and Motivation

Problem Background

1. Failure of Skumanich's Law: In 1972, Skumanich established the t^(-1/2) relationship between stellar rotation rate and age, which was widely verified in the 1-4.5 Gyr age range and became the foundation of stellar gyrochronology.
2. Observed Anomalies: Recent asteroseismic age determinations have revealed that some solar-type stars slightly older than the Sun rotate significantly faster than predicted by Skumanich's law, indicating a sharp decrease in angular momentum loss rate after solar age.
3. Limitations of Existing Explanations: Current mainstream explanations attribute this to the shutdown or mode change of the large-scale dynamo, but this implies that the current solar dynamo operates near a critical state, presenting theoretical difficulties.

Research Motivation

This paper proposes an alternative mechanism: explaining the sharp decrease in angular momentum loss through thermal collapse of coronal winds caused by reduced coronal heating energy input, providing a new physical perspective for understanding rotational evolution in older stars.

Core Contributions

1. Revealed complex dynamics of low-temperature coronal winds: Discovered that at low coronal temperatures approaching the hydrostatic limit, magnetic centrifugal effects become important, and angular momentum loss rate exhibits a non-monotonic relationship with coronal temperature.
2. Established linear relationship between mass loss and energy input: Demonstrated a strict linear relationship between mass loss rate and wind energy input across all coronal temperatures.
3. Quantified conditions for breaking Skumanich's law: By introducing a power-law relationship T₀ ∝ B_^σ, found that σ ≳ 1.5 is required to reproduce the observed breakdown of stellar gyrochronology.
4. Provided energetic constraints: Found that the required decrease in coronal energy input exceeds three orders of magnitude, which appears observationally unrealistic.

Methodology Details

Theoretical Framework

Weber-Davis Magnetohydrodynamic Wind Model

This paper employs the axisymmetric magnetohydrodynamic coronal wind model of Weber & Davis (1967):

Basic Assumptions:

• Axisymmetric (∂/∂φ = 0), equatorial plane (θ = π/2), steady-state (∂/∂t = 0) ideal MHD
• Radial magnetic field B_r(r) ∝ 1/r²
• Polytropic relation p ∝ ρ^α (1 ≤ α ≤ 5/3)

Key Physical Quantities:

Angular momentum conservation integral:

L = ru_φ - rB_φB_r/(μ₀ρu_r)

Alfvén radius condition:

L = Ωr²_A

Angular momentum loss rate:

J̇ = -8π/(3μ₀) × r₀⁴B²_{r0}ΩA⁻¹_{rA}

MacGregor-Brenner Rotational Evolution Model

A two-zone model describing core-envelope coupling in stellar interiors:

I_c(dΩ_c/dt) = -ΔJ/τ_c
I_e(dΩ_e/dt) = ΔJ/τ_c + J̇

where τ_c is the core-envelope coupling timescale.

Numerical Solution Methods

Numerical Solution of Weber-Davis Equations

The WD model is reduced to a six-dimensional nonlinear system:

• Regularity conditions at slow and fast magnetoacoustic critical points (4 equations)
• Energy conservation condition (2 equations)

Computational efficiency is enhanced through pre-computed four-dimensional parameter space grids and interpolation methods.

Hydrostatic Limit

The hydrostatic limit temperature for polytropic coronae:

T₀ ≤ (α-1)/α × GM⊙μm_p/(k_Br₀)

For α = 1.1, the limit temperature is T₀ = 0.91×10⁶ K.

Experimental Setup

Model Parameters

• Reference radius: r₀ = 1.15R⊙
• Polytropic index: α ∈ 1.05, 1.15
• Coronal temperature: T₀ ∈ 0.81, 4.50 MK
• Rotation rate: Ω/Ω⊙ ∈ 0.6, 25
• Magnetic field strength: B_ ∈ 1.67, 69.8 G

Dynamo Relation Assumption

A "standard" dynamo relation is adopted: B_ ∝ Ω

Temperature-Magnetic Field Relation

A power-law relationship is introduced:

T₀(Ω)/T⊙ = (Ω/Ω⊙)^σ

Experimental Results

Dynamical Characteristics of Low-Temperature Coronal Winds

Non-monotonic Angular Momentum Loss

At fixed solar rotation rate and magnetic field strength, the angular momentum loss rate J̇ exhibits non-monotonic behavior with coronal temperature T₀:

• For T₀ > 1.45×10⁶ K, J̇ increases as T₀ decreases
• Peak is reached at T₀ ≈ 1.35×10⁶ K
• Sharp decline for T₀ < 1.35×10⁶ K

Importance of Magnetic Centrifugal Effects

At low coronal temperatures (T₀ ≈ 1.23×10⁶ K):

• Magnetic centrifugal driving becomes significant, with asymptotic wind speeds reaching 858 km/s
• In contrast, pure thermal Parker wind speeds are only 162 km/s
• However, base flow speeds are essentially identical, with little difference in mass loss rates

Conditions for Breaking Skumanich's Law

Required Power-Law Index

To significantly deviate from Skumanich's law, σ ≳ 1.5 is required:

• At σ = 1.44, the rotation period at 8 Gyr increases by only ≈10%
• Even when mass loss rates decrease by several orders of magnitude, angular momentum loss rate changes remain limited

Sharp Decrease in Energy Input

The case σ ≥ 1.5 corresponds to:

• Coronal energy input decreasing by more than three orders of magnitude during the 5-8 Gyr period
• Declining from ≈10² W/m² to ≈10⁻¹ W/m²

Wind Energetics Analysis

Relationship Between Mass Loss and Energy Input

A strict linear relationship is found: Ṁ ∝ S (power input)

• This relationship holds across a broad range of (α, T₀) parameters
• Different (α, T₀) combinations can produce identical energy input and mass loss rates

Volumetric Heating Distribution

The reconstructed volumetric heating rate s(r) exhibits a steep power-law distribution: s(r) ∝ r⁻⁴

• Heating is concentrated at the coronal base
• Heating profiles from different parameter combinations can differ by more than two orders of magnitude

Related Work

Stellar Rotational Evolution Theory

• Skumanich (1972): Established the t⁻¹/² rotation-age relationship
• Durney (1972): Derived Skumanich's law from coronal wind theory
• MacGregor & Brenner (1991): Developed core-envelope coupling models

Magnetized Wind Theory

• Weber & Davis (1967): Laid the foundation for magnetized coronal wind theory
• Belcher & MacGregor (1976): Refined numerical solutions of the WD model
• Matt et al. (2012): Modern semi-empirical angular momentum loss models

Stellar Gyrochronology Observations

• van Saders et al. (2016): First report of rapid rotation in older stars
• Metcalfe et al. (2023): Further confirmation of gyrochronology breakdown
• Saunders et al. (2024): Latest observational evidence

Conclusions and Discussion

Main Conclusions

1. Complexity of Low-Temperature Coronal Winds: Near the hydrostatic limit, magnetic centrifugal effects become important, leading to non-monotonic behavior in angular momentum loss rate.
2. Decisive Role of Energy Input: Mass loss rate is entirely determined by wind energy input, exhibiting a strict linear relationship.
3. Limitations of Thermal Collapse Mechanism: Temperature decrease alone cannot fully explain the observed breakdown of stellar gyrochronology.
4. Necessity of Multiple Physical Mechanisms: Additional mechanisms such as dynamo shutdown or topological reconstruction are required.

Limitations

1. Simplified Magnetic Field Geometry: The WD model employs a dipole magnetic field, while actual stellar fields are more complex
2. Polytropic Approximation: Neglects detailed coronal heating physics
3. Empirical Parameter Relations: The T₀ ∝ B_^σ relationship lacks theoretical foundation
4. Insufficient Observational Constraints: Limited observational data on coronal parameters of older stars

Future Directions

1. More Realistic Magnetic Field Models: Consider multipole fields and active region effects
2. Coronal Heating Physics: Incorporate specific mechanisms such as magnetic reconnection and wave heating
3. Multidimensional MHD Simulations: Verify the applicability of axisymmetric models
4. Observational Verification: Seek direct evidence for decreased coronal activity in older stars

In-Depth Evaluation

Strengths

1. Deep Physical Insights: Reveals counterintuitive behavior of low-temperature magnetized winds, providing new perspectives for understanding stellar wind dynamics
2. Advanced Numerical Methods: Establishes an efficient framework for solving WD equations, laying the foundation for subsequent research
3. Comprehensive Energetic Analysis: First systematic analysis of the relationship between coronal wind energetics and stellar rotational evolution
4. Objective Results Presentation: Honestly reports the insufficiency of the thermal collapse mechanism, avoiding overinterpretation

Weaknesses

1. High Degree of Model Simplification: The axisymmetric assumption of the WD model may not apply to complex stellar magnetic fields
2. Lack of Physical Basis for Parameter Relations: The T₀ ∝ B_^σ relationship is purely phenomenological, lacking theoretical support
3. Weak Observational Constraints: Insufficient discussion of observational evidence for coronal temperature changes in older stars
4. Limited Discussion of Alternative Mechanisms: Insufficient exploration of other possible physical mechanisms

Impact

1. Theoretical Contribution: Provides a new physical mechanism for gyrochronology breakdown, enriching the theoretical framework
2. Methodological Value: The established numerical methods can be applied to broader stellar wind research
3. Observational Guidance: Provides theoretical expectations and constraints for future stellar coronal observations
4. Cross-disciplinary Significance: Connects multiple fields including stellar interior structure, magnetic field generation, and atmospheric physics

Applicable Scenarios

This research is applicable to:

• Rotational evolution modeling of solar-type stars (G-type main-sequence stars)
• Improvement and correction of gyrochronology methods
• Theoretical research on stellar wind physics
• Observational interpretation of the relationship between stellar magnetic activity and age

References

The paper cites numerous important works, including:

• Skumanich (1972): Foundational work in stellar gyrochronology
• Weber & Davis (1967): Theoretical foundation of magnetized wind theory
• MacGregor & Brenner (1991): Stellar interior coupling models
• van Saders et al. (2016): Observational discovery of gyrochronology breakdown
• Metcalfe et al. (2023): Latest observational confirmation

This paper provides important theoretical contributions to understanding the anomalous rotational evolution of older stars. While a single thermal collapse mechanism cannot fully explain the observed phenomena, the physical processes it reveals lay the foundation for constructing a more complete theoretical framework.