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.

Galaxy clusters provide a multitude of observational data across wavelengths and their structure and morphology are of considerable interest in cosmology as well as astrophysics. We develop a framework that allows the combination of lensing and non-lensing observations in a free-form and mesh-free approach to infer the projected mass distribution of individual galaxy clusters. This method can be used to test common assumptions on the morphology of clusters in parametric models. We make use of the lensing reconstruction code SaWLens2 and expand its capabilities by incorporating an estimate of the projected gravitational potential based on X-ray data that are deprojected using the local Richardson-Lucy method and used to infer the Newtonian potential of the cluster and we discuss how potentially arising numerical artefacts can be treated. We demonstrate the feasibility of our method on a simplified mock NFW halo and on a cluster from a realistic hydrodynamical simulation and show how the combination of X-ray and weak lensing data can affect a free-form reconstruction, improving the accuracy in the central region in some cases by a factor of two.

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.

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)

The hot electrons in the intra-cluster medium produce a spectral distorsion of the cosmic microwave background (CMB) black body emission, the thermal Sunyaev-Zel'dovich effect (tSZ). This characteristic spectral distorsion is now commonly used to detect and characterize the properties of galaxy clusters. The tSZ effect spectral distorsion does not depends on the redshift, and is only slightly affected by the galaxy cluster properties via the relativistic corrections, when the electrons reach relativistic velocities. The present work proposes a linear component separation approach to extract the tSZ effect Compton parameter and relativistic corrections for next-generation CMB experiments. We demonstrated that relativistic corrections, if neglected, would induce a significant bias on galaxy cluster Compton parameter, tSZ scaling relation slope, and tSZ angular power spectrum shape measurements. We showed that tSZ relativistic corrections mapping can be achieved at high signal-to-noise ratio with a low level of contamination up to $\ell=3000$ for next generation CMB experiments. At smaller angular scales the contamination produced by infra-red emission will be a significant source of bias. Such tSZ relativistic corrections mapping enables the study of galaxy cluster temperature profile via the tSZ effect only.

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.

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).

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.

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.

We present a self-consistent interpretation of the of very-high-energy neutrino signal from the direction of the inner Galaxy, which is a part of the astronomical neutrino signal reported by IceCube. We demonstrate that an estimate of the neutrino flux in the E>100 TeV energy range lies at the high-energy power-law extrapolation of the spectrum of diffuse gamma-ray emission from the Galactic Ridge, as observed by Fermi telescope. This suggests that IceCube neutrino and Fermi/LAT gamma-ray fluxes are both produced in interactions of cosmic rays with the interstellar medium in the Norma arm and/or in the Galactic Bar. Cosmic rays responsible for the gamma-ray and neutrino flux are characterised by hard spectrum with the slope harder than -2.4 and cut-off energy higher than 10 PeV.

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.

Context. The imaging and timing properties of {\gamma}-ray emission from electromagnetic cascades initiated by very-high-energy (VHE) {\gamma}-rays in the intergalactic medium depend on the strength B and correlation length {\lambda}B of intergalactic magnetic fields (IGMF). Aims. We study the possibility of measuring both B and {\lambda}B via observations of the cascade emission with {\gamma}-ray telescopes. Methods. For each measurement method, we find two characteristics of the cascade signal, which are sensitive to the IGMF B and {\lambda}B values in different combinations. For the case of IGMF measurement using the observation of extended emission around extragalactic VHE {\gamma}-ray sources, the two characteristics are the slope of the surface brightness profile and the overall size of the cascade source. For the case of IGMF measurement from the time delayed emission, these two characteristics are the initial slope of the cascade emission light curve and the overall duration of the cascade signal. Results. We show that measurement of the slope of the cascade induced extended emission and/or light curve can both potentially provide measure of the IGMF correlation length, provided it lies within the range 10 kpc< {\lambda}B <1 Mpc. For correlation lengths outside this range, gamma-ray observations can provide upper or lower bound on {\lambda}B. The latter of the two methods holds great promise in the near future for providing a measurement/constraint using measurements from present/next-generation {\gamma}-ray-telescopes. Conclusions. Measurement of the IGMF correlation length will provide an important constraint on its origin. In particular, it will enable to distinguish between an IGMF of galactic wind origin from an IGMF of cosmological origin.

We explore the detectability of the neutrino flux from the entire Galactic Plane or from a part of it with IceCube. We calculate the normalization and the spectral index of the neutrino power law spectrum from different regions of the Galactic plane, based on the observed spectral characteristics of the pion decay gamma-ray diffuse emission observed by the Fermi/LAT telescope in the energy band above 100 GeV. We compare the neutrino flux calculated in this way with the sensitivity of IceCube for the detection of extended sources. Assuming a binned extended source analysis method, we find that the only possible evidence for neutrino emission for sources located in the Northern hemisphere is from the Cygnus region after 20 years of exposure. For other parts of the Galactic Plane even a 20 years exposure with IceCube is not sufficient for the detection. Taking into account marginal significance of the detectable source in the Cygnus region, we find a precise position and size of the source region which optimizes the signal-to-noise ratio for neutrinos. We also calculate the low-energy threshold above which the neutrino signal could be detected with the highest signal-to-noise ratio. This calculation of precise source position, size and energy range, based on the gamma-ray data, could be used to remove the 'trial factor' in the analysis of the real neutrino data of IceCube. We notice that the diffuse neutrino emission from the inner Galactic Plane in the Southern Hemisphere is much brighter. A neutrino detector with characteristics equivalent to IceCube, but placed at the Northern Hemisphere (such as KM3NeT), would detect several isolated neutrino sources in the Galactic Plane within just 5 years exposure at 5{\sigma} level. These isolated sources of ~TeV neutrinos would unambiguously localize sources of cosmic rays which operated over the last 10 thousand years in the Galaxy.[abridged]

Context: Hadronic models, involving matter (proton or nuclei) acceleration in blazar jets, imply high energy photon and neutrino emissions due to interactions of high-energy protons with matter and/or radiation in the source environment. Aims: This paper shows that the sensitivity of the IceCube neutrino telescope in its 40-string configuration (IC-40) is already at the level of constraining the parameter space of purely hadronic scenarios of activity of blazars. Methods: Assuming that the entire source power originates from hadronic interactions, and assuming that the models describe the data, we estimate the expected neutrino flux from blazars based on the observed gamma-ray flux by Fermi, simultaneously with IC-40 observations. We consider two cases separately to keep the number of constrainable parameters at an acceptable level: proton-proton or proton-gamma interactions are dominant. Comparing the IC-40 sensitivity to the neutrino flux expected from some of the brightest blazars, we constrain model parameters characterizing the parent high-energy proton spectrum. Results: We find that when pp interactions dominate, some constraints on the primary proton spectrum can be imposed. For instance, for the tightest constrained source 3C 454.3, the very high energy part of the spectra of blazars is constrained to be harder than E^-2 with cut-off energies in the range of Ecut >10^18 eV. When interactions of high-energy protons on soft photon fields dominate, we can find similarly tight constraints on the proton spectrum parameters. [abridged]