Search

A multitude of JWST studies reveal a surprising over-abundance of over-massive accreting super-massive blackholes (SMBHs) -- leading to a deepening tension between theory and observation in the first billion years of cosmic time. Across X-ray to infrared wavelengths, models built off of pre-JWST predictions fail to easily reproduce observed AGN signatures (or lack thereof), driving uncertainty around the true nature of these sources. Using a sample of JWST AGN identified via their broadened Halpha emission and covered by the deepest X-ray surveys, we find neither any measurable X-ray emission nor any detection of high-ionization emission lines frequently associated with accreting SMBHs. We propose that these sources are accreting at or beyond the Eddington limit, which reduces the need for efficient production of heavy SMBH seeds at cosmic dawn. Using a theoretical model of super-Eddington accretion, we can produce the observed relative dearth of both X-ray and ultraviolet emission, as well as the high Balmer decrements, without the need for significant dust attenuation. This work indicates that super-Eddington accretion is easily achieved through-out the early Universe, and further study is required to determine what environments are required to trigger this mode of black hole growth., Comment: Submitted

We present a multi-wavelength analysis of MCG+11-11-032, a nearby AGN with the unique classification of both a binary and a dual AGN candidate. With new Chandra observations we aim to resolve any dual AGN system via imaging data, and search for signs of a binary AGN via analysis of the X-ray spectrum. Analyzing the Chandra spectrum, we find no evidence of previously suggested double-peaked Fe K$\alpha$ lines; the spectrum is instead best fit by an absorbed powerlaw with a single Fe K$\alpha$ line, as well as an additional line centered at $\approx$7.5 keV. The Chandra observation reveals faint, soft, and extended X-ray emission, possibly linked to low-level nuclear outflows. Further analysis shows evidence for a compact, hard source -- MCG+11-11-032 X2 -- located 3.27'' from the primary AGN. Modeling MCG+11-11-032 X2 as a compact source, we find that it is relatively luminous ($L_{\text{2$-$10 keV}} = 1.52_{-0.48}^{+0.96}\times 10^{41}$ erg s$^{-1}$), and the location is coincident with an compact and off-nuclear source resolved in Hubble Space Telescope infrared (F105W) and ultraviolet (F621M, F547M) bands. Pairing our X-ray results with a 144 MHz radio detection at the host galaxy location, we observe X-ray and radio properties similar to those of ESO 243-49 HLX-1, suggesting that MCG+11-11-032 X2 may be a hyper-luminous X-ray source. This detection with Chandra highlights the importance of a high-resolution X-ray imager, and how previous binary AGN candidates detected with large-aperture instruments benefit from high-resolution follow-up. Future spatially resolved optical spectra, and deeper X-ray observations, can better constrain the origin of MCG+11-11-032 X2., Comment: 16 pages, 7 figures, submitted to ApJ

We present a comprehensive X-ray analysis and spectral energy distribution (SED) fitting of WISEA J171419.96+602724.6, an extremely luminous type 2 quasar at $z$ = 2.99. The source was suggested as a candidate Compton-thick (column density N$_{\rm H}>$1.5 $\times$ 10$^{24}$ cm$^{-2}$) quasar by a short XMM-Newton observation in 2011. We recently observed the source with deep NuSTAR and XMM-Newton exposures in 2021 and found that the source has a lower obscuration of N$_{\rm H}\sim$5 $\times$ 10$^{22}$ cm$^{-2}$ with an about four times lower flux. The two epochs of observations suggested that the source was significantly variable in X-ray obscuration, flux, and intrinsic luminosity at 2-3~$\sigma$ in less than 2.5 years (in the source rest frame). We performed SED fitting of this source using CIGALE thanks to its great availability of multiwavelength data (from hard X-rays to radio). The source is very luminous with a bolometric luminosity of $L_{\rm BOL}\sim$ 2.5 $\times$ 10$^{47}$ erg s$^{-1}$. Its host galaxy has a huge star formation rate (SFR) of $\sim$1280 Solar mass yr$^{-1}$ and a huge stellar mass of $\sim$1.1 $\times$ 10$^{12}$ Solar mass. The correlation between the SFR and stellar mass of this source is consistent with what was measured in the high-$z$ quasars. It is also consistent with what was measured in the main-sequence star-forming galaxies, suggesting that the presence of the active nucleus in our target does not enhance or suppress the SFR of its host galaxy. The source is an Infrared hyper-luminous, obscured galaxy with significant amount of hot dust in its torus and shares many similar properties with hot, dust obscured galaxies., Comment: 19 pages, 10 figure, 5 tables, accepted to ApJ

The Chandra Source Catalog (CSC) is a virtual X-ray astrophysics facility that enables both detailed individual source studies and statistical studies of large samples of X-ray sources detected in ACIS and HRC-I imaging observations obtained by the Chandra X-ray Observatory. The catalog provides carefully-curated, high-quality, and uniformly calibrated and analyzed tabulated positional, spatial, photometric, spectral, and temporal source properties, as well as science-ready X-ray data products. The latter includes multiple types of source- and field-based FITS format products that can be used as a basis for further research, significantly simplifying followup analysis of scientifically meaningful source samples. We discuss in detail the algorithms used for the CSC Release 2 Series, including CSC 2.0, which includes 317,167 unique X-ray sources on the sky identified in observations released publicly through the end of 2014, and CSC 2.1, which adds Chandra data released through the end of 2021 and expands the catalog to 407,806 sources. Besides adding more recent observations, the CSC Release 2 Series includes multiple algorithmic enhancements that provide significant improvements over earlier releases. The compact source sensitivity limit for most observations is ~5 photons over most of the field of view, which is ~2x fainter than Release 1, achieved by co-adding observations and using an optimized source detection approach. A Bayesian X-ray aperture photometry code produces robust fluxes even in crowded fields and for low count sources. The current release, CSC 2.1, is tied to the Gaia-CRF3 astrometric reference frame for the best sky positions for catalog sources., Comment: 66 pages, 17 figures, 16 tables, accepted for publication in The Astrophysical Journal Supplement Series

Recent James Webb Space Telescope (JWST) observations have revealed a surprisingly abundant population of faint, dusty active galactic nuclei (AGNs) at z~4-7. Together with the presence of supermassive black holes (SMBHs) at z>6, this raises questions about the formation and growth histories of early black holes. Current theories for the formation of seed black holes from the death of the first stars (i.e. light seeds) and/or the direct collapse of primordial gas clouds (i.e. heavy seeds) still lack observational confirmation. Here, we present LID-568, a low-mass (7.2e6Msun) black hole hosting powerful outflows that is observed in an extreme phase of rapid growth at z~4. This object is similar to other JWST-discovered faint AGN populations, but is bright in X-ray emission and accreting at more than 4000% of the limit at which radiation pressure exceeds the force of gravitational attraction of the black hole (i.e. super-Eddington accretion). Analysis of JWST NIRSpec/IFU data reveals spatially extended Ha emission with velocities of ~ -600 - -500 km/s relative to the central black hole, indicative of robust nuclear-driven outflows. LID-568 represents an elusive low-mass black hole experiencing super-Eddington accretion as invoked by models of early black hole formation. This discovery showcases a previously undiscovered key parameter space and offers crucial insights into rapid black hole growth mechanisms in the early universe., Comment: Submitted

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., Comment: Accepted for publication in ApJ:S. 57 pages, 32 figures

We present the third catalog release of the wide-area (31.3 deg$^2$) Stripe 82 X-ray survey. This catalog combines previously published X-ray source properties with multiwavelength counterparts and photometric redshifts, presents 343 new spectroscopic redshifts, and provides black hole masses for 1297 Type 1 Active Galactic Nuclei (AGN). With spectroscopic redshifts for 3457 out of 6181 Stripe 82X sources, the survey has a spectroscopic completeness of 56%. This completeness rises to 90% when considering the contiguous portions of the Stripe 82X survey with homogeneous X-ray coverage at an optical magnitude limit of $r<22$. Within that portion of the survey, 23% of AGN can be considered obscured by being either a Type 2 AGN, reddened ($R-K > 4$, Vega), or X-ray obscured with a column density $N_{\rm H} > 10^{22}$ cm$^{-2}$. Unlike other surveys, there is only a 18% overlap between Type 2 and X-ray obscured AGN. We calculated black hole masses for Type 1 AGN that have SDSS spectra using virial mass estimators calibrated on the H$\beta$,MgII, H$\alpha$, and CIV emission lines. We find wide scatter in these black hole mass estimates, indicating that statiscal analyses should use black hole masses calculated from the same formula to minimize bias. We find that the AGN with the highest X-ray luminosities are accreting at the highest Eddington ratios, consistent with the picture that most black hole mass accretion happens in the phase when the AGN is luminous ($L_{\rm 2-10 keV} > 10^{45}$ erg s$^{-1}$)., Comment: 30 pages, 8 figures, accepted for publication to ApJ

We present the second NuSTAR and XMM-Newton extragalactic survey of the JWST North Ecliptic Pole (NEP) Time-Domain Field (TDF). The first NuSTAR NEP-TDF survey (Zhao et al. 2021) had 681 ks total exposure time executed in NuSTAR cycle 5, in 2019 and 2020. This second survey, acquired from 2020 to 2022 in cycle 6, adds 880 ks of NuSTAR exposure time. The overall NuSTAR NEP-TDF survey is the most sensitive NuSTAR extragalactic survey to date, and a total of 60 sources were detected above the 95% reliability threshold. We constrain the hard X-ray number counts, logN-log S, down to 1.7 x 10$^{-14}$ erg cm$^{-2}$ s$^{-1}$ at 8-24 keV and detect an excess of hard X-ray sources at the faint end. About 47% of the NuSTAR-detected sources are heavily obscured (NH > 10$^{23}$ cm$^{-2}$), and 18+20% of the NuSTAR-detected sources are Compton-thick (N>10$^{24}$ cm$^{-2}$). These fractions are consistent with those measured in other NuSTAR surveys. Four sources presented >2$\sigma$ variability in the 3-year survey. In addition to NuSTAR, a total of 62 ks of XMM-Newton observations were taken during NuSTAR cycle 6. The XMM-Newton observations provide soft X-ray (0.5-10keV) coverage in the same field and enable more robust identification of the visible and infrared counterparts of the NuSTAR-detected sources. A total of 286 soft X-ray sources were detected, out of which 214 XMM-Newton sources have secure counterparts from multiwavelength catalogs., Comment: 37 pages, 27 figures, Accepted by ApJ

A fundamental goal of modern-day astrophysics is to understand the connection between supermassive black hole (SMBH) growth and galaxy evolution. Merging galaxies offer one of the most dramatic channels for galaxy evolution known, capable of driving inflows of gas into galactic nuclei, potentially fueling both star formation and central SMBH activity. Dual active galactic nuclei (dual AGNs) in late-stage mergers with nuclear pair separations $<10$ kpc are thus ideal candidates to study SMBH growth along the merger sequence since they coincide with the most transformative period for galaxies. However, dual AGNs can be extremely difficult to confirm and study. Hard X-ray ($>10$ keV) studies offer a relatively contamination-free tool for probing the dense obscuring environments predicted to surround the majority of dual AGN in late-stage mergers. To date, only a handful of the brightest and closest systems have been studied at these energies due to the demanding instrumental requirements involved. We demonstrate the unique capabilities of HEX-P to spatially resolve the soft and - for the first time - hard X-ray counterparts of closely-separated ($\sim2''-5''$) dual AGNs in the local Universe. By incorporating state-of-the-art physical torus models, we reproduce realistic broadband X-ray spectra expected for deeply embedded accreting SMBHs. Hard X-ray spatially resolved observations of dual AGNs - accessible only to HEX-P - will hence transform our understanding of dual AGN in the nearby Universe., Comment: Under review for Frontiers in Astronomy and Space Sciences

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., Comment: 28 pages, 13 Figures, 2 tables; submitted to Frontiers in Astronomy and Space Sciences; minor edits in the text, abstract and figure labels

Spectral energy distributions (SEDs) from X-ray to far-infrared (FIR) wavelengths are presented for a sample of 1246 X-ray luminous active galactic nuclei (AGN; $L_{0.5-10\rm{keV}}>10^{43}$ erg s$^{-1}$), with $z_{\rm{spec}}<1.2$, selected from Stripe 82X, COSMOS, and GOODS-N/S. The rest-frame SEDs show a wide spread ($\sim2.5$ dex) in the relative strengths of broad continuum features at X-ray, ultraviolet (UV), mid-infrared (MIR), and FIR wavelengths. A linear correlation (log-log slope of 0.7$\pm0.04$) is found between $L_{\rm{MIR}}$ and $L_{\rm{X}}$. There is significant scatter in the relation between the $L_{\rm{UV}}$ and $L_{\rm{X}}$ due to heavy obscuration, however the most luminous and unobscured AGN show a linear correlation (log-log slope of 0.8$\pm0.06$) in the relation above this scatter. The relation between $L_{\rm{FIR}}$ and $L_{\rm{X}}$ is predominantly flat, but with decreasing dispersion at $L_{\rm{X}}>10^{44}$ erg s$^{-1}$. The ratio between the "galaxy subtracted" bolometric luminosity and the intrinsic $L_{\rm{X}}$ increases from a factor of $\sim$$10-70$ from log $L_{\rm{bol}}/{\rm(erg\; s}^{-1})=44.5-46.5$. Characteristic SED shapes have been determined by grouping AGN based on relative strengths of the UV and MIR emission. The average $L_{1\mu\rm{m}}$ is constant for the majority of these SED shapes, while AGN with the strongest UV and MIR emission have elevated $L_{1\mu\rm{m}}$, consistent with the AGN emission dominating their SEDs at optical and NIR wavelengths. A strong correlation is found between the SED shape and both the $L_{\rm{X}}$ and $L_{\rm{bol}}$, such that $L_{\rm{bol}}/L_{\rm{X}}=20.4\pm1.8$, independent of the SED shape. This is consistent with an evolutionary scenario of increasing $L_{\rm{bol}}$ with decreasing obscuration as the AGN blows away circumnuclear gas., Comment: 24 pages, 20 figures, 4 appendix tables, Accepted to ApJ

Direct-collapse black holes (DCBHs) of mass $\sim 10^4$-$10^5 {M}_\odot$ that form in HI-cooling halos in the early Universe are promising progenitors of the $\gtrsim 10^9 {M}_\odot$ supermassive black holes that fuel observed $z \gtrsim 7$ quasars. Efficient accretion of the surrounding gas onto such DCBH seeds may render them sufficiently bright for detection with the JWST up to $z\approx 20$. Additionally, the very steep and red spectral slope predicted across the $\approx 1$-5 $\mu$m wavelength range of the JWST/NIRSpec instrument during their initial growth phase should make them photometrically identifiable up to very high redshifts. In this work, we present a search for such DCBH candidates across the 34 arcmin$^{2}$ in the first two spokes of the JWST cycle-1 PEARLS survey of the north ecliptic pole time-domain field covering eight NIRCam filters down to a maximum depth of $\sim$ 29 AB mag. We identify two objects with spectral energy distributions consistent with the Pacucci et al. (2016) DCBH models. However, we also note that even with data in eight NIRCam filters, objects of this type remain degenerate with dusty galaxies and obscured active galactic nuclei over a wide range of redshifts. Follow-up spectroscopy would be required to pin down the nature of these objects. Based on our sample of DCBH candidates and assumptions on the typical duration of the DCBH steep-slope state, we set a conservative upper limit of $\lesssim 5\times 10^{-4}$ comoving Mpc$^{-3}$ (cMpc$^{-3}$) on the comoving density of host halos capable of hosting DCBHs with spectral energy distributions similar to the Pacucci et al. (2016) models at $z\approx 6$-14., Comment: 10 pages, 6 figures, 1 table, accepted for publication in A&A

A large fraction of the accreting supermassive black hole population is shrouded by copious amounts of gas and dust, particularly in the distant ($z\gtrsim1$) Universe. While much of the obscuration is attributed to a parsec-scale torus, there is a known contribution from the larger-scale host galaxy. Using JWST/NIRCam imaging from the COSMOS-Web survey, we probe the galaxy-wide dust distribution in X-ray selected AGN up to $z\sim2$. Here, we focus on a sample of three AGNs with their host galaxies exhibiting prominent dust lanes, potentially due to their edge-on alignment. These represent 27% (3 out of 11 with early NIRCam data) of the heavily obscured ($N_H>10^{23}$ cm$^{-2}$) AGN population. With limited signs of a central AGN in the optical and near-infrared, the NIRCam images are used to produce reddening maps $E(B-V)$ of the host galaxies. We compare the mean central value of $E(B-V)$ to the X-ray obscuring column density along the line-of-sight to the AGN ($N_H\sim10^{23-23.5}$ cm$^{-2}$). We find that the extinction due to the host galaxy is present ($0.6\lesssim E(B-V) \lesssim 0.9$; $1.9 \lesssim A_V \lesssim 2.8$) and significantly contributes to the X-ray obscuration at a level of $N_H\sim10^{22.5}$ cm$^{-2}$ assuming an SMC gas-to-dust ratio which amounts to $\lesssim$30% of the total obscuring column density. These early results, including three additional cases from CEERS, demonstrate the ability to resolve such dust structures with JWST and separate the different circumnuclear and galaxy-scale obscuring structures., Comment: 12 pages, 8 figures, Accepted for publication in ApJL

We present the analysis of the rest frame ultraviolet and optical spectra of 30 bright blue quasars at $z\sim3$, selected to examine the suitability of AGN as cosmological probes. In our previous works, we found an unexpectedly high fraction ($\approx 25 \%$) of X-ray weak quasars in the sample. The latter sources also display a flatter UV continuum and a broader and fainter CIV profile in the archival UV data with respect to their X-ray normal counterparts. Here we present new observations with the LBT in both the $zJ$ (rest-frame $\simeq$2300-3100 $\rm \mathring{A}$) and the $K_S$ ($\simeq$4750-5350 $\rm \mathring{A}$) bands. We estimated black hole masses ($M_{\rm BH}$) and Eddington ratios ($\lambda_{\rm Edd}$) from the from the H$\beta$ and MgII emission lines, finding that our $z\sim3$ quasars are on average highly accreting ($\langle \lambda_{\rm Edd} \rangle\simeq 1.2$ and $\langle M_{\rm BH} \rangle\simeq 10^{9.7}M_\odot$), with no difference in $\lambda_{\rm Edd}$ or $M_{\rm BH}$ between X-ray weak and X-ray normal quasars. From the $zJ$ spectra, we derive flux and equivalent width of MgII and FeII, finding that X-ray weak quasars display higher FeII/MgII ratios with respect to typical quasars. FeII/MgII ratios of X-ray normal quasars are instead consistent with other estimates up to $z\simeq6.5$, corroborating the idea of already chemically mature BLRs at early cosmic time. From the $K_S$ spectra, we find that all the X-ray weak quasars present generally weaker [OIII] emission (EW<10 $\rm \mathring{A}$) than the normal ones. The sample as a whole, however, abides by the known X-ray/[OIII] luminosity correlation, hence the different [OIII] properties are likely due to an intrinsically weaker [OIII] emission in X-ray weak objects, associated to the shape of the spectral energy distribution. We interpret these results in the framework of accretion-disc winds., Comment: 25 pages, 17 figures, Accepted for publication on Astronomy & Astrophysics

We report the detection of CO(6-5) and CO(7-6) and their underlying continua from the host galaxy of quasar J100758.264+211529.207 (P\=oniu\=a'ena) at z=7.5419, obtained with the NOrthern Extended Millimeter Array (NOEMA). P\=oniu\=a'ena belongs to the HYPerluminous quasars at the Epoch of ReionizatION (HYPERION) sample of 17 $z>6$ quasars selected to be powered by supermassive black holes (SMBH) which experienced the fastest mass growth in the first Gyr of the Universe. The one reported here is the highest-redshift measurement of the cold and dense molecular gas to date. The host galaxy is unresolved and the line luminosity implies a molecular reservoir of $\rm M(H_2)=(2.2\pm0.2)\times 10^{10}$ $\rm M_\odot$, assuming a CO spectral line energy distribution typical of high-redshift quasars and a conversion factor $\alpha=0.8$ $\rm M_{\odot} (K\,km \, s^{-1} \,pc^{2})^{-1} $. We model the cold dust spectral energy distribution (SED) to derive a dust mass of M$_{\rm dust} =(2.1\pm 0.7)\times 10^8$ $\rm M_\odot$, and thus a gas to dust ratio $\sim100$. Both the gas and dust mass are not dissimilar from the reservoir found for luminous quasars at $z\sim6$. We use the CO detection to derive an estimate of the cosmic mass density of $\rm H_2$, $\Omega_{H_2} \simeq 1.31 \times 10^{-5}$. This value is in line with the general trend suggested by literature estimates at $ z < 7 $ and agrees fairly well with the latest theoretical expectations of non-equilibrium molecular-chemistry cosmological simulations of cold gas at early times., Comment: Submitted to ApJ Letters

By characterizing the contribution of stray light to large datasets from the NuSTAR X-ray observatory collected over 2012--2017, we report a measurement of the cosmic X-ray background in the 3--20 keV energy range. These data represent $\sim20\%$ sky coverage while avoiding Galactic Ridge X-ray emission and are less weighted by deep, survey fields than previous measurements with NuSTAR. Images in narrow energy bands are stacked in detector space and spatially fit with a model representing the stray light and uniform pattern expected from the cosmic X-ray background and the instrumental background, respectively. We establish baseline flux values from Earth-occulted data and validate the fitting method on stray light observations of the Crab, which further serve to calibrate the resulting spectra. We present independent spectra of the cosmic X-ray background with the FPMA and FPMB detector arrays, which are in excellent agreement with the canonical characterization by HEAO 1 and are $10\%$ lower than most subsequent measurements; $F_{\rm{3-20~keV}}^{FPMA} = 2.63 \times 10^{-11}~\rm{erg~s^{-1}~cm^{-2}~deg^{-2}}$ and $F_{\rm{3-20~keV}}^{FPMB} = 2.58 \times 10^{-11}~\rm{erg~s^{-1}~cm^{-2}~deg^{-2}}$. We discuss these results in light of previous measurements of the cosmic X-ray background and consider the impact of systematic uncertainties on our spectra., Comment: 23 pages, 15 figures, 4 tables

We present ALMA Band 9 continuum observation of the ultraluminous quasi-stellar object (QSO) SDSS J0100+2802, providing a $\sim 10\sigma$ detection at $\sim 670$ GHz. SDSS J0100+2802 is the brightest QSO with the most massive super massive black hole (SMBH) known at $z>6$, and we study its dust spectral energy distribution in order to determine the dust properties and the star formation rate (SFR) of its host-galaxy. We obtain the most accurate estimate so far of the temperature, mass and emissivity index of the dust, having $T_{\rm dust}=48.4\pm2.3$ K, $M_{\rm dust}=(2.29\pm0.83)\times 10^7$ M$_\odot$, $\beta=2.63\pm 0.23$. This allows us to measure the SFR with the smallest statistical error for this QSO, SFR$=265\pm 32\ \rm M_\odot yr^{-1}$. Our results enable us to evaluate the relative growth of the SMBH and host galaxy of J0100+2802, finding that the SMBH is dominating the process of BH-galaxy growth in this QSO at $z=6.327$, when the Universe was $865$ Myr old. Such unprecedented constraints on the host galaxy SFR and dust temperature can only be obtained through high frequency observations, and highlight the importance of ALMA Band 9 to obtain a robust overview of the build-up of the first quasars' host galaxies at $z>6$., Comment: 10 pages, 4 figures, 1 table. Accepted for publication in ApJL

Extremely variable quasars can also show strong changes in broad-line emission strength and are known as changing-look quasars (CLQs). To study the CLQ transition mechanism, we present a pilot sample of CLQs with X-ray observations in both the bright and faint states. From a sample of quasars with bright-state archival SDSS spectra and (Chandra or XMM-Newton) X-ray data, we identified five new CLQs via optical spectroscopic follow-up, and then obtained new target-of-opportunity X-ray observations with Chandra. No strong absorption is detected in either the bright- or the faint-state X-ray spectra. The intrinsic X-ray flux generally changes along with the optical variability, and the X-ray power-law slope becomes harder in the faint state. Large amplitude mid-infrared variability is detected in all five CLQs, and the MIR variability echoes the variability in the optical with a time lag expected from the light-crossing time of the dusty torus for CLQs with robust lag measurements. The changing-obscuration model is not consistent with the observed X-ray spectra and spectral energy distribution changes seen in these CLQs. It is highly likely that the observed changes are due to the changing accretion rate of the supermassive black hole, so the multiwavelength emission varies accordingly, with promising analogies to the accretion states of X-ray binaries., Comment: Accepted for publication in ApJ; 30 pages, 12 figures

We give an overview and describe the rationale, methods, and first results from NIRCam images of the JWST "Prime Extragalactic Areas for Reionization and Lensing Science" ("PEARLS") project. PEARLS uses up to eight NIRCam filters to survey several prime extragalactic survey areas: two fields at the North Ecliptic Pole (NEP); seven gravitationally lensing clusters; two high redshift proto-clusters; and the iconic backlit VV 191 galaxy system to map its dust attenuation. PEARLS also includes NIRISS spectra for one of the NEP fields and NIRSpec spectra of two high-redshift quasars. The main goal of PEARLS is to study the epoch of galaxy assembly, AGN growth, and First Light. Five fields, the JWST NEP Time-Domain Field (TDF), IRAC Dark Field (IDF), and three lensing clusters, will be observed in up to four epochs over a year. The cadence and sensitivity of the imaging data are ideally suited to find faint variable objects such as weak AGN, high-redshift supernovae, and cluster caustic transits. Both NEP fields have sightlines through our Galaxy, providing significant numbers of very faint brown dwarfs whose proper motions can be studied. Observations from the first spoke in the NEP TDF are public. This paper presents our first PEARLS observations, their NIRCam data reduction and analysis, our first object catalogs, the 0.9-4.5 $\mu$m galaxy counts and Integrated Galaxy Light. We assess the JWST sky brightness in 13 NIRCam filters, yielding our first constraints to diffuse light at 0.9-4.5 {\mu}m. PEARLS is designed to be of lasting benefit to the community., Comment: Accepted to AJ, comments welcome. We ask anyone who uses our public PEARLS (NEP TDF) data to refer to this overview paper

We have computed obscured AGN redshifts using the XZ method, adopting a broad treatment in which we employed a wide-ranging data set and worked primarily at the XZ counts sensitivity threshold, culminating with a redshift catalog containing 121 sources that lack documented redshifts. We considered 363 obscured AGN from the Chandra Source Catalog Release 2.0, 59 of which were selected using multiwavelength criteria while 304 were X-ray selected. One-third of the data set had cross-matched spectroscopic or photometric redshifts. These sources, dominated by low-$z$ and low-$N_H$ AGN, were supplemented by 1000 simulations to form a data set for testing the XZ method. We used a multi-layer perceptron neural network to examine and predict cases in which XZ fails to reproduce the known redshift, yielding a classifier that can identify and discard poor redshift estimates. This classifier demonstrated a statistically significant $\sim$3$\sigma$ improvement over the existing XZ redshift information gain filter. We applied the machine learning model to sources with no documented redshifts, resulting in the 121-source new redshift catalog, all of which were X-ray selected. Our neural network's performance suggests that nearly 90% of these redshift estimates are consistent with hypothetical spectroscopic or photometric measurements, strengthening the notion that redshifts can be reliably estimated using only X-rays, which is valuable to current and future missions such as Athena. We have also identified a possible Compton-thick candidate that warrants further investigation., Comment: 34 pages, 20 figures, accepted by ApJ

We present a study of the relation between X-rays and ultraviolet emission in quasars for a sample of broad-line, radio-quiet objects obtained from the cross-match of the Sloan Digital Sky Survey DR14 with the latest Chandra Source Catalog 2.0 (2,332 quasars) and the Chandra COSMOS Legacy survey (273 quasars). The non-linear relation between the ultraviolet (at 2500 A, $L_{O}$) and the X-ray (at 2 keV, $L_{X}$) emission in quasars has been proved to be characterised by a smaller intrinsic dispersion than the observed one, as long as a homogeneous selection, aimed at preventing the inclusion of contaminants in the sample, is fulfilled. By leveraging on the low background of Chandra, we performed a complete spectral analysis of all the data available for the SDSS-CSC2.0 quasar sample (i.e. 3,430 X-ray observations), with the main goal of reducing the uncertainties on the source properties (e.g. flux, spectral slope). We analysed whether any evolution of the $L_{X}-L_{O}$ relation exists by dividing the sample in narrow redshift intervals across the redshift range spanned by our sample, $z \simeq 0.5-4$. We find that the slope of the relation does not evolve with redshift and it is consistent with the literature value of $0.6$ over the explored redshift range, implying that the mechanism underlying the coupling of the accretion disc and hot corona is the same at the different cosmic epochs. We also find that the dispersion decreases when examining the highest redshifts, where only pointed observations are available. These results further confirm that quasars are `standardisable candles', that is we can reliably measure cosmological distances at high redshifts where very few cosmological probes are available., Comment: 20 pages, 14 figures, accepted for publication in A&A

We analyse the properties of the CIV broad emission line in connection with the X-ray emission of 30 bright SDSS quasars at z~3.0-3.3 with pointed XMM-Newton observations, which were selected to test the suitability of AGN as cosmological tools. In our previous work, we found that a large fraction (~25%) of the quasars in this sample are X-ray underluminous by factors of >3-10. As absorbing columns of >10$^{23}$ cm$^{-2}$ can be safely ruled out, their weakness is most likely intrinsic. Here we explore possible correlations between the UV and X-ray features of these sources to investigate the origin of X-ray weakness. We fit their UV SDSS spectra and analyse their CIV properties (e.g., equivalent width, EW; line peak velocity, $\upsilon_{\rm peak}$) as a function of the X-ray photon index and 2-10 keV flux. We confirm the trends of CIV $\upsilon_{\rm peak}$ and EW with UV luminosity at 2500 angstrom for both X-ray weak and X-ray normal quasars, as well as the correlation between X-ray weakness and CIV EW. In contrast to some recent work, we do not observe any clear relation between the 2-10 keV luminosity and $\upsilon_{\rm peak}$. We find a correlation between the hard X-ray flux and the integrated CIV flux for X-ray normal quasars, whilst X-ray weak quasars deviate from the main trend by more than 0.5 dex. We argue that X-ray weakness might be interpreted in a starved X-ray corona picture associated with an ongoing disc-wind phase. If the wind is ejected in the vicinity of the black hole, the extreme-UV radiation that reaches the corona will be depleted, depriving the corona of seeds photons and generating an X-ray weak quasar. Yet, at the largest UV luminosities (>10$^{47}$ erg s$^{-1}$), there will still be an ample reservoir of ionising photons that can explain the excess CIV emission observed in the X-ray weak quasars with respect to normal sources of similar X-ray luminosities., Comment: 22 pages, 15 figures (with 3 more figures in the Appendix), abstract abridged. Accepted for publication in A&A

Many X-ray bright active galactic nuclei (AGN) are predicted to follow an extended stage of obscured black hole growth. In support of this picture we examine the X-ray undetected AGNs in the COSMOS field and compare their host galaxies with X-ray bright AGNs. We examine galaxies with M_\ast>10^{9.5}M_\odot for the presence of AGNs at redshifts $z=0.5-3$. We select AGNs in the infrared using \textit{Spitzer} and \textit{Herschel} detections and use color selection techniques to select AGNs within strongly star forming hosts. We stack \textit{Chandra} X-ray data of galaxies with an IR detection but lacking an X-ray detection to obtain soft and hard fluxes, allowing us to measure the energetics of these AGNs. We find a clear correlation between X-ray luminosity and IR AGN luminosity in the stacked galaxies. We also find that X-ray undetected AGNs all lie on the main sequence -- the tight correlation between SFR and $M_\ast$ that holds for the majority of galaxies, regardless of mass or redshift. This work demonstrates that there is a higher population of obscured AGNs than previously thought.

The discovery of hundreds of QSOs in the first Gyr of the Universe powered by already grown SMBHs challenges our knowledge of SMBH formation. In particular, investigations of $z>6$ QSOs presenting notable properties can provide unique information on the physics of fast SMBH growth in the early universe. We present the results of follow-up observations of the $z=6.515$ radio-quiet QSO PSO167-13, which is interacting with a close companion galaxy. The PSO167-13 system has been recently proposed to host the first heavily obscured X-ray source at high redshift. We observed PSO167-13 with Chandra/ACIS-S (177 ks), and obtained new spectroscopic observations (7.2 h) with Magellan/FIRE. No significant X-ray emission is detected from the PSO167-13 system, suggesting that the obscured X-ray source previously tentatively detected was either due to a strong background fluctuation or is highly variable. The upper limit (90% confidence level) on the X-ray emission of PSO167-13 ($L_{2-10\,\mathrm{keV}}<8.3\times10^{43}\,\mathrm{erg s^{-1}}$) is the lowest available for a $z>6$ QSO. The ratio between the X-ray and UV luminosity of $\alpha_{ox}<-1.95$ makes PSO167-13 a strong outlier from the $\alpha_{ox}-L_{UV}$ and $L_X-L_{\mathrm{bol}}$ relations. In particular, its X-ray emission is $>6$ times weaker than the expectation based on its UV luminosity. The new Magellan/FIRE spectrum of PSO167-13 is strongly affected by the unfavorable sky conditions, but the tentatively detected C IV and Mg II emission lines appear strongly blueshifted. The most plausible explanations for the X-ray weakness of PSO167-13 are intrinsic weakness or small-scale absorption by Compton-thick material. The possible strong blueshift of its emission lines hints at the presence of nuclear winds, which could be related to its X-ray weakness., Comment: Accepted for publication on A&A

The sensitivity of X-ray facilities and our ability to detect fainter active galactic nuclei (AGNs) will increase with the upcoming Athena mission and the AXIS and Lynx concept missions, thus improving our understanding of supermassive black holes (BHs) in a luminosity regime that can be dominated by X-ray binaries. We analyze the population of faint AGN (L_x (2-10 keV) < 10^42 erg/s) in the Illustris, TNG100, EAGLE, and SIMBA cosmological simulations, and find that the properties of their host galaxies vary from one simulation to another. In Illustris and EAGLE, faint AGN are powered by low-mass BHs located in low-mass star-forming galaxies. In TNG100 and SIMBA, they are mostly associated with more massive BHs in quenched massive galaxies. We model the X-ray binary populations (XRB) of the simulated galaxies, and find that AGN often dominate the galaxy AGN+XRB hard X-ray luminosity at z>2, while XRBs dominate in some simulations at z<2. Whether the AGN or XRB emission dominates in star-forming and quenched galaxies depends on the simulations. These differences in simulations can be used to discriminate between galaxy formation models with future high-resolution X-ray observations. We compare the luminosity of simulated faint AGN host galaxies to observations of stacked galaxies from Chandra. Our comparison indicates that the simulations post-processed with our X-ray modeling tend to overestimate the AGN+XRB X-ray luminosity; luminosity that can be strongly affected by AGN obscuration. Some simulations reveal clear AGN trends as a function of stellar mass (e.g., galaxy luminosity drop in massive galaxies), which are not apparent in the observations., Comment: Published version - MNRAS

Spin measurements of supermassive black holes (SMBHs) provide crucial constraints on the accretion processes that power active galactic nuclei (AGN), fuel outflows, and trigger black hole growth. However, spin measurements are mainly limited to a few dozen nearby sources for which high quality, high S/N spectra (e.g., from Chandra, XMM-Newton, Suzaku, NuSTAR) are available. Here we measure the average SMBH spin of $\sim$1900 AGN in the Chandra COSMOS-Legacy survey using spectral stacking analysis. We find broad Fe K$\alpha$ line emission in the average COSMOS spectrum (Gaussian width $\sigma=0.27\pm0.05$ keV), and by fitting this emission line profile with relativistic line models, we measure the average black hole spin parameter $a=0.62~\substack{+0.07 \\ -0.17}$. The sample size, availability of multiwavelength data, and spatial resolution of the COSMOS Legacy field also provide a unique environment to investigate the average SMBH spin as a function of other observables (e.g., redshift, luminosity) up to $z\sim5.3$. We find that optically classified Type 1 sources have broader Fe K$\alpha$ line emission than Type 2 sources. X-ray unobscured and obscured sources, as defined by their column densities, have widths that are consistent with the optically defined unobscured and obscured sources, respectively. There is some evidence for evolution of the Fe K$\alpha$ width and black hole spin parameter with luminosity, but not conclusively with redshift. The results of this work provide insights into the average spins of SMBHs in AGN, shedding light on their growth mechanisms and observed co-evolution with their host galaxies., Comment: 24 pages, 23 figures, submitted to ApJ

We present a comprehensive search for the 3.5 keV line, using $\sim$51 Ms of archival Chandra observations peering through the Milky Way's Dark Matter Halo from across the entirety of the sky, gathered via the Chandra Source Catalog Release 2.0. We consider the data's radial distribution, organizing observations into four data subsets based on angular distance from the Galactic Center. All data is modeled using both background-subtracted and background-modeled approaches to account for the particle instrument background, demonstrating statistical limitations of the currently-available $\sim$1 Ms of particle background data. A non-detection is reported in the total data set, allowing us to set an upper-limit on 3.5 keV line flux and constrain the sterile neutrino dark matter mixing angle. The upper-limit on sin$^2$(2$\theta$) is $2.58 \times 10^{-11}$ (though systematic uncertainty may increase this by a factor of $\sim$2), corresponding to the upper-limit on 3.5 keV line flux of $2.34 \times 10^{-7}$ ph s$^{-1}$ cm$^{-2}$. These limits show consistency with recent constraints and several prior detections. Non-detections are reported in all radial data subsets, allowing us to constrain the spatial profile of 3.5 keV line intensity, which does not conclusively differ from Navarro-Frenk-White predictions. Thus, while offering heavy constraints, we do not entirely rule out the sterile neutrino dark matter scenario or the more general decaying dark matter hypothesis for the 3.5 keV line. We have also used the non-detection of any unidentified emission lines across our continuum to further constrain the sterile neutrino parameter space., Comment: 37 pages, 27 figures, accepted by ApJ

The integrated X-ray luminosity ($L_{\mathrm{X}}$) of high-mass X-ray binaries (HMXBs) in a galaxy is correlated with its star formation rate (SFR), and the normalization of this correlation increases with redshift. Population synthesis models suggest that the redshift evolution of $L_{\mathrm{X}}$/SFR is driven by the metallicity ($Z$) dependence of HMXBs, and the first direct evidence of this connection was recently presented using galaxies at $z\sim2$. To confirm this result with more robust measurements and better constrain the $L_{\mathrm{X}}$-SFR-$Z$ relation, we have studied the $Z$ dependence of $L_{\mathrm{X}}$/SFR at lower redshifts. Using samples of star-forming galaxies at $z=0.1-0.9$ with optical spectra from the hCOSMOS and zCOSMOS surveys, we stacked \textit{Chandra} data from the COSMOS Legacy survey to measure the average $L_{\mathrm{X}}$/SFR as a function of $Z$ in three redshift ranges: $z=0.1-0.25$, $0.25-0.4$, and $0.5-0.9$. We find no significant variation of the $L_{\mathrm{X}}$-SFR-$Z$ relation with redshift. Our results provide further evidence that the $Z$ dependence of HMXBs is responsible for the redshift evolution of $L_{\mathrm{X}}$/SFR. Combining all available $z>0$ measurements together, we derive a best-fitting $L_{\mathrm{X}}$-SFR-$Z$ relation and assess how different population synthesis models describe the data. These results provide the strongest constraints to date on the $L_{\mathrm{X}}$-SFR-$Z$ relation in the range of $8.0<$12+log(O/H)$<9.0$., Comment: 13 pages, 7 figures, 2 tables

The masses of supermassive black holes at the centres of local galaxies appear to be tightly correlated with the mass and velocity dispersions of their galactic hosts. However, the local Mbh-Mstar relation inferred from dynamically measured inactive black holes is up to an order-of-magnitude higher than some estimates from active black holes, and recent work suggests that this discrepancy arises from selection bias on the sample of dynamical black hole mass measurements. In this work we combine X-ray measurements of the mean black hole accretion luminosity as a function of stellar mass and redshift with empirical models of galaxy stellar mass growth, integrating over time to predict the evolving Mbh-Mstar relation. The implied relation is nearly independent of redshift, indicating that stellar and black hole masses grow, on average, at similar rates. Matching the de-biased local Mbh-Mstar relation requires a mean radiative efficiency ~0.15, in line with theoretical expectations for accretion onto spinning black holes. However, matching the "raw" observed relation for inactive black holes requires a mean radiative efficiency around 0.02, far below theoretical expectations. This result provides independent evidence for selection bias in dynamically estimated black hole masses, a conclusion that is robust to uncertainties in bolometric corrections, obscured active black hole fractions, and kinetic accretion efficiency. For our fiducial assumptions, they favour moderate-to-rapid spins of typical supermassive black holes, to achieve a mean radiative efficiency ~0.12-0.20. Our approach has similarities to the classic Soltan analysis, but by using galaxy-based data instead of integrated quantities we are able to focus on regimes where observational uncertainties are minimized., Comment: 12 pages, 8 Figures. MNRAS, accepted. A discussion around the points raised by arXiv:1909.10821 is included in the Appendix

We investigate the cosmic evolution of the ratio between black hole mass (MBH) and host galaxy total stellar mass (Mstellar) out to z~2.5 for a sample of 100 X-ray-selected moderate-luminosity, broad-line active galactic nuclei (AGNs) in the Chandra-COSMOS Legacy Survey. By taking advantage of the deep multi-wavelength photometry and spectroscopy in the COSMOS field, we measure in a uniform way the galaxy total stellar mass using a SED decomposition technique and the black hole mass based on broad emission line measurements and single-epoch virial estimates. Our sample of AGN host galaxies has total stellar masses of 10^10-12Msun, and black hole masses of 10^7.0-9.5Msun. Combining our sample with the relatively bright AGN samples from the literature, we find no significant evolution of the MBH-Mstellar relation with black hole-to-host total stellar mass ratio of MBH/Mstellar~0.3% at all redshifts probed. We conclude that the average black hole-to-host stellar mass ratio appears to be consistent with the local value within the uncertainties, suggesting a lack of evolution of the MBH-Mstellar relation up to z~2.5., Comment: 16 pages, 7 figures, accepted for publication in ApJ

We present rest-frame far-infrared (FIR) and optical size measurements of AGN hosts and star-forming galaxies in the COSMOS field, enabled by high-resolution ALMA/1 mm (0.1 arcsec - 0.4 arcsec) and HST/F814W imaging (~ 0.1 arcsec). Our sample includes 27 galaxies at z<2.5, classified as infrared-selected AGN (3 sources), X-ray selected AGN (4 sources), and non-AGN star-forming galaxies (20 sources), for which high-resolution Band 6/7 ALMA images are available at 1 mm from our own observing program as well as archival observations. The sizes and SFR surface densities measured from both ALMA/1 mm and HST/F814W images show that obscured AGN host galaxies are more compact than non-AGN star-forming galaxies at similar redshift and stellar mass. This result suggests that the obscured accretion phase may be related to galaxies experiencing a compaction of their gaseous component, which could be associated with enhanced central star formation before a subsequent quenching driving the formation of compact passive galaxies. Moreover, most of the detected and stacked rest-frame FIR sizes of AGNs in our sample are similar or more compact than their rest-frame optical sizes, which is consistent with recent results of ALMA detected sources. This might be explained by the fact that the dusty starbursts take place in the compact regions, and suggests that the star formation mechanisms in the compact regions of AGN hosts are similar to those observed in star-forming galaxies observed with ALMA., Comment: accepted for publication in ApJ

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

X-ray emission from QSOs has been used to assess SMBH accretion properties up to $z$~6. However, at $z>6$ only ~15 QSOs are covered by sensitive X-ray observations, preventing a statistically significant investigation of the X-ray properties of QSOs in the first Gyr of the Universe. We present new Chandra observations of 10 $z>6$ QSOs, selected to have virial black-hole mass estimates from Mg II line spectroscopy. Adding archival X-ray data for an additional 15 $z>6$ QSOs, we investigate the X-ray properties of the QSO population in the first Gyr of the Universe, focusing in particular on the $L_{UV}-L_{X}$ relation, which is traced by the $\alpha_{ox}$ parameter, and the shape of their X-ray spectra. We performed photometric analyses to derive estimates of the X-ray luminosities, and thus the $\alpha_{ox}$ values and bolometric corrections ($K_{bol}=L_{bol}/L_{X}$). We compared the resulting $\alpha_{ox}$ and $K_{bol}$ distributions with the results found for QSO samples at lower redshift. Finally, we performed a basic X-ray spectral analysis of the brightest $z>6$ QSOs to derive their individual photon indices, and joint spectral analysis of the whole sample to estimate the average photon index. We confirm a lack of significant evolution of $\alpha_{ox}$ with redshift, extending the results from previous works up to $z>6$, and the trend of an increasing bolometric correction with increasing luminosity found for QSOs at lower redshifts. The average power-law photon index of our sample ($\Gamma=2.20_{-0.34}^{+0.39}$ and $\Gamma=2.13_{-0.13}^{+0.13}$ for sources with $<30$ and $>30$ net counts, respectively) is slightly steeper than, but still consistent with, typical QSOs at $z=1-6$. All these results point toward a lack of substantial evolution of the inner accretion-disk/hot-corona structure in QSOs from low redshift to $z>6$. Our data hint at generally high Eddington ratios at $z>6$., Comment: 15 pages. 10 figures. 7 tables. Accepted for publication in A&A

Active Galactic Nuclei (AGNs) hosting disk water megamasers are well known to be obscured by large amounts of gas, likely due to the presence along the line of sight of an almost edge-on disky structure orbiting the supermassive black hole. Correcting for the high obscuration is crucial to infer parameters intrinsic to the source, like its luminosity. We present a broadband X-ray spectral analysis of a water megamaser AGN in an early merger (NGC 5765B), combining Chandra and NuSTAR data. NGC 5765B is highly Compton-thick and reflection-dominated, following the general trend among disk megamasers. Combining the exquisite black hole mass from masers with our X-ray spectroscopy, the Eddington ratio of the megamaser is estimated to be in the $2-14\%$ range, and its robustness is confirmed through SED fitting., Comment: 10 pages, 8 figures. Accepted for publication in ApJ

While theoretical arguments predict that most of the early growth of supermassive black holes (SMBHs) happened during heavily obscured phases of accretion, current methods used for selecting $z>6$ quasars (QSOs) are strongly biased against obscured QSOs, thus considerably limiting our understanding of accreting SMBHs during the first Gyr of the Universe from an observational point of view. We report the $Chandra$ discovery of the first heavily obscured QSO candidate in the early universe, hosted by a close ($\approx5$ kpc) galaxy pair at $z=6.515$. One of the members is an optically classified type 1 QSO, PSO167-13. The companion galaxy was first detected as a [C II] emitter by ALMA. An X-ray source is significantly ($P=0.9996$) detected by $Chandra$ in the 2-5 keV band, with $<1.14$ net counts in the 0.5-2 keV band, although the current positional uncertainty does not allow a conclusive association with either PSO167-13 or its companion galaxy. From X-ray photometry and hardness-ratio arguments, we estimated an obscuring column density of $N_H>2\times10^{24}\,\mathrm{cm^{-2}}$ and $N_H>6\times10^{23}\,\mathrm{cm^{-2}}$ at $68\%$ and $90\%$ confidence levels, respectively. Thus, regardless of which of the two galaxies is associated with the X-ray emission, this source is the first heavily obscured QSO candidate at $z>6$., Comment: A&A Letters, accepted

We present the X-ray properties of 108 Dust-Obscured Galaxies (DOGs; F$_{24 \mu m}$/F$_{R} >$ 1000) in the COSMOS field, all of which detected in at least three far-infrared bands with the Herschel Observatory. Out of the entire sample, 22 are individually detected in the hard 2-8 keV X-ray band by the Chandra COSMOS Legacy survey, allowing us to classify them as AGN. Of them, 6 (27%) are Compton Thick AGN candidates with column densities N$_{H}$$>$10$^{24}$ cm$^{-2}$ while 15 are moderately obscured AGNs with 10$^{22}$ $<$ N$_{H}$ $<$ 10$^{24}$ cm$^{-2}$. Additionally, we estimate AGN contributions to the IR luminosity (8-1000$\mu$m rest-frame) greater than 20% for 19 DOGs based on SED decomposition using Spitzer/MIPS 24$\mu$m and the five Herschel bands (100-500 $\mu$m). Only 7 of these are detected in X-rays individually. We performed a X-ray stacking analysis for the 86 undetected DOGs. We find that the AGN fraction in DOGs increases with 24$\mu$m flux and that it is higher than that of the general 24$\mu$m population. However, no significant difference is found when considering only X-ray detections. This strongly motivates the combined use of X-ray and far-IR surveys to successfully probe a wider population of AGNs, particularly for the most obscured ones., Comment: To be published in ApJ

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 ($\S$ 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 ($\S$2.2 and $\S$2.3), enabling us to resolve the connection between SMBH growth and their environment., Comment: 5 pages, 3 figures; Science white paper submitted to the Astro2020 Decadal Survey

Gravitational wave (GW) and gravitational slingshot recoil kicks, which are natural products of SMBH evolution in merging galaxies, can produce active galactic "nuclei" that are offset from the centers of their host galaxies. Detections of offset AGN would provide key constraints on SMBH binary mass and spin evolution and on GW event rates. Although numerous offset AGN candidates have been identified, none have been definitively confirmed. Multi-wavelength observations with next-generation telescopes, including systematic large-area surveys, will provide unprecedented opportunities to identify and confirm candidate offset AGN from sub-parsec to kiloparsec scales. We highlight ways in which these observations will open a new avenue for multi-messenger studies in the dawn of low-frequency (~ nHz - mHz) GW astronomy., Comment: Science White Paper submitted to the Astro2020 Decadal Survey. arXiv admin note: substantial text overlap with arXiv:1810.06609

We investigate the multi-wavelength properties of host galaxies of 3701 X-ray-selected active galactic nuclei (AGNs) out to z~5 in the Chandra-COSMOS Legacy Survey. Thanks to the extensive multi-wavelength photometry available in the COSMOS field, we derive AGN luminosities, host stellar masses, and star formation rates (SFRs) via a multi-component SED fitting technique. Type 1 and Type 2 AGNs follow the same intrinsic Lx-L6um relation, suggesting that mid-infrared emission is a reasonably good measure of the AGN accretion power regardless of obscuration. We find that there is a strong increase in Type 1 AGN fraction toward higher AGN luminosity, possibly due to the fact that Type 1 AGNs tend to be hosted by more massive galaxies. The AGN luminosity and SFR are consistent with an increase toward high stellar mass, while both the Mstellar-dependence is weaker towards the high-mass end, which could be interpreted as a consequence of quenching both star formation and AGN activity in massive galaxies. AGN host galaxies tend to have SFRs that are consistent with normal star-forming galaxies, independent of AGN luminosities. We confirm that black hole accretion rate and SFR are correlated up to z~5, when forming stars. The majority (~73%) of our AGN sample are faint in the far-infrared, implying that the moderate-luminosity AGNs seem to be still active after the star formation is suppressed. It is not certain whether AGN activity plays a role in quenching the star formation. We conclude that both AGN activity and star formation might be more fundamentally related to host stellar mass., Comment: 30 pages, 18 figures, accepted for publication in ApJ

We present the results of a search for galaxy clusters and groups in the $\sim2$ square degree of the COSMOS field using all available X-ray observations from the XMM-Newton and Chandra observatories. We reach an X-ray flux limit of $3\times10^{-16}\;ergs\;cm^{-2}\;s^{-1}$ in 0.5--2 keV range, and identify 247 X-ray groups with $M_{200c}=8\times10^{12}-3\times10^{14}\;M_{\odot}$ at a redshift range of $0.08\leq z<1.53$, using the multiband photometric redshift and the master spectroscopic redshift catalogues of the COSMOS. The X-ray centres of groups are determined using high-resolution Chandra imaging. We investigate the relations between the offset of the brightest group galaxies (BGGs) from halo X-ray centre and group properties and compare with predictions from semi-analytic models and hydrodynamical simulations. We find that BGG offset decreases with both increasing halo mass and decreasing redshift with no strong dependence on the X-ray flux and SNR. We show that the BGG offset decreases as a function of increasing magnitude gap with no considerable redshift dependent trend. The stellar mass of BGGs in observations extends over a wider dynamic range compared to model predictions. At $z<0.5$, the central dominant BGGs become more massive than those with large offsets by up to 0.3dex, in agreement with model prediction. The observed and predicted lognormal scatter in the stellar mass of both low- and large-offset BGGs at fixed halo mass is $\sim0.3$dex., Comment: Accepted 2018 November 16 (MNRAS)

A primary aim of the Nuclear Spectroscopic Telescope Array (NuSTAR) mission is to find and characterize heavily obscured Active Galactic Nuclei (AGNs). Based on mid-infrared photometry from the Wide-Field Infrared Survey Explorer (WISE) and optical photometry from the Sloan Digital Sky Surveys, we have selected a large population of luminous obscured AGN (i.e., obscured quasars). Here we report NuSTAR observations of four WISE-selected heavily obscured quasars for which we have optical spectroscopy from the Southern African Large Telescope and W. M. Keck Observatory. Optical diagnostics confirm that all four targets are AGNs. With NuSTAR hard X-ray observations, three of the four objects are undetected, while the fourth has a marginal detection. We confirm that these objects have observed hard X-ray (10-40 keV) luminosities at or below ~ 10^43 erg s^-1. We compare X-ray and IR luminosities to obtain estimates of the hydrogen column densities (N_H) based on the suppression of the hard X-ray emission. We estimate N_H of these quasars to be at or larger than 10^25 cm^-2, confirming that WISE and optical selection can identify very heavily obscured quasars that may be missed in X-ray surveys, and do not contribute significantly to the cosmic X-ray background. From the optical Balmer decrements, we found that our three extreme obscured targets lie in highly reddened host environments. This galactic extinction is not adequate to explain the more obscured AGN, but it may imply a different scale of obscuration in the galaxy., Comment: Accepted for publication in the Astrophysical Journal, 11 pages, 10 figures

Gravitational-wave (GW) and gravitational slingshot recoil kicks, which are natural products of SMBH evolution in merging galaxies, can produce active galactic "nuclei" that are offset from the centers of their host galaxies. Detections of offset AGN would provide key constraints on SMBH binary mass and spin evolution and on GW event rates. Although numerous offset AGN candidates have been identified, none have been definitively confirmed. The ngVLA offers unparalleled capabilities to identify and confirm candidate offset AGN from sub-parsec to kiloparsec scales, opening a new avenue for multi-messenger studies in the dawn of low-frequency GW astronomy., Comment: To be published in the ASP Monograph Series, "Science with a Next-Generation VLA", ed. E. J. Murphy (ASP, San Francisco, CA)

We report on the results of active galactic nuclei (AGNs) detection by NuSTAR performed in three extragalactic survey fields (COSMOS, UDS, ECDFS) in three hard bands, namely H1 (8-16 keV), H2 (16-24 keV) and VH (35-55 keV). The aggregated area of the surveys is $\sim 2.7$ deg$^2$. While a large number of sources is detected in the H1 band (72 at the $97\%$ level of reliability), the H2 band directly probing close to the peak of the Cosmic X-ray Background (CXB) returns four significant detections, and two tentative, although not significant, detections are found in the VH band. All the sources detected above 16 keV are also detected at lower energies. We compute the integral number counts for sources in such bands, which show broad consistency with population synthesis models of the CXB. We furthermore identify two Compton-thick AGNs, one in the COSMOS field, associated with a hard and faint Chandra source, and one in the UDS field, never detected in the X-ray band before. Both sources are at the same redshift $z \sim 1.25$, which shifts their Compton-hump into the H1 band, and were previously missed in the usually employed NuSTAR bands, confirming the potential of using the H1 band to discover obscured AGNs at $z > 1$ in deep surveys., Comment: 12 pages, 5 figures, accepted for publication in The Astrophysical Journal

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., Comment: 15 pages, 11 figures, published in the Astrophysics Journal Supplement Series

We present an X-ray stacking analysis of $\sim$75,000 star-forming galaxies between $0.11.3$ exhibit a significant X-ray excess indicating nuclear emission. The AGN emission is strongly correlated with stellar mass but does not exhibit an additional correlation with SFR. The hardness ratios of the high-redshift stacks indicate that the AGN are substantially obscured ($N_{\mathrm{H}}\sim10^{23}$ cm$^{-2}$). These obscured AGN are not identified by IRAC color selection and have $L_X\sim10^{41}-10^{43}$ erg s$^{-1}$, consistent with accretion at an Eddington rate of $\sim10^{-3}$ onto $10^7-10^8 M_{\odot}$ black holes. Combining our results with other X-ray studies suggests that AGN obscuration depends on stellar mass and an additional variable, possibly the Eddington rate., Comment: 24 pages, 12 figures, 3 tables (2 tables available in full online or by request)

SPHEREx is a proposed NASA MIDEX mission selected for Phase A study. SPHEREx would carry out the first all-sky spectral survey in the near infrared. At the end of its two-year mission, SPHEREx would obtain 0.75-to-5$\mu$m spectra of every 6.2 arcsec pixel on the sky, with spectral resolution R>35 and a 5-$\sigma$ sensitivity AB$>$19 per spectral/spatial resolution element. More details concerning SPHEREx are available at http://spherex.caltech.edu. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. Though these three themes are undoubtedly compelling, they are far from exhausting the scientific output of SPHEREx. Indeed, SPHEREx would create a unique all-sky spectral database including spectra of very large numbers of astronomical and solar system targets, including both extended and diffuse sources. These spectra would enable a wide variety of investigations, and the SPHEREx team is dedicated to making the data available to the community to enable these investigations, which we refer to as Legacy Science. To that end, we have sponsored two workshops for the general scientific community to identify the most interesting Legacy Science themes and to ensure that the SPHEREx data products are responsive to their needs. In February of 2016, some 50 scientists from all fields met in Pasadena to develop these themes and to understand their implications for the SPHEREx mission. The 2016 workshop highlighted many synergies between SPHEREx and other contemporaneous astronomical missions, facilities, and databases. Consequently, in January 2018 we convened a second workshop at the Center for Astrophysics in Cambridge to focus specifically on these synergies. This white paper reports on the results of the 2018 SPHEREx workshop., Comment: 50 pages, 24 figures, more details at http://spherex.caltech.edu

We present a study of the connection between black hole accretion, star formation, and galaxy morphology at z<=2.5. We focus on active galactic nuclei (AGNs) selected by their mid-IR power-law emission. By fitting optical to far-IR photometry with state-of-the-art spectral energy distribution (SED) techniques, we derive stellar masses, star formation rates, dust properties, and AGN contributions in galaxies over the whole COSMOS field. We find that obscured AGNs lie within or slightly above the star-forming sequence. We confirm our previous finding about compact host galaxies of obscured AGNs at z~1, and find that galaxies with 20-50% AGN contributions tend to have smaller sizes, by ~25-50%, compared to galaxies without AGNs. Furthermore, we find that a high merger fraction of up to 0.5 is appropriate for the most luminous (log (LIR/Lsun) ~ 12.5) AGN hosts and non-AGN galaxies, but not for the whole obscured AGN sample. Moreover, merger fraction depends on the total and star-forming infrared luminosity, rather than the decomposed AGN infrared luminosity. Our results suggest that major mergers are not the main driver of AGN activity, and therefore obscured AGNs might be triggered by internal mechanisms, such as secular processes, disk instabilities, and compaction in a particular evolutionary stage. We make the SED modeling results publicly available., Comment: Accepted for publication in ApJS