The SRG/eROSITA All-Sky Survey. Detection of shock-heated gas beyond the halo boundary into the accretion region

Using SRG/eROSITA data from 680 galaxy clusters, this study detects shock-heated gas extending to 4.5 Mpc, identifies the virial radius as the transition point where cosmic filaments connect to clusters, and finds that observed gas distribution at large radii is more efficient than predicted by IllustrisTNG simulations.

The Gist

Imagine the universe not as an empty void, but as a giant, bustling city made of invisible threads and glowing fog. In this cosmic city, galaxy clusters are the massive downtown skyscrapers where hundreds of galaxies live together. But what about the space between these skyscrapers? That's where this paper takes us on a journey.

Here is the story of what astronomers found, explained simply:

1. The Mystery of the "Fog"

Galaxy clusters are mostly made of invisible dark matter (the "skeleton" of the city) and hot gas (the "fog" or "atmosphere"). We know a lot about the gas right inside the cluster, like the air inside a building. But what happens outside? Does the fog just stop at the door, or does it spill out into the neighborhood?

For a long time, this "outer fog" was a mystery. It's so faint and spread out that looking at just one cluster is like trying to see a single drop of mist in a hurricane. You can't see it.

2. The "Stacking" Trick: Building a Giant Cloud

To solve this, the astronomers didn't look at one cluster. They looked at 680 clusters at once.

Think of it like this: If you take a photo of a single firefly in the dark, it's hard to see. But if you take 680 photos of fireflies and stack them perfectly on top of each other, suddenly you see a giant, glowing cloud. This is called stacking.

Using data from the eROSITA telescope (a giant X-ray eye in space), the team stacked the images of these 680 clusters. The result? They finally saw the "fog" extending far beyond the cluster's edge, all the way out to 4.5 million light-years.

3. The "Splash" and the "Shock"

The paper explains two different types of gas behavior, using a great analogy:

• The Dark Matter (The Ghosts): Dark matter particles are like ghosts. They don't bump into each other. When they fall toward a cluster, they swing around and shoot back out, creating a "splash" at the edge. This is called the Splashback Radius.
• The Hot Gas (The Cars): The hot gas is different. It's like a stream of cars crashing into a wall. As the gas falls toward the cluster, it hits an invisible wall of other gas and shocks. This crash heats the gas up to millions of degrees, making it glow in X-rays. This is the Accretion Shock.

The astronomers found that this glowing shock wave extends much further out than we thought, filling the space between the clusters and the cosmic "highways" (filaments) that connect them.

4. The Cosmic Highway vs. The Empty Desert

The team also used supercomputer simulations (like a video game of the universe called IllustrisTNG) to understand the shape of this gas.

They discovered that the gas isn't spread out evenly like a smooth blanket. It's more like a river:

• Along the Filaments (The River): If you look toward the cosmic "highways" (filaments) connecting clusters, the gas is thick, hot, and flowing in like a river.
• Away from Filaments (The Desert): If you look toward the empty voids between highways, the gas is much thinner and stops abruptly at a shock wave.

The paper found that the transition point where the "river" of gas starts to dominate over the "desert" is right around the edge of the cluster ($r_{200m}$). This suggests that $r_{200m}$ is actually the spot where the cosmic highways physically connect to the cluster.

5. The "Feedback" Surprise: The Universe is Louder Than We Thought

Here is the biggest surprise. The astronomers compared their real observations with the computer simulations.

• The Simulation: The computer model predicted that the gas should be very concentrated near the center of the cluster, like a dense fog bank right around the skyscraper.
• The Reality: The real data showed the gas is much more spread out, extending far into the suburbs.

What does this mean? It means the "heating systems" inside galaxy clusters (called feedback) are working harder than the computer models thought. Imagine a giant fan in the center of the cluster blowing the gas outward, pushing it further away than the models predicted. The universe is more turbulent and energetic than our current simulations give it credit for.

6. The Big Picture

This paper is a milestone because it's the first time we've clearly seen the "atmosphere" of galaxy clusters spilling out into the deep universe.

• Before: We thought the gas stopped at the cluster's edge.
• Now: We know it flows out, connects to cosmic highways, and creates shock waves millions of light-years away.

It's like realizing that the smoke from a campfire doesn't just stay in the camp; it travels miles, connects to other fires, and shapes the weather of the whole forest. This discovery helps us understand how the universe builds its structure, one cluster at a time.

Technical Summary

Here is a detailed technical summary of the paper "Detection of shock-heated gas beyond the halo boundary into the accretion region" by Zhang et al. (2026).

1. Problem Statement

Galaxy clusters are the most massive collapsed structures in the Universe, positioned at the nodes of the cosmic web. While the inner regions of these clusters are well-studied, the outskirts (regions beyond $r_{500c}$ and extending beyond the virial radius $r_{200m}$) remain poorly explored. This region contains critical baryonic components:

• Shock-heated gas: Gas heated by accretion shocks as it falls into the cluster potential.
• Cosmic filaments: The large-scale structure connecting clusters to the rest of the universe.
• The "Missing Baryons": A significant fraction of the universe's baryons is expected to reside in the Warm-Hot Intergalactic Medium (WHIM) in these regions, but they are difficult to detect due to low surface brightness.

Previous studies have been limited by the rapid decline of X-ray surface brightness and the Sunyaev-Zeldovich (SZ) signal at large radii, often restricting analysis to within $r_{200c}$. Furthermore, distinguishing between the gas bound to the halo (one-halo term) and gas in infalling filaments or nearby halos (two-halo term) has been challenging.

2. Methodology

Observational Data

• Survey: The authors utilized data from the first two years of the SRG/eROSITA All-Sky Survey (eRASS).
• Sample: They selected 680 galaxy clusters from the western Galactic hemisphere.
  • Selection Criteria: X-ray selected from the eRASS1 catalog, luminosity $L_{0.5-2keV} > 2 \times 10^{43}$ erg s$^{-1}$, redshift $0.03 < z < 0.2$, and optical richness$\lambda > 20$.
  • Mass Range: $M_{500c} \approx 1.3 - 11.6 \times 10^{14} M_\odot$ (median $2.6 \times 10^{14} M_\odot$).
• Data Reduction:
  • Used the first four scans (eRASS:4) in the 0.2–2.3 keV band.
  • Applied rigorous masking for point sources, other clusters, and stray light (using a kernel convolution method).
  • Corrected for Galactic foreground and instrumental background.

Stacking Analysis

• Technique: A weighted stacking of surface brightness profiles was performed to boost the signal-to-noise ratio (S/N) in the low-surface-brightness outskirts.
• Weighting: Weights were inversely proportional to the projected sky solid angle ($w \propto D_A^2 r_{500c}^{-2}$) to minimize bias from massive or nearby objects.
• Uncertainty Estimation: Used a HEALPix-oriented bootstrap resampling method (500 samples) to derive covariance matrices.

Modeling Framework

The stacked surface brightness profile was modeled using a combination of three components:

1. One-Halo Term: Modeled using a Generalized NFW (gNFW) profile to describe the virialized intracluster medium (ICM). A cutoff was applied where the gas density drops to the mean cosmic baryon density ($\langle \rho_b \rangle$), representing the accretion shock boundary.
2. Two-Halo Term: Modeled theoretically to account for X-ray emission from correlated nearby halos and unvirialized gas in cosmic filaments. This utilized the halo bias, halo mass function (HMF), and luminosity-mass ($L_X-M$) scaling relations.
3. Constant Background: Instrumental background, Galactic foreground, and cosmic X-ray background.

Simulation Validation

• Simulation: Used the IllustrisTNG (TNG300-1) cosmological hydrodynamical simulation.
• Purpose: To validate the stacking/modeling framework and investigate the 3D thermodynamic structure of the gas.
• Analysis: Extracted gas properties (density, temperature, entropy, pressure) out to 5 $r_{200m}$. Crucially, they separated Line-of-Sight (LOS) directions into "into-filament" and "off-filament" (void) directions to study anisotropy.

3. Key Contributions & Results

A. Detection of Extended Gas

• The stacked X-ray surface brightness profile reveals a statistically significant signal (12$\sigma$) extending out to **$2 \times r_{200m}$** ($\sim 4.5$ Mpc).
• This represents the first robust detection of X-ray emission from galaxy clusters extending deep into the accretion region beyond the virial radius.

B. Density Profile and Transition Radius

• Best-fit Density: At $r_{200m}$, the gas number density is $n \approx 2.5 \times 10^{-5}$ cm$^{-3}$, corresponding to a baryon overdensity of $\sim 30$.
• One-Halo vs. Two-Halo: The transition between the one-halo (virialized) and two-halo (infalling/filament) dominated regimes occurs at approximately $r_{200m}$.
  • Inside $r_{200m}$: The profile is dominated by the gNFW term.
  • Outside $r_{200m}$: The two-halo term becomes significant, indicating the onset of gas connected to cosmic filaments and infalling sub-halos.
• Accretion Shock Estimate: By extrapolating the one-halo density profile to the mean cosmic baryon density, the authors estimate the accretion shock radius ($r_{shock}$) to be approximately $3 \times r_{200m}$.

C. Anisotropy and Filament Connection (Simulation Results)

• Analysis of TNG300-1 reveals strong directional dependence:
  • Off-filament directions: Hot gas is confined by the accretion shock, showing a sharp drop in temperature and density beyond $r_{200m}$.
  • Into-filament directions: Gas density and temperature remain elevated far beyond $r_{200m}$ due to continuous inflow.
• The radius where these profiles diverge aligns with $r_{200m}$, suggesting this radius marks the halo-filament connection point.

D. Comparison with Simulations (Feedback Efficiency)

• Discrepancy: The observed gas density profile is more extended (less centrally peaked) than predicted by the TNG300-1 simulation.
• Gas Fraction:
  • Without correction, the observed gas fraction at $r_{200m}$ exceeds the universal baryon fraction ($\Omega_b/\Omega_m \approx 0.158$), implying overestimation due to clumping.
  • After correcting for gas clumping (using a factor derived from simulations), the gas fraction at $r_{200m}$ drops to $\sim 80\%$ of the cosmic mean.
  • Even with this correction, the observed gas fraction is higher than in TNG300-1, suggesting that feedback processes in real clusters are more efficient at distributing gas to large radii than in the IllustrisTNG model.

4. Significance

1. Mapping the Cosmic Web: This study provides direct observational evidence of the transition from the virialized cluster core to the large-scale accretion region, effectively mapping the "bridge" between clusters and cosmic filaments.
2. Defining Halo Boundaries: It supports the use of $r_{200m}$ as a characteristic radius marking the connection to filaments, distinct from the splashback radius (which is a dark matter feature). The accretion shock is identified as the true physical boundary for hot gas in off-filament directions.
3. Baryon Census: The detection of shock-heated gas in the outskirts helps account for the "missing baryons" in the local universe, although clumping effects must be carefully accounted for to avoid overestimating gas mass.
4. Feedback Constraints: The mismatch between observations and TNG simulations regarding gas distribution suggests that current feedback models (AGN and stellar feedback) may be underestimating the efficiency of gas redistribution to the cluster outskirts.
5. Future Prospects: The work establishes a baseline for future high-sensitivity missions like NewAthena, which will be required to resolve individual systems and probe the detailed physics of the accretion shock and filamentary gas.

Conclusion

Zhang et al. successfully detected and modeled the shock-heated gas extending from galaxy cluster outskirts into the accretion region using eROSITA data. By combining a robust stacking analysis with IllustrisTNG simulations, they identified $r_{200m}$ as the transition radius to filamentary gas and found that real clusters exhibit more extended gas distributions than current simulations predict, pointing to the need for stronger feedback mechanisms in cosmological models.