arxiv: v2 [astro-ph.co] 5 Apr 2013

Transcription

1 Draft version July 25, 2018 Preprint typeset using L A TEX style emulateapj v. 8/13/10 arxiv: v2 [astro-ph.co] 5 Apr 2013 THE CLUSTERING OF GALAXIES IN THE SDSS-III BARYON OSCILLATION SPECTROSCOPIC SURVEY: LUMINOSITY AND COLOR DEPENDENCE AND REDSHIFT EVOLUTION Hong Guo 1,2, Idit Zehavi 1, Zheng Zheng 2, David H. Weinberg 3, Andreas A. Berlind 4, Michael Blanton 5, Yanmei Chen 6, Daniel J. Eisenstein 7, Shirley Ho 8,9, Eyal Kazin 10, Marc Manera 11, Claudia Maraston 11,12, Cameron K. McBride 7, Sebastián E. Nuza 13, Nikhil Padmanabhan 14, John K. Parejko 14, Will J. Percival 11, Ashley J. Ross 11, Nicholas P. Ross 8, Lado Samushia 11,15, Ariel G. Sánchez 16, David J. Schlegel 8, Donald P. Schneider 17,18, Ramin A. Skibba 19, Molly E. C. Swanson 7, Jeremy L. Tinker 5, Rita Tojeiro 11, David A. Wake 20,21, Martin White 8,22,23, Neta A. Bahcall 24, Dmitry Bizyaev 25, Howard Brewington 25, Kevin Bundy 26, Luiz N. A. da Costa 27,28, Garrett Ebelke 25, Elena Malanushenko 25, Viktor Malanushenko 25, Daniel Oravetz 25, Graziano Rossi 29,30, Audrey Simmons 25, Stephanie Snedden 25, Alina Streblyanska 31,32, and Daniel Thomas 11,12 1 Department of Astronomy, Case Western Reserve University, OH 44106, USA 2 Department of Physics and Astronomy, University of Utah, UT 84112, USA 3 Department of Astronomy and CCAPP, Ohio State University, Columbus, OH 43210, USA 4 Department of Physics and Astronomy, Vanderbilt University, Nashville, TN 37235, USA 5 Center for Cosmology and Particle Physics, New York University, New York, NY 10003, USA 6 Department of Astronomy, Nanjing University, Nanjing , China 7 Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138, USA 8 Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA 9 Department of Physics, Carnegie Mellon University, 5000 Forbes Avenue, Pittsburgh, PA 15213, USA 10 Center for Astrophysics and Supercomputing, Swinburne University of Technology, P.O. Box 218, Hawthorn, Victoria 3122, Australia 11 Institute of Cosmology & Gravitation, Dennis Sciama Building, University of Portsmouth, Portsmouth, PO1 3FX, UK 12 SEPnet, South East Physics Network, 13 Leibniz-Institut für Astrophysik Potsdam (AIP), An der Sternwarte 16, D Potsdam, Germany 14 Department of Physics, Yale University, 260 Whitney Ave, New Haven, CT 06520, USA 15 National Abastumani Astrophysical Observatory, Ilia State University, 2A Kazbegi Ave., GE-1060 Tbilisi, Georgia 16 Max-Planck-Institut für extraterrestrische Physik, Postfach 1312, Giessenbachstr., D Garching, Germany 17 Department of Astronomy and Astrophysics, The Pennsylvania State University, University Park, PA 16802, USA 18 Institute for Gravitation and the Cosmos, The Pennsylvania State University, University Park, PA 16802, USA 19 Center for Astrophysics and Space Sciences, University of California, 9500 Gilman Drive, San Diego, CA 92093, USA 20 Department of Astronomy, Yale University, New Haven, CT 06520, USA 21 Department of Astronomy, University of Wisconsin-Madison, 475 N. Charter Street, Madison, WI, 53706, USA 22 Department of Physics, University of California, 366 LeConte Hall, Berkeley, CA 94720, USA 23 Department of Astronomy, 601 Campbell Hall, University of California at Berkeley, Berkeley, CA 94720, USA 24 Princeton University Observatory, Princeton, NJ 08544, USA 25 Apache Point Observatory, P.O. Box 59, Sunspot, NM , USA 26 Kavli Institute for the Physics and Mathematics of the Universe, University of Tokyo, Kashiwa, , Japan 27 Observatório Nacional, Rua Gal. José Cristino 77, Rio de Janeiro, RJ , Brazil 28 Laboratório Interinstitucional de e-astronomia - LineA, Rua Gal. José Cristino 77, Rio de Janeiro, RJ , Brazil 29 CEA, Centre de Saclay, Irfu/SPP, F Gif-sur-Yvette, France 30 Paris Center for Cosmological Physics (PCCP) and Laboratoire APC, Université Paris 7, F Paris, France 31 Instituto de Astrofisica de Canarias (IAC), E La Laguna, Tenerife, Spain and 32 Dept. Astrofisica, Universidad de La Laguna (ULL), E La Laguna, Tenerife, Spain Draft version July 25, 2018 Abstract We measure the luminosity and color dependence and the redshift evolution of galaxy clustering in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Ninth Data Release. We focus on the projected two-point correlation function (2PCF) of subsets of its CMASS sample, which includes about 260,000 galaxies over 3,300deg 2 in the redshift range 0.43 < z < 0.7. To minimize the selection effect on galaxy clustering, we construct well-defined luminosity and color subsamples by carefully accounting for the CMASS galaxy selection cuts. The 2PCF of the whole CMASS sample, if approximated by a power-law, has a correlation length of r 0 = 7.93±0.06h 1 Mpc and an index of γ = 1.85 ± Clear dependences on galaxy luminosity and color are found for the projected 2PCF in all redshift bins, with more luminous and redder galaxies generally exhibiting stronger clustering and steeper 2PCF. The color dependence is also clearly seen for galaxies within the red sequence, consistent with the behavior of SDSS-II main sample galaxies at lower redshifts. At a given luminosity (k+e corrected), no significant evolution of the projected 2PCFs with redshift is detected for red sequence galaxies. We also construct galaxy samples of fixed number density at different redshifts, using redshift-dependent magnitude thresholds. The clustering of these galaxies in the CMASS redshift range is found to be consistent with that predicted by passive evolution. Our measurements of the luminosity and color dependence and redshift evolution of galaxy clustering will allow for detailed modeling of the relation between galaxies and dark matter halos and new constraints on galaxy formation and evolution. Subject headings: cosmology: observations cosmology: theory galaxies: distances and redshifts

2 2 Guo et al. galaxies: halos galaxies: statistics large-scale structure of universe 1. INTRODUCTION Galaxies in the universe are observed to display a wide range of properties, such as luminosity, color, stellar mass, age, morphology, and spectral type. These properties encode information about galaxy formation and evolution and are related to the environment hosting the galaxies. Different populations of galaxies are thus expected to trace the underlying dark matter distribution in different ways. Contemporary galaxy redshift surveys, most notably the Sloan Digital Sky Survey (SDSS; York et al. 2000), have transformed the study of large-scale structure, enabling pristine measurements and detailed studies of galaxy clustering. Galaxy clustering provides a powerful approach to probe the complex relation between galaxies and the underlying dark matter distribution (e.g., Kaiser 1984). The dependence of galaxy clustering on galaxy properties has been observed in numerous galaxy surveys (e.g., Davis & Geller 1976; Davis et al. 1988; Hamilton 1988; Loveday et al. 1995; Benoist et al. 1996; Guzzo et al. 1997; Norberg et al. 2001, 2002; Zehavi et al. 2002, 2005b, 2011; Budavári et al. 2003; Madgwick et al. 2003; Li et al. 2006; Coil et al. 2006, 2008; Meneux et al. 2006, 2008; Wang et al. 2007; Wake et al. 2008; Swanson et al. 2008; Meneux et al. 2009; Ross & Brunner 2009; Skibba et al. 2009; Loh et al. 2010; Ross et al. 2010, 2011a; Wake et al. 2011; Christodoulou et al. 2012; Mostek et al. 2012). In general, more luminous and redder galaxies are found to be more strongly clustered than their fainter and bluer counterparts. Similarly, early-type (elliptical) galaxies exhibit stronger clustering than latetype (spiral) ones. Galaxy luminosity and color are perhaps the two major readily observed properties, which also facilitate comparison between different surveys and are less dependent on the stellar evolution models. Moreover, they have proven to be the two properties most predictive of galaxy environment, such that any residual dependence on morphology or surface brightness is weak (Blanton et al. 2005). In this paper, we measure the dependence of galaxy clustering on color and luminosity for galaxies in the SDSS- III Baryon Oscillation Spectroscopic Survey (BOSS; Eisenstein et al.2011; Dawson et al.2013) and study the implications. Galaxy clustering measurements, once properly interpreted, can provide key information about galaxy formation and evolution. In particular, the theoretical understanding of galaxy clustering has been greatly advanced with the development of dark matter halo models (see Cooray & Sheth 2002 and references therein). In the cosmological constant + cold dark matter (ΛCDM) paradigm, galaxies form and evolve in dark matter halos. While the properties and clustering of dark matter halos are well understood with the help of analytic models and numerical simulations (e.g., Mo & White 1996; Springel et al. 2005), the properties of galaxies are hard to predict because of the complex baryonic processes and the lack of complete theory of galaxy formation. Galaxy clustering offers an opportunity to connect galaxies to dark matter halos, providing a new direction in studying galaxy formation and evolution. The halo occupation distribution (HOD) framework (see e.g., Peacock & Smith 2000; Seljak 2000; Scoccimarro et al. 2001; Berlind & Weinberg 2002; Berlind et al. 2003; Zheng et al. 2005) or the closely related conditional luminosity function (CLF) method (Yang et al. 2003, 2005) describe the number of galaxies as a function of halo mass, and galaxy clustering is used to constrain the HOD or CLF parameters. The subhalo abundance matching method makes use of subhalos in high resolution N- body simulations and connects them to galaxies to interpret galaxy clustering (see, e.g., Kravtsov et al. 2004; Conroy et al. 2006; Guo et al. 2010; Nuza et al. 2012). Such models essentially convert galaxy clustering measurements to the relation between galaxies and halos, which provides strong tests of galaxy formation models. The galaxy-halo connections inferred at different redshifts, together with the theoretically known halo evolution, can lead to empirical constraints on galaxy evolution. For example, Zheng et al. (2007) compare the HODs for z 1 DEEP2 and z 0 SDSS galaxies and find a halo mass dependent growth of stellar mass of central galaxies, separating into contributions from star formation and mergers. Tinker & Wetzel (2010) analyze four samples of galaxies from redshift z = 0.4 to z = 2.0. They find that more than 75% of the red satellite galaxies move onto the red sequence because of halo mergers, while the mechanism for central galaxies to move to the red sequence evolves from z = 0.5 to z = 0. To advance our understanding of galaxy evolution, improved measurements of galaxy clustering in large galaxy surveys at different redshifts are necessary. The color and luminosity dependence of clustering has been studied in detail at z 1 for DEEP2 galaxies (Coil et al. 2008) and at z 0 for SDSS galaxies (Zehavi et al. 2011; denoted Z11 hereafter). The trends with color and luminosity are generally similar. However, while Coil et al. (2008) find no changes of the clustering within the red sequence, Z11 find a continuous trend (in both amplitude and slope) of stronger clustering with color in SDSS galaxies. This discrepancy may be caused by different sample selections, but it may also be a signature of galaxy evolution (e.g., related to the buildup of red sequence). Studying a sample from an intermediate redshift range can provide new important information and a better understanding of galaxy evolution. In this paper, we use the recently released CMASS sample of the BOSS Data Release 9 (DR9; Ahn et al. 2012) to measure the clustering of galaxies in the redshift range 0.43 < z < 0.7 and study the dependence on galaxy luminosity and color. The sample is constructed to contain a roughly volume-limited set of massive and luminous galaxies (with a typical stellar mass of h 1 M ; Maraston et al. 2012) in this redshift range. This sample has recently been used to accurately measure the baryon acoustic oscillation (BAO) signature (see Weinberg et al for a recent comprehensive review) on large scales ( 100h 1 Mpc) by Anderson et al. (2012). The sample has been thoroughly vetted, and the robustness of the results and cosmological constraints

3 Luminosity and Color Dependent Clustering 3 from the BAO and redshift-space distortions are explored in a series of papers (Reid et al. 2012; Ross et al. 2012, 2013; Sánchez et al. 2012; Samushia et al. 2013; Scóccola et al. 2012; Tojeiro et al. 2012a,b). The smaller scale clustering measurements and HOD fits are first presented by White et al. (2011) (for an earlier smaller sample) and Nuza et al. (2012). Here we study the small to intermediate scale ( h 1 Mpc) two-point correlation functions (2PCFs) of CMASS galaxies, focusing especially on the dependence on luminosity and color and the implications on galaxy evolution from simple models. We will study the evolution of CMASS galaxies based on HOD modeling in a companion paper. The structure of the paper is organized as follows. In Section 2, we briefly describe the CMASS sample and our method of measuring the 2PCFs. The division to specific subsamples, the clustering measurements and detailed dependence on luminosity and color, and the implications for galaxy evolution are presented in Section 3. We summarize our results in Section 4. Appendix A discusses the effect of different stellar evolution models on our results, and Appendix B explores the robustness of the jackknife error estimates used. Throughout the paper, we assume a spatially flat ΛCDM cosmology as in Anderson et al. (2012), with Ω m = 0.274, h = 0.7, Ω b h 2 = , n s = 0.95, and σ 8 = 0.8, consistent with the best-fit model from the Wilkinson Microwave Anisotropy Probe 7-year data (Komatsu et al. 2011). 2. OBSERVATIONS AND METHODS 2.1. Data As part of the SDSS-III survey, BOSS selects luminous galaxies from the multiple-band SDSS imaging (Fukugita et al.1996; Gunn et al.1998,2006; York et al. 2000) for spectroscopic observation to probe the largescale BAO signals (Anderson et al. 2012). Dawson et al. (2013) provide a comprehensive overview of BOSS, while thetechnicaldetailsofbossarepresentedinsmee et al. (2012) and Bolton et al. (2012). The selection of BOSS galaxies is a union of targets in two different redshift intervals. One is an extension of the SDSS-I/II Luminous Red Galaxy (LRG) sample (Eisenstein et al. 2001), referred to as LOWZ with 0.2 < z < 0.43 (Parejko et al. 2013). The other, denoted as CMASS (White et al. 2011; Anderson et al. 2012), includes 260, 000 galaxies in DR9, and is approximately stellar-mass limited at higherredshifts(0.43 < z < 0.7)withaneffectivevolume of 2.2Gpc 3 andaneffectiveareaofabout3,300deg 2. In thispaper, ourstudyfocusesonthecmasssample. The target selection cuts for CMASS galaxies are fainter and bluer than the LRG sample in order to achieve a higher number density of about h 3 Mpc 3 and sample a wider range of galaxies. The detailed target selection of the CMASS sample is described in Padmanabhan et al. (2013), and summarized in Eisenstein et al. (2011) and Anderson et al. (2012). We briefly describe here the major selection cuts that will affect our analysis of luminosity and color dependence of the 2PCF. The CMASS sample aims at selecting galaxies following a roughly constant stellar mass cut at redshift z > 0.4. The selection criteria of the CMASS galaxies are defined by, 17.5 < i cmod <19.9 (1) d >0.55, (2) i cmod < (d 0.8) (3) i fib2 <21.5 (4) r mod i mod <2.0 (5) where d is defined as d = r mod i mod (g mod r mod )/8 (6) and all magnitudes are extinction corrected (Schlegel et al. 1998) and are in the observed frame. While the magnitudes are calculated using cmodel magnitudes (denoted by the subscript cmod ), the colors are computed using model magnitudes (denoted by the subscript mod ; see Anderson et al for details). The magnitude i fib2 corresponds to the i-band flux within the 2 fiber size. CMASS objects also pass specific star-galaxy separation cuts, as described in Anderson et al. (2012) Selection Cuts and Sample Completeness The CMASS galaxies are chosen by applying complex target selection cuts, and it is thus difficult to construct exact volume-limited samples. Since we intend to investigate the luminosity and color dependence of galaxy clustering in this paper, we must pay particular attention to the completeness in luminosity and color. With this goal in mind, we first investigate the target selection cuts projected to the color-luminosity plane at different redshifts, which will provide the appropriate boundaries in constructing our samples. Figure 1 shows the color-magnitude diagram (CMD) at six narrow redshift ranges with z = The contours represent the number density distribution of CMASS galaxies in the CMD. The absolute magnitude M i and r i color are k+e corrected to z = 0.55 throughout the paper using a global Flexible Stellar Population Synthesis (FSPS) model (Conroy et al. 2009; Conroy & Gunn 2010; Tojeiro et al. 2011, 2012b). Using other stellar evolution models slightly changes the resulting magnitudes and colors, but does not affect our analysis of galaxyclustering, as we discuss in Appendix A. In each panel of Figure 1, only the approximate positions of the selection cuts are shown as dashed lines, since the selection cuts are made in apparent magnitudes and in three bands (g,r, and, i) while our CMD is shown for r i color and M i magnitude. The slopes of the dashed lines follow the boundaries seen in the CMD contours in each redshift bin. We focus on the three main cuts (Equations 1 3). Following Zehavi et al. (2011), we adopt a luminosity-dependent color cut to separate red and blue galaxies (discussed in detail in Section 3.3), (r i) cut = (M i +20), (7) represented by the green lines in Figure 1. The horizontal cut in each panel corresponds to the i-band faint-end flux limit (Equation 1, i < 19.9), which selects galaxies brighter than the corresponding absolute magnitude at each redshift. The slightly tilted vertical dashed line on the left of the distribution is the d cut (Equation 2), which removes galaxies with

4 4 Guo et al. Fig. 1. Color-magnitude diagram (CMD) of CMASS galaxies in different redshift intervals. Both magnitude and color are k+e corrected to z = The contours represent the number density distribution of CMASS galaxies in the CMD. The approximate positions of the selection cuts are shown as dashed lines in each panel, labeled with the corresponding equation numbers of the selection cuts (see the text). The green solid lines are our color cut (Equation 7) for red and blue galaxies. (A color version of this figure is available in the online journal.) bluer colors and with lower redshifts(cannon et al. 2006; Padmanabhan et al. 2013). The bottom-left dashed line is the i-band sliding cut (Equation 3), which excludes the fainter and bluer (thus lower stellar mass) galaxies from the sample. These three main cuts evolve with redshift. A consequence of the d and sliding cuts is that there are more blue galaxies at higher redshifts. At z > 0.55, while the CMASS sample is less affected by the d cut, the sliding cut must be carefully taken into account at all redshifts in constructing complete galaxy samples. At z < 0.55, blue galaxiesarehighlyincomplete (see anestimation of the fraction of star-forming galaxies in CMASS in Figure 16 of Chen et al. 2012). At z < 0.45, the d cut leads to incompleteness even for the most luminous red galaxies. A more sophisticated method to study the sample completeness would be to simulate the galaxy properties (stellar mass, luminosity and color) at different redshifts by assuming certain galaxy stellar evolution models and to apply the selection cuts to simulated galaxies (Swanson et al. 2013). While it would generally depend on the assumptions in the galaxy evolution models, such a study can provide a quantitative estimate of the sample completeness as a function of color and luminosity. In this paper, however, we do not intend to evoke such a model in our clustering analysis. Instead, we empirically use the CMD and selection cuts at each redshift to construct approximately complete galaxy samples. As we proceed in our analysis, we keep in mind the boundaries of completeness defined by the selection cuts when constructing our galaxy samples. One advantage of such an empirical method is that it is largely model-independent Clustering Measurements In this paper, we focus our discussion on the galaxy 2PCF. We use the Landy-Szalay estimator (Landy & Szalay 1993) to measure the 2PCF of galaxies, ξ(r) = DD 2DR+RR (8) RR where DD, DR and RR are the data data, data random, and random random pair counts measured from the data of N galaxies and random samples consisting of N R random points. These pair counts are appropriately normalized by N(N 1)/2, NN R, and N R (N R 1)/2, respectively. We measure the three-dimensional (3D) 2PCF ξ(r p,r π ), where r p and r π are the separations of galaxy pairs perpendicular and parallel to the line of sight. The redshift-space 2PCF ξ(r p,r π ) differs from the realspace one because of redshift distortions induced by galaxy peculiar velocities. The redshift distortions can be mitigated by projecting the 2PCF along the line-of-sight direction, with the projected 2PCF w p (r p ) (Davis & Peebles 1983) defined and measured as w p (r p ) = 2 0 ξ(r p,r π )dr π = 2 i ξ(r p,r π,i ) r π,i (9) where r π,i and r π,i are the i th bin of the line-of-sight

5 Luminosity and Color Dependent Clustering 5 separation and its corresponding bin size. In practice, we sum ξ(r p,r π ) along the line-of-sight direction up to r π,max = 80h 1 Mpc to include most of the correlated pairs. As our analysis focuses on 2PCF measurements up to r p 25h 1 Mpc, the clustering measurements do not depend significantly on the assumed r π,max once it is sufficiently larger. For example, if we integrate the lineof-sight direction to 200h 1 Mpc, the contribution from 80h 1 Mpc < r π < 200h 1 Mpc to w p (r p ) only introduces noisy fluctuations of about 2%. The projected 2PCF can be related to the real-space correlation function, ξ(r), by w p (r p ) = 2 r p rdrξ(r)(r 2 r 2 p ) 1/2 (10) (Davis & Peebles 1983). It is common practice to characterize ξ(r) by a power law on small scales, ξ(r) = (r/r 0 ) γ, and in such a case w p (r p ) can be expressed as ( ) γ r0 w p (r p ) = r p Γ r p ( ) 1 Γ 2 ( γ 1 2 ) / ( γ Γ (11) 2) The covariance matrix of the correlation function is estimated using the jackknife resampling method (following Zehavi et al. 2005b, 2011) Cov(ξ i,ξ j ) = N 1 N N (ξi l ξ i )(ξj l ξ j ) (12) l=1 where N 100 is the total number of jackknife samples, ξi l is the 2PCF in the i th pair separation bin measured from the l th jackknife sample, and ξ i is the average over all samples. We expand on the tests performed by Zehavi et al. (2005b) and investigate the accuracy of the jackknife error estimates using CMASS mock catalogs and find good agreement between the jackknife estimates and those from multiple mocks(see Appendix B for more details). We focus on analyzing our measurement of the galaxy 2PCF on the small to intermediate scales (0.05h 1 Mpc < r p < 25h 1 Mpc), aiming at probing the relation between galaxies and their host halos. However, the small-scale measurement of the 2PCF is limited by the fiber collision effect in the spectrograph system of SDSS-III, where two fibers on the same plate cannot be placed closer than an angular separation of 62 (Dawson et al. 2013; Anderson et al. 2012). As a result, about 5.5% of CMASS galaxies do not have redshifts, strongly affecting small-scale clustering measurements. We correct for this effect using the method proposed and tested by Guo et al. (2012), which divides the galaxy sample into two distinct populations, one free of fiber collisions (referred to as D 1 ) and the other consisting of potentially collided galaxies (denoted as D 2 ). The total clustering signal is a combination of the contributions from these two populations, where the contribution of the collided population is estimated from the resolved galaxies in tile-overlap regions. As discussed in detail in Guo et al. (2012), there are two main systematics that could impact the accuracy of such a method to treat fiber collisions. One is possible density variations between the tile overlap and nonoverlap regions, which is found to be insignificant from tests with mock galaxy catalogs. Another effect is that galaxies in collided triplets (or even higher-order colliding groups) can only be fully recovered in regions covered by at least three tiles, making the estimation from close pairs in two-tile regions not accurate (because of the lack of D 2 D 2 close pairs). This effect is alleviated by an additional correction term using the measured close pairs in the two-tile regions. After full application of the Guo et al.(2012) method, we estimate the remaining systematic errors to be less than 3 5%. It is generally difficult to have an unbiased correction on all scales, but this approach provides the best estimate of the galaxy 2PCF on small scales, compared with other possible methods. When counting pairs for the 2PCF, each galaxy is assigned a series of weights to reduce variance in the measurements and take into account different effects. Following Anderson et al. (2012), we apply a scaleindependent weight to optimize the clustering measurements (Feldman et al. 1994) w FKP = 1 1+ n(z)p 0, (13) where n(z) is the mean density at redshift z, and P 0 = h 3 Mpc 3. This equation provides similar results to the minimum variance weight used by Hamilton (1993) and Zehavi et al. (2002). Another employed weight, w rf, accounts for the fact that not every galaxy with a spectrum taken has a reliable redshift measurement. The redshift failures are dependent on the positions of fibers on the plates (Ross et al. 2012) and are corrected by up-weighting the nearest galaxy that has an accurate redshift. The final weight, w sys, is caused by the scarcity of galaxies detected due to foreground bright stars (Ross et al. 2011a). Ross et al. (2012) present a comprehensive study of potential systematic effects in the 2PCF analysis, and compute a set of weights w sys based on the stellar density and i fib2 magnitude. Therefore, the total weight applied to each galaxy is w tot = w FKP w sys w rf. (14) The quantity w tot is applied to both the D 1 and D 2 populations in the fiber collision correction. The systematic weight w sys only has a small effect on small and intermediate scales, but significantly changes the clustering measurements on BAO scales. We construct the random catalogs according to the detailed angular selection of the DR9 galaxy sample. The radial selection function for each sample is taken into account by assigning the shuffled galaxy redshifts to the random objects. The shuffling method provides a better representation of the true distribution compared with a smooth spline fit to the observed galaxy redshift distribution, as detailed in Ross et al. (2012). This process is done separately for the northern and southern Galactic Caps to account for the different number density distributions (see Anderson et al. 2012, for details). To apply the fiber collision correction, separate random catalogs for the D 1 and D 2 populations are used. Denoting the fraction of recovered D 2 galaxies as N 2 /N 2, for the D 2 random catalogs we apply an additional angular mask N 2/N 2 in each sector (see Guo et al for details).

6 6 Guo et al. 3. RESULTS PCF of the Full CMASS Sample Before presenting the luminosity and color dependence of the 2PCFs for CMASS galaxies, we show in Figure 2 the projected and redshift-space 2PCFs for the entire CMASS sample in the redshift range 0.43 < z < 0.7 (solid lines). For comparison, we also display the recent CMASS DR9 measurements of Nuza et al. (2012, open symbols), which are limited to slightly larger scales ( 0.5h 1 Mpc) and show good agreement on all measured scales. The vertical dashed lines indicate the maximal fiber collision scale, 0.53h 1 Mpc at redshift z = 0.7. By applying the fiber-collision correction method of Guo et al. (2012), we are able to robustly measure the small-scale clustering (note the small error bars). This should enable better constraints on the spatial distribution of galaxies inside dark matter halos and a better constraint on the fraction of satellite galaxies, which we will address in future work. The dotted line in the left panel of Figure 2 is the power-law fit to w p for 0.05h 1 Mpc < r p < 25h 1 Mpc, based on Equation 11, using the full error covariance matrix. Although w p clearly deviates from a power-law (also shown from the χ 2 /dof of the fitting), a power-law fit provides a simple characterization of the clustering and allows for easy comparisons among the 2PCFs of different galaxy samples. The correlation length for the CMASS sample is r 0 = 7.93±0.06h 1 Mpc and the slope γ = 1.85±0.01 (note that given the large bestfit χ 2, the error bars here should be interpreted with care). These values are similar to clustering measurements of LRGs at z 0.55 from the 2dF-SDSS LRG and QSO (2SLAQ) survey(cannon et al. 2006; Ross et al. 2007; Wake et al. 2008), as expected, since the CMASS galaxy color selections were based on the 2SLAQ LRG selection. In order to study the large-scale galaxy bias, we focus on scales larger than typical dark matter halo sizes. By fitting the ratio between the measured galaxy w p and the theoretically computed dark matter w p on scales 3h 1 Mpc < r p < 25h 1 Mpc (detailed in Section 3.5), we obtain a large-scale galaxy bias factor of b = 2.16± 0.01, consistent with the measurements of other authors (White et al.2011; Nuza et al.2012;sánchez et al.2012; Shen et al. 2012). The galaxy bias determined from the projected 2PCF has less scale dependence than that from the redshift-space 2PCF. The exact value of the implied bias can depend on the fitting scales and fitting methods. We find that if the minimal fitting scale is changed from 3h 1 Mpc to 5h 1 Mpc (or larger), the resulting bias only varies slightly at the 2σ level. In the right panel of Figure 2 we present the redshift-space 2PCF ξ(s) on scales above s > 0.1h 1 Mpc, where it can be reliably measured. The fiber collisions in that case impact larger scales than indicated by the dashed line, since s includes the contribution from the line-of-sight separations, s 2 = rp 2 + r2 π (see more discussion in Guo et al. 2012) Luminosity Dependence Luminosity Cuts We now investigate the luminosity dependence of CMASS galaxy clustering. To minimize the influence of sample incompleteness, we carefully construct samples of different luminosities by accounting for the selection cuts as a function of redshifts discussed in Section 2.2. Wedividegalaxiesintotworedshiftbins, 0.43 < z < 0.55 and 0.55 < z < 0.7. The color-magnitude distributions in these two redshift bins, the overall magnitude-redshift distribution, and the cuts used to define our samples are shown in Figure 3. In the right-most panel, some galaxies lie below the faint flux limit (denoted by the lower green line), reflecting the change between the photometry at the time of targeting and that from the final processing. To keep a uniform criterion, we construct our luminosity samples based on the targeting photometry (see details in Anderson et al. 2012). We avoid the impact of the sliding cut by only selecting galaxies brighter than the intersections between the d and sliding cuts. The incompleteness caused by the d cut at z < 0.55 cannot be avoided. Such a limitation means that the blue galaxies are incomplete, especially for the low-redshift samples, while the red galaxies are close to complete for z > 0.5, which is a caveat to remember when interpreting the results. We will thus also study the luminosity dependence limited to the more complete red galaxies. The sliding cut also impacts our ability to study fainter galaxies at high redshift, resulting in the unsampled triangle region above z = 0.55 in the right panel of Figure 3. We have three and two luminosity bin samples at lower and higher redshifts, respectively, each with a bin width of 0.3 magnitude, as shown in the figure. The total numbers of galaxies, N tot, and the comoving volume, V z, in each luminosity bin are shown in Table 1. We also provide the number of blue and red galaxies using the color cut (Equation 7). Powerlaw fits of Equation 11 to the projected 2PCF w p for 0.1h 1 Mpc < r p < 25h 1 Mpc are also given in the table. It is evident that at lower redshifts, red galaxies dominate the samples, and blue galaxies contribute less than 10%. At higher redshifts, approximately 32% of the sample are blue galaxies. This could still be partly related to the CMASS selection effects. If the change in blue galaxy fraction is caused by evolution, i.e., blue galaxies turning red with time, we would expect the luminosity dependence of the 2PCF for red galaxies to also evolve with redshift The Dependence of Galaxy 2PCF on Luminosity The projected 2PCFs of different luminosity samples at the two redshift bins are shown in the top panels of Figure 4. At both redshifts, the luminosity dependence of w p (r p ) is evident, with more luminous galaxies exhibiting stronger clustering, consistent with the results of the SDSS-I/II main sample (Zehavi et al. 2005b, 2011; Li et al. 2006). At z < 0.55, since red galaxiescontribute 90% of the CMASS galaxy population, the luminosity dependence mostly reflects the clustering environment of the red galaxies. At z > 0.55, our measurements of the 2PCF become noisier because of the lower number of galaxies. Even after accounting for the uncertainties in the measurements, the luminosity dependence of clustering is significant in the redshift range of CMASS galaxies. According to the HOD modeling results in SDSS-I/II (Zehavi et al. 2005b, 2011), the increase in the clustering amplitude for more luminous samples reflects the shift in the host halo mass toward the high mass end. The mod-

7 Luminosity and Color Dependent Clustering 7 Fig. 2. Projected (left panel) and redshift-space(right panel) 2PCFs for the entire CMASS sample in the range of 0.43 < z < 0.7. The solid lines present our measurements. The open symbols are the recent measurements of Nuza et al. (2012), which are in good agreement. The vertical dashed lines correspond to the maximal fiber collision scale ( 0.53h 1 Mpc) for z = 0.7. The dotted line in the left panel is a power-law fit to w p for 0.05h 1 Mpc < r p < 25h 1 Mpc with the corresponding parameters as labeled. Fig. 3. Color-magnitude diagrams of CMASS galaxies in the two redshift bins we use, as well as overall distribution of galaxies in i-band absolute magnitude and redshift. The red lines delineate the luminosity bin samples we study. The two green lines in the right panel represent the i-band flux limits of Equation 1, which is also k+e corrected to z = (A color version of this figure is available in the online journal.) TABLE 1 Samples of different luminosities M i range z range N tot N red N blue V z[gpc/h] 3 r 0 γ χ 2 /dof -21.9, , ± ± / , , ± ± / , , ± ± / , , ± ± / , , ± ± /10 Note. r 0 and γ are obtained from fitting a power-law to w p(r p) using the full error covariance matrices for 0.1h 1 Mpc < r p < 25h 1 Mpc. The ratios between χ 2 and degrees-of-freedom (dof) of the fits are also shown.

8 8 Guo et al. Fig. 4. Projected correlation functions, w p(r p), for the various luminosity subsamples at low (top left) and high redshift (top right). The bottom panels present the redshift evolution of w p(r p) in the luminosity interval 22.8 < M i < 22.2 for all the galaxies in the sample (left) and only for the red galaxies (right). Error bars shown are from the diagonal elements of the jackknife covariance matrices. The dotted lines in the top left panel are the power-law fits to the w p in the range of 0.1h 1 Mpc < r p < 2h 1 Mpc to provide a guide of the small-scale slope. (A color version of this figure is available in the online journal.)

9 Luminosity and Color Dependent Clustering 9 eling results in Zehavi et al. (2005b, 2011) alsoshow that the satellite fraction drops as the luminosity of galaxies increases. Our measurements naively seem to support such a result although the 2PCFs becomes generally noisier for higher luminosity samples, the uncertainties in the 2PCFs on small scales increase faster. At such scales, satellite galaxies have a significant contribution to the clustering signal. Therefore, the increase in the measurement errors could be a reflection of the lack of satellites in more luminous samples. The shapesofw p forallluminosity samplesaresimilar. The deviation from a power-law in w p is somewhat more apparent for brighter galaxies, consistent with the results from main galaxies (Zehavi et al. 2004, 2005b, 2011). In thehalomodel, theslopeofw p hasarapidchangearound a few Mpc, indicating the transition from intra-halo galaxy pairs (1-halo term) to inter-halo galaxy pairs (2- halo term). In our measurements, we see that this transition scale increases with increasing luminosity, in agreement with the interpretation that more luminous galaxies reside in more massive (hence larger) halos. The dotted lines in the top left panel of Figure 4 show the power-law fits to w p in the range of 0.1h 1 Mpc < r p < 2h 1 Mpc. There is an apparent weak trend that brighter galaxies have a steeper slope in w p on small scales, in line with the result in Z11. The dotted lines in the top left panel show the power-law fits to the w p in the range of 0.1h 1 Mpc < r p < 2h 1 Mpc, the slope varies from 1.86±0.04(the faintest sample) to 2.13±0.11(the brightest sample). The trend canalsobe interpretedas aresult of the change in the host halo mass scale (see Figure 7 and Appendix A in Zheng et al. 2009). The bottom panels of Figure 4 show the redshift evolution of the 2PCFs of galaxies in a fixed luminosity bin 22.5 < M i < The left panel is for all the galaxies (both blue and red). Galaxies at a lower redshift appear to have a higher clustering strength. This result may be due to the incompleteness of blue galaxies at lower redshifts, and the inclusion of more (less clustered) blue galaxies at higher redshifts. We therefore also examine the redshift evolution of the 2PCF in the 22.5 < M i < 22.2 sample for red galaxies only, as shown in the bottom-right panel of Figure 4, where the red galaxies are defined by the color cut in Equation 7. We find only slight evolution with redshift for the red galaxies (at most 18% in w p for r p > 1h 1 Mpc), not significant within the measurement errors, implying that the differences in the sample of all galaxies (bottom-left panel) are mostly induced by the blue galaxies. In the redshift range of 0.16 < z < 0.44, Zehavi et al. (2005a) find no strong evolution trend in w p (r p ) in the SDSS LRG sample. Their sample of 23.2 < M g < 21.8at 0.16 < z < 0.44 has a comoving number density of n h 3 Mpc 3 (their Figure 2), similar to the number density of our sample of 22.5 < M i < 22.2 at 0.43 < z < 0.7 shown in the bottom right panel of Figure 4. Combining their results with ours, we would inferthat there isnostrongredshift evolutionin w p ofluminousredgalaxiesintheredshiftrangeof0.1 < z < 0.7, consistent with the results of Wake et al. (2008), implying that the effect of structure growth roughly cancels that of evolution of galaxy bias. As will be discussed in Section 3.5, within the errorbars, the trend is also roughly consistent with passive evolution Color Dependence Color Cuts The various target selection cuts make it difficult to discern the red sequence and blue cloud in the CMD of Figure 1. Motivated by the bimodal color distribution of galaxies(e.g., Strateva et al. 2001; Baldry et al. 2004), we construct the division between red and blue galaxies by fitting the galaxy color distribution with two Gaussian distributions. Figure 5 shows the probability distribution functionofr icolor(k+ecorrectedtoz=0.55)insixredshift slices as in Figure 1. The distribution at each redshift is computed from galaxy samples that are as complete in luminosity as possible, i.e., we only use galaxies more luminous than the luminosity given by the intersectionofline(2)andline(3)infigure1, whichcorresponds to an i-band apparent magnitude of i cmod = The CMASS sample shows a clear bimodal distribution in color, similar to the findings in SDSS-I/II(Strateva et al. 2001; Baldry et al. 2004; Skibba & Sheth 2009). We can use the intersection of the best-fit two Gaussian distributions to divide galaxies into blue and red samples. We do not perform the two-gaussian fit to the galaxy distribution in redshift interval 0.4 < z < 0.45, as these galaxy samples are far from complete. The blue galaxies are essentially missing from this sample, and even the red galaxy colors at this redshift are not well described by a Gaussian distribution. As shown in Figure 1, galaxies at this redshift suffer from the d selection cut, which eliminates the blue galaxies and a fraction of red ones. In the 0.45 < z < 0.5 redshift bin, we still miss galaxies with r i < 0.9, and the distribution is dominated by the Gaussian profile from the red population. At 0.5 < z < 0.55, luminous blue galaxies are excluded from the sample by the d cut (see Figure 1), and the contribution to the blue Gaussian profile is mainly from faint blue galaxies. Therefore, at z < 0.55, blue galaxies in the CMASS sample are far from complete. For red galaxies, we find that the centers of the color distribution do not significantly change with redshift. Thus for the analysis of the whole CMASS sample, we use the redshiftindependent color cut for simplicity (see Equation 7). The red/blue color division cut shows a mild redshift dependence, becoming bluer at higher redshift. Since both the color and magnitude used in this paper have been k+e corrected (i.e., the evolution effects are removed), such a mild evolution might indicate that the global evolution correction is not accurate. On the other hand, the photometric errors increase for larger redshift (see below), making the two Gaussian profiles broader, which can lead to a shift of the red-blue division cut towards the blue end even if there is no change in the blue and red populations. Moreover, the blue sample is generally incomplete due to the selection effects, which may also introduce additional change of the color cut. Therefore, the weak dependence of the red-blue division cut on redshift may not reveal much about the evolution. We further investigate the dependence of the color distribution on luminosity for the whole CMASS sample at 0.43 < z < 0.7. As shown in Figure 6, the peak of the red sequence, as well as the intersection of the two Gaussian profiles, become slightly redder as

10 10 Guo et al. 0.4<z< <z< <z< <z< <z< <z<0.7 Fig. 5. Probability distribution function of r i color at different redshift slices, for CMASS samples of high completeness (see text). The black lines are the histogram of r i for all galaxies. The red and blue lines are the bimodality fitting using two Gaussian distributions, with the green curves as their combination. We do not fit the distribution of r i for 0.4 < z < 0.45 because both the red and blue galaxies are far from complete in this redshift interval. (A color version of this figure is available in the online journal.) Fig. 6. Similar to Figure 5, but for different luminosity intervals at 0.43 < z < 0.7. Note that blue galaxies in the left-most panel are not complete. (A color version of this figure is available in the online journal.) the luminosity increases, reflecting the well-known tilt of the red sequence. The tilted red sequence is likely a reflection of differences in the chemical composition, where the more luminous galaxies are richer in metals while the smaller galaxies suffer from the loss of metalenriched gas (Kodama & Arimoto 1997; Gallazzi et al. 2006). The tilted red sequence may also reflect the role of dry mergers(i.e., of gas-poor galaxies) in the evolution of red galaxies (e.g., Skelton et al. 2009), with galaxies increasing their mass (luminosity) from mergers and becoming older (redder) as a result of stellar evolution (see Faber et al for a comprehensive review). Note that in the left-most panel, blue galaxies are not complete in the redshift range of 0.43 < z < 0.7 (see Figures 1 and 3), which leads to the non-monotonic behavior across the three panels. In order to find a reasonable color cut for red and blue galaxies in CMASS, we fit the bimodal color distribution as a function of luminosity for galaxies in the range of 0.5 < z < 0.7, where the samples are less affected by incompleteness. The resulting color cut is the one already presented in Equation 7. With such a color cut, if we

11 Luminosity and Color Dependent Clustering 11 naively count the CMASS galaxies disregarding the incompleteness at lower redshifts, we find that only about 13%ofthegalaxiesin CMASSarebluegalaxies,and80% of these blue galaxies are from z > We emphasize that our color cut is based on the k+e corrected colors. Masters et al. (2011) proposed an observer-frame color cut of g i = 2.35 for CMASS galaxies, motivated by the color and morphology distribution for a matched sample between CMASS and Hubble Space Telescope imaging of the Cosmic Evolution Survey (COSMOS; Scoville et al. 2007). They demonstrated that the g i cut can be used to separate the elliptical and spiral galaxies in CMASS. However, there is still some fraction of LRG progenitors with g i < 2.35 due to small amounts of star formation (Tojeiro et al. 2012b). For our purpose of studying the color dependence of clustering, an observer-frame color cut would mix galaxies from different populations, since galaxies of the same intrinsic colors will appear to have different observed colors at different redshifts. In addition, the g-band magnitude in CMASS is especially faint and thus large measurement errors will make g-band based colors prone to spurious fluctuations. We therefore prefer to use the intrinsic (i.e., k+e corrected) r i color in our study. We note that the z = 0.55 observed r- and i-bands are close to rest-frame g- and r-bands, respectively. So our analysis based on z = 0.55 M i and r i approximates that of rest-frame M r and g r, as adopted in Z11 for z 0 SDSS galaxies. To study the color dependence in different luminosity and redshift intervals in more detail, we further decompose the sample into finer color subsamples. The Gaussian fittings provide the centers and 1σ widths of the blue cloud and red sequence, which are used in defining the fine color cuts (r i) bc = blue center (15) (r i) br = red center 0.5 red width (16) (r i) rr = red center+0.5 red width (17) With the three cuts, we can form blue (below the bc cut), green (between the bc and br cuts), redseq (between the br and rr cuts), and reddest (above the rr cut) samples. In each redshift interval, the luminosity-dependent color cuts are fitted with a straight line, r i = a j M i +b j (18) where j = 1, 2, 3 for the bc, br, and rr cuts, respectively. The linear fits for these cuts are listed in Table 2. For clarity, we show again the CMD in Figure 7, with the fine color cuts superimposed. The redshift interval of 0.4 < z < 0.45 is omitted because of the high sample incompleteness. Examining Figures 5-7, there is a noticeable trend that the width of the red sequence appears narrower for more luminous galaxies and at lower redshifts. We highlight this effect in Figure 8, which presents the 1σ width of the red sequence as a function of luminosity and redshift. Since the photometric errors become larger for fainter galaxies and at higher redshifts, we subtract their contribution in quadrature to obtain the intrinsic color width of the red sequence galaxies. The r i color photometric errors are estimated by simply combining in quadrature the errors in r- and i-band magnitudes (neglecting any additional errors in the k+e corrections and in the correlation between r- and i-band photometric errors). As shownin Figure 8, the abovetrend persistsfor the intrinsic color scatter, suggesting an evolutionary effect which we discuss further in Appendix A The Dependence of Galaxy 2PCF on Color With the color cuts defined in the previous subsection, we investigate the dependence of galaxy 2PCFs on the r i color. First, we examine the 2PCFs for blue and red galaxies in the whole CMASS sample. The red and blue samples here are defined using the color cut in Equation 7. The samples are flux-limited (in addition to other selection cuts) and are by no means complete. The purpose of this exercise is simply to have an overall view of the difference in red and blue galaxy clustering. The ξ(r p,r π ) measurements for blue and red galaxies are shown in Figure 9. For reference, the dotted circles in both panels are the angle-averaged redshift-space correlation function ξ(s). Red galaxies are more strongly clustered. The Fingers-of-God feature (Jackson 1972) on small scales caused by random motions of galaxies in virializedstructurescanbeclearlyseenforbothredandblue galaxies. Red galaxies have a stronger Fingers-of-God effect, reflecting their stronger motions within halos. On large scales (e.g., above r p = 10h 1 Mpc, the outmost contours), the contours for both blue and red galaxies show the flattening trend caused by coherent large-scale infall (Kaiser 1987). On these scales, the Kaiser squashing effect appears to be stronger for blue galaxies, since the effect is determined by Ωm 0.55 /b and blue galaxies have a smaller galaxy bias factor b. We now investigate the color-dependent 2PCFs as a function of luminosity and redshift from the fine color samples. In order to minimize the effect of incompleteness, the luminosity and redshift bins are selected using Figure 1 to make sure that red galaxies are not affected by the selection cuts. The blue galaxies are generally not complete at most redshifts and the results of the blue samples need to be interpreted with care. Nevertheless, the blue samples are still useful in comparison with the red galaxies. The main results of the color-dependent 2PCFs are summarized in Figure 10. The top panels display the dependence of w p on color in the magnitude range 22.5 < M i < 21.5at twodifferent redshift intervals. The trend with color is obvious at both redshifts there is a continuous increase in the clustering amplitude as galaxy color goes from blue to red. This result is consistent with the behavior observed in the SDSS-I/II main galaxy sample (Z11). On small scales (below the inflection scale of 1 2h 1 Mpc), there appears to be a trend that redder galaxies have a steeper slope in w p, which is weaker than that measured by Z11. According to the HOD modeling resultinz11,forgalaxiesinafixedluminosityrange,redder galaxies generally have a higher fraction of satellites residing in massive halos. Our results therefore implies that a larger fraction of redder galaxies reside in more massive halos, giving rise to a larger clustering amplitude. The steepening of w p on smaller scales may also indicate a halo mass scale shift with color, leading to a relative increase in the contribution from the 1-halo central-satellite pairs with respect to the 1-halo satellitesatellite pairs (see the Appendix A of Zheng et al. 2009). The bottom panels of Figure 10 present the 2PCFs