Extensive surveys with the \textit{Hubble Space Telescope} (HST) over the past decade, targeting some of the most massive clusters in the sky, have uncovered dozens of galaxy-cluster strong lenses. The massive cluster strong-lens scale is typically $\theta_{E}\sim10\arcsec$ to $\sim30-35\arcsec$, with only a handful of clusters known with Einstein radii $\theta_{E}\sim40\arcsec$ or above (for $z_{source}=2$, nominally). Here we report another very large cluster lens, RXC J0032.1+1808 ($z=0.3956$), the second richest cluster in the redMapper cluster catalog and the 85th most massive cluster in the Planck Sunyaev-Zel'dovich catalog. With our Light-Traces-Mass and fully parametric (dPIEeNFW) approaches, we construct strong lensing models based on 18 multiple images of 5 background galaxies newly identified in the \textit{Hubble} data mainly from the \textit{Reionization Lensing Cluster Survey} (RELICS), in addition to a known sextuply imaged system in this cluster. Furthermore, we compare these models to Lenstool and GLAFIC models that were produced independently as part of the RELICS program. All models reveal a large effective Einstein radius of $\theta_{E}\simeq40\arcsec$ ($z_{source}=2$), owing to the obvious concentration of substructures near the cluster center. Although RXC J0032.1+1808 has a very large critical area and high lensing strength, only three magnified high-redshift candidates are found within the field targeted by RELICS. Nevertheless, we expect many more high-redshift candidates will be seen in wider and deeper observations with \textit{Hubble} or \emph{JWST}. Finally, the comparison between several algorithms demonstrates that the total error budget is largely dominated by systematic uncertainties.

The isotropy of the late Universe is an assumption greatly used in astronomy. However, many studies have recently reported deviations from isotropy with a definitive conclusion yet to be made. New, independent methods to robustly test the cosmic isotropy are of crucial importance. In this work, we investigate the directional behavior of the X-ray luminosity-temperature ($L_X-T$) relation of galaxy clusters. A tight correlation exists between the luminosity and temperature of the X-ray-emitting intracluster medium. While the measured luminosity depends on the underlying cosmology, the temperature can be determined without any cosmological assumptions. By exploiting this property one can effectively test the isotropy of cosmological parameters over the full extragalactic sky. Here, we used 313 homogeneously selected X-ray galaxy clusters from the MCXC catalog and obtained core-excised temperatures for all of them. We find that the behavior of the $L_X-T$ relation heavily depends on the direction of the sky. Strong anisotropies are detected at a $\sim 4\sigma$ level toward $(l,b)\sim (280^{\circ}, -20^{\circ})$. Several X-ray and cluster-related effects that could potentially explain these anisotropies were examined, but none did so. Interestingly, two other available cluster samples appear to have a similar behavior throughout the sky, while being fully independent of each other and our sample. Performing a joint analysis of the three samples, the final anisotropy is further intensified ($\sim 5\sigma$), toward $(l,b)\sim (303^{\circ}, -27^{\circ})$, which is in good agreement with other cosmological probes. This result demonstrates that X-ray galaxy cluster studies that assume perfect isotropy can produce strongly biased results whether the underlying reason is cosmological or related to X-rays. The identification of the exact nature of these anisotropies is therefore crucial.

We report the scaling relations derived by fitting the X-ray parameters determined from analyzing the XMM-Newton observations of 120 galaxy clusters in the Planck Early Sunyaev-Zel'dovich sample spanning the redshift range of 0.059$<$$z$$<$0.546. We find that the slopes of all the investigated scaling relations significantly deviate from the self-similar predictions, if self-similar redshift evolution is assumed. When the redshift evolution is left free to vary, the derived slopes are more in agreement with the self-similar predictions. Relaxed clusters have on average $\sim$30$\%$ higher X-ray luminosity than disturbed clusters at a given mass, a difference that, depending on the relative fraction of relaxed and disturbed clusters in the samples (e.g. SZ vs X-ray selected), have a strong impact in the normalization obtained in different studies. Using the core-excised cluster luminosities reduces the scatter and brings into better agreement the $L$-$M_{tot}$ and $L$-$T$ relations determined for different samples. $M_{tot}$-$T$, $M_{tot}$-$Y_X$, and $M_{tot}$-$M_{gas}$ relations show little dependence on the dynamical state of the clusters, but the normalizations of these relations may depend on the mass range investigated. Although most of the clusters investigated in this work reside at relatively low redshift, the fits prefer values of $\gamma$, the parameter accounting for the redshift evolution, different from the self-similar predictions. This suggests an evolution ($<$2$\sigma$ level, with the exception of the $M_{tot}$-$T$ relation) of the scaling relations. For the first time, we find significant evolution ($>$3$\sigma$) of the $M_{tot}$-$T$ relation, pointing to an increase of the kinetic-to-thermal energy ratio with redshift. This is consistent with a scenario in which higher redshift clusters are on average more disturbed than their lower redshift counterparts.

Massive merging galaxy clusters often host diffuse Mpc-scale radio synchrotron emission. This emission originates from relativistic electrons in the ionized intracluster medium (ICM). An important question is how these synchrotron emitting relativistic electrons are accelerated. Our aim is to search for diffuse emission in the Frontier Fields clusters Abell S1063 and Abell 370 and characterize its properties. While these clusters are very massive and well studied at some other wavelengths, no diffuse emission has been reported for these clusters so far. We obtained 325 MHz Giant Metrewave Radio Telescope (GMRT) and 1--4 GHz Jansky Very Large Array (VLA) observations of Abell S1063 and Abell 370. We complement these data with Chandra and XMM-Newton X-ray observations. In our sensitive images, we discover radio halos in both clusters. In Abell S1063, a giant radio halo is found with a size of $\sim 1.2$ Mpc. The integrated spectral index between 325 MHz and 1.5 GHz is $-0.94\pm0.08$ and it steepens to $-1.77 \pm 0.20$ between 1.5 and 3.0 GHz. This spectral steepening provides support for the turbulent re-acceleration model for radio halo formation. Abell 370 hosts a faint radio halo mostly centred on the southern part of this binary merging cluster, with a size of $\sim 500-700$ kpc. The spectral index between 325 MHz and 1.5 GHz is $-1.10\pm0.09$. Both radio halos follow the known scaling relation between the cluster mass proxy $Y_{500}$ and radio power, consistent with the idea that they are related to ongoing cluster merger events.

We present a machine learning approach for estimating galaxy cluster masses, trained using both Chandra and eROSITA mock X-ray observations of 2,041 clusters from the Magneticum simulations. We train a random forest regressor, an ensemble learning method based on decision tree regression, to predict cluster masses using an input feature set. The feature set uses core-excised X-ray luminosity and a variety of morphological parameters, including surface brightness concentration, smoothness, asymmetry, power ratios, and ellipticity. The regressor is cross-validated and calibrated on a training sample of 1,615 clusters (80% of sample), and then results are reported as applied to a test sample of 426 clusters (20% of sample). This procedure is performed for two different mock observation series in an effort to bracket the potential enhancement in mass predictions that can be made possible by including dynamical state information. The first series is computed from idealized Chandra-like mock cluster observations, with high spatial resolution, long exposure time (1 Ms), and the absence of background. The second series is computed from realistic-condition eROSITA mocks with lower spatial resolution, short exposures (2 ks), instrument effects, and background photons modeled. We report a 20% reduction in the mass estimation scatter when either series is used in our random forest model compared to a standard regression model that only employs core-excised luminosity. The morphological parameters that hold the highest feature importance are smoothness, asymmetry, and surface brightness concentration. Hence, these parameters, which encode the dynamical state of the cluster, can be used to make more accurate predictions of cluster masses in upcoming surveys, offering a crucial step forward for cosmological analyses.

We study the dynamical state and the integrated total mass profiles of 75 massive (M500 > 5 e+14 M_sun) SZ-selected clusters at 0.08<z< 1.1. The sample is built from the Planck catalogue, with the addition of 4 SPT clusters at z>0.9. Using XMM imaging observations, we characterise the dynamical state with the centroid shift, the concentration, and their combination, M, which simultaneously probes the core and the large scale gas morphology. Using spatially-resolved spectroscopy and assuming hydrostatic equilibrium, we derive the total integrated mass profiles. The mass profile shape is quantified by the sparsity, the ratio of M500 to M2500, the masses at density contrast 500 and 2500, respectively. We study the correlations between the various parameters and their dependence on redshift. We confirm that SZ-selected samples, thought to reflect most closely the underlying cluster population, are dominated by disturbed and non-cool core objects at all z. There is no significant evolution or mass dependence of either the cool core fraction or the centroid shift parameter. The M parameter evolves slightly with z, having a correlation coefficient of rho= -0.2 $\pm$ 0.1 and a null hypothesis p-value of 0.01. In the high mass regime considered here, the sparsity evolves minimally with redshift, increasing by 10% between z<0.2 and z>0.55, an effect significant at less than 2 sigma. In contrast, the dependence of the sparsity on dynamical state is much stronger, increasing by a factor of $\sim$60% from the 1/3 most relaxed to the 1/3 most disturbed objects, an effect significant at more than 3 sigma. This is the first observational evidence that the shape of the integrated total mass profile in massive clusters is principally governed by the dynamical state, and is only mildly dependent on redshift. We discuss the consequences for the comparison between observations and theoretical predictions.

The M49 group, resident outside the virial radius of the Virgo cluster, is falling onto the cluster from the south. We report results from deep {\sl XMM-Newton} mosaic observations of M49. Its hot gas temperature is 0.8\,keV at the group center and rises to 1.5\,keV beyond the brightest group galaxy (BGG). The group gas extends to radii of $\sim300$\,kpc to the north and south. The observations reveal a cold front $\sim20$\,kpc north of the BGG center and an X-ray bright stripped tail 70\,kpc long and 10\,kpc wide to the southwest of the BGG. We argue that the atmosphere of the infalling group was slowed by its encounter with the Virgo cluster gas, causing the BGG to move forward subsonically relative to the group gas. We measure declining temperature and metallicity gradients along the stripped tail. The tail gas can be traced back to the cooler and enriched gas uplifted from the BGG center by buoyant bubbles, implying that AGN outbursts may have intensified the stripping process. We extrapolate to a virial radius of 740\,kpc and derive a virial mass of $4.6\times10^{13}\,M_\odot$ for the M49 group. Its group atmosphere appears truncated and deficient when compared with isolated galaxy groups of similar temperatures. If M49 is on its first infall to Virgo, the infall region of a cluster could have profound impacts on galaxies and groups that are being accreted onto galaxy clusters. Alternatively, M49 may have already passed through Virgo once.

Up to now, the largest sample of galaxy clusters selected in X-rays comes from the ROSAT All-Sky Survey (RASS). Although there have been many interesting clusters discovered with the RASS data, the broad point spread function (PSF) of the ROSAT satellite limits the amount of spatial information of the detected objects. This leads to the discovery of new cluster features when a re-observation is performed with higher resolution X-ray satellites. Here we present the results from XMM-Newton observations of three clusters: RXCJ2306.6-1319, ZwCl1665 and RXCJ0034.6-0208, for which the observations reveal a double or triple system of extended components. These clusters belong to the extremely expanded HIghest X-ray FLUx Galaxy Cluster Sample (eeHIFLUGCS), which is a flux-limited cluster sample ($f_\textrm{X,500}\geq 5\times10^{-12}$ erg s$^{-1}$ cm$^{-2}$ in the $0.1-2.4$ keV energy band). For each structure in each cluster, we determine the redshift with the X-ray spectrum and find that the components are not part of the same cluster. This is confirmed by an optical spectroscopic analysis of the galaxy members. Therefore, the total number of clusters is actually 7 and not 3. We derive global cluster properties of each extended component. We compare the measured properties to lower-redshift group samples, and find a good agreement. Our flux measurements reveal that only one component of the ZwCl1665 cluster has a flux above the eeHIFLUGCS limit, while the other clusters will no longer be part of the sample. These examples demonstrate that cluster-cluster projections can bias X-ray cluster catalogues and that with high-resolution X-ray follow-up this bias can be corrected.

Supermassive Black Holes (BHs) residing in brightest cluster galaxies (BCGs) are overly massive when considering the local relationships between the BH mass and stellar bulge mass or velocity dispersion. Due to the location of these BHs within the cluster, large-scale cluster processes may aid the growth of BHs in BCGs. In this work, we study a sample of 71 galaxy clusters to explore the relationship between the BH mass, stellar bulge mass of the BCG, and the total gravitating mass of the host clusters. Due to difficulties in obtaining dynamically measured BH masses in distant galaxies, we use the Fundamental Plane relationship of BHs to infer their masses. We utilize X-ray observations taken by $Chandra$ to measure the temperature of the intra-cluster medium (ICM), which is a proxy for the total mass of the cluster. We analyze the $\rm M_{BH}-kT$ and $\rm M_{BH}-M_{Bulge}$ relationships and establish the best-fitting power laws:$\log_{10}(M_{\rm BH} /10^9 M_{\odot})=-0.35+2.08 \log_{10}(kT / 1 \rm keV)$ and $\log_{10}(\rm M_{BH}/10^9M_{\odot})= -1.09+ 1.92 \log_{10}(M_{\rm bulge}/10^{11}M_{\odot})$. Both relations are comparable with that established earlier for a sample of brightest group/cluster galaxies with dynamically measured BH masses. Although both the $\rm M_{BH}-kT$ and the $\rm M_{BH}-M_{Bulge}$ relationships exhibit large intrinsic scatter, based on Monte Carlo simulations we conclude that dominant fraction of the scatter originates from the Fundamental Plane relationship. We split the sample into cool core and non-cool core resembling clusters, but do not find statistically significant differences in the $\rm M_{BH}-kT$ relation. We speculate that the overly massive BHs in BCGs may be due to frequent mergers and cool gas inflows onto the cluster center.

Strong gravitational lensing by clusters of galaxies probes the mass distribution at the core of each cluster and magnifies the universe behind it. MACS J0417.5-1154 at z=0.443 is one of the most massive clusters known based on weak lensing, X-ray, and Sunyaev-Zel'dovich analyses. Here we compute a strong lens model of MACS J0417 based on Hubble Space Telescope imaging observations collected, in part, by the Reionization Lensing Cluster Survey (RELICS), and recently reported spectroscopic redshifts from the MUSE instrument on the Very Large Telescope (VLT). We measure an Einstein radius of ThetaE=36'' at z = 9 and a mass projected within 200 kpc of M(200 kpc) = 1.78+0.01-0.03x10**14Msol. Using this model, we measure a ratio between the mass attributed to cluster-member galaxy halos and the main cluster halo of order 1:100. We assess the probability to detect magnified high-redshift galaxies in the field of this cluster, both for comparison with RELICS HST results and as a prediction for the James Webb Space Telescope (JWST) Guaranteed Time Observations upcoming for this cluster. Our lensing analysis indicates that this cluster has similar lensing strength to other clusters in the RELICS program. Our lensing analysis predicts a detection of at least a few z~6-8 galaxies behind this cluster, at odds with a recent analysis that yielded no such candidates in this field. Reliable strong lensing models are crucial for accurately predicting the intrinsic properties of lensed galaxies. As part of the RELICS program, our strong lensing model produced with the Lenstool parametric method is publicly available through the Mikulski Archive for Space Telescopes (MAST).

Large surveys of galaxy clusters with the Hubble and Spitzer Space Telescopes, including CLASH and the Frontier Fields, have demonstrated the power of strong gravitational lensing to efficiently deliver large samples of high-redshift galaxies. We extend this strategy through a wider, shallower survey named RELICS, the Reionization Lensing Cluster Survey. This survey, described here, was designed primarily to deliver the best and brightest high-redshift candidates from the first billion years after the Big Bang. RELICS observed 41 massive galaxy clusters with Hubble and Spitzer at 0.4-1.7um and 3.0-5.0um, respectively. We selected 21 clusters based on Planck PSZ2 mass estimates and the other 20 based on observed or inferred lensing strength. Our 188-orbit Hubble Treasury Program obtained the first high-resolution near-infrared images of these clusters to efficiently search for lensed high-redshift galaxies. We observed 46 WFC3/IR pointings (~200 arcmin^2) with two orbits divided among four filters (F105W, F125W, F140W, and F160W) and ACS imaging as needed to achieve single-orbit depth in each of three filters (F435W, F606W, and F814W). As previously reported by Salmon et al., we discovered 322 z ~ 6 - 10 candidates, including the brightest known at z ~ 6, and the most distant spatially-resolved lensed arc known at z ~ 10. Spitzer IRAC imaging (945 hours awarded, plus 100 archival) has crucially enabled us to distinguish z ~ 10 candidates from z ~ 2 interlopers. For each cluster, two HST observing epochs were staggered by about a month, enabling us to discover 11 supernovae, including 3 lensed supernovae, which we followed up with 20 orbits from our program. We delivered reduced HST images and catalogs of all clusters to the public via MAST and reduced Spitzer images via IRSA. We have also begun delivering lens models of all clusters, to be completed before the JWST GO call for proposals.

Strong gravitational lensing (SL) is a powerful means to map the distribution of dark matter. In this work, we perform a SL analysis of the prominent X-ray cluster RXJ0152.7-1357 (z=0.83, also known as CL 0152.7-1357) in \textit{Hubble Space Telescope} images, taken in the framework of the Reionization Lensing Cluster Survey (RELICS). On top of a previously known $z=3.93$ galaxy multiply imaged by RXJ0152.7-1357, for which we identify an additional multiple image, guided by a light-traces-mass approach we identify seven new sets of multiply imaged background sources lensed by this cluster, spanning the redshift range [1.79-3.93]. A total of 25 multiple images are seen over a small area of ~0.4 $arcmin^2$, allowing us to put relatively high-resolution constraints on the inner matter distribution. Although modestly massive, the high degree of substructure together with its very elongated shape make RXJ0152.7-1357 a very efficient lens for its size. This cluster also comprises the third-largest sample of z~6-7 candidates in the RELICS survey. Finally, we present a comparison of our resulting mass distribution and magnification estimates with those from a Lenstool model. These models are made publicly available through the MAST archive.

We study a sample of 207 nearby galaxy groups and clusters observed with XMM-Newton. Key aspects of this sample include the large size, the high data quality, and the large diversity of cluster dynamical states. We determine the overall metallicity within 0.3R$_{500}$ and the radial distribution of the metals. On average, we find a mild dependence of the core metallicity with the average temperature of the system in agreement with previous results. However, we identify the cause of this mild dependence to be due to relaxed systems only; disturbed systems do not show this trend, on average. The large scatter observed in this relation is strongly associated with the dynamical state of the systems: relaxed systems have on average a higher metallicity in the core than disturbed objects. The radial profiles of relaxed systems are centrally peaked and show a steep decrease with radius, flattening beyond 0.3-0.4R$_{500}$. The metallicity of disturbed systems is also higher in the center but at much lower values than what is observed for relaxed objects. This finding is consistent with the picture that cluster mergers mix the abundance distribution by inducing large scale motions. The scatter of the radial profiles is quite large, but while for relaxed systems it decreases almost monotonically as function of the radius, for disturbed systems it shows a significant boost at large radii. Systems with a central radio source have a flatter profile indicating that central AGNs are an efficient mechanism to uplift and redistribute the metals in the ICM.

Brightest Cluster Galaxies (BCGs) residing in the centers of galaxy clusters are typically quenched giant ellipticals. A recent study hinted that star-forming galaxies with large disks, so-called superluminous spirals and lenticulars, are the BCGs of a subset of galaxy clusters. Based on the existing optical data it was not possible to constrain whether the superluminous disk galaxies reside at the center of galaxy clusters. In this work, we utilize XMM-Newton X-ray observations of five galaxy clusters to map the morphology of the intracluster medium (ICM), characterize the galaxy clusters, determine the position of the cluster center, and measure the offset between the cluster center and the superluminous disk galaxies. We demonstrate that one superluminous lenticular galaxy, 2MASX J10405643-0103584, resides at the center of a low-mass ($M_{\rm 500} = 10^{14} \ \rm{M_{\odot}}$) galaxy cluster. This represents the first conclusive evidence that a superluminous disk galaxy is the central BCG of a galaxy cluster. We speculate that the progenitor of 2MASX J10405643-0103584 was an elliptical galaxy, whose extended disk was re-formed due to the merger of galaxies. We exclude the possibility that the other four superluminous disk galaxies reside at the center of galaxy clusters, as their projected distance from the cluster center is $150-1070$ kpc, which corresponds to $(0.27-1.18)R_{\rm 500}$. We conclude that these clusters host quiescent massive elliptical galaxies at their center.

We present a multi-wavelength analysis based on archival radio, optical and X-ray data of the complex radio source 3C 196.1, whose host is the brightest cluster galaxy of a $z=0.198$ cluster. HST data show H$\alpha$+[N II] emission aligned with the jet 8.4 GHz radio emission. An H$\alpha$+[N II] filament coincides with the brightest X-ray emission, the northern hotspot. Analysis of the X-ray and radio images reveals cavities located at galactic- and cluster- scales. The galactic-scale cavity is almost devoid of 8.4 GHz radio emission and the south-western H$\alpha$+[N II] emission is bounded (in projection) by this cavity. The outer cavity is co-spatial with the peak of 147 MHz radio emission, and hence we interpret this depression in X-ray surface brightness as being caused by a buoyantly rising bubble originating from an AGN outburst $\sim$280 Myrs ago. A \textit{Chandra} snapshot observation allowed us to constrain the physical parameters of the cluster, which has a cool core with a low central temperature $\sim$2.8 keV, low central entropy index $\sim$13 keV cm$^2$ and a short cooling time of $\sim$500 Myr, which is $<0.05$ of the age of the Universe at this redshift. By fitting jumps in the X-ray density we found Mach numbers between 1.4 and 1.6, consistent with a shock origin. We also found compelling evidence of a past merger, indicated by a morphology reminiscent of gas sloshing in the X-ray residual image. Finally, we computed the pressures, enthalpies $E_{cav}$ and jet powers $P_{jet}$ associated with the cavities: $E_{cav}\sim7\times10^{58}$ erg, $P_{jet}\sim1.9\times10^{44}$ erg s$^{-1}$ for the inner cavity and $E_{cav}\sim3\times10^{60}$ erg, $P_{jet}\sim3.4\times10^{44}$ erg s$^{-1}$ for the outer cavity.

We present a strong-lensing analysis of four massive galaxy clusters imaged with the Hubble Space Telescope in the Reionization Lensing Cluster Survey. We use a Light-Traces-Mass technique to uncover sets of multiply images and constrain the mass distribution of the clusters. These mass models are the first published for Abell S295 and MACS J0159.8-0849, and are improvements over previous models for Abell 697 and MACS J0025.4-1222. Our analysis for MACS J0025.4-1222 and Abell S295 shows a bimodal mass distribution supporting the merger scenarios proposed for these clusters. The updated model for MACS J0025.4-1222 suggests a substantially smaller critical area than previously estimated. For MACS J0159.8-0849 and Abell 697 we find a single peak and relatively regular morphology, revealing fairly relaxed clusters. Despite being less prominent lenses, three of these clusters seem to have lensing strengths, i.e. cumulative area above certain magnification, similar to the Hubble Frontier Fields clusters (e.g., A($\mu>5$) $\sim 1-3$ arcmin$^2$, A($\mu>10$) $\sim 0.5-1.5$ arcmin$^2$), which in part can be attributed to their merging configurations. We make our lens models publicly available through the Mikulski Archive for Space Telescopes. Finally, using Gemini-N/GMOS spectroscopic observations we detect a single emission line from a high-redshift $J_{125}\simeq25.7$ galaxy candidate lensed by Abell 697. While we cannot rule out a lower-redshift solution, we interpret the line as Ly$\alpha$ at $z=5.800\pm 0.001$, in agreement with its photometric redshift and dropout nature. Within this scenario we measure a Ly$\alpha$ rest-frame equivalent width of $52\pm22$ \AA, and an observed Gaussian width of $117\pm 15$ km/s.

We present a lens model for the cluster SPT-CLJ0615$-$5746, which is the highest redshift ($z=0.972$) system in the Reionization of Lensing Clusters Survey (RELICS), making it the highest redshift cluster for which a full strong lens model is published. We identify three systems of multiply-imaged lensed galaxies, two of which we spectroscopically confirm at $z=1.358$ and $z=4.013$, which we use as constraints for the model. We find a foreground structure at $z\sim0.4$, which we include as a second cluster-sized halo in one of our models; however two different statistical tests find the best-fit model consists of one cluster-sized halo combined with three individually optimized galaxy-sized halos, as well as contributions from the cluster galaxies themselves. We find the total projected mass density within $r=26.7"$ (the region where the strong lensing constraints exist) to be $M=2.51^{+0.15}_{-0.09}\times 10^{14}$~M$_{\odot}$. If we extrapolate out to $r_{500}$, our projected mass density is consistent with the mass inferred from weak lensing and from the Sunyaev-Zel'dovich effect ($M\sim10^{15}$~M$_{\odot}$). This cluster is lensing a previously reported $z\sim10$ galaxy, which, if spectroscopically confirmed, will be the highest-redshift strongly lensed galaxy known.

Strong gravitational lensing by galaxy clusters has become a powerful tool for probing the high-redshift Universe, magnifying distant and faint background galaxies. Reliable strong lensing (SL) models are crucial for determining the intrinsic properties of distant, magnified sources and for constructing their luminosity function. We present here the first SL analysis of MACS J0308.9+2645 and PLCK G171.9-40.7, two massive galaxy clusters imaged with the Hubble Space Telescope in the framework of the Reionization Lensing Cluster Survey (RELICS). We use the Light-Traces-Mass modeling technique to uncover sets of multiply imaged galaxies and constrain the mass distribution of the clusters. Our SL analysis reveals that both clusters have particularly large Einstein radii ($\theta_E>30"$ for a source redshift of $z_s=2$), providing fairly large areas with high magnifications, useful for high-redshift galaxy searches ($\sim2$ arcmin$^{2}$ with $\mu>5$ to $\sim1$ arcmin$^{2}$ with $\mu>10$, similar to a typical \textit{Hubble Frontier Fields} cluster). We also find that MACS J0308.9+2645 hosts a promising, apparently bright (J$\sim23.2-24.6$ AB), multiply imaged high-redshift candidate at $z\sim6.4$. These images are amongst the brightest high-redshift candidates found in RELICS. Our mass models, including magnification maps, are made publicly available for the community through the Mikulski Archive for Space Telescopes.

Massive foreground galaxy clusters magnify and distort the light of objects behind them, permitting a view into both the extremely distant and intrinsically faint galaxy populations. We present here the z ~ 6 - 8 candidate high-redshift galaxies from the Reionization Lensing Cluster Survey (RELICS), a Hubble and Spitzer Space Telescope survey of 41 massive galaxy clusters spanning an area of ~200 arcmin^2. These clusters were selected to be excellent lenses and we find similar high-redshift sample sizes and magnitude distributions as CLASH. We discover 321 candidate galaxies with photometric redshifts between z ~ 6 to z ~ 8, including extremely bright objects with H-band magnitudes of m_AB ~ 23 mag. As a sample, the observed (lensed) magnitudes of these galaxies are among the brightest known at z> 6, comparable to much wider, blank-field surveys. RELICS demonstrates the efficiency of using strong gravitational lenses to produce high-redshift samples in the epoch of reionization. These brightly observed galaxies are excellent targets for follow-up study with current and future observatories, including the James Webb Space Telescope.

Strong gravitational lensing by galaxy clusters magnifies background galaxies, enhancing our ability to discover statistically significant samples of galaxies at z>6, in order to constrain the high-redshift galaxy luminosity functions. Here, we present the first five lens models out of the Reionization Lensing Cluster Survey (RELICS) Hubble Treasury Program, based on new HST WFC3/IR and ACS imaging of the clusters RXC J0142.9+4438, Abell 2537, Abell 2163, RXC J2211.7-0349, and ACT-CLJ0102-49151. The derived lensing magnification is essential for estimating the intrinsic properties of high-redshift galaxy candidates, and properly accounting for the survey volume. We report on new spectroscopic redshifts of multiply imaged lensed galaxies behind these clusters, which are used as constraints, and detail our strategy to reduce systematic uncertainties due to lack of spectroscopic information. In addition, we quantify the uncertainty on the lensing magnification due to statistical and systematic errors related to the lens modeling process, and find that in all but one cluster, the magnification is constrained to better than 20% in at least 80% of the field of view, including statistical and systematic uncertainties. The five clusters presented in this paper span the range of masses and redshifts of the clusters in the RELICS program. We find that they exhibit similar strong lensing efficiencies to the clusters targeted by the Hubble Frontier Fields within the WFC3/IR field of view. Outputs of the lens models are made available to the community through the Mikulski Archive for Space Telescopes

Supermassive black holes (BHs) residing in the brightest cluster galaxies are over-massive relative to the stellar bulge mass or central stellar velocity dispersion of their host galaxies. As BHs residing at the bottom of the galaxy cluster's potential well may undergo physical processes that are driven by the large-scale characteristics of the galaxy clusters, it is possible that the growth of these BHs is (indirectly) governed by the properties of their host clusters. In this work, we explore the connection between the mass of BHs residing in the brightest group/cluster galaxies (BGGs/BCGs) and the virial temperature, and hence total gravitating mass, of galaxy groups/clusters. To this end, we investigate a sample of 17 BGGs/BCGs with dynamical BH mass measurements and utilize XMM-Newton X-ray observations to measure the virial temperatures and infer the $M_{\rm 500}$ mass of the galaxy groups/clusters. We find that the $M_{\rm BH} - kT$ relation is significantly tighter and exhibits smaller scatter than the $M_{\rm BH} - M_{\rm bulge}$ relations. The best-fitting power-law relations are $ \log_{10} (M_{\rm BH}/10^{9} \ \rm{M_{\odot}}) = 0.20 + 1.74 \log_{10} (kT/1 \ \rm{keV}) $ and $ \log_{10} (M_{\rm BH}/10^{9} \ \rm{M_{\odot}}) = -0.80 + 1.72 \log_{10} (M_{\rm bulge}/10^{11} \ M_{\odot})$. Thus, the BH mass of BGGs/BCGs may be set by physical processes that are governed by the properties of the host galaxy group/cluster. These results are confronted with the Horizon-AGN simulation, which reproduces the observed relations well, albeit the simulated relations exhibit notably smaller scatter.

Observations of relaxed, massive and distant clusters can provide important tests of standard cosmological models e.g. using the gas mass fraction. We study the very luminous, high redshift ($z=0.902$) galaxy cluster ClJ120958.9+495352 using XMM-Newton data and measure the temperature profile and cooling time to investigate the dynamical status with respect to the presence of a cool core as well as global cluster properties. We use HST weak lensing data to estimate its total mass and determine the gas mass fraction. We perform a spectral analysis using an XMM-Newton observation of 15ks cleaned exposure time. As the treatment of the background is crucial, we use two different approaches to account for the background emission to verify our results. We account for point-spread-function effects and deproject our results to estimate the gas mass fraction of the cluster. We measure weak lensing galaxy shapes from mosaic HST imaging and select background galaxies photometrically in combination with imaging data from the William Herschel Telescope. The X-ray luminosity of ClJ120958.9+495352 in the 0.1-2.4keV band estimated from our XMM-Newton data is $L_X = (13.4_{-1.0}^{+1.2})\times10^{44}$erg/s and thus it is one of the most X-ray luminous clusters known at similarly high redshift. We find clear indications for the presence of a cool core from the temperature profile and the central cooling time, which is very rare at such high redshifts. Based on the weak lensing analysis we estimate a cluster mass of $M_\mathrm{500}/10^{14}M_\odot=4.4^{+2.2}_{-2.0}(\mathrm{stat.})\pm0.6(\mathrm{sys.})$ and a gas mass fraction of $f_\mathrm{gas,2500} = 0.11_{-0.03}^{+0.06}$ in good agreement with previous findings for high redshift and local clusters.

The progenitor systems of the class of "Ca-rich transients" is a key open issue in time domain astrophysics. These intriguing objects exhibit unusually strong calcium line emissions months after explosion, fall within an intermediate luminosity range, are often found at large projected distances from their host galaxies, and may play a vital role in enriching galaxies and the intergalactic medium. Here we present multi-wavelength observations of iPTF15eqv in NGC 3430, which exhibits a unique combination of properties that bridge those observed in Ca-rich transients and Type Ib/c supernovae. iPTF15eqv has among the highest [Ca II]/[O I] emission line ratios observed to date, yet is more luminous and decays more slowly than other Ca-rich transients. Optical and near-infrared photometry and spectroscopy reveal signatures consistent with the supernova explosion of a < 10 solar mass star that was stripped of its H-rich envelope via binary interaction. Distinct chemical abundances and ejecta kinematics suggest that the core collapse occurred through electron capture processes. Deep limits on possible radio emission made with the Jansky Very Large Array imply a clean environment ($n <$ 0.1 cm$^{-3}$) within a radius of $\sim 10^{17}$ cm. Chandra X-ray Observatory observations rule out alternative scenarios involving tidal disruption of a white dwarf by a black hole, for masses > 100 solar masses). Our results challenge the notion that spectroscopically classified Ca-rich transients only originate from white dwarf progenitor systems, complicate the view that they are all associated with large ejection velocities, and indicate that their chemical abundances may vary widely between events.

X-ray observations show that galaxy clusters have a very large range of morphologies. The most disturbed systems which are good to study how clusters form and grow and to test physical models, may potentially complicate cosmological studies because the cluster mass determination becomes more challenging. Thus, we need to understand the cluster properties of our samples to reduce possible biases. This is complicated by the fact that different experiments may detect different cluster populations. For example, SZ selected cluster samples have been found to include a greater fraction of disturbed systems than X-ray selected samples. In this paper we determined eight morphological parameters for the Planck Early Sunyaev-Zeldovich (ESZ) objects observed with XMM-Newton. We found that two parameters, concentration and centroid-shift, are the best to distinguish between relaxed and disturbed systems. For each parameter we provide the values that allow one to select the most relaxed or most disturbed objects from a sample. We found that there is no mass dependence on the cluster dynamical state. By comparing our results with what was obtained with REXCESS clusters, we also confirm that indeed the ESZ clusters tend to be more disturbed, as found by previous studies.

The chemical yields of supernovae and the metal enrichment of the hot intra-cluster medium (ICM) are not well understood. This paper introduces the CHEmical Enrichment RGS Sample (CHEERS), which is a sample of 44 bright local giant ellipticals, groups and clusters of galaxies observed with XMM-Newton. This paper focuses on the abundance measurements of O and Fe using the reflection grating spectrometer (RGS). The deep exposures and the size of the sample allow us to quantify the intrinsic scatter and the systematic uncertainties in the abundances using spectral modeling techniques. We report the oxygen and iron abundances as measured with RGS in the core regions of all objects in the sample. We do not find a significant trend of O/Fe as a function of cluster temperature, but we do find an intrinsic scatter in the O and Fe abundances from cluster to cluster. The level of systematic uncertainties in the O/Fe ratio is estimated to be around 20-30%, while the systematic uncertainties in the absolute O and Fe abundances can be as high as 50% in extreme cases. We were able to identify and correct a systematic bias in the oxygen abundance determination, which was due to an inaccuracy in the spectral model. The lack of dependence of O/Fe on temperature suggests that the enrichment of the ICM does not depend on cluster mass and that most of the enrichment likely took place before the ICM was formed. We find that the observed scatter in the O/Fe ratio is due to a combination of intrinsic scatter in the source and systematic uncertainties in the spectral fitting, which we are unable to disentangle. The astrophysical source of intrinsic scatter could be due to differences in AGN activity and ongoing star formation in the BCG. The systematic scatter is due to uncertainties in the spatial line broadening, absorption column, multi-temperature structure and the thermal plasma models. (Abbreviated).

We derive and compare the fractions of cool-core clusters in the {\em Planck} Early Sunyaev-Zel'dovich sample of 164 clusters with $z \leq 0.35$ and in a flux-limited X-ray sample of 100 clusters with $z \leq 0.30$, using {\em Chandra} observations. We use four metrics to identify cool-core clusters: 1) the concentration parameter: the ratio of the integrated emissivity profile within 0.15 $r_{500}$ to that within $r_{500}$, and 2) the ratio of the integrated emissivity profile within 40 kpc to that within 400 kpc, 3) the cuspiness of the gas density profile: the negative of the logarithmic derivative of the gas density with respect to the radius, measured at 0.04 $r_{500}$, and 4) the central gas density, measured at 0.01 $r_{500}$. We find that the sample of X-ray selected clusters, as characterized by each of these metrics, contains a significantly larger fraction of cool-core clusters compared to the sample of SZ selected clusters (44$\pm$7\% vs. 28$\pm$4\% using the concentration parameter in the 0.15--1.0 $r_{500}$ range, 61$\pm$8\% vs. 36$\pm$5\% using the concentration parameter in the 40--400 kpc range, 64$\pm$8\% vs. 38$\pm$5\% using the cuspiness, and 53$\pm$7\% vs. 39$\pm$5\% using the central gas density). Qualitatively, cool-core clusters are more X-ray luminous at fixed mass. Hence, our X-ray flux-limited sample, compared to the approximately mass-limited SZ sample, is over-represented with cool-core clusters. We describe a simple quantitative model that uses the excess luminosity of cool-core clusters compared to non-cool-core clusters at fixed mass to successfully predict the observed fraction of cool-core clusters in X-ray selected samples.

We report Chandra X-ray observations and optical weak-lensing measurements from Subaru/Suprime-Cam images of the double galaxy cluster Abell 2465 (z=0.245). The X-ray brightness data are fit to a beta-model to obtain the radial gas density profiles of the northeast (NE) and southwest (SW) sub-components, which are seen to differ in structure. We determine core radii, central temperatures, the gas masses within $r_{500c}$, and the total masses for the broader NE and sharper SW components assuming hydrostatic equilibrium. The central entropy of the NE clump is about two times higher than the SW. Along with its structural properties, this suggests that it has undergone merging on its own. The weak-lensing analysis gives virial masses for each substructure, which compare well with earlier dynamical results. The derived outer mass contours of the SW sub-component from weak lensing are more irregular and extended than those of the NE. Although there is a weak enhancement and small offsets between X-ray gas and mass centers from weak lensing, the lack of large amounts of gas between the two sub-clusters indicates that Abell 2465 is in a pre-merger state. A dynamical model that is consistent with the observed cluster data, based on the FLASH program and the radial infall model, is constructed, where the subclusters currently separated by ~1.2Mpc are approaching each other at ~2000km/s and will meet in ~0.4Gyr.

We investigated the gas properties of the galaxy group UGC03957 up to $1.4R_{200} \approx 1.4$Mpc in four azimuthal directions with the Suzaku satellit and performed a spectral analysis of five Suzaku observations of UGC03957 with 138ks good exposure time in total and five Chandra snapshot observations for point source detection. We found a temperature drop of a factor of ${\sim} 3$ from the center to the outskirts that is consistent with previous results for galaxy clusters. The metal abundance profile shows a flat behavior towards large radii, which is a hint for galactic winds as the primary ICM enrichment process. The entropy profile is consistent with numerical simulations after applying a gas mass fraction correction. Previous analyses for clusters and groups often showed an entropy flattening or even a drop around ${\sim} R_{200}$, which can be an indication of clumping or non-equilibrium states in the outskirts. Such entropy behavior is absent in UGC03957. The gas mass fraction is well below the cosmic mean but rises above this value beyond $R_{200}$, which could be a hint for deviations from hydrostatic equilibrium at these large radii. We determined the relative number of different supernovae types and found that the abundance pattern can be described by a relative contribution of 80% - 100% of core-collapse supernovae. This result is in agreement with previous measurements for galaxy groups.

Using the EPIC and RGS data from a deep (~200 ks) XMM-Newton observation, we investigate the temperature structure (kT and sigma_T ) and the abundances of 9 elements (O, Ne, Mg, Si, S, Ar, Ca, Fe and Ni) of the intra-cluster medium (ICM) in the nearby (z=0.046) cool-core galaxy cluster Abell 4059. Next to a deep analysis of the cluster core, a careful modelling of the EPIC background allows us to build radial profiles up to 12' (~650 kpc) from the core. Probably because of projection effects, the temperature ICM is found not to be in single phase, even in the outer parts of the cluster. The abundances of Ne, Si, S, Ar, Ca and Fe, but also O are peaked towards the core. Fe and O are still significantly detected in the outermost annuli; suggesting that the enrichment by both type Ia and core-collapse SNe started in the early stages of the cluster formation. However, the particularly high Ca/Fe ratio that we find in the core is not well reproduced by the standard SNe yield models. Finally, 2-D maps of temperature and Fe abundance are presented and confirm the existence of a denser, colder, and Fe-rich ridge southwest of the core, previously observed by Chandra. The origin of this asymmetry in the hot gas of the cluster core is still unclear, but might be explained by a past intense ram-pressure stripping event near the central cD galaxy.

Cosmological constraints from clusters rely on accurate gravitational mass estimates, which strongly depend on cluster gas temperature measurements. Therefore, systematic calibration differences may result in biased, instrument-dependent cosmological constraints. This is of special interest in the light of the tension between the Planck results of the primary temperature anisotropies of the CMB and Sunyaev-Zel'dovich plus X-ray cluster counts analyses. We quantify in detail the systematics and uncertainties of the cross-calibration of the effective area between five X-ray instruments, EPIC-MOS1/MOS2/PN onboard XMM-Newton and ACIS-I/S onboard Chandra, and the influence on temperature measurements. Furthermore, we assess the impact of the cross calibration uncertainties on cosmology. Using the HIFLUGCS sample, consisting of the 64 X-ray brightest galaxy clusters, we constrain the ICM temperatures through spectral fitting in the same, mostly isothermal, regions and compare them. Our work is an extension to a previous one using X-ray clusters by the IACHEC. Performing spectral fitting in the full energy band we find that best-fit temperatures determined with XMM-Newton/EPIC are significantly lower than Chandra/ACIS temperatures. We demonstrate that effects like multitemperature structure and different relative sensitivities of the instruments at certain energy bands cannot explain the observed differences. We conclude that using XMM-Newton/EPIC, instead of Chandra/ACIS to derive full energy band temperature profiles for cluster mass determination results in an 8% shift towards lower OmegaM values and <1% shift towards higher sigma8 values in a cosmological analysis of a complete sample of galaxy clusters. Such a shift is insufficient to significantly alleviate the tension between Planck CMB anisotropies and SZ plus XMM-Newton cosmological constraints.

The galaxy cluster Abell 2256 hosts one of the most intriguing examples in the class of radio relics. It has been found that this radio relic has a rather flat integrated spectrum at low frequencies that would imply an injection spectral index for the electrons that is inconsistent with the flattest allowed by the test particle diffusive shock acceleration (DSA). We performed new high-frequency observations at 2273, 2640, and 4850 MHz. Combining these new observations with images available in the literature, we constrain the radio integrated spectrum of the radio relic in Abell 2256 over the widest sampled frequency range collected so far for this class of objects (63 -10450 MHz). Moreover, we used X-ray observations of the cluster to check the temperature structure in the regions around the radio relic. We find that the relic keeps an unusually flat behavior up to high frequencies. Although the relic integrated spectrum between 63 and 10450 MHz is not inconsistent with a single power law with $\alpha_{63}^{10450}= 0.92\pm 0.02$, we find hints of a steepening at frequencies > 1400 MHz. The two frequency ranges 63-1369 MHz and 1369-10450 MHz are, indeed, best represented by two different power laws, with $\alpha_{63}^{1369}= 0.85\pm 0.01$ and $\alpha_{1369}^{10450}= 1.00\pm 0.02$. This broken power law would require special conditions to be explained in terms of test-particle DSA, e.g., non-stationarity of the spectrum and/or non-stationarity of the shock. On the other hand, the single power law would make of this relic the one with the flattest integrated spectrum known so far, even flatter than what allowed in the test-particle approach to DSA. We find a rather low temperature ratio of $T_2/T_1 \sim 1.7$ across the G region of the radio relic and no temperature jump across the H region.

We used X-ray 2D spectrally resolved maps to resolve structure in temperature and metal abundance. To perform stellar population analysis we applied the spectral fitting technique with STARLIGHT to the optical spectrum of the central galaxy. We simulated the chemical evolution of the central galaxy. While the temperature, pseudo-pressure, and pseudo-entropy maps showed no inhomogeneities, the iron spatial distribution shows a filamentary structure in the core of this group, which is spatially correlated with the central galaxy, suggesting a connection between the two. The analysis of the optical spectrum of the central galaxy showed no contribution by any recent AGN activity. Using the star formation history as input to chemical evolution models, we predicted the iron and oxygen mass released by supernovae (SNe) winds in the central galaxy up to the present time. Comparing the predicted amount of mass released by the NGC4325 galaxy to the ones derived through X-ray analysis we conclude that the winds from the central galaxy alone play a minor role in the IGM metal enrichment of this group inside r2500. The SNe winds are responsible for no more than 3% of it and of the iron mass and 21% of the oxygen mass enclosed within r2500. Our results suggest that oxygen has been produced in the early stages of the group formation, becoming well mixed and leading to an almost flat profile. Instead, the iron distribution is centrally peaked, indicating that this element is still being added to the IGM specifically in the core by the SNIa. A possible scenario to explain the elongated iron-rich structure in the core of the NGC4325 is a past AGN activity, in which our results suggest an episode older than ~10^7-10^8 yrs and younger than 5x10^8.

Upcoming X-ray surveys like eROSITA require precise calibration between X-ray observables and mass down to the low mass regime to set tight constraints on the fundamental cosmological parameters. Since an individual mass measurement is only possible for a relatively small number of objects it is crucial to have robust and well understood scaling relations that relate the total mass to easily observable quantities. The main goal of this work is to constrain the galaxy group scaling relations corrected for selection effects, and to quantify the influence of non-gravitational physics at the low-mass regime. We analyzed XMM-Newton observations for a complete sample of galaxy groups selected from the ROSAT All-Sky Survey and we compared the derived scaling properties with a galaxy cluster sample. To investigate the role played by the different non-gravitational processes we then compared the observational data with the predictions of hydrodynamical simulations. After applying the correction for selection effects (e.g. Malmquist bias) the L-M relation is steeper than the observed one. Its slope is also steeper than the value obtained by using the more massive systems of the HIFLUGCS sample. This behaviour can be explained by a gradual change of the true L-M relation which should be taken into account when converting the observational parameters into masses. The other observed scaling relations (not corrected for selection biases) do not show any break although the comparison with the simulations suggests that feedback processes play an important role in the formation and evolution of galaxy groups. Thanks to our master sample of 82 objects spanning two order of magnitude in mass we tightly constrain the dependence of the gas mass fraction on the total mass, finding almost a factor of two difference between groups and clusters.

We aim to investigate the bolometric $L_{\mathrm{X}}-T$ relation for galaxy groups, and study the impact of gas cooling, feedback from supermassive black holes, and selection effects on it. With a sample of 26 galaxy groups we obtained the best fit $L_{\mathrm{X}}-T$ relation for five different cases depending on the ICM core properties and central AGN radio emission, and determined the slopes, normalisations, intrinsic and statistical scatters for both temperature and luminosity. Simulations were undertaken to correct for selection effects (e.g. Malmquist bias) and the bias corrected relations for groups and clusters were compared. The slope of the bias corrected $L_{\mathrm{X}}-T$ relation is marginally steeper but consistent with clusters ($\sim 3$). Groups with a central cooling time less than 1 Gyr (SCC groups) show indications of having the steepest slope and the highest normalisation. For the groups, the bias corrected intrinsic scatter in $L_{\mathrm{X}}$ is larger than the observed scatter for most cases, which is reported here for the first time. Lastly, we see indications that the groups with an extended central radio source have a much steeper slope than those groups which have a CRS with only core emission. Additionally, we also see indications that the more powerful radio AGN are preferentially located in NSCC groups rather than SCC groups.

We present a simulation setup for studying the dynamical and chemical evolution of the intracluster medium (ICM) and analyze a sample of 12 galaxy clusters that are diverse both kinetically (pre-merger, merging, virialized) and in total mass (M vir = 1.17 x 10^14 - 1.06 x 10^15 M). We analyzed the metal mass fraction in the ICM as a function of redshift and discuss radial trends as well as projected 2D metallicity maps. The setup combines high mass resolution N-body simulations with the semi-analytical galaxy formation model Galacticus for consistent treatment of the subgrid physics (such as galactic winds and ram-pressure stripping) in the cosmological hydrodynamical simulations. The interface between Galacticus and the hydro simulation of the ICM with FLASH is discussed with respect to observations of star formation rate histories, radial star formation trends in galaxy clusters, and the metallicity at different redshifts. As a test for the robustness of the wind model, we compare three prescriptions from different approaches. For the wind model directly taken from Galacticus, we find mean ICM metallicities between 0.2 - 0.8Z within the inner 1Mpc at z = 0. The main contribution to the metal mass fraction comes from galactic winds. The outflows are efficiently mixed in the ICM, leading to a steady homogenization of metallicities until ram-pressure stripping becomes effective at low redshifts. We find a very peculiar and yet common drop in metal mass fractions within the inner ~200kpc of the cool cores, which is due to a combination of wind suppression by outer pressure within our model and a lack of mixing after the formation of these dense regions.

The nature of dark energy is imprinted in the large-scale structure of the Universe and thus in the mass and redshift distribution of galaxy clusters. The upcoming eROSITA mission will exploit this method of probing dark energy by detecting roughly 100,000 clusters of galaxies in X-rays. For a precise cosmological analysis the various galaxy cluster properties need to be measured with high precision and accuracy. To predict these characteristics of eROSITA galaxy clusters and to optimise optical follow-up observations, we estimate the precision and the accuracy with which eROSITA will be able to determine galaxy cluster temperatures and redshifts from X-ray spectra. Additionally, we present the total number of clusters for which these two properties will be available from the eROSITA survey directly. During its four years of all-sky surveys, eROSITA will determine cluster temperatures with relative uncertainties of Delta(T)/T<10% at the 68%-confidence level for clusters up to redshifts of z~0.16 which corresponds to ~1,670 new clusters with precise properties. Redshift information itself will become available with a precision of Delta(z)/(1+z)<10% for clusters up to z~0.45. Additionally, we estimate how the number of clusters with precise properties increases with a deepening of the exposure. Furthermore, the biases in the best-fit temperatures as well as in the estimated uncertainties are quantified and shown to be negligible in the relevant parameter range in general. For the remaining parameter sets, we provide correction functions and factors. The eROSITA survey will increase the number of galaxy clusters with precise temperature measurements by a factor of 5-10. Thus the instrument presents itself as a powerful tool for the determination of tight constraints on the cosmological parameters.

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

Major astrophysical questions related to the formation and evolution of structures, and more specifically of galaxy groups and clusters, will still be open in the coming decade and beyond: what is the interplay of galaxy, supermassive black hole, and intergalactic gas evolution in the most massive objects in the Universe - galaxy groups and clusters? What are the processes driving the evolution of chemical enrichment of the hot diffuse gas in large-scale structures? How and when did the first galaxy groups in the Universe, massive enough to bind more than 10^7 K gas, form? Focussing on the period when groups and clusters assembled (0.5<z<2.5), we show that, due to the continuum and line emission of this hot intergalactic gas at X-ray wavelengths, Athena+, combining high sensitivity with excellent spectral and spatial resolution, will deliver breakthrough observations in view of the aforementioned issues. Indeed, the physical and chemical properties of the hot intra-cluster gas, and their evolution across time, are a key to understand the co-evolution of galaxy and supermassive black hole within their environments.

Well-calibrated scaling relations between the observable properties and the total masses of clusters of galaxies are important for understanding the physical processes that give rise to these relations. They are also a critical ingredient for studies that aim to constrain cosmological parameters using galaxy clusters. For this reason much effort has been spent during the last decade to better understand and interpret relations of the properties of the intra-cluster medium. Improved X-ray data have expanded the mass range down to galaxy groups, whereas SZ surveys have openened a new observational window on the intracluster medium. In addition,continued progress in the performance of cosmological simulations has allowed a better understanding of the physical processes and selection effects affecting the observed scaling relations. Here we review the recent literature on various scaling relations, focussing on the latest observational measurements and the progress in our understanding of the deviations from self similarity.

We review the methods adopted to reconstruct the mass profiles in X-ray luminous galaxy clusters. We discuss the limitations and the biases affecting these measurements and how these mass profiles can be used as cosmological proxies.

Until recently, only about 10% of the total intracluster gas volume had been studied with high accuracy, leaving a vast region essentially unexplored. This is now changing and a wide area of hot gas physics and chemistry awaits discovery in galaxy cluster outskirts. Also, robust large-scale total mass profiles and maps are within reach. First observational and theoretical results in this emerging field have been achieved in recent years with sometimes surprising findings. Here, we summarize and illustrate the relevant underlying physical and chemical processes and review the recent progress in X-ray, Sunyaev--Zel'dovich, and weak gravitational lensing observations of cluster outskirts, including also brief discussions of technical challenges and possible future improvements.

We investigate the origin and nature of the multiple sloshing cold fronts in the core of Abell 496 by direct comparison between observations and dedicated hydrodynamical simulations. Our simulations model a minor merger with a 4{\times}10^13M{\circ} subcluster crossing A496 from the south-west to the north-north-east, passing the cluster core in the south-east at a pericentre distance 100 to a few 100 kpc about 0.6 to 0.8 Gyr ago. The gas sloshing triggered by the merger can reproduce almost all observed features, e.g. the characteristic spiral-like brightness residual distribution in the cluster centre and its asymmetry out to 500 kpc, also the positions of and contrasts across the cold fronts. If the subcluster passes close (100 kpc) to the cluster core, the resulting shear flows are strong enough to trigger Kelvin-Helmholtz instabilities that in projection resemble the peculiar kinks in the cold fronts of Abell 496. Finally, we show that sloshing does not lead to a significant modification of the global ICM profiles but a mild oscillation around the initial profiles.

The hot gas that fills the space between galaxies in clusters is rich in metals. In their large potential wells, galaxy clusters accumulate metals over the whole cluster history and hence they retain important information on cluster formation and evolution. We use a sample of 5 cool core clusters to study the distribution of metals in the ICM. We investigate whether the X-ray observations yield good estimates for the metal mass and whether the heavy elements abundances are consistent with a certain relative fraction of SN Ia to SNCC. We derive detailed metallicity maps of the clusters from XMM - Newton observations and we use them as a measure for the metal mass in the ICM. We determine radial profiles for several elements and using population synthesis and chemical enrichment models, we study the agreement between the measured abundances and the theoretical yields. We show that even in relaxed clusters the distribution of metals show a lot of inhomogeneities. Using metal maps usually gives a metal mass 10-30% higher than the metal mass computed using a single extraction region, hence it is expected that most previous metal mass determination have underestimated metal mass. The abundance ratio of {\alpha}-elements to Fe, even in the central parts of clusters, are consistent with an enrichment due to the combination of SN Ia and SNCC.

We use XMM-Newton data of the merging cluster Abell 3667 to analyze its metallicity distribution. A detailed abundance map of the central 1.1x1.1 Mpc region indicates that metals are inhomogeneously distributed in the cluster showing a non-uniform and very complex metal pattern. The highest peak in the map corresponds to a cold region, slightly offset South of the X-ray center. This could be interpreted as stripped gas due to a merger between a group moving from NW towards the SE and the main cluster. We note several clumps of high metallicity also in the opposite direction with respect to the X-ray peak. Furthermore we determined abundances for 5 elements (O, Si, S, Ar, Fe) in four different regions of the cluster. Comparisons between these observed abundances and theoretical supernovae yields allow to get constraints on the relative number of SN Ia and II contributing to the enrichment of the intra-cluster medium. To reproduce the observed abundances of the best determined elements (Fe, O and Si) in a region of 7 arcmin around the X-ray center, 65-80% of SN II are needed. The comparison between the metal map, a galaxy density map obtained using 550 spectroscopically confirmed cluster members and our simulations suggest a recent merger between the main cluster and the group in the SE.

(Abriged) We present the analysis of the baryonic content of 52 X-ray luminous galaxy clusters observed with Chandra in the redshift range 0.3-1.273. We use the deprojected X-ray surface brightness profiles and the measured values of the gas temperature to recover the gas and total mass profiles. By assuming that galaxy clusters are representative of the cosmic baryon budget, the distribution of the cluster baryon fraction in the hottest (T> 4 keV) systems as a function of redshift is used to constrain the cosmological parameters. We discuss how our constraints are affected by several systematics, namely the isothermality, the assumed baryon fraction in stars, the depletion parameter and the sample selection. By using only the cluster baryon fraction as a proxy for the cosmological parameters, we obtain that Omega is very well constrained at the value of 0.35 with a relative statistical uncertainty of 11% (1 sigma level; w=-1) and a further systematic error of about (-6,+7)%. On the other hand, constraints on Lambda (without the prior of flat geometry) and w (using the prior of flat geometry) are definitely weaker due to the presence of larger statistical and systematic uncertainties (of the order of 40 per cent on Lambda and larger than 50 per cent on w). If the WMAP 5-year best-fit results are assumed to fix the cosmological parameters, we limit the contributions expected from non-thermal pressure support and ICM clumpiness to be lower than about 10 per cent, leaving also room to accommodate baryons not accounted for either in the X-ray emitting plasma or in stars of the order of 18 per cent of the total cluster baryon budget.