arxiv:astro-ph/ v2 18 Aug 2003

Transcription

1 ASTROPHYSICAL JOURNAL SUPPLEMENT SERIES, IN PRESS: AUGUST 4, 2003 Preprint typeset using L A TEX style emulateapj v. 26/01/00 THE OPTICAL AND NEAR-INFRARED PROPERTIES OF GALAXIES: I. LUMINOSITY AND STELLAR MASS FUNCTIONS ERIC F. BELL Max-Planck-Institut für Astronomie, Königstuhl 17, D Heidelberg, Germany; bell@mpia.de DANIEL H. MCINTOSH, NEAL KATZ, AND MARTIN D. WEINBERG Department of Astronomy, University of Massachusetts, 710 North Pleasant Street, Amherst, MA ; dmac@hamerkop.astro.umass.edu, nsk@kaka.astro.umass.edu, weinberg@astro.umass.edu ASTROPHYSICAL JOURNAL SUPPLEMENT SERIES, IN PRESS: August 4, 2003 arxiv:astro-ph/ v2 18 Aug 2003 ABSTRACT We use a large sample of galaxies from the Two Micron All Sky Survey (2MASS) and the Sloan Digital Sky Survey (SDSS) to calculate galaxy luminosity and stellar mass functions in the local Universe. We estimate corrections for passband shifting and galaxy evolution, as well as present-day stellar mass-to-light (M/L) ratios, by fitting the optical near-infrared galaxy data with simple models. Accounting for the 8% galaxy overdensity in the SDSS early data release region, the optical and near-infrared luminosity functions we construct for this sample agree with most recent literature optical and near-infrared determinations within the uncertainties. We argue that 2MASS is biased against low surface brightness galaxies, and use SDSS plus our knowledge of stellar populations to estimate the true K-band luminosity function. This has a steeper faint end slope and a slightly higher overall luminosity density than the direct estimate. Furthermore, assuming a universally-applicable stellar initial mass function (IMF), we find good agreement between the stellar mass function we derive from the 2MASS/SDSS data and that derived by Cole et al. (2001; MNRAS, 326, 255). The faint end slope slope for the stellar mass function is steeper than 1.1, reflecting the low stellar M/L ratios characteristic of low-mass galaxies. We estimate an upper limit to the stellar mass density in the local Universe Ω h = 2.0 ± by assuming an IMF as rich in low-mass stars as allowed by observations of galaxy dynamics in the local Universe. The stellar mass density may be lower than this value if a different IMF with fewer low-mass stars is assumed. Finally, we examine typedependence in the optical and near-infrared luminosity functions and the stellar mass function. In agreement with previous work, we find that the characteristic luminosity or mass of early-type galaxies is larger than for later types, and the faint end slope is steeper for later types than for earlier types. Accounting for typing uncertainties, we estimate that at least half, and perhaps as much as 3/4, of the stellar mass in the Universe is in early-type galaxies. As an aid to workers in the field, we present in an appendix the relationship between model stellar M/L ratios and colors in SDSS/2MASS passbands, an updated discussion of near-infrared stellar M/L ratio estimates, and the volume-corrected distribution of g and K-band stellar M/L ratios as a function of stellar mass. Subject headings: galaxies: luminosity function, mass function galaxies: general galaxies: evolution galaxies: stellar content 1. INTRODUCTION 1 The distribution of galaxy luminosities and stellar masses in the present-day Universe is of fundamental importance for studying the assembly of galaxies over cosmic time, both observationally and theoretically (e.g., Lilly et al. 1995; Lin et al. 1999; Cole et al. 2000; Brinchmann & Ellis 2000; Somerville, Primack, & Faber 2001; Wolf et al. 2003). In addition to providing the zero redshift baseline for luminosity function (LF) evolution, the local LF constrains powerfully much of the important physics affecting the assembly of baryons in dark matter halos. For example, gas accretion and cooling dominates the bright end of the LF, whereas feedback and photoionization affect primarily fainter galaxies (e.g., Cole et al. 2000; Benson et al. 2002). Near-infrared (NIR) luminosities of galaxies are particularly useful as the mass-to-light (M/L) ratios in the NIR vary only by a factor of two or less across a wide range of star formation (SF) histories (Bell & de Jong 2001, see also the Appendix), contrasting with a factor of ten change in M/L ratio at the blue end of the optical regime. Therefore, NIR luminosities provide a cleaner estimate of galaxy stellar masses, which are more robustly predicted by the theoretical models (e.g., Gardner et al. 1997; Cole et al. 2001; Kochanek et al. 2001). The goal of this paper is to use the NIR Two Micron All Sky Survey (2MASS; Skrutskie et al. 1997) in conjunction with optical data and redshifts from the Sloan Digital Sky Survey (SDSS; York et al. 2000) to explore the distribution of galaxy luminosities in the optical and NIR, and to use these data to estimate the distribution of stellar masses in the local Universe. There have been a number of recent studies that have estimated LFs and mass functions (MFs), based on a number of recent large surveys. Around the knee of the LF, which represents the dominant contribution to the overall luminosity density, the agreement between the LFs from different surveys is good. In the optical, luminosity densities agree at typically the 20% level or better, accounting for differences in filter bandpasses and median redshift (e.g., Norberg et al. 2002; Liske et al. 2003; Blanton et al. 2003c). A similar conclusion is found for the NIR K-band (e.g., Gardner et al. 1997; Cole et al. 2001; Kochanek et al. 2001). There are some indications that the behavior of the difficult-to-measure fainter galaxies may depend on environment (Tully et al. 2002), although these galaxies do not exist in sufficient numbers to contribute significantly to the luminosity density of the local Universe (e.g., Zabludoff & Mulchaey 2000; Trentham & Tully 2002). Three notable exceptions to this concordance of recent LF measurements are the Las Campanas Redshift Survey (LCRS; Lin et al. 1996), the early SDSS LF from Blanton et al. (2001),

2 2 Bell et al. and the K-band LF estimate of Huang et al. (2003). The LCRS estimates are relatively consistent with more recent estimates of the optical LFs (e.g., Blanton et al. 2003c), but because of two offsetting effects: (i) the neglect of evolution, which biases the luminosity density to higher values; and (ii) the use of isophotal magnitudes, which biases the luminosity density back down to lower values (Blanton et al. 2001, 2003c). Blanton et al. (2001) find 50% more luminosity density in the local Universe than more recent SDSS or Two Degree Field Galaxy Redshift Survey (2dFGRS; Colless et al. 2001) estimates. This offset is due mostly to the neglect of galaxy evolution and partially to the use of crude k-corrections (Blanton et al. 2003c). The difference between the K-band LF of Huang et al. (2003) and other local estimates is less well-understood, but could stem from the neglect of evolution corrections, LF fitting uncertainties and/or large-scale structure (their LF estimate comes from an area of sky 50 times smaller than the area studied in this work; we discuss this issue in more detail in 4.2). Furthermore, it is unclear if the optical and NIR LFs are mutually consistent. Cole et al. (2001) compared the optical z- band LF from Blanton et al. (2001) with their hybrid J/K-band LF, finding poor agreement. Wright (2001) finds over a factor of two offset between extrapolations from the optical LFs of Blanton et al. (2001) and 2MASS-derived K-band LFs (Cole et al. 2001; Kochanek et al. 2001). Given the above argument that luminosity densities in the optical and NIR are basically known to within 20%, it is unclear whether this discrepancy can be simply accounted for by the neglect of evolution corrections by Blanton et al. (2001), or whether, for example, this is an indication of gross global incompleteness in K-band LFs. Furthermore, the landmark stellar MFs derived by Cole et al. (2001) have not been, as yet, tested systematically. In this paper, the first in a series of papers focusing on the optical and NIR properties of galaxies in the local Universe, we use the NIR 2MASS in conjunction with optical data and redshifts from SDSS to explore in detail the LFs of galaxies over a factor of 6 in wavelength from the u-band (0.35µm) to the K-band (2.15µm). We then, following the methodology of Bell & de Jong (2000, 2001), use the constraints on the optical-nir spectral energy distributions (SEDs) in conjunction with stateof-the-art stellar population synthesis (SPS) models to investigate in detail the stellar MF of galaxies over a factor of 1000 in stellar mass, assuming a universally-applicable stellar initial mass function (IMF). We used these stellar mass estimates in conjunction with a statistically-estimated cold gas masses (HI and H 2 ) to construct the cold baryonic MF in the local Universe and the efficiency of galaxy formation (Bell et al. 2003a, see also, e.g., Salucci & Persic 1999). In subsequent papers we will examine, e.g., the K-band size distribution of galaxies, the photometric properties of a K-selected sample, the dust contents and SF histories of disk galaxies, and the K-band LF of bulges and disks separately, amongst other goals. This paper is arranged as follows. In 2, we discuss the data, focusing on the most important sources of error and incompleteness. In 3, we discuss our method for deriving k- corrections, evolution corrections and stellar M/L ratios. In 4, we construct and discuss optical and NIR LFs for our sample of galaxies. We construct stellar MFs in 5 and discuss these further in 6. We summarize in 7. In the Appendix, we present the distribution of color-derived stellar M/L ratio estimates as a function of galaxy mass and fits to the color- M/L ratio correlations in the SDSS/2MASS passbands as aids to workers in the field. We assume Ω matter = 0.3, Ω Λ = 0.7, and H 0 = 100h km s 1 Mpc 1. For estimating evolution corrections, we assume h = 0.7. Sections 2, 3, and the Appendix go into considerable detail regarding the uncertainties and stellar M/L ratios; thus, readers interested mainly in the results should read 2.1 and then skip directly to THE DATA, DATA QUALITY, AND SELECTION EFFECTS 2.1. Overview We use the SDSS Early Data Release (EDR; Stoughton et al. 2002) to provide a nearly complete 13 r 17.5 sample of galaxies over 414 square degrees with accurate ugriz fluxes and magnitudes. We match these SDSS spectroscopic sample galaxies with the 2MASS extended source catalog (XSC; Jarrett et al. 2000) and point source catalog (PSC) 1. To match the catalogs, we choose the closest galaxy within 2 as the best match (for reference, the random and systematic positional uncertainties of 2MASS and SDSS are 200 and 50 milliarcseconds, respectively; Pier et al. 2003). In this way, we have a reasonably complete 13 r 17.5 sample of galaxies with redshifts, of which have a match in the 2MASS XSC (and therefore have ugrizk fluxes, half-light radii and concentrations in r and K-bands), 6629 of which have a match in the 2MASS PSC (and therefore have ugrizk fluxes, and halflight radii and concentration parameters in r-band), and 3965 of which have no match in either 2MASS catalog (and thus have the optical data only). We choose to use only the 2MASS K-band at the present time. A complete description of these catalogs is far beyond the scope of this paper (see e.g., Jarrett et al. 2000; Blanton et al. 2001; Cole et al. 2001; Stoughton et al. 2002, for more details). Here, we discuss the most important aspects for our purposes: the accuracy of the magnitudes, concentrations and surface brightnesses, and the completeness of the catalogs Magnitude Accuracy An important focus of this paper is the discussion of the offset between optical and 2MASS K-band LFs (e.g., Wright 2001). Because the NIR luminosity density seemed a factor of two below expectations, one of the principal concerns was a large shortfall in either the magnitudes or numbers of galaxies in 2MASS. In this section, we discuss the accuracy of the 2MASS K-band magnitudes in detail, and briefly summarize the expected accuracy of magnitudes in the other passbands. In many respects, K-band data from 2MASS is the ideal tool for constraining galaxy LFs and the stellar MF. 2 K-band galaxy luminosities are five to ten times less sensitive to dust and stellar population effects than optical luminosities, allowing an accurate census of stellar mass in the local Universe (e.g., Bell & de Jong 2001). Furthermore, 2MASS covers the entire sky homogeneously, with 1% systematic variations in zero point (Nikolaev et al. 2000). However, in the NIR the sky background is roughly a factor of 100 times brighter than the mean surface brightness of luminous galaxies, and the exposure time of 2MASS is short (7.8 seconds with a 1.3-m telescope; Skrutskie et al. 1997). Thus, low surface brightness (LSB) features, such as LSB galaxies or the outer regions of normal galaxies, may be missed by 2MASS We note that, strictly speaking, 2MASS adopts a K s-band that peaks at rather shorter wavelengths than the standard K-band, but we we will call it K-band for brevity in this paper.

3 Galaxy Luminosities and Stellar Masses 3 FIG. 1. Comparison of 2MASS K-band Kron magnitudes with deep K- band magnitudes for 223 galaxies from Loveday (2000). The lower panel shows the difference (2MASS Loveday) in K-band magnitudes. The bold line represents the mean magnitude difference 0.01 ± 0.04 (0.10 ± 0.02) for K < 11 (K 11). The scatter is 0.2 mag. FIG. 2. Comparison of 2MASS K-band Kron magnitudes with deep K- band magnitudes (transformed from H-band using the 2MASS measured H K color) for 1017 galaxies from papers by Gavazzi and coworkers. The lower panel shows the difference (2MASS Gavazzi) in K-band magnitudes. The bold line represents the mean magnitude difference 0.06 ± 0.02 (0.17 ± 0.01) for K < 10 (K 10). The scatter is 0.2 mag. To test how much light 2MASS misses in the LSB outer parts of galaxies, we compare 2MASS K-band magnitudes from the XSC with K-band magnitudes from deeper imaging data. Following Cole et al. (2001), we show Loveday s K-band Kron 3 magnitude from relatively deep data (10 min- 3 Kron (1980) magnitudes are measured in apertures that are related to the utes on the Cerro Tololo International Observatory 1.5-m telescope) against 2MASS K-band Kron magnitudes (Fig. 1). At K < 11, 2MASS Kron magnitudes seem quite accurate, with a systematic offset of 0.01±0.04 mag. At fainter magnitudes, 2MASS Kron magnitudes underestimate the true magnitude by 0.10±0.02 mag (scatter 0.2 mag). Cole et al. (2001) found a larger offset between Second Incremental Data Release 2MASS K-band and total magnitudes; since the Second Incremental Data Release there have been improvements to the reduction pipeline that have improved the quality of 2MASS K-band Kron magnitudes. We check this offset by comparison with a larger sample of galaxies imaged in the H-band by Gavazzi et al. (1996a,b, 2000) and Boselli et al. (2000) 4. We adopt the 2MASS H K color to estimate the total K-band magnitude; the typical value is H K 0.25, almost independent of galaxy type. The average offset brighter (fainter) than K = 10 mag is 0.06 (0.17) mag, in the sense that Gavazzi s magnitudes are slightly brighter than the 2MASS Kron magnitudes (see Fig. 2). We do not adopt this correction in this paper owing to uncertainties in transforming H-band data into K-band. We do, nevertheless, choose to adopt an offset of 0.1 mag for all galaxies (not just galaxies with K 11) to better match Gavazzi s offset. We note that magnitudes corrected in this way will be within 0.1 mag of total, independent of whether one compares them to Gavazzi s or Loveday s total magnitudes. We adopt a 0.1 mag uncertainty in the correction to total K-band fluxes, added in quadrature with the 2MASS random magnitude error. We tested whether the correction to total is a function of K-band surface brightness, as one could imagine that the fraction of light lost may be larger for lower surface brightness galaxies. We found no correlation between the correction to total flux and K-band surface brightness within the errors, supporting our use of a blanket 0.1 mag offset. We also choose to match to the 2MASS PSC. There are very few matches to the comparison samples: 8 from Loveday (2000) and 8 from the sample from Gavazzi and coworkers. We find mean offsets of 0.9 mag (0.3 mag RMS) and 0.8 mag (0.3 mag RMS) for the two samples. We disregard two outliers (with no offset and a 3.1 mag offset) from the eight of the Gavazzi sample. We account for the large PSC offsets by subtracting 0.85 mag from the 2MASS PSC K-band magnitudes, and setting their K-band errors to 0.5 mag. These magnitudes are clearly of very limited use. We use them primarily to constrain only roughly the k-correction, evolution correction, and stellar M/L ratio estimates. In particular, our choice of K- band magnitude limit (extinction-corrected K-band Kron magnitude of 13.57, with the offset included after galaxy selection) includes only 66 galaxies from the PSC, or just over 1% of our K-band selected sample. Because of its high signal-to-noise, SDSS Petrosian 5 ugriz magnitudes are expected to be accurate to better than 0.05 mag in a random and systematic sense (Strauss et al. 2002; Blanton et al. 2003c). Sloan papers typically make the distinction between preliminary magnitudes presented by the EDR in the natural Sloan 2.5-m telescope system, denoted u g r i z, and galaxy radius (for 2MASS, not less than 5 ). 4 This sample was used to test the circular isophotal magnitudes used by Kochanek et al. (2001) for their K-band derived LF. 5 SDSS Petrosian magnitudes are estimated within an aperture that is twice the radius at which the local surface brightness is 1/5 of the mean surface brightness within that radius (Strauss et al. 2002).

4 4 Bell et al. the true Sloan magnitudes ugriz. We denote the EDR Petrosian magnitudes ugriz for brevity. Petrosian magnitudes of well-resolved early-type galaxies (with close to R 1/4 law luminosity profiles) underestimate the total flux by 0.1 mag because their surface brightness profiles fall off very slowly at large radii (Strauss et al. 2002; Blanton et al. 2003c). In this paper, we crudely correct for this effect by subtracting 0.1 mag from the magnitude of any galaxy with an r-band concentration parameter of c r > 2.6 (defined in the next section). While simplistic, it allows us to estimate the total fluxes for early-type galaxies to within 0.05 mag. We adopt a magnitude error of 0.05 mag for all galaxies, added in quadrature to the (tiny) internal SDSS random magnitude errors. Note that we apply all magnitude offsets after galaxy selection Concentration Parameters and Surface Brightnesses In this paper, we study primarily the overall luminosities and stellar masses of galaxies, choosing not to focus on their structural parameters, such as concentration parameter or surface brightness. Nevertheless, we do use concentration parameter as a crude discriminant between early and late-type galaxies and surface brightnesses when examining the completeness of the galaxy samples. We adopt as our primary morphological classifier the r-band concentration parameter, c r = r 90 /r 50, where r 90 and r 50 are the circular aperture radii within which 90% and 50% of the Petrosian flux are contained, respectively. The concentration parameter has been extensively used within the SDSS collaboration to separate between early and late-type galaxies in a rudimentary fashion; early-type galaxies have higher c r than later types. This is motivated by the work of Strateva et al. (2001) and Shimasaku et al. (2001), who find a scattered but reasonable correlation between qualitative morphological classifications and c r. Strateva et al. (2001) suggest a c r 2.6 selection for early-type galaxies and this cut has been adopted by Kauffmann et al. (2003a). We also adopt this criterion, primarily because it is easily reproducible, facilitating easy comparison with our results by other workers. Blanton et al. (2003b) note the sensitivity of the the concentration parameter to seeing; more heavily smoothed early-type galaxies appear less concentrated than they would be either if they were observed with better seeing or were closer. Because of this, the early type definition is conservative; intrinsically smaller or more distant early types maybe misclassified as later types owing to seeing effects. 2MASS also gives a concentration parameter c K = r 75 /r 25, where r 75 and r 25 are the elliptical aperture within which 75% and 25% of the flux are contained. In Fig. 3, we compare the SDSS c r to the 2MASS c K in an effort to explore systematic uncertainties in the use of concentration parameters as a morphological typing tool (we compare c r to color selection later also in 4.4). We restrict our comparison to K < 12 galaxies, which have sufficient S/N to estimate c K. It is clear that there are systematic differences between the two definitions, which are manifested by zero point shifts, a non-unity slope, and a substantial scatter. Nevertheless, making the crude approximation that c K c r + 1 (the solid line), we can compare the fractions classified as early-type with both definitions. Using c r 2.6, we find that 401/603 K < 12 EDR galaxies are classified as early-type. Using c K 3.6, we find that 439/603 galaxies are classified as early-type. Furthermore, 355/401 galaxies (89%±7%) of r-band classified early-types are classified as early type using the K-band classification. FIG. 3. K-band concentration parameter c K = r 75 /r 25 against r-band concentration parameter c r = r 90 /r 50. The solid line denotes c K c r + 1, and the dashed lines denote the two rough early-type cuts at c r 2.6 and c K 3.6. We use g, r, and K-band surface brightnesses only as a rough check on the completeness properties of the 2MASS-matched sample and on the K-band LF. These surface brightnesses are defined to be the average surface brightness within the halflight radii. Since the magnitudes for SDSS and 2MASS XSC are accurate to at worst 20%, we expect that the half-light surface brightnesses will be accurate to 30%, given the 20% error in total magnitude, added in quadrature with the effect of a 20% scale size error, which is typical of scale-size comparisons between different authors (see, e.g., Bell & de Jong 2000). Accuracies of this order are more than sufficient for our present purposes Completeness The homogeneity and completeness of SDSS and 2MASS make them powerful tools for understanding the characteristics of galaxies in the local Universe. To construct meaningful LFs from these datasets, we must understand the completeness characteristics of each survey. We choose SDSS EDR spectroscopic sample galaxies with Galactic foreground extinction-corrected 13 r 17.5 (Schlegel, Finkbeiner, & Davis 1998), following Stoughton et al. (2002). This galaxy sample is nearly complete, as discussed in much more detail by e.g. Blanton et al. (2001), Stoughton et al. (2002) or Strauss et al. (2002). We find that the area covered by the SDSS spectroscopic sample is 414 square degrees, 90% of the 462 square degrees covered by the EDR imaging data (Stoughton et al. 2002). This is quite consistent with the the statement by Stoughton et al. (2002) that only 93% of the spectroscopic tiles were attempted; we adopt the 3% difference as our systematic error in determining the sky coverage of this sample. We estimate a total completeness within this area by querying SDSS EDR photometric catalog galaxies satisfying the spectroscopic catalog inclusion criteria as outlined by Strauss et al. (2002). This value is 78%, which is consistent with a 2% loss of galaxies due to bright stars, a > 99% redshift success rate, and between 80% and 90% targeting efficiency

5 Galaxy Luminosities and Stellar Masses 5 (e.g., Blanton et al. 2001; Stoughton et al. 2002). A value of 85% was recently found by Nakamura et al. (2003) for bright SDSS galaxies (r 15.9) in the EDR: we adopt the difference between our and Nakamura et al. s measurements as the systematic error in the completeness, which is propagated through into the φ and j estimates later. We do not take account of the detailed, position-dependent completeness of the sample. While a detailed accounting for the completeness as a function of position is pivotal for estimating galaxy clustering properties, it is of only minor importance for estimating the overall LF. Finally, we note that there is little systematic bias within SDSS against galaxies within the selection limits. Of order 0.1% of the lowest surface brightness galaxies are not targeted because a spectrum would be impossible to obtain, and because over 3/4 of the lowest surface brightness features in the SDSS imaging survey are artifacts (Strauss et al. 2002). Also, 5% of bright galaxies are rejected because they overlap a bright, saturated star, or because they have a very bright fiber magnitude and are not targeted to avoid severe cross-talk between the fiber spectra. Neither of these biases will significantly affect our analysis. In addition to estimating the completeness of SDSS internally, we determine whether the SDSS EDR area is overdense using the full coverage of 2MASS. We estimate overdensities by comparing the number of 2MASS extended sources with 10 < K < 13.5 in the sky outside of the Galactic Plane ( b 30 ) with the number of similar sources in the SDSS EDR region. We use an area that is slightly less than the 414 square degrees that we calculate for the spectroscopic EDR coverage because we choose rectangular areas that are fully enclosed by the SDSS EDR boundaries. We show in Fig. 4 that the EDR is overdense over the entire magnitude range 10 < K < We include the estimated density in the 2dFGRS region for comparison and give the number counts for each region in Table 1. The SDSS EDR is 8% overdense (with a 1% Poisson uncertainty), and the 2dFGRS region (used by Cole et al. 2001) is 2% underdense, compared to the whole sky. Although this estimate is admittedly rough because 10 < K < 13.5 galaxies are a somewhat different set of galaxies than those with 13 r 17.5, the overdensity estimate is accurate given that we compare to half of the entire sky ( b 30 ). Furthermore, our estimate is insensitive to Galactic foreground extinction. We account for the EDR region overdensity in our analysis by multiplying the effective survey area by 1.08 when constructing our LFs. We focus our study on the 2MASS matches to the 13 r 17.5 SDSS catalog. As stated earlier, out of the r 17.5 galaxies in the EDR spectroscopic sample, we match galaxies in the 2MASS XSC and 6629 galaxies in the 2MASS PSC. In Fig. 5, we explore the properties of the 2MASS matched and unmatched galaxies in more detail. In the upper panels, we show the distribution of galaxy g r color (left) and r-band surface brightness µ r (right). The solid his- TABLE 1 10 < K < 13.5 GALAXY NUMBER COUNTS Region N Area (deg 2 ) n gal (deg 2 ) XSC b dFGRS Sloan EDR tograms show galaxies with 2MASS counterparts, the dashed histograms show those without. We give the K-band surface brightness µ K (left) and apparent magnitude K (right) in the lower panels. We estimate µ K and K for galaxies that have no 2MASS data using the SDSS µ r and r-band apparent magnitude in conjunction with the r K color of the best-fit SED model (as described in 3). We test this procedure by using the optical data only to predict the K-band magnitudes of the galaxies with K-band XSC data. We find that this procedure is accurate to 0.4 mag RMS. We see that the galaxies that are unmatched in 2MASS are preferentially blue and LSB in the optical and NIR. There are 84 LSB galaxies ( µ K 19.1 mag arcsec 2 ) with estimated K < 13.57, thus, there may be a small population of LSB galaxies missed by 2MASS. Faint, LSB galaxies are visible only in the very nearest parts of an apparent magnitude-limited survey (such as 2MASS), and therefore carry a large weight 1/V max. Therefore, this small bias ( 1%) may translate into a larger bias when considering the LF or luminosity density. This bias would affect all published 2MASS LFs (e.g., Cole et al. 2001; Kochanek et al. 2001), as well as our own. We show later that this bias affects the faint end of the LF, as one would expect given the surface brightness dependence of the LF (see, e.g., de Jong & Lacey 2000; Cross & Driver 2002). We also estimate the degree of incompleteness using the optical data in conjunction with our knowledge of stellar populations to push the K-band LF and stellar MF down to lower galaxy masses. We select samples for estimating LFs in different passbands using passband-dependent magnitude limits, following Blanton et al. (2003c). Specifically, when constructing ugizk LFs, we select the magnitude limit in ugizk so that the V max for each galaxy is constrained by the ugizk limit for 98% of the sample, and is defined by the r 17.5 limit for the other 2% of the galaxies. Functionally, these limits are u =18.50, g =17.74, i =16.94, z =16.59, and K = METHODOLOGY: K-CORRECTIONS, EVOLUTION CORRECTIONS, AND STELLAR M/L RATIOS 3.1. The Method To estimate LFs and stellar MFs using the redshift and ugrizk data for the SDSS EDR galaxies, we must estimate k-corrections and stellar M/L ratios. Furthermore, Blanton et al. (2003c) and Norberg et al. (2002) stress the need to include the effects of galaxy evolution. We estimate k-corrections, evolution corrections, and galaxy stellar M/L ratios by comparing the ugrizk galaxy fluxes with state-of-the art stellar population synthesis (SPS) models. For each galaxy, we construct a grid of stellar populations with a range of metallicities and star formation histories (SFHs) at both the real galaxy redshift and at redshift zero. We use the PÉGASE model (see Fioc & Rocca-Volmerange 1997, for a description of an earlier version of the model), choosing ten galaxy metallicities from 0.5% to 250% solar. The SFHs vary exponentially with time t: ψ = [τ 1 (1 e T0/τ ) 1 ]e t/τ, where ψ is the star formation rate (SFR), τ is the exponential e-folding time of the SFR, and T 0 is the age of the galaxy (the time since SF commenced). The term in the square brackets is simply a normalization to keep the total mass of stars formed by the present day at one solar mass. We choose a grid of 29 τ values between 0 (single burst) and (continuous), continuing through to and then to 1 Gyr (strongly increasing to the present day). Our grid covers color space relatively uniformly.

6 6 Bell et al. FIG. 4. Number counts of galaxies per square degree per magnitude as a function of 2MASS K-band Kron apparent magnitude for the whole sky with b 30 (solid line), the SDSS EDR (dot-dashed) and the 2dFGRS (dashed line). Linear (left) and logarithmic (right) scales are shown. The points with error bars show the error in the EDR galaxy number density (vertical error bars) and the magnitude range used for each bin (horizontal error bars). This grid is produced at both redshift zero and at the real galaxy redshift. The galaxy age is 12 Gyr for the redshift zero model, and is younger for the non-zero redshift model assuming h = 0.7 (in essence, we choose a formation redshift of 4). For example, this gives an age of 10.7 Gyr for a galaxy at z = 0.1. We least-squares fit the model galaxies at the real galaxy redshift to the observed galaxy colors to choose the best model galaxy template. We then estimate the evolution correction, k- correction, and present-day stellar M/L ratio by comparing the non-zero redshift model with the evolved redshift zero model. Thus, in essence, we correct for evolution by assuming that the SFH indicated by the colors of the galaxy at the observed redshift continues smoothly to the present day. To estimate stellar masses, we adopt the z = 0 model galaxy M/L ratios in each passband, assuming solar absolute magnitudes of (6.41, 5.15, 4.67, 4.56, 4.53, 3.32) in ugrizk respectively, estimated using the PÉGASE SPS model. Those wishing to convert our luminosity densities into physical units or SDSS or 2MASS-calibrated absolute magnitudes per cubic Mpc can easily use the above solar absolute magnitudes for conversion without loss of accuracy. We estimate uncertainties in k-corrections, evolutionary corrections, and stellar M/L ratio values via three methods: (i) omitting one passband at a time from the SED fit (the jackknife method; 6 fits); (ii) uniformly weighting all data points in the fit for each galaxy (1 fit); and (iii) adding random magnitude offsets with sizes corresponding to the magnitude error to all the galaxy photometry and re-doing the fits (5 times). We then compute the errors from the RMS difference between these 12 different fits to the SED and our original SED fit. Typical k-correction and evolution correction random errors derived in this way for the g-band selected sample are (0.06, 0.03, 0.03, 0.03, 0.03, 0.02) mag in ugrizk respectively. Typical random errors in stellar M/L ratio are (0.08, 0.06, 0.05, 0.04, 0.04, 0.05) dex, again in ugrizk. We adopt a diet Salpeter (1955) stellar IMF (following Bell & de Jong 2001) that has the same colors and luminosity as a normal Salpeter IMF, but with only 70% of the mass due to a lower number of faint low-mass stars. This yields stellar M/L ratios 30% lower at a given color than a Salpeter IMF. Bell & de Jong (2001) show that this IMF is maximum disk, inasmuch as IMFs richer in low-mass stars over-predict the rotation velocity of Ursa Major Cluster galaxies with K-band photometry and well-resolved HI rotation curves. This prescription thus gives the maximum possible stellar M/L ratio. Naturally, a different choice of IMF allows lower M/L ratios. For example, the popular Kennicutt (1983) or Kroupa, Tout & Gilmore (1993) IMFs have roughly 30% lower M/L ratios than this IMF, and are thus submaximal. We discuss this issue in more detail in Comparison with Other Constraints Our k-corrections and evolution corrections are quite robust. K-band k-corrections are insensitive to galaxy spectral type. In particular, we find k(z) 2.1 ± 0.3z, which is in good agreement with k(z) 2.25z from Glazebrook et al. (1995). We test the optical k-corrections by comparing with a simple powerlaw interpolation, including the effects of bandpass widening. Blanton et al. (2003a) find that this approximation is good to around 0.1 mag in all passbands, but better in riz as the spectral shapes are simpler there. We find also that our optical and NIR k-corrections are consistent with the simple powerlaw recipe to within 0.1 mag in all passbands. These offsets decrease to 0.05 mag in riz. This agreement is more than adequate, bearing in mind our mag k-correction errors. We quantify our evolution corrections by comparing the mean k+evolution correction with the mean k-correction for our g- band selected galaxy sample. The mean evolution corrections are (2.3,1.6,1.3,1.1,1.0,0.8)z in ugrizk, in the sense that galaxies are fainter at the present day, owing mostly to passive evolution. This can be compared to the Q values derived by Blanton et al. (2003c), who estimate the evolution by fitting for it explicitly in their LF estimation. They find an evolution of (4.2±0.9,2.0±0.5,1.6±0.3,1.6±0.4,0.8±0.3)z in ugriz. Therefore, we find satisfactory agreement between our color-based evolutionary corrections and direct estimates from the LF evolution by Blanton et al. (2003c), except perhaps in

7 Galaxy Luminosities and Stellar Masses 7 FIG. 5. The g r colors, r and K-band surface brightnesses, and K-band apparent magnitudes of SDSS 13 r 17.5 galaxies with (solid) and without (dashed) 2MASS matches. K-band surface brightnesses and apparent magnitudes for the 3965 galaxies without 2MASS data are estimated using the r-band derived quantities and the estimated r K color of the best-fit SED model, as described in 3. The dotted line in the upper right panel shows the surface brightness limit of SDSS, and the dotted line in the lower right panel shows the K < magnitude limit that we adopt for our LF analysis. the u-band, where the photometric and k-correction uncertainties are largest, and our assumption of smoothly-varying SFHs could easily prove inadequate. Independently, Bernardi et al. (2003a) find evolution of (1.2,0.9,0.8,0.6)z in griz for earlytype galaxies using a similar (but totally independent) technique to Blanton et al. (2003c), again within 0.05 mag of our corrections over the redshift interval of interest. Bell & de Jong (2001) demonstrate that for galaxies with relatively smooth SFHs, stellar M/L ratio and optical color should correlate quite tightly. We present a test of our stellar M/L ratio estimates in Fig. 6. Using an independent method that accounts for bursts of SF based on the strengths of the 4000Å break and the Hδ line, Kauffmann et al. (2003a) construct stellar M/L ratios for over 120,000 SDSS galaxies. In Fig. 19 of that paper, they compare their M/L ratios in g-band with the g r color, estimated at z = 0.1, and find a strong correlation. To compare to the Kauffmann et al. (2003a) correlation we estimate a color correction (g r) z= (g r), assuming a power-law k-correction. Moreover, we account for the IMF difference; our diet Salpeter IMF is 0.15 dex heavier at a given color, because of its larger number of low-mass stars, than the Kennicutt IMF that Kauffmann et al. (2003a) adopt. To within 20% random scatter our multi-color method gives results consistent with their spectral method (comparing the points with the solid line in Fig. 6). This is particularly impressive given the very different methodologies and the different stellar population models used 6. With the low scatter in the g r versus M/L g ratio correlation it is possible to predict stellar mass to within 20% using g- and r-band data alone, compared to the maximum-likelihood SED fits of up to 6 optical/nir passbands. Kauffmann et al. (2003a) find a scatter closer to 50%; this is likely due to the 6 Lee, McCall & Richer (2003) show a comparison of their color-derived two-population B-band stellar M/L ratio estimates (Lee et al. 2003) with optical color in the Appendix to that paper. They find excellent agreement with color-derived stellar M/L ratios from Bell & de Jong (2001), showing again that different methodologies yield consistent estimates of stellar M/L ratio.

8 8 Bell et al. different methods adopted to derive stellar M/L ratios by our group and Kauffmann et al. (2003a). Kauffmann et al. (2003a) use 3 aperture spectra, plus an r i color that is emission-line sensitive for dust estimation. Thus, they are sensitive to the aperture mismatch between the spectra and colors, and model mismatches between color and spectral features. In contrast, we minimize the residuals explicitly between our galaxy model colors and the observed colors; therefore, we explicitly minimize the spread in the color M/L ratio correlation with our method. Either way, it is clear that we can use SDSS color data alone, plus a redshift, to estimate the stellar mass of galaxies to between 20% and 50%, relative to the answer that one obtains using K-band data or spectra. In particular, this allows us to use the SDSS data to fill in areas of parameter space not covered as completely by 2MASS, such as blue LSB galaxies. We can independently check these M/L ratios using recent results from Bernardi et al. (2003b). They construct estimates of total M/L ratio (including the potentially non-negligible contribution of dark matter) using kinematic constraints, by multiplying the half-light radius by the velocity dispersion squared, and then dividing by half the luminosity. Assuming the Hubble Space Telescope Key Project distance scale (Freedman et al. 2001), they then compare these M/L ratios with g r color, finding a strong correlation (the lower right-hand panel of their Fig. 5). Correcting for their application of a 0.08 mag bluewards offset in g r color, we estimate their log 10 (M/L r ) (g r). Over the g r range of interest (0.3 g r 1), this is within 25% at the blue end and 5% of the red end of our maximum-disk tuned stellar population model expectation (see the Appendix for details): log 10 (M/L r ) = (g r). Furthermore, their total scatter (including contributions from observational error) is 0.15 dex, or 40% in terms of M/L ratio, in agreement with our earlier estimate of 20% to 50%. The agreement between these two totally independent methodologies, each with their own sources of systematic and random error, is astonishing; both predict roughly a factor of five change in stellar M/L ratio from the blue to the red end of the galaxy population, and both have the same absolute stellar M/L ratio scale. This agreement is another, powerful argument in favor of a color-based stellar M/L ratio of the type discussed in the Appendix, or by Bell & de Jong (2001). FIG. 6. Our g-band stellar M/L ratio estimate (from the maximumlikelihood fit to the galaxy SED) against g r color for the g-band selected sample of galaxies. The dashed line is the bi-weight least-squares fit to the data. The solid line is a rough fit to the relationship in Fig. 19 of Kauffmann et al. (2003a), accounting for the 0.15 dex offset between a Kennicutt IMF and our diet Salpeter IMF. In addition, we transform the Kauffmann et al. (2003a) z = 0.1 color to a z = 0 color using (g r) z= (g r) (see text for full details). The arrows show the average (short arrow) and maximum (long arrow) effect of dust on this relationship, as discussed in Systematic uncertainties The above prescription for estimating k-corrections, evolution corrections, and stellar M/L ratios assumes that the colors of a stellar population are driven by star formation history (SFH) and metallicity alone. What are the systematic uncertainties introduced by neglecting the effects of dust and more complex SFHs? Our k-corrections are robust, inasmuch as we simply use a physically-motivated model to interpolate between the observations (a simple power-law interpolation suffices also to roughly 0.1 mag). Furthermore, our evolutionary corrections agree with independent estimates, and since they are a relatively small correction ( 0.2 mag typically), small errors in the evolutionary correction will not substantially affect our results. However, there may be a significant uncertainty in stellar M/L ratio estimates that is not accounted for by our prescription. Overall galaxy age (i.e. the time since SF started) can change the stellar M/L ratio at a given color in a systematic sense by a small but non-negligible amount, e.g. ±0.05 dex for an age difference of ±3 Gyr. Furthermore, it is not a priori clear what effects dust may have on the stellar M/L ratios. Bell & de Jong (2001) show that, to first order, the effects of dust cancel out to within dex when estimating color-derived stellar M/L ratios. This cancellation occurs because the stellar populations and dust each predict roughly the same amount of reddening per unit fading in most passbands. However, the random uncertainties of this technique are only 20% in terms of M/L ratio (see above earlier in 3), so the second order difference between the effects of dust and stellar populations could be significant. We explore the possible effects of dust on our results in a simple way, following Tully et al. (1998). Tully et al. (1998) estimate the luminosity-dependent dust content of disk galaxies by minimizing the scatter in the color-magnitude relation (CMR) in BRIK passbands. They find that luminous galaxies suffer from a 1.7 (0.3) mag dimming in their B-band (K-band) flux when going from face-on to nearly edge-on, while faint galaxies show very little evidence for dust. We adopt a rough dimming of (1.6, 1.3, 0.3) mag from face-on to edge-on ( 80 ) in grk passbands for massive disk galaxies (masses > h 2 M and c r < 2.6). We allow this dimming to decrease linearly with logarithmic mass to zero for stellar masses below h 2 M. We assume a simple slab model, with an optical depth at arbitrary inclination of τ(i) = τ(80 )cos80 /cosi, where τ(80 ) is the optical depth at edge-on derived from the above quoted difference between edge-on and face-on. Concentrated, i.e. early-type, galaxies are assigned a factor of 3 less dust. We assume a random distribution of orientations. In Fig. 6, we show schematically the effect of the average (short arrow) and maximum possible (long arrow) dust contents, according to our admittedly ad hoc description. The arrows show the bluing of color and reduction of M/L ratio when dust is taken into account. It is clear that the effects of dust and

9 Galaxy Luminosities and Stellar Masses 9 stellar population are mostly degenerate in agreement with Bell & de Jong (2001). Nevertheless, there is a slight systematic difference between the two effects. For this dust prescription, we overestimate the average stellar M/L ratio in g-band by 0.06 dex when we fit dust-reddened colors with pure stellar populations. Another source of systematic uncertainty is from bursts of SF. Bell & de Jong (2001) find that large bursts of SF can cause an over-estimate of the true M/L ratio of the stellar population, if the stellar population is interpreted in terms of smoothly varying SFHs. We attempt to constrain the magnitude of this error for our purposes using a simple model. We choose two solar metallicity stellar populations, one with a decreasing SFR to the present day (τ = 4 Gyr), and one with constant SF (τ = Gyr). We then apply random variations in SFR over timescales of 10 8 yr, distributed in a log-normal fashion with a dispersion σ of a dex, i.e. the SFR can easily change by more than an order of magnitude from its baseline rate. We then examine the offset from the g r color and M/L g ratio correlation of these bursty models, compared to smooth SFH models. For both SF models we find that SF bursts generate a 25% scatter about the color M/L ratio relation, and a 10% offset to slightly lower M/L ratio at a given color. A full order of magnitude variation in SFH over 10 8 yr timescales is likely to be an upper limit for all but the strongest present-day star-bursting galaxies; therefore, we demonstrate that the bias we impose by assuming such simplistic SFHs is 10%. To summarize, the random uncertainties of color-based stellar M/L ratio estimation are 20%. Systematic uncertainties from galaxy age, dust, and bursts of SF are 0.1 dex, or 25%. These systematic uncertainties will not cancel out with larger galaxy samples, and will dominate, along with stellar IMF, the error budget of the stellar mass density of the Universe. 4. LUMINOSITY FUNCTIONS 4.1. LF estimation We estimate LFs using the simple and intuitive V/V max formalism of, e.g., Felten (1977). This method has the disadvantage that it is somewhat sensitive to galaxy density variations. For example, if the near part of the survey is rather overdense, where a magnitude-limited survey is most sensitive to low-luminosity galaxies, then the V/V max estimator will yield a somewhat larger number of low-luminosity galaxies than it should. Both the Step-Wise Maximum Likelihood (SWML) method of Efstathiou et al. (1988) and the parametric method of Sandage et al. (1979) are insensitive to density fluctuations of this type (although both methods are sensitive to density fluctuations when calculating the overall LF normalization). Nevertheless, both the SWML and parametric method make the assumption that the shape of the LF is independent of environment, yet there is impressive evidence against this assumption, at least in the optical (De Propris et al. 2003; Hütsi et al. 2003). In contrast, the V/V max method does not make this assumption. Furthermore, it does not make any a priori assumptions regarding the form of the LF, unlike Sandage et al. (1979) parametric method. Therefore, we use the V/V max method and note, that among others, Cole et al. (2001) find with a similar dataset that LFs derived using V/V max and SWML are identical within the errors. For an unbiased estimate of V max, we estimate the maximum distance that a galaxy of a given absolute magnitude would be visible, accounting for Galactic foreground extinction and k- and evolution-corrections, not including the early-type and K-band-to-total corrections. Due to slight inaccuracies in the k- and evolution corrections, we find a few galaxies with distances that are larger by a few observational sigma than expected, which gives a few galaxies with V/V max > 1. Not including these galaxies does not affect any of the results in this paper. Our formal error estimates include Poisson, Monte- Carlo magnitude, evolution, k-correction, and V/V max bootstrap uncertainties, plus random stellar M/L ratio errors for the MFs. There are also systematic sources of error: e.g., the 25% systematic uncertainty from dust and bursts of SF, and the 5% systematic uncertainty in absolute magnitude calibration (see, e.g., Fukugita et al. 1996). We have to first order avoided uncertainties from galaxy clustering because we have renormalized the luminosity functions to account for the EDR s 8% overdensity ( 2.4). However, we neither sample all of the EDR region, nor can we securely extrapolate to different galaxy populations (e.g., the 13 r 17.5 population as opposed to the 10 K 13.5 population), so we attach a 10% percent systematic uncertainty to the normalization φ and luminosity density j to account for clustering. We summarize the systematic error budget in Table 2. We calculate LFs using pseudo-ugrizk-limited samples (where 98% of the galaxies are limited in the passband of interest, and only 2% of the galaxies are limited by their r-band flux). We present our results in Table 3 and include some relevant comparisons from the literature. We discuss two sets of LFs in more detail in this paper: a joint r and K-band selected sample, and a joint g and r-band selected sample such that only 2% of galaxies are r-band limited, and 98% are limited by the magnitude limit in the other passband. We do not calculate the LF or stellar MF for magnitude bins with less than 5 galaxies The K-band limited sample In Fig. 7, we plot V/V max versus K-band absolute magnitude. In an unbiased sample, an average value of V/V max = 0.5 is expected, as galaxies uniformly fill the volume. For our sample, the average value is ± Excluding the 66/6282 (1%) of galaxies with V/V max > 1 gives V/V max = This indicates a slight tendency for galaxies to be in the more distant half of the sample, perhaps reflecting uncertainty in the evolution correction or small amounts of large-scale structure. Nevertheless, any bias in the sample is weak; for example, Cole et al. (2001) find V/V max 0.52 for their sample of 2MASS/2dFGRS galaxies, yet obtain excellent agreement between V/V max and SWML estimates of the LF. We plot the LFs derived using the K-band limited sample in Figs. 8 and 9. In Fig. 8, we show the K-band LF for samples using different K-band surface brightness µ K cuts. All galaxies in our sample (solid line) have µ K < 20 mag arcsec 2. We fit the Schechter (1976) function to the V/V max data points: ( ) α ( L φ(l)dl = φ L exp L ) dl L L, (1) where φ is the LF normalization, L is the characteristic luminosity at the knee of the LF where the form changes from exponential to power law, and α is the faint end slope. In common with other work (e.g., de Jong & Lacey 2000; Cross & Driver 2002), we find that a fainter surface brightness limit increases φ somewhat and substantially affects α. The K- band luminosity density for our galaxy sample is 5.77± hl Mpc 3 (formal error only). As discussed in Table 2, we

10 10 Bell et al. TABLE 2 SYSTEMATIC ERROR BUDGET Quantity Error Source Ref. (1) (2) (3) (4) Luminosity Function φ 10% Uncertainty in exact sky coverage (3%), completeness (7%), Poisson error in normalization 2.4, 4.1 (1%), and differences between behavior of the 10 < K < 13.5 sample and our EDR sample M 5% Uncertainty in absolute calibration of ugrizk system (1) 10% K only: Extrapolation to total 2.2 α 0.1? Optical: from departures from a Schechter function NIR: from strong departures from a Schechter function, and LSB galaxy incompleteness 2.4, 4.2 j 15% Optical: from φ and M uncertainty above +35% 15% NIR: from φ, M and α uncertainty above & 4.2 Stellar Mass Function M & ρ 30% Dust, bursts of SF, galaxy age, and absolute calibration uncertainty above & % 60% Stellar IMF 6.1 References. (1) Fukugita et al. (1996) Note. Column (1) describes the quantity, (2) the contribution to the systematic error budget, (3) describes the error in more detail, and (4) gives any relevant references (section number or literature citation). TABLE 3 GALAXY LUMINOSITY FUNCTION FITS Band m lim N gal V/V max z φ M 5log 10 h α j j literature Ref. (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) u ± (8) 18.13(3) 0.95(3) a g ± (5) 19.73(3) 1.03(3) a r ± (7) 20.57(3) 1.07(3) a i ± (4) 21.00(3) 1.11(3) a z ± (4) 21.34(2) 1.06(2) a K ± (7) 23.29(5) 0.77(4) , b,c Note. The passband (1), corresponding limiting magnitude (2), and number of galaxies (3). The mean V/V max is in (4), and the mean redshift in (5). Each passband LF is fit with a Schechter function and described by three parameters the normalization φ in h 3 Mpc 3 mag 1 (6), the characteristic luminosity L (7), and the faint end slope α (8). Our luminosity density estimate j (9), compared to an estimate from the literature (10), both in units of hl Mpc 3. The literature references are as follows: (a) Blanton et al. (2003c); (b) Cole et al. (2001); (c) Kochanek et al. (2001). The formal error estimates for quantities are given in parentheses. Table 2 gives a complete summary of the systematic error sources, in addition to the formal errors calculated above. The optical luminosity densities give the formal error as the positive error bar, and the influence of the correction of early-type galaxy magnitudes by 0.1 mag in SDSS as the negative error bar. The K-band error estimate includes a substantial uncertainty from 2MASS s bias against LSB galaxies. See Fig. 15 for a graphical representation of the luminosity density literature comparison.

11 Galaxy Luminosities and Stellar Masses 11 FIG. 7. V/V max versus K-band absolute magnitude. The median (thick solid line), and upper and lower quartiles (shaded area), are shown as a function of K-band absolute magnitude. The average value for the whole sample is ± 0.004, which is reasonably consistent with the expected value of 0.5 (thin solid line). FIG. 9. K-band LF split by morphological type. The solid line represents the total LF. The dotted and dashed lines represent the LF for late and earlytype galaxies, separated using c r = 2.6. The solid grey line denotes the 2MASS LF of (Kochanek et al. 2001), while the grey thick dashed and dotted lines denote the early/late types from that paper, respectively. FIG. 8. K-band V/V max LFs using different surface brightness cuts. The black solid line with data points represents the total sample (µ K < 20 mag arcsec 2 ) LF. The dotted and dashed lines represent the LF for µ K < 17 mag arcsec 2 and µ K < 18 mag arcsec 2 subsamples, which shows the LF steepening at the faint end as the surface brightness limit gets fainter. The thick grey dashed line denotes the predicted K-band LF, in the absence of selection bias (see the text for more details). The thin grey line shows the hybrid Schechter+power-law fit to the predicted K-band LF. The grey solid line with error bars denotes the 2MASS+2dFGRS LF of (Cole et al. 2001), and the thin solid line is the Schechter fit to our total LF, described in Table 3. For reference, the Schechter function fit to the total Huang et al. (2003) K-band LF is shown as a dash-dotted curve. estimate a ±15% systematic uncertainty from our extrapolation to total flux, absolute magnitude calibration, and sky coverage uncertainty. Thus, our raw estimate of K-band luminosity density is 5.8± hl Mpc 3, including the sources of random and systematic error. Earlier, we expressed concern regarding incompleteness in 2MASS for LSB galaxies. It is interesting to use the full SDSS+2MASS dataset to estimate what the K-band LF should look like, in the absence of selection bias. We use the g-band selected galaxy sample to construct a K-band LF using 9307 real K-band magnitudes and 2541 synthesized K-band magnitudes (estimated to 0.4 mag accuracy using ugriz as a constraint). This is denoted in Fig. 8 as a thick grey dashed line. The agreement is excellent at the bright end; nevertheless, the faint end slope of the predicted K-band LF is substantially steeper. A Schechter function is a poor fit to this LF owing to the kink at M K 5log 10 h 21, thus we fit a power law to the faint end between 21 M K 5log 10 h 18, which has a slope of α This bias against faint, LSB galaxies affects all 2MASS-derived estimates of not just the faint end slope, but also the total K-band luminosity density. Using this rough hybrid Schechter+power law, which has (φ,m,α) = (0.0149, 23.33, 0.88) brightwards of M K 5log 10 h = 21 and continues with power law slope 1.33 faintwards of this limit (the solid grey line in Fig. 8), we estimate that the total K- band luminosity density may be as high as hl Mpc 3. Thus, we see that 2MASS s bias against LSB galaxies may bias the faint end slope downwards, and the luminosity density estimates downwards by 25%. This conclusion is qualitatively and quantitatively consistent with a more direct assessment of light missed by 2MASS s relatively shallow exposures by Andreon (2002) 7. 7 There is another argument that suggests that 2MASS misses LSB galaxies. The luminosity-density weighted r K color (AB-Vega) of the galaxy population is 2.75±0.05 (when either r or K-band luminosity weighted), and is dominated by luminous, L galaxies. SDSS should not miss large numbers of LSB galaxies, therefore we can use the r-band luminosity density (which should be

12 12 Bell et al. How do our luminosity functions and luminosity density estimates of hl Mpc 3 compare with the literature? In Fig. 8 (solid grey line) we compare with the LF estimate of Cole et al. (2001) and find excellent agreement in both the shape of the LF and the overall normalization. Cole et al. (2001) find 5.9± hl Mpc 3 but do not account for the bias against LSB galaxies inherent in the 2MASS data. We compare our LF with Kochanek et al. (2001) in Fig. 9 (shown as a solid grey line), finding that their LF is well within our error bars. 8 They find a somewhat steeper faint end slope than we do, leading them to a slightly high luminosity density of 7 ± hl Mpc 3. Therefore, accounting for all of the sources of error, it is clear that we are consistent with both Cole et al. (2001) and Kochanek et al. (2001). Our determination has the advantage, however, that we have considerably reduced the large-scale structure uncertainty by renormalizing our LF to the whole b > 30 sky, and that we understand in detail 2MASS s bias against LSB galaxies. It will indeed be interesting to see if large, deep K-band surveys will converge towards the steeper faint end slope predicted by our analysis. In this context, comparison to the rather deeper survey of Huang et al. (2003) is particularly interesting. They have a sample of 1000 galaxies over an area of sky 50 times smaller than our area, and are unable to normalize their luminosity function to the whole b > 30 sky. Thus, their luminosity function is highly susceptible to the effects of large-scale structure. We disagree with their luminosity function (see Fig. 8) and total luminosity density of hl Mpc 3. We attribute much of this mismatch to large-scale structure. Two other effects may also contribute. First, the knee of their LF is 0.2 mag brighter than ours, which is likely caused by the uncorrected galaxy evolution in their sample. Given a median redshift of z 0.15, a 0.12 ± 0.05 mag offset is needed to correct for this evolution. Furthermore, ignoring even modest evolution can cause faint end slope over-estimation with the maximum-likelihood SWML and STY methods (Blanton et al. 2003c). Second, a Schechter fit poorly represents their dataset; an improved fit would have a sharper knee and a shallower faint-end slope, giving a substantially lower luminosity density (see Fig. 2 of Huang et al. 2003). However, it is intriguing that with deeper data they find a steep faint end slope, roughly parallel with our predicted K-band LF. We predict that further work will show that Huang et al. (2003) indeed found roughly the right faint end slope, but that they were adversely affected by large-scale structure and a small offset from ignoring evolution corrections. In Fig. 9, we show the LF split crudely by morphological type using the SDSS r-band concentration parameter ( 2.3). Recall that Strateva et al. (2001) and Hogg et al. (2002) show that most concentrated c r 2.6 galaxies are early-type (earlier than Sa), although c r could be affected by seeing (see Blanton et al. 2003b, and 2.3). We recover the classic result that the LF for early types has a flat or decreasing faint end slope and has a brighter L than late types, which have a somewhat steeper LF quite complete) plus the luminosity density-weighted r K color (which reflects the behavior of L galaxies) to estimate the K-band luminosity density. This estimate is j K hl Mpc 3, somewhat higher than the uncorrected value of hl Mpc 3. The virtue of this r (r K)-based estimate is that the faint end slope is determined by the well-constrained r-band LF, whereas in the direct approach we are forced to estimate the faint end slope from the K-band data directly. 8 They also do not account for 2MASS s bias against LSB galaxies but still find a relatively steep faint end slope to the LF, making it hard to estimate the effect of the LSB bias in their case. FIG. 10. V/V max against g-band absolute magnitude. The median (thick solid line), and upper and lower quartiles (shaded area), are shown as a function of g-band absolute magnitude. The average value for the whole sample is ± 0.003, reasonably consistent with the expected value of 0.5 (thin solid line). (e.g., Bromley et al. 1998; Blanton et al. 2001). Our type-split LFs agree qualitatively with Kochanek et al. (2001), who find a larger L for early types, although their overall LF is slightly offset from ours. We explore the role of morphological selection in The g-band limited sample One strength of our combined SDSS and 2MASS sample is that we can construct LFs in the optical ugriz passbands to accompany our NIR K-band LF. In this section, we derive a g-band limited LF. The Schechter fits for other passbands are given in Table 3 for reference, and are discussed further in 4.5. Analogous to the K-band limited sample, we show the distribution of galaxy V/V max with absolute g-band magnitude (Fig. 10), the g-band LF derived using different surface brightness limits (Fig. 11), and the g-band LF split by morphological type (Fig. 12). The g-band limited sample has V/V max = ± 0.003, which is slightly higher (by 3σ) than the purely random distribution expectation of 0.5. Nevertheless, as with the K sample, this departure is small and should not affect our results at more than the few percent level. 9 In Fig. 11, our g-band LF (solid line) compares well with that of Blanton et al. (2003c), shown as the dash-dotted line. Furthermore, our g-band luminosity density (see Table 3) is 7% larger than Blanton et al. s value. Blanton et al. did not include light lost from the low surface brightness wings of early-type galaxies in their luminosity density estimate, however. When we account for the differences 9 We note that there is substantial structure in the V/V max distribution at g > 16. This structure is likely caused by large scale structure, owing to the small V max characteristic of galaxies with faint absolute magnitudes in apparent magnitude-limited samples. Furthermore, this structure is the probable origin of the fluctuations in the g-band LF at g > 16 in Figs. 11 and 12.