We develop semi-empirical models of the supermassive black hole and active galactic nucleus (AGN) populations, which incorporate the black hole growth implied by the observed AGN luminosity function assuming a radiative efficiency \epsilon, and a distribution of Eddington ratios \lambda. By generalizing these continuity-equation models to allow a distribution P(\lambda|mbh,z) we are able to draw on constraints from observationally estimated P(\lambda) distributions and active galaxy fractions while accounting for the luminosity thresholds of observational samples. We consider models with a Gaussian distribution of log \lambda, and Gaussians augmented with a power-law tail to low \lambda. Within our framework, reproducing the high observed AGN fractions at low redshift requires a characteristic Eddington ratio \lambda_c that declines at late times, and matching observed Eddington ratio distributions requires P(\lambda) that broadens at low redshift. To reproduce the observed increase of AGN fraction with black hole or galaxy mass, we also require a \lambda_c that decreases with increasing black hole mass, reducing the AGN luminosity associated with the most massive black holes. Finally, achieving a good match to the high mass end of the local black hole mass function requires an increased radiative efficiency at high black hole mass. We discuss the potential impact of black hole mergers or a \lambda-dependent bolometric correction, and we compute evolutionary predictions for black hole and galaxy specific accretion rates. Despite the flexibility of our framework, no one model provides a good fit to all the data we consider.

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)

We present Hubble Space Telescope imaging and spectroscopic observations of three Brightest Cluster Galaxies, Abell 1836-BCG, Abell 2052-BCG, and Abell 3565-BCG, obtained with the Wide Field and Planetary Camera 2, the Advanced Camera for Surveys and the Space Telescope Imaging Spectrograph. The data provide detailed information on the structure and mass profile of the stellar component, the dust optical depth, and the spatial distribution and kinematics of the ionized gas within the innermost region of each galaxy. Dynamical models, which account for the observed stellar mass profile and include the contribution of a central supermassive black hole (SBH), are constructed to reproduce the kinematics derived from the Halpha and [N II](lambda 6548,6583) emission lines. Secure SBH detection with M_bh=3.61(+0.41,-0.50)x10^9 M_sun and M_bh=1.34(+0.21,-0.19)x10^9 M_sun, respectively, are obtained for Abell 1836-BCG and Abell 3565-BCG, which show regular rotation curves and strong central velocity gradients. In the case of Abell 2052-BCG, the lack of an orderly rotational motion prevents a secure determination, although an upper limit of M_bh < 4.60x10^9 M_sun can be placed on the mass of the central SBH. These measurements represent an important step forward in the characterization of the high-mass end of the SBH mass function.

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.

We outline the physical picture of the high-redshift Ly-alpha forest that has emerged from cosmological simulations, describe statistical characteristics of the forest that can be used to test theories of structure formation, present a preliminary comparison between simulation results and measurements from Keck HIRES spectra, and discuss a recent determination of the slope and amplitude of the linear mass power spectrum P(k) at z=2.5 from moderate resolution spectra. The physical picture is simple if each QSO spectrum is viewed as a continuous non-linear map of the line-of-sight density field rather than a collection of discrete absorption lines. The distribution of flux decrements depends mainly on the amplitude and PDF (Gaussian vs. non-Gaussian) of the primordial density fluctuations. The threshold crossing frequency, analogous to the 3-d genus curve, responds to the shape and amplitude of P(k) and to the values of Omega and Lambda. Open and Lambda-CDM models agree well with the measured flux decrement distribution at smoothing lengths of 25 km/s and 100 km/s and with the threshold crossing frequency at 100 km/s. Discrepancy with the observed threshold crossing frequency at 25 km/s may reflect the combined effects of noise in the data and limited mass resolution of the simulations. The slope of the measured P(k) agrees with the predictions of inflation+CDM models. Combining the amplitude with COBE normalization imposes a constraint on these models of the form Omega h^x n^y Omega_b^z = constant. Assuming Gaussian primordial fluctuations and a power spectrum shape parameter Gamma~0.2, consistency of the measured P(k) with the observed cluster mass function at z=0 requires Omega=0.46^{+0.12}_{-0.10} for an open universe and Omega=0.34^{+0.13}_{-0.09} for a flat universe (1-sigma errors). (Shortened)