Showing 1–50 of 103 results for author: Aird, J

The NuSTAR Serendipitous Survey: the 80-month catalog and source properties of the high-energy emitting AGN and quasar population

Authors: Claire L. Greenwell, Lizelke Klindt, George B. Lansbury, David J. Rosario, David M. Alexander, James Aird, Daniel Stern, Karl Forster, Michael J. Koss, Franz E. Bauer, Claudio Ricci, John Tomsick, William N. Brandt, Thomas Connor, Peter G. Boorman, Adlyka Annuar, David R. Ballantyne, Chien-Ting Chen, Francesca Civano, Andrea Comastri, Victoria A. Fawcett, Francesca M. Fornasini, Poshak Gandhi, Fiona Harrison, Marianne Heida , et al. (10 additional authors not shown)

Abstract: We present a catalog of hard X-ray serendipitous sources detected in the first 80 months of observations by the Nuclear Spectroscopic Telescope Array (NuSTAR). The NuSTAR serendipitous survey 80-month (NSS80) catalog has an unprecedented $\sim$ 62 Ms of effective exposure time over 894 unique fields (a factor of three increase over the 40-month catalog), with an areal coverage of $\sim $36 deg… ▽ More We present a catalog of hard X-ray serendipitous sources detected in the first 80 months of observations by the Nuclear Spectroscopic Telescope Array (NuSTAR). The NuSTAR serendipitous survey 80-month (NSS80) catalog has an unprecedented $\sim$ 62 Ms of effective exposure time over 894 unique fields (a factor of three increase over the 40-month catalog), with an areal coverage of $\sim $36 deg$^2$, larger than all NuSTAR extragalactic surveys. NSS80 provides 1274 hard X-ray sources in the $3-24$ keV band (822 new detections compared to the previous 40-month catalog). Approximately 76% of the NuSTAR sources have lower-energy ($<10$ keV) X-ray counterparts from Chandra, XMM-Newton, and Swift-XRT. We have undertaken an extensive campaign of ground-based spectroscopic follow-up to obtain new source redshifts and classifications for 427 sources. Combining these with existing archival spectroscopy provides redshifts for 550 NSS80 sources, of which 547 are classified. The sample is primarily composed of active galactic nuclei (AGN), detected over a large range in redshift ($z$ = 0.012-3.43), but also includes 58 spectroscopically confirmed Galactic sources. In addition, five AGN/galaxy pairs, one dual AGN system, one BL Lac candidate, and a hotspot of 4C 74.26 (radio quasar) have been identified. The median rest-frame $10-40$ keV luminosity and redshift of the NSS80 are $\langle{L_\mathrm{10-40 keV}}\rangle$ = 1.2 $\times$ 10$^{44}$ erg s$^{-1}$ and $\langle z \rangle = 0.56$. We investigate the optical properties and construct composite optical spectra to search for subtle signatures not present in the individual spectra, finding an excess of redder BL AGN compared to optical quasar surveys predominantly due to the presence of the host-galaxy and, at least in part, due to dust obscuration. △ Less

Submitted 26 April, 2024; originally announced April 2024.

Comments: Accepted for publication in ApJ:S. 57 pages, 32 figures

Accretion properties of X-ray AGN: Evidence for radiation-regulated obscuration with redshift-dependent host galaxy contribution

Authors: Brivael Laloux, Antonis Georgakakis, David M. Alexander, Johannes Buchner, Carolina Andonie, Nischal Acharya, James Aird, Alba V. Alonso-Tetilla, Angela Bongiorno, Ryan C. Hickox, Andrea Lapi, Blessing Musiimenta, Cristina Ramos Almeida, Carolin Vellforth, Francesco Shankar

Abstract: We adopt a Bayesian X-ray spectral approach to investigate the accretion properties of unobscured ($20<\log(N_{\rm H}/{\rm cm}^{-2}<22$) and obscured ($22< \log(N_{\rm H}/{\rm cm}^{-2}<24$) active galactic nuclei (AGN) to shed light on the orientation vs evolution scenarios for the origin of the obscuring material. For a sample of 3882 X-ray-selected AGN from the Chandra COSMOS Legacy, AEGIS and C… ▽ More We adopt a Bayesian X-ray spectral approach to investigate the accretion properties of unobscured ($20<\log(N_{\rm H}/{\rm cm}^{-2}<22$) and obscured ($22< \log(N_{\rm H}/{\rm cm}^{-2}<24$) active galactic nuclei (AGN) to shed light on the orientation vs evolution scenarios for the origin of the obscuring material. For a sample of 3882 X-ray-selected AGN from the Chandra COSMOS Legacy, AEGIS and CDFS extragalactic surveys, we constrain their stellar masses, $M_\star$, intrinsic X-ray luminosities, $L_{\rm X}$, obscuring column densities, $N_{\rm H}$, and specific accretion rates $λ\propto L_{\rm X}/M_\star$. By combining these observables within a Bayesian non-parametric approach, we infer, for the first time, the specific accretion rate distribution (SARD) of obscured and unobscured AGN to $z\approx3$, i.e. the probability of a galaxy with mass $M_\star$ at redshift $z$ hosting an AGN with column density $N_{\rm H}$ and specific accretion rate $λ$. Our findings indicate that (1) both obscured and unobscured SARDs share similar shapes, shifting towards higher accretion rates with redshift, (2) unobscured SARDs exhibit a systematic offset towards higher $λ$ compared to obscured SARD for all redshift intervals, (3) the obscured AGN fraction declines sharply at $\logλ_{\rm break} \sim-2$ for $z <0.5$, but shifts to higher $λ$ values with increasing redshift, (4) the incidence of AGN within the theoretically unstable blow-out region of the $λ-N_{\rm H}$ plane increases with redshift. These observations provide compelling evidence for AGN "downsizing" and radiation-regulated nuclear-scale obscuration with an increasing host galaxy contribution towards higher redshifts. △ Less

Submitted 11 March, 2024; originally announced March 2024.

Comments: Submitted to MNRAS

The High Energy X-ray Probe (HEX-P): Bringing the Cosmic X-ray Background into Focus

Authors: Francesca Civano, Xiurui Zhao, Peter Boorman, Stefano Marchesi, Tonima Ananna, Samantha Creech, Chien-Ting Chen, Ryan Hickox, Daniel Stern, Kristin Madsen, Javier A. García, Ross Silver, James Aird, David M. Alexander, Mislav Balokovíc, William N. Brandt, Johannes Buchner, Poshak Gandhi, Elias Kammoun, Stephanie LaMassa, Giorgio Lanzuisi, Andrea Merloni, Alberto Moretti, Kirpal Nandra, Emanuele Nardini , et al. (7 additional authors not shown)

Abstract: Since the discovery of the cosmic X-ray background (CXB), astronomers have strived to understand the accreting supermassive black holes (SMBHs) contributing to its peak in the 10-40 keV band. Existing soft X-ray telescopes could study this population up to only 10 keV, and, while NuSTAR (focusing on 3--24 keV) made great progress, it also left significant uncertainties in characterizing the hard X… ▽ More Since the discovery of the cosmic X-ray background (CXB), astronomers have strived to understand the accreting supermassive black holes (SMBHs) contributing to its peak in the 10-40 keV band. Existing soft X-ray telescopes could study this population up to only 10 keV, and, while NuSTAR (focusing on 3--24 keV) made great progress, it also left significant uncertainties in characterizing the hard X-ray population, crucial for calibrating current population synthesis models. This paper presents an in-depth analysis of simulations of two extragalactic surveys (deep and wide) with the High-Energy X-ray Probe (HEX-P), each observed for 2 Ms. Applying established source detection techniques, we show that HEX-P surveys will reach a flux of $\sim$10$^{-15}$ erg s$^{-1}$ cm$^{-2}$ in the 10-40 keV band, an order of magnitude fainter than current NuSTAR surveys. With the large sample of new hard X-ray detected sources ($\sim2000$), we showcase HEX-P's ability to resolve more than 80% of the CXB up to 40 keV into individual sources. The expected precision of HEX-P's resolved background measurement will allow us to distinguish between population synthesis models of SMBH growth. HEX-P leverages accurate broadband (0.5-40 keV) spectral analysis and the combination of soft and hard X-ray colors to provide obscuration constraints even for the fainter sources, with the overall objective of measuring the Compton-thick fraction. With unprecedented sensitivity in the 10--40 keV band, HEX-P will explore the hard X-ray emission from AGN to flux limits never reached before, thus expanding the parameter space for serendipitous discoveries. Consequently, it is plausible that new models will be needed to capture the population HEX-P will unveil. △ Less

Submitted 16 November, 2023; v1 submitted 8 November, 2023; originally announced November 2023.

Comments: 28 pages, 13 Figures, 2 tables; submitted to Frontiers in Astronomy and Space Sciences; minor edits in the text, abstract and figure labels

TRINITY IV: Predictions for Supermassive Black Holes at $z \gtrsim 7$

Authors: Haowen Zhang, Peter Behroozi, Marta Volonteri, Joseph Silk, Xiaohui Fan, James Aird, Jinyi Yang, Feige Wang, Philip F. Hopkins

Abstract: We present predictions for the high-redshift halo-galaxy-supermassive black hole (SMBH) connection from the TRINITY model. Constrained by a comprehensive compilation of galaxy ($0\leq z \leq 10$) and SMBH datasets ($0\leq z \leq 6.5$), TRINITY finds: 1) The number of SMBHs with $M_\bullet > 10^9 M_\odot$ in the observable Universe increases by six orders of magnitude from $z\sim10$ to $z\sim2$, an… ▽ More We present predictions for the high-redshift halo-galaxy-supermassive black hole (SMBH) connection from the TRINITY model. Constrained by a comprehensive compilation of galaxy ($0\leq z \leq 10$) and SMBH datasets ($0\leq z \leq 6.5$), TRINITY finds: 1) The number of SMBHs with $M_\bullet > 10^9 M_\odot$ in the observable Universe increases by six orders of magnitude from $z\sim10$ to $z\sim2$, and by another factor of $\sim 3$ from $z\sim2$ to $z=0$; 2) The $M_\bullet > 10^9/10^{10} M_\odot$ SMBHs at $z\sim 6$ live in haloes with $\sim (2-3)/(3-5) \times 10^{12} M_\odot$; 3) the new JWST AGNs at $7\lesssim z \lesssim 11$ are broadly consistent with the median SMBH mass-galaxy mass relation for AGNs from TRINITY; 4) Seeds from runaway mergers in nuclear star clusters are viable progenitors for the SMBHs in GN-z11 ($z=10.6$) and CEERS_1019 ($z=8.7$); 5) $z=6-10$ quasar luminosity functions from wide area surveys by, e.g., Roman and Euclid, will reduce uncertainties in the $z=6-10$ SMBH mass-galaxy mass relation by up to $\sim 0.5$ dex. △ Less

Submitted 13 September, 2023; originally announced September 2023.

Comments: 15 pages, 12 figures, submitted to MNRAS. Questions and comments are welcome!

The intrinsic X-ray luminosity distribution of an optically-selected SDSS quasar population

Authors: Amy L. Rankine, James Aird, Angel Ruiz, Antonis Georgakakis

Abstract: In active galactic nuclei, the relationship between UV and X-ray luminosity is well studied (often characterised by $α_\text{ox}$) but often with heterogeneous samples. We have parametrized the intrinsic distribution of X-ray luminosity, $L_\text{X}$, for the optically-selected sample of SDSS quasars in the Stripe 82 and XXL fields across redshifts 0.5-3.5. We make use of the available XMM observa… ▽ More In active galactic nuclei, the relationship between UV and X-ray luminosity is well studied (often characterised by $α_\text{ox}$) but often with heterogeneous samples. We have parametrized the intrinsic distribution of X-ray luminosity, $L_\text{X}$, for the optically-selected sample of SDSS quasars in the Stripe 82 and XXL fields across redshifts 0.5-3.5. We make use of the available XMM observations and a custom pipeline to produce Bayesian sensitivity curves that are used to derive the intrinsic X-ray distribution in a hierarchical Bayesian framework. We find that the X-ray luminosity distribution is well described by a Gaussian function in ${\log_{10}}L_\text{X}$ space with a mean that is dependent on the monochromatic 2500A UV luminosity, $L_{2500}$. We also observe some redshift dependence of the distribution. The mean of the $L_\text{X}$ distribution increases with redshift while the width decreases. This weak but significant redshift dependence leads to $L_{2500}$-$L_\text{X}$ and $L_{2500}$-$α_\text{ox}$ relations that evolve with redshift, and we produce a redshift- and $L_{2500}$-dependent $α_\text{ox}$ equation. Neither black hole mass nor Eddington ratio appear to be potential drivers of the redshift evolution. △ Less

Submitted 28 November, 2023; v1 submitted 31 July, 2023; originally announced July 2023.

Comments: 18 pages, 16 figures (including appendices). Accepted for publication in MNRAS

The incidence of AGN in galaxies with different stellar population ages

Authors: Q. Ni, J. Aird, A. Merloni, K. L. Birchall, J. Buchner, M. Salvato, G. Yang

Abstract: It has been argued that recycled gas from stellar mass loss in galaxies might serve as an important fuelling source for black holes (BHs) in their centers. Utilizing spectroscopic samples of galaxies from the Sloan Digital Sky Survey (SDSS) at $z = 0-0.35$ and the Large Early Galaxy Astrophysics Census (LEGA-C) survey at $z = 0.6-1$ that have X-ray coverage from XMM-Newton or Chandra, we test this… ▽ More It has been argued that recycled gas from stellar mass loss in galaxies might serve as an important fuelling source for black holes (BHs) in their centers. Utilizing spectroscopic samples of galaxies from the Sloan Digital Sky Survey (SDSS) at $z = 0-0.35$ and the Large Early Galaxy Astrophysics Census (LEGA-C) survey at $z = 0.6-1$ that have X-ray coverage from XMM-Newton or Chandra, we test this stellar mass loss fuelling scenario by investigating how AGN activity and BH growth vary with the break strength at 4000 $Å$, $\rm D_{n}4000$ (which is closely related to the age of stellar populations), as younger galaxies are considered to have higher stellar mass loss rates. We found that when controlling for host-galaxy properties, the fraction of log $L_{\rm X}$/$M_\star$ > 32 (which roughly corresponds to Eddington ratios $\gtrsim 1$%) AGN and sample-averaged black hole accretion rate ($\rm \overline{BHAR}$) decrease with $\rm D_{n}4000$ among $\rm D_{n}4000$ $\lesssim$ 1.9 galaxies, suggesting a higher level of AGN activity among younger galaxies, which supports the stellar mass loss fuelling scenario. For the oldest and most massive galaxies at $z = 0-0.35$, this decreasing trend is not present anymore. We found that, among these most massive galaxies at low redshift, the fraction of low specific-accretion-rate (31 $<$ log $L_{\rm X}$/$M_\star$ $<$ 32) AGNs increases with $\rm D_{n}4000$, which may be associated with additional fuelling from hot halo gas and/or enhanced accretion capability. △ Less

Submitted 30 June, 2023; originally announced July 2023.

Comments: 24 pages, 28 figures. Accepted for publication in MNRAS

The relationship between the incidence of X-ray selected AGN in nearby galaxies and star-formation rate

Authors: Keir L. Birchall, M. G. Watson, J. Aird, R. L. C. Starling

Abstract: We present the identification and analysis of an X-ray selected AGN sample that lie within the local ($z < 0.35$) galaxy population. From a parent sample of 22,079 MPA-JHU (based on SDSS DR8) galaxies, we identified 917 galaxies with central, excess X-ray emission (from 3XMM-DR7) likely originating from an AGN. We measured the host galaxies' star formation rates and classified them as either star-… ▽ More We present the identification and analysis of an X-ray selected AGN sample that lie within the local ($z < 0.35$) galaxy population. From a parent sample of 22,079 MPA-JHU (based on SDSS DR8) galaxies, we identified 917 galaxies with central, excess X-ray emission (from 3XMM-DR7) likely originating from an AGN. We measured the host galaxies' star formation rates and classified them as either star-forming or quiescent based on their position relative to main sequence of star formation. Only 72% of the X-ray selected sample were identified as AGN using BPT selection; this technique is much less effective in quiescent hosts, only identifying 50% of the X-ray AGN. We also calculated the growth rates of the black holes powering these AGN in terms of their specific accretion rate ($\propto \mathrm{L_X/M_*}$) and found quiescent galaxies, on average, accrete at a lower rate than star-forming galaxies. Finally, we measured the sensitivity function of 3XMM so we could correct for observational bias and construct probability distributions as a function of accretion rate. AGN were found in galaxies across the full range of star formation rates ($\log_{10} \mathrm{SFR/M_\odot\ yr^{-1}} = -3\ \mathrm{to}\ 2$) in both star-forming and quiescent galaxies. The incidence of AGN was enhanced by a factor 2 (at a 3.5$σ$ significance) in star-forming galaxies compared to quiescent galaxies of equivalent stellar mass and redshift, but we also found a significant population of AGN hosted by quiescent galaxies. △ Less

Submitted 19 June, 2023; originally announced June 2023.

Comments: 14 pages, 8 figures, 2 appendices; published in MNRAS

TRINITY III: Quasar Luminosity Functions Decomposed by Halo, Galaxy, and Black Hole Masses and Eddington Ratios from z=0-10

Authors: Haowen Zhang, Peter Behroozi, Marta Volonteri, Joseph Silk, Xiaohui Fan, James Aird, Jinyi Yang, Philip F. Hopkins

Abstract: We present the redshift evolution of quasar luminosity functions decomposed by halo mass, galaxy mass, supermassive black hole (SMBH) mass, and Eddington ratio, as well as SMBH kinetic/radiative energy output ratios from TRINITY, a flexible empirical model that self-consistently infers the halo--galaxy--SMBH connection that match observational data. Key findings include: 1) The normalization of QL… ▽ More We present the redshift evolution of quasar luminosity functions decomposed by halo mass, galaxy mass, supermassive black hole (SMBH) mass, and Eddington ratio, as well as SMBH kinetic/radiative energy output ratios from TRINITY, a flexible empirical model that self-consistently infers the halo--galaxy--SMBH connection that match observational data. Key findings include: 1) The normalization of QLF increases by ~3-4 dex from z~10 to z~4, due to the fast mass build-up of different SMBH populations; 2) From z~4 to z~1, less massive galaxies and SMBHs make up bigger and bigger fractions of QLFs, due to the AGN downsizing effect; 3) At z~0, massive haloes/galaxies/SMBHs are responsible for most bright quasars due to low Eddington ratios among all SMBHs; 4) The bright ends of quasar luminosity functions (QLFs) are dominated by SMBHs that are at least 0.3 dex over-massive relative to the median SMBH mass-galaxy mass relation; 5) QLFs at z~6-7 are dominated by SMBHs accreting at Eddington ratios 0.1 < $η_\mathrm{rad}$ < 1, but super-Eddington AGNs contribute more significantly to QLFs towards z~9-10. △ Less

Submitted 1 March, 2024; v1 submitted 30 May, 2023; originally announced May 2023.

Comments: 18 pages, 14 figures. Accepted by MNRAS. Comments welcome!

TRINITY II: The Luminosity-dependent Bias of the Supermassive Black Hole Mass--Galaxy Mass Relation for Bright Quasars at $z=6$

Authors: Haowen Zhang, Peter Behroozi, Marta Volonteri, Joseph Silk, Xiaohui Fan, Philip F. Hopkins, Jinyi Yang, James Aird

Abstract: Using recent empirical constraints on the dark matter halo--galaxy--supermassive black hole (SMBH) connection from $z=0-7$, we infer how undermassive, typical, and overmassive SMBHs contribute to the quasar luminosity function (QLF) at $z=6$. We find that beyond $L_\mathrm{bol} = 5 \times 10^{46}$ erg/s, the $z=6$ QLF is dominated by SMBHs that are at least 0.3 dex above the $z=6$ median… ▽ More Using recent empirical constraints on the dark matter halo--galaxy--supermassive black hole (SMBH) connection from $z=0-7$, we infer how undermassive, typical, and overmassive SMBHs contribute to the quasar luminosity function (QLF) at $z=6$. We find that beyond $L_\mathrm{bol} = 5 \times 10^{46}$ erg/s, the $z=6$ QLF is dominated by SMBHs that are at least 0.3 dex above the $z=6$ median $M_\bullet-M_*$ relation. The QLF is dominated by typical SMBHs (i.e., within $\pm 0.3$ dex around the $M_\bullet-M_*$ relation) at $L_\mathrm{bol} \lesssim 10^{45}$ erg/s. At $z\sim 6$, the intrinsic $M_\bullet-M_*$ relation for all SMBHs is slightly steeper than the $z=0$ scaling, with a similar normalization at $M_* \sim 10^{11} M_\odot$. We also predict the $M_\bullet-M_*$ relation for $z=6$ bright quasars selected by different bolometric luminosity thresholds, finding very good agreement with observations. For quasars with $L_\mathrm{bol} > 3 \times 10^{46}$ ($10^{48}$) erg/s, the scaling relation is shifted upwards by $\sim0.35$ (1.0) dex for $10^{11} M_\odot$ galaxies. To accurately measure the intrinsic $M_\bullet-M_*$ relation, it is essential to include fainter quasars with $L_\mathrm{bol} \lesssim 10^{45}$ erg/s. At high redshifts, low-luminosity quasars are thus the best targets for understanding typical formation paths for SMBHs in galaxies. △ Less

Submitted 14 March, 2023; originally announced March 2023.

Comments: 5 pages, 3 figures. Submitted to MNRAS Letters. Comments welcome!

The Extragalactic Serendipitous Swift Survey (ExSeSS) -- I. Survey definition and measurements of the X-ray number counts

Authors: Jack N. Delaney, James Aird, Phil A. Evans, Cassandra Barlow-Hall, Julian P. Osborne, Michael G. Watson

Abstract: We present the Extragalactic Serendipitous Swift Survey (ExSeSS), providing a new well-defined sample constructed from the observations performed using the Swift X-ray Telescope. The ExSeSS sample consists of 79,342 sources detected in the medium (1-2 keV), hard (2-10 keV) or total (0.3-10 keV) energy bands, covering 2086.6 deg$^{2}$ of sky across a flux range of… ▽ More We present the Extragalactic Serendipitous Swift Survey (ExSeSS), providing a new well-defined sample constructed from the observations performed using the Swift X-ray Telescope. The ExSeSS sample consists of 79,342 sources detected in the medium (1-2 keV), hard (2-10 keV) or total (0.3-10 keV) energy bands, covering 2086.6 deg$^{2}$ of sky across a flux range of $f_\mathrm{0.3-10keV}\sim10^{-15}-10^{-10}$ erg s$^{-1}$ cm$^{-2}$. Using the new ExSeSS sample we present measurements of the differential number counts of X-ray sources as a function of 2-10 keV flux that trace the population of Active Galactic Nuclei (AGN) in a previously unexplored regime. We find that taking the line-of-sight absorption column density into account has an effect on the differential number count measurements and is vital to obtain agreement with previous results. In the hard band, we obtain a good agreement between the ExSeSS measurements and previous, higher energy data from NuSTAR and Swift/BAT when taking into account the varying column density of the ExSeSS sample as well as the X-ray spectral parameters of each of the samples we are comparing to. We also find discrepancies between the ExSeSS measurements and AGN population synthesis models, indicating a change in the properties of the AGN population over this flux range that is not fully described by current models at these energies, hinting at a larger, moderately obscured population at low redshifts ($z\lesssim0.2$) that the models are not currently taking into account. △ Less

Submitted 14 December, 2022; originally announced December 2022.

Comments: 14 pages, 11 figures, accepted for publication by MNRAS. Supplementary data provided at https://www.swift.ac.uk/2SXPS/exsess/

Chasing Super-Massive Black Hole merging events with $Athena$ and LISA

Authors: L. Piro, M. Colpi, J. Aird, A. Mangiagli, A. C. Fabian, M. Guainazzi, S. Marsat, A. Sesana, P. McNamara, M. Bonetti, E. M. Rossi, N. R. Tanvir, J. G. Baker, G. Belanger, T. Dal Canton, O. Jennrich, M. L. Katz, N. Luetzgendorf

Abstract: The European Space Agency is studying two large-class missions bound to operate in the decade of the 30s, and aiming at investigating the most energetic and violent phenomena in the Universe. $Athena$ is poised to study the physical conditions of baryons locked in large-scale structures from the epoch of their formation, as well as to yield an accurate census of accreting super-massive black holes… ▽ More The European Space Agency is studying two large-class missions bound to operate in the decade of the 30s, and aiming at investigating the most energetic and violent phenomena in the Universe. $Athena$ is poised to study the physical conditions of baryons locked in large-scale structures from the epoch of their formation, as well as to yield an accurate census of accreting super-massive black holes down to the epoch of reionization; LISA will extend the hunt for Gravitational Wave (GW) events to the hitherto unexplored mHz regime. We discuss in this paper the science that their concurrent operation could yield, and present possible $Athena$ observational strategies. We focus on Super-Massive (M$\lesssim10^7\rm M_{\odot}$) Black Hole Mergers (SMBHMs), potentially accessible to $Athena$ up to $z\sim2$. The simultaneous measurement of their electro-magnetic (EM) and GW signals may enable unique experiments in the domains of astrophysics, fundamental physics, and cosmography, such as the magneto-hydrodynamics of fluid flows in a rapidly variable space-time, the formation of coronae and jets in Active Galactic Nuclei, and the measurement of the speed of GW, among others. Key to achieve these breakthrough results will be the LISA capability of locating a SMBHM event with an error box comparable to, or better than the field-of-view of the $Athena$ Wide Field Imager ($\simeq0.4\,$deg$^2$) and $Athena$ capability to slew fast to detect the source during the inspiral phase and the post-merger phase. Together, the two observatories will open in principle the exciting possibility of truly concurrent EM and GW studies of the SMBHMs △ Less

Submitted 28 February, 2023; v1 submitted 24 November, 2022; originally announced November 2022.

Comments: 17 pages, 8 figures. Accepted for publication in MNRAS

The demographics of obscured AGN from X-ray spectroscopy guided by multiwavelength information

Authors: Brivael Laloux, Antonis Georgakakis, Carolina Andonie, David M. Alexander, Angel Ruiz, David J. Rosario, James Aird, Johannes Buchner, Francisco J. Carrera, Andrea Lapi, Cristina Ramos Almeida, Mara Salvato, Francesco Shankar

Abstract: A complete census of Active Galactic Nuclei (AGN) is a prerequisite for understanding the growth of supermassive black holes across cosmic time. A significant challenge toward this goal is the whereabouts of heavily obscured AGN that remain uncertain. This paper sets new constraints on the demographics of this population by developing a methodology that combines X-ray spectral information with pri… ▽ More A complete census of Active Galactic Nuclei (AGN) is a prerequisite for understanding the growth of supermassive black holes across cosmic time. A significant challenge toward this goal is the whereabouts of heavily obscured AGN that remain uncertain. This paper sets new constraints on the demographics of this population by developing a methodology that combines X-ray spectral information with priors derived from multiwavelength observations. We select X-ray AGN in the Chandra COSMOS Legacy survey and fit their $2.2-500μm$ spectral energy distributions with galaxy and AGN templates to determine the mid-infrared ($6μm$) luminosity of the AGN component. Empirical correlations between X-ray and $6μm$ luminosities are then adopted to infer the intrinsic accretion luminosity at X-rays for individual AGN. This is used as prior information in our Bayesian X-ray spectral analysis to estimate physical properties, such as line-of-sight obscuration. Our approach breaks the degeneracies between accretion luminosity and obscuration that affect X-ray spectral analysis, particularly for the most heavily obscured (Compton-Thick) AGN with low photon counts X-ray spectra. The X-ray spectral results are then combined with the selection function of the Chandra COSMOS Legacy survey to derive the AGN space density and a Compton-Thick fraction of $21.0^{+16.1}_{-9.9}\%$ at redshifts $z<0.5$. At higher redshift, our analysis suggests upper limits to the Compton-Thick AGN fraction of $\le 40\%$. These estimates are at the low end of the range of values determined in the literature and underline the importance of multiwavelength approaches for tackling the challenge of heavily obscured AGN demographics. △ Less

Submitted 7 November, 2022; originally announced November 2022.

Comments: Accepted for publication at MNRAS, 21 pages, 20 figures, 6 tables

A super-linear "radio-AGN main sequence'' links mean radio-AGN power and galaxy stellar mass since z$\sim$3

Authors: I. Delvecchio, E. Daddi, M. T. Sargent, J. Aird, J. R. Mullaney, B. Magnelli, D. Elbaz, L. Bisigello, L. Ceraj, S. Jin, B. S. Kalita, D. Liu, M. Novak, I. Prandoni, J. F. Radcliffe, C. Spingola, G. Zamorani, V. Allevato, G. Rodighiero, V. Smolcic

Abstract: Mapping the average AGN luminosity across galaxy populations and over time encapsulates important clues on the interplay between supermassive black hole (SMBH) and galaxy growth. This paper presents the demography, mean power and cosmic evolution of radio AGN across star-forming galaxies (SFGs) of different stellar masses (${M_{*}}$). We exploit deep VLA-COSMOS 3 GHz data to build the rest-frame 1… ▽ More Mapping the average AGN luminosity across galaxy populations and over time encapsulates important clues on the interplay between supermassive black hole (SMBH) and galaxy growth. This paper presents the demography, mean power and cosmic evolution of radio AGN across star-forming galaxies (SFGs) of different stellar masses (${M_{*}}$). We exploit deep VLA-COSMOS 3 GHz data to build the rest-frame 1.4 GHz AGN luminosity functions at 0.1$\leq$$z$$\leq$4.5 hosted in SFGs. Splitting the AGN luminosity function into different ${M_{*}}$ bins reveals that, at all redshifts, radio AGN are both more frequent and more luminous in higher ${M_*}$ than in lower ${M_*}$ galaxies. The cumulative kinetic luminosity density exerted by radio AGN in SFGs peaks at $z$$\sim$2, and it is mostly driven by galaxies with 10.5$\leq$$\log$(${M_{*}}$/${M_{\odot}}$)$<$11. Averaging the cumulative radio AGN activity across all SFGs at each (${M_{*}}$,$z$) results in a "radio-AGN main sequence" that links the time-averaged radio-AGN power $\langle$$L_{1.4}^{AGN}$$\rangle$ and galaxy stellar mass, in the form: $\log$$\langle$[$L_{1.4}^{AGN}$/ W Hz$^{-1}]\rangle$ = (20.97$\pm$0.16) + (2.51$\pm$0.34)$\cdot$$\log$(1+$z$) + (1.41$\pm$0.09)$\cdot$($\log$[${M_{*}}$/${M_{\odot}}$] -10). The super-linear dependence on ${M_{*}}$, at fixed redshift, suggests enhanced radio-AGN activity in more massive SFGs, as compared to star formation. We ascribe this enhancement to both a higher radio AGN duty cycle and a brighter radio-AGN phase in more massive SFGs. A remarkably consistent ${M_{*}}$ dependence is seen for the evolving X-ray AGN population in SFGs. This similarity is interpreted as possibly driven by secular cold gas accretion fueling both radio and X-ray AGN activity in a similar fashion over the galaxy's lifetime. △ Less

Submitted 26 September, 2022; originally announced September 2022.

Comments: Accepted for publication in A&A. 20 pages + Appendices

Journal ref: A&A 668, A81 (2022)

The MOSDEF Survey: A New View of a Remarkable z=1.89 Merger

Authors: Jordan N. Runco, Alice E. Shapley, Mariska Kriek, Michele Cappellari, Michael W. Topping, Ryan L. Sanders, Vasily I. Kokorev, Sedona H. Price, Naveen A. Reddy, Alison L. Coil, Bahram Mobasher, Brian Siana, Tom Zick, Georgios E. Magdis, Gabriel Brammer, James Aird

Abstract: We present a detailed study of a galaxy merger taking place at $z=1.89$ in the GOODS-S field. Here we analyze Keck/MOSFIRE spectroscopic observations from the MOSFIRE Deep Evolution Field (MOSDEF) survey along with multi-wavelength photometry assembled by the 3D-HST survey. The combined dataset is modeled to infer the past star-formation histories (SFHs) of both merging galaxies. They are found to… ▽ More We present a detailed study of a galaxy merger taking place at $z=1.89$ in the GOODS-S field. Here we analyze Keck/MOSFIRE spectroscopic observations from the MOSFIRE Deep Evolution Field (MOSDEF) survey along with multi-wavelength photometry assembled by the 3D-HST survey. The combined dataset is modeled to infer the past star-formation histories (SFHs) of both merging galaxies. They are found to be massive, with log$_{10}(M_{\ast}/M_{\odot}) > 11$, with a close mass ratio satisfying the typical major-merger definition. Additionally, in the context of delayed-$τ$ models, GOODS-S 43114 and GOODS-S 43683 have similar SFHs and low star-formation rates (log$_{10}$(SFR(SED)/$M_{\odot}/\rm{yr}^{-1}$) $<$ 1.0) compared to their past averages. The best-fit model SEDs show elevated H$δ_{\rm{A}}$ values for both galaxies, indicating that their stellar spectra are dominated by A-type stars, and that star formation peaked $\sim0.5-1$ Gyr ago and has recently declined. Additionally, based on SED fitting both merging galaxies turned on and shut off star formation within a few hundred Myr of each other, suggesting that their bursts of star formation may be linked. Combining the SFHs and H$δ_{\rm{A}}$ results with recent galaxy merger simulations, we infer that these galaxies have recently completed their first pericentric passage and are moving apart. Finally, the relatively low second velocity moment of GOODS-S 43114 given its stellar mass, suggests a disk-like structure. However, including the geometry of the galaxy in the modeling does not completely resolve the discrepancy between the dynamical and stellar masses. Future work is needed to resolve this inconsistency in mass. △ Less

Submitted 6 October, 2022; v1 submitted 10 June, 2022; originally announced June 2022.

Comments: 12 pages, 7 figures, 1 table, Accepted for publication in MNRAS

AGN accretion and black hole growth across compact and extended galaxy evolution phases

Authors: James Aird, Alison L. Coil, Dale D. Kocevski

Abstract: The extent of black hole growth during different galaxy evolution phases and the connection between galaxy compactness and AGN activity remain poorly understood. We use Hubble Space Telescope imaging of the CANDELS fields to identify star-forming and quiescent galaxies at z=0.5-3 in both compact and extended phases and use Chandra X-ray imaging to measure the distribution of AGN accretion rates an… ▽ More The extent of black hole growth during different galaxy evolution phases and the connection between galaxy compactness and AGN activity remain poorly understood. We use Hubble Space Telescope imaging of the CANDELS fields to identify star-forming and quiescent galaxies at z=0.5-3 in both compact and extended phases and use Chandra X-ray imaging to measure the distribution of AGN accretion rates and track black hole growth within these galaxies. Accounting for the impact of AGN light changes ~20% of the X-ray sources from compact to extended galaxy classifications. We find that ~10-25% of compact star-forming galaxies host an AGN, a mild enhancement (by a factor ~2) compared to extended star-forming galaxies or compact quiescent galaxies of equivalent stellar mass and redshift. However, AGN are not ubiquitous in compact star-forming galaxies and this is not the evolutionary phase, given its relatively short timescale, where the bulk of black hole mass growth takes place. Conversely, we measure the highest AGN fractions (~10-30%) within the relatively rare population of extended quiescent galaxies. For massive galaxies that quench at early cosmic epochs, substantial black hole growth in this extended phase is crucial to produce the elevated black hole mass-to-galaxy stellar mass scaling relation observed for quiescent galaxies at z~0. We also show that AGN fraction increases with compactness in star-forming galaxies and decreases in quiescent galaxies within both the compact and extended sub-populations, demonstrating that AGN activity depends closely on the structural properties of galaxies. △ Less

Submitted 25 July, 2022; v1 submitted 27 January, 2022; originally announced January 2022.

Comments: 30 pages, 18 figures, accepted for publication in MNRAS. Primary results are shown in Fig 7 and summarised by Fig 11. See Fig 15 and 16 for key interpretation/conclusions

Constraints on the X-ray Luminosity Function of AGN at z=5.7-6.4 with the Extragalactic Serendipitous Swift Survey

Authors: Cassandra L. Barlow-Hall, Jack Delaney, James Aird, Philip A. Evans, Julian P. Osborne, Michael G. Watson

Abstract: X-ray luminosity functions (XLFs) of Active Galactic Nuclei (AGN) 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 \redit{place limits} on the bright-end of the XLF at z=5.7-6.4 using high-redshift AGN ide… ▽ More X-ray luminosity functions (XLFs) of Active Galactic Nuclei (AGN) 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 \redit{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}\times10^{44}\mathrm{erg.s^{-1}}$ and $z=6.31\pm0.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 toward 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\gtrsim6$, although they require a large X-ray to bolometric correction factor at these high luminosities. △ Less

Submitted 15 December, 2022; v1 submitted 26 January, 2022; originally announced January 2022.

Comments: 10pages, 6 figures. Resubmitted to MNRAS, following referee comments

The eROSITA Final Equatorial Depth Survey (eFEDS): X-ray emission around star-forming and quiescent galaxies at $0.05<z<0.3$

Authors: Johan Comparat, Nhut Truong, Andrea Merloni, Annalisa Pillepich, Gabriele Ponti, Simon Driver, Sabine Bellstedt, Joe Liske, James Aird, Marcus Brüggen, Esra Bulbul, Luke Davies, Justo Antonio González Villalba, Antonis Georgakakis, Frank Haberl, Teng Liu, Chandreyee Maitra, Kirpal Nandra, Paola Popesso, Peter Predehl, Aaron Robotham, Mara Salvato, Jessica E. Thorne, Yi Zhang

Abstract: We aim at characterizing the hot phase of the Circum-Galactic Medium in a large sample of galaxies. We stack X-ray events from the SRG/eROSITA eFEDS survey around central galaxies in the GAMA 9hr field to construct radially projected soft X-ray luminosity profiles as a function of their stellar mass and specific star formation rate. We consider samples of quiescent (star-forming) galaxies in the s… ▽ More We aim at characterizing the hot phase of the Circum-Galactic Medium in a large sample of galaxies. We stack X-ray events from the SRG/eROSITA eFEDS survey around central galaxies in the GAMA 9hr field to construct radially projected soft X-ray luminosity profiles as a function of their stellar mass and specific star formation rate. We consider samples of quiescent (star-forming) galaxies in the stellar mass range $2\times 10^{10}$ -- $10^{12}$ M$_\odot$ ($3\times 10^9$ -- $6\times10^{11}$ M$_\odot$). For quiescent galaxies, the X-ray profiles are clearly extended throughout the available mass range; however, the measured profile is likely biased high due to projection effects, as these galaxies tend to live in dense and hot environments. For the most massive star forming samples ($\geq10^{11}$ M$_\odot$), there is a hint of detection of extended emission. For star-forming galaxies with $< 10^{11}$ M$_\odot$ the X-ray stacked profiles are compatible with unresolved sources and consistent with the expected emission from faint Active Galactic Nuclei and X-ray binaries. We measure for the first time the mean relation between average X-ray luminosity and stellar mass separately for quiescent and star-forming galaxies. High-mass ($\geq 10^{11}$ M$_\odot$) star-forming or quiescent galaxies follow the expected scaling of virialized hot haloes, while lower mass star-forming galaxies show a less prominent luminosity and a weaker dependence on stellar mass, consistent with empirical models of the weak AGN population. When comparing our results with state-of-the art numerical simulations, we find an overall consistency on large ($>80$ kpc) scales at masses $\geq 10^{11}$ M$_\odot$, but disagreement on the small scales, where brighter than observed compact cores are predicted. Simulations also do not predict the clear differentiation that we observe between quiescent and star-forming galaxies. △ Less

Submitted 17 August, 2022; v1 submitted 13 January, 2022; originally announced January 2022.

Comments: 23 pages, 11 figures, 4 tables, accepted in A&A

Journal ref: A&A 666, A156 (2022)

The Incidence of X-ray selected AGN in Nearby Galaxies

Authors: Keir L. Birchall, M. G. Watson, J. Aird, R. L. C. Starling

Abstract: We present the identification and analysis of an unbiased sample of 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 \leq 0.33$). After confirming that there was strictly no AGN light contaminating stellar mass and star-formation rate calculations, we identified 917 galaxies with cen… ▽ More We present the identification and analysis of an unbiased sample of 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 \leq 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% agreement in the lowest mass bin to 93% in the highest. We then calculated the growth rates of the black holes powering these AGN in terms of their specific accretion rates ($\propto L_X/M_*$). Our sample exhibits a wide range of accretion rates, with the majority accreting at rates $\leq 0.5\%$ 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\%$) over stellar masses of $8 < \log \mathrm{M_*/M_\odot} < 12$ and increases (from $\approx 1\%$ to $10\%$) with redshift. △ Less

Submitted 14 December, 2021; v1 submitted 6 December, 2021; originally announced December 2021.

Comments: 18 pages, 10 figures, 2 appendices. Accepted for publication in MNRAS

Supermassive black holes in cosmological simulations II: the AGN population and predictions for upcoming X-ray missions

Authors: Melanie Habouzit, Rachel S. Somerville, Yuan Li, Shy Genel, James Aird, Daniel Anglés-Alcázar, Romeel Davé, Iskren Y. Georgiev, Stuart McAlpine, Yetli Rosas-Guevara, Yohan Dubois, Dylan Nelson, Eduardo Bañados, Lars Hernquist, Sébastien Peirani, Mark Vogelsberger

Abstract: In large-scale hydrodynamical cosmological simulations, the fate of massive galaxies is mainly dictated by the modeling of feedback from active galactic nuclei (AGN). The amount of energy released by AGN feedback is proportional to the mass that has been accreted onto the BHs, but the exact sub-grid modeling of AGN feedback differs in all simulations. Whilst modern simulations reliably produce pop… ▽ More In large-scale hydrodynamical cosmological simulations, the fate of massive galaxies is mainly dictated by the modeling of feedback from active galactic nuclei (AGN). The amount of energy released by AGN feedback is proportional to the mass that has been accreted onto the BHs, but the exact sub-grid modeling of AGN feedback differs in all simulations. Whilst modern simulations reliably produce populations of quiescent massive galaxies at z<2, it is also crucial to assess the similarities and differences of the responsible AGN populations. Here, we compare the AGN population of the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA simulations. The AGN luminosity function (LF) varies significantly between simulations. Although in agreement with current observational constraints at z=0, at higher redshift the agreement of the LFs deteriorates with most simulations producing too many AGN of L_{x, 2-10 keV}~10^43-10^44 erg/s. AGN feedback in some simulations prevents the existence of any bright AGN with L_{x, 2-10 keV}>=10^45 erg/s (although this is sensitive to AGN variability), and leads to smaller fractions of AGN in massive galaxies than in the observations at z<=2. We find that all the simulations fail at producing a number density of AGN in good agreement with observational constraints for both luminous (L_{x, 2-10 keV}~10^43-10^45 erg/s) and fainter (L_{x, 2-10 keV}~10^42-10^43 erg/s) AGN, and at both low and high redshift. These differences can aid us in improving future BH and galaxy subgrid modeling in simulations. Upcoming X-ray missions (e.g., Athena, AXIS, and LynX) will bring faint AGN to light and new powerful constraints. After accounting for AGN obscuration, we find that the predicted number density of detectable AGN in future surveys spans at least one order of magnitude across the simulations, at any redshift. △ Less

Submitted 2 November, 2021; originally announced November 2021.

Comments: Published in MNRAS

Multi-messenger-Athena Synergy White Paper

Authors: L. Piro, M. Ahlers, A. Coleiro, M. Colpi, E. de Oña Wilhelmi, M. Guainazzi, P. G. Jonker, P. Mc Namara, D. A. Nichols, P. O'Brien, E. Troja, J. Vink, J. Aird, L. Amati, S. Anand, E. Bozzo, F. J. Carrera, A. C. Fabian, C. Fryer, E. Hall, O. Korobkin, V. Korol, A. Mangiagli, S. Martínez-Núñez, S. Nissanke , et al. (8 additional authors not shown)

Abstract: In this paper we explore the scientific synergies between Athena and some of the key multi-messenger facilities that should be operative concurrently with Athena. These facilities include LIGO A+, Advanced Virgo+ and future detectors for ground-based observation of gravitational waves (GW), LISA for space-based observations of GW, IceCube and KM3NeT for neutrino observations, and CTA for very high… ▽ More In this paper we explore the scientific synergies between Athena and some of the key multi-messenger facilities that should be operative concurrently with Athena. These facilities include LIGO A+, Advanced Virgo+ and future detectors for ground-based observation of gravitational waves (GW), LISA for space-based observations of GW, IceCube and KM3NeT for neutrino observations, and CTA for very high energy observations. These science themes encompass pressing issues in astrophysics, cosmology and fundamental physics such as: the central engine and jet physics in compact binary mergers, accretion processes and jet physics in Super-Massive Binary Black Holes (SMBBHs) and in compact stellar binaries, the equation of state of neutron stars, cosmic accelerators and the origin of Cosmic Rays (CRs), the origin of intermediate and high-Z elements in the Universe, the Cosmic distance scale and tests of General Relativity and the Standard Model. Observational strategies for implementing the identified science topics are also discussed. A significant part of the sources targeted by multi-messenger facilities is of transient nature. We have thus also discussed the synergy of \textsl{Athena} with wide-field high-energy facilities, taking THESEUS as a case study for transient discovery. This discussion covers all the Athena science goals that rely on follow-up observations of high-energy transients identified by external observatories, and includes also topics that are not based on multi-messenger observations, such as the search for missing baryons or the observation of early star populations and metal enrichment at the cosmic dawn with Gamma-Ray Bursts (GRBs). △ Less

Submitted 29 October, 2021; originally announced October 2021.

Comments: White paper of the multi-messenger Athena synergy excersise

The connection between star formation and supermassive Black Hole activity in the local Universe

Authors: Olena Torbaniuk, Maurizio Paolillo, Francisco Carrera, Stefano Cavuoti, Cristian Vignali, Giuseppe Longo, James Aird

Abstract: We present a study of the active galactic nucleus (AGN) activity in the local Universe (z < 0.33) and its correlation with the host galaxy properties, derived from a Sloan Digital Sky Survey (SDSS DR8) sample with spectroscopic star-formation rate (SFR) and stellar mass ($\mathcal{M}_{\ast}$) determination. To quantify the level of AGN activity we used X-ray information from the XMM-Newton Serendi… ▽ More We present a study of the active galactic nucleus (AGN) activity in the local Universe (z < 0.33) and its correlation with the host galaxy properties, derived from a Sloan Digital Sky Survey (SDSS DR8) sample with spectroscopic star-formation rate (SFR) and stellar mass ($\mathcal{M}_{\ast}$) determination. To quantify the level of AGN activity we used X-ray information from the XMM-Newton Serendipitous Source Catalogue (3XMM DR8). Applying multiwavelength AGN selection criteria (optical BPT-diagrams, X-ray/optical ratio etc) we found that 24% of the detected sources are efficiently-accreting AGN with moderate-to-high X-ray luminosity, which are twice as likely to be hosted by star-forming galaxies than by quiescent ones. The distribution of the specific Black Hole accretion rate (sBHAR, $λ_{\mathrm{sBHAR}}$) shows that nuclear activity in local, non-AGN dominated galaxies peaks at very low accretion rates ($-4 \lesssim \logλ_{\mathrm{sBHAR}} \lesssim -3$) in all stellar mass ranges. However, we observe systematically larger values of sBHAR for galaxies with active star-formation than for quiescent ones, as well as an increase of the mean $λ_{\mathrm{sBHAR}}$ with SFR for both star-forming and quiescent galaxies. These findings confirm the decreased level of AGN activity with cosmic time and are consistent with a scenario where both star-formation and AGN activity are fuelled by a common gas reservoir. △ Less

Submitted 21 June, 2021; originally announced June 2021.

Comments: 20 pages, 20 figures, accepted for publication in MNRAS

Journal ref: Volume 506 (2021), pp. 2619-2637

The XMM-SERVS survey: XMM-Newton point-source catalogs for the W-CDF-S and ELAIS-S1 fields

Authors: Q. Ni, W. N. Brandt, C. -T. Chen, B. Luo, K. Nyland, G. Yang, F. Zou, J. Aird, D. M. Alexander, F. E. Bauer, M. Lacy, B. D. Lehmer, L. Mallick, M. Salvato, D. P. Schneider, P. Tozzi, I. Traulsen, M. Vaccari, C. Vignali, F. Vito, Y. Xue, M. Banerji, K. Chow, A. Comastri, A. Del Moro , et al. (8 additional authors not shown)

Abstract: We present the X-ray point-source catalogs in two of the XMM-Spitzer Extragalactic Representative Volume Survey (XMM-SERVS) fields, W-CDF-S (4.6 deg$^2$) and ELAIS-S1 (3.2 deg$^2$), aiming to fill the gap between deep pencil-beam X-ray surveys and shallow X-ray surveys over large areas. The W-CDF-S and ELAIS-S1 regions were targeted with 2.3 Ms and 1.0 Ms of XMM-Newton observations, respectively;… ▽ More We present the X-ray point-source catalogs in two of the XMM-Spitzer Extragalactic Representative Volume Survey (XMM-SERVS) fields, W-CDF-S (4.6 deg$^2$) and ELAIS-S1 (3.2 deg$^2$), aiming to fill the gap between deep pencil-beam X-ray surveys and shallow X-ray surveys over large areas. The W-CDF-S and ELAIS-S1 regions were targeted with 2.3 Ms and 1.0 Ms of XMM-Newton observations, respectively; 1.8 Ms and 0.9 Ms exposures remain after flare filtering. The survey in W-CDF-S has a flux limit of 1.0 $\times$ 10$^{-14}$ erg cm$^{-2}$ s$^{-1}$ over 90% of its area in the 0.5-10 keV band; 4053 sources are detected in total. The survey in ELAIS-S1 has a flux limit of 1.3 $\times$ 10$^{-14}$ erg cm$^{-2}$ s$^{-1}$ over 90% of its area in the 0.5-10 keV band; 2630 sources are detected in total. Reliable optical-to-IR multiwavelength counterpart candidates are identified for $\approx$ 89% of the sources in W-CDF-S and $\approx$ 87% of the sources in ELAIS-S1. 3186 sources in W-CDF-S and 1985 sources in ELAIS-S1 are classified as AGNs. We also provide photometric redshifts for X-ray sources; $\approx$ 84% of the 3319/2001 sources in W-CDF-S/ELAIS-S1 with optical-to-NIR forced photometry available have either spectroscopic redshifts or high-quality photometric redshifts. The completion of the XMM-Newton observations in the W-CDF-S and ELAIS-S1 fields marks the end of the XMM-SERVS survey data gathering. The $\approx$ 12,000 point-like X-ray sources detected in the whole $\approx$ 13 deg$^2$ XMM-SERVS survey will benefit future large-sample AGN studies. △ Less

Submitted 19 June, 2021; originally announced June 2021.

Comments: 40 pages, 35 figures. Accepted for publication in ApJS. Data products available at: https://personal.psu.edu/wnb3/xmmservs/xmmservs.html

Trinity I: Self-Consistently Modeling the Dark Matter Halo-Galaxy-Supermassive Black Hole Connection from $z=0-10$

Authors: Haowen Zhang, Peter Behroozi, Marta Volonteri, Joseph Silk, Xiaohui Fan, Philip F. Hopkins, Jinyi Yang, James Aird

Abstract: We present Trinity, a flexible empirical model that self-consistently infers the statistical connection between dark matter haloes, galaxies, and supermassive black holes (SMBHs). Trinity is constrained by galaxy observables from $0<z<10$ (galaxies' stellar mass functions, specific and cosmic SFRs, quenched fractions, and UV luminosity functions) and SMBH observables from $0<z<6.5$ (quasar luminos… ▽ More We present Trinity, a flexible empirical model that self-consistently infers the statistical connection between dark matter haloes, galaxies, and supermassive black holes (SMBHs). Trinity is constrained by galaxy observables from $0<z<10$ (galaxies' stellar mass functions, specific and cosmic SFRs, quenched fractions, and UV luminosity functions) and SMBH observables from $0<z<6.5$ (quasar luminosity functions, quasar probability distribution functions, active black hole mass functions, local SMBH mass--bulge mass relations, and the observed SMBH mass distributions of high redshift bright quasars). The model includes full treatment of observational systematics (e.g., AGN obscuration and errors in stellar masses). From these data, Trinity infers the average SMBH mass, SMBH accretion rate, merger rate, and Eddington ratio distribution as functions of halo mass, galaxy stellar mass, and redshift. Key findings include: 1) the normalization and the slope of the SMBH mass--bulge mass relation increases mildly from $z=0$ to $z=10$; 2) The best-fitting AGN radiative$+$kinetic efficiency is $\sim 0.05-0.06$, but can range from $\sim 0.035-0.07$ with alternative input assumptions; 3) AGNs show downsizing, i.e., the Eddington ratios of more massive SMBHs start to decrease earlier than those of lower-mass objects; 4) The average ratio between average SMBH accretion rate and SFR is $\sim 10^{-3}$ for low-mass galaxies, which are primarily star-forming. This ratio increases to $\sim 10^{-1}$ for the most massive haloes below $z\sim 1$, where star formation is quenched but SMBHs continue to accrete. △ Less

Submitted 9 September, 2022; v1 submitted 21 May, 2021; originally announced May 2021.

Comments: Accepted by MNRAS. 42 pages, 27 figures. Trinity code available at: https://github.com/HaowenZhang/TRINITY Comments welcome!

The AGN-galaxy-halo connection: The distribution of AGN host halo masses to z=2.5

Authors: James Aird, Alison L. Coil

Abstract: It is widely reported, based on clustering measurements of observed active galactic nuclei (AGN) samples, that AGN reside in similar mass host dark matter halos across the bulk of cosmic time, with log $M/M_\odot$~12.5-13.0 to z~2.5. We show that this is due in part to the AGN fraction in galaxies rising with increasing stellar mass, combined with AGN observational selection effects that exacerbat… ▽ More It is widely reported, based on clustering measurements of observed active galactic nuclei (AGN) samples, that AGN reside in similar mass host dark matter halos across the bulk of cosmic time, with log $M/M_\odot$~12.5-13.0 to z~2.5. We show that this is due in part to the AGN fraction in galaxies rising with increasing stellar mass, combined with AGN observational selection effects that exacerbate this trend. Here, we use AGN specific accretion rate distribution functions determined as a function of stellar mass and redshift for star-forming and quiescent galaxies separately, combined with the latest galaxy-halo connection models, to determine the parent and sub-halo mass distribution function of AGN to various observational limits. We find that while the median (sub-)halo mass of AGN, $\approx10^{12}M_\odot$, is fairly constant with luminosity, specific accretion rate, and redshift, the full halo mass distribution function is broad, spanning several orders of magnitude. We show that widely used methods to infer a typical dark matter halo mass based on an observed AGN clustering amplitude can result in biased, systematically high host halo masses. While the AGN satellite fraction rises with increasing parent halo mass, we find that the central galaxy is often not an AGN. Our results elucidate the physical causes for the apparent uniformity of AGN host halos across cosmic time and underscore the importance of accounting for AGN selection biases when interpreting observational AGN clustering results. We further show that AGN clustering is most easily interpreted in terms of the relative bias to galaxy samples, not from absolute bias measurements alone. △ Less

Submitted 2 February, 2021; v1 submitted 6 October, 2020; originally announced October 2020.

Comments: 19 pages, 20 figures, accepted for publication in MNRAS. Updated to show median halo masses following referee's helpful comments

The Chandra Deep Wide-Field Survey: A New Chandra Legacy Survey in the Boötes Field I. X-ray Point Source Catalog, Number Counts and Multi-Wavelength Counterparts

Authors: A. Masini, R. C. Hickox, C. M. Carroll, J. Aird, D. M. Alexander, R. J. Assef, R. Bower, M. Brodwin, M. J. I. Brown, S. Chatterjee, C. -T. J. Chen, A. Dey, M. A. DiPompeo, K. J. Duncan, P. R. M. Eisenhardt, W. R. Forman, A. H. Gonzalez, A. D. Goulding, K. N. Hainline, B. T. Jannuzi, C. Jones, C. S. Kochanek, R. Kraft, K. -S. Lee, E. D. Miller , et al. (8 additional authors not shown)

Abstract: We present a new, ambitious survey performed with the Chandra X-ray Observatory of the 9.3 deg$^2$ Boötes field of the NOAO Deep Wide-Field Survey. The wide field probes a statistically representative volume of the Universe at high redshift. The Chandra Deep Wide-Field Survey exploits the excellent sensitivity and angular resolution of Chandra over a wide area, combining 281 observations spanning… ▽ More We present a new, ambitious survey performed with the Chandra X-ray Observatory of the 9.3 deg$^2$ Boötes field of the NOAO Deep Wide-Field Survey. The wide field probes a statistically representative volume of the Universe at high redshift. The Chandra Deep Wide-Field Survey exploits the excellent sensitivity and angular resolution of Chandra over a wide area, combining 281 observations spanning 15 years, for a total exposure time of 3.4 Ms, and detects 6891 X-ray point sources down to limiting fluxes of $4.7\times10^{-16}$, $1.5\times10^{-16}$, and $9\times10^{-16}$ erg cm$^{-2}$ s$^{-1}$, in the $0.5-7$ keV, $0.5-2$ keV, and $2-7$ keV bands, respectively. The robustness and reliability of the detection strategy is validated through extensive, state-of-the-art simulations of the whole field. Accurate number counts, in good agreement with previous X-ray surveys, are derived thanks to the uniquely large number of point sources detected, which resolve $65.0 \pm 12.8\%$ of the cosmic X-ray background between $0.5-2$ keV and $81.0 \pm 11.5\%$ between $2-7$ keV. Exploiting the wealth of multi-wavelength data available on the field, we assign redshifts to $\sim 94\%$ of the X-ray sources, estimate their obscuration and derive absorption-corrected luminosities. We provide an electronic catalog containing all the relevant quantities needed for future investigations. △ Less

Submitted 7 September, 2020; originally announced September 2020.

Comments: Accepted for publication in ApJS. The catalogs associated to this paper can be accessed at http://www.dartmouth.edu/~hickox/cdwfs.php

Local AGN Survey (LASr): I. Galaxy sample, infrared colour selection and predictions for AGN within 100 Mpc

Authors: D. Asmus, C. L. Greenwell, P. Gandhi, P. G. Boorman, J. Aird, D. M. Alexander, R. J. Assef, R. D. Baldi, R. I. Davies, S. F. Hönig, C. Ricci, D. J. Rosario, M. Salvato, F. Shankar, D. Stern

Abstract: In order to answer some of the major open questions in the fields of supermassive black hole (SMBH) and galaxy evolution, a complete census of SMBH growth, i.e., active galactic nuclei (AGN), is required. Thanks to deep all-sky surveys, such as those by the Wide-field Infrared Survey Explorer (WISE) and the Spectrum-Roentgen-Gamma (SRG) missions, this task is now becoming feasible in the nearby Un… ▽ More In order to answer some of the major open questions in the fields of supermassive black hole (SMBH) and galaxy evolution, a complete census of SMBH growth, i.e., active galactic nuclei (AGN), is required. Thanks to deep all-sky surveys, such as those by the Wide-field Infrared Survey Explorer (WISE) and the Spectrum-Roentgen-Gamma (SRG) missions, this task is now becoming feasible in the nearby Universe. Here, we present a new survey, the Local AGN Survey (LASr), with the goal of identifying AGN unbiased against obscuration and determining the intrinsic Compton-thick (CT) fraction. First, we construct the most complete all-sky sample of galaxies within 100 Mpc from astronomical databases (90% completeness for log M_star/M_sun ~9.4), four times deeper than the current local galaxy reference, the Two Micron All-Sky Survey Redshift Survey (2MRS), which turns out to miss 20% of known luminous AGN. These 49k galaxies serve as parent sample for LASr, called LASr-GPS. It contains 4.3k already known AGN, >=82% of these are estimated to have L(12um) < 10^42.3 erg/s, i.e., are low-luminosity AGN. As a first method for identifying Seyfert-like AGN, we use WISE-based infrared colours, finding 221 galaxies at L(12um) >= 10^42.3 erg/s to host an AGN at 90% reliability, This includes 61 new AGN candidates and implies and optical type 2 fraction of 50 to 71%. We quantify the efficiency of this technique and estimate the total number of AGN with L_int(2-10 keV) >= 10^42 erg/s in the volume to be 362^+145_-116 (8.6^+3.5_-2.8 10^-5 Mpc^-3). X-ray brightness estimates indicate the CT fraction to be 40-55% to explain the Swift non-detections of the infrared selected objects. One third of the AGN within 100 Mpc remain to be identified and we discuss the prospects for the eROSITA all-sky survey. △ Less

Submitted 12 March, 2020; originally announced March 2020.

Comments: Accepted for publication in MNRAS. 30 pages, 19 figures. For full, digital versions of Table 1, 2 and 3, see: https://my.pcloud.com/publink/show?code=kZiqWnkZgyLixKFSgruXOsYYglos20Y3kB6V

The evolving AGN duty cycle in galaxies since z$\sim$3 as encoded in the X-ray luminosity function

Authors: I. Delvecchio, E. Daddi, J. Aird, J. R. Mullaney, E. Bernhard, L. P. Grimmett, R. Carraro, A. Cimatti, G. Zamorani, N. Caplar, F. Vito, D. Elbaz, G. Rodighiero

Abstract: We present a new modeling of the X-ray luminosity function (XLF) of Active Galactic Nuclei (AGN) out to z$\sim$3, dissecting the contribution of main-sequence (MS) and starburst (SB) galaxies. For each galaxy population, we convolved the observed galaxy stellar mass (M$_{\star}$) function with a grid of M$_{\star}$-independent Eddington ratio ($λ_{\rm EDD}$) distributions, normalised via empirical… ▽ More We present a new modeling of the X-ray luminosity function (XLF) of Active Galactic Nuclei (AGN) out to z$\sim$3, dissecting the contribution of main-sequence (MS) and starburst (SB) galaxies. For each galaxy population, we convolved the observed galaxy stellar mass (M$_{\star}$) function with a grid of M$_{\star}$-independent Eddington ratio ($λ_{\rm EDD}$) distributions, normalised via empirical black hole accretion rate (BHAR) to star formation rate (SFR) relations. Our simple approach yields an excellent agreement with the observed XLF since z$\sim$3. We find that the redshift evolution of the observed XLF can only be reproduced through an intrinsic flattening of the $λ_{\rm EDD}$ distribution, and with a positive shift of the break $λ^{*}$, consistent with an anti-hierarchical behavior. The AGN accretion history is predominantly made by massive (10$^{10}<$M$_{\star}<$10$^{11}$ M$_{\odot}$) MS galaxies, while SB-driven BH accretion, possibly associated with galaxy mergers, becomes dominant only in bright quasars, at $\log$(L$_{\rm X}$/erg s$^{-1}$)$>$44.36 + 1.28$\cdot$(1+z). We infer that the probability of finding highly-accreting ($λ_{\rm EDD}>$ 10%) AGN significantly increases with redshift, from 0.4% (3.0%) at z=0.5 to 6.5% (15.3%) at z=3 for MS (SB) galaxies, implying a longer AGN duty cycle in the early Universe. Our results strongly favor a M$_{\star}$-dependent ratio between BHAR and SFR, as BHAR/SFR $\propto$ M$_{\star}^{0.73[+0.22,-0.29]}$, supporting a non-linear BH buildup relative to the host. Finally, this framework opens potential questions on super-Eddington BH accretion and different $λ_{\rm EDD}$ prescriptions for understanding the cosmic BH mass assembly. △ Less

Submitted 20 February, 2020; originally announced February 2020.

Comments: Accepted for publication in ApJ. 25 pages, 15 figures, 2 tables. The best X-ray luminosity functions are available at https://github.com/idelvecchio/XLF_Delvecchio2020.git

X-ray Detected AGN in SDSS Dwarf Galaxies

Authors: Keir L Birchall, M G Watson, J Aird

Abstract: In this work we present a robust quantification of X-ray selected AGN in local ($z \leq 0.25$) dwarf galaxies ($M_\mathrm{*} \leq 3 \times 10^9 \mathrm{M_\odot}$). We define a parent sample of 4,331 dwarf galaxies found within the footprint of both the MPA-JHU galaxy catalogue (based on SDSS DR8) and 3XMM DR7, performed a careful review of the data to remove misidentifications and produced a sampl… ▽ More In this work we present a robust quantification of X-ray selected AGN in local ($z \leq 0.25$) dwarf galaxies ($M_\mathrm{*} \leq 3 \times 10^9 \mathrm{M_\odot}$). We define a parent sample of 4,331 dwarf galaxies found within the footprint of both the MPA-JHU galaxy catalogue (based on SDSS DR8) and 3XMM DR7, performed a careful review of the data to remove misidentifications and produced a sample of 61 dwarf galaxies that exhibit nuclear X-ray activity indicative of an AGN. We examine the optical emission line ratios of our X-ray selected sample and find that optical AGN diagnostics fail to identify 85% of the sources. We then calculated the growth rates of the black holes powering our AGN in terms of their specific accretion rates ($\propto L_\mathrm{X}/M_\mathrm{*}$, an approximate tracer of the Eddington ratio). Within our observed sample, we found a wide range of specific accretion rates. After correcting the observed sample for the varying sensitivity of 3XMM, we found further evidence for a wide range of X-ray luminosities and specific accretion rates, described by a power law. Using this corrected AGN sample we also define an AGN fraction describing their relative incidence within the parent sample. We found the AGN fraction increases with host galaxy mass (up to $\approx$ 6%) for galaxies with X-ray luminosities between $10^{39} \mathrm{erg/s}$ and $10^{42} \mathrm{erg/s}$, and by extrapolating the power law to higher luminosities, we found evidence to suggest the fraction of luminous AGN ($L_\mathrm{X} \geq 10^{42.4} \mathrm{erg/s}$) is constant out to $z \approx 0.7$. △ Less

Submitted 28 January, 2020; v1 submitted 9 January, 2020; originally announced January 2020.

Comments: 19 pages, 12 figures, 1 table. Accepted for publication in MNRAS

Journal ref: MNRAS, Volume 492, Issue 2, February 2020, Pages 2268 to 2284

The galaxy's gas content regulated by the dark matter halo mass results in a super-linear M$_{\rm BH}$-M$_{\star}$ relation

Authors: I. Delvecchio, E. Daddi, F. Shankar, J. R. Mullaney, G. Zamorani, J. Aird, E. Bernhard, A. Cimatti, D. Elbaz, M. Giavalisco, L. P. Grimmett

Abstract: Supermassive black holes (SMBHs) are tightly correlated with their hosts but the origin of such connection remains elusive. To explore the cosmic build-up of this scaling relation, we present an empirically-motivated model that tracks galaxy and SMBH growth down to z=0. Starting from a random mass seed distribution at z=10, we assume that each galaxy evolves on the star-forming "main sequence" (MS… ▽ More Supermassive black holes (SMBHs) are tightly correlated with their hosts but the origin of such connection remains elusive. To explore the cosmic build-up of this scaling relation, we present an empirically-motivated model that tracks galaxy and SMBH growth down to z=0. Starting from a random mass seed distribution at z=10, we assume that each galaxy evolves on the star-forming "main sequence" (MS) and each BH follows the recently-derived stellar mass (M$_{\star}$) dependent ratio between BH accretion rate and star formation rate, going as BHAR/SFR$\propto$M$_{\star}^{0.73[+0.22,-0.29]}$. Our simple recipe naturally describes the BH-galaxy build-up in two stages. At first, the SMBH lags behind the host that evolves along the MS. Later, as the galaxy grows in M$_{\star}$, our M$_{\star}$-dependent BHAR/SFR induces a super-linear BH growth, as M$_{\rm BH}$$\propto$M$_{\star}^{1.7}$. According to this formalism, smaller BH seeds increase their relative mass faster and earlier than bigger BH seeds, at fixed M$_{\star}$, thus setting along a gradually tighter M$_{\rm BH}$-M$_{\star}$ locus towards higher M$_{\star}$. Assuming reasonable values of the radiative efficiency $ε\sim$0.1, our empirical trend agrees with both high-redshift model predictions and intrinsic M$_{\rm BH}$-M$_{\star}$ relations of local BHs. We speculate that the observed non-linear BH-galaxy build-up is reflected in a twofold behavior with dark matter halo mass (M$_{\rm DM}$), displaying a clear turnover at M$_{\rm DM}\sim$2$\times$10$^{12}$M$_{\odot}$. While Supernovae-driven feedback suppresses BH growth in smaller halos (BHAR/SFR$\propto$M$_{\rm DM}^{1.6}$), above the M$_{\rm DM}$ threshold cold gas inflows possibly fuel both BH accretion and star formation in a similar fashion (BHAR/SFR$\propto$M$_{\rm DM}^{0.3}$). △ Less

Submitted 17 October, 2019; originally announced October 2019.

Comments: Accepted for publication in ApJ Letter. 8 pages, 5 figures

The MOSDEF Survey: The Metallicity Dependence of X-ray Binary Populations at $z\sim2$

Authors: Francesca M. Fornasini, Mariska Kriek, Ryan L. Sanders, Irene Shivaei, Francesca Civano, Naveen A. Reddy, Alice E. Shapley, Alison L. Coil, Bahram Mobasher, Brian Siana, James Aird, Mojegan Azadi, William R. Freeman, Gene C. K. Leung, Sedona H. Price, Tara Fetherolf, Tom Zick, Guillermo Barro

Abstract: Population synthesis models predict that high-mass X-ray binary (HMXB) populations produced in low metallicity environments should be more X-ray luminous, a trend supported by studies of nearby galaxies. This trend may be responsible for the observed increase of the X-ray luminosity ($L_{\mathrm{X}}$) per star formation rate (SFR) with redshift due to the decrease of metallicity ($Z$) at fixed ste… ▽ More Population synthesis models predict that high-mass X-ray binary (HMXB) populations produced in low metallicity environments should be more X-ray luminous, a trend supported by studies of nearby galaxies. This trend may be responsible for the observed increase of the X-ray luminosity ($L_{\mathrm{X}}$) per star formation rate (SFR) with redshift due to the decrease of metallicity ($Z$) at fixed stellar mass as a function of redshift. To test this hypothesis, we use a sample of 79 $z\sim2$ star-forming galaxies with oxygen abundance measurements from the MOSDEF survey, which obtained rest-frame optical spectra for $\sim1500$ galaxies in the CANDELS fields at $1.37<z<3.80$. Using Chandra data from the AEGIS-X Deep, Deep Field North, and Deep Field South surveys, we stack the X-ray data at the galaxy locations in bins of redshift and $Z$ because the galaxies are too faint to be individually detected. In agreement with previous studies, the average $L_{\mathrm{X}}$/SFR of our $z\sim2$ galaxy sample is enhanced by $\approx0.4-0.8$ dex relative to local HMXB $L_{\mathrm{X}}$-SFR scaling relations. Splitting our sample by $Z$, we find that $L_{\mathrm{X}}$/SFR and $Z$ are anti-correlated with 97% confidence. This observed $Z$ dependence for HMXB-dominated galaxies is consistent both with the local $L_{\mathrm{X}}$-SFR-$Z$ relation and a subset of population synthesis models. Although the statistical significance of the observed trends is weak due to the low X-ray statistics, these results constitute the first direct evidence connecting the redshift evolution of $L_{\mathrm{X}}$/SFR and the $Z$ dependence of HMXBs. △ Less

Submitted 18 September, 2019; originally announced September 2019.

Comments: 22 pages, 8 figures, 2 tables

The MOSDEF survey: a census of AGN-driven ionized outflows at $z = 1.4-3.8$

Authors: Gene C. K. Leung, Alison L. Coil, James Aird, Mojegan Azadi, Mariska Kriek, Bahram Mobasher, Naveen Reddy, Alice Shapley, Brian Siana, Tara Fetherolf, Francesca M. Fornasini, William R. Freeman, Sedona H. Price, Ryan L. Sanders, Irene Shivaei, Tom Zick

Abstract: Using data from the MOSFIRE Deep Evolution Field (MOSDEF) survey, we present a census of AGN-driven ionized outflows in a sample of 159 AGNs at $1.4 \le z \le 3.8$. The sample spans AGN bolometric luminosities of $10^{44-47} \mathrm{~erg~s}^{-1}$ and includes both quiescent and star-forming galaxies extending across three orders of magnitude in stellar mass. We identify and characterize outflows f… ▽ More Using data from the MOSFIRE Deep Evolution Field (MOSDEF) survey, we present a census of AGN-driven ionized outflows in a sample of 159 AGNs at $1.4 \le z \le 3.8$. The sample spans AGN bolometric luminosities of $10^{44-47} \mathrm{~erg~s}^{-1}$ and includes both quiescent and star-forming galaxies extending across three orders of magnitude in stellar mass. We identify and characterize outflows from the \hbeta, [OIII], \halpha ~and [NII] emission line spectra. We detect outflows in $17\%$ of the AGNs, seven times more often than in a mass-matched sample of inactive galaxies in MOSDEF. The outflows are fast and galaxy-wide, with velocities of $\sim 400-3500 ~\mathrm{km~s}^{-1}$ and spatial extents of $0.3-11.0$ kpc. The incidence of outflows among AGNs is independent of the stellar mass of the host galaxy, with outflows detected in both star-forming and quiescent galaxies. This suggests that outflows exist across different phases in galaxy evolution. We investigate relations between outflow kinematic, spatial, and energetic properties and both AGN and host galaxy properties. Our results show that AGN-driven outflows are widespread in galaxies along the star-forming main sequence. The mass-loading factors of the outflows are typically $0.1-1$ and increase with AGN luminosity, capable of exceeding unity at $L_\mathrm{AGN} \gtrsim 10^{46.3} \mathrm{~erg~s}^{-1}$. In these more luminous sources the ionized outflow alone is likely sufficient to regulate star formation, and when combined with outflowing neutral and molecular gas may be able to quench star formation in their host galaxies. △ Less

Submitted 25 September, 2019; v1 submitted 30 May, 2019; originally announced May 2019.

Comments: 24 pages, 12 figures, accepted for publication in the The Astrophysical Journal

Cosmic evolution of supermassive black holes: A view into the next two decades

Authors: Francesca Civano, Nico Cappelluti, Ryan Hickox, Rebecca Canning, James Aird, Marco Ajello, Steve Allen, Eduardo Bañados, Laura Blecha, William N. Brandt, Marcella Brusa, Francisco Carrera, Massimo Cappi, Andrea Comastri, Klaus Dolag, Megan Donahue, Martin Elvis, Giuseppina Fabbiano, Francesca Fornasini, Poshak Gandhi, Antonis Georgakakis, Kelly Holley-Bockelmann, Anton Koekemoer, Andrew Goulding, Mackenzie Jones , et al. (16 additional authors not shown)

Abstract: The discoveries made over the past 20 years by Chandra and XMM-Newton surveys in conjunction with multiwavelength imaging and spectroscopic data available in the same fields have significantly changed the view of the supermassive black hole (SMBH) and galaxy connection. These discoveries have opened up several exciting questions that are beyond the capabilities of current X-ray telescopes and will… ▽ More The discoveries made over the past 20 years by Chandra and XMM-Newton surveys in conjunction with multiwavelength imaging and spectroscopic data available in the same fields have significantly changed the view of the supermassive black hole (SMBH) and galaxy connection. These discoveries have opened up several exciting questions that are beyond the capabilities of current X-ray telescopes and will need to be addressed by observatories in the next two decades. As new observatories peer into the early Universe, we will begin to understand the physics and demographics of SMBH infancy (at $z>6$) and investigate the influence of their accretion on the formation of the first galaxies ($§$ 2.1). We will also be able to understand the accretion and evolution over the cosmic history (at $z\sim$1-6) of the full population of black holes in galaxies, including low accretion rate, heavily obscured AGNs at luminosities beyond the reach of current X-ray surveys ($§$2.2 and $§$2.3), enabling us to resolve the connection between SMBH growth and their environment. △ Less

Submitted 26 March, 2019; originally announced March 2019.

Comments: 5 pages, 3 figures; Science white paper submitted to the Astro2020 Decadal Survey

Deep ugrizY Imaging and DEEP2/3 Spectroscopy: A Photometric Redshift Testbed for LSST and Public Release of Data from the DEEP3 Galaxy Redshift Survey

Authors: Rongpu Zhou, Michael C. Cooper, Jeffrey A. Newman, Matthew L. N. Ashby, James Aird, Christopher J. Conselice, Marc Davis, Aaron A. Dutton, S. M. Faber, Jerome J. Fang, G. G. Fazio, Puragra Guhathakurta, Dale Kocevski, David C. Koo, Kirpal Nandra, Andrew C. Phillips, David J. Rosario, Edward F. Schlafly, Jonathan R. Trump, Benjamin Weiner, Christopher N. A. Willmer, Renbin Yan

Abstract: We present catalogs of calibrated photometry and spectroscopic redshifts in the Extended Groth Strip, intended for studies of photometric redshifts (photo-z's). The data includes ugriz photometry from CFHTLS and Y-band photometry from the Subaru Suprime camera, as well as spectroscopic redshifts from the DEEP2, DEEP3 and 3D-HST surveys. These catalogs incorporate corrections to produce effectively… ▽ More We present catalogs of calibrated photometry and spectroscopic redshifts in the Extended Groth Strip, intended for studies of photometric redshifts (photo-z's). The data includes ugriz photometry from CFHTLS and Y-band photometry from the Subaru Suprime camera, as well as spectroscopic redshifts from the DEEP2, DEEP3 and 3D-HST surveys. These catalogs incorporate corrections to produce effectively matched-aperture photometry across all bands, based upon object size information available in the catalog and Moffat profile point spread function fits. We test this catalog with a simple machine learning-based photometric redshift algorithm based upon Random Forest regression, and find that the corrected aperture photometry leads to significant improvement in photo-z accuracy compared to the original SExtractor catalogs from CFHTLS and Subaru. The deep ugrizY photometry and spectroscopic redshifts are well-suited for empirical tests of photometric redshift algorithms for LSST. The resulting catalogs are publicly available. We include a basic summary of the strategy of the DEEP3 Galaxy Redshift Survey to accompany the recent public release of DEEP3 data. △ Less

Submitted 12 August, 2019; v1 submitted 19 March, 2019; originally announced March 2019.

Comments: 20 pages, 9 figures; published in MNRAS; for associated catalogs, see http://d-scholarship.pitt.edu/36064/

Journal ref: MNRAS 488, 4565-4584 (2019)

The Advanced X-ray Imaging Satellite

Authors: Richard F. Mushotzky, James Aird, Amy J. Barger, Nico Cappelluti, George Chartas, Lia Corrales, Rafael Eufrasio, Andrew C. Fabian, Abraham D. Falcone, Elena Gallo, Roberto Gilli, Catherine E. Grant, Martin Hardcastle, Edmund Hodges-Kluck, Erin Kara, Michael Koss, Hui Li, Carey M. Lisse, Michael Loewenstein, Maxim Markevitch, Eileen T. Meyer, Eric D. Miller, John Mulchaey, Robert Petre, Andrew J. Ptak , et al. (19 additional authors not shown)

Abstract: Much of the baryonic matter in the Universe, including the most active and luminous sources, are best studied in the X-ray band. Key advances in X-ray optics and detectors have paved the way for the Advanced X-ray Imaging Satellite (AXIS), a Probe-class mission that is a major improvement over Chandra, which has generated a steady stream of important discoveries for the past 2 decades. AXIS can be… ▽ More Much of the baryonic matter in the Universe, including the most active and luminous sources, are best studied in the X-ray band. Key advances in X-ray optics and detectors have paved the way for the Advanced X-ray Imaging Satellite (AXIS), a Probe-class mission that is a major improvement over Chandra, which has generated a steady stream of important discoveries for the past 2 decades. AXIS can be launched in the late 2020s and will transform our understanding in several major areas of astrophysics, including the growth and fueling of supermassive black holes, galaxy formation and evolution, the microphysics of cosmic plasmas, the time-variable universe, and a wide variety of cutting-edge studies. Relative to Chandra, the AXIS PSF is nearly twice as sharp on-axis; its field of view for subarcsecond imaging 70 times larger by area; its effective area at 1 keV is 10 times larger. The low-Earth orbit ensures a low and stable detector background, resulting in 50 times greater sensitivity than Chandra for extended sources. AXIS has a rapid repointing response with operations similar to Swift, but is 100 times more sensitive for time-domain science. These capabilities open up a vast discovery space and complement the next generation of astronomical observatories. A high-spectral-resolution mission (Athena) operating at the same time as a high-angular-resolution mission (AXIS) greatly increases the range of scientific discovery. AXIS will use lightweight X-ray optics made of thin single-crystal silicon mirrors developed at NASA Goddard. The detector array builds on a long legacy of X-ray CCD and provides improved photon localization, much faster readout time, and broader energy band. The estimated mission costs are consistent with the $1B Probe mission cost guideline. △ Less

Submitted 14 March, 2019; v1 submitted 10 March, 2019; originally announced March 2019.

Comments: A Probe-class mission study commissioned by NASA for the NAS Astro2020 Decadal Survey. Cost section redacted. 66 pages, 41 figures. v2: minor fixes

4MOST: Project overview and information for the First Call for Proposals

Authors: R. S. de Jong, O. Agertz, A. Agudo Berbel, J. Aird, D. A. Alexander, A. Amarsi, F. Anders, R. Andrae, B. Ansarinejad, W. Ansorge, P. Antilogus, H. Anwand-Heerwart, A. Arentsen, A. Arnadottir, M. Asplund, M. Auger, N. Azais, D. Baade, G. Baker, S. Baker, E. Balbinot, I. K. Baldry, M. Banerji, S. Barden, P. Barklem , et al. (313 additional authors not shown)

Abstract: We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolut… ▽ More We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = λ/Δλ\sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations. △ Less

Submitted 1 April, 2019; v1 submitted 6 March, 2019; originally announced March 2019.

Comments: Part of the 4MOST issue of The Messenger, published in preparation of 4MOST Community Workshop, see http://www.eso.org/sci/meetings/2019/4MOST.html

Journal ref: The Messenger, 2019, 175, 3

Exploring the halo occupation of AGN using dark-matter cosmological simulations

Authors: A. Georgakakis, J. Comparat, A. Merloni, L. Ciesla, J. Aird, A. Finoguenov

Abstract: A semi-empirical model is presented that describes the distribution of Active Galactic Nuclei (AGN) on the cosmic web. It populates dark-matter halos in N-body simulations (MultiDark) with galaxy stellar masses using empirical relations based on abundance matching techniques, and then paints accretion events on these galaxies using state-of-the-art measurements of the AGN occupation of galaxies. T… ▽ More A semi-empirical model is presented that describes the distribution of Active Galactic Nuclei (AGN) on the cosmic web. It populates dark-matter halos in N-body simulations (MultiDark) with galaxy stellar masses using empirical relations based on abundance matching techniques, and then paints accretion events on these galaxies using state-of-the-art measurements of the AGN occupation of galaxies. The explicit assumption is that the large-scale distribution of AGN is independent of the physics of black-hole fueling. The model is shown to be consistent with current measurements of the two-point correlation function of AGN samples. It is then used to make inferences on the halo occupation of the AGN population. Mock AGN are found in halos with a broad distribution of masses with a mode of $\approx 10^{12}\,h^{-1} \, M_{\odot}$ and a tail extending to cluster-size halos. The clustering properties of the model AGN depend only weakly on accretion luminosity and redshift. The fraction of satellite AGN in the model increases steeply toward more massive halos, in contrast with some recent observational results. This discrepancy, if confirmed, could point to a dependence of the halo occupation of AGN on the physics of black-hole fueling. △ Less

Submitted 10 December, 2018; originally announced December 2018.

Comments: 22 pages, submitted to MNRAS, light-cones and correlation functions available at https://doi.org/10.5281/zenodo.2114283

Chandra Observations of NuSTAR Serendipitous Sources near the Galactic Plane

Authors: John A. Tomsick, George B. Lansbury, Farid Rahoui, James Aird, David M. Alexander, Maica Clavel, AnaSofija Cuturilo, Francesca M. Fornasini, JaeSub Hong, Lizelke Klindt, Daniel Stern

Abstract: The NuSTAR serendipitous survey has already uncovered a large number of Active Galactic Nuclei (AGN), providing new information about the composition of the Cosmic X-ray Background. For the AGN off the Galactic plane, it has been possible to use the existing X-ray archival data to improve source localizations, identify optical counterparts, and classify the AGN with optical spectroscopy. However,… ▽ More The NuSTAR serendipitous survey has already uncovered a large number of Active Galactic Nuclei (AGN), providing new information about the composition of the Cosmic X-ray Background. For the AGN off the Galactic plane, it has been possible to use the existing X-ray archival data to improve source localizations, identify optical counterparts, and classify the AGN with optical spectroscopy. However, near the Galactic Plane, better X-ray positions are necessary to achieve optical or near-IR identifications due to the higher levels of source crowding. Thus, we have used observations with the Chandra X-ray Observatory to obtain the best possible X-ray positions. With eight observations, we have obtained coverage for 19 NuSTAR serendips within 12 deg of the plane. One or two Chandra sources are detected within the error circle of 15 of the serendips, and we report on these sources and search for optical counterparts. For one source (NuSTAR J202421+3350.9), we obtained a new optical spectrum and detected the presence of hydrogen emission lines. The source is Galactic, and we argue that it is likely a Cataclysmic Variable. For the other sources, the Chandra positions will enable future classifications in order to place limits on faint Galactic populations, including high-mass X-ray binaries and magnetars. △ Less

Submitted 13 November, 2018; originally announced November 2018.

Comments: 16 pages, accepted by ApJ

X-rays across the galaxy population - III. The incidence of AGN as a function of star formation rate

Authors: James Aird, Alison L. Coil, Antonis Georgakakis

Abstract: We map the co-eval growth of galaxies and their central supermassive black holes in detail by measuring the incidence of Active Galactic Nuclei (AGN) in galaxies as a function of star formation rate (SFR) and redshift (to z~4). We combine large galaxy samples with deep Chandra X-ray imaging to measure the probability distribution of specific black hole accretion rates (LX relative to stellar mass)… ▽ More We map the co-eval growth of galaxies and their central supermassive black holes in detail by measuring the incidence of Active Galactic Nuclei (AGN) in galaxies as a function of star formation rate (SFR) and redshift (to z~4). We combine large galaxy samples with deep Chandra X-ray imaging to measure the probability distribution of specific black hole accretion rates (LX relative to stellar mass) and derive robust AGN fractions and average specific accretion rates. First, we consider galaxies along the main sequence of star formation. We find a linear correlation between the average SFR and both the AGN fraction and average specific accretion rate across a wide range in stellar mass ($M_* \sim 10^{8.5-11.5}M_\odot$) and to at least z~2.5, indicating that AGN in main-sequence galaxies are driven by the stochastic accretion of cold gas. We also consider quiescent galaxies and find significantly higher AGN fractions than predicted, given their low SFRs, indicating that AGN in quiescent galaxies are fuelled by additional mechanisms (e.g. stellar winds). Next, we bin galaxies according to their SFRs relative to the main sequence. We find that the AGN fraction is significantly elevated for galaxies that are still star-forming but with SFRs below the main sequence, indicating further triggering mechanisms enhance AGN activity within these sub-main-sequence galaxies. We also find that the incidence of high-accretion-rate AGN is enhanced in starburst galaxies and evolves more mildly with redshift than within the rest of the galaxy population, suggesting mergers play a role in driving AGN activity in such high-SFR galaxies. △ Less

Submitted 9 January, 2019; v1 submitted 10 October, 2018; originally announced October 2018.

Comments: 19 pages, 15 figures. This is a pre-copyedited, author-produced version of an article accepted for publication in MNRAS following peer review

X-UDS: The Chandra Legacy Survey of the UKIDSS Ultra Deep Survey Field

Authors: Dale D. Kocevski, Guenther Hasinger, Murray Brightman, Kirpal Nandra, Antonis Georgakakis, Nico Cappelluti, Francesca Civano, Yuxuan Li, Yanxia Li, James Aird, David M. Alexander, Omar Almaini, Marcella Brusa, Johannes Buchner, Andrea Comastri, Christopher J. Conselice, Mark A. Dickinson, Alexis Finoguenov, Roberto Gilli, Anton M. Koekemoer, Takamitsu Miyaji, James R. Mullaney, Casey Papovich, David Rosario, Mara Salvato , et al. (3 additional authors not shown)

Abstract: We present the X-UDS survey, a set of wide and deep Chandra observations of the Subaru-XMM Deep/UKIDSS Ultra Deep Survey (SXDS/UDS) field. The survey consists of 25 observations that cover a total area of 0.33 deg$^{2}$. The observations are combined to provide a nominal depth of ~600 ksec in the central 100 arcmin$^{2}$ region of the field that has been imaged with Hubble/WFC3 by the CANDELS surv… ▽ More We present the X-UDS survey, a set of wide and deep Chandra observations of the Subaru-XMM Deep/UKIDSS Ultra Deep Survey (SXDS/UDS) field. The survey consists of 25 observations that cover a total area of 0.33 deg$^{2}$. The observations are combined to provide a nominal depth of ~600 ksec in the central 100 arcmin$^{2}$ region of the field that has been imaged with Hubble/WFC3 by the CANDELS survey and $\sim$200 ksec in the remainder of the field. In this paper, we outline the survey's scientific goals, describe our observing strategy, and detail our data reduction and point source detection algorithms. Our analysis has resulted in a total of 868 band-merged point sources detected with a false-positive Poisson probability of $<1\times10^{-4}$. In addition, we present the results of an X-ray spectral analysis and provide best-fitting neutral hydrogen column densities, $N_{\rm H}$, as well as a sample of 51 Compton-thick active galactic nucleus candidates. Using this sample, we find the intrinsic Compton-thick fraction to be 30-35% over a wide range in redshift ($z=0.1-3$), suggesting the obscured fraction does not evolve very strongly with epoch. However, if we assume that the Compton-thick fraction is dependent on luminosity, as is seen for Compton-thin sources, then our results are consistent with a rise in the obscured fraction out to $z\sim3$. Finally, an examination of the host morphologies of our Compton-thick candidates shows a high fraction of morphological disturbances, in agreement with our previous results. All data products described in this paper are made available via a public website. △ Less

Submitted 21 September, 2018; originally announced September 2018.

Comments: 15 pages, 11 figures, published in the Astrophysics Journal Supplement Series

Journal ref: Bibliographic Code: 2018ApJS..236...48K

THE MOSDEF SURVEY: The Nature of Mid-Infrared Excess Galaxies and a Comparison of IR and UV Star Formation Tracers at z~2

Authors: Mojegan Azadi, Alison Coil, James Aird, Irene Shivaei, Naveen Reddy, Alice Shapley, Mariska Kriek, William R. Freeman, Gene C. K. Leung, Bahram Mobasher, Sedona H. Price, Ryan L. Sanders, Brian Siana, Tom Zick

Abstract: We present an analysis using the MOSFIRE Deep Evolution Field (MOSDEF) survey on the nature of "MIR-excess" galaxies, which have star formation rates (SFR) inferred from mid-infrared (MIR) data that is substantially elevated relative to that estimated from dust-corrected UV data. We use a sample of $\sim$200 galaxies and AGN at $1.40<z<2.61$ with 24 $μ$m detections (rest-frame 8$μ$m) from MIPS/\te… ▽ More We present an analysis using the MOSFIRE Deep Evolution Field (MOSDEF) survey on the nature of "MIR-excess" galaxies, which have star formation rates (SFR) inferred from mid-infrared (MIR) data that is substantially elevated relative to that estimated from dust-corrected UV data. We use a sample of $\sim$200 galaxies and AGN at $1.40<z<2.61$ with 24 $μ$m detections (rest-frame 8$μ$m) from MIPS/\textit{Spitzer}. We find that the identification of MIR-excess galaxies strongly depends on the methodologies used to estimate IR luminosity ($\rm L_{IR}$) and to correct the UV light for dust attenuation. We find that extrapolations of the SFR from the observed 24 $μ$m flux, using luminosity-dependent templates based on local galaxies, substantially overestimate $\rm L_{IR}$ in $z\sim2$ galaxies. By including \textit{Herschel} observations and using a stellar mass-dependent, luminosity-independent $\rm L_{IR}$, we obtain more reliable estimates of the SFR and a lower fraction of MIR-excess galaxies. Once stellar mass selection biases are taken into account, we identify $\sim24\%$ of our galaxies as MIR-excess. However, $\rm SFR_{Hα}$ is not elevated in MIR-excess galaxies compared to MIR-normal galaxies, indicating that the intrinsic fraction of MIR-excess may be lower. Using X-ray, IR, and optically-selected AGN in MOSDEF, we do not find a higher prevalence for AGN in MIR-excess galaxies relative to MIR-normal galaxies. A stacking analysis of X-ray undetected galaxies does not reveal a harder spectrum in MIR-excess galaxies relative to MIR-normal galaxies. Our analysis indicates that AGN activity does not contribute substantially to the MIR excess and instead implies that it is likely due to the enhanced PAH emission. △ Less

Submitted 23 June, 2018; originally announced June 2018.

Evidence for a mass-dependent AGN Eddington ratio distribution via the flat relationship between SFR and AGN luminosity

Authors: E. Bernhard, J. R. Mullaney, J. Aird, R. C. Hickox, M. L. Jones, F. Stanley, L. P. Grimmett, E. Daddi

Abstract: The lack of a strong correlation between AGN X-ray luminosity ($L_X$; a proxy for AGN power) and the star formation rate (SFR) of their host galaxies has recently been attributed to stochastic AGN variability. Studies using population synthesis models have incorporated this by assuming a broad, universal (i.e. does not depend on the host galaxy properties) probability distribution for AGN specific… ▽ More The lack of a strong correlation between AGN X-ray luminosity ($L_X$; a proxy for AGN power) and the star formation rate (SFR) of their host galaxies has recently been attributed to stochastic AGN variability. Studies using population synthesis models have incorporated this by assuming a broad, universal (i.e. does not depend on the host galaxy properties) probability distribution for AGN specific X-ray luminosities (i.e. the ratio of $L_X$ to host stellar mass; a common proxy for Eddington ratio). However, recent studies have demonstrated that this universal Eddington ratio distribution fails to reproduce the observed X-ray luminosity functions beyond z$\sim$1.2. Furthermore, empirical studies have recently shown that the Eddington ratio distribution may instead depend upon host galaxy properties, such as SFR and/or stellar mass. To investigate this further we develop a population synthesis model in which the Eddington ratio distribution is different for star-forming and quiescent host galaxies. We show that, although this model is able to reproduce the observed X-ray luminosity functions out to z$\sim$2, it fails to simultaneously reproduce the observed flat relationship between SFR and X-ray luminosity. We can solve this, however, by incorporating a mass dependency in the AGN Eddington ratio distribution for star-forming host galaxies. Overall, our models indicate that a relative suppression of low Eddington ratios ($λ_{\rm Edd}\lesssim$0.1) in lower mass galaxies (M<$10^{10-11}$Msun) is required to reproduce both the observed X-ray luminosity functions and the observed flat SFR/X-ray relationship. △ Less

Submitted 19 February, 2018; v1 submitted 23 January, 2018; originally announced January 2018.

Comments: 17 pages, 12 figures, accepted for publication in MNRAS (fixed: issues with some figures not up-to-date)

The NuSTAR Extragalactic Surveys: X-ray spectroscopic analysis of the bright hard-band selected sample

Authors: L. Zappacosta, A. Comastri, F. Civano, S. Puccetti, F. Fiore, J. Aird, A. Del Moro, G. B. Lansbury, G. Lanzuisi, A. Goulding, J. R. Mullaney, D. Stern, M. Ajello, D. M. Alexander, D. R. Ballantyne, F. E. Bauer, W. N. Brandt, C. -T. J. Chen, D. Farrah, F. A. Harrison, P. Gandhi, L. Lanz, A. Masini, S. Marchesi, C. Ricci , et al. (1 additional authors not shown)

Abstract: We discuss the spectral analysis of a sample of 63 Active Galactic Nuclei (AGN) detected above a limiting flux of $S(8-24$ keV$)=7\times10^{-14}$ erg/s/cm$^2$ in the multi-tiered NuSTAR Extragalactic Survey program. The sources span a redshift range z=0-2.1 (median $\langle$z$\rangle=$0.58). The spectral analysis is performed over the broad 0.5-24 keV energy range, combining NuSTAR with Chandra an… ▽ More We discuss the spectral analysis of a sample of 63 Active Galactic Nuclei (AGN) detected above a limiting flux of $S(8-24$ keV$)=7\times10^{-14}$ erg/s/cm$^2$ in the multi-tiered NuSTAR Extragalactic Survey program. The sources span a redshift range z=0-2.1 (median $\langle$z$\rangle=$0.58). The spectral analysis is performed over the broad 0.5-24 keV energy range, combining NuSTAR with Chandra and/or XMM-Newton data and employing empirical and physically motivated models. This constitutes the largest sample of AGN selected at $ >10$ keV to be homogeneously spectrally analyzed at these flux levels. We study the distribution of spectral parameters such as photon index, column density ($N_{\rm H}$), reflection parameter ($R$) and 10-40 keV luminosity ($L_{X}$). Heavily obscured ($\log [N_{\rm H}/cm^{-2}]\ge23$) and Compton Thick (CT; $\log [N_{\rm H}/cm^{-2}]\ge24$) AGN constitute $\sim$25% (15-17 sources) and $\sim$2-3% ( 1-2 sources) of the sample, respectively. The observed $N_{\rm H}$ distribution fairly agrees with predictions of Cosmic X-ray Background population synthesis models (CXBPSM). We estimate the intrinsic fraction of AGN as a function of $N_{H}$, accounting for the bias against obscured AGN in a flux-selected sample. The fraction of CT AGN relative to $\log [N_{\rm H}/cm^{-2}]=20-24$ AGN is poorly constrainted, formally in the range 2-56% (90% upper limit of 66%). We derived a fraction ($f_{abs}$) of obscured AGN ($\log [N_{\rm H}/cm^{-2}]=22-24$) as a function of $L_{X}$ in agreement with CXBPSM and previous z < 1 X-ray determinations. Furthermore $f_{abs}$ at z=0.1-0.5 and $\log (L_{X}/erg$ $s^{-1})\approx43.6-44.3$ agrees with observational measurements/trends obtained over larger redshift intervals. We report a significant anti-correlation of $R$ with $L_{X}$ (confirmed by our companion paper on stacked spectra) with considerable scatter around the median $R$ values. △ Less

Submitted 12 January, 2018; originally announced January 2018.

Comments: Comments: 38 pages, 17 figures, 8 tables. Accepted for publication in ApJ

The NuSTAR Extragalactic Survey: Average broad-band X-ray spectral properties of the NuSTAR detected AGN

Authors: A. Del Moro, D. M. Alexander, J. A. Aird, F. E. Bauer, F. Civano, J. R. Mullaney, D. R. Ballantyne, W. N. Brandt, A. Comastri, P. Gandhi, F. A. Harrison, G. B. Lansbury, L. Lanz, B. Luo, S. Marchesi, S. Puccetti, C. Ricci, C. Saez, D. Stern, E. Treister, L. Zappacosta

Abstract: We present a study of the average X-ray spectral properties of the sources detected by the NuSTAR extragalactic survey, comprising observations of the E-CDFS, EGS and COSMOS fields. The sample includes 182 NuSTAR sources (64 detected at 8-24 keV), with 3-24 keV fluxes ranging between $f_{\rm 3-24 keV}\approx10^{-14}$ and $6\times10^{-13}$ erg/cm$^2$/s (… ▽ More We present a study of the average X-ray spectral properties of the sources detected by the NuSTAR extragalactic survey, comprising observations of the E-CDFS, EGS and COSMOS fields. The sample includes 182 NuSTAR sources (64 detected at 8-24 keV), with 3-24 keV fluxes ranging between $f_{\rm 3-24 keV}\approx10^{-14}$ and $6\times10^{-13}$ erg/cm$^2$/s ($f_{\rm 8-24 keV}\approx3\times10^{-14}-3\times10^{-13}$ erg/cm$^2$/s) and redshifts of $z=0.04-3.21$. We produce composite spectra from the Chandra+NuSTAR data ($E\approx2-40$ keV, rest frame) for all the sources with redshift identifications (95%) and investigate the intrinsic, average spectra of the sources, divided into broad-line (BL) and narrow-line (NL) AGN, and also in different bins of X-ray column density and luminosity. The average power-law photon index for the whole sample is $Γ=1.65_{-0.03}^{+0.03}$, flatter than $Γ\approx1.8$ typically found for AGN. While the spectral slope of BL and X-ray unabsorbed AGN is consistent with typical values ($Γ=1.79_{-0.01}^{+0.01}$), a significant flattening is seen in NL AGN and heavily-absorbed sources ($Γ=1.60_{-0.05}^{+0.08}$ and $Γ=1.38_{-0.12}^{+0.12}$, respectively), likely due to the effect of absorption and to the contribution from Compton reflection to the high-energy flux (E>10 keV). We find that the typical reflection fraction in our spectra is $R\approx0.5$ (for $Γ=1.8$), with a tentative indication of an increase of the reflection strength with column density. While there is no significant evidence for a dependence of the photon index with X-ray luminosity in our sample, we find that $R$ decreases with luminosity, with relatively high levels of reflection ($R\approx1.2$) for $L_{\rm 10-40 keV}<10^{44}$ erg/s and $R\approx0.3$ for $L_{\rm 10-40 keV}>10^{44}$ erg/s AGN, assuming $Γ=1.8$. △ Less

Submitted 3 October, 2017; originally announced October 2017.

Comments: 19 pages, 14 figures, 4 tables. Accepted for publication in ApJ

The NuSTAR Serendipitous Survey: Hunting for The Most Extreme Obscured AGN at >10 keV

Authors: G. B. Lansbury, D. M. Alexander, J. Aird, P. Gandhi, D. Stern, M. Koss, I. Lamperti, M. Ajello, A. Annuar, R. J. Assef, D. R. Ballantyne, M. Balokovic, F. E. Bauer, N. Brandt, M. Brightman, C. -T. J. Chen, F. Civano, A. Comastri, A. D. Moro, C. Fuentes, F. A. Harrison, S. Marchesi, A. Masini, J. R. Mullaney, C. Ricci , et al. (5 additional authors not shown)

Abstract: We identify sources with extremely hard X-ray spectra (i.e., with photon indices of Gamma<0.6 in the 13 sq. deg. NuSTAR serendipitous survey, to search for the most highly obscured AGNs detected at >10 keV. Eight extreme NuSTAR sources are identified, and we use the NuSTAR data in combination with lower energy X-ray observations (from Chandra, Swift XRT, and XMM-Newton) to characterize the broad-b… ▽ More We identify sources with extremely hard X-ray spectra (i.e., with photon indices of Gamma<0.6 in the 13 sq. deg. NuSTAR serendipitous survey, to search for the most highly obscured AGNs detected at >10 keV. Eight extreme NuSTAR sources are identified, and we use the NuSTAR data in combination with lower energy X-ray observations (from Chandra, Swift XRT, and XMM-Newton) to characterize the broad-band (0.5-24 keV) X-ray spectra. We find that all of the extreme sources are highly obscured AGNs, including three robust Compton-thick (CT; N_H > 1.5e24 cm^-2) AGNs at low redshift (z<0.1), and a likely-CT AGN at higher redshift (z=0.16). Most of the extreme sources would not have been identified as highly obscured based on the low energy (<10 keV) X-ray coverage alone. The multiwavelength properties (e.g., optical spectra and X-ray/MIR luminosity ratios) provide further support for the eight sources being significantly obscured. Correcting for absorption, the intrinsic rest-frame 10-40 keV luminosities of the extreme sources cover a broad range, from ~ 5 x 10^42 to 10^45 erg s^-1. The estimated number counts of CT AGNs in the NuSTAR serendipitous survey are in broad agreement with model expectations based on previous X-ray surveys, except for the lowest redshifts (z<0.07) where we measure a high CT fraction of f_CT^obs = 30 (+16 -12) %. For the small sample of CT AGNs, we find a high fraction of galaxy major mergers (50 +/- 33%) compared to control samples of "normal" AGNs. △ Less

Submitted 16 August, 2017; v1 submitted 20 July, 2017; originally announced July 2017.

Comments: 19 pages, 15 figures, 5 tables; Accepted for publication in ApJ; Author list updated

The mean star formation rates of unobscured QSOs: searching for evidence of suppressed or enhanced star formation

Authors: F. Stanley, D. M. Alexander, C. M. Harrison, D. J. Rosario, L. Wang, J. A. Aird, N. Bourne, L. Dunne, S. Dye, S. Eales, K. K. Knudsen, M. J. Michalowski, E. Valiante, G. De Zotti, C. Furlanetto, R. Ivison, S. Maddox, M. W. L. Smith

Abstract: We investigate the mean star formation rates (SFRs) in the host galaxies of ~3000 optically selected QSOs from the SDSS survey within the Herschel-ATLAS fields, and a radio-luminous sub-sample, covering the redshift range of z = 0.2-2.5. Using WISE & Herschel photometry (12 - 500μm) we construct composite SEDs in bins of redshift and AGN luminosity. We perform SED fitting to measure the mean infra… ▽ More We investigate the mean star formation rates (SFRs) in the host galaxies of ~3000 optically selected QSOs from the SDSS survey within the Herschel-ATLAS fields, and a radio-luminous sub-sample, covering the redshift range of z = 0.2-2.5. Using WISE & Herschel photometry (12 - 500μm) we construct composite SEDs in bins of redshift and AGN luminosity. We perform SED fitting to measure the mean infrared luminosity due to star formation, removing the contamination from AGN emission. We find that the mean SFRs show a weak positive trend with increasing AGN luminosity. However, we demonstrate that the observed trend could be due to an increase in black hole (BH) mass (and a consequent increase of inferred stellar mass) with increasing AGN luminosity. We compare to a sample of X-ray selected AGN and find that the two populations have consistent mean SFRs when matched in AGN luminosity and redshift. On the basis of the available virial BH masses, and the evolving BH mass to stellar mass relationship, we find that the mean SFRs of our QSO sample are consistent with those of main sequence star-forming galaxies. Similarly, the radio-luminous QSOs have mean SFRs that are consistent with both the overall QSO sample and with star-forming galaxies on the main sequence. In conclusion, on average QSOs reside on the main sequence of star-forming galaxies, and the observed positive trend between the mean SFRs and AGN luminosity can be attributed to BH mass and redshift dependencies. △ Less

Submitted 17 July, 2017; originally announced July 2017.

Comments: 22 pages (incl. appendix), 17 Figures, re-submitted to MNRAS with referee comments addressed

The MOSDEF Survey: First Measurement of Nebular Oxygen Abundance at $z>4$

Authors: Alice E. Shapley, Ryan L. Sanders, Naveen A. Reddy, Mariska Kriek, William R. Freeman, Bahram Mobasher, Brian Siana, Alison L. Coil, Gene C. K. Leung, Laura deGroot, Irene Shivaei, Sedona H. Price, Mojegan Azadi, James Aird

Abstract: We present the first spectroscopic measurement of multiple rest-frame optical emission lines at $z>4$. During the MOSFIRE Deep Evolution Field (MOSDEF) survey, we observed the galaxy GOODSN-17940 with the Keck I/MOSFIRE spectrograph. The K-band spectrum of GOODSN-17940 includes significant detections of the [OII]$λλ3726,3729$, [NeIII]$\lambda3869$, and H$γ$ emission lines and a tentative detection… ▽ More We present the first spectroscopic measurement of multiple rest-frame optical emission lines at $z>4$. During the MOSFIRE Deep Evolution Field (MOSDEF) survey, we observed the galaxy GOODSN-17940 with the Keck I/MOSFIRE spectrograph. The K-band spectrum of GOODSN-17940 includes significant detections of the [OII]$λλ3726,3729$, [NeIII]$\lambda3869$, and H$γ$ emission lines and a tentative detection of H$δ$, indicating $z_{\rm{spec}}=4.4121$. GOODSN-17940 is an actively star-forming $z>4$ galaxy based on its K-band spectrum and broadband spectral energy distribution. A significant excess relative to the surrounding continuum is present in the Spitzer/IRAC channel 1 photometry of GOODSN-17940, due primarily to strong H$α$ emission with a rest-frame equivalent width of $\mbox{EW(H}α)=1200$ Å. Based on the assumption of $0.5 Z_{\odot}$ models and the Calzetti attenuation curve, GOODSN-17940 is characterized by $M_*=5.0^{+4.3}_{-0.2}\times 10^9 M_{\odot}$. The Balmer decrement inferred from H$α$/H$γ$ is used to dust correct the H$α$ emission, yielding $\mbox{SFR(H}α)=320^{+190}_{-140} M_{\odot}\mbox{ yr}^{-1}$. These $M_*$ and SFR values place GOODSN-17940 an order of magnitude in SFR above the $z\sim 4$ star-forming "main sequence." Finally, we use the observed ratio of [NeIII]/[OII] to estimate the nebular oxygen abundance in GOODSN-17940, finding $\mbox{O/H}\sim 0.2 \mbox{ (O/H)}_{\odot}$. Combining our new [NeIII]/[OII] measurement with those from stacked spectra at $z\sim 0, 2, \mbox{ and } 3$, we show that GOODSN-17940 represents an extension to $z>4$ of the evolution towards higher [NeIII]/[OII] (i.e., lower $\mbox{O/H}$) at fixed stellar mass. It will be possible to perform the measurements presented here out to $z\sim 10$ using the James Webb Space Telescope. △ Less

Submitted 23 August, 2017; v1 submitted 17 July, 2017; originally announced July 2017.

Comments: 7 pages, 3 figures, accepted to ApJ Letters

Galactic Sources Detected in the NuSTAR Serendipitous Survey

Authors: John A. Tomsick, George B. Lansbury, Farid Rahoui, Maica Clavel, Francesca M. Fornasini, JaeSub Hong, James Aird, David M. Alexander, Arash Bodaghee, Jeng-Lun Chiu, Jonathan E. Grindlay, Charles J. Hailey, Fiona A. Harrison, Roman A. Krivonos, Kaya Mori, Daniel Stern

Abstract: The Nuclear Spectroscopic Telescope Array (NuSTAR) provides an improvement in sensitivity at energies above 10 keV by two orders of magnitude over non-focusing satellites, making it possible to probe deeper into the Galaxy and Universe. Lansbury and collaborators recently completed a catalog of 497 sources serendipitously detected in the 3-24 keV band using 13 deg2 of NuSTAR coverage. Here, we rep… ▽ More The Nuclear Spectroscopic Telescope Array (NuSTAR) provides an improvement in sensitivity at energies above 10 keV by two orders of magnitude over non-focusing satellites, making it possible to probe deeper into the Galaxy and Universe. Lansbury and collaborators recently completed a catalog of 497 sources serendipitously detected in the 3-24 keV band using 13 deg2 of NuSTAR coverage. Here, we report on an optical and X-ray study of 16 Galactic sources in the catalog. We identify eight of them as stars (but some or all could have binary companions), and use information from Gaia to report distances and X-ray luminosities for three of them. There are four CVs or CV candidates, and we argue that NuSTAR J233426-2343.9 is a relatively strong CV candidate based partly on an X-ray spectrum from XMM-Newton. NuSTAR J092418-3142.2, which is the brightest serendipitous source in the Lansbury catalog, and NuSTAR J073959-3147.8 are LMXB candidates, but it is also possible that these two sources are CVs. One of the sources is a known HMXB, and NuSTAR J105008-5958.8 is a new HMXB candidate, which has strong Balmer emission lines in its optical spectrum and a hard X-ray spectrum. We discuss the implications of finding these HMXBs for the surface density (logN-logS) and luminosity function of Galactic HMXBs. We conclude that, with the large fraction of unclassified sources in the Galactic plane detected by NuSTAR in the 8-24 keV band, there could be a significant population of low luminosity HMXBs. △ Less

Submitted 23 May, 2017; originally announced May 2017.

Comments: 16 pages, accepted by ApJ

ESO-Athena Synergy White Paper

Authors: P. Padovani, F. Combes, M. Diaz Trigo, S. Ettori, E. Hatziminaoglou, P. Jonker, M. Salvato, S. Viti, C. Adami, J. Aird, D. Alexander, P. Casella, C. Ceccarelli, E. Churazov, M. Cirasuolo, E. Daddi, A. Edge, C. Feruglio, V. Mainieri, S. Markoff, A. Merloni, F. Nicastro, P. O'Brien, L. Oskinova, F. Panessa , et al. (7 additional authors not shown)

Abstract: The Advanced Telescope for High ENergy Astrophysics (Athena) is the X-ray observatory mission selected by ESA within its Cosmic Vision 2015-2025 programme to address the Hot and Energetic Universe scientific theme. The ESO-Athena Synergy Team (EAST) has been tasked to single out the potential scientific synergies between Athena and optical/near-infrared (NIR) and sub/mm ground based facilities, in… ▽ More The Advanced Telescope for High ENergy Astrophysics (Athena) is the X-ray observatory mission selected by ESA within its Cosmic Vision 2015-2025 programme to address the Hot and Energetic Universe scientific theme. The ESO-Athena Synergy Team (EAST) has been tasked to single out the potential scientific synergies between Athena and optical/near-infrared (NIR) and sub/mm ground based facilities, in particular those of ESO (i.e., the VLT and ELT, ALMA and APEX), by producing a White Paper to identify and develop the: 1. needs to access ESO ground-based facilities to achieve the formulated Athena science objectives; 2. needs to access Athena to achieve the formulated science objectives of ESO facilities contemporary to Athena; 3. science areas where the synergetic use of Athena and ESO facilities in the late 2020s will result in scientific added value. Community input to the process happened primarily via a dedicated ESO - Athena Synergy Workshop that took place on Sept. 14 - 16, 2016 at ESO, Garching. This White Paper presents the results of the EAST's work, sorted by synergy area, and deals with the following topics: 1. the Hot Universe: Early groups and clusters and their evolution, Physics of the Intracluster medium, Missing baryons in cosmic filaments; 2. the Energetic Universe: Supermassive black hole (SMBH) history, SMBH accretion disks, Active Galactic Nuclei feedback - Molecular outflows, Ultra-fast outflows, Accretion Physics, Transient Science; 3. Observatory Science: Star Formation, Stars. It then discusses the optical-NIR-sub-mm perspective by providing details on VLT/MOONS, the E-ELT instruments, in particular the MOS, VISTA/4MOST, the ESO and ALMA archives, future ALMA and ESO developments, and finally the (likely) ESO - Athena astronomical scene in the 2020s. (abridged) △ Less

Submitted 17 May, 2017; originally announced May 2017.

Comments: 70 pages, 16 figures

Observational constraints on the specific accretion-rate distribution of X-ray selected AGN

Authors: A. Georgakakis, J. Aird, A. Schulze, T. Dwelly, M. Salvato, K. Nandra, A. Merloni, Donald P. Schneider

Abstract: This paper estimates the specific accretion-rate distribution of AGN using a sample of 4821 X-ray sources from both deep and shallow surveys. The specific accretion-rate distribution is defined as the probability of a galaxy with a given stellar mass and redshift hosting an active nucleus with a certain specific accretion rate. We find that the probability of a galaxy hosting an AGN increases with… ▽ More This paper estimates the specific accretion-rate distribution of AGN using a sample of 4821 X-ray sources from both deep and shallow surveys. The specific accretion-rate distribution is defined as the probability of a galaxy with a given stellar mass and redshift hosting an active nucleus with a certain specific accretion rate. We find that the probability of a galaxy hosting an AGN increases with decreasing specific accretion rate. There is evidence that this trend reverses at low specific accretion rates, $λ<10^{-4}-10^{-3}$ (in Eddington units). There is also a break close to the Eddington limit, above which the probability of an accretion event decreases steeply. The specific accretion-rate distribution evolves such that the fraction of AGN among galaxies drops toward lower redshifts. This decrease in the AGN duty cycle is responsible for the strong evolution of the accretion density of the Universe from redshift $z\approx1-1.5$ to the present day. Our analysis also suggests that this evolution is accompanied by a decoupling of accretion events onto black holes from the formation of stars in galaxies. There is also evidence that at earlier times the relative probability of high vs low specific accretion-rate events among galaxies increases. We argue that this differential redshift evolution of the AGN duty cycle with respect to $λ$ produces the AGN downsizing trend, whereby luminous sources peak at earlier epochs compared to less luminous ones. Finally, we also find a stellar-mass dependence of the specific accretion-rate distribution, with more massive galaxies avoiding high specific accretion-rate events. △ Less

Submitted 2 May, 2017; originally announced May 2017.

Comments: MNRAS submitted

X-rays across the galaxy population - II. The distribution of AGN accretion rates as a function of stellar mass and redshift

Authors: James Aird, Alison L. Coil, Antonis Georgakakis

Abstract: We use deep Chandra X-ray imaging to measure the distribution of specific black hole accretion rates ($L_X$ relative to the stellar mass of the galaxy) and thus trace AGN activity within star-forming and quiescent galaxies, as a function of stellar mass (from $10^{8.5}-10^{11.5} M_\odot$) and redshift (to $z \sim 4$). We adopt near-infrared selected samples of galaxies from the CANDELS and UltraVI… ▽ More We use deep Chandra X-ray imaging to measure the distribution of specific black hole accretion rates ($L_X$ relative to the stellar mass of the galaxy) and thus trace AGN activity within star-forming and quiescent galaxies, as a function of stellar mass (from $10^{8.5}-10^{11.5} M_\odot$) and redshift (to $z \sim 4$). We adopt near-infrared selected samples of galaxies from the CANDELS and UltraVISTA surveys, extract X-ray data for every galaxy, and use a flexible Bayesian method to combine these data and to measure the probability distribution function of specific black hole accretion rates, $λ_{sBHAR}$. We identify a broad distribution of $λ_{sBHAR}$ in both star-forming and quiescent galaxies---likely reflecting the stochastic nature of AGN fuelling---with a roughly power-law shape that rises toward lower $λ_{sBHAR}$, a steep cutoff at $λ_{sBHAR} \gtrsim 0.1-1$ (in Eddington equivalent units), and a turnover or flattening at $λ_{sBHAR} \lesssim 10^{-3}-10^{-2}$. We find that the probability of a star-forming galaxy hosting a moderate $λ_{sBHAR}$ AGN depends on stellar mass and evolves with redshift, shifting toward higher $λ_{sBHAR}$ at higher redshifts. This evolution is truncated at a point corresponding to the Eddington limit, indicating black holes may self-regulate their growth at high redshifts when copious gas is available. The probability of a quiescent galaxy hosting an AGN is generally lower than that of a star-forming galaxy, shows signs of suppression at the highest stellar masses, and evolves strongly with redshift. The AGN duty cycle in high-redshift ($z\gtrsim2$) quiescent galaxies thus reaches $\sim$20 per cent, comparable to the duty cycle in star-forming galaxies of equivalent stellar mass and redshift. △ Less

Submitted 13 October, 2017; v1 submitted 2 May, 2017; originally announced May 2017.

Comments: 20 pages + 3 appendices (27 pages total) with 15 figures. The sketch shown in Figure 9 summarizes our findings. Data tables available at http://doi.org/10.5281/zenodo.1009605 . Updated FAST SED fitting code available at https://github.com/jamesaird/FAST . This is a pre-copyedited, author-produced version of an article accepted for publication in MNRAS following peer review