Black hole accretion versus star formation rate: theory confronts observations

We use a suite of hydrodynamical simulations of galaxy mergers to compare star formation rate (SFR) and black hole accretion rate (BHAR) for galaxies before the interaction (stochastic phase), during the `merger proper, lasting ~0.2-0.3 Gyr, and in the `remnant phase. We calculate the bi-variate distribution of SFR and BHAR and define the regions in the SFR-BHAR plane that the three phases occupy. No strong correlation between BHAR and galaxy-wide SFR is found. A possible exception are galaxies with the highest SFR and the highest BHAR. We also bin the data in the same way used in several observational studies, by either measuring the mean SFR for AGN in different luminosity bins, or the mean BHAR for galaxies in bins of SFR. We find that the apparent contradiction or SFR versus BHAR for observed samples of AGN and star forming galaxies is actually caused by binning effects. The two types of samples use different projections of the full bi-variate distribution, and the full information would lead to unambiguous interpretation. We also find that a galaxy can be classified as AGN-dominated up to 1.5 Gyr after the merger-driven starburst took place. Our study is consistent with the suggestion that most low-luminosity AGN hosts do not show morphological disturbances.

Growing black holes and galaxies: black hole accretion versus star formation rate

We present a new suite of hydrodynamical simulations and use it to study, in detail, black hole and galaxy properties. The high time, spatial and mass resolution, and realistic orbits and mass ratios, down to 1:6 and 1:10, enable us to meaningfully compare star formation rate (SFR) and BH accretion rate (BHAR) timescales, temporal behaviour and relative magnitude. We find that (i) BHAR and galaxy-wide SFR are typically temporally uncorrelated, and have different variability timescales, except during the merger proper, lasting ~0.2-0.3 Gyr. BHAR and nuclear (<100 pc) SFR are better correlated, and their variability are similar. Averaging over time, the merger phase leads typically to an increase by a factor of a few in the BHAR/SFR ratio. (ii) BHAR and nuclear SFR are intrinsically proportional, but the correlation lessens if the long-term SFR is measured. (iii) Galaxies in the remnant phase are the ones most likely to be selected as systems dominated by an active galactic nucleus (AGN), because of the long time spent in this phase. (iv) The timescale over which a given diagnostic probes the SFR has a profound impact on the recovered correlations with BHAR, and on the interpretation of observational data.

Active galactic nuclei at z ~ 1.5: I. Spectral energy distribution and accretion discs

The physics of active super massive black holes (BHs) is governed by their mass (M_BH), spin (a*) and accretion rate ($dot{M}$). This work is the first in a series of papers with the aim of testing how these parameters determine the observable attributes of active galactic nuclei (AGN). We have selected a sample in a narrow redshift range, centered on z~1.55, that covers a wide range in M_BH and $dot{M}$, and are observing them with X-shooter, covering rest wavelengths ~1200-9800 AA. The current work covers 30 such objects and focuses on the origin of the AGN spectral energy distribution (SED). After estimating M_BH and $dot{M}$ based on each observed SED, we use thin AD models and a Bayesian analysis to fit the observed SEDs in our sample. We are able to fit 22/30 of the SEDs. Out of the remaining 8 SEDs, 3 can be fit by the thin AD model by correcting the observed SED for reddening within the host galaxy and 4 can be fit by adding a disc wind to the model. In four of these 8 sources, Milky Way-type extinction, with the strong 2175AA feature, provides the best reddening correction. The distribution in spin parameter covers the entire range, from -1 to 0.998, and the most massive BHs have spin parameters greater than 0.7. This is consistent with the spin-up model of BH evolution. Altogether, these results indicate that thin ADs are indeed the main power houses of AGN, and earlier claims to the contrary are likely affected by variability and a limited observed wavelength range.

Broad-line region size and black hole mass in high-z AGN

In this contribution we briefly review the reverberation mapping technique and its results for low and intermidiate luminosity AGNs. Then we present a monitoring campaign of high-luminosity high-redshift quasars which will extend these results by two orders of magnitude, probing the broad-line region size and black hole (BH) mass of luminous AGN at redshift ~2-3.

The Hard X-ray Spectrum as a Probe for Black-Hole Growth in Radio-Quiet Active Galactic Nuclei

We study the hard-X-ray spectral properties of ten highly luminous radio-quiet (RQ) active galactic nuclei (AGNs) at z=1.3-3.2, including new XMM-Newton observations of four of these sources. We find a significant correlation between the normalized accretion rate (L/L_Edd) and the hard-X-ray photon index (Gamma) for 35 moderate-high luminosity RQ AGNs including our ten highly luminous sources. Within the limits of our sample, we show that a measurement of Gamma and L_X can provide an estimate of L/L_Edd and black-hole (BH) mass (M_BH) with a mean uncertainty of a factor of <~3 on the predicted values of these properties. This may provide a useful probe for tracing the history of BH growth in the Universe, utilizing samples of X-ray-selected AGNs for which L/L_Edd and M_BH have not yet been determined systematically. It may prove to be a useful way to probe BH growth in distant Compton-thin type 2 AGNs. We also find that the optical-X-ray spectral slope (a_ox) depends primarily on optical-UV luminosity rather than on L/L_Edd in a sample of RQ AGNs spanning five orders of magnitude in luminosity and over two orders of magnitude in L/L_Edd. We detect a significant Compton-reflection continuum in two of our highly luminous sources, and in the stacked X-ray spectrum of seven other sources with similar luminosities, we obtain a mean relative Compton reflection of R=0.9^{+0.6}_{-0.5} and an upper limit on the rest-frame equivalent width of a neutral Fe Ka line of 105 eV. We do not detect a significant steepening of the X-ray power-law spectrum below rest-frame 2 keV in any of our highly luminous sources, suggesting that a soft-excess feature, commonly observed in local AGNs, either does not depend strongly on L/L_Edd, or is not accessible at high redshifts using current X-ray detectors. [Abridged]

Black-Hole Mass and Growth Rate at High Redshift

We present new H and K bands spectroscopy of 15 high luminosity active galactic nuclei (AGNs) at redshifts 2.3-3.4 obtained on Gemini South. We combined the data with spectra of additional 29 high-luminosity sources to obtain a sample with 10^{45.2}<lambda L_{lambda}(5100A)<10^{47.3} ergs/sec and black hole (BH) mass range, using reverberation mapping relationships based on the H_beta method, of 10^{8.8}-10^{10.7} M_sun. We do not find a correlation of L/L_Edd with M_BH but find a correlation with lambda L_{lambda}(5100A) which might be due to selection effects. The L/L_Edd distribution is broad and covers the range ~0.07-1.6, similar to what is observed in lower redshift, lower luminosity AGNs. We suggest that this consistently measured and calibrated sample gives the best representation of L/L_Edd at those redshifts and note potential discrepancies with recent theoretical and observational studies. The lower accretion rates are not in accord with growth scenarios for BHs at such redshifts and the growth times of many of the sources are longer than the age of the universe at the corresponding epochs. This suggests earlier episodes of faster growth at z>~3 for those sources. The use of the C IV method gives considerably different results and a larger scatter; this method seems to be a poor M_BH and L/L_Edd estimator at very high luminosity.

Spitzer Quasar and ULIRG Evolution Study (QUEST): II. The Spectral Energy Distributions of Palomar-Green Quasars

This is the second paper studying the QSOs in the spitzer QUEST sample. Previously we presented new PAH measurements and argued that most of the observed far infrared (FIR) radiation is due to star-forming activity. Here we present spectral energy distributions (SEDs) by supplementing our data with optical, NIR and FIR observations. We define two sub-groups of ``weak FIR and ``strong FIR QSOs, and a third group of FIR non-detections. Assuming a starburst origin for the FIR, we obtain ``intrinsic AGN SEDs by subtracting a starburst template from the mean SEDs. The resulting SEDs are remarkably similar for all groups. They show three distinct peaks corresponding to two silicate emission features and a 3mic bump that we interpret as the signature of the hottest AGN dust. They also display drops beyond 20mic that we interpret as the signature of the minimum temperature (about 200K) dust. This component must be optically thin to explain the silicate emission and the slope of the long wavelength continuum. We discuss the merits of an alternative model where most of the FIR emission is due to AGN heating. Such models are unlikely to explain the properties of our QSOs but they cannot be ruled out for more luminous objects. We also find correlations between the luminosity at 5100A and two infrared starburst indicators: L(60mic) and L(PAH 7.7mic). The correlation of L(5100A) with L(60mic) can be used to measure the relative growth rates and lifetimes of the black hole and the new stars.