X-ray Evolution of Young Stars: Early Dimming and Coronal Softening in Solar-Mass Stars with Implications for Planetary Atmospheres

TL;DR

The paper addresses how X-ray luminosity and coronal temperatures evolve in solar-mass and slightly lower-mass stars up to about 750 Myr and how this high-energy radiation shapes early planetary atmospheres. It uses new Chandra observations of five ~45–100 Myr clusters, Gaia-based memberships, and archival ROSAT/Chandra data for older clusters, applying joint X-ray–Gaia analyses and Kaplan–Meier statistics to derive mass- and age-stratified L_X trends and coronal temperatures. The key finding is a mass-dependent decay with solar-mass stars dimming faster and softening their coronal spectra by ~100 Myr, while sub-solar stars retain hotter coronae longer; these trends differ from prior XUV-rotation models, suggesting lower atmospheric mass loss and altered photochemistry for close-in planets. The results have significant implications for planetary evolution, habitability windows, and预biotic chemistry, indicating a need to recalibrate models of early atmospheric erosion and high-energy ionization histories.

Abstract

X-ray and ultraviolet (XUV) emission from young stars plays a critical role in shaping the evolution of planetary atmospheres and the conditions for habitability. To assess the long-term impact of high-energy stellar radiation, it is essential to empirically trace how X-ray luminosities and spectral hardness evolve during the first ~<1 Gyr, when atmospheric loss and chemical processing are most active. This study extends the X-ray activity-mass-age analysis of <25 Myr stars by Getman et al. (2022) to ages up to 750 Myr, using Gaia-based cluster memberships, new Chandra observations of five rich open clusters (~45--100 Myr), and archival ROSAT and Chandra data for three older clusters (~220--750 Myr). We find a mass-dependent decay in X-ray luminosity: solar-mass stars undergo a far more rapid and sustained decline, accompanied by coronal softening and the disappearance of hot plasma by ~100 Myr, compared to their lower-mass siblings. These trends in solar-mass stars are likely linked to reduced magnetic dynamo efficiency and diminished ability to sustain large-scale, high-temperature coronal structures. The trends are significantly stronger than predicted by widely used XUV-rotation-age relations. The revised trends imply systematically lower rates of atmospheric mass loss and water photolysis, as well as altered ionization environments and chemical pathways relevant to the formation of prebiotic molecules, for planets in close orbits around solar analogs. These effects persist throughout at least the ~<750 Myr interval probed in this study.

Figures (11)

• Figure 1: Adaptively smoothed, low-resolution mosaic images of the five open clusters observed with Chandra ACIS-I, shown in the $(0.5-8)$ keV energy band, revealing a few thousand of the brighter X-ray sources across the fields.
• Figure 2: Summary of member identification and properties for NGC 6242. (a,b) Final likely members: non-X-ray (Groups 1 and 2, red), X-ray (Groups 5 and 6, green), and additional members outside the Chandra field from Cantat-Gaudin2020 ( purple), compared with likely non-members (grey). Shown in 3D space defined by proper motion distance from the cluster's motion centroid, parallax, and photometric offset from the empirical CMD sequence. Insets show zoomed-in member views. (c) CMD for Groups 1, 2, 5, and 6 plus Cantat-Gaudin2020 members (black points with $G_{BP} - G_{RP}$ errors in grey). Overlaid are the best-fit PARSEC isochrone (red), bounding isochrones (cyan/blue), and the smoothed empirical sequence from locfit (dashed green). (d) IMF histogram for Groups 1, 2, 5, and 6 across the Chandra field (black line with 95% Poisson error bars). Theoretical IMF from Maschberger2013 shown in grey. The dashed vertical line marks $M = 0.75$ M$_{\odot}$; the legend lists the observed $M > 0.75$ M$_{\odot}$ count and extrapolated population down to 0.1 M$_{\odot}$. (e) Stacked Chandra spectrum of Group 5 and 6 stars with best-fit two-temperature plasma model. (f) X-ray luminosity vs. angular distance from cluster center. Group 5 (green), Group 6 (blue), and likely contaminants (Groups 7–10, grey) are shown. A red local regression to the lower envelope is used to estimate X-ray upper limits for undetected (Groups 1 and 2) cluster members. Analogous figures for other clusters appear in the Appendix (Figures \ref{['fig:mem_TRUMPLER3']}-\ref{['fig:mem_NGC2301']}).
• Figure 3: (a, d, g) Gaia CMDs of all cluster members with $Prob > 0.5$ in the $>150$ Myr clusters NGC 6475, M37, and Praesepe HuntReffert2023. The best-fit PARSEC isochrone to the MSTO is shown in green, with bounding isochrones surrounding the MSTO in red and blue. (d, e, h) IMF histograms of the same cluster members (black), with the theoretical IMF from Maschberger2013 shown in grey. Dashed vertical lines mark $0.75$ and $1.2$ M$_{\odot}$, and the legend lists both the observed number of stars with $M>0.75$ M$_\odot$ and the extrapolated population down to $0.1$ M$_{\odot}$. (c, f, i) Spatial distributions of Gaia cluster members with $Prob > 0.5$ in the mass ranges $0.75-0.9$ M$_{\odot}$ (green diamonds) and $0.9-1.2$ M$_{\odot}$ (red boxes), overlaid with X-ray detections (blue $\times$) from ROSAT (NGC 6475 and Praesepe) and Chandra (M37). For Praesepe, the ROSAT PSPC catalog of Randich1995 (blue $\times$), limited to stars considered members at the time, is supplemented with additional ROSAT sources from the Second ROSAT PSPC catalog (black $+$) following Gaia-based membership updates. Only X-ray and non-X-ray members located within the orange contours are used in the X-ray temporal evolution analysis (Figure \ref{['fig:lx_vs_t_main']}) and are recorded in Table \ref{['tab:tab_ngc6475_m37_praesepe']}. The demarcating orange circles for NGC 6475 and M37 have radii of $R=18\arcmin$ and $5\arcmin$, respectively, while the orange square for Praesepe measures $70 \times 70$ arcmin$^2$.
• Figure 4: Temporal evolution of X-ray luminosity (in the $0.5-8$ keV band) across four stellar mass strata. The corresponding mass bins are indicated in the panel legends. Younger clusters from Getman2022 are shown as brown points ($t<5$ Myr) and black points ($7-22$ Myr). The five new ($30-150$ Myr) Chandra clusters (Table \ref{['tab:cluster_props']}) are marked by dark green, cyan, purple, orange, and blue points, while three older clusters ($>150$ Myr) from an extended sample with previously published X-ray detections (Table \ref{['tab:extra3cluster_props']}) are represented by large pink points. The latter appear only in the first two panels. The legends list the cluster names and the total number of stars (X-ray detections and non-detections) used in the Kaplan-Meier (KM) estimation of $L_X$, as well as the best-fit parameters from weighted least squares regression fits to the relation $L_X = a \times t^b$, performed over the $7-100$ Myr baseline using the R lm function Sheather09. The regression fits are shown as black lines with 68% confidence intervals shaded in gray; and the small $p$-values (except for panel (d)) reported in the legends indicate strong evidence against the null hypothesis of a zero slope. In addition, for panel (b) only, the locfit best-fit (pink curve) shows a declining trend of X-ray luminosity over the extended $7-750$ Myr baseline. (a,b) Panels show the evolution of $L_X$ for lower-mass stars. Median $L_X$ values (50% quartiles) are shown with 68% bootstrap confidence intervals: circles or diamonds for open clusters and boxes for younger MYStIX/SFiNCs clusters. The 25% and 75% quartiles, plotted only for samples with more than 15 stars and where permitted by the KM estimator, are marked by green triangles. For comparison, the model $L_X$ predictions in the $0.5-7$ keV band from Johnstone2021 are shown: red points with error bars for medians and their 68% confidence intervals, and red lines for the interquartile range (25%-75%). (c,d) Panels show analogous information for higher-mass stars, but only the 85% $L_X$ quartiles (as permitted by the KM estimator) are plotted, with their associated 68% bootstrap confidence intervals.
• Figure 5: (a-f) Empirical cumulative distribution functions (ECDFs) of source' X-ray median photon energies, stratified by stellar mass and stacked for clusters of similar ages and observation epochs. Stellar mass groups are color-coded: solar-mass (red), sub-solar (green), and lower-mass (gray). Figure legends indicate the number of stars in each group, along with $p$-values from Anderson-Darling two-sample tests (using the ad.test function from the R package kSamples) for pairwise comparisons between sub-solar and solar-mass samples (black text), and sub-solar and lower-mass samples (gray text). (g) For the oldest clusters, NGC 2353 and NGC 2301 ($\sim 100$ Myr), ECDFs of the observed, background-subtracted photon energies are shown for solar-mass (red) and sub-solar (green) stars. These distributions are overlaid with simulated CDFs from one-temperature coronal plasma models with varying plasma temperatures. The figure legend specifies the plasma temperatures corresponding to each model curve, as well as the total number of net counts used to construct the observed photon energy distributions.