X-ray Luminosity–Age Relation in Stars

• X-ray luminosity/age relationship is a correlation where stellar magnetic activity and rotation drive the X-ray output decay over time.
• Empirical studies using eROSITA, Chandra, and XMM-Newton data reveal a power-law decay with significant scatter influenced by mass, metallicity, and binary evolution.
• Precise age-dating of stars requires multi-epoch observations to account for intrinsic variability and systematic uncertainties in X-ray measurements.

The relationship between X-ray luminosity ($L_X$) and stellar age is a fundamental diagnostic of magnetic activity evolution in stars, populations, and galaxies. The X-ray/age relation encapsulates the effects of stellar rotation, magnetic dynamo efficiency, binary evolution, and high-energy feedback processes, and serves as a clock for age-dating individual stars, stellar populations, and environments hosting exoplanets. Empirical calibration of this relation has advanced with precise ages from asteroseismology and cluster membership, as well as large X-ray surveys (eROSITA, Chandra, XMM-Newton), revealing both broad correlations and substantial astrophysical scatter.

1. Physical Foundations and Formalism

Stellar X-ray emission originates primarily from magnetically heated coronae in single stars and accretion onto compact objects (primarily in binaries) in populations. For single field stars and cluster members, magnetic activity is regulated by rotation and interior structure. Major formalisms include:

• Surface-area normalized relation: log $L_{X,n} = m\,\log \tau + b$ with $L_{X,n} \equiv L_X/(R_*/R_\odot)^2$, age $\tau$ in Gyr, intercept $b$, and slope $m$ ($\equiv \beta$).
• Fractional luminosity: $R_X = L_X / L_{\rm bol}$, with $R_{X,\rm sat}$ indicating a saturated "plateau" (typically $10^{-3.1}$ to $L_{X,n} = m\,\log \tau + b$0 across spectral types for young rapid rotators (Jackson et al., 2011)).
• Population scaling: $L_{X,n} = m\,\log \tau + b$1 for integrated binaries and field populations, linked via $L_{X,n} = m\,\log \tau + b$2 or broken power-law luminosity functions (XLFs).

In X-ray binary populations, two principal regimes are recognized:

• High-mass X-ray binaries (HMXBs) dominate at $L_{X,n} = m\,\log \tau + b$3 Myr, tracking recent star formation.
• Low-mass X-ray binaries (LMXBs) emerge at $L_{X,n} = m\,\log \tau + b$4 Gyr, with normalization and slope steepening as donor-mass declines (Lehmer et al., 2024, Lehmer et al., 2017, Lehmer et al., 2014).

2. Empirical Results for Single and Main-Sequence Stars

Age and Mass Dependence

Aldarondo Quiñones et al. (Quiñones et al., 10 Nov 2025) provide updated eROSITA-based calibrations for main-sequence F–M stars older than 1 Gyr:

• Best-fit slope: $L_{X,n} = m\,\log \tau + b$5, intercept $L_{X,n} = m\,\log \tau + b$6 (surface-area normalized).
• Intrinsic scatter: $L_{X,n} = m\,\log \tau + b$7 dex (factor 2–3 in $L_{X,n} = m\,\log \tau + b$8).
• With mass term: log $L_{X,n} = m\,\log \tau + b$9, $L_{X,n} \equiv L_X/(R_*/R_\odot)^2$0, and $L_{X,n} \equiv L_X/(R_*/R_\odot)^2$1 steepens to $L_{X,n} \equiv L_X/(R_*/R_\odot)^2$2.

A prior study by Booth et al. (Booth et al., 2017) found a steeper decay in old stars ($L_{X,n} \equiv L_X/(R_*/R_\odot)^2$3), but eROSITA data indicate that accounting for variability ("jitter") the decay slope matches that of younger stars ($L_{X,n} \equiv L_X/(R_*/R_\odot)^2$4 to $L_{X,n} \equiv L_X/(R_*/R_\odot)^2$5, Jackson et al. 2012 (Jackson et al., 2011)).

Astrophysical implications:

• X-ray luminosity decays with age, but intrinsic scatter and variability dominate, limiting age precision to ≳30–50%.
• Shallow decay plus scatter implies $L_{X,n} \equiv L_X/(R_*/R_\odot)^2$6 should not be used in isolation for field star ages beyond 1 Gyr (Quiñones et al., 10 Nov 2025).

Activity "Plateau" and Dynamo Evolution

Young ($L_{X,n} \equiv L_X/(R_*/R_\odot)^2$71 Gyr) and pre-main-sequence (PMS) stars show a saturated X-ray regime with $L_{X,n} \equiv L_X/(R_*/R_\odot)^2$8 and $L_{X,n} \equiv L_X/(R_*/R_\odot)^2$9 nearly constant:

• Saturation values range from $\tau$0 (late K) to $\tau$1 (early F) (Jackson et al., 2011).
• Saturation timescale $\tau$2–200 Myr, not monotonic with spectral type.
• After saturation, fractional luminosity decays as $\tau$3, $\tau$4.

For PMS stars, X-ray output is constant for $\tau$5 Myr; decay ($\tau$6) steepens with mass and radiative-core development (Getman et al., 2022):

• $\tau$7 for 0.75–1 M$\tau$8; $\tau$91.8 for 1–3.5 M$b$0; up to $b$13–4 for 3.5–7 M$b$2.

3. Population and Binary Scaling Relations

Galaxy and Cluster Populations

Star-forming galaxies and open clusters reveal a pronounced decline in integrated X-ray output per unit mass with stellar population age, traced both by individual source XLFs and by aggregate $b$3:

• Empirical frameworks (e.g., (Lehmer et al., 2024)): XLF normalization per mass drops by 2–3 dex from 10 Myr to 10 Gyr; metallicity modulates this decline (slower at low $b$4).
• Quantitatively, (Gilbertson et al., 2021) find $b$5 drops from $b$632.4 (young; $b$710 Myr) to $b$829 for $b$93 Gyr—a $m$0 decline.
• M51 spatially-resolved analysis (Lehmer et al., 2017): $m$1 declines by $m$23 dex over 10 Myr–10 Gyr, and the bright-end XLF slope steepens from $m$31.4 (HMXB-like) to 3 (LMXB-dominated).

LMXBs in Early-Type Galaxies

Zhang et al. (Zhang et al., 2012) systematically show $m$450% more LMXBs per unit mass in older ($m$56 Gyr) compared to younger galaxies, with a two-parameter scaling: $m$6 where $m$7 is age (Gyr) and $m$8 is globular cluster frequency. The cumulative field data (Lehmer et al., 2014) reveal field LMXB $m$9 declines as $\equiv \beta$0; excesses in young galaxies are factors $\equiv \beta$1–$\equiv \beta$2 above old ones, in agreement with population synthesis models (Fragos et al.).

4. Spectral Type, Mass, and Metallicity Dependencies

• Spectral type: Faster $\equiv \beta$3 decay for higher-mass stars, longer "saturation" for late-Ms (e.g. $\equiv \beta$4 drops with $\equiv \beta$5 post 2.1 Gyr for M2.5–6.5 vs. $\equiv \beta$6 for early Ms (Engle, 2023)).
• Metallicity: Lower $\equiv \beta$7 enhances normalization and delays the decline in $\equiv \beta$8 due to reduced stellar wind mass-loss and higher donor mass in binaries (Lehmer et al., 2024, Gilbertson et al., 2021).
• Mass: Weak $\equiv \beta$9 evidence for steeper $R_X = L_X / L_{\rm bol}$0 decay in high-mass stars ($R_X = L_X / L_{\rm bol}$1 in Aldarondo Quiñones et al. (Quiñones et al., 10 Nov 2025)).

5. Astrophysical Variability and Systematic Uncertainties

Intrinsic variability ("jitter") in $R_X = L_X / L_{\rm bol}$2 (magnetic cycles, rotational modulation, flares) introduces an irreducible scatter (0.3–0.5 dex) above formal errors (Quiñones et al., 10 Nov 2025):

• This scatter exceeds the age-driven trend in old stars, limiting age inference for individuals to $R_X = L_X / L_{\rm bol}$3–$R_X = L_X / L_{\rm bol}$4 uncertainty.
• Systematics: Instrument cross-calibration; age estimation errors (asteroseismology, WD companions); unresolved binaries and contaminating flux.
• For population studies: completeness corrections, cosmic variance, and globular cluster LMXB contributions affect normalization.

6. Applications, Limitations, and Future Prospects

• X-ray luminosity is a robust ensemble age diagnostic, notably for population studies (e.g., field galaxies, clusters), but is too variable for precise single-star ages at $R_X = L_X / L_{\rm bol}$51 Gyr except when multi-epoch data are available (Quiñones et al., 10 Nov 2025).
• Calibration of exoplanet evaporation histories must include the early saturated phase, which delivers the bulk (≥75%) of irradiation-driven mass loss within the first Gyr (Jackson et al., 2011).
• Multi-dimensional frameworks that incorporate metallicity, mass, and SFH (e.g. (Lehmer et al., 2024)) are now available and required for accurate modeling.
• Larger time-resolved X-ray samples, especially multi-epoch monitoring, are needed to deconvolve variability from intrinsic age trends and tighten precision.
• Combination with complementary indicators (rotation, chromospheric lines, asteroseismology) yields the strongest constraints for stellar dating and environmental effects.

7. Summary Table: Key Empirical Age–$R_X = L_X / L_{\rm bol}$6 Relations

Regime	Best-fit Relation	Slope/Exponent (β)	Typical Scatter (dex)
PMS plateau (0–5 Myr)	$R_X = L_X / L_{\rm bol}$7 (0.75–1 M$R_X = L_X / L_{\rm bol}$8)	~0	$R_X = L_X / L_{\rm bol}$90.3
Early MS (6–625 Myr, G/K/M)	$R_{X,\rm sat}$0	G: –0.61; K: –0.82; M: –0.4	0.2–0.4
Field stars ($R_{X,\rm sat}$11 Gyr, F–M)	$R_{X,\rm sat}$2	$R_{X,\rm sat}$3	$R_{X,\rm sat}$4
LMXBs in field galaxies	$R_{X,\rm sat}$5	–0.9	$R_{X,\rm sat}$60.2–0.3
Integrated galaxy $R_{X,\rm sat}$7	Decline by $R_{X,\rm sat}$83 dex, 10 Myr–10 Gyr

The X-ray luminosity/age relationship is thus characterized by a plateau at early times, a subsequent power-law decay with an exponent near –1 for field stars and steeper in select regimes, and substantial astrophysical jitter. Population effects, metallicity, mass, and intrinsic variability all modulate this decay and must be accounted for in precise applications. The relation remains a core input for stellar astrophysics, galactic evolution modeling, and exoplanet habitability studies.