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

Recent observations have revealed the presence of strong CIII] emission (EW$_{\rm{CIII]}}>20$ \r{A}) in $z>6$ galaxies, the origin of which remains unclear. In an effort to understand the nature of these line emitters, we have initiated a survey targeting CIII] emission in gravitationally-lensed reionization era analogs identified in HST imaging of clusters from the RELICS survey. Here we report initial results on four galaxies selected to have low stellar masses (2-8$\times$10$^7$ M$_\odot$) and J$_{125}$-band flux excesses indicative of intense [OIII]+H$\beta$ emission (EW$_{\rm{[OIII]+H\beta}}$=500-2000 \r{A}), similar to what has been observed at $z>6$. We detect CIII] emission in three of the four sources, with the CIII] EW reaching values seen in the reionization era (EW$_{\rm{CIII]}}\simeq 17-22$ \r{A}) in the two sources with the strongest optical line emission (EW$_{\rm{[OIII]+H\beta}}\simeq 2000$ \r{A}). We have obtained a Magellan/FIRE near-infrared spectrum of the strongest CIII] emitter in our sample, revealing gas that is both metal poor and highly ionized. Using photoionization models, we are able to simultaneously reproduce the intense CIII] and optical line emission for extremely young (2-3 Myr) and metal poor (0.06-0.08 Z$_\odot$) stellar populations, as would be expected after a substantial upturn in the SFR of a low mass galaxy. The sources in this survey are among the first for which CIII] has been used as the primary means of redshift confirmation. We suggest that it should be possible to extend this approach to $z>6$ with current facilities, using CIII] to measure redshifts of objects with IRAC excesses indicating EW$_{\rm{[OIII]+H\beta}}\simeq 2000$ \r{A}, providing a method of spectroscopic confirmation independent of Ly$\alpha$.

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

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

We present deep Chandra observations of A3411-12, a remarkable merging cluster that hosts the most compelling evidence for electron re-acceleration at cluster shocks to date. Using the $Y_X-M$ scaling relation, we find $r_{500} \sim 1.3$ Mpc, $M_{500} = (7.1 \pm 0.7) \times 10^{14} \ M_{\rm{\odot}}$, $kT=6.5\pm 0.1$ keV, and a gas mass of $M_{\rm g,500} = (9.7 \pm 0.1) \times 10^{13} M_\odot$. The gas mass fraction within $r_{500}$ is $f_{\rm g} = 0.14 \pm 0.01$. We compute the shock strength using density jumps to conclude that the Mach number of the merging subcluster is small ($M \leq 1.15_{-0.09}^{+0.14}$). We also present pseudo-density, projected temperature, pseudo-pressure, and pseudo-entropy maps. Based on the pseudo-entropy map we conclude that the cluster is undergoing a mild merger, consistent with the small Mach number. On the other hand, radio relics extend over Mpc scale in the A3411-12 system, which strongly suggests that a population of energetic electrons already existed over extended regions of the cluster.

The largest galaxy clusters are observed still to be forming through major cluster-cluster mergers, often showing observational signatures such as radio relics and giant radio haloes. Using LOFAR Two-meter Sky Survey data, we present new detections of both a radio halo (with a spectral index of $\alpha_{143}^{1400}=1.48^{+0.06}_{-0.23}$) and a likely radio relic in Abell 959, a massive cluster at a redshift of z=0.288. Using a sample of clusters with giant radio haloes from the literature (80 in total), we show that the radio halo in A959 lies reasonably well on the scaling relations between the thermal and non-thermal power of the system. Additionally, we find evidence that steep-spectrum haloes tend to reside in clusters with high X-ray luminosities relative to those expected from cluster LM scaling relations, indicating that such systems may preferentially lie at an earlier stage of the merger, consistent with the theory that some steep-spectrum haloes result from low-turbulence mergers. Lastly, we find that halo systems containing radio relics tend to lie at lower X-ray luminosities, relative to those expected from cluster LM scaling relations, for a given halo radio power than those without relics, suggesting that the presence of relics indicates a later stage of the merger, in line with simulations.

Measurements of stellar properties of galaxies when the universe was less than one billion years old yield some of the only observational constraints of the onset of star formation. We present here the inclusion of \textit{Spitzer}/IRAC imaging in the spectral energy distribution fitting of the seven highest-redshift galaxy candidates selected from the \emph{Hubble Space Telescope} imaging of the Reionization Lensing Cluster Survey (RELICS). We find that for 6/8 \textit{HST}-selected $z\gtrsim8$ sources, the $z\gtrsim8$ solutions are still strongly preferred over $z\sim$1-2 solutions after the inclusion of \textit{Spitzer} fluxes, and two prefer a $z\sim 7$ solution, which we defer to a later analysis. We find a wide range of intrinsic stellar masses ($5\times10^6 M_{\odot}$ -- $4\times10^9$ $M_{\odot}$), star formation rates (0.2-14 $M_{\odot}\rm yr^{-1}$), and ages (30-600 Myr) among our sample. Of particular interest is Abell1763-1434, which shows evidence of an evolved stellar population at $z\sim8$, implying its first generation of star formation occurred just $< 100$ Myr after the Big Bang. SPT0615-JD, a spatially resolved $z\sim10$ candidate, remains at its high redshift, supported by deep \textit{Spitzer}/IRAC data, and also shows some evidence for an evolved stellar population. Even with the lensed, bright apparent magnitudes of these $z \gtrsim 8$ candidates (H = 26.1-27.8 AB mag), only the \textit{James Webb Space Telescope} will be able further confirm the presence of evolved stellar populations early in the universe.

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

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

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

We present the results from new deep {\it Chandra} ($\sim410$~ks) and {\it Suzaku} ($\sim180$ ks) observations of the merging galaxy cluster ZwCl\,0008.8+5215 ($z=0.104$). Previous radio observations revealed the presence of a double radio relic located diametrically west and east of the cluster center. Using our new {\it Chandra} data, we find evidence for the presence of a shock at the location of the western relic, RW, with a Mach number $\mathcal{M}_{S_X}=1.48^{+0.50}_{-0.32}$ from the density jump. We also measure $\mathcal{M}_{T_X}=2.35^{+0.74}_{-0.55}$ and $\mathcal{M}_{T_X}=2.02^{+0.74}_{-0.47}$ from the temperature jump, with {\it Chandra} and {\it Suzaku} respectively. These values are consistent with the Mach number estimate from a previous study of the radio spectral index, under the assumption of diffusive shock acceleration ($\mathcal{M}_{\rm RW}=2.4^{+0.4}_{-0.2}$). Interestingly, the western radio relic does not entirely trace the X-ray shock. A possible explanation is that the relic traces fossil plasma from nearby radio galaxies which is re-accelerated at the shock. For the eastern relic we do not detect an X-ray surface brightness discontinuity, despite the fact that radio observations suggest a shock with $\mathcal{M}_{\rm RE}=2.2^{+0.2}_{-0.1}$. The low surface brightness and reduced integration time for this region might have prevented the detection. {\it Chandra} surface brightness profile suggests $\mathcal{M}\lesssim1.5$, while {\it Suzaku} temperature measurements found $\mathcal{M}_{T_X}=1.54^{+0.65}_{-0.47}$. Finally, we also detect a merger induced cold front on the western side of the cluster, behind the shock that traces the western relic.

We examine the possible acceleration mechanisms of the relativistic particles responsible for the extended radio emission in Abell 520. We used new LOFAR 145 MHz, archival GMRT 323 MHz and VLA 1.5 GHz data to study the morphological and spectral properties of extended cluster emission. The observational properties are discussed in the framework of particle acceleration models associated with cluster merger turbulence and shocks. In Abell 520, we confirm the presence of extended synchrotron radio emission that has been classified as a radio halo. The comparison between the radio and X-ray brightness suggests that the halo might originate in a cocoon rather than from the central X-ray bright regions of the cluster. The halo spectrum is roughly uniform on the scale of 66 kpc. There is a hint of spectral steepening from the SW edge towards the cluster centre. Assuming DSA, the radio data are suggestive of a shock of $\mathcal{M}_{SW}=2.6_{-0.2}^{+0.3}$ that is consistent with the X-ray derived estimates. This is in line with the scenario in which relativistic electrons in the SW radio edge gain their energies at the shock front via acceleration of either thermal or fossil electrons. We do not detect extended radio emission ahead of the SW shock that is predicted if the emission is the result of adiabatic compression. An X-ray surface brightness discontinuity is detected towards the NE region that may be a counter shock of $\mathcal{M}_{NE}^{X}=1.52\pm0.05$. This is lower than the value predicted from the radio emission ($\mathcal{M}_{NE}=2.1\pm0.2$). Our observations indicate that the SW radio emission in Abell 520 is likely effected by the prominent X-ray detected shock in which radio emitting particles are (re-)accelerated through the Fermi-I mechanism. The NE X-ray discontinuity that is approximately collocated with an edge in the radio emission hints at the presence of a counter shock.

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

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

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

Golovich et al. 2017b presents an optical imaging and spectroscopic survey of 29 radio relic merging galaxy clusters. In this paper, we study this survey to identify substructure and quantify the dynamics of the mergers. Using a combined photometric and spectroscopic approach, we identify the minimum number of substructures in each system to describe the galaxy populations and estimate the line of sight velocity difference between likely merging subclusters. We find that the line-of-sight velocity components of the mergers are typically small compared with the maximum three dimensional relative velocity (usually $<1000$ km s$^{-1}$ and often consistent with zero). This suggests that the merger axes of these systems are generally in or near the plane of the sky matching findings in magneto-hydrodynamical simulations. In 28 of the 29 systems we identify substructures in the galaxy population aligned with the radio relic(s) and presumed associated merger induced shock. From this ensemble, we identify eight systems to include in a `gold' sample that is prime for further observation, modeling, and simulation study. Additional papers will present weak lensing mass maps and dynamical modeling for each merging system, ultimately leading to new insight into a wide range of astrophysical phenomena at some of the largest scales in the universe.

We present a study of X-ray source populations in M87, the cD galaxy of the Virgo cluster, using 12 archival {\it Chandra} observations with a total exposure of $\sim$680 ks spanning about a decade. A total of 346 point-like sources are detected down to a limiting 0.5--8 keV luminosity of $4\times10^{37}{\rm~erg~s^{-1}}$ and out to a galactocentric radius of $\sim$40 kpc. We cross-correlate the X-ray sources with published catalogs of globular clusters (GCs), derived from the ACS Virgo Cluster Survey and the Next Generation Virgo Cluster Survey. This results in 122 matches, making it one of the largest samples of GC-hosting X-ray sources in an external galaxy. These sources, most likely low-mass X-ray binaries (LMXBs), correspond to $\sim$5\% of all known GCs within the {\it Chandra} field-of-view. Conversely, $\sim$50\% of the detected X-ray sources are found in a GC. Moreover, red (metal-rich) GCs are $\sim$2.2 times more likely to host an X-ray source than blue (metal-poor) GCs. We also examine 76 currently known ultra-compact dwarf galaxies around M87, but find no significant X-ray counterparts. After statistically accounting for the cosmic X-ray background, we identify $\sim$110 field-LMXBs. The GC-LMXBs and field-LMXBs differ in their luminosity function and radial distribution, which indicates that the latter cannot be primarily originated from GCs. Using another set of deep {\it Chandra} observations toward $\sim$100 kpc northwest of the M87 center, we statistically constrain the abundance of field-LMXBs in the stellar halo, which is consistent with that found in the central region. We also identify 40 variable X-ray sources, among which one source is likely a black hole binary residing in a GC.

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

We present LOFAR $120-168$ MHz images of the merging galaxy cluster Abell 1240 that hosts double radio relics. In combination with the GMRT $595-629$ MHz and VLA $2-4$ GHz data, we characterised the spectral and polarimetric properties of the radio emission. The spectral indices for the relics steepen from their outer edges towards the cluster centre and the electric field vectors are approximately perpendicular to the major axes of the relics. The results are consistent with the picture that these relics trace large-scale shocks propagating outwards during the merger. Assuming diffusive shock acceleration (DSA), we obtain shock Mach numbers of $\mathcal{M}=2.4$ and $2.3$ for the northern and southern shocks, respectively. For $\mathcal{M}\lesssim3$ shocks, a pre-existing population of mildly relativistic electrons is required to explain the brightness of the relics due to the high ($>10$ per cent) particle acceleration efficiency required. However, for $\mathcal{M}\gtrsim4$ shocks the required efficiency is $\gtrsim1\%$ and $\gtrsim0.5\%$, respectively, which is low enough for shock acceleration directly from the thermal pool. We used the fractional polarization to constrain the viewing angle to $\geqslant(53\pm3)^\circ$ and $\geqslant(39\pm5)^\circ$ for the northern and southern shocks, respectively. We found no evidence for diffuse emission in the cluster central region. If the halo spans the entire region between the relics ($\sim1.8\,\text{Mpc}$) our upper limit on the power is $P_\text{1.4 GHz}=(1.4\pm0.6)\times10^{23}\,\text{W}\text{Hz}^{-1}$ which is approximately equal to the anticipated flux from a cluster of this mass. However, if the halo is smaller than this, our constraints on the power imply that the halo is underluminous.

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

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

The most distant galaxies known are at z~10-11, observed 400-500 Myr after the Big Bang. The few z~10-11 candidates discovered to date have been exceptionally small- barely resolved, if at all, by the Hubble Space Telescope. Here we present the discovery of SPT0615-JD, a fortuitous z~10 (z_phot=9.9+/-0.6) galaxy candidate stretched into an arc over ~2.5" by the effects of strong gravitational lensing. Discovered in the Reionization Lensing Cluster Survey (RELICS) Hubble Treasury program and companion S-RELICS Spitzer program, this candidate has a lensed H-band magnitude of 25.7+/-0.1 AB mag. With a magnification of \mu~4-7 estimated from our lens models, the de-lensed intrinsic magnitude is 27.6+/-0.3 AB mag, and the half-light radius is r_e<0.8 kpc, both consistent with other z>9 candidates. The inferred stellar mass (log [M* /M_Sun]=9.7^{+0.7}_{-0.5}) and star formation rate (\log [SFR/M_Sun yr^{-1}]=1.3^{+0.2}_{-0.3}) indicate that this candidate is a typical star-forming galaxy on the z>6 SFR-M* relation. We note that three independent lens models predict two counterimages, at least one of which should be of a similar magnitude to the arc, but these counterimages are not yet detected. Counterimages would not be expected if the arc were at lower redshift. However, the only spectral energy distributions capable of fitting the Hubble and Spitzer photometry well at lower redshifts require unphysical combinations of z~2 galaxy properties. The unprecedented lensed size of this z~10 candidate offers the potential for the James Webb Space Telescope to study the geometric and kinematic properties of a galaxy observed 500 Myr after the Big Bang.

We present Chandra observations of hot gas structures, characteristic of gas stripping during infall, in the Virgo cluster elliptical galaxy M60 (NGC4649) located $1$ Mpc east of M87. $0.5-2$ keV Chandra X-ray images show a sharp leading edge in the surface brightness $12.4 \pm 0.1$ kpc north and west of the galaxy center in the direction of M87 characteristic of a merger cold front due to M60's motion through the Virgo ICM. We measured a temperature of $1.00 \pm 0.02$ keV for abundance $0.5 Z_\odot$ inside the edge and $1.37^{+0.35}_{-0.19}$ keV for abundance $0.1 Z_\odot$ in the Virgo ICM free stream region. We find that the observed jump in surface brightness yields a density ratio of $6.44^{+1.04}_{-0.67}$ between gas inside the edge and in the cluster free stream region. If the edge is a cold front due solely to the infall of M60 in the direction of M87, we find a pressure ratio of $4.7^{+1.7}_{-1.4}$ and Mach number $1.7 \pm 0.3$. For 1.37 keV Virgo gas we find a total infall velocity for M60 of $1030 \pm 180$ kms$^{-1}$. We calculate the motion in the plane of the sky to be $1012^{+183}_{-192}$ km$^{-1}$ implying an inclination angle $\xi = 11 \pm 3$ degrees. Surface brightness profiles show the presence of a faint diffuse gaseous tail. We identify filamentary, gaseous wing structures caused by the galaxy's motion through the ICM. The structure and dimensions of these wings are consistent with simulations of Kelvin-Helmholtz instabilities as expected if the gas stripping is close to inviscid.

We report L-band VLA observations of 1RXS J0603.3+4214, a cluster that hosts a bright radio relic, known as the Toothbrush, and an elongated giant radio halo. These new observations allow us to study the surface brightness distribution down to one arcsec resolution with very high sensitivity. Our images provide an unprecedented detailed view of the Toothbrush, revealing enigmatic filamentary structures. To study the spectral index distribution, we complement our analysis with published LOFAR and GMRT observations. The bright `brush' of the Toothbrush shows a prominent narrow ridge to its north with a sharp outer edge. The spectral index at the ridge is in the range $-0.70\leq\alpha\leq-0.80$. We suggest that the ridge is caused by projection along the line of sight. With a simple toy model for the smallest region of the ridge, we conclude that the magnetic field is below $5\,\rm\mu G$ and varies significantly across the shock front. Our model indicates that the actual Mach number is higher than that obtained from the injection index and agrees well with the one derived from the overall spectrum, namely ${\cal M}=3.78^{+0.3}_{-0.2}$. The radio halo shows an average spectral index of $\alpha=-1.16\pm0.05$ and a slight gradient from north to south. The southernmost part of the halo is steeper and possibly related to a shock front. Excluding the southernmost part, the halo morphology agrees very well with the X-ray morphology. A power-law correlation is found between the radio and X-ray surface brightness

The evolution of the X-ray emitting gas mass fraction ($f_{gas}$) in massive galaxy clusters can be used as an independent cosmological tool to probe the expansion history of the Universe. Its use, however, depends upon a crucial quantity, i.e., the depletion factor $\gamma$, which corresponds to the ratio by which $f_{gas}$ is depleted with respect to the universal baryonic mean. This quantity is not directly observed and hydrodynamical simulations performed in a specific cosmological model (e.g., a flat $\Lambda$CDM cosmology) have been used to calibrate it. In this work, we obtain for the first time self-consistent observational constraints on the gas depletion factor combining 40 X-ray emitting gas mass fraction measurements and luminosity distance measurements from type Ia supernovae. Using Gaussian Processes to reconstruct a possible redshift evolution of $\gamma$, we find no evidence for such evolution, which confirms the current results from hydrodynamical simulations. Moreover, our constraints on $\gamma$ can be seen as a data prior for cosmological analyses on different cosmological models. The current measurements are systematic limited, so future improvements will depend heavily on a better mass calibration of galaxy clusters and their measured density profiles.

Multi-band photometric and multi-object spectroscopic surveys of merging galaxy clusters allow for the characterization of the distributions of constituent dark matter and galaxy populations, constraints on the dynamics of the merging subclusters, and an understanding of galaxy evolution of member galaxies. We present deep photometric observations from Subaru/SuprimeCam and a catalog of $\sim$5400 spectroscopic cluster members from Keck/DEIMOS across 29 merging galaxy clusters ranging in redshift from $z=0.07$ to $0.55$. The ensemble is compiled based on the presence of radio relics, which highlight cluster scale collisionless shocks in the intra-cluster medium. Together with the spectroscopic and photometric information, the velocities, timescales, and geometries of the respective merging events may be tightly constrained. In this preliminary analysis, the velocity distributions of 28 of the 29 clusters are shown to be well fit by single Gaussians. This indicates that radio relic mergers largely occur transverse to the line of sight and/or near apocenter. In this paper, we present our optical and spectroscopic surveys, preliminary results, and a discussion of the value of radio relic mergers for developing accurate dynamical models of each system.

LOFAR observations at 144 MHz have revealed large-scale radio sources in the unrelaxed galaxy cluster Abell 1132. The cluster hosts diffuse radio emission on scales of $\sim$650 kpc near the cluster center and a head-tail (HT) radio galaxy, extending up to 1 Mpc, South of the cluster center. The central diffuse radio emission is not seen in NVSS, FIRST, WENSS, nor in C & D array VLA observations at 1.4 GHz, but is detected in our follow-up GMRT observations at 325 MHz. Using LOFAR and GMRT data, we determine the spectral index of the central diffuse emission to be $\alpha=-1.75\pm0.19$ ($S\propto\nu^{\alpha}$). We classify this emission as an ultra-steep spectrum radio halo and discuss the possible implications for the physical origin of radio halos. The HT radio galaxy shows narrow, collimated emission extending up to 1 Mpc and another 300 kpc of more diffuse, disturbed emission, giving a full projected linear size of 1.3 Mpc - classifying it as a giant radio galaxy (GRG) and making it the longest HT found to date. The head of the GRG coincides with an elliptical galaxy (SDSS J105851.01$+$564308.5) belonging to Abell 1132. In our LOFAR image, there appears to be a connection between the radio halo and the GRG. The turbulence that may have produced the halo may have also affected the tail of the GRG. In turn, the GRG may have provided seed electrons for the radio halo.

We study the prototypical Seyfert 2 galaxy, Markarian 3, based on imaging and high-resolution spectroscopy observations taken by the Chandra X-ray Observatory. We construct a deconvolved X-ray image, which reveals the S-shaped morphology of the hot gas in the narrow line region (NLR). While this morphology is similar to the radio and [O III] emission, the distribution of the X-ray gas is broader than that obtained at these other wavelengths. By mapping the density and temperature distribution of the hot gas in the NLR, we demonstrate the presence of shocks towards the west ($M=2.5^{+1.0}_{-0.6}$) and east ($M=1.5^{+1.0}_{-0.5}$). Moreover, we compute the flux ratios between the [O III] and $0.5-2$ keV band X-ray luminosity and show that it is non-uniform in the NLR with the western side of the NLR being more highly ionized. In addition, based on the Chandra grating data we investigate the line ratios of the Si XIII triplet, which are not consistent with pure photoionization. Based on these results, we suggest that in the NLR of Mrk 3 both photoionization and collisional ionization act as excitation mechanisms. We conclude that the canonical picture, in which photoionization is solely responsible for exciting the interstellar medium in the NLR of Seyfert galaxies, may be overly simplistic. Given that weak and small-scale radio jets are commonly detected in Seyfert galaxies, it is possible that shock heating plays a non-negligible role in the NLR of these galaxies.

We examine the reconstruction of galaxy cluster radial density profiles obtained from Chandra and XMM X-ray observations, using high quality data for a sample of twelve objects covering a range of morphologies and redshifts. By comparing the results obtained from the two observatories and by varying key aspects of the analysis procedure, we examine the impact of instrumental effects and of differences in the methodology used in the recovery of the density profiles. We find that the final density profile shape is particularly robust. We adapt the photon weighting vignetting correction method developed for XMM for use with Chandra data, and confirm that the resulting Chandra profiles are consistent with those corrected a posteriori for vignetting effects. Profiles obtained from direct deprojection and those derived using parametric models are consistent at the 1% level. At radii larger than $\sim$6", the agreement between Chandra and XMM is better than 1%, confirming an excellent understanding of the XMM PSF. We find no significant energy dependence. The impact of the well-known offset between Chandra and XMM gas temperature determinations on the density profiles is found to be negligible. However, we find an overall normalisation offset in density profiles of the order of $\sim$2.5%, which is linked to absolute flux cross-calibration issues. As a final result, the weighted ratios of Chandra to XMM gas masses computed at R2500 and R500 are r=1.03$\pm$0.01 and r=1.03$\pm$0.03, respectively. Our study confirms that the radial density profiles are robustly recovered, and that any differences between Chandra and XMM can be constrained to the $\sim$ 2.5% level, regardless of the exact data analysis details. These encouraging results open the way for the true combination of X-ray observations of galaxy clusters, fully leveraging the high resolution of Chandra and the high throughput of XMM.

Cluster mergers leave distinct signatures in the ICM in the form of shocks and diffuse cluster radio sources that provide evidence for the acceleration of relativistic particles. However, the physics of particle acceleration in the ICM is still not fully understood. Here we present new 1-4 GHz Jansky Very Large Array (VLA) and archival Chandra observations of the HST Frontier Fields Cluster Abell 2744. In our new VLA images, we detect the previously known $\sim2.1$ Mpc radio halo and $\sim1.5$ Mpc radio relic. We carry out a radio spectral analysis from which we determine the relic's injection spectral index to be $\alpha_{\rm{inj}} = -1.12 \pm 0.19$. This corresponds to a shock Mach number of $\mathcal{M}$ = 2.05$^{+0.31}_{-0.19}$ under the assumption of diffusive shock acceleration. We also find evidence for spectral steepening in the post-shock region. We do not find evidence for a significant correlation between the radio halo's spectral index and ICM temperature. In addition, we observe three new polarized diffuse sources and determine two of these to be newly discovered giant radio relics. These two relics are located in the southeastern and northwestern outskirts of the cluster. The corresponding integrated spectral indices measure $-1.81 \pm 0.26$ and $-0.63 \pm 0.21$ for the SE and NW relics, respectively. From an X-ray surface brightness profile we also detect a possible density jump of $R=1.39^{+0.34}_{-0.22}$ co-located with the newly discovered SE relic. This density jump would correspond to a shock front Mach number of $\mathcal{M}=1.26^{+0.25}_{-0.15}$.

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

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

The rich galaxy cluster Abell 2204 exhibits edges in its X-ray surface brightness at $\sim 65$ and $35 {\rm~ kpc}$ west and east of its center, respectively. The presence of these edges, which were interpreted as sloshing cold fronts, implies that the intracluster medium was recently disturbed. We analyze the properties of the intracluster medium using multiple Chandra observations of Abell 2204. We find a density ratio $n_{\rm in}/n_{\rm out} = 2.05\pm0.05$ and a temperature ratio $T_{\rm out}/T_{\rm in} = 1.91\pm0.27$ (projected, or $1.87\pm0.56$ deprojected) across the western edge, and correspondingly $n_{\rm in}/n_{\rm out} = 1.96\pm0.05$ and $T_{\rm out}/T_{\rm in} =1.45\pm0.15$ (projected, or $1.25\pm0.26$ deprojected) across the eastern edge. These values are typical of cold fronts in galaxy clusters. This, together with the spiral pattern observed in the cluster core, supports the sloshing scenario for Abell 2204. No Kelvin-Helmholtz eddies are observed along the cold front surfaces, indicating that they are effectively suppressed by some physical mechanism. We argue that the suppression is likely facilitated by the magnetic fields amplified in the sloshing motion, and deduce from the measured gas properties that the magnetic field strength should be greater than $24\pm6$ $\mu$G and $32\pm8$ $\mu$G along the west and east cold fronts, respectively.

To investigate the relationship between thermal and non-thermal components in merger galaxy clusters, we present deep JVLA and Chandra observations of the HST Frontier Fields cluster MACS J0717.5+3745. The Chandra image shows a complex merger event, with at least four components belonging to different merging subclusters. NW of the cluster, $\sim 0.7$ Mpc from the center, there is a ram-pressure-stripped core that appears to have traversed the densest parts of the cluster after entering the ICM from the direction of a galaxy filament to the SE. We detect a density discontinuity NNE of this core which we speculate is associated with a cold front. Our radio images reveal new details for the complex radio relic and radio halo in this cluster. In addition, we discover several new filamentary radio sources with sizes of 100-300 kpc. A few of these seem to be connected to the main radio relic, while others are either embedded within the radio halo or projected onto it. A narrow-angled-tailed (NAT) radio galaxy, a cluster member, is located at the center of the radio relic. The steep spectrum tails of this AGN leads into the large radio relic where the radio spectrum flattens again. This morphological connection between the NAT radio galaxy and relic provides evidence for re-acceleration (revival) of fossil electrons. The presence of hot $\gtrsim 20$ keV ICM gas detected by Chandra near the relic location provides additional support for this re-acceleration scenario.

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

On the largest scales, the Universe consists of voids and filaments making up the cosmic web. Galaxy clusters are located at the knots in this web, at the intersection of filaments. Clusters grow through accretion from these large-scale filaments and by mergers with other clusters and groups. In a growing number of galaxy clusters, elongated Mpc-size radio sources have been found, so-called radio relics. These relics are thought to trace relativistic electrons in the intracluster plasma accelerated by low-Mach number collisionless shocks generated by cluster-cluster merger events. A long-standing problem is how low-Mach number shocks can accelerate electrons so efficiently to explain the observed radio relics. Here we report on the discovery of a direct connection between a radio relic and a radio galaxy in the merging galaxy cluster Abell 3411-3412. This discovery indicates that fossil relativistic electrons from active galactic nuclei are re-accelerated at cluster shocks. It also implies that radio galaxies play an important role in governing the non-thermal component of the intracluster medium in merging clusters.

We present results of deep 153 ks Chandra observations of the hot, 11 keV, galaxy cluster associated with the radio galaxy 3C 438. By mapping the morphology of the hot gas and analyzing its surface brightness and temperature distributions, we demonstrate the presence of a merger bow shock. We identify the presence of two jumps in surface brightness and in density located at $\sim$400 kpc and $\sim$800 kpc from the cluster's core. At the position of the inner jump, we detect a factor of $2.3\pm 0.2$ density jump, while at the location of the outer jump, we detect a density drop of a factor of $3.5 \pm 0.7$. Combining this with the temperature distribution within the cluster, we establish that the pressure of the hot gas is continuous at the 400 kpc jump, while there is a factor of $6.2 \pm 2.8$ pressure discontinuity at 800 kpc jump. From the magnitude of the outer pressure discontinuity, using the Rankine-Hugoniot jump conditions, we determine that the sub-cluster is moving at $M = 2.3\pm 0.5$, or approximately $2600\pm 565$ km/s through the surrounding intracluster medium, creating the conditions for a bow shock. Based on these findings, we conclude that the pressure discontinuity is likely the result of an ongoing major merger between two massive clusters. Since few observations of bow shocks in clusters have been made, this detection can contribute to the study of the dynamics of cluster mergers, which offers insight on how the most massive clusters may have formed.

We report on high-resolution JVLA and Chandra observations of the HST Frontier Cluster MACS J0717.5+3745. MACS J0717.5+3745 offers the largest contiguous magnified area of any known cluster, making it a promising target to search for lensed radio and X-ray sources. With the high-resolution 1.0-6.5 GHz JVLA imaging in A and B configuration, we detect a total of 51 compact radio sources within the area covered by the HST imaging. Within this sample we find 7 lensed sources with amplification factors larger than $2$. None of these sources are identified as multiply-lensed. Based on the radio luminosities, the majority of these sources are likely star forming galaxies with star formation rates of 10-50 M$_\odot$ yr$^{-1}$ located at $1 \lesssim z \lesssim 2$. Two of the lensed radio sources are also detected in the Chandra image of the cluster. These two sources are likely AGN, given their $2-10$ keV X-ray luminosities of $\sim 10^{43-44}$ erg s$^{-1}$. From the derived radio luminosity function, we find evidence for an increase in the number density of radio sources at $0.6<z<2.0$, compared to a $z < 0.3$ sample. Our observations indicate that deep radio imaging of lensing clusters can be used to study star forming galaxies, with star formation rates as low as $\sim10$ M$_{\odot}$ yr$^{-1}$, at the peak of cosmic star formation history.

We analyzed deep $Chandra$ ACIS-I exposures of the cluster-scale X-ray halo surrounding the radio source 4C+37.11. This remarkable system hosts the closest resolved pair of super-massive black hole and an exceptionally luminous elliptical galaxy, the likely product of a series of past mergers. We characterize the halo with $r_{500} = 0.95$ Mpc, $M_{500} = (2.5 \pm 0.2) \times 10^{14} \ M_{\rm{\odot}}$, $ kT = 4.6\pm 0.2$ keV, and a gas mass of $M_{\rm g,500} = (2.2 \pm 0.1) \times 10^{13} M_\odot$. The gas mass fraction within $r_{500}$ is $f_{\rm g} = 0.09 \pm 0.01$. The entropy profile shows large non-gravitational heating in the central regions. We see several surface brightness jumps, associated with substantial temperature and density changes, but approximate pressure equilibrium, implying that these are sloshing structures driven by a recent merger. A residual intensity image shows core spiral structure closely matching that seen for the Perseus cluster, although at $z=0.055$ the spiral pattern is less distinct. We infer the most recent merger occurred $1-2$ Gyr ago and that the event that brought the two observed super-massive black holes to the system core is even older. Under that interpretation, this black hole binary pair has, unusually, remained at pc-scale separation for more than 2 Gyr.

We present results from recent Suzaku and Chandra X-ray, and MMT optical observations of the strongly merging "double cluster" A1750 out to its virial radius, both along and perpendicular to a putative large-scale structure filament. Some previous studies of individual clusters have found evidence for ICM entropy profiles that flatten at large cluster radii, as compared with the self-similar prediction based on purely gravitational models of hierarchical cluster formation, and gas fractions that rise above the mean cosmic value. Weakening accretion shocks and the presence of unresolved cool gas clumps, both of which are expected to correlate with large scale structure filaments, have been invoked to explain these results. In the outskirts of A1750, we find entropy profiles that are consistent with self-similar expectations, and gas fractions that are consistent with the mean cosmic value, both along and perpendicular to the putative large scale filament. Thus, we find no evidence for gas clumping in the outskirts of A1750, in either direction. This may indicate that gas clumping is less common in lower temperature (kT~4keV), less massive systems, consistent with some (but not all) previous studies of low mass clusters and groups. Cluster mass may therefore play a more important role in gas clumping than dynamical state. Finally, we find evidence for diffuse, cool (<1 keV) gas at large cluster radii (R200) along the filament, which is consistent with the expected properties of the denser, hotter phase of the WHIM.

The HST Frontier Fields cluster MACS J1149.6+2223 is one of the most complex merging clusters, believed to consist of four dark matter halos. We present results from deep (365 ks) Chandra observations of the cluster, which reveal the most distant cold front (z=0.544) discovered to date. In the cluster outskirts, we also detect hints of a surface brightness edge that could be the bow shock preceding the cold front. The substructure analysis of the cluster identified several components with large relative radial velocities, thus indicating that at least some collisions occur almost along the line of sight. The inclination of the mergers with respect to the plane of the sky poses significant observational challenges at X-ray wavelengths. MACS J1149.6+2223 possibly hosts a steep-spectrum radio halo. If the steepness of the radio halo is confirmed, then the radio spectrum, combined with the relatively regular ICM morphology, could indicate that MACS J1149.6+2223 is an old merging cluster.

We present deep LOFAR observations between 120-181 MHz of the "Toothbrush" (RX J0603.3+4214), a cluster that contains one of the brightest radio relic sources known. Our LOFAR observations exploit a new and novel calibration scheme to probe 10 times deeper than any previous study in this relatively unexplored part of the spectrum. The LOFAR observations, when combined with VLA, GMRT, and Chandra X-ray data, provide new information about the nature of cluster merger shocks and their role in re-accelerating relativistic particles. We derive a spectral index of $\alpha = -0.8 \pm 0.1$ at the northern edge of the main radio relic, steepening towards the south to $\alpha \approx - 2$. The spectral index of the radio halo is remarkably uniform ($\alpha = -1.16$, with an intrinsic scatter of $\leq 0.04$). The observed radio relic spectral index gives a Mach number of $\mathcal{M} = 2.8^{+0.5}_{-0.3}$, assuming diffusive shock acceleration (DSA). However, the gas density jump at the northern edge of the large radio relic implies a much weaker shock ($\mathcal{M} \approx 1.2$, with an upper limit of $\mathcal{M} \approx 1.5$). The discrepancy between the Mach numbers calculated from the radio and X-rays can be explained if either (i) the relic traces a complex shock surface along the line of sight, or (ii) if the radio relic emission is produced by a re-accelerated population of fossil particles from a radio galaxy. Our results highlight the need for additional theoretical work and numerical simulations of particle acceleration and re-acceleration at cluster merger shocks.

LOFAR, the Low-Frequency Array, is a powerful new radio telescope operating between 10 and 240 MHz. LOFAR allows detailed sensitive high-resolution studies of the low-frequency radio sky. At the same time LOFAR also provides excellent short baseline coverage to map diffuse extended emission. However, producing high-quality deep images is challenging due to the presence of direction dependent calibration errors, caused by imperfect knowledge of the station beam shapes and the ionosphere. Furthermore, the large data volume and presence of station clock errors present additional difficulties. In this paper we present a new calibration scheme, which we name facet calibration, to obtain deep high-resolution LOFAR High Band Antenna images using the Dutch part of the array. This scheme solves and corrects the direction dependent errors in a number of facets that cover the observed field of view. Facet calibration provides close to thermal noise limited images for a typical 8 hr observing run at $\sim$ 5arcsec resolution, meeting the specifications of the LOFAR Tier-1 northern survey.

Merging galaxy clusters leave long-lasting signatures on the baryonic and non-baryonic cluster constituents, including shock fronts, cold fronts, X-ray substructure, radio halos, and offsets between the dark matter and the gas components. Using observations from Chandra, the Jansky Very Large Array, the Giant Metrewave Radio Telescope, and the Hubble Space Telescope, we present a multiwavelength analysis of the merging Frontier Fields cluster MACS J0416.1-2403 (z=0.396), which consists of a NE and a SW subclusters whose cores are separated on the sky by ~250 kpc. We find that the NE subcluster has a compact core and hosts an X-ray cavity, yet it is not a cool core. Approximately 450 kpc south-south west of the SW subcluster, we detect a density discontinuity that corresponds to a compression factor of ~1.5. The discontinuity was most likely caused by the interaction of the SW subcluster with a less massive structure detected in the lensing maps SW of the subcluster's center. For both the NE and the SW subclusters, the dark matter and the gas components are well-aligned, suggesting that MACS J0416.1-2403 is a pre-merging system. The cluster also hosts a radio halo, which is unusual for a pre-merging system. The halo has a 1.4 GHz power of (1.06 +/- 0.09) x 10^{24} W Hz^{-1}, which is somewhat lower than expected based on the X-ray luminosity of the cluster. We suggest that we are either witnessing the birth of a radio halo, or have discovered a rare ultra-steep spectrum halo.

We present Chandra and XMM-Newton observations of PLCK G036.7+14.9 from the Chandra-Planck Legacy Program. The high resolution X-ray observations reveal two close subclusters, G036N and G036S, which were not resolved by previous ROSAT, optical, or recent Planck observations. We perform detailed imaging and spectral analyses and use a simplified model to study the kinematics of this system. The basic picture is that PLCK G036.7+14.9 is undergoing a major merger (mass ratio close to unity) between the two massive subclusters, with the merger largely along the line-of-sight and probably at an early stage. G036N hosts a small, moderate cool-core, while G036S has at most a very weak cool-core in the central 40 kpc region. The difference in core cooling times is unlikely to be caused by the ongoing merger disrupting a pre-existing cool-core in G036S. G036N also hosts an unresolved radio source in the center, which may be heating the gas if the radio source is extended. The Planck derived mass is higher than the X-ray measured mass of either subcluster, but is lower than the X-ray measured mass of the whole cluster, due to the fact that Planck does not resolve PLCK G036.7+14.9 into subclusters and interprets it as a single cluster. This mass discrepancy could induce significant bias to the mass function if such previously unresolved systems are common in the Planck cluster sample. High resolution X-ray observations are necessary to identify the fraction of such systems and correct such a bias for the purpose of precision cosmological studies.

We present results from Chandra, XMM-Newton, and ROSAT observations of the Planck SZ-detected cluster A3716 (PLCKG345.40-39.34 - G345). We show that G345 is, in fact, two subclusters separated on the sky by 400 kpc. We measure the subclusters' gas temperatures (~ 2-3 keV), total (~ 1-2 x 10^14 solar masses) and gas (~ 1-2 x 10^13 solar masses) masses, gas mass fraction within r500, entropy profiles, and X-ray luminosities (~ 10^43 erg/s). Using the gas density and temperature profiles for both subclusters, we show that there is good (0.8 sigma) agreement between the expected Sunyaev-Zel'dovich signal predicted from the X-ray data and that measured from the Planck mission, and better agreement within 0.6 sigma when we re-computed the Planck value assuming a two component cluster model, with relative amplitudes fixed based on the X-ray data. Dynamical analysis shows that the two galaxy subclusters are very likely (> 97% probability) gravitationally bound, and in the most likely scenario, the subclusters will undergo core passage in 500 +- 200 Myr. The northern subcluster is centrally peaked and has a low entropy core, while the southern subcluster has a high central entropy. The high central entropy in the southern subcluster can be explained either by the mergers of several groups, as suggested by the presence of five giant ellipticals or by AGN energy injection, as suggested by the presence of a strong radio source in one of its massive elliptical galaxies, or by a combination of both processes.

We present the X-ray and optical properties of the galaxy groups selected in the Chandra X-Bo\"otes survey. Our final sample comprises 32 systems at \textbf{$z<1.75$}, with 14 below $z = 0.35$. For these 14 systems we estimate velocity dispersions ($\sigma_{gr}$) and perform a virial analysis to obtain the radii ($R_{200}$ and $R_{500}$) and total masses ($M_{200}$ and $M_{500}$) for groups with at least five galaxy members. We use the Chandra X-ray observations to derive the X-ray luminosity ($L_X$). We examine the performance of the group properties $\sigma_{gr}$, $L_{opt}$ and $L_X$, as proxies for the group mass. Understanding how well these observables measure the total mass is important to estimate how precisely the cluster/group mass function is determined. Exploring the scaling relations built with the X-Bo\"otes sample and comparing these with samples from the literature, we find a break in the $L_X$-$M_{500}$ relation at approximately $M_{500} = 5\times10^{13}$ M$_\odot$ (for $M_{500} > 5\times10^{13}$ M$_\odot$, $M_{500} \propto L_X^{0.61\pm0.02}$, while for $M_{500} \leq 5\times10^{13}$ M$_\odot$, $M_{500} \propto L_X^{0.44\pm0.05}$). Thus, the mass-luminosity relation for galaxy groups cannot be described by the same power law as galaxy clusters. A possible explanation for this break is the dynamical friction, tidal interactions and projection effects which reduce the velocity dispersion values of the galaxy groups. By extending the cluster luminosity function to the group regime, we predict the number of groups that new X-ray surveys, particularly eROSITA, will detect. Based on our cluster/group luminosity function estimates, eROSITA will identify $\sim$1800 groups ($L_X = 10^{41}-10^{43}$ ergs s$^{-1}$) within a distance of 200 Mpc. Since groups lie in large scale filaments, this group sample will map the large scale structure of the local universe.

We present results from Chandra and XMM-Newton observations of Abell 98 (A98), a galaxy cluster with three major components: a relatively bright subcluster to the north (A98N), a disturbed subcluster to the south (A98S), and a fainter subcluster to the far south (A98SS). We find evidence for surface brightness and temperature asymmetries in A98N consistent with a shock-heated region to the south, which could be created by an early stage merger between A98N and A98S. Deeper observations are required to confirm this result. We also find that A98S has an asymmetric core temperature structure, likely due to a separate ongoing merger. Evidence for this is also seen in optical data. A98S hosts a wide-angle tail (WAT) radio source powered by a central active galactic nucleus (AGN). We find evidence for a cavity in the intracluster medium (ICM) that has been evacuated by one of the radio lobes, suggesting that AGN feedback is operating in this system. Examples of cavities in non-cool core clusters are relatively rare. The three subclusters lie along a line in projection, suggesting the presence of a large-scale filament. We observe emission along the filament between A98N and A98S, and a surface brightness profile shows emission consistent with the overlap of the subcluster extended gas haloes. We find the temperature of this region is consistent with the temperature of the gas at similar radii outside this bridge region. Lastly, we examine the cluster dynamics using optical data. We conclude A98N and A98S are likely bound to one another, with a 67% probability, while A98S and A98SS are not bound at a high level of significance.

We will discuss here how structures observed in clusters of galaxies can provide us insight on the formation and evolution of these objects. We will focus primarily on X-ray observations and results from hydrodynamical $N$-body simulations. This paper is based on a talk given at the School of Cosmology Jose Plinio Baptista -- `Cosmological perturbations and Structure Formation' in Ubu/ES, Brazil.

We report the discovery of extended radio emission in the Phoenix cluster (SPT-CL J2344-4243, z=0.596) with the GMRT at 610 MHz. The diffuse emission extends over a region of at least 400-500 kpc and surrounds the central radio source of the Brightest Cluster Galaxy, but does not appear to be directly associated with it. We classify the diffuse emission as a radio mini-halo, making it the currently most distant mini-halo known. Radio mini-halos have been explained by synchrotron emitting particles re-accelerated via turbulence, possibly induced by gas sloshing generated from a minor merger event. Chandra observations show a non-concentric X-ray surface brightness distribution, which is consistent with this sloshing interpretation. The mini-halo has a flux density of $17\pm5$ mJy, resulting in a 1.4 GHz radio power of ($10.4\pm3.5) \times 10^{24}$ W Hz$^{-1}$. The combined cluster emission, which includes the central compact radio source, is also detected in a shallow GMRT 156 MHz observation and together with the 610 MHz data we compute a spectral index of $-0.84\pm0.12$ for the overall cluster radio emission. Given that mini-halos typically have steeper radio spectra than cluster radio galaxies, this spectral index should be taken as an upper limit for the mini-halo.

Structure formation models predict that clusters of galaxies contain numerous massive subhalos. The gravity of a subhalo in a cluster compresses the surrounding intracluster gas and enhances its X-ray emission. We present a simple model, which treats subhalos as slow moving and gasless, for computing this effect. Recent weak lensing measurements by Okabe et al. have determined masses of ~ 10^13 solar masses for three mass concentrations projected within 300 kpc of the center of the Coma Cluster, two of which are centered on the giant elliptical galaxies NGC 4889 and NGC 4874. Adopting a smooth spheroidal beta-model for the gas distribution in the unperturbed cluster, we model the effect of these subhalos on the X-ray morphology of the Coma Cluster, comparing our results to Chandra and XMM-Newton X-ray data. The agreement between the models and the X-ray morphology of the central Coma Cluster is striking. With subhalo parameters from the lensing measurements, the distances of the three subhalos from the Coma Cluster midplane along our line of sight are all tightly constrained. Using the model to fit the subhalo masses for NGC 4889 and NGC 4874 gives 9.1 x 10^12 and 7.6 x 10^12 solar masses, respectively, in good agreement with the lensing masses. These results lend strong support to the argument that NGC 4889 and NGC 4874 are each associated with a subhalo that resides near the center of the Coma Cluster. In addition to constraining the masses and 3-d location of subhalos, the X-ray data show promise as a means of probing the structure of central subhalos.

We present a new method to quantify substructures in clusters of galaxies, based on the analysis of the intensity of structures. This analysis is done in a residual image that is the result of the subtraction of a surface brightness model, obtained by fitting a two-dimensional analytical model (beta-model or S\'ersic profile) with elliptical symmetry, from the X-ray image. Our method is applied to 34 clusters observed by the Chandra Space Telescope that are in the redshift range 0.02<z<0.2 and have a signal-to-noise ratio greater than 100. We present the calibration of the method and the relations between the substructure level with physical quantities, such as the mass, X-ray luminosity, temperature, and cluster redshift. We use our method to separate the clusters in two sub-samples of high and low substructure levels. We conclude, using Monte Carlo simulations, that the method recuperates very well the true amount of substructure for small angular core radii clusters (with respect to the whole image size) and good signal-to-noise observations. We find no evidence of correlation between the substructure level and physical properties of the clusters such as mass, gas temperature, X-ray luminosity and redshift. The scaling relations for the two sub-samples (high and low substructure level clusters) are different (they present an off-set, i.e., given a fixed mass or temperature, low substructure clusters tend to be more X-ray luminous), which is an important result for cosmological tests using the mass-luminosity relation to obtain the cluster mass function, since they rely on the assumption that clusters do not present different scaling relations according to their dynamical state.

X-ray data analysis have found that fairly complex structures at cluster centres are more common than expected. Many of these structures have similar morphologies, which exhibit spiral-like substructure. It is not yet well known how these structures formed or are maintained. Understanding the origin of these spiral-like features at the centre of some clusters is the major motivation behind this work. We analyse deep \textit{Chandra} observations of 15 nearby galaxy clusters (0.01 $ < z < $ 0.06), and use X-ray temperature and substructure maps to detect small features at the cores of the clusters. We detect spiral-like features at the centre of 7 clusters: A85, A426, A496, Hydra A cluster, Centaurus, Ophiuchus, and A4059. These patterns are similar to those found in numerical hydrodynamic simulations of cluster mergers with non-zero impact parameter. In some clusters of our sample, a strong radio source also occupies the inner region of the cluster, which indicates a possible connection between the two. Our investigation implies that these spiral-like structures may be caused by off-axis minor mergers. Since these features occur in regions of high density, they may confine radio emission from the central galaxy producing, in some cases, unusual radio morphology.

The luminous material in clusters of galaxies falls primarily into two forms: the visible galaxies and the X-ray emitting intra-cluster medium. The hot intra-cluster gas is the major observed baryonic component of clusters, about six times more massive than the stellar component. The mass contained within visible galaxies amounts to approximately 3% of the dynamical mass. Our aim was to analyze both baryonic components, combining X-ray and optical data of a sample of five galaxy clusters (Abell 496, 1689, 2050, 2631 and 2667), within the redshift range 0.03 < z < 0.3. We determined the contribution of stars in galaxies and the intra-cluster medium to the total baryon budget. We used public XMM-Newton data to determine the gas mass and to obtain the X-ray substructures. Using the optical counterparts from SDSS or CFHT we determined the stellar contribution. We examine the relative contribution of galaxies, intra-cluster light and intra-cluster medium to baryon budget in clusters through the stellar-to-gas mass ratio, estimated with use of recent data. We find that the stellar-to-gas mass ratio within r_500 (the radius which the mean cluster density exceeds the critical density by a factor of 500), is anti-correlated with the ICM temperature, ranging from 24% to 6% whereas the temperature ranges from 4.0 to 8.3 keV. This indicates that less massive cold clusters are more prolific star forming environments than massive hot clusters.