Search

Strong gravitational lensing offers a powerful probe of the detailed distribution of matter in lenses, while magnifying and bringing faint background sources into view. Observed strong lensing by massive galaxy clusters, which are often in complex dynamical states, has also been used to map their dark matter substructures on smaller scales. Deep high resolution imaging has revealed the presence of strong lensing events associated with these substructures, namely galaxy-scale sub-halos. However, an inventory of these observed galaxy-galaxy strong lensing (GGSL) events is noted to be discrepant with state-of-the-art {\Lambda}CDM simulations. Cluster sub-halos appear to be over-concentrated compared to their simulated counterparts yielding an order of magnitude higher value of GGSL. In this paper, we explore the possibility of resolving this observed discrepancy by redistributing the mass within observed cluster sub-halos in ways that are consistent within the {\Lambda}CDM paradigm of structure formation. Lensing mass reconstructions from data provide constraints on the mass enclosed within apertures and are agnostic to the detailed mass profile within them. Therefore, as the detailed density profile within cluster sub-halos currently remains unconstrained by data, we are afforded the freedom to redistribute the enclosed mass. We investigate if rearranging the mass to a more centrally concentrated density profile helps alleviate the GGSL discrepancy. We report that refitting cluster sub-halos to the ubiquitous {\Lambda}CDM-motivated Navarro-Frenk-White profile, and further modifying them to include significant baryonic components, does not resolve this tension. A resolution to this persisting GGSL discrepancy may require more careful exploration of alternative dark matter models., Comment: 20 pages, 5 figures

This work aims at assessing the impact of DM self-interactions on the properties of galaxy clusters. In particular, the goal is to study the angular dependence of the cross section by testing rare (large angle scattering) and frequent (small angle scattering) SIDM models with velocity-dependent cross sections. We re-simulate six galaxy cluster zoom-in initial conditions with a dark matter only run and with a full-physics setup simulations that includes a self-consistent treatment of baryon physics. We test the dark matter only setup and the full physics setup with either collisionless cold dark matter, rare self-interacting dark matter, and frequent self-interacting dark matter models. We then study their matter density profiles as well as their subhalo population. Our dark matter only SIDM simlations agree with theoretical models, and when baryons are included in simulations, our SIDM models substantially increase the central density of galaxy cluster cores compared to full-physics simulations using collisionless dark matter. SIDM subhalo suppression in full-physics simulations is milder compared to the one found in dark matter only simulations, because of the cuspier baryionic potential that prevent subhalo disruption. Moreover SIDM with small-angle scattering significantly suppress a larger number of subhaloes compared to large angle scattering SIDM models. Additionally, SIDM models generate a broader range of subhalo concentration values, including a tail of more diffuse subhaloes in the outskirts of galaxy clusters and a population of more compact subhaloes in the cluster cores., Comment: 14 pages, 15 figures, abstract abridged to fit arXiv character limit, submitted to Astronomy and Astrophysics

The Local Volume Complete Cluster Survey (LoVoCCS) is an on-going program to observe nearly a hundred low-redshift X-ray-luminous galaxy clusters (redshifts $0.0310^{44}$ erg/s) with the Dark Energy Camera (DECam), capturing data in $u,g,r,i,z$ bands with a $5\sigma$ point source depth of approximately 25-26th AB magnitudes. Here, we map the aperture masses in 58 galaxy cluster fields using weak gravitational lensing. These clusters span a variety of dynamical states, from nearly relaxed to merging systems, and approximately half of them have not been subject to detailed weak lensing analysis before. In each cluster field, we analyze the alignment between the 2D mass distribution described by the aperture mass map, the 2D red-sequence (RS) galaxy distribution, and the brightest cluster galaxy (BCG). We find that the orientations of the BCG and the RS distribution are strongly aligned throughout the interiors of the clusters: the median misalignment angle is 19 deg within 2 Mpc. We also observe the alignment between the orientations of the RS distribution and the overall cluster mass distribution (by a median difference of 32 deg within 1 Mpc), although this is constrained by galaxy shape noise and the limitations of our cluster sample size. These types of alignment suggest long-term dynamical evolution within the clusters over cosmic timescales., Comment: 40 pages, 16 figures, 5 tables; revised and accepted for publication in ApJ

We aim at a direct measurement of the compactness of three galaxy-scale lenses in massive clusters, testing the accuracy of the scaling laws that describe the members in strong lensing (SL) models of galaxy clusters. We selected the multiply imaged sources MACS J0416.1$-$2403 ID14 ($z=3.221$), MACS J0416.1$-$2403 ID16 ($z=2.095$), and MACS J1206.2$-$0847 ID14 ($z=3.753$). Eight images were observed for the first SL system, and six for the latter two. We focused on the main deflector of each galaxy-scale SL system (identified as members 8971, 8785, and 3910, respectively), and modelled its total mass distribution with a truncated isothermal sphere. We accounted for the lensing effects of the remaining cluster components, and included the uncertainty on the cluster-scale mass distribution through a bootstrapping procedure. We measured a truncation radius value of $6.1^{+2.3}_{-1.1} \, \mathrm{kpc}$, $4.0^{+0.6}_{-0.4} \, \mathrm{kpc}$, and $5.2^{+1.3}_{-1.1} \, \mathrm{kpc}$ for members 8971, 8785, and 3910, respectively. Alternative non-truncated models with a higher number of free parameters do not lead to an improved description of the SL system. We measured the stellar-to-total mass fraction within the effective radius $R_e$ for the three members, finding $0.51\pm0.21$, $1.0\pm0.4$, and $0.39\pm0.16$, respectively. We find that a parameterisation of the properties of cluster galaxies in SL models based on power-law scaling relations with respect to the total luminosity cannot accurately describe their compactness over their full total mass range. Our results agree with modelling of the cluster members based on the Fundamental Plane relation. Finally, we report good agreement between our values of the stellar-to-total mass fraction within $R_e$ and those of early-type galaxies from the SLACS Survey. Our work significantly extends the regime of the current samples of lens galaxies., Comment: Accepted for publication in Astronomy and Astrophysics, 15 pages, 12 figures, 8 tables

The degree of mass loss, i.e. the fraction of stars lost by globular clusters, and specifically by their different populations, is still poorly understood. Many scenarios of the formation of multiple stellar populations, especially the ones involving self-enrichment, assume that the first generation (FG) was more massive at birth than now to reproduce the current mass of the second generation (SG). This assumption implies that, during their long-term evolution, clusters lose around $90\%$ of the FG. We have tested whether such strong mass loss could take place in a massive globular cluster orbiting the Milky Way at $4\ {\rm kpc}$ from the centre and composed of two generations. We perform a series of $N$-body simulations for ${12\ \rm Gyr}$ to probe the parameter space of internal cluster properties. We have derived that, for an extended FG and a low-mass second one, the cluster loses almost $98\%$ of its initial FG mass and the cluster mass can be as much as 20 times lower after a Hubble time. Furthermore, under these conditions, the derived fraction of SG stars, $f_{\rm enriched}$, falls in the range occupied by observed clusters of similar mass ($\sim 0.6-0.8$). In general, the parameters that affect the most the degree of mass loss are the presence or not of primordial segregation, the depth of the central potential, $W_{0,FG}$, the initial mass of the SG, $M^{ini}_{SG}$, and the initial half-mass radius of the SG, $r_{h,SG}$. Higher $M^{ini}_{SG}$ have not been found to imply higher final $f_{\rm enriched}$ due to the deeper cluster potential well which slows down mass loss., Comment: 16 pages, 9 figures. Accepted for publication on Astronomy & Astrophysics

In this study, we investigate a recent finding based on strong lensing observations, which suggests that the sub-halos observed in clusters exhibit greater compactness compared to those predicted by $\Lambda$CDM simulations. To address this discrepancy, we performed a comparative analysis by comparing the cumulative mass function of sub-halos and the $M_{\text{sub}}$-$V_{\text{circ}}$ relation between observed clusters and 324 simulated clusters from The Three Hundred project, focusing on re-simulations using GADGET-X and GIZMO-SIMBA baryonic models. The sub-halos' cumulative mass function of the GIZMO-SIMBA simulated clusters agrees with observations, while the GADGET-X simulations exhibit discrepancies in the lower sub-halo mass range possibly due to its strong SuperNova feedback. Both GADGET-X and GIZMO-SIMBA simulations demonstrate a redshift evolution of the sub-halo mass function and the $V_{max}$ function, with slightly fewer sub-halos observed at lower redshifts. Neither the GADGET-X nor GIZMO-SIMBA(albeit a little closer) simulated clusters' predictions for the $M_{\text{sub}}$-$V_{\text{circ}}$ relation align with the observational result. Further investigations on the correlation between sub-halo/halo properties and the discrepancy in the $M_{\text{sub}}$-$V_{\text{circ}}$ relation reveals that the sub-halo's half mass radius and galaxy stellar age, the baryon fraction and sub-halo distance from the cluster's centre, as well as the halo relaxation state play important roles on this relation. Nevertheless, we think it is still challenging in accurately reproducing the observed $M_{\text{sub}}$-$V_{\text{circ}}$ relation in our current hydrodynamic cluster simulation under the standard $\Lambda$CDM cosmology., Comment: see the companion paper, Meneghetti et al. (2023), on today's arxiv list

We introduce a new, multi-zone chemical evolution model of the DustPedia galaxy M74, calibrated by means of MCMC methods. We take into account the observed stellar and gas density profiles and use Bayesian analysis to constrain two fundamental parameters characterising the gas accretion and star formation timescale, i.e. the infall timescale tau and the SF efficiency nu, respectively, as a function of galactocentric radius R. Our analysis supports an infall timescale increasing with R and a star formation efficiency decreasing with R, thus supporting an 'Inside-Out' formation for M74. For both tau and nu, we find a weaker radial dependence than in the Milky Way. We also investigate the dust content of M74, comparing the observed dust density profile with the results of our chemical evolution models. Various prescriptions have been considered for two key parameters, i.e. the typical dust accretion timescale and the mass of gas cleared out of the dust by a supernova remnant, regulating the dust growth and destruction rate, respectively. Two models with a different current balance between destruction and accretion, i.e. with equilibrium and dominion of accretion over destruction, can equally reproduce the observed dust profile of M74. This outlines the degeneracy between these parameters in shaping the interstellar dust content in galaxies. Our methods will be extended to more DustPedia galaxies to shed more light on the relative roles of dust production and destruction., Comment: MNRAS, accepted for publication, 19 pages, 14 figures

We present a new high-precision, JWST-based, strong lensing model for the galaxy cluster Abell 2744 at $z=0.3072$. By combining the deep, high-resolution JWST imaging from the GLASS-JWST and UNCOVER programs and a Director's Discretionary Time program, with newly obtained VLT/MUSE data, we identify 32 multiple images from 11 background sources lensed by two external sub-clusters at distances of ~160" from the main cluster. The new MUSE observations enable the first spectroscopic confirmation of a multiple image system in the external clumps. Moreover, the re-analysis of the spectro-photometric archival and JWST data yields 27 additional multiple images in the main cluster. The new lens model is constrained by 149 multiple images ($\sim66\%$ more than in our previous Bergamini et al. 2023 model) covering an extended redshift range between 1.03 and 9.76. The subhalo mass component of the cluster includes 177 member galaxies down to $m_{\rm F160W}=21$, 163 of which are spectroscopically confirmed. Internal velocity dispersions are measured for 85 members. The new lens model is characterized by a remarkably low scatter between predicted and observed positions of the multiple images (0.43"). This precision is unprecedented given the large multiple image sample, the complexity of the cluster mass distribution, and the large modeled area. The improved accuracy and resolution of the cluster total mass distribution provides a robust magnification map over a $\sim\!45$ arcmin$^2$ area, which is critical for inferring the intrinsic physical properties of the highly magnified, high-$z$ sources. The lens model and the new MUSE redshift catalog are released with this publication., Comment: Submitted to ApJ., 23 pages, 10 figures, 3 tables

We report the detection of a high density of redshift $z\approx 10$ galaxies behind the foreground cluster Abell 2744, selected from imaging data obtained recently with NIRCam onboard {\it JWST} by three programs -- GLASS-JWST, UNCOVER, and DDT\#2756. To ensure robust estimates of the lensing magnification $\mu$, we use an improved version of our model that exploits the first epoch of NIRCam images and newly obtained MUSE spectra, and avoids regions with $\mu>5$ where the uncertainty may be higher. We detect seven bright $z\approx 10$ galaxies with demagnified rest-frame $-22 \lesssim M_{\rm UV}\lesssim -19$ mag, over an area of $\sim37$ sq. arcmin. Taking into account photometric incompleteness and the effects of lensing on luminosity and cosmological volume, we find that the density of $z\approx 10$ galaxies in the field is about $10\times$ ($3\times$) larger than the average at $M_{UV}\approx -21~ (-20)$ mag reported so far. The density is even higher when considering only the GLASS-JWST data, which are the deepest and the least affected by magnification and incompleteness. The GLASS-JWST field contains 5 out of 7 galaxies, distributed along an apparent filamentary structure of 2 Mpc in projected length, and includes a close pair of candidates with $M_{\rm UV}< -20$ mag having a projected separation of only 16 kpc. These findings suggest the presence of a $z\approx 10$ overdensity in the field. In addition to providing excellent targets for efficient spectroscopic follow-up observations, our study confirms the high density of bright galaxies observed in early {\it JWST} observations, but calls for multiple surveys along independent lines of sight to achieve an unbiased estimate of their average density and a first estimate of their clustering., Comment: Accepted for publication in ApJL, 13 pages, 4 figures

We use hydrodynamical cosmological simulations to test the differences between cold and self-interacting dark matter models (CDM and SIDM) in the mass range of massive galaxies ($10^{12}M_{\odot}h^{-1}10^{13}M_{\odot}h^{-1}$) show cored profiles, while the least massive ones ($M<10^{12.5}M_{\odot}h^{-1}$) have cuspier profiles. Finally, we compare the predictions of our simulations to observational results, by looking at the dark matter fractions and the distribution of strong lensing Einstein radii. We find that in SIDM the DM-fractions decrease more rapidly with increasing stellar mass than in CDM, leading to lower fractions at $M_{*}>10^{11}M_{\odot}$, a distinctive signature of SIDM. At the same time, the distribution of Einstein radii, derived from both CDM and SIDM hydro runs, is comparable to observed samples of strong lenses with $M>10^{13}M_{\odot}h^{-1}$. We conclude that the interplay between self-interaction and baryons can greatly reduce the expected differences between CDM and SIDM models at this mass scale, and that techniques able to separate the dark and luminous mass in the inner regions of galaxies are needed to constrain self-interactions., Comment: 13 pages, 12 figures

In this paper, we assess the impact of numerical resolution and of the implementation of energy input from AGN feedback models on the inner structure of cluster sub-haloes in hydrodynamic simulations. We compare several zoom-in re-simulations of a sub-sample of the cluster-sized haloes studied in Meneghetti et al. (2020), obtained by varying mass resolution, softening length and AGN energy feedback scheme. We study the impact of these different setups on the subhalo abundances, their radial distribution, their density and mass profiles and the relation between the maximum circular velocity, which is a proxy for subhalo compactness. Regardless of the adopted numerical resolution and feedback model, subhaloes with masses Msub < 1e11Msun/h, the most relevant mass-range for galaxy-galaxy strong lensing, have maximum circular velocities ~30% smaller than those measured from strong lensing observations of Bergamini et al. (2019). We also find that simulations with less effective AGN energy feedback produce massive subhaloes (Msub> 1e11 Msun/h ) with higher maximum circular velocity and that their Vmax - Msub relation approaches the observed one. However the stellar-mass number count of these objects exceeds the one found in observations and we find that the compactness of these simulated subhaloes is the result of an extremely over-efficient star formation in their cores, also leading to larger-than-observed subhalo stellar mass. We conclude that simulations are unable to simultaneously reproduce the observed stellar masses and compactness (or maximum circular velocities) of cluster galaxies. Thus, the discrepancy between theory and observations that emerged from the analysis of Meneghetti et al. (2020) persists. It remains an open question as to whether such a discrepancy reflects limitations of the current implementation of galaxy formation models or the LCDM paradigm., Comment: 11 pages, 10 figures, abstract is redacted to fit arXiv character count limit

Meneghetti et al. (2020) recently reported an excess of galaxy-galaxy strong lensing (GGSL) in galaxy clusters compared to expectations from the LCDM cosmological model. Theoretical estimates of the GGSL probability are based on the analysis of numerical hydrodynamical simulations in the LCDM cosmology. We quantify the impact of the numerical resolution and AGN feedback scheme adopted in cosmological simulations on the predicted GGSL probability and determine if varying these simulation properties can alleviate the gap with observations. We repeat the analysis of Meneghetti et al. (2020) on cluster-size halos simulated with different mass and force resolutions and implementing several independent AGN feedback schemes. We find that improving the mass resolution by a factor of ten and twenty-five, while using the same galaxy formation model that includes AGN feedback, does not affect the GGSL probability. We find similar results regarding the choice of gravitational softening. On the contrary, adopting an AGN feedback scheme that is less efficient at suppressing gas cooling and star formation leads to an increase in the GGSL probability by a factor between three and six. However, we notice that such simulations form overly massive subhalos whose contribution to the lensing cross-section would be significant while their Einstein radii are too large to be consistent with the observations. The primary contributors to the observed GGSL cross-sections are subhalos with smaller masses, that are compact enough to become critical for lensing. The population with these required characteristics appears to be absent in simulations., Comment: 13 pages, 11 figures. Submitted for publication on Astronomy and Astrophysics

We present a strong lensing analysis of the galaxy cluster SDSS J1029+2623 at $z=0.588$, one of the few currently known lens clusters with multiple images of a background ($z=2.1992$) quasar with a measured time delay. We use archival Hubble Space Telescope multi-band imaging and new Multi Unit Spectroscopic Explorer follow-up spectroscopy to build an accurate lens mass model, a crucial step towards future cosmological applications. The spectroscopic data enable the secure identification of 57 cluster members and of two nearby perturbers along the line-of-sight. We estimate the inner kinematics of a sub-set of 20 cluster galaxies to calibrate the scaling relations parametrizing the sub-halo mass component. We also reliably determine the redshift of 4 multiply imaged sources, provide a tentative measurement for one system, and report the discovery of a new four-image system. The final catalog comprises 26 multiple images from 7 background sources, spanning a wide redshift range, from 1.02 to 5.06. We present two parametric lens models, with slightly different cluster mass parametrizations. The observed positions of the multiple images are accurately reproduced within approximately $0''.2$, the three image positions of the quasar within only $\sim0''.1$. We estimate a cluster projected total mass of $M(<300~ {\rm kpc}) \sim 2.1 \times 10^{14}~ M_{\odot}$, with a statistical uncertainty of a few percent. Both models, that include a small galaxy close to one of the quasar images, predict magnitude differences and time delays between the quasar images that are consistent with the observations., Comment: 24 pages, 10 figures, 6 Tables. Submitted to ApJ. The full MUSE spectroscopic catalog of SDSS J1029+2623 presented here is made publicly available at https://www.fe.infn.it/astro/lensing

We present a comparison of the nitrogen-to-oxygen ratio (N/O) in 37 high-redshift galaxies at $z\sim$2 taken from the KMOS Lensed Emission Lines and VElocity Review (KLEVER) Survey with a comparison sample of local galaxies, taken from the Sloan Digital Sky Survey (SDSS). The KLEVER sample shows only a mild enrichment in N/O of $+$0.1 dex when compared to local galaxies at a given gas-phase metallicity (O/H), but shows a depletion in N/O of $-$0.36 dex when compared at a fixed stellar mass (M$_*$). We find a strong anti-correlation in local galaxies between N/O and SFR in the M$_*$-N/O plane, similar to the anti-correlation between O/H and SFR found in the mass-metallicity relation (MZR). We use this anti-correlation to construct a fundamental nitrogen relation (FNR), analogous to the fundamental metallicity relation (FMR). We find that KLEVER galaxies are consistent with both the FMR and the FNR. This suggests that the depletion of N/O in high-$z$ galaxies when considered at a fixed M$_*$ is driven by the redshift-evolution of the mass-metallicity relation in combination with a near redshift-invariant N/O-O/H relation. Furthermore, the existence of an fundamental nitrogen relation suggests that the mechanisms governing the fundamental metallicity relation must be probed by not only O/H, but also N/O, suggesting pure-pristine gas inflows are not the primary driver of the FMR, and other properties such as variations in galaxy age and star formation efficiency must be important., Comment: 23 pages, 11 figures, Accepted for publication in Monthly Notices of the Royal Astronomical Society

Galaxy clusters have a triaxial matter distribution. The weak-lensing signal, an important part in cosmological studies, measures the projected mass of all matter along the line-of-sight, and therefore changes with the orientation of the cluster. Studies suggest that the shape of the brightest cluster galaxy (BCG) in the centre of the cluster traces the underlying halo shape, enabling a method to account for projection effects. We use 324 simulated clusters at four redshifts between 0.1 and 0.6 from `The Three Hundred Project' to quantify correlations between the orientation and shape of the BCG and the halo. We find that haloes and their embedded BCGs are aligned, with an average $\sim$20 degree angle between their major axes. The bias in weak lensing cluster mass estimates correlates with the orientation of both the halo and the BCG. Mimicking observations, we compute the projected shape of the BCG, as a measure of the BCG orientation, and find that it is most strongly correlated to the weak-lensing mass for relaxed clusters. We also test a 2-dimensional cluster relaxation proxy measured from BCG mass isocontours. The concentration of stellar mass in the projected BCG core compared to the total stellar mass provides an alternative proxy for the BCG orientation. We find that the concentration does not correlate to the weak-lensing mass bias, but does correlate with the true halo mass. These results indicate that the BCG shape and orientation for large samples of relaxed clusters can provide information to improve weak-lensing mass estimates., Comment: 15 pages, 9 figures, Figure 7 is key plot, Updated to match version accepted by MNRAS

We perform a systematic study of the recently discovered fundamental plane of galaxy clusters (CFP) using ~250 simulated clusters from The Three Hundred project, focusing on the stability of the plane against different temperature definitions and its dependence on the dynamical relaxation state of clusters. The CFP is characterised by $T\propto M_\mathrm{s}^\alpha\,r_\mathrm{s}^\beta$, defined with the gas temperature ($T$) and the characteristic halo scale radius and mass ($r_\mathrm{s}$ and $M_\mathrm{s}$) assuming a Navarro-Frenk-White halo description. We explore two definitions of weighted temperatures, namely mass-weighted and spectroscopic-like temperatures, in three radial ranges. The Three Hundred project clusters at $z=0$ lie on a thin plane whose parameters ($\alpha,\beta$) and dispersion ($0.015-0.030$ dex) depend on the gas temperature definition. The CFP for mass-weighted temperatures is closer to the virial equilibrium expectation ($\alpha=1,\beta=-1$) with a smaller dispersion. For gas temperatures measured within $500h^{-1}$ kpc, the resulting CFP deviates the most from the virial expectation and shifts towards the similarity solution for a secondary infall model ($\alpha=1.5,\beta=-2$). Independently of the temperature definition, we find that clusters at $z=1$ and relaxed clusters form a CFP similar to the virial expectation, unlike disturbed clusters exhibiting stronger evolution. Only systems formed over the last 4 Gyr present a CFP that is closer to the similarity solution. All these findings are compatible with the CFP obtained for a CLASH subsample excluding the hottest clusters with $T_\mathrm{X}>12$ keV., Comment: 19 pages, 10 figures, 9 tables. Accepted for publication in MNRAS

From Hubble Frontier Fields photometry, and data from the Multi Unit Spectroscopic Explorer on the Very Large Telescope, we build the Fundamental Plane (FP) relation for the early-type galaxies of the cluster Abell S1063. We use this relation to develop an improved strong lensing model of the total mass distribution of the cluster, determining the velocity dispersions of all 222 cluster members included in the model from their measured structural parameters. Fixing the hot gas component from X-ray data, the mass density distributions of the diffuse dark matter haloes are optimised by comparing the observed and model-predicted positions of 55 multiple images of 20 background sources, distributed over the redshift range $0.73-6.11$. We determine the uncertainties on the model parameters with Monte Carlo Markov chains. Compared to previous works, our model allows for the inclusion of a scatter on the relation between the total mass and the velocity dispersion of cluster members, which also shows a shallower slope. We notice a lower statistical uncertainty on the value of some parameters, such as the core radius, of the diffuse mass component of the cluster. Thanks to a new estimate of the stellar mass of all members, we measure the cumulative projected mass profiles out to a radius of 350 kpc, for all baryonic and dark matter components of the cluster. At the outermost radius, we find a baryon fraction of $0.147 \pm 0.002$. We compare the sub-haloes as described by our model with recent hydrodynamical cosmological simulations. We find good agreement in terms of the stellar over total mass fraction. On the other hand, we report some discrepancies in terms of the maximum circular velocity, which is an indication of their compactness, and the sub-halo mass function in the central cluster regions., Comment: 15 pages, 13 figures, 5 tables, Astronomy and Astrophysics, in press

We search for unresolved X-ray emission from lensed sources in the FOV of 11 CLASH clusters with Chandra data. We consider the solid angle in the lens plane corresponding to a magnification $\mu>1.5$, that amounts to a total of ~100 arcmin$^2$. Our main goal is to assess the efficiency of massive clusters as cosmic telescopes to explore the faint end of X-ray extragalactic source population. We search for X-ray emission from strongly lensed sources identified in the optical, and perform an untargeted detection of lensed X-ray sources. We detect X-ray emission only in 9 out of 849 lensed/background optical sources. The stacked emission of the sources without detection does not reveal any signal in any band. Based on the untargeted detection, we find 66 additional X-ray sources that are consistent with being lensed sources. After accounting for completeness and sky coverage, we measure for the first time the soft- and hard-band number counts of lensed X-ray sources. The results are consistent with current modelization of the AGN population distribution. The distribution of de-lensed fluxes of the sources identified in moderately deep CLASH fields reaches a flux limit of ~$10^{-16}$ and ~$10^{-15}$ erg/s/cm$^{2}$ in the soft and hard bands, respectively. We conclude that, in order to match the depth of the CDFS exploiting massive clusters as cosmic telescopes, the required number of cluster fields is about two orders of magnitude larger than that offered by the 20 years Chandra archive. A significant step forward will be made when future X-ray facilities, with ~1' angular resolution and large effective area, will allow the serendipitous discovery of rare, strongly lensed high-$z$ X-ray sources, enabling the study of faint AGN activity in early Universe and the measurement of gravitational time delays in the X-ray variability of multiply imaged AGN., Comment: Accepted for publication in A&A

We present a new high-precision parametric strong lensing model of the galaxy cluster MACS J0416.1-2403, at z=0.396, which takes advantage of the MUSE Deep Lensed Field (MDLF), with 17.1h integration in the northeast region of the cluster, and Hubble Frontier Fields data. We spectroscopically identify 182 multiple images from 48 background sources at 0.9

Cold dark matter (CDM) constitutes most of the matter in the Universe. The interplay between dark and luminous matter in dense cosmic environments like galaxy clusters is studied theoretically using cosmological simulations. Observed gravitational lensing is used to test and characterize the properties of substructures - the small-scale distribution of dark matter - in clusters. An apt metric, the probability of strong lensing events produced by dark matter substructure, is devised and computed for 11 galaxy clusters. We report that observed cluster substructures are more efficient lenses than predicted by CDM simulations, by more than an order of magnitude. We suggest that hitherto undiagnosed systematic issues with simulations or incorrect assumptions about the properties of dark matter could explain our results., Comment: 56 pages, 15 figures. Published in Science

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

We present near-infrared observations of 42 gravitationally lensed galaxies obtained in the framework of the KMOS Lensed Emission Lines and VElocity Review (KLEVER) Survey, a program aimed at investigating the spatially resolved properties of the ionised gas in 1.2 3$\sigma$) "inverted" gradients are also found, showing an anti-correlation between metallicity and star formation rate density on local scales, possibly suggesting recent episodes of pristine gas accretion or strong radial flows in place. Nevertheless, the individual metallicity maps are characterised by a variety of different morphologies, with flat radial gradients sometimes hiding non-axisymmetric variations on kpc scales which are washed out by azimuthal averages, especially in interacting systems or in those undergoing local episodes of recent star formation., Comment: Main text: 23 pages, 11 figures. Appendix: 12 pages, 11 figures. Accepted for publication on MNRAS

According to the standard cold dark matter (CDM) cosmology, the structure of dark halos including those of galaxy clusters reflects their mass accretion history. Older clusters tend to be more concentrated than younger clusters. Their structure, represented by the characteristic radius $r_s$ and mass $M_s$ of the Navarro--Frenk--White (NFW) density profile, is related to their formation time. In~this study, we showed that $r_s$, $M_s$, and the X-ray temperature of the intracluster medium (ICM), $T_X$, form a thin plane in the space of $(\log r_s, \log M_s, \log T_X)$. This tight correlation indicates that the ICM temperature is also determined by the formation time of individual clusters. Numerical simulations showed that clusters move along the fundamental plane as they evolve. The plane and the cluster evolution within the plane could be explained by a similarity solution of structure formation of the universe. The angle of the plane shows that clusters have not achieved "virial equilibrium" in the sense that mass/size growth and pressure at the boundaries cannot be ignored. The distribution of clusters on the plane was related to the intrinsic scatter in the halo concentration--mass relation, which originated from the variety of cluster ages. The well-known mass--temperature relation of clusters ($M_\Delta\propto T_X^{3/2}$) can be explained by the fundamental plane and the mass dependence of the halo concentration without the assumption of virial equilibrium. The fundamental plane could also be used for calibration of cluster masses., Comment: Invited review article, to be published in "From Dark Haloes to Visible Galaxies", special issue of Galaxies

Observations and numerical simulations have shown that the relation between the mass scaled with the critical density of the universe and the X-ray temperature of galaxy clusters is approximately represented by $M_\Delta \propto T_X^{3/2}$ (e.g. $\Delta=500$). This relation is often interpreted as evidence that clusters are in virial equilibrium. However, the recently discovered fundamental plane (FP) of clusters indicates that the temperature of clusters primarily depends on a combination of the characteristic mass $M_s$ and radius $r_s$ of the Navarro-Frenk-White profile rather than $M_\Delta$. Moreover, the angle of the FP revealed that clusters are not in virial equilibrium because of continuous mass accretion from the surrounding matter. By considering both the FP and the mass dependence of the cluster concentration parameter, we show that this paradox can be solved and the relation $M_\Delta \propto T_X^{3/2}$ actually reflects the central structure of clusters. We also find that the intrinsic scatter in the halo concentration-mass relation can largely account for the spread of clusters on the FP. We also show that X-ray data alone form the FP and the angle and the position are consistent with those of the FP constructed from gravitational lensing data. We demonstrate that a possible shift between the two FPs can be used to calibrate cluster masses obtained via X-ray observations., Comment: Published on ApJ. Matched to published version

We reconstruct the two-dimensional (2D) matter distributions in 20 high-mass galaxy clusters selected from the CLASH survey by using the new approach of performing a joint weak lensing analysis of 2D shear and azimuthally averaged magnification measurements. This combination allows for a complete analysis of the field, effectively breaking the mass-sheet degeneracy. In a Bayesian framework, we simultaneously constrain the mass profile and morphology of each individual cluster assuming an elliptical Navarro-Frenk-White halo characterized by the mass, concentration, projected axis ratio, and position angle of the projected major axis.. We find that spherical mass estimates of the clusters from azimuthally averaged weak-lensing measurements in previous work are in excellent agreement with our results from a full 2D analysis. Combining all 20 clusters in our sample, we detect the elliptical shape of weak-lensing halos at the $5\sigma$ significance level within a scale of 2Mpc$/h$. The median projected axis ratio is $0.67\pm 0.07$ at a virial mass of $M_\mathrm{vir}=(15.2\pm 2.8) \times 10^{14} M_\odot$, which is in agreement with theoretical predictions of the standard collisionless cold dark matter model. We also study misalignment statistics of the brightest cluster galaxy, X-ray, thermal Sunyaev-Zel'dovich effect, and strong-lensing morphologies with respect to the weak-lensing signal. Among the three baryonic tracers studied here, we find that the X-ray morphology is best aligned with the weak-lensing mass distribution, with a median misalignment angle of $21\pm 7$ degrees. We also conduct a stacked quadrupole shear analysis assuming that the X-ray major axis is aligned with that of the projected mass distribution. This yields a consistent axis ratio of $0.67\pm 0.10$, suggesting again a tight alignment between the intracluster gas and dark matter., Comment: Minor changes to match the version published in ApJ. One of three new companion papers of the CLUMP-3D project (I-Non Chiu et al., arXiv:1804.00676; Mauro Sereno et al., arXiv:1804.00667)

In cold dark matter (CDM) cosmology, objects in the Universe have grown under the effect of gravity of dark matter. The intracluster gas in a galaxy cluster was heated when the dark-matter halo formed through gravitational collapse. The potential energy of the gas was converted to thermal energy through this process. However, this process and the thermodynamic history of the gas have not been clearly characterized in connection with with the formation and evolution of the internal structure of dark-matter halos. Here, we show that observational CLASH data of high-mass galaxy clusters lie on a plane in the three-dimensional logarithmic space of their characteristic radius $r_s$, mass $M_s$, and X-ray temperature $T_X$ with a very small orthogonal scatter. The tight correlation indicates that the gas temperature was determined at a specific cluster formation time, which is encoded in $r_s$ and $M_s$. The plane is tilted with respect to $T_X \propto M_s/r_s$, which is the plane expected in case of simplified virial equilibrium. We show that this tilt can be explained by a similarity solution, which indicates that clusters are not isolated but continuously growing through matter accretion from their outer environments. Numerical simulations reproduce the observed plane and its angle. This result holds independently of the gas physics implemented in the code, revealing the fundamental origin of this plane., Comment: Replaced with a revised version to match the ApJ accepted version

We map the lensing-inferred substructure in the first three clusters observed by the Hubble Space Telescope Frontier Fields Initiative (HSTFF): Abell 2744 (z = 0.308), MACSJ0416, (z = 0.396) and MACSJ1149 (z = 0.543). Statistically resolving dark-matter subhaloes down to ~10^{9.5} solar masses, we compare the derived subhalo mass functions (SHMFs) to theoretical predictions from analytical models and with numerical simulations in a Lambda Cold Dark Matter (LCDM) cosmology. Mimicking our observational cluster member selection criteria in the HSTFF, we report excellent agreement in both amplitude and shape of the SHMF over four decades in subhalo mass (10^{9-13} solar masses). Projection effects do not appear to introduce significant errors in the determination of SHMFs from simulations. We do not find evidence for a substructure crisis, analogous to the missing satellite problem in the Local Group, on cluster scales, but rather excellent agreement of the count-matched HSTFF SHMF down to M_{sub halo}/M_{halo} ~ 10^{-5}. However, we do find discrepancies in the radial distribution of sub haloes inferred from HSTFF cluster lenses compared to determinations from simulated clusters. This suggests that although the selected simulated clusters match the HSTFF sample in mass, they do not adequately capture the dynamical properties and complex merging morphologies of these observed cluster lenses. Therefore, HSTFF clusters are likely observed in a transient evolutionary stage that is presently insufficiently sampled in cosmological simulations. The abundance and mass function of dark matter substructure in cluster lenses continues to offer an important test of the LCDM paradigm, and at present we find no tension between model predictions and observations., Comment: 20 pages, 29 figures, Accepted for publication, in press, MNRAS, replaced with added author affiliations and added reference

Morphology is often used to infer the state of relaxation of galaxy clusters. The regularity, symmetry, and degree to which a cluster is centrally concentrated inform quantitative measures of cluster morphology. The Cluster Lensing and Supernova survey with Hubble Space Telescope (CLASH) used weak and strong lensing to measure the distribution of matter within a sample of 25 clusters, 20 of which were deemed to be relaxed based on their X-ray morphology and alignment of the X-ray emission with the BCG. Towards a quantitative characterization of this important sample of clusters, we present uniformly estimated X-ray morphological statistics for all 25 CLASH clusters. We compare X-ray morphologies of CLASH clusters with those identically measured for a large sample of simulated clusters from the MUSIC-2 simulations, selected by mass. We confirm a threshold in X-ray surface brightness concentration of C>0.4 for cool-core clusters, where C is the ratio of X-ray emission inside 100 kpc/h70 compared to inside 500 kpc/h70. We report and compare morphologies of these clusters inferred from Sunyaev-Zeldovich Effect (SZE) maps of the hot gas and in from projected mass maps based on strong and weak lensing. We find a strong agreement in alignments of the orientation of major axes for the lensing, X-ray, and SZE maps of nearly all of the CLASH clusters at radii of 500 kpc (approximately 0.5R500 for these clusters). We also find a striking alignment of clusters shapes at the 500 kpc scale, as measured with X-ray, SZE, and lensing, with that of the near-infrared stellar light at 10 kpc scales for the 20 "relaxed" clusters. This strong alignment indicates a powerful coupling between the cluster- and galaxy-scale galaxy formation processes., Comment: 13 pages. Accepted to ApJ 11 Jan 2016

We present profiles of temperature (Tx), gas mass, and hydrostatic mass estimated from new and archival X-ray observations of CLASH clusters. We compare measurements derived from XMM and Chandra observations with one another and compare both to gravitational lensing mass profiles derived with CLASH HST and ground-based lensing data. Radial profiles of Chandra and XMM electron density and enclosed gas mass are nearly identical, indicating that differences in hydrostatic masses inferred from X-ray observations arise from differences in Tx measurements. Encouragingly, cluster Txs are consistent with one another at ~100-200 kpc radii but XMM Tx systematically decline relative to Chandra Tx at larger radii. The angular dependence of the discrepancy suggests additional investigation on systematics such as the XMM point spread function correction, vignetting and off-axis responses. We present the CLASH-X mass-profile comparisons in the form of cosmology-independent and redshift-independent circular-velocity profiles. Ratios of Chandra HSE mass profiles to CLASH lensing profiles show no obvious radial dependence in the 0.3-0.8 Mpc range. However, the mean mass biases inferred from the WL and SaWLens data are different. e.g., the weighted-mean value at 0.5 Mpc is = 0.12 for the WL comparison and = -0.11 for the SaWLens comparison. The ratios of XMM HSE mass profiles to CLASH lensing profiles show a pronounced radial dependence in the 0.3-1.0 Mpc range, with a weighted mean mass bias of value rising to ~0.3 at ~1 Mpc for the WL comparison and of 0.25 for SaWLens comparison. The enclosed gas mass profiles from both Chandra and XMM rise to a value 1/8 times the total-mass profiles inferred from lensing at 0.5 Mpc and remain constant outside of that radius, suggesting that [8xMgas] profiles may be an excellent proxy for total-mass profiles at >0.5 Mpc in massive galaxy clusters., Comment: Accepted to ApJ; 24 pages; scheduled to appear in the Oct 10, 2014 issue. This version corrects the typographical error in the superscripts for Equation (2) to include the square of (r/r_core). The correct version of this equation was used in the analysis

We present a joint shear-and-magnification weak-lensing analysis of a sample of 16 X-ray-regular and 4 high-magnification galaxy clusters at 0.19

We study how well halo properties of galaxy clusters, like mass and concentration, are recovered using lensing data. In order to generate a large sample of systems at different redshifts we use the code MOKA. We measure halo mass and concentration using weak lensing data alone (WL), fitting to an NFW profile the reduced tangential shear profile, or by combining weak and strong lensing data, by adding information about the size of the Einstein radius (WL+SL). For different redshifts, we measure the mass and the concentration biases and find that these are mainly caused by the random orientation of the halo ellipsoid with respect to the line-of-sight. Since our simulations account for the presence of a bright central galaxy, we perform mass and concentration measurements using a generalized NFW profile which allows for a free inner slope. This reduces both the mass and the concentration biases. We discuss how the mass function and the concentration mass relation change when using WL and WL+SL estimates. We investigate how selection effects impact the measured concentration-mass relation showing that strong lens clusters may have a concentration 20-30% higher than the average, at fixed mass, considering also the particular case of strong lensing selected samples of relaxed clusters. Finally, we notice that selecting a sample of relaxed galaxy clusters, as is done in some cluster surveys, explain the concentration-mass relation biases., Comment: (1) DIFA-UniBO, (2) INAF-OABo, (3) INFN-BO, (4) JPL-Pasadena 18 pages, 19 figures - accepted for publication by MNRAS, two figures added for comparison with SGAS-SDSS and LoCuSS clusters

The Cluster Lensing And Supernova survey with Hubble (CLASH) is a 524-orbit multi-cycle treasury program to use the gravitational lensing properties of 25 galaxy clusters to accurately constrain their mass distributions. The survey, described in detail in this paper, will definitively establish the degree of concentration of dark matter in the cluster cores, a key prediction of CDM. The CLASH cluster sample is larger and less biased than current samples of space-based imaging studies of clusters to similar depth, as we have minimized lensing-based selection that favors systems with overly dense cores. Specifically, twenty CLASH clusters are solely X-ray selected. The X-ray selected clusters are massive (kT > 5 keV; 5 - 30 x 10^14 M_solar) and, in most cases, dynamically relaxed. Five additional clusters are included for their lensing strength (Einstein radii > 35 arcsec at z_source = 2) to further quantify the lensing bias on concentration, to yield high resolution dark matter maps, and to optimize the likelihood of finding highly magnified high-redshift (z > 7) galaxies. The high magnification, in some cases, provides angular resolutions unobtainable with any current UVOIR facility and can yield z > 7 candidates bright enough for spectroscopic follow-up. A total of 16 broadband filters, spanning the near-UV to near-IR, are employed for each 20-orbit campaign on each cluster. These data are used to measure precise (sigma_phz < 0.02(1+z)) photometric redshifts for dozens of newly discovered multiply-lensed images per cluster. Observations of each cluster are spread over 8 epochs to enable a search, primarily in the parallel fields, for Type Ia supernovae at z > 1 to improve constraints on the time dependence of the dark energy equation of state and the evolution of such supernovae in an epoch when the universe is matter dominated., Comment: Accepted for publication in the Astrophysical Journal Supplements, 22 pages, 16 figures. Updated Tables 3,4,8 and figures 6 and 8 to reflect replacement of Abell 963 with Abell 1423 in CLASH survey. A963 cannot be observed with WFC3 due to the lack of usable guide stars