skip to main content


Title: A new step forward in realistic cluster lens mass modelling: analysis of Hubble Frontier Field Cluster Abell S1063 from joint lensing, X-ray, and galaxy kinematics data
ABSTRACT

We present a new method to simultaneously and self-consistently model the mass distribution of galaxy clusters that combines constraints from strong lensing features, X-ray emission, and galaxy kinematics measurements. We are able to successfully decompose clusters into their collisionless and collisional mass components thanks to the X-ray surface brightness, as well as use the dynamics of cluster members, to obtain more accurate masses exploiting the fundamental plane of elliptical galaxies. Knowledge from all observables is included through a consistent Bayesian approach in the likelihood or in physically motivated priors. We apply this method to the galaxy cluster Abell S1063 and produce a mass model that we publicly release with this paper. The resulting mass distribution presents different ellipticities for the intra-cluster gas and the other large-scale mass components as well as deviation from elliptical symmetry in the main halo. We assess the ability of our method to recover the masses of the different elements of the cluster using a mock cluster based on a simplified version of our Abell S1063 model. Thanks to the wealth of mutliwavelength information provided by the mass model and the detected X-ray emission, we also found evidence for an ongoing merger event with gas sloshing from a smaller infalling structure into the main cluster. In agreement with previous findings, the total mass, gas profile, and gas mass fraction are all consistent with small deviations from the hydrostatic equilibrium. This new mass model for Abell S1063 is publicly available, as the lenstool extension used to construct it.

 
more » « less
NSF-PAR ID:
10475132
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
527
Issue:
2
ISSN:
0035-8711
Format(s):
Medium: X Size: p. 3246-3275
Size(s):
["p. 3246-3275"]
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Context. Inferences about dark matter, dark energy, and the missing baryons all depend on the accuracy of our model of large-scale structure evolution. In particular, with cosmological simulations in our model of the Universe, we trace the growth of structure, and visualize the build-up of bigger structures from smaller ones and of gaseous filaments connecting galaxy clusters. Aims. Here we aim to reveal the complexity of the large-scale structure assembly process in great detail and on scales from tens of kiloparsecs up to more than 10 Mpc with new sensitive large-scale observations from the latest generation of instruments. We also aim to compare our findings with expectations from our cosmological model. Methods. We used dedicated SRG/eROSITA performance verification (PV) X-ray, ASKAP/EMU Early Science radio, and DECam optical observations of a ~15 deg 2 region around the nearby interacting galaxy cluster system A3391/95 to study the warm-hot gas in cluster outskirts and filaments, the surrounding large-scale structure and its formation process, the morphological complexity in the inner parts of the clusters, and the (re-)acceleration of plasma. We also used complementary Sunyaev-Zeldovich (SZ) effect data from the Planck survey and custom-made Galactic total (neutral plus molecular) hydrogen column density maps based on the HI4PI and IRAS surveys. We relate the observations to expectations from cosmological hydrodynamic simulations from the Magneticum suite. Results. We trace the irregular morphology of warm and hot gas of the main clusters from their centers out to well beyond their characteristic radii, r 200 . Between the two main cluster systems, we observe an emission bridge on large scale and with good spatial resolution. This bridge includes a known galaxy group but this can only partially explain the emission. Most gas in the bridge appears hot, but thanks to eROSITA’s unique soft response and large field of view, we discover some tantalizing hints for warm, truly primordial filamentary gas connecting the clusters. Several matter clumps physically surrounding the system are detected. For the “Northern Clump,” we provide evidence that it is falling towards A3391 from the X-ray hot gas morphology and radio lobe structure of its central AGN. Moreover, the shapes of these X-ray and radio structures appear to be formed by gas well beyond the virial radius, r 100 , of A3391, thereby providing an indirect way of probing the gas in this elusive environment. Many of the extended sources in the field detected by eROSITA are also known clusters or new clusters in the background, including a known SZ cluster at redshift z = 1. We find roughly an order of magnitude more cluster candidates than the SPT and ACT surveys together in the same area. We discover an emission filament north of the virial radius of A3391 connecting to the Northern Clump. Furthermore, the absorption-corrected eROSITA surface brightness map shows that this emission filament extends south of A3395 and beyond an extended X-ray-emitting object (the “Little Southern Clump”) towards another galaxy cluster, all at the same redshift. The total projected length of this continuous warm-hot emission filament is 15 Mpc, running almost 4 degrees across the entire eROSITA PV observation field. The Northern and Southern Filament are each detected at >4 σ . The Planck SZ map additionally appears to support the presence of both new filaments. Furthermore, the DECam galaxy density map shows galaxy overdensities in the same regions. Overall, the new datasets provide impressive confirmation of the theoretically expected structure formation processes on the individual system level, including the surrounding warm-hot intergalactic medium distribution; the similarities of features found in a similar system in the Magneticum simulation are striking. Our spatially resolved findings show that baryons indeed reside in large-scale warm-hot gas filaments with a clumpy structure. 
    more » « less
  2. null (Ed.)
    The presence of relativistic electrons within the diffuse gas phase of galaxy clusters is now well established, thanks to deep radio observations obtained over the last decade, but their detailed origin remains unclear. Cosmic ray protons are also expected to accumulate during the formation of clusters. They may explain part of the radio signal and would lead to γ -ray emission through hadronic interactions within the thermal gas. Recently, the detection of γ -ray emission has been reported toward the Coma cluster with Fermi -LAT. Assuming that this γ -ray emission arises essentially from pion decay produced in proton-proton collisions within the intracluster medium (ICM), we aim at exploring the implication of this signal on the cosmic ray content of the Coma cluster and comparing it to observations at other wavelengths. We use the MINOT software to build a physical model of the Coma cluster, which includes the thermal target gas, the magnetic field strength, and the cosmic rays, to compute the corresponding expected γ -ray signal. We apply this model to the Fermi -LAT data using a binned likelihood approach, together with constraints from X-ray and Sunyaev-Zel’dovich observations. We also consider contamination from compact sources and the impact of various systematic effects on the results. We confirm that a significant γ -ray signal is observed within the characteristic radius θ 500 of the Coma cluster, with a test statistic TS ≃ 27 for our baseline model. The presence of a possible point source (4FGL J1256.9+2736) may account for most of the observed signal. However, this source could also correspond to the peak of the diffuse emission of the cluster itself as it is strongly degenerate with the expected ICM emission, and extended models match the data better. Given the Fermi -LAT angular resolution and the faintness of the signal, it is not possible to strongly constrain the shape of the cosmic ray proton spatial distribution when assuming an ICM origin of the signal, but preference is found in a relatively flat distribution elongated toward the southwest, which, based on data at other wavelengths, matches the spatial distribution of the other cluster components well. Assuming that the whole γ -ray signal is associated with hadronic interactions in the ICM, we constrain the cosmic ray to thermal energy ratio within R 500 to X CRp = 1.79 −0.30 +1.11 % and the slope of the energy spectrum of cosmic rays to α = 2.80 −0.13 +0.67 ( X CRp = 1.06 −0.22 +0.96 % and α = 2.58 −0.09 +1.12 when including both the cluster and 4FGL J1256.9+2736 in our model). Finally, we compute the synchrotron emission associated with the secondary electrons produced in hadronic interactions assuming steady state. This emission is about four times lower than the overall observed radio signal (six times lower when including 4FGL J1256.9+2736), so that primary cosmic ray electrons or reacceleration of secondary electrons is necessary to explain the total emission. We constrain the amplitude of the primary to secondary electrons, or the required boost from reacceleration with respect to the steady state hadronic case, depending on the scenario, as a function of radius. Our results confirm that γ -ray emission is detected in the direction of the Coma cluster. Assuming that the emission is due to hadronic interactions in the intracluster gas, they provide the first quantitative measurement of the cosmic ray proton content in a galaxy cluster and its implication for the cosmic ray electron populations. 
    more » « less
  3. ABSTRACT

    We present a multiwavelength observation of a cool core that does not appear to be associated with any galaxy, in a nearby cluster, Abell 1142. Its X-ray surface brightness peak of ≲2 keV is cooler than the ambient intracluster gas of ≳3 keV, and is offset from its brightest cluster galaxy (BCG) by 80 kpc in projection, representing the largest known cool core – BCG separation. This BCG-less cool core allows us to measure the metallicity of a cluster centre with a much-reduced contribution from the interstellar medium (ISM) of the BCG. XMM–Newton observation reveals a prominent Fe abundance peak of $1.07^{+0.16}_{-0.15}$ Z⊙ and an α/Fe abundance ratio close to the solar ratio, fully consistent with those found at the centres of typical cool core clusters. This finding hints that BCGs play a limited role in enriching the cluster centres. However, the discussion remains open, given that the α/Fe abundance ratios of the orphan cool core and the BCG ISM are not significantly different. Abell 1142 may have experienced a major merger more than 100 Myr ago, which has dissociated its cool core from the BCG. This implies that the Fe abundance peak in cool core clusters can be resilient to cluster mergers. Our recent Institut de Radio Astronomie Millimétrique 30-m observation did not detect any CO emission at its X-ray peak and we find no evidence for massive runaway cooling in the absence of recent active galactic nucleus feedback. The lack of a galaxy may contribute to an inefficient conversion of the ionized warm gas to the cold molecular gas.

     
    more » « less
  4. null (Ed.)
    Context. The excitation of the filamentary gas structures surrounding giant elliptical galaxies at the center of cool-core clusters, also known as brightest cluster galaxies (BCGs), is key to our understanding of active galactic nucleus (AGN) feedback, and of the impact of environmental and local effects on star formation. Aims. We investigate the contribution of thermal radiation from the cooling flow surrounding BCGs to the excitation of the filaments. We explore the effects of small levels of extra heating (turbulence), and of metallicity, on the optical and infrared lines. Methods. Using the C LOUDY code, we modeled the photoionization and photodissociation of a slab of gas of optical depth A V  ≤ 30 mag at constant pressure in order to calculate self-consistently all of the gas phases, from ionized gas to molecular gas. The ionizing source is the extreme ultraviolet (EUV) and soft X-ray radiation emitted by the cooling gas. We tested these models comparing their predictions to the rich multi-wavelength observations from optical to submillimeter, now achieved in cool core clusters. Results. Such models of self-irradiated clouds, when reaching sufficiently large A V , lead to a cloud structure with ionized, atomic, and molecular gas phases. These models reproduce most of the multi-wavelength spectra observed in the nebulae surrounding the BCGs, not only the low-ionization nuclear emission region like optical diagnostics, [O  III ] λ 5007 Å/H β , [N  II ] λ 6583 Å/H α , and ([S  II ] λ 6716 Å+[S  II ] λ 6731 Å)/H α , but also the infrared emission lines from the atomic gas. [O  I ] λ 6300 Å/H α , instead, is overestimated across the full parameter space, except for very low A V . The modeled ro-vibrational H 2 lines also match observations, which indicates that near- and mid-infrared H 2 lines are mostly excited by collisions between H 2 molecules and secondary electrons produced naturally inside the cloud by the interaction between the X-rays and the cold gas in the filament. However, there is still some tension between ionized and molecular line tracers (i.e., CO), which requires optimization of the cloud structure and the density of the molecular zone. The limited range of parameters over which predictions match observations allows us to constrain, in spite of degeneracies in the parameter space, the intensity of X-ray radiation bathing filaments, as well as some of their physical properties like A V or the level of turbulent heating rate. Conclusions. The reprocessing of the EUV and X-ray radiation from the plasma cooling is an important powering source of line emission from filaments surrounding BCGs. C LOUDY self-irradiated X-ray excitation models coupled with a small level of turbulent heating manage to simultaneously reproduce a large number of optical-to-infrared line ratios when all the gas phases (from ionized to molecular) are modeled self-consistently. Releasing some of the simplifications of our model, like the constant pressure, or adding the radiation fields from the AGN and stars, as well as a combination of matter- and radiation-bounded cloud distribution, should improve the predictions of line emission from the different gas phases. 
    more » « less
  5. ABSTRACT

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

     
    more » « less