We present the first metal abundance profiles for a representative sample of massive clusters. Our measures extend to $R_{500}$ and are corrected for a systematic error plaguing previous outskirt estimates. Our profiles flatten out at large radii, admittedly not a new result, however the radial range and representative nature of our sample extends its import well beyond previous findings. We find no evidence of segregation between cool-core and non-cool-core systems beyond $\sim 0.3 R_{500}$, implying that, as was found for thermodynamic properties (Ghirardini et al, 2019), the physical state of the core does not affect global cluster properties. Our mean abundance within $R_{500}$ shows a very modest scatter, $< $15%, suggesting the enrichment process must be quite similar in all these massive systems. This is a new finding and has significant implications on feedback processes. Together with results from thermodynamic properties presented in a previous X-COP paper, it affords a coherent picture where feedback effects do not vary significantly from one system to another. By combing ICM with stellar measurements we have found the amount of Fe diffused in the ICM to be about ten times higher than that locked in stars. Although our estimates suggest, with some strength, that the measured iron mass in clusters is well in excess of the predicted one, systematic errors prevent us from making a definitive statement. Further advancements will only be possible when systematic uncertainties, principally those associated to stellar masses, both within and beyond $R_{500}$, can be reduced.

Using the largest mosaic of Hubble Space Telescope images around a galaxy cluster, we map the distribution of dark matter throughout a $\sim$$6\times6$ Mpc$^2$ area centred on the cluster MS 0451-03 ($z=0.54$, $M_{200}=1.65\times10^{15} \rm{M}_\odot$). Our joint strong- and weak-lensing analysis shows three possible filaments extending from the cluster, encompassing six group-scale substructures. The dark-matter distribution in the cluster core is elongated, consists of two distinct components, and is characterized by a concentration parameter of $c_{200}=3.79\pm0.36$. By contrast, XMM-Newton observations show the gas distribution to be more spherical, with excess entropy near the core, and a lower concentration of $c_{200}=2.35^{+0.89}_{-0.70}$ (assuming hydrostatic equilibrium). Such a configuration is predicted in simulations of major mergers 2-7Gyr after the first core passage, when the two dark-matter halos approach second turnaround, and before their gas has relaxed. This post-merger scenario finds further support in optical spectroscopy of the cluster's member galaxies, which shows that star formation was abruptly quenched 5 Gyr ago. MS 0451-03 will be an ideal target for future studies of the growth of structure along filaments, star-formation processes after a major merger, and the late-stage evolution of cluster collisions.

We present the results of deep Chandra and XMM-Newton observations of a complex merging galaxy cluster Abell 2256 (A2256) that hosts a spectacular radio relic. The temperature and metallicity maps show clear evidence of a merger between the western subcluster and the primary cluster. We detect five X-ray surface brightness edges. Three of them near the cluster center are cold fronts (CFs): CF1 is associated with the infalling subcluster; CF2 is located in the east of the primary cluster and CF3 is to the west of the primary cluster core. The other two edges at cluster outskirts are shock fronts (SFs): SF1 near the radio relic in the NW has Mach numbers derived from the temperature and the density jumps respectively of $M_T=1.62\pm0.12$ and $M_\rho=1.23\pm0.06$; SF2 in the SE has $M_T=1.54\pm0.05$ and $M_\rho=1.16\pm0.13$. In the region of the radio relic there is no evidence for the correlation between X-ray and radio substructures, from which we estimate an upper limit for the inverse Compton emission, and therefore set a lower limit on the magnetic field ($\sim$ 450 kpc from primary cluster center) of $B>1.0\ \mu$G for a single power-law electron spectrum or $B>0.4\ \mu$G for a broken power-law electron spectrum. We propose a merger scenario including a primary cluster, a subcluster, and a group. Our merger scenario accounts for the X-ray edges, diffuse radio features, and galaxy kinematics, as well as projection effects.

In this paper, we present resolved observations of the Sunyaev-Zel'dovich (SZ) effect, obtained with the NIKA2 camera, towards the cluster of galaxies XLSSC102, a relatively low-mass system ($M_{500} \sim 2 \times 10^{14}$ M$_{\odot}$) at $z = 0.97$ detected from the XXL survey. We combine NIKA2 SZ data, XMM-Newton X-ray data, and Megacam optical data to explore, respectively, the spatial distribution of the gas electron pressure, the gas density, and the galaxies themselves. We find significant offsets between the X-ray peak, the SZ peak, the brightest cluster galaxy, and the peak of galaxy density. Additionally, the galaxy distribution and the gas present elongated morphologies. This is interpreted as the sign of a recent major merging event, which induced a local boost of the gas pressure towards the north of XLSSC102 and stripped the gas out of the galaxy group. The NIKA2 data are also combined with XXL data to construct the thermodynamic profiles of XLSSC102, obtaining relatively tight constraints up to about $\sim r_{500}$, and revealing properties that are typical of disturbed systems. We also explore the impact of the cluster centre definition and the implication of local pressure substructure on the recovered profiles. Finally, we derive the global properties of XLSSC102 and compare them to those of high-mass-and-low-redshift systems, finding no strong evidence for non-standard evolution. We also use scaling relations to obtain alternative mass estimates from our profiles. The variation between these different mass estimates reflects the difficulty to accurately measure the mass of low-mass clusters at z$\sim$1, especially with low signal-to-noise ratio (S/N) data and for a disturbed system. However, it also highlights the strength of resolved SZ observations alone and in combination with survey-like X-ray data.

The thermodynamic properties of the hot plasma in galaxy clusters retains information on the processes leading to the formation and evolution of the gas in their deep, dark matter potential wells. These processes are dictated not only by gravity but also by gas physics, e.g. AGN feedback and turbulence. In this work, we study the thermodynamic properties, e.g. density, temperature, pressure, and entropy, of the most massive and the most distant ($z > 1.2$) SPT-selected clusters, and compare them with those of the nearby clusters ($z<0.1$) to constrain their evolution as a function of time and radius. We find that thermodynamic properties in the outskirts of high redshift clusters are remarkably similar to the low redshift clusters, and their evolution follows the prediction of the self-similar model. Their intrinsic scatter is larger, indicating that the physical properties that lead to the formation and virialization of cluster outskirts show evolving variance. On the other hand, thermodynamic properties in the cluster cores deviates significantly from self-similarity indicating that the processes that regulate the core are already in place in these very high redshift clusters. This result is supported by the unevolving physical scatter of all thermodynamic quantities in cluster cores.

The degree of turbulent pressure support by residual gas motions in galaxy clusters is not well known. Mass modelling of combined X-ray and Sunyaev Zel'dovich observations provides an estimate of turbulent pressure support in the outer regions of several galaxy clusters. Here, we test two different filtering techniques to disentangle bulk from turbulent motions in non-radiative high-resolution cosmological simulations of galaxy clusters using the cosmological hydro code ENZO. We find that the radial behavior of the ratio of non-thermal pressure to total gas pressure as a function of cluster-centric distance can be described by a simple polynomial function. The typical non-thermal pressure support in the centre of clusters is $\sim$5%, increasing to $\sim$15% in the outskirts, in line with the pressure excess found in recent X-ray observations. While the complex dynamics of the ICM makes it impossible to reconstruct a simple correlation between turbulent motions and hydrostatic bias, we find that a relation between them can be established using the median properties of a sample of objects. Moreover, we estimate the contribution of radial accelerations to the non-thermal pressure support and conclude that it decreases moving outwards from 40% (in the core) to 15% (in the cluster's outskirts). Adding this contribution to one provided by turbulence, we show that it might account for the entire observed hydrostatic bias in the innermost regions of the clusters, and for less than 80% of it at $r > 0.8 r_{200, m}$.

An accurate modelling of baryonic feedback effects is required to exploit the full potential of future weak-lensing surveys such as Euclid or LSST. In this second paper in a series of two, we combine Euclid-like mock data of the cosmic shear power spectrum with an eROSITA X-ray mock of the cluster gas fraction to run a combined likelihood analysis including both cosmological and baryonic parameters. Following the first paper of this series, the baryonic effects (based on the baryonic correction model of Schneider et al. 2019) are included in both the tomographic power spectrum and the covariance matrix. However, this time we assume the more realistic case of a $\Lambda$CDM cosmology with massive neutrinos, and we consider several extensions of the currently favoured cosmological model. For the standard $\Lambda$CDM case, we show that including X-ray data reduces the uncertainties on the sum of the neutrino mass by $\sim30$ percent, while there is only a mild improvement on other parameters such as $\Omega_m$ and $\sigma_8$. As extensions of $\Lambda$CDM, we consider the cases of a dynamical dark energy model (wCDM), a $f(R)$ gravity model (fRCDM), and a mixed dark matter model ($\Lambda$MDM) with both a cold and a warm/hot dark matter component. We find that combining weak lensing with X-ray data only leads to a mild improvement of the constraints on the additional parameters of wCDM, while the improvement is more substantial for both fRCDM and $\Lambda$MDM. Ignoring baryonic effects in the analysis pipeline leads to significant false-detections of either phantom dark energy or a light subdominant dark matter component. Overall we conclude that for all cosmologies considered, a general parametrisation of baryonic effects is both necessary and sufficient to obtain tight constraints on cosmological parameters.

We present a weak-lensing analysis of X-ray galaxy groups and clusters selected from the XMM-XXL survey using the first-year data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program. Our joint weak-lensing and X-ray analysis focuses on 136 spectroscopically confirmed X-ray-selected systems at 0.031 < z < 1.033 detected in the 25sqdeg XXL-N region. We characterize the mass distributions of individual clusters and establish the concentration-mass (c-M) relation for the XXL sample, by accounting for selection bias and statistical effects, and marginalizing over the remaining mass calibration uncertainty. We find the mass-trend parameter of the c-M relation to be \beta = -0.07 \pm 0.28 and the normalization to be c200 = 4.8 \pm 1.0 (stat) \pm 0.8 (syst) at M200=10^{14}Msun/h and z = 0.3. We find no statistical evidence for redshift evolution. Our weak-lensing results are in excellent agreement with dark-matter-only c-M relations calibrated for recent LCDM cosmologies. The level of intrinsic scatter in c200 is constrained as \sigma(\ln[c200]) < 24% (99.7% CL), which is smaller than predicted for the full population of LCDM halos. This is likely caused in part by the X-ray selection bias in terms of the relaxation state. We determine the temperature-mass (Tx-M500) relation for a subset of 105 XXL clusters that have both measured HSC lensing masses and X-ray temperatures. The resulting Tx-M500 relation is consistent with the self-similar prediction. Our Tx-M500 relation agrees with the XXL DR1 results at group scales, but has a slightly steeper mass trend, implying a smaller mass scale in the cluster regime. The overall offset in the Tx-M500 relation is at the $1.5\sigma$ level, corresponding to a mean mass offset of (34\pm 20)%. We also provide bias-corrected, weak-lensing-calibrated M200 and M500 mass estimates of individual XXL clusters based on their measured X-ray temperatures.

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

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.

One key ingredient in using galaxy clusters (GCs) as a precision cosmological probe in large X-ray surveys is to understand selection effects. The dependence of the X-ray emission on the square of the gas density leads to a predominant role of cool cores in the detection of GCs. The contribution of cool cores to the X-ray luminosity does not scale with GC mass and cosmology and therefore affects the use of X-ray GCs in producing cosmological constraints. One of the main science goals of the eROSITA mission is to constrain cosmology with a wide X-ray survey. We propose an eROSITA GC detection scheme that avoids the use of X-ray GC centers in detection. We calculate theoretical expectations and characterize the performance of this scheme by simulations. Performing realistic simulations of point sources (PSs) in survey mode we search for spatial scales where the extended signal is uncontaminated by the PS flux. We derive a combination of scales and thresholds, which result in a clean extended source catalog. We design the output of the GC detection which enables calibrating the core-excised luminosity using external mass measurements. We provide a way to incorporate the results of this calibration in the production of final core-excised luminosity. Similarly to other GC detection pipelines, we sample the flux - core radius detection space of our method and find many similarities with the pipeline used in the 400d survey. Both detection methods require large statistics on compact GCs, in order to reduce the contamination from PSs. The benefit of our pipeline consists in the sensitivity to the outer GC shapes, which are characterized by large core sizes with little GC to GC variation at a fixed total mass. GC detection through cluster outskirts improves the GC characterization using eROSITA survey data and is expected to yield well characterized GC catalogs having simple selection functions.

In the context of the BUFFALO (Beyond Ultra-deep Frontier Fields And Legacy Observations) survey, we present a new analysis of the merging galaxy cluster MACS\,J0416.1-2403 ($z = 0.397$) and its parallel field using the data collected by the Hubble Frontier Fields (HFF) campaign. In this work, we measure the surface mass density from a weak-lensing analysis, and characterise the overall matter distribution in both the cluster and parallel fields. The surface mass distribution derived for the parallel field shows clumpy overdensities connected by filament-like structures elongated in the direction of the cluster core. We also characterise the X-ray emission of the cluster, and compare it with the lensing mass distribution. We identify five substructures at the $>5\sigma$ level over the two fields, four of them being in the cluster one. Furthermore, three of them are located close to the edges of the field of view, and border issues can significantly hamper the determination of their physical parameters. Finally, we compare our results with the predicted subhalo distribution of one of the Hydrangea/C-EAGLE simulated cluster. Significant differences are obtained suggesting the simulated cluster is at a more advanced evolutionary state than MACS\,J0416.1-2403. Our results anticipate the upcoming BUFFALO observations that will link the two HFF fields, extending further the \emph{HST} coverage, and thus allowing a better characterisation of the reported substructures.

Scaling relations trace the formation and evolution of galaxy clusters. We exploited multi-wavelength surveys -- the XXL survey at \emph{XMM-Newton} in the X-ray band, and the Hyper Suprime-Cam Subaru Strategic Program for optical weak lensing -- to study an X-ray selected, complete sample of clusters and groups. The scalings of gas mass, temperature, and soft-band X-ray luminosity with the weak lensing mass show imprints of radiative cooling and AGN feedback in groups. From the multi-variate analysis, we found some evidence for steeper than self-similar slopes for gas mass ($\beta_{m_\text{g}|m}=1.73 \pm0.80$) and luminosity ($\beta_{l|m}=1.91\pm0.94$) and a nearly self-similar slope for the temperature ($\beta_{t|m}=0.78\pm0.43$). Intrinsic scatters of X-ray properties appear to be positively correlated at a fixed mass (median correlation factor $\rho_{X_1X_2|m}\sim0.34$) due to dynamical state and merger history of the halos. Positive correlations with the weak lensing mass (median correlation factor $\rho_{m_\text{wl}X|m}\sim0.35$) can be connected to triaxiality and orientation. Comparison of weak lensing and hydrostatic masses suggests a small role played by non-thermal pressure support ($9\pm17\%$).

We present the constraints on the helium abundance in 12 X-ray luminous galaxy clusters that have been mapped in their X-ray and Sunyaev-Zeldovich (SZ) signals out to $R_{200}$ for the XMM-Newton Cluster Outskirts Project (X-COP). The unprecedented precision available for the estimate of $H_0$ allows us to investigate how much the reconstructed X-ray and SZ signals are consistent with the expected ratio $x$ between helium and proton densities of 0.08-0.1. We find that a $H_0$ around 70 km/s/Mpc is preferred from our measurements, with lower values of $H_0$ as requested from the Planck collaboration (67 km/s/Mpc) requiring a 34% higher value of $x$. On the other hand, higher values of $H_0$, as obtained by measurements in the local universe, impose $x$, from the primordial nucleosynthesis calculations and current solar abundances, reduced by 37--44\%.

We demonstrate a novel technique for calibrating the energy scale of the XMM EPIC-pn detector, which allows us to measure bulk flows in the intracluster medium (ICM) of the Perseus and Coma clusters. The procedure uses the instrumental lines present in all observations, in particular, Cu-Ka. By studying their spatial and temporal variations, in addition to incorporating calibration observations, we refined the absolute energy scale to better than 150 km/s at the Fe-K line, a large improvement over the nominal accuracy of 550 km/s. We then mapped the bulk motions over much of the central 1200 and 800 kpc of Perseus and Coma, respectively, in spatial regions down to 65 and 140 kpc size. We cross-checked our procedure by comparing our measurements with those found in Perseus by Hitomi for an overlapping region, finding consistent results. For Perseus, there is a LoS velocity increase of 480+-210 km/s (1sigma) 250 kpc east of the nucleus. This region is associated with a cold front, providing direct evidence of the ICM sloshing in the potential well. Assuming the intrinsic distribution of bulk motions is Gaussian, its width is 214+-85 km/s, excluding systematics. Removing the sloshing region, this is reduced to 20-150 km/s, which is similar in magnitude to the Hitomi line width measurements in undisturbed regions. In Coma, the line-of-sight velocity of the ICM varies between the velocities of the two central galaxies. Maps of the gas velocity and metallicity provide clues about the merger history of the Coma, with material to the north and east of the cluster core having a velocity similar to NGC 4874, while that to the south and west has velocities close to NGC 4889. Our results highlight the difference between a merging system, such as Coma, where we observe a ~1000 km/s range in velocity, and a relatively relaxed system, such as Perseus, with much weaker bulk motions. [abridged]

Context. Scaling relations between cluster properties embody the formation and evolution of cosmic structure. Intrinsic scatters and correlations between X-ray properties are determined from merger history, baryonic processes, and dynamical state. Aims. We look for an unbiased measurement of the scatter covariance matrix between the three main X-ray observable quantities attainable in large X-ray surveys -- temperature, luminosity, and gas mass. This also gives us the cluster property with the lowest conditional intrinsic scatter at fixed mass. Methods. Intrinsic scatters and correlations can be measured under the assumption that the observable properties of the intra-cluster medium hosted in clusters are log-normally distributed around power-law scaling relations. The proposed method is self-consistent, based on minimal assumptions, and requires neither the external calibration by weak lensing, dynamical, or hydrostatic masses nor the knowledge of the mass completeness. Results. We analyzed the 100 brightest clusters detected in the XXL Survey and their X-ray properties measured within a fixed radius of 300 kpc. The gas mass is the less scattered proxy (~8%). The temperature (~20%) is intrinsically less scattered than the luminosity (~30%) but it is measured with a larger observational uncertainty. We found some evidence that gas mass, temperature and luminosity are positively correlated. Time-evolutions are in agreement with the self-similar scenario, but the luminosity-temperature and the gas mass-temperature relations are steeper. Conclusions. Positive correlations between X-ray properties can be determined by the dynamical state and the merger history of the halos. The slopes of the scaling relations are affected by radiative processes.

We present results from simultaneous modeling of high angular resolution GBT/MUSTANG-2 90 GHz Sunyaev-Zel'dovich effect (SZE) measurements and XMM-XXL X-ray images of three rich galaxy clusters selected from the HSC-SSP Survey. The combination of high angular resolution SZE and X-ray imaging enables a spatially resolved multi-component analysis, which is crucial to understand complex distributions of cluster gas properties. The targeted clusters have similar optical richnesses and redshifts, but exhibit different dynamical states in their member galaxy distributions; a single-peaked cluster, a double-peaked cluster, and a cluster belonging to a supercluster. A large-scale residual pattern in both regular Compton-parameter $y$ and X-ray surface brightness distributions is found in the single-peaked cluster, indicating a sloshing mode. The double-peaked cluster shows an X-ray remnant cool core between two SZE peaks associated with galaxy concentrations. The temperatures of the two peaks reach $\sim20-30$ keV in contrast to the cool core component of $\sim2$ keV, indicating a violent merger. The highest SZE signal for the supercluster is elongated along a direction perpendicular to the major axis of the X-ray core, suggesting a minor merger before core passage. The $S_X$ and $y$ distributions are thus perturbed at some level, regardless of the optical properties. We find that an integrated Compton $y$ parameter and a temperature for the major merger are significantly boosted from those expected by the weak-lensing mass and those for the other two clusters show no significant deviations, which is consistent with predictions of numerical simulations.

We carry out a comprehensive Bayesian correlation analysis between hot halos and direct masses of supermassive black holes (SMBHs), by retrieving the X-ray plasma properties (temperature, luminosity, density, pressure, masses) over galactic to cluster scales for 85 diverse systems. We find new key scalings, with the tightest relation being the $M_\bullet-T_{\rm x}$, followed by $M_\bullet-L_{\rm x}$. The tighter scatter (down to 0.2 dex) and stronger correlation coefficient of all the X-ray halo scalings compared with the optical counterparts (as the $M_\bullet-\sigma_{\rm e}$) suggest that plasma halos play a more central role than stars in tracing and growing SMBHs (especially those that are ultramassive). Moreover, $M_\bullet$ correlates better with the gas mass than dark matter mass. We show the important role of the environment, morphology, and relic galaxies/coronae, as well as the main departures from virialization/self-similarity via the optical/X-ray fundamental planes. We test the three major channels for SMBH growth: hot/Bondi-like models have inconsistent anti-correlation with X-ray halos and too low feeding; cosmological simulations find SMBH mergers as sub-dominant over most of the cosmic time and too rare to induce a central-limit-theorem effect; the scalings are consistent with chaotic cold accretion (CCA), the rain of matter condensing out of the turbulent X-ray halos that sustains a long-term self-regulated feedback loop. The new correlations are major observational constraints for models of SMBH feeding/feedback in galaxies, groups, and clusters (e.g., to test cosmological hydrodynamical simulations), and enable the study of SMBHs not only through X-rays, but also via the Sunyaev-Zel'dovich effect (Compton parameter), lensing (total masses), and cosmology (gas fractions).

The formation of cosmic structure culminates with the assembly of galaxy clusters, a process quite different from cluster to cluster. We present the study of the structure and dynamics of the Lyra complex formed of the two clusters RXC J1825.3+3026 and CIZA J1824.1+3029, very recently studied using both X-ray and radio data. This is the first analysis based on kinematics of member galaxies. New spectroscopic data for 285 galaxies were acquired at the Italian Telescopio Nazionale Galileo and used in combination with PanSTARRS photometry. The result of our member selection is a sample of 198 galaxies. For RXCJ1825 and CIZAJ1824 we report the redshifts, z=0.0645 and z=0.0708, the first estimates of velocity dispersion, sigmav=995 and 700 km/s, and dynamical mass, M200=1.1E15 and 4E14 Msun. The past assembly of RXCJ1825 is traced by the two dominant galaxies, both aligned with the major axis of the galaxy distribution along the East-West direction, and by a minor North-East substructure. We also detect a quite peculiar high velocity field in the South-West region of the Lyra complex. This feature is likely related to a high velocity, very luminous galaxy, suggested to be the central galaxy of a group in interaction with RXCJ1825 by very recent studies based on X-ray and radio data. The redshift of the whole Lyra complex is z=0.067. Assuming that the redshift difference between RXCJ1825 and CIZAJ1824 is due to the relative kinematics, the projected distance between the cluster centers is 1.3 Mpc and the los velocity difference is 1750 km/s. A dynamical analysis of the system shows that the two clusters are likely to be gravitationally bound, in a pre-merger phase, with CIZAJ1824 in front of RXCJ1825 and going toward it. Our results corroborate a picture where the Lyra region is the place of a very complex scenario of cluster assembly.

Nearby clusters of galaxies, z<0.1, are cosmic structures still under formation. Understanding the thermodynamic properties of merging clusters can provide crucial information on how they grow in the local universe. A detailed study of the intra-cluster medium (ICM) properties of un-relaxed systems is essential to understand the fate of in-falling structures and, more generally, the virialization process. We analyzed a mosaic of XMM-Newton observations (240 ks) of the Lyra system (z=0.067) that shows a complex dynamical state. We find the main cluster RXC J1825.3+3026 to be in a late merger phase, whereas its companion CIZA J1824.1+3029 is a relaxed cool-core cluster. We estimate a mass ratio of ~1:2 for the pair. No diffuse X-ray emission is found in the region between them, indicating that these clusters are in a pre-merger phase. We found evidence of a galaxy group infalling on RXC J1825.3+3026 in an advanced state of disruption. The Southern Galaxy, one of the brightest galaxies in the Lyra complex, was very likely at the center of the infalling group. This galaxy has a gaseous corona indicating that it was able to retain some of its gas after the ram-pressure stripping of the intra-group medium. In this scenario the diffuse emission excess observed southwest of RXC J1825.3+3026 could be due to gas once belonging to the group and/or to cluster ICM dislocated by the passage of the group. Finally, we identified three high-velocity galaxies aligned between RXC J1825.3+3026 and the SG, two of these showing evidence of gas stripped from them during infall. We estimate them to be currently falling onto the main cluster at an infall velocity of ~ 3000 km/s. Our study of the Lyra complex provides important clues about the processes presiding over the virialization of massive clusters in the local Universe.

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.

The imaging of galaxy clusters through the Sunyaev--Zel'dovich effect is a valuable tool to probe the thermal pressure of the intra-cluster gas, especially in the outermost regions where X-ray observations suffer from photon statistics. For the first time, we produce maps of the Comptonization parameter by applying a locally parametric algorithm for sparse component separation to the latest frequency maps released by Planck. The algorithm takes into account properties of real cluster data through the two-component modelling of the spectral energy density of thermal dust, and the masking of bright point sources. Its robustness has been improved in the low signal-to-noise regime, thanks to the implementation of a deconvolution of Planck beams in the chi-square minimisation of each wavelet coefficient. We applied this procedure to twelve low-redshift galaxy clusters detected by Planck with the highest signal-to-noise ratio, considered in the XMM Cluster Oustkirts Project (X-COP). Our images show the presence of anisotropic features, such as small-scale blobs and filamentary substructures that are located in the outskirts of a number of clusters in the sample. The significance of their detection is established via a bootstrap-based procedure we propose here for the first time. In particular, we present a qualitative comparison with X-ray data for two interesting systems, namely A2029 and RXCJ1825. Our results are in agreement with the features detected in the outskirts of the clusters in the two bands.

Diffuse radio emission associated with the intra-cluster medium (ICM) is observed in a number of merging galaxy clusters. It is currently believed that in mergers a fraction of the kinetic energy is channeled into non-thermal components, such as turbulence, cosmic rays and magnetic fields, that may lead to the formation of giant synchrotron sources in the ICM. Studying merging galaxy clusters in different evolutionary phases is fundamental to understanding the origin of radio emission in the ICM. We observed the nearby galaxy cluster pair RXC J1825.3+3026 ($z\sim0.065$) and CIZA J1824.1+3029 ($z\sim0.071$) at 120-168 MHz with the LOw Frequency ARray (LOFAR) and made use of a deep (240 ks) XMM-Newton dataset to study the non-thermal and thermal properties of the system. RXC J1825.3+3026 is in a complex dynamical state, with a primary on-going merger in the E-W direction and a secondary later stage merger with a group of galaxies in the SW, while CIZA J1824.1+3029 is dynamically relaxed. These two clusters are in a pre-merger phase. We report the discovery of a Mpc-scale radio halo with a low surface brightness extension in RXC J1825.3+3026 that follows the X-ray emission from the cluster center to the remnant of a galaxy group in the SW. This is among the least massive systems and the faintest giant radio halo known to date. Contrary to this, no diffuse radio emission is observed in CIZA J1824.1+3029 nor in the region between the pre-merger cluster pair. The power spectra of the X-ray surface brightness fluctuations of RXC J1825.3+3026 and CIZA J1824.1+3029 are in agreement with the findings for clusters exhibiting a radio halo and the ones where no radio emission has been detected, respectively. We provide quantitative support to the idea that cluster mergers play a crucial role in the generation of non-thermal components in the ICM.

The majority of the ordinary matter in the local Universe has been heated by strong structure formation shocks and resides in a largely unexplored hot, diffuse, X-ray emitting plasma that permeates the halos of galaxies, galaxy groups and clusters, and the cosmic web. We propose a next-generation "Cosmic Web Explorer" that will permit a complete and exhaustive understanding of these unseen baryons. This will be the first mission capable to reach the accretion shocks located several times farther than the virial radii of galaxy clusters, and reveal the out-of-equilibrium parts of the intra-cluster medium which are live witnesses to the physics of cosmic accretion. It will also enable a view of the thermal, kinematic, and chemical composition of the circum-galactic medium in galaxies with masses similar to the Milky Way, at the same level of detail that $Athena$ will unravel for the virialized regions of massive galaxy clusters, delivering a transformative understanding of the evolution of those galaxies in which most of the stars and metals in the Universe were formed. Finally, the proposed X-ray satellite will connect the dots of the large-scale structure by mapping, at high spectral resolution, as much as 100% of the diffuse gas hotter than $10^6$ K that fills the filaments of the cosmic web at low redshifts, down to an over-density of 1, both in emission and in absorption against the ubiquitous cosmic X-ray background, surveying up to half of the sky over 5 years in orbit. This requires a large effective area (~10 m$^2$ at 1 keV) over a large field of view (~1 deg$^2$), a megapixel cryogenic microcalorimeter array providing integral field spectroscopy with a resolving power $E/\Delta E$ = 2000 at 0.6 keV and a spatial resolution of 5" in the soft X-ray band, and a Low-Earth orbit that ensures a stable instrumental background and high sensitivity to faint, extended emission.

In the framework of the hierarchical model the intra-cluster medium properties of galaxy clusters are tightly linked to structure formation, which makes X-ray surveys well suited for cosmological studies. To constrain cosmological parameters accurately by use of galaxy clusters in X-ray surveys, a better understanding of selection effects related to the detection method is needed. We aim at a better understanding of galaxy cluster morphologies to include corrections between the different core types and covariances with X-ray luminosities in selection functions. We stress the morphological deviations between a newly described surface brightness (SB) profile characterization and a commonly used single $\beta$-model. We investigate a novel approach to describe SB profiles, where the excess cool-core emission in the galaxy cluster centres is modelled using wavelet decomposition. Morphological parameters and the residuals are compared to classical single $\beta$-models. Using single $\beta$-models to describe the ensemble of overall SB profiles leads on average to a non-zero bias ($0.032 \pm 0.003$) in the outer part of the clusters, i.e. a $\sim 3\%$ systematic difference in the SB at large radii. In addition $\beta$-models show a general trend towards underestimating the flux in the outskirts for smaller core radii. Fixing the $\beta$ parameter to $2/3$ doubles the bias and increases the residuals from a single $\beta$-model up to more than $40\%$. Modelling the core region in the fitting procedure reduces the impact of these two effects significantly. We find a positive scaling between shape parameters and temperature, as well as a negative correlation ($\sim-0.4$) between extent and luminosity. Our non-parametric analysis of the self-similar scaled emission measure profiles indicates no systematic core-type differences of median profiles in the galaxy clusters outskirts.

We present the results of a joint XMM-Newton and NuSTAR observation (200 ks) of the galaxy cluster Abell 523 at $z=0.104$. The peculiar morphology of the cluster radio halo and its outlier position in the radio power P(1.4 GHz) - X-ray luminosity plane make it an ideal candidate for the study of radio-X-ray correlations and for the search of inverse Compton (IC) emission. We constructed thermodynamic maps derived from the XMM observations to describe in detail the physical and dynamical state of the ICM. We performed a point-to-point comparison in terms of surface brightness between the X-ray and radio emissions, to quantify their morphological discrepancies. Making use of NuSTAR's hard X-ray focusing capability, we looked for IC emission both globally and locally, after modeling the purely thermal component with a multi-temperature description. The thermodynamic maps suggest the presence of a secondary merging process that could be responsible for the peculiar radio halo morphology. This hypothesis is supported by the comparison between the X-ray and radio surface brightnesses, which shows a broad intrinsic scatter and a series of outliers. The global NuSTAR spectrum can be explained by purely thermal gas emission, and there is no convincing evidence that an IC component is needed. The $3\sigma$ upper limit on the IC flux in the 20-80 keV band is in the range $\left[2.2 - 4.0\right] \times 10^{-13} \, \mathrm{erg} \, \mathrm{s}^{-1} \, \mathrm{cm}^{-2}$, implying a lower limit on the magnetic field strength in the range $B > [0.23 - 0.31] \, \mu G$. Locally, we looked for IC emission in the central region of the cluster radio halo finding a $3\sigma$ upper limit on the 20-80 keV non-thermal flux of $3.17 \times 10^{-14} \, \mathrm{erg} \, \mathrm{s}^{-1} \, \mathrm{cm}^{-2}$, corresponding to a lower limit on the magnetic field strength of $B \gtrsim 0.81 \, \mu G$.

In this work, we investigate the relation between the radially-resolved thermodynamic quantities of the intracluster medium in the X-COP cluster sample, aiming to assess the stratification properties of the ICM. We model the relations between radius, gas temperature, density and pressure using a combination of power-laws, also evaluating the intrinsic scatter in these relations. We show that the gas pressure is remarkably well correlated to the density, with very small scatter. Also, the temperature correlates with gas density with similar scatter. The slopes of these relations have values that show a clear transition from the inner cluster regions to the outskirts. This transition occurs at the radius $r_t = 0.19(\pm0.04)R_{500}$ and electron density $n_t = (1.91\pm0.21)\cdot10^{-3} cm^{-3} E^2 (z)$. We find that above 0.2 $R_{500}$ the radial thermodynamic profiles are accurately reproduced by a well defined and physically motivated framework, where the dark matter follows the NFW potential and the gas is represented by a polytropic equation of state. By modeling the gas temperature dependence upon both the gas density and radius, we propose a new method to reconstruct the hydrostatic mass profile based only on the quite inexpensive measurement of the gas density profile.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

Clusters of galaxies, the largest collapsed structures in the Universe, are located at the intersection of extended filaments of baryons and dark matter. Cosmological accretion onto clusters through large scale filaments adds material at cluster outskirts. Kinetic energy in the form of bulk motions and turbulence due to this accretion provides a form of pressure support against gravity, supplemental to thermal pressure. Significant amount of non-thermal pressure support could bias cluster masses derived assuming hydrostatic equilibrium, the primary proxy for cluster cosmology studies. Sensitive measurements of Doppler broadening and shift of astrophysical lines, and the relative fluctuations in thermodynamical quantities (e.g., density, pressure, and entropy) are primary diagnostic tools. Forthcoming planned and proposed X-ray (with large etendue, throughput, and high spectral resolution) and SZ observatories will provide crucial information on the assembly and virialisation processes of clusters, involving turbulent eddies cascading at various spatial scales and larger gas bulk motions in their external regions to the depth or their potential wells.

In recent years, the outskirts of galaxy clusters have emerged as one of the new frontiers and unique laboratories for studying the growth of large scale structure in the universe. Modern cosmological hydrodynamical simulations make firm and testable predictions of the thermodynamic and chemical evolution of the X-ray emitting intracluster medium. However, recent X-ray and Sunyaev-Zeldovich effect observations have revealed enigmatic disagreements with theoretical predictions, which have motivated deeper investigations of a plethora of astrophysical processes operating in the virialization region in the cluster outskirts. Much of the physics of cluster outskirts is fundamentally different from that of cluster cores, which has been the main focus of X-ray cluster science over the past several decades. A next-generation X-ray telescope, equipped with sub-arcsecond spatial resolution over a large field of view along with a low and stable instrumental background, is required in order to reveal the full story of the growth of galaxy clusters and the cosmic web and their applications for cosmology.

We show that there is a new class of gas tails - slingshot tails - which form as a subhalo (i.e. a subcluster or early-type cluster galaxy) moves away from the cluster center towards the apocenter of its orbit. These tails can point perpendicular or even opposite to the subhalo direction of motion, not tracing the recent orbital path. Thus, the observed tail direction can be misleading, and we caution against naive conclusions regarding the subhalo's direction of motion based on the tail direction. A head-tail morphology of a galaxy's or subcluster's gaseous atmosphere is usually attributed to ram pressure stripping and the widely applied conclusion is that gas stripped tail traces the most recent orbit. However, during the slingshot tail stage, the subhalo is not being ram pressure stripped (RPS) and the tail is shaped by tidal forces more than just the ram pressure. Thus, applying a classic RPS scenario to a slingshot tail leads not only to an incorrect conclusion regarding the direction of motion, but also to incorrect conclusions in regard to the subhalo velocity, expected locations of shear flows, instabilities and mixing. We describe the genesis and morphology of slingshot tails using data from binary cluster merger simulations, discuss their observable features and how to distinguish them from classic RPS tails. We identify three examples from the literature that are not RPS tails but slingshot tails and discuss other potential candidates.

The total mass of a galaxy cluster is one of its most fundamental properties. Together with the redshift, the mass links observation and theory, allowing us to use the cluster population to test models of structure formation and to constrain cosmological parameters. Building on the rich heritage from X-ray surveys, new results from Sunyaev-Zeldovich and optical surveys have stimulated a resurgence of interest in cluster cosmology. These studies have generally found fewer clusters than predicted by the baseline Planck LCDM model, prompting a renewed effort on the part of the community to obtain a definitive measure of the true cluster mass scale. Here we review recent progress on this front. Our theoretical understanding continues to advance, with numerical simulations being the cornerstone of this effort. On the observational side, new, sophisticated techniques are being deployed in individual mass measurements and to account for selection biases in cluster surveys. We summarise the state of the art in cluster mass estimation methods and the systematic uncertainties and biases inherent in each approach, which are now well identified and understood, and explore how current uncertainties propagate into the cosmological parameter analysis. We discuss the prospects for improvements to the measurement of the mass scale using upcoming multi-wavelength data, and the future use of the cluster population as a cosmological probe.

The detailed velocity structure of the diffuse X-ray emitting intra-cluster medium (ICM) remains one of the last missing key ingredients in understanding the microphysical properties of these hot baryons and constraining our models of the growth and evolution of structure on the largest scales in the Universe. Direct measurements of the gas velocities from the widths and shifts of X-ray emission lines were recently provided for the central region of the Perseus Cluster of galaxies by $Hitomi$, and upcoming high-resolution X-ray microcalorimeters onboard $XRISM$ and $Athena$ are expected to extend these studies to many more systems. In the mean time, several other direct and indirect methods have been proposed for estimating the velocity structure in the ICM, ranging from resonant scattering to X-ray surface brightness fluctuation analysis, the kinematic Sunyaev-Zeldovich effect, or using optical line emitting nebulae in the brightest cluster galaxies as tracers of the motions of the ambient plasma. Here, we review and compare the existing estimates of the velocities of the hot baryons, as well as the various overlapping physical processes that drive motions in the ICM, and discuss the implications of these measurements for constraining the viscosity and identifying the source of turbulence in clusters of galaxies.

As the largest virialized structures in the universe, galaxy clusters continue to grow and accrete matter from the cosmic web. Due to the low gas density in the outskirts of clusters, measurements are very challenging, requiring extremely sensitive telescopes across the entire electromagnetic spectrum. Observations using X-rays, the Sunyaev-Zeldovich effect, and weak lensing and galaxy distributions from the optical band, have over the last decade helped to unravel this exciting new frontier of cluster astrophysics, where the infall and virialization of matter takes place. Here, we review the current state of the art in our observational and theoretical understanding of cluster outskirts, and discuss future prospects for exploration using newly planned and proposed observatories.

A number of merging galaxy clusters show the presence of large-scale radio emission associated with the intra-cluster medium (ICM). These synchrotron sources are generally classified as radio haloes and radio relics. Whilst it is commonly accepted that mergers play a crucial role in the formation of radio haloes and relics, not all the merging clusters show the presence of giant diffuse radio sources and this provides important information concerning current models. The Abell 781 complex is a spectacular system composed of an apparent chain of clusters on the sky. Its main component is undergoing a merger and hosts peripheral emission that is classified as a candidate radio relic and a disputed radio halo. We used new LOw Frequency ARay (LOFAR) observations at 143 MHz and archival Giant Metrewave Radio Telescope (GMRT) observations at 325 and 610 MHz to study radio emission from non-thermal components in the ICM of Abell 781. Complementary information came from XMM-Newton data, which allowed us to investigate the connection with the thermal emission and its complex morphology. The origin of the peripheral emission is still uncertain. We speculate that it is related to the interaction between a head tail radio galaxy and shock. However, the current data allow us only to set an upper limit of $\mathcal{M} < 1.4$ on the Mach number of this putative shock. Instead, we successfully characterise the surface brightness and temperature jumps of a shock and two cold fronts in the main cluster component of Abell 781. Their positions suggest that the merger is involving three substructures. We do not find any evidence for a radio halo either at the centre of this system or in the other clusters of the chain. We place an upper limit to the diffuse radio emission in the main cluster of Abell 781 that is a factor of 2 below the current radio power-mass relation for giant radio haloes.

We present the reconstruction of hydrostatic mass profiles in 13 X-ray luminous galaxy clusters that have been mapped in their X-ray and SZ signal out to $R_{200}$ for the XMM-Newton Cluster Outskirts Project (X-COP). Using profiles of the gas temperature, density and pressure that have been spatially resolved out to (median value) 0.9 $R_{500}$, 1.8 $R_{500}$, and 2.3 $R_{500}$, respectively, we are able to recover the hydrostatic gravitating mass profile with several methods and using different mass models. The hydrostatic masses are recovered with a relative (statistical) median error of 3% at $R_{500}$ and 6% at $R_{200}$. By using several different methods to solve the equation of the hydrostatic equilibrium, we evaluate some of the systematic uncertainties to be of the order of 5% at both $R_{500}$ and $R_{200}$. A Navarro-Frenk-White profile provides the best-fit in nine cases out of 13, with the remaining four cases that do not show a statistically significant tension with it. The distribution of the mass concentration follows the correlations with the total mass predicted from numerical simulations with a scatter of 0.18 dex, with an intrinsic scatter on the hydrostatic masses of 0.15 dex. We compare them with the estimates of the total gravitational mass obtained through X-ray scaling relations applied to $Y_X$, gas fraction and $Y_{SZ}$, and from weak lensing and galaxy dynamics techniques, and measure a substantial agreement with the results from scaling laws, from WL at both $R_{500}$ and $R_{200}$ (with differences below 15%), from cluster velocity dispersions, but a significant tension with the caustic masses that tend to underestimate the hydrostatic masses by 40% at $R_{200}$. We also compare these measurements with predictions from alternative models to the Cold Dark Matter, like the "Emergent Gravity" and MOND scenarios.

The hot plasma in galaxy clusters is expected to be heated to high temperatures through shocks and adiabatic compression. The thermodynamical properties of the gas encode information on the processes leading to the thermalization of the gas in the cluster's potential well as well as non-gravitational processes such as gas cooling, AGN feedback and kinetic energy. In this work we present the radial profiles of the thermodynamic properties of the intracluster medium (ICM) out to the virial radius for a sample of 12 galaxy clusters selected from the Planck all-sky survey. We determine the universal profiles of gas density, temperature, pressure, and entropy over more than two decades in radius. We exploit jointly X-ray information from XMM and Sunyaev-Zel'dovich constraints from Planck to recover thermodynamic properties out to 2 R500. We provide average functional forms for the radial dependence of the main quantities and quantify the slope and intrinsic scatter of the population as a function of radius. We find that gas density and pressure profiles steepen steadily with radius, in excellent agreement with previous observational results. Entropy profiles beyond R500 closely follow the predictions for the gravitational collapse of structures. The scatter in all thermodynamical quantities reaches a minimum in the range [0.2-0.8] R500 and increases outwards. Somewhat surprisingly, we find that pressure is substantially more scattered than temperature and density. Our results indicate that once accreting substructures are properly excised, the properties of the ICM beyond the cooling region R > 0.3 R500) follow remarkably well the predictions of simple gravitational collapse and require little non-gravitational corrections.

In the currently debated context of using clusters of galaxies as cosmological probes, the need for well-defined cluster samples is critical. The XXL Survey has been specifically designed to provide a well characterised sample of some 500 X-ray detected clusters suitable for cosmological studies. The main goal of present article is to make public and describe the properties of the cluster catalogue in its present state, as well as of associated catalogues as super-clusters and fossil groups. We release a sample containing 365 clusters in total. We give the details of the follow-up observations and explain the procedure adopted to validate the cluster spectroscopic redshifts. Considering the whole XXL cluster sample, we have provided two types of selection, both complete in a particular sense: one based on flux-morphology criteria, and an alternative based on the [0.5-2] keV flux within one arcmin of the cluster centre. We have also provided X-ray temperature measurements for 80$\%$ of the clusters having a flux larger than 9$\times$10$^{-15}$$\rm \thinspace erg \, s^{-1} \, cm^{-2}$. Our cluster sample extends from z$\sim$0 to z$\sim$1.2, with one cluster at z$\sim$2. Clusters were identified through a mean number of six spectroscopically confirmed cluster members. Our updated luminosity function and luminosity-temperature relation are compatible with our previous determinations based on the 100 brightest clusters, but show smaller uncertainties. We also present an enlarged list of super-clusters and a sample of 18 possible fossil groups. This intermediate publication is the last before the final release of the complete XXL cluster catalogue when the ongoing C2 cluster spectroscopic follow-up is complete. It provides a unique inventory of medium-mass clusters over a 50~\dd\ area out to z$\sim$1.

Galaxy clusters are the endpoints of structure formation and are continuously growing through the merging and accretion of smaller structures. Numerical simulations predict that a fraction of their energy content is not yet thermalized, mainly in the form of kinetic motions (turbulence, bulk motions). Measuring the level of non-thermal pressure support is necessary to understand the processes leading to the virialization of the gas within the potential well of the main halo and to calibrate the biases in hydrostatic mass estimates. We present high-quality measurements of hydrostatic masses and intracluster gas fraction out to the virial radius for a sample of 12 nearby clusters with available XMM-Newton and Planck data. We compare our hydrostatic gas fractions with the expected universal gas fraction to constrain the level of non-thermal pressure support. We find that hydrostatic masses require little correction and infer a median non-thermal pressure fraction of $\sim6\%$ and $\sim10\%$ at $R_{500}$ and $R_{200}$, respectively. Our values are lower than the expectations of hydrodynamical simulations, possibly implying a faster thermalization of the gas. If instead we use the mass calibration adopted by the Planck team, we find that the gas fraction of massive local systems implies a mass bias $1-b=0.85\pm0.05$ for SZ-derived masses, with some evidence for a mass-dependent bias. Conversely, the high bias required to match Planck CMB and cluster count cosmology is excluded by the data at high significance, unless the most massive halos are missing a substantial fraction of their baryons.

A recent comparison of the massive galaxy cluster Abell 2744 with the Millennium XXL (MXXL) N-body simulation has hinted at a tension between the observed substructure distribution and the predictions of LambdaCDM. Follow-up investigations indicated that this could be due to the contribution from the host halo and the subhalo finding algorithm used. To be independent of any subhalo finding algorithm, we therefore investigate the particle data of the MXXL simulation directly. We propose a new method to find substructures in 2D mass maps using a wavelet transform, which treats the simulation and observations equally. Using the same criteria to define a subhalo in observations and simulated data, we find three Abell 2744 analogues in the MXXL simulation. Thus the observations in Abell 2744 are in agreement with the predictions of LambdaCDM. We investigate the reasons for the discrepancy between the results obtained from the SUBFIND and full particle data analyses. We find that this is due to incompatible substructure definitions in observations and SUBFIND.

Due to their late formation in cosmic history, clusters of galaxies are not fully in hydrostatic equilibrium and the gravitational pull of their mass at a given radius is expected not to be entirely balanced by the thermal gas pressure. Turbulence may supply additional pressure, and recent (X-ray and SZ) hydrostatic mass reconstructions claim a pressure support of $\sim 5-15\%$ of the total pressure at $R_{\rm 200}$. In this work we show that, after carefully disentangling bulk from small-scale turbulent motions in high-resolution simulations of galaxy clusters, we can constrain which fraction of the gas kinetic energy effectively provides pressure support in the cluster's gravitational potential. While the ubiquitous presence of radial inflows in the cluster can lead to significant bias in the estimate of the non-thermal pressure support, we report that only a part of this energy effectively acts as a source of pressure, providing a support of the order of $\sim 10\%$ of the total pressure at $R_{\rm 200}$.

We present a gravitational lensing and X-ray analysis of a massive galaxy cluster and its surroundings. The core of MACS\,J0717.5+3745 ($M(R<1\,{\rm Mpc})\sim$\,$2$$\times$$10^{15}\,\msun$, $z$=$0.54$) is already known to contain four merging components. We show that this is surrounded by at least seven additional substructures with masses ranging from $3.8-6.5\times10^{13}\,\msun$, at projected radii $1.6$ to $4.9$\,Mpc. We compare MACS\,J0717 to mock lensing and X-ray observations of similarly rich clusters in cosmological simulations. The low gas fraction of substructures predicted by simulations turns out to match our observed values of $1$--$4\%$. Comparing our data to three similar simulated halos, we infer a typical growth rate and substructure infall velocity. That suggests MACS\,J0717 could evolve into a system similar to, but more massive than, Abell\,2744 by $z=0.31$, and into a $\sim$\,$10^{16}\,\msun$ supercluster by $z=0$. The radial distribution of infalling substructure suggests that merger events are strongly episodic; however we find that the smooth accretion of surrounding material remains the main source of mass growth even for such massive clusters.

We present a new cosmological probe for galaxy clusters, the halo sparsity. This characterises halos in terms of the ratio of halo masses measured at two different radii and carries cosmological information encoded in the halo mass profile. Building upon the work of Balmes et al. (2014) we test the properties of the sparsity using halo catalogs from a numerical N-body simulation of ($2.6$ Gpc/h)$^3$ volume with $4096^3$ particles. We show that at a given redshift the average sparsity can be predicted from prior knowledge of the halo mass function. This provides a quantitative framework to infer cosmological parameter constraints using measurements of the sparsity of galaxy clusters. We show this point by performing a likelihood analysis of synthetic datasets with no systematics, from which we recover the input fiducial cosmology. We also perform a preliminary analysis of potential systematic errors and provide an estimate of the impact of baryonic effects on sparsity measurements. We evaluate the sparsity for a sample of 104 clusters with hydrostatic masses from X-ray observations and derive constraints on the cosmic matter density $\Omega_m$ and the normalisation amplitude of density fluctuations at the $8$ Mpc h$^{-1}$ scale, $\sigma_8$. Assuming no systematics, we find $\Omega_m=0.42\pm 0.17$ and $\sigma_8=0.80\pm 0.31$ at $1\sigma$, corresponding to $S_8\equiv \sigma_8\sqrt{\Omega_m}=0.48\pm 0.11$. Future cluster surveys may provide opportunities for precise measurements of the sparsity. A sample of a few hundreds clusters with mass estimate errors at a few percent level can provide competitive cosmological parameter constraints complementary to those inferred from other cosmic probes.

The mass of galaxy clusters is not a direct observable, nonetheless it is commonly used to probe cosmological models. Based on the combination of all main cluster observables, that is, the X-ray emission, the thermal Sunyaev-Zel'dovich (SZ) signal, the velocity dispersion of the cluster galaxies, and gravitational lensing, the gravitational potential of galaxy clusters can be jointly reconstructed. We derive the two main ingredients required for this joint reconstruction: the potentials individually reconstructed from the observables and their covariance matrices, which act as a weight in the joint reconstruction. We show here the method to derive these quantities. The result of the joint reconstruction applied to a real cluster will be discussed in a forthcoming paper. We apply the Richardson-Lucy deprojection algorithm to data on a two-dimensional (2D) grid. We first test the 2D deprojection algorithm on a $\beta$-profile. Assuming hydrostatic equilibrium, we further reconstruct the gravitational potential of a simulated galaxy cluster based on synthetic SZ and X-ray data. We then reconstruct the projected gravitational potential of the massive and dynamically active cluster Abell 2142, based on the X-ray observations collected with XMM-Newton and the SZ observations from the Planck satellite. Finally, we compute the covariance matrix of the projected reconstructed potential of the cluster Abell 2142 based on the X-ray measurements collected with XMM-Newton. The gravitational potentials of the simulated cluster recovered from synthetic X-ray and SZ data are consistent, even though the potential reconstructed from X-rays shows larger deviations from the true potential. Regarding Abell 2142, the projected gravitational cluster potentials recovered from SZ and X-ray data reproduce well the projected potential inferred from gravitational-lensing observations. (abridged)

We propose a novel method to constrain turbulence and bulk motions in massive galaxies, groups and clusters, exploring both simulations and observations. As emerged in the recent picture of the top-down multiphase condensation, the hot gaseous halos are tightly linked to all other phases in terms of cospatiality and thermodynamics. While hot halos (10^7 K) are perturbed by subsonic turbulence, warm (10^4 K) ionized and neutral filaments condense out of the turbulent eddies. The peaks condense into cold molecular clouds (< 100 K) raining in the core via chaotic cold accretion (CCA). We show all phases are tightly linked via the ensemble (wide-aperture) velocity dispersion along the line of sight. The correlation arises in complementary long-term AGN feedback simulations and high-resolution CCA runs, and is corroborated by the combined Hitomi and new IFU measurements in Perseus cluster. The ensemble multiphase gas distributions are characterized by substantial spectral line broadening (100-200 km/s) with mild line shift. On the other hand, pencil-beam detections sample the small-scale clouds displaying smaller broadening and significant line shift up to several 100 km/s, with increased scatter due to the turbulence intermittency. We present new ensemble sigma_v of the warm Halpha+[NII] gas in 72 observed cluster/group cores: the constraints are consistent with the simulations and can be used as robust proxies for the turbulent velocities, in particular for the challenging hot plasma (otherwise requiring extremely long X-ray exposures). We show the physically motivated criterion C = t_cool/t_eddy ~ 1 best traces the condensation extent region and presence of multiphase gas in observed clusters/groups. The ensemble method can be applied to many available datasets and can substantially advance our understanding of multiphase halos in light of the next-generation multiwavelength missions.

We present the joint analysis of the X-ray and SZ signals in A2319, the galaxy cluster with the highest signal-to-noise ratio in Planck maps and that has been surveyed within our XMM Cluster Outskirts Project (X-COP). We recover the thermodynamical profiles by the geometrical deprojection of the X-ray surface brightness, of the SZ comptonization parameter, and an accurate and robust spectroscopic measurements of the temperature. We resolve the clumpiness of the density to be below 20 per cent demonstrating that most of this clumpiness originates from the ongoing merger and can be associated to large-scale inhomogeneities. This analysis is done in azimuthally averaged radial bins and in eight independent angular sectors, enabling us to study in details the azimuthal variance of the recovered properties. Given the exquisite quality of the X-ray and SZ datasets, we constrain at $R_{200}$ the total hydrostatic mass, modelled with a NFW profile, with very high precision ($M_{200} = 9.76 \pm 0.16^{stat.} \pm 0.31^{syst.} \times 10^{14} M_\odot$). We identify the ongoing merger and how it is affecting differently the gas properties in the resolved azimuthal sectors. We have several indications that the merger has injected a high level of non-thermal pressure in this system: the clumping free density profile is above the average profile obtained by stacking Rosat observations; the gas mass fraction exceeds the expected cosmic gas fraction beyond $R_{500}$; the pressure profile is flatter than the fit obtained by the Planck collaboration; the entropy profile is flatter than the mean one predicted from non-radiative simulations; the analysis in azimuthal sectors has revealed that these deviations occur in a preferred region of the cluster. All these tensions are resolved by requiring a relative support of about 40 per cent from non-thermal to the total pressure at $R_{200}$.

We estimate the "non-gravitational" entropy-injection profiles, $\Delta K$, and the resultant energy feedback profiles, $\Delta E$, of the intracluster medium for 17 clusters using their Planck SZ and ROSAT X-Ray observations, spanning a large radial range from $0.2r_{500}$ up to $r_{200}$. The feedback profiles are estimated by comparing the observed entropy, at fixed gas mass shells, with theoretical entropy profiles predicted from non-radiative hydrodynamic simulations. We include non-thermal pressure and gas clumping in our analysis. The inclusion of non-thermal pressure and clumping results in changing the estimates for $r_{500}$ and $r_{200}$ by 10\%-20\%. When clumpiness is not considered it leads to an under-estimation of $\Delta K\approx300$ keV cm$^2$ at $r_{500}$ and $\Delta K\approx1100$ keV cm$^2$ at $r_{200}$. On the other hand, neglecting non-thermal pressure results in an over-estimation of $\Delta K\approx 100$ keV cm$^2$ at $r_{500}$ and under-estimation of $\Delta K\approx450$ keV cm$^2$ at $r_{200}$. For the estimated feedback energy, we find that ignoring clumping leads to an under-estimation of energy per particle $\Delta E\approx1$ keV at $r_{500}$ and $\Delta E\approx1.5$ keV at $r_{200}$. Similarly, neglect of the non-thermal pressure results in an over-estimation of $\Delta E\approx0.5$ keV at $r_{500}$ and under-estimation of $\Delta E\approx0.25$ keV at $r_{200}$. We find entropy floor of $\Delta K\approx300$ keV cm$^2$ is ruled out at $\approx3\sigma$ throughout the entire radial range and $\Delta E\approx1$ keV at more than 3$\sigma$ beyond $r_{500}$, strongly constraining ICM pre-heating scenarios. We also demonstrate robustness of results w.r.t sample selection, X-Ray analysis procedures, entropy modeling etc.

The XXL survey currently covers two 25 sq. deg. patches with XMM observations of ~10ks. We summarise the scientific results associated with the first release of the XXL data set, that occurred mid 2016. We review several arguments for increasing the survey depth to 40 ks during the next decade of XMM operations. X-ray (z<2) cluster, (z<4) AGN and cosmic background survey science will then benefit from an extraordinary data reservoir. This, combined with deep multi-$\lambda$ observations, will lead to solid standalone cosmological constraints and provide a wealth of information on the formation and evolution of AGN, clusters and the X-ray background. In particular, it will offer a unique opportunity to pinpoint the z>1 cluster density. It will eventually constitute a reference study and an ideal calibration field for the upcoming eROSITA and Euclid missions.

In the local Universe, the growth of massive galaxy clusters mainly operates through the continuous accretion of group-scale systems. The infalling group in Abell 2142 is the poster child of such an accreting group, and as such, it is an ideal target to study the astrophysical processes induced by structure formation. We present the results of a deep (200 ks) observation of this structure with Chandra, which highlights the complexity of this system in exquisite detail. In the core of the group, the spatial resolution of Chandra reveals the presence of a leading edge and a complex AGN-induced activity. The morphology of the stripped gas tail appears straight in the innermost 250 kpc, suggesting that magnetic draping efficiently shields the gas from its surroundings. However, beyond $\sim300$ kpc from the core, the tail flares and the morphology becomes strongly irregular, which could be explained by a breaking of the drape, e.g. because of turbulent motions. The power spectrum of surface-brightness fluctuations is relatively flat ($P_{2D}\propto k^{-2.3}$), which indicates that thermal conduction is strongly inhibited even beyond the region where magnetic draping is effective. The amplitude of density fluctuations in the tail is consistent with a mild level of turbulence with a Mach number $M_{3D}\sim0.1-0.25$. Overall, our results show that the processes leading to the thermalization and mixing of the infalling gas are slow and relatively inefficient.

We present the radial distribution of the dark matter in two massive, X-ray luminous galaxy clusters, Abell~2142 and Abell~2319, and compare it with the quantity predicted as apparent manifestation of the baryonic mass in the context of the "Emergent Gravity" scenario, recently suggested from Verlinde (2016). Thanks to the observational strategy of the \xmm\ Cluster Outskirt Programme (X-COP), using the X-ray emission mapped with \xmm\ and the SZ signal in the Planck survey, we recover the gas density, temperature and thermal pressure profiles up to $\sim R_{200}$, allowing to constrain at unprecedented level the total mass through the hydrostatic equilibrium equation. We show that, also including systematic uncertainties related to the X-ray based mass modelling, the apparent "dark" matter shows a radial profile that has a shape different from the traditional dark matter distribution, with larger discrepancies (by a factor 2--3) in the inner ($r<200$ kpc) cluster's regions and a remarkable agreement only across $R_{500}$.

Giant radio halos are Mpc-scale diffuse radio sources associated with the central regions of galaxy clusters. The most promising scenario to explain the origin of these sources is that of turbulent re-acceleration, in which MeV electrons injected throughout the formation history of galaxy clusters are accelerated to higher energies by turbulent motions mostly induced by cluster mergers. In this Letter, we use the amplitude of density fluctuations in the intracluster medium as a proxy for the turbulent velocity and apply this technique to a sample of 51 clusters with available radio data. Our results indicate a segregation in the turbulent velocity of radio halo and radio quiet clusters, with the turbulent velocity of the former being on average higher by about a factor of two. The velocity dispersion recovered with this technique correlates with the measured radio power through the relation $P_{\rm radio}\propto\sigma_v^{3.3\pm0.7}$, which implies that the radio power is nearly proportional to the turbulent energy rate. Our results provide an observational confirmation of a key prediction of the turbulent re-acceleration model and possibly shed light on the origin of radio halos.

We study the intra-cluster magnetic field in the poor galaxy cluster Abell 194 by complementing radio data, at different frequencies, with data in the optical and X-ray bands. We analyze new total intensity and polarization observations of Abell 194 obtained with the Sardinia Radio Telescope (SRT). We use the SRT data in combination with archival Very Large Array observations to derive both the spectral aging and Rotation Measure (RM) images of the radio galaxies 3C40A and 3C40B embedded in Abell 194. The optical analysis indicates that Abell 194 does not show a major and recent cluster merger, but rather agrees with a scenario of accretion of small groups. Under the minimum energy assumption, the lifetimes of synchrotron electrons in 3C40B measured from the spectral break are found to be 157 Myrs. The break frequency image and the electron density profile inferred from the X-ray emission are used in combination with the RM data to constrain the intra-cluster magnetic field power spectrum. By assuming a Kolmogorov power law power spectrum, we find that the RM data in Abell 194 are well described by a magnetic field with a maximum scale of fluctuations of Lambda_max=64 kpc and a central magnetic field strength of <B0>=1.5 microG. Further out, the field decreases with the radius following the gas density to the power of eta=1.1. Comparing Abell 194 with a small sample of galaxy clusters, there is a hint of a trend between central electron densities and magnetic field strengths.

We report on a spectral study at radio frequencies of the giant radio halo in A2142 (z=0.0909), which we performed to explore its nature and origin. A2142 is not a major merger and the presence of a giant radio halo is somewhat surprising. We performed deep radio observations with the GMRT at 608 MHz, 322 MHz, and 234 MHz and with the VLA in the 1-2 GHz band. We obtained high-quality images at all frequencies in a wide range of resolutions. The radio halo is well detected at all frequencies and extends out to the most distant cold front in A2142. We studied the spectral index in two regions: the central part of the halo and a second region in the direction of the most distant south-eastern cold front, selected to follow the bright part of the halo and X-ray emission. We complemented our observations with a preliminary LOFAR image at 118 MHz and with the re-analysis of archival VLA data at 1.4 GHz. The two components of the radio halo show different observational properties. The central brightest part has higher surface brightess and a spectrum whose steepness is similar to those of the known radio halos, i.e. $\alpha^{\rm 1.78~GHz}_{\rm 118~MHz}=1.33\pm 0.08$. The ridge, which fades into the larger scale emission, is broader in size and has considerably lower surface brightess and a moderately steeper spectrum, i.e. $\alpha^{\rm 1.78~GHz}_{\rm 118~MHz}\sim 1.5$. We propose that the brightest part of the radio halo is powered by the central sloshing in A2142, similar to what has been suggested for mini-halos, or by secondary electrons generated by hadronic collisions in the ICM. On the other hand, the steeper ridge may probe particle re-acceleration by turbulence generated either by stirring the gas and magnetic fields on a larger scale or by less energetic mechanisms, such as continuous infall of galaxy groups or an off-axis merger.

We characterized the population of galaxy clusters detected with the SZ effect with Planck, by measuring the cool core state of the objects in a well-defined subsample of the Planck catalogue. We used as indicator the concentration parameter Santos et al. (2008). The fraction of cool core clusters is $29 \pm 4 \%$ and does not show significant indications of evolution in the redshift range covered by our sample. We compare the distribution of the concentration parameter in the Planck sample with the one of the X-ray selected sample MACS (Mann & Ebeling, 2011): the distributions are significantly different and the cool core fraction in MACS is much higher ($59 \pm 5 \%$). Since X-ray selected samples are known to be biased towards cool cores due to the presence of their prominent surface brightness peak, we simulated the impact of the "cool core bias" following Eckert et al. (2011). We found that it plays a large role in the difference between the fractions of cool cores in the two samples. We examined other selection effects that could in principle bias SZ-surveys against cool cores but we found that their impact is not sufficient to explain the difference between Planck and MACS. The population of X-ray under-luminous objects, which are found in SZ-surveys but missing in X-ray samples (Planck Collaboration 2016), could possibly contribute to the difference, as we found most of them to be non cool cores, but this hypothesis deserves further investigation.

The massive substructures found in Abell 2744 by Jauzac et al. (2016) present a challenge to the cold dark matter paradigm due to their number and proximity to the cluster centre. We use one of the biggest N-body simulations, the Millennium XXL, to investigate the substructure in a large sample of massive dark matter haloes. A range of effects which influence the comparison with the observations is considered, extending the preliminary evaluation carried out by Jauzac et al. (2016). There are many tens of haloes in the simulation with a total mass comparable with or larger than that of Abell 2744. However, we find no haloes with a number and distribution of massive substructures (> 5 10^13 Msun) that is close to that inferred from the observations of Abell 2744. The application of extreme value statistics suggests that we would need a simulation of at least ten times the volume of the Millennium XXL to find a single dark matter halo with a similar internal structure to Abell 2744. Explaining the distribution of massive substructures in clusters is a new hurdle for hierarchical models to negotiate, which is not weakened by appeals to baryonic physics or uncertainty over the nature of the dark matter particle.

Galaxy cluster merger shocks are the main agent for the thermalization of the intracluster medium and the energization of cosmic ray particles in it. Shock propagation changes the state of the tenuous intracluster plasma, and the corresponding signal variations are measurable with the current generation of X-ray and Sunyaev-Zel'dovich (SZ) effect instruments. Additionally, non-thermal electrons (re-)energized by the shocks sometimes give rise to extended and luminous synchrotron sources known as radio relics, which are prominent indicators of shocks propagating roughly in the plane of the sky. In this short review, we discuss how the joint modeling of the non-thermal and thermal signal variations across radio relic shock fronts is helping to advance our knowledge of the gas thermodynamical properties and magnetic field strengths in the cluster outskirts. We describe the first use of the SZ effect to measure the Mach numbers of relic shocks, for both the nearest (Coma) and the farthest (El Gordo) clusters with known radio relics.

We present ALMA measurements of a merger shock using the thermal Sunyaev-Zel'dovich (SZ) effect signal, at the location of a radio relic in the famous El Gordo galaxy cluster at $z \approx 0.9$. Multi-wavelength analysis in combination with the archival Chandra data and a high-resolution radio image provides a consistent picture of the thermal and non-thermal signal variation across the shock front and helps to put robust constraints on the shock Mach number as well as the relic magnetic field. We employ a Bayesian analysis technique for modeling the SZ and X-ray data self-consistently, illustrating respective parameter degeneracies. Combined results indicate a shock with Mach number ${\cal M} = 2.4^{+1.3}_{-0.6}$, which in turn suggests a high value of the magnetic field (of the order of $4-10 ~\mu$G) to account for the observed relic width at 2 GHz. At roughly half the current age of the universe, this is the highest-redshift direct detection of a cluster shock to date, and one of the first instances of an ALMA-SZ observation in a galaxy cluster. It shows the tremendous potential for future ALMA-SZ observations to detect merger shocks and other cluster substructures out to the highest redshifts.

Galaxy clusters are thought to grow hierarchically through the continuous merging and accretion of smaller structures across cosmic time. In the local Universe, these phenomena are still active in the outer regions of massive clusters ($R>R_{500}$), where the matter distribution is expected to become clumpy and asymmetric because of the presence of accreting structures. We present the XMM-Newton Cluster Outskirts Project (X-COP), which targets the outer regions of a sample of 13 massive clusters ($M_{500}>3\times10^{14}M_\odot$) in the redshift range 0.04-0.1 at uniform depth. The sample was selected based on the signal-to-noise ratio in the Planck Sunyaev-Zeldovich (SZ) survey with the aim of combining high-quality X-ray and SZ constraints throughout the entire cluster volume. Our observing strategy allows us to reach a sensitivity of $3\times10^{-16}$ ergs cm$^{-2}$ s$^{-1}$ arcmin$^{-2}$ in the [0.5-2.0] keV range thanks to a good control of systematic uncertainties. The combination of depth and field of view achieved in X-COP will allow us to pursue the following main goals: i) measure the distribution of entropy and thermal energy to an unprecedented level of precision; ii) assess the presence of non-thermal pressure support in cluster outskirts; iii) study the occurrence and mass distribution of infalling gas clumps. We illustrate the capabilities of the program with a pilot study on the cluster Abell 2142.

Non-gravitational feedback affects the nature of the intra-cluster medium (ICM). X-ray cooling of the ICM and in situ energy feedback from AGN's and SNe as well as {\it preheating} of the gas at epochs preceding the formation of clusters are proposed mechanisms for such feedback. While cooling and AGN feedbacks are dominant in cluster cores, the signatures of a preheated ICM are expected to be present even at large radii. To estimate the degree of preheating, with minimum confusion from AGN feedback/cooling, we study the excess entropy and non-gravitational energy profiles upto $r_{200}$ for a sample of 17 galaxy clusters using joint data sets of {\it Planck} SZ pressure and {\it ROSAT/PSPC} gas density profiles. The canonical value of preheating entropy floor of $\gtrsim 300$ keV cm$^2$, needed in order to match cluster scalings, is ruled out at $\approx 3\sigma$. We also show that the feedback energy of 1 keV/particle is ruled out at 5.2$\sigma$ beyond $r_{500}$. Our analysis takes both non-thermal pressure and clumping into account which can be important in outer regions. Our results based on the direct probe of the ICM in the outermost regions do not support any significant preheating.

The growth of brightest cluster galaxies is closely related to the properties of their host cluster. We present evidence for dry mergers as the dominant source of BCG mass growth at $z\lesssim1$ in the XXL 100 brightest cluster sample. We use the global red sequence, H$\alpha$ emission and mean star formation history to show that BCGs in the sample possess star formation levels comparable to field ellipticals of similar stellar mass and redshift. XXL 100 brightest clusters are less massive on average than those in other X-ray selected samples such as LoCuSS or HIFLUGCS. Few clusters in the sample display high central gas concentration, rendering inefficient the growth of BCGs via star formation resulting from the accretion of cool gas. Using measures of the relaxation state of their host clusters, we show that BCGs grow as relaxation proceeds. We find that the BCG stellar mass corresponds to a relatively constant fraction 1\%\ of the total cluster mass in relaxed systems. We also show that, following a cluster scale merger event, the BCG stellar mass lags behind the expected value from the M$_{cluster}$ - M$_{BCG}$ relation but subsequently accretes stellar mass via dry mergers as the BCG and cluster evolve towards a relaxed state.

We present the results of the search for decaying dark matter with particle mass in the 6-40 keV range with NuSTAR deep observations of COSMOS and ECDFS empty sky fields. We show that main contribution to the decaying dark matter signal from the Milky Way galaxy comes through the aperture of the NuSTAR detector, rather than through the focusing optics. High sensitivity of the NuSTAR detector, combined with the large aperture and large exposure times of the two observation fields allow us to improve previously existing constraints on the dark matter decay time by up to an order of magnitude in the mass range 10-30 keV. In the particular case of the nuMSM sterile neutrino dark matter, our constraints impose an upper bound m<20 keV on the dark matter particle mass. We report detection of four unidentified spectral lines in our data set. These line detections are either due to the systematic effects (uncertainties of calibrations of the NuSTAR detectors) or have an astrophysical origin. We discuss different possibilities for testing the nature of the detected lines.

Radio relics are Mpc-scale diffuse radio sources at the peripheries of galaxy clusters which are thought to trace outgoing merger shocks. We present XMM-Newton and Suzaku observations of the galaxy cluster Abell 2744 (z=0.306), which reveal the presence of a shock front 1.5 Mpc East of the cluster core. The surface-brightness jump coincides with the position of a known radio relic. Although the surface-brightness jump indicates a weak shock with a Mach number $\mathcal{M}=1.7_{-0.3}^{+0.5}$, the plasma in the post-shock region has been heated to a very high temperature ($\sim13$ keV) by the passage of the shock wave. The low acceleration efficiency expected from such a weak shock suggests that mildly relativistic electrons have been re-accelerated by the passage of the shock front.

Context. Galaxy clusters are continuously growing through the accretion of matter in their outskirts. This process induces inhomogeneities in the gas density distribution (clumping) which need to be taken into account to recover the physical properties of the intracluster medium (ICM) at large radii. Aims. We studied the thermodynamic properties in the outskirts (R > R500) of the massive galaxy cluster Abell 2142 by combining the Sunyaev Zel'dovich (SZ) effect with the X-ray signal. Methods. We combined the SZ pressure profile measured by Planck with the XMM-Newton gas density profile to recover radial profiles of temperature, entropy and hydrostatic mass out to 2R500. We used a method that is insensitive to clumping to recover the gas density, and we compared the results with traditional X-ray measurement techniques. Results. When taking clumping into account, our joint SZ/X-ray entropy profile is consistent with the predictions from pure gravitational collapse, whereas a significant entropy flattening is found when the effect of clumping is neglected. The hydrostatic mass profile recovered using joint X-ray/SZ data agrees with that obtained from spectroscopic X-ray measurements and with mass reconstructions obtained through weak lensing and galaxy kinematics. Conclusions. We found that clumping can explain the entropy flattening observed by Suzaku in the outskirts of several clusters. When using a method insensitive to clumping for the reconstruction of the gas density, the thermodynamic properties of Abell 2142 are compatible with the assumption that the thermal gas pressure sustains gravity and that the entropy is injected at accretion shocks, with no need to evoke more exotic physics. Our results highlight the need for X-ray observations with sufficient spatial resolution, and large collecting area, to understand the processes at work in cluster outer regions.

We present a joint optical/X-ray analysis of the massive galaxy cluster Abell 2744 (z=0.308). Our strong- and weak-lensing analysis within the central region of the cluster, i.e., at R<1Mpc from the brightest cluster galaxy, reveals eight substructures, including the main core. All of these dark-matter halos are detected with a significance of at least 5sigma and feature masses ranging from 0.5 to 1.4x10^{14}Msun within R<150kpc. Merten et al. (2011) and Medezinski et al. (2016) substructures are also detected by us. We measure a slightly higher mass for the main core component than reported previously and attribute the discrepancy to the inclusion of our tightly constrained strong-lensing mass model built on Hubble Frontier Fields data. X-ray data obtained by XMM-Newton reveal four remnant cores, one of them a new detection, and three shocks. Unlike Merten et al. (2011), we find all cores to have both dark and luminous counterparts. A comparison with clusters of similar mass in the MXXL simulations yields no objects with as many massive substructures as observed in Abell 2744, confirming that Abell 2744 is an extreme system. We stress that these properties still do not constitute a challenge to $\Lambda$CDM, as caveats apply to both the simulation and the observations: for instance, the projected mass measurements from gravitational lensing and the limited resolution of the sub-haloes finders. We discuss implications of Abell 2744 for the plausibility of different dark-matter candidates and, finally, measure a new upper limit on the self-interaction cross-section of dark matter of sigma_{DM}<1.28cm2/g(68\% CL), in good agreement with previous results from Harvey et al. (2015).

We present results of a search for the 3.5 keV emission line in our recent very long (~ 1.4 Ms) XMM-Newton observation of the Draco dwarf spheroidal galaxy. The astrophysical X-ray emission from such dark matter-dominated galaxies is faint, thus they provide a test for the dark matter origin of the 3.5 keV line previously detected in other massive, but X-ray bright objects, such as galaxies and galaxy clusters. We do not detect a statistically significant emission line from Draco; this constrains the lifetime of a decaying dark matter particle to tau > (7-9) x 10^27 s at 95% CL (combining all three XMM-Newton cameras; the interval corresponds to the uncertainty of the dark matter column density in the direction of Draco). The PN camera, which has the highest sensitivity of the three, does show a positive spectral residual (above the carefully modeled continuum) at E = 3.54 +/- 0.06 keV with a 2.3 sigma significance. The two MOS cameras show less-significant or no positive deviations, consistently within 1 sigma with PN. Our Draco limit on tau is consistent with previous detections in the stacked galaxy clusters, M31 and the Galactic Center within their 1-2 sigma uncertainties, but is inconsistent with the high signal from the core of the Perseus cluster (which has itself been inconsistent with the rest of the detections). We conclude that this Draco observation does not exclude the dark matter interpretation of the 3.5 keV line in those objects.

X-ray extragalactic surveys are ideal laboratories for the study of the evolution and clustering of active galactic nuclei (AGN). The XXL Survey spans two fields of a combined 50 $deg^2$ observed for more than 6Ms with XMM-Newton, occupying the parameter space between deep surveys and very wide area surveys; at the same time it benefits from a wealth of ancillary data. This paper marks the first release of the XXL point source catalogue selected in the 2-10 keV energy band with limiting flux $F_{2-10keV}=4.8\cdot10^{-14}\rm{erg\,s^{-1}\,cm^{-2}}$. We use both public and proprietary data sets to identify the counterparts of the X-ray point-like sources and improved upon the photometric redshift determination for AGN by applying a Random Forest classification trained to identify for each object the optimal photometric redshift model library. We also assign a probability to each source to be a star or an outlier. We model with Bayesian analysis the X-ray spectra assuming a power-law model with the presence of an absorbing medium. We find an average unabsorbed photon index of $\Gamma=1.85$ and average hydrogen column density $\log{N_{H}}=21.07 cm^{-2}$. We find no trend of $\Gamma$ or $N_H$ with redshift and a fraction of 26% absorbed sources ($\log N_{H}>22$). We show that the XXL-1000-AGN number counts extended the number counts of the COSMOS survey to higher fluxes and are fully consistent with the Euclidean expectation. We constrain the intrinsic luminosity function of AGN in the 2-10 keV energy band where the unabsorbed X-ray flux is estimated from the X-ray spectral fit up to z=3. Finally, we demonstrate the presence of a supercluster size structure at redshift 0.14, identified by means of percolation analysis of the XXL-1000-AGN sample. The XXL survey, reaching a medium flux limit and covering a wide area is a stepping stone between current deep fields and planned wide area surveys.

We study the acceleration of cosmic rays by collisionless structure formation shocks with ENZO grid simulations. Data from the FERMI satellite enable the use of galaxy clusters as a testbed for particle acceleration models. Based on advanced cosmological simulations that include different prescriptions for gas and cosmic rays physics, we use the predicted {\gamma}-ray emission to constrain the shock acceleration efficiency. We infer that the efficiency must be on average <0.1% for cosmic shocks, particularly for the 2<M<5 merger shocks that are mostly responsible for the thermalisation of the intracluster medium. These results emerge, both, from non-radiative and radiative runs including feedback from active galactic nuclei, as well as from zoomed resimulations of a cluster resembling MACSJ1752.0+0440. The limit on the acceleration efficiency we report is lower than what has been assumed in the literature so far. Combined with the information from radio emission in clusters, it appears that a revision of the present understanding of shock acceleration in the ICM is unavoidable.

In the current epoch, one of the main mechanisms driving the growth of galaxy clusters is the continuous accretion of group-scale halos. In this process, the ram pressure applied by the hot intracluster medium on the gas content of the infalling group is responsible for stripping the gas from its dark-matter halo, which gradually leads to the virialization of the infalling gas in the potential well of the main cluster. Using deep wide-field observations of the poor cluster Hydra A/A780 with XMM-Newton and Suzaku, we report the discovery of an infalling galaxy group 1.1 Mpc south of the cluster core. The presence of a substructure is confirmed by a dynamical study of the galaxies in this region. A wake of stripped gas is trailing behind the group over a projected scale of 760 kpc. The temperature of the gas along the wake is constant at kT ~ 1.3 keV, which is about a factor of two less than the temperature of the surrounding plasma. We observe a cold front pointing westwards compared to the peak of the group, which indicates that the group is currently not moving in the direction of the main cluster, but is moving along an almost circular orbit. The overall morphology of the group bears remarkable similarities with high-resolution numerical simulations of such structures, which greatly strengthens our understanding of the ram-pressure stripping process.

We present a high-precision mass model of the galaxy cluster MACSJ1149.6+2223, based on a strong-gravitational-lensing analysis of Hubble Space Telescope Frontier Fields (HFF) imaging data and spectroscopic follow-up with Gemini/GMOS and VLT/MUSE. Our model includes 12 new multiply imaged galaxies, bringing the total to 22, comprised of 65 individual lensed images. Unlike the first two HFF clusters, Abell 2744 and MACSJ0416.1-2403, MACSJ1149 does not reveal as many multiple images in the HFF data. Using the Lenstool software package and the new sets of multiple images, we model the cluster with several cluster-scale dark-matter halos and additional galaxy-scale halos for the cluster members. Consistent with previous analyses, we find the system to be complex, composed of five cluster-scale halos. Their spatial distribution and lower mass, however, makes MACSJ1149 a less powerful lens. Our best-fit model predicts image positions with an RMS of 0.91". We measure the total projected mass inside a 200~kpc aperture as ($1.840\pm 0.006$)$\times 10^{14}$M$_{\odot}$, thus reaching again 1% precision, following our previous HFF analyses of MACSJ0416.1-2403 and Abell 2744. In light of the discovery of the first resolved quadruply lensed supernova, SN Refsdal, in one of the multiply imaged galaxies identified in MACSJ1149, we use our revised mass model to investigate the time delays and predict the rise of the next image between November 2015 and January 2016.

Context. The XXL Survey is the largest survey carried out by the XMM-Newton satellite and covers a total area of 50 square degrees distributed over two fields. It primarily aims at investigating the large-scale structures of the Universe using the distribution of galaxy clusters and active galactic nuclei as tracers of the matter distribution. Aims. This article presents the XXL bright cluster sample, a subsample of 100 galaxy clusters selected from the full XXL catalogue by setting a lower limit of $3\times 10^{-14}\,\mathrm{erg \,s^{-1}cm^{-2}}$ on the source flux within a 1$^{\prime}$ aperture. Methods. The selection function was estimated using a mixture of Monte Carlo simulations and analytical recipes that closely reproduce the source selection process. An extensive spectroscopic follow-up provided redshifts for 97 of the 100 clusters. We derived accurate X-ray parameters for all the sources. Scaling relations were self-consistently derived from the same sample in other publications of the series. On this basis, we study the number density, luminosity function, and spatial distribution of the sample. Results. The bright cluster sample consists of systems with masses between $M_{500}=7\times 10^{13}$ and $3\times 10^{14} M_\odot$, mostly located between $z=0.1$ and 0.5. The observed sky density of clusters is slightly below the predictions from the WMAP9 model, and significantly below the predictions from the Planck 2015 cosmology. In general, within the current uncertainties of the cluster mass calibration, models with higher values of $\sigma_8$ and/or $\Omega_m$ appear more difficult to accommodate. We provide tight constraints on the cluster differential luminosity function and find no hint of evolution out to $z\sim1$. We also find strong evidence for the presence of large-scale structures in the XXL bright cluster sample and identify five new superclusters.

We present the XXL Survey, the largest XMM programme totaling some 6.9 Ms to date and involving an international consortium of roughly 100 members. The XXL Survey covers two extragalactic areas of 25 deg2 each at a point-source sensitivity of ~ 5E-15 erg/sec/cm2 in the [0.5-2] keV band (completeness limit). The survey's main goals are to provide constraints on the dark energy equation of state from the space-time distribution of clusters of galaxies and to serve as a pathfinder for future, wide-area X-ray missions. We review science objectives, including cluster studies, AGN evolution, and large-scale structure, that are being conducted with the support of approximately 30 follow-up programmes. We describe the 542 XMM observations along with the associated multi-lambda and numerical simulation programmes. We give a detailed account of the X-ray processing steps and describe innovative tools being developed for the cosmological analysis. The paper provides a thorough evaluation of the X-ray data, including quality controls, photon statistics, exposure and background maps, and sky coverage. Source catalogue construction and multi-lambda associations are briefly described. This material will be the basis for the calculation of the cluster and AGN selection functions, critical elements of the cosmological and science analyses. The XXL multi-lambda data set will have a unique lasting legacy value for cosmological and extragalactic studies and will serve as a calibration resource for future dark energy studies with clusters and other X-ray selected sources. With the present article, we release the XMM XXL photon and smoothed images along with the corresponding exposure maps. The XMM XXL observation list (Table B.1) is available in electronic form at the CDS. The present paper is the first in a series reporting results of the XXL-XMM survey.

The XXL Survey is the largest homogeneous and contiguous survey carried out with XMM-Newton. Covering an area of 50 square degrees distributed over two fields, it primarily investigates the large-scale structures of the Universe using the distribution of galaxy clusters and active galactic nuclei as tracers of the matter distribution. Given its depth and sky coverage, XXL is particularly suited to systematically unveiling the clustering of X-ray clusters and to identifying superstructures in a homogeneous X-ray sample down to the typical mass scale of a local massive cluster. A friends-of-friends algorithm in three-dimensional physical space was run to identify large-scale structures. In this paper we report the discovery of the highest redshift supercluster of galaxies found in the XXL Survey. We describe the X-ray properties of the clusters members of the structure and the optical follow-up. The newly discovered supercluster is composed of six clusters of galaxies at a median redshift z around 0.43 and distributed across approximately 30 by 15 arc minutes (10 by 5 Mpc on sky) on the sky. This structure is very compact with all the clusters residing in one XMM pointing; for this reason this is the first supercluster discovered with the XXL Survey. Spectroscopic follow-up with WHT (William Herschel Telescope) and NTT (New Technology Telescope) confirmed a median redshift of z = 0.43. An estimate of the X-ray mass and luminosity of this supercluster and of its total gas mass put XLSSC-e at the average mass range of superclusters; its appearance, with two members of equal size, is quite unusual with respect to other superclusters and provides a unique view of the formation process of a massive structure.

Traditionally, galaxy clusters have been expected to retain all the material accreted since their formation epoch. For this reason, their matter content should be representative of the Universe as a whole, and thus their baryon fraction should be close to the Universal baryon fraction. We make use of the sample of the 100 brightest galaxy clusters discovered in the XXL Survey to investigate the fraction of baryons in the form of hot gas and stars in the cluster population. We measure the gas masses of the detected halos and use a mass--temperature relation directly calibrated using weak-lensing measurements for a subset of XXL clusters to estimate the halo mass. We find that the weak-lensing calibrated gas fraction of XXL-100-GC clusters is substantially lower than was found in previous studies using hydrostatic masses. Our best-fit relation between gas fraction and mass reads $f_{\rm gas,500}=0.055_{-0.006}^{+0.007}\left(M_{\rm 500}/10^{14}M_\odot\right)^{0.21_{-0.10}^{+0.11}}$. The baryon budget of galaxy clusters therefore falls short of the Universal baryon fraction by about a factor of two at $r_{\rm 500}$. Our measurements require a hydrostatic bias $1-b=M_X/M_{\rm WL}=0.72_{-0.07}^{+0.08}$ to match the gas fraction obtained using lensing and hydrostatic equilibrium. Comparing our gas fraction measurements with the expectations from numerical simulations, our results favour an extreme feedback scheme in which a significant fraction of the baryons are expelled from the cores of halos. This model is, however, in contrast with the thermodynamical properties of observed halos, which might suggest that weak-lensing masses are overestimated. We note that a mass bias $1-b=0.58$ as required to reconcile Planck CMB and cluster counts should translate into an even lower baryon fraction, which poses a major challenge to our current understanding of galaxy clusters. [Abridged]

The XXL survey is the largest survey carried out by XMM-Newton. Covering an area of 50deg$^2$, the survey contains $\sim450$ galaxy clusters out to a redshift $\sim$2 and to an X-ray flux limit of $\sim5\times10^{-15}erg\,s^{-1}cm^{-2}$. This paper is part of the first release of XXL results focussed on the bright cluster sample. We investigate the scaling relation between weak-lensing mass and X-ray temperature for the brightest clusters in XXL. The scaling relation is used to estimate the mass of all 100 clusters in XXL-100-GC. Based on a subsample of 38 objects that lie within the intersection of the northern XXL field and the publicly available CFHTLenS catalog, we derive the $M_{WL}$ of each system with careful considerations of the systematics. The clusters lie at $0.1<z<0.6$ and span a range of $ T\simeq1-5keV$. We combine our sample with 58 clusters from the literature, increasing the range out to 10keV. To date, this is the largest sample of clusters with $M_{WL}$ measurements that has been used to study the mass-temperature relation. The fit ($M\propto T^b$) to the XXL clusters returns a slope $b=1.78^{+0.37}_{-0.32}$ and intrinsic scatter $\sigma_{\ln M|T}\simeq0.53$; the scatter is dominated by disturbed clusters. The fit to the combined sample of 96 clusters is in tension with self-similarity, $b=1.67\pm0.12$ and $\sigma_{\ln M|T}\simeq0.41$. Overall our results demonstrate the feasibility of ground-based weak-lensing scaling relation studies down to cool systems of $\sim1keV$ temperature and highlight that the current data and samples are a limit to our statistical precision. As such we are unable to determine whether the validity of hydrostatic equilibrium is a function of halo mass. An enlarged sample of cool systems, deeper weak-lensing data, and robust modelling of the selection function will help to explore these issues further.

Observations of the cosmic microwave background indicate that baryons account for 5% of the Universe's total energy content. In the local Universe, the census of all observed baryons falls short of this estimate by a factor of two. Cosmological simulations indicate that the missing baryons might not have condensed into virialized haloes, but reside throughout the filaments of the cosmic web (where matter density is larger than average) as a low-density plasma at temperatures of $10^5-10^7$ kelvin, known as the warm-hot intergalactic medium. There have been previous claims of the detection of warm baryons along the line of sight to distant blazars and of hot gas between interacting clusters. These observations were, however, unable to trace the large-scale filamentary structure, or to estimate the total amount of warm baryons in a representative volume of the Universe. Here we report X-ray observations of filamentary structures of gas at $10^7$ kelvin associated with the galaxy cluster Abell 2744. Previous observations of this cluster were unable to resolve and remove coincidental X-ray point sources. After subtracting these, we reveal hot gas structures that are coherent over scales of 8 mergaparsecs. The filaments coincide with over-densities of galaxies and dark matter, with 5-10% of their mass in baryonic gas. This gas has been heated up by the cluster's gravitational pull and is now feeding its core. Our findings strengthen evidence for a picture of the Universe in which a large fraction of the missing baryons reside in the filaments of the cosmic web.

We want to characterize the dynamical state of galaxy clusters detected with the Sunyaev-Zeldovich (SZ) effect by Planck and compare them with the dynamical state of clusters selected in X-rays survey. We analyzed a representative subsample of the Planck SZ catalogue, containing the 132 clusters with the highest signal to noise ratio and characterize their dynamical state using as indicator the projected offset between the peak of the X-ray emission and the position of the Brightest cluster galaxy. We study the distribution of our indicator in our sample and compare it to its distribution in X-ray selected samples (HIFLUGCS, MACS and REXCESS). The distributions are significantly different and the fraction of relaxed objects is smaller in the Planck sample ($52 \pm 4 \%$) than in X-ray samples ($\simeq 74\%$) We interpret this result as an indication of different selection effects affecting X-rays (e.g. "cool core bias") and SZ surveys of galaxy clusters.

Our goal is to provide a robust estimate of the metal content of the ICM in massive clusters. We make use of published abundance profiles for a sample of ~60 nearby systems, we include in our estimate uncertainties associated to the measurement process and to the almost total lack of measures in cluster outskirts. We perform a first, albeit rough, census of metals finding that the mean abundance of the ICM within r_180 is very poorly constrained, 0.06Z_sol < Z < 0.26Z_sol, and presents no tension with expectations. Similarly, the question of if and how the bulk of the metal content in clusters varies with cosmic time, is very much an open one. A solid estimate of abundances in cluster outskirts could be achieved by combining observations of the two experiments which will operate on board Athena, the XIFU and the WFI, provided they do not fall victim to the de-scoping process that has afflicted several space observatories over the last decade.

The complex X-ray system IGR J11014-6103 (a.k.a. the Lighthouse nebula) is composed of a bow-shock pulsar wind nebula (PWN) as well as large-scale jet-like features, all launched by IGR J11014-6103 which is moving supersonically in the interstellar medium. Previous observations suggested that the jet features stem from a ballistic jet of relativistic particles. In order to confirm the nature of the jet and the marginally detected counter-jet, we obtained a new deep 250 ks Chandra observation of the Lighthouse nebula. We performed detailed spatial and spectral analysis of all X-ray components of the system. The X-ray PWN is now better resolved and shows a clear bi-modal morphology. The overall helical pattern of the main jet is confirmed. However, there are large deviations from a simple helical model at small and large scales. Significant extended emission is now detected, encompassing the main jet all along its length. The brightness dip of the main jet at ~50" distance from the pulsar is confirmed, the extended emission however prevents conclusions about the coherence at this position of the jet. The counter-jet is now detected at high statistical significance. In addition, we found two small-scale "arcs" departing from the pulsar towards the jets. We also looked for possible bow-shock emission in front of the pulsar, with a short VLT/FORS2 H-alpha observation. No clear emission is found, most likely because of contamination from a diffuse surrounding nebulosity. The results of our X-ray analysis show that current expectations from a ballistic nature of the jets can explain satisfactorily some of the observational evidences but cannot fully reproduce the observations. The alternative scenario (diffusion of particles along preexisting magnetic field lines in the surrounding medium) however also continues to suffer from conflicts with the observations.

We measure the thermodynamic properties of cluster outer regions to provide constraints on the processes that rule the formation of large scale structures. We derived the thermodynamic properties of the intracluster gas (temperature, entropy) by combining the SZ thermal pressure from Planck and the X-ray gas density from ROSAT. This method allowed us to reconstruct for the first time temperature and entropy profiles out to the virial radius and beyond in a large sample of objects. At variance with several recent Suzaku studies, we find that the entropy rises steadily with radius, albeit at at a somewhat lower rate than predicted by self-similar expectations. We note significant differences between relaxed, cool-core systems and unrelaxed clusters in the outer regions. Relaxed systems appear to follow the self-similar expectations more closely than perturbed objects. Our results indicate that the well-known entropy excess observed in cluster cores extends well beyond the central regions. When correcting for the gas depletion, the observed entropy profiles agree with the prediction from gravitational collapse only, especially for cool-core clusters.

Several recent studies used the hot gas fraction of galaxy clusters as a standard ruler to constrain dark energy, which provides competitive results compared to other techniques. This method, however, relies on the assumption that the baryon fraction in clusters agrees with the cosmic value Omega_b/Omega_m, and does not differ from one system to another. We test this hypothesis by measuring the gas mass fraction over the entire cluster volume in a sample of local clusters. Combining the SZ thermal pressure from Planck and the X-ray gas density from ROSAT, we measured for the first time the average gas fraction (fgas) out to the virial radius and beyond in a large sample of clusters. We also obtained azimuthally-averaged measurements of the gas fraction for 18 individual systems, which we used to compute the scatter of fgas around the mean value at different radii and its dependence on the cluster's temperature. The gas mass fraction increases with radius and reaches the cosmic baryon fraction close to R200. At R200, we measure fgas,200=0.176+/-0.009. We find significant differences between the baryon fraction of relaxed, cool-core (CC) systems and unrelaxed, non-cool core (NCC) clusters in the outer regions. In average, the gas fraction in NCC clusters slightly exceeds the cosmic baryon fraction, while in CC systems the gas fraction converges to the expected value when accounting for the stellar content, without any evidence for variations from one system to another. We find that fgas estimates in NCC systems slightly disagree with the cosmic value approaching R200. This result could be explained either by a violation of the assumption of hydrostatic equilibrium or by an inhomogeneous distribution of the gas mass. Conversely, cool-core clusters are found to provide reliable constraints on fgas at overdensities >200, which makes them suitable for cosmological studies.

We present the analysis of a local (z = 0.04 - 0.2) sample of 31 galaxy clusters with the aim of measuring the density of the X-ray emitting gas in cluster outskirts. We compare our results with numerical simulations to set constraints on the azimuthal symmetry and gas clumping in the outer regions of galaxy clusters. We exploit the large field-of-view and low instrumental background of ROSAT/PSPC to trace the density of the intracluster gas out to the virial radius. We perform a stacking of the density profiles to detect a signal beyond r200 and measure the typical density and scatter in cluster outskirts. We also compute the azimuthal scatter of the profiles with respect to the mean value to look for deviations from spherical symmetry. Finally, we compare our average density and scatter profiles with the results of numerical simulations. As opposed to some recent Suzaku results, and confirming previous evidence from ROSAT and Chandra, we observe a steepening of the density profiles beyond \sim r500. Comparing our density profiles with simulations, we find that non-radiative runs predict too steep density profiles, whereas runs including additional physics and/or treating gas clumping are in better agreement with the observed gas distribution. We report for the first time the high-confidence detection of a systematic difference between cool-core and non-cool core clusters beyond \sim 0.3r200, which we explain by a different distribution of the gas in the two classes. Beyond \sim r500, galaxy clusters deviate significantly from spherical symmetry, with only little differences between relaxed and disturbed systems. We find good agreement between the observed and predicted scatter profiles, but only when the 1% densest clumps are filtered out in the simulations. [Abridged]

Several recent works have reported the detection of an unidentified X-ray line at 3.55~keV, which could possibly be attributed to the decay of dark matter (DM) particles in the halos of galaxy clusters and in the M31 galaxy. We analyze all publicly-available XMM-Newton data of dwarf spheroidal galaxies to test the possible DM origin of the line. Dwarf spheroidal galaxies have high mass-to-light ratios and their interstellar medium is not a source of diffuse X-ray emission, thus they are expected to provide the cleanest DM decay line signal. Our analysis shows no evidence for the presence of the line in the stacked spectra of the dwarf galaxies. It excludes the sterile neutrino DM decay origin of the 3.5 keV line reported by Bulbul et al. (2014) at the level of 4.1sigma under standard assumptions about the Galactic DM column density in the direction of selected dwarf galaxies and at the level of 3.2sigma assuming minimal Galactic DM column density. Our analysis is still consistent with the estimate of sterile neutrino DM parameters by Boyarsky et al. (2014), because of its larger uncertainty. However, the central value of their estimate of the mixing angle is inconsistent with our dwarf spheroidals data at 3.4sigma (2.5sigma) level assuming the mean (minimal) Galactic DM column density. As a by-product of our analysis, we provide updated upper limits to the mixing angle of sterile neutrino DM in the mass range between 2 and 20 keV.

We present a high-precision mass model of galaxy cluster Abell 2744, based on a strong-gravitational-lensing analysis of the \emph{Hubble Space Telescope Frontier Fields} (HFF) imaging data, which now include both \emph{Advanced Camera for Surveys} and \emph{Wide-Field Camera 3} observations to the final depth. Taking advantage of the unprecedented depth of the visible and near-infrared data, we identify 34 new multiply imaged galaxies, bringing the total to 61, comprising 181 individual lensed images. In the process, we correct previous erroneous identifications and positions of multiple systems in the northern part of the cluster core. With the \textsc{Lenstool} software and the new sets of multiple images, we model the cluster using two cluster-scale dark matter halos plus galaxy-scale halos for the cluster members. Our best-fit model predicts image positions with an \emph{RMS} error of 0.69$\arcsec$, which constitutes an improvement by almost a factor of two over previous parametric models of this cluster. We measure the total projected mass inside a 200~kpc aperture as ($2.162\pm 0.005$)$\times 10^{14}M_{\odot}$, thus reaching 1\% level precision for the second time, following the recent HFF measurement of MACSJ0416.1-2403. Importantly, the higher quality of the mass model translates into an overall improvement by a factor of 4 of the derived magnification factor. % for the high-redshift lensed background galaxies. Together with our previous HFF gravitational lensing analysis, this work demonstrates that the HFF data enables high-precision mass measurements for massive galaxy clusters and the derivation of robust magnification maps to probe the early Universe.

Radio relics in galaxy clusters are associated with powerful shocks that (re)accelerate relativistic electrons. It is widely believed that the acceleration proceeds via diffusive shock acceleration. In the framework of thermal leakage, the ratio of the energy in relativistic electrons to the energy in relativistic protons should should be smaller than $K_{\rm e/p} \sim 10^{-2}$. The relativistic protons interact with the thermal gas to produce $\gamma$-rays in hadronic interactions. Combining observations of radio relics with upper limits from $\gamma$-ray observatories can constrain the ratio $K_{\rm e/p}$. In this work we selected 10 galaxy clusters that contain double radio relics, and derive new upper limits from the stacking of $\gamma$-ray observations by FERMI. We modelled the propagation of shocks using a semi-analytical model, where we assumed a simple geometry for shocks and that cosmic ray protons are trapped in the intracluster medium. Our analysis shows that diffusive shock acceleration has difficulties in matching simultaneously the observed radio emission and the constraints imposed by FERMI, unless the magnetic field in relics is unrealistically large ($\gg 10 ~\rm \mu G$). In all investigated cases (also including realistic variations of our basic model and the effect of re-acceleration) the mean emission of the sample is of the order of the stacking limit by FERMI, or larger. These findings put tension on the commonly adopted model for the powering of radio relics, and imply that the relative acceleration efficiency of electrons and protons is at odds with predictions of diffusive shock acceleration, requiring $K_{\rm e/p} \geq 10-10^{-2}$.

Context. Galaxy clusters can be used as cosmological probes, but to this end, they need to be thoroughly understood. Combining all cluster observables in a consistent way will help us to understand their global properties and their internal structure. Aims. We provide proof of the concept that the projected gravitational potential of galaxy clusters can directly be reconstructed from X-ray observations. We also show that this joint analysis can be used to locally test the validity of the equilibrium assumptions in galaxy clusters. Methods. We used a newly developed reconstruction method, based on Richardson-Lucy deprojection, that allows reconstructing projected gravitational potentials of galaxy clusters directly from X-ray observations. We applied this algorithm to the well-studied cluster Abell 1689 and compared the gravitational potential reconstructed from X-ray observables to the potential obtained from gravitational lensing measurements. [...] Results. Assuming spherical symmetry and hydrostatic equilibrium, the potentials recovered from gravitational lensing and from X-ray emission agree very well beyond 500 kpc. Owing to the fact that the Richardson-Lucy deprojection algorithm allows deprojecting each line of sight independently, this result may indicate that non-gravitational effects and/or asphericity are strong in the central regions of the clusters. Conclusions. We demonstrate the robustness of the potential reconstruction method based on the Richardson-Lucy deprojection algorithm and show that gravitational lensing and X-ray emission lead to consistent gravitational potentials. Our results illustrate the power of combining galaxy-cluster observables in a single, non-parametric, joint reconstruction of consistent cluster potentials that can be used to locally constrain the physical state of the gas.

The reconstruction of galaxy cluster's gas density profiles is usually performed by assuming spherical symmetry and averaging the observed X-ray emission in circular annuli. In the case of a very inhomogeneous and asymmetric gas distribution, this method has been shown to return biased results in numerical simulations because of the $n^2$ dependence of the X-ray emissivity. We propose a method to recover the true density profiles in the presence of inhomogeneities, based on the derivation of the azimuthal median of the surface brightness in concentric annuli. We demonstrate the performance of this method with numerical simulations, and apply it to a sample of 31 galaxy clusters in the redshift range 0.04-0.2 observed with ROSAT/PSPC. The clumping factors recovered by comparing the mean and the median are mild and show a slight trend of increasing bias with radius. For $R<R_{500}$, we measure a clumping factor $\sqrt{C}<1.1$, which indicates that the thermodynamic properties and hydrostatic masses measured in this radial range are only mildly affected by this effect. Comparing our results with three sets of hydrodynamical numerical simulations, we found that non-radiative simulations significantly overestimate the level of inhomogeneities in the ICM, while the runs including cooling, star formation, and AGN feedback reproduce the observed trends closely. Our results indicate that most of the accretion of X-ray emitting gas is taking place in the diffuse, large-scale accretion patterns rather than in compact structures.

We use a joint optical/X-ray analysis to constrain the geometry and history of the ongoing merging event in the massive galaxy cluster MACSJ0416.1-2403 (z=0.397). Our investigation of cluster substructure rests primarily on a combined strong- and weak-lensing mass reconstruction based on the deep, high-resolution images obtained for the Hubble Frontier Fields initiative. To reveal the system's dynamics, we complement this lensing analysis with a study of the intra-cluster gas using shallow Chandra data, and a three-dimensional model of the distribution and motions of cluster galaxies derived from over 100 spectroscopic redshifts. The multi-scale grid model obtained from our combined lensing analysis extends the high-precision strong-lensing mass reconstruction recently performed to cluster-centric distances of almost 1 Mpc. Our analysis detects the two well known mass concentrations in the cluster core. A pronounced offset between collisional and collisionless matter is only observed for the SW cluster component, while excellent alignment is found for the NE cluster. Both the lensing analysis and the distribution of cluster light strongly suggest the presence of a third massive structure, almost 2 arcmin SW of the cluster centre. Since no X-ray emission is detected in this region, we conclude that this structure is non-virialised and speculate that it might be part of a large-scale filament almost aligned with our line of sight. Combining all evidence from the distribution of dark and luminous matter, we propose two alternative scenarios for the trajectories of the components of MACSJ0416.1-2403. Upcoming deep X-ray observations that allow the detection of shock fronts, cold cores, and sloshing gas (all key diagnostics for studies of cluster collisions) will allow us to test, and distinguish between these two scenarios.

Structure formation in the current Universe operates through the accretion of group-scale systems onto massive clusters. The detection and study of such accreting systems is crucial to understand the build-up of the most massive virialized structures we see today. We report the discovery with XMM-Newton of an irregular X-ray substructure in the outskirts of the massive galaxy cluster Abell 2142. The tip of the X-ray emission coincides with a concentration of galaxies. The bulk of the X-ray emission of this substructure appears to be lagging behind the galaxies and extends over a projected scale of at least 800 kpc. The temperature of the gas in this region is 1.4 keV, which is a factor of ~4 lower than the surrounding medium and is typical of the virialized plasma of a galaxy group with a mass of a few 10^13M_sun. For this reason, we interpret this structure as a galaxy group in the process of being accreted onto the main dark-matter halo. The X-ray structure trailing behind the group is due to gas stripped from its original dark-matter halo as it moves through the intracluster medium (ICM). This is the longest X-ray trail reported to date. For an infall velocity of ~1,200 km s-1 we estimate that the stripped gas has been surviving in the presence of the hot ICM for at least 600 Myr, which exceeds the Spitzer conduction timescale in the medium by a factor of >~400. Such a strong suppression of conductivity is likely related to a tangled magnetic field with small coherence length and to plasma microinstabilities. The long survival time of the low-entropy intragroup medium suggests that the infalling material can eventually settle within the core of the main cluster.

We present a high-precision mass model of the galaxy cluster MACSJ0416.1-2403, based on a strong-gravitational-lensing analysis of the recently acquired Hubble Space Telescope Frontier Fields (HFF) imaging data. Taking advantage of the unprecedented depth provided by HST/ACS observations in three passbands, we identify 51 new multiply imaged galaxies, quadrupling the previous census and bringing the grand total to 68, comprising 194 individual lensed images. Having selected a subset of the 57 most securely identified multiply imaged galaxies, we use the Lenstool software package to constrain a lens model comprised of two cluster-scale dark-matter halos and 98 galaxy-scale halos. Our best-fit model predicts image positions with an $RMS$ error of 0.68'', which constitutes an improvement of almost a factor of two over previous, pre-HFF models of this cluster. We find the total projected mass inside a 200~kpc aperture to be $(1.60\pm0.01)\times 10^{14}\ M_\odot$, a measurement that offers a three-fold improvement in precision, reaching the percent level for the first time in any cluster. Finally, we quantify the increase in precision of the derived gravitational magnification of high-redshift galaxies and find an improvement by a factor of $\sim$2.5 in the statistical uncertainty. Our findings impressively confirm that HFF imaging has indeed opened the domain of high-precision mass measurements for massive clusters of galaxies.

We report on the X-ray observation of a strong lensing selected group, SL2S J08544-0121, with a total mass of $2.4 \pm 0.6 \times 10^{14}$ $\rm{M_\odot}$ which revealed a separation of $124\pm20$ kpc between the X-ray emitting collisional gas and the collisionless galaxies and dark matter (DM), traced by strong lensing. This source allows to put an order of magnitude estimate to the upper limit to the interaction cross section of DM of 10 cm$^2$ g$^{-1}$. It is the lowest mass object found to date showing a DM-baryons separation and it reveals that the detection of bullet-like objects is not rare and confined to mergers of massive objects opening the possibility of a statistical detection of DM-baryons separation with future surveys.

Observations in radio have shown that galaxy clusters are giant reservoirs of cosmic rays (CR). Although a gamma- ray signal from the cluster volume is expected to arise through interactions of CR protons with the ambient plasma, a confirming observation is still missing. We search for a cumulative gamma-ray emission in direction of galaxy clusters by analysing a collection of stacked Fermi-LAT count maps. Additionally, we investigate possible systematic differences in the emission between cool-core and non-cool core cluster populations. Making use of a sample of 53 clusters selected from the HIFLUGCS catalog, we do not detect a significant signal from the stacked sample. The upper limit on the average flux per cluster derived for the total stacked sample is at the level of a few 1e-11 ph cm-2 s-1 at 95% confidence level in the 1-300 GeV band, assuming power-law spectra with photon indices 2.0, 2.4, 2.8 and 3.2. Separate stacking of the cool core and non-cool core clusters in the sample lead to similar values of around 5e-11 ph cm-2 s-1 and 2e-11 ph cm-2 s-1, respectively. Under the assumption that decaying pi0, produced in collisions between CRs and the ambient thermal gas, are responsible for the gamma-ray emission, we set upper limits on the average CR content in galaxy clusters. For the entire cluster population, our upper limit on the gamma-ray flux translates into an upper limit on the average CR-to-thermal energy ratio of 4.6% for a photon index of 2.4, although it is possible for individual systems to exceed this limit. Our 95% upper limits are at the level expected from numerical simulations, which likely suggests that the injection of CR at cosmological shocks is less efficient than previously assumed.

As the nodes of the cosmic web, clusters of galaxies trace the large-scale distribution of matter in the Universe. They are thus privileged sites in which to investigate the complex physics of structure formation. However, the complete story of how these structures grow, and how they dissipate the gravitational and non-thermal components of their energy budget over cosmic time, is still beyond our grasp. Fundamental questions such as How do hot diffuse baryons accrete and dynamically evolve in dark matter potentials? How and when was the energy that we observe in the ICM generated and distributed? Where and when are heavy elements produced and how are they circulated? are still unanswered. Most of the cluster baryons exists in the form of a diffuse, hot, metal-enriched plasma that radiates primarily in the X-ray band (the intracluster medium, ICM), allowing the X-ray observations of the evolving cluster population to provide a unique opportunity to address these topics. Athena+ with its large collecting area and unprecedented combination of high spectral and angular resolution offers the only way to make major advances in answering these questions. Athena+ will show how the baryonic gas evolves in the dark matter potential wells by studying the motions and turbulence in the ICM. Athena+ will be able to resolve the accreting region both spatially and spectroscopically, probing the true nature and physical state of the X-ray emitting plasma. Athena+ has the capabilities to permit a definitive understanding of the formation and evolution of large-scale cosmic structure through the study of the cluster population.

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.

We present results obtained with a new XMM-Newton observation of A2142, a famous textbook example of cluster with multiple cold fronts, which has been studied in detail with Chandra but whose large scale properties are presented here for the first time. We report the discovery of a a new cold front, the most distant one ever detected in a galaxy cluster, at about one Mpc from the center to the SE. Residual images, thermodynamics and metal abundance maps are qualitatively in agreement with predictions from numerical simulations of the sloshing phenomenon. However, the scales involved are much larger, similarly to what recently observed in the Perseus cluster. These results show that sloshing is a cluster-wide phenomenon, not confined in the cores, which extends well beyond the cooling region involving a large fraction of the ICM up to almost half of the virial radius. The absence of a cool core and a newly discovered giant radio halo in A2142, in spite of its relaxed X-ray morphology, suggest that large scale sloshing, or the intermediate merger which caused it, may trigger Mpc-scale radio emission and may lead to the disruption of the cluster cool core

We present a combined X-ray, optical, and radio analysis of the galaxy group IC 1860 using the currently available Chandra and XMM data, literature multi-object spectroscopy data and GMRT data. The Chandra and XMM imaging and spectroscopy reveal two surface brightness discontinuities at 45 and 76 kpc shown to be consistent with a pair of cold fronts. These features are interpreted as due to sloshing of the central gas induced by an off-axis minor merger with a perturber. This scenario is further supported by the presence of a peculiar velocity of the central galaxy IC 1860 and the identification of a possible perturber in the optically disturbed spiral galaxy IC 1859. The identification of the perturber is consistent with the comparison with numerical simulations of sloshing. The GMRT observation at 325 MHz shows faint, extended radio emission contained within the inner cold front, as seen in some galaxy clusters hosting diffuse radio mini-halos. However, unlike mini-halos, no particle reacceleration is needed to explain the extended radio emission, which is consistent with aged radio plasma redistributed by the sloshing. There is strong analogy of the X-ray and optical phenomenology of the IC 1860 group with two other groups, NGC 5044 and NGC 5846, showing cold fronts. The evidence presented in this paper is among the strongest supporting the currently favored model of cold-front formation in relaxed objects and establishes the group scale as a chief environment to study this phenomenon.

The spatial distribution of gas matter inside galaxy clusters is not completely smooth, but may host gas clumps associated with substructures. These overdense gas substructures are generally a source of unresolved bias of X-ray observations towards high density gas, but their bright luminosity peaks may be resolved sources within the ICM, that deep X-ray exposures may be (already) capable to detect. In this paper we aim at investigating both features, using a set of high-resolution cosmological simulations with ENZO. First, we monitor how the bias by unresolved gas clumping may yield incorrect estimates of global cluster parameters and affects the measurements of baryon fractions by X-ray observations. We find that based on X-ray observations of narrow radial strips, it is difficult to recover the real baryon fraction to better than 10 - 20 percent uncertainty. Second, we investigated the possibility of observing bright X-ray clumps in the nearby Universe (z<=0.3). We produced simple mock X-ray observations for several instruments (XMM, Suzaku and ROSAT) and extracted the statistics of potentially detectable bright clumps. Some of the brightest clumps predicted by simulations may already have been already detected in X- ray images with a large field of view. However, their small projected size makes it difficult to prove their existence based on X-ray morphology only. Preheating, AGN feedback and cosmic rays are found to have little impact on the statistical properties of gas clumps.

X-ray astronomers often divide galaxy clusters into two classes: "cool core" (CC) and "non-cool core" (NCC) objects. The origin of this dichotomy has been the subject of debate in recent years, between "evolutionary" models (where clusters can evolve from CC to NCC, mainly through mergers) and "primordial" models (where the state of the cluster is fixed "ab initio" by early mergers or pre-heating). We found that in a well-defined sample (clusters in the GMRT Radio halo survey with available Chandra or XMM-Newton data), none of the objects hosting a giant radio halo can be classified as a cool core. This result suggests that the main mechanisms which can start a large scale synchrotron emission (most likely mergers) are the same that can destroy CC and therefore strongly supports "evolutionary" models of the CC-NCC dichotomy. Moreover combining the number of objects in the CC and NCC state with the number of objects with and without a radio-halo, we estimated that the time scale over which a NCC cluster relaxes to the CC state, should be larger than the typical life-time of radio-halos and likely shorter than about 3 Gyr. This suggests that NCC transform into CC more rapidly than predicted from the cooling time, which is about 10 Gyr in NCC systems, allowing the possibility of a cyclical evolution between the CC and NCC states.

We wish to reassess the properties of the ICM at large radii in the galaxy cluster PKS 0745-191 in light of the recent Suzaku measurements. We analyze an archival 10.5 ksec ROSAT/PSPC observation to extract the surface-brightness profile of PKS 0745-191 and infer the deprojected density profile. We then compare the ROSAT surface-brightness profile with the Suzaku result. We perform a mass analysis combining the ROSAT density profile and the published temperature profiles from different instruments. We find that the ROSAT surface-brightness profile is statistically inconsistent (7.7 sigma) with the Suzaku result around and beyond the value of r200 estimated by Suzaku. We argue that, thanks to its large field of view and low background, ROSAT/PSPC is to the present day the most sensitive instrument to low surface-brightness X-ray emission in the 0.4-2.0 keV band. We also note that the Suzaku temperature and mass profiles are at odds with the results from at least two other satellites (XMM-Newton and Swift). The difference in surface brightness between ROSAT and Suzaku is most likely explained by the existence of additional foreground components at the low Galactic latitude of the source, which were not taken into account in the Suzaku background modeling. In light of our mass analysis, we conclude that any estimate of the fraction of the virial radius reached by X-ray measures is affected by systematic errors of the order of 25%. As a result, the properties of the ICM at the virial radius are still uncertain, and the Suzaku results should be considered with caution.

Diffuse radio emission in galaxy clusters, and their connection with cluster mergers, are still debated. We seek to explore the internal dynamics of the radio halo cluster Abell 545. This cluster is also peculiar for hosting in its center a very bright, red, diffuse intracluster light due to an old, stellar population, so bright to be named as "star pile". Our analysis is based on redshift data for 110 galaxies. We identify 95 cluster members and analyze the cluster internal dynamics by combining galaxy velocities and positions. We also use both photometric and X-ray data. We estimate the cluster redshift, z=0.1580, a velocity dispersion of 1200 km/s, and ICM temperature kT_X~8 keV. Our optical and X-ray analyses detect substructures. Optical data reveal three main galaxy clumps (center, NNW, and NE), and possibly a fourth clump at South. There is not a dominant galaxy and the four brightest galaxies avoid the cluster core (>~0.4h distant from the cluster center) and are >~1500 km/s far from the mean cluster velocity. The analysis of the X-ray surface brightness distribution provides us evidence of a disturbed dynamical phase. Located in the star pile region there is the brightest galaxies of the cluster core (CBCG) and a very compact elliptical galaxy. We show that the star pile has a similar redshift to that of the CBCG. Both the star pile and the CBCG are at rest in the cluster rest frame. The emerging picture of Abell 545 is that of a massive, M(R<1.6 h_70^-1 Mpc)=1.1-1.8x10^15 h_70^-1 Msun, very complex cluster with merging occurring along two directions. A545 gives another proof in the favor of the connection between cluster merger and extended, diffuse radio emission. The star pile, likely due to the process of a brightest galaxy forming in the cluster core. A545 represents a textbook cluster where to study the simultaneous formation of a galaxy system and its brightest galaxy.

From the point of view of X-ray astronomers, galaxy clusters are usually divided into two classes: "cool core"(CC) and "non-cool core" (NCC) objects. The origin of this dichotomy has been subject of debate in recent years, between "evolutionary" models (where clusters can evolve from CC to NCC, mainly through mergers) and "primordial" models (where the state of the cluster is fixed "ab initio" by early mergers or pre-heating). We found that in a representative sample (clusters in the GMRT Radio halo survey with available X-ray data), none of the objects hosting a giant radio halo can be classified as a cool core. This result suggests that the main mechanisms which can produce the ingredients to start a large scale synchrotron emission (most likely mergers) are the same that can destroy CC and therefore strongly supports "evolutionary" models of the CC-NCC dichotomy.

When selecting flux-limited cluster samples, the detection efficiency of X-ray instruments is not the same for centrally-peaked and flat objects, which introduces a bias in flux-limited cluster samples. We quantify this effect in the case of a well-known cluster sample, HIFLUGCS. We simulate a population of X-ray clusters with various surface-brightness profiles, and use the instrumental characteristics of the ROSAT All-Sky Survey (RASS) to select flux-limited samples similar to the HIFLUGCS sample and predict the expected bias. For comparison, we also estimate observationally the bias in the HIFLUGCS sample using XMM-Newton and ROSAT data. We find that the selection of X-ray cluster samples is significantly biased ($\sim29%$) in favor of the peaked, Cool-Core (CC) objects, with respect to Non-Cool-Core (NCC) systems. Interestingly, we find that the bias affects the low-mass, nearby objects (groups, poor clusters) much more than the more luminous objects (i.e massive clusters). We also note a moderate increase of the bias for the more distant systems. Observationally, we propose to select the objects according to their flux in a well-defined physical range excluding the cores, $0.2r_{500}-r_{500}$, to get rid of the bias. From the fluxes in this range, we reject 13 clusters out of the 64 in the HIFLUGCS sample, none of which appears to be NCC. As a result, we estimate that less than half (35-37%) of the galaxy clusters in the local Universe are strong CC. In the paradigm where the CC objects trace relaxed clusters as opposed to unrelaxed, merging objects, this implies that to the present day the majority of the objects are not in a relaxed state. From this result, we estimate a rate of heating events of $\sim1/3$ Gyr$^{-1}$ per dark-matter halo.

Extragalactic nuclear activity is best explored with observations at high energies, where the most extreme flux and spectral variations are expected to occur, witnessing changes in the accretion flow or in the kinematics of the plasma. In active galactic nuclei of blazar type, these variations are the most dramatic. By following blazar outbursts from their onset and by correlating the observed variations at many different wavelengths we can reconstruct the behavior of the plasma and map out the development of the flare within the jet. The advent of the Fermi satellite has allowed the start of a systematic and intensive monitoring program of blazars. Blazar outbursts are very effectively detected by the LAT instrument in the MeV-GeV domain, and these can be promptly followed up with other facilities. Based on a Fermi LAT detection of a high MeV-GeV state, we have observed the blazar PKS 1502+106 with the INTEGRAL satellite between 9 and 11 August 2008. Simultaneous Swift observations have been also accomplished, as well as optical follow-up at the Nordic Optical Telescope. The IBIS instrument onboard INTEGRAL detected a source at a position inconsistent with the optical coordinates of PKS 1502+106, but consistent with those of the Seyfert 1 galaxy Mkn 841, located at 6.8 arcmin south-west from the blazar, which is therefore responsible for all the hard X-ray flux detected by IBIS. At the location of the blazar, IBIS sets an upper limit of ~10^{-11} erg/s/cm2 on the 15-60 keV flux, that turns out to be consistent with a model of inverse Compton scattering accounting for the soft X-ray and gamma-ray spectra measured by Swift XRT and Fermi LAT, respectively. The gamma-ray spectrum during the outburst indicates substantial variability of the characteristic energy of the inverse Compton component in this blazar.(abridged).

Context: In 6 years of operation, INTEGRAL/ISGRI revealed more than 500 sources. Many of these sources are variable. Taking into account that nearly half of INTEGRAL/ISGRI sources are new and many of them are still unidentified, the variability properties of the sources can serve as additional parameters that may help to classify and identify the unknown sources. Aims: In order to study the variability properties of the sources detected by INTEGRAL/ISGRI we develop a method to quantify the variability of a source. We describe here our techniques and compile a catalog of the sources that fit our criteria of variability. Methods: We use the natural time binning of INTEGRAL observations called Science Window ($\approx 2000$ seconds) and test the hypothesis that the detected sources are constant using a $\chi^2$ all-sky map in three energy bands (20-40, 40-100, 100-200 keV). We calculate an intrinsic variance of the flux in individual pixels and use it to define the fractional variability of a source. The method is sensitive to the source variability on time scales of one Science Window and higher. We concentrate only on the sources which were already reported to be detected by INTEGRAL. Results: We present a catalog of 202 sources which are found to be significantly variable. For the catalog sources we give the measure of variability and fluxes with corresponding errors in 20-40, 40-100, 100-200 keV energy bands, and we present some statistics about the population of variable sources. The description of the physical properties of the variable sources will be given in a forthcoming paper.

VLA observations at 1477 MHz revealed the presence of a radio mini-halo surrounding the faint central point-like radio source in the Ophiuchus cluster of galaxies. In this work we present a study of the radio emission from this cluster of galaxies at lower radio frequencies. We observed the Ophiuchus cluster at 153, 240, and 614 MHz with the GMRT. The mini-halo is clearly detected at 153 and 240 MHz while it is not detected at 610 MHz. The most prominent feature at low frequencies is a patch of diffuse steep spectrum emission located at about 5' south-east from the cluster center. By combining these images with that at 1477 MHz, we derived the spectral index of the mini-halo. Globally, the mini-halo has a low-frequency spectral index of alpha_240^153 ~1.4 +/- 0.3 and an high-frequency spectral index of alpha_1477^240 ~ 1.60 +/- 0.05. Moreover, we measure a systematic increase of the high-frequency spectral index with radius: the azimuthal radial average of alpha_1477^240 increases from about 1.3, at the cluster center, up to about 2.0 in the mini-halo outskirts. The observed radio spectral index is in agreement with that obtained by modeling the non-thermal hard X-ray emission in this cluster of galaxies. We assume that the X-ray component arises from inverse Compton scattering between the photons of the cosmic microwave background and a population of non-thermal electrons which are isotropically distributed and whose energy spectrum is a power law with index p. We derive that the electrons energy spectrum should extend from a minimum Lorentz factor of gamma_min < 700 up to a maximum Lorentz factor of gamma_max =3.8 x 10^4 with an index p=3.8 +/- 0.4. The volume-averaged strength for a completely disordered intra-cluster magnetic field is B_V ~0.3 +/- 0.1 micro-G.

Aims. We use the IBIS/ISGRI telescope on-board INTEGRAL to measure the position of the centroid of the 20-200 keV emission from the Crab region. Methods. We find that the astrometry of the IBIS telescope is affected by the temperature of the IBIS mask during the observation. After correcting for this effect, we show that the systematic errors in the astrometry of the telescope are of the order of 0.5 arcsec. In the case of the Crab nebula and several other bright sources, the very large number of photons renders the level of statistical uncertainty in the centroid smaller or comparable to this value. Results. We find that the centroid of the Crab nebula in hard X-rays (20-40 keV) is shifted by 8.0 arcsec with respect to the Crab pulsar in the direction of the X-ray centroid of the nebula. A similar shift is also found at higher energies (40-100 and 100-200 keV). We observe a trend of decreasing shift with energy, which can be explained by an increase in the pulsed fraction. To differentiate between the contribution of the pulsar and the nebula, we divide our data into an on-pulse and off-pulse sample. Surprisingly, the nebular emission (i.e., off-pulse) is located significantly away from the X-ray centroid of the nebula. Conclusions. In all 3 energy bands (20-40, 40-100, and 100-200 keV), we find that the centroid of the nebula is significantly offset from the predicted position. We interpret this shift in terms of a cut-off in the electron spectrum in the outer regions of the nebula, which is probably the origin of the observed spectral break around 100 keV. From a simple spherically-symmetric model for the nebula, we estimate that the electrons in the external regions of the torus (d ~ 0.35 pc from the pulsar) reach a maximal energy slightly below 10^14 eV.

We investigated the non-thermal hard X-ray emission in the Ophiuchus cluster of galaxies. Our aim was to characterise the physical properties of the non-thermal component and its interaction with the cosmic microwave background. We performed spatially resolved spectroscopy and imaging using XMM-Newton data to model the thermal emission. Combining this with INTEGRAL ISGRI data, we modeled the 0.6-140 keV band total emission in the central 7 arcmin region. The models that best describe both PN and ISGRI data contain a power-law component with a photon index in a range 2.2-2.5. This component produces ~10% of the total flux in the 1-10 keV band. The pressure of the non-thermal electrons is ~1% of that of the thermal electrons. Our results support the scenario whereby a relativistic electron population, which produces the recently detected radio mini-halo in Ophiuchus, also produces the hard X-rays via inverse compton scattering of the CMB photons. The best-fit models imply a differential momentum spectrum of the relativistic electrons with a slope of 3.4-4.0 and a magnetic field strength B=0.05-0.15 microG. The lack of evidence for a recent major merger in the Ophiuchus center allows the possibility that the relativistic electrons are produced by turbulence or hadronic collisions.

We present the results of a 500 ksec observation of the Perseus cluster with INTEGRAL, with the aim of investigating the possible diffuse non-thermal component detected in a previous Chandra observation. In the 3-20 keV band with the JEM-X instrument, we detect the source with high significance and resolve it spatially. Above 20 keV with IBIS/ISGRI, we find that the source is point-like, and the cluster could be detected up to 120 keV. From the broad-band ISGRI/JEM-X spectrum, although we detect a non-thermal component, we find that the high-energy flux is variable and is consistent with the extrapolation of the 2-10 keV flux of the central AGN, NGC 1275. The extrapolation of the non-thermal component claimed from Chandra data exceeds the INTEGRAL spectrum by a factor of 3.

We present the results of deep observations of the Ophiuchus cluster of galaxies with INTEGRAL in the 3-80 keV band. We analyse 3 Ms of INTEGRAL data on the Ophiuchus cluster with the IBIS/ISGRI hard X-ray imager and the JEM-X X-ray monitor. In the X-ray band using JEM-X, we show that the source is extended, and that the morphology is compatible with the results found by previous missions. Above 20 keV, we show that the size of the source is slightly larger than the PSF of the instrument, and is consistent with the soft X-ray morphology found with JEM-X and ASCA. Thanks to the constraints on the temperature provided by JEM-X, we show that the spectrum of the cluster is not well fitted by a single-temperature thermal Bremsstrahlung model, and that another spectral component is needed to explain the high energy data. We detect the high energy tail with a higher detection significance (6.4 sigma) than the BeppoSAX claim (2 sigma). Because of the imaging capabilities of JEM-X and ISGRI, we are able to exclude the possibility that the excess emission comes from very hot regions or absorbed AGN, which proves that the excess emission is indeed of non-thermal origin. Using the available radio data together with the non-thermal hard X-ray flux, we estimate a magnetic field B ~ 0.1-0.2 mu G.

We present the results of a deep (1.1 Ms) observation of the Coma cluster of galaxies in the 18-30 keV band with the IBIS/ISGRI imager on board the INTEGRAL satellite. We show that the source extension in the North-East to South-West (SW) direction ($\sim 17'$) significantly exceeds the size of the point spread function of ISGRI, and that the centroid of the image of the source in the 18-30 keV band is displaced in the SW direction compared to the centroid in the 1-10 keV band. To test the nature of the SW extension we fit the data assuming different models of source morphology. The best fit is achieved with a diffuse source of elliptical shape, although an acceptable fit can be achieved assuming an additional point source SW of the cluster core. In the case of an elliptical source, the direction of extension of the source coincides with the direction toward the subcluster falling onto the Coma cluster. If the SW excess is due to the presence of a point source with a hard spectrum, we show that there is no obvious X-ray counterpart for this additional source, and that the closest X-ray source is the quasar EXO 1256+281, which is located $6.1'$ from the centroid of the excess. Finally, we show that the hard X-ray emission coincides with the 1.4 GHz radio emission, which suggests that the hard X-ray emission comes from the same population of electrons that is responsible for radio haloes through synchrotron emission.

We explore the morphology of hard (18-30 keV) X-ray emission from the Coma cluster of galaxies. We analyze a deep (1.1 Ms) observation of the Coma cluster with the ISGRI imager on board the \emph{INTEGRAL} satellite. We show that the source extension in the North-East to South-West (SW) direction ($\sim 17'$) significantly exceeds the size of the point spread function of ISGRI, and that the centroid of the image of the source in the 18-30 keV band is displaced in the SW direction compared to the centroid in the 1-10 keV band. To test the nature of the SW extension we fit the data assuming different models of source morphology. The best fit is achieved with a diffuse source of elliptical shape, although an acceptable fit can be achieved assuming an additional point source SW of the cluster core. In the case of an elliptical source, the direction of extension of the source coincides with the direction toward the subcluster falling onto the Coma cluster. If the SW excess is due to the presence of a point source with a hard spectrum, we show that there is no obvious X-ray counterpart for this additional source, and that the closest X-ray source is the quasar EXO 1256+281, which is located $6.1'$ from the centroid of the excess. The observed morphology of the hard X-ray emission clarifies the nature of the hard X-ray "excess" emission from the Coma cluster, which is due to the presence of an extended hard X-ray source SW of the cluster core.

We report the discovery of a high-energy transient with the IBIS/ISGRI detector on board the INTEGRAL observatory. The source, namely IGR J00291+5934, was first detected on 2nd December 2004 in the routine monitoring of the IBIS/ISGRI 20--60 keV images. The observations were conducted during Galactic Plane Scans, which are a key part of the INTEGRAL Core Programme observations. After verifying the basic source behaviour, the discovery was announced on 3rd December. The transient shows a hard Comptonised spectrum, with peak energy release at about 20 keV and a total luminosity of ~ 0.9E36 erg/s in the 5--100 keV range, assuming a distance of 3 kpc. Following the INTEGRAL announcement of the discovery of IGR J00291+5934, a number of observations were made by other instruments. We summarise the results of those observations and, together with the INTEGRAL data, identifiy IGR J00291+5934 as the 6th member of a class of accreting X-ray millisecond pulsars.