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 structure formation models. 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) and, in most cases, dynamically relaxed. Five additional clusters are included for their lensing strength (θ Ein > 35 at z s = 2) to optimize the likelihood of finding highly magnified high-z (z > 7) galaxies. 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 (σ z ∼ 0.02(1+z)) photometric redshifts for newly discovered arcs. Observations of each cluster are spread over 8 epochs to enable a search 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 supernovae. We present newly re-derived X-ray luminosities, temperatures, and Fe abundances for the CLASH clusters as well as a representative source list for MACS1149.6+2223 (z = 0.544).

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 < ∼ z < ∼ 0.69 selected from the Cluster Lensing And Supernova survey with Hubble (CLASH). Our analysis uses wide-field multi-color imaging, taken primarily with Suprime-Cam on the Subaru Telescope. From a stacked shear-only analysis of the X-ray-selected subsample, we detect the ensemble-averaged lensing signal with a total signal-to-noise ratio of 25 in the radial range of 200 to 3500 kpc h −1 , providing integrated constraints on the halo profile shape and concentration-mass relation. The stacked tangential-shear signal is well described by a family of standard density profiles predicted for dark-matter-dominated halos in gravitational equilibrium, namely the Navarro-Frenk-White (NFW), truncated variants of NFW, and Einasto models. For the NFW model, we measure a mean concentration of c 200c = 4.01 +0.35 −0.32 at an effective halo mass of M 200c = 1.34 +0.10 −0.09 × 10 15 M . We show this is in excellent agreement with Λ cold-dark-matter (ΛCDM) predictions when the CLASH X-ray selection function and projection effects are taken into account. The best-fit Einasto shape parameter is α E = 0.191 +0.071 −0.068 , which is consistent with the NFWequivalent Einasto parameter of ∼ 0.18. We reconstruct projected mass density profiles of all CLASH clusters from a joint likelihood analysis of shear-and-magnification data, and measure cluster masses at several characteristic radii assuming an NFW density profile. We also derive an ensemble-averaged total projected mass profile of the X-ray-selected subsample by stacking their individual mass profiles. The stacked total mass profile, constrained by the shear+magnification data, is shown to be consistent with our shear-based halo-model predictions including the effects of surrounding large-scale structure as a two-halo term, establishing further consistency in the context of the ΛCDM model.

Abstract.We study the concentration parameters, their mass dependence and redshift evolution, of dark-matter halos in different dark-energy cosmologies with constant and time-variable equation of state, and compare them with "standard" ΛCDM and OCDM models. We find that previously proposed algorithms for predicting halo concentrations can be well adapted to dark-energy models. When centred on the analytically expected values, halo concentrations show a log-normal distribution with a uniform standard deviation of ∼0.2. The dependence of averaged halo concentrations on mass and redshift permits a simple fit of the form (1 + z) c = c 0 (M/M 0 ) α , with α ≈ −0.1 throughout. We find that the cluster concentration depends on the dark energy equation of state at the cluster formation redshift z coll through the linear growth factor D + (z coll ). As a simple correction accounting for dark-energy cosmologies, we propose scaling c 0 from ΛCDM with the ratio of linear growth factors, c 0 → c 0 D + (z coll )/D +, ΛCDM (z coll ).

We present results from a comprehensive lensing analysis in HST data, of the complete Cluster Lensing And Supernova survey with Hubble (CLASH) cluster sample. We identify new multiple-images previously undiscovered, allowing improved or first constraints on the cluster inner mass distributions and profiles. We combine these strong-lensing constraints with weak-lensing shape measurements within the HST FOV to jointly constrain the mass distributions. The analysis is performed in two different common parameterizations (one adopts light-traces-mass for both galaxies and dark matter while the other adopts an analytical, elliptical NFW form for the dark matter), to provide a better assessment of the underlying systematics -which is most important for deep, cluster-lensing surveys, especially when studying magnified high-redshift objects. We find that the typical (median), relative systematic differences throughout the central FOV are ∼ 40% in the (dimensionless) mass density, κ, and ∼ 20% in the magnification, µ. We show maps of these differences for each cluster, as well as the mass distributions, critical curves, and 2D integrated mass profiles. For the Einstein radii (z s = 2) we find that all typically agree within 10% between the two models, and Einstein masses agree, typically, within ∼ 15%. At larger radii, the total projected, 2D integrated mass profiles of the two models, within r ∼ 2 , differ by ∼ 30%. Stacking the surface-density profiles of the sample from the two methods together, we obtain an average slope of d log(Σ)/d log(r) ∼ −0.64 ± 0.1, in the radial range [5,350] kpc. Lastly, we also characterize the behavior of the average magnification, surface density, and shear differences between the two models, as a function of both the radius from the center, and the best-fit values of these quantities. All mass models and magnification maps are made publicly available for the community.

Context. Among the methods employed to measure the mass of galaxy clusters, the techniques based on lensing and X-ray analyses are perhaps the most widely used; however, the comparison between these mass estimates is often difficult and, in several clusters, the results apparently inconsistent. Aims. We aim at investigating potential biases in lensing and X-ray methods to measure the cluster mass profiles. Methods. We performed realistic simulations of lensing and X-ray observations that were subsequently analyzed using observational techniques. The resulting mass estimates were compared with the input models. Three clusters obtained from state-of-the-art hydrodynamical simulations, each of which projected along three independent lines-of-sight, were used for this analysis. Results. We find that strong lensing models can be trusted over a limited region around the cluster core. Extrapolating the strong lensing mass models to outside the Einstein ring can lead to significant biases in the mass estimates, if the BCG is not modeled properly, for example. Weak-lensing mass measurements can be strongly affected by substructures, depending on the method implemented to convert the shear into a mass estimate. Using nonparametric methods which combine weak and strong lensing data, the projected masses within R 200 can be constrained with a precision of ∼10%. Deprojection of lensing masses increases the scatter around the true masses by more than a factor of two because of cluster triaxiality. X-ray mass measurements have much smaller scatter (about a factor of two less than the lensing masses), but they are generally biased toward low values between 5 and 10%. This bias is entirely ascribable to bulk motions in the gas of our simulated clusters. Using the lensing and the X-ray masses as proxies for the true and the hydrostatic equilibrium masses of the simulated clusters and by averaging over the cluster sample, we are able to measure the lack of hydrostatic equilibrium in the systems we have investigated. Conclusions. Although the comparison between lensing and X-ray masses may be difficult in individual systems due to triaxiality and substructures, using a large number of clusters with both lensing and X-ray observations may lead to important information about their gas physics and allow use of lensing masses to calibrate the X-ray scaling relations.

We used near-infrared integral field spectroscopic observations from the AMAZE and LSD ESO programs to constrain the metallicity in a sample of 40 star-forming galaxies at 3 < z < 5 (most of which are at z ∼ 3.4). We measured metallicities by exploiting strong emission-line diagnostics. We found that a significant fraction of star-forming galaxies at z ∼ 3.4 deviate from the fundamental metallicity relation (FMR), with a metallicity of up to a factor of ten lower than expected according to the FMR. This deviation does not correlate with the dynamical properties of the galaxy or with the presence of interactions. To investigate the origin of the metallicity deviation in more detail, we also inferred information on the gas content by inverting the Schmidt-Kennicutt relation, assuming that the latter does not evolve out to z ∼ 3.4. In agreement with recent CO observational data, we found that in contrast with the steeply rising trend at 0 < z < 2, the gas fraction in massive galaxies remains constant, with an indication of a marginal decline at 2 < z < 3.5. When combined with the metallicity information, we infer that to explain the low metallicity and gas content in z ∼ 3.4 galaxies, both prominent outflows and massive pristine gas inflows are needed. In ten galaxies we can also spatially resolve the metallicity distribution. We found that the metallicity generally anticorrelates with the distribution of star formation and with the gas surface density. We discuss these findings in terms of pristine gas inflows toward the center, and outflows of metal-rich gas from the center toward the external regions.

The early Universe at redshift z∼6-11 marks the reionization of the intergalactic medium, following the formation of the first generation of stars. However, those young galaxies at a cosmic age of < ∼ 500 million years (Myr, at z > ∼ 10) remain largely unexplored as they are at or beyond the sensitivity limits of current large telescopes. Gravitational lensing by galaxy clusters enables the detection of high-redshift galaxies that are fainter than what otherwise could be found in the deepest images of the sky. We report the discovery of an object found in the multi-band observations of the cluster MACS1149+22 that has a high probability of being a gravitationally magnified object from the early universe. The object is firmly detected (12σ) in the two reddest bands of HST/WFC3, and not detected below 1.2 µm, matching the characteristics of z∼9 objects. We derive a robust photometric redshift of z = 9.6 ± 0.2, corresponding to a cosmic age of 490 ± 15 Myr (i.e., 3.6% of the age of the Universe).The large number of bands used to derive the redshift estimate make it one of the most accurate estimates ever obtained for such a distant object. The significant magnification by cluster lensing (a factor of ∼15) allows us to analyze the object's ultra-violet and optical luminosity in its rest-2 frame, thus enabling us to constrain on its stellar mass, star-formation rate and age. If the galaxy is indeed at such a large redshift, then its age is less than 200 Myr (at the 95% confidence level), implying a formation redshift of z f < ∼ 14. The object is the first z>9 candidate that is bright enough for detailed spectroscopic studies with JWST, demonstrating the unique potential of galaxy cluster fields for finding highly magnified, intrinsically faint galaxies at the highest redshifts.Observational cosmology has established that the age of the Universe is 13.7 billion years, and the reionization of the vast intergalactic medium (IGM) started around redshift z ∼ 11, 1 as the result of radiation from the first generation of stars. The task of probing the most distant Universe is progressively challenging: While more than 10 5 quasars have been found, only one is at z > 7; 2 while thousands of gamma-ray burst events have been recorded, only one 3 is confirmed at z=8.3; and while thousands of galaxy candidates have been found at z ∼ 6, only one has been reported at z ∼ 10, 4 which is based on a single-band detection. Galaxies at z ∼ 10 are expected to be at a magnitude of ∼ 29 (in the AB system, used hereafter) 4, 5 , near the detection limits of the deepest fields observed by Hubble Space Telescope (HST), and beyond the spectroscopic capability of even the next generation of large telescopes.In this Letter we report the discovery of a gravitationally lensed source whose most likely redshift is z ∼ 9.6. The source, hereafter called MACS1149-JD1, is selected from a near-infrared detection image at significance of 22σ. MACS1149-JD1 has a unique flux distribution characterized by a) no detection at Galaxy clusters are the largest r...

We present the results of a numerical study based on the analysis of the MUSIC-2 N-body/hydrodynamical simulations, aimed at estimating the expected concentration-mass relation for the CLASH cluster sample. We study nearly 1400 halos simulated at high spatial and mass resolution, which were projected along many linesof-sight each. We study the shape of both their density and surface-density profiles and fit them with a variety of radial functions, including the Navarro-Frenk-White, the generalised Navarro-Frenk-White, and the Einasto density profiles. We derive concentrations and masses from these fits and investigate their distributions as a function of redshift and halo relaxation. We use the X-ray image simulator X-MAS to produce simulated Chandra observations of the halos and we use them to identify objects resembling the X-ray morphologies and masses of the clusters in the CLASH X-ray selected sample. We also derive a concentration-mass relation for strong-lensing clusters. We find that the sample of simulated halos which resemble the X-ray morphology of the CLASH clusters is composed mainly by relaxed halos, but it also contains a significant fraction of unrelaxed systems. For such a heterogeneous sample we measure an average 2D concentration which is ∼ 11% higher than found for the full sample of simulated halos. After accounting for projection and selection effects, the average NFW concentrations of CLASH clusters are expected to be intermediate between those predicted in 3D for relaxed and super-relaxed halos. Matching the simulations to the individual CLASH clusters on the basis of the X-ray morphology, we expect that the NFW concentrations recovered from the lensing analysis of the CLASH clusters are in the range [3 − 6], with an average value of 3.87 and a standard deviation of 0.61. Simulated halos with X-ray morphologies similar to those of the CLASH clusters are affected by a modest orientation bias.