Source author record

Keiichi Umetsu

Keiichi Umetsu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

99works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

99 published item(s)

preprint2025arXiv

The SRG/eROSITA All-Sky Survey. The Weak-Lensing Mass Calibration and the Stellar Mass-to-Halo Mass Relation from the Hyper Suprime-Cam Subaru Strategic Program

We present the weak-lensing mass calibration and constrain the BCG (brightest cluster galaxy) stellar-mass-to-halo-mass-and-redshift ($M_{\star,\mathrm{BCG}}-M-z$) relation for a sample of $124$ galaxy clusters and groups at redshift $0.1<z<0.8$ from the first Data Release of the $eROSITA$ All-Sky Survey (eRASS1), using data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program. The cluster survey is conducted by the $eROSITA$ X-ray telescope aboard the Spectrum-Roentgen-Gamma (SRG) space observatory. The cluster sample is X-ray-selected and optically confirmed with a negligibly low contamination rate ($\approx5%$). On a basis of individual clusters, the shear profiles of $96$ clusters are derived using the HSC Three-Year (HSC-Y3) weak-lensing data, while the BCG stellar masses of $101$ clusters are estimated using the SED template fitting to the HSC five-band ($grizY$) photometry. The observed X-ray photon count rate is used as the mass proxy, based on which individual halo masses are obtained at the given count rate in a population modelling while accounting for systematic uncertainties in the weak-lensing modelling through a simulation-calibrated weak-lensing mass-to-halo-mass relation. The count rate and BCG stellar mass relations are simultaneously constrained in a forward and population modelling. In agreement with the results based on the weak-lensing data from the DES and KiDS surveys, we obtain a count rate relation with a self-similar redshift scaling and a mass trend that is steeper than the self-similar prediction. Our results suggest that the BCG stellar mass at a fixed halo mass has remained stable with a moderate increase at a level of $\left(20\pm8\right)%$ since redshift $z\approx0.8$. This finding supports the picture of the ``rapid-then-slow'' BCG formation, where the majority of the stellar mass must have been assembled at much earlier cosmic time.

preprint2023arXiv

Exploring the stellar populations of backsplash galaxies

Backsplash galaxies are those that traverse and overshoot cluster cores as they fall into these structures. They are affected by environment, and should stand out in contrast to the infalling population. We target galaxies in the vicinity of clusters (R>R200) and select a sample in projected phase space (PPS), from the compilation of Sampaio et al. based on SDSS data. We present a statistical analysis, comparing two regions in PPS, with the same projected distance to the cluster but different velocity. The analysis relies on the presence of variations in the stellar population content of backsplash galaxies. We find a lower limit in the fractional contribution of ~5% with respect to the general sample of infalling galaxies at similar group-centric distance when using single line strength analysis, or ~15-30% when adopting bivariate distributions. The stellar populations show a subtle but significant difference towards older ages, and a higher fraction of quiescent galaxies. We also compare this set with a general field sample, where a substantially larger difference in galaxy properties is found, with the field sample being consistently younger, metal poorer and with a lower fraction of quiescent galaxies. Noting that our "cluster" sample is located outside of the virial radius, we expect this difference to be caused by pre-processing of the infalling galaxies in the overall higher density regions.

preprint2022arXiv

HSC-XXL : Baryon budget of the 136 XXL Groups and Clusters

We present our determination of the baryon budget for an X-ray-selected XXL sample of 136 galaxy groups and clusters spanning nearly two orders of magnitude in mass ($M_{500}\sim 10^{13}-10^{15}M_\odot$) and the redshift range $0< z < 1$. Our joint analysis is based on the combination of HSC-SSP weak-lensing mass measurements, XXL X-ray gas mass measurements, and HSC and SDSS multiband photometry. We carry out a Bayesian analysis of multivariate mass-scaling relations of gas mass, galaxy stellar mass, stellar mass of brightest cluster galaxies (BCGs), and soft-band X-ray luminosity, by taking into account the intrinsic covariance between cluster properties, selection effect, weak-lensing mass calibration, and observational error covariance matrix. The mass-dependent slope of the gas mass--total mass ($M_{500}$) relation is found to be $1.29_{-0.10}^{+0.16}$, which is steeper than the self-similar prediction of unity, whereas the slope of the stellar mass--total mass relation is shallower than unity, $0.85_{-0.09}^{+0.12}$. The BCG stellar mass weakly depends on cluster mass with a slope of $0.49_{-0.10}^{+0.11}$. The baryon, gas mass, and stellar mass fractions as a function of $M_{500}$ agree with the results from numerical simulations and previous observations. We successfully constrain the full intrinsic covariance of the baryonic contents. The BCG stellar mass shows the larger intrinsic scatter at a given halo total mass, followed in order by stellar mass and gas mass. We find a significant positive intrinsic correlation coefficient between total (and satellite) stellar mass and BCG stellar mass and no evidence for intrinsic correlation between gas mass and stellar mass. All the baryonic components show no redshift evolution.

preprint2022arXiv

Line-of-sight Elongation and Hydrostatic Mass Bias of the Frontier Fields Galaxy Cluster Abell 370

We present a detailed weak-lensing and X-ray study of the Frontier Fields galaxy cluster Abell 370, one of the most massive known lenses on the sky, using wide-field BRz Subaru/Sprime-Cam and Chandra X-ray observations. By combining two-dimensional (2D) shear and azimuthally averaged magnification constraints derived from Subaru data, we perform a lensing mass reconstruction in a free-form manner, which allows us to determine both radial structure and 2D morphology of the cluster mass distribution. In a triaxial framework assuming a Navarro-Frenk-White density profile, we constrain the intrinsic structure and geometry of the cluster halo by forward modeling the reconstructed mass map. We obtain a halo mass $M_{200}=(1.54 \pm 0.29)\times 10^{15}h^{-1}M_\odot$, a halo concentration $c_{200}=5.27 \pm 1.28$, and a minor-major axis ratio $q_a=0.62 \pm 0.23$ with uninformative priors. Using a prior on the line-of-sight alignment of the halo major axis derived from binary merger simulations constrained by multi-probe observations, we find that the data favor a more prolate geometry with lower mass and lower concentration. From triaxial lens modeling with the line-of-sight prior, we find a spherically enclosed gas mass fraction of $f_\mathrm{gas}=(8.4 \pm 1.0)\%$ at $0.7h^{-1}$ Mpc $\sim 0.7r_{500}$. When compared to the hydrostatic mass estimate from Chandra observations, our triaxial weak-lensing analysis yields spherically enclosed mass ratios of $1-b \equiv M_\mathrm{HE}/M_\mathrm{WL} = 0.56 \pm 0.09$ and $0.51 \pm 0.09$ at $0.7h^{-1}$ Mpc with and without using the line-of-sight prior, respectively. Since the cluster is in a highly disturbed dynamical state, this represents the likely maximum level of hydrostatic bias in galaxy clusters.

preprint2022arXiv

MUSUBI (MegaCam Ultra-deep Survey: $u^\ast$-Band Imaging)-Data for the COSMOS and SXDS Fields

The Subaru Hyper Suprime-Cam (HSC) Strategic Survey is the latest-generation multi-band optical imaging survey for galaxy evolution and structure formation. The "Ultra-Deep" component of the HSC survey provides $grizy$ broad-band images over $\sim3.4$ deg$^2$ to detection limits of $\sim26$-28 AB, along with narrow-band images, in the COSMOS and the SXDS fields. These images provide an unprecedented combination of depths and area coverage, for the studies galaxies up to $z\sim7$. However, the lack of coverage at $<4000$ Ang implies incomplete sampling of the rest-frame UV at $z\lesssim 3$, which is critically needed for understanding the buildup of stellar mass in the later cosmic time. We conducted a multi-year CFHT $u^\ast$-band imaging campaign in the two HSC Ultra-Deep fields with CFHT MegaCam. By including shallower archival data, we reach 5-$σ$ depths of $u^\ast=28.1$ and 28.4 (AB) at the centers of the COSMOS and SXDS fields, respectively, and $u^\ast=27.7$ and 27.8 in the central 1 deg$^2$ fields. The image quality is $\gtrsim0.90$ arcsec, fairly good for the $u^\ast$ band. Both the photometric and astrometric quality of our data are excellent. We show that the combination of our $u^\ast$-band and HSC data can lead to high-quality photometric redshifts at $z=0$-3, and robust measurements of rest-frame UV on galaxies at $0.4<z<0.6$ for distinguishing green-valley galaxies from star-forming and quiescent galaxies. We publicly release our reduced $u^\ast$-band images and reference catalogs that can be used readily for scientific studies.

preprint2022arXiv

Star Formation Properties of Sloan Digital Sky Survey BOSS Void Galaxies in the Hyper Suprime-Cam Survey

We utilize the Hyper Suprime-Cam (HSC) Wide Survey to explore the properties of galaxies located in the voids identified from the Baryon Oscillation Spectroscopic Survey (BOSS) up to z~0.7. The HSC reaches i~25, allowing us to characterize the void galaxies down to 10$^{9.2}$ solar mass. We find that the revised void galaxy densities, when including faint galaxies in voids defined by bright galaxies, are still underdense compared to the mean density from the entire field. In addition, we classify galaxies into star-forming, quiescent, and green valley populations, and find that void galaxies tend to have slightly higher fractions of star-forming galaxies under the mass and redshift control, although the significance of this result is only moderate (2$σ$). However, when we focus on the star-forming population, the distribution of the specific star formation rate (sSFR) of void galaxies shows little difference from that of the control galaxies. Similarly, the median sSFR of star-forming void galaxies is also in good agreement with that of the star-forming control galaxies. Moreover, the effective green valley fraction of void galaxies, defined as the number of green valley galaxies over the number of nonquiescent galaxies, is comparable to that of the control ones, supporting the suggestion that void and control galaxies evolve under similar physical processes and quenching frequencies. Our results thus favor a scenario of the galaxy assembly bias.

preprint2022arXiv

Testing non-local gravity by clusters of galaxies

Extended theories of gravity have been extensively investigated during the last thirty years, aiming at fixing infrared and ultraviolet shortcomings of General Relativity and of the associated $Λ$CDM cosmological model. Recently, non-local theories of gravity have drawn increasing attention due to their potential to ameliorate both the ultraviolet and infrared behavior of gravitational interaction. In particular, Integral Kernel theories of Gravity provide a viable mechanism to explain the late time cosmic acceleration so as to avoid the introduction of any form of unknown dark energy. On the other hand, these models represent a natural link towards quantum gravity. Here, we study a scalar-tensor equivalent model of General Relativity corrected with non-local terms, where corrections are selected by the existence of Noether symmetries. After performing the weak field limit and generalizing the results to extended mass distributions, we analyse the non-local model at galaxy cluster scales, by comparing the theoretical predictions with gravitational lensing observations from the CLASH program. We obtain agreement with data at the same level of statistical significance as General Relativity. We also provide constraints for the Navarro--Frenk--White parameters and lower bounds for the non-local length scales. The results are finally compared with those from the literature.

preprint2022arXiv

The Three Hundred Project: dissecting the fundamental plane of galaxy clusters up to $z=1$

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

preprint2021arXiv

Likelihood-free Forward Modeling for Cluster Weak Lensing and Cosmology

Likelihood-free inference provides a rigorous approach to preform Bayesian analysis using forward simulations only. The main advantage of likelihood-free methods is its ability to account for complex physical processes and observational effects in forward simulations. Here we explore the potential of likelihood-free forward modeling for Bayesian cosmological inference using the redshift evolution of the cluster abundance combined with weak-lensing mass calibration. We use two complementary likelihood-free methods, namely Approximate Bayesian Computation (ABC) and Density-Estimation Likelihood-Free Inference (DELFI), to develop an analysis procedure for inference of the cosmological parameters $(Ω_\mathrm{m},σ_8)$ and the mass scale of the survey sample. Adopting an eROSITA-like selection function and a 10-percent scatter in the observable-mass relation in a flat $Λ$CDM cosmology with $Ω_\mathrm{m}=0.286$ and $σ_8=0.82$, we create a synthetic catalog of observable-selected NFW clusters in a survey area of 50 deg$^2$. The stacked tangential shear profile and the number counts in redshift bins are used as summary statistics for both methods. By performing a series of forward simulations, we obtain convergent solutions for the posterior distribution from both methods. We find that ABC recovers broader posteriors than DELFI, especially for the $Ω_\mathrm{m}$ parameter. For a weak-lensing survey with a source density of $n_\mathrm{g}=20$ arcmin$^{-2}$, we obtain posterior constraints on $S_8=σ_8(Ω_\mathrm{m}/0.3)^{0.3}$ of $0.836 \pm 0.032$ and $0.810 \pm 0.019$ from ABC and DELFI, respectively. The analysis framework developed in this study will be particularly powerful for cosmological inference with ongoing cluster cosmology programs, such as the XMM-XXL survey and the eROSITA all-sky survey, in combination with wide-field weak-lensing surveys.

preprint2021arXiv

Multi-component DHOST analysis in galaxy clusters

Extended Theories of Gravity with additional scalar degrees of freedom have recently acquired increasing interest due to the presence of a screening mechanism that allows suppressing at small scales (e.g., the Solar System scale) every modification restoring General Relativity. In this work, we consider a second-order Extended Theory of Gravity belonging to the family of Degenerate High Order Scalar Tensor theories (DHOST) characterized by a partial breaking of the Vainshtein screening mechanism. We study this model in two different scenarios as a description of dark energy only and as a description of both dark matter and dark energy. Such scenarios have been tested here by analysing a sample of 16 high-mass galaxy clusters targeted by the Cluster Lensing and Supernova survey with Hubble (CLASH) program using two complementary probes, namely X-ray and strong-and-weak gravitational lensing observations. In mass modelling, we adopt a multi-component approach including hot gas and galactic stellar contributions. For the majority of the clusters in our sample, results show mild Bayesian evidence in favour of the DHOST model as a description of dark energy over General Relativity. This model also appears to alleviate the discrepancy present in General Relativity between X-ray hydrostatic and lensing mass estimates. For the second scenario where gravity acts as both dark energy and dark matter due to the partial breaking of the Vainshtein screening mechanism at cluster scales, the model is statistically disfavoured compared to General Relativity.

preprint2021arXiv

Probing Vainsthein-screening gravity with galaxy clusters using internal kinematics and strong and weak lensing

We use high-precision combined strong/weak lensing and kinematics measurements of the total mass profiles of the observed galaxy clusters MACS~J1206.2-0847 and Abell~S1063, to constrain the relativistic sector of the general DHOST dark energy theories, which exhibit a partial breaking of the so called Vainsthein screening mechanism, on the linear level of scalar fluctuations around a cosmological background. In particular, by using the \textsc{MG-MAMMPOSSt} framework developed in Pizzuti et al. 2021, for the kinematics analysis of member galaxies in clusters along with lensing mass profile reconstructions, we provide new constraints on the coupling $Y_2$ which governs the theory's relativistic contribution to the lensing potential. The new bound from the combination of kinematics and lensing measurements of MACS 1206, $Y_2=-0.12^{+0.66}_{-0.67}$ at $2σ$, provides about a 2-fold improvement on previous constraints. In the case of Abell~S1063 a $>2σ$ tension with the GR expectation arises. We discuss this in some detail, and we investigate the possible sources of systematics which can explain the tension. We further discuss why the combination of kinematics of member galaxies with lensing is capable of providing much tighter bounds compared to kinematics or lensing alone, and we explain how the number density profile of tracers, as well as the choice of the velocity anisotropy profile affects the final results.

preprint2021arXiv

The eROSITA Final Equatorial-Depth Survey (eFEDS): X-ray Observable-to-Mass-and-Redshift Relations of Galaxy Clusters and Groups with Weak-Lensing Mass Calibration from the Hyper Suprime-Cam Subaru Strategic Program Survey

We present the first weak-lensing mass calibration and X-ray scaling relations of galaxy clusters and groups selected in the $eROSITA$ Final Equatorial Depth Survey (eFEDS) observed by Spectrum Roentgen Gamma/$eROSITA$ over a contiguous footprint with an area of $\approx140$ deg$^2$, using the three-year (S19A) weak-lensing data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. In this work, a sample of $434$ optically confirmed galaxy clusters (and groups) at redshift $0.01\lesssim z \lesssim1.3$ with a median of $0.35$ is studied, of which $313$ systems are uniformly covered by the HSC survey to enable the extraction of the weak-lensing shear observable. In a Bayesian population modelling, we perform a blind analysis for the weak-lensing mass calibration by simultaneously modelling the observed count rate $η$ and the shear profile $g$ of individual clusters through the count rate-to-mass-and-redshift ($η$--$M_{500}$--$z$) and weak-lensing mass-to-mass-and-redshift ($M_{\mathrm{WL}}$--$M_{500}$--$z$) relations, respectively, while accounting for the bias in these observables using simulation-based calibrations. As a result, the count rate-inferred and lensing-calibrated cluster mass is obtained from the joint modelling of the scaling relations, as the ensemble mass spanning a range of $10^{13}h^{-1}M_{\odot}\lesssim M_{500}\lesssim10^{15} h^{-1}M_{\odot}$ with a median of $\approx10^{14} h^{-1}M_{\odot}$ for the eFEDS sample. With the mass calibration, we further model the X-ray observable-to-mass-and-redshift relations, including the rest-frame soft-band and bolometric luminosity ($L_{\mathrm{X}}$ and $L_{\mathrm{b}}$), the emission-weighted temperature $T_{\mathrm{X}}$, the mass of intra-cluster medium $M_{\mathrm{g}}$, and the mass proxy $Y_{\mathrm{X}}$, which is the product of $T_{\mathrm{X}}$ and $M_{\mathrm{g}}$. (abridged)

preprint2021arXiv

The XXL survey: XLVI. Forward cosmological analysis of the C1 cluster sample

We present the forward cosmological analysis of an $XMM$ selected sample of galaxy clusters out to a redshift of unity. Following our previous 2018 study based on the dn/dz quantity alone, we perform an upgraded cosmological analysis of the same XXL C1 cluster catalogue (178 objects), with a detailed account of the systematic errors. We follow the ASpiX methodology: the distribution of the observed X-ray properties of the cluster population is analysed in a 3D observable space (count rate, hardness ratio, redshift) and modelled as a function of cosmology. Compared to more traditional methods, ASpiX allows the inclusion of clusters down to a few tens of photons. We obtain an improvement by a factor of 2 compared to the previous analysis by letting the normalisation of the M-T relation and the evolution of the L-T relation free. Adding constraints from the XXL cluster 2-point correlation function and the BAO from various surveys decreases the uncertainties by 23 and 53 % respectively, and 62% when adding both. Switching to the scaling relations from the Subaru analysis, and letting free more parameters, our final constraints are $σ_8$ = $0.99^{+0.14}_{-0.23}$, $Ω_m$ = 0.296 $\pm$ 0.034 ($S_8 = 0.98^{+0.11}_{-0.21}$) for the XXL sample alone. Finally, we combine XXL ASpiX, the XXL cluster 2-point correlation function and the BAO, with 11 free parameters, allowing for the cosmological dependence of the scaling relations in the fit. We find $σ_8$ = $0.793^{+0.063}_{-0.12}$, $Ω_m$ = 0.364 $\pm$ 0.015 ($S_8 = 0.872^{+0.068}_{-0.12}$), but still compatible with Planck CMB at 2.2$σ$. The results obtained by the ASpiX method are promising; further improvement is expected from the final XXL cosmological analysis involving a cluster sample twice as large. Such a study paves the way for the analysis of the eROSITA and future Athena surveys.

preprint2020arXiv

A clustering-based self-calibration of the richness-to-mass relation of CAMIRA galaxy clusters out to $z\approx1.1$ in the Hyper Suprime-Cam survey

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. These correlation functions are $ξ_{\mathrm{cc}}$ of CAMIRA clusters, the auto-correlation function $ξ_{\mathrm{gg}}$ of the CMASS galaxies spectroscopically observed in the BOSS survey, and the cross-correlation function $ξ_{\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 $ξ_{\mathrm{cc}}$, $ξ_{\mathrm{cc}}+ξ_{\mathrm{cg}}$, and $ξ_{\mathrm{cc}} + ξ_{\mathrm{cg}} + ξ_{\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σ$. This implies that the absolute mass scale of CAMIRA clusters inferred from clustering is mildly higher than those from the independent methods. [abridged]

preprint2020arXiv

Active gas features in three HSC-SSP CAMIRA clusters revealed by high angular resolution analysis of MUSTANG-2 SZE and XXL X-ray observations

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 main 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 the integrated Compton $y$ parameter and the temperature for the major merger are 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.

preprint2020arXiv

Gas Density Perturbations in Cool Cores of CLASH Galaxy Clusters

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.

preprint2020arXiv

Redshift Evolution of Green Valley Galaxies in Different Environments from the Hyper Suprime-Cam Survey

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.

preprint2020arXiv

RELICS: A Very Large ($θ_{E}\sim40"$) Cluster Lens -- RXC J0032.1+1808

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 $θ_{E}\sim10\arcsec$ to $\sim30-35\arcsec$, with only a handful of clusters known with Einstein radii $θ_{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 $θ_{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.

preprint2020arXiv

The BUFFALO HST Survey

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.

preprint2020arXiv

The Dynamical State of the Frontier Fields Galaxy Cluster Abell 370

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. Moreover, our best model reproduces the observed velocity dispersion of cluster member galaxies, which supports the large total mass of A370 derived from weak lensing. 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.

preprint2020arXiv

The Radial Acceleration Relation in CLASH Galaxy Clusters

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 $Λ$CDM framework. Our results also predict the presence of a baryonic Faber--Jackson relation ($σ_v^4\propto M_\mathrm{bar}$) on cluster scales.

preprint2020arXiv

The Richness-to-Mass Relation of CAMIRA Galaxy Clusters from Weak-lensing Magnification in the Subaru Hyper Suprime-Cam Survey

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

preprint2020arXiv

Weak lensing Analysis of X-Ray-selected XXL Galaxy Groups and Clusters with Subaru HSC Data

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 β= -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 σ(\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σ$ 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.

preprint2019arXiv

XXL Survey groups and clusters in the Hyper Suprime-Cam Survey. Scaling relations between X-ray properties and weak lensing mass

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 ($β_{m_\text{g}|m}=1.73 \pm0.80$) and luminosity ($β_{l|m}=1.91\pm0.94$) and a nearly self-similar slope for the temperature ($β_{t|m}=0.78\pm0.43$). Intrinsic scatters of X-ray properties appear to be positively correlated at a fixed mass (median correlation factor $ρ_{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 $ρ_{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\%$).

preprint2018arXiv

The Reionization Lensing Cluster Survey (RELICS) and the Brightest High-z Galaxies

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.

preprint2016arXiv

AMiBA: Cluster Sunyaev-Zel'dovich Effect Observations with the Expanded 13-Element Array

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.

preprint2016arXiv

CLASH: Joint Analysis of Strong-Lensing, Weak-Lensing Shear and Magnification Data for 20 Galaxy Clusters

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 $σ(\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.

preprint2016arXiv

Probing WHIM around Galaxy Clusters with Fast Radio Bursts and the Sunyaev-Zel'dovich effect

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.

preprint2016arXiv

The Detailed Science Case for the Maunakea Spectroscopic Explorer: the Composition and Dynamics of the Faint Universe

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)

preprint2015arXiv

CLASH-X: A Comparison of Lensing and X-ray Techniques for Measuring the Mass Profiles of Galaxy Clusters

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.

preprint2015arXiv

Frontier Fields: Subaru Weak-Lensing Analysis of the Merging Galaxy Cluster A2744

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.

preprint2015arXiv

Hubble Space Telescope Combined Strong and Weak Lensing Analysis of the CLASH Sample: Mass and Magnification Models and Systematic Uncertainties

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, $κ$, and $\sim20\%$ in the magnification, $μ$. We show maps of these differences for each cluster, as well as the mass distributions, critical curves, and 2D integrated mass profiles. For the Einstein radii ($z_{s}=2$) we find that all typically agree within $10\%$ between the two models, and Einstein masses agree, typically, within $\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 (Σ)/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.

preprint2015arXiv

Radio Relics Tracing the Projected Mass Distribution in CIZA J2242.8+5301

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

preprint2015arXiv

The Detection and Statistics of Giant Arcs Behind CLASH Clusters

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.

preprint2015arXiv

Three-dimensional Multi-probe Analysis of the Galaxy Cluster A1689

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

preprint2014arXiv

CLASH-VLT: Constraints on the Dark Matter Equation of State from Accurate Measurements of Galaxy Cluster Mass Profiles

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ρ)=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 $ρ$ is the density. We plan to further improve our constraint by applying the same procedure to all clusters from the ongoing CLASH-VLT program.

preprint2014arXiv

CLASH: Extreme Emission Line Galaxies and Their Implication on Selection of High-Redshift Galaxies

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.

preprint2014arXiv

CLASH: Weak-Lensing Shear-and-Magnification Analysis of 20 Galaxy Clusters

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

preprint2014arXiv

Comparing gravitational redshifts of SDSS galaxy clusters with the magnification redshift enhancement of background BOSS galaxies

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.

preprint2014arXiv

Measurements of the Sunyaev-Zel'dovich Effect in MACS J0647.7+7015 and MACS J1206.2-0847 at High Angular Resolution with MUSTANG

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 $γ= 0.9$ fits the MUSTANG image with $χ^{2}_{red}=1.005$ and probability to exceed (PTE) = 0.34. For MACS J1206.2-0847, we find $γ=0.7$, $χ^{2}_{red}=0.993$, and PTE = 0.70. In addition, we find a significant ($>$3-$σ$) 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.

preprint2014arXiv

Three Gravitationally Lensed Supernovae Behind CLASH Galaxy Clusters

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.

preprint2013arXiv

A Measurement of the Kinetic Sunyaev-Zel'dovich Signal Toward MACS J0717.5+3745

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

preprint2013arXiv

CLASH: Complete Lensing Analysis of the Largest Cosmic Lens MACS J0717.5+3745 and Surrounding Structures

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σ_κ. 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σ departure above the maximum simulated halo mass at this redshift.

preprint2013arXiv

Cluster Lensing Profiles Derived from a Redshift Enhancement of Magnified BOSS-Survey Galaxies

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.

preprint2013arXiv

Galaxy halo truncation and Giant Arc Surface Brightness Reconstruction in the Cluster MACSJ1206.2-0847

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

preprint2013arXiv

Model-Free Multi-Probe Lensing Reconstruction of Cluster Mass Profiles

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.

preprint2013arXiv

Platform Deformation Phase Correction for the AMiBA-13 Co-planar Interferometer

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.

preprint2013arXiv

Sunyaev-Zel'dovich-Measured Pressure Profiles from the Bolocam X-ray/SZ Galaxy Cluster Sample

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.

preprint2013arXiv

Suzaku observations of Abell 1835 outskirts: Deviation from hydrostatic equilibrium

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.

preprint2013arXiv

Tangential Velocity of the Dark Matter in the Bullet Cluster from Precise Lensed Image Redshifts

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.

preprint2013arXiv

The Contribution of Halos with Different Mass Ratios to the Overall Growth of Cluster-Sized Halos

We provide a new observational test for a key prediction of the Λ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 $σ$ agreement with Λ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.

preprint2013arXiv

The Contribution of Radio Galaxy Contamination to Measurements of the Sunyaev-Zel'dovich Decrement in Massive Galaxy Clusters at 140 GHz with Bolocam

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.

preprint2012arXiv

A highly magnified candidate for a young galaxy seen when the Universe was 500 Myrs old

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

preprint2012arXiv

Bolocam Observations of Two Unconfirmed Galaxy Cluster Candidates from the Planck Early SZ Sample

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

preprint2012arXiv

CLASH: Mass Distribution in and around MACS J1206.2-0847 from a Full Cluster Lensing Analysis

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

preprint2012arXiv

CLASH: Precise New Constraints on the Mass Profile of Abell 2261

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 Λ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 Λ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%.

preprint2012arXiv

Cluster Lensing And Supernova survey with Hubble (CLASH): An Overview

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.

preprint2012arXiv

Cluster Mass Profiles from a Bayesian Analysis of Weak Lensing Distortion and Magnification Measurements: Applications to Subaru Data

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

preprint2012arXiv

Mass, shape and thermal properties of A1689 by a multi-wavelength X-ray, lensing and Sunyaev-Zel'dovich analysis

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.

preprint2012arXiv

Miscentring in Galaxy Clusters: Dark Matter to Brightest Cluster Galaxy Offsets in 10,000 SDSS Clusters

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, Δ, are distributed equally around zero with no preferred direction, and are well described by a log-normal distribution with <log_{10}(Δ[h^{-1} Mpc])>=-1.895^{+0.003}_{-0.004}, and σ=0.501\pm0.004 (95% confidence levels), or <log_{10}(Δ[\arcsec])>=0.564\pm0.005, and σ=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]

preprint2012arXiv

The Universal Einstein Radius Distribution from 10,000 SDSS Clusters

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(θ_{e}\arcsec)>=0.73^{+0.02}_{-0.03}$, $σ=0.316^{+0.004}_{-0.002}$, and with higher abundance of large $θ_{e}$ clusters than predicted by $Λ$CDM. We visually inspect each of the clusters with $θ_{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 $θ_{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.

preprint2012arXiv

The Weight of Emptiness: The Gravitational Lensing Signal of Stacked Voids

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.

preprint2011arXiv

A Precise Cluster Mass Profile Averaged from the Highest-Quality Lensing Data

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

preprint2011arXiv

A Weak Lensing Detection of the Cosmological Distance-Redshift Relation Behind Three Massive Clusters

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.

preprint2011arXiv

Cluster-Cluster Lensing and the Case of Abell 383

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

preprint2011arXiv

LoCuSS: Subaru Weak Lensing Study of 30 Galaxy Clusters

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

preprint2011arXiv

Star Formation and AGN activity in Galaxies classified using the 1.6 μm Bump and PAH features at $z = 0.4-2$

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.

preprint2011arXiv

Triaxiality and non-thermal gas pressure in Abell 1689

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.

preprint2011arXiv

Weak and strong lensing analysis of the triaxial matter distribution of A1689

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.

preprint2010arXiv

1.2 Meter Shielded Cassegrain Antenna for Close-Packed Radio Interferometer

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.

preprint2010arXiv

AMiBA Wideband Analog Correlator

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.

preprint2010arXiv

AMiBA: Sunyaev-Zel'dovich effect derived properties and scaling relations of massive galaxy clusters

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 β-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) β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.

preprint2010arXiv

Constraining Intra-cluster Gas Models with AMiBA13

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.

preprint2010arXiv

Contamination of the Central Sunyaev--Zel'dovich Decrements in AMiBA Galaxy Cluster Observations

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.

preprint2010arXiv

Detailed Cluster Mass and Light profiles of A1703, A370 and RXJ1347-11 from Deep Subaru Imaging

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, α~= -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.

preprint2010arXiv

Full Lensing Analysis of Abell 1703: Comparison of Independent Lens-Modelling Techniques

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 $Λ$CDM model, similar to other massive relaxed clusters with accurately determined lensing-based profiles.

preprint2010arXiv

Strong-Lensing Analysis of MS 1358.4+6245: New Multiple Images and Implications for the Well-Resolved z=4.92 Galaxy

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

preprint2010arXiv

Suzaku Observation of Abell 1689: Anisotropic Temperature and Entropy Distributions Associated with the Large-Scale Structure

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

preprint2010arXiv

The Mass Structure of the Galaxy Cluster Cl0024+1654 from a Full Lensing Analysis of Joint Subaru and ACS/NIC3 Observations

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.

preprint2009arXiv

AMiBA: Broadband Heterodyne CMB Interferometry

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.

preprint2009arXiv

AMiBA: scaling relations between the integrated Compton-y and X-ray derived temperature, mass, and luminosity

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.

preprint2009arXiv

AMiBA: System Performance

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.

preprint2009arXiv

Dynamical Study of A1689 from Wide-Field VLT/VIMOS Spectroscopy: Mass Profile, Concentration Parameter, and Velocity Anisotropy

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

preprint2009arXiv

Mass and Hot Baryons in Massive Galaxy Clusters from Subaru Weak Lensing and AMiBA SZE Observations

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 = Ω_b/Ω_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.

preprint2009arXiv

New Multiply-Lensed Galaxies Identified in ACS/NIC3 Observations of Cl0024+1654 Using an Improved Mass Model

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

preprint2009arXiv

The AMiBA Hexapod Telescope Mount

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.

preprint2009arXiv

The Yuan-Tseh Lee Array for Microwave Background Anisotropy

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.

preprint2008arXiv

Comparison of Cluster Lensing Profiles with Lambda CDM Predictions

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.

preprint2008arXiv

Subaru Weak Lensing Study of Seven Merging Clusters: Distributions of Mass and Baryons

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.

preprint2008arXiv

Tests of AMiBA Data Integrity

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.

preprint2007arXiv

A New Measure for Weak Lensing Flexion

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.

preprint2007arXiv

Observational Constraints on the ICM Temperature Enhancement by Cluster Mergers

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.

preprint2007arXiv

Using Weak Lensing Dilution to Improve Measurements of the Luminous and Dark Matter in A1689

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.

preprint2005arXiv

Discovery of a Strongly Lensed Galaxy at z=3.9 behind a z=0.83 Galaxy Cluster

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.

preprint2005arXiv

Simulation of a Combined SZE and Weak Lensing Cluster Survey for AMiBA Experiment

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.

preprint2004arXiv

The Surprisingly Steep Mass Profile of Abell 1689, from a Lensing Analysis of Subaru Images

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.

preprint2002arXiv

The topology of weak lensing field in the neighborhood of MS1054-03

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.

preprint2000arXiv

Cluster Mass Reconstruction by a Weak Shear Field

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.