We present a statistical weak-lensing magnification analysis on an optically selected sample of 3029 \texttt{CAMIRA} galaxy clusters with richness $N>15$ at redshift $0.2\leq z <1.1$ in the Subaru Hyper Suprime-Cam (HSC) survey. We use two distinct populations of color-selected, flux-limited background galaxies, namely the low-$z$ and high-$z$ samples at mean redshifts of $\approx1.1$ and $\approx1.4$, respectively, from which to measure the weak-lensing magnification signal by accounting for cluster contamination as well as masking effects. Our magnification bias measurements are found to be uncontaminated according to validation tests against the "null-test" samples for which the net magnification bias is expected to vanish. The magnification bias for the full \texttt{CAMIRA} sample is detected at a significance level of $9.51\sigma$, which is dominated by the high-$z$ background. We forward-model the observed magnification data to constrain the normalization of the richness-to-mass ($N$--$M$) relation for the \texttt{CAMIRA} sample with informative priors on other parameters. The resulting scaling relation is $N\propto {M_{500}}^{0.92\pm0.13} (1 + z)^{-0.48\pm0.69}$, with a characteristic richness of $N=\left(17.72\pm2.60\right)$ and intrinsic log-normal scatter of $0.15\pm0.07$ at $M_{500} = 10^{14}h^{-1}M_{\odot}$. With the derived $N$--$M$ relation, we provide magnification-calibrated mass estimates of individual \texttt{CAMIRA} clusters, with the typical uncertainty of $\approx39\%$ and $\approx32\%$ at richness$\approx20$ and $\approx40$, respectively. We further compare our magnification-inferred $N$--$M$ relation with those from the shear-based results in the literature, finding good agreement.

We present a weak-lensing analysis of X-ray galaxy groups and clusters selected from the XMM-XXL survey using the first-year data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program. Our joint weak-lensing and X-ray analysis focuses on 136 spectroscopically confirmed X-ray-selected systems at 0.031 < z < 1.033 detected in the 25sqdeg XXL-N region. We characterize the mass distributions of individual clusters and establish the concentration-mass (c-M) relation for the XXL sample, by accounting for selection bias and statistical effects, and marginalizing over the remaining mass calibration uncertainty. We find the mass-trend parameter of the c-M relation to be \beta = -0.07 \pm 0.28 and the normalization to be c200 = 4.8 \pm 1.0 (stat) \pm 0.8 (syst) at M200=10^{14}Msun/h and z = 0.3. We find no statistical evidence for redshift evolution. Our weak-lensing results are in excellent agreement with dark-matter-only c-M relations calibrated for recent LCDM cosmologies. The level of intrinsic scatter in c200 is constrained as \sigma(\ln[c200]) < 24% (99.7% CL), which is smaller than predicted for the full population of LCDM halos. This is likely caused in part by the X-ray selection bias in terms of the relaxation state. We determine the temperature-mass (Tx-M500) relation for a subset of 105 XXL clusters that have both measured HSC lensing masses and X-ray temperatures. The resulting Tx-M500 relation is consistent with the self-similar prediction. Our Tx-M500 relation agrees with the XXL DR1 results at group scales, but has a slightly steeper mass trend, implying a smaller mass scale in the cluster regime. The overall offset in the Tx-M500 relation is at the $1.5\sigma$ level, corresponding to a mean mass offset of (34\pm 20)%. We also provide bias-corrected, weak-lensing-calibrated M200 and M500 mass estimates of individual XXL clusters based on their measured X-ray temperatures.

We present results from simultaneous modeling of high angular resolution GBT/MUSTANG-2 90 GHz Sunyaev-Zel'dovich effect (SZE) measurements and XMM-XXL X-ray images of three rich galaxy clusters selected from the HSC-SSP Survey. The combination of high angular resolution SZE and X-ray imaging enables a spatially resolved multi-component analysis, which is crucial to understand complex distributions of cluster gas properties. The targeted clusters have similar optical richnesses and redshifts, but exhibit different dynamical states in their member galaxy distributions; a single-peaked cluster, a double-peaked cluster, and a cluster belonging to a supercluster. A large-scale residual pattern in both regular Compton-parameter $y$ and X-ray surface brightness distributions is found in the single-peaked cluster, indicating a sloshing mode. The double-peaked cluster shows an X-ray remnant cool core between two SZE peaks associated with galaxy concentrations. The temperatures of the two peaks reach $\sim20-30$ keV in contrast to the cool core component of $\sim2$ keV, indicating a violent merger. The highest SZE signal for the supercluster is elongated along a direction perpendicular to the major axis of the X-ray core, suggesting a minor merger before core passage. The $S_X$ and $y$ distributions are thus perturbed at some level, regardless of the optical properties. We find that an integrated Compton $y$ parameter and a temperature for the major merger are significantly boosted from those expected by the weak-lensing mass and those for the other two clusters show no significant deviations, which is consistent with predictions of numerical simulations.

Abell 2744, a massive Hubble Frontier Fields merging galaxy cluster with many multiple images in the core has been the subject of many lens inversions using different methods. While most existing studies compare various inversion methods, we focus on a comparison of reconstructions that use different input lensing data. Since the quantity and quality of lensing data is constantly improving, it makes sense to ask if the estimated uncertainties are robust against changes in the data. We address this question using free-form Grale, which takes only image information as input, and nothing pertaining to cluster galaxies. We reconstruct Abell 2744 using two sets of strong lensing data from the Hubble Frontier Fields community. Our first and second reconstructions use 55 and 91 images, respectively, and only 10 of the 91 images have the same positions and redshifts as in the first reconstruction. Comparison of the two mass maps shows that Grale uncertainties are robust against these changes, as well as small modifications in the inversion routine. Additionally, applying the methods used in Sebesta et al. (2016) for MACS J0416, we conclude that, in a statistical sense, light follows mass in Abell 2744, with brighter galaxies clustering stronger with the recovered mass than the fainter ones. We also show that the faintest galaxies are anti-correlated with mass, which is likely the result of light contamination from bright galaxies, and lensing magnification bias acting on galaxies background to the cluster.

We present multi-wavelength studies of optically-defined merging clusters, based on the Hyper Suprime-Cam Subaru Strategic Program. Luminous red galaxies, tracing cluster mass distributions, enable to identify cluster subhalos at various merging stages, and thus make a homogeneous sample of cluster mergers, which is unbiased with respect to the merger boost of the intracluster medium (ICM). We define, using a peak-finding method, merging clusters with multiple-peaks and single clusters with single-peaks from the CAMIRA cluster catalog. Stacked weak-lensing analysis indicates that our sample of the merging clusters is categorized into major mergers. The average halo concentration for the merging clusters is $\sim70\%$ smaller than that of the single-peak clusters, which agrees well with predictions of numerical simulations. The spatial distribution of subhalos is less centrally concentrated than the mass distribution of the main halo. The fractions of red galaxies in the merging clusters are not higher than those of the single-peak clusters. We find a signature of the merger boost of the ICM from stacked Planck Sunyaev-Zeldovich effect and ROSAT X-ray luminosity, but not in optical richness. The stacked X-ray surface brightness distribution, aligned with the main-subhalo pairs of low redshift and massive clusters, shows that the central gas core is elongated along the merger axis and overall gas distribution is misaligned by $\sim60$ deg. The homogeneous, unbiased sample of cluster mergers and multi-wavelength follow-up studies provide a unique opportunity to make a complete picture of merger physics over the whole process.

Community access to deep (i ~ 25), highly-multiplexed optical and near-infrared multi-object spectroscopy (MOS) on 8-40m telescopes would greatly improve measurements of cosmological parameters from LSST. The largest gain would come from improvements to LSST photometric redshifts, which are employed directly or indirectly for every major LSST cosmological probe; deep spectroscopic datasets will enable reduced uncertainties in the redshifts of individual objects via optimized training. Such spectroscopy will also determine the relationship of galaxy SEDs to their environments, key observables for studies of galaxy evolution. The resulting data will also constrain the impact of blending on photo-z's. Focused spectroscopic campaigns can also improve weak lensing cosmology by constraining the intrinsic alignments between the orientations of galaxies. Galaxy cluster studies can be enhanced by measuring motions of galaxies in and around clusters and by testing photo-z performance in regions of high density. Photometric redshift and intrinsic alignment studies are best-suited to instruments on large-aperture telescopes with wider fields of view (e.g., Subaru/PFS, MSE, or GMT/MANIFEST) but cluster investigations can be pursued with smaller-field instruments (e.g., Gemini/GMOS, Keck/DEIMOS, or TMT/WFOS), so deep MOS work can be distributed amongst a variety of telescopes. However, community access to large amounts of nights for surveys will still be needed to accomplish this work. In two companion white papers we present gains from shallower, wide-area MOS and from single-target imaging and spectroscopy.

Do cool-core (CC) and noncool-core (NCC) clusters live in different environments? We make novel use of H$\alpha$ emission lines in the central galaxies of redMaPPer clusters as proxies to construct large (1,000's) samples of CC and NCC clusters, and measure their relative assembly bias using both clustering and weak lensing. We increase the statistical significance of the bias measurements from clustering by cross-correlating the clusters with an external galaxy redshift catalog from the Sloan Digital Sky Survey III, the LOWZ sample. Our cross-correlations can constrain assembly bias up to a statistical uncertainty of 6%. Given our H$\alpha$ criteria for CC and NCC, we find no significant differences in their clustering amplitude. Interpreting this difference as the absence of halo assembly bias, our results rule out the possibility of having different large-scale (tens of Mpc) environments as the source of diversity observed in cluster cores. Combined with recent observations of the overall mild evolution of CC and NCC properties, such as central density and CC fraction, this would suggest that either the cooling properties of the cluster core are determined early on solely by the local (<200 kpc) gas properties at formation or that local merging leads to stochastic CC relaxation and disruption in a periodic way, preserving the average population properties over time. Studying the small-scale clustering in clusters at high redshift would help shed light on the exact scenario.

We measure cosmic weak lensing shear power spectra with the Subaru Hyper Suprime-Cam (HSC) survey first-year shear catalog covering 137deg$^2$ of the sky. Thanks to the high effective galaxy number density of $\sim$17 arcmin$^{-2}$ even after conservative cuts such as magnitude cut of $i<24.5$ and photometric redshift cut of $0.3\leq z \leq 1.5$, we obtain a high significance measurement of the cosmic shear power spectra in 4 tomographic redshift bins, achieving a total signal-to-noise ratio of 16 in the multipole range $300 \leq \ell \leq 1900$. We carefully account for various uncertainties in our analysis including the intrinsic alignment of galaxies, scatters and biases in photometric redshifts, residual uncertainties in the shear measurement, and modeling of the matter power spectrum. The accuracy of our power spectrum measurement method as well as our analytic model of the covariance matrix are tested against realistic mock shear catalogs. For a flat $\Lambda$ cold dark matter ($\Lambda$CDM) model, we find $S_8\equiv \sigma_8(\Omega_{\rm m}/0.3)^\alpha=0.800^{+0.029}_{-0.028}$ for $\alpha=0.45$ ($S_8=0.780^{+0.030}_{-0.033}$ for $\alpha=0.5$) from our HSC tomographic cosmic shear analysis alone. In comparison with Planck cosmic microwave background constraints, our results prefer slightly lower values of $S_8$, although metrics such as the Bayesian evidence ratio test do not show significant evidence for discordance between these results. We study the effect of possible additional systematic errors that are unaccounted in our fiducial cosmic shear analysis, and find that they can shift the best-fit values of $S_8$ by up to $\sim 0.6\sigma$ in both directions. The full HSC survey data will contain several times more area, and will lead to significantly improved cosmological constraints.

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.

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

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.

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.

We present weak-lensing measurements using the first-year data from the Hyper Suprime-Cam Strategic Survey Program on the Subaru telescope for eight galaxy clusters selected through their thermal Sunyaev-Zel'dovich (SZ) signal measured at 148 GHz with the Atacama Cosmology Telescope Polarimeter experiment. The overlap between the two surveys in this work is 33.8 square degrees, before masking bright stars. The signal-to-noise ratio of individual cluster lensing measurements ranges from 2.2 to 8.7, with a total of 11.1 for the stacked cluster weak-lensing signal. We fit for an average weak-lensing mass distribution using three different profiles, a Navarro-Frenk-White profile, a dark-matter-only emulated profile, and a full cosmological hydrodynamic emulated profile. We interpret the differences among the masses inferred by these models as a systematic error of 10\%, which is currently smaller than the statistical error. We obtain the ratio of the SZ-estimated mass to the lensing-estimated mass (the so-called hydrostatic mass bias $1-b$) of $0.74^{+0.13}_{-0.12}$, which is comparable to previous SZ-selected clusters from the Atacama Cosmology Telescope and from the {\sl Planck} Satellite. We conclude with a discussion of the implications for cosmological parameters inferred from cluster abundances compared to cosmic microwave background primary anisotropy measurements.

Hyper Suprime-Cam (HSC) is a wide-field imaging camera on the prime focus of the 8.2m Subaru telescope on the summit of Maunakea in Hawaii. A team of scientists from Japan, Taiwan and Princeton University is using HSC to carry out a 300-night multi-band imaging survey of the high-latitude sky. The survey includes three layers: the Wide layer will cover 1400 deg$^2$ in five broad bands ($grizy$), with a $5\,\sigma$ point-source depth of $r \approx 26$. The Deep layer covers a total of 26~deg$^2$ in four fields, going roughly a magnitude fainter, while the UltraDeep layer goes almost a magnitude fainter still in two pointings of HSC (a total of 3.5 deg$^2$). Here we describe the instrument, the science goals of the survey, and the survey strategy and data processing. This paper serves as an introduction to a special issue of the Publications of the Astronomical Society of Japan, which includes a large number of technical and scientific papers describing results from the early phases of this survey.

We have selected a sample of eleven massive clusters of galaxies observed by the Hubble Space Telescope in order to study the impact of the dynamical state on the IntraCluster Light (ICL) fraction, the ratio of total integrated ICL to the total galaxy member light. With the exception of the Bullet cluster, the sample is drawn from the Cluster Lensing and Supernova Survey and the Frontier Fields program, containing five relaxed and six merging clusters. The ICL fraction is calculated in three optical filters using the CHEFs IntraCluster Light Estimator, a robust and accurate algorithm free of a priori assumptions. We find that the ICL fraction in the three bands is, on average, higher for the merging clusters, ranging between $\sim7-23\%$, compared with the $\sim 2-11\%$ found for the relaxed systems. We observe a nearly constant value (within the error bars) in the ICL fraction of the regular clusters at the three wavelengths considered, which would indicate that the colors of the ICL and the cluster galaxies are, on average, coincident and, thus, their stellar populations. However, we find a higher ICL fraction in the F606W filter for the merging clusters, consistent with an excess of lower-metallicity/younger stars in the ICL, which could have migrated violently from the outskirts of the infalling galaxies during the merger event.

We present the result of searching for clusters of galaxies based on weak gravitational lensing analysis of the $\sim 160$~deg$^2$ area surveyed by Hyper Suprime-Cam (HSC) as a Subaru Strategic Program. HSC is a new prime focus optical imager with a 1.5 diameter field of view on the 8.2-meter Subaru telescope. The superb median seeing on the HSC $i$-band images of $0.56$ arcsec allows the reconstruction of high angular resolution mass maps via weak lensing, which is crucial for the weak lensing cluster search. We identify 65 mass map peaks with signal-to-noise (SN) ratio larger than 4.7, and carefully examine their properties by cross-matching the clusters with optical and X-ray cluster catalogs. We find that all the 39 peaks with SN$>5.1$ have counterparts in the optical cluster catalogs, and only 2 out of the 65 peaks are probably false positives. The upper limits of X-ray luminosities from ROSAT All Sky Survey (RASS) imply the existence of an X-ray under-luminous cluster population. We show that the X-rays from the shear selected clusters can be statistically detected by stacking the RASS images. The inferred average X-ray luminosity is about half that of the X-ray selected clusters of the same mass. The radial profile of the dark matter distribution derived from the stacking analysis is well modeled by the Navarro-Frenk-White profile with a small concentration parameter value of $c_{500}\sim 2.5$, which suggests that the selection bias on the orientation or the internal structure for our shear selected cluster sample is not strong.

Using $\sim$140 deg$^2$ Subaru Hyper Suprime-Cam (HSC) survey data, we stack the weak lensing (WL) signal around five Planck clusters found within the footprint. This yields a 15$\sigma$ detection of the mean Planck cluster mass density profile. The five Planck clusters span a relatively wide mass range, $M_{\rm WL,500c} = (2-30)\times10^{14}\,M_\odot/h$ with a mean mass of $M_{\rm WL,500c} = (4.15\pm0.61)\times10^{14}\,M_\odot/h$. The ratio of the stacked Planck Sunyaev-Zel'dovich (SZ) mass to the stacked WL mass is $ \langle M_{\rm SZ}\rangle/\langle M_{\rm WL}\rangle = 1-b = 0.80\pm0.14$. This mass bias is consistent with previous WL mass calibrations of Planck clusters within the errors. We discuss the implications of our findings for the calibration of SZ cluster counts and the much discussed tension between Planck SZ cluster counts and Planck $\Lambda$CDM cosmology.

We present optimized source galaxy selection schemes for measuring cluster weak lensing (WL) mass profiles unaffected by cluster member dilution from the Subaru Hyper Suprime-Cam Strategic Survey Program (HSC-SSP). The ongoing HSC-SSP survey will uncover thousands of galaxy clusters to $z\lesssim1.5$. In deriving cluster masses via WL, a critical source of systematics is contamination and dilution of the lensing signal by cluster {members, and by foreground galaxies whose photometric redshifts are biased}. Using the first-year CAMIRA catalog of $\sim$900 clusters with richness larger than 20 found in $\sim$140 deg$^2$ of HSC-SSP data, we devise and compare several source selection methods, including selection in color-color space (CC-cut), and selection of robust photometric redshifts by applying constraints on their cumulative probability distribution function (PDF; P-cut). We examine the dependence of the contamination on the chosen limits adopted for each method. Using the proper limits, these methods give mass profiles with minimal dilution in agreement with one another. We find that not adopting either the CC-cut or P-cut methods results in an underestimation of the total cluster mass ($13\pm4\%$) and the concentration of the profile ($24\pm11\%$). The level of cluster contamination can reach as high as $\sim10\%$ at $R\approx 0.24$ Mpc/$h$ for low-z clusters without cuts, while employing either the P-cut or CC-cut results in cluster contamination consistent with zero to within the 0.5% uncertainties. Our robust methods yield a $\sim60\sigma$ detection of the stacked CAMIRA surface mass density profile, with a mean mass of $M_\mathrm{200c} = (1.67\pm0.05({\rm {stat}}))\times 10^{14}\,M_\odot/h$.

We present and characterize the catalog of galaxy shape measurements that will be used for cosmological weak lensing measurements in the Wide layer of the first year of the Hyper Suprime-Cam (HSC) survey. The catalog covers an area of 136.9 deg$^2$ split into six fields, with a mean $i$-band seeing of $0.58$ arcsec and $5\sigma$ point-source depth of $i\sim 26$. Given conservative galaxy selection criteria for first year science, the depth and excellent image quality results in unweighted and weighted source number densities of 24.6 and 21.8 arcmin$^{-2}$, respectively. We define the requirements for cosmological weak lensing science with this catalog, then focus on characterizing potential systematics in the catalog using a series of internal null tests for problems with point-spread function (PSF) modeling, shear estimation, and other aspects of the image processing. We find that the PSF models narrowly meet requirements for weak lensing science with this catalog, with fractional PSF model size residuals of approximately $0.003$ (requirement: 0.004) and the PSF model shape correlation function $\rho_1<3\times 10^{-7}$ (requirement: $4\times 10^{-7}$) at 0.5$^\circ$ scales. A variety of galaxy shape-related null tests are statistically consistent with zero, but star-galaxy shape correlations reveal additive systematics on $>1^\circ$ scales that are sufficiently large as to require mitigation in cosmic shear measurements. Finally, we discuss the dominant systematics and the planned algorithmic changes to reduce them in future data reductions.

Collisions and interactions between gas-rich galaxies are thought to be pivotal stages in their formation and evolution, causing the rapid production of new stars, and possibly serving as a mechanism for fueling supermassive black holes (BH). Harnessing the exquisite spatial resolution (~0.5 arcsec) afforded by the first ~170 deg^2 of the Hyper Suprime-Cam (HSC) Survey, we present our new constraints on the importance of galaxy-galaxy major mergers (1:4) in growing BHs throughout the last ~8 Gyrs. Utilizing mid-infrared observations in the WISE All-Sky survey, we robustly select active galactic nuclei (AGN) and mass-matched control galaxy samples, totaling ~140,000 spectroscopically confirmed systems at i<22 mag. We identify galaxy interaction signatures using a novel machine-learning random forest decision tree technique allowing us to select statistically significant samples of major-mergers, minor-mergers/irregular-systems, and non-interacting galaxies. We use these samples to show that galaxies undergoing mergers are a factor ~2-7 more likely to contain luminous obscured AGN than non-interacting galaxies, and this is independent of both stellar mass and redshift to z < 0.9. Furthermore, based on our comparison of AGN fractions in mass-matched samples, we determine that the most luminous AGN population (L_AGN > 10^45 erg/s) systematically reside in merging systems over non-interacting galaxies. Our findings show that galaxy-galaxy interactions do, on average, trigger luminous AGN activity substantially more often than in secularly evolving non-interacting galaxies, and we further suggest that the BH growth rate may be closely tied to the dynamical time of the merger system.

We present a joint X-ray, optical and weak-lensing analysis for X-ray luminous galaxy clusters selected from the MCXC (Meta-Catalog of X-Ray Detected Clusters of Galaxies) cluster catalog in the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) survey field with S16A data, As a pilot study of our planned series papers, we measure hydrostatic equilibrium (H.E.) masses using XMM-Newton data for four clusters in the current coverage area out of a sample of 22 MCXC clusters. We additionally analyze a non-MCXC cluster associated with one MCXC cluster. We show that H.E. masses for the MCXC clusters are correlated with cluster richness from the CAMIRA catalog (Oguri et al. 2017), while that for the non-MCXC cluster deviates from the scaling relation. The mass normalization of the relationship between the cluster richness and H.E. mass is compatible with one inferred by matching CAMIRA cluster abundance with a theoretical halo mass function. The mean gas mass fraction based on H.E. masses for the MCXC clusters is $\langle f_{\rm gas} \rangle = 0.125\pm0.012$ at spherical overdensity $\Delta=500$, which is $\sim80-90$ percent of the cosmic mean baryon fraction, $\Omega_b/\Omega_m$, measured by cosmic microwave background experiments. We find that the mean baryon fraction estimated from X-ray and HSC-SSP optical data is comparable to $\Omega_b/\Omega_m$. A weak-lensing shear catalog of background galaxies, combined with photometric redshifts, is currently available only for three clusters in our sample. Hydrostatic equilibrium masses roughly agree with weak-lensing masses, albeit with large uncertainty. This study demonstrates that further multiwavelength study for a large sample of clusters using X-ray, HSC-SSP optical and weak lensing data will enable us to understand cluster physics and utilize cluster-based cosmology.

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

The Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) is a three-layered imaging survey aimed at addressing some of the most outstanding questions in astronomy today, including the nature of dark matter and dark energy. The survey has been awarded 300 nights of observing time at the Subaru Telescope and it started in March 2014. This paper presents the first public data release of HSC-SSP. This release includes data taken in the first 1.7 years of observations (61.5 nights) and each of the Wide, Deep, and UltraDeep layers covers about 108, 26, and 4 square degrees down to depths of i~26.4, ~26.5, and ~27.0 mag, respectively (5sigma for point sources). All the layers are observed in five broad bands (grizy), and the Deep and UltraDeep layers are observed in narrow bands as well. We achieve an impressive image quality of 0.6 arcsec in the i-band in the Wide layer. We show that we achieve 1-2 per cent PSF photometry (rms) both internally and externally (against Pan-STARRS1), and ~10 mas and 40 mas internal and external astrometric accuracy, respectively. Both the calibrated images and catalogs are made available to the community through dedicated user interfaces and database servers. In addition to the pipeline products, we also provide value-added products such as photometric redshifts and a collection of public spectroscopic redshifts. Detailed descriptions of all the data can be found online. The data release website is https://hsc-release.mtk.nao.ac.jp/.

We present accurate photometric redshifts for galaxies observed by the Cluster Lensing and Supernova survey with Hubble (CLASH). CLASH observed 25 massive galaxy cluster cores with the Hubble Space Telescope in 16 filters spanning 0.2 - 1.7 $\mu$m. Photometry in such crowded fields is challenging. Compared to our previously released catalogs, we make several improvements to the photometry, including smaller apertures, ICL subtraction, PSF matching, and empirically measured uncertainties. We further improve the Bayesian Photometric Redshift (BPZ) estimates by adding a redder elliptical template and by inflating the photometric uncertainties of the brightest galaxies. The resulting photometric redshift accuracies are dz/(1+z) $\sim$ 0.8\%, 1.0\%, and 2.0\% for galaxies with I-band F814W AB magnitudes $<$ 18, 20, and 23, respectively. These results are consistent with our expectations. They improve on our previously reported accuracies by a factor of 4 at the bright end and a factor of 2 at the faint end. Our new catalog includes 1257 spectroscopic redshifts, including 382 confirmed cluster members. We also provide stellar mass estimates. Finally, we include lensing magnification estimates of background galaxies based on our public lens models. Our new catalog of all 25 CLASH clusters is available via MAST. The analysis techniques developed here will be useful in other surveys of crowded fields, including the Frontier Fields and surveys carried out with J-PAS and JWST.

We present an optically-selected cluster catalog from the Hyper Suprime-Cam (HSC) Subaru Strategic Program. The HSC images are sufficiently deep to detect cluster member galaxies down to $M_*\sim 10^{10.2}M_\odot$ even at $z\sim 1$, allowing a reliable cluster detection at such high redshifts. We apply the CAMIRA algorithm to the HSC Wide S16A dataset covering $\sim 232$ deg$^2$ to construct a catalog of 1921 clusters at redshift $0.1<z<1.1$ and richness $\hat{N}_{\rm mem}>15$ that roughly corresponds to $M_{\rm 200m}\gtrsim 10^{14}h^{-1}M_\odot$. We confirm good cluster photometric redshift performance, with the bias and scatter in $\Delta z/(1+z)$ being better than 0.005 and 0.01 over most of the redshift range, respectively. We compare our cluster catalog with large X-ray cluster catalogs from XXL and XMM-LSS surveys and find good correlation between richness and X-ray properties. We also study the miscentering effect from the distribution of offsets between optical and X-ray cluster centers. We confirm the high ($>0.9$) completeness and purity for high mass clusters by analyzing mock galaxy catalogs.

We present the procedure to build and validate the bright-star masks for the Hyper-Suprime-Cam Strategic Subaru Proposal (HSC-SSP) survey. To identify and mask the saturated stars in the full HSC-SSP footprint, we rely on the Gaia and Tycho-2 star catalogues. We first assemble a pure star catalogue down to $G_{\rm Gaia} < 18$ after removing $\sim1.5\%$ of sources that appear extended in the Sloan Digital Sky Survey (SDSS). We perform visual inspection on the early data from the S16A internal release of HSC-SSP, finding that our star catalogue is $99.2\%$ pure down to $G_{\rm Gaia} < 18$. Second, we build the mask regions in an automated way using stacked detected source measurements around bright stars binned per $G_{\rm Gaia}$ magnitude. Finally, we validate those masks from visual inspection and comparison with the literature of galaxy number counts and angular two-point correlation functions. This version (Arcturus) supersedes the previous version (Sirius) used in the S16A internal and DR1 public releases. We publicly release the full masks and tools to flag objects in the entire footprint of the planned HSC-SSP observations at this address: ftp://obsftp.unige.ch/pub/coupon/brightStarMasks/HSC-SSP/.

There are well-observed differences between cool-core (CC) and non-cool-core (NCC) clusters, but the origin of this distinction is still largely unknown. Competing theories can be divided into internal (inside-out), in which internal physical processes transform or maintain the NCC phase, and external (outside-in), in which the cluster type is determined by its initial conditions, which in turn lead to different formation histories (i.e., assembly bias). We propose a new method that uses the relative assembly bias of CC to NCC clusters, as determined via the two-point cluster-galaxy cross-correlation function (CCF), to test whether formation history plays a role in determining their nature. We apply our method to 48 ACCEPT clusters, which have well resolved central entropies, and cross-correlate with the SDSS-III/BOSS LOWZ galaxy catalog. We find that the relative bias of NCC over CC clusters is $b = 1.42 \pm 0.35$ ($1.6\sigma$ different from unity). Our measurement is limited by the small number of clusters with core entropy information within the BOSS footprint, 14 CC and 34 NCC. Future compilations of X-ray cluster samples, combined with deep all-sky redshift surveys, will be able to better constrain the relative assembly bias of CC and NCC clusters and determine the origin of the bimodality.

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

We present VIMOS-VLT spectroscopy of the Frontier Fields cluster MACS~J0416.1-2403 (z=0.397). Taken as part of the CLASH-VLT survey, the large spectroscopic campaign provided more than 4000 reliable redshifts, including ~800 cluster member galaxies. The unprecedented sample of cluster members at this redshift allows us to perform a highly detailed dynamical and structural analysis of the cluster out to ~2.2$r_{200}$ (~4Mpc). Our analysis of substructures reveals a complex system composed of a main massive cluster ($M_{200}$~0.9$\times 10^{15} M_{\odot}$) presenting two major features: i) a bimodal velocity distribution, showing two central peaks separated by $\Delta V_{rf}$~1100 km s$^{-1}$ with comparable galaxy content and velocity dispersion, ii) a projected elongation of the main substructures along the NE-SW direction, with a prominent subclump ~600 kpc SW of the center and an isolated BCG approximately halfway between the center and the SW clump. We also detect a low mass structure at z~0.390, ~10' S of the cluster center, projected at ~3Mpc, with a relative line-of-sight velocity of $\Delta V_{rf}$~-1700 km s$^{-1}$. The cluster mass profile that we obtain through our dynamical analysis deviates significantly from the "universal" NFW, being best fit by a Softened Isothermal Sphere model instead. The mass profile measured from the galaxy dynamics is found to be in relatively good agreement with those obtained from strong and weak lensing, as well as with that from the X-rays, despite the clearly unrelaxed nature of the cluster. Our results reveal overall a complex dynamical state of this massive cluster and support the hypothesis that the two main subclusters are being observed in a pre-collisional phase, in line with recent findings from radio and deep X-ray data. With this article we also release the entire redshift catalog of 4386 sources in the field of this cluster.

We present a weak-lensing analysis of the merging {\em Frontier Fields} (FF) cluster Abell~2744 using new Subaru/Suprime-Cam imaging. The wide-field lensing mass distribution reveals this cluster is comprised of four distinct substructures. Simultaneously modeling the two-dimensional reduced shear field using a combination of a Navarro--Frenk--White (NFW) model for the main core and truncated NFW models for the subhalos, we determine their masses and locations. The total mass of the system is constrained as $M_\mathrm{200c} = (2.06\pm0.42)\times10^{15}\,M_\odot$. The most massive clump is the southern component with $M_\mathrm{200c} = (7.7\pm3.4)\times10^{14}\,M_\odot$, followed by the western substructure ($M_\mathrm{200c} = (4.5\pm2.0)\times10^{14}\,M_\odot$) and two smaller substructures to the northeast ($M_\mathrm{200c} = (2.8\pm1.6)\times10^{14}\,M_\odot$) and northwest ($M_\mathrm{200c} = (1.9\pm1.2)\times10^{14}\,M_\odot$). The presence of the four substructures supports the picture of multiple mergers. Using a composite of hydrodynamical binary simulations we explain this complicated system without the need for a "slingshot" effect to produce the northwest X-ray interloper, as previously proposed. The locations of the substructures appear to be offset from both the gas ($87^{+34}_{-28}$ arcsec, 90\% CL) and the galaxies ($72^{+34}_{-53}$ arcsec, 90\% CL) in the case of the northwestern and western subhalos. To confirm or refute these findings, high resolution space-based observations extending beyond the current FF limited coverage to the west and northwestern area are essential.

We investigate the dust-obscured star formation properties of the massive, X-ray selected galaxy cluster MACS J1931.8-2634 at $z$=0.352. Using far-infrared (FIR) imaging in the range 100-500$\mu$m obtained with the \textit{Herschel} telescope, we extract 31 sources (2$\sigma$) within $r\sim$1 Mpc from the brightest cluster galaxy (BCG). Among these sources we identify six cluster members for which we perform an analysis of their spectral energy distributions (SEDs). We measure total infrared luminosity (L$_{IR}$), star formation rate (SFR) and dust temperature. The BCG, with L$_{IR}$=1.4$\times$10$^{12}$L$_\odot$ is an Ultra Luminous Infrared Galaxy and hosts a type II AGN. We decompose its FIR SED into AGN and starburst components and find equal contributions from AGN and starburst. We also recompute the SFR of the BCG finding SFR=150$\pm$15 M$_\odot$yr$^{-1}$. We search for an isobaric cooling flow in the cool core using {\sl Chandra} X-ray data, and find no evidence for gas colder than 1.8 keV in the inner 30 kpc, for an upper limit to the istantaneous mass-deposition rate of 58 M$_\odot$yr$^{-1}$ at 95 % c.l. This value is $3\times$ lower than the SFR in the BCG, suggesting that the on-going SF episode lasts longer than the ICM cooling events.

SN HFF14Tom is a Type Ia Supernova (SN) discovered at z = 1.3457 +- 0.0001 behind the galaxy cluster Abell 2744 (z = 0.308). In a cosmology-independent analysis, we find that HFF14Tom is 0.77 +- 0.15 magnitudes brighter than unlensed Type Ia SNe at similar redshift, implying a lensing magnification of mu_obs = 2.03 +- 0.29. This observed magnification provides a rare opportunity for a direct empirical test of galaxy cluster lens models. Here we test 17 lens models, 13 of which were generated before the SN magnification was known, qualifying as pure "blind tests". The models are collectively fairly accurate: 8 of the models deliver median magnifications that are consistent with the measured mu to within 1-sigma. However, there is a subtle systematic bias: the significant disagreements all involve models overpredicting the magnification. We evaluate possible causes for this mild bias, and find no single physical or methodological explanation to account for it. We do find that model accuracy can be improved to some extent with stringent quality cuts on multiply-imaged systems, such as requiring that a large fraction have spectroscopic redshifts. In addition to testing model accuracies as we have done here, Type Ia SN magnifications could also be used as inputs for future lens models of Abell 2744 and other clusters, providing valuable constraints in regions where traditional strong- and weak-lensing information is unavailable.

We perform a 3D multi-probe analysis of the rich galaxy cluster A1689 by combining improved weak-lensing data from new BVRi'z' Subaru/Suprime-Cam observations with strong-lensing, X-ray, and Sunyaev-Zel'dovich effect (SZE) data sets. We reconstruct the projected matter distribution from a joint weak-lensing analysis of 2D shear and azimuthally integrated magnification constraints, the combination of which allows us to break the mass-sheet degeneracy. The resulting mass distribution reveals elongation with axis ratio ~0.7 in projection. When assuming a spherical halo, our full weak-lensing analysis yields a projected concentration of $c_{200c}^{2D}=8.9\pm 1.1$ ($c_{vir}^{2D}\sim 11$), consistent with and improved from earlier weak-lensing work. We find excellent consistency between weak and strong lensing in the region of overlap. In a parametric triaxial framework, we constrain the intrinsic structure and geometry of the matter and gas distributions, by combining weak/strong lensing and X-ray/SZE data with minimal geometric assumptions. We show that the data favor a triaxial geometry with minor-major axis ratio 0.39+/-0.15 and major axis closely aligned with the line of sight (22+/-10 deg). We obtain $M_{200c}=(1.2\pm 0.2)\times 10^{15} M_{\odot}/h$ and $c_{200c}=8.4\pm 1.3$, which overlaps with the $>1\sigma$ tail of the predicted distribution. The shape of the gas is rounder than the underlying matter but quite elongated with minor-major axis ratio 0.60+/-0.14. The gas mass fraction within 0.9Mpc is 10^{+3}_{-2}%. The thermal gas pressure contributes to ~60% of the equilibrium pressure, indicating a significant level of non-thermal pressure support. When compared to Planck's hydrostatic mass estimate, our lensing measurements yield a spherical mass ratio of $M_{Planck}/M_{GL}=0.70\pm 0.15$ and $0.58\pm 0.10$ with and without corrections for lensing projection effects, respectively.

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 <b> = 0.12 for the WL comparison and <b> = -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 <b>~0.3 at ~1 Mpc for the WL comparison and <b> 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.

In the effort to understand the link between the structure of galaxy clusters and their galaxy populations, we focus on MACSJ1206.2-0847 at z~0.44 and probe its substructure in the projected phase space through the spectrophotometric properties of a large number of galaxies from the CLASH-VLT survey. Our analysis is mainly based on an extensive spectroscopic dataset of 445 member galaxies, mostly acquired with VIMOS@VLT as part of our ESO Large Programme, sampling the cluster out to a radius ~2R200 (4 Mpc). We classify 412 galaxies as passive, with strong Hdelta absorption (red and blue galaxies, and with emission lines from weak to very strong. A number of tests for substructure detection are applied to analyze the galaxy distribution in the velocity space, in 2D space, and in 3D projected phase-space. Studied in its entirety, the cluster appears as a large-scale relaxed system with a few secondary, minor overdensities in 2D distribution. We detect no velocity gradients or evidence of deviations in local mean velocities. The main feature is the WNW-ESE elongation. The analysis of galaxy populations per spectral class highlights a more complex scenario. The passive galaxies and red strong Hdelta galaxies trace the cluster center and the WNW-ESE elongated structure. The red strong Hdelta galaxies also mark a secondary, dense peak ~2 Mpc at ESE. The emission line galaxies cluster in several loose structures, mostly outside R200. The observational scenario agrees with MACS J1206.2-0847 having WNW-ESE as the direction of the main cluster accretion, traced by passive galaxies and red strong Hdelta galaxies. The red strong Hdelta galaxies, interpreted as poststarburst galaxies, date a likely important event 1-2 Gyr before the epoch of observation. The emission line galaxies trace a secondary, ongoing infall where groups are accreted along several directions.

We utilize the CLASH (Cluster Lensing And Supernova survey with Hubble) observations of 25 clusters to search for extreme emission-line galaxies (EELGs). The selections are carried out in two central bands: F105W (Y105) and F125W (J125), as the flux of the central bands could be enhanced by the presence of [O III] 4959, 5007 at redshift of about 0.93-1.14 and 1.57-1.79, respectively. The multi-band observations help to constrain the equivalent widths of emission lines. Thanks to cluster lensing, we are able to identify 52 candidates down to an intrinsic limiting magnitude of 28.5 and to a rest-frame [O III] 4959,5007 equivalent width of about 3737 angstrom. Our samples include a number of EELGs at lower luminosities that are missed in other surveys, and the extremely high equivalent width can be only found in such faint galaxies. These EELGs can mimic the dropout feature similar to that of high redshift galaxies and contaminate the color-color selection of high redshift galaxies when the S/N ratio is limited or the band coverage is incomplete. We predict that the fraction of EELGs in the future high redshift galaxy selections cannot be neglected.

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 3500kpc/h. 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 $M_{200c}=1.34^{+0.10}_{-0.09} 10^{15}M_{\odot}$. We show this is in excellent agreement with Lambda cold-dark-matter (LCDM) predictions when the CLASH X-ray selection function and projection effects are taken into account. The best-fit Einasto shape parameter is $\alpha_E=0.191^{+0.071}_{-0.068}$, which is consistent with the NFW-equivalent Einasto parameter of $\sim 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. 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 LCDM model.

We utilise a two-color Lyman-Break selection criterion to search for z~9-10 galaxies over the first 19 clusters in the CLASH program. A systematic search yields three z~9-10 candidates. While we have already reported the most robust of these candidates, MACS1149-JD, two additional z~9 candidates are also found and have H_{160}-band magnitudes of ~26.2-26.9. A careful assessment of various sources of contamination suggests <~1 contaminants for our z~9-10 selection. To determine the implications of these search results for the LF and SFR density at z~9, we introduce a new differential approach to deriving these quantities in lensing fields. Our procedure is to derive the evolution by comparing the number of z~9-10 galaxy candidates found in CLASH with the number of galaxies in a slightly lower redshift sample (after correcting for the differences in selection volumes), here taken to be z~8. This procedure takes advantage of the fact that the relative volumes available for the z~8 and z~9-10 selections behind lensing clusters are not greatly dependent on the details of the lensing models. We find that the normalization of the UV LF at z~9 is just 0.28_{-0.20}^{+0.39}\times that at z~8, ~1.4_{-0.8}^{+3.0}x lower than extrapolating z~4-8 LF results. While consistent with the evolution in the UV LF seen at z~4-8, these results marginally favor a more rapid evolution at z>8. Compared to similar evolutionary findings from the HUDF, our result is less insensitive to large-scale structure uncertainties, given our many independent sightlines on the high-redshift universe.

Context. The study of the galaxy stellar mass function (SMF) in relation to the galaxy environment and the stellar mass density profile, rho(r), is a powerful tool to constrain models of galaxy evolution. Aims. We determine the SMF of the z=0.44 cluster of galaxies MACS J1206.2-0847 separately for passive and star-forming (SF) galaxies, in different regions of the cluster, from the center out to approximately 2 virial radii. We also determine rho(r) to compare it to the number density and total mass density profiles. Methods. We use the dataset from the CLASH-VLT survey. Stellar masses are obtained by SED fitting on 5-band photometric data obtained at the Subaru telescope. We identify 1363 cluster members down to a stellar mass of 10^9.5 Msolar. Results. The whole cluster SMF is well fitted by a double Schechter function. The SMFs of cluster SF and passive galaxies are statistically different. The SMF of the SF cluster galaxies does not depend on the environment. The SMF of the passive population has a significantly smaller slope (in absolute value) in the innermost (<0.50 Mpc), highest density cluster region, than in more external, lower density regions. The number ratio of giant/subgiant galaxies is maximum in this innermost region and minimum in the adjacent region, but then gently increases again toward the cluster outskirts. This is also reflected in a decreasing radial trend of the average stellar mass per cluster galaxy. On the other hand, the stellar mass fraction, i.e., the ratio of stellar to total cluster mass, does not show any significant radial trend. Conclusions. Our results appear consistent with a scenario in which SF galaxies evolve into passive galaxies due to density-dependent environmental processes, and eventually get destroyed very near the cluster center to become part of a diffuse intracluster medium.

We utilize 16 band Hubble Space Telescope (HST) observations of 18 lensing clusters obtained as part of the Cluster Lensing And Supernova survey with Hubble (CLASH) Multi-Cycle Treasury program to search for $z\sim6-8$ galaxies. We report the discovery of 204, 45, and 13 Lyman-break galaxy candidates at $z\sim6$, $z\sim7$, and $z\sim8$, respectively, identified from purely photometric redshift selections. This large sample, representing nearly an order of magnitude increase in the number of magnified star-forming galaxies at $z\sim 6-8$ presented to date, is unique in that we have observations in four WFC3/UVIS UV, seven ACS/WFC optical, and all five WFC3/IR broadband filters, which enable very accurate photometric redshift selections. We construct detailed lensing models for 17 of the 18 clusters to estimate object magnifications and to identify two new multiply lensed $z \gtrsim 6$ candidates. The median magnifications over the 17 clusters are 4, 4, and 5 for the $z\sim6$, $z\sim7$, and $z\sim8$ samples, respectively, over an average area of 4.5 arcmin$^2$ per cluster. We compare our observed number counts with expectations based on convolving "blank" field UV luminosity functions through our cluster lens models and find rough agreement down to $\sim27$ mag, where we begin to suffer significant incompleteness. In all three redshift bins, we find a higher number density at brighter observed magnitudes than the field predictions, empirically demonstrating for the first time the enhanced efficiency of lensing clusters over field surveys. Our number counts also are in general agreement with the lensed expectations from the cluster models, especially at $z\sim6$, where we have the best statistics.

We present a new determination of the concentration-mass relation for galaxy clusters based on our comprehensive lensing analysis of 19 X-ray selected galaxy clusters from the Cluster Lensing and Supernova Survey with Hubble (CLASH). Our sample spans a redshift range between 0.19 and 0.89. We combine weak lensing constraints from the Hubble Space Telescope (HST) and from ground-based wide field data with strong lensing constraints from HST. The result are reconstructions of the surface-mass density for all CLASH clusters on multi-scale grids. Our derivation of NFW parameters yields virial masses between 0.53 x 10^15 and 1.76 x 10^15 M_sol/h and the halo concentrations are distributed around c_200c ~ 3.7 with a 1-sigma significant negative trend with cluster mass. We find an excellent 4% agreement between our measured concentrations and the expectation from numerical simulations after accounting for the CLASH selection function based on X-ray morphology. The simulations are analyzed in 2D to account for possible biases in the lensing reconstructions due to projection effects. The theoretical concentration-mass (c-M) relation from our X-ray selected set of simulated clusters and the c-M relation derived directly from the CLASH data agree at the 90% confidence level.

We present the results of a numerical study based on the analysis of the MUSIC-2 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 lines-of-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 un-relaxed systems. For such a 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.

We aim at constraining the assembly history of clusters by studying the intra cluster light (ICL) properties, estimating its contribution to the fraction of baryons in stars, f*, and understanding possible systematics/bias using different ICL detection techniques. We developed an automated method, GALtoICL, based on the software GALAPAGOS to obtain a refined version of typical BCG+ICL maps. We applied this method to our test case MACS J1206.2-0847, a massive cluster located at z=0.44, that is part of the CLASH sample. Using deep multi-band SUBARU images, we extracted the surface brightness (SB) profile of the BCG+ICL and we studied the ICL morphology, color, and contribution to f* out to R500. We repeated the same analysis using a different definition of the ICL, SBlimit method, i.e. a SB cut-off level, to compare the results. The most peculiar feature of the ICL in MACS1206 is its asymmetric radial distribution, with an excess in the SE direction and extending towards the 2nd brightest cluster galaxy which is a Post Starburst galaxy. This suggests an interaction between the BCG and this galaxy that dates back to t <= 1.5 Gyr. The BCG+ICL stellar content is 8% of M_(*,500) and the (de-) projected baryon fraction in stars is f*=0.0177 (0.0116), in excellent agreement with recent results. The SBlimit method provides systematically higher ICL fractions and this effect is larger at lower SB limits. This is due to the light from the outer envelopes of member galaxies that contaminate the ICL. Though more time consuming, the GALtoICL method provides safer ICL detections that are almost free of this contamination. This is one of the few ICL study at redshift z > 0.3. At completion, the CLASH/VLT program will allow us to extend this analysis to a statistically significant cluster sample spanning a wide redshift range: 0.2<z<0.6.

We report observations of three gravitationally lensed supernovae (SNe) in the Cluster Lensing And Supernova survey with Hubble (CLASH) Multi-Cycle Treasury program. These objects, SN CLO12Car (z = 1.28), SN CLN12Did (z = 0.85), and SN CLA11Tib (z = 1.14), are located behind three different clusters, MACSJ1720.2+3536 (z = 0.391), RXJ1532.9+3021 (z = 0.345), and Abell 383 (z = 0.187), respectively. Each SN was detected in Hubble Space Telescope (HST) optical and infrared images. Based on photometric classification, we find that SNe CLO12Car and CLN12Did are likely to be Type Ia supernovae (SNe Ia), while the classification of SN CLA11Tib is inconclusive. Using multi-color light-curve fits to determine a standardized SN Ia luminosity distance, we infer that SN CLO12Car was approximately 1.0 +/- 0.2 mag brighter than field SNe Ia at a similar redshift and ascribe this to gravitational lens magnification. Similarly, SN CLN12Did is approximately 0.2 +/- 0.2 mag brighter than field SNe Ia. We derive independent estimates of the predicted magnification from CLASH strong+weak lensing maps of the clusters: 0.83 +/- 0.16 mag for SN CLO12Car, 0.28 +/- 0.08 mag for SN CLN12Did, and 0.43 +/- 0.11 mag for SN CLA11Tib. The two SNe Ia provide a new test of the cluster lens model predictions: we find that the magnifications based on the SN Ia brightness and those predicted by the lens maps are consistent. Our results herald the promise of future observations of samples of cluster-lensed SNe Ia (from the ground or space) to help illuminate the dark-matter distribution in clusters of galaxies, through the direct determination of absolute magnifications.

We present a strong lensing system in which a double source is imaged 5 times by 2 early-type galaxies. We take advantage in this target of the multi-band photometry obtained as part of the CLASH program, complemented by the spectroscopic data of the VLT/VIMOS and FORS2 follow-up campaign. We use a photometric redshift of 3.7 for the source and confirm spectroscopically the membership of the 2 lenses to the galaxy cluster MACS J1206.2-0847 at redshift 0.44. We exploit the excellent angular resolution of the HST/ACS images to model the 2 lenses in terms of singular isothermal sphere profiles and derive robust effective velocity dispersions of (97 +/- 3) and (240 +/- 6) km/s. The total mass distribution of the cluster is also well characterized by using only the local information contained in this lensing system, that is located at a projected distance of more than 300 kpc from the cluster luminosity center. According to our best-fitting lensing and composite stellar population models, the source is magnified by a total factor of 50 and has a luminous mass of about (1.0 +/- 0.5) x 10^{9} M_{Sun}. By combining the total and luminous mass estimates of the 2 lenses, we measure luminous over total mass fractions projected within the effective radii of 0.51 +/- 0.21 and 0.80 +/- 0.32. With these lenses we can extend the analysis of the mass properties of lens early-type galaxies by factors that are about 2 and 3 times smaller than previously done with regard to, respectively, velocity dispersion and luminous mass. The comparison of the total and luminous quantities of our lenses with those of astrophysical objects with different physical scales reveals the potential of studies of this kind for investigating the internal structure of galaxies. These studies, made possible thanks to the CLASH survey, will allow us to go beyond the current limits posed by the available lens samples in the field.

A pressureless scenario for the Dark Matter (DM) fluid is a widely adopted hypothesis, despite the absence of a direct observational evidence. According to general relativity, the total mass-energy content of a system shapes the gravitational potential well, but different test particles perceive this potential in different ways depending on their properties. Cluster galaxy velocities, being $\ll$c, depend solely on the gravitational potential, whereas photon trajectories reflect the contributions from the gravitational potential plus a relativistic-pressure term that depends on the cluster mass. We exploit this phenomenon to constrain the Equation of State (EoS) parameter of the fluid, primarily DM, contained in galaxy clusters. We use the complementary information provided by the kinematic and lensing mass profiles of the galaxy cluster MACS 1206.2-0847 at $z=0.44$, as obtained in an extensive imaging and spectroscopic campaign within the CLASH survey. The unprecedented high quality of our data-set and the properties of this cluster are well suited to determine the EoS parameter of the cluster fluid. Since baryons contribute at most $15\%$ to the total mass in clusters and their pressure is negligible, the EoS parameter we derive describes the behavior of the DM fluid. We obtain the most stringent constraint on the DM EoS parameter to date, $w=(p_r+2\,p_t)/(3\,c^2\rho)=0.00\pm0.15\mathrm{(stat)}\pm0.08\mathrm{(syst)}$, averaged over the radial range $0.5\,\mathrm{Mpc}\leq$$r$$\leq$$r_{200}$, where $p_r$ and $p_t$ are the radial and tangential pressure, and $\rho$ is the density. We plan to further improve our constraint by applying the same procedure to all clusters from the ongoing CLASH-VLT program.

We present the supernova (SN) sample and Type-Ia SN (SN Ia) rates from the Cluster Lensing And Supernova survey with Hubble (CLASH). Using the Advanced Camera for Surveys and the Wide Field Camera 3 on the Hubble Space Telescope (HST), we have imaged 25 galaxy-cluster fields and parallel fields of non-cluster galaxies. We report a sample of 27 SNe discovered in the parallel fields. Of these SNe, ~13 are classified as SN Ia candidates, including four SN Ia candidates at redshifts z > 1.2. We measure volumetric SN Ia rates to redshift 1.8 and add the first upper limit on the SN Ia rate in the range 1.8 < z < 2.4. The results are consistent with the rates measured by the HST/GOODS and Subaru Deep Field SN surveys. We model these results together with previous measurements at z < 1 from the literature. The best-fitting SN Ia delay-time distribution (DTD; the distribution of times that elapse between a short burst of star formation and subsequent SN Ia explosions) is a power law with an index of -1.00 +0.06(0.09) -0.06(0.10) (statistical) +0.12 -0.08 (systematic), where the statistical uncertainty is a result of the 68% and 95% (in parentheses) statistical uncertainties reported for the various SN Ia rates (from this work and from the literature), and the systematic uncertainty reflects the range of possible cosmic star-formation histories. We also test DTD models produced by an assortment of published binary population synthesis (BPS) simulations. The shapes of all BPS double-degenerate DTDs are consistent with the volumetric SN Ia measurements, when the DTD models are scaled up by factors of 3-9. In contrast, all BPS single-degenerate DTDs are ruled out by the measurements at a >99% significance level.

The Cluster Lensing And Supernovae survey with Hubble (CLASH) is an Hubble Space Telescope (HST) Multi-Cycle Treasury program observing 25 massive galaxy clusters. CLASH observations are carried out in 16 bands from UV to NIR to derive accurate and reliable estimates of photometric redshifts. We present the CLASH photometric redshifts and study the photometric redshift accuracy of the arcs in more detail for the case of MACS1206.2-0847. We use the publicly available Le Phare and BPZ photometric redshift codes on 17 CLASH galaxy clusters. Using Le Phare code for objects with StoN>=10, we reach a precision of 3%(1+z) for the strong lensing arcs, which is reduced to 2.4%(1+z) after removing outliers. For galaxies in the cluster field the corresponding values are 4%(1+z) and 3%(1+z). Using mock galaxy catalogues, we show that 3%(1+z) precision is what one would expect from the CLASH photometry when taking into account extinction from dust, emission lines and the finite range of SEDs included in the photo-z template library. We study photo-z results for different aperture photometry and find that the SExtractor isophotal photometry works best. Le Phare and BPZ give similar photo-z results for the strong lensing arcs as well as galaxies of the cluster field. Results are improved when optimizing the photometric aperture shape showing an optimal aperture size around 1" radius giving results which are equivalent to isophotal photometry. Tailored photometry of the arcs improve the photo-z results.

We present a quintuply lensed z ~ 6 candidate discovered in the field of the galaxy cluster RXC J2248.7-4431 (z ~ 0.348) targeted within the Cluster Lensing and Supernova survey with Hubble (CLASH) and selected in the deep HST Frontier Fields survey. Thanks to the CLASH 16-band HST imaging, we identify the quintuply lensed z ~ 6 candidate as an optical dropout in the inner region of the cluster, the brightest image having magAB=24.81+-0.02 in the f105w filter. We perform a detailed photometric analysis to verify its high-z and lensed nature. We get as photometric redshift z_phot ~ 5.9, and given the extended nature and NIR colours of the lensed images, we rule out low-z early type and galactic star contaminants. We perform a strong lensing analysis of the cluster, using 13 families of multiple lensed images identified in the HST images. Our final best model predicts the high-z quintuply lensed system with a position accuracy of 0.8''. The magnifications of the five images are between 2.2 and 8.3, which leads to a delensed UV luminosity of L_1600 ~ 0.5L*_1600 at z=6. We also estimate the UV slope from the observed NIR colours, finding a steep beta=-2.89+-0.38. We use singular and composite stellar population SEDs to fit the photometry of the hiz candidate, and we conclude that it is a young (age <300 Myr) galaxy with mass of M ~ 10^8Msol, subsolar metallicity (Z<0.2Zsol) and low dust content (AV ~ 0.2-0.4).

We present VIsible Multi-Object Spectrograph (VIMOS) observations of a z 6 galaxy quintuply imaged by the Frontier Fields galaxy cluster RXC J2248.7-4431 (z=0.348). This sub-L^*, high-z galaxy has been recently discovered by Monna et al. (2013) using dropout techniques with the 16-band HST photometry acquired as part of the Cluster Lensing And Supernova survey with Hubble (CLASH). Obtained as part of the CLASH-VLT survey, the VIMOS medium-resolution spectra of this source show a very faint continuum between ~8700A and ~9300A and a prominent emission line at 8643A, which can be readily identified with Lyman-alpha at z=6.110. The emission line exhibits an asymmetric profile, with a more pronounced red wing. The rest-frame equivalent width of the line is EW=79+-10A. After correcting for magnification, the star-formation rate (SFR) estimated from the Lya line is SFR(Lya)=11 M_{sol}/yr and that estimated from the UV data is SFR(UV)=3 M_{sol}/yr. We estimate that the effective radius of the source is R_e<~0.4 kpc, which implies a star formation surface mass density Sigma_{SFR}>6 M_{sol}/yr/kpc^2 and, using the Kennicutt-Schmidt relation, a gas surface mass density Sigma_{gas}>10^3 M_{sol}/pc^2. Our results support the idea that this magnified, distant galaxy is a young and compact object with 0.4 L^* at z=6, with comparable amount of mass in gas and stars. Future follow-up observations with ALMA will provide valuable insight into the SFR and molecular gas content of this source. In the spirit of the Frontier Fields initiative, we also publish the redshifts of several multiply imaged sources and other background objects which will help improving the strong lensing model of this galaxy cluster.

The galaxy cluster MACS J0717.5+3745 (z=0.55) is the largest known cosmic lens, with complex internal structures seen in deep X-ray, Sunyaev-Zel'dovich effect and dynamical observations. We perform a combined weak and strong lensing analysis with wide-field BVRi'z' Subaru/Suprime-Cam observations and 16-band Hubble Space Telescope observations taken as part of the Cluster Lensing And Supernova survey with Hubble (CLASH). We find consistent weak distortion and magnification measurements of background galaxies, and combine these signals to construct an optimally estimated radial mass profile of the cluster and its surrounding large-scale structure out to 5 Mpc/h. We find consistency between strong-lensing and weak-lensing in the region where these independent data overlap, <500 kpc/h. The two-dimensional weak-lensing map reveals a clear filamentary structure traced by distinct mass halos. We model the lensing shear field with 9 halos, including the main cluster, corresponding to mass peaks detected above 2.5\sigma_\kappa. The total mass of the cluster as determined by the different methods is M_{vir}=(2.8\pm0.4) \times 10^15 M_sun. Although this is the most massive cluster known at z>0.5, in terms of extreme value statistics we conclude that the mass of MACS J0717.5+3745 by itself is not in serious tension with LambdaCDM, representing only a ~2{\sigma} departure above the maximum simulated halo mass at this redshift.

We provide a new observational test for a key prediction of the \Lambda CDM cosmological model: the contributions of mergers with different halo-to-main-cluster mass ratios to cluster-sized halo growth. We perform this test by dynamically analyzing seven galaxy clusters, spanning the redshift range $0.13 < z_c < 0.45$ and caustic mass range $0.4-1.5$ $10^{15} h_{0.73}^{-1}$ M$_{\odot}$, with an average of 293 spectroscopically-confirmed bound galaxies to each cluster. The large radial coverage (a few virial radii), which covers the whole infall region, with a high number of spectroscopically identified galaxies enables this new study. For each cluster, we identify bound galaxies. Out of these galaxies, we identify infalling and accreted halos and estimate their masses and their dynamical states. Using the estimated masses, we derive the contribution of different mass ratios to cluster-sized halo growth. For mass ratios between ~0.2 and ~0.7, we find a ~1 $\sigma$ agreement with \Lambda CDM expectations based on the Millennium simulations I and II. At low mass ratios, $\lesssim 0.2$, our derived contribution is underestimated since the detection efficiency decreases at low masses, $\sim 2 \times 10^{14}$ $h_{0.73}^{-1}$ M$_{\odot}$. At large mass ratios, $\gtrsim 0.7$, we do not detect halos probably because our sample, which was chosen to be quite X-ray relaxed, is biased against large mass ratios. Therefore, at large mass ratios, the derived contribution is also underestimated.

We use an unprecedented data-set of about 600 redshifts for cluster members, obtained as part of a VLT/VIMOS large programme, to constrain the mass profile of the z=0.44 cluster MACS J1206.2-0847 over the radial range 0-5 Mpc (0-2.5 virial radii) using the MAMPOSSt and Caustic methods. We then add external constraints from our previous gravitational lensing analysis. We invert the Jeans equation to obtain the velocity-anisotropy profiles of cluster members. With the mass-density and velocity-anisotropy profiles we then obtain the first determination of a cluster pseudo-phase-space density profile. The kinematics and lensing determinations of the cluster mass profile are in excellent agreement. This is very well fitted by a NFW model with mass M200=(1.4 +- 0.2) 10^15 Msun and concentration c200=6 +- 1, only slightly higher than theoretical expectations. Other mass profile models also provide acceptable fits to our data, of (slightly) lower (Burkert, Hernquist, and Softened Isothermal Sphere) or comparable (Einasto) quality than NFW. The velocity anisotropy profiles of the passive and star-forming cluster members are similar, close to isotropic near the center and increasingly radial outside. Passive cluster members follow extremely well the theoretical expectations for the pseudo-phase-space density profile and the relation between the slope of the mass-density profile and the velocity anisotropy. Star-forming cluster members show marginal deviations from theoretical expectations. This is the most accurate determination of a cluster mass profile out to a radius of 5 Mpc, and the only determination of the velocity-anisotropy and pseudo-phase-space density profiles of both passive and star-forming galaxies for an individual cluster [abridged]

Growing observational evidence now indicates that nebular line emission has a significant impact on the rest-frame optical fluxes of z~5-7 galaxies observed with Spitzer. This line emission makes z~5-7 galaxies appear more massive, with lower specific star formation rates. However, corrections for this line emission have been very difficult to perform reliably due to huge uncertainties on the overall strength of such emission at z>~5.5. Here, we present the most direct observational evidence yet for ubiquitous high-EW [OIII]+Hbeta line emission in Lyman-break galaxies at z~7, while also presenting a strategy for an improved measurement of the sSFR at z~7. We accomplish this through the selection of bright galaxies in the narrow redshift window z~6.6-7.0 where the IRAC 4.5 micron flux provides a clean measurement of the stellar continuum light. Observed 4.5 micron fluxes in this window contrast with the 3.6 micron fluxes which are contaminated by the prominent [OIII]+Hbeta lines. To ensure a high S/N for our IRAC flux measurements, we consider only the brightest (H_{160}<26 mag) magnified galaxies we have identified in CLASH and other programs targeting galaxy clusters. Remarkably, the mean rest-frame optical color for our bright seven-source sample is very blue, [3.6]-[4.5]=-0.9+/-0.3. Such blue colors cannot be explained by the stellar continuum light and require that the rest-frame EW of [OIII]+Hbeta be greater than 637 Angstroms for the average source. The bluest four sources from our seven-source sample require an even more extreme EW of 1582 Angstroms. Our derived lower limit for the mean [OIII]+Hbeta EW could underestimate the true EW by ~2x based on a simple modeling of the redshift distribution of our sources. We can also set a robust lower limit of >~4 Gyr^-1 on the specific star formation rates based on the mean SED for our seven-source sample. (abridged)

In this work we analyze the mass distribution of MACSJ1206.2-0847, especially focusing on the halo properties of its cluster members. The cluster appears relaxed in its X-ray emission, but has significant amounts of intracluster light which is not centrally concentrated, suggesting that galaxy-scale interactions are still ongoing despite the overall relaxed state. The cluster lenses 12 background galaxies into multiple images and one galaxy at $z=1.033$ into a giant arc and its counterimage. The multiple image positions and the surface brightness distribution (SFB) of the arc which is bent around several cluster members are sensitive to the cluster galaxy halo properties. We model the cluster mass distribution with a NFW profile and the galaxy halos with two parameters for the mass normalization and extent of a reference halo assuming scalings with their observed NIR--light. We match the multiple image positions at an r.m.s. level of $0.85\arcsec$ and can reconstruct the SFB distribution of the arc in several filters to a remarkable accuracy based on this cluster model. The length scale where the enclosed galaxy halo mass is best constrained is about 5 effective radii -- a scale in between those accessible to dynamical and field strong lensing mass estimates on one hand and galaxy--galaxy weak lensing results on the other hand. The velocity dispersion and halo size of a galaxy with $m_{\rm 160W,AB}=19.2$ or $M_{\rm B,Vega}=-20.7$ are $\sigma=150 \rm kms^{-1}$ and $r\approx 26\pm 6 \rm kpc$, indicating that the halos of the cluster galaxies are tidally stripped. We also reconstruct the unlensed source (which is smaller by a factor of $\sim5.8$ in area), demonstrating the increase of morphological information due to lensing and conclude that this galaxy has likely star--forming spiral arms with a red (older) central component.

We perform a strong-lensing analysis of the merging galaxy cluster MACS J0416.1-2403 (M0416; z=0.42) in recent CLASH/HST observations. We identify 70 new multiple images and candidates of 23 background sources in the range 0.7<z_{phot}<6.14 including two probable high-redshift dropouts, revealing a highly elongated lens with axis ratio ~5:1, and a major axis of ~100\arcsec (z_{s}~2). Compared to other well-studied clusters, M0416 shows an enhanced lensing efficiency. Although the critical area is not particularly large (~0.6 \square\arcmin; z_{s}~2), the number of multiple images, per critical area, is anomalously high. We calculate that the observed elongation boosts the number of multiple images, \emph{per critical area}, by a factor of ~2.5\times, due to the increased ratio of the caustic area relative to the critical area. Additionally, we find that the observed separation between the two main mass components enlarges the critical area by a factor of ~2. These geometrical effects can account for the high number (density) of multiple images observed. We find in numerical simulations, that only ~4% of the clusters (with M_{vir}>6 x 10^{14} h^{-1}M_{\odot}) exhibit as elongated critical curves as M0416.

Hubble Space Telescope images of the galaxy cluster Abell 2261, obtained as part of the Cluster Lensing And Supernova survey with Hubble, show that the brightest galaxy in the cluster, A2261-BCG, has the largest core yet detected in any galaxy. The cusp radius of A2261-BCG is 3.2 kpc, twice as big as the next largest core known, and ~3x bigger than those typically seen in the most luminous BCGs. The morphology of the core in A2261-BCG is also unusual, having a flat or even slightly-depressed interior surface brightness profile, rather than the typical shallow cusp. This implies that the galaxy has a core with constant or even centrally decreasing stellar density. Interpretation of the core as an end product of the "scouring" action of a binary supermassive black hole implies a total black hole mass ~1E+10 M_sun from the extrapolation of most relationships between core structure and black hole mass. The core falls 1-sigma above the cusp-radius versus galaxy luminosity relation. Its large size in real terms, and the extremely large black hole mass required to generate it, raise the possibility that the core has been enlarged by additional processes, such as the ejection of the black holes that originally generated the core. The flat central stellar density profile is consistent with this hypothesis. The core is also displaced by 0.7 kpc from the center of the surrounding envelope, consistent with a local dynamical perturbation of the core.

We examine the inner mass distribution of the relaxed galaxy cluster Abell 383 in deep 16-band HST/ACS+WFC3 imaging taken as part of the CLASH multi-cycle treasury program. Our program is designed to study the dark matter distribution in 25 massive clusters, and balances depth with a wide wavelength coverage to better identify lensed systems and generate precise photometric redshifts. This information together with the predictive strength of our strong-lensing analysis method identifies 13 new multiply-lensed images and candidates, so that a total of 27 multiple-images of 9 systems are used to tightly constrain the inner mass profile, $d\log \Sigma/d\log r\simeq -0.6\pm 0.1$ (r<160 kpc). We find consistency with the standard distance-redshift relation for the full range spanned by the lensed images, 1.01<z<6.03, with the higher redshift sources deflected through larger angles as expected. The inner mass profile derived here is consistent with the results of our independent weak-lensing analysis of wide-field Subaru images, with good agreement in the region of overlap. The overall mass profile is well fitted by an NFW profile with M_{vir}=(5.37^{+0.70}_{-0.63}\pm 0.26) x 10^{14}M_{\odot}/h and a relatively high concentration, c_{vir}=8.77^{+0.44}_{-0.42}\pm 0.23, which lies above the standard c-M relation similar to other well-studied clusters. The critical radius of Abell 383 is modest by the standards of other lensing clusters, r_{E}\simeq16\pm2\arcsec (for z_s=2.55), so the relatively large number of lensed images uncovered here with precise photometric redshifts validates our imaging strategy for the CLASH survey. In total we aim to provide similarly high-quality lensing data for 25 clusters, 20 of which are X-ray selected relaxed clusters, enabling a precise determination of the representative mass profile free from lensing bias. (ABRIDGED)

We directly construct model-independent mass profiles of galaxy clusters from combined weak-lensing distortion and magnification measurements within a Bayesian statistical framework,which allows for a full parameter-space extraction of the underlying signal. This method applies to the full range of radius outside the Einstein radius, and recovers the absolute mass normalization. We apply our method to deep Subaru imaging of five high-mass (>10^{15}M_{sun}) clusters, A1689, A1703, A370, Cl0024+17, and RXJ1347-11, to obtain accurate profiles to beyond the virial radius (r_{vir}). For each cluster the lens distortion and magnification data are shown to be consistent with each other, and the total signal-to-noise ratio of the combined measurements ranges from 13 to 24 per cluster. We form a model-independent mass profile from stacking the clusters, which is detected at 37{\sigma} out to R ~ 1.7r_{vir}. The projected logarithmic slope steepens from -1.01 \pm 0.09 at R ~ 0.1r_{vir} to -1.92 \pm 0.51 at R ~ 0.9r_{vir}. We also derive for each cluster inner strong-lensing based mass profiles from deep HST/ACS observations, which we show overlap well with the outer Subaru-based profiles and together are well described by a generalized form of the Navarro-Frenk-White profile, except for the ongoing merger RXJ1347-11, with modest variations in the central cusp slope (-dln{\rho}/dlnr < 0.9). The improvement here from adding the magnification measurements is significant, ~30% in terms of cluster mass profile measurements, compared with the lensing distortion signal.

We derive an accurate mass distribution of the galaxy cluster MACS J1206.2-0847 (z=0.439) from a combined weak-lensing distortion, magnification, and strong-lensing analysis of wide-field Subaru BVRIz' imaging and our recent 16-band Hubble Space Telescope observations taken as part of the Cluster Lensing And Supernova survey with Hubble (CLASH) program. We find good agreement in the regions of overlap between several weak and strong lensing mass reconstructions using a wide variety of modeling methods, ensuring consistency. The Subaru data reveal the presence of a surrounding large scale structure with the major axis running approximately north-west south-east (NW-SE), aligned with the cluster and its brightest galaxy shapes, showing elongation with a \sim 2:1 axis ratio in the plane of the sky. Our full-lensing mass profile exhibits a shallow profile slope dln\Sigma/dlnR\sim -1 at cluster outskirts (R>1Mpc/h), whereas the mass distribution excluding the NW-SE excess regions steepens further out, well described by the Navarro-Frenk-White form. Assuming a spherical halo, we obtain a virial mass M_{vir}=(1.1\pm 0.2\pm 0.1)\times 10^{15} M_{sun}/h and a halo concentration c_{vir} = 6.9\pm 1.0\pm 1.2 (\sim 5.7 when the central 50kpc/h is excluded), which falls in the range 4< <c> <7 of average c(M,z) predictions for relaxed clusters from recent Lambda cold dark matter simulations. Our full lensing results are found to be in agreement with X-ray mass measurements where the data overlap, and when combined with Chandra gas mass measurements, yield a cumulative gas mass fraction of 13.7^{+4.5}_{-3.0}% at 0.7Mpc/h (\approx 1.7r_{2500}), a typical value observed for high mass clusters.

The early Universe at redshift z\sim6-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 \lesssim 500 million years (Myr, at z \gtrsim 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 sigma) in the two reddest bands of HST/WFC3, and not detected below 1.2 {\mu}m, matching the characteristics of z\sim9 objects. We derive a robust photometric redshift of z = 9.6 \pm 0.2, corresponding to a cosmic age of 490 \pm 15Myr (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 \sim15) allows us to analyze the object's ultra-violet and optical luminosity in its rest-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 zf \lesssim 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.

We report the discovery of a z_{phot}=6.18^{+0.05}_{-0.07} (95% confidence level) dwarf galaxy, lensed into four images by the galaxy cluster MACS J0329.6-0211 (z_{l}=0.45). The galaxy is observed as a high-redshift dropout in HST/ACS/WFC3 CLASH and Spitzer/IRAC imaging. Its redshift is securely determined due to a clear detection of the Lyman-break in the 18-band photometry, making this galaxy one of the highest-redshift multiply-lensed objects known to date with an observed magnitude of F125W=24.00\pm0.04 AB mag for its highest-magnified image. We also present the first strong-lensing analysis of this cluster uncovering 15 additional multiply-imaged candidates of five lower-redshift sources spanning the range z_{s}~2-4. The mass model independently supports the high photometric redshift and reveals magnifications of 11.6^{+8.9}_{-4.1}, 17.6^{+6.2}_{-3.9}, 3.9^{+3.0}_{-1.7}, and 3.7^{+1.3}_{-0.2}, respectively, for the four images of the high-redshift galaxy. With this we construct a source image with a physical resolution of ~200 pc when the universe was ~0.9 Gyr old, where the z~6.2 galaxy occupies a source-plane area of approximately 2.2 kpc^{2}. Modeling the observed spectral energy distribution using population synthesis models, we find a demagnified stellar mass of ~10^{9} {M}_{sun}, subsolar metallicity (Z/Z_{sun}~0.5), low dust content (A_{V}~0.1 mag), a demagnified star formation rate (SFR) of ~3.2 {M}_{sun} yr^{-1}, and a specific SFR of ~3.4 Gyr^{-1}, all consistent with the properties of local dwarf galaxies.

We precisely constrain the inner mass profile of Abell 2261 (z=0.225) for the first time and determine this cluster is not "over-concentrated" as found previously, implying a formation time in agreement with {\Lambda}CDM expectations. These results are based on strong lensing analyses of new 16-band HST imaging obtained as part of the Cluster Lensing and Supernova survey with Hubble (CLASH). Combining this with revised weak lensing analyses of Subaru wide field imaging with 5-band Subaru + KPNO photometry, we place tight new constraints on the halo virial mass M_vir = 2.2\pm0.2\times10^15 M\odot/h70 (within r \approx 3 Mpc/h70) and concentration c = 6.2 \pm 0.3 when assuming a spherical halo. This agrees broadly with average c(M,z) predictions from recent {\Lambda}CDM simulations which span 5 <~ <c> <~ 8. Our most significant systematic uncertainty is halo elongation along the line of sight. To estimate this, we also derive a mass profile based on archival Chandra X-ray observations and find it to be ~35% lower than our lensing-derived profile at r2500 ~ 600 kpc. Agreement can be achieved by a halo elongated with a ~2:1 axis ratio along our line of sight. For this elongated halo model, we find M_vir = 1.7\pm0.2\times10^15 M\odot/h70 and c_vir = 4.6\pm0.2, placing rough lower limits on these values. The need for halo elongation can be partially obviated by non-thermal pressure support and, perhaps entirely, by systematic errors in the X-ray mass measurements. We estimate the effect of background structures based on MMT/Hectospec spectroscopic redshifts and find these tend to lower Mvir further by ~7% and increase cvir by ~5%.

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.

We present a strong-lensing analysis of the galaxy cluster MACS J1206.2-0847 ($z$=0.44) using UV, Optical, and IR, HST/ACS/WFC3 data taken as part of the CLASH multi-cycle treasury program, with VLT/VIMOS spectroscopy for some of the multiply-lensed arcs. The CLASH observations, combined with our mass-model, allow us to identify 47 new multiply-lensed images of 12 distant sources. These images, along with the previously known arc, span the redshift range $1\la z\la5.5$, and thus enable us to derive a detailed mass distribution and to accurately constrain, for the first time, the inner mass-profile of this cluster. We find an inner profile slope of $d\log \Sigma/d\log \theta\simeq -0.55\pm 0.1$ (in the range [1\arcsec, 53\arcsec], or $5\la r \la300$ kpc), as commonly found for relaxed and well-concentrated clusters. Using the many systems uncovered here we derive credible critical curves and Einstein radii for different source redshifts. For a source at $z_{s}\simeq2.5$, the critical curve encloses a large area with an effective Einstein radius of $\theta_{E}=28\pm3\arcsec$, and a projected mass of $1.34\pm0.15\times10^{14} M_{\odot}$. From the current understanding of structure formation in concordance cosmology, these values are relatively high for clusters at $z\sim0.5$, so that detailed studies of the inner mass distribution of clusters such as MACS J1206.2-0847 can provide stringent tests of the $\Lambda$CDM paradigm.

Extensive surveys of galaxy clusters motivate us to assess the likelihood of cluster-cluster lensing (CCL), namely, gravitational-lensing of a background cluster by a foreground cluster. We briefly describe the characteristics of CCLs in optical, X-ray and SZ measurements, and calculate their predicted numbers for $\Lambda$CDM parameters and a viable range of cluster mass functions and their uncertainties. The predicted number of CCLs in the strong-lensing regime varies from several ($<10$) to as high as a few dozen, depending mainly on whether lensing triaxiality bias is accounted for, through the c-M relation. A much larger number is predicted when taking into account also CCL in the weak-lensing regime. In addition to few previously suggested CCLs, we report a detection of a possible CCL in A383, where background candidate high-$z$ structures are magnified, as seen in deep Subaru observations.

We outline our methods for obtaining high precision mass profiles, combining independent weak-lensing distortion, magnification, and strong-lensing measurements. For massive clusters the strong and weak lensing regimes contribute equal logarithmic coverage of the radial profile. The utility of high-quality data is limited by the cosmic noise from large scale structure along the line of sight. This noise is overcome when stacking clusters, as too are the effects of cluster asphericity and substructure, permitting a stringent test of theoretical models. We derive a mean radial mass profile of four similar mass clusters of high-quality HST and Subaru images, in the range R=40kpc/h to 2800kpc/h, where the inner radial boundary is sufficiently large to avoid smoothing from miscentering effects. The stacked mass profile is detected at 58-sigma significance over the entire radial range, with the contribution from the cosmic noise included. We show that the projected mass profile has a continuously steepening gradient out to beyond the virial radius, in remarkably good agreement with the standard Navarro-Frenk-White form predicted for the family of CDM-dominated halos in gravitational equilibrium. The central slope is constrained to lie in the range, -dln{\rho}/dln{r}=0.89^{+0.27}_{-0.39}. The mean concentration is c_{vir}=7.68^{+0.42}_{-0.40} (at a mean virial mass 1.54^{+0.11}_{-0.10}\times 10^{15} M_{sun}/h), which is high for relaxed, high-mass clusters, but consistent with LCDM when a sizable projection bias estimated from N-body simulations is considered. This possible tension will be more definitively explored with new cluster surveys, such as CLASH, LoCuSS, Subaru HSC, and XXM-XXL, to construct the c-M relation over a wider mass range.

The amplitude of weak lensing should increase with source distance, rising steeply behind a lens and saturating at high redshift, providing a model-independent means of measuring cosmic geometry. We measure the amplitude of weak lensing with redshift for three massive clusters, A370 (z=0.375), ZwCl0024+17 (z=0.395) and RXJ1347-11 (z=0.451), using deep, three-colour Subaru imaging. We define the depth of lensed populations with reference to the COSMOS and GOODS fields, providing a consistency check of photo-z estimates over a wide range of redshift and magnitude. The predicted distance-redshift relation is followed well for the deepest dataset, A370, for a wide range of cosmologies, and is consistent with less accurate data for the other two clusters. Scaling this result to a new survey of ~25 massive clusters should provide a useful cosmological constraint on w, complementing existing techniques, with distance measurements covering the untested redshift range, 1<z<5.

The inner mass-profile of the relaxed cluster Abell 1703 is analysed by two very different strong-lensing techniques applied to deep ACS and WFC3 imaging. Our parametric method has the accuracy required to reproduce the many sets of multiple images, based on the assumption that mass approximately traces light. We test this assumption with a fully non-parametric, adaptive grid method, with no knowledge of the galaxy distribution. Differences between the methods are seen on fine scales due to member galaxies which must be included in models designed to search for lensed images, but on the larger scale the general distribution of dark matter is in good agreement, with very similar radial mass profiles. We add undiluted weak-lensing measurements from deep multi-colour Subaru imaging to obtain a fully model-independent mass profile out to the virial radius and beyond. Consistency is found in the region of overlap between the weak and strong lensing, and the full mass profile is well-described by an NFW model of a concentration parameter, $c_{\rm vir}\simeq 7.15\pm0.5$ (and $M_{vir}\simeq 1.22\pm0.15 \times 10^{15}M_{\odot}/h$). Abell 1703 lies above the standard $c$--$M$ relation predicted for the standard $\Lambda$CDM model, similar to other massive relaxed clusters with accurately determined lensing-based profiles.

We derive an accurate mass distribution of the rich galaxy cluster Cl0024+1654 (z=0.395) based on deep Subaru BR_{c}z' imaging and our recent comprehensive strong lensing analysis of HST/ACS/NIC3 observations. We obtain the weak lensing distortion and magnification of undilted samples of red and blue background galaxies by carefully combining all color and positional information. Unlike previous work, the weak and strong lensing are in excellent agreement where the data overlap. The joint mass profile continuously steepens out to the virial radius with only a minor contribution \sim 10% in the mass from known subcluster at a projected distance of \sim 700kpc/h. The projected mass distribution for the entire cluster is well fitted with a single Navarro-Frenk-White model with a virial mass, M_{vir} = (1.2 \pm 0.2) \times 10^{15} M_{sun}/h, and a concentration, c_{vir} = 9.2^{+1.4}_{-1.2}. This model fit is fully consistent with the depletion of the red background counts, providing independent confirmation. Careful examination and interpretation of X-ray and dynamical data strongly suggest that this cluster system is in a post collision state, which we show is consistent with our well-defined mass profile for a major merger occurring along the line of sight, viewed approximately 2-3Gyr after impact when the gravitational potential has had time to relax in the center, before the gas has recovered and before the outskirts are fully virialized. Finally, our full lensing analysis provides a model-independent constraint of M_{2D}(<r_{vir}) = (1.4 \pm 0.3) \times 10^{15} M_{sun}/h for the projected mass of the whole system, including any currently unbound material beyond the virial radius, which can constrain the sum of the two pre-merger cluster masses when designing simulations to explore this system.

Weak lensing work can be badly compromised by unlensed foreground and cluster members which dilute the true lensing signal. We show how the lensing amplitude in multi-colour space can be harnessed to securely separate cluster members from the foreground and background populations for three massive clusters, A1703 (z=0.258), A370 (z=0.375) and RXJ1347-11 (z=0.451) imaged with Subaru. The luminosity functions of these clusters when corrected for dilution, show similar faint-end slopes, \alpha ~= -1.0, with no marked faint-end upturn to our limit of M_R ~= -15.0, and only a mild radial gradient. In each case, the radial profile of the M/L ratio peaks at intermediate radius, ~=0.2r_{vir}, at a level of 300-500(M/L_R)_\odot, and then falls steadily towards ~100(M/L_R)_{\odot} at the virial radius, similar to the mean field level. This behaviour is likely due to the relative paucity of central late-type galaxies, whereas for the E/S0-sequence only a mild radial decline in M/L is found for each cluster. We discuss this behaviour in the context of detailed simulations where predictions for tidal stripping may now be tested accurately with observations.

We present an improved strong-lensing analysis of Cl0024+1654 ($z$=0.39) using deep HST/ACS/NIC3 images, based on 33 multiply-lensed images of 11 background galaxies. These are found with a model that assumes mass approximately traces light, with a low order expansion to allow for flexibility on large scales. The model is constrained initially by the well known 5-image system ($z$=1.675) and refined as new multiply-lensed systems are identified using the model. Photometric redshifts of these new systems are then used to constrain better the mass profile by adopting the standard cosmological relation between redshift and lensing distance. Our model requires only 6 free parameters to describe well all positional and redshift data. The resulting inner mass profile has a slope of $d\log M/d\log r\simeq -0.55$, consistent with new weak-lensing measurements where the data overlap, at $r\simeq200$ kpc/$h_{70}$. The combined profile is well fitted by a high concentration NFW mass profile, $C_{\rm vir}\sim 8.6\pm1.6$, similar to other well studied clusters, but larger than predicted with standard $\Lambda$CDM. A well defined radial critical curve is generated by the model and is clearly observed at $r \simeq 12\arcsec$, outlined by elongated images pointing towards the centre of mass. The relative fluxes of the multiply-lensed images are found to agree well with the modelled magnifications, providing an independent consistency check.

We present a multiwavelength analysis of a sample of four hot (T_X>8keV) X-ray galaxy clusters (A1689, A2261, A2142, and A2390) using joint AMiBA Sunyaev-Zel'dovich effect (SZE) and Subaru weak lensing observations, combined with published X-ray temperatures, to examine the distribution of mass and the intracluster medium (ICM) in massive cluster environments. Our observations show that A2261 is very similar to A1689 in terms of lensing properties. Many tangential arcs are visible around A2261, with an effective Einstein radius \sim 40 arcsec (at z \sim 1.5), which when combined with our weak lensing measurements implies a mass profile well fitted by an NFW model with a high concentration c_{vir} \sim 10, similar to A1689 and to other massive clusters. The cluster A2142 shows complex mass substructure, and displays a shallower profile (c_{vir} \sim 5), consistent with detailed X-ray observations which imply recent interaction. The AMiBA map of A2142 exhibits an SZE feature associated with mass substructure lying ahead of the sharp north-west edge of the X-ray core suggesting a pressure increase in the ICM. For A2390 we obtain highly elliptical mass and ICM distributions at all radii, consistent with other X-ray and strong lensing work. Our cluster gas fraction measurements, free from the hydrostatic equilibrium assumption, are overall in good agreement with published X-ray and SZE observations, with the sample-averaged gas fraction of <f_{gas}> = 0.133 \pm 0.027, for our sample <M_{vir}> = (1.2 \pm 0.1) \times 10^{15} M_{sun} h^{-1}. When compared to the cosmic baryon fraction f_b = \Omega_b/\Omega_m constrained by the WMAP 5-year data, this indicates <f_{gas}>/f_b = 0.78 \pm 0.16, i.e., (22 \pm 16)% of the baryons are missing from the hot phase of clusters.

We derive lens distortion and magnification profiles of four well known clusters observed with Subaru. Each cluster is very well fitted by the general form predicted for Cold Dark Matter (CDM) dominated halos, with good consistency found between the independent distortion and magnification measurements. The inferred level of mass concentration is surprisingly high, 8<c_{vir}<15 (<c_{vir}> = 10.4 \pm 0.9), compared to the relatively shallow profiles predicted by the Lambda CDM model, c_{vir}=5.1 \pm 1.1 (for <M_{vir}> =1.25\times 10^{15}M_{\odot}/h). This represents a 4sigma discrepancy, and includes the relatively modest effects of projection bias and profile evolution derived from N-body simulations, which oppose each other with little residual effect. In the context of CDM based cosmologies, this discrepancy implies clusters collapse earlier (z\geq 1) than predicted (z<0.5), when the Universe was correspondingly denser.

The E/SO sequence of a cluster defines a boundary redward of which a reliable weak lensing signal can be obtained from background galaxies, uncontaminated by cluster members. For bluer colors, both background and cluster members are present, reducing the distortion signal by the proportion of unlensed cluster members. In deep Subaru and HST/ACS images of A1689 the tangential distortion of galaxies with bluer colors falls rapidly toward the cluster center relative to the lensing signal of the red background. We use this dilution effect to derive the cluster light profile and luminosity function to large radius, with the advantage that no subtraction of far-field background counts is required. The light profile declines smoothly to the limit of the data, r<2Mpc/h, with a constant slope, dlog(L)/dlog(r)=-1.12+-0.06, unlike the lensing mass profile which steepens continuously with radius, so that M/L peaks at an intermediate radius, ~100kpc/h. A flatter behavior is found for the more physically meaningful ratio of dark-matter to stellar-matter, when accounting for the color-mass relation of cluster members. The cluster luminosity function has a flat slope, alpha=-1.05+-0.07, independent of radius and with no faint upturn to M_i'<-12. We establish that the very bluest objects are negligibly contaminated by the cluster V-i'<0.2, because their distortion profile rises towards the center following the red background, but offset higher by ~20%. This larger amplitude is consistent with the greater estimated depth of the faint blue galaxies, z~=2.0 compared to z~=0.85 for the red background, a purely geometric effect related to cosmological parameters. Finally, we improve upon our earlier mass profile by combining both the red and blue background populations, clearly excluding low concentration CDM profiles.