Andreas A. Berlind

The Three-Dimensional Power Spectrum of Galaxies from the Sloan Digital Sky Survey

We measure the large-scale real-space power spectrum P(k) using a sample of 205,443 galaxies from the Sloan Digital Sky Survey, covering 2417 square degrees with mean redshift z~0.1. We employ a matrix-based method using pseudo-Karhunen-Loeve eigenmodes, producing uncorrelated minimum-variance measurements in 22 k-bands of both the clustering power and its anisotropy due to redshift-space distortions, with narrow and well-behaved window functions in the range 0.02 h/Mpc < k < 0.3h/Mpc. We pay particular attention to modeling, quantifying and correcting for potential systematic errors, nonlinear redshift distortions and the artificial red-tilt caused by luminosity-dependent bias. Our final result is a measurement of the real-space matter power spectrum P(k) up to an unknown overall multiplicative bias factor. Our calculations suggest that this bias factor is independent of scale to better than a few percent for k<0.1h/Mpc, thereby making our results useful for precision measurements of cosmological parameters in conjunction with data from other experiments such as the WMAP satellite. As a simple characterization of the data, our measurements are well fit by a flat scale-invariant adiabatic cosmological model with h Omega_m =0.201+/- 0.017 and L* galaxy sigma_8=0.89 +/- 0.02 when fixing the baryon fraction Omega_b/Omega_m=0.17 and the Hubble parameter h=0.72; cosmological interpretation is given in a companion paper.We measure the large-scale real-space power spectrum P(k) by using a sample of 205,443 galaxies from the Sloan Digital Sky Survey, covering 2417 effective square degrees with mean redshift z ≈ 0.1. We employ a matrix-based method using pseudo-Karhunen-Loeve eigenmodes, producing uncorrelated minimum-variance measurements in 22 k-bands of both the clustering power and its anisotropy due to redshift-space distortions, with narrow and well-behaved window functions in the range 0.02 h Mpc-1 < k < 0.3 h Mpc-1. We pay particular attention to modeling, quantifying, and correcting for potential systematic errors, nonlinear redshift distortions, and the artificial red-tilt caused by luminosity-dependent bias. Our results are robust to omitting angular and radial density fluctuations and are consistent between different parts of the sky. Our final result is a measurement of the real-space matter power spectrum P(k) up to an unknown overall multiplicative bias factor. Our calculations suggest that this bias factor is independent of scale to better than a few percent for k < 0.1 h Mpc-1, thereby making our results useful for precision measurements of cosmological parameters in conjunction with data from other experiments such as the Wilkinson Microwave Anisotropy Probe satellite. The power spectrum is not well-characterized by a single power law but unambiguously shows curvature. As a simple characterization of the data, our measurements are well fitted by a flat scale-invariant adiabatic cosmological model with h Ωm = 0.213 ± 0.023 and σ8 = 0.89 ± 0.02 for L* galaxies, when fixing the baryon fraction Ωb/Ωm = 0.17 and the Hubble parameter h = 0.72; cosmological interpretation is given in a companion paper.

Cosmological constraints from the SDSS luminous red galaxies

We measure the large-scale real-space power spectrum P(k) using luminous red galaxies (LRGs) in the Sloan Digital Sky Survey (SDSS) and use this measurement to sharpen constraints on cosmological parameters from the Wilkinson Microwave Anisotropy Probe (WMAP). We employ a matrix-based power spectrum estimation method using Pseudo-Karhunen-Loeve eigenmodes, producing uncorrelated minimum-variance measurements in 20 k-bands of both the clustering power and its anisotropy due to redshift-space distortions, with narrow and well-behaved window functions in the range 0.01h/Mpc 0.1h/Mpc and associated nonlinear complications, yet agree well with more aggressive published analyses where nonlinear modeling is crucial.

The Halo Occupation Distribution: Toward an Empirical Determination of the Relation between Galaxies and Mass

We investigate galaxy bias in the framework of the halo occupation distribution (HOD), which defines the bias of a population of galaxies by the conditional probability P(N|M) that a dark matter halo of virial mass M contains N galaxies, together with prescriptions that specify the relative spatial and velocity distributions of galaxies and dark matter within halos. By populating the halos of a cosmological N-body simulation using a variety of HOD models, we examine the sensitivity of different galaxy clustering statistics to properties of the HOD. The galaxy correlation function responds to different aspects of P(N|M) on different scales. Obtaining the observed power-law form of ξg(r) requires rather specific combinations of HOD parameters, implying a strong constraint on the physics of galaxy formation; the success of numerical and semianalytic models in reproducing this form is entirely nontrivial. Other clustering statistics such as the galaxy-mass correlation function, the bispectrum, the void probability function, the pairwise velocity dispersion, and the group multiplicity function are sensitive to different combinations of HOD parameters and thus provide complementary information about galaxy bias. We outline a strategy for determining the HOD empirically from redshift survey data. This method starts from an assumed cosmological model, but we argue that cosmological and HOD parameters will have nondegenerate effects on galaxy clustering, so that a substantially incorrect cosmological model will not reproduce the observations for any choice of HOD. Empirical determinations of the HOD as a function of galaxy type from the Two-Degree Field (2dF) and Sloan Digital Sky Survey (SDSS) redshift surveys will provide a detailed target for theories of galaxy formation, insight into the origin of galaxy properties, and sharper tests of cosmological models.

THE LUMINOSITY AND COLOR DEPENDENCE OF THE GALAXY CORRELATION FUNCTION

Westudytheluminosityandcolordependenceofthegalaxytwo-pointcorrelationfunctionintheSloanDigitalSky Survey, starting from a sample of � 200,000 galaxies over 2500 deg 2 . We concentrate our analysis on volume-limited subsamples of specified luminosity ranges, for which we measure the projected correlation function wp(rp), which is directly related to the real-space correlation function � (r). The amplitude of wp(rp) rises continuously with luminosity from Mr �� 17: 5t oMr �� 22:5, with the most rapid increase occurring above the characteristic luminosity L� (Mr �� 20:5). Over the scales 0:1 h � 1 Mpc � 22 can be approximated, imperfectly, by power-law three-dimensional correlation functions � (r) ¼ (r/r0) � � with � � 1:8 and r0(L� ) � 5:0 h � 1 Mpc. The brightest subsample, � 23 < Mr < � 22, has a significantly steeper � (r). When we divide samples by color, redder galaxies exhibit a higher amplitude and steeper correlation function at all luminosities. The correlation amplitude of blue galaxies increases continuously with luminosity, but the luminosity dependence for red galaxies is less regular, with bright red galaxies exhibiting the strongest clustering at large scales and faint red galaxies exhibiting the strongest clustering at small scales. We interpret these results using halo occupation distribution (HOD) models assuming concordance cosmological parameters. For most samples, an HOD model with two adjustable parameters fits the wp(rp) data better than a power law, explaining inflections at rp � 1 3 h � 1 Mpc as the transition between the one-halo and two-halo regimes of � (r). The implied minimum mass for a halo hosting a central galaxy more luminous than L grows steadily, with Mmin / L at low luminosities and a steeper dependence above L� . The mass at which a halo has, on average, one satellite galaxy brighter than L is M1 � 23Mmin(L), at all luminosities. These results imply a conditional luminosity function (at fixed halo mass) in which central galaxies lie far above a Schechter function extrapolation of the satellite population. The HOD model fits nicely explain the color dependence of wp(rp) and the cross correlation between red and blue galaxies. For galaxies with Mr < � 21, halos slightly above Mmin have blue central galaxies, while more massive halos have red central galaxies and predominantly red satellite populations. The fraction of blue central galaxies increases steadily with decreasing luminosity and host halo mass. The strong clustering offaint red galaxies follows from the fact that nearly all of them are satellite systems in high-mass halos. The HOD fitting results are in good qualitative agreement with the predictions of numerical and semianalytic models of galaxy formation. Subject headingg cosmology: observations — cosmology: theory — galaxies: distances and redshifts — galaxies: halos — galaxies: statistics — large-scale structure of universe

The halo occupation distribution and the physics of galaxy formation

The halo occupation distribution (HOD) describes the bias between galaxies and dark matter by specifying (1) the probability P(N|M) that a halo of virial mass M contains N galaxies of a particular class and (2) the relative spatial and velocity distributions of galaxies and dark matter within halos. We calculate and compare the HODs predicted by a smoothed particle hydrodynamics (SPH) simulation of a ΛCDM cosmological model (cold dark matter with a cosmological constant) and by a semianalytic galaxy formation model applied to the same cosmology. Although the two methods predict different galaxy mass functions, their HOD predictions for samples of the same space density agree remarkably well. In a sample defined by a baryonic mass threshold, the mean occupation function NM exhibits a sharp cutoff at low halo masses, a slowly rising plateau in which N climbs from 1 to 2 over nearly a decade in halo mass, and a more steeply rising high-occupancy regime at high halo mass. In the low-occupancy regime, the factorial moments N(N - 1) and N(N - 1)(N - 2) are well below the values of N2 and N3 expected for Poisson statistics, with important consequences for the small-scale behavior of the two- and three-point correlation functions. The HOD depends strongly on galaxy age, with high-mass halos populated mainly by old galaxies and low-mass halos by young galaxies. The distribution of galaxies within SPH halos supports the assumptions usually made in semianalytic calculations: the most massive galaxy lies close to the halo center and moves near the halos mean velocity, while the remaining, satellite galaxies have the same radial profile and velocity dispersion as the dark matter. The mean occupation at fixed halo mass in the SPH simulation is independent of the halos larger scale environment, supporting both the merger tree approach of the semianalytic method and the claim that the HOD provides a complete statistical characterization of galaxy bias. We discuss the connections between the predicted HODs and the galaxy formation physics incorporated in the SPH and semianalytic approaches. These predictions offer useful guidance to theoretical models of galaxy clustering, and they will be tested empirically by ongoing analyses of galaxy redshift surveys. By applying the HODs to a large-volume N-body simulation, we show that both methods predict slight departures from a power-law galaxy correlation function, similar to features detected in recent observational analyses.

On departures from a power law in the galaxy correlation function

We measure the projected correlation function wp from the Sloan Digital Sky Survey for a flux-limited sample of 118,000 galaxies and a volume-limited subset of 22,000 galaxies with absolute magnitude Mr M1 = 4.74 ? 1013 h-1 M? is M = 0.89, with 75% of the galaxies residing in less massive, single-galaxy halos and simple auxiliary assumptions about the spatial distribution of galaxies within halos and the fluctuations about the mean occupation. This physically motivated model has the same number of free parameters as a power law, and it fits the wp data better, with a ?2/dof = 0.93, compared to 6.12 (for 10 degrees of freedom, incorporating the covariance of the correlation function errors). Departures from a power-law correlation function encode information about the relation between galaxies and dark matter halos. Higher precision measurements of these departures for multiple classes of galaxies will constrain galaxy bias and provide new tests of the theory of galaxy formation.

Percolation Galaxy Groups and Clusters in the SDSS Redshift Survey: Identification, Catalogs, and the Multiplicity Function

We identify galaxy groups and clusters in volume-limited samples of the Sloan Digital Sky Survey (SDSS) redshift survey, using a redshift-space friends-of-friends algorithm. We optimize the friends-of-friends linking lengths to recover galaxy systems that occupy the same dark matter halos, using a set of mock catalogs created by populating halos of N-body simulations with galaxies. Extensive tests with these mock catalogs show that no combination of perpendicular and line-of-sight linking lengths is able to yield groups and clusters that simultaneously recover the true halo multiplicity function, projected size distribution, and velocity dispersion. We adopt a linking length combination that yields, for galaxy groups with 10 or more members: a group multiplicity function that is unbiased with respect to the true halo multiplicity function; an unbiased median relation between the multiplicities of groups and their associated halos; a spurious group fraction of less than ~1%; a halo completeness of more than ~97%; the correct projected size distribution as a function of multiplicity; and a velocity dispersion distribution that is ~20% too low at all multiplicities. These results hold over a range of mock catalogs that use different input recipes of populating halos with galaxies. We apply our group-finding algorithm to the SDSS data and obtain three group and cluster catalogs for three volume-limited samples that cover 3495.1 deg2 on the sky, go out to redshifts of 0.1, 0.068, and 0.045, and contain 57,138, 37,820, and 18,895 galaxies, respectively. We correct for incompleteness caused by fiber collisions and survey edges and obtain measurements of the group multiplicity function, with errors calculated from realistic mock catalogs. These multiplicity function measurements provide a key constraint on the relation between galaxy populations and dark matter halos.

THE CLUSTERING OF MASSIVE GALAXIES AT z ∼ 0.5 FROM THE FIRST SEMESTER OF BOSS DATA

We calculate the real- and redshift-space clustering of massive galaxies at z ∼ 0.5 using the first semester of data by the Baryon Oscillation Spectroscopic Survey (BOSS). We study the correlation functions of a sample of 44,000 massive galaxies in the redshift range 0.4 <z< 0.7. We present a halo-occupation distribution modeling of the clustering results and discuss the implications for the manner in which massive galaxies at z ∼ 0.5 occupy dark matter halos. The majority of our galaxies are central galaxies living in halos of mass 10 13 h −1 M� ,b ut 10% are satellites living in halos 10 times more massive. These results are broadly in agreement with earlier investigations of massive galaxies at z ∼ 0.5. The inferred large-scale bias (b � 2) and relatively high number density ( ¯ n = 3 × 10 −4 h 3 Mpc −3 ) imply that BOSS galaxies are excellent tracers of large-scale structure, suggesting BOSS will enable a wide range of investigations on the distance scale, the growth of large-scale structure, massive galaxy evolution, and other topics.

The baryonic acoustic feature and large-scale clustering in the Sloan Digital Sky Survey luminous red galaxy sample

We examine the correlation function \xi of the Sloan Digital Sky Survey (SDSS) Luminous Red Galaxy sample (LRG) at large scales (60<s<400 Mpc/h) using the final data release (DR7; 105,831 LRGs between 0.16<z<0.47). Using mock catalogs, we demonstrate that the observed baryonic acoustic peak and larger scale signal are consistent with LCDM at the 1.5\sigma level. The signal at 155<s<200 Mpc/h tends to be high relative to theoretical expectations; this slight deviation can be attributed to a bright subsample of the LRGs. Fitting data to a non-linear, redshift-space, template based-model, we constrain the peak position at s_p=103.6+3.6-2.4 Mpc/h when fitting the range 60<s<150 Mpc/h (1\sigma uncertainties measured from the mocks. This redshift-space distance s_p is related to the comoving sound horizon scale r_s after taking into account matter clustering non-linearities, redshift distortions and galaxy clustering bias. Mock catalogs show that the probability that a DR7-sized sample would not have an identifiable peak is at least 10%. As a consistency check of a fiducial cosmology, we use the observed s_p to obtain the distance D_V=[(1+z)^2D_A^2cz/H(z)]^(1/3) relative to the acoustic scale. We find r_s/D_V(z=0.278)=0.1394+-0.0049. This result is in excellent agreement with Percival et. al (2009), who examine roughly the same data set, but using the power spectrum. Comparison with other determinations in the literature are also in very good agreement. We have tested our results against a battery of possible systematic effects, finding all effects are smaller than our estimated sample variance.We examine the correlation function ? of the Sloan Digital Sky Survey Luminous Red Galaxy sample at large scales (60 h ?1 Mpc < s < 400?h ?1?Mpc) using the final data release (DR7). Focusing on a quasi-volume-limited (0.16 < z < 0.36) subsample and utilizing mock galaxy catalogs, we demonstrate that the observed baryonic acoustic peak and larger scale signal are consistent with ?CDM at 70%-95% confidence. Fitting data to a non-linear, redshift-space, template-based model, we constrain the peak position at s p = 101.7? 3.0?h ?1?Mpc when fitting the range 60 h ?1 Mpc < s < 150?h ?1?Mpc (1? uncertainties). This redshift-space distance s p is related to the comoving sound horizon scale rs after taking into account matter-clustering non-linearities, redshift distortions, and galaxy-clustering bias. Mock catalogs show that the probability that a DR7-sized sample would not have an identifiable peak is at least ~10%. As a consistency check of a fiducial cosmology, we use the observed s p to obtain the distance relative to the acoustic scale. We find rs /DV (z = 0.278) = 0.1389?? 0.0043. This result is in excellent agreement with Percival et al., who examine roughly the same data set, but use the power spectrum. Comparison with other determinations in the literature are also in very good agreement. The signal of the full sample at 125 h ?1 Mpc < s < 200?h ?1?Mpc tends to be high relative to theoretical expectations; this slight deviation can probably be attributed to sample variance. We have tested our results against a battery of possible systematic effects, finding all effects are smaller than our estimated sample variance.

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: a large sample of mock galaxy catalogues

We present a fast method for producing mock galaxy catalogues that can be used to compute the covariance of large-scale clustering measurements and test analysis techniques. Our method populates a second-order Lagrangian perturbation theory (2LPT) matter field, where we calibrate masses of dark matter haloes by detailed comparisons with N-body simulations. We demonstrate that the clustering of haloes is recovered at ∼10 per cent accuracy. We populate haloes with mock galaxies using a halo occupation distribution (HOD) prescription, which has been calibrated to reproduce the clustering measurements on scales between 30 and 80 h−1 Mpc. We compare the sample covariance matrix from our mocks with analytic estimates, and discuss differences. We have used this method to make catalogues corresponding to Data Release 9 of the Baryon Oscillation Spectroscopic Survey (BOSS), producing 600 mock catalogues of the ‘CMASS’ galaxy sample. These mocks have enabled detailed tests of methods and errors, and have formed an integral part of companion analyses of these galaxy data.