Hard X-ray Luminosity Function Research Articles

X-Ray Sources in the SSA22 Chandra Field

The Hawaii Survey Field SSA22 is the fourth deepest Chandra X-ray field. To allow for the fullest exploration of this field, we present new optical spectroscopy from Keck/DEIMOS and Keck/LRIS, which, in combination with the literature, brings the spectroscopic completeness of the 2–8 keV sample to 62%. We also make optical spectral classifications and estimate photometric redshifts for the sources without spectroscopic redshifts. We then determine hard X-ray luminosity functions (XLFs) for the full sample of active galactic nuclei (AGNs), as well as for the broad-line AGNs (BLAGNs) and the non-BLAGNs separately. Our XLF for the full sample is in good agreement with the literature, showing relatively strong evolution over the redshift range 0.25 ≤ z < 4. The XLFs for the BLAGNs and the non-BLAGNs imply distinct evolution with redshift, with BLAGNs becoming increasingly dominant at higher redshifts and X-ray luminosities.

The quasar luminosity function at redshift 4 with the Hyper Suprime-Cam Wide Survey

Abstract We present the luminosity function of z ∼ 4 quasars based on the Hyper Suprime-Cam Subaru Strategic Program Wide layer imaging data in the g, r, i, z, and y bands covering 339.8 deg2. From stellar objects, 1666 z ∼ 4 quasar candidates are selected via the g-dropout selection down to i = 24.0 mag. Their photometric redshifts cover the redshift range between 3.6 and 4.3, with an average of 3.9. In combination with the quasar sample from the Sloan Digital Sky Survey in the same redshift range, a quasar luminosity function covering the wide luminosity range of M1450 = −22 to −29 mag is constructed. The quasar luminosity function is well described by a double power-law model with a knee at M1450 = −25.36 ± 0.13 mag and a flat faint-end slope with a power-law index of −1.30 ± 0.05. The knee and faint-end slope show no clear evidence of redshift evolution from those seen at z ∼ 2. The flat slope implies that the UV luminosity density of the quasar population is dominated by the quasars around the knee, and does not support the steeper faint-end slope at higher redshifts reported at z &amp;gt; 5. If we convert the M1450 luminosity function to the hard X-ray 2–10 keV luminosity function using the relation between the UV and X-ray luminosity of quasars and its scatter, the number density of UV-selected quasars matches well with that of the X-ray-selected active galactic nuclei (AGNs) above the knee of the luminosity function. Below the knee, the UV-selected quasars show a deficiency compared to the hard X-ray luminosity function. The deficiency can be explained by the lack of obscured AGNs among the UV-selected quasars.

Supermassive black holes in the EAGLE Universe. Revealing the observables of their growth

We investigate the evolution of supermassive black holes in the `Evolution and Assembly of GaLaxies and their Environments' (EAGLE) cosmological hydrodynamic simulations. The largest of the EAGLE volumes covers a $(100 \,\rm cMpc)^3$ and includes state-of-the-art physical models for star formation and black hole growth that depend only on local gas properties. We focus on the black hole mass function, Eddington ratio distribution and the implied duty cycle of nuclear activity. The simulation is broadly consistent with observational constraints on these quantities. In order to make a more direct comparison with observational data, we calculate the soft and hard X-ray luminosity functions of the active galactic nuclei (AGN). Between redshifts $0$ and $1$, the simulation is in agreement with data. At higher redshifts, the simulation tends to underpredict the luminosities of the brightest observed AGN. This may be due to the limited volume of the simulation, or a fundamental deficiency of the underlying model. It seems unlikely that additional unresolved variability can account for this difference. The simulation shows a similar `downsizing' of the AGN population as seen in observational surveys.

Hard X-ray luminosity function of tidal disruption events: First results from the MAXI extragalactic survey

Abstract We derive the first hard X-ray luminosity function (XLF) of stellar tidal disruption events (TDEs) by supermassive black holes (SMBHs), which gives an occurrence rate of TDEs per unit volume as a function of peak luminosity and redshift, utilizing an unbiased sample observed by the Monitor of All-sky X-ray Image (MAXI). On the basis of the light curves characterized by a power-law decay with an index of −5/3, a systematic search using the MAXI data detected four TDEs in the first 37 months of observations, all of which have been found in the literature. To formulate the TDE XLF, we consider the mass function of SMBHs, that of disrupted stars, the specific TDE rate as a function of SMBH mass, and the fraction of TDEs with relativistic jets. We perform an unbinned maximum likelihood fit to the MAXI TDE list and check the consistency with the observed TDE rate in the ROSAT all-sky survey. The results suggest that the intrinsic fraction of the jet-accompanying events is 0.0007%–34%. We confirm that at z ≲ 1.5 the contamination of the hard X-ray luminosity functions of active galactic nuclei by TDEs is not significant and hence that their contribution to the growth of SMBHs is negligible at the redshifts.

The Illustris simulation: the evolving population of black holes across cosmic time

We study the properties of black holes and their host galaxies across cosmic time in the Illustris simulation. Illustris is a large scale cosmological hydrodynamical simulation which resolves a (106.5 Mpc)^3 volume with more than 12 billion resolution elements and includes state-of-the-art physical models relevant for galaxy formation. We find that the black hole mass density for redshifts z = 0 - 5 and the black hole mass function at z = 0 predicted by Illustris are in very good agreement with the most recent observational constraints. We show that the bolometric and hard X-ray luminosity functions of AGN at z = 0 and 1 reproduce observational data very well over the full dynamic range probed. Unless the bolometric corrections are largely underestimated, this requires radiative efficiencies to be on average low, epsilon_r <= 0.1, noting however that in our model radiative efficiencies are degenerate with black hole feedback efficiencies. Cosmic downsizing of the AGN population is in broad agreement with the findings from X-ray surveys, but we predict a larger number density of faint AGN at high redshifts than currently inferred. We also study black hole -- host galaxy scaling relations as a function of galaxy morphology, colour and specific star formation rate. We find that black holes and galaxies co-evolve at the massive end, but for low mass, blue and star-forming galaxies there is no tight relation with either their central black hole masses or the nuclear AGN activity.

Cosmological simulations of black hole growth: AGN luminosities and downsizing

In this study, we present a detailed, statistical analysis of black hole growth and the evolution of active galactic nuclei (AGN) using cosmological hydrodynamic simulations run down to $z=0$. The simulations self-consistently follow radiative cooling, star formation, metal enrichment, black hole growth and associated feedback processes from both supernovae typeII/Ia and AGN. We consider two simulation runs, one with a large co-moving volume of $(500\ \mathrm{Mpc})^3$ and one with a smaller volume of $(68\ \mathrm{Mpc})^3$ but with a by a factor of almost 20 higher mass resolution. Consistently with previous results, our simulations can widely match observed black hole properties of the local Universe. Furthermore, our simulations can successfully reproduce the evolution of the bolometric AGN luminosity function for both the low-luminosity and the high-luminosity end up to $z=3.0$. In addition, the smaller but higher resolution run is able to match the observational data of the low bolometric luminosity end at higher redshifts $z=3-4$. We also perform a direct comparison with the observed soft and hard X-ray luminosity functions of AGN, including an empirical correction for a torus-level obscuration, and find a similarly good agreement. These results nicely demonstrate that the observed "anti-hierarchical" trend in the AGN number density evolution (i.e. the number densities of luminous AGN peak at higher redshifts than those of faint AGN) is self-consistently predicted by our simulations. Implications of this downsizing behaviour on active black holes, their masses and Eddington-ratios are discussed. Overall, the downsizing behaviour in the AGN number density as a function of redshift can be mainly attributed to the evolution of the gas density in the resolved vicinity of a (massive) black hole. (shortened)

On the contribution of active galactic nuclei to reionization

The electron scattering optical depth constraints on reionization suggest that there may be other sources that contribute to the ionization of hydrogen aside from observable star forming galaxies. Often the calculated value of the electron scattering optical depth, \tau_{es}, falls below the measurements derived from observations of the CMB or an assumption about non-observable sources must be made in order to reach agreement. Here, we calculate the hydrogen ionization fraction as a function of redshift and the electron scattering optical depth from both galaxies and active galactic nuclei (AGN) factoring in the secondary collisional ionizations from the AGN X-ray emission. In this paper we use the most current determination of the evolving hard X-ray luminosity function and extrapolate its evolution beyond z = 6. The AGN spectral energy distributions (SEDs) include both UV and X-ray ionizing photons. To search for the largest possible effect, all AGN are assumed to have \lambda_{Edd} = 1.0 and be completely unobscured. The results show that AGNs produce a perturbative effect on the reionization of hydrogen and remains in agreement with current constraints. Our calculations find the epoch of reionization still ends at z \approx 6 and only increases the electron scattering optical depth by ~2.3% under the most optimal conditions. This can only be moderately increased by assuming a constant black hole mass of M_{BH} = 10^(5) solar masses. As a result, we conclude that there is a need for other sources beyond observable galaxies and AGNs that contribute to the reionization of hydrogen at z > 6.

THE SDSS-III BARYON OSCILLATION SPECTROSCOPIC SURVEY: THE QUASAR LUMINOSITY FUNCTION FROM DATA RELEASE NINE

We present a new measurement of the optical Quasar Luminosity Function (QLF), using data from the Sloan Digital Sky Survey-III: Baryon Oscillation Spectroscopic Survey (SDSS-III: BOSS). From the SDSS-III Data Release Nine (DR9), we select a uniform sample of 22,301 i<=21.8 quasars over an area of 2236 sq. deg with confirmed spectroscopic redshifts between 2.2<z<3.5, filling in a key part of the luminosity-redshift plane for optical quasar studies. We derive the completeness of the survey through simulated quasar photometry, and check this completeness estimate using a sample of quasars selected by their photometric variability within the BOSS footprint. We investigate the level of systematics associated with our quasar sample using the simulations, in the process generating color-redshift relations and a new quasar k-correction. We probe the faint end of the QLF to M_i(z=2.2) = -24.5 and see a clear break in the QLF at all redshifts up to z=3.5. We find that a log-linear relation (in log[Phi*] - M*) for a luminosity and density evolution (LEDE) model adequately describes our data within the range 2.2<z<3.5; across this interval the break luminosity increases by a factor of ~2.3 while Phi* declines by a factor of ~6. At z<2.2 our data is reasonably well fit by a pure luminosity evolution (PLE) model. We see only a weak signature of "AGN downsizing", in line with recent studies of the hard X-ray luminosity function. We compare our measured QLF to a number of theoretical models and find that models making a variety of assumptions about quasar triggering and halo occupation can fit our data over a wide range of redshifts and luminosities.

Origin of the antihierarchical growth of black holes

Observational studies have revealed a "downsizing" trend in black hole (BH) growth: the number densities of luminous AGN peak at higher redshifts than those of faint AGN. This would seem to imply that massive black holes formed before low mass black holes, in apparent contradiction to hierarchical clustering scenarios. We investigate whether this observed "downsizing" in BH growth is reproduced in a semi-analytic model for the formation and evolution of galaxies and black holes, set within the hierarchical paradigm for structure formation (Somerville et al. 2008; S08). In this model, black holes evolve from light seeds (\sim100M\odot) and their growth is merger-driven. The original S08 model (baseline model) reproduces the number density of AGN at intermediate redshifts and luminosities, but underproduces luminous AGN at very high redshift (z > 3) and overproduces them at low redshift (z < 1). In addition, the baseline model underproduces low-luminosity AGN at low redshift (z < 1). To solve these problems we consider several modifications to the physical processes in the model: (1) a 'heavy' black hole seeding scenario (2) a sub-Eddington accretion rate ceiling that depends on the cold gas fraction, and (3) an additional black hole accretion mode due to disk instabilities. With these three modifications, the models can explain the observed downsizing, successfully reproduce the bolometric AGN luminosity function and simultaneously reproduce galaxy and black hole properties in the local Universe. We also perform a comparison with the observed soft and hard X-ray luminosity functions of AGN, including an empirical correction for torus-level obscuration, and reach similar conclusions. Our best-fit model suggests a scenario in which disk instabilities are the main driver for moderately luminous Seyfert galaxies at low redshift, while major mergers are the main trigger for luminous AGN.

Luminosity function of high-mass X-ray binaries and anisotropy in the distribution of active galactic nuclei toward the Large Magellanic Cloud

In 2003-2012, the INTEGRAL observatory has performed long-term observations of the Large Magellanic Cloud (LMC). At present, this is one of the deepest hard X-ray (20-60 keV) surveys of extragalactic fields in which more than 20 sources of different natures have been detected. We present the results of a statistical analysis of the population of high-mass X-ray binaries in the LMC and active galactic nuclei (AGNs) observed in its direction. The hard X-ray luminosity function of high-mass X-ray binaries is shown to be described by a power law with a slope alpha~1.8, that in agreement with the luminosity function measurements both in the LMC itself, but made in the soft X-ray energy band, and in other galaxies. At the same time, the number of detected AGNs toward the LMC turns out to be considerably smaller than the number of AGNs registered in other directions, in particular, toward the source 3C 273. The latter confirms the previously made assumption that the distribution of matter in the local Universe is nonuniform.

A TALE OF TWO POPULATIONS: THE CONTRIBUTION OF MERGER AND SECULAR PROCESSES TO THE EVOLUTION OF ACTIVE GALACTIC NUCLEI

Due to the co-evolution of supermassive black holes and their host galaxies, understanding the mechanisms that trigger active galactic nuclei (AGN) are imperative to understanding galaxy evolution and the formation of massive galaxies. It is observationally difficult to determine the trigger of a given AGN due to the difference between the AGN lifetime and triggering timescales. Here, we utilize AGN population synthesis modeling to determine the importance of different AGN triggering mechanisms. An AGN population model is computed by combining an observationally motivated AGN triggering rate and a theoretical AGN light curve. The free parameters of the AGN light curve are constrained by minimizing a \chi squared test with respect to the observed AGN hard X-ray luminosity function. The observed black hole space density, AGN number counts, and X-ray background spectrum are also considered as observational constraints. It is found that major mergers are not able to account for the entire AGN population. Therefore, non-merger processes, such as secular mechanisms, must also trigger AGN. Indeed, non-merger processes are the dominant AGN triggering mechanism at z \lesssim 1--1.5. Furthermore, the shape and evolution of the black hole mass function of AGN triggered by major mergers is intrinsically different from the shape and evolution of the black hole mass function of AGN triggered by secular processes.

Evolution of the luminosity function and obscuration of active galactic nuclei: comparison between X-ray and infrared

We present a detailed comparison between the 2-10 keV hard X-ray and infrared (IR) luminosity function (LF) of active galactic nuclei (AGN). The composite X-ray to IR spectral energy distributions (SEDs) of AGN used for connecting the hard X-ray LF (HXLF) and IR LF (IRLF) are modeled with a simple but well tested torus model based on the radiative transfer and photoionization code CLOUDY. Four observational determinations of the evolution of 2-10 keV HXLF and six evolution models of the obscured type-2 AGN fraction ($f_2$) have been considered. The 8.0 and 15 \micron LFs for the total, unobscured type-1 and obscured type-2 AGN are predicted from the HXLFs, and then compared with the measurements currently available. We find that the IRLFs predicted from HXLFs tend to underestimate the number of the most IR-luminous AGN. This is independent of the choices of HXLF and $f_2$, and even more obvious for the HXLFs recently measured. We show that the discrepancy between the HXLFs and IRLFs can be largely resolved when the anticorrelation between the UV to X-ray slope $\alpha_{\mathrm{ox}}$ and UV luminosity $L_{\rm UV}$ is appropriately considered. We also discuss other possible explanations for the discrepancy, such as the missing population of Compton-thick AGN and possible contribution of star-formation in the host to the mid-IR. Meanwhile, we find that the HXLFs and IRLFs of AGN can be more consistent with each other if the obscuration mechanisms of quasars and Seyferts are assumed to be different, corresponding to their different triggering and fueling mechanisms. More accurate measurements of the IRLFs of AGN, especially that determined at smaller redshift bins and more accurately separated to that for type-1 and type-2, are very helpful for clarifying these interesting issues.

LIFTING THE VEIL ON OBSCURED ACCRETION: ACTIVE GALACTIC NUCLEI NUMBER COUNTS AND SURVEY STRATEGIES FOR IMAGING HARD X-RAY MISSIONS

Finding and characterizing the population of active galactic nuclei (AGNs) that produces the X-ray background (XRB) is necessary to connect the history of accretion to observations of galaxy evolution at longer wavelengths. The year 2012 will see the deployment of the first hard X-ray imaging telescope that, through deep extragalactic surveys, will be able to measure the AGN population at the energies where the XRB peaks (~20-30 keV). Here, we present predictions of AGN number counts in three hard X-ray bandpasses: 6-10 keV, 10-30 keV and 30-60 keV. Separate predictions are presented for the number counts of Compton thick AGNs, the most heavily obscured active galaxies. The number counts are calculated for five different models of the XRB that differ in the assumed hard X-ray luminosity function, the evolution of the Compton thick AGNs, and the underlying AGN spectral model. The majority of the hard X-ray number counts will be Compton thin AGNs, but there is a >10x increase in the Compton thick number counts from the 6-10 keV to the 10-30 keV band. The Compton thick population show enough variation that a hard X-ray number counts measurement will constrain the models. The computed number counts are used to consider various survey strategies for the NuSTAR mission, assuming a total exposure time of 6.2 Ms. We find that multiple surveys will allow a measurement of Compton thick evolution. The predictions presented here should be useful for all future imaging hard X-ray missions.

THE OPTX PROJECT. IV. HOW RELIABLE IS [O III] AS A MEASURE OF AGN ACTIVITY?

We compare optical and hard X-ray identifications of AGNs using a uniformly selected (above a flux limit of f_2-8 keV = 3.5e-15 erg/cm2/s) and highly optically spectroscopically complete ( > 80% for f_2-8 keV > 1e-14 erg/cm2/s and > 60% below) 2-8 keV sample observed in three Chandra fields (CLANS, CLASXS, and the CDF-N). We find that empirical emission-line ratio diagnostic diagrams misidentify 20-50% of the X-ray selected AGNs that can be put on these diagrams as star formers, depending on which division is used. We confirm that there is a large (2 orders in magnitude) dispersion in the log ratio of the [OIII]5007A to hard X-ray luminosities for the non-broad line AGNs, even after applying reddening corrections to the [OIII] luminosities. We find that the dispersion is similar for the broad-line AGNs, where there is not expected to be much X-ray absorption from an obscuring torus around the AGN nor much obscuration from the galaxy along the line-of-sight if the AGN is aligned with the galaxy. We postulate that the X-ray selected AGNs that are misidentified by the diagnostic diagrams have low [OIII] luminosities due to the complexity of the structure of the narrow-line region, which causes many ionizing photons from the AGN not to be absorbed. This would mean that the [OIII] luminosity can only be used to predict the X-ray luminosity to within a factor of ~3 (one sigma). Despite selection effects, we show that the shapes and normalizations of the [OIII] and transformed hard X-ray luminosity functions show reasonable agreement, suggesting that the [OIII] samples are not finding substantially more AGNs at low redshifts than hard X-ray samples.

The evolution of the hard X-ray luminosity function of AGN

We present new observational determinations of the evolution of the 2-10keV X-ray luminosity function (XLF) of AGN. We utilise data from a number of surveys including both the 2Ms Chandra Deep Fields and the AEGIS-X 200ks survey, enabling accurate measurements of the evolution of the faint end of the XLF. We combine direct, hard X-ray selection and spectroscopic follow-up or photometric redshift estimates at z<1.2 with a rest-frame UV colour pre-selection approach at higher redshifts to avoid biases associated with catastrophic failure of the photometric redshifts. Only robust optical counterparts to X-ray sources are considered using a likelihood ratio matching technique. A Bayesian methodology is developed that considers redshift probability distributions, incorporates selection functions for our high redshift samples, and allows robust comparison of different evolutionary models. We find that the XLF retains the same shape at all redshifts, but undergoes strong luminosity evolution out to z~1, and an overall negative density evolution with increasing redshift, which thus dominates the evolution at earlier times. We do not find evidence that a Luminosity-Dependent Density Evolution, and the associated flattening of the faint-end slope, is required to describe the evolution of the XLF. We find significantly higher space densities of low-luminosity, high-redshift AGN than in prior studies, and a smaller shift in the peak of the number density to lower redshifts with decreasing luminosity. The total luminosity density of AGN peaks at z=1.2+/-0.1, but there is a mild decline to higher redshifts. We find >50% of black hole growth takes place at z>1, with around half in Lx<10^44 erg/s AGN.

Properties of the Galactic population of cataclysmic variables in hard X-rays

We measure the spatial distribution and hard X-ray luminosity function of cataclysmic variables (CVs) using the INTEGRAL all-sky survey in the 17-60 keV energy band. The vast majority of the INTEGRAL detected CVs are intermediate polars with luminosities in the range 10^{32}-10^{34} erg/sec. The scale height of the Galactic disk population of CVs is found to be 130{+90}{-50} pc. The CV luminosity function measured with INTEGRAL in hard X-rays is compatible with that previously determined at lower energies (3--20 keV) using a largely independent sample of sources detected by RXTE (located at |b|>10deg as opposed to the INTEGRAL sample, strongly concentrated to the Galactic plane). The cumulative 17-60 keV luminosity density of CVs per unit stellar mass is found to be (1.3+/-0.3)x10^{27} erg/sec/Msun and is thus comparable to that of low-mass X-ray binaries in this energy band. Therefore, faint but numerous CVs are expected to provide an important contribution to the cumulative hard X-ray emission of galaxies.

A synthesis model for AGN evolution: supermassive black holes growth and feedback modes

We present a comprehensive synthesis model for the active galactic nuclei (AGN) evolution and the growth of supermassive black hole (SMBH) in the Universe. We assume that black holes accrete in just three distinct physical states, or 'modes': at low Eddington ratio, only a radiatively inefficient, kinetically dominated mode is allowed [low kinetic (LK)]; at high Eddington ratio, instead, AGN may display both a purely radiative [radio quiet, high radiative (HR)] and a kinetic [radio loud, high kinetic (HK)] mode. We solve the continuity equation for the black hole mass function using the locally determined one as a boundary condition, and the hard X-ray luminosity function as tracer of the AGN growth rate distribution, supplemented with a luminosity-dependent bolometric correction and an absorbing column distribution. Differently from most previous semi-analytic and numerical models for black hole growth, we do not assume any specific distribution of Eddington ratios, rather we determine it empirically by coupling the mass and luminosity functions and a set of fundamental relations between observables in the three accretion modes. SMBH always show a very broad accretion rate distribution, and we discuss the profound consequences of this fact for our understanding of observed AGN fractions in galaxies, as well as for the empirical determination of SMBH mass functions with large surveys. We confirm previous results and clearly demonstrate that, at least for z≤ 1.5, SMBH mass function evolves antihierarchically, i.e. the most massive holes grew earlier and faster than less massive ones. For the first time, we find hints of a reversal of such a downsizing behaviour at redshifts above the peak of the black hole accretion rate density (z ≈2). We also derive tight constraints on the (mass-weighted) average radiative efficiency of AGN: under the simplifying assumption that the mass density of both high redshift (z ∼5) and 'wandering' black holes ejected from galactic nuclei after merger events are negligible compared to the local mass density, we find that 0.065 0.07, where ξ 0 is the local SMBH mass density in units of 4.3 x 10 5 M ⊙ Mpc -3 . We trace the cosmological evolution of the kinetic luminosity function of AGN, and find that the overall efficiency of SMBH in converting accreted rest mass energy into kinetic power, e kin , ranges between e kin ≃ (3-5) x 10 -3 , depending on the choice of the radio core luminosity function. Such a 'kinetic efficiency' varies however strongly with SMBH mass and redshift, being maximal for very massive holes at late times, as required for the AGN feedback by many galaxy formation models in cosmological contexts.

The Luminosity Function of X‐Ray–selected Active Galactic Nuclei: Evolution of Supermassive Black Holes at High Redshift

We present a measure of the hard (2Y8 keV) X-ray luminosity function (XLF) of AGNs up to z � 5. At high redshifts, the wide area coverage of the Chandra Multiwavength Project is crucial to detect rare and luminous (LX > 10 44 ergs s � 1 ) AGNs. The inclusion of samples from deeper published surveys, such as the Chandra Deep Fields, allows us to span the lower LX range of the XLF. Our sample is selected from both the hard (z 6:3 ; 10 � 16 ergs cm � 2 s � 1 ) and soft (z > 3, f0:5Y2:0 keV > 1:0 ; 10 � 16 ergs cm � 2 s � 1 ) energy band detections. Within our optical magnitude limits (r 0 ; i 0 50%) regarding X-ray source identification (i.e., redshift). Wefind that the luminosity function is similar to that found in previous X-ray surveys up to z � 3 with an evolution dependent on both luminosity and redshift. At z > 3, there is a significant decline in the numbers of AGNs with an evolution rate similar to that found by studies of optically selected QSOs. Based on our XLF, we assess the resolved fraction of the cosmic X-ray background, the cumulative mass density of SMBHs, and the comparison of the mean accretion rate onto SMBHs and the star formation history of galaxies as a function of redshift. A coevolution scenario up to z � 2 is plausible, although at higher redshifts the accretion rate onto SMBHs drops more rapidly. Finally, we highlight the need for better statistics of high-redshift AGNs at zk3, which is achievable with the upcoming Chandra surveys. Subject headingg galaxies: active — quasars: general — surveys — X-rays: galaxies

First Measurement of a Rapid Increase in the AGN Fraction in High-Redshift Clusters of Galaxies

We present the first measurement of the AGN fraction in high-redshift clusters of galaxies (z~0.6) with spectroscopy of one cluster and archival data for three additional clusters. We identify 8 AGN in all four of these clusters from the Chandra data, which are sensitive to AGN with hard X-ray (2-10keV) luminosity L_{X,H} > 10^43 erg/s in host galaxies more luminous than a rest frame M_R 10^43 erg/s we discovered in our earlier study of eight low-redshift clusters with z=0.06-0.31 (average z~0.2). Three of the four high-redshift cluster datasets are sensitive to nearly L_{X,H} > 10^42 erg/s and we identify seven AGN above this luminosity limit, compared to two in eight, low-redshift clusters. Based on membership estimates for each cluster, we determine that the AGN fraction at z~0.6 is f_A(L_X>10^42;M_R 10^43;M_R<-20) = 0.020 (+0.012/-0.008). These values are approximately a factor of 20 greater than the AGN fractions in lower-redshift (average z~0.2) clusters of galaxies and represent a substantial increase over the factors of 1.5 and 3.3 increase, respectively, in the measured space density evolution of the hard X-ray luminosity function over this redshift range. Potential systematic errors would only increase the significance of our result. The cluster AGN fraction increases more rapidly with redshift than the field and the increase in cluster AGN indicates the presence of an AGN Butcher-Oemler Effect.

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

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 and dispersion is assumed. This increases the contribution of the Compton-thin AGN population to the 30 keV XRB intensity by ~ with respect to the case of null dispersion (i.e. a single primary AGN powerlaw with ) 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.