The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) will soon start to scan thousands of square degrees of the northern extragalactic sky with a unique set of $56$ optical filters from a dedicated $2.55$m telescope, JST, at the Javalambre Astrophysical Observatory. Before the arrival of the final instrument (a 1.2 Gpixels, 4.2deg$^2$ field-of-view camera), the JST was equipped with an interim camera (JPAS-Pathfinder), composed of one CCD with a 0.3deg$^2$ field-of-view and resolution of 0.23 arcsec pixel$^{-1}$. To demonstrate the scientific potential of J-PAS, with the JPAS-Pathfinder camera we carried out a survey on the AEGIS field (along the Extended Groth Strip), dubbed miniJPAS. We observed a total of $\sim 1$ deg$^2$, with the $56$ J-PAS filters, which include $54$ narrow band (NB, $\rm{FWHM} \sim 145$Angstrom) and two broader filters extending to the UV and the near-infrared, complemented by the $u,g,r,i$ SDSS broad band (BB) filters. In this paper we present the miniJPAS data set, the details of the catalogues and data access, and illustrate the scientific potential of our multi-band data. The data surpass the target depths originally planned for J-PAS, reaching $\rm{mag}_{\rm {AB}}$ between $\sim 22$ and $23.5$ for the NB filters and up to $24$ for the BB filters ($5\sigma$ in a $3$~arcsec aperture). The miniJPAS primary catalogue contains more than $64,000$ sources extracted in the $r$ detection band with forced photometry in all other bands. We estimate the catalogue to be complete up to $r=23.6$ for point-like sources and up to $r=22.7$ for extended sources. Photometric redshifts reach subpercent precision for all sources up to $r=22.5$, and a precision of $\sim 0.3$% for about half of the sample. (Abridged)

Weak gravitational lensing of background galaxies provides a direct probe of the projected matter distribution in and around galaxy clusters. Here we present a self-contained pedagogical review of cluster--galaxy weak lensing, covering a range of topics relevant to its cosmological and astrophysical applications. We begin by reviewing the theoretical foundations of gravitational lensing from the first principles, with special attention to the basics and advanced techniques of weak gravitational lensing. We summarize and discuss key findings from recent cluster--galaxy weak-lensing studies on both observational and theoretical grounds, with a focus on cluster mass profiles, the concentration--mass relation, the splashback radius, and implications from extensive mass calibration efforts for cluster cosmology.

We perform a self-calibration of the richness-to-mass ($N$-$M$) relation of CAMIRA galaxy clusters with richness $N\geq15$ at redshift $0.2\leq z<1.1$ by modeling redshift-space two-point correlation functions, namely, the auto-correlation function $\xi_{\mathrm{cc}}$ of CAMIRA clusters, the auto-correlation function $\xi_{\mathrm{gg}}$ of the CMASS galaxies spectroscopically observed in the BOSS survey, and the cross-correlation function $\xi_{\mathrm{cg}}$ between these two samples. We focus on constraining the normalization $A_{\mathrm{N}}$ of the $N$-$M$ relation in a forward-modeling approach, carefully accounting for the redshift-space distortion, the Finger-of-God effect, and the uncertainty in photometric redshifts of CAMIRA clusters. The modeling also takes into account the projection effect on the halo bias of CAMIRA clusters. The parameter constraints are shown to be unbiased according to validation tests using a large set of mock catalogs constructed from N-body simulations. At the pivotal mass $M_{500}=10^{14}h^{-1}M_{\odot}$ and the pivotal redshift $z_{\mathrm{piv}} = 0.6$, the resulting normalization $A_{\mathrm{N}}$ is constrained as $13.8^{+5.8}_{-4.2}$, $13.2^{+3.4}_{-2.7}$, and $11.9^{+3.0}_{-1.9}$ by modeling $\xi_{\mathrm{cc}}$, $\xi_{\mathrm{cc}}+\xi_{\mathrm{cg}}$, and $\xi_{\mathrm{cc}} + \xi_{\mathrm{cg}} + \xi_{\mathrm{gg}}$, with average uncertainties at levels of $36\%$, $23\%$, and $21\%$, respectively. We find that the resulting $A_{\mathrm{N}}$ is statistically consistent with those independently obtained from weak-lensing magnification and from a joint analysis of shear and cluster abundance, with a preference for a lower value at a level of $\lesssim1.9\sigma$. This implies that the absolute mass scale of CAMIRA clusters inferred from clustering is mildly higher than those from the independent methods. [abridged]

Green valley galaxies represent the population that is likely to transition from the star-forming to the quiescent phases. To investigate the role of the environment in quenching star formation, we use the wide-field data from the Hyper Suprime-Cam Strategic Subaru Proposal survey to quantify the frequency of green valley galaxies in different environments and their redshift evolution. We find that the green valley fraction, in general, is less than 20% in any redshift and environment. The green valley fraction, when normalized to the total population, is higher in the field than that in groups or clusters and decreases with a decreasing redshift and increasing mass. The lower fraction of transitional galaxies in denser environments could be a consequence of the lack of star-forming galaxies, which could be the progenitors of green valley galaxies. To assess the effect of the environment on star formation quenching, we define the effective green valley fraction as the ratio of the number of green valley galaxies to that of nonquiescent galaxies only. The effective green valley fraction for field galaxies is lower than that for group or cluster galaxies, which reveals a strong positive mass dependence and mild redshift evolution. Moreover, the specific star formation rate (sSFR) is reduced by 0.1-0.3 dex in groups or clusters. Our results thus imply that an ongoing slow quenching process has been acting in the dense environment since z~1.

Extensive surveys with the \textit{Hubble Space Telescope} (HST) over the past decade, targeting some of the most massive clusters in the sky, have uncovered dozens of galaxy-cluster strong lenses. The massive cluster strong-lens scale is typically $\theta_{E}\sim10\arcsec$ to $\sim30-35\arcsec$, with only a handful of clusters known with Einstein radii $\theta_{E}\sim40\arcsec$ or above (for $z_{source}=2$, nominally). Here we report another very large cluster lens, RXC J0032.1+1808 ($z=0.3956$), the second richest cluster in the redMapper cluster catalog and the 85th most massive cluster in the Planck Sunyaev-Zel'dovich catalog. With our Light-Traces-Mass and fully parametric (dPIEeNFW) approaches, we construct strong lensing models based on 18 multiple images of 5 background galaxies newly identified in the \textit{Hubble} data mainly from the \textit{Reionization Lensing Cluster Survey} (RELICS), in addition to a known sextuply imaged system in this cluster. Furthermore, we compare these models to Lenstool and GLAFIC models that were produced independently as part of the RELICS program. All models reveal a large effective Einstein radius of $\theta_{E}\simeq40\arcsec$ ($z_{source}=2$), owing to the obvious concentration of substructures near the cluster center. Although RXC J0032.1+1808 has a very large critical area and high lensing strength, only three magnified high-redshift candidates are found within the field targeted by RELICS. Nevertheless, we expect many more high-redshift candidates will be seen in wider and deeper observations with \textit{Hubble} or \emph{JWST}. Finally, the comparison between several algorithms demonstrates that the total error budget is largely dominated by systematic uncertainties.

We study the dynamics of Abell 370 (A370), a highly massive Hubble Frontier Fields galaxy cluster, using self-consistent three-dimensional N-body/hydrodynamical simulations. Our simulations are constrained by X-ray, optical spectroscopic and gravitational lensing, and Sunyaev-Zel'dovich (SZ) effect observations. Analyzing archival Chandra observations of A370 and comparing the X-ray morphology to the latest gravitational lensing mass reconstruction, we find offsets of ~30 kpc and ~100 kpc between the two X-ray surface brightness peaks and their nearest mass surface density peaks, suggesting that it is a merging system, in agreement with previous studies. Based on our dedicated binary cluster merger simulations, we find that initial conditions of the two progenitors with virial masses of 1.7 x 10^(15) M_sun and 1.6 x 10^(15) M_sun, an infall velocity of 3500 km/s, and an impact parameter of 100 kpc can explain the positions and the offsets between the peaks of the X-ray emission and mass surface density, the amplitude of the integrated SZ signal, and the observed relative line-of-sight velocity. Our simulations suggest that A370 is a major merger after the second core passage in the infalling phase, just before the third core passage. In this phase, the gas has not settled down in the gravitational potential well of the cluster, which explains why A370 does not follow closely the galaxy cluster scaling relations.

The radial acceleration relation (RAR) in galaxies describes a tight empirical scaling law between the total acceleration $g_\mathrm{tot}(r)=GM_\mathrm{tot}(<r)/r^2$ observed in galaxies and that expected from their baryonic mass $g_\mathrm{bar}(r)=GM_\mathrm{bar}(<r)/r^2$, with a characteristic acceleration scale of $g_\dagger\simeq 1.2\times 10^{-10}$ms$^{-2}$. Here, we examine if such a correlation exists in galaxy clusters using weak-lensing, strong-lensing, and X-ray data sets available for 20 high-mass clusters targeted by the CLASH survey. By combining our CLASH data with stellar mass estimates for the brightest cluster galaxies (BCGs) and accounting for the stellar baryonic component in clusters, we determine, for the first time, an RAR on BCG--cluster scales. The resulting RAR is well described by a tight power-law relation, $g_\mathrm{tot}\propto g_\mathrm{bar}^{0.51^{+0.04}_{-0.05}}$, with lognormal intrinsic scatter of $14.7^{+2.9}_{-2.8}\%$. The slope is consistent with the low acceleration limit of the RAR in galaxies, $g_\mathrm{tot}=\sqrt{g_\dagger\,g_\mathrm{bar}}$, whereas the intercept implies a much higher acceleration scale of $g_\ddagger = (2.02\pm0.11)\times 10^{-9}$ms$^{-2}$, indicating that there is no universal RAR that holds on all scales from galaxies to clusters. We find that the observed RAR in CLASH clusters is consistent with predictions from a semi-analytical model developed in the standard $\Lambda$CDM framework. Our results also predict the presence of a baryonic Faber--Jackson relation ($\sigma_v^4\propto M_\mathrm{bar}$) on cluster scales.

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.

In this paper, we present resolved observations of the Sunyaev-Zel'dovich (SZ) effect, obtained with the NIKA2 camera, towards the cluster of galaxies XLSSC102, a relatively low-mass system ($M_{500} \sim 2 \times 10^{14}$ M$_{\odot}$) at $z = 0.97$ detected from the XXL survey. We combine NIKA2 SZ data, XMM-Newton X-ray data, and Megacam optical data to explore, respectively, the spatial distribution of the gas electron pressure, the gas density, and the galaxies themselves. We find significant offsets between the X-ray peak, the SZ peak, the brightest cluster galaxy, and the peak of galaxy density. Additionally, the galaxy distribution and the gas present elongated morphologies. This is interpreted as the sign of a recent major merging event, which induced a local boost of the gas pressure towards the north of XLSSC102 and stripped the gas out of the galaxy group. The NIKA2 data are also combined with XXL data to construct the thermodynamic profiles of XLSSC102, obtaining relatively tight constraints up to about $\sim r_{500}$, and revealing properties that are typical of disturbed systems. We also explore the impact of the cluster centre definition and the implication of local pressure substructure on the recovered profiles. Finally, we derive the global properties of XLSSC102 and compare them to those of high-mass-and-low-redshift systems, finding no strong evidence for non-standard evolution. We also use scaling relations to obtain alternative mass estimates from our profiles. The variation between these different mass estimates reflects the difficulty to accurately measure the mass of low-mass clusters at z$\sim$1, especially with low signal-to-noise ratio (S/N) data and for a disturbed system. However, it also highlights the strength of resolved SZ observations alone and in combination with survey-like X-ray data.

Context: The shape of the mass density profiles of cosmological halos informs us of the nature of DM and DM-baryons interactions. Previous estimates of the inner slope of the mass density profiles of clusters of galaxies are in opposition to predictions derived from numerical simulations of CDM. Aims: We determine the inner slope of the DM density profile of a massive cluster of galaxies, Abell S1063 (RXC J2248.7$-$4431) at $z=0.35$, with a dynamical analysis based on an extensive spectroscopic campaign carried out with the VIMOS and MUSE spectrographs at the ESO VLT. This new data set provides an unprecedented sample of 1234 spectroscopic members, 104 of which are located in the cluster core (R < 200 kpc), extracted from the MUSE integral field spectroscopy. The latter also allows the stellar velocity dispersion profile of the brightest cluster galaxy (BCG) to be measured out to 40 kpc. Methods: We used an upgraded version of the MAMPOSSt technique to perform a joint maximum likelihood fit to the velocity dispersion profile of the BCG and to the velocity distribution of cluster member galaxies over a radial range from 1 kpc to the virial radius (~ 2.7 Mpc). Results: We find a value of $\gamma_{DM} =0.99 \pm 0.04$ for the inner logarithmic slope of the DM density profile after marginalizing over all the other parameters. The newly determined dynamical mass profile is found to be in excellent agreement with the mass density profiles obtained from the independent X-ray hydrostatic analysis based on deep Chandra data, as well as the strong and weak lensing analyses. Our value of \gamma_{DM} is in very good agreement with predictions from cosmological CDM simulations. We will extend our analysis to more clusters in future works. If confirmed on a larger cluster sample, our result makes this DM model more appealing than alternative models.

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 a systematic study of gas density perturbations in cool cores of high-mass galaxy clusters. We select 12 relaxed clusters from the Cluster Lensing And Supernova survey with Hubble (CLASH) sample and analyze their cool core features observed with the Chandra X-ray Observatory. We focus on the X-ray residual image characteristics after subtracting their global profile of the X-ray surface brightness distribution. We find that all the galaxy clusters in our sample have, at least, both one positive and one negative excess regions in the X-ray residual image, indicating the presence of gas density perturbations. We identify and characterize the locally perturbed regions using our detection algorithm, and extract X-ray spectra of the intracluster medium (ICM). The ICM temperature in the positive excess region is lower than that in the negative excess region, whereas the ICM in both regions is in pressure equilibrium in a systematic manner. These results indicate that gas sloshing in cool cores takes place in more than 80% of relaxed clusters (95% CL). We confirm this physical picture by analyzing synthetic X-ray observations of a cool-core cluster from a hydrodynamic simulation, finding that our detection algorithm can accurately extract both the positive and negative excess regions and can reproduce the temperature difference between them. Our findings support the picture that the gas density perturbations are induced by gas sloshing, and a large fraction of cool-core clusters have undergone gas sloshing, indicating that gas sloshing may be capable of suppressing runaway cooling of the ICM.

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

Scaling relations trace the formation and evolution of galaxy clusters. We exploited multi-wavelength surveys -- the XXL survey at \emph{XMM-Newton} in the X-ray band, and the Hyper Suprime-Cam Subaru Strategic Program for optical weak lensing -- to study an X-ray selected, complete sample of clusters and groups. The scalings of gas mass, temperature, and soft-band X-ray luminosity with the weak lensing mass show imprints of radiative cooling and AGN feedback in groups. From the multi-variate analysis, we found some evidence for steeper than self-similar slopes for gas mass ($\beta_{m_\text{g}|m}=1.73 \pm0.80$) and luminosity ($\beta_{l|m}=1.91\pm0.94$) and a nearly self-similar slope for the temperature ($\beta_{t|m}=0.78\pm0.43$). Intrinsic scatters of X-ray properties appear to be positively correlated at a fixed mass (median correlation factor $\rho_{X_1X_2|m}\sim0.34$) due to dynamical state and merger history of the halos. Positive correlations with the weak lensing mass (median correlation factor $\rho_{m_\text{wl}X|m}\sim0.35$) can be connected to triaxiality and orientation. Comparison of weak lensing and hydrostatic masses suggests a small role played by non-thermal pressure support ($9\pm17\%$).

We present 90, 140, and 268GHz sub-arcminute resolution imaging of the Sunyaev-Zel'dovich effect (SZE) in MACSJ0717.5+3745. Our 90GHz SZE data result in a sensitive, 34uJy/bm map at 13" resolution using MUSTANG. Our 140 and 268GHz SZE imaging, with resolutions of 58" and 31" and sensitivities of 1.8 and 3.3mJy/beam respectively, was obtained using Bolocam. We compare these maps to a 2-dimensional pressure map derived from Chandra X-ray observations. Our MUSTANG data confirm previous indications from Chandra of a pressure enhancement due to shock-heated, >20keV gas immediately adjacent to extended radio emission seen in low-frequency radio maps. The MUSTANG data also detect pressure substructure that is not well-constrained by the X-ray data in the remnant core of a merging subcluster. We find that the small-scale pressure enhancements in the MUSTANG data amount to ~2% of the total pressure measured in the 140GHz Bolocam observations. The X-ray template also fails on larger scales to accurately describe the Bolocam data, particularly at the location of a subcluster known to have a high line of sight optical velocity (~3200km/s). Our Bolocam data are adequately described when we add an additional component - not described by a thermal SZE spectrum - coincident with this subcluster. Using flux densities extracted from our model fits, and marginalizing over the temperature constraints for the region, we fit a thermal+kinetic SZE spectrum to our data and find the subcluster has a best-fit line of sight proper velocity of 3600+3440/-2160km/s. This agrees with the optical velocity estimates for the subcluster. The probability of velocity<0 given our measurements is 2.1%. Repeating this analysis using flux densities measured non-parametrically results in a 3.4% probability of a velocity<=0. We note that this tantalizing result for the kinetic SZE is on resolved, subcluster scales.

We carry out a detailed strong lensing analysis of a sub-sample of eight galaxy clusters of the CLASH survey, in the redshift range of $ z_{\rm cluster} = [0.23-0.59]$, using extensive spectroscopic information, primarily MUSE archival data complemented with CLASH-VLT redshift measurements. Different models are tested in each cluster depending on the complexity of its mass distribution and on the number of detected multiple images. Four clusters show more than five spectroscopically confirmed multiple image families. In this sample, we do not make use of families that are only photometrically identified, in order to reduce model degeneracies and systematics due to the potential misidentifications of some multiple images. We present spectroscopic confirmation of 27 multiply lensed sources, with no previous spectroscopic measurements, spanning over the redshift range of $z_{\rm src}=[0.7-6.1]$. Moreover, we confirm an average of $48$ galaxy members in the core of each cluster, thanks to the high efficiency and large field of view of MUSE. Despite having different properties (i.e., number of mass components, total mass, redshift, etc), the projected total mass and mass density profiles of all clusters have very similar shapes, when rescaled by independent measurements of $M_{200c}$ and $R_{200c}$. Specifically, we measure the mean value of the projected total mass of our cluster sample within 10 (20)% of $R_{200c}$ to be 0.13 (0.32) of $M_{200c}$, with a remarkably small scatter of 5 (6)%. Furthermore, the large number of high-z sources and the precise magnification maps derived in this work for four clusters add up to the sample of high-quality gravitational telescopes to be used to study the faint and distant Universe. The strong lensing models and the full redshift catalogues from MUSE are publicly available.

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.

Strong gravitational lensing by clusters of galaxies probes the mass distribution at the core of each cluster and magnifies the universe behind it. MACS J0417.5-1154 at z=0.443 is one of the most massive clusters known based on weak lensing, X-ray, and Sunyaev-Zel'dovich analyses. Here we compute a strong lens model of MACS J0417 based on Hubble Space Telescope imaging observations collected, in part, by the Reionization Lensing Cluster Survey (RELICS), and recently reported spectroscopic redshifts from the MUSE instrument on the Very Large Telescope (VLT). We measure an Einstein radius of ThetaE=36'' at z = 9 and a mass projected within 200 kpc of M(200 kpc) = 1.78+0.01-0.03x10**14Msol. Using this model, we measure a ratio between the mass attributed to cluster-member galaxy halos and the main cluster halo of order 1:100. We assess the probability to detect magnified high-redshift galaxies in the field of this cluster, both for comparison with RELICS HST results and as a prediction for the James Webb Space Telescope (JWST) Guaranteed Time Observations upcoming for this cluster. Our lensing analysis indicates that this cluster has similar lensing strength to other clusters in the RELICS program. Our lensing analysis predicts a detection of at least a few z~6-8 galaxies behind this cluster, at odds with a recent analysis that yielded no such candidates in this field. Reliable strong lensing models are crucial for accurately predicting the intrinsic properties of lensed galaxies. As part of the RELICS program, our strong lensing model produced with the Lenstool parametric method is publicly available through the Mikulski Archive for Space Telescopes (MAST).

Large surveys of galaxy clusters with the Hubble and Spitzer Space Telescopes, including CLASH and the Frontier Fields, have demonstrated the power of strong gravitational lensing to efficiently deliver large samples of high-redshift galaxies. We extend this strategy through a wider, shallower survey named RELICS, the Reionization Lensing Cluster Survey. This survey, described here, was designed primarily to deliver the best and brightest high-redshift candidates from the first billion years after the Big Bang. RELICS observed 41 massive galaxy clusters with Hubble and Spitzer at 0.4-1.7um and 3.0-5.0um, respectively. We selected 21 clusters based on Planck PSZ2 mass estimates and the other 20 based on observed or inferred lensing strength. Our 188-orbit Hubble Treasury Program obtained the first high-resolution near-infrared images of these clusters to efficiently search for lensed high-redshift galaxies. We observed 46 WFC3/IR pointings (~200 arcmin^2) with two orbits divided among four filters (F105W, F125W, F140W, and F160W) and ACS imaging as needed to achieve single-orbit depth in each of three filters (F435W, F606W, and F814W). As previously reported by Salmon et al., we discovered 322 z ~ 6 - 10 candidates, including the brightest known at z ~ 6, and the most distant spatially-resolved lensed arc known at z ~ 10. Spitzer IRAC imaging (945 hours awarded, plus 100 archival) has crucially enabled us to distinguish z ~ 10 candidates from z ~ 2 interlopers. For each cluster, two HST observing epochs were staggered by about a month, enabling us to discover 11 supernovae, including 3 lensed supernovae, which we followed up with 20 orbits from our program. We delivered reduced HST images and catalogs of all clusters to the public via MAST and reduced Spitzer images via IRSA. We have also begun delivering lens models of all clusters, to be completed before the JWST GO call for proposals.

Strong gravitational lensing (SL) is a powerful means to map the distribution of dark matter. In this work, we perform a SL analysis of the prominent X-ray cluster RXJ0152.7-1357 (z=0.83, also known as CL 0152.7-1357) in \textit{Hubble Space Telescope} images, taken in the framework of the Reionization Lensing Cluster Survey (RELICS). On top of a previously known $z=3.93$ galaxy multiply imaged by RXJ0152.7-1357, for which we identify an additional multiple image, guided by a light-traces-mass approach we identify seven new sets of multiply imaged background sources lensed by this cluster, spanning the redshift range [1.79-3.93]. A total of 25 multiple images are seen over a small area of ~0.4 $arcmin^2$, allowing us to put relatively high-resolution constraints on the inner matter distribution. Although modestly massive, the high degree of substructure together with its very elongated shape make RXJ0152.7-1357 a very efficient lens for its size. This cluster also comprises the third-largest sample of z~6-7 candidates in the RELICS survey. Finally, we present a comparison of our resulting mass distribution and magnification estimates with those from a Lenstool model. These models are made publicly available through the MAST archive.

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 line-of-sight gas sloshing first found in a cool core in a galaxy cluster. The galaxy cluster Abell 907 is identified as a relaxed cluster owing to its global X-ray surface brightness taken by the Chandra X-ray Observatory. The X-ray residual image after removing the global emission of the intracluster medium (ICM), however, shows an arc-like positive excess and a negative excess surrounding the central positive excess in the cluster core, which in turn indicates a disturbance of the ICM. We analyze the X-ray spectra extracted from both regions and find that (1) the ICM temperature and the metal abundance in the positive excess are lower and higher than those in the negative excess, respectively, and (2) the ICM is nearly in pressure equilibrium. We also find a slight redshift difference between the positive and the negative excesses, which corresponds to the velocity shear of $1680^{+1300}_{-920}$ km s$^{-1}$ ($1\sigma$). The X-ray residual image and the ICM properties are consistent with those expected by line-of-sight gas sloshing. Assuming that the gas is moving toward inverse-parallel to each other along the line-of-sight, the shear velocity is expected to be $\sim 800$ km s$^{-1}$. The velocity field of this level is able to provide non-thermal pressure support by $\sim 34\%$ relative to the thermal one. The total kinetic energy inferred from the shear velocity corresponds to $\sim 30\%$ of the bolometric luminosity of the sloshing ICM. Abell 907 is therefore complementary to galaxy clusters in which gas sloshing takes place in the plane of the sky, and is important for understanding gas dynamics driven by sloshing and its influence on the heating to prevent runaway cooling.

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.

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

To investigate the dynamics of the galaxy cluster A2142, we compile an extended catalog of 2239 spectroscopic redshifts of sources, including newly measured 237 redshifts, within 30 arcmin from the cluster center. With the $\sigma$-plateau algorithm from the caustic method, we identify 868 members and a number of substructures in the galaxy distribution both in the outskirts, out to $\sim$3.5 Mpc from the cluster center, and in the central region. In the outskirts, one substructure overlaps a falling clump of gas previously identified in the X-ray band. These substructures suggests the presence of multiple minor mergers, which are responsible for the complex dynamics of A2142, and the absence of recent or ongoing major mergers. We show that the distribution of the galaxies in the cluster core and in several substructures are consistent with the mass distribution inferred from the weak lensing signal. Moreover, we use spatially-resolved X-ray spectroscopy to measure the redshift of different regions of the intracluster medium within $\sim$3 arcmin from the cluster center. We find a ring of gas near the two X-ray cold fronts identified in previous analyses and measure a velocity of this ring of $810 \pm 330 {\rm km s^{-1}}$ larger than the cluster mean velocity. Our analysis suggests the presence of another ring surrounding the core, whose velocity is $660 \pm 300 {\rm km s^{-1}}$ larger than the cluster velocity. These X-ray features are not associated to any optical substructures, and support the core-sloshing scenario suggested in previous work.

We investigate the splashback features of dark-matter halos based on cosmic density and velocity fields. Besides the density correlation function binned by the halo orientation angle which was used in the literature, we introduce, for the first time, the corresponding velocity statistic, alignment momentum correlation function, to take into account the asphericity of halos. Using large-volume, high-resolution N-body simulations, we measure the alignment statistics of density and velocity. On halo scales, $x\sim R_\mathrm{200m} \sim 1Mpc/h$, we detect a sharp steepening in the momentum correlation associated with the physical halo boundary, or the splashback feature, which is found more prominent than in the density correlation. We also find that the splashback radius determined from the density correlation becomes $\sim 3.5\%$ smaller than that from the momentum correlation, with their correlation coefficient being 0.605. Moreover, the orientation-dependent splashback feature due to halo asphericity is measured when the density profile is determined by dark-matter particles, which can be used as a test of collisional CDM since the halo shape is predicted to be rounder in such a model.

We present a strong-lensing analysis of four massive galaxy clusters imaged with the Hubble Space Telescope in the Reionization Lensing Cluster Survey. We use a Light-Traces-Mass technique to uncover sets of multiply images and constrain the mass distribution of the clusters. These mass models are the first published for Abell S295 and MACS J0159.8-0849, and are improvements over previous models for Abell 697 and MACS J0025.4-1222. Our analysis for MACS J0025.4-1222 and Abell S295 shows a bimodal mass distribution supporting the merger scenarios proposed for these clusters. The updated model for MACS J0025.4-1222 suggests a substantially smaller critical area than previously estimated. For MACS J0159.8-0849 and Abell 697 we find a single peak and relatively regular morphology, revealing fairly relaxed clusters. Despite being less prominent lenses, three of these clusters seem to have lensing strengths, i.e. cumulative area above certain magnification, similar to the Hubble Frontier Fields clusters (e.g., A($\mu>5$) $\sim 1-3$ arcmin$^2$, A($\mu>10$) $\sim 0.5-1.5$ arcmin$^2$), which in part can be attributed to their merging configurations. We make our lens models publicly available through the Mikulski Archive for Space Telescopes. Finally, using Gemini-N/GMOS spectroscopic observations we detect a single emission line from a high-redshift $J_{125}\simeq25.7$ galaxy candidate lensed by Abell 697. While we cannot rule out a lower-redshift solution, we interpret the line as Ly$\alpha$ at $z=5.800\pm 0.001$, in agreement with its photometric redshift and dropout nature. Within this scenario we measure a Ly$\alpha$ rest-frame equivalent width of $52\pm22$ \AA, and an observed Gaussian width of $117\pm 15$ km/s.

We present a lens model for the cluster SPT-CLJ0615$-$5746, which is the highest redshift ($z=0.972$) system in the Reionization of Lensing Clusters Survey (RELICS), making it the highest redshift cluster for which a full strong lens model is published. We identify three systems of multiply-imaged lensed galaxies, two of which we spectroscopically confirm at $z=1.358$ and $z=4.013$, which we use as constraints for the model. We find a foreground structure at $z\sim0.4$, which we include as a second cluster-sized halo in one of our models; however two different statistical tests find the best-fit model consists of one cluster-sized halo combined with three individually optimized galaxy-sized halos, as well as contributions from the cluster galaxies themselves. We find the total projected mass density within $r=26.7"$ (the region where the strong lensing constraints exist) to be $M=2.51^{+0.15}_{-0.09}\times 10^{14}$~M$_{\odot}$. If we extrapolate out to $r_{500}$, our projected mass density is consistent with the mass inferred from weak lensing and from the Sunyaev-Zel'dovich effect ($M\sim10^{15}$~M$_{\odot}$). This cluster is lensing a previously reported $z\sim10$ galaxy, which, if spectroscopically confirmed, will be the highest-redshift strongly lensed galaxy known.

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.

The $\Lambda$CDM model of structure formation makes strong predictions on concentration and shape of DM (dark matter) halos, which are determined by mass accretion processes. Comparison between predicted shapes and observations provides a geometric test of the $\Lambda$CDM model. Accurate and precise measurements needs a full three-dimensional analysis of the cluster mass distribution. We accomplish this with a multi-probe 3D analysis of the X-ray regular CLASH (Cluster Lensing And Supernova survey with Hubble) clusters combining strong and weak lensing, X-ray photometry and spectroscopy, and the Sunyaev-Zel'dovich effect. The cluster shapes and concentrations are consistent with $\Lambda$CDM predictions. The CLASH clusters are randomly oriented, as expected given the sample selection criteria. Shapes agree with numerical results for DM-only halos, which hints at baryonic physics being not so effective in making halos rounder.

We perform a three-dimensional triaxial analysis of 16 X-ray regular and 4 high-magnification galaxy clusters selected from the CLASH survey by combining two-dimensional weak-lensing and central strong-lensing constraints. In a Bayesian framework, we constrain the intrinsic structure and geometry of each individual cluster assuming a triaxial Navarro-Frenk-White halo with arbitrary orientations, characterized by the mass $M_{200\mathrm{c}}$, halo concentration $C_{200\mathrm{c}}$, and triaxial axis ratios ($q_{\mathrm{a}} \le q_{\mathrm{b}}$), and investigate scaling relations between these halo structural parameters. From triaxial modeling of the X-ray-selected subsample, we find that the halo concentration decreases with increasing cluster mass, with a mean concentration of $C_{200\mathrm{c}} = 4.82\pm0.30$ at the pivot mass $M_{200\mathrm{c}}=10^{15}M_{\odot}h^{-1}$. This is consistent with the result from spherical modeling, $C_{200\mathrm{c}}=4.51\pm 0.14$. Independently of the priors, the minor-to-major axis ratio $q_{\mathrm{a}}$ of our full sample exhibits a clear deviation from the spherical configuration ($q_{\mathrm{a}}=0.52 \pm 0.04$ at $10^{15}M_{\odot}h^{-1}$ with uniform priors), with a weak dependence on the cluster mass. Combining all 20 clusters, we obtain a joint ensemble constraint on the minor-to-major axis ratio of $q_{\mathrm{a}}=0.652^{+0.162}_{-0.078}$ and a lower bound on the intermediate-to-major axis ratio of $q_{\mathrm{b}}>0.63$ at the $2\sigma$ level from an analysis with uniform priors. Assuming priors on the axis ratios derived from numerical simulations, we constrain the degree of triaxiality for the full sample to be $\mathcal{T}=0.79 \pm 0.03$ at $10^{15}M_{\odot}h^{-1}$, indicating a preference for a prolate geometry of cluster halos. We find no statistical evidence for an orientation bias ($f_{\mathrm{geo}}=0.93 \pm 0.07$) (abridged)

Strong gravitational lensing by galaxy clusters has become a powerful tool for probing the high-redshift Universe, magnifying distant and faint background galaxies. Reliable strong lensing (SL) models are crucial for determining the intrinsic properties of distant, magnified sources and for constructing their luminosity function. We present here the first SL analysis of MACS J0308.9+2645 and PLCK G171.9-40.7, two massive galaxy clusters imaged with the Hubble Space Telescope in the framework of the Reionization Lensing Cluster Survey (RELICS). We use the Light-Traces-Mass modeling technique to uncover sets of multiply imaged galaxies and constrain the mass distribution of the clusters. Our SL analysis reveals that both clusters have particularly large Einstein radii ($\theta_E>30"$ for a source redshift of $z_s=2$), providing fairly large areas with high magnifications, useful for high-redshift galaxy searches ($\sim2$ arcmin$^{2}$ with $\mu>5$ to $\sim1$ arcmin$^{2}$ with $\mu>10$, similar to a typical \textit{Hubble Frontier Fields} cluster). We also find that MACS J0308.9+2645 hosts a promising, apparently bright (J$\sim23.2-24.6$ AB), multiply imaged high-redshift candidate at $z\sim6.4$. These images are amongst the brightest high-redshift candidates found in RELICS. Our mass models, including magnification maps, are made publicly available for the community through the Mikulski Archive for Space Telescopes.

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.

Massive foreground galaxy clusters magnify and distort the light of objects behind them, permitting a view into both the extremely distant and intrinsically faint galaxy populations. We present here the z ~ 6 - 8 candidate high-redshift galaxies from the Reionization Lensing Cluster Survey (RELICS), a Hubble and Spitzer Space Telescope survey of 41 massive galaxy clusters spanning an area of ~200 arcmin^2. These clusters were selected to be excellent lenses and we find similar high-redshift sample sizes and magnitude distributions as CLASH. We discover 321 candidate galaxies with photometric redshifts between z ~ 6 to z ~ 8, including extremely bright objects with H-band magnitudes of m_AB ~ 23 mag. As a sample, the observed (lensed) magnitudes of these galaxies are among the brightest known at z> 6, comparable to much wider, blank-field surveys. RELICS demonstrates the efficiency of using strong gravitational lenses to produce high-redshift samples in the epoch of reionization. These brightly observed galaxies are excellent targets for follow-up study with current and future observatories, including the James Webb Space Telescope.

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.

Strong gravitational lensing by galaxy clusters magnifies background galaxies, enhancing our ability to discover statistically significant samples of galaxies at z>6, in order to constrain the high-redshift galaxy luminosity functions. Here, we present the first five lens models out of the Reionization Lensing Cluster Survey (RELICS) Hubble Treasury Program, based on new HST WFC3/IR and ACS imaging of the clusters RXC J0142.9+4438, Abell 2537, Abell 2163, RXC J2211.7-0349, and ACT-CLJ0102-49151. The derived lensing magnification is essential for estimating the intrinsic properties of high-redshift galaxy candidates, and properly accounting for the survey volume. We report on new spectroscopic redshifts of multiply imaged lensed galaxies behind these clusters, which are used as constraints, and detail our strategy to reduce systematic uncertainties due to lack of spectroscopic information. In addition, we quantify the uncertainty on the lensing magnification due to statistical and systematic errors related to the lens modeling process, and find that in all but one cluster, the magnification is constrained to better than 20% in at least 80% of the field of view, including statistical and systematic uncertainties. The five clusters presented in this paper span the range of masses and redshifts of the clusters in the RELICS program. We find that they exhibit similar strong lensing efficiencies to the clusters targeted by the Hubble Frontier Fields within the WFC3/IR field of view. Outputs of the lens models are made available to the community through the Mikulski Archive for Space Telescopes

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$.

The most distant galaxies known are at z~10-11, observed 400-500 Myr after the Big Bang. The few z~10-11 candidates discovered to date have been exceptionally small- barely resolved, if at all, by the Hubble Space Telescope. Here we present the discovery of SPT0615-JD, a fortuitous z~10 (z_phot=9.9+/-0.6) galaxy candidate stretched into an arc over ~2.5" by the effects of strong gravitational lensing. Discovered in the Reionization Lensing Cluster Survey (RELICS) Hubble Treasury program and companion S-RELICS Spitzer program, this candidate has a lensed H-band magnitude of 25.7+/-0.1 AB mag. With a magnification of \mu~4-7 estimated from our lens models, the de-lensed intrinsic magnitude is 27.6+/-0.3 AB mag, and the half-light radius is r_e<0.8 kpc, both consistent with other z>9 candidates. The inferred stellar mass (log [M* /M_Sun]=9.7^{+0.7}_{-0.5}) and star formation rate (\log [SFR/M_Sun yr^{-1}]=1.3^{+0.2}_{-0.3}) indicate that this candidate is a typical star-forming galaxy on the z>6 SFR-M* relation. We note that three independent lens models predict two counterimages, at least one of which should be of a similar magnitude to the arc, but these counterimages are not yet detected. Counterimages would not be expected if the arc were at lower redshift. However, the only spectral energy distributions capable of fitting the Hubble and Spitzer photometry well at lower redshifts require unphysical combinations of z~2 galaxy properties. The unprecedented lensed size of this z~10 candidate offers the potential for the James Webb Space Telescope to study the geometric and kinematic properties of a galaxy observed 500 Myr after the Big Bang.

We present a new method for photometering objects in galaxy clusters. We introduce a mode-filtering technique for removing spatially variable backgrounds, improving both detection and photometric accuracy (roughly halving the scatter in the red sequence compared to previous catalogs of the same clusters). This method is based on robustly determining the distribution of background pixel values and should provide comparable improvement in photometric analysis of any crowded fields. We produce new multiwavelength catalogs for the 25 CLASH cluster fields in all 16 bandpasses from the UV through the near IR, as well as rest-frame magnitudes. A comparison with spectroscopic values from the literature finds a ~30% decrease in the redshift deviation from previously-released CLASH photometry. This improvement in redshift precision, in combination with a detection scheme designed to maximize purity, yields a substantial upgrade in cluster member identification over the previous CLASH galaxy catalog. We construct luminosity functions for each cluster, reliably reaching depths of at least 4.5 mag below M* in every case, and deeper still in several clusters. We measure M* , $\alpha$, and their redshift evolution, assuming the cluster populations are coeval, and find little to no evolution of $\alpha$, $-0.9\lesssim\langle\alpha\rangle\lesssim -0.8$, and M* values consistent with passive evolution. We present a catalog of galaxy photometry, photometric and spectroscopic redshifts, and rest-frame photometry for the full fields of view of all 25 CLASH clusters. Not only will our new photometric catalogs enable new studies of the properties of CLASH clusters, but mode-filtering techniques, such as those presented here, should greatly enhance the data quality of future photometric surveys of crowded fields.

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.

We present a new strong lensing analysis of the galaxy cluster MACS J1206.2-0847 (MACS 1206), at z=0.44, using deep spectroscopy from CLASH-VLT and VLT/MUSE archival data in combination with imaging from the Cluster Lensing and Supernova survey with Hubble. MUSE observations enable the spectroscopic identification of 23 new multiply imaged sources, extending the previous compilations by a factor of approximately five. In total, we use the positional measurements of 82 spectroscopic multiple images belonging to 27 families at z=1.0-6.1 to reconstruct the projected total mass distribution of MACS 1206. Remarkably, 11 multiple images are found within 50 kpc of the brightest cluster galaxy, making this an unprecedented set of constraints for the innermost projected mass distribution of a galaxy cluster. We thus find that, although dynamically relaxed, the smooth matter component (dark matter plus hot gas) of MACS 1206 shows a significant asymmetry, which closely follows the asymmetric distribution of the stellar component (galaxy members and intracluster light). We determine the value of the innermost logarithmic slope of the projected total mass density profile and find it to be close to the canonical Navarro-Frenk-White value. We demonstrate that this quantity is very robust against different parametrizations of the diffuse mass component; however, this is not the case when only one central image is used in the mass reconstruction. We also show that the mass density profile from our new strong lensing model is in very good agreement with dynamical and X-ray measurements at larger radii, where they overlap.

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}$.

We perform a maximum likelihood kinematic analysis of the two dynamically relaxed galaxy clusters MACS J1206.2-0847 at $z=0.44$ and RXC J2248.7-4431 at $z=0.35$ to determine the total mass profile in modified gravity models, using a modified version of the MAMPOSSt code of Mamon, Biviano and Bou\'e. Our work is based on the kinematic and lensing mass profiles derived using the data from the Cluster Lensing And Supernova survey with Hubble (hereafter CLASH) and the spectroscopic follow-up with the Very Large Telescope (hereafter CLASH-VLT). We assume a spherical Navarro-Frenk-White (NFW hereafter) profile in order to obtain a constraint on the fifth force interaction range $\lambda$ for models in which the dependence of this parameter on the enviroment is negligible at the scale considered (i.e. $\lambda=const$) and fixing the fifth force strength to the value predicted in $f(R)$ gravity. We then use information from lensing analysis to put a prior on the other NFW free parameters. In the case of MACSJ 1206 the joint kinematic+lensing analysis leads to an upper limit on the effective interaction range $\lambda\le1.61\,\mbox{Mpc}$ at $\Delta\chi^{2}=2.71$ on the marginalized distribution. For RXJ 2248 instead a possible tension with the $\Lambda$CDM model appears when adding lensing information, with a lower limit $\lambda\ge0.14\,\mbox{Mpc}$ at $\Delta\chi^{2}=2.71$. This is consequence of the slight difference between the lensing and kinematic data, appearing in GR for this cluster, that could be explained in terms of modifications of gravity. We discuss the impact of systematics and the limits of our analysis as well as future improvements of the results obtained.[...]

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/.

The XXL survey currently covers two 25 sq. deg. patches with XMM observations of ~10ks. We summarise the scientific results associated with the first release of the XXL data set, that occurred mid 2016. We review several arguments for increasing the survey depth to 40 ks during the next decade of XMM operations. X-ray (z<2) cluster, (z<4) AGN and cosmic background survey science will then benefit from an extraordinary data reservoir. This, combined with deep multi-$\lambda$ observations, will lead to solid standalone cosmological constraints and provide a wealth of information on the formation and evolution of AGN, clusters and the X-ray background. In particular, it will offer a unique opportunity to pinpoint the z>1 cluster density. It will eventually constitute a reference study and an ideal calibration field for the upcoming eROSITA and Euclid missions.

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 investigate the effects of intrinsic alignments (IA) of dark-matter halo shapes on cosmic density and velocity fields from cluster to cosmic scales beyond 100 Mpc/h. Besides the density correlation function binned by the halo orientation angle which was used in the literature, we introduce, for the first time, the corresponding two velocity statistics, the angle-binned pairwise infall momentum and momentum correlation function. Using large-volume, high-resolution N-body simulations, we measure the alignment statistics of density and velocity, both in real and redshift space. We find that the alignment signal is not amplified by redshift-space distortions at linear scales. Behaviors of IA in the velocity statistics are similar to those in the density statistics, except that the halo orientations are aligned with the velocity field up to a scale larger than those with the density field, x>100 Mpc/h. On halo scales, x~ R_{200m} ~ 1 Mpc/h, we detect a sharp steepening in the momentum correlation associated with the physical halo boundary, or the splashback feature, which is found more prominent than in the density correlation. Our results indicate that observations of IA with the velocity field can provide additional information on cosmological models from large scales and on physical sizes of halos from small scales.

We present a combined strong and weak lensing analysis of the J085007.6+360428 (J0850) field, which was selected by its high projected concentration of luminous red galaxies and contains the massive cluster Zwicky 1953. Using Subaru/Suprime-Cam $BVR_{c}I_{c}i^{\prime}z^{\prime}$ imaging and MMT/Hectospec spectroscopy, we first perform a weak lensing shear analysis to constrain the mass distribution in this field, including the cluster at $z = 0.3774$ and a smaller foreground halo at $z = 0.2713$. We then add a strong lensing constraint from a multiply-imaged galaxy in the imaging data with a photometric redshift of $z \approx 5.03$. Unlike previous cluster-scale lens analyses, our technique accounts for the full three-dimensional mass structure in the beam, including galaxies along the line of sight. In contrast with past cluster analyses that use only lensed image positions as constraints, we use the full surface brightness distribution of the images. This method predicts that the source galaxy crosses a lensing caustic such that one image is a highly-magnified "fold arc", which could be used to probe the source galaxy's structure at ultra-high spatial resolution ($< 30$ pc). We calculate the mass of the primary cluster to be $\mathrm{M_{vir}} = 2.93_{-0.65}^{+0.71} \times 10^{15}~\mathrm{M_{\odot}}$ with a concentration of $\mathrm{c_{vir}} = 3.46_{-0.59}^{+0.70}$, consistent with the mass-concentration relation of massive clusters at a similar redshift. The large mass of this cluster makes J0850 an excellent field for leveraging lensing magnification to search for high-redshift galaxies, competitive with and complementary to that of well-studied clusters such as the HST Frontier Fields.

We construct a large, redshift complete sample of distant galaxy clusters by correlating Sloan Digital Sky Survey (SDSS) Data Release 12 (DR12) redshifts with clusters identified with the red-sequence Matched-filter Probabilistic Percolation (redMaPPer) algorithm. Our spectroscopic completeness is 97% for ~ 7,000 clusters within the redMaPPer selection limit, $z \leqslant$ 0.325, so that our cluster correlation functions are much more precise than earlier work and not suppressed by photometric redshifts. We derive an accurate power-law mass-richness relation from the observed abundance with respect to the mass function from Millennium XXL (MXXL) simulation, adjusted to the Planck weighted cosmology. The number density of clusters is found to decline by 20% over the range 0.1 $< z <$ 0.3, in good agreement with the evolution predicted by MXXL. Our projected three-dimensional correlation function scales with richness, $\lambda$, rising from $r_0=$ 14 h$^{-1}$ Mpc at $\lambda\simeq$ 25, to $r_0=$ 22 h$^{-1}$ Mpc at $\lambda\simeq$ 60, with a gradient that matches MXXL when applying our mass-richness relation, whereas the observed amplitude of the correlation function at $\left<z\right>=$ 0.24 exceeds the MXXL prediction by 20% at the $\simeq$ 2.5$\sigma$ level. This tension cannot be blamed on spurious, randomly located clusters as this would reduce the correlation amplitude. Full consistency between the correlation function and the abundances is achievable for the pre-Planck values of $\sigma_8=$ 0.9, $\Omega_m=$ 0.25, and h = 0.73, matching the improved distance ladder estimate of the Hubble constant.

The lensing signal around galaxy clusters can, in principle, be used to test detailed predictions for their average mass profile from numerical simulations. However, the intrinsic shape of the profiles can be smeared out when a sample that spans a wide range of cluster masses is averaged in physical length units. This effect especially conceals rapid changes in gradient such as the steep drop associated with the splashback radius, a sharp edge corresponding to the outermost caustic in accreting halos. We optimize the extraction of such local features by scaling individual halo profiles to a number of spherical overdensity radii, and apply this method to 16 X-ray-selected high-mass clusters targeted in the Cluster Lensing And Supernova survey with Hubble. By forward-modeling the weak- and strong-lensing data presented by Umetsu et al., we show that, regardless of the scaling overdensity, the projected ensemble density profile is remarkably well described by an NFW or Einasto profile out to $R \sim 2.5h^{-1}$Mpc, beyond which the profiles flatten. We constrain the NFW concentration to $c_{200c} = 3.66 \pm 0.11$ at $M_{200c} \simeq 1.0 \times 10^{15}h^{-1}M_\odot$, consistent with and improved from previous work that used conventionally stacked lensing profiles, and in excellent agreement with theoretical expectations. Assuming the profile form of Diemer & Kravtsov and generic priors calibrated from numerical simulations, we place a lower limit on the splashback radius of the cluster halos, if it exists, of $R_{sp}/r_{200m} > 0.89$ ($R_{sp} > 1.83h^{-1}$Mpc) at 68% confidence. The corresponding density feature is most pronounced when the cluster profiles are scaled by $r_{200m}$, and smeared out when scaled to higher overdensities.

Multi-wavelength techniques can probe the distribution and the physical properties of baryons and dark matter in galaxy clusters from the inner regions out to the peripheries. We present a full three-dimensional analysis combining strong and weak lensing, X-ray surface brightness and temperature, and the Sunyaev-Zel'dovich effect. The method is applied to MACS J1206.2-0847, a remarkably regular, face-on, massive, M_{200}=(1.1+-0.2)*10^{15}M_Sun/h, cluster at z=0.44. The measured concentration, c_{200}=6.3+-1.2, and the triaxial shape are common to halos formed in a LCDM scenario. The gas has settled in and follows the shape of the gravitational potential, which is evidence of pressure equilibrium via the shape theorem. There is no evidence for significant non-thermal pressure and the equilibrium is hydrostatic.

To investigate the relationship between thermal and non-thermal components in merger galaxy clusters, we present deep JVLA and Chandra observations of the HST Frontier Fields cluster MACS J0717.5+3745. The Chandra image shows a complex merger event, with at least four components belonging to different merging subclusters. NW of the cluster, $\sim 0.7$ Mpc from the center, there is a ram-pressure-stripped core that appears to have traversed the densest parts of the cluster after entering the ICM from the direction of a galaxy filament to the SE. We detect a density discontinuity NNE of this core which we speculate is associated with a cold front. Our radio images reveal new details for the complex radio relic and radio halo in this cluster. In addition, we discover several new filamentary radio sources with sizes of 100-300 kpc. A few of these seem to be connected to the main radio relic, while others are either embedded within the radio halo or projected onto it. A narrow-angled-tailed (NAT) radio galaxy, a cluster member, is located at the center of the radio relic. The steep spectrum tails of this AGN leads into the large radio relic where the radio spectrum flattens again. This morphological connection between the NAT radio galaxy and relic provides evidence for re-acceleration (revival) of fossil electrons. The presence of hot $\gtrsim 20$ keV ICM gas detected by Chandra near the relic location provides additional support for this re-acceleration scenario.

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 a joint analysis of Chandra X-ray observations, Bolocam thermal Sunyaev-Zel'dovich (SZ) effect observations, Hubble Space Telescope (HST) strong lensing data, and HST and Subaru Suprime-Cam weak lensing data. The multiwavelength dataset is used to constrain parametric models for the distribution of dark and baryonic matter in a sample of six massive galaxy clusters selected from the Cluster Lensing And Supernova survey with Hubble (CLASH). For five of the six clusters, the multiwavelength dataset is well described by a relatively simple model that assumes spherical symmetry, hydrostatic equilibrium, and entirely thermal pressure support. The joint analysis yields considerably better constraints on the total mass and concentration of the cluster compared to analysis of any one dataset individually. The subsample of five galaxy clusters is used to place an upper limit on the fraction of pressure support in the intracluster medium (ICM) due to nonthermal processes, such as turbulence and bulk flow of the gas. We constrain the nonthermal pressure fraction at r500c to be less than 0.11 at 95 percent confidence. This is in tension with state-of-the-art hydrodynamical simulations, which predict a nonthermal pressure fraction of approximately 0.25 at r500c for clusters of similar mass and redshift. This tension may be explained by the sample selection and/or our assumption of spherical symmetry.

We propose a new method to probe the Warm Hot Intergalactic Medium (WHIM) beyond the virial radius (R_200) of a cluster of galaxies, where X-ray observations are not easily achievable. In this method, we use dispersion measures (DMs) of Fast Radio Bursts (FRBs) that appear behind the cluster and the Sunyaev-Zel'dovich (SZ) effect towards the cluster. The DMs reflect the density of the intracluster medium (ICM) including the WHIM. If we observe a sufficient number of FRBs in the direction of the cluster, we can derive the density profile from the DMs. Similarly, we can derive the pressure profile from the SZ effect. By combining the density and the pressure profiles, the temperature profile can be obtained. Based on mock observations of nearby clusters, we find that the density of the WHIM can be determined even at > 2 R_200 from the cluster center when FRB observations with the Square Kilometre Array (SKA) become available. The temperature can be derived out to r~ 1.5 R_200, and the radius is limited by the current sensitivity of SZ observations.

The Yuan-Tseh Lee Array for Microwave Background Anisotropy (AMiBA) is a co-planar interferometer array operating at a wavelength of 3mm to measure the Sunyaev-Zeldovich effect (SZE) of galaxy clusters. In the first phase of operation -- with a compact 7-element array with 0.6m antennas (AMiBA-7) -- we observed six clusters at angular scales from 5\arcmin to 23\arcmin. Here, we describe the expansion of AMiBA to a 13-element array with 1.2m antennas (AMiBA-13), its subsequent commissioning, and our cluster SZE observing program. The most important changes compared to AMiBA-7 are (1) array re-configuration with baselines ranging from 1.4m to 4.8m covering angular scales from 2\arcmin to 11.5\arcmin, (2) thirteen new lightweight carbon-fiber-reinforced plastic (CFRP) 1.2m reflectors, and (3) additional correlators and six new receivers. From the AMiBA-13 SZE observing program, we present here maps of a subset of twelve clusters. In highlights, we combine AMiBA-7 and AMiBA-13 observations of Abell 1689 and perform a joint fitting assuming a generalized NFW pressure profile. Our cylindrically integrated Compton-y values for this cluster are consistent with the BIMA/OVRA, SZA, and Planck results. We report the first targeted SZE detection towards the optically selected galaxy cluster RCS J1447+0828, and we demonstrate the ability of AMiBA SZE data to serve as a proxy for the total cluster mass. Finally, we show that our AMiBA-SZE derived cluster masses are consistent with recent lensing mass measurements in the literature.

We report on high-resolution JVLA and Chandra observations of the HST Frontier Cluster MACS J0717.5+3745. MACS J0717.5+3745 offers the largest contiguous magnified area of any known cluster, making it a promising target to search for lensed radio and X-ray sources. With the high-resolution 1.0-6.5 GHz JVLA imaging in A and B configuration, we detect a total of 51 compact radio sources within the area covered by the HST imaging. Within this sample we find 7 lensed sources with amplification factors larger than $2$. None of these sources are identified as multiply-lensed. Based on the radio luminosities, the majority of these sources are likely star forming galaxies with star formation rates of 10-50 M$_\odot$ yr$^{-1}$ located at $1 \lesssim z \lesssim 2$. Two of the lensed radio sources are also detected in the Chandra image of the cluster. These two sources are likely AGN, given their $2-10$ keV X-ray luminosities of $\sim 10^{43-44}$ erg s$^{-1}$. From the derived radio luminosity function, we find evidence for an increase in the number density of radio sources at $0.6<z<2.0$, compared to a $z < 0.3$ sample. Our observations indicate that deep radio imaging of lensing clusters can be used to study star forming galaxies, with star formation rates as low as $\sim10$ M$_{\odot}$ yr$^{-1}$, at the peak of cosmic star formation history.

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

We present a comprehensive analysis of strong-lensing, weak-lensing shear and magnification data for a sample of 16 X-ray-regular and 4 high-magnification galaxy clusters selected from the CLASH survey. Our analysis combines constraints from 16-band HST observations and wide-field multi-color imaging taken primarily with Subaru/Suprime-Cam. We reconstruct surface mass density profiles of individual clusters from a joint analysis of the full lensing constraints, and determine masses and concentrations for all clusters. We find internal consistency of the ensemble mass calibration to be $\le 5\% \pm 6\%$ by comparison with the CLASH weak-lensing-only measurements of Umetsu et al. For the X-ray-selected subsample, we examine the concentration-mass relation and its intrinsic scatter using a Bayesian regression approach. Our model yields a mean concentration of $c|_{z=0.34} = 3.95 \pm 0.35$ at $M_{200c} \simeq 14\times 10^{14}M_\odot$ and an intrinsic scatter of $\sigma(\ln c_{200c}) = 0.13 \pm 0.06$, in excellent agreement with LCDM predictions when the CLASH selection function based on X-ray morphological regularity and the projection effects are taken into account. We also derive an ensemble-averaged surface mass density profile for the X-ray-selected subsample by stacking their individual profiles. The stacked mass profile is well described by a family of density profiles predicted for cuspy dark-matter-dominated halos, namely, the NFW, Einasto, and DARKexp models, whereas the single power-law, cored isothermal and Burkert density profiles are disfavored by the data. We show that cuspy halo models that include the two-halo term provide improved agreement with the data. For the NFW halo model, we measure a mean concentration of $c_{200c} = 3.79^{+0.30}_{-0.28}$ at $M_{200c} = 14.1^{+1.0}_{-1.0}\times 10^{14}M_\odot$, demonstrating consistency between complementary analysis methods.

The HST Frontier Fields cluster MACS J1149.6+2223 is one of the most complex merging clusters, believed to consist of four dark matter halos. We present results from deep (365 ks) Chandra observations of the cluster, which reveal the most distant cold front (z=0.544) discovered to date. In the cluster outskirts, we also detect hints of a surface brightness edge that could be the bow shock preceding the cold front. The substructure analysis of the cluster identified several components with large relative radial velocities, thus indicating that at least some collisions occur almost along the line of sight. The inclination of the mergers with respect to the plane of the sky poses significant observational challenges at X-ray wavelengths. MACS J1149.6+2223 possibly hosts a steep-spectrum radio halo. If the steepness of the radio halo is confirmed, then the radio spectrum, combined with the relatively regular ICM morphology, could indicate that MACS J1149.6+2223 is an old merging cluster.

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 use high-precision kinematic and lensing measurements of the total mass profile of the dynamically relaxed galaxy cluster MACS J1206.2-0847 at $z=0.44$ to estimate the value of the ratio $\eta=\Psi/\Phi$ between the two scalar potentials in the linear perturbed Friedmann-Lemaitre-Robertson-Walker metric.[...] Complementary kinematic and lensing mass profiles were derived from exhaustive analyses using the data from the Cluster Lensing And Supernova survey with Hubble (CLASH) and the spectroscopic follow-up with the Very Large Telescope (CLASH-VLT). Whereas the kinematic mass profile tracks only the time-time part of the perturbed metric (i.e. only $\Phi$), the lensing mass profile reflects the contribution of both time-time and space-space components (i.e. the sum $\Phi+\Psi$). We thus express $\eta$ as a function of the mass profiles and perform our analysis over the radial range $0.5\,Mpc\le r\le r_{200}=1.96\,Mpc$. Using a spherical Navarro-Frenk-White mass profile, which well fits the data, we obtain $\eta(r_{200})=1.01\,_{-0.28}^{+0.31}$ at the 68\% C.L. We discuss the effect of assuming different functional forms for mass profiles and of the orbit anisotropy in the kinematic reconstruction. Interpreting this result within the well-studied $f(R)$ modified gravity model, the constraint on $\eta$ translates into an upper bound to the interaction length (inverse of the scalaron mass) smaller than 2 Mpc. This tight constraint on the $f(R)$ interaction range is however substantially relaxed when systematic uncertainties in the analysis are considered. Our analysis highlights the potential of this method to detect deviations from general relativity, while calling for the need of further high-quality data on the total mass distribution of clusters and improved control on systematic effects.

We present a strong lensing system, composed of 4 multiple images of a source at z = 2.387, created by two lens galaxies, G1 and G2, belonging to the galaxy cluster MACS J1115.9+0129 at z = 0.353. We use observations taken as part of the Cluster Lensing and Supernova survey with Hubble, CLASH, and its spectroscopic follow-up programme at the Very Large Telescope, CLASH-VLT, to estimate the total mass distributions of the two galaxies and the cluster through strong gravitational lensing models. We find that the total projected mass values within the half-light radii, R_{e}, of the two lens galaxies are M_{T,G1}(< R_{e,G1}) = (3.6 +/- 0.4) x 10^{10}M_{Sun} and M_{T,G2}(< R_{e,G2}) = (4.2 +/- 1.6) x 10^{10}M_{Sun}. The effective velocity dispersion values of G1 and G2 are (122 +/- 7) km/s and (137 +/- 27) km/s, respectively. We remark that these values are relatively low when compared to those of ~200-300 km/s, typical of lens galaxies found in the field by previous surveys. By fitting the spectral energy distributions of G1 and G2, we measure projected luminous over total mass fractions within R_{e} of 0.11 +/- 0.03, for G1, and 0.73 +/- 0.32, for G2. The fact that the less massive galaxy, G1, is dark-matter dominated in its inner regions raises the question of whether the dark matter fraction in the core of early-type galaxies depends on their mass. Further investigating strong lensing systems will help us understand the influence that dark matter has on the structure and evolution of the inner regions of galaxies.

We perform a comprehensive study of the total mass distribution of the galaxy cluster RXCJ2248 ($z=0.348$) with a set of high-precision strong lensing models, which take advantage of extensive spectroscopic information on many multiply lensed systems. In the effort to understand and quantify inherent systematics in parametric strong lensing modelling, we explore a collection of 22 models where we use different samples of multiple image families, parametrizations of the mass distribution and cosmological parameters. As input information for the strong lensing models, we use the CLASH HST imaging data and spectroscopic follow-up observations, carried out with the VIMOS and MUSE spectrographs, to identify bona-fide multiple images. A total of 16 background sources, over the redshift range $1.0-6.1$, are multiply lensed into 47 images, 24 of which are spectroscopically confirmed and belong to 10 individual sources. The cluster total mass distribution and underlying cosmology in the models are optimized by matching the observed positions of the multiple images on the lens plane. We show that with a careful selection of a sample of spectroscopically confirmed multiple images, the best-fit model reproduces their observed positions with a rms of $0.3$ in a fixed flat $\Lambda$CDM cosmology, whereas the lack of spectroscopic information lead to biases in the values of the model parameters. Allowing cosmological parameters to vary together with the cluster parameters, we find (at $68\%$ confidence level) $\Omega_m=0.25^{+0.13}_{-0.16}$ and $w=-1.07^{+0.16}_{-0.42}$ for a flat $\Lambda$CDM model, and $\Omega_m=0.31^{+0.12}_{-0.13}$ and $\Omega_\Lambda=0.38^{+0.38}_{-0.27}$ for a universe with $w=-1$ and free curvature. Using toy models mimicking the overall configuration of RXCJ2248, we estimate the impact of the line of sight mass structure on the positional rms to be $0.3\pm 0.1$.(ABRIDGED)

We perform a weak-lensing study of the nearby cool-core galaxy clusters, Hydra A ($z=0.0538$) and A478 ($z=0.0881$), of which brightest cluster galaxies (BCGs) host powerful activities of active galactic nuclei (AGNs). For each cluster, the observed tangential shear profile is well described either by a single Navarro--Frenk--White model or a two-component model including the BCG as an unresolved point mass. For A478, we determine the BCG and its host-halo masses from a joint fit to weak-lensing and stellar photometry measurements. We find that the choice of initial mass functions (IMFs) can introduce a factor of two uncertainty in the BCG mass, whereas the BCG host halo mass is well constrained by data. We perform a joint analysis of weak-lensing and stellar kinematics data available for the Hydra A cluster, which allows us to constrain the central mass profile without assuming specific IMFs. We find that the central mass profile ($r<300$kpc) determined from the joint analysis is in excellent agreement with those from independent measurements, including dynamical masses estimated from the cold gas disk component, X-ray hydrostatic total mass estimates, and the central stellar mass estimated based on the Salpeter IMF. The observed dark-matter fraction around the BCG for Hydra A is found to be smaller than those predicted by adiabatic contraction models, suggesting the importance of other physical processes, such as the the AGN feedback and/or dissipationless mergers.

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 developed an algorithm to find and characterize gravitationally lensed galaxies (arcs) to perform a comparison of the observed and simulated arc abundance. Observations are from the Cluster Lensing And Supernova survey with Hubble (CLASH). Simulated CLASH images are created using the MOKA package and also clusters selected from the high resolution, hydrodynamical simulations, MUSIC, over the same mass and redshift range as the CLASH sample. The algorithm' s arc elongation accuracy, completeness and false positive rate are determined and used to compute an estimate of the true arc abundance. We derive a lensing efficiency of $4 \pm 1$ arcs (with length $\ge 6"$ and length-to-width ratio $\ge 7$) per cluster for the X-ray selected CLASH sample, $4 \pm 1$ arcs per cluster for the MOKA simulated sample and $3 \pm 1$ arcs per cluster for the MUSIC simulated sample. The observed and simulated arc statistics are in full agreement. We measure the photometric redshifts of all detected arcs and find a median redshift $z_s = 1.9$ with 33% of the detected arcs having $z_s > 3$. We find that the arc abundance does not depend strongly on the source redshift distribution but is sensitive to the mass distribution of the dark matter halos (e.g. the $c-M$ relation). Our results show that consistency between the observed and simulated distributions of lensed arc sizes and axial ratios can be achieved by using cluster-lensing simulations that are carefully matched to the selection criteria used in the observations.

We test the assumption of hydrostatic equilibrium in an X-ray luminosity selected sample of 50 galaxy clusters at $0.15<z<0.3$ from the Local Cluster Substructure Survey (LoCuSS). Our weak-lensing measurements of $M_{500}$ control systematic biases to sub-4 per cent, and our hydrostatic measurements of the same achieve excellent agreement between XMM-Newton and Chandra. The mean ratio of X-ray to lensing mass for these 50 clusters is $\beta_{\rm X}=0.95\pm0.05$, and for the 44 clusters also detected by Planck, the mean ratio of Planck mass estimate to LoCuSS lensing mass is $\beta_{\rm P}=0.95\pm0.04$. Based on a careful like-for-like analysis, we find that LoCuSS, the Canadian Cluster Comparison Project (CCCP), and Weighing the Giants (WtG) agree on $\beta_{\rm P}\simeq0.9-0.95$ at $0.15<z<0.3$. This small level of hydrostatic bias disagrees at $\sim5\sigma$ with the level required to reconcile Planck cosmology results from the cosmic microwave background and galaxy cluster counts.

We present a weak-lensing analysis for the merging galaxy cluster, CIZA J2242.8+5301, hosting double radio relics, using three-band Subaru/Suprime-Cam imaging ($Br'z'$). Since the lifetime of dark matter halos colliding into clusters is longer than that of X-ray emitting gas halos, weak-lensing analysis is a powerful method to constrain a merger dynamics. Two-dimensional shear fitting using a clean background catalog suggests that the cluster undergoes a merger with a mass ratio of about 2:1. The main halo is located around the gas core in the southern region, while no concentrated gas core is associated with the northern sub halo. We find that the projected cluster mass distribution resulting from an unequal-mass merger is in excellent agreement with the curved shapes of the two radio relics and the overall X-ray morphology except for the lack of the northern gas core. The lack of a prominent radio halo enables us to constrain an upper limit of the fractional energy of magneto-hydrodynamics turbulence of $(\delta B/B)^2<\mathcal{O}(10^{-6})$ at a resonant wavenumber, by a balance between the acceleration time and the time after the core passage or the cooling time, with an assumption of resonant acceleration by second-order Fermi process.

Merging galaxy clusters leave long-lasting signatures on the baryonic and non-baryonic cluster constituents, including shock fronts, cold fronts, X-ray substructure, radio halos, and offsets between the dark matter and the gas components. Using observations from Chandra, the Jansky Very Large Array, the Giant Metrewave Radio Telescope, and the Hubble Space Telescope, we present a multiwavelength analysis of the merging Frontier Fields cluster MACS J0416.1-2403 (z=0.396), which consists of a NE and a SW subclusters whose cores are separated on the sky by ~250 kpc. We find that the NE subcluster has a compact core and hosts an X-ray cavity, yet it is not a cool core. Approximately 450 kpc south-south west of the SW subcluster, we detect a density discontinuity that corresponds to a compression factor of ~1.5. The discontinuity was most likely caused by the interaction of the SW subcluster with a less massive structure detected in the lensing maps SW of the subcluster's center. For both the NE and the SW subclusters, the dark matter and the gas components are well-aligned, suggesting that MACS J0416.1-2403 is a pre-merging system. The cluster also hosts a radio halo, which is unusual for a pre-merging system. The halo has a 1.4 GHz power of (1.06 +/- 0.09) x 10^{24} W Hz^{-1}, which is somewhat lower than expected based on the X-ray luminosity of the cluster. We suggest that we are either witnessing the birth of a radio halo, or have discovered a rare ultra-steep spectrum halo.

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 scaling relations between the integrated Sunyaev-Zel'dovich Effect (SZE) signal, $Y_{\rm SZ}$, its X-ray analogue, $Y_{\rm X}\equiv M_{\rm gas}T_{\rm X}$, and total mass, $M_{\rm tot}$, for the 45 galaxy clusters in the Bolocam X-ray-SZ (BOXSZ) sample. All parameters are integrated within $r_{2500}$. $Y_{2500}$ values are measured using SZE data collected with Bolocam, operating at 140 GHz at the Caltech Submillimeter Observatory (CSO). The temperature, $T_{\rm X}$, and mass, $M_{\rm gas,2500}$, of the intracluster medium are determined using X-ray data collected with Chandra, and $M_{\rm tot}$ is derived from $M_{\rm gas}$ assuming a constant gas mass fraction. Our analysis accounts for several potential sources of bias, including: selection effects, contamination from radio point sources, and the loss of SZE signal due to noise filtering and beam-smoothing effects. We measure the $Y_{2500}$--$Y_{\rm X}$ scaling to have a power-law index of $0.84\pm0.07$, and a fractional intrinsic scatter in $Y_{2500}$ of $(21\pm7)\%$ at fixed $Y_{\rm X}$, both of which are consistent with previous analyses. We also measure the scaling between $Y_{2500}$ and $M_{2500}$, finding a power-law index of $1.06\pm0.12$ and a fractional intrinsic scatter in $Y_{2500}$ at fixed mass of $(25\pm9)\%$. While recent SZE scaling relations using X-ray mass proxies have found power-law indices consistent with the self-similar prediction of 5/3, our measurement stands apart by differing from the self-similar prediction by approximately 5$\sigma$. Given the good agreement between the measured $Y_{2500}$--$Y_{\rm X}$ scalings, much of this discrepancy appears to be caused by differences in the calibration of the X-ray mass proxies adopted for each particular analysis.

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 present results from a comprehensive lensing analysis in HST data, of the complete CLASH cluster sample. We identify new multiple-images previously undiscovered allowing improved or first constraints on the cluster inner mass distributions and profiles. We combine these strong-lensing constraints with weak-lensing shape measurements within the HST FOV to jointly constrain the mass distributions. The analysis is performed in two different common parameterizations (one adopts light-traces-mass for both galaxies and dark matter while the other adopts an analytical, elliptical NFW form for the dark matter), to provide a better assessment of the underlying systematics - which is most important for deep, cluster-lensing surveys, especially when studying magnified high-redshift objects. We find that the typical (median), relative systematic differences throughout the central FOV are $\sim40\%$ in the (dimensionless) mass density, $\kappa$, and $\sim20\%$ in the magnification, $\mu$. We show maps of these differences for each cluster, as well as the mass distributions, critical curves, and 2D integrated mass profiles. For the Einstein radii ($z_{s}=2$) we find that all typically agree within $10\%$ between the two models, and Einstein masses agree, typically, within $\sim15\%$. At larger radii, the total projected, 2D integrated mass profiles of the two models, within $r\sim2\arcmin$, differ by $\sim30\%$. Stacking the surface-density profiles of the sample from the two methods together, we obtain an average slope of $d\log (\Sigma)/d\log(r)\sim-0.64\pm0.1$, in the radial range [5,350] kpc. Lastly, we also characterize the behavior of the average magnification, surface density, and shear differences between the two models, as a function of both the radius from the center, and the best-fit values of these quantities.

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 detailed mass reconstruction and a novel study on the substructure properties in the core of the CLASH and Frontier Fields galaxy cluster MACS J0416.1-2403. We show and employ our extensive spectroscopic data set taken with the VIMOS instrument as part of our CLASH-VLT program, to confirm spectroscopically 10 strong lensing systems and to select a sample of 175 plausible cluster members to a limiting stellar mass of log(M_*/M_Sun) ~ 8.6. We reproduce the measured positions of 30 multiple images with a remarkable median offset of only 0.3" by means of a comprehensive strong lensing model comprised of 2 cluster dark-matter halos, represented by cored elliptical pseudo-isothermal mass distributions, and the cluster member components. The latter have total mass-to-light ratios increasing with the galaxy HST/WFC3 near-IR (F160W) luminosities. The measurement of the total enclosed mass within the Einstein radius is accurate to ~5%, including systematic uncertainties. We emphasize that the use of multiple-image systems with spectroscopic redshifts and knowledge of cluster membership based on extensive spectroscopic information is key to constructing robust high-resolution mass maps. We also produce magnification maps over the central area that is covered with HST observations. We investigate the galaxy contribution, both in terms of total and stellar mass, to the total mass budget of the cluster. When compared with the outcomes of cosmological $N$-body simulations, our results point to a lack of massive subhalos in the inner regions of simulated clusters with total masses similar to that of MACS J0416.1-2403. Our findings of the location and shape of the cluster dark-matter halo density profiles and on the cluster substructures provide intriguing tests of the assumed collisionless, cold nature of dark matter and of the role played by baryons in the process of structure formation.

We present high resolution (9$^{\prime \prime}$) imaging of the Sunyaev-Zel'dovich Effect (SZE) toward two massive galaxy clusters, MACS J0647.7+7015 ($z=0.591$) and MACS J1206.2-0847 ($z=0.439$). We compare these 90 GHz measurements, taken with the MUSTANG receiver on the Green Bank Telescope, with generalized Navarro-Frenk-White (gNFW) models derived from Bolocam 140 GHz SZE data as well as maps of the thermal gas derived from {\it Chandra} X-ray observations. For MACS J0647.7+7015, we find a gNFW profile with core slope parameter $\gamma= 0.9$ fits the MUSTANG image with $\chi^{2}_{red}=1.005$ and probability to exceed (PTE) = 0.34. For MACS J1206.2-0847, we find $\gamma=0.7$, $\chi^{2}_{red}=0.993$, and PTE = 0.70. In addition, we find a significant ($>$3-$\sigma$) residual SZE feature in MACS J1206.2-0847 coincident with a group of galaxies identified in VLT data and filamentary structure found in a weak-lensing mass reconstruction. We suggest the detected sub-structure may be the SZE decrement from a low mass foreground group or an infalling group. GMRT measurements at 610 MHz reveal diffuse extended radio emission to the west, which we posit is either an AGN-driven radio lobe, a bubble expanding away from disturbed gas associated with the SZE signal, or a bubble detached and perhaps re-accelerated by sloshing within the cluster. Using the spectroscopic redshifts available, we find evidence for a foreground ($z=0.423$) or infalling group, coincident with the residual SZE feature.

A clean measurement of the evolution of the galaxy cluster mass function can significantly improve our understanding of cosmology from the rapid growth of cluster masses below z < 0.5. Here we examine the consistency of cluster catalogues selected from the SDSS by applying two independent gravity-based methods using all available spectroscopic redshifts from the DR10 release. First, we detect a gravitational redshift related signal for 20,119 and 13,128 clusters with spectroscopic redshifts contained in the GMBCG and redMaPPer catalogues, respectively, at a level of $\sim - 10$ km s$^{-1}$. This we show is consistent with the magnitude expected using the richness-mass relations provided by the literature and after applying recently clarified relativistic and flux bias corrections. This signal is also consistent with the richest clusters in the larger catalogue of Wen et al. (2012), corresponding to $M_{200m} \gtrsim 2 \times 10^{14}\,\mathrm{M}_\odot\,h^{-1}$, however we find no significant detection of gravitational redshift signal for less riched clusters, which may be related to bulk motions from substructure and spurious cluster detections. Second, we find all three catalogues generate mass-dependent levels of lensing magnification bias, which enhances the mean redshift of flux-selected background galaxies from the BOSS survey. The magnitude of this lensing effect is generally consistent with the corresponding richness-mass relations advocated for the surveys. We conclude that all catalogues comprise a high proportion of reliable clusters, and that the GMBCG and redMaPPer cluster finder algorithms favor more relaxed clusters with a meaningful gravitational redshift signal, as anticipated by the red-sequence colour selection of the GMBCG and redMaPPer samples.

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 conduct a joint X-ray and weak-lensing study of four relaxed galaxy clusters (Hydra A, A478, A1689 and A1835) observed by both Suzaku and Subaru out to virial radii, with an aim to understand recently-discovered unexpected feature of the ICM in cluster outskirts. We show that the average hydrostatic-to-lensing total mass ratio for the four clusters decreases from \sim 70% to \sim 40% as the overdensity contrast decreases from 500 to the virial value.The average gas mass fraction from lensing total mass estimates increases with cluster radius and agrees with the cosmic mean baryon fraction within the virial radius, whereas the X-ray-based gas fraction considerably exceeds the cosmic values due to underestimation of the hydrostatic mass. We also develop a new advanced method for determining normalized cluster radial profiles for multiple X-ray observables by simultaneously taking into account both their radial dependence and multivariate scaling relations with weak-lensing masses. Although the four clusters span a range of halo mass, concentration, X-ray luminosity and redshift, we find that the gas entropy, pressure, temperature and density profiles are all remarkably self-similar when scaled with the lensing M_200 mass and r_200 radius.The entropy monotonically increases out to \sim 0.5r_200 following the accretion shock heating model K(r)\propto r^1.1, and flattens at \simgt 0.5r_200.The universality of the scaled entropy profiles indicates that the thermalization mechanism over the entire cluster region (>0.1r_200) is controlled by gravitation in a common to all clusters, although the heating efficiency in the outskirts needs to be modified from the standard law.The bivariate scaling functions of the gas density and temperature reveal that the flattening of the outskirts entropy profile is caused by the steepening of the temperature, rather than the flattening of the gas density.

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.

The Javalambre-Physics of the Accelerated Universe Astrophysical Survey (J-PAS) is a narrow band, very wide field Cosmological Survey to be carried out from the Javalambre Observatory in Spain with a purpose-built, dedicated 2.5m telescope and a 4.7 sq.deg. camera with 1.2Gpix. Starting in late 2015, J-PAS will observe 8500sq.deg. of Northern Sky and measure $0.003(1+z)$ photo-z for $9\times10^7$ LRG and ELG galaxies plus several million QSOs, sampling an effective volume of $\sim 14$ Gpc$^3$ up to $z=1.3$ and becoming the first radial BAO experiment to reach Stage IV. J-PAS will detect $7\times 10^5$ galaxy clusters and groups, setting constrains on Dark Energy which rival those obtained from its BAO measurements. Thanks to the superb characteristics of the site (seeing ~0.7 arcsec), J-PAS is expected to obtain a deep, sub-arcsec image of the Northern sky, which combined with its unique photo-z precision will produce one of the most powerful cosmological lensing surveys before the arrival of Euclid. J-PAS unprecedented spectral time domain information will enable a self-contained SN survey that, without the need for external spectroscopic follow-up, will detect, classify and measure $\sigma_z\sim 0.5\%$ redshifts for $\sim 4000$ SNeIa and $\sim 900$ core-collapse SNe. The key to the J-PAS potential is its innovative approach: a contiguous system of 54 filters with $145\AA$ width, placed $100\AA$ apart over a multi-degree FoV is a powerful "redshift machine", with the survey speed of a 4000 multiplexing low resolution spectrograph, but many times cheaper and much faster to build. The J-PAS camera is equivalent to a 4.7 sq.deg. "IFU" and it will produce a time-resolved, 3D image of the Northern Sky with a very wide range of Astrophysical applications in Galaxy Evolution, the nearby Universe and the study of resolved stellar populations.

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.

Hubble Space Telescope imaging of the galaxy cluster Abell 1689 has revealed an exceptional number of strongly lensed multiply-imaged galaxies, including high-redshift candidates. Previous studies have used this data to obtain the most detailed dark matter reconstructions of any galaxy cluster to date, resolving substructures ~25 kpc across. We examine Abell 1689 (hereafter, A1689) non-parametrically, combining strongly lensed images and weak distortions from wider field Subaru imaging, and we incorporate member galaxies to improve the lens solution. Strongly lensed galaxies are often locally affected by member galaxies, however, these perturbations cannot be recovered in grid based reconstructions because the lensing information is too sparse to resolve member galaxies. By adding luminosity-scaled member galaxy deflections to our smooth grid we can derive meaningful solutions with sufficient accuracy to permit the identification of our own strongly lensed images, so our model becomes self consistent. We identify 11 new multiply lensed system candidates and clarify previously ambiguous cases, in the deepest optical and NIR data to date from Hubble and Subaru. Our improved spatial resolution brings up new features not seen when the weak and strong lensing effects are used separately, including clumps and filamentary dark matter around the main halo. Our treatment means we can obtain an objective mass ratio between the cluster and galaxy components, for examining the extent of tidal stripping of the luminous member galaxies. We find a typical mass-to-light ratios of M/L_B = 21 inside the r<1 arcminute region that drops to M/L_B = 17 inside the r<40 arcsecond region. Our model independence means we can objectively evaluate the competitiveness of stacking cluster lenses for defining the geometric lensing-distance-redshift relation in a model independent way.

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.

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 report our analysis of MACS J0717.5+3745 using 140 and 268 GHz Bolocam data collected at the Caltech Submillimeter Observatory. We detect extended Sunyaev-Zel'dovich (SZ) effect signal at high significance in both Bolocam bands, and we employ Herschel-SPIRE observations to subtract the signal from dusty background galaxies in the 268 GHz data. We constrain the two-band SZ surface brightness toward two of the sub-clusters of MACS J0717.5+3745: the main sub-cluster (named C), and a sub-cluster identified in spectroscopic optical data to have a line-of-sight velocity of +3200 km/s (named B). We determine the surface brightness in two separate ways: via fits of parametric models and via direct integration of the images. For both sub-clusters, we find consistent surface brightnesses from both analysis methods. We constrain spectral templates consisting of relativistically corrected thermal and kinetic SZ signals, using a jointly-derived electron temperature from Chandra and XMM-Newton under the assumption that each sub-cluster is isothermal. The data show no evidence for a kinetic SZ signal toward sub-cluster C, but they do indicate a significant kinetic SZ signal toward sub-cluster B. The model-derived surface brightnesses for sub-cluster B yield a best-fit line-of-sight velocity of v_z = +3450 +- 900 km/s, with (1 - Prob[v_z > 0]) = 1.3 x 10^-5 (4.2 sigma away from 0 for a Gaussian distribution). The directly integrated sub-cluster B SZ surface brightnesses provide a best-fit v_z = +2550 +- 1050 km/s, with (1 - Prob[v_z > 0]) = 2.2 x 10^-3 (2.9 sigma).

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)

We show that the fast moving component of the "bullet cluster" (1E0657-56) can induce potentially resolvable redshift differences between multiply-lensed images of background galaxies. The moving cluster effect can be expressed as the scalar product of the lensing deflection angle with the tangential velocity of the mass components, and it is maximal for clusters colliding in the plane of the sky with velocities boosted by their mutual gravity. The bullet cluster is likely to be the best candidate for the first measurement of this effect due to the large collision velocity and because the lensing deflection and the cluster fields can be calculated in advance. We derive the deflection field using multiply-lensed background galaxies detected with the Hubble Space Telescope. The velocity field is modeled using self-consistent N-body/hydrodynamical simulations constrained by the observed X-ray and gravitational lensing features of this system. We predict that the triply-lensed images of systems "G" and "H" straddling the critical curve of the bullet component will show the largest frequency shifts up to ~0.5 km/sec. This is within the range of the Atacama Large Millimeter/sub-millimeter Array (ALMA) for molecular emission, and is near the resolution limit of the new generation high-throughput optical-IR spectrographs. A detection of this effect measures the tangential motion of the subclusters directly, thereby clarifying the tension with LCDM, which is inferred from gas motion less directly. This method may be extended to smaller redshift differences using the Ly-alpha forest towards QSOs lensed by more typical clusters of galaxies. More generally, the tangential component of the peculiar velocities of clusters derived by our method complements the radial component determined by the kinematic SZ effect, providing a full 3-dimensional description of velocities.

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 describe Sunyaev-Zel'dovich (SZ) effect measurements and analysis of the intracluster medium (ICM) pressure profiles of a set of 45 massive galaxy clusters imaged using Bolocam at the Caltech Submillimeter Observatory. We have used masses determined from Chandra X-ray observations to scale each cluster's profile by the overdensity radius R500 and the mass-and-redshift-dependent normalization factor P500. We deproject the average pressure profile of our sample into 13 logarithmically spaced radial bins between 0.07R500 and 3.5R500. We find that a generalized Navarro, Frenk, and White (gNFW) profile describes our data with sufficient goodness-of-fit and best-fit parameters (C500, alpha, beta, gamma, P0 = 1.18, 0.86, 3.67, 0.67, 4.29). We also use the X-ray data to define cool-core and disturbed subsamples of clusters, and we constrain the average pressure profiles of each of these subsamples. We find that given the precision of our data the average pressure profiles of disturbed and cool-core clusters are consistent with one another at R>~0.15R500, with cool-core systems showing indications of higher pressure at R<~0.15R500. In addition, for the first time, we place simultaneous constraints on the mass scaling of cluster pressure profiles, their ensemble mean profile, and their radius-dependent intrinsic scatter between 0.1R500 and 2.0R500. The scatter among profiles is minimized at radii between ~0.2R500 and ~0.5R500, with a value of ~20%. The best-fit mass scaling has a power-law slope of 0.49, which is shallower than the nominal prediction of 2/3 from self-similar hydrostatic equilibrium models. These results for the intrinsic scatter and mass scaling are largely consistent with previous analyses, most of which have relied heavily on X-ray derived pressures of clusters at significantly lower masses and redshifts compared to our sample.

We present a new way to solve the platform deformation problem of co-planar interferometers. The platform of a co-planar interferometer can be deformed due to driving forces and gravity. A deformed platform will induce extra components into the geometric delay of each baseline, and change the phases of observed visibilities. The reconstructed images will also be diluted due to the errors of the phases. The platform deformations of The Yuan-Tseh Lee Array for Microwave Background Anisotropy (AMiBA) were modelled based on photogrammetry data with about 20 mount pointing positions. We then used the differential optical pointing error between two optical telescopes to fit the model parameters in the entire horizontal coordinate space. With the platform deformation model, we can predict the errors of the geometric phase delays due to platform deformation with given azimuth and elevation of the targets and calibrators. After correcting the phases of the radio point sources in the AMiBA interferometric data, we recover 50% - 70% flux loss due to phase errors. This allows us to restore more than 90% of a source flux. The method outlined in this work is not only applicable to the correction of deformation for other co-planar telescopes but also to single dish telescopes with deformation problems. This work also forms the basis of the upcoming science results of AMiBA-13.

Lens magnification by galaxy clusters induces characteristic spatial variations in the number counts of background sources, amplifying their observed fluxes and expanding the area of sky, the net effect of which, known as magnification bias, depends on the intrinsic faint-end slope of the source luminosity function. The bias is strongly negative for red galaxies, dominated by the geometric area distortion, whereas it is mildly positive for blue galaxies, enhancing the blue counts toward the cluster center. We generalize the Bayesian approach of Umetsu et al. for reconstructing projected cluster mass profiles, by incorporating multiple populations of background sources for magnification bias measurements and combining them with complementary lens distortion measurements, effectively breaking the mass-sheet degeneracy and improving the statistical precision of cluster mass measurements. The approach can be further extended to include strong-lensing projected mass estimates, thus allowing for non-parametric absolute mass determinations in both the weak and strong regimes. We apply this method to our recent CLASH lensing measurements of MACS J1206.2-0847, and demonstrate how combining multi-probe lensing constraints can improve the reconstruction of cluster mass profiles. This method will also be useful for a stacked lensing analysis, combining all lensing-related effects in the cluster regime, for a definitive determination of the averaged mass profile.

We present a stacked weak-lensing analysis of an approximately mass-selected sample of 50 galaxy clusters at 0.15<z<0.3, based on observations with Suprime-Cam on the Subaru Telescope. We develop a new method for selecting lensed background galaxies from which we estimate that our sample of red background galaxies suffers just 1% contamination. We detect the stacked tangential shear signal from the full sample of 50 clusters, based on this red sample of background galaxies, at a total signal-to-noise ratio of S/N=32.7. The Navarro-Frenk-White model is an excellent fit to the data, yielding sub-10% statistical precision on mass and concentration: Mvir=7.19^{+0.53}_{-0.50}\times10^{14}h^{-1}Msol, cvir=5.41^{+0.49}_{-0.45} (c_{200}=4.22^{+0.40}_{-0.36}). Tests of a range of possible systematic errors, including shear calibration and stacking-related issues, indicate that they are sub-dominant to the statistical errors. The concentration parameter obtained from stacking our approximately mass-selected cluster sample is broadly in line with theoretical predictions. Moreover, the uncertainty on our measurement is comparable with the differences between the different predictions in the literature. Overall our results highlight the potential for stacked weak-lensing methods to probe the mean mass density profile of cluster-scale dark matter halos with upcoming surveys, including Hyper-Suprime-Cam, Dark Energy Survey, and KIDS.

We report the first detection of a redshift-depth enhancement of background galaxies magnified by foreground clusters. Using 300,000 BOSS-Survey galaxies with accurate spectroscopic redshifts, we measure their mean redshift depth behind four large samples of optically selected clusters from the SDSS surveys, totalling 5,000-15,000 clusters. A clear trend of increasing mean redshift towards the cluster centers is found, averaged over each of the four cluster samples. In addition we find similar but noisier behaviour for an independent X-ray sample of 158 clusters lying in the foreground of the current BOSS sky area. By adopting the mass-richness relationships appropriate for each survey we compare our results with theoretical predictions for each of the four SDSS cluster catalogs. The radial form of this redshift enhancement is well fitted by a richness-to-mass weighted composite Navarro-Frenk-White profile with an effective mass ranging between M_200 ~ 1.4-1.8 10^14 M_sun for the optically detected cluster samples, and M_200 ~ 5.0 10^14 M_sun for the X-ray sample. This lensing detection helps to establish the credibility of these SDSS cluster surveys, and provides a normalization for their respective mass-richness relations. In the context of the upcoming bigBOSS, Subaru-PFS, and EUCLID-NISP spectroscopic surveys, this method represents an independent means of deriving the masses of cluster samples for examining the cosmological evolution, and provides a relatively clean consistency check of weak-lensing measurements, free from the systematic limitations of shear calibration.

We describe in detail our characterization of the compact radio source population in 140 GHz Bolocam observations of a set of 45 massive galaxy clusters. We use a combination of 1.4 and 30 GHz data to select a total of 28 probable cluster-member radio galaxies and also to predict their 140 GHz flux densities. All of these galaxies are steep-spectrum radio sources and they are found preferentially in the cool-core clusters within our sample. In particular, 11 of the 12 brightest cluster member radio sources are associated with cool-core systems. Although none of the individual galaxies are robustly detected in the Bolocam data, the ensemble-average flux density at 140 GHz is consistent with, but slightly lower than, the extrapolation from lower frequencies assuming a constant spectral index. In addition, our data indicate an intrinsic scatter of 30 percent around the power-law extrapolated flux densities at 140 GHz, although our data do not tightly constrain this scatter. For our cluster sample, which is composed of high-mass and moderate-redshift systems, we find that the maximum fractional change in the Sunyaev-Zel'dovich signal integrated over any single cluster due to the presence of these radio sources is 20 percent, and only 1/4 of the clusters show a fractional change of more than 1 percent. The amount of contamination is strongly dependent on cluster morphology, and nearly all of the clusters with more than 1 percent contamination are cool-core systems. This result indicates that radio contamination is not significant compared to current noise levels in 140 GHz images of massive clusters and is in good agreement with the level of radio contamination found in previous results based on lower frequency data or simulations.

We present results of four-pointing Suzaku X-ray observations (total ~200 ks) of the intracluster medium (ICM) in the Abell 1835 galaxy cluster (kT ~ 8 keV, z = 0.253) out to the virial radius (r_vir ~ 2.9 Mpc) and beyond. Faint X-ray emission from the ICM out to r_vir is detected. The temperature gradually decreases with radius from ~8 keV in the inner region to ~2 keV at r_vir. The entropy profile is shown to flatten beyond r_500, in disagreement with the r_1.1 dependence predicted from the accretion shock heating model. The thermal pressure profile in the range 0.3r_500 < r < r_vir agrees well with that obtained from the stacked Sunyaev-Zel'dovich effect observations with the Planck satellite. The hydrostatic mass profile in the cluster outskirts (r_500 < r < r_vir) falls well short of the weak lensing one derived from Subaru/Suprime-Cam observations, showing an unphysical decrease with radius. The gas mass fraction at r_vir defined with the lensing total mass agrees with the cosmic baryon fraction from the WMAP 7-year data. All these results indicate, rather than the gas-clumping effect, that the bulk of the ICM in the cluster outskirts is far from hydrostatic equilibrium and infalling matter retained some of its kinetic energy. Finally, combining with our recent Suzaku and lensing analysis of Abell 1689, a cluster of similar mass, temperature, and redshift, we show that the cluster temperature distribution in the outskirts is significantly correlated with the galaxy density field in the surrounding large-scale environment at (1-2)r_vir.

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.

The upcoming new generation of spectroscopic galaxy redshift surveys will provide large samples of cosmic voids, the distinct, large underdense structures in the universe. Combining these with future galaxy imaging surveys, we study the prospects of probing the underlying matter distribution in and around cosmic voids via the weak gravitational lensing effects of stacked voids, utilizing both shear and magnification information. The statistical precision is greatly improved by stacking together a large number of voids along different lines of sight, even when taking into account the impact of inherent miscentering and projection effects. We show that Dark Energy Task Force Stage IV surveys, such as the Euclid satellite and the Large Synoptic Survey Telescope, should be able to detect the void lensing signal with sufficient precision from stacking abundant medium-sized voids, thus providing direct constraints on the matter density profile of voids independent of assumptions on galaxy bias.

Knowledge of mass and concentration of galaxy clusters is crucial to understand their formation and evolution. Unbiased estimates require the understanding of the shape and orientation of the halo as well as its equilibrium status. We propose a novel method to determine the intrinsic properties of galaxy clusters from a multi-wavelength data set spanning from X-ray spectroscopic and photometric data to gravitational lensing to the Sunyaev-Zel'dovich effect (SZe). The method relies on two quite non informative geometrical assumptions: the distributions of total matter or gas are approximately ellipsoidal and co-aligned; they have different, constant axial ratios but share the same degree of triaxiality. Weak and strong lensing probe the features of the total mass distribution in the plane of the sky. X-ray data measure size and orientation of the gas in the plane of the sky. Comparison with the SZ amplitude fixes the elongation of the gas along the line of sight. These constraints are deprojected thanks to Bayesian inference. The mass distribution is described as a Navarro-Frenk-White halo with arbitrary orientation, gas density and temperature are modelled with parametric profiles. We applied the method to Abell 1689. Independently of the priors, the cluster is massive, M_{200}=(1.3+-0.2)*10^{15}M_sun, and over-concentrated, c_{200}=8+-1, but still consistent with theoretical predictions. The total matter is triaxial (minor to major axis ratio ~0.5+-0.1 exploiting priors from N-body simulations) with the major axis nearly orientated along the line of sight. The gas is rounder (minor to major axis ratio ~0.6+-0.1) and deviates from hydrostatic equilibrium. The contribution of non-thermal pressure is ~20-50 per cent in inner regions, <~ 300 kpc, and ~25+-5 per cent at ~1.5 Mpc.

We characterise the typical offset between the Dark Matter (DM) projected centre and the Brightest Cluster Galaxy (BCG) in 10,000 SDSS clusters. To place constraints on the centre of DM, we use an automated strong-lensing analysis, mass-modelling technique which is based on the well-tested assumption that light traces mass. The cluster galaxies are modelled with a steep power-law, and the DM component is obtained by smoothing the galaxy distribution fitting a low-order 2D polynomial (via spline interpolation), while probing a whole range of polynomial degrees and galaxy power laws. We find that the offsets between the BCG and the peak of the smoothed light map representing the DM, \Delta, are distributed equally around zero with no preferred direction, and are well described by a log-normal distribution with <log_{10}(\Delta [h^{-1} Mpc])>=-1.895^{+0.003}_{-0.004}, and \sigma=0.501\pm0.004 (95% confidence levels), or <log_{10}(\Delta [\arcsec])>=0.564\pm0.005, and \sigma=0.475\pm0.007. Some of the offsets originate in prior misidentifications of the BCG or other bright cluster members by the cluster finding algorithm, whose level we make an additional effort to assess, finding that ~10% of the clusters in the probed catalogue are likely to be misidentified, contributing to higher-end offsets in general agreement with previous studies. Our results constitute the first statistically-significant high-resolution distributions of DM-to-BCG offsets obtained in an observational analysis, and importantly show that there exists such a typical non-zero offset in the probed catalogue. The offsets show a weak positive correlation with redshift, so that higher separations are generally found for higher-z clusters in agreement with the hierarchical growth of structure, which in turn could help characterise the merger, relaxation and evolution history of clusters, in future studies. [ABRIDGED]

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 present results from strong-lens modelling of 10,000 SDSS clusters, to establish the universal distribution of Einstein radii. Detailed lensing analyses have shown that the inner mass distribution of clusters can be accurately modelled by assuming light traces mass, successfully uncovering large numbers of multiple-images. Approximate critical curves and the effective Einstein radius of each cluster can therefore be readily calculated, from the distribution of member galaxies and scaled by their luminosities. We use a subsample of 10 well-studied clusters covered by both SDSS and HST to calibrate and test this method, and show that an accurate determination of the Einstein radius and mass can be achieved by this approach "blindly", in an automated way, and without requiring multiple images as input. We present the results of the first 10,000 clusters analysed in the range $0.1<z<0.55$, and compare them to theoretical expectations. We find that for this all-sky representative sample the Einstein radius distribution is log-normal in shape, with $< Log(\theta_{e}\arcsec)>=0.73^{+0.02}_{-0.03}$, $\sigma=0.316^{+0.004}_{-0.002}$, and with higher abundance of large $\theta_{e}$ clusters than predicted by $\Lambda$CDM. We visually inspect each of the clusters with $\theta_{e}>40 \arcsec$ ($z_{s}=2$) and find that $\sim20%$ are boosted by various projection effects detailed here, remaining with $\sim40$ real giant-lens candidates, with a maximum of $\theta_{e}=69\pm12 \arcsec$ ($z_{s}=2$) for the most massive candidate, in agreement with semi-analytic calculations. The results of this work should be verified further when an extended calibration sample is available.

We present Bolocam observations of two galaxy cluster candidates reported as unconfirmed in the Planck early Sunyaev-Zel'dovich (eSZ) sample, PLCKESZ G115.71+17.52 and PLCKESZ G189.84-37.24. We observed each of these candidates with Bolocam at 140 GHz from the Caltech Submm Observatory in October 2011. The resulting images have white noise levels of ~30 {\mu}KCMB-arcmin in their central regions. We find a significant SZ decrement towards PLCKESZ G115.71. This decrement has a false detection probability of 5.3\times10-5, and we therefore confirm PLCKESZ G115.71 as a cluster. The maximum SZ decrement towards PLCKESZ G189.84 corresponds to a false detection probability of 0.027, and it therefore remains as an unconfirmed cluster candidate. In order to make our SZ-derived results more robust, we have also analyzed data from the Wide-field Infrared Survey Explorer (WISE) at the location of each cluster candidate. We find an overdensity of WISE sources consistent with other clusters in the eSZ at the location of PLCKESZ G115.71, providing further evidence that it is a cluster. We do not find a significant overdensity of WISE sources at the location of PLCKESZ G189.84.

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 have studied the star-formation and AGN activity of massive galaxies in the redshift range $z=0.4-2$, which are detected in a deep survey field using the AKARI InfraRed (IR) astronomical satellite and {\em Subaru} telescope toward the North Ecliptic Pole (NEP). The AKARI/IRC Mid-InfraRed (MIR) multiband photometry is used to trace their star-forming activities with the Polycyclic-Aromatic Hydrocarbon (PAH) emissions, which is also used to distinguish star-forming populations from AGN dominated ones and to estimate the Star Formation Rate (SFR) derived from their total emitting IR (TIR) luminosities. In combination with analyses of their stellar components, we have studied the MIR SED features of star-forming and AGN-harboring galaxies.

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.

Halos formed in the standard Lambda cold dark matter framework should follow an universal mass density profile and fit a well defined mass-concentration relation. Lensing analyses of clusters with a large Einstein radius seem to contradict this scenario, with the massive cluster Abell 1689 being often claimed as a notable example of a highly over-concentrated halo. Shape and orientation biases in lensing studies might be at the basis of this disagreement between theory and observations. We developed a method for a full three-dimensional analysis of strong and weak lensing data. Surface density maps estimated from lensing are de-projected to infer the actual triaxial structure of the cluster, whose mass distribution is approximated as an ellipsoidal Navarro-Frenk-White halo with arbitrary orientation. Inversion is performed under competing a priori assumptions, integrated in the method thanks to Bayesian statistics. We applied the method to A1689. Whatever the considered priors on shape and orientation, both weak and strong lensing analyses found the halo to be slightly over-concentrated but still consistent with theoretical predictions. We found some evidence for a mildly triaxial lens (minor to major axis ratio ~ 0.5 +- 0.2) with the major axis orientated along the line of sight. Exploiting priors from N-body simulations, we found mass M_{200} = (1.3 +- 0.4)x 10^{15}M_Sun and concentration c_{200} =10+-3 for the weak lensing analysis of Subaru data, M_{200} = (1.7 +- 0.3)x 10^{15}M_Sun and c_{200} =6.1+-0.9 for the strong lensing analysis of multiple image systems, and M_{200} = (1.3 +- 0.2)\times 10^{15}M_Sun and c_{200} =7.3+-0.8 for the combined weak plus strong analysis.

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.

Clusters of galaxies are uniquely important cosmological probes of the evolution of the large scale structure, whose diagnostic power depends quite significantly on the ability to reliably determine their masses. Clusters are typically modeled as spherical systems whose intracluster gas is in strict hydrostatic equilibrium (i.e., the equilibrium gas pressure is provided entirely by thermal pressure), with the gravitational field dominated by dark matter, assumptions that are only rough approximations. In fact, numerical simulations indicate that galaxy clusters are typically triaxial, rather than spherical, and that turbulent gas motions (induced during hierarchical merger events) provide an appreciable pressure component. Extending our previous work, we present results of a joint analysis of X-ray, weak and strong lensing measurements of Abell 1689. The quality of the data allows us to determine both the triaxial shape of the cluster and the level of non-thermal pressure that is required if the intracluster gas is in hydrostatic equilibrium. We find that the dark matter axis ratios are 1.24 +/- 0.13 and 2.02 +/- 0.01 on the plane of the sky and along the line of sight, respectively, and that about 20% of the pressure is non-thermal. Our treatment demonstrates that the dynamical properties of clusters can be determined in a (mostly) bias-free way, enhancing the use of clusters as more precise cosmological probes.

(Abridged) We use Subaru data to conduct a detailed weak-lensing study of the dark matter distribution in a sample of 30 X-ray luminous galaxy clusters at 0.15<z<0.3. A weak-lensing signal is detected at high statistical significance in each cluster, the total detection S/N ranging from 5 to 13. In this paper we concentrate on fitting spherical models to the tangential distortion profiles of the clusters. When the models are fitted to the clusters individually, we are unable to discriminate statistically between SIS and NFW models. However when the tangential distortion profiles of the individual clusters are combined, and models fitted to the stacked profile, the SIS model is rejected at 6- and 11-sigma, respectively, for low- and high-mass bins. We also use the individual cluster NFW model fits to investigate the relationship between cluster mass (M_vir) and concentration (c_vir), finding an anti-correlation of c_vir and M_vir. The best-fit c_vir-M_vir relation is: c_vir(M_vir) propto M_vir^{-alpha} with alpha=0.41+/-0.19 -- i.e. a non-zero slope is detected at 2sigma significance. We then investigate the optimal radius within which to measure cluster mass, finding that the typical fractional errors are improved to sigma(M_Delta)/M_Delta ~ 0.1-0.2 for cluster masses at higher over-densities Delta=500-2000, from 0.2-0.3 for the virial over-density (~110). Further comparisons between mass measurements based on spherical model fitting and the model-independent aperture mass method reveal that the 2D aperture mass enclosed within a cylinder of a given aperture radius is systematically greater than the 3D spherical mass obtained from NFW model fitting: M_2D/M_3D= 1.34 and 1.40 for Delta=500 and 110, respectively. The amplitude of this effect agrees well with that predicted by integrating the NFW model along the line-of-sight.

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.

Interferometric millimeter observations of the cosmic microwave background and clusters of galaxies with arcmin resolutions require antenna arrays with short spacings. Having all antennas co-mounted on a single steerable platform sets limits to the overall weight. A 25 kg lightweight novel carbon-fiber design for a 1.2 m diameter Cassegrain antenna is presented. The finite element analysis predicts excellent structural behavior under gravity, wind and thermal load. The primary and secondary mirror surfaces are aluminum coated with a thin TiO$_2$ top layer for protection. A low beam sidelobe level is achieved with a Gaussian feed illumination pattern with edge taper, designed based on feedhorn antenna simulations and verified in a far field beam pattern measurement. A shielding baffle reduces inter-antenna coupling to below $\sim$ -135 dB. The overall antenna efficiency, including a series of efficiency factors, is estimated to be around 60%, with major losses coming from the feed spillover and secondary blocking. With this new antenna, a detection rate of about 50 clusters per year is anticipated in a 13-element array operation.

Clusters of galaxies have been used extensively to determine cosmological parameters. A major difficulty in making best use of Sunyaev-Zel'dovich (SZ) and X-ray observations of clusters for cosmology is that using X-ray observations it is difficult to measure the temperature distribution and therefore determine the density distribution in individual clusters of galaxies out to the virial radius. Observations with the new generation of SZ instruments are a promising alternative approach. We use clusters of galaxies drawn from high-resolution adaptive mesh refinement (AMR) cosmological simulations to study how well we should be able to constrain the large-scale distribution of the intra-cluster gas (ICG) in individual massive relaxed clusters using AMiBA in its configuration with 13 1.2-m diameter dishes (AMiBA13) along with X-ray observations. We show that non-isothermal beta models provide a good description of the ICG in our simulated relaxed clusters. We use simulated X-ray observations to estimate the quality of constraints on the distribution of gas density, and simulated SZ visibilities (AMiBA13 observations) for constraints on the large-scale temperature distribution of the ICG. We find that AMiBA13 visibilities should constrain the scale radius of the temperature distribution to about 50% accuracy. We conclude that the upgraded AMiBA, AMiBA13, should be a powerful instrument to constrain the large-scale distribution of the ICG.

Accurate mass determination of clusters of galaxies is crucial if they are to be used as cosmological probes. However, there are some discrepancies between cluster masses determined based on gravitational lensing, and X-ray observations assuming strict hydrostatic equilibirium (i.e., the equilibrium gas pressure is provided entirely by thermal pressure). Cosmological simulations suggest that turbulent gas motions remaining from hierarchical structure formation may provide a significant contribution to the equilibrium pressure in clusters. We analyze a sample of massive clusters of galaxies drawn from high resolution cosmological simulations, and find a significant contribution (20%-45%) from non-thermal pressure near the center of relaxed clusters, and, in accord with previous studies, a minimum contribution at about 0.1 Rvir, growing to about 30%-45% at the virial radius, Rvir. Our results strongly suggest that relaxed clusters should have significant non-thermal support in their core region. As an example, we test the validity of strict hydrostatic equilibirium in the well-studied massive galaxy cluster Abell 1689 using the latest high resolution gravitational lensing and X-ray observations. We find a contribution of about 40% from non-thermal pressure within the core region of A1689, suggesting an alternate explanation for the mass discrepancy: the strict hydrostatic equilibirium is not valid in this region.

We investigate the contamination of the Sunyaev--Zel'dovich (SZ) effect for six galaxy clusters, A1689, A1995, A2142, A2163, A2261, and A2390, observed by the Y. T. Lee Array for Microwave Background Anisotropy during 2007. With the range of baselines used, we find that the largest effect (of order 13%-50% of the central SZ flux density) comes from primary anisotropies in the cosmic microwave background and exceeds the thermal noise in all six cases. Contamination from discrete radio sources is estimated to be at a level of (3%-60%) of the central SZ flux density. We use the statistics of these contaminating sources to estimate and correct the errors in the measured SZ effects of these clusters.

We present a strong-lensing analysis of the galaxy cluster MS 1358.4+6245 ($z=0.33$), in deep 6-band ACS/HST imaging. In addition to the well-studied system at $z=4.92$, our modelling method uncovers 19 new multiply-lensed images so that a total of 23 images and their redshifts are used to accurately constrain the inner mass distribution. We derive a relatively shallow inner mass profile, $d\log \Sigma/d\log r\simeq -0.33 \pm0.05$ ($r<200$ kpc), with a much higher magnification than estimated previously by models constrained only by the $z=4.92$ system. Using these many new images we can apply a non-parametric adaptive-grid method, which also yields a shallow mass profile without prior assumptions, strengthening our conclusions. The total magnification of the $z_s=4.92$ galaxy is high, about a $\sim100\times$ over its four images, so that the inferred source size, luminosity and star-formation rate are about $\sim5\times$ smaller than previous estimates, corresponding to a dwarf-sized galaxy of radius $\simeq1$ kpc. A detailed image of the interior morphology of the source is generated with a high effective resolution of only $\simeq$50 pc, thanks to the high magnification and to the declining angular diameter distance above $z\sim1.5$ for the standard cosmology, so that this image apparently represents the best resolved object known at high redshift.

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 present a joint weak-lensing/X-ray study of galaxy cluster mass-observable scaling relations, motivated by the critical importance of accurate calibration of mass proxies for future X-ray missions, including eROSITA. We use a sample of 12 clusters at z\simeq0.2 that we have observed with Subaru and XMM-Newton to construct relationships between the weak-lensing mass (M), and three X-ray observables: gas temperature (T), gas mass (Mgas), and quasi-integrated gas pressure (Yx) at overdensities of \Delta=2500, 1000, and 500 with respect to the critical density. We find that Mgas at \Delta\le1000 appears to be the most promising mass proxy of the three, because it has the lowest intrinsic scatter in mass at fixed observable: \sigma_lnM\simeq0.1, independent of cluster dynamical state. The scatter in mass at fixed T and Yx is a factor of \sim2-3 larger than at fixed Mgas, which are indicative of the structural segregation that we find in the M-T and M-Yx relationships. Undisturbed clusters are found to be \sim40% and \sim20% more massive than disturbed clusters at fixed T and Yx respectively at \sim2\sigma significance. In particular, A1914 - a well-known merging cluster - significantly increases the scatter and lowers the the normalization of the relation for disturbed clusters. We also investigated the covariance between intrinsic scatter in M-Mgas and M-T relations, finding that they are positively correlated. This contradicts the adaptive mesh refinement simulations that motivated the idea that Yx may be a low scatter mass proxy, and agrees with more recent smoothed particle hydrodynamic simulations based on the Millennium Simulation. We also propose a method to identify a robust mass proxy based on principal component analysis. The statistical precision of our results are limited by the small sample size and the presence of the extreme merging cluster in our sample.

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 gravitational lensing of background galaxies is a unique, direct probe of the distribution of matter in clusters of galaxies. We review several important aspects of cluster weak gravitational lensing together with recent advances in weak lensing techniques for measuring cluster lensing profiles and constraining cluster structure parameters.

A wideband analog correlator has been constructed for the Yuan-Tseh Lee Array for Microwave Background Anisotropy. Lag correlators using analog multipliers provide large bandwidth and moderate frequency resolution. Broadband IF distribution, backend signal processing and control are described. Operating conditions for optimum sensitivity and linearity are discussed. From observations, a large effective bandwidth of around 10 GHz has been shown to provide sufficient sensitivity for detecting cosmic microwave background variations.

The Sunyaev-Zel'dovich Effect (SZE) has been observed toward six massive galaxy clusters, at redshifts 0.091 \leq z \leq 0.322 in the 86-102 GHz band with the Y. T. Lee Array for Microwave Background Anisotropy (AMiBA). We modify an iterative method, based on the isothermal \beta-models, to derive the electron temperature T_e, total mass M_t, gas mass M_g, and integrated Compton Y within r_2500, from the AMiBA SZE data. Non-isothermal universal temperature profile (UTP) \beta models are also considered in this paper. These results are in good agreement with those deduced from other observations. We also investigate the embedded scaling relations, due to the assumptions that have been made in the method we adopted, between these purely SZE-deduced T_e, M_t, M_g and Y. Our results suggest that cluster properties may be measurable with SZE observations alone. However, the assumptions built into the pure-SZE method bias the results of scaling relation estimations and need further study.

(Abridged) We present results of Suzaku observations of the intracluster medium (ICM) in Abell 1689, combined with complementary analysis of the SDSS data and weak and strong lensing analysis of Subaru/Suprime-Cam and HST/ACS observations. Faint X-ray emission from the ICM around the virial radius is detected at 4.0 sigma significance. We find anisotropic gas temperature and entropy distributions in cluster outskirts correlated with large-scale structure of galaxies. The high temperature and entropy region in the northeastern (NE) outskirts is connected to an overdense filamentary structure. The outskirt regions in contact with low density void environments have low gas temperatures and entropies, deviating from hydrostatic equilibrium. These results suggest that thermalization of the ICM occurs faster along the filamentary structures than the void regions. A joint X-ray and lensing analysis shows that the hydrostatic mass is $\sim60-90%$ of spherical lensing one but comparable to a triaxial halo mass within errors in $0.6r_{2500} \simlt r \simlt 0.8r_{500}$, and that it is significantly biased as low as $\simlt60%$ within $0.4r_{2500}$, irrespective of mass models. The thermal gas pressure within $r_{500}$ is, at most, $\sim50$--60% of the total pressure to balance fully the gravity of the spherical lensing mass, and $\sim30$--40% around the virial radius. Although these constitute lower limits when one considers the possible halo triaxiality, these small relative contributions of thermal pressure would require additional sources of pressure, such as bulk and/or turbulent motions.

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 investigate the scaling relations between the X-ray and the thermal Sunyaev-Zel'dovich Effect (SZE) properties of clusters of galaxies, using data taken during 2007 by the Y.T. Lee Array for Microwave Background Anisotropy (AMiBA) at 94 GHz for the six clusters A1689, A1995, A2142, A2163, A2261, and A2390. The scaling relations relate the integrated Compton-y parameter Y_{2500} to the X-ray derived gas temperature T_{e}, total mass M_{2500}, and bolometric luminosity L_X within r_{2500}. Our results for the power-law index and normalization are both consistent with the self-similar model and other studies in the literature except for the Y_{2500}-L_X relation, for which a physical explanation is given though further investigation may be still needed. Our results not only provide confidence for the AMiBA project but also support our understanding of galaxy clusters.

The Array for Microwave Background Anisotropy (AMiBA) is a radio interferometer for research in cosmology, currently operating 7 0.6m diameter antennas co-mounted on a 6m diameter platform driven by a hexapod mount. AMiBA is currently the largest hexapod telescope. We briefly summarize the hexapod operation with the current pointing error model. We then focus on the upcoming 13-element expansion with its potential difficulties and solutions. Photogrammetry measurements of the platform reveal deformations at a level which can affect the optical pointing and the receiver radio phase. In order to prepare for the 13-element upgrade, two optical telescopes are installed on the platform to correlate optical pointing tests. Being mounted on different locations, the residuals of the two sets of pointing errors show a characteristic phase and amplitude difference as a function of the platform deformation pattern. These results depend on the telescope's azimuth, elevation and polarization position. An analytical model for the deformation is derived in order to separate the local deformation induced error from the real hexapod pointing error. Similarly, we demonstrate that the deformation induced radio phase error can be reliably modeled and calibrated, which allows us to recover the ideal synthesized beam in amplitude and shape of up to 90% or more. The resulting array efficiency and its limits are discussed based on the derived errors.

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 examine the dynamics structure of the rich cluster A1689, combining VLT/VIMOS spectroscopy with Subaru/Suprime-Cam imaging. The radial velocity distribution of $\sim 500$ cluster members is bounded by a pair of clearly defined velocity caustics, with a maximum amplitude of $\sim|4000|$ km/s at $\simeq$ 300 h$^{-1}$ kpc, beyond which the amplitude steadily declines, approaching zero velocity at a limiting radius of $\sim$ 2 h$^{-1}$ Mpc. We derive the 3D velocity anisotropy and galaxy number density profiles using a model-independent method to solve the Jeans equation, simultaneously incorporating the observed velocity dispersion profile, the galaxy counts from deep Subaru imaging, and our previously derived cluster mass profile from a joint lensing and X-ray analysis. The velocity anisotropy is found to be predominantly radial at large radius, becoming increasingly tangential towards the center, in accord with expectations. We also analyze the galaxy data independently of our previous analysis using two different methods: The first is based on a solution of the Jeans equation assuming an NFW form for the mass distribution, whereas in the second method the caustic amplitude is used to determine the escape velocity. The cluster virial mass derived by both of these dynamical methods is in good agreement with results from our earlier lensing and X-ray analysis. We also confirm the high NFW concentration parameter, with results from both methods combined to yield $c_{\rm vir}>13$ (1$\sigma$). The inferred virial radius is consistent with the limiting radius where the caustics approach zero velocity and where the counts of cluster members drop off, suggesting that infall onto A1689 is currently not significant.

We present the first galaxy evolution results from the Local Cluster Substructure Survey (LoCuSS), a multi-wavelength survey of 100 X-ray selected galaxy clusters at 0.15<z<0.30. LoCuSS combines far-UV through far-IR observations of cluster galaxies with gravitational lensing analysis and X-ray data to investigate the interplay between the hierarchical assembly of clusters and the evolution of cluster galaxies. Here we present new panoramic Spitzer/MIPS 24micron observations of the merging cluster Abell 1758 at z=0.279 spanning 6.5x6.5Mpc and reaching a 90% completeness limit of 400uJy. We estimate a global cluster SFR of 910\pm320 M_sun/yr within 3 Mpc of the cluster centre, originating from 42 galaxies with L_IR > 5x10^10 L_sun. The obscured activity in A1758 is therefore comparable with that in Cl 0024+1654, the most active cluster previously studied at 24um. The obscured galaxies faithfully trace the cluster potential as revealed by the weak-lensing mass map of the cluster, including numerous mass peaks at R~2-3Mpc that are likely associated with infalling galaxy groups and filamentary structures. However the core (R<500kpc) of A1758N is 2x more active in the IR than that of A1758S, likely reflecting differences in the recent dynamical history of the two clusters. The 24micron results from A1758 therefore suggest that dust-obscured cluster galaxies are common in merging clusters and suggests that obscured activity in clusters is triggered by both the details of cluster-cluster mergers and processes that operate at larger radii including those within in-falling groups. Our ongoing far-UV through far-IR observations of a large sample of clusters should allow us to disentangle the different physical processes responsible for triggering obscured star formation in clusters.

The Y. T. Lee Array for Microwave Background (AMiBA) has reported the first science results on the detection of galaxy clusters via the Sunyaev Zel'dovich effect. The science objectives required small reflectors in order to sample large scale structures (20') while interferometry provided modest resolutions (2'). With these constraints, we designed for the best sensitivity by utilizing the maximum possible continuum bandwidth matched to the atmospheric window at 86-102GHz, with dual polarizations. A novel wide-band analog correlator was designed that is easily expandable for more interferometer elements. MMIC technology was used throughout as much as possible in order to miniaturize the components and to enhance mass production. These designs will find application in other upcoming astronomy projects. AMiBA is now in operations since 2006, and we are in the process to expand the array from 7 to 13 elements.

The Y.T. Lee Array for Microwave Background Anisotropy (AMiBA) started scientific operation in early 2007. This work describes the optimization of the system performance for the measurements of the Sunyaev-Zel'dovich effect for six massive galaxy clusters at redshifts $0.09 - 0.32$. We achieved a point source sensitivity of $63\pm 7$ mJy with the seven 0.6m dishes in 1 hour of on-source integration in 2-patch differencing observations. We measured and compensated for the delays between the antennas of our platform-mounted interferometer. Beam switching was used to cancel instrumental instabilities and ground pick up. Total power and phase stability were good on time scales of hours, and the system was shown to integrate down on equivalent timescales of 300 hours per baseline/correlation, or about 10 hours for the entire array. While the broadband correlator leads to good sensitivity, the small number of lags in the correlator resulted in poorly measured bandpass response. We corrected for this by using external calibrators (Jupiter and Saturn). Using Jupiter as the flux standard, we measured the disk brightness temperature of Saturn to be $149^{+5}_{-12}$ K.

AMiBA is the largest hexapod astronomical telescope in current operation. We present a description of this novel hexapod mount with its main mechanical components -- the support cone, universal joints, jack screws, and platform -- and outline the control system with the pointing model and the operating modes that are supported. The AMiBA hexapod mount performance is verified based on optical pointing tests and platform photogrammetry measurements. The photogrammetry results show that the deformations in the inner part of the platform are less than 120 micron rms. This is negligible for optical pointing corrections, radio alignment and radio phase errors for the currently operational 7-element compact configuration. The optical pointing error in azimuth and elevation is successively reduced by a series of corrections to about 0.4 arcmin rms which meets our goal for the 7-element target specifications.

We present observations, analysis and results for the first-year operation of AMiBA, an interferometric experiment designed to study cosmology via the measurement of Cosmic Microwave Background (CMB). AMiBA is the first CMB interferometer operating at 3 mm to have reported successful results, currently with seven close-packed antennas of 60-cm diameter giving a synthesized resolution of around 6'. During 2007 AMiBA detected the Sunyaev-Zel'dovich effects (SZE) of six galaxy clusters at redshift 0.091 <= z <= 0.322. An observing strategy with on-off-source switching is used to minimize the effects from electronic offset and ground pickup. Planets were used to test the observational capability of AMiBA and to calibrate the conversion from correlator time-lag data to visibilities. The detailed formalism for data analysis is given. We summarize our early tests including observations of planets and quasars, and present images, visibility profiles, the estimated central coordinates, sizes, and SZE amplitudes of the galaxy clusters. Science results are summarized. We also discuss possible systematic effects in the results.

The Yuan-Tseh Lee Array for Microwave Background Anisotropy (AMiBA) is the first interferometer dedicated to studying the cosmic microwave background (CMB) radiation at 3mm wavelength. The choice of 3mm was made to minimize the contributions from foreground synchrotron radiation and Galactic dust emission. The initial configuration of seven 0.6m telescopes mounted on a 6-m hexapod platform was dedicated in October 2006 on Mauna Loa, Hawaii. Scientific operations began with the detection of a number of clusters of galaxies via the thermal Sunyaev-Zel'dovich effect. We compare our data with Subaru weak lensing data in order to study the structure of dark matter. We also compare our data with X-ray data in order to derive the Hubble constant.

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 describe methods used to validate data from the Y.T. Lee Array for Microwave Background Anisotropy (AMiBA), an interferometric array designed to measure the Sunyaev-Zel'dovich effect and the anisotropy of the Cosmic Microwave Background (CMB). We perform several statistical tests on data from pointed galaxy cluster observations taken in 2007 and noise data from long-term blank sky observations and measurements with the feeds covered by the absorbers. We apply power spectrum analysis, cross power spectrum analysis among different outputs with different time lags in our analog correlator, and sample variance law tests to noise data. We find that (1) there is no time variation of electronic offsets on the time scale of our two-patch observations (~10 minutes); (2) noise is correlated by less than 10% between different lags; and (3) the variance of noise scales with the inverse of time. To test the Gaussianity of the data, we apply Kolmogorov-Smirnov (K-S) tests to cluster data, and find that a 5% significance level efficiently detects data sets with known hardware problems without rejecting an excess of acceptable data. We also calculate third- and fourth-order moments and cumulants for the noise residual visibilities and find that about 95% of our data are within the 99% confidence regions of Gaussianity.

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.

We present and compare projected distributions of mass, galaxies, and the intracluster medium (ICM) for a sample of merging clusters of galaxies based on the joint weak-lensing, optical photometric, and X-ray analysis. Our sample comprises seven nearby Abell clusters, for which we have conducted systematic, deep imaging observations with Suprime-Cam on Subaru telescope. Our seven target clusters, representing various merging stages and conditions, allow us to investigate in details the physical interplay between dark matter, ICM, and galaxies associated with cluster formation and evolution. A1750 and A1758 are binary systems consisting of two cluster-sized components, A520, A754, A1758N, A1758S, and A1914 are on-going cluster mergers, and A2034 and A2142 are cold-front clusters. In the binary clusters, the projected mass, optical light, and X-ray distributions are overall similar and regular without significant substructures. On-going and cold-front merging clusters, on the other hand, reveal highly irregular mass distributions. Overall the mass distribution appears to be similar to the galaxy luminosity distribution, whereas their distributions are quite different from the ICM distribution in a various ways. We also measured for individual targets the global cluster parameters such as the cluster mass,the mass-to-light ratio, and the ICM temperature. A comparison of the ICM and virial temperatures of merging clusters from X-ray and weak-lensing analyses, respectively, shows that the ICM temperature of on-going and cold-front clusters is significantly higher than the cluster virial temperature by a factor of $\sim 2$. This temperature excess in the ICM could be explained by the effects of merger boosts.

[Abridged] We derive a projected 2D mass map of the well studied galaxy cluster A1689 based on an entropy-regularized maximum-likelihood combination of the lens magnification and distortion of red background galaxies registered in deep Subaru images. The known strong lensing information is also readily incorporated in this approach, represented as a central pixel with a mean surface density close to the critical value. We also utilize the distortion measurements to locally downweight the intrinsic clustering noise, which otherwise perturbs the depletion signal. The projected mass profile continuously steepens with radius and is well fitted by the NFW model, but with a surprising large concentration c_vir=13.4^{+5.3}_{-3.3}, lying far from the predicted value of c_vir ~ 5, corresponding to the measured virial mass, M_vir=(2.1\pm 0.2)\times 10^15M_{\odot}, posing a challenge to the standard assumptions defining the LCDM model. We examine the consistency of our results with estimates derived with the standard weak lensing estimators and by comparison with the inner mass profile obtained from strong lensing. All the reconstructions tested here imply a virial mass in the range, M_vir=(1.5-2.1)\times 10^15M_{\odot}, and the combined ACS and Subaru-2D mass reconstruction yields a tight constraint on the concentration parameter, c_vir=12.7 \pm 1(stat.) \pm 2.8(systematic) (c_200 ~ 10), improving upon the statistical accuracy of our earlier 1D analysis. Importantly, our best fitting profile properly reproduces the observed Einstein radius of 45 arcsec (z_s=1), in contrast to other weak lensing work, reporting lower concentration profiles, which underestimate the observed Einstein radius.

We have developed a method for measuring higher-order weak lensing distortions of faint background galaxies, namely the weak gravitational flexion, by fully extending the Kaiser, Squires & Broadhurst method to include higher-order lensing image characteristics (HOLICs) introduced by Okura, Umetsu, & Futamase. We take into account explicitly the weight function in calculations of noisy shape moments and the effect of higher-order PSF anisotropy, as well as isotropic PSF smearing. Our HOLICs formalism allows accurate measurements of flexion from practical observational data in the presence of non-circular, anisotropic PSF. We test our method using mock observations of simulated galaxy images and actual, ground-based Subaru observations of the massive galaxy cluster A1689 ($z=0.183$). From the high-precision measurements of spin-1 first flexion, we obtain a high-resolution mass map in the central region of A1689. The reconstructed mass map shows a bimodal feature in the central $4'\times 4'$ region of the cluster. The major, pronounced peak is associated with the brightest cluster galaxy and central cluster members, while the secondary mass peak is associated with a local concentration of bright galaxies. The refined, high-resolution mass map of A1689 demonstrates the power of the generalized weak lensing analysis techniques for quantitative and accurate measurements of the weak gravitational lensing signal.

We present results from a weak lensing analysis of the galaxy cluster A1689 (z=0.183) based on deep wide-field imaging data taken with Suprime-Cam on Subaru telescope. A maximum entropy method has been used to reconstruct directly the projected mass distribution of A1689 from combined lensing distortion and magnification measurements of red background galaxies.The resulting mass distribution is clearly concentrated around the cD galaxy, and mass and light in the cluster are similarly distributed in terms of shape and orientation. The azimuthally-averaged mass profile from the two-dimensional reconstruction is in good agreement with the earlier results from the Subaru one-dimensional analysis of the weak lensing data, supporting the assumption of quasi-circular symmetry in the projected mass distribution of the cluster.

We present results from a combined weak lensing and X-ray analysis of the merging cluster A1914 at a redshift of $z=0.1712$ based on $R_{\rm}$-band imaging data with Subaru/Suprime-Cam and archival Chandra X-ray data. Based on the weak-lensing and X-ray data we explore the relationships between cluster global properties, namely the gravitational mass, the bolometric X-lay luminosity and temperature, the gas mass and the gas mass fraction, as a function of radius. We found that the gas mass fractions within $r_{2500}$ and $r_{{\rm vir}}$ are consistent with the results of earlier X-ray cluster studies and of cosmic microwave background studies based on the WMAP observations, respectively. However, the observed temperature, $k_B T_{\rm ave}=9.6\pm0.3 {\rm keV}$, is significantly higher than the virial temperature, $k_B T_{\rm vir}=4.7\pm0.3 {\rm keV}$ derived from the weak lensing distortion measurement. The X-ray bolometric luminosity and temperature ($L_X-T$) relation is consistent with the $L_X-T$ relation derived by previous statistical X-ray studies of galaxy clusters. Such correlations among the global cluster properties are invaluable observational tools for studying the cluster merger physics. Our results demonstrate that the combination of X-ray and weak-lensing observations is a promising, powerful probe of the physical processes associated with cluster mergers as well as of their mass properties.

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.

We study a possibility to use the octopole moment of gravitationally lensed images as a direct measure of the third-order weak gravitational lensing effect, or the gravitational flexion. It turns out that there is a natural relation between flexion and certain combinations of octopole/higher-multipole moments which we call the Higher Order Lensing Image's Characteristics (HOLICs). This will allow one to measure directly flexion from observable octopole and higher-multipole moments of background images. We show based on simulated observations how the use of HOLICs can improve the accuracy and resolution of a reconstructed mass map, in which we assume Gaussian uncertainties in the shape measurements estimated using deep i'-band data of blank fields observed with Suprime-Cam on the Subaru telescope.

The technique of weak-lensing aperture mass densitometry, so called the zeta-statistic, has recently been popular in actual observations for measurement of individual cluster mass. It has however been anticipated that the line-of-sight projection by foreground and background matter can adversely affect the cluster mass determination with not only substantial error dispersion but also a sizable positive systematic bias. Additionally, the finite number of background galaxies even at a reasonable observing depth can also introduce Poisson noise to the mass estimate. In this paper, we quantitatively investigate the degree of errors separately contributed by the two sources to the mass determination of those galaxy clusters with M_{200}>10^{14}M_{\odot}. We find that the aperture mass of zeta-statistic turns out to be a mass estimator of much reduced systematic bias, due to the cancellation by the positively biased local background mass sheet. However, the error dispersion of M_{200} arising from both projection effect and Poisson noise is found to be still sizable (40%-90%), even for the shear-selected, clean sample where multiple clusters located within a suitable projected aperture are removed. We also investigate how to remedy this large-error problem in weak lensing measurements, and propose a plausible alternative mass estimator, M(<theta_{1000}), an aperture mass measured within about half the virial radius. The aperture mass M(<theta_{1000}) is free of bias and has a substantially reduced error dispersion, 39% for the worst case of high-z, low-mass clusters, that can be smaller than the error dispersion of M_{200} as much as a factor 3.

We present a weak-lensing analysis of a region around the galaxy cluster Cl 1604+4304 (z=0.897) on the basis of the deep observations with the HST/WFPC2. We apply a variant of Schneider's aperture mass technique to the observed WFPC2 field and obtain the distribution of weak-lensing signal-to-noise (S/N) ratio within the field. The resulting S/N map reveals a clear pronounced peak located about 1.7 arcmin (850h_{50}^{-1} kpc at z=0.897) southwest of the second peak associated with the optical cluster center determined from the dynamical analysis of Postman et al. A non-linear finite-field inversion method has been used to reconstruct the projected mass distribution from the observed shear field. The reconstructed mass map shows a super-critical feature at the location of the S/N peak as well as in the cluster central region. Assuming the redshift distribution of field galaxies, we obtain the total mass in the observed field to be 1.0 h_{50}^{-1} 10^{15} M_sun for <z>=1.0. The estimated mass within a circular aperture of radius 280h_{50}^{-1} kpc centered on the dark clump is 2.4h_{50}^{-1} 10^{14} M_sun. We have confirmed the existence of the ` dark ' mass concentration from another deep HST observation with a slightly different ~20 arcsec pointing.

We present simulations of interferometric Sunyaev-Zel'dovich effect (SZE) and optical weak lenisng observations for the forthcoming AMiBA experiment, aiming at searching for high-redshift clusters of galaxies. On the basis of simulated sky maps, we have derived theoretical halo number counts and redshift distributions of selected halo samples for an AMiBA SZE survey and a weak lensing follow-up survey. By utilizing the conditional number counts of weak lensing halos with the faint SZE detection, we show that a combined SZE and weak lensing survey can gain an additional fainter halo sample at a given false positive rate, which cannot be obtained from either survey alone.

We report the discovery and spectroscopic confirmation of three gravitationally-lensed images of a galaxy at z=3.9 in the background of a distant, rich cluster of galaxies at z=0.83, on the basis of observations with Faint Object Camera And Spectrograph (FOCAS) on the Subaru telescope. We construct a simple lens model of the cluster mass distribution based on Jee et al.'s weak lensing mass estimates from deep, high-resolution images by Advanced Camera for Surveys (ACS) on the Hubble Space Telescope. The lens model can account simultaneously for the observed image configuration and the flux ratio of the closer pair located close to the critical curve. The parities of the three images are also consistent with the lensing hypothesis. Since this galaxy is apparently bright (i'_{AB} \sim 23.7) for its redshift due to the magnification, it serves as a good high redshift target on which we can make extensive and detailed studies based on multi-wavelength imaging and spectroscopy.

The steep mass profile of A1689 derived from recent detailed lensing observations is not readily reconciled with the low concentration halos predicted by the standard CDM model. However, halo triaxiality may act to bias the profile constraints derived assuming a spherically symmetric mass distribution, since lensing relates only to the projected mass distribution. A degree of halo triaxiality is inherent to the CDM structure formation, arising from the collision-less nature of the dark matter. Here we compare the CDM-based model predictions of triaxial halo with the precise lensing measurements of A1689 based on the ACS/HST and Subaru data, over a wide range of 10kpc<r<2Mpc. The model lensing profiles cover the intrinsic spread of halo mass and shape, and are projected over all inclinations when comparing with the data. We show that the model parameters are only weakly constrained and strongly degenerate mainly because of the lack of information along the line of sight. In particular, the limits on the concentration parameter become less restrictive with increasing triaxiality. Yet, by comparing the obtained constraints with expected probability distributions for the axis ratio and concentration parameters computed from numerical simulations, we find that ~6% of cluster-size halos in the CDM model can match the A1689 lensing observations at the 2-sigma level, corresponding to cases where the major-axis of the halo is aligned with the line of sight. Thus halo triaxiality could reduce the apparent discrepancy between theory and observation. This hypothesis needs to be further explored by a statistical lensing study for other clusters as well as by complementary three-dimensional information derived using X-ray, kinematics, and SZ effect observations.

We discuss the environmental dependence of galaxy evolution based on deep panoramic imaging of two distant clusters taken with Suprime-Cam as part of the PISCES project. By combining with the SDSS data as a local counterpart for comparison, we construct a large sample of galaxies that spans wide ranges in environment, time, and stellar mass (or luminosity). We find that colours of galaxies, especially those of faint galaxies ($M_V>M_V^*+1$), change from blue to red at a break density as we go to denser regions. Based on local and global densities of galaxies, we classify three environments: field, groups, and clusters. We show that the cluster colour-magnitude relation is already built at $z=0.83$. In contrast to this, the bright-end of the field colour-magnitude relation has been vigorously built all the way down to the present-day and the build-up at the faint-end has not started yet. A possible interpretation of these results is that galaxies evolve in the 'down-sizing' fashion. That is, massive galaxies complete their star formation first and the truncation of star formation is propagated to smaller objects as time progresses. This trend is likely to depend on environment since the build-up of the colour-magnitude relation is delayed in lower-density environments. Therefore, we may suggest that the evolution of galaxies took place earliest in massive galaxies and in high density regions, and it is delayed in less massive galaxies and in lower density regions.

The cross-correlation of Sunyaev-Zel'dovich effect (SZ) and weak-lensing imaging surveys can be used to test how well hot baryons trace dark matter in clusters of galaxies. We examine this concept using mock SZ and weak-lensing surveys based on the forthcoming AMiBA experiment and generated from a pre-heated cosmological N-body/hydrodynamic simulation. A cross-correlation diagram between matched lensing convergence and Compton y maps exhibits butterfly-wing-like structures, corresponding to individual clusters, that encode rich information about the distributions of hot gas and dark matter. When the cluster redshift and temperature are available the slope of a wing reveals the cluster gas fraction and the width of the wing indicates how badly the hot gas traces dark matter. On the basis of simulated data we discuss systematic errors in the projected gas fraction estimates that would be obtained from such survey comparisons.

Subaru observations of A1689 (z=0.183) are used to derive an accurate, model-independent mass profile for the entire cluster, r<2 Mpc/h, by combining magnification bias and distortion measurements. The projected mass profile steepens quickly with increasing radius, falling away to zero at r~1.0 Mpc/h, well short of the anticipated virial radius. Our profile accurately matches onto the inner profile, r<200 kpc/h, derived from deep HST/ACS images. The combined ACS and Subaru information is well fitted by an NFW profile with virial mass, (1.93 \pm 0.20)10^15 M_sun, and surprisingly high concentration, c_vir=13.7^{+1.4}_{-1.1}, significantly larger than theoretically expected (c_vir~4), corresponding to a relatively steep overall profile. A slightly better fit is achieved with a steep power-law model that has its 2D logarithmic slope -3 and core radius theta_c~1.7' (r_c~210 kpc/h), whereas an isothermal profile is strongly rejected. These results are based on a reliable sample of background galaxies selected to be redder than the cluster E/S0 sequence. By including the faint blue galaxy population a much smaller distortion signal is found, demonstrating that blue cluster members significantly dilute the true signal for r~400 kpc/h. This contamination is likely to affect most weak lensing results to date.

We have simulated the interferometric observation of the Cosmic Microwave Background (CMB) temperature and polarization fluctuations. We have constructed data pipelines from the time-ordered raw visibility samples to the CMB power spectra which utilize the methods of data compression, maximum likelihood analysis, and optimal subspace filtering. They are customized for three observational strategies, such as the single pointing, the mosaicking, and the drift-scanning. For each strategy, derived are the optimal strategy parameters that yield band power estimates with minimum uncertainty. The results are general and can be applied to any close-packed array on a single platform such as the CBI and the forthcoming AMiBA experiments. We have also studied the effect of rotation of the array platform on the band power correlation by simulating the CBI single pointing observation. It is found that the band power anti-correlations can be reduced by rotating the platform and thus densely sampling the visibility plane. This enables us to increase the resolution of the power spectrum in the l-space down to the limit of the sampling theorem (Delta l = 226 = pi / theta), which is narrower by a factor of about sqrt{2} than the resolution limit (Delta l = 300) used in the recent CBI single pointing observation. The validity of this idea is demonstrated for a two-element interferometer that samples visibilities uniformly in the uv-annulus. From the fact that the visibilities are the Fourier modes of the CMB field convolved with the beam, a fast unbiased estimator (FUE) of the CMB power spectra is developed and tested. It is shown that the FUE gives results very close to those from the quadratic estimator method without requiring large computer resources even though uncertainties in the results increase.

We report the first measurement of genus curves for the two-dimensional mass map in the neighborhood of rich, X-ray luminous galaxy cluster MS1054-03 at z=0.83, reconstructed from weak lensing data obtained by Suprime-Cam on the prime focus of 8.2m Subaru telescope. We find that the genus curve measured in the whole survey field deviates from that expected from a random Gaussian field. We show that this non-Gaussianity is induced by the rich cluster in this region, and that the genus curve for the region without the cluster is consistent with the prediction for a random Gaussian field. The results demonstrate the feasibility of the genus statistics to examine the non-Gaussianity due to the large scale structures and to probe the statistical properties of the large scale clustering.

The tidal gravitational field of galaxy clusters causes a coherent distortion of the images of background sources. Since the distribution of image distortions, namely the shear field, traces the local gravitational potential of a deflector, it can be used to reconstruct the two-dimensional mass distribution of clusters of galaxies. Moreover, the shear field can provide unique information on the redshift distribution of high-redshift galaxies. In this review we summarize recently-developed parameter-free methods of cluster-mass reconstruction based on the shear field, and we apply a mass-reconstruction method to the cluster Abell 370 at redshift 0.375.

We have shown that there exists low-frequency growing modes driven by a global temperature gradient in electron and ion plasmas, by linear perturbation analysis within the frame work of plasma Kinetic theory. The driving force of the instability is the local deviation of the distribution function from the Maxwell-Boltzmann due to global temperature gradient. Application to the intracluster medium shows that scattering of the particles due to waves excited by the instability is possible to reduce mean free paths of electron and ion down to five to seven order of magnitude than the mean free paths due to Coulomb collisions. This may provide a hint to explain why hot and cool gas can co-exist in the intracluster medium in spite of the very short evaporation time scale due to thermal conduction if the conductivity is the classical Spitzer value. Our results suggest that the realization of the global thermal equilibrium is postponed by the local instability which is induced for quicker realization of local thermal equilibrium state in plasmas. The instability provides a new possibility to create and grow cosmic magnetic fields without any seed magnetic field.