Rachel Mandelbaum

Galaxy halo masses and satellite fractions from galaxy–galaxy lensing in the Sloan Digital Sky Survey: stellar mass, luminosity, morphology and environment dependencies

The relationship between galaxies and dark matter (DM) can be characterized by the halo mass of the central galaxy and the fraction of galaxies that are satellites. Here, we present observational constraints from the Sloan Digital Sky Survey on these quantities as a function of r-band luminosity and stellar mass using galaxy-galaxy weak lensing, with a total of 351 507 lenses. We use stellar masses derived from spectroscopy and virial halo masses derived from weak gravitational lensing to determine the efficiency with which baryons in the halo of the central galaxy have been converted into stars. We find that an L* galaxy with a stellar mass of 6 x 1010 M ⊙ is hosted by a halo with mass of 1.4 x 10 12 h -1 M ⊙ , independent of morphology, yielding baryon conversion efficiencies of 17 +10 -5 per cent (early-types) and 16 +15 -6 per cent (late-types) at the 95 per cent confidence level (statistical, not including systematic uncertainty due to assumption of a universal initial mass function). We find that for a given stellar mass, the halo mass is independent of morphology below M stellar = 10 11 M ⊙ , in contrast to typically a factor of 2 difference in halo mass between ellipticals and spirals at a fixed luminosity. This suggests that stellar mass is a good proxy for halo mass in this range and should be used preferentially whenever a halo mass selected sample is needed. For higher stellar masses, the conversion efficiency is a declining function of stellar mass, and the differences in halo mass between early- and late-types become larger, reflecting the fact that most group and cluster haloes with masses above 10 13 M ⊙ host ellipticals at the centre, while even the brightest central spirals are hosted by haloes of mass below 10 13 M ⊙ We find that the fraction of spirals that are satellites is roughly 10-15 per cent independent of stellar mass or luminosity, while for ellipticals this fraction decreases with stellar mass from 50 per cent at 10 10 M ⊙ to 10 per cent at 3 x 10 11 M ⊙ or 20 per cent at the maximum luminosity considered. We split the elliptical sample by local density, and find that at a given luminosity there is no difference in the signal on scales below 100 h -1 kpc between high- and low-density regions, suggesting that tidal stripping inside large haloes does not remove most of the DM from the early-type satellites. This result is dominated by haloes in the mass range 10 13 -10 14 h -1 M ⊙ , and is an average over all separations from the group or cluster centre.

The shear testing programme 2 : factors affecting high-precision weak-lensing analyses.

The Shear Testing Programme (STEP) is a collaborative project to improve the accuracy and reliability of weak-lensing measurement, in preparation for the next generation of wide-field surveys. We review 16 current and emerging shear-measurement methods in a common language, and assess their performance by running them (blindly) on simulated images that contain a known shear signal. We determine the common features of algorithms that most successfully recover the input parameters. A desirable goal would be the combination of their best elements into one ultimate shear-measurement method. In this analysis, we achieve previously unattained discriminatory precision via a combination of more extensive simulations and pairs of galaxy images that have been rotated with respect to each other. That removes the otherwise overwhelming noise from their intrinsic ellipticities. Finally, the robustness of our simulation approach is confirmed by testing the relative calibration of methods on real data. Weak-lensing measurements have improved since the first STEP paper. Several methods now consistently achieve better than 2 per cent precision, and are still being developed. However, we can now distinguish all methods from perfect performance. Our main concern continues to be the potential for a multiplicative shear calibration bias: not least because this cannot be internally calibrated with real data. We determine which galaxy populations are responsible for bias and, by adjusting the simulated observing conditions, we also investigate the effects of instrumental and atmospheric parameters. The simulated point spread functions are not allowed to vary spatially, to avoid additional confusion from interpolation errors. We have isolated several previously unrecognized aspects of galaxy shape measurement, in which focused development could provide further progress towards the sub-per cent level of precision desired for future surveys. These areas include the suitable treatment of image pixellization and galaxy morphology evolution. Ignoring the former effect affects the measurement of shear in different directions, leading to an overall underestimation of shear and hence the amplitude of the matter power spectrum. Ignoring the second effect could affect the calibration of shear estimators as a function of galaxy redshift, and the evolution of the lensing signal, which will be vital to measure parameters including the dark energy equation of state.

Confirmation of general relativity on large scales from weak lensing and galaxy velocities

Although general relativity underlies modern cosmology, its applicability on cosmological length scales has yet to be stringently tested. Such a test has recently been proposed, using a quantity, EG, that combines measures of large-scale gravitational lensing, galaxy clustering and structure growth rate. The combination is insensitive to ‘galaxy bias’ (the difference between the clustering of visible galaxies and invisible dark matter) and is thus robust to the uncertainty in this parameter. Modified theories of gravity generally predict values of EG different from the general relativistic prediction because, in these theories, the ‘gravitational slip’ (the difference between the two potentials that describe perturbations in the gravitational metric) is non-zero, which leads to changes in the growth of structure and the strength of the gravitational lensing effect. Here we report that EG = 0.39 ± 0.06 on length scales of tens of megaparsecs, in agreement with the general relativistic prediction of EG ≈ 0.4. The measured value excludes a model within the tensor–vector–scalar gravity theory, which modifies both Newtonian and Einstein gravity. However, the relatively large uncertainty still permits models within f() theory, which is an extension of general relativity. A fivefold decrease in uncertainty is needed to rule out these models.

Systematic errors in weak lensing: application to SDSS galaxy–galaxy weak lensing

Weak lensing is emerging as a powerful observational tool to constrain cosmological models, but is at present limited by an incomplete understanding of many sources of systematic error. Many of these errors are multiplicative and depend on the population of background galaxies. We show how the commonly cited geometric test, which is rather insensitive to cosmology, can be used as a ratio test of systematics in the lensing signal at the 1 per cent level. We apply this test to the galaxy-galaxy lensing analysis of the Sloan Digital Sky Survey (SDSS), which at present is the sample with the highest weak lensing signal-to-noise ratio and has the additional advantage of spectroscopic redshifts for lenses. This allows one to perform meaningful geometric tests of systematics for different subsamples of galaxies at different mean redshifts, such as brighter galaxies, fainter galaxies and high-redshift luminous red galaxies, both with and without photometric redshift estimates. We use overlapping objects between SDSS and the DEEP2 and 2df-Sloan LRG and Quasar (2SLAQ) spectroscopic surveys to establish accurate calibration of photometric redshifts and to determine the redshift distributions for SDSS. We use these redshift results to compute the projected surface density contrast ΔΣ around 259 609 spectroscopic galaxies in the SDSS; by measuring ΔΣ with different source samples we establish consistency of the results at the 10 per cent level (1σ). We also use the ratio test to constrain shear calibration biases and other systematics in the SDSS survey data to determine the overall galaxy-galaxy weak lensing signal calibration uncertainty. We find no evidence of any inconsistency among many subsamples of the data.

Cosmological parameter constraints from galaxy–galaxy lensing and galaxy clustering with the SDSS DR7

Recent studies have shown that the cross-correlation coefficient between galaxies and dark matter is very close to unity on scales outside a few virial radii of galaxy haloes, independent of the details of how galaxies populate dark matter haloes. This finding makes it possible to determine the dark matter clustering from measurements of galaxy–galaxy weak lensing and galaxy clustering. We present new cosmological parameter constraints based on large-scale measurements of spectroscopic galaxy samples from the Sloan Digital Sky Survey (SDSS) data release 7. We generalize the approach of Baldauf et al. to remove small-scale information (below 2 and 4 h^(−1) Mpc for lensing and clustering measurements, respectively), where the cross-correlation coefficient differs from unity. We derive constraints for three galaxy samples covering 7131 deg^2, containing 69 150, 62 150 and 35 088 galaxies with mean redshifts of 0.11, 0.28 and 0.40. We clearly detect scale-dependent galaxy bias for the more luminous galaxy samples, at a level consistent with theoretical expectations. When we vary both σ_8 and Ω_m (and marginalize over non-linear galaxy bias) in a flat Λ cold dark matter model, the best-constrained quantity is σ_8(Ω_m/0.25)^(0.57) = 0.80 ± 0.05 (1σ, stat. + sys.), where statistical and systematic errors (photometric redshift and shear calibration) have comparable contributions, and we have fixed n_s = 0.96 and h = 0.7. These strong constraints on the matter clustering suggest that this method is competitive with cosmic shear in current data, while having very complementary and in some ways less serious systematics. We therefore expect that this method will play a prominent role in future weak lensing surveys. When we combine these data with Wilkinson Microwave Anisotropy Probe 7-year (WMAP7) cosmic microwave background (CMB) data, constraints on σ_8, Ω_m, H_0, w_(de) and ∑m_ν become 30–80 per cent tighter than with CMB data alone, since our data break several parameter degeneracies.

Detection of large scale intrinsic ellipticity-density correlation from the sloan digital sky survey and implications for weak lensing surveys

The power spectrum of weak lensing shear caused by large-scale structure is an emerging tool for precision cosmology, in particular for measuring the effects of dark energy on the growth of structure at low redshift. One potential source of systematic error is intrinsic alignments of ellipticities of neighbouring galaxies [the intrinsic ellipticity-intrinsic ellipticity (II) correlation] that could mimic the correlations due to lensing. A related possibility pointed out by Hirata & Seljak is correlation between the intrinsic ellipticities of galaxies and the density field responsible for gravitational lensing shear [the gravitational shear-intrinsic ellipticity (GI) correlation]. We present constraints on both the II and GI correlations using 265 908 spectroscopic galaxies from the Sloan Digital Sky Survey (SDSS) and using galaxies as tracers of the mass in the case of the GI analysis. The availability of redshifts in the SDSS allows us to select galaxies at small radial separations, which both reduces noise in the intrinsic alignment measurement and suppresses galaxy-galaxy lensing (which otherwise swamps the GI correlation). While we find no detection of the II correlation, our results are none the less statistically consistent with recent detections found using the SuperCOSMOS survey. Extrapolation of these limits to cosmic shear surveys at z ∼ 1 suggests that the II correlation is unlikely to have been a significant source of error for current measurements of σ 8 with ∼ 10 per cent accuracy, but may still be an issue for future surveys with projected statistical errors below the 1 per cent level unless eliminated using photometric redshifts. In contrast, we have a clear detection of GI correlation in galaxies brighter than L. that persists to the largest scales probed (60 h -1 Mpc) and with a sign predicted by theoretical models. This correlation could cause the existing lensing surveys at z ∼ 1 to underestimate the linear amplitude of fluctuations by as much as 20 per cent depending on the source sample used, while for surveys at z ∼ 0.5 the underestimation may reach 30 per cent. The GI contamination is dominated by the brightest galaxies, possibly due to the alignment of brightest cluster galaxies (BCGs) with the cluster ellipticity due to anisotropic infall along filaments, although other sources of contamination cannot be excluded at this point. We propose that cosmic shear surveys should consider rejection of BCGs from their source catalogues as a test for GI contamination. Future high-precision weak lensing surveys must develop methods to search for and remove this contamination if they are to achieve their promise.

Results of the GREAT08 Challenge: an image analysis competition for cosmological lensing

We present the results of the Gravitational LEnsing Accuracy Testing 2008 (GREAT08) Challenge, a blind analysis challenge to infer weak gravitational lensing shear distortions from images. The primary goal was to stimulate new ideas by presenting the problem to researchers outside the shear measurement community. Six GREAT08 Team methods were presented at the launch of the Challenge and five additional groups submitted results during the 6-month competition. Participants analyzed 30 million simulated galaxies with a range in signal-to-noise ratio, point spread function ellipticity, galaxy size and galaxy type. The large quantity of simulations allowed shear measurement methods to be assessed at a level of accuracy suitable for currently planned future cosmic shear observations for the first time. Different methods perform well in different parts of simulation parameter space and come close to the target level of accuracy in several of these. A number of fresh ideas have emerged as a result of the Challenge including a re-examination of the process of combining information from different galaxies, which reduces the dependence on realistic galaxy modelling. The image simulations will become increasingly sophisticated in future GREAT Challenges, meanwhile the GREAT08 simulations remain as a benchmark for additional developments in shear measurement algorithms.

A halo mass—concentration relation from weak lensing

We perform a statistical weak lensing analysis of dark matter profiles around tracers of halo mass from galaxy-size to cluster-size halos. In this analysis we use 170 640 isolated ∼L∗ galaxies split into ellipticals and spirals, 38 236 groups traced via isolated spectroscopic luminous red galaxies and 13 823 maxBCG clusters from the Sloan Digital Sky Survey covering a wide range of richness. Together these three samples allow a determination of the density profiles of dark matter halos over three orders of magnitude in mass, from 10^12M☉ to 10^15M☉. The resulting lensing signal is consistent with a Navarro–Frenk–White (NFW) or Einasto profile on scales outside the central region. In the inner regions, uncertainty in modeling of the proper identification of the halo center and inclusion of baryonic effects from the central galaxy make the comparison less reliable. We find that the NFW concentration parameter c200b decreases with halo mass, from around 10 for galactic halos to 4 for cluster halos. Assuming its dependence on halo mass in the form of c200b = c0(M/10^14h^−1 M☉)−β we find c0 = 4.6 ± 0.7 (at z = 0.22) and β = 0.13 ± 0.07, with very similar results for the Einasto profile. The slope (β) is in agreement with theoretical predictions, while the amplitude is about two standard deviations below the predictions for this mass and redshift, but we note that the published values in the literature differ at a level of 10–20% and that for a proper comparison our analysis should be repeated in simulations. We compare our results to other recent determinations, some of which find significantly higher concentrations. We discuss the implications of our results for the baryonic effects on the shear power spectrum: since these are expected to increase the halo concentration, the fact that we see no evidence of high concentrations on scales above 20% of the virial radius suggests that baryonic effects are limited to small scales, and are not a significant source of uncertainty for the current weak lensing measurements of the dark matter power spectrum.

Intrinsic galaxy alignments from the 2SLAQ and SDSS surveys: Luminosity and redshift scalings and implications for weak lensing surveys

Correlations between intrinsic shear and the density field on large scales, a potentially important contaminant for cosmic shear surveys, have been robustly detected at low redshifts with bright galaxies in Sloan Digital Sky Survey (SDSS) data. Here we present a more detailed characterization of this effect, which can cause anticorrelations between gravitational lensing shear and intrinsic ellipticity (GI correlations). This measurement uses 36 278 luminous red galaxies (LRGs) from the SDSS spectroscopic sample with 0.15 3σ detections of the effect on large scales (up to 60 h−1 Mpc) for all galaxy subsamples within the SDSS LRG sample; for the 2SLAQ sample, we find a 2σ detection for a bright subsample, and no detection for a fainter subsample. Fitting formulae are provided for the scaling of the GI correlations with luminosity, transverse separation and redshift (for which the 2SLAQ sample, while small, provides crucial constraints due to its longer baseline in redshift). We estimate contamination in the measurement of σ8 for future cosmic shear surveys on the basis of the fitted dependence of GI correlations on galaxy properties. We find contamination to the power spectrum ranging from −1.5 per cent (optimistic) to −33 per cent (pessimistic) for a toy cosmic shear survey using all galaxies to a depth of R = 24 using scales l ≈ 500, though the central value of predicted contamination is −6.5 per cent. This corresponds to a bias in σ8 of Δσ8 = −0.004 (optimistic), −0.02 (central) or −0.10 (pessimistic). We provide a prescription for inclusion of this error in cosmological parameter estimation codes. The principal uncertainty is in the treatment of the L ≤ L blue galaxies, for which we have no detection of the GI signal, but which could dominate the GI contamination if their GI amplitude is near our upper limits. Characterization of the tidal alignments of these galaxies, especially at redshifts relevant for cosmic shear, should be a high priority for the cosmic shear community.

Cross-correlation Weak Lensing of SDSS Galaxy Clusters. I. Measurements

This is the first in a series of papers on the weak lensing effect caused by clusters of galaxies in Sloan Digital Sky Survey. The photometrically selected cluster sample, known as MaxBCG, includes ~130,000 objects between redshift 0.1 and 0.3, ranging in size from small groups to massive clusters. We split the clusters into bins of richness and luminosity and stack the surface density contrast to produce mean radial profiles. The mean profiles are detected over a range of scales, from the inner halo (25 kpc h^(–1)) well into the surrounding large-scale structure (30 Mpc h^(–1)), with a significance of 15 to 20 in each bin. The signal over this large range of scales is best interpreted in terms of the cluster-mass cross-correlation function. We pay careful attention to sources of systematic error, correcting for them where possible. The resulting signals are calibrated to the ~10% level, with the dominant remaining uncertainty being the redshift distribution of the background sources. We find that the profiles scale strongly with richness and luminosity. We find that the signal within a given richness bin depends upon luminosity, suggesting that luminosity is more closely correlated with mass than galaxy counts. We split the samples by redshift but detect no significant evolution. The profiles are not well described by power laws. In a subsequent series of papers, we invert the profiles to three-dimensional mass profiles, show that they are well fit by a halo model description, measure mass-to-light ratios, and provide a cosmological interpretation.