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: Under review at the Astrophysical Journal, 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. Submitted to ApJ. Comments are welcome and appreciated

The gravitationally lensed quasar J014516.6-094517 at z=2.719 has been observed with the ESPRESSO instrument at the ESO VLT to obtain high-fidelity spectra of the two images A and B with a resolving power R=70000. At the redshifts under investigation (2.1 < z < 2.7), the Lyman forests along the two sightlines are separated by sub-kiloparsec physical distances and exhibit a strong correlation. We find that the two forests are indistinguishable at the present level of signal-to-noise ratio and do not show any global velocity shift, with the cross-correlation peaking at $\Delta v = 12 \pm 48$ m/s. The distribution of the difference in velocity of individual Lyman-$\alpha$ features is compatible with a null average and a mean absolute deviation of 930 m/s. Significant differences in NHI column density are not detected, putting a limit to the RMS fluctuation in the baryon density on $\leq 1$ proper kpc scales of $\Delta \rho / \rho < 3$%. On the other hand, metal lines show significant differences both in velocity structure and in column density. A toy model shows that the difference in velocity of the metal features between the two sightlines is compatible with the the motions of the baryonic component associated to dark matter halos of typical mass $M\simeq 2\times 10^{10} M_\odot$, also compatible with the observed incidence of the metal systems. The present observations confirm the feasibility of the Sandage test of the cosmic redshift drift with high-fidelity spectroscopy of the Lyman forest of distant, bright quasars, but also provide an element of caution about the intrinsic noise associated to the usage of metal features for the same purpose., Comment: 17 pages, 10 figures, 5 tables, accepted for publication in MNRAS

We present a new high-precision strong lensing model of PLCK G287.0$+$32.9, a massive lens galaxy cluster at $z=0.383$, with the aim to get an accurate estimation of its effective Einstein radius and total mass distribution. We also present a spectroscopic catalog containing accurate redshift measurements for 490 objects, including multiply-lensed sources and cluster member galaxies. We exploit high-quality spectroscopic data from three pointings of the VLT Multi Unit Spectroscopic Explorer, covering a central $3~\rm{arcmin}^2$ region of the cluster. We complete the spectroscopic catalog by including redshift measurements from VLT-VIMOS and KECK-DEIMOS. We identify 129 spectroscopic cluster member galaxies, with redshift values $0.360 \leq z \leq 0.405$ and $m_{\rm{F160W}} \leq 21$, and 24 photometric ones identified with a Convolutional Neural Network from ancillary HST imaging. We also identify 114 multiple images from 28 background sources, of which 84 images from 16 sources are new and the remaining ones were identified in previous work. The best-fitting lens model shows a root mean square separation value between the predicted and observed positions of the multiple images of $0.75''$, corresponding to an improvement in reconstructing the observed positions of the multiple images of a factor of $2.5$ with respect to previous models. Using the predictive power of our new lens model we find 3 new multiple images and we confirm the configuration of three systems of multiple images that were not used for the optimization of the model. The derived total mass distribution confirms this cluster to be a very prominent gravitational lens with an effective Einstein $\theta_{E} = 43.4'' \pm 0.1''$, that is in agreement with previous estimates and corresponds to a total mass enclosed in the critical curve of $M_E = {3.33}_{-0.07}^{+0.02} \times{ 10^{14} M_\odot}$., Comment: submitted to "Astronomy & Astrophysics" and accepted the on January 29th 2024

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

Previous studies have revealed that the estimated probability of galaxy-galaxy strong lensing in observed galaxy clusters exceeds the expectations from the $\Lambda$ Cold Dark Matter cosmological model by one order of magnitude. We aim to understand the origin of this excess by analyzing a larger set of simulated galaxy clusters and investigating how the theoretical expectations vary under different adopted prescriptions and numerical implementations of star formation and feedback in simulations. We perform a ray-tracing analysis of 324 galaxy clusters from the Three Hundred project, comparing the Gadget-X and Gizmo-Simba runs. These simulations, which start from the same initial conditions, are performed with different implementations of hydrodynamics and galaxy formation models tailored to match different observational properties of the Intra-Cluster-Medium and cluster galaxies. We find that galaxies in the Gizmo-Simba simulations develop denser stellar cores than their Gadget-X counterparts. Consequently, their probability for galaxy-galaxy strong lensing is higher by a factor of $\sim 3$. This increment is still insufficient to fill the gap with observations, as a discrepancy by a factor $\sim 4$ still persists. In addition, we find that several simulated galaxies have Einstein radii that are too large compared to observations. We conclude that a persistent excess of galaxy-galaxy strong lensing exists in observed galaxy clusters. The origin of this discrepancy with theoretical predictions is still unexplained in the framework of the cosmological hydrodynamical simulations. This might signal a hitherto unknown issue with either the simulation methods or our assumptions regarding the standard cosmological model., Comment: 6 pages, 2 figures, accepted for publication on A&A Letters. See companion paper Srivastava et al. (2023)

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

The detection of the accelerated expansion of the Universe has been one of the major breakthroughs in modern cosmology. Several cosmological probes (CMB, SNe Ia, BAO) have been studied in depth to better understand the nature of the mechanism driving this acceleration, and they are being currently pushed to their limits, obtaining remarkable constraints that allowed us to shape the standard cosmological model. In parallel to that, however, the percent precision achieved has recently revealed apparent tensions between measurements obtained from different methods. These are either indicating some unaccounted systematic effects, or are pointing toward new physics. Following the development of CMB, SNe, and BAO cosmology, it is critical to extend our selection of cosmological probes. Novel probes can be exploited to validate results, control or mitigate systematic effects, and, most importantly, to increase the accuracy and robustness of our results. This review is meant to provide a state-of-art benchmark of the latest advances in emerging beyond-standard cosmological probes. We present how several different methods can become a key resource for observational cosmology. In particular, we review cosmic chronometers, quasars, gamma-ray bursts, standard sirens, lensing time-delay with galaxies and clusters, cosmic voids, neutral hydrogen intensity mapping, surface brightness fluctuations, stellar ages of the oldest objects, secular redshift drift, and clustering of standard candles. The review describes the method, systematics, and results of each probe in a homogeneous way, giving the reader a clear picture of the available innovative methods that have been introduced in recent years and how to apply them. The review also discusses the potential synergies and complementarities between the various probes, exploring how they will contribute to the future of modern cosmology., Comment: Invited review article for Living Reviews in Relativity. Accepted version. 191 pages, 54 figures, 11 tables

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 generation of simulated convergence maps is of key importance in fully exploiting weak lensing by Large Scale Structure (LSS) from which cosmological parameters can be derived. In this paper we present an extension of the PINOCCHIO code which produces catalogues of dark matter haloes so that it is capable of simulating weak lensing by LSS. Like WL-MOKA, the method starts with a random realisation of cosmological initial conditions, creates a halo catalogue and projects it onto the past-light-cone, and paints in haloes assuming parametric models for the mass density distribution within them. Large scale modes that are not accounted for by the haloes are constructed using linear theory. We discuss the systematic errors affecting the convergence power spectra when Lagrangian Perturbation Theory at increasing order is used to displace the haloes within PINOCCHIO, and how they depend on the grid resolution. Our approximate method is shown to be very fast when compared to full ray-tracing simulations from an N-Body run and able to recover the weak lensing signal, at different redshifts, with a few percent accuracy. It also allows for quickly constructing weak lensing covariance matrices, complementing PINOCCHIO's ability of generating the cluster mass function and galaxy clustering covariances and thus paving the way for calculating cross covariances between the different probes. This work advances these approximate methods as tools for simulating and analysing surveys data for cosmological purposes., Comment: Replaced to match the accepted version for publication in MNRAS. Added: comparison of diagonal and off-diagonal terms of the covariance matrices; cosmological constraints. Acceptance date: 2020-06-02

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

We present a convolutional neural network to classify distinct cosmological scenarios based on the statistically similar weak-lensing maps they generate. Modified gravity (MG) models that include massive neutrinos can mimic the standard concordance model ($\Lambda$CDM) in terms of Gaussian weak-lensing observables. An inability to distinguish viable models that are based on different physics potentially limits a deeper understanding of the fundamental nature of cosmic acceleration. For a fixed redshift of sources, we demonstrate that a machine learning network trained on simulated convergence maps can discriminate between such models better than conventional higher-order statistics. Results improve further when multiple source redshifts are combined. To accelerate training, we implement a novel data compression strategy that incorporates our prior knowledge of the morphology of typical convergence map features. Our method fully distinguishes $\Lambda$CDM from its most similar MG model on noise-free data, and it correctly identifies among the MG models with at least 80% accuracy when using the full redshift information. Adding noise lowers the correct classification rate of all models, but the neural network still significantly outperforms the peak statistics used in a previous analysis., Comment: 12 pages, 10 figures, 1 table. Text updated and expanded throughout

Based on the DUSTGRAIN-pathfinder suite of simulations, we investigate observational degeneracies between nine models of modified gravity and massive neutrinos. Three types of machine learning techniques are tested for their ability to discriminate lensing convergence maps by extracting dimensional reduced representations of the data. Classical map descriptors such as the power spectrum, peak counts and Minkowski functionals are combined into a joint feature vector and compared to the descriptors and statistics that are common to the field of digital image processing. To learn new features directly from the data we use a Convolutional Neural Network (CNN). For the mapping between feature vectors and the predictions of their underlying model, we implement two different classifiers; one based on a nearest-neighbour search and one that is based on a fully connected neural network. We find that the neural network provides a much more robust classification than the nearest-neighbour approach and that the CNN provides the most discriminating representation of the data. It achieves the cleanest separation between the different models and the highest classification success rate of 59% for a single source redshift. Once we perform a tomographic CNN analysis, the total classification accuracy increases significantly to 76% with no observational degeneracies remaining. Visualising the filter responses of the CNN at different network depths provides us with the unique opportunity to learn from very complex models and to understand better why they perform so well., Comment: 21 pages, 14 figures, 10 tables. Associated code and data respository at https://www.bitbucket.org/jmerten82/mydnn . Accepted for publication by the MNRAS