We present an evolutionary model for starbursts, quasars, and spheroidal galaxies in which mergers between gas-rich galaxies drive nuclear inflows of gas, producing starbursts and feeding the buried growth of supermassive black holes (BHs) until feedback expels gas and renders a briefly visible optical quasar. The quasar lifetime and obscuring column density depend on both the instantaneous and peak quasar luminosity, and we determine this dependence using a large set of galaxy merger simulations varying galaxy properties, orbital geometry, and gas physics. We use these fits to deconvolve observed quasar luminosity functions and obtain the evolution of the formation rate of quasars with peak luminosity, (Lpeak, z). Quasars spend extended periods at luminosities well below peak, so (Lpeak) has a maximum corresponding to the break in the observed luminosity function. From (Lpeak) and our simulations, we obtain self-consistent hard and soft X-ray and optical luminosity functions and predict many observables at multiple redshifts, including column density distributions of optical and X-ray samples, the luminosity function of broad-line quasars in X-ray samples and broad-line fraction versus luminosity, active BH mass functions, the distribution of Eddington ratios, the mass function of relic BHs and total BH mass density, and the cosmic X-ray background. In every case, our predictions agree well with observed estimates, without invoking ad hoc assumptions about source properties or distributions. We provide a library of Monte Carlo realizations of our models for comparison with observations.

/pdf/a-unified-merger-driven-model-of-the-origin-of-starbursts-1n8df6ysis.pdf

A Unified, Merger-driven Model of the Origin of Starbursts, Quasars, the Cosmic X-Ray Background, Supermassive Black Holes, and Galaxy Spheroids

We combine a large set of quasar luminosity function (QLF) measurements from the rest-frame optical, soft and hard X-ray, and near- and mid-IR bands to determine the bolometric QLF in the redshift interval z = 0-6. Accounting for the observed distributions of quasar column densities and variation of SED shapes, as well as their dependence on luminosity, makes it possible to integrate the observations in a reliable manner and provides a baseline in redshift and luminosity larger than that of any individual survey. We infer the QLF break luminosity and faint-end slope out to z ~ 4.5 and confirm at high significance (10 σ) previous claims of a flattening in both the faint- and bright-end slopes with redshift. With the best-fit estimates of the column density distribution and quasar SED, which both depend on luminosity, a single bolometric QLF self-consistently reproduces the observed QLFs in all bands and at all redshifts for which we compile measurements. Ignoring this luminosity dependence does not yield a self-consistent bolometric QLF and there is no evidence for any additional dependence on redshift. We calculate the expected relic black hole mass function and mass density, cosmic X-ray background, and ionization rate as a function of redshift and find that they are consistent with existing measurements. The peak in the total quasar luminosity density is well constrained at z = 2.15 ± 0.05. We provide a number of fitting functions to the bolometric QLF and its manifestations in various bands, as well as a script to return the QLF at arbitrary frequency and redshift from these fits.

https://iopscience.iop.org/article/10.1086/509629/pdf

An Observational Determination of the Bolometric Quasar Luminosity Function

We present a detailed and self-consistent modeling of the cosmic X-ray background (XRB) based on the most up-to-date X-ray luminosity functions (XLF) and evolution of Active Galactic Nuclei (AGN). The large body of observational results collected by soft (0.5-2 keV) and hard (2-10 keV) X-ray surveys are used to constrain at best the properties of the Compton-thin AGN population and its contribution to the XRB emission. The number ratio R between moderately obscured (Compton-thin) AGN and unobscured AGN is fixed by the comparison between the soft and hard XLFs, which suggests that R decreases from 4 at low luminosities to 1 at high luminosities. From the same comparison there is no clear evidence of an evolution of the obscured AGN fraction with redshift. The distribution of the absorbing column densities in obscured AGN is determined by matching the soft and hard source counts. A distribution rising towards larger column densities is able to reproduce the soft and hard AGN counts over about 6 dex in flux. The model also reproduces with excellent accuracy the fraction of obscured objects in AGN samples selected at different X-ray fluxes. The integrated emission of the Compton-thin AGN population is found to underestimate the XRB flux at about 30 keV, calling for an additional population of extremely obscured (Compton-thick) AGN. Since the number of Compton-thick sources required to fit the 30 keV XRB emission strongly depends on the spectral templates assumed for unobscured and moderately obscured AGN, we explored the effects of varying the spectral templates. In particular, in addition to the column density distribution, we also considered a distribution in the intrinsic powerlaw spectral indices of variable width. In our baseline model a Gaussian distribution of photon indices with mean (F) = 1.9 and dispersion or = 0.2 is assumed. This increases the contribution of the Compton-thin AGN population to the 30 keV XRB intensity by ∼30% with respect to the case of null dispersion (i.e. a single primary AGN powerlaw with I= 1.9) but is not sufficient to match the 30 keV XRB emission. Therefore a population of heavily obscured -Compton-thick- AGN, as large as that of moderately obscured AGN, is required to fit the residual background emission. Remarkably, the fractions of Compton-thick AGN observed in the Chandra Deep Field South and in the first INTEGRAL and Swift catalogs of AGN selected above 10 keV are in excellent agreement with the model predictions.

/pdf/the-synthesis-of-the-cosmic-x-ray-background-in-the-chandra-318u61qsjy.pdf

The synthesis of the cosmic X-ray background in the Chandra and XMM-Newton era

What drives the growth of black holes

We present partial-correlation analyses that examine the strengths of the relationships between l2500 A, l2 keV, αOX, and redshift for optically selected active galactic nuclei (AGNs). We extend the work of Strateva and coworkers, which analyzed optically selected AGNs from the Sloan Digital Sky Survey (SDSS), by including 52 moderate-luminosity, optically selected AGNs from the COMBO-17 survey with corresponding deep (≈250 ks to 1 Ms) X-ray observations from the Extended Chandra Deep Field-South. The COMBO-17 survey extends ~3 mag deeper than the SDSS and probes the moderate-luminosity AGNs that numerically dominate the AGN population in the universe. We also include recently published observations of 19 high-redshift, optically selected AGNs and 46 luminous, low-redshift AGNs from the Bright Quasar Survey. The full sample used in our analysis consists of 333 AGNs, extending out to z ~ 6, with 293 (88%) having X-ray detections. The sample spans five decades in UV luminosity and four decades in X-ray luminosity. We confirm that αOX is strongly anticorrelated with l2500 A (13.6 σ), the highest significance found for this relation to date, and find evidence suggesting that the slope of this relation may be dependent on l2500 A. We find that no significant correlation exists between αOX and redshift (1.3 σ) and constrain the maximum evolution of AGN UV-to-X-ray flux ratios to be less than 30% (1 σ) out to z = 5. Using our sample's high X-ray detection fraction, we also find a significant anticorrelation (3.0 σ) between αOX and l2 keV. We make comparisons to earlier studies on this topic and discuss implications for X-ray and optical AGN luminosity functions.

https://iopscience.iop.org/article/10.1086/503627/pdf

The X-Ray-to-Optical Properties of Optically Selected Active Galaxies over Wide Luminosity and Redshift Ranges

We have determined the cosmological evolution of the density of active galactic nuclei (AGNs) and of their NH distribution as a function of the unabsorbed 2–10 keV luminosity up to redshift 4. We used the HELLAS2XMM sample combined with other published catalogs, yielding a total of 508 AGNs. Our best fit is obtained with a luminosity-dependent density evolution (LDDE) model where low-luminosity (LX � 10 43 ergs s � 1 ) AGNs peak at z � 0:7, while high-luminosity AGNs (LX > 10 45 ergs s � 1 )p eak atz � 2:0. A pure luminosity evolution model (PLE) can instead be rejected. There is evidence that the fraction of absorbed (NH > 10 22 cm � 2 ) AGNs decreases with the intrinsic X-ray luminosity and increases with the redshift. Our best-fit solution provides a good fit to the observed counts, the cosmic X-ray background, and to the observed fraction of absorbed AGNs as a function of the flux in the 10 � 15 ergs cm � 2 s � 1 10 44 ergs s � 1 ) AGNs have a density of 267 deg � 2 at fluxesS2� 10 > 10 � 15 ergs cm � 2 s � 1 . Using these results, "# (

A. Lamastra

Papers

The HELLAS2XMM Survey. VII. The Hard X-Ray Luminosity Function of AGNs up to z = 4: More Absorbed AGNs at Low Luminosities and High Redshifts