CIZA J2242.8+5301 is a merging system with a prominent (~2 Mpc long) radio relic, which together with the morphology of the X-ray emission provides strong evidence for a violent collision along the N-S axis. We present our constraints on the dark mat
ter distribution of this unusual system using Subaru and CFHT imaging data. Measuring a high S/N lensing signal from this cluster is potentially a challenging task because of its proximity to the Milky Way plane (|b|~5 deg). We overcome this challenge with careful observation planning and systematics control, which enables us to successfully map the dark matter distribution of the cluster with high fidelity. The resulting mass map shows that the mass distribution is highly elongated along the N-S merger axis inferred from the orientation of the radio relics. Based on our mass reconstruction, we identify two sub-clusters, which coincide with the cluster galaxy distributions. We determine their masses using MCMC analysis by simultaneously fitting two NFW halos without fixing their centroids. The resulting masses of the northern and southern systems are $M_{200}=11.0_{-3.2}^{+3.7}times10^{14} M_{sun}$ and $9.8_{-2.5}^{+3.8}times10^{14} M_{sun}$, respectively, indicating that we are witnessing a post-collision of two giant systems of nearly equal mass. When the mass and galaxy centroids are compared in detail, we detect ~ 1 (~190 kpc) offsets in both northern and southern sub-clusters. We find that the galaxy luminosity-mass offset for the northern clump is statistically significant at the ~2 sigma level whereas the detection is only marginal for the southern sub-cluster in part because of a relatively large mass centroid error. We conclude that it is yet premature to uniquely attribute the galaxy-mass misalignment to SIDM and discuss caveats.
We present the results from a weak gravitational lensing study of the merging cluster A520 based on the analysis of Hubble Space Telescope/Advanced Camera for Surveys (ACS) data. The excellent data quality allows us to reach a mean number density of
source galaxies of ~109 per sq. arcmin, which improves both resolution and significance of the mass reconstruction compared to a previous study based on Wide Field Planetary Camera 2 (WFPC2) images. We take care in removing instrumental effects such as the trailing of charge due to radiation damage of the ACS detector and the position-dependent point spread function (PSF). This new ACS analysis confirms the previous claims that a substantial amount of dark mass is present between two luminous subclusters. We examine the distribution of cluster galaxies and observe very little light at this location. We find that the centroid of the dark peak in the current ACS analysis is offset to the southwest by ~1 arcmin with respect to the centroid from the WFPC2 analysis. Interestingly, this new centroid is in better spatial agreement with the location where the X-ray emission is strongest, and the mass-to-light ratio estimated with this centroid is much higher 813+-78 M_sun/L_Rsun than the previous value; the aperture mass based on the WFPC2 centroid provides a slightly lower, but consistent mass. Although we cannot provide a definite explanation for the presence of the dark peak, we discuss a revised scenario, wherein dark matter with a more conventional range sigma_DM/m_DM < 1 cm^2/g of self-interacting cross-section can lead to the detection of this dark substructure. If supported by detailed numerical simulations, this hypothesis opens up the possibility that the A520 system can be used to establish a lower limit of the self-interacting cross-section of dark matter.
(Abridged) We present a HST weak-lensing study of the merging galaxy cluster El Gordo (ACT-CL J0102-4915) at z=0.87 discovered by the Atacama Cosmology Telescope collaboration as the strongest SZ decrement in its ~1000 sq. deg survey. Our weak-lensin
g analysis confirms that ACT-CL J0102-4915 is indeed an extreme system consisting of two massive (~10^15 Msun each) subclusters with a projected separation of ~0.7 Mpc. This binary mass structure revealed by our lensing study is consistent with the cluster galaxy distribution and the dynamical study carried out with 89 spectroscopic members. We estimate the mass of ACT-CL J0102-4915 by simultaneously fitting two axisymmetric NFW profiles allowing their centers to vary. Our MCMC analysis shows that the masses of the northwestern (NW) and the southeastern (SE) components are M200c=(1.38+-0.22) x 10^15 Msun and (0.78+-0.20) x 10^15 Msun, respectively. The lensing-based velocity dispersions are consistent with their spectroscopic measurements. The centroids of both components are tightly constrained (~4) and close to the optical luminosity centers. The X-ray and mass peaks are spatially offset by ~8 (~62 kpc), which is significant at the ~2 sigma confidence level and confirms that the baryonic and dark matter in this cluster are disassociated. The dark matter peak, however, does not lead the gas peak in the direction expected if we are viewing the cluster soon after first core passage during a high speed merger. Under the assumption that the merger is happening in the plane of the sky, extrapolation of the two NFW halos to a radius r200a=2.4 Mpc yields a combined mass of M200a=(3.13+-0.56) x 10^15 Msun. This extrapolated total mass is consistent with our two-component-based dynamical analysis and previous X-ray measurements, projecting ACT-CL J0102-4915 to be the most massive cluster at z>0.6 known to date.
We present a cosmic shear study from the Deep Lens Survey (DLS), a deep BVRz multi-band imaging survey of five 4 sq. degree fields with two National Optical Astronomy Observatory (NOAO) 4-meter telescopes at Kitt Peak and Cerro Tololo. For both teles
copes, the change of the point-spread-function (PSF) shape across the focal plane is complicated, and the exposure-to-exposure variation of this position-dependent PSF change is significant. We overcome this challenge by modeling the PSF separately for individual exposures and CCDs with principal component analysis (PCA). We find that stacking these PSFs reproduces the final PSF pattern on the mosaic image with high fidelity, and the method successfully separates PSF-induced systematics from gravitational lensing effects. We calibrate our shears and estimate the errors, utilizing an image simulator, which generates sheared ground-based galaxy images from deep Hubble Space Telescope archival data with a realistic atmospheric turbulence model. For cosmological parameter constraints, we marginalize over shear calibration error, photometric redshift uncertainty, and the Hubble constant. We use cosmology-dependent covariances for the Markov Chain Monte Carlo analysis and find that the role of this varying covariance is critical in our parameter estimation. Our current non-tomographic analysis alone constrains the Omega_M-sigma_8 likelihood contour tightly, providing a joint constraint of Omega_M=0.262+-0.051 and sigma_8=0.868+-0.071. We expect that a future DLS weak-lensing tomographic study will further tighten these constraints because explicit treatment of the redshift dependence of cosmic shear more efficiently breaks the Omega_M-sigma_8 degeneracy. Combining the current results with the Wilkinson Microwave Anisotropy Probe 7-year (WMAP7) likelihood data, we obtain Omega_M=0.278+-0.018 and sigma_8=0.815+-0.020.
We present a Hubble Space Telescope/Wide Field Planetary Camera 2 weak-lensing study of A520, where a previous analysis of ground-based data suggested the presence of a dark mass concentration. We map the complex mass structure in much greater detail
leveraging more than a factor of three increase in the number density of source galaxies available for lensing analysis. The dark core that is coincident with the X-ray gas peak, but not with any stellar luminosity peak is now detected with more than 10 sigma significance. The ~1.5 Mpc filamentary structure elongated in the NE-SW direction is also clearly visible. Taken at face value, the comparison among the centroids of dark matter, intracluster medium, and galaxy luminosity is at odds with what has been observed in other merging clusters with a similar geometric configuration. To date, the most remarkable counter-example might be the Bullet Cluster, which shows a distinct bow-shock feature as in A520, but no significant weak-lensing mass concentration around the X-ray gas. With the most up-to-date data, we consider several possible explanations that might lead to the detection of this peculiar feature in A520. However, we conclude that none of these scenarios can be singled out yet as the definite explanation for this puzzle.
We present weak gravitational lensing analysis of 22 high-redshift (z >~1) clusters based on Hubble Space Telescope images. Most clusters in our sample provide significant lensing signals and are well detected in their reconstructed two-dimensional m
ass maps. Combining the current results and our previous weak-lensing studies of five other high-z clusters, we compare gravitational lensing masses of these clusters with other observables. We revisit the question whether the presence of the most massive clusters in our sample is in tension with the current LambdaCDM structure formation paradigm. We find that the lensing masses are tightly correlated with the gas temperatures and establish, for the first time, the lensing mass-temperature relation at z >~ 1. For the power law slope of the M-TX relation (M propto T^{alpha}), we obtain alpha=1.54 +/- 0.23. This is consistent with the theoretical self-similar prediction alpha=3/2 and with the results previously reported in the literature for much lower redshift samples. However, our normalization is lower than the previous results by 20-30%, indicating that the normalization in the M-TX relation might evolve. After correcting for Eddington bias and updating the discovery area with a more conservative choice, we find that the existence of the most massive clusters in our sample still provides a tension with the current Lambda CDM model. The combined probability of finding the four most massive clusters in this sample after marginalization over current cosmological parameters is less than 1%.
The weak-lensing science of the LSST project drives the need to carefully model and separate the instrumental artifacts from the intrinsic lensing signal. The dominant source of the systematics for all ground based telescopes is the spatial correlati
on of the PSF modulated by both atmospheric turbulence and optical aberrations. In this paper, we present a full FOV simulation of the LSST images by modeling both the atmosphere and the telescope optics with the most current data for the telescope specifications and the environment. To simulate the effects of atmospheric turbulence, we generated six-layer phase screens with the parameters estimated from the on-site measurements. For the optics, we combined the ray-tracing tool ZEMAX and our simulated focal plane data to introduce realistic aberrations and focal plane height fluctuations. Although this expected flatness deviation for LSST is small compared with that of other existing cameras, the fast f-ratio of the LSST optics makes this focal plane flatness variation and the resulting PSF discontinuities across the CCD boundaries significant challenges in our removal of the systematics. We resolve this complication by performing PCA CCD-by-CCD, and interpolating the basis functions using conventional polynomials. We demonstrate that this PSF correction scheme reduces the residual PSF ellipticity correlation below 10^-7 over the cosmologically interesting scale. From a null test using HST/UDF galaxy images without input shear, we verify that the amplitude of the galaxy ellipticity correlation function, after the PSF correction, is consistent with the shot noise set by the finite number of objects. Therefore, we conclude that the current optical design and specification for the accuracy in the focal plane assembly are sufficient to enable the control of the PSF systematics required for weak-lensing science with the LSST.
We present a weak-lensing analysis of the z=1.4 galaxy cluster XMMU J2235.3-2557, based on deep Advanced Camera for Surveys images. Despite the observational challenge set by the high redshift of the lens, we detect a substantial lensing signal at th
e >~ 8 sigma level. This clear detection is enabled in part by the high mass of the cluster, which is verified by our both parametric and non-parametric estimation of the cluster mass. Assuming that the cluster follows a Navarro-Frenk-White mass profile, we estimate that the projected mass of the cluster within r=1 Mpc is (8.5+-1.7) x 10^14 solar mass, where the error bar includes the statistical uncertainty of the shear profile, the effect of possible interloping background structures, the scatter in concentration parameter, and the error in our estimation of the mean redshift of the background galaxies. The high X-ray temperature 8.6_{-1.2}^{+1.3} keV of the cluster recently measured with Chandra is consistent with this high lensing mass. When we adopt the 1-sigma lower limit as a mass threshold and use the cosmological parameters favored by the Wilkinson Microwave Anisotropy Probe 5-year (WMAP5) result, the expected number of similarly massive clusters at z >~ 1.4 in the 11 square degree survey is N ~ 0.005. Therefore, the discovery of the cluster within the survey volume is a rare event with a probability < 1%, and may open new scenarios in our current understanding of cluster formation within the standard cosmological model.
We present a weak-lensing analysis of the galaxy cluster CL J1226+3332 at z=0.89 using Hubble Space Telescope Advanced Camera for Surveys images. The cluster is the hottest (>10 keV), most X-ray luminous system at z>0.6 known to date. The relaxed X-r
ay morphology, as well as its high temperature, is unusual at such a high redshift. Our mass reconstruction shows that on a large scale the dark matter distribution is consistent with a relaxed system with no significant substructures. However, on a small scale the cluster core is resolved into two mass clumps highly correlated with the cluster galaxy distribution. The dominant mass clump lies close to the brightest cluster galaxy whereas the other less massive clump is located ~40 (~310 kpc) to the southwest. Although this secondary mass clump does not show an excess in the X-ray surface brightness, the gas temperature of the region is much higher (12~18 keV) than those of the rest. We propose a scenario in which the less massive system has already passed through the main cluster and the X-ray gas has been stripped during this passage. The elongation of the X-ray peak toward the southwestern mass clump is also supportive of this possibility. We measure significant tangential shears out to the field boundary (~1.5 Mpc), which are well described by an Navarro-Frenk-White profile with a concentration parameter of c200=2.7+-0.3 and a scale length of rs=78+-19 (~600 kpc) with chi^2/d.o.f=1.11. Within the spherical volume r200=1.6 Mpc, the total mass of the cluster becomes M(r<r200)=(1.4+-0.2) x 10^15 solar mass. Our weak-lensing analysis confirms that CL1226+3332 is indeed the most massive cluster known to date at z>0.6.
We describe the time- and position-dependent point spread function (PSF) variation of the Wide Field Channel (WFC) of the Advanced Camera for Surveys (ACS) with the principal component analysis (PCA) technique. The time-dependent change is caused by
the temporal variation of the $HST$ focus whereas the position-dependent PSF variation in ACS/WFC at a given focus is mainly the result of changes in aberrations and charge diffusion across the detector, which appear as position-dependent changes in elongation of the astigmatic core and blurring of the PSF, respectively. Using >400 archival images of star cluster fields, we construct a ACS PSF library covering diverse environments of the $HST$ observations (e.g., focus values). We find that interpolation of a small number ($sim20$) of principal components or ``eigen-PSFs per exposure can robustly reproduce the observed variation of the ellipticity and size of the PSF. Our primary interest in this investigation is the application of this PSF library to precision weak-lensing analyses, where accurate knowledge of the instruments PSF is crucial. However, the high-fidelity of the model judged from the nice agreement with observed PSFs suggests that the model is potentially also useful in other applications such as crowded field stellar photometry, galaxy profile fitting, AGN studies, etc., which similarly demand a fair knowledge of the PSFs at objects locations. Our PSF models, applicable to any WFC image rectified with the Lanczos3 kernel, are publicly available.
