Outskirts Of Galaxy Clusters Research Articles

Exploring the spectral properties of radio relics – I: integrated spectral index and Mach number

ABSTRACT Radio relics are the manifestation of electrons presumably being shock (re-)accelerated to high energies in the outskirts of galaxy clusters. However, estimates of the shocks’ strength yield different results when measured with radio or X-ray observations. In general, Mach numbers obtained from radio observations are larger than the corresponding X-ray measurements. In this work, we investigate this Mach number discrepancy. For this purpose, we used the cosmological code enzo to simulate a sample of galaxy clusters that host bright radio relics. For each relic, we computed the radio Mach number from the integrated radio spectrum and the X-ray Mach number from the X-ray surface brightness and temperature jumps. Our analysis suggests that the differences in the Mach number estimates follow from the way in which different observables are related to different parts of the underlying Mach number distribution: radio observations are more sensistive to the high Mach numbers present only in a small fraction of a shock’s surface, while X-ray measurements reflect the average of the Mach number distribution. Moreover, X-ray measurements are very sensitive to the relic’s orientation. If the same relic is observed from different sides, the measured X-ray Mach number varies significantly. On the other hand, the radio measurements are more robust, as they are unaffected by the relic’s orientation.

Evolution of splashback boundaries and gaseous outskirts: insights from mergers of self-similar galaxy clusters

ABSTRACT A self-similar spherical collapse model predicts a dark matter (DM) splashback and accretion shock in the outskirts of galaxy clusters while missing a key ingredient of structure formation – processes associated with mergers. To fill this gap, we perform simulations of merging self-similar clusters and investigate their DM and gas evolution in an idealized cosmological context. Our simulations show that the cluster rapidly contracts during the major merger and the splashback radius rsp decreases, approaching the virial radius rvir. While in the self-similar model rsp depends on a smooth mass accretion rate parameter Γs, our simulations show that in the presence of mergers, rsp responds to the changes in the total mass accretion rate Γvir, which accounts for both mergers and smooth accretion. The scatter of the Γvir − rsp/rvir relation indicates a generally low Γs ∼ 1 in clusters in cosmological simulations. In contrast to the DM, the hot gaseous atmospheres significantly expand by the merger-accelerated (MA-) shocks formed when the runaway merger shocks overtake the outer accretion shock. After a major merger, the MA-shock radius is larger than rsp by a factor of up to ∼1.7 for Γs ≲ 1 and is ∼rsp for Γs ≳ 3. This implies that (1) mergers could easily generate the MA-shock-splashback offset measured in cosmological simulations, and (2) the smooth mass accretion rate is small in regions away from filaments where MA-shocks reside. We further discuss the shapes of the DM haloes, various shocks, and contact discontinuities formed at different epochs of the merger, and the ram-pressure stripping in cluster outskirts.

The Colors of Bulges and Disks in the Core and Outskirts of Galaxy Clusters

Abstract The role of the environment on the formation of S0 galaxies is still not well understood, specifically in the outskirts of galaxy clusters. We study eight low-redshift clusters, analyzing galaxy members up to cluster-centric distances of ∼2.5 R 200. We perform 2D photometric bulge–disk decomposition in the g, r, and i bands from which we identify 469 double-component galaxies. We separately analyze the colors of the bulges and disks and their dependence on the projected cluster-centric distance and local galaxy density. For our sample of cluster S0 galaxies, we find that bulges are redder than their surrounding disks, show a significant color–magnitude trend, and have colors that do not correlate with environment metrics. On the other hand, the disks associated with our cluster S0s become significantly bluer with increasing cluster-centric radius but show no evidence for a color–magnitude relation. The disk color–radius relation is mainly driven by galaxies in the cluster core at 0 ≤ R/R 200 < 0.5. No significant difference is found for the disk colors of backsplash and infalling galaxies in the projected phase space (PPS). Beyond R 200, the disk colors do not change with the local galaxy density, indicating that the colors of double-component galaxies are not affected by preprocessing. A significant color–density relation is observed for single-component disk-dominated galaxies beyond R 200. We conclude that the formation of cluster S0 galaxies is primarily driven by cluster core processes acting on the disks, while evidence of preprocessing is found for single-component disk-dominated galaxies. We publicly release the data from the bulge–disk decomposition.

The threehundred: the structure and properties of cosmic filaments in the outskirts of galaxy clusters

ABSTRACT Galaxy cluster outskirts are described by complex velocity fields induced by diffuse material collapsing towards filaments, gas, and galaxies falling into clusters, and gas shock processes triggered by substructures. A simple scenario that describes the large-scale tidal fields of the cosmic web is not able to fully account for this variety, nor for the differences between gas and collisionless dark matter. We have studied the filamentary structure in zoom-in resimulations centred on 324 clusters from the threehundred project, focusing on differences between dark and baryonic matter. This paper describes the properties of filaments around clusters out to five R200, based on the diffuse filament medium where haloes had been removed. For this, we stack the remaining particles of all simulated volumes to calculate the average profiles of dark matter and gas filaments. We find that filaments increase their thickness closer to nodes and detect signatures of gas turbulence at a distance of ${\sim}2 \rm {{{~h^{-1}\,{\rm Mpc}}}}$ from the cluster. These are absent in dark matter. Both gas and dark matter collapse towards filament spines at a rate of ${\sim}200 \,\rm {km ~ s^{-1}\, h^{-1}}$. We see that gas preferentially enters the cluster as part of filaments, and leaves the cluster centre outside filaments. We further see evidence for an accretion shock just outside the cluster. For dark matter, this preference is less obvious. We argue that this difference is related to the turbulent environment. This indicates that filaments act as highways to fuel the inner regions of clusters with gas and galaxies.

Stellar splashback: the edge of the intracluster light

ABSTRACT We examine the outskirts of galaxy clusters in the C-EAGLE simulations to quantify the ‘edges’ of the stellar and dark matter distribution. The radius of the steepest slope in the dark matter, commonly used as a proxy for the splashback radius, is located at $\sim \, r_{200 \rm m}$; the strength and location of this feature depends on the recent mass accretion rate, in good agreement with previous work. Interestingly, the stellar distribution (or intracluster light, ICL) also has a well-defined edge, which is directly related to the splashback radius of the halo. Thus, detecting the edge of the ICL can provide an independent measure of the physical boundary of the halo, and the recent mass accretion rate. We show that these caustics can also be seen in the projected density profiles, but care must be taken to account for the influence of substructures and other non-diffuse material, which can bias and/or weaken the signal of the steepest slope. This is particularly important for the stellar material, which has a higher fraction bound in subhaloes than the dark matter. Finally, we show that the ‘stellar splashback’ feature is located beyond current observational constraints on the ICL, but these large projected distances (≫1 Mpc) and low surface brightnesses (μ ≫ 32 mag arcsec−2) can be reached with upcoming observational facilities such as the Vera C. Rubin Observatory, the Nancy Grace Roman Space Telescope, and Euclid.

Anisotropic infall in the outskirts of OmegaWINGS galaxy clusters

ABSTRACT We study the effects of the environment on galaxy quenching in the outskirts of clusters at 0.04 < z < 0.08. We use a subsample of 14 WINGS and OmegaWINGS clusters that are linked to other groups/clusters by filaments and study separately galaxies located in two regions in the outskirts of these clusters according to whether they are located towards the filaments’ directions or not. We also use samples of galaxies in clusters and field as a comparison. Filamentary structures linking galaxy groups/clusters were identified over the Six Degree Field Galaxy Redshift Survey Data Release 3. We find a fraction of passive galaxies in the outskirts of clusters intermediate between that of the clusters and the field’s. We find evidence of a more effective quenching in the direction of the filaments. We also analyse the abundance of post-starburst (PS) galaxies in the outskirts of clusters focusing our study on two extreme sets of galaxies according to their phase-space position: backsplash and true infallers. We find that up to $\sim 70{{\ \rm per\ cent}}$ of PS galaxies in the direction of filaments are likely backsplash, while this number drops to $\sim 40{{\ \rm per\ cent}}$ in the isotropic infall region. The presence of this small fraction of galaxies in filaments that are falling into clusters for the first time and have been recently quenched, supports a scenario in which a significant number of filament galaxies have been quenched long time ago.

Analytical model for cluster radio relics

ABSTRACT Radio relics are vast synchrotron sources that sit on the outskirts of merging galaxy clusters. In this work we model their formation using a Press–Schechter formalism to simulate merger rates, analytical models for the intracluster medium and the shock dynamics, as well as a simple model for the cosmic ray electrons at the merger shocks. We show that the statistical properties of the population of radio relics are strongly dependent on key physical parameters, such as the acceleration efficiency, the magnetic field strength at the relic, the geometry of the relic and the duration of the electron acceleration at merger shocks. It turns out that the flux distribution as well as the power–mass relation can constrain key parameters of the intracluster medium. With the advent of new large-area radio surveys, statistical analyses of radio relics will complement what we have learned from observations of individual objects.

Spectral imaging of X-COP galaxy clusters with the Sunyaev–Zel’dovich effect

The Sunyaev–Zel’dovich effect is the ideal probe for investigating the outskirts of galaxy clusters. To map this signal, we apply a spectral imaging technique which combines parametric component separation and sparse representations. Our procedure is an improved version of an existing algorithm, which now features a better treatment of astrophysical contaminants, and the implementation of a new beam deconvolution. We use the most recent frequency maps delivered byPlanck, and we consider the clusters analysed in theXMMCluster Outskirts Project (X-COP). In particular, we focus on the images of two clusters which may be possibly interacting with neighbouring objects, namely A2029 and RXCJ1825. We also highlight the advantages of the new beam deconvolution method, through a comparison with the original version of the imaging algorithm.

Detection of metal-rich, cool-warm gas in the outskirts of galaxy clusters

ABSTRACT We present an ultraviolet quasar absorption line analysis of metal lines associated with three strong intervening H i absorbers (with $N(\rm {{H}\,{\small I}})$ > 1016.5 cm−2) detected in the outskirts of Sunyaev–Zel’dovich (SZ) effect-selected galaxy clusters (zcl ∼ 0.4–0.5), within clustocentric impact parameters of ρcl ∼ (1.6–4.7)r500. Discovered in a recent set of targeted far-UV HST/COS spectroscopic observations, these absorbers have among the highest H i column densities ever observed in the outskirts of galaxy clusters, and are also rich in metal absorption lines. Photoionization models yield single phase solutions for the three absorbers with gas densities of nH ∼ 10−3–10−4 cm−3 and metallicities of [X/H] > −1.0 (from one-tenth solar to near-solar). The widths of detected absorption lines suggest gas temperatures of T ∼ 104 K. The inferred densities (temperatures) are significantly higher (lower) compared to the X-ray emitting intracluster medium in cluster cores. The absorbers are tracing a cool phase of the intracluster gas in the cluster outskirts, either associated with gas stripped from cluster galaxies via outflows, tidal streams or ram-pressure forces, or denser regions within the intracluster medium that were uniformly chemically enriched from an earlier epoch of enhanced supernova and Active Galactic Nucleus (AGN) feedback.

Electron Preacceleration in Weak Quasi-perpendicular Shocks in High-beta Intracluster Medium

Abstract Giant radio relics in the outskirts of galaxy clusters are known to be lit up by the relativistic electrons produced via diffusive shock acceleration (DSA) in shocks with low sonic Mach numbers, M s ≲ 3. The particle acceleration at these collisionless shocks critically depends on the kinetic plasma processes that govern the injection to DSA. Here, we study the preacceleration of suprathermal electrons in weak, quasi-perpendicular (Q ⊥) shocks in the hot, high-β (β = P gas/P B) intracluster medium (ICM) through two-dimensional particle-in-cell simulations. Guo et al. showed that, in high-β Q ⊥-shocks, some of the incoming electrons could be reflected upstream and gain energy via shock drift acceleration (SDA). The temperature anisotropy due to the SDA-energized electrons then induces the electron firehose instability (EFI), and oblique waves are generated, leading to a Fermi-like process and multiple cycles of SDA in the preshock region. We find that such electron preacceleration is effective only in shocks above a critical Mach number . This means that, in ICM plasmas, Q ⊥-shocks with M s ≲ 2.3 may not efficiently accelerate electrons. We also find that, even in Q ⊥-shocks with M s ≳ 2.3, electrons may not reach high enough energies to be injected to the full Fermi-I process of DSA, because long-wavelength waves are not developed via the EFI alone. Our results indicate that additional electron preaccelerations are required for DSA in ICM shocks, and the presence of fossil relativistic electrons in the shock upstream region may be necessary to explain observed radio relics.

High density of active galactic nuclei in the outskirts of distant galaxy clusters

We present a study of the distribution of X-ray detected active galactic nuclei (AGN) in the five most massive, M500SZ > 1014 M⊙, and distant, z ∼ 1, galaxy clusters in the Planck and South Pole Telescope (SPT) surveys. The spatial and thermodynamic individual properties of each cluster have been defined with unprecedented accuracy at this redshift using deep X-ray observations. This is an essential property of our sample in order to precisely determine the R500Yx radius of the clusters. For our purposes, we computed the X-ray point-like source surface density in 0.5 R500Yx wide annuli up to a clustercentric distance of 4 R500Yx, statistically subtracting the background and accounting for the respective average density of optical galaxies. We found a significant excess of X-ray point sources between 2 and 2.5 R500Yx at the 99.9% confidence level. The results clearly display for the first time strong observational evidence of AGN triggering in the outskirts of high-redshift massive clusters with such a high statistical significance. We argue that the particular conditions at this distance from the cluster centre increase the galaxy merging rate, which is probably the dominant mechanism of AGN triggering in the outskirts of massive clusters.

First detection of a virial shock with SZ data: implication for the mass accretion rate of Abell 2319

Shocks produced by the accretion of infalling gas in the outskirts of galaxy clusters are expected in the hierarchical structure formation scenario, as found in cosmological hydrodynamical simulations. Here, we report the detection of a shock front at a large radius in the pressure profile of the galaxy cluster A2319 at a significance of 8.6σ, using Planck thermal Sunyaev-Zel’dovich data. The shock is located at (2.93 ± 0.05) × R500 and is not dominated by any preferential radial direction. Using a parametric model of the pressure profile, we derive a lower limit on the Mach number of the infalling gas, ℳ > 3.25 at 95% confidence level. These results are consistent with expectations derived from hydrodynamical simulations. Finally, we use the shock location to constrain the accretion rate of A2319 to Ṁ ≃ (1.4 ± 0.4) × 1014 M⊙ Gyr−1 for a total mass of M200 ≃ 1015 M⊙.

Astrophysics with the Spatially and Spectrally Resolved Sunyaev-Zeldovich Effects

In recent years, observations of the Sunyaev-Zeldovich (SZ) effect have had significant cosmological implications and have begun to serve as a powerful and independent probe of the warm and hot gas that pervades the Universe. As a few pioneering studies have already shown, SZ observations both complement X-ray observations -- the traditional tool for studying the intra-cluster medium -- and bring unique capabilities for probing astrophysical processes at high redshifts and out to the low-density regions in the outskirts of galaxy clusters. Advances in SZ observations have largely been driven by developments in centimetre-, millimetre-, and submillimetre-wave instrumentation on ground-based facilities, with notable exceptions including results from the Planck satellite. Here we review the utility of the thermal, kinematic, relativistic, non-thermal, and polarised SZ effects for studies of galaxy clusters and other large scale structures, incorporating the many advances over the past two decades that have impacted SZ theory, simulations, and observations. We also discuss observational results, techniques, and challenges, and aim to give an overview and perspective on emerging opportunities, with the goal of highlighting some of the exciting new directions in this field.

X-ray and Radio Observations of the Radio Relic Galaxy Clusters 1RXS J0603.3+4214 and RXC J1053.7+5453

We study two galaxy clusters with radio relics, 1RXS J0603.3+4214 and RXC J1053.7+5453, through X-ray and radio observations. Radio relics are diffuse non-thermal radio sources found in outskirts of galaxy clusters. Because of their shape and location, they are thought to be related to cluster merger shocks. The galaxy cluster 1RXS J0603.3+4214 has a well-known linear-shape “toothbrush” radio relic. We investigate the temperature structure across the relic to constrain the Mach number of the associated shock. The results are compared with radio spectral results, which suggest that a simple diffusive shock acceleration model does not hold for this relic. The RXC J1053.7+5453 harbors a standard arc-like relic. We also get the Mach number from the temperature profile. In addition, we found an edge-like structure in the X-ray image between the X-ray peak and relic. We investigate the density and temperature profiles across the edge and found that the structure is likely relevant to not a shock but a contact discontinuity.

Three dynamically distinct stellar populations in the halo of M49

Context. M49 (NGC 4472) is the dominant galaxy in subcluster B of the Virgo Cluster, and a benchmark for studying the build-up of the extended halos of brightest group galaxies in the outskirts of galaxy clusters. Aims. We investigate the kinematics in the outer halo of M49, look for substructures, and describe the transition to the surrounding intra-group light. Methods. As kinematic tracers, we use planetary nebulae (PNe), combining kinematics from the extended Planetary Nebula Spectrograph (PN.S) early-type galaxy survey with our recent deep photometric sample. We study the position-velocity-plane for bright and faint PN populations out to 95 kpc radius, and employ a multi-Gaussian model for the velocity distribution to identify stellar populations with distinct kinematics and histories. Results. We report the detection of stellar-kinematic substructure associated with the interaction of M49 with the dwarf irregular galaxy VCC 1249. We find two kinematically distinct PN populations associated with the main M49 halo and the extended intra-group light (IGL). These have velocity dispersions σhalo ≃ 170 km s−1 and σIGL ≃ 400 km s−1 at 10–80 kpc radii. The overall luminosity profile and velocity dispersion at ~80 kpc are consistent with a flat circular velocity curve extrapolated from X-ray observations. The dispersion of the PNe associated with the IGL joins onto that of the satellite galaxies in subcluster B at ~100 kpc radius. This is the first time that the transition from halo to IGL is observed based on the velocities of individual stars. Conclusions. Therefore the halo of M49, consisting of at least three distinct components, has undergone an extended accretion history within its parent group potential. The blue colours of the IGL component are consistent with a population of stars formed in low-mass galaxies at redshift ~0.5 that has since evolved passively, as suggested by other data.

In-flight performance of the soft x-ray spectrometer detector system on Astro-H

The soft x-ray spectrometer (SXS) instrument was launched aboard the Astro-H (Hitomi) observatory on February 17, 2016. The SXS is based on a high-sensitivity x-ray calorimeter detector system that has been successfully deployed in many ground and suborbital spectrometers. The instrument was to provide essential diagnostics for nearly every class of x-ray emitting objects from the atmosphere of Jupiter to the outskirts of galaxy clusters, without degradation for spatially extended objects. The SXS detector system consisted of a 36-pixel cryogenic microcalorimeter array operated at a heat sink temperature of 50 mK. In preflight testing, the detector system demonstrated a resolving power of better than 1300 at 6 keV with a simultaneous bandpass from below 0.3 keV to above 12 keV with a timing precision better than 100 μs. In addition, a solid-state anticoincidence detector was placed directly behind the detector array for background suppression. The detector error budget included the measured interference from the SXS cooling system and the spacecraft. Additional margin for on-orbit gain stability and on-orbit spacecraft interference were also included predicting an on-orbit performance that meets or exceeds the 7-eV FWHM at 6-keV requirement. The actual on-orbit spectral resolution was better than 5 eV FWHM at 6 keV, easily satisfying the instrument requirement. Here, we discuss the actual on-orbit performance of the SXS detector system and compare this to performance in preflight testing and the on-orbit predictions. We will also discuss the on-orbit gain stability, additional on-orbit interference, and measurements of the on-orbit background.

The origin of ICM enrichment in the outskirts of present-day galaxy clusters from cosmological hydrodynamical simulations

The uniformity of the intra-cluster medium (ICM) enrichment level in the outskirts of nearby galaxy clusters suggests that chemical elements were deposited and widely spread into the intergalactic medium before the cluster formation. This observational evidence is supported by numerical findings from cosmological hydrodynamical simulations, as presented in Biffi et al. (2017), including the effect of thermal feedback from active galactic nuclei. Here, we further investigate this picture, by tracing back in time the spatial origin and metallicity evolution of the gas residing at z=0 in the outskirts of simulated galaxy clusters. In these regions, we find a large distribution of iron abundances, including a component of highly-enriched gas, already present at z=2. At z>1, the gas in the present-day outskirts was distributed over tens of virial radii from the the main cluster and had been already enriched within high-redshift haloes. At z=2, about 40% of the most Fe-rich gas at z=0 was not residing in any halo more massive than 1e11 Msun/h in the region and yet its average iron abundance was already 0.4, w.r.t. the solar value by Anders & Grevesse (1989). This confirms that the in situ enrichment of the ICM in the outskirts of present-day clusters does not play a significant role, and its uniform metal abundance is rather the consequence of the accretion of both low-metallicity and pre-enriched (at z>2) gas, from the diffuse component and through merging substructures. These findings do not depend on the mass of the cluster nor on its core properties.

Quenching of satellite galaxies at the outskirts of galaxy clusters

We find, using cosmological simulations of galaxy clusters, that the hot X-ray emitting intra-cluster medium (ICM) enclosed within the outer accretion shock extends out to $R_{\rm shock}\sim(2 - 3) R_{\rm vir}$, where $R_{\rm vir}$ is the standard virial radius of the halo. Using a simple analytic model for satellite galaxies in the cluster, we evaluate the effect of ram pressure stripping on the gas in the inner discs and in the haloes at different distances from the cluster centre. We find that significant removal of star-forming disc gas occurs only at $r \lesssim 0.5 R_{\rm vir}$, while gas removal from the satellite halo is more effective and can occur when the satellite is found between $ R_{\rm vir}$ and $R_{\rm shock}$. Removal of halo gas sets the stage for quenching of the star formation by starvation over $2\textrm{--}3\,\mathrm{Gyr}$, prior to the satellite entry to the inner cluster halo. This scenario explains the presence of quenched galaxies, preferentially discs, at the outskirts of galaxy clusters, and the delayed quenching of satellites compared to central galaxies.

Non-thermal pressure in the outskirts of Abell 2142

Clumping and turbulence are expected to affect the matter accreted onto the outskirts of galaxy clusters. To determine their impact on the thermodynamic properties of Abell 2142 we perform an analysis of the X-ray temperature data from XMM-Newton via our SuperModel, a state-of-the-art tool for investigating the astrophysics of the intracluster medium already tested on many individual clusters (since Cavaliere et al. 2009). Using the gas density profile corrected for clumpiness derived by Tchernin et al. (2016), we find evidence for the presence of a nonthermal pressure component required to sustain gravity in the cluster outskirts of Abell 2142, that amounts to about 30% of the total pressure at the virial radius. The presence of the nonthermal component implies the gas fraction to be consistent with the universal value at the virial radius and the electron thermal pressure profile to be in good agreement with that inferred from the SZ data. Our results indicate that the presence of gas clumping and of a nonthermal pressure component are both necessary to recover the observed physical properties in the cluster outskirts. Moreover, we stress that an alternative method often exploited in the literature (included Abell 2142) to determine the temperature profile k_BT = P_e/n_e basing on a combination of the Sunyaev-Zel'dovich (SZ) pressure P_e and of the X-ray electron density n_e does not allow to highlight the presence of nonthermal pressure support in the cluster outskirts.

SHOCK ACCELERATION MODEL WITH POSTSHOCK TURBULENCE FOR GIANT RADIO RELICS

We explore the shock acceleration model for giant radio relics, in which relativistic electrons are accelerated via diffusive shock acceleration (DSA) by merger-driven shocks in the outskirts of galaxy clusters. In addition to DSA, turbulent acceleration by compressive MHD mode downstream of the shock is included as well as energy losses of postshock electrons by Coulomb scattering, synchrotron emission, and inverse Compton scattering off the cosmic background radiation. Considering that only a small fraction of merging clusters host radio relics, we favor the reacceleration scenario in which radio relics are generated preferentially when shocks encounter the regions containing low-energy ($\gamma_{\rm e} \lesssim 300$) cosmic ray electrons (CRe). We perform time-dependent DSA simulations of spherically expanding shocks with physical parameters relevant for the Sausage radio relic, and calculate the radio synchrotron emission from the accelerated CRe. We find that significant level of postshock turbulent acceleration is required in order to reproduce broad profiles of the observed radio flux densities of the Sausage relic. Moreover, the spectral curvature in the observed integrated radio spectrum can be explained, if the putative shock should have swept up and exited out of the preshock region of fossil CRe about 10~Myr ago.