Source author record

Eduardo Rozo

Eduardo Rozo 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

49works
6topics
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

49 published item(s)

preprint2022arXiv

KaRMMa -- Kappa Reconstruction for Mass Mapping

We present KaRMMa, a novel method for performing mass map reconstruction from weak-lensing surveys. We employ a fully Bayesian approach with a physically motivated lognormal prior to sample from the posterior distribution of convergence maps. We test KaRMMa on a suite of dark matter N-body simulations with simulated DES Y1-like shear observations. We show that KaRMMa outperforms the basic Kaiser-Squires mass map reconstruction in two key ways: 1) our best map point estimate has lower residuals compared to Kaiser-Squires; and 2) unlike the Kaiser-Squires reconstruction, the posterior distribution of KaRMMa maps are nearly unbiased in all summary statistics we considered, namely: one-point and two-point functions, and peak/void counts. In particular, KaRMMa successfully captures the non-Gaussian nature of the distribution of $κ$ values in the simulated maps. We further demonstrate that the KaRMMa posteriors correctly characterize the uncertainty in all summary statistics we considered.

preprint2022arXiv

Map-based cosmology inference with lognormal cosmic shear maps

Most cosmic shear analyses to date have relied on summary statistics (e.g. $ξ_+$ and $ξ_-$). These types of analyses are necessarily sub-optimal, as the use of summary statistics is lossy. In this paper, we forward-model the convergence field of the Universe as a lognormal random field conditioned on the observed shear data. This new map-based inference framework enables us to recover the joint posterior of the cosmological parameters and the convergence field of the Universe. Our analysis properly accounts for the covariance in the mass maps across tomographic bins, which significantly improves the fidelity of the maps relative to single-bin reconstructions. We verify that applying our inference pipeline to Gaussian random fields recovers posteriors that are in excellent agreement with their analytical counterparts. At the resolution of our maps -- and to the extent that the convergence field can be described by the lognormal model -- our map posteriors allow us to reconstruct \it all \rm summary statistics (including non-Gaussian statistics). We forecast that a map-based inference analysis of LSST-Y10 data can improve cosmological constraints in the $σ_8$--$Ω_{\rm m}$ plane by $\approx 30\%$ relative to the currently standard cosmic shear analysis. This improvement happens almost entirely along the $S_8=σ_8Ω_{\rm m}^{1/2}$ directions, meaning map-based inference fails to significantly improve constraints on $S_8$.

preprint2022arXiv

Optical selection bias and projection effects in stacked galaxy cluster weak lensing

Cosmological constraints from current and upcoming galaxy cluster surveys are limited by the accuracy of cluster mass calibration. In particular, optically identified galaxy clusters are prone to selection effects that can bias the weak lensing mass calibration. We investigate the selection bias of the stacked cluster lensing signal associated with optically selected clusters, using clusters identified by the redMaPPer algorithm in the Buzzard simulations as a case study. We find that at a given cluster halo mass, the residuals of redMaPPer richness and weak lensing signal are positively correlated. As a result, for a given richness selection, the stacked lensing signal is biased high compared with what we would expect from the underlying halo mass probability distribution. The cluster lensing selection bias can thus lead to overestimated mean cluster mass and biased cosmology results. We show that the lensing selection bias exhibits a strong scale-dependence and is approximately 20 to 60 percent for $ΔΣ$ at large scales. This selection bias largely originates from spurious member galaxies within +/- 20 to 60 Mpc/h along the line of sight, highlighting the importance of quantifying projection effects associated with the broad redshift distribution of member galaxies in photometric cluster surveys. While our results qualitatively agree with those in the literature, accurate quantitative modelling of the selection bias is needed to achieve the goals of cluster lensing cosmology and will require synthetic catalogues covering a wide range of galaxy-halo connection models.

preprint2021arXiv

The Phase Space Structure of Dark Matter Halos

The phase space structure of dark matter halos can be used to measure the mass of the halo, infer mass accretion rates, and probe the effects of modified gravity. Previous studies showed that the splashback radius can be measured in position space using the slope of the density profile. Using N-body simulations, we show that the phase space structure of the dark matter halo does not end at this splashback radius. Instead, there exists a region where infalling, splashback, and virialized halos are mixed spatially. We model the distribution of the three kinematically distinct populations and show that there exists an "edge radius" beyond which a dark matter halo has no orbiting substructures. This radius is a fixed multiple of the splashback radius as defined in previous works, and can be interpreted as a radius which contains a fixed fraction of the apocenters of dark matter particles. Our results provide a firm theoretical foundation to the satellite galaxy model adopted in the companion paper by Tomooka et al., where we analyzed the phase space distribution of SDSS redMaPPer clusters.

preprint2020arXiv

Dark Energy Survey Year 1 Results: Cosmological Constraints from Cluster Abundances and Weak Lensing

We perform a joint analysis of the counts and weak lensing signal of redMaPPer clusters selected from the Dark Energy Survey (DES) Year 1 dataset. Our analysis uses the same shear and source photometric redshifts estimates as were used in the DES combined probes analysis. Our analysis results in surprisingly low values for $S_8 =σ_8(Ω_{\rm m}/0.3)^{0.5}= 0.65\pm 0.04$, driven by a low matter density parameter, $Ω_{\rm m}=0.179^{+0.031}_{-0.038}$, with $σ_8-Ω_{\rm m}$ posteriors in $2.4σ$ tension with the DES Y1 3x2pt results, and in $5.6σ$ with the Planck CMB analysis. These results include the impact of post-unblinding changes to the analysis, which did not improve the level of consistency with other data sets compared to the results obtained at the unblinding. The fact that multiple cosmological probes (supernovae, baryon acoustic oscillations, cosmic shear, galaxy clustering and CMB anisotropies), and other galaxy cluster analyses all favor significantly higher matter densities suggests the presence of systematic errors in the data or an incomplete modeling of the relevant physics. Cross checks with X-ray and microwave data, as well as independent constraints on the observable--mass relation from SZ selected clusters, suggest that the discrepancy resides in our modeling of the weak lensing signal rather than the cluster abundance. Repeating our analysis using a higher richness threshold ($λ\ge 30$) significantly reduces the tension with other probes, and points to one or more richness-dependent effects not captured by our model.

preprint2020arXiv

Localizing Transformations of the Galaxy-Galaxy Lensing Observable

Modern cosmological analyses of galaxy-galaxy lensing face a theoretical systematic effect arising from the non-locality of the observed galaxy-galaxy lensing signal. Because the predicted tangential shear signal at a given separation depends on the physical modeling on all scales internal to that separation, systematic uncertainties in the modeling of non-linear small scales are propagated outwards to larger scales. Even in the absence of other limiting factors, this systematic effect alone can necessitate conservative small-scale cuts, resulting in significant losses of information in the tangential shear data vector. We construct a simple linear transformation of the observable that removes this non-locality, enabling more aggressive small-scale cuts for a given theoretical model. Our modified galaxy-galaxy lensing observable makes it possible to include observations on significantly smaller scales than those under the standard approach in cosmological analyses. More importantly, it ensures that the cosmological signal contained within the observable is exclusively drawn from well-understood physical scales.

preprint2019arXiv

Concordance Cosmology?

We propose a new intuitive metric for evaluating the tension between two experiments, and apply it to several data sets. While our metric is non-optimal, if evidence of tension is detected, this evidence is robust and easy to interpret. Assuming a flat $Λ$CDM cosmological model, we find that there is a modest $2.2σ$ tension between the DES Year 1 results and the ${\it Planck}$ measurements of the Cosmic Microwave Background (CMB). This tension is driven by the difference between the amount of structure observed in the late-time Universe and that predicted from fitting the ${\it Planck}$ data, and appears to be unrelated to the tension between ${\it Planck}$ and local esitmates of the Hubble rate. In particular, combining DES, Baryon Acoustic Oscillations (BAO), Big-Bang Nucleosynthesis (BBN), and supernovae (SNe) measurements recovers a Hubble constant and sound horizon consistent with ${\it Planck}$, and in tension with local distance-ladder measurements. If the tension between these various data sets persists, it is likely that reconciling ${\it all}$ current data will require breaking the flat $Λ$CDM model in at least two different ways: one involving new physics in the early Universe, and one involving new late-time Universe physics.

preprint2016arXiv

Constraining the Mass-Richness Relationship of redMaPPer Clusters with Angular Clustering

The potential of using cluster clustering for calibrating the mass-observable relation of galaxy clusters has been recognized theoretically for over a decade. Here, we demonstrate the feasibility of this technique to achieve high precision mass calibration using redMaPPer clusters in the Sloan Digital Sky Survey North Galactic Cap. By including cross-correlations between several richness bins in our analysis we significantly improve the statistical precision of our mass constraints. The amplitude of the mass-richness relation is constrained to 7% statistical precision. However, the error budget is systematics dominated, reaching an 18% total error that is dominated by theoretical uncertainty in the bias-mass relation for dark matter halos. We perform a detailed treatment of the effects of assembly bias on our analysis, finding that the contribution of such effects to our parameter uncertainties is somewhat greater than that of measurement noise. We confirm the results from Miyatake et al. (2015) that the clustering amplitude of redMaPPer clusters depends on galaxy concentration, and provide additional evidence in support of this effect being due to some form of assembly bias. The results presented here demonstrate the power of cluster clustering for mass calibration and cosmology provided the current theoretical systematics can be ameliorated.

preprint2016arXiv

Cosmology from large scale galaxy clustering and galaxy-galaxy lensing with Dark Energy Survey Science Verification data

We present cosmological constraints from the Dark Energy Survey (DES) using a combined analysis of angular clustering of red galaxies and their cross-correlation with weak gravitational lensing of background galaxies. We use a 139 square degree contiguous patch of DES data from the Science Verification (SV) period of observations. Using large scale measurements, we constrain the matter density of the Universe as Omega_m = 0.31 +/- 0.09 and the clustering amplitude of the matter power spectrum as sigma_8 = 0.74 +/- 0.13 after marginalizing over seven nuisance parameters and three additional cosmological parameters. This translates into S_8 = sigma_8(Omega_m/0.3)^{0.16} = 0.74 +/- 0.12 for our fiducial lens redshift bin at 0.35 <z< 0.5, while S_8 = 0.78 +/- 0.09 using two bins over the range 0.2 <z< 0.5. We study the robustness of the results under changes in the data vectors, modelling and systematics treatment, including photometric redshift and shear calibration uncertainties, and find consistency in the derived cosmological parameters. We show that our results are consistent with previous cosmological analyses from DES and other data sets and conclude with a joint analysis of DES angular clustering and galaxy-galaxy lensing with Planck CMB data, Baryon Accoustic Oscillations and Supernova type Ia measurements.

preprint2016arXiv

Detection of the Splashback Radius and Halo Assembly bias of Massive Galaxy Clusters

We show that the projected number density profiles of SDSS photometric galaxies around galaxy clusters displays strong evidence for the splashback radius, a sharp halo edge corresponding to the location of the first orbital apocenter of satellite galaxies after their infall. We split the clusters into two subsamples with different mean projected radial distances of their members, $\langle R_{\rm mem}\rangle$, at fixed richness and redshift, and show that the sample with smaller $\langle R_{\rm mem}\rangle$ has a smaller ratio of the splashback radius to the traditional halo boundary $R_{\rm 200m}$, than the subsample with larger $\langle R_{\rm mem}\rangle$, indicative of different mass accretion rates for the two subsamples. The same cluster samples were recently used by Miyatake et al. to show that their large-scale clustering differs despite their similar weak lensing masses, demonstrating strong evidence for halo assembly bias. We expand on this result by presenting a 6.6-$σ$ detection of halo assembly bias using the cluster-photometric galaxy cross-correlations. Our measured splashback radii are smaller, while the strength of the assembly bias signal is stronger, than expectations from N-body simulations based on the $Λ$-dominated, cold dark matter structure formation model. Dynamical friction or cluster-finding systematics such as miscentering or projection effects are not likely to be the sole source of these discrepancies.

preprint2016arXiv

Evidence of Halo Assembly Bias in Massive Clusters

We present significant evidence of halo assembly bias for SDSS redMaPPer galaxy clusters in the redshift range $[0.1, 0.33]$. By dividing the 8,648 clusters into two subsamples based on the average member galaxy separation from the cluster center, we first show that the two subsamples have very similar halo mass of $M_{\rm 200m}\simeq 1.9\times 10^{14}~h^{-1}M_\odot$ based on the weak lensing signals at small radii $R<\sim 10~h^{-1}{\rm Mpc}$. However, their halo bias inferred from both the large-scale weak lensing and the projected auto-correlation functions differs by a factor of $\sim$1.5, which is a signature of assembly bias. The same bias hypothesis for the two subsamples is excluded at 2.5$σ$ in the weak lensing and 4.4$σ$ in the auto-correlation data, respectively. This result could bring a significant impact on both galaxy evolution and precision cosmology.

preprint2016arXiv

Galaxy Cluster Mass Estimation from Stacked Spectroscopic Analysis

We use simulated galaxy surveys to study: i) how galaxy membership in redMaPPer clusters maps to the underlying halo population, and ii) the accuracy of a mean dynamical cluster mass, $M_σ(λ)$, derived from stacked pairwise spectroscopy of clusters with richness $λ$. Using $\sim\! 130,000$ galaxy pairs patterned after the SDSS redMaPPer cluster sample study of Rozo et al. (2015 RMIV), we show that the pairwise velocity PDF of central--satellite pairs with $m_i < 19$ in the simulation matches the form seen in RMIV. Through joint membership matching, we deconstruct the main Gaussian velocity component into its halo contributions, finding that the top-ranked halo contributes $\sim 60\%$ of the stacked signal. The halo mass scale inferred by applying the virial scaling of Evrard et al. (2008) to the velocity normalization matches, to within a few percent, the log-mean halo mass derived through galaxy membership matching. We apply this approach, along with mis-centering and galaxy velocity bias corrections, to estimate the log-mean matched halo mass at $z=0.2$ of SDSS redMaPPer clusters. Employing the velocity bias constraints of Guo et al. (2015), we find $\langle \ln(M_{200c})|λ\rangle = \ln(M_{30}) + α_m \ln(λ/30)$ with $M_{30} = 1.56 \pm 0.35 \times 10^{14} M_\odot$ and $α_m = 1.31 \pm 0.06_{stat} \pm 0.13_{sys}$. Systematic uncertainty in the velocity bias of satellite galaxies overwhelmingly dominates the error budget.

preprint2016arXiv

Intrinsic alignments in redMaPPer clusters -- I. Central galaxy alignments and angular segregation of satellites

The shapes of cluster central galaxies are not randomly oriented, but rather exhibit coherent alignments with the shapes of their parent clusters as well as with the surrounding large-scale structures. In this work, we aim to identify the galaxy and cluster quantities that most strongly predict the central galaxy alignment phenomenon among a large parameter space with a sample of 8237 clusters and 94817 members within 0.1<z<0.35, based on the redMaPPer cluster catalog constructed from the Sloan Digital Sky Survey. We first quantify the alignment between the projected central galaxy shapes and the distribution of member satellites, to understand what central galaxy and cluster properties most strongly correlate with these alignments. Next, we investigate the angular segregation of satellites with respect to their central galaxy major axis directions, to identify the satellite properties that most strongly predict their angular segregation. We find that central galaxies are more aligned with their member galaxy distributions in clusters that are more elongated and have higher richness, and for central galaxies with larger physical size, higher luminosity and centering probability, and redder color. Satellites with redder color, higher luminosity, located closer to the central galaxy, and with smaller ellipticity show a stronger angular segregation toward their central galaxy major axes. Finally, we provide physical explanations for some of the identified correlations, and discuss the connection to theories of central galaxy alignments, the impact of primordial alignments with tidal fields, and the importance of anisotropic accretion.

preprint2016arXiv

Measuring subhalo mass in redMaPPer clusters with CFHT Stripe 82 Survey

We use the shear catalog from the CFHT Stripe-82 Survey to measure the subhalo masses of satellite galaxies in redMaPPer clusters. Assuming a Chabrier Initial Mass Function (IMF) and a truncated NFW model for the subhalo mass distribution, we find that the sub-halo mass to galaxy stellar mass ratio increases as a function of projected halo-centric radius $r_p$, from $M_{\rm sub}/M_{\rm star}=4.43^{+ 6.63}_{- 2.23}$ at $r_p \in [0.1,0.3]$ $h^{-1}Mpc$ to $M_{\rm sub}/M_{\rm star}=75.40^{+ 19.73}_{- 19.09}$ at $r_p \in [0.6,0.9]$ $h^{-1}Mpc$. We also investigate the dependence of subhalo masses on stellar mass by splitting satellite galaxies into two stellar mass bins: $10<\log(M_{\rm star}/M_{\rm sun})<10.5$ and $11<\log(M_{\rm star}/M_{\rm sun})<12$. The best-fit subhalo mass of the more massive satellite galaxy bin is larger than that of the less massive satellites: $\log(M_{\rm sub}/M_{\rm sun})=11.14 ^{+ 0.66 }_{- 0.73}$ ($M_{\rm sub}/M_{\rm star}=19.5^{+19.8}_{-17.9}$) versus $\log(M_{\rm sub}/M_{\rm sun})=12.38 ^{+ 0.16 }_{- 0.16}$ ($M_{\rm sub}/M_{\rm star}=21.1^{+7.4}_{-7.7}$).

preprint2016arXiv

SPIDERS: the spectroscopic follow-up of X-ray selected clusters of galaxies in SDSS-IV

SPIDERS (The SPectroscopic IDentification of eROSITA Sources) is a program dedicated to the homogeneous and complete spectroscopic follow-up of X-ray AGN and galaxy clusters over a large area ($\sim$7500 deg$^2$) of the extragalactic sky. SPIDERS is part of the SDSS-IV project, together with the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) and the Time-Domain Spectroscopic Survey (TDSS). This paper describes the largest project within SPIDERS before the launch of eROSITA: an optical spectroscopic survey of X-ray selected, massive ($\sim 10^{14}$ to $10^{15}~M_{\odot}$) galaxy clusters discovered in ROSAT and XMM-Newton imaging. The immediate aim is to determine precise ($Δ_z \sim 0.001$) redshifts for 4,000-5,000 of these systems out to $z \sim 0.6$. The scientific goal of the program is precision cosmology, using clusters as probes of large-scale structure in the expanding Universe. We present the cluster samples, target selection algorithms and observation strategies. We demonstrate the efficiency of selecting targets using a combination of SDSS imaging data, a robust red-sequence finder and a dedicated prioritization scheme. We describe a set of algorithms and work-flow developed to collate spectra and assign cluster membership, and to deliver catalogues of spectroscopically confirmed clusters. We discuss the relevance of line-of-sight velocity dispersion estimators for the richer systems. We illustrate our techniques by constructing a catalogue of 230 spectroscopically validated clusters ($0.031 < z < 0.658$), found in pilot observations. We discuss two potential science applications of the SPIDERS sample: the study of the X-ray luminosity-velocity dispersion ($L_X-σ$) relation and the building of stacked phase-space diagrams.

preprint2016arXiv

The DESI Experiment Part I: Science,Targeting, and Survey Design

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$α$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

preprint2016arXiv

The DESI Experiment Part II: Instrument Design

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= λ/Δλ$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

preprint2016arXiv

Tidal stripping as a test of satellite quenching in redMaPPer clusters

When dark matter halos are accreted by massive host clusters, strong gravitational tidal forces begin stripping mass from the accreted subhalos. This stripping eventually removes all mass beyond a subhalo's tidal radius, but the unbound mass remains in the vicinity of the satellite for at least a dynamical time t_dyn. The N-body subhalo study of Chamberlain et al. verified this picture and pointed out a useful observational consequence: measurements of subhalo correlations beyond the tidal radius are sensitive to the infall time, t_infall, of the subhalo onto its host. We perform this cross-correlation measurement using ~ 160,000 red satellite galaxies in SDSS redMaPPer clusters and find evidence that subhalo correlations do persist well beyond the tidal radius, suggesting that many of the observed satellites fell into their current host less than a dynamical time ago, t_infall < t_dyn. Combined with estimated dynamical times t_dyn ~ 3-5 Gyr and SED fitting results for the time at which satellites stopped forming stars, t_quench ~ 6 Gyr, we infer that for a significant fraction of the satellites, star formation quenched before those satellites entered their current hosts. The result holds for red satellites over a large range of cluster-centric distances 0.1 - 0.6 Mpc/h. We discuss the implications of this result for models of galaxy formation.

preprint2015arXiv

Environment-based selection effects of Planck clusters

We investigate whether the large scale structure environment of galaxy clusters imprints a selection bias on Sunyaev Zel'dovich (SZ) catalogs. Such a selection effect might be caused by line of sight (LoS) structures that add to the SZ signal or contain point sources that disturb the signal extraction in the SZ survey. We use the Planck PSZ1 union catalog (Planck Collab- oration et al. 2013a) in the SDSS region as our sample of SZ selected clusters. We calculate the angular two-point correlation function (2pcf) for physically correlated, foreground and background structure in the RedMaPPer SDSS DR8 catalog with respect to each cluster. We compare our results with an optically selected comparison cluster sample and with theoretical predictions. In contrast to the hypothesis of no environment-based selection, we find a mean 2pcf for background structures of -0.049 on scales of $\lesssim 40'$, significantly non-zero at $\sim 4 σ$, which means that Planck clusters are more likely to be detected in regions of low background density. We hypothesize this effect arises either from background estimation in the SZ survey or from radio sources in the background. We estimate the defect in SZ signal caused by this effect to be negligibly small, of the order of $\sim 10^{-4}$ of the signal of a typical Planck detection. Analogously, there are no implications on X-ray mass measurements. However, the environ- mental dependence has important consequences for weak lensing follow up of Planck galaxy clusters: we predict that projection effects account for half of the mass contained within a 15' radius of Planck galaxy clusters. We did not detect a background underdensity of CMASS LRGs, which also leaves a spatially varying redshift dependence of the Planck SZ selection function as a possible cause for our findings.

preprint2015arXiv

Luminous Red Galaxies in Clusters: Central Occupation, Spatial Distributions, and Mis-centering

Luminous Red Galaxies (LRG) from the Sloan Digital Sky Survey are considered among the best understood samples of galaxies, and they are employed in a broad range of cosmological studies. Because they form a relatively homogeneous population, with high stellar masses and red colors, they are expected to occupy halos in a relatively simple way. In this paper, we study how LRGs occupy massive halos via direct counts in clusters and we reveal several unexpected trends suggesting that the connection between LRGs and dark matter halos may not be straightforward. Using the redMaPPer cluster catalog, we derive the central occupation of LRGs as a function richness, Ncen(λ). Assuming no correlation between cluster mass and central galaxy luminosity at fixed richness, we show that clusters contain a significantly lower fraction of central LRGs than predicted from the two-point correlation function. At halo masses of 10^14.5 Msun, we find Ncen=0.73, compared to Ncen of 0.89 from correlation studies. Our central occupation function for LRGs converges to 0.95 at large halo masses. A strong anti-correlation between central luminosity and cluster mass at fixed richness is required to reconcile our results with those based on clustering studies. We also derive P_BNC, the probability that the brightest cluster member is not the central galaxy. We find P_BNC ~ 20-30% which is a factor of ~2 lower than the value found by Skibba et al. 2011. Finally, we study the radial offsets of bright non-central LRGs from cluster centers and show that bright non-central LRGs follow a different radial distribution compared to red cluster members, which follow a Navarro-Frank-White profile. This work demonstrates that even the most massive clusters do not always have an LRG at the center, and that the brightest galaxy in a cluster is not always the central galaxy.

preprint2014arXiv

Orientation Bias of Optically Selected Galaxy Clusters and Its Impact on Stacked Weak Lensing Analyses

Weak-lensing measurements of the averaged shear profiles of galaxy clusters binned by some proxy for cluster mass are commonly converted to cluster mass estimates under the assumption that these cluster stacks have spherical symmetry. In this paper we test whether this assumption holds for optically selected clusters binned by estimated optical richness. Using mock catalogues created from N-body simulations populated realistically with galaxies, we ran a suite of optical cluster finders and estimated their optical richness. We binned galaxy clusters by true cluster mass and estimated optical richness and measure the ellipticity of these stacks. We find that the processes of optical cluster selection and richness estimation are biased, leading to stacked structures that are elongated along the line-of-sight. We show that weak-lensing alone cannot measure the size of this orientation bias. Weak lensing masses of stacked optically selected clusters are overestimated by up to 3-6 per cent when clusters can be uniquely associated with haloes. This effect is large enough to lead to significant biases in the cosmological parameters derived from large surveys like the Dark Energy Survey, if not calibrated via simulations or fitted simultaneously. This bias probably also contributes to the observed discrepancy between the observed and predicted Sunyaev-Zel'dovich signal of optically-selected clusters.

preprint2014arXiv

redMaPPer III: A Detailed Comparison of the Planck 2013 and SDSS DR8 RedMaPPer Cluster Catalogs

We compare the Planck Sunyaev-Zeldovich (SZ) cluster sample (PSZ1) to the Sloan Digital Sky Survey (SDSS) redMaPPer catalog, finding that all Planck clusters within the redMaPPer mask and within the redshift range probed by redMaPPer are contained in the redMaPPer cluster catalog. These common clusters define a tight scaling relation in the richness-SZ mass ($λ$--$M_{SZ}$) plane, with an intrinsic scatter in richness of $σ_{λ|M_{SZ}} = 0.266 \pm 0.017$. The corresponding intrinsic scatter in true cluster halo mass at fixed richness is $\approx 21\%$. The regularity of this scaling relation is used to identify failures in both the redMaPPer and Planck cluster catalogs. Of the 245 galaxy clusters in common, we identify three failures in redMaPPer and 36 failures in the PSZ1. Of these, at least 12 are due to clusters whose optical counterpart was correctly identified in the PSZ1, but where the quoted redshift for the optical counterpart in the external data base used in the PSZ1 was incorrect. The failure rates for redMaPPer and the PSZ1 are $1.2\%$ and $14.7\%$ respectively, or 9.8% in the PSZ1 after subtracting the external data base errors. We have further identified 5 PSZ1 sources that suffer from projection effects (multiple rich systems along the line-of-sight of the SZ detection) and 17 new high redshift ($z\gtrsim 0.6$) cluster candidates of varying degrees of confidence. Should all of the high-redshift cluster candidates identified here be confirmed, we will have tripled the number of high redshift Planck clusters in the SDSS region. Our results highlight the power of multi-wavelength observations to identify and characterize systematic errors in galaxy cluster data sets, and clearly establish photometric data both as a robust cluster finding method, and as an important part of defining clean galaxy cluster samples.

preprint2014arXiv

redMaPPer IV: Photometric Membership Identification of Cluster Galaxies with 1% Precision

In order to study the galaxy population of galaxy clusters with photometric data one must be able to accurately discriminate between cluster members and non-members. The redMaPPer cluster finding algorithm treats this problem probabilistically. Here, we utilize SDSS and GAMA spectroscopic membership rates to validate the redMaPPer membership probability estimates for clusters with $z\in[0.1,0.3]$. We find small - but correctable - biases, sourced by three different systematics. The first two were expected a priori, namely blue cluster galaxies and correlated structure along the line of sight. The third systematic is new: the redMaPPer template fitting exhibits a non-trivial dependence on photometric noise, which biases the original redMaPPer probabilities when utilizing noisy data. After correcting for these effects, we find exquisite agreement ($\approx 1\%$) between the photometric probability estimates and the spectroscopic membership rates, demonstrating that we can robustly recover cluster membership estimates from photometric data alone. As a byproduct of our analysis we find that on average unavoidable projection effects from correlated structure contribute $\approx 6\%$ of the richness of a redMaPPer galaxy cluster. This work also marks the second public release of the SDSS redMaPPer cluster catalog.

preprint2014arXiv

The X-CLASS - redMaPPer galaxy cluster comparison: I. Identification procedures

We performed a detailed and, for a large part interactive, analysis of the matching output between the X-CLASS and redMaPPer cluster catalogues. The overlap between the two catalogues has been accurately determined and possible cluster positional errors were manually recovered. The final samples comprise 270 and 355 redMaPPer and X-CLASS clusters respectively. X-ray cluster matching rates were analysed as a function of optical richness. In a second step, the redMaPPer clusters were correlated with the entire X-ray catalogue, containing point and uncharacterised sources (down to a few 10^{-15} erg s^{-1} cm^{-2} in the [0.5-2] keV band). A stacking analysis was performed for the remaining undetected optical clusters. Main results show that neither of the wavebands misses any massive cluster (as coded by X-ray luminosity or optical richness). After correcting for obvious pipeline short-comings (about 10% of the cases both in optical and X-ray), ~50% of the redMaPPer (down to a richness of 20) are found to coincide with an X-CLASS cluster; when considering X-ray sources of any type, this fraction increases to ~ 80%; for the remaining objects, the stacking analysis finds a weak signal within 0.5 Mpc around the cluster optical centers. The fraction of clusters totally dominated by AGN-type emission appears to be of the order of a few percent. Conversely ~ 40% of the X-CLASS clusters are identified with a redMaPPer (down to a richness of 20) - part of the non-matches being due to the fact that the X-CLASS sample extends further out than redMaPPer (z<1 vs z<0.6); extending the correlation down to a richness of 5, raises the matching rate to ~ 65%.

preprint2013arXiv

A Maximum Likelihood Approach to Estimating Correlation Functions

We define a Maximum Likelihood (ML for short) estimator for the correlation function, ξ, that uses the same pair counting observables (D, R, DD, DR, RR) as the standard Landy and Szalay (1993, LS for short) estimator. The ML estimator outperforms the LS estimator in that it results in smaller measurement errors at any fixed random point density. Put another way, the ML estimator can reach the same precision as the LS estimator with a significantly smaller random point catalog. Moreover, these gains are achieved without significantly increasing the computational requirements for estimating ξ. We quantify the relative improvement of the ML estimator over the LS estimator, and discuss the regimes under which these improvements are most significant. We present a short guide on how to implement the ML estimator, and emphasize that the code alterations required to switch from a LS to a ML estimator are minimal.

preprint2013arXiv

Cluster Cosmology at a Crossroads: Neutrino Masses

Galaxy clusters --- in combination with CMB and BAO data --- can provide precise constraints on the sum of neutrino masses. However, these constraints depend on the calibration of the mass--observable relation. For instance, the mass calibration employed in Planck Collaboration (2011a,b) rules out the minimal 6-parameter ΛCDM model at 3.7σ, and implies a sum of neutrino masses \sum m_ν= 0.39 \pm 0.10. By contrast, the mass calibration favored by Rozo et al. (2012b) from a self-consistent analysis of X-ray, SZ, and optical scaling relations is consistent with a minimal flat ΛCDM model with no massive neutrinos (1.7σ), and is a better fit to additional data (e.g. H0). We discuss these results in light of the recent SPT and ACT analysis, noting that the Rozo et al. (2012b) model suggests the current mild "tension" (<2σ) between CMB and BAO+H_0 data will decrease as the uncertainties in these measurements decrease.

preprint2013arXiv

Growth of Cosmic Structure: Probing Dark Energy Beyond Expansion

The quantity and quality of cosmic structure observations have greatly accelerated in recent years. Further leaps forward will be facilitated by imminent projects, which will enable us to map the evolution of dark and baryonic matter density fluctuations over cosmic history. The way that these fluctuations vary over space and time is sensitive to the nature of dark matter and dark energy. Dark energy and gravity both affect how rapidly structure grows; the greater the acceleration, the more suppressed the growth of structure, while the greater the gravity, the more enhanced the growth. While distance measurements also constrain dark energy, the comparison of growth and distance data tests whether General Relativity describes the laws of physics accurately on large scales. Modified gravity models are able to reproduce the distance measurements but at the cost of altering the growth of structure (these signatures are described in more detail in the accompanying paper on Novel Probes of Gravity and Dark Energy). Upcoming surveys will exploit these differences to determine whether the acceleration of the Universe is due to dark energy or to modified gravity. To realize this potential, both wide field imaging and spectroscopic redshift surveys play crucial roles. Projects including DES, eBOSS, DESI, PFS, LSST, Euclid, and WFIRST are in line to map more than a 1000 cubic-billion-light-year volume of the Universe. These will map the cosmic structure growth rate to 1% in the redshift range 0<z<2, over the last 3/4 of the age of the Universe.

preprint2013arXiv

Observational Probes of Cosmic Acceleration

The accelerating expansion of the universe is the most surprising cosmological discovery in many decades, implying that the universe is dominated by some form of "dark energy" with exotic physical properties, or that Einstein's theory of gravity breaks down on cosmological scales. The profound implications of cosmic acceleration have inspired ambitious experimental efforts to measure the history of expansion and growth of structure with percent-level precision or higher. We review in detail the four most well established methods for making such measurements: Type Ia supernovae, baryon acoustic oscillations (BAO), weak gravitational lensing, and galaxy clusters. We pay particular attention to the systematic uncertainties in these techniques and to strategies for controlling them at the level needed to exploit "Stage IV" dark energy facilities such as BigBOSS, LSST, Euclid, and WFIRST. We briefly review a number of other approaches including redshift-space distortions, the Alcock-Paczynski test, and direct measurements of H_0. We present extensive forecasts for constraints on the dark energy equation of state and parameterized deviations from GR, achievable with Stage III and Stage IV experimental programs that incorporate supernovae, BAO, weak lensing, and CMB data. We also show the level of precision required for other methods to provide constraints competitive with those of these fiducial programs. We emphasize the value of a balanced program that employs several of the most powerful methods in combination, both to cross-check systematic uncertainties and to take advantage of complementary information. Surveys to probe cosmic acceleration produce data sets with broad applications, and they continue the longstanding astronomical tradition of mapping the universe in ever greater detail over ever larger scales.

preprint2013arXiv

redMaPPer II: X-ray and SZ Performance Benchmarks for the SDSS Catalog

We evaluate the performance of the SDSS DR8 redMaPPer photometric cluster catalog by comparing it to overlapping X-ray and SZ selected catalogs from the literature. We confirm the redMaPPer photometric redshifts are nearly unbiased (<Δz> < 0.005), have low scatter (σ_z ~ 0.006-0.02, depending on redshift), and have a low catastrophic failure rate (~ 1%). Both the T_X-λ and Mgas-λ scaling relations are consistent with a mass scatter of σ_{\ln M|λ} ~ 25%, albeit with a ~ 1% outlier rate due to projection effects. This failure rate is somewhat lower than that expected for the full cluster sample, but is consistent with the additional selection effects introduced by our reliance on X-ray and SZ selected reference cluster samples. Where the redMaPPer DR8 catalog is volume limited (z < 0.35), the catalog is 100% complete above T_X > 3.5 keV, and L_X > 2\times 10^{44} erg/s, decreasing to 90% completeness at L_X ~ 10^{43} erg/s. All rich (λ> 100), low redshift (z < 0.25) redMaPPer clusters are X-ray detected in the ROSAT All Sky Survey (RASS), and 86% of the clusters are correctly centered. Compared to other SDSS photometric cluster catalogs, redMaPPer has the highest completeness and purity, and the best photometric redshift performance, though some algorithms do achieve comparable performance to redMaPPer in subsets of the above categories and/or in limited redshift ranges. The redMaPPer richness is clearly the one that best correlates with X-ray temperature and gas mass. Most algorithms (including redMaPPer) have very similar centering performance, with only one exception which performs worse.

preprint2012arXiv

A Comparative Study of Local Galaxy Clusters: I. Derived X-ray Observables

We examine systematic differences in the derived X-ray properties of galaxy clusters as reported by three different groups: Vikhlinin et al. (2009a), Mantz et al. (2010b), and Planck Collaboration (2011b). The sample overlap between any two pairs of works ranges between 16 to 28 galaxy clusters in common. We find systematic differences in most reported properties, including the total cluster mass, M500. The most extreme case is an average 45% \pm 5% difference in cluster mass between the Planck Collaboration (2011b) and Mantz et al. (2010b), for clusters at z > 0.13 (averaged over 16 clusters). These mass differences induce differences in cluster observables defined within an R500 aperture. After accounting for aperture differences, we find very good agreement in gas mass estimates between the different groups. However, the soft-band X-ray luminosity, LX, core-excised spectroscopic temperature, TX, and gas thermal energy, YX = MgasTX display mean differences at the 5%-15% level. We also find that the low (z \leq 0.13) and high (z \geq 0.13) galaxy cluster samples in Planck Collaboration (2011b) appear to be systematically different: the YSZ/YX ratio for these two sub- samples is ln(YSZ/YX) = -0.06 \pm 0.04 and ln(YSZ/YX) = 0.08 \pm 0.04 for the low and high redshift sub-samples respectively.

preprint2012arXiv

A Comparative Study of Local Galaxy Clusters: II: X-ray and SZ Scaling Relations

We compare cluster scaling relations published for three different samples selected via X-ray and Sunyaev-Zel'dovich (SZ) signatures. We find tensions driven mainly by two factors: i) systematic differences in the X-ray cluster observables used to derive the scaling relations, and ii) uncertainty in the modeling of how the gas mass of galaxy clusters scales with total mass. All scaling relations are in agreement after accounting for these two effects. We describe a multivariate scaling model that enables a fully self-consistent treatment of multiple observational catalogs in the presence of property covariance, and apply this formalism when interpreting published results. The corrections due to scatter and observable covariance can be significant. For instance, our predicted Ysz-Lx scaling relation differs from that derived using the naive "plug in" method by \approx 25%. Finally, we test the mass normalization for each of the X-ray data sets we consider by applying a space density consistency test: we compare the observed REFLEX luminosity function to expectations from published Lx-M relations convolved with the mass function for a WMAP7 flat ΛCDM model.

preprint2012arXiv

Closing the Loop: A Self-Consistent Model of Optical, X-ray, and SZ Scaling Relations for Clusters of Galaxies

We demonstrate that optical data from SDSS, X-ray data from ROSAT and Chandra, and SZ data from Planck, can be modeled in a fully self-consistent manner. After accounting for systematic errors and allowing for property covariance, we find that scaling relations derived from optical and X-ray selected cluster samples are consistent with one another. Moreover, these clusters scaling relations satisfy several non-trivial spatial abundance constraints and closure relations. Given the good agreement between optical and X-ray samples, we combine the two and derive a joint set of LX-M and YSZ-M relations. Our best fit YSZ-M relation is in good agreement with the observed amplitude of the thermal SZ power spectrum for a WMAP7 cosmology, and is consistent with the masses for the two CLASH galaxy clusters published thus far. We predict the halo masses of the remaining z \leq 0.4 CLASH clusters, and use our scaling relations to compare our results with a variety of X-ray and weak lensing cluster masses from the literature.

preprint2012arXiv

Cosmological Constraints from the Large Scale Weak Lensing of SDSS MaxBCG Clusters

We derive constraints on the matter density \Om and the amplitude of matter clustering \sig8 from measurements of large scale weak lensing (projected separation R=5-30\hmpc) by clusters in the Sloan Digital Sky Survey MaxBCG catalog. The weak lensing signal is proportional to the product of \Om and the cluster-mass correlation function \xicm. With the relation between optical richness and cluster mass constrained by the observed cluster number counts, the predicted lensing signal increases with increasing \Om or \sig8, with mild additional dependence on the assumed scatter between richness and mass. The dependence of the signal on scale and richness partly breaks the degeneracies among these parameters. We incorporate external priors on the richness-mass scatter from comparisons to X-ray data and on the shape of the matter power spectrum from galaxy clustering, and we test our adopted model for \xicm against N-body simulations. Using a Bayesian approach with minimal restrictive priors, we find \sig8(\Om/0.325)^{0.501}=0.828 +/- 0.049, with marginalized constraints of \Om=0.325_{-0.067}^{+0.086} and \sig8=0.828_{-0.097}^{+0.111}, consistent with constraints from other MaxBCG studies that use weak lensing measurements on small scales (R<=2\hmpc). The (\Om,\sig8) constraint is consistent with and orthogonal to the one inferred from WMAP CMB data, reflecting agreement with the structure growth predicted by GR for an LCDM cosmological model. A joint constraint assuming LCDM yields \Om=0.298 +/- 0.020 and \sig8=0.831 +/- 0.020. Our cosmological parameter errors are dominated by the statistical uncertainties of the large scale weak lensing measurements, which should shrink sharply with current and future imaging surveys.

preprint2012arXiv

Spatial Anisotropy of Galaxy Kinematics in Sloan Digital Sky Survey Galaxy Clusters

Measurements of galaxy cluster kinematics are important in understanding the dynamical state and evolution of clusters of galaxies, as well as constraining cosmological models. While it is well established that clusters exhibit non-spherical geometries, evident in the distribution of galaxies on the sky, azimuthal variations of galaxy kinematics within clusters have yet to be observed. Here we measure the azimuthal dependence of the line-of-sight velocity dispersion profile in a stacked sample of 1743 galaxy clusters from the Sloan Digital Sky Survey (SDSS). The clusters are drawn from the SDSS DR8 redMaPPer catalog. We find that the line-of-sight velocity dispersion of galaxies lying along the major axis of the central galaxy is larger than those that lie along the minor axis. This is the first observational detection of anisotropic kinematics of galaxies in clusters. We show that the result is consistent with predictions from numerical simulations. Furthermore we find that the degree of projected anisotropy is strongly dependent on the line-of-sight orientation of the galaxy cluster, opening new possibilities for assessing systematics in optical cluster finding.

preprint2012arXiv

The Ysz--Yx Scaling Relation as Determined from Planck and Chandra

SZ clusters surveys like Planck, the South Pole Telescope, and the Atacama Cosmology Telescope, will soon be publishing several hundred SZ-selected systems. The key ingredient required to transport the mass calibration from current X-ray selected cluster samples to these SZ systems is the Ysz--Yx scaling relation. We constrain the amplitude, slope, and scatter of the Ysz--Yx scaling relation using SZ data from Planck, and X-ray data from Chandra. We find a best fit amplitude of \ln (D_A^2\Ysz/CY_X) = -0.202 \pm 0.024 at the pivot point CY_X=8\times 10^{-5} Mpc^2. This corresponds to a Ysz/Yx-ratio of 0.82\pm 0.024, in good agreement with X-ray expectations after including the effects of gas clumping. The slope of the relation is α=0.916\pm 0.032, consistent with unity at \approx 2.3σ. We are unable to detect intrinsic scatter, and find no evidence that the scaling relation depends on cluster dynamical state.

preprint2011arXiv

Cosmological Constraints from Galaxy Clustering and the Mass-to-Number Ratio of Galaxy Clusters

We place constraints on the average density (Omega_m) and clustering amplitude (sigma_8) of matter using a combination of two measurements from the Sloan Digital Sky Survey: the galaxy two-point correlation function, w_p, and the mass-to-galaxy-number ratio within galaxy clusters, M/N, analogous to cluster M/L ratios. Our w_p measurements are obtained from DR7 while the sample of clusters is the maxBCG sample, with cluster masses derived from weak gravitational lensing. We construct non-linear galaxy bias models using the Halo Occupation Distribution (HOD) to fit both w_p and M/N for different cosmological parameters. HOD models that match the same two-point clustering predict different numbers of galaxies in massive halos when Omega_m or sigma_8 is varied, thereby breaking the degeneracy between cosmology and bias. We demonstrate that this technique yields constraints that are consistent and competitive with current results from cluster abundance studies, even though this technique does not use abundance information. Using w_p and M/N alone, we find Omega_m^0.5*sigma_8=0.465+/-0.026, with individual constraints of Omega_m=0.29+/-0.03 and sigma_8=0.85+/-0.06. Combined with current CMB data, these constraints are Omega_m=0.290+/-0.016 and sigma_8=0.826+/-0.020. All errors are 1-sigma. The systematic uncertainties that the M/N technique are most sensitive to are the amplitude of the bias function of dark matter halos and the possibility of redshift evolution between the SDSS Main sample and the maxBCG sample. Our derived constraints are insensitive to the current level of uncertainties in the halo mass function and in the mass-richness relation of clusters and its scatter, making the M/N technique complementary to cluster abundances as a method for constraining cosmology with future galaxy surveys.

preprint2011arXiv

Extrinsic Sources of Scatter in the Richness-Mass Relation of Galaxy Clusters

Maximizing the utility of upcoming photometric cluster surveys requires a thorough understanding of the richness-mass relation of galaxy clusters. We use Monte Carlo simulations to study the impact of various sources of observational scatter on this relation. Cluster ellipticity, photometric errors, photometric redshift errors, and cluster-to-cluster variations in the properties of red-sequence galaxies contribute negligible noise. Miscentering, however, can be important, and likely contributes to the scatter in the richness-mass relation of galaxy maxBCG clusters at the low mass end, where centering is more difficult. We also investigate the impact of projection effects under several empirically motivated assumptions about cluster environments. Using SDSS data and the maxBCG cluster catalog, we demonstrate that variations in cluster environments can rarely (\approx 1% - 5% of the time) result in significant richness boosts. Due to the steepness of the mass/richness function, the corresponding fraction of optically selected clusters that suffer from these projection effects is \approx 5% - 15%. We expect these numbers to be generic in magnitude, but a precise determination requires detailed, survey-specific modeling.

preprint2011arXiv

The Sunyaev-Zeldovich Signal of the maxBCG SDSS Galaxy Clusters in WMAP

The Planck Collaboration measured the Sunyaev-Zel'dovich (SZ) decrement of optically selected clusters from the Sloan Digital Sky Survey, finding that it falls significantly below expectations based on existing mass calibration of the maxBCG galaxy clusters. Resolving this tension requires either the data to go up, or the theoretical expectations to come down. Here, we use data from the Wilkinson Microwave Anisotropy Probe (WMAP) to perform an independent estimate of the SZ decrement of maxBCG clusters. The recovered signal is consistent with that obtained using Planck, though with larger error bars due to WMAP's larger beam size and smaller frequency range. Nevertheless, this detection serves as an independent confirmation of the magnitude of the effect, and demonstrates that the observed discrepancy must be theoretical in origin.

preprint2010arXiv

A GMBCG Galaxy Cluster Catalog of 55,424 Rich Clusters from SDSS DR7

We present a large catalog of optically selected galaxy clusters from the application of a new Gaussian Mixture Brightest Cluster Galaxy (GMBCG) algorithm to SDSS Data Release 7 data. The algorithm detects clusters by identifying the red sequence plus Brightest Cluster Galaxy (BCG) feature, which is unique for galaxy clusters and does not exist among field galaxies. Red sequence clustering in color space is detected using an Error Corrected Gaussian Mixture Model. We run GMBCG on 8240 square degrees of photometric data from SDSS DR7 to assemble the largest ever optical galaxy cluster catalog, consisting of over 55,000 rich clusters across the redshift range from 0.1 < z < 0.55. We present Monte Carlo tests of completeness and purity and perform cross-matching with X-ray clusters and with the maxBCG sample at low redshift. These tests indicate high completeness and purity across the full redshift range for clusters with 15 or more members.

preprint2010arXiv

Annealing a Follow-up Program: Improvement of the Dark Energy Figure of Merit for Optical Galaxy Cluster Surveys

The precision of cosmological parameters derived from galaxy cluster surveys is limited by uncertainty in relating observable signals to cluster mass. We demonstrate that a small mass-calibration follow-up program can significantly reduce this uncertainty and improve parameter constraints, particularly when the follow-up targets are judiciously chosen. To this end, we apply a simulated annealing algorithm to maximize the dark energy information at fixed observational cost, and find that optimal follow-up strategies can reduce the observational cost required to achieve a specified precision by up to an order of magnitude. Considering clusters selected from optical imaging in the Dark Energy Survey, we find that approximately 200 low-redshift X-ray clusters or massive Sunyaev-Zel'dovich clusters can improve the dark energy figure of merit by 50%, provided that the follow-up mass measurements involve no systematic error. In practice, the actual improvement depends on (1) the uncertainty in the systematic error in follow-up mass measurements, which needs to be controlled at the 5% level to avoid severe degradation of the results; and (2) the scatter in the optical richness-mass distribution, which needs to be made as tight as possible to improve the efficacy of follow-up observations.

preprint2010arXiv

Precision Measurements of the Cluster Red Sequence using an Error Corrected Gaussian Mixture Model

The red sequence is an important feature of galaxy clusters and plays a crucial role in optical cluster detection. Measurement of the slope and scatter of the red sequence are affected both by selection of red sequence galaxies and measurement errors. In this paper, we describe a new error corrected Gaussian Mixture Model for red sequence galaxy identification. Using this technique, we can remove the effects of measurement error and extract unbiased information about the intrinsic properties of the red sequence. We use this method to select red sequence galaxies in each of the 13,823 clusters in the maxBCG catalog, and measure the red sequence ridgeline location and scatter of each. These measurements provide precise constraints on the variation of the average red galaxy populations in the observed frame with redshift. We find that the scatter of the red sequence ridgeline increases mildly with redshift, and that the slope decreases with redshift. We also observe that the slope does not strongly depend on cluster richness. Using similar methods, we show that this behavior is mirrored in a spectroscopic sample of field galaxies, further emphasizing that ridgeline properties are independent of environment.

preprint2010arXiv

Stacked Weak Lensing Mass Calibration: Estimators, Systematics, and Impact on Cosmological Parameter Constraints

When extracting the weak lensing shear signal, one may employ either locally normalized or globally normalized shear estimators. The former is the standard approach when estimating cluster masses, while the latter is the more common method among peak finding efforts. While both approaches have identical signal-to-noise in the weak lensing limit, it is possible that higher order corrections or systematics considerations make one estimator preferable over the other. In this paper, we consider the efficacy of both estimators within the context of stacked weak lensing mass estimation in the Dark Energy Survey (DES). We find the two estimators have nearly identical statistical precision, even after including higher order corrections, but that these corrections must be incorporated into the analysis to avoid observationally relevant biases in the recovered masses. We also demonstrate that finite bin-width effects may be significant if not properly accounted for, and that the two estimators exhibit different systematics, particularly with respect to contamination of the source catalog by foreground galaxies. Thus, the two estimators may be employed as a systematics cross-check of each other. Stacked weak lensing in the DES should allow for the mean mass of galaxy clusters to be calibrated to about 2% precision (statistical only), which can improve the figure of merit of the DES cluster abundance experiment by a factor of ~3 relative to the self-calibration expectation. A companion paper (Schmidt & Rozo, 2010) investigates how the two types of estimators considered here impact weak lensing peak finding efforts.

preprint2010arXiv

Weak Lensing Mass Calibration with Shear and Magnification

We study how joint shear and magnification measurements improve the statistical precision of weak lensing mass calibration experiments, relative to standard shear-only analysis. For our magnification measurements, we consider not only the impact of lensing on the source density of galaxies, but also the apparent increase in source sizes. The combination of all three lensing probes - density, size, and shear - can improve the statistical precision of mass estimates by as much as 40% - 50%. This number is insensitive to survey assumptions, though it depends on the degree of knowledge of the parameters controlling the magnification measurements, and on the value of the parameter q that characterizes the response of the source population to magnification. Furthermore, the combination of magnification and shear allows for powerful cross-checks on residual systematics (such as point spread function corrections) at the few percent level.

preprint2010arXiv

Weak Lensing Peak Finding: Estimators, Filters, and Biases

Large catalogs of shear-selected peaks have recently become a reality. In order to properly interpret the abundance and properties of these peaks, it is necessary to take into account the effects of the clustering of source galaxies, among themselves and with the lens. In addition, the preferred selection of lensed galaxies in a flux- and size-limited sample leads to fluctuations in the apparent source density which correlate with the lensing field (lensing bias). In this paper, we investigate these issues for two different choices of shear estimators which are commonly in use today: globally-normalized and locally-normalized estimators. While in principle equivalent, in practice these estimators respond differently to systematic effects such as lensing bias and cluster member dilution. Furthermore, we find that which estimator is statistically superior depends on the specific shape of the filter employed for peak finding; suboptimal choices of the estimator+filter combination can result in a suppression of the number of high peaks by orders of magnitude. Lensing bias generally acts to increase the signal-to-noise νof shear peaks; for high peaks the boost can be as large as Δν~ 1-2. Due to the steepness of the peak abundance function, these boosts can result in a significant increase in the abundance of shear peaks. A companion paper (Rozo et al., 2010) investigates these same issues within the context of stacked weak lensing mass estimates.

preprint2009arXiv

Cosmological Constraints from the SDSS maxBCG Cluster Catalog

We use the abundance and weak lensing mass measurements of the SDSS maxBCG cluster catalog to simultaneously constrain cosmology and the richness--mass relation of the clusters. Assuming a flat \LambdaCDM cosmology, we find σ_8(Ω_m/0.25)^{0.41} = 0.832\pm 0.033 after marginalization over all systematics. In common with previous studies, our error budget is dominated by systematic uncertainties, the primary two being the absolute mass scale of the weak lensing masses of the maxBCG clusters, and uncertainty in the scatter of the richness--mass relation. Our constraints are fully consistent with the WMAP five-year data, and in a joint analysis we find σ_8=0.807\pm 0.020 and Ω_m=0.265\pm 0.016, an improvement of nearly a factor of two relative to WMAP5 alone. Our results are also in excellent agreement with and comparable in precision to the latest cosmological constraints from X-ray cluster abundances. The remarkable consistency among these results demonstrates that cluster abundance constraints are not only tight but also robust, and highlight the power of optically-selected cluster samples to produce precision constraints on cosmological parameters.

preprint2007arXiv

Astrometric Perturbations in Substructure Lensing

In recent years, gravitational lensing has been used as a means to detect substructure in galaxy-sized halos, via anomalous flux ratios in quadruply-imaged lenses. In addition to causing anomalous flux ratios, substructure may also perturb the positions of lensed images at observable levels. In this paper, we numerically investigate the scale of such astrometric perturbations using realistic models of substructure distributions. Substructure distributions that project clumps near the Einstein radius of the lens result in perturbations that are the least degenerate with the best-fit smooth macromodel, with residuals at the milliarcsecond scale. Degeneracies between the center of the lens potential and astrometric perturbations suggest that milliarcsecond constraints on the center of the lensing potential boost the observed astrometric perturbations by an order of magnitude compared to leaving the center of the lens as a free parameter. In addition, we discuss methods of substructure detection via astrometric perturbations that avoid full lens modeling in favor of local image observables and also discuss modeling of systems with luminous satellites to constrain the masses of those satellites.

preprint2007arXiv

MaxBCG: A Red Sequence Galaxy Cluster Finder

Measurements of galaxy cluster abundances, clustering properties, and mass to- light ratios in current and future surveys can provide important cosmological constraints. Digital wide-field imaging surveys, the recently-demonstrated fidelity of red-sequence cluster detection techniques, and a new generation of realistic mock galaxy surveys provide the means for construction of large, cosmologicallyinteresting cluster samples, whose selection and properties can be understood in unprecedented depth. We present the details of the "maxBCG" algorithm, a cluster-detection technique tailored to multi-band CCD-imaging data. MaxBCG primarily relies on an observational cornerstone of massive galaxy clusters: they are marked by an overdensity of bright, uniformly red galaxies. This detection scheme also exploits classical brightest cluster galaxies (BCGs), which are often found at the center of these same massive clusters. (ABRIDGED)

preprint2006arXiv

The Impact of Baryonic Cooling on Giant Arc Abundances

Using ray tracing for simple analytic profiles, we demonstrate that the lensing cross section for producing giant arcs has distinct contributions due to arcs formed through image distortion only, and arcs form from the merging of two or three images. We investigate the dependence of each of these contributions on halo ellipticity and on the slope of the density profile, and demonstrate that at fixed Einstein radius, the lensing cross section increases as the halo profile becomes steeper. We then compare simulations with and without baryonic cooling of the same cluster for a sample of six clusters, and demonstrate that cooling can increase the overall abundance of giant arcs by factors of a few. The net boost to the lensing probability for individual clusters is mass dependent, and can lower the effective low mass limit of lensing clusters. This last effect can potentially increase the number of lensing clusters by an extra 50%. While the magnitude of these effects may be overestimated due to the well known overcooling problem in simulations, it is evident that baryonic cooling has a non-negligible impact on the expected abundance of giant arcs, and hence cosmological constraints from giant arc abundances may be subject to large systematic errors.

preprint2004arXiv

Halo Model Analysis of Cluster Statistics

We use the halo model formalism to provide expressions for cluster abundances and bias, as well as estimates for the correlation matrix between these observables. Off-diagonal elements due to scatter in the mass tracer scaling with mass are included, as are observational effects such as biases/scatter in the data, detection rates (completeness), and false detections (purity). We apply the formalism to a hypothetical volume limited optical survey where the cluster mass tracer is chosen to be the number of member galaxies assigned to a cluster. Such a survey can strongly constrain $σ_8$ ($Δσ_8\approx 0.05$), the power law index $α$ where $<N_{gal}|m>= 1+(m/M_1)^α$ ($Δα\approx0.03$), and perhaps even the Hubble parameter ($Δh\approx 0.07$). We find cluster abundances and bias not well suited for constraining $Ω_m$ or the amplitude $M_1$. We also find that without bias information $σ_8$ and $α$ are degenerate, implying constraints on the former are strongly dependent on priors used for the latter and vice-versa. The degeneracy stems from an intrinsic scaling relation of the halo mass function, and hence it should be present regardless of the mass tracer used in the survey.