X-ray Luminosity Function Research Articles

On the Cosmic Evolution of AGN Obscuration and the X-Ray Luminosity Function: XMM-Newton and Chandra Spectral Analysis of the 31.3 deg2 Stripe 82X

We present X-ray spectral analysis of XMM-Newton and Chandra observations in the 31.3 deg2 Stripe-82X (S82X) field. Of the 6181 unique X-ray sources in this field, we analyze a sample of 2937 candidate active galactic nuclei (AGNs) with solid redshifts and sufficient counts determined by simulations. Our results show an observed population with median values of spectral index , column density log and intrinsic, de-absorbed, 2–10 keV luminosity log, in the redshift range 0–4. We derive the intrinsic, model-independent, fraction of AGNs that are obscured (), finding a significant increase in the obscured AGN fraction with redshift and a decline with increasing luminosity. The average obscured AGN fraction is 57% ± 4% for log L X/erg s−1 > 43. This work constrains the AGN obscuration and spectral shape of the still uncertain high-luminosity and high-redshift regimes (log L X/erg s−1 > 45.5, z > 3), where the obscured AGN fraction rises to 64% ± 12%. We report a luminosity and density evolution of the X-ray luminosity function, with obscured AGNs dominating at all luminosities at z > 2, and unobscured sources prevailing at log L X/erg s−1 > 45 at lower redshifts. Our results agree with the evolutionary models in which the bulk of AGN activity is triggered by gas-rich environments and in a downsizing scenario. Moreover, the black hole accretion density (BHAD) is found to evolve similarly to the star formation rate density, confirming the coevolution between AGN and host galaxy, but suggesting different timescales in their growing history. The derived BHAD evolution shows that Compton-thick AGNs contribute to the accretion history of AGNs as much as all other AGN populations combined.

Constraints on the X-ray luminosity function of AGN at z = 5.7–6.4 with the Extragalactic Serendipitous Swift Survey

ABSTRACT X-ray luminosity functions (XLFs) of active galactic nuclei (AGNs) trace the growth and evolution of supermassive black hole populations across cosmic time. However, current XLF models are poorly constrained at redshifts of z > 6, with a lack of spectroscopic constraints at these high redshifts. In this work, we place limits on the bright-end of the XLF at z = 5.7–6.4 using high-redshift AGN identified within the Extragalactic Serendipitous Swift Survey (ExSeSS) catalogue. Within ExSeSS, we find one serendipitously X-ray detected z > 6 AGN, ATLAS J025.6821-33.4627, with an X-ray luminosity of $L_\mathrm{X}=8.47^{+3.40}_{-3.13}\times 10^{44}\mathrm{erg.s^{-1}}$ and z = 6.31 ± 0.03, making it the highest redshift, spectroscopically confirmed, serendipitously X-ray detected quasar known to date. We also calculate an upper limit on the space density at higher luminosities where no additional sources are found, enabling us to place constraints on the shape of the XLF. Our results are consistent with the rapid decline in the space densities of high-luminosity AGN towards high redshift as predicted by extrapolations of existing parametric models of the XLF. We also find that our X-ray based measurements are consistent with estimates of the bolometric quasar luminosity function based on UV measurements at z ≳ 6, although they require a large X-ray to bolometric correction factor at these high luminosities.

A new method to determine X-ray luminosity functions of AGN and their evolution with redshift

ABSTRACTAlmost all massive galaxies today are understood to contain supermassive black holes (SMBH) at their centres. SMBHs grew by accreting material from their surroundings, emitting X-rays as they did so. X-ray luminosity functions (XLFs) of active galactic nuclei (AGN) have been extensively studied in order to understand the AGN population’s cosmological properties and evolution. We present a new fixed rest-frame method to achieve a more accurate study of the AGN XLF evolution over cosmic time. Normally, XLFs are constructed in a fixed observer-frame energy band, which can be problematic because it probes different rest-frame energies at different redshifts. In the new method, we construct XLFs in the fixed rest-frame band instead, by varying the observed energy band with redshift. We target a rest-frame 2–8 keV band using XMM-Newton and HEAO 1 X-ray data, with seven observer-frame energy bands that vary with redshift for 0 < z < 3. We produce the XLFs using two techniques; one to construct a binned XLF, and one using a maximum likelihood (ML) fit, which makes use of the full unbinned source sample. We find that our ML best-fitting pure luminosity evolution results for both methods are consistent with each other, suggesting that performing XLF evolution studies with the high-redshift data limited to high-luminosity AGN is not very sensitive to the choice of fixed observer-frame or rest-frame energy band, which is consistent with our expectation that high-luminosity AGN typically show little absorption. We have demonstrated the viability of the new method in measuring the XLF evolution.

Forecasts for cosmological measurements based on the angular power spectra of AGN and clusters of galaxies in the SRG/eROSITA all-sky survey

Context. The eROSITA X-ray telescope aboard the Spectrum-Roentgen-Gamma (SRG) orbital observatory, in the course of its all-sky survey, is expected to detect about three million active galactic nuclei (AGNs) and approximately one hundred thousand clusters and groups of galaxies. Such a sample, clean and uniform, complemented with redshift information, will open a new window into the studies of the large-scale structure (LSS) of the Universe and the determination of its cosmological parameters. Aims. The purpose of this work is to assess the prospects of cosmological measurements with the eROSITA sample of AGNs and clusters of galaxies. We assumed the availability of photometric redshift measurements for eROSITA sources and explored the impact of their quality on our forecasts. Methods. As the LSS probe, we use the redshift-resolved angular power spectrum of the density fluctuations of objects. We employed a Fisher-matrix formalism and assumed flat ΛCDM cosmology to forecast the constraining power of eROSITA samples of AGNs and clusters of galaxies. We computed the LSS-relevant characteristics of AGNs and clusters in the framework of the halo model and their X-ray luminosity functions. As the baseline scenario, we considered the full four-year all-sky survey and investigated the impact of reducing the survey length to two years. Results. We find that the accuracy of photometric redshift estimates has a more profound effect on cosmological measurements than the fraction of catastrophic errors. Under realistic assumptions about the photometric redshift quality, the marginalised errors on the cosmological parameters achieve 1 − 10% accuracy depending on the cosmological priors used from other experiments. The statistical significance of Baryon acoustic oscillation detection in angular power spectra of AGNs and clusters of galaxies considered individually achieves 5 − 6σ. Our results demonstrate that the eROSITA sample of AGNs and clusters of galaxies used in combination with currently available photometric redshift estimates will provide cosmological constraints on a par with dedicated optical LSS surveys.

X-ray emission from a rapidly accreting narrow-line Seyfert 1 galaxy at z = 6.56

Context. The space density of X-ray-luminous, blindly selected active galactic nuclei (AGN) traces the population of rapidly accreting super-massive black holes through cosmic time. It is encoded in the X-ray luminosity function, whose bright end remains poorly constrained in the first billion years after the Big Bang as X-ray surveys have thus far lacked the required cosmological volume. With the eROSITA Final Equatorial-Depth Survey (eFEDS), the largest contiguous and homogeneous X-ray survey to date, X-ray AGN population studies can now be extended to new regions of the luminosity–redshift space (L2 − 10 keV > 1045 erg s−1 and z > 6). Aims. The current study aims at identifying luminous quasars at z > 5.7 among X-ray-selected sources in the eFEDS field in order to place a lower limit on black hole accretion well into the epoch of re-ionisation. A secondary goal is the characterisation of the physical properties of these extreme coronal emitters at high redshifts. Methods. Cross-matching eFEDS catalogue sources to optical counterparts from the DESI Legacy Imaging Surveys, we confirm the low significance detection with eROSITA of a previously known, optically faint z = 6.56 quasar from the Subaru High-z Exploration of Low-luminosity Quasars (SHELLQs) survey. We obtained a pointed follow-up observation of the source with the Chandra X-ray telescope in order to confirm the low-significance eROSITA detection. Using new near-infrared spectroscopy, we derived the physical properties of the super-massive black hole. Finally, we used this detection to infer a lower limit on the black hole accretion density rate at z > 6. Results. The Chandra observation confirms the eFEDS source as the most distant blind X-ray detection to date. The derived X-ray luminosity is high with respect to the rest-frame optical emission of the quasar. With a narrow MgII line, low derived black hole mass, and high Eddington ratio, as well as its steep photon index, the source shows properties that are similar to local narrow-line Seyfert 1 galaxies, which are thought to be powered by young super-massive black holes. In combination with a previous high-redshift quasar detection in the field, we show that quasars with L2 − 10 keV > 1045 erg s−1 dominate accretion onto super-massive black holes at z ∼ 6.

Decomposition of galactic X-ray emission with PHOX

Context.X-ray observations of galaxies with high spatial resolution instruments such asChandrahave revealed that major contributions to their diffuse emission originate from X-ray-bright point sources in the galactic stellar field. It has been established that these point sources, called X-ray binaries, are accreting compact objects with stellar donors in a binary configuration. They are classified according to the predominant accretion process: wind-fed in the case of high-mass donors and Roche-lobe mass transfer in the case of low-mass donors. Observationally, it is challenging to reliably disentangle these two populations from each other because of their similar spectra.Aims.We provide a numerical framework with which spatially and spectrally accurate representations of X-ray binary populations can be studied from hydrodynamical cosmological simulations. We construct average spectra, accounting for a hot gas component, and verify the emergence of observed scaling relations between galaxy-wide X-ray luminosity (LX) and stellar mass (M*) and betweenLXand the star-formation rate (SFR).Methods.Using simulated galaxy halos extracted from the (48 h−1 cMpc)3volume of the Magneticum Pathfinder cosmological simulations atz = 0.07, we generate mock spectra with the X-ray photon-simulator PHOX. We extend the PHOXcode to account for the stellar component in the simulation and study the resulting contribution in composite galactic spectra.Results.Well-known X-ray binary scaling relations with galactic SFR andM*emerge self-consistently, verifying our numerical approach. Average X-ray luminosity functions are perfectly reproduced up to the one-photon luminosity limit. Comparing our resultingLX − SFR − M*relation for X-ray binaries with recent observations of field galaxies in the Virgo galaxy cluster, we find significant overlap. Invoking a metallicity-dependent model for high-mass X-ray binaries yields an anticorrelation between mass-weighted stellar metallicity and SFR-normalized luminosity. The spatial distribution of high-mass X-ray binaries coincides with star-formation regions of simulated galaxies, while low-mass X-ray binaries follow the stellar mass surface density. X-ray binary emission is the dominant contribution in the hard X-ray band (2–10 keV) in the absence of an actively accreting central super-massive black hole, and it provides a ∼50% contribution in the soft X-ray band (0.5–2 keV), rivaling the hot gas component.Conclusions.We conclude that our modeling remains consistent with observations despite the uncertainties connected to our approach. The predictive power and easily extendable framework hold great value for future investigations of galactic X-ray spectra.

The Obscured Fraction of Quasars at Cosmic Noon

Statistical studies of X-ray selected active galactic nuclei (AGN) indicate that the fraction of obscured AGN increases with increasing redshift, and the results suggest that a significant part of the accretion growth occurs behind obscuring material in the early universe. We investigate the obscured fraction of highly accreting X-ray AGN at around the peak epoch of supermassive black hole growth utilizing the wide and deep X-ray and optical/IR imaging data sets. A unique sample of luminous X-ray selected AGNs above z > 2 was constructed by matching the XMM-SERVS X-ray point-source catalog with a point-spread function convolved photometric catalog covering the u*–4.5 μm bands. Photometric redshift, hydrogen column density, and 2–10 keV AGN luminosity of the X-ray selected AGN candidates were estimated. Using the sample of 306 2–10 keV detected AGN at above redshift 2, we estimate the fraction of AGN with , assuming parametric X-ray luminosity and absorption functions. The results suggest that of luminous quasars () above redshift 2 are obscured. The fraction indicates an increased contribution of obscured accretion at high redshift than that in the local universe. We discuss the implications of the increasing obscured fraction with increasing redshift based on the AGN obscuration scenarios, which describe obscuration properties in the local universe. Both the obscured and unobscured z > 2 AGN show a broad range of SEDs and morphology, which may reflect the broad variety of host galaxy properties and physical processes associated with the obscuration.

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.

BASS. XXXIII. Swift-BAT Blazars and Their Jets through Cosmic Time

We derive the most up-to-date Swift-Burst Alert Telescope (BAT) blazar luminosity function in the 14–195 keV range, making use of a clean sample of 118 blazars detected in the BAT 105 month survey catalog, with newly obtained redshifts from the BAT Active Galatic Nucleus Spectroscopic Survey. We determine the best-fit X-ray luminosity function for the whole blazar population, as well as for flat-spectrum radio quasars (FSRQs) alone. The main results are: (1) at any redshift, BAT detects the most luminous blazars, above any possible break in their luminosity distribution, which means we cannot differentiate between density and luminosity evolution; (2) the whole blazar population, dominated by FSRQs, evolves positively up to redshift z ∼ 4.3, confirming earlier results and implying lower number densities of blazars at higher redshifts than previously estimated. The contribution of this source class to the cosmic X-ray background at 14–195 keV can range from 5%–18%, while possibly accounting for 100% of the MeV background. We also derived the average 14 keV–10 GeV spectral energy distribution for BAT blazars, which allows us to predict the number counts of sources in the MeV range, as well as the expected number of high-energy (>100 TeV) neutrinos. A mission like COSI will detect 40 MeV blazars, of which two may have coincident neutrino detections. Finally, taking into account beaming selection effects, the distribution and properties of the parent population of these extragalactic jets are derived. We find that the distribution of viewing angles is quite narrow, with most sources aligned within <5° of the line of sight. Moreover, the average Lorentz factor, 〈Γ〉 = 8–12, is lower than previously suggested for these powerful sources.

Properties of intra-cluster low-mass X-ray binaries in Fornax globular clusters

Aims. We present a study of the intra-cluster population of low-mass X-ray binaries (LMXB) residing in globular clusters (GC) in the central 1 deg2 of the Fornax galaxy cluster. Differently from previous studies, which were restricted to the innermost regions of individual galaxies, this work is aimed at comparing the properties of the intra-cluster population of GC-LMXBs with those of the host galaxy. Methods. The data used in this work are a combination of the VLT Survey Telescope (VST) and Chandra observations. We performed a cross-match between the optical and the X-ray catalogue in order to identify the LMXBs residing in GCs. We divided the GC-LMXBs into host-galaxy and intra-cluster objects based on their distance from the nearest galaxy in terms of effective radius (Reff). We found 82 intra-cluster GC-LMXBs and 86 objects that are hosted in galaxies. As the formation of LMXBs also depends on the host GC colour, we performed a Gaussian mixture model to divide the population into red and blue GCs. Results. As has been found for the innermost regions of galaxies, LMXBs tend to form in red and bright GCs in intra-cluster space as well. We find, however, that the likelihood of a red GC to host an LMXB decreases with galactocentric distance, but it remains approximately constant for the blue GC population. Investigating the X-ray properties of the LMXBs residing in GCs, we find a difference in the X-ray luminosity function between the intra-cluster and host-galaxy sample: both follow a power-law down to ∼8.5 × 1037 erg s−1, which is consistent with field LMXBs for the intra-cluster sample, while the latter agree with previous estimates for LMXBs in GCs. We observe a deficiency of bright LMXBs in blue intra-cluster GCs, however. This might indicate a lack of black hole binaries in metal-poor systems. We further investigated the spectral properties of the GC-LMXBs through their hardness-ratio. We detect a tentative difference in the hardness ratio of two populations, where the intra-cluster GC-LMXBs appear to have harder spectra than the host-galaxy objects. We find the same trend when we compare red and blue GC-LMXBs: the spectra of the blue sample are harder spectra than those of the red sample. This result could suggest a relation between the spectral properties of LMXBs and the host GC colour and therefore its metallicity. We discuss the possibilities of spatial biases due to uncertainties in the X-ray spectral response correction and due to contamination by background active galactic nuclei.

BASS. XXX. Distribution Functions of DR2 Eddington Ratios, Black Hole Masses, and X-Ray Luminosities

We determine the low-redshift X-ray luminosity function, active black hole mass function (BHMF), and Eddington ratio distribution function (ERDF) for both unobscured (Type 1) and obscured (Type 2) active galactic nuclei (AGNs), using the unprecedented spectroscopic completeness of the BAT AGN Spectroscopic Survey (BASS) data release 2. In addition to a straightforward 1/V max approach, we also compute the intrinsic distributions, accounting for sample truncation by employing a forward-modeling approach to recover the observed BHMF and ERDF. As previous BHMFs and ERDFs have been robustly determined only for samples of bright, broad-line (Type 1) AGNs and/or quasars, ours are the first directly observationally constrained BHMF and ERDF of Type 2 AGNs. We find that after accounting for all observational biases, the intrinsic ERDF of Type 2 AGNs is significantly more skewed toward lower Eddington ratios than the intrinsic ERDF of Type 1 AGNs. This result supports the radiation-regulated unification scenario, in which radiation pressure dictates the geometry of the dusty obscuring structure around an AGN. Calculating the ERDFs in two separate mass bins, we verify that the derived shape is consistent, validating the assumption that the ERDF (shape) is mass-independent. We report the local AGN duty cycle as a function of mass and Eddington ratio, by comparing the BASS active BHMF with the local mass function for all supermassive black holes. We also present the logN−logS of the Swift/BAT 70 month sources.

The X-Ray Luminosity Function Evolution of Quasars and the Correlation between the X-Ray and Ultraviolet Luminosities

We explore the evolution of the X-ray luminosity function of quasars and the intrinsic correlation between the X-ray and 2500 Å ultraviolet luminosities, utilizing techniques verified in previous works and a sample of over 4000 quasars detected with Chandra and XMM-Newton in the range 0 < z < 5. We find that quasars have undergone significantly less evolution with redshift in their total X-ray luminosity than in other wave bands. We then determine that the best-fit intrinsic power-law correlation between the X-ray and ultraviolet luminosities, of the form , is γ = 0.28 ± 0.03, and we derive the luminosity function and density evolution in the X-ray band. We discuss the implications of these results for models of quasar systems.

The 700 ks Chandra Spiderweb Field

Aims. We present an analysis of the 700 ks Chandra ACIS-S observation of the field around the radio galaxy J1140-2629 (the Spiderweb Galaxy) at z = 2.156, focusing on the nuclear activity in the associated large-scale environment. Methods. We identified unresolved X-ray sources in the field down to flux limits of 1.3 × 10−16 and 3.9 × 10−16 erg s−1 cm−2 in the soft (0.5–2.0 keV) and hard (2–10 keV) band, respectively. We searched for counterparts in the optical, near-infrared, and submillimeter catalogs available in the literature to identify X-ray sources belonging to the protocluster and derived their X-ray properties. Results. We detect 107 X-ray unresolved sources within 5 arcmin (corresponding to 2.5 Mpc) of J1140-2629, among which 13 have optical counterparts with spectroscopic redshift 2.11 &lt; z &lt; 2.20, and 1 source has a photometric redshift consistent with this range. The X-ray-emitting protocluster members are distributed approximately over a ∼3.2 × 1.3 Mpc2 rectangular region. An X-ray spectral analysis for all the sources within the protocluster shows that their intrinsic spectral slope is consistent with an average ⟨Γ⟩∼1.84 ± 0.04. Excluding the Spiderweb Galaxy, the best-fit intrinsic absorption for five protocluster X-ray members is NH &gt; 1023 cm−2, while another six have upper limits of the order of a few times 1022 cm−2. Two sources can only be fitted with very flat Γ ≤ 1, and are therefore considered Compton-thick candidates. The 0.5–10 keV rest-frame luminosities of the 11 Compton-thin protocluster members corrected for intrinsic absorption are greater than 2 × 1043 erg s−1. These values are typical for the bright end of a Seyfert-like distribution and significantly greater than X-ray luminosities expected from star formation activity. The X-ray luminosity function of the AGN in the volume associated to the Spiderweb protocluster in the range 1043 &lt; LX &lt; 1044.5 erg s−1 is at least ten times higher than that in the field at the same redshift and significantly flatter, implying an increasing excess at the bright end. The X-ray AGN fraction is measured to be 25.5 ± 4.5% of the spectroscopically confirmed members in the stellar mass range log(M*/M⊙) &gt; 10.5. This value corresponds to an enhancement factor of 6.0−3.0+9.0 for the nuclear activity with L0.5 − 10 keV &gt; 4 × 1043 erg s−1 with respect to the COSMOS field at comparable redshifts and stellar mass range. Conclusions. We conclude that the galaxy population in the Spiderweb protocluster is characterized by enhanced X-ray nuclear activity triggered by environmental effects on megaparsec scales.

EROSITA study of the globular cluster 47 Tucanae

Aims. We present the results of the analysis of five observations of the globular clutser 47 Tucanae (47 Tuc) with the extended Roentgen Survey with an Imaging Telescope Array (eROSITA) on board the Spektrum-Roentgen-Gamma (Spektr-RG, SRG). We study the X-ray population in the field of one of the most massive globular clusters in our Milky Way. We focused on the classification of point-like sources in the field of 47 Tuc. The unresolved dense core of 47 Tuc (1.7 radius) and also sources that show extended emission are excluded from this study. Methods. We applied different methods of X-ray spectral and timing analysis together with multi-wavelength studies to classify the X-rays sources in the field of 47 Tuc. Results. We detected 888 point-like sources in the energy range of 0.2–5.0 keV. We identified 126 background active galactic nuclei and 25 foreground stars. One of the foreground stars is classified as a variable M dwarf. We also classified 14 X-ray sources as members of 47 Tuc, including one symbiotic star, two quiescent low-mass X-ray binaries, and four cataclysmic variables. There are also five X-ray sources that can either be a cataclysmic variable or a contact binary, and also one X-ray source can be an active binary (type RS CVn). We identified one X-ray binary that belongs to the Small Magellanic Cloud. Moreover, we calculated the X-ray luminosity function of 47 Tuc. No significant population that seems to belong to the globular cluster is observed in the energy range of 0.5–2.0 keV using eROSITA observations.

The eROSITA Final Equatorial-Depth Survey (eFEDS)

Aims.The eROSITA Final Equatorial-Depth Survey has been carried out during the performance verification phase of the Spectrum-Roentgen-Gamma/eROSITA telescope and was completed in November 2019. This survey is designed to provide the first eROSITA-selected sample of clusters and groups and to test the predictions for the all-sky survey in the context of cosmological studies with clusters of galaxies.Methods.In the area of ~140 square degrees covered by eFEDS, 542 candidate clusters and groups of galaxies were detected as extended X-ray sources with theeSASSsource detection algorithm. We performed imaging and spectral analysis of the 542 cluster candidates with eROSITA X-ray data and studied the properties of the sample.Results.We provide the catalog of candidate galaxy clusters and groups detected by eROSITA in the eFEDS field down to a flux of ~10–14erg s–1cm–2in the soft band (0.5–2 keV) within 1’. The clusters are distributed in the redshift range ɀ=[0.01, 1.3] with a median redshift ɀmedian= 0.35. With eROSITA X-ray data, we measured the temperature of the intracluster medium within two radii, 300 kpc and 500 kpc, and constrained the temperature with &gt;2σconfidence level for ~1/5 (102 out of 542) of the sample. The average temperature of these clusters is ~2 keV. Radial profiles of flux, luminosity, electron density, and gas mass were measured from the precise modeling of the imaging data. The selection function, the purity, and the completeness of the catalog are examined and discussed in detail. The contamination fraction is ~1/5 in this sample and is dominated by misidentified point sources. The X-ray luminosity function of the clusters agrees well with the results obtained from other recent X-ray surveys. We also find 19 supercluster candidates in this field, most of which are located at redshifts between 0.1 and 0.5, including one cluster at ɀ ~ 0.36 that was presented previously.Conclusions.The eFEDS cluster and group catalog at the final eRASS equatorial depth provides a benchmark proof of concept for the eROSITA All-Sky Survey extended source detection and characterization. We confirm the excellent performance of eROSITA for cluster science and expect no significant deviations from our pre-launch expectations for the final all-sky survey.

Consistent Analysis of the AGN LF in X-Ray and MIR in the XMM-LSS Field

Abstract The luminosity function of active galactic nuclei (AGN) probes the history of supermassive black hole assembly and growth across cosmic time. To mitigate selection biases, we present a consistent analysis of the AGN luminosity functions (LFs) derived for both X-ray and mid-infrared (MIR) selected AGN in the XMM-Large Scale Structure field. There are 4268 AGN used to construct the MIR luminosity function (IRLF) and 3427 AGN used to construct the X-ray luminosity function (XLF), providing the largest census of the AGN population out to z = 4 in both bands with significant reduction in uncertainties. We are able for the first time to see the knee of the IRLF at z &gt; 2 and observe a flattening of the faint-end slope as redshift increases. The bolometric luminosity density, a proxy for the cosmic black hole accretion history, computed from our LFs, shows a peak at z ≈ 2.25, consistent with recent estimates of the peak in the star formation rate density (SFRD). However, at earlier epochs, the AGN luminosity density is flatter than the SFRD. If confirmed, this result suggests that the build up of black hole mass outpaces the growth of stellar mass in high-mass systems at z ≳ 2.5. This is consistent with observations of redshift z ∼ 6 quasars that lie above the local M − σ relationship. The luminosity density derived from the IRLF is higher than that from the XLF at all redshifts. This is consistent with the dominant role of obscured AGN activity in the cosmic growth of supermassive black holes.

The incidence of X-ray selected AGN in nearby galaxies

ABSTRACT We present the identification and analysis of an unbiased sample of active galactic nuclei (AGN) that lie within the local galaxy population. Using the MPA-JHU catalogue (based on SDSS DR8) and 3XMM DR7 we define a parent sample of 25 949 local galaxies (z ≤ 0.33). After confirming that there was strictly no AGN light contaminating stellar mass and star-formation rate calculations, we identified 917 galaxies with central, excess X-ray emission likely originating from an AGN. We analysed their optical emission lines using the BPT diagnostic and confirmed that such techniques are more effective at reliably identifying sources as AGN in higher mass galaxies: rising from 30 per cent agreement in the lowest mass bin to 93 per cent in the highest. We then calculated the growth rates of the black holes powering these AGN in terms of their specific accretion rates (∝LX/M*). Our sample exhibits a wide range of accretion rates, with the majority accreting at rates $\le 0.5\ \mathrm{ per \, cent}$ of their Eddington luminosity. Finally, we used our sample to calculate the incidence of AGN as a function of stellar mass and redshift. After correcting for the varying sensitivity of 3XMM, we split the galaxy sample by stellar mass and redshift and investigated the AGN fraction as a function of X-ray luminosity and specific black hole accretion rate. From this we found the fraction of galaxies hosting AGN above a fixed specific accretion rate limit of 10−3.5 is constant (at $\approx 1\ \mathrm{ per \, cent}$) over stellar masses of 8 &amp;lt; log M*/M⊙ &amp;lt; 12 and increases (from $\approx 1\ \mathrm{ per \, cent}$ to 10 per cent) with redshift.

X-Ray Binaries in M51 I: Catalog and Statistics

Abstract We used archival data from the Chandra X-ray Observatory (Chandra) and the Hubble Space Telescope, to identify 334 candidate X-ray binary systems and their potential optical counterparts in the interacting galaxy pair NGC 5194/5195 (M51). We present the catalog and data analysis of X-ray and optical properties for those sources, from the deep 892 ks Chandra observations, along with the magnitudes of candidate optical sources as measured in the 8.16 ks Hubble Space Telescope observations. The X-ray luminosity function of the X-ray sources above a few times 1036 erg s−1 follows a power law N ( &gt; L X , b ) ∝ L X , b 1 − α with α = 1.65 ± 0.03. Approximately 80% of sources are variable over a 30 day window. Nearly half of the X-ray sources (173/334) have an optical counterpart within 0.″5.

LeMMINGs – IV. The X-ray properties of a statistically complete sample of the nuclei in active and inactive galaxies from the Palomar sample

ABSTRACT All 280 of the statistically complete Palomar sample of nearby (&amp;lt;120 Mpc) galaxies δ &amp;gt; 20° have been observed at 1.5 GHz as part of the LeMMINGs e-MERLIN legacy survey. Here, we present Chandra X-ray observations of the nuclei of 213 of these galaxies, including a statistically complete sub-set of 113 galaxies in the declination range 40° &amp;lt;δ &amp;lt; 65°. We observed galaxies of all optical spectral types, including ‘active’ galaxies [e.g. low-ionization nuclear emission line regions (LINERs) and Seyferts] and ‘inactive’ galaxies like ${\rm H\, \small {II}}$ galaxies and absorption line galaxies (ALG). The X-ray flux limit of our survey is 1.65 × 10−14 erg s−1 cm−2 (0.3−10 keV). We detect X-ray emission coincident within 2 arcsec of the nucleus in 150/213 galaxies, including 13/14 Seyferts, 68/77 LINERs, 13/22 ALGs and 56/100 ${\rm H\, \small {II}}$ galaxies, but cannot completely rule out contamination from non-AGN processes in sources with nuclear luminosities $\lesssim 10^{39}$ erg s−1. We construct an X-ray Luminosity function (XLF) and find that the local galaxy XLF, when including all active galactic nucleus (AGN) types, can be represented as a single power law of slope −0.54 ± 0.06. The Eddington ratio of the Seyferts is usually 2−4 decades higher than that of the LINERs, ALGs, and ${\rm H\, \small {II}}$ galaxies, which are mostly detected with Eddington ratios $\lesssim 10^{-3}$. Using [${\rm O\, \small {III}}$] line measurements and black hole masses from the literature, we show that LINERs, ${\rm H\, \small {II}}$ galaxies and ALGs follow similar correlations to low luminosities, suggesting that some ‘inactive’ galaxies may harbour AGN.

First tidal disruption events discovered by SRG/eROSITA: X-ray/optical properties and X-ray luminosity function at z &amp;lt; 0.6

ABSTRACT We present the first sample of tidal disruption events (TDEs) discovered during the SRG all-sky survey. These 13 events were selected among X-ray transients detected in the 0° &amp;lt; l &amp;lt; 180° hemisphere by eROSITA during its second sky survey (2020 June 10 to December 14) and confirmed by optical follow-up observations. The most distant event occurred at z = 0.581. One TDE continued to brighten at least 6 months. The X-ray spectra are consistent with nearly critical accretion on to black holes of a few ×103 to $10^8\, \mathrm{ M}_\odot$, although supercritical accretion is possibly taking place. In two TDEs, a spectral hardening is observed 6 months after the discovery. Four TDEs showed an optical brightening apart from the X-ray outburst. The other nine TDEs demonstrate no optical activity. All 13 TDEs are optically faint, with Lg/LX &amp;lt; 0.3 (Lg and LX being the g band and 0.2–6 keV luminosity, respectively). We have constructed a TDE X-ray luminosity function, which can be fit by a power law with a slope of −0.6 ± 0.2, similar to the trend observed for optically selected TDEs. The total rate is estimated at (1.1 ± 0.5) × 10−5 TDEs per galaxy per year, an order of magnitude lower than inferred from optical studies. This suggests that X-ray bright events constitute a minority of TDEs, consistent with models predicting that X-rays can only be observed from directions close to the axis of a thick accretion disc formed from the stellar debris. Our TDE detection threshold can be lowered by a factor of ∼2, which should allow a detection of ∼700 TDEs by the end of the SRG survey.