Search

We investigate the characteristic radiative efficiency \epsilon, Eddington ratio \lambda, and duty cycle P_0 of high-redshift active galactic nuclei (AGN), drawing on measurements of the AGN luminosity function at z=3-6 and, especially, on recent measurements of quasar clustering at z=3-4.5 from the Sloan Digital Sky Survey. The free parameters of our models are \epsilon, \lambda, and the normalization, scatter, and redshift evolution of the relation between black hole mass \mbh and halo virial velocity V_vir. We compute the luminosity function from the implied growth of the black hole mass function and the quasar correlation length from the bias of the host halos. We test our adopted formulae for the halo mass function and halo bias against measurements from the large N-body simulation developed by the MICE collaboration. The strong clustering of AGNs observed at z=3 and, especially, at z=4 implies that massive black holes reside in rare, massive dark matter halos. Reproducing the observed luminosity function then requires high efficiency \epsilon and/or low Eddington ratio \lambda, with a lower limit (based on 2\sigma agreement with the measured z=4 correlation length) \epsilon> 0.7\lambda/(1+0.7\lambda), implying \epsilon > 0.17 for \lambda > 0.25. Successful models predict high duty cycles, P_0~0.2, 0.5, and 0.9 at z=3.1, 4.5 and 6, respectively, and they require that the fraction of halo baryons locked in the central black hole is much larger than the locally observed value. The rapid drop in the abundance of the massive and rare host halos at z>7 implies a proportionally rapid decline in the number density of luminous quasars, much stronger than simple extrapolations of the z=3-6 luminosity function would predict. (abridged), Comment: Replaced with version accepted by ApJ. More detailed analysis including black hole mergers. Results unchanged

Using a homogenous sample of 38,208 quasars with a sky coverage of $4000 {\rm deg^2}$ drawn from the SDSS Data Release Five quasar catalog, we study the dependence of quasar clustering on luminosity, virial black hole mass, quasar color, and radio loudness. At $z<2.5$, quasar clustering depends weakly on luminosity and virial black hole mass, with typical uncertainty levels $\sim 10%$ for the measured correlation lengths. These weak dependences are consistent with models in which substantial scatter between quasar luminosity, virial black hole mass and the host dark matter halo mass has diluted any clustering difference, where halo mass is assumed to be the relevant quantity that best correlates with clustering strength. However, the most luminous and most massive quasars are more strongly clustered (at the $\sim 2\sigma$ level) than the remainder of the sample, which we attribute to the rapid increase of the bias factor at the high-mass end of host halos. We do not observe a strong dependence of clustering strength on quasar colors within our sample. On the other hand, radio-loud quasars are more strongly clustered than are radio-quiet quasars matched in redshift and optical luminosity (or virial black hole mass), consistent with local observations of radio galaxies and radio-loud type 2 AGN. Thus radio-loud quasars reside in more massive and denser environments in the biased halo clustering picture. Using the Sheth et al.(2001) formula for the linear halo bias, the estimated host halo mass for radio-loud quasars is $\sim 10^{13} h^{-1}M_\odot$, compared to $\sim 2\times 10^{12} h^{-1}M_\odot$ for radio-quiet quasar hosts at $z\sim 1.5$., Comment: Updated version; accepted for publication in ApJ

We examine the accretion and merger histories of central and satellite galaxies in a smoothed particle hydrodynamics (SPH) cosmological simulation that resolves galaxies down to 7e9 M[Sun]. Most friends-of-friends halos in the simulation have a distinct central galaxy, typically two to five times more massive than the most massive satellite. As expected, satellites have systematically higher assembly redshifts than central galaxies of the same baryonic mass, and satellites in more massive halos form earlier. However, contrary to the simplest expectations, satellite galaxies continue to accrete gas and convert it to stars; the gas accretion declines steadily over a period of 0.5-1 Gyr after the satellite halo merges with a larger parent halo. Satellites in a cluster mass halo eventually begin to lose baryonic mass. Since z=1, 27% of central galaxies (above 3e10 M[Sun]) and 22% of present-day satellite galaxies have merged with a smaller system above a 1:4 mass ratio; about half of the satellite mergers occurred after the galaxy became a satellite and half before. In effect, satellite galaxies can remain "central" objects of halo substructures, with continuing accretion and mergers, making the transition in assembly histories and physical properties a gradual one. Implementing such a gradual transformation in semi-analytic models would improve their agreement with the observed colour distributions of satellite galaxies in groups and with the observed colour dependence of galaxy clustering., Comment: 14 pages, 14 figures, Submitted to MNRAS

We perform Halo Occupation Distribution (HOD) modeling to interpret small-scale and intermediate-scale clustering of 35,000 luminous early-type galaxies and their cross-correlation with a reference imaging sample of normal L* galaxies in the Sloan Digital Sky Survey. The modeling results show that most of these luminous red galaxies (LRGs) are central galaxies residing in massive halos of typical mass M ~ a few times 10^13 to 10^14 Msun/h, while a few percent of them have to be satellites within halos in order to produce the strong auto-correlations exhibited on smaller scales. The mean luminosity Lc of central LRGs increases with the host halo mass, with a rough scaling relation of Lc \propto M^0.5. The halo mass required to host on average one satellite LRG above a luminosity threshold is found to be about 10 times higher than that required to host a central LRG above the same threshold. We find that in massive halos the distribution of L* galaxies roughly follows that of the dark matter and their mean occupation number scales with halo mass as M^1.5. The HOD modeling results also allows for an intuitive understanding of the scale-dependent luminosity dependence of the cross-correlation between LRGs and L_* galaxies. Constraints on the LRG HOD provide tests to models of formation and evolution of massive galaxies, and they are also useful for cosmological parameter investigations. In one of the appendices, we provide LRG HOD parameters with dependence on cosmology inferred from modeling the two-point auto-correlation functions of LRGs., Comment: 21 pages, 10 figures, accepted for publication in ApJ

We study the formation of galaxies in a (50 Mpc/h)^3 cosmological simulation (2x288^3 particles), evolved using the entropy conserving SPH code Gadget-2. Most of the baryonic mass in galaxies of all masses is originally acquired through filamentary "cold mode" accretion of gas that was never shock heated to its halo virial temperature, confirming the key feature of our earlier results obtained with a different SPH code (Keres et al. 2005). Atmospheres of hot, virialized gas develop in halos above ~2.5e11 Msun, a transition mass that is nearly constant from z=3 to z=0. Cold accretion persists in halos above the transition mass, especially at z>=2. It dominates the growth of galaxies in low mass halos at all times, and it is the main driver of the cosmic star formation history. Satellite galaxies have accretion rates similar to central galaxies of the same baryonic mass at high redshifts, but they have less accretion than comparable central galaxies at low redshift. Relative to our earlier results, the Gadget-2 simulations predict much lower rates of "hot mode" accretion from the virialized gas component of massive halos. At z<=1, typical hot accretion rates in halos above 5e12 Msun are below 1 Msun/yr, even though our simulation does not include AGN heating or other forms of "preventive" feedback. The inner density profiles of hot gas in these halos are shallow, with long associated cooling times. The cooling recipes typically used in semi-analytic models can overestimate the accretion rates in these halos by orders of magnitude, so such models may overemphasize the role of preventive feedback in producing observed galaxy masses and colors. A fraction of the massive halos develop cuspy profiles and significant cooling rates between z=1 and z=0, a redshift trend similar to the observed trend in the frequency of cooling flow clusters., Comment: MNRAS in press, 24 pages, 11 figures. Minor changes and updated references. A copy with high resolution figures (Highly Recommended) can be found at http://www.cfa.harvard.edu/~dkeres/coldmode/ms_keres_rv1.ps.gz

The shape of the primordial matter power spectrum Plin(k) encodes critical information on cosmological parameters. At large scales, the observable galaxy power spectrum Pobs(k) is expected to follow the shape of Plin(k), but on smaller scales the effects of nonlinearity and galaxy bias make the ratio Pobs(k)/Plin(k) scale-dependent. We develop a method that can extend the dynamic range of the Plin(k) recovery by incorporating constraints on the galaxy halo occupation distribution (HOD) from the projected galaxy correlation function wp. We devise an analytic model to calculate Pobs(k) in real-space and redshift-space. Once HOD parameters are determined by matching wp for a given cosmological model, galaxy bias is completely specified, and our analytic model predicts both the shape and normalization of Pobs(k). Applying our method to SDSS main galaxy samples, we find that the real-space Pobs(k) follows the shape of the nonlinear matter power spectrum at the 1-2% level up to k=0.2 h/Mpc. When we apply our method to SDSS LRG samples, the linear bias approximation is accurate to 5% at k<0.08 h/Mpc, but the scale-dependence of LRG bias prevents the use of linear theory at k>0.08 h/Mpc. Our HOD model prediction is in good agreement with the recent SDSS LRG Pobs(k) measurements at all measured scales (k<0.2 h/Mpc), naturally explaining the shape of Pobs(k). The "Q-model" prescription is a poor description of galaxy bias for the LRG samples, and it can lead to biased cosmological parameter estimates when measurements at k>0.1 h/Mpc are included in the analysis. We quantify the potential bias and constraints on cosmological parameters that arise from applying linear theory and Q-model fitting, and we demonstrate the utility of HOD modeling of future high precision measurements of Pobs(k) on quasi-linear scales., Comment: 19 pages, 15 figures, submitted to The Astrophysical Journal

We construct evolutionary models of the populations of AGN and supermassive black holes, in which the black hole mass function grows at the rate implied by the observed luminosity function, given assumptions about the radiative efficiency and the Eddington ratio. We draw on a variety of recent X-ray and optical measurements to estimate the bolometric AGN luminosity function and compare to X-ray background data and the independent estimate of Hopkins et al. (2007) to assess remaining systematic uncertainties. The integrated AGN emissivity closely tracks the cosmic star formation history, suggesting that star formation and black hole growth are closely linked at all redshifts. Observational uncertainties in the local black hole mass function remain substantial, with estimates of the integrated black hole mass density \rho_BH spanning the range 3-5.5x10^5 Msun/Mpc^3. We find good agreement with estimates of the local mass function for a reference model where all active black holes have efficiency \eps=0.065 and L_bol/L_Edd~0.4. In this model, the duty cycle of 10^9 Msun black holes declines from 0.07 at z=3 to 0.004 at z=1 and 0.0001 at z=0. The decline is shallower for less massive black holes, a signature of "downsizing" evolution in which more massive black holes build their mass earlier. The predicted duty cycles and AGN clustering bias in this model are in reasonable accord with observational estimates. If the typical Eddington ratio declines at z<2, then the "downsizing" of black hole growth is less pronounced. Matching the integrated AGN emissivity to the local black hole mass density implies \eps=0.075 (\rho_BH/4.5x10^5 Msun/Mpc^3)^{-1} for our standard luminosity function estimate (25% higher for Hopkins et al.'s), lower than the values \eps=0.16-0.20 predicted by MHD simulations of disk accretion., Comment: replaced with version accepted by ApJ. Minor revisions

We use measurements of the projected galaxy correlation function w_p and galaxy void statistics to test whether the galaxy content of halos of fixed mass is systematically different in low density environments. We present new measurements of the void probability function (VPF) and underdensity probability function (UPF) from Data Release Four of the Sloan Digital Sky Survey, as well as new measurements of the VPF from the full data release of the Two-Degree Field Galaxy Redshift Survey. We compare these measurements to predictions calculated from models of the Halo Occupation Distribution (HOD) that are constrained to match both w_p and the space density of galaxies. The standard implementation of the HOD assumes that galaxy occupation depends on halo mass only, and is independent of local environment. For luminosity-defined samples, we find that the standard HOD prediction is a good match to the observations, and the data exclude models in which galaxy formation efficiency is reduced in low-density environments. More remarkably, we find that the void statistics of red and blue galaxies (at L ~ 0.4L_*) are perfectly predicted by standard HOD models matched to the correlation function of these samples, ruling out "assembly bias" models in which galaxy color is correlated with large-scale environment at fixed halo mass. We conclude that the luminosity and color of field galaxies are determined predominantly by the mass of the halo in which they reside and have little direct dependence on the environment in which the host halo formed. In broader terms, our results show that the sizes and emptiness of voids found in the distribution of L > 0.2L_* galaxies are in excellent agreement with the predictions of a standard cosmological model with a simple connection between galaxies and dark matter halos. (abridged), Comment: 20 emulateapj pages, 9 figures. submitted to ApJ

We present the fourth edition of the Sloan Digital Sky Survey (SDSS) Quasar Catalog. The catalog contains 77,429 objects; this is an increase of over 30,000 entries since the previous edition. The catalog consists of the objects in the SDSS Fifth Data Release that have luminosities larger than M_i = -22.0 (in a cosmology with H_0 = 70 km/s/Mpc, Omega_M = 0.3, and Omega_Lambda = 0.7) have at least one emission line with FWHM larger than 1000 km/s, or have interesting/complex absorption features, are fainter than i=15.0, and have highly reliable redshifts. The area covered by the catalog is 5740 sq. deg. The quasar redshifts range from 0.08 to 5.41, with a median value of 1.48; the catalog includes 891 quasars at redshifts greater than four, of which 36 are at redshifts greater than five. Approximately half of the catalog quasars have i < 19; nearly all have i < 21. For each object the catalog presents positions accurate to better than 0.2 arcsec. rms per coordinate, five-band (ugriz) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains basic radio, near-infrared, and X-ray emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra cover the wavelength region 3800--9200A at a spectral resolution of ~2000. The spectra can be retrieved from the public database using the information provided in the catalog. The average SDSS colors of quasars as a function of redshift, derived from the catalog entries, are presented in tabular form. Approximately 96% of the objects in the catalog were discovered by the SDSS., Comment: 37 pages, Accepted for publication in AJ

Recent observations support the idea that nuclear black holes grew by gas accretion while shining as luminous quasars at high redshift, and they establish a relation of the black hole mass with the host galaxy's spheroidal stellar system. We develop an analytic model to calculate the expected impact of mergers on the masses of black holes in massive clusters of galaxies. We use the extended Press-Schechter formalism to generate Monte Carlo merger histories of halos with a mass 10^{15} h^{-1} Msun. We assume that the black hole mass function at z=2 is similar to that inferred from observations at z=0 (since quasar activity declines markedly at z<2), and we assign black holes to the progenitor halos assuming a monotonic relation between halo mass and black hole mass. We follow the dynamical evolution of subhalos within larger halos, allowing for tidal stripping, the loss of orbital energy by dynamical friction, and random orbital perturbations in gravitational encounters with subhalos, and we assume that mergers of subhalos are followed by mergers of their central black holes. Our analytic model reproduces numerical estimates of the subhalo mass function. We find that the most massive black holes in massive clusters typically grow by a factor ~ 2 by mergers after gas accretion has stopped. In our ten realizations of 10^{15} h^{-1} Msun clusters, the highest initial (z=2) black hole masses are 5-7 x 10^9 Msun, but four of the clusters contain black holes in the range 1-1.5 x 10^{10} Msun at z=0. Satellite galaxies may host black holes whose mass is comparable to, or even greater than, that of the central galaxy. Thus, black hole mergers can significantly extend the very high end of the black hole mass function., Comment: 13 pages, 7 figures, accepted for publication in The Astrophysical Journal

We present a Fourier analysis of the clustering of galaxies in the combined Main galaxy and Luminous Red Galaxy (LRG) Sloan Digital Sky Survey (SDSS) Data Release 5 (DR5) sample. The aim of our analysis is to consider how well we can measure the cosmological matter density using the signature of the horizon at matter-radiation equality embedded in the large-scale power spectrum. The new data constrains the power spectrum on scales 100--600h^-1Mpc with significantly higher precision than previous analyses of just the SDSS Main galaxies, due to our larger sample and the inclusion of the LRGs. This improvement means that we can now reveal a discrepancy between the shape of the measured power and linear CDM models on scales 0.01

We measure the cosmological matter density by observing the positions of baryon acoustic oscillations in the clustering of galaxies in the Sloan Digital Sky Survey (SDSS). We jointly analyse the main galaxies and LRGs in the SDSS DR5 sample, using over half a million galaxies in total. The oscillations are detected with 99.74% confidence (3.0sigma assuming Gaussianity) compared to a smooth power spectrum. When combined with the observed scale of the peaks within the CMB, we find a best-fit value of Omega_m=0.256+0.029-0.024 (68% confidence interval), for a flat Lambda cosmology when marginalising over the Hubble parameter and the baryon density. This value of the matter density is derived from the locations of the baryon oscillations in the galaxy power spectrum and in the CMB, and does not include any information from the overall shape of the power spectra. This is an extremely clean cosmological measurement as the physics of the baryon acoustic oscillation production is well understood, and the positions of the oscillations are expected to be independent of systematics such as galaxy bias., Comment: 5 pages, 3 figures, updated to match version accepted by ApJ

We measure the relation between galaxy luminosity and disk circular velocity (the Tully-Fisher [TF] relation), in the g, r, i, and z-bands, for a broadly selected sample of galaxies from the Sloan Digital Sky Survey, with the goal of providing well defined observational constraints for theoretical models of galaxy formation. The input sample of 234 galaxies has a roughly flat distribution of absolute magnitudes in the range -18.5 > Mr > -22, and our only morphological selection is an axis-ratio cut b/a < 0.6 to allow accurate inclination corrections. Long-slit spectroscopy yields usable H-alpha rotation curves for 162 galaxies. Observational errors, including distance errors due to peculiar velocities, are small compared to the intrinsic scatter of the TF relation. The slope of the forward TF relation steepens from -5.5 +/- 0.2 mag/log(km/s) in the g-band to -6.6 +/- 0.2 mag/log(km/s) in the z-band. The intrinsic scatter is approximately 0.4 mag in all bands. The scatter is not dominated by rare outliers or by any particular class of galaxies, though it drops slightly, to 0.36 mag, if we restrict the sample to nearly bulgeless systems. Correlations of TF residuals with other galaxy properties are weak: bluer galaxies are significantly brighter than average in the g-band but only marginally brighter in the i-band; more concentrated galaxies are slightly fainter than average; and the TF residual is virtually independent of half-light radius, contrary to the trend expected for gravitationally dominant disks. The observed residual correlations do not account for most of the intrinsic scatter, implying that this scatter is instead driven largely by variations in the ratio of dark to luminous matter within the disk galaxy population., Comment: 23 figures, accepted by AJ, includes cosmological corrections to the data

We investigate structural properties of dark matter halos of disk galaxies in LCDM cosmology, using a well-defined sample of 81 disk-dominated galaxies from the SDSS redshift survey. We model the mass-velocity and fundamental plane relations of these galaxies, which are constructed from the galaxy stellar mass, disk scale length, and optical Halpha rotation velocity at 2.2 scale lengths. We calculate a sequence of model galaxy populations, defined by the distribution of the stellar disk-to-total mass fraction, m_d. We include the effect of adiabatic contraction of dark matter halos in response to condensation of baryons. We find that models with constant m_d underpredict the intrinsic scatter of the TF and FP relations and predict an (unobserved) strong correlation between TF residuals. Introducing a scatter of disk mass fractions and allowing the mean value m_d to scale with the stellar surface density significantly improves observational match of both the slope and intercept of the TF relation and reduces the predicted residual correlation enough to be consistent with the data. Our best-fit models with a Kroupa stellar IMF over-produce the galaxy stellar mass function and predict the virial r-band mass-to-light ratios, M_vir/L_r, systematically lower than those inferred from galaxy-galaxy weak lensing and satellite dynamics. We investigate three possible solutions to these problems: (1) ignoring the effects of adiabatic contraction, (2) adopting a ``light'' stellar IMF, or (3) considering the lower halo concentrations predicted for a low cosmological power spectrum normalization. Any of these solutions yields acceptable residual correlations and relieves most of the observational tension between the TF relation and the galaxy stellar mass function (abridged)., Comment: 20 pages, submitted to ApJ

We compare the galaxy population of an SPH simulation to those predicted by the GalICS semi-analytic model and a stripped down version without supernova and AGN feedback. The SPH simulation and the no-feedback GalICS model make similar predictions for the baryonic mass functions of galaxies and for the dependence of these mass functions on environment and redshift. The two methods also make similar predictions for the galaxy content of dark matter haloes as a function of halo mass and for the gas accretion history of galaxies. Both the SPH and no-feedback GalICS models predict a bimodal galaxy population at z=0. The "red'' sequence of gas poor, old galaxies is populated mainly by satellite systems while, contrary to observations, the central galaxies of massive haloes lie on the "blue'' star-forming sequence as a result of continuing hot gas accretion at late times. Furthermore, both models overpredict the observed baryonic mass function, especially at the high mass end. In the full GalICS model, supernova-driven outflows reduce the masses of low and intermediate mass galaxies by about a factor of two. AGN feedback suppresses gas cooling in large haloes, producing a sharp cut-off in the baryonic mass function and moving the central galaxies of these massive haloes to the red sequence. Our results imply that the observational failings of the SPH simulation and the no-feedback GalICS model are a consequence of missing input physics rather than computational inaccuracies, that truncating gas accretion by satellite galaxies automatically produces a bimodal galaxy distribution with a red sequence, but that explaining the red colours of the most massive galaxies requires a mechanism like AGN feedback that suppresses the accretion onto central galaxies in large haloes., Comment: 17 pages, 11 figures, submitted to MNRAS

The total extragalactic background radiation can be an important test of the global star formation history (SFH). Using direct observational estimates of the SFH, along with standard assumptions about the initial mass function (IMF), we calculate the total extragalactic background radiation and the observed stellar density today. We show that plausible SFHs allow a significant range in each quantity, but that their ratio is very tightly constrained. Current estimates of the stellar mass and extragalactic background are difficult to reconcile, as long as the IMF is fixed to the Salpeter slope above 1 Msun. The joint confidence interval of these two quantities only agrees with that determined from the allowed range of SFH fits at the 3-sigma level, and for our best-fit values the discrepancy is about a factor of two. Alternative energy sources that contribute to the background, such as active galactic nuclei (AGN), Population III stars, or decaying particles, appear unlikely to resolve the discrepancy. However, changes to the IMF allow plausible solutions to the background problem. The simplest is an average IMF with an increased contribution from stars around 1.5--4 Msun. A ``paunchy'' IMF of this sort could emerge as a global average if low mass star formation is suppressed in galaxies experiencing rapid starbursts. Such an IMF is consistent with observations of star-forming regions, and would help to reconcile the fossil record of star formation with the directly observed SFH., Comment: 21 pages, 7 figures, 3 tables; submitted to Monthly Notices

By comparing a collisionless cosmological N-body simulation (DM) to an SPH simulation with the same initial conditions, we investigate the correspondence between the dark matter subhalos produced by collisionless dynamics and the galaxies produced by dissipative gas dynamics in a dark matter background. When galaxies in the SPH simulation become satellites in larger groups, they retain local dark matter concentrations (SPH subhalos) whose mass is typically five times their baryonic mass. The more massive subhalos of the SPH simulation have corresponding subhalos of similar mass and position in the DM simulation; at lower masses, there is fairly good correspondence, but some DM subhalos are in different spatial positions and some suffer tidal stripping or disruption. The halo occupation statistics of DM subhalos -- the mean number of subhalos, pairs, and triples as a function of host halo mass -- are very similar to those of SPH subhalos and SPH galaxies. Gravity of the dissipative baryon component amplifies the density contrast of subhalos in the SPH simulation, making them more resistant to tidal disruption. Relative to SPH galaxies and SPH subhalos, the DM subhalo population is depleted in the densest regions of the most massive halos. The good agreement of halo occupation statistics between the DM subhalo and SPH galaxy populations leads to good agreement of their two-point correlation functions and higher order moments on large scales. The depletion of DM subhalos in dense regions depresses their clustering at R<1 Mpc/h. In these simulations, the "conversation" between dark matter and baryons is mostly one-way, with dark matter dynamics telling galaxies where to form and how to cluster, but the "back talk" of the baryons influences small scale clustering by enhancing the survival of substructure in the densest environments., Comment: 32 pages including 16 figs. Submitted to ApJ. PDF file with higher quality versions of Figs 2 and 3 available at http://www.astronomy.ohio-state.edu/~dhw/Preprints/subhalo.pdf

(Abridged) We make predictions for the pairwise velocity dispersion (PVD) of galaxies with models that are constrained to match the projected correlation function and luminosity function of galaxies in the Two-Degree Field Galaxy Redshift Survey (2dFGRS). We use these data to constrain the halo occupation distribution (HOD), then calculate the PVD by populating the halos of a high resolution N-body simulation. We examine the luminosity and scale dependence of the predicted PVD. At r<1 Mpc/h and r>3 Mpc/h, we find that the PVD decreases with increasing galaxy luminosity. This result is mostly driven by the fraction of satellite galaxies f_sat, which is well-constrained by the correlation function. We find f_sat~25% for galaxies fainter than L_star, while for brighter galaxies the satellite fraction rapidly declines, creating the decrease in the PVD with luminosity. At r=1 Mpc/h, the PVD has no dependence on luminosity because satellite galaxies dominate the statistics for all objects. Recent measurements of the PVD in Fourier space using the "dispersion model" have reported a strong decline in PVD with increasing luminosity at k=1 h/Mpc. We test this method with our HOD models, finding that there is no consistent comparison between the PVD at a given k and the true dispersion at a given value of r. This results in a luminosity dependence in k-space that is stronger than in configuration space. The luminosity dependence of the HOD results in Fourier space are consistent with those measured at k=1 h/Mpc; thus the recent measurements of the PVD are fully explainable in the context of halo occupation models. The normalization of the PVD is lower than predicted by our fiducial model, and reproducing it requires a lower value of Omega_m (~0.2 instead of 0.3), a lower value of sigma_8 (~0.7 instead of 0.9), or strong velocity bias., Comment: 15 emulateapj pages, 8 figures, submitted to ApJ

(Abridged) We investigate the power of void statistics to constrain galaxy bias and the amplitude of dark matter fluctuations. We use the halo occupation distribution (HOD) framework to describe the relation between galaxies and dark matter. After choosing HOD parameters that reproduce the mean space density n_gal and projected correlation function w_p measured for galaxy samples with M_r<-19 and M_r<-21 from the Sloan Digital Sky Survey (SDSS), we predict the void probability function (VPF) and underdensity probability function (UPF) of these samples by populating the halos of a large, high-resolution N-body simulation. If we make the conventional assumption that the HOD is independent of large scale environment at fixed halo mass, then models constrained to match n_gal and w_p predict nearly identical void statistics, independent of the scatter between halo mass and central galaxy luminosity or uncertainties in HOD parameters. Models with sigma_8=0.7 and sigma_8=0.9 also predict very similar void statistics. However, the VPF and UPF are sensitive to environmental variations of the HOD in a regime where these variations have little impact on w_p. For example, doubling the minimum host halo mass in regions with large scale (5 Mpc/h) density contrast delta<-0.65 has a readily detectable impact on void probabilities of M_r<-19 galaxies, and a similar change for delta<-0.2 alters the void probabilities of M_r<-21 galaxies at a detectable level. The VPF and UPF provide complementary information about the onset and magnitude of density- dependence in the HOD. By detecting or ruling out HOD changes in low density regions, void statistics can reduce systematic uncertainties in the cosmological constraints derived from HOD modeling, and, more importantly, reveal connections between halo formation history and galaxy properties., Comment: emulateapj, 16 pages, 13 figures

We identify galaxy groups and clusters in volume-limited samples of the SDSS redshift survey, using a redshift-space friends-of-friends algorithm. We optimize the friends-of-friends linking lengths to recover galaxy systems that occupy the same dark matter halos, using a set of mock catalogs created by populating halos of N-body simulations with galaxies. Extensive tests with these mock catalogs show that no combination of perpendicular and line-of-sight linking lengths is able to yield groups and clusters that simultaneously recover the true halo multiplicity function, projected size distribution, and velocity dispersion. We adopt a linking length combination that yields, for galaxy groups with ten or more members: a group multiplicity function that is unbiased with respect to the true halo multiplicity function; an unbiased median relation between the multiplicities of groups and their associated halos; a spurious group fraction of less than ~1%; a halo completeness of more than ~97%; the correct projected size distribution as a function of multiplicity; and a velocity dispersion distribution that is ~20% too low at all multiplicities. These results hold over a range of mock catalogs that use different input recipes of populating halos with galaxies. We apply our group-finding algorithm to the SDSS data and obtain three group and cluster catalogs for three volume-limited samples that cover 3495.1 square degrees on the sky. We correct for incompleteness caused by fiber collisions and survey edges, and obtain measurements of the group multiplicity function, with errors calculated from realistic mock catalogs. These multiplicity function measurements provide a key constraint on the relation between galaxy populations and dark matter halos., Comment: 26 emulateapj pages including 19 figures. Replaced with final ApJ version

Adopting the framework of the Halo Occupation Distribution (HOD), we investigate the ability of galaxy clustering measurements to simultaneously constrain cosmological parameters and galaxy bias. Starting with a fiducial cosmological model and galaxy HOD, we calculate spatial clustering observables on a range of length and mass scales, dynamical clustering observables that depend on galaxy peculiar velocities, and the galaxy-matter cross-correlation measurable by weak lensing. We then change one or more cosmological parameters and use chi^2-minimization to find the galaxy HOD that best reproduces the original clustering. Our parameterization of the HOD incorporates a flexible relation between galaxy occupation numbers and halo mass and allows spatial and velocity bias of galaxies within dark matter halos. Despite this flexibility, we find that changes to the HOD cannot mask substantial changes to the matter density Omega_m, the matter clustering amplitude sigma_8, or the shape parameter Gamma of the linear matter power spectrum -- cosmology and bias are not degenerate. With the conservative assumption of 10% fractional errors, the set of observables considered here can provide ~10% (1 sigma) constraints on sigma_8, Omega_m, and Gamma, using galaxy clustering data ALONE. The combination sigma_8*Omega_m^0.75 is constrained to ~5%. In combination with traditional methods that focus on large scale structure in the "perturbative" regime, HOD modeling can greatly amplify the cosmological power of galaxy redshift surveys by taking advantage of high-precision clustering measurements at small and intermediate scales (from sub-Mpc to ~20Mpc/h). At the same time, the inferred constraints on the galaxy HOD provide valuable tests of galaxy formation theory., Comment: 31 pages, 16 figures. Accepted for publication in ApJ. Minor changes to match the accepted version

Galaxy-galaxy lensing measures the mean excess surface density DS(r) around a sample of lensing galaxies. We develop a method for combining DS(r) with the galaxy correlation function xi_gg(r) to constrain Omega_m and sigma_8, going beyond the linear bias model to reach the level of accuracy demanded by current and future measurements. We adopt the halo occupation distribution (HOD) framework, and we test its applicability to this problem by examining the effects of replacing satellite galaxies in the halos of an SPH simulation with randomly selected dark matter particles from the same halos. The difference between dark matter and satellite galaxy radial profiles has a ~10% effect on DS(r) at r<1 Mpc/h. However, if radial profiles are matched, the remaining impact of individual subhalos around satellite galaxies and environmental dependence of the HOD at fixed halo mass is <5% in DS(r) for 0.12 Mpc/h, but it fails at the 30-50% level on smaller scales. The scaling of DS(r) ~ Omega_m^a(r) sigma_8^b(r) approaches the linear bias expectation a=b=1 at r>10 Mpc/h, but a(r) and b(r) vary from 0.8 to 1.6 at smaller r. We calculate a fiducial DS(r) and scaling indices a(r) and b(r) for two SDSS galaxy samples; galaxy-galaxy lensing measurements for these samples can be combined with our predictions to constrain Omega_m and sigma_8., Comment: 18 pages, 10 figures, accepted for publication in The Astrophysical Journal

We present an overview of some of the issues surrounding current models of galaxy formation, highlighting recent insights obtained from cosmological hydrodynamic simulations. Detailed examination of gas accretion processes show a hot mode of gas cooling from near the halo's virial temperature, and a previously underappreciated cold mode where gas flows in along filaments on dynamical timescales, emitting its energy in line radiation. Cold mode dominates in systems with halo masses slightly smaller than the Milky Way and below, and hence dominates the global accretion during the heydey of galaxy formation. This rapid accretion path enables prompt assembly of massive galaxies in the early universe, and results in $z\sim 4$ galaxy properties in broad agreement with observations, with the most massive galaxies being the most rapid star formers. Massive galaxies today are forming stars at a much reduced rate, a trend called downsizing. The trend of downsizing is naturally reproduced in simulations, owing to a transition from cold mode accretion in the early growth phase to slower hot mode accretion once their halos grow large. However, massive galaxies at the present epoch are still observed to have considerably redder colors than simulations suggest, suggesting that star formation is not sufficiently truncated in models by the transition to hot mode, and that another process not included in current simulations is required to suppress star formation., Comment: 10 pages, to appear in proceedings of "The Fabulous Destiny of Galaxies: Bridging Past and Present", Marseille, 2005

The ability to test the nature of dark mass-energy components in the universe through large-scale structure studies hinges on accurate predictions of sky survey expectations within a given world model. Numerical simulations predict key survey signatures with varying degrees of confidence, limited mainly by the complex astrophysics of galaxy formation. As surveys grow in size and scale, systematic uncertainties in theoretical modeling can become dominant. Dark energy studies will challenge the computational cosmology community to critically assess current techniques, develop new approaches to maximize accuracy, and establish new tools and practices to efficiently employ globally networked computing resources., Comment: 5 pages, theory white paper submitted to the Dark Energy Task Force

Nearly all proposed tests for the nature of dark energy measure some combination of four fundamental observables: the Hubble parameter H(z), the distance-redshift relation d(z), the age-redshift relation t(z), or the linear growth factor D_1(z). I discuss the sensitivity of these observables to the value and redshift history of the equation of state parameter w, emphasizing where these different observables are and are not complementary. Demonstrating time-variability of w is difficult in most cases because dark energy is dynamically insignificant at high redshift. Time-variability in which dark energy tracks the matter density at high redshift and changes to a cosmological constant at low redshift is {\it relatively} easy to detect. However, even a sharp transition of this sort at z_c=1 produces only percent-level differences in d(z) or D_1(z) over the redshift range 0.4 < z < 1.8$, relative to the closest constant-w model. Estimates of D_1(z) or H(z) at higher redshift, potentially achievable with the Ly-alpha forest, galaxy redshift surveys, and the CMB power spectrum, can add substantial leverage on such models, given precise distance constraints at z < 2. The most promising routes to obtaining sub-percent precision on dark energy observables are space-based studies of Type Ia supernovae, which measure d(z) directly, and of weak gravitational lensing, which is sensitive to d(z), D_1(z), and H(z)., Comment: 13 pages, to appear in proceedings of Wide Field Imaging From Space, New Astronomy Reviews, eds. T. McKay, A. Fruchter, and E. Linder

In a large scale view of the universe, galaxies are the basic unit of structure. A typical bright galaxy may contain 100 billion stars and span tens of thousands of light years, but the empty expanses between the galaxies are much larger still. Galaxies are not randomly distributed in space, but instead reside in groups and clusters, which are themselves arranged in an intricate lattice of filaments and walls, threaded by tunnels and pocked with bubbles. Two ambitious new surveys, the Two-Degree Field Galaxy Redshift Survey (2dFGRS) and the Sloan Digital Sky Survey (SDSS), have mapped the three-dimensional distribution of galaxies over an unprecedented range of scales. Astronomers are using these maps to learn about conditions in the early universe, the matter and energy contents of the cosmos, and the physics of galaxy formation., Comment: 4 pages, 1 figure. A brief, semi-popular review of large scale structure studies with the 2dFGRS and the SDSS

The Dark Energy Survey (DES) will use a new imaging camera on the Blanco 4-m telescope at CTIO to image 5000 square degrees of sky in the South Galactic Cap in four optical bands, and to carry out repeat imaging over a smaller area to identify and measure lightcurves of Type Ia supernovae. The main imaging area overlaps the planned Sunyaev-Zel'dovich survey of the South Pole Telescope. The idea behind DES is to use four distinct and largely independent methods to probe the properties of dark energy: baryon oscillations of the power spectrum, abundance and spatial distribution of clusters, weak gravitational lensing, and Type Ia supernovae. This white paper outlines, in broad terms, some of the theoretical issues associated with the first three of these probes (the issues for supernovae are mostly different in character), and with the general task of characterizing dark energy and distinguishing it from alternative explanations for cosmic acceleration. A companion white paper discusses the kind of numerical simulations and other theoretical tools that will be needed to address the these issues and to create mock catalogs that allow end-to-end tests of analysis procedures. Although we have been thinking about these problems in the specific context of DES, many of them are also relevant to other planned dark energy studies., Comment: 5 pages. White paper submitted to the Dark Energy Task Force

We construct merger trees for galaxies identified in a cosmological hydrodynamical simulation and use them to characterize predicted merger rates as a function of redshift, galaxy mass, and merger mass ratio. At z=0.3, we find a mean rate of 0.054 mergers per galaxy per Gyr above a 1:2 mass ratio threshold for massive galaxies (baryonic mass above 6.4E10 \Msun), but only 0.018 / Gyr for lower mass galaxies. The mass ratio distribution is proportional to R_merg^-1.2 for the massive galaxy sample, so high mass mergers dominate the total merger growth rate. The predicted rates increase rapidly with increasing redshift, and they agree reasonably well with observational estimates. A substantial fraction of galaxies do not experience any resolved mergers during the course of the simulation, and even for the high mass sample only 50% of galaxies experience a greater than 1:4 merger since z=1. Typical galaxies thus have fairly quiescent merger histories. We assign bulge-to-disk ratios to simulated galaxies by assuming that mergers above a mass ratio threshold R_major convert stellar disks into spheroids. With R_major of 1:4, we obtain a fairly good match to the observed dependence of early-type fraction on galaxy mass. However, the predicted fraction of truly bulge-dominated systems (f_bulge} > 0.8) is small, and producing a substantial population of bulge-dominated galaxies may require a mechanism that shuts off gas accretion at late times and/or additional processes (besides major mergers) for producing bulges., Comment: submitted to MNRAS, 11 pages, 11 figures, a version of the paper with a high resolution version of Figure 2 at http://www.astro.umass.edu/~ari/pub.html

Using high resolution cosmological simulations, we study hydrogen and helium gravitational cooling radiation. We focus on the HeII cooling lines, which arise from gas with a different temperature history (T_max ~ 10^5K) than HI line emitting gas. We examine whether three major atomic cooling lines, HI 1216A, HeII 1640A and HeII 304A are observable, finding that HI 1216A and HeII 1640A cooling emission at z=2-3 are potentially detectable with deep narrow band (R>100) imaging and/or spectroscopy from the ground. While the expected strength of HI 1216A cooling emission depends strongly on the treatment of the self-shielded phase of the IGM in the simulations, our predictions for the HeII 1640A line are more robust because the HeII 1640A emissivity is negligible below T~10^4.5 K and less sensitive to the UV background. Although HeII 1640A cooling emission is fainter than HI 1216A by at least a factor of 10 and, unlike HI 1216A, might not be resolved spatially with current observational facilities, it is more suitable to study gas accretion in the galaxy formation process because it is optically thin and less contaminated by the recombination lines from star-forming galaxies. The HeII 1640A line can be used to distinguish among mechanisms for powering the so-called "Lyman alpha blobs" -- including gravitational cooling radiation, photoionization by stellar populations, and starburst-driven superwinds -- because (1) HeII 1640A emission is limited to very low metallicity (log(Z/Z_sun) < -5.3) and Population III stars, and (2) the blob's kinematics are probed unambiguously through the HeII 1640A line width, which, for cooling radiation, is narrower (sigma < 400 km/s) than typical wind speeds., Comment: 15 pages, 9 figures. Accepted for publication in ApJ

We study the distribution of Eddington luminosity ratios, L_bol/L_edd, of active galactic nuclei (AGNs) discovered in the AGN and Galaxy Evolution Survey (AGES). We combine H-beta, MgII, and CIV line widths with continuum luminosities to estimate black hole (BH) masses in 407 AGNs, covering the redshift range z~0.3-4 and the bolometric luminosity range L_bol~10^45-10^47 erg/s. The sample consists of X-ray or mid-infrared (24 micron) point sources with optical magnitude R<=21.5 mag and optical emission line spectra characteristic of AGNs. For the range of luminosity and redshift probed by AGES, the distribution of estimated Eddington ratios is well described as log-normal with a peak at L_bol/L_edd ~ 1/4 and a dispersion of 0.3 dex. Since additional sources of scatter are minimal, this dispersion must account for contributions from the scatter between estimated and true BH mass and the scatter between estimated and true bolometric luminosity. Therefore, we conclude that: (1) neither of these sources of error can contribute more than \~0.3 dex rms; and (2) the true Eddington ratios of optically luminous AGNs are even more sharply peaked. Because the mass estimation errors must be smaller than ~0.3 dex, we can also investigate the distribution of Eddington ratios at fixed BH mass. We show for the first time that the distribution of Eddington ratios at fixed BH mass is peaked, and that the dearth of AGNs at a factor ~10 below Eddington is real and not an artifact of sample selection. These results provide strong evidence that supermassive BHs gain most of their mass while radiating close to the Eddington limit, and they suggest that the fueling rates in luminous AGNs are ultimately determined by BH self-regulation of the accretion flow rather than galactic scale dynamical disturbances., Comment: 34 pages including 12 figures. Incorporates referee's comments. Accepted for publication in ApJ

We measure the topology of volume-limited galaxy samples selected from a parent sample of 314,050 galaxies in the Sloan Digital Sky Survey(SDSS), which is now complete enough to describe the fully three-dimensional topology and its dependence on galaxy properties. We compare the observed genus statistic G\nu_f) to predictions for a Gaussian random field and to the genus measured for mock surveys constructed from new large-volume simulations of the LambdaCDM cosmology. In this analysis we carefully examine the dependence of the observed genus statistic on the Gaussian smoothing scale R_G from 3.5 to 11 h^-1Mpc and on the luminosity of galaxies over the range -22.50

We investigate the correlations among stellar mass (M_*), disk scale length (R_d), and rotation velocity at 2.2 disk scale lengths (V_2.2) for a sample of 81 disk-dominated galaxies (disk/total >= 0.9) selected from the SDSS. We measure V_2.2 from long-slit H-alpha rotation curves and infer M_* from galaxy i-band luminosities (L_i) and g-r colors. We find logarithmic slopes of 2.60+/-0.13 and 3.05+/-0.12 for the L_i-V_2.2 and M_*-V_2.2 relations, somewhat shallower than most previous studies, with intrinsic scatter of 0.13 dex and 0.16 dex. Our direct estimates of the total-to-stellar mass ratio within 2.2R_d, assuming a Kroupa IMF, yield a median ratio of 2.4 for M_*>10^10 Msun and 4.4 for M_*=10^9-10^10 Msun, with large scatter at a given M_* and R_d. The typical ratio of the rotation speed predicted for the stellar disk alone to the observed rotation speed at 2.2R_d is ~0.65. The distribution of R_d at fixed M_* is broad, but we find no correlation between disk size and the residual from the M_*-V_2.2 relation, implying that this relation is an approximately edge-on view of the disk galaxy fundamental plane. Independent of the assumed IMF, this result implies that stellar disks do not, on average, dominate the mass within 2.2R_d. We discuss our results in the context of infall models of disk formation in cold dark matter halos. A model with a disk-to-halo mass ratio m_d=0.05 provides a reasonable match to the R_d-M_* distribution for spin parameters \lambda ranging from ~0.04-0.08, and it yields a reasonable match to the mean M_*-V_2.2 relation. A model with m_d=0.1 predicts overly strong correlations between disk size and M_*-V_2.2 residual. Explaining the wide range of halo-to-disk mass ratios within 2.2R_d requires significant scatter in m_d values, with systematically lower m_d for galaxies with lower $M_*$., Comment: 18 pages, 2 tables, 7 figures, Accepted to ApJ, Table 1 updated, otherwise minor changes

[Abridged] We test the reliability of a method to measure the mean halo mass of Damped Ly-alpha absorbers (DLAs). The method is based on measuring the ratio of the cross-correlation between DLAs and galaxies to the auto-correlation of the galaxies themselves ($w_{\rm dg}/w_{\rm gg}$), which is (in linear theory) the ratio of their bias factor. This is shown to be true irrespective of the galaxy redshift distribution, provided that one uses the same galaxies for the two correlation functions. The method is applicable to all redshifts. Here, we focus on z=3 DLAs and we demonstrate that the method robustly constrains the mean DLA halo mass using smoothed particle hydrodynamics (SPH) cosmological simulations. If we use the bias formalism of Mo & White with the DLA and galaxy mass distributions of these simulations, we predict a bias ratio of 0.771. Direct measurement from the simulations of $w_{\rm dg}/w_{\rm gg}$ st yields a ratio of 0.73+/-0.08, in excellent agreement with that prediction. Equivalently, inverting the measured correlation ratio to infer a mean DLA halo mass yields (log. averaging, in solar units) =11.13+/-013, in excellent agreement with the true value in the simulations: 11.16. The cross- correlation method thus appears to yield a robust estimate of the average host halo mass even though the DLAs and the galaxies occupy a broad mass spectrum of halos, and massive halos contain multiple galaxies with DLAs. We show that the inferred mean DLA halo mass is independent of the galaxy sub-sample used, i.e. the cross-correlation technique is also reliable. Our results imply that the cross-correlation length between DLAs and LBGs is predicted to be, at most, 2.85 Mpc. Future observations will soon distinguish models in which DLAs are in low mass halos from those in which DLAs are in massive halos., Comment: 15 pages, 7 figures, to be published in ApJ 2005 July 20th (Full resolution of Fig.2 at http://www.mpe.mpg.de/~nbouche/papers/Xcorr/f2-orig.eps); minor changes to match the published text

We show how redshift-space distortions of the galaxy correlation function or power spectrum can constrain the matter density parameter Omega_m and the linear matter fluctuation amplitude sigma_8. We improve on previous treatments by adopting a fully non-linear description of galaxy clustering and bias, which allows us to break parameter degeneracies by combining large-scale and small- scale distortions. We consider different combinations of Omega_m and sigma_8 and find parameters of the galaxy halo occupation distribution (HOD) that yield nearly identical galaxy correlation functions in real space. We use these HOD parameters to populate the dark matter halos of large N-body simulations, from which we measure redshift-space distortions on small and large scales. We include a velocity bias parameter alpha_v that allows the velocity dispersions of satellite galaxies in halos to be systematically higher or lower than those of dark matter. Large-scale distortions are determined by the parameter combination beta = Omega_m^{0.6}/b_g, where b_g is the galaxy bias, in agreement with linear theory. However, linear theory does not accurately describe the distortions themselves on scales accessible to our simulations. We provide fitting formulas to estimate beta from the redshift-space correlation function or power spectrum, and we show that these formulas are significantly more accurate than those in the existing literature. On small scales, the ``finger-of-god'' distortions at projected separations ~0.1 Mpc/h depend on Omega_m*alpha_v^2 but are independent of sigma_8, while at intermediate separations they depend on sigma_8 as well. One can thus use redshift-space distortions over a wide range of scales to separately determine Omega_m, sigma_8, and alpha_v. (Abridged), Comment: 25 pages, submitted to Monthly Notices

We examine the dependence of the mass-to-light (M/L) ratio of large-scale structure on cosmological parameters, in models that are constrained to match observations of the projected galaxy correlation function w(rp). For a sequence of cosmological models with a fixed P(k) shape and increasing normalization \sig8, we find parameters of the galaxy halo occupation distribution (HOD) that reproduce SDSS measurements of w(rp) as a function of luminosity. Using these HOD models we calculate mean M/L ratios as a function of halo mass and populate halos of N-body simulations to compute M/L in larger scale environments, including cluster infall regions. For all cosmological models, the M/L ratio in high mass halos or high density regions is approximately independent of halo mass or smoothing scale. However, the "plateau" value of M/L depends on \sig8 as well as \Omega_m, and it represents the universal mass-to-light ratio only for models in which the galaxy correlation function is approximately unbiased, i.e., with \sig8 ~ \sig8_gal. Our results for cluster mass halos follow the trend M/L = 577(\Omega_m/0.3)(\sig8/0.9)^{1.7} h Msun/Lsun. Combined with Carlberg et al.'s (1996) mean M/L ratio of CNOC galaxy clusters, this relation implies (\sig8/0.9)(\Omega_m/0.3)^{0.6} = 0.75 +/- 0.06. M/L ratios of clusters from the SDSS and CAIRNS surveys yield similar results. This constraint is inconsistent with parameter values \Omega_m ~ 0.3, \sig8 ~ 0.9 favored by recent joint analyses of CMB measurements and other large-scale structure data. We discuss possible resolutions, none of which seems entirely satisfactory. Appendices present an improved formula for halo bias factors and an improved analytic technique for calculating the galaxy correlation function from a given cosmological model and HOD. (Abridged), Comment: Accepted to ApJ (v 630, no 2). Replaced with accepted version

We present H-alpha rotation curves of two low luminosity field galaxies with r-band absolute magnitudes M_r=-13.9^{+0.8}_{-0.5} and M_r=-14.7^{+0.3}_{-0.2} (for h=0.7; the large error bars reflect distance uncertainties). Most previously studied galaxies in this luminosity range are members of groups defined by brighter galaxies, but these two systems, selected from Blanton et al.'s (2004) sample of low luminosity galaxies in the Sloan Digital Sky Survey (SDSS), appear to have no bright companions. The measured rotation speeds at the outer extent of the H-alpha rotation curves, 34.8+/-3.8 km/s and 30.9+/-7.2 km/s, are larger than the values of 16.0^{+6.1}_{-5.4} km/s and 20.9^{+6.2}_{-5.2} km/s predicted by extrapolating the inverse Tully-Fisher relation of luminous SDSS galaxies to these faint luminosities. However, a previous HI measurement of the first galaxy shows that it has a gas mass similar to its stellar mass, and the total baryonic mass is consistent with that predicted by McGaugh et al.'s (2000) "baryonic Tully-Fisher relation." We find r-band dynamical mass-to-light ratios within the radii of the last H-alpha data points (about 1.8 disk scale lengths in each case) of 12.6^{+4.7}_{-4.5} Msun/Lsun and 4.8^{+2.5}_{-2.1} Msun/Lsun, much higher than the values ~ 1 Msun/Lsun expected for the stellar populations. The dynamical properties of these galaxies, including the rotation speeds and evidence for high gas fractions and dark matter domination within the luminous extent of the galaxy, are consistent with those of previously studied faint galaxies in nearby groups. Further studies of the SDSS sample will allow characterization of low luminosity galaxies over the full range of environments in which they reside., Comment: 15 pages, 3 figures, Submitted to the ApJL

The halo occupation distribution (HOD) describes the relation between galaxies and dark matter at the level of individual dark matter halos. The properties of galaxies residing at the centers of halos differ from those of satellite galaxies because of differences in their formation histories. Using a smoothed particle hydrodynamics (SPH) simulation and a semi-analytic (SA) galaxy formation model, we examine the separate contributions of central and satellite galaxies to the HOD, more specifically to the probability P(N|M) that a halo of virial mass M contains N galaxies of a particular class. In agreement with earlier results for dark matter subhalos, we find that the mean occupation function for galaxies above a baryonic mass threshold can be approximated by a step function for central galaxies plus a power law for satellites, and that the distribution of satellite numbers is close to Poisson at fixed halo mass. For galaxy samples defined by different baryonic mass thresholds, there is a nearly linear relation between the minimum halo mass Mmin required to host a central galaxy and the mass M1 at which an average halo hosts one satellite, with M1 ~ 14 Mmin (SPH) or M1 ~ 18 Mmin (SA). The mean occupation number of young galaxies exhibits a local minimum at M ~ 10 Mmin where halos are too massive to host a young central galaxy but not massive enough to host satellites. We show that the conditional galaxy mass function at fixed halo mass cannot be described by a Schechter function because central galaxies produce a "bump" at high masses. We suggest parameterizations for the HOD and the conditional luminosity function that can be used to model observed galaxy clustering. Many of our predictions are in good agreement with recent results inferred from clustering in the Sloan Digital Sky Survey., Comment: 38 pages, 13 figures. Figures are added to show model predictions on CLF and the conditional mass distribution of halos hosting galaxies at fixed luminosity. Accepted for publication in ApJ

We use the projected correlation function w_p(r_p) of a volume-limited subsample of the Sloan Digital Sky Survey (SDSS) main galaxy redshift catalogue to measure the halo occupation distribution (HOD) of the galaxies of the sample. Simultaneously, we allow the cosmology to vary within cosmological constraints imposed by cosmic microwave background experiments in a Lambda-CDM model. We find that combining w_p(r_p) for this sample alone with the observations by WMAP, ACBAR, CBI and VSA can provide one of the most precise techniques available to measure cosmological parameters. For a minimal flat six-parameter Lambda-CDM model with an HOD with three free parameters, we find Omega_m=0.278^{+0.027}_{-0.027}, sigma_8=0.812^{+0.028}_{-0.027}, and H_0=69.8^{+2.6}_{-2.6}km s^{-1} Mpc^{-1}; these errors are significantly smaller than from CMB alone and similar to those obtained by combining CMB with the large-scale galaxy power spectrum assuming scale-independent bias. The corresponding HOD parameters describing the minimum halo mass and the normalization and cut-off of the satellite mean occupation are M_min=(3.03^{+0.36}_{-0.36})x 10^{12} h^{-1} M_sun, M_1 = (4.58^{+0.60}_{-0.60})x 10^{13} h^{-1} M_sun, and kappa=4.44^{+0.51}_{-0.69}. When more parameters are added to the HOD model, the error bars on the HOD parameters increase because of degeneracies, but the error bars on the cosmological parameters do not increase greatly. Similar modeling for other galaxy samples could reduce the statistical errors on these results, while more thorough investigations of the cosmology dependence of nonlinear halo bias and halo mass functions are needed to eliminate remaining systematic uncertainties, which may be comparable to statistical uncertainties., Comment: v2: includes an expanded modeling discussion, a test for halo spatial bias model dependence, and a small but non-negligible contribution from sample variance

We examine the temperature history of gas accreted by forming galaxies in SPH simulations. About half the gas shock heats to roughly the virial temperature of the galaxy potential well before cooling, condensing, and forming stars, but the other half radiates its acquired gravitational energy at much lower temperatures, typically T<10^5 K, and the histogram of maximum gas temperatures is clearly bimodal. The "cold mode" of gas accretion dominates for low mass galaxies (M_baryon < 10^{10.3}Msun or M_halo < 10^{11.4}Msun), while the conventional "hot mode" dominates the growth of high mass systems. Cold accretion is often directed along filaments, allowing galaxies to efficiently draw gas from large distances, while hot accretion is quasi-spherical. The galaxy and halo mass dependence leads to redshift and environment dependence of cold and hot accretion rates, with cold mode dominating at high redshift and in low density regions today, and hot mode dominating in group and cluster environments at low redshift. Star formation rates closely track accretion rates, and we discuss the physics behind the observed environment and redshift dependence of galactic scale star formation. If we allowed hot accretion to be suppressed by conduction or AGN feedback, then the simulation predictions would change in interesting ways, perhaps resolving conflicts with the colors of ellipticals and the cutoff of the galaxy luminosity function. The transition between cold and hot accretion at M_h ~ 10^{11.4}Msun is similar to that found by Birnboim & Dekel (2003) using 1-d simulations and analytic arguments. The corresponding baryonic mass is tantalizingly close to the scale at which Kauffmann et al. (2003) find a marked shift in galaxy properties. We speculate on connections between these theoretical and observational transitions., Comment: 1 figure added, Appendix discussing SAMs added, some text changes. Matches the version accepted by MNRAS. 31 pages (MNRAS style), 21 figures,For high resolution version of the paper (highly recommended) follow http://www.astro.umass.edu/~keres/paper/ms2.ps.gz

We study the relationship between galaxy luminosity, color, and environment in a cosmological simulation of galaxy formation. We compare the predicted relationship with that found for SDSS galaxies and find that the model successfully predicts most of the qualitative features seen in the data, but also shows some interesting differences. Specifically, the simulation predicts that the local density around bright red galaxies is a strong increasing function of luminosity, but does not depend much on color at fixed luminosity. Moreover, we show that this is due to central galaxies in dark matter halos whose baryonic masses correlate strongly with halo mass. The simulation also predicts that the local density around blue galaxies is a strong increasing function of color, but does not depend much on luminosity at fixed color. We show that this is due to satellite galaxies in halos whose stellar ages correlate with halo mass. Finally, the simulation fails to predict the luminosity dependence of environment observed around low luminosity red galaxies. However, we show that this is most likely due to the simulation's limited resolution. A study of a higher resolution, smaller volume simulation suggests that this dependence is caused by the fact that all low luminosity red galaxies are satellites in massive halos, whereas intermediate luminosity red galaxies are a mixture of satellites in massive halos and central galaxies in less massive halos., Comment: 8 pages including 5 figures (emulateapj). Submitted to ApJ

We measure the power spectrum, P_F(k,z), of the transmitted flux in the Ly-alpha forest using 3035 high redshift quasar spectra from the Sloan Digital Sky Survey. This sample is almost two orders of magnitude larger than any previously available data set, yielding statistical errors of ~0.6% and ~0.005 on, respectively, the overall amplitude and logarithmic slope of P_F(k,z). This unprecedented statistical power requires a correspondingly careful analysis of the data and of possible systematic contaminations in it. For this purpose we reanalyze the raw spectra to make use of information not preserved by the standard pipeline. We investigate the details of the noise in the data, resolution of the spectrograph, sky subtraction, quasar continuum, and metal absorption. We find that background sources such as metals contribute significantly to the total power and have to be subtracted properly. We also find clear evidence for SiIII correlations with the Ly-alpha forest and suggest a simple model to account for this contribution to the power. While it is likely that our newly developed analysis technique does not eliminate all systematic errors in the P_F(k,z) measurement below the level of the statistical errors, our tests indicate that any residual systematics in the analysis are unlikely to affect the inference of cosmological parameters from P_F(k,z). These results should provide an essential ingredient for all future attempts to constrain modeling of structure formation, cosmological parameters, and theories for the origin of primordial fluctuations., Comment: 92 pages, 45 of them figures, submitted to ApJ, data available at http://feynman.princeton.edu/~pmcdonal/LyaF/sdss.html

We have developed a flexible framework for constructing physical models of quasar evolution that can incorporate a wide variety of observational constraints, such as multi-wavelength quasar luminosity functions (QLFs), estimated masses and accretion rates of active black holes, space densities of quasar host galaxies, clustering measurements, and the mass function of black holes in the local universe. In this brief contribution we focus on the observed decline in the QLF break luminosity at $z<2$, which can be explained either by a shift toward lower characteristic accretion rates at low $z$ or by preferential suppression of activity in higher mass black holes., Comment: 2 pages, 1 figure, to be published in the Proceedings of "Multiwavelength AGN Surveys", Cozumel, Dec 8 - 12, 2003

We present estimates of cosmological parameters from the application of the Karhunen-Loeve transform to the analysis of the 3D power spectrum of density fluctuations using Sloan Digital Sky Survey galaxy redshifts. We use Omega_m*h and f_b = Omega_b/Omega_m to describe the shape of the power spectrum, sigma8 for the (linearly extrapolated) normalization, and beta to parametrize linear theory redshift space distortions. On scales k < 0.16 h/Mpc, our maximum likelihood values are Omega_m*h = 0.264 +/-0.043, f_b = 0.286 +/- 0.065, sigma8 = 0.966 +/- 0.048, and beta = 0.45 +/- 0.12. When we take a prior on Omega_b from WMAP, we find Omega_m*h = 0.207 +/- 0.030, which is in excellent agreement with WMAP and 2dF. This indicates that we have reasonably measured the gross shape of the power spectrum but we have difficulty breaking the degeneracy between Omega_m*h and f_b because the baryon oscillations are not resolved in the current spectroscopic survey window function., Comment: LaTeX, 8 pages, 3 figures; figures updated, references added, minor changes to text; accepted by ApJ

We present a flexible framework for constructing physical models of quasar evolution that can incorporate a variety of observational constraints, such as multi-wavelength luminosity functions, estimated masses and accretion rates of active black holes, space densities of quasar hosts, and the local black hole mass function. We focus on the accretion rate distribution p(mdot|M,z), the probability that a black hole of mass M at redshift z accretes at a rate mdot in Eddington units. Given the radiative efficiency as a function of mdot, the quasar luminosity function (QLF) is determined by a convolution of p(mdot|M,z) with the black hole mass function n(M,z). In the absence of mergers p(mdot|M,z) also determines the full evolution of n(M,z), given a "boundary value" of n(M) at some redshift. Matching the observed decline of the QLF break luminosity at z<2 requires either a shift towards lower characteristic accretion rates or an evolving mass dependence of p(mdot) that preferentially shuts off accretion onto high mass black holes. These two scenarios make different predictions for the masses and accretion rates of active black holes. If the first mechanism dominates, then the QLF changes character between z=2 and z=0, shifting from a sequence of black hole mass towards a sequence of L/Ledd. We construct and compare five models that illustrate different assumptions about the quasar population: short and long lifetime models dominated by unobscured thin-disk accretion, a model with a high fraction of obscured quasars, and models in which mergers or ADAF accretion produce substantial black hole growth at low redshift. We discuss the observational advances that would be most valuable for distinguishing such models and for pinning down the physics that drives black hole and quasar evolution. (Abridged), Comment: submitted to ApJ, 79pp with 17 figs

{Abridged} We report the discovery of six compact, starburst galaxy candidates with redshifts 2.3 < z < 2.8 and r-band magnitudes 19.8-20.5 in the Quasar Catalog of the Sloan Digital Sky Survey First Data Release (SDSS DR1). The SDSS spectra of these objects resemble the composite spectrum of Lyman Break Galaxies (LBGs) at z ~ 3 (albeit with some differences and broader spectral lines), but the systems are 4-5 magnitudes brighter than an ``L*'' LBG and 2-3 magnitudes brighter than the most luminous objects in typical LBG spectroscopic surveys. These objects could be extremely luminous LBGs, normal LBGs amplified by gravitational lensing (like MS 1512-cB58), or a rare class of unusual BAL quasars. Under the first hypothesis, star formation rates inferred from the UV continuum luminosities are about 300-1000 M_sun yr^-1 with no correction for dust extinction. The SDSS images show no evidence of multiple imaging or foreground lensing structures, but lensing amplification cannot be ruled out with the present data. The spectra have fairly broad absorption features and prominent high-ionization absorption, but they do not have high ionization emission lines, nor do they resemble spectra of known BAL quasars. Assuming that lensing and AGN contributions are unimportant, we estimate the high end of the z ~ 2.5 galaxy luminosity function, making the crude but conservative assumption of unit detection efficiency over our observed redshift and magnitude range. The estimated comoving space density is above a Schechter function extrapolation of the LBG luminosity function from fainter magnitudes. Improved optical spectra and measurements at X-ray, IR, and sub-mm wavelengths could help determine the nature of these remarkable objects., Comment: 12 pages, 2 tables, 2 figures, submitted to ApJ Letters; replaced figure 1 and revised text in response to comments on the original version; final refereed version

The cosmic ultraviolet (UV) ionizing background impacts the formation of dwarf galaxies in the low-redshift universe (z=3) by suppressing gas infall into galactic halos with circular velocities up to v(circ)=75 km/s. Using a one-dimensional, spherically symmetric hydrodynamics code (Thoul & Weinberg 1995), we examine the effect of an ionizing background on low-mass galaxies forming at high redshifts (z>10). We find that the importance of photoionization feedback is greatly reduced, because (1) at high redshift, dwarf-galaxy sized objects can self-shield against the ionizing background, (2) collisional cooling processes at high redshift are more efficient, (3) the amplitude of the ionizing background at high redshift is lower, and (4) the ionizing radiation turns on when the perturbation that will become the dwarf galaxy has already grown to a substantial overdensity. We find that because of these reasons, gas can collect inside halos with circular velocities as low as v(circ)=10 km/s. This result has important implications for the reionization history of the universe., Comment: submitted to ApJ, 10 emulateapj pages with 4 figures included

Many open questions in X-ray astronomy are limited by the relatively small number of objects in uniform optically-identified samples, especially when rare subclasses are considered, or subsets isolated to search for evolution or correlations between wavebands. We describe initial results of a program aimed to ultimately yield 10^4 X-ray source identifications--a sample about an order of magnitude larger than earlier efforts. The technique employs X-ray data from the ROSAT All-Sky Survey (RASS), and optical imaging and spectroscopic followup from the Sloan Digital Sky Survey (SDSS). Optical objects in the SDSS catalogs are automatically cross-correlated with RASS X-ray source positions; then priorities for follow-on SDSS optical spectra of candidate counterparts are automatically assigned using an algorithm based on the known fx/fopt ratios for various classes of X-ray emitters. SDSS parameters for optical morphology, magnitude, colors, plus FIRST radio data, serve as proxies for object class. Initial application of this approach to 1400 deg^2 of sky provides a catalog of 1200 spectroscopically confirmed quasars/AGN that are probable RASS identifications. Most of the IDs are new, and only a few percent of the AGN are likely to be random superpositions. The magnitude and redshift ranges of the counterparts extend over 15

We measure the projected correlation function w_p(r_p) from the Sloan Digital Sky Survey for a flux-limited sample of 118,000 galaxies and for a volume limited subset of 22,000 galaxies with absolute magnitude M_r<-21. Both correlation functions show subtle but systematic departures from the best-fit power law, in particular a change in slope at r_p~1-2 Mpc/h. These departures are stronger for the volume-limited sample, which is restricted to relatively luminous galaxies. We show that the inflection point in w_p(r_p) can be naturally explained by contemporary models of galaxy clustering, according to which it marks the transition from a large scale regime dominated by galaxy pairs in separate dark matter halos to a small scale regime dominated by galaxy pairs in the same dark matter halo. For example, given the dark halo population predicted by an inflationary cold dark matter scenario, the projected correlation function of the volume-limited sample can be well reproduced by a model in which the mean number of M_r<-21 galaxies in a halo of mass M>M_1=4.74 X 10^{13}\msun/h is _M=(M/M_1)^{0.89}, with 75% of the galaxies residing in less massive, single-galaxy halos. This physically motivated model has the same number of free parameters as a power law, and it fits the w_p(r_p) data better, with a \chi^2/d.o.f.=0.93 compared to 6.12 (for 10 degrees of freedom, incorporating the covariance of the correlation function errors). Departures from a power-law correlation function encode information about the relation between galaxies and dark matter halos. Higher precision measurements of these departures for multiple classes of galaxies will constrain galaxy bias and provide new tests of the theory of galaxy formation., Comment: 23 pages, 5 figures. Minor changes to match accepted ApJ version

We review recent developments in the theory of the Lyman-alpha forest and their implications for the role of the forest as a test of cosmological models. Simulations predict a relatively tight correlation between the local Lya optical depth and the local gas or dark matter density. Statistical properties of the transmitted flux can constrain the amplitude and shape of the matter power spectrum at high redshift, test the assumption of Gaussian initial conditions, and probe the evolution of dark energy by measuring the Hubble parameter H(z). Simulations predict increased Lya absorption in the vicinity of galaxies, but observations show a Lya deficit within Delta_r ~ 0.5 Mpc/h (comoving). We investigate idealized models of "winds" and find that they must eliminate neutral hydrogen out to comoving radii ~1.5 Mpc/h to marginally explain the data. Winds of this magnitude suppress the flux power spectrum by \~0.1 dex but have little effect on the distribution function or threshold crossing frequency. In light of the stringent demands on winds, we consider the alternative possibility that extended Lya emission from target galaxies replaces absorbed flux, but we conclude that this explanation is unlikely. Taking full advantage of the data coming from large telescopes and from the Sloan Digital Sky Survey will require more complete understanding of the galaxy proximity effect, careful attention to continuum determination, and more accurate numerical predictions, with the goal of reaching 5-10% precision on key cosmological quantities., Comment: 13 pages, to appear in "The Emergence of Cosmic Structure," Proceedings of the 13th Annual Astrophysics Conference in Maryland, eds. S. Holt and C. Reynolds, AIP Press, 2003

Lyman alpha absorber correlations contain information about the underlying density distribution associated with a particular class of absorbers. As such, they provide an opportunity to independently measure the "bias" of the Lyman alpha forest, i.e. the relationship between HI column density and underlying dark matter density. In these proceedings we use hydrodynamic simulations to investigate whether the evolution of this bias is measurable from observable correlations. Unfortunately, the increasingly complex physics in the IGM at z<1 makes a direct measurement of the bias difficult. Nevertheless, current simulations do make predictions for HI absorber correlations that are in broad agreement with observations at both high and low redshift, thus reinforcing the bias evolution predictions given by these models., Comment: 6 pages, to appear in "The IGM/Galaxy Connection: The Distribution of Baryons at z=0", eds. J. L. Rosenberg and M. E. Putman