Probing habitable regions with SRG/eROSITA

Stellar high-energy radiation is a key driver of atmospheric erosion and evolution in exoplanets, directly affecting their long-term habitability. We present a comprehensive study on stellar high-energy radiation and its impact on exoplanetary atmospheres, leveraging data from the \textit{SRG/eROSITA} all-sky survey. Our sample consists of 3750 main-sequence stars identified by cross-matching with \textit{Gaia} DR3. Utilizing X-ray spectral fits from the \textit{eROSITA} catalog, we computed X-ray ($L_X$) and combined extreme-ultraviolet (EUV) luminosities ($L_{\mathrm{EUV}}$), which we used to derive XUV fluxes at the habitable zone ($F_{\mathrm{XUV,HZ}}$). We find that the majority of stars in our sample are significantly more XUV-active than the Sun, with habitable zone fluxes ranging from $10^0$ to $10^5$ erg~cm$^{-2}$~s$^{-1}$. The ratio of $L_{\mathrm{XUV}}/L_{\mathrm{bol}}$ is found to be higher for cooler, magnetically active stars, highlighting their potentially hazardous nature for planetary atmospheres. Applying the energy-limited escape model, we computed atmospheric mass-loss rates for hypothetical earth-like planets located at the habitable zone of each star. We also present local maps for distances up to $500$~pc of the average XUV flux, revealing ``hazard zones’’ where stellar radiation could significantly influence planetary atmospheric evolution. This work demonstrates the power of X-ray surveys in constraining the high-energy environments of exoplanets and underscores the critical role of stellar activity in planetary habitability.

💡 Research Summary

This paper presents a comprehensive assessment of the high‑energy radiation environments around main‑sequence stars and evaluates the implications for planetary habitability, using data from the SRG/eROSITA all‑sky X‑ray survey combined with Gaia DR3 astrometry. The authors begin by constructing a clean sample of 3,750 main‑sequence stars. They start from the eRASS1 coronal source catalog, remove objects flagged for optical loading or non‑coronal emission, and then identify main‑sequence members by comparing Gaia absolute magnitudes (derived from apparent G‑band magnitudes and parallactic distances) with the empirical dwarf sequence of Pecaut & Mamajek, allowing a tolerance of one magnitude. This yields a homogeneous set of stars spanning a wide range of spectral types (F through M) and bolometric luminosities.

X‑ray spectral analysis follows the methodology of Gatuzz et al. (2024). Each source is fitted with a suite of phenomenological models (absorbed power‑law, single‑temperature APEC, black‑body, double‑temperature APEC, and a double‑temperature model with free abundances). The double‑thermal models (D and E) provide the best statistical fits and physically realistic coronal parameters. Unabsorbed fluxes (F_X) are extracted from the preferred model, and X‑ray luminosities are computed via L_X = 4π d² F_X, where distances d come from Gaia parallaxes. To ensure realistic coronal outputs, a luminosity cut of 10²⁶–10³² erg s⁻¹ is applied, which brackets the known range from solar‑like quiet stars to saturated, highly active coronae.

Because the extreme‑ultraviolet (EUV) band (0.013–0.1 keV) cannot be observed directly, the authors adopt the empirical scaling relation of Sanz‑Forcada et al. (2011): log L_EUV = (4.80 ± 1.99) + (0.860 ± 0.073) log L_X (with L_X defined in the 0.1–2.4 keV band). This yields EUV luminosities for every star, which are then summed with the X‑ray luminosities to obtain the combined high‑energy output L_XUV = L_X + L_EUV. The distribution of L_XUV spans 10²⁷–10³³ erg s⁻¹, and the authors explore its dependence on Gaia color (BP–RP), absolute magnitude, and stellar mass, confirming that more massive, hotter stars tend to have higher absolute L_XUV, whereas cooler K‑ and M‑type dwarfs exhibit larger L_XUV/L_bol ratios, reflecting a higher fraction of their total energy emitted in the high‑energy bands.

To translate these stellar outputs into the irradiation experienced by a planet in the habitable zone (HZ), the paper employs the effective stellar flux prescription of Kopparapu et al. (2014). The effective flux S_eff is expressed as a fourth‑order polynomial in the deviation of the stellar effective temperature from the solar value (T_* = T_eff – 5780 K), with separate coefficient sets for the inner and outer HZ boundaries. The HZ distance is then d_HZ = √(L_bol / (L_⊙ S_eff)) AU. Using this distance, the XUV flux at the HZ is calculated as F_XUV,HZ = L_XUV / (4π d_HZ²). The resulting fluxes range from 10⁰ to 10⁵ erg cm⁻² s⁻¹, substantially exceeding the present solar XUV flux at Earth (≈4 erg cm⁻² s⁻¹). Low‑mass stars (M ≲ 0.6 M_⊙) deliver the highest F_XUV,HZ because their HZs lie very close to the star, while more massive stars, despite higher intrinsic L_XUV, have more distant HZs and thus lower fluxes.

The authors then estimate atmospheric mass‑loss rates for a hypothetical Earth‑analog (M_p = M_⊕, R_p = R_⊕) using the energy‑limited escape formula: Ṁ = ε π R_p³ F_XUV,HZ / (G M_p K). Here ε is the heating efficiency (adopted 0.1–0.3), and K accounts for the reduction of the gravitational potential at the Roche lobe. The calculated Ṁ values typically lie between 10⁹ and 10¹¹ g s⁻¹ for planets around M‑type stars, indicating that sustained high‑energy irradiation could strip a terrestrial atmosphere on timescales of a few hundred million years, potentially precluding long‑term habitability. For planets around G‑type stars, the mass‑loss rates are an order of magnitude lower but still often exceed the modern Earth’s escape rate.

To visualize the spatial distribution of hazardous radiation environments, the study constructs three‑dimensional maps of the average F_XUV,HZ within a 500 pc sphere centered on the Sun. By binning stars into cubic cells (≈5 pc per side) and averaging their HZ fluxes, the authors produce “hazard zone” maps that highlight regions of elevated XUV irradiation, typically associated with young stellar associations or dense clusters where active, fast‑rotating stars are abundant. These maps serve as a practical tool for prioritizing target selection in future exoplanet surveys, allowing observers to avoid regions where stellar high‑energy radiation would likely dominate atmospheric evolution.

In the discussion, the paper emphasizes several key points: (1) The majority of nearby main‑sequence stars are more XUV‑active than the contemporary Sun, especially low‑mass, magnetically active dwarfs; (2) The XUV‑to‑bolometric luminosity ratio is a strong function of effective temperature, underscoring the importance of stellar magnetic activity in shaping planetary environments; (3) Energy‑limited escape calculations suggest that many Earth‑like planets in the HZs of M dwarfs could experience rapid atmospheric erosion, challenging their habitability; (4) Large‑scale X‑ray surveys like eROSITA provide a powerful, uniform dataset for quantifying stellar high‑energy output across the Galaxy, complementing UV and optical observations; (5) Future work should refine EUV scaling relations with simultaneous X‑ray/UV observations, incorporate realistic atmospheric chemistry, and explore the role of stellar flares and coronal mass ejections, which can episodically boost XUV fluxes by orders of magnitude.

In conclusion, the authors demonstrate that SRG/eROSITA’s all‑sky X‑ray data, when combined with Gaia distances and bolometric corrections, enable a systematic, galaxy‑wide assessment of the high‑energy radiation environments that exoplanets experience. Their analysis reveals extensive “hazard zones” where stellar XUV radiation is sufficiently intense to drive significant atmospheric loss, especially around low‑mass stars. This work provides a valuable framework for evaluating planetary habitability and for guiding the selection of promising exoplanet targets for detailed atmospheric characterization with upcoming facilities such as JWST, ELT, and ARIEL.