Ross J. McLure

Quasars, their host galaxies and their central black holes

We present the final results from our deep Hubble Space Telescope (HST) imaging study of the host galaxies of radio-quiet quasars (RQQs), radio-loud quasars (RLQs) and radio galaxies (RGs). We describe and analyse new Wide Field & Planetary Camera 2 (WFPC2) R-band observations for 14 objects, which when combined with the first tranche of HST imaging reported in McLure et al., provide a complete and consistent set of deep, red, line-free images for statistically matched samples of 13 RQQs, 10 RLQs and 10 RGs in the redshift band 0.1 < z < 0.25. We also report the results of new deep VLA imaging that has yielded a 5-GHz detection of all but one of the 33 active galactic nuclei (AGN) in our sample. Careful modelling of our images, aided by a high dynamic-range point spread function, has allowed us to determine accurately the morphology, luminosity, scalelength and axial ratio of every host galaxy in our sample. Armed with this information we have undertaken a detailed comparison of the properties of the hosts of these three types of powerful AGN, both internally and with the galaxy population in general. We find that spheroidal hosts become more prevalent with increasing nuclear luminosity such that, for nuclear luminosities MV < −23.5, the hosts of both radio-loud and radio-quiet AGN are virtually all massive ellipticals. Moreover, we demonstrate that the basic properties of these hosts are indistinguishable from those of quiescent, evolved, low-redshift ellipticals of comparable mass. This result rules out the possibility that radio-loudness is determined by host-galaxy morphology, and also sets severe constraints on evolutionary schemes that attempt to link low-z ultraluminous infrared galaxies with RQQs. Instead, we show that our results are as expected given the relationship between black hole and spheroid mass established for nearby galaxies, and apply this relation to estimate the mass of the black hole in each object. The results agree remarkably well with completely independent estimates based on nuclear emission-line widths; all the quasars in our sample have Mbh > 5 × 108 M⊙, while the radio-loud objects are confined to Mbh > 109 M⊙. This apparent mass-threshold difference, which provides a natural explanation for why RQQs outnumber RLQs by a factor of 10, appears to reflect the existence of a minimum and a maximum level of black hole radio output, which is a strong function of black hole mass (∝M2−2.5bh). Finally, we use our results to estimate the fraction of massive spheroids/black holes that produce quasar-level activity. This fraction is ≃0.1 per cent at the present day, rising to >10 per cent at z≃ 2–3.

The cosmological evolution of quasar black hole masses

Virial black hole mass estimates are presented for 12 698 quasars in the redshift interval 0.1 ≤ z ≤ 2.1, based on modelling of spectra from the Sloan Digital Sky Survey (SDSS) first data release. The black hole masses of the SDSS quasars are found to lie between ≃10 7 M ○. and an upper limit of ≃3 x 10 9 M ○. , entirely consistent with the largest black hole masses found to date in the local Universe. The estimated Eddington ratios of the broad-line quasars (full width at half-maximum ≥ 2000 km s -1 ) show a clear upper boundary at L bol /L Edd ≃ 1, suggesting that the Eddington luminosity is still a relevant physical limit to the accretion rate of luminous broad-line quasars at z ≤ 2. By combining the black hole mass distribution of the SDSS quasars with the two degree field (2dF) quasar luminosity function, the number density of active black holes at z ≃ 2 is estimated as a function of mass. In addition, we independently estimate the local black hole mass function for early-type galaxies using the M bh -σ and M bh -L bulge correlations. Based on the SDSS velocity dispersion function and the Two Micron All Sky Survey (2MASS) K-band luminosity function, both estimates are found to be consistent at the high-mass end (M bh ≥ 10 8 M ○. ). By comparing the estimated number density of active black holes at z ≃ 2 with the local mass density of dormant black holes, we set lower limits on the quasar lifetimes and find that the majority of black holes with mass ≥10 8.5 M ○. are in place by ≃2.

NEW CONSTRAINTS ON COSMIC REIONIZATION FROM THE 2012 HUBBLE ULTRA DEEP FIELD CAMPAIGN

Understanding cosmic reionization requires the identification and characterization of early sources of hydrogen-ionizing photons. The 2012 Hubble Ultra Deep Field (UDF12) campaign has acquired the deepest infrared images with the Wide Field Camera 3 aboard Hubble Space Telescope and, for the first time, systematically explored the galaxy population deep into the era when cosmic microwave background (CMB) data indicate reionization was underway. The UDF12 campaign thus provides the best constraints to date on the abundance, luminosity distribution, and spectral properties of early star-forming galaxies. We synthesize the new UDF12 results with the most recent constraints from CMB observations to infer redshift-dependent ultraviolet (UV) luminosity densities, reionization histories, and electron scattering optical depth evolution consistent with the available data. Under reasonable assumptions about the escape fraction of hydrogen-ionizing photons and the intergalactic medium clumping factor, we find that to fully reionize the universe by redshift z ~ 6 the population of star-forming galaxies at redshifts z ~ 7-9 likely must extend in luminosity below the UDF12 limits to absolute UV magnitudes of M UV ~ –13 or fainter. Moreover, low levels of star formation extending to redshifts z ~ 15-25, as suggested by the normal UV colors of z ≃ 7-8 galaxies and the smooth decline in abundance with redshift observed by UDF12 to z ≃ 10, are additionally likely required to reproduce the optical depth to electron scattering inferred from CMB observations.

Measuring the black hole masses of high-redshift quasars

A new technique is presented for determining the black-hole masses of high-redshift quasars from optical spectroscopy. The new method utilizes the full-width half maximum (FWHM) of the low-ionization Mgii emission line and the correlation between broad-line region (BLR) radius and continuum luminosity at 3000u Using archival UV spectra it is found that the correlation between BLR radius and 3000u luminosity is tighter than the established correlation with 5100u luminosity. Furthermore, it is found that the correlation between BLR radius and 3000u continuum luminosity is consistent with a relation of the form RBLR / �L 0.5 � , as expected for a constant ionization parameter. Using a sample of objects with broad-line radii determined from

Keck Spectroscopy of 3 < z < 7 Faint Lyman Break Galaxies: The Importance of Nebular Emission in Understanding the Specific Star Formation Rate and Stellar Mass Density

The physical properties inferred from the spectral energy distributions (SEDs) of z > 3 galaxies have been influential in shaping our understanding of early galaxy formation and the role galaxies may play in cosmic reionization. Of particular importance is the stellar mass density at early times, which represents the integral of earlier star formation. An important puzzle arising from the measurements so far reported is that the specific star formation rates (sSFRs) evolve far less rapidly than expected in most theoretical models. Yet the observations underpinning these results remain very uncertain, owing in part to the possible contamination of rest-optical broadband light from strong nebular emission lines. To quantify the contribution of nebular emission to broadband fluxes, we investigate the SEDs of 92 spectroscopically confirmed galaxies in the redshift range 3.8 4 than previously thought, supporting up to a 5× increase between z ≃ 2 and 7. Such a trend is much closer to theoretical expectations. Given our findings, we discuss the prospects for verifying quantitatively the nebular emission line strengths prior to the launch of the James Webb Space Telescope.

THE CANADA-FRANCE HIGH-z QUASAR SURVEY: NINE NEW QUASARS AND THE LUMINOSITY FUNCTION AT REDSHIFT 6

We present discovery imaging and spectroscopy for nine new z ∼ 6 quasars found in the Canada–France High-z Quasar Survey (CFHQS) bringing the total number of CFHQS quasars to 19. By combining the CFHQS with the more luminous Sloan Digital Sky Survey sample, we are able to derive the quasar luminosity function from a sample of 40 quasars at redshifts 5.74 <z< 6.42. Our binned luminosity function shows a slightly lower normalization and flatter slope than found in previous work. The binned data also suggest a break in the luminosity function at M1450 ≈− 25. A double power-law maximum likelihood fit to the data is consistent with the binned results. The luminosity function is strongly constrained (1σ uncertainty < 0.1 dex) over the range −27.5 <M 1450 < −24.7. The best-fit parameters are Φ(M ∗ 1450 ) = 1.14 × 10 −8 Mpc −3 mag −1 , break magnitude M ∗ 1450 =− 25.13, and bright end slope β =− 2.81. However, the covariance between β and M ∗ 1450 prevents strong constraints being placed on either parameter. For a break magnitude in the range −26 <M ∗ 1450 < −24, we find −3.8 <β <−2.3 at 95% confidence. We calculate the z = 6 quasar intergalactic ionizing flux and show it is between 20 and 100 times lower than that necessary for reionization. Finally, we use the luminosity function to predict how many higher redshift quasars may be discovered in future near-IR imaging surveys.

A new multifield determination of the galaxy luminosity function at z = 7-9 incorporating the 2012 Hubble Ultra-Deep Field imaging

We present a new determination of the ultraviolet (UV) galaxy luminosity function (LF) at redshift z ≃ 7 and 8, and a first estimate at z ≃ 9. An accurate determination of the form and evolution of the galaxy LF during this era is of key importance for improving our knowledge of the earliest phases of galaxy evolution and the process of cosmic reionization. Our analysis exploits to the full the new, deepest Wide Field Camera 3/infrared imaging from our Hubble Space Telescope (HST) Ultra-Deep Field 2012 (UDF12) campaign, with dynamic range provided by including a new and consistent analysis of all appropriate, shallower/wider area HST survey data. Our new measurement of the evolving LF at z ≃ 7 to 8 is based on a final catalogue of ≃600 galaxies, and involves a step-wise maximum-likelihood determination based on the photometric redshift probability distribution for each object; this approach makes full use of the 11-band imaging now available in the Hubble Ultra-Deep Field (HUDF), including the new UDF12 F140W data, and the latest Spitzer IRAC imaging. The final result is a determination of the z ≃ 7 LF extending down to UV absolute magnitudes M_1500 = −16.75 (AB mag) and the z ≃ 8 LF down to M_1500 = −17.00. Fitting a Schechter function, we find M*_1500 = −19.90^(+0.23)_(−0.28), log ϕ* = −2.96^(+0.18)_(−0.23) and a faint-end slope α = −1.90^(+0.14)_(−0.15) at z ≃ 7, and M*_1500 = −20.12^(+0.37)_(−0.48), log ϕ* = −3.35^(+0.28)_(−0.47) and α = −2.02^(+0.22)_(-0.23) at z ≃ 8. These results strengthen previous suggestions that the evolution at z > 7 appears more akin to ‘density evolution’ than the apparent ‘luminosity evolution’ seen at z ≃ 5 − 7. We also provide the first meaningful information on the LF at z ≃ 9, explore alternative extrapolations to higher redshifts, and consider the implications for the early evolution of UV luminosity density. Finally, we provide catalogues (including derived z_phot, M_1500 and photometry) for the most robust z ∼ 6.5-11.9 galaxies used in this analysis. We briefly discuss our results in the context of earlier work and the results derived from an independent analysis of the UDF12 data based on colour–colour selection.

The Evolution of the Galaxy Rest-Frame Ultraviolet Luminosity Function Over the First Two Billion Years

We present a robust measurement and analysis of the rest-frame ultraviolet (UV) luminosity functions at z = 4-8. We use deep Hubble Space Telescope imaging over the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey/GOODS fields, the Hubble Ultra Deep Field, and the Hubble Frontier Field deep parallel observations near the Abell 2744 and MACS J0416.1-2403 clusters. The combination of these surveys provides an effective volume of 0.6-1.2 x 10(6) Mpc(3) over this epoch, allowing us to perform a robust search for faint (M-UV = -18) and bright (Muv \textless -21) high-redshift galaxies. We select candidate galaxies using a well-tested photometric redshift technique with careful screening of contaminants, finding a sample of 7446 candidate galaxies at 3.5 \textless z \textless 8.5, with \textgreater1000 galaxies at z approximate to 6-8. We measure both a stepwise luminosity function for candidate galaxies in our redshift samples, and a Schechter function, using a Markov Chain Monte Carlo analysis to measure robust uncertainties. At the faint end, our UV luminosity functions agree with previous studies, yet we find a higher abundance of UV-bright candidate galaxies at z \textgreater= 6. Our best-fit value of the characteristic magnitude MN is consistent with -21 at z \textgreater=, 5, which is different than that inferred based on previous trends at lower redshift, and brighter at similar to 2 sigma significance than previous measures at z = 6 and 7. At z = 8, a single power law provides an equally good fit to the UV luminosity function, while at z = 6 and 7 an exponential cutoff at the bright end is moderately preferred. We compare our luminosity functions to semi-analytical models, and find that the lack of evolution in M-UV(*) is consistent with models where the impact of dust attenuation on the bright end of the luminosity function decreases at higher redshift, although a decreasing impact of feedback may also be possible. We measure the evolution of the cosmic star-formation rate (SFR) density by integrating our observed luminosity functions to M-UV = -17, correcting for dust attenuation, and find that the SFR density declines proportionally to (1 +z)(-4.3 +/- 0 5) at z \textgreater 4, which is consistent with observations at z \textgreater= 9. Our observed luminosity functions are consistent with a reionization history that starts at z greater than or similar to 10, completes at z \textgreater 6, and reaches a midpoint (xH = 0.5) at 6.7 \textless z \textless9.4. Finally, using a constant cumulative number density selection and an empirically derived rising star-formation history, our observations predict that the abundance of bright z = 9 galaxies is likely higher than previous constraints, although consistent with recent estimates of bright z similar to 10 galaxies.

KECK SPECTROSCOPY OF FAINT 3 < z < 8 LYMAN BREAK GALAXIES: EVIDENCE FOR A DECLINING FRACTION OF EMISSION LINE SOURCES IN THE REDSHIFT RANGE 6 < z < 8

Using deep Keck spectroscopy of Lyman break galaxies selected from infrared imaging data taken with the Wide Field Camera 3 on board the Hubble Space Telescope, we present new evidence for a reversal in the redshift-dependent fraction of star-forming galaxies with detectable Lyman alpha (Lyα) emission in the redshift range 6.3 < z < 8.8. Our earlier surveys with the DEIMOS spectrograph demonstrated a significant increase with redshift in the fraction of line emitting galaxies over the interval 4 < z < 6, particularly for intrinsically faint systems which dominate the luminosity density. Using the longer wavelength sensitivities of Low Resolution Imaging Spectrometer and NIRSPEC, we have targeted 19 Lyman break galaxies selected using recent WFC3/IR data whose photometric redshifts are in the range 6.3 < z < 8.8 and which span a wide range of intrinsic luminosities. Our spectroscopic exposures typically reach a 5σ sensitivity of <50 A for the rest-frame equivalent width (EW) of Lyα emission. Despite the high fraction of emitters seen only a few hundred million years later, we find only two convincing and one possible line emitter in our more distant sample. Combining with published data on a further seven sources obtained using FORS2 on the ESO Very Large Telescope, and assuming continuity in the trends found at lower redshift, we discuss the significance of this apparent reversal in the redshift-dependent Lyα fraction in the context of our range in continuum luminosity. Assuming all the targeted sources are at their photometric redshift and our assumptions about the Lyα EW distribution are correct, we would expect to find so few emitters in less than 1% of the realizations drawn from our lower redshift samples. Our new results provide further support for the suggestion that, at the redshifts now being probed spectroscopically, we are entering the era where the intergalactic medium is partially neutral. With the arrival of more sensitive multi-slit infrared spectrographs, the prospects for improving the statistical validity of this result are promising.

The black hole masses of Seyfert galaxies and quasars

The central black hole masses of a sample of 30 luminous quasars are estimated using Hβ full width half-maximum (FWHM) measurements from a combination of new and previously published nuclear spectra. The quasar black hole mass estimates are combined with reverberation-mapping measurements for a sample of Seyfert galaxies in order to study active galatic nucleus (AGN) black hole masses over a wide range in nuclear luminosity. The link between bulge luminosity and black hole mass is investigated using two-dimensional disc/bulge decompositions of the host galaxy images, the vast majority of which are high-resolution Hubble Space Telescope (HST) observations. It is found that black hole mass and bulge luminosity are well correlated and follow a relation consistent with that expected if black hole and bulge mass are directly proportional. Contrary to the recent results of Wandel, no evidence that Seyfert galaxies follow a different relation to quasars is found. However, the black hole mass distributions of the radio-loud and radio-quiet quasar subsamples are found to be significantly different, with the median black hole mass of the radio-loud quasars a factor of three larger than their radio-quiet counterparts. Finally, utilizing the elliptical galaxy fundamental plane to provide stellar velocity dispersion estimates, a comparison is performed between the virial Hβ black hole mass estimates and those of the correlations of Gebhardt et al. and Merritt & Ferrarese. With the disc geometry of the broad-line region adopted in this paper, the virial Hβ black hole masses indicate that the correct normalization of the black hole versus bulge mass relation is , while the standard assumption of purely random broad-line velocities leads to . The normalization of provided by the disc model is in remarkably good agreement with that inferred for our quasar sample using the (completely independent) correlations.