The eROSITA view of the Abell 3391/95 field: Cluster Outskirts and Filaments

Aims. We aim to characterize the gas properties in the cluster outskirts ($R_{500}R_{200}$) of the A3391/95 system and to compare them to predictions. Methods. We performed X-ray image and spectral analyses using the eROSITA PV data to assess the gas morphology and properties in the outskirts and the filaments in the directions of the previously detected Northern and Southern Filament of the A3391/95 system. We took particular care of the foreground. Results. In the filament-facing outskirts of A3391 and the Northern Clump, we find higher temperatures than typical cluster outskirts profiles, with a significance of $1.6-2.8\sigma$, suggesting heating due to their connections with the filaments. We confirm SB excess in the profiles of the Northern, Eastern, and Southern Filaments. From spectral analysis, we detect hot gas of ~1 keV for the Northern and Southern Filaments. The filament metallicities are below 10% solar metallicity and the $n_e$ range between 2.6 and $6.3\times10^{-5}~\mathrm{cm^{-3}}$. The characteristic properties of the Little Southern Clump (~1.5$R_{200}$ from A3395S in the Southern Filament) suggest that it is a small galaxy group. Excluding the LSC from the analysis of the Southern Filament decreases the gas density by 30%. This shows the importance of taking into account any clumps to avoid overestimation of the gas measurement in the outskirts and filament regions. Conclusions. The $n_e$ of the filaments are consistent with the WHIM properties as predicted by cosmological simulations, but the temperatures are higher, close to the upper WHIM temperature limit. As both filaments are short and located in a denser environment, stronger gravitational heating may be responsible for this temperature enhancement. The metallicities are low, but still within the expected range from the simulations.


Introduction
The missing baryon problem states that about 30% − 40% of the total baryons in the Universe are still unaccounted for (e.g., Shull et al. 2012).According to the large-scale cosmological simulations (e.g., Cen & Ostriker 1999;Davé et al. 2001), they should be found in the warm-hot intergalactic medium (WHIM) located in the so-called cosmic web filaments (Bond et al. 1996).
These filaments connect galaxy clusters and play an important role in transferring matter onto the clusters, thus supporting their growth (West et al. 1995;Tanaka et al. 2007;Bond et al. 2010).
Due to the faint nature of the WHIM emission (the electron densities and temperature ranging between n e ≈ 10 −6 − 10 −4 cm −3 and T = 10 5 − 10 7 K (Nicastro et al. 2017), respectively) detecting the WHIM is challenging.However, with the X-ray emissivity being proportional to the square of the density, Fig. 1: eROSITA PIBsubtracted, exposure-and Galactic-absorption-corrected image in the 0.3 − 2.0 keV band.The image has been adaptively smoothed with S/N set to 12.The point sources have been removed and refilled with their surrounding background values.The defined regions are used for spectral analyses (see Sec. 2.2).The colorbar is in the units of counts per second (as well as the following X-ray images).
The binary galaxy cluster system Abell 3391 and Abell 3395 (A3391/95) is also one of the systems where a search for the WHIM has been carried out, and the system has therefore been observed extensively using various instruments of different wavelengths, including several X-ray telescopes;for example by ROSAT and ASCA and more recently by XMM-Newton, Chandra, and Suzaku.Earlier studies showed that the gas between the A3391 and the double-peaked A3395S/N (hereafter referred to as the bridge) is consistent with a filament (Tittley & Henriksen 2001).While latter studies concluded that the gas properties in these regions are more typical of tidally stripped intracluster medium (ICM) gas from A3391 and A3395, a sign of an early cluster merger phase (Sugawara et al. 2017;Alvarez et al. 2018).A study of the A3391/95 system was also carried out using Planck data, and a significant excess of signal near the A3395S/N system in the thermal Sunyaev-Zel'dovich (tSZ) residual map was reported (Planck Collaboration et al. 2013).However, the individual components of this cluster cannot be resolved by the data due to the modest resolution of Planck.
The extended ROentgen Survey with an Imaging Telescope Array (eROSITA, Predehl et al. 2021) is the newly launched German X-ray telescope on board the Spectrum-Roentgen-Gamma (SRG) mission.It has a wide field-of-view (FoV) of 1.03 degrees in diameter.In the soft energy band of 0.5 − 2.0 keV, the effective area of the seven combined eROSITA telescope modules (TMs) is slightly higher than the combined XMM-Newton cameras (MOS1+MOS2+pn) (Predehl et al. 2021).Additionally, eROSITA operates in the scan and survey modes, which allows us to easily reach the virial radii of clusters and to constrain foreand background emission.The FoV-averaged PSF (∼ 26 ′′ ) is significantly better than that of the Suzaku satellite, allowing us to remove unrelated background sources to much lower flux levels.The energy resolution at soft energies is better than for any other operating X-ray satellite CCD, allowing us to better characterize soft band emission lines.And the short focal length of 1.6m minimizes the particle-induced background (PIB) levels.These properties make eROSITA an excellent instrument for studying the outskirts of galaxy clusters and finding large-scale structures.
In October 2019, the A3391/95 system was observed as one of the eROSITA performance verification (PV) targets.The field was observed four times during this phase, with three raster-scan observations and one pointing observation.The combination of all four observations covers an area of about 15 square degrees, which gives us the deepest large-scale X-ray view of the system to date.Using these observations, eROSITA revealed indications of the presence of warm gas in the bridge region, on top of the tidally stripped hot ICM gas (Reiprich et al. 2021).This finding is also seen in the analogous A3391/95 cluster system from Magneticum simulations (Biffi et al. 2022).Furthermore, a continuous soft emission extending from the north to the south was dis- covered.This ∼15 Mpc filamentary structure seems to connect at least five galaxy groups and clusters.Those are A3391, A3395S, A3395N, the Northern Clump1 (NC; MCXC J0621.7-5242cluster or MS 0620.6-5239cluster, Tittley & Henriksen 2001;Piffaretti et al. 2011;Veronica et al. 2022), the Little Southern Clump (LSC), and the MCXC J0631.3-5610cluster.The LSC is embedded in the Southern Filament and has fully entered the 2R 200 of A3395, while the center of the MCXC J0631.3-5610cluster is just outside the FoV of these PV observations.
In this work, we utilize the eROSITA A3391/95 PV data to investigate the gas properties of the outskirts and the detected filaments.These filaments are the Northern and Eastern Filaments extending north and northeastward from the A3391 cluster, and the Southern Filament in the south of the A3395S/N cluster.The present paper is structured as follows: in Sect.2, we describe our observations, the data reduction steps, and our analysis strategy.In Sect.3, we present the results of our imaging analysis (Sect.3.1) and those of our spectral analysis (Sect.3.2).In Sect.4, we discuss these results.Finally, in Sect.5, we summarize our findings and conclude the results.

Data reduction and analysis
All eROSITA A3391/95 PV observations are listed in Table 1.We used eROSITA data processing in configuration c001.The data reduction steps of all 16 eROSITA data sets were performed using the extended Science Analysis Software (eSASS, Brunner et al. 2022) version eSASSusers_201009.The data reduction steps and the image-correction steps, were done following the procedure described in Sect.2.1 and 3.3 of Reiprich et al. (2021), which include the particle-induced background (PIB) subtraction, exposure correction (including vignetting), and Galactic absorption correction across the FoV.

Imaging analysis
To focus on the soft emission in the outskirts and from the filaments, we used the energy band of 0.3 − 2.0 keV for the imaging analysis.The lower energy limit used for the telescope modules (TMs) with an on-chip filter (TM1, 2, 3, 4, 6; the combination of these TMs is referred to as TM8) was set to 0.3 keV, while for the TMs without an on-chip filter (TM5 and 7; the combination of these TMs is referred to as TM9) due to the optical light leak contamination (Predehl et al. 2021), the lower energy limit was set to 0.5 keV.The count rates of the final image (all observations combined, and fully corrected) correspond to an effective area given by one TM with an on-chip filter in the energy band 0.3 − 2.0 keV.
The Galactic absorption across the FoV was corrected by calculating the ratio of the expected eROSITA count rates with the total hydrogen column density (N H ) at a given position to the count rate of the median N H across the FoV.The count rates of different N H values were estimated assuming a typical X-ray fore-and background model, that is one unabsorbed with one absorbed from diffuse thermal components and one power law component.These components represent the Local Hot Bubble (LHB), the Milky Way halo (MWH), and unresolved sources, respectively.The computed N H values are from the A3391/95 total N H map (for more details, see Sect.2016).A ratio map relative to the median N H was generated for each type of eROSITA TMs.We divided the exposure map by their corresponding ratio map for the Galactic absorption correction.
The point source catalog was generated using the method described in Pacaud et al. (2006) and Ramos-Ceja et al. (2019), that is, the point sources were detected on a wavelet-filtered image by the Source Extractor software (Bertin & Arnouts 1996).Using the obtained catalog, we excluded them from our analyses.The so-called ghosting procedure was also performed to generate the surface brightness image.Through this procedure, the areas from which the point sources were excised were refilled with surrounding background photons.Any ghosted images generated in this work were only used for visualization purposes.To avoid any biases introduced by using the artificially refilled pixels in the analysis, we excluded these pixels from any calculations and the missing area was properly taken into account.The final PIB-subtracted, point-sources-subtracted and refilled, exposure-corrected, and Galactic-absorption-corrected image is shown in Fig. 1.The image has been adaptively smoothed with a signal-to-noise ratio (S/N) set to 12 in order to enhance lowsurface-brightness emission and large-scale features.
The X-ray surface brightness profiles were calculated using the final eROSITA products in the 0.3 − 2.0 keV energy band, including the photon image, the fully corrected exposure map, and the PIB map, as well as the point source catalog.The large eROSITA FoV and the SRG scan mode in which most of the observations were taken, allow us to model the cosmic X-ray background (CXB) from regions within the FoV.To account for the variation of the CXB, we placed ten boxes across the field.To ensure minimum emission from the clusters and the filaments, which could cause overestimation in the CXB level (and therefore underestimation of the surface brightness values), they were placed beyond the R 2002 of each cluster and the defined area of the filaments.The configuration of these CXB boxes can be found in Fig. F.1. of Reiprich et al. (2021).The average of the surface brightness values of these ten boxes was taken to be the CXB level in the field, SB CXB , while the root-mean-square deviation of these values was taken as the standard deviation, σ CXB , which is larger than the statistical uncertainty.We report that the CXB surface brightness value in the A3391/95 system is of the typical value.The CXB value we obtained is of the same order of magnitude as, for example, to the CXB surface brightness value calculated for the A3158 surface brightness profile using the eROSITA calibration data (Whelan et al. 2022), which is not affected by additional foreground components, such as the eROSITA Bubble or other extended Galactic sources.
In their review, Reiprich et al. (2013) define the outskirts of galaxy clusters as the regions found within the range of R 500 < r < 3R 200 .In the present work, we also define the lower boundary of the outskirts region to be R 500 , but as the 3R 200 of the A3391/95 clusters extends to the entire FoV of the eROSITA PV data, we chose an arbitrary upper boundary of R 200 .The filament was defined as the inter-cluster emission beyond the R 200 of a cluster.The used cluster radii, as well as their corresponding emission peak coordinates and their references, are listed in Table 2.The R 500 of A3391 and A3395N/S were calculated from ROSAT and ASCA observations (Reiprich & Böhringer 2002); the R 500 of the NC was estimated from a dedicated XMM-Newton observation (Veronica et al. 2022); and the R 500 of the LSC is a product of the present study (eROSITA PV data; see Sec. 3.2.1).

Spectral analysis
We utilized the eSASS task srctool to extract all eROSITA spectra, the Ancillary Response Files (ARFs), and the Response Matrix Files (RMFs) for the source and background regions.The spectral fitting was performed with XSPEC (Arnaud 1996) version: 12.10.1f.
For the spectral analysis, we used the third eROSITA scan observation (ObsID: 300016), where all seven eROSITA TMs are available.The configuration of the regions used in the spectral analysis is shown in Fig. 1.Any source regions related to the Northern (Eastern) Filament are indicated with dashed white (green) color.For instance, the A3391 outskirts region (R 500 − R 200 ) in the direction of the Northern (Eastern) Filament is plotted as a white (green) sector, labeled A3391-N (A3391-E).The Northern (Eastern) Filament source regions are represented by two white (green) boxes, labeled NF BOX1+2 (EF BOX1+2).In the south, the source regions are depicted with magenta boxes, labeled A3395S and SF BOX1+2 for the outskirts and filament regions, respectively.The cyan circle marks the R 500 of the LSC.For each filament, we placed two identical boxes in line with one another outside the R 200 of the corresponding parent clusters.For the analyses of the Northern and the Eastern Filament, we took the average redshift between the A3391 cluster and the NC, z N/E = 0.0533.While for the Southern Filament, we used the average redshift between the A3395N, A3395S, and the MCXC J0631.3-5610clusters, z S = 0.0525.The height (h) and width (2r) of each box are 14 ′ (0.90 Mpc) and 50 ′ (3.23 Mpc) for the Northern Filament, 25 ′ (1.62 Mpc) and 33 ′ (2.13 Mpc) for the Eastern Filament, and 20.83 ′ (1.33 Mpc) and 41.67 ′ (2.66 Mpc) for the Southern Filament, respectively.
The X-ray spectral fitting model, including the CXB and the source emission, is described in the following equation: where the first term represents the CXB components.The absorption parametrized by the hydrogen column density along the line of sight, N H , is represented by TBabs (Wilms et al. 2000).
The N H values used in this work are taken from the A3391/95 total N H map (see Sect. 2.1).The X-ray foreground emission from the LHB and the MWH are described by apec 1 and apec 2 , respectively.Lastly, the cosmic X-ray background from the unresolved sources (e.g., Luo et al. 2017) is represented by the power law component, powerlaw.
The second term, TBabs×apec 3 , represents the absorbed cluster or filament emission.The apec parameters include temperature k B T , metallicity Z, redshift z, and normalization norm.The parameter norm is defined as where D A is the angular diameter distance to the source in centimeters (cm), while n e and n H are the electron and hydrogen densities in cm −3 , respectively.
To take into account the N H variations in the FoV, two different CXB regions were used.These were defined based on how closely the average N H (see Sect. (5.63 +0.16 −0.15 ) × 10 −7 N/E from CXB used for the Northern or Eastern Filament analysis, S from CXB used for the Southern Filament analysis, † [10 22 atoms cm −2 ], * [cm −5 /arcmin 2 ], ‡ [photons/keV/cm 2 /s/arcmin 2 at 1 keV] Fig. 2: eROSITA spectrum of the northern outskirts sector of A3391 (A3391-N).The spectra and the corresponding response files of TM3 and TM4 are merged for better visualization.The black points are the spectral data, while the red, green, and blue lines represent the total model, the cluster and CXB model, and the instrumental background model, respectively.
The third term of Eq. 1, PIB, represents the instrumental background model.We utilized the results of the eROSITA EDR Filter Wheel Closed (FWC) 3 data analysis to model the instrumental background of TM8 and the modified version for TM9, where the parameters were constrained including the lowerenergy FWC spectra.The components of the instrumental background consist of a combination of a power law and an exponential cut-off to model the signal above ∼1 keV and two power laws to model the background increase due to the detector noise at low energy.Additionally, 14 Gaussian lines are included to model the fluorescence lines caused by the interaction between the cosmic 3 https://erosita.mpe.mpg.de/edr/eROSITAObservations/EDRFWC/particles with the detector components.Among these, the brightest lines are Al-K α at 1.486 keV and Fe-K α at 6.391 keV, and the weaker lines are Ti-K α at 4.504 keV, Co-K α at 6.915 keV, Ni-K α at 7.461, Cu-K α at 8.027 keV, and Zn-K α at 8.615 keV.Further details of the eROSITA in-flight background are discussed in Freyberg et al. (2022).
We performed the fitting in the energy band of 0.3 − 9.0 keV for the TMs with an on-chip filter.Below 0.3 keV, the detector noise becomes stronger, while no strong source signal is expected above 9.0 keV given eROSITA's effective area.We find that the light leak contamination in the A3391/95 eROSITA PV observations is lower than in the other Cal/PV observations.Based on the light-leak diagnostic of this field, we can use a low energy limit of 0.5 keV for TM5 and TM7.
The spectral fitting procedure started with fitting the CXB spectra of all TMs in order to obtain best-fit values and the errors on the normalizations of the CXB components.The information from this first fitting in the respective CXB regions, including the CXB components and the resulting best-fit values are listed in Table 3. Subsequently, the source and CXB spectra of all TMs were fitted simultaneously.In this second fitting, we freed the apec 3 parameters, k B T , Z, and norm, of the source spectra whenever the statistics allowed.Otherwise, we fixed Z to 0.3 of Asplund et al. (2009) for the cluster outskirts (e.g., Reiprich et al. 2013;Urban et al. 2017), and 0.1 and 0.2 for the filaments (e.g., Tanimura et al. 2020;Biffi et al. 2022).The normalizations of the instrumental background components were also freed.The normalizations of the CXB components in this second fitting were also initially freed.However, we found one case where the fit was seen to be driven by one of the CXB components, that is the Eastern Filament fit with metallicity fixed to 0.1 solar.In this case, the obtained temperature for the Eastern Filament is 0.23 keV, which is close to the MWH temperature of 0.25 keV.This temperature drop is accompanied by an increase in source normalization, while the MWH normalization decreases.Translating the Eastern Filament normalization value in this case to the gas overdensity (δ b ), we find δ b ≈ 356 (3.5 times higher than the value obtained when the CXB normalizations are restricted; see Table 4).This value is not physical for cluster outskirts nor filament regions, but more common for the hot and denser ICM gas.The MWH normalization in this case drops by 18%.Therefore, to avoid any source parameters being driven by the CXB components, we restricted the CXB normalizations of the source spectra to the minimum and maximum values obtained from the CXB spectra from the first fit.We note that the resulting gas properties of the outskirts and filaments between freeing and restricting the CXB normalizations are consistent within the 1σ uncertainties, except for the Eastern Filament case described above.In this work, we present only the results from the restricted source CXB normalization fit.The C-statistic (Cash 1979) was adopted.An example of the eROSITA spectrum and its fitted model are shown in Fig. 2.

Configuration
We present the configuration of the source region used for the surface brightness analysis around A3391 and A3395S clusters in the left and right panels of Fig. 3, respectively.The corresponding PIB-subtracted surface brightness profiles in the 0.3 − 2.0 keV band are shown in Fig. 4a and b where SB r>R 200 and σ r>R 200 are the surface brightness value of the entire r > R 200 region and its statistical standard deviation, and SB CXB and σ CXB are the sky background surface brightness and its standard deviation, which includes a systematic uncertainty (see Sect. 2.1).
In the directions of the Northern and Eastern Filament, we calculated A3391 sectorized surface brightness profiles out to R 200 , namely A3391-N and A3391-E (cyan and red sectors in Fig. 3, left).The surface brightness values of the Northern (Eastern) Filament were then calculated from the four cyan (red) boxes of the same size, which were placed outside the R 200 of the A3391 cluster and the NC.Each Northern Filament box has a size of 7.25 ′ × 50 ′ and is angled at 234.71 • from the R 200 of the A3391 in the NC direction.For the Eastern Filament, each box has the size of 12.5 ′ × 33.33 ′ and is angled at 300 • from the R 200 of the A3391.In the Northern Filament direction, we continued with the surface brightness profile of NC, which was calculated from annuli around the center of the cluster.The data points related to the Northern Filament (filN) and Eastern Filament (filE) are shown respectively as blue upward triangles and red squares in Fig. 4a.To compare with the surface brightness profiles from the filament regions, we also computed surface brightness profiles in the western (green) and eastern (orange) directions out to 57.33 ′ (∼ 2R 200 ) from the cluster center.The western and eastern sectors have position angles of 316 • to 378 • and 149 • to 209 • , respectively.The two apparent extended sources at the edge of the eastern sector are background sources located at z ∼ 0.1, and therefore we excluded them from our analysis.Due to the influence of the A3395S/N cluster, we did not calculate any surface brightness profiles in the southern direction.The black data points in Fig. 4a are the average surface brightness profile out to the R 200 of A3391 (excluding the southern area).By comparing the sectorized profiles with this average profile, we checked for any deviations (enhancement or depression) from the spherically symmetric surface brightness profile.
As the A3395S/N cluster comprises two components, a different approach was taken to place the surface brightness regions compared to that used to place those around the A3391 cluster, such that we avoided treating the cluster as a single entity.To minimize the amount of emission from A3395N and to get the coverage of the Southern Filament, we placed five sectors covering the lower half of A3395S (Fig. 3, right).These sectors were centered at the A3395S (see Table 2) and have an opening angle of 36 • each, with the position angle from Sector1 to Sector5 spanning from 215 • to 395 • .From each sector, we calculated the surface brightness profile out to 72.17 ′ (∼ 2R 200 ) from the center.The black data points in Fig. 4b are the average surface brightness profile of all sectors.

Inner region (r ≤ R 500 )
As seen in Fig. 4a (top plot), within 1 − 11 ′ , the surface brightness values of the A3391-N (blue downward-pointing triangles) are significantly lower than the other directions.In this radial range, compared with the full annuli profile (black data points), the A3391-N surface brightness values are lower, with relative differences ranging 8% − 36%.The lowest surface brightness value in this radial range is located at 6.3 ′ with a significance of  11.5σ to the full annulus value.This shows that A3391 has an elliptical morphology, where the minor axis is in the north-south direction while the major axis is in the east-west direction.
Similarly, for A3395S, we observe significantly higher surface brightness in Sector5 (Fig. 4b, cyan upward-pointing triangles).Within 2.7 − 21.6 ′ , the relative difference values with respect to the average surface brightness profile (black data points) range between 16% and 93%, with the peak located at 7.1 ′ with a 13.7σ significance.This bump shows the apparent cluster elongation in this direction.From the center to ∼ 15 ′ , the profiles of Sector1 (blue squares) and Sector2 (red diamonds) drop rapidly, and then flatten outward, denoting the pronounced shorter edge of the cluster in the southern direction.
Within R 500 < r < R 200 , the enhancement in the western profile of A3391 (Fig. 4a, green diamonds) is apparent.This is not unexpected given that the elongation of the cluster is in this direction.The average enhancement in this radial range with respect to the annuli profile (black data points) is 14 ± 2%.For the other sectors, the surface brightness profiles are on average lower than the annuli with relative differences of 4 ± 2% (A3391-N), 3 ± 2% (A3391-E), and 12 ± 2% (E), meaning that the A3391-N and A3391-E are still brighter than the E sector.While ICM emission is still expected in this regime, other sources of emission, for example from the filaments, may contribute to this outskirts excess.
For the A3395 profiles (Fig. 4b) in the R 500 < r < R 200 range, an enhancement is observed for Sector5 (cyan) with 15 ± 2% average excess with respect to the average profile in the same regime (black points).We attribute this enhancement to the emission from the cluster in this direction.Decrements with an average of 5 ± 3% and 9 ± 2% are observed for Sector3 (green) and Sector4 (orange).Similar to the case for A3391, insignificant decrements with respect to the average profile of 4 ± 2% and 2 ± 3% are seen for Sector1 (blue) and Sector2 (red).

Filament region (r > R 200 )
Beyond R 200 , on average we observe higher surface brightness values in the filN and filE profiles for A3391 (blue upwardpointing triangles and red squares) than in the western (green) and eastern profiles (orange), indicating an enhancement in the intercluster regime.The overall significance values (S/N all ; Eq. 3) for these r > R 200 regions show that the Northern and Eastern Filament are > 3.2σ away from the CXB level, while the Western and Eastern Sectors are 1.5 − 1.7σ away.We note that the first r > R 200 bins of all sectors are greater than 3.2σ.In general, the enhancement near R 200 is expected to originate from the ICM.As can already be observed from the image (e.g., Fig. 3, left), there is an apparent elongation of the cluster to the west.Moreover, the surface brightness elevation of the eastern profile could result from additional residual emission from the excluded background bright galaxies, seen at ∼ 28 ′ and 50 ′ .Removing these first bins in the r > R 200 regime results in S/N all = 3.5σ (filN), 2.5σ (filE), 0.5σ (W), and 1.2σ (E).This shows that there is still significant excess emission in the Northern and Eastern Filament regions in comparison to the other directions.
For A3395S's r > R 200 , a striking surface brightness peak is seen at around 54 ′ (∼ 1.5R 200 ) of Sector1.This peak is attributed to the LSC.For Sector1 and Sector2, the S/N all values in this regime are 2.7 and 2.3σ, respectively.Combining Sector1 and Sector2, we obtain an excess with an overall significance of 2.5σ.Meanwhile, for Sector3-5, the S/N all values are consistent with the CXB level, such that the S/N all values lie around zero (Fig. 4b, bottom plot).

The Little Southern Clump
The LSC is an extended source that resides in the Southern Filament and has completely entered the 2R 200 of the A3395 cluster, as shown in Fig. 1.We show a zoomed-in image around the LSC in Fig. 5.We overlaid the image with the A3391/95 galaxy density map (white contour; also see Sect.2016)).The white circles mark the galaxies retrieved from the NED4 with known redshifts and further selected with 0.0424 ≤ z ≤ 0.0636.We labeled their redshifts on top of each circle.The cyan circle indicates the calculated R 500 of the LSC and the white dashed lines mark the R 200 and 2R 200 of the A3395S cluster.
To characterize the LSC and to exclude it during the analysis of the Southern Filament, we need to first determine some characteristic properties of this object.Through spectroscopy, we acquired the k B T 500 , which was then used to infer the M 500 and consequently, the R 500 .The k B T 500 was calculated through an iterative procedure, such that we extracted and fit spectra from an annulus centered at the WISEA J063004.80-551751.5 galaxy at (α, δ) = (6 : 30 : 04.8, −55 : 17 : 51.504).The in-ner and outer radii were varied until they correspond roughly to the 0.2 − 0.5R 500 of the source region, where the R 500 values are listed in Table 2. Based on this, we obtain k B T 500 = 0.94 +0.06 −0.07 keV and Z 500 = 0.10 +0.05 −0.04 Z ⊙ from an annulus with inner and outer radii of 1.3 ′ and 3.2 ′ , respectively.We input this value into the mass−temperature (M − T ) scaling relation by Lovisari et al. (2015), where a = 1.65 ± 0.07, b = 0.19 ± 0.02, C 1 = 5 × 10 13 h −1 70 M ⊙ , and C 2 = 2.0 keV.Assuming spherical symmetry and taking 500 times the critical density of the Universe at the WISEA J063004.80-551751.5 redshift as ρ 500 (z = 0.0562) = 4.51×10 −27 g cm −3 , we obtain M 500 = (2.33 +0.31  −0.34 ) × 10 13 M ⊙ and R 500 = (6.44 +0.29  −0.32 ) ′ ≈ 437.23 +19.38 −21.47 kpc.The calculated mass is lower than the one calculated using the eFEDS scaling relation (Chiu et al. 2022), (4.30 ± 2.8) × 10 13 M ⊙ , but is nevertheless in agreement within their 1σ uncertainties.As the studied group is closer to the cluster population considered in Lovisari et al. (2015), which are average-and low-mass systems at low-z, we use the M 500 = (2.33 +0.31  −0.34 ) × 10 13 M ⊙ as the default.Based on the results, we classify the LSC as a small group of galaxies.This is further supported by the DECam/DES galaxy density map (Fig. 5, white contours), where galaxy overdensity was observed within the R 500 and extending towards the southwest.We also notice an overdensity in the northeast direction that coincides with six known galaxies.As has been discovered in simulations, galaxy groups that are within the vicinity of a cluster can accrete some field galaxies and grow through this process (Vijayaraghavan & Ricker 2013;Kuchner et al. 2022).However, these galaxies, together with the LSC, may also simply be on their way to being accreted towards the A3395 cluster with the filament.Deeper pointed observations, for example by XMM-Newton and Chandra, will allow us to study the LSC and its surroundings in more detail.With the newly determined R 500 of the LSC, we performed two analyses of the Southern Filament, namely including and omitting the emission from the LSC within this radius (Sect.3.2.3).

The outskirts (R
We list the results of the eROSITA spectral analysis in Table 4. Starting from the northern outskirts of A3391 (A3391-N; Fig. 1, white dashed lines sector), we report an upper limit of the metallicity of 0.31Z ⊙ .Meanwhile, due to a lack of statistics, we find that it is necessary to fix the metallicity to 0.3 during the fit of the eastern sector (A3391-E; Fig. 1, green sector).This is a typical metallicity value adopted for cluster outskirts (e.g., Reiprich et al. 2013;Urban et al. 2017).Overall, the normalizations and temperatures obtained from A3391-N and A3391-E are in good agreement with each other.The temperature values range between 2.2 and 4.2 keV.
The temperature of the outskirts of the NC (Fig. 1, white dashed lines sector) is k B T = 0.96 +0.18  −0.19 keV, which is lower than the temperature measured for the outskirts of the NC by Veronica et al. (2022) based on data from XMM-Newton, of namely k B T = 1.71 +0.66 −0.37 keV.Although there is overlap between the definitions of the source regions, they are not identical: the eROSITA region is placed in the south of the cluster at 10.62 − 16.34 ′ , while the XMM-Newton source region is placed at 10.0 − 12.5 ′ .This temperature difference might also be the product of some other systematic differences, such  2013) using k B T 0.2−0.5R500 as mean temperature of each cluster (see Table 5).The black vertical solid (dotted) lines mark the R 500 (R 200 ) of the parent clusters.
Article number, page 9 of 17 as background treatment, point source selection, and/or instrumental calibration.Furthermore, the yielded metallicity is Z = 0.05 +0.06 −0.04 Z ⊙ , whereas for XMM-Newton the value has to be frozen at 0.3Z ⊙ for statistical reasons.For the A3395 outskirts region (Fig. 1, magenta box), the temperature and metallicity values are 1.68 +0.31  −0.15 keV and 0.72 +0.60 −0.33 Z ⊙ , respectively.Furthermore, we fit eROSITA spectra extracted over 0.2 − 0.5R 500 of the A3391, A3395S/N, and NC clusters.The fitted eROSITA temperatures constrained using the R 500 (k B T 0.2−0.5R500 ) from the literature and the eROSITA newly defined radii (k B T 0.2−0.5R500,eRO ), as well as other derived quantities (M 500,eRO and R 500,eRO are presented in Table 5.We report that the eROSITA temperatures for A3391 and A3395N are lower than ASCA values listed in Reiprich & Böhringer (2002), with significances of 2.4σ and 4.9σ, respectively.The eROSITA best-fit temperature for A3395S is higher than that reported for A3391.We argue that since A3395S and A3395N are interacting, this hotter gas could be the result of some merging processes, for instance, shock heating or adiabatic compression.For NC, comparing with the XMM-Newton temperature reported in Veronica et al. (2022), the eROSITA result is 25% lower with 2.7σ significance.We note that this could either originate from some differences in the analyses, such as the source catalog used, or different background treatments or a calibration issue between the two instruments.
To obtain the eROSITA temperatures (k B T eRO ) that were calculated within 0.2−0.5R500,eRO , we repeated the M−T scaling relation procedure (see Sect. 3.2.1),starting with the extraction of source spectra from an annulus region.We iterated the steps until the inner and outer radii of the source region corresponded to 0.2 and 0.5 of the calculated R 500,eRO .We acquire lower k B T eRO values, which consequently mean smaller eROSITA cluster radii (R 500,eRO ).If we consider these newly calculated eROSITA radii, for example R 500 = 1.25R 500,eRO for A3391 cluster, and take into account our definition of inter-cluster emission, that is, that found at r > R 200 , we find the filaments to be longer.However, for consistency, the characteristic radii used to present the results and discussion throughout this work (and other eROSITA A3391/95 papers) are those listed in Table 2, unless otherwise stated.These eROSITA cluster properties of the A3391/95 clusters are presented in Table 5.
We employed the eROSITA k B T values constrained from within apertures of 0.2 − 0.5R 500 to compute the universal temperature profile from Reiprich et al. (2013).The profile describes the cluster temperatures within 0.3R 200 < r < 1.15R 200 , such that k B T = (1.19− 0.84r/R 200 )⟨k B T ⟩.We took ⟨k B T ⟩ to be the eROSITA temperatures (Table 5, second row).We display the temperature profiles of the outskirts and filament regions related to the Northern, Eastern, and Southern Filaments in the top, middle, and bottom panels of Fig. 6, respectively (data points).In each plot, the universal temperature profile from the corresponding cluster is plotted as the pink-shaded area.As shown, the acquired temperatures for the outskirts of A3391-N, A3391-E, and NC (Fig. 6, top and middle plots) are higher than the predicted values with a significance of 2.4σ, 2.8σ, and 1.6σ, respectively.

The filaments
We used the obtained apec normalizations to derive the electron densities, n e , of the filament regions.We assumed a simple geometry for the filament, that is a cylinder with its axis in the plane of the sky.The n e was calculated by rewriting Eq. 2: where, for instance, D A (z = 0.0533) = 222.22Mpc.The radius and height of the cylinder are represented by r and h.The values of r and h of the filament boxes were introduced in Sect.2.2.The hydrogen density is taken to be n H ≈ n e /1.17.Given that the radii of our filament cylinders ranges from 1.07 to 1.62 Mpc, and that the gas in filaments is isothermal up to ∼ 1.5 Mpc from their cores, as reported by Galárraga-Espinosa et al. ( 2021), we deem our assumptions for calculating n e to be reasonable.
Using n e , we derived the filament gas density, ρ gas , that is where µ e ≈ 1.14 is the mean molecular weight per electron, and m p = 1.6726×10 −24 g is the proton mass.Then, we calculated the gas overdensity To account for the systematic errors introduced by the assumed geometrical shape, we also calculated the n e and δ b using two alternative shape assumptions.The first shape is an elliptical cylinder, where the major axis and the height are assumed to be the same as the radius and the height of the previously assumed cylinder, while the minor axis is taken to be half of the radius.The second assumption is a cuboid shape, where we take the width of the base and the height to be the diameter and the height of the cylinder.The density contained in the elliptical cylinder shape is 41% higher than that of the cylinder, while the density of the cuboid shape is 11% lower.Using either of the geometrical shape assumptions does not change the classification of the A3391/95 filaments based on the overdensity criteria stated in Shull et al. (2012), namely that they fall between the WHIM criteria.In Table 4, we therefore only list the n e and δ b values calculated whilst assuming the common cylindrical shape.
We compared our results to the WHIM properties predicted by cosmological hydrodynamic simulations (e.g., Davé et al. 2001;Martizzi et al. 2019.For observational constraints of WHIM properties, see, e.g., Shull et al. 2012;Nicastro et al. 2018;Kovács et al. 2019).Here, we adopted temperature cut of T = 10 5 − 10 7 K (k B T = 0.009 − 0.862 keV) and the gas overdensity δ b = 0.001 − 316.22 for the WHIM (Shull et al. 2012).Whereas, gas with k B T > 0.862 keV and any number density n H was considered as the hot medium (HM), which is the shock-heated gas in and near the most massive dark matter haloes (Martizzi et al. 2019).
We present the results of the spectral analysis of the different filaments in Table 4. Firstly, for the Northern Filament, we obtain a best-fit temperature of ∼ 1 keV and an upper metallicity limit of 0.01Z ⊙ .To investigate whether this low metallicity is caused by a bias of some kind, for example, by multitemperature components, we performed a two-temperature (2T) component fit.We obtain a cool and a warm component of 0.1 and 1.3 keV, respectively.The cool component temperature is well within the range of the expected WHIM temperatures.The temperature of the warm component, on the other hand, agrees Article number, page 10 of 17 Angie Veronica et al.: The eROSITA view of the Abell 3391/95 field: Table 4: eROSITA spectral analysis results, the derived electron density (n e ), and the filament gas overdensity (δ b ).

Regions
Fitting  well with the 1T fit temperature and gives a higher metallicity upper limit of 0.12Z ⊙ .However, through the Bayesian information criterion (BIC, Schwarz 1978) and the Akaike information criterion (AIC, Akaike 1974) tests, we find that adding a second thermal emission component to the model does not improve the fit.We note that deeper data are required to differentiate between the two models.We list the result of the two-temperature fitting in Table B.1 of Appendix B. The derived gas overdensity, δ b , of the Northern Filament is 225 ± 13, which is within the predicted WHIM overdensities.
For the Eastern Filament, we needed to fix the metallicity during the fit to typical values found for filaments from observations and simulations (e.g., Tanimura et al. 2020;Biffi et al. 2022), that is 0.1 and 0.2Z ⊙ .The average temperature of both fits is 1.9 keV and the best-fit values of these fits are always within 1σ of one another.We note that fixing the metallicity has an impact on the normalization, such that we observe an anticorrelation between the value fixed for the metallicity and the resulting normalization.The δ b from both metallicity fits ranges between 169 and 202.
Article number, page 11 of 17 A&A proofs: manuscript no.A339195_outskirts From the Southern Filament boxes, we obtain a 1.1 keV gas with a metallicity value of ∼ 0.10Z ⊙ .The eROSITA data enable us to identify the LSC as a clump within the filament, therefore allowing us to assess the influence of the clump on the Southern Filament (Sect.3.2.1).We performed spectral fittings excluding the LSC emission within its R 500 (hereafter, woLSC).The obtained best-fit temperature and metallicity are 1.2 keV and ∼ 0.11Z ⊙ , which means that we observe no significant changes in these parameters when including or excluding the LSC.However, the normalization per unit area drops by 47%, and the gas overdensity decreases by 30% (from δ b = 166 ± 14 to 117 +17 −19 ) with a significance of 2.2σ.Through this analysis, we directly demonstrate the extent to which the filament density is overestimated when ignoring the presence of a dense clump (e.g., Nagai & Lau 2011;Eckert et al. 2015b;Mirakhor & Walker 2021).

Discussion
4.1.Inner (r ≤ R 500 ) Fig. 4 shows the surface brightness profiles of the A3391 (top) and A3395S (bottom) clusters from their center out to ∼ 2R 200 .From the inner region of the profiles (r < R 500 ), we observe that the clusters are not spherically symmetric, but are rather elliptical.The elongation of the A3391 cluster is apparent in the east−west direction, while the A3395S cluster is extended in the northwest direction.A recent study employing galaxy groups and clusters at z = 0 from large hydrodynamical simulations shows a correlation between their mass and morphology, as well as a correlation between their morphology and the number of filaments connected (connectivity) to them (Gouin et al. 2021).In these simulations, the unrelaxed groups and clusters of high connectivity, such as the A3391 cluster, are more elliptical than their relaxed and weakly connected counterparts.

Outskirts
For the outskirt regions (R 500 < r < R 200 ), we observe higher surface brightness values in the elongation directions of the clusters (Fig. 3 and 4).Despite this, the filament-facing sectors are still brighter than their non-filament-facing counterparts (see Sect. 3.1.3).For instance, A3391-N and A3391-E of the A3391 cluster are brighter than the E Sector, and Sectors 1 and 2 of the A3395S cluster are brighter than Sectors 3 and 4.
The spectral analysis (Sect.3.2.2 and Table 4) suggests that the outskirts temperatures of A3391-N, A3391-E, and NC are higher than the universal temperature profile (Reiprich et al. (2013); Fig. 6, pink shaded areas) with a significance of 2.4σ, 2.8σ, and 1.6σ, respectively.A similar comparison to the outskirts of the A3391 and NC clusters reveals a temperature enhancement in the northeastern outskirts of the A1689 cluster (Kawaharada et al. 2010).These latter authors speculated that this is due to the thermalization process induced by the filament in the mentioned direction.
The constrained metallicity values in the outskirts of A3391 and NC (excluding A3391-E, where the metallicity is frozen) range between 0.1 and 0.3Z ⊙ .This is in good agreement with other measurements; for example, Suzaku measurements of the outskirts of the Perseus cluster (Simionescu et al. 2011), the A3112 cluster (Ezer et al. 2017), and ten other nearby galaxy clusters (Urban et al. 2017), as well as the cosmological hydrodynamical simulations (e.g., Biffi et al. 2018).The metallicity measurements from the observational data reported in the literature are constrained from various radial bins out to R 200 .In the regime defined as outskirts in our work (R 500 < r < R 200 ), the reported metallicities in the literature are consistent with ∼ 0.3Z ⊙ , with large uncertainties ranging between 0.1 and 0.6Z ⊙ .The metallicity of the outskirts of A3395 is high (0.72 +0.60  −0.33 Z ⊙ ), but is still broadly in agreement with the literature values.The rather high metal content found in the outskirts is supposedly the result of the accretion of already enriched gas together with more pristine gas during the cluster assembly (Ezer et al. 2017;Biffi et al. 2018).Metal production occurs at higher redshifts (z ≈ 2 − 3) through some supernova (SN) events (SN Type II and SN Type Ia) and asymptotic giant branch (AGB) stars (e.g., Werner et al. 2008a;Nomoto et al. 2013) and this material is later mixed and distributed throughout the ICM by various processes, such as AGN feedback, galactic winds, or ram-pressure stripping (e.g., Vazza et al. 2010;Ettori et al. 2013).

Filaments (r > R 200 )
From the surface brightness analysis (Sect.3.1.4),we calculated the significance of the emission of all directions in the r > R 200 regime with respect to the CXB level (S/N all ; Eq. 3).We confirmed the surface brightness excess of the detected filaments (Reiprich et al. 2021).For the A3391 cluster, we obtain S/N all = 3.5σ and 2.5σ for the Northern and Eastern Filament, respectively, and 0.5σ and 1.2σ for the western and eastern sectors (non-filament-facing directions).For A3395S, excesses with a significance of S/N all = 2.7σ and 2.3σ are acquired from the filament sectors, while the surface brightness of the other sectors is around the CXB level (S/N all ≈ 0σ).While the significance of the excess emission of the Northern and Eastern Filaments is in good agreement with that reported in Reiprich et al. (2021, see their Table 5), the significance of the excess emission of the Southern Filament is ∼ 1.5σ lower.As our goal here is to construct surface brightness profiles, we do not have identical surface brightness regions for the filaments as defined in Reiprich et al. (2021) and therefore small differences are expected.
We characterized the properties of the filament gas through spectral analysis (Sect.3.2.3).While the gas overdensity of our filaments (98 ≤ δ b ≤ 238) falls within the predicted WHIM gas overdensity, the temperatures obtained for both the Northern and Southern Filament (k B T ≈ 1.0 keV) are more consistent with the HM phase.Nevertheless, these temperatures are close to the upper limit of the WHIM temperature (∼ 0.9 keV, as adopted in the simulations of Shull et al. 2012).Our findings for temperature and density are in accordance with the hot gas phase of filaments in the IllustrisTNG simulation (lower right panel of Fig. 6 and 8 of Galárraga-Espinosa et al. 2021).These latter authors found that the profiles of the hot gas phase of the filaments are flat with an average temperature of below 2 × 10 7 K (∼ 1.7 keV), while the density profiles start at ∼ 10 −4 cm −3 and drop to ∼ 10 −5 cm −3 at ∼ 1 Mpc from the filament spines.We recall that the radial widths of the A3391/95 filaments in this work are r = 0.45 Mpc for Northern Filament and 1.33 Mpc for Southern Filament, and their lengths are 1.8 Mpc and 2.7 Mpc, respectively, which means they are considered to be short filaments 5 .The shorter filaments usually trace a denser large-scale environment and are therefore more likely to be associated with the over-dense structures rather than the skeleton of the cosmic web (Galárraga-Espinosa et al. 2020, 2021, 2022;Vurm et al. 2023).Furthermore, they are in deeper potential wells and experience stronger gravitational heating.Our temperatures and gas overdensities are also within the ranges found for gas at r > R 200 in the EAGLE hydrodynamical simulation (Tuominen et al. 2021).From comparison with other observations, we see that the temperatures of our filaments are consistent with the filament temperatures reported in Tanimura et al. (2020) from the ROSAT stacked analysis of 15 165 filaments (0.9 +1.0 −0.6 keV) and with those reported in Tanimura et al. (2022) from the eFEDS stacked analysis of 463 filaments (1.0 +0.3 −0.2 keV).Tanimura et al. (2020Tanimura et al. ( , 2022) ) restricted their analyses to a radial range of r < 2 Mpc, and the radial widths of the A3391/95 filaments in this work fall within this range.
We find very good agreement between our findings and the gas properties of the simulated A3391/95 analog clusters (outside of their R 200 ) of Biffi et al. (2022).As depicted in their gas phase-space diagram (Fig. 9b), these authors find the bridge gas to be in the warm-hot phase with a typical temperature of ∼ 1 keV and a median overdensity of ∼ 100 (Biffi et al. 2022).
Furthermore, we find that the Northern and Southern Filaments are poorly enriched by metals; this is expected in these low-density regions as they are further away from any metal production sites, such as star-forming regions (Biffi et al. 2022).While the metallicity values for these filaments are unexpectedly low, these values are not completely excluded by the simulations (e.g., Fig. 9b of Biffi et al. 2022).Deeper data is required to confirm the low metallicity value.
The gas temperature of the Eastern Filament ranges between 1.5 and 2.7 keV.Nevertheless, as reported by Reiprich et al. (2021), the significance of the detection of the Eastern Filament is less than 3σ, which is tentative.We note that the Northern and Southern Filament are confirmed in the Planck-SZ and DECam galaxy density maps, as shown in Reiprich et al. (2021).From the spectral modeling, we obtain 3.6 − 4.4σ significance for the Northern and Southern Filament and < 3σ significance for the Eastern and Southern Filament without LSC.As the source regions used and the number of the free parameters in the imaging and spectral analyses differ, small differences in the values obtained from both analyses are not unexpected.
We performed additional tests on the outskirts regions of A3391 and A3395, the Northern Filament, and Southern Filaments.These tests involved treating the filament boxes as a possible sky foreground.The results and discussion of these tests can be found in Appendix B.

Summary and conclusions
We investigated the outskirts (R 500 < r < R 200 ) and the detected inter-cluster filaments (r > R 200 ) of the A3391/95 system using the eROSITA PV data.We focused our analysis on the northern and northeastern directions of A3391, and the southern direction of A3395, which are in the directions of the Northern, Eastern, and Southern Filaments.We created images to generate a PIBsubtracted, exposure-, and Galactic-absorption-corrected image in the soft energy band of 0.3 − 2.0 keV.Using the final corrected products, we calculated surface brightness profiles out to ∼ 2R 200 of the A3391 and A3395 clusters in various directions.We performed comprehensive spectral analyses in the outskirts and filament regions.We constrained the gas properties, including the normalizations, temperatures, and metallicities.Utilizing the acquired normalizations and assuming a cylinder shape with its axis in the plane of the sky, we derived other quantities, namely the electron densities (n e ) and gas overdensities (δ b ).We compared our results with simulated filament properties.We summarize our findings below: -The X-ray surface brightness profiles at r < R 500 in various directions emphasize the morphology of A3391 and A3395S.We notice significant surface brightness decrements below 10 ′ in the northern direction of A3391 and higher values below 22 ′ in the northwestern profile of A3395S.These show the apparent ellipticity of A3391 in the east−west direction and the extension of A3395S in the northwestern direction, respectively.-The temperatures in the filament-facing outskirts of A3391-N, A3391-E, and NC are higher than the temperatures expected for a typical cluster outskirt profile predicted using the universal temperature profile of Reiprich et al. (2013), with significances of 2.4σ, 2.8σ, and 1.6σ, respectively.These enhancements may be related to heating processes induced by the filaments.-We confirm surface brightness excesses in the profiles of the Northern, Eastern, and Southern Filaments.In the surface brightness profile of the Southern Filament, we observe a peak at ∼ 1.5R 200 due to the LSC.-We detect hot gas beyond the R 200 of the cluster (filament regions).The Northern Filament has a best-fit temperature of 0.96 +0.17 −0.14 keV, while the Southern Filament has a best-fit temperature of 1.09 +0.06 −0.13 keV.The gas overdensities of the Northern and Southern Filaments are found to be within the ranges of 212 < δ b < 237 and 152 < δ b < 180, respectively.The density of these filaments is within the expected range of the simulated WHIM properties.The enhanced temperatures of both filaments may be due to the fact that they are shorter and are located in a denser environment; consequently, they may experience stronger gravitational heating (Galárraga-Espinosa et al. 2020).However, these values are close to the upper limit of WHIM temperatures (∼ 0.9 keV, as adopted in the simulations of Shull et al. 2012).More detailed comparisons to the simulated filaments of similar lengths and environments are required, as well as comparisons to other systems in observations.-An upper metallicity value of < 0.01Z ⊙ is obtained for the Northern Filament and 0.10 +0.05 −0.04 Z ⊙ for the Southern Filament.
-The eROSITA data allowed us to identify the LSC, a clump within a filament.We characterized the properties of the LSC, such that k B T 500 = 0.94 +0.06 −0.07 keV, Z 500 = 0.10 +0.05 −0.04 Z ⊙ , M 500 = 2.33 +0.31  −0.34 × 10 13 M ⊙ , and R 500 = 6.44 +0.29 −0.32 ′ ≈ 437.23 +19.38  −21.47 kpc.From the optical DECam/DES galaxy density map, we observe some galaxy overdensities around the LSC.These galaxies may be being accreted towards the LSC, or they are simply a part of the filament and are on the way to the A3395 cluster.Deeper observations will enable us to study the LSC and its surroundings.
-In this work, we directly show the extent to which the filament gas density is overestimated when ignoring the presence of a clump.For example, when excluding the LSC from the analysis of the Southern Filament, we observe a decrease of 47% in the normalization per unit area and 30% in the gas overdensity with a significance of 2.2σ.No significant changes to the temperature and metallicity of the gas were observed.
The outskirts of galaxy clusters are of paramount importance not only for locating large-scale structures but also for studying various accretion physics phenomena, such as shocks, mergers, and clumping.However, due to the faint nature and their location, studying the outskirts is most often hindered by a lack of sensitivity in the soft energy band and/or the limited FoV of the Article number, page 13 of 17 instruments.eROSITA, equipped with wide FoV, scan, and survey observation modes, and superior soft energy response, is an outstanding instrument with which to probe these regions.Indeed, through the A3391/95 PV observations, eROSITA has already demonstrated its capability to capturing these nearby clusters beyond their virial radii.Moreover, these observations reveal the filamentary structure that connects at least five galaxy groups and clusters.Here, we were able to go one step further and characterize the filament properties in terms of density, temperature, and metallicity.With the eROSITA all-sky survey (eRASS), further such systems will be found and this will improve the statistics of observational studies of inter-cluster filaments.
2.5 of Reiprich et al. 2021) generated from the IRAS 100 µm (Miville-Deschênes & Lagache 2005) and the HI4PI data (HI4PI Collaboration et al. 2.1) values within them resemble the average N H values of the sources.The first region (Fig. 1, yellow circle) was used for the Northern and Eastern Filament, and the second (magenta circle) was used for the Southern Filament.We fixed the N H values of the source and CXB spectra to their respective average values in the X-ray spectral fitting.The placement of the source and CXB regions on the N H map is shown in Fig. A.1 of Appendix A. In this figure, we also labeled the used average N H value of each region.

Fig. 3 :
Fig. 3: Configuration of surface brightness profiles related to the A3391 cluster, NC, Northern and Eastern Filament (left), and the A3395 cluster and Southern Filament (right) overlaid on the adaptively-smoothed eROSITA PIB subtracted, exposure corrected, and Galactic absorption corrected image in the 0.3 − 2.0 keV band.

Fig. 4 :
Fig.4: PIB-subtracted SB profiles in the 0.3 − 2.0 keV band (top plots) and the significance values with respect to CXB for r > R 200 data points (bottom plots).The overall significance, S/N all , of the excess emission in this regime is listed in the legends.a: Profiles from A3391-N, NC, and the Northern Filament in blue; A3391-E and Eastern Filament in red; and west and east of A3391 in green and orange (see Fig.3, left).b: Profiles from sector1-5 of A3395S.Sector1 (blue) and 2 (red) are related to the Southern Filament (seeFig.3, right).The CXB level and its 1σ range are plotted as the magenta horizontal dashed-lines and magenta shaded area.The R 500 (R 200 ) of each cluster is plotted as vertical solid (dashed) lines.

Fig. 5 :
Fig. 5: eROSITA fully corrected count rate image in the 0.3 − 2.0 keV band (smoothed with a 64 ′′ Gaussian kernel).The cyan circle indicates the R 500 of the LSC, centered at the coordinates of the WISEA J063004.80-551751.5 galaxy.The small white circles are known galaxies within 0.0424 ≤ z ≤ 0.0636 with their redshifts labeled.The overlaid white contours are from the optical DECam/DES galaxy density map.The white dashed lines mark the R 200 and 2R 200 of the A3395S cluster.
2.2 of Reiprich et al. 2021) obtained from the Dark Energy Camera (DECam, Flaugher et al. (2015)) as part of the Dark Energy Survey (DES, Dark Energy Survey Collaboration et al. ( Fig.6: eROSITA temperature profiles of the outskirts and filaments for the northern (top), eastern (middle), and southern (bottom) regions.The pink shaded areas are the 1σ temperature profiles fromReiprich et al. (2013) using k B T 0.2−0.5R500 as mean temperature of each cluster (see Table5).The black vertical solid (dotted) lines mark the R 500 (R 200 ) of the parent clusters.

Table 1 :
eROSITA A3391/95 PV observations used in this work.There are 16 datasets in total among the observations.The exposure times listed are the average filtered on-axis exposure time across the available TMs of each observation. *

Table 2 :
Information of the A3391/95 groups and clusters.

Table 3 :
Information on the parameters of the eROSITA CXB components and their best-fit normalization values from the first fitting step (see Sect. 2.2).