Showing 1–50 of 190 results for author: Hickox, R

Optical Spectropolarimetric Variability Properties in Blazars PKS 0637-75 and PKS 1510-089

Authors: Stephanie A. Podjed, Ryan C. Hickox, Jedidah C. Isler, Markus Böttcher, Hester M. Schutte

Abstract: Spectropolarimetry is a powerful tool to investigate the central regions of active galactic nuclei (AGNs) as polarization signatures are key to probing magnetic field structure, evolution, and the physics of particle acceleration in jets. Optical linear polarization of blazars is typically greater than a few percent, indicating the emission is dominated by nonthermal synchrotron radiation, while p… ▽ More Spectropolarimetry is a powerful tool to investigate the central regions of active galactic nuclei (AGNs) as polarization signatures are key to probing magnetic field structure, evolution, and the physics of particle acceleration in jets. Optical linear polarization of blazars is typically greater than a few percent, indicating the emission is dominated by nonthermal synchrotron radiation, while polarization less than a few percent is common for other type 1 AGNs. We present a spectropolarimetric study of PKS 0637-75 and PKS 1510-089 to determine how the head-on orientation of a jet and dominant emission processes influence polarimetric variations in the broad lines and continuum. Observations were obtained biweekly from the Robert Stobie Spectrograph on the Southern African Large Telescope. Variability in the continuum polarization is detected for both PKS 0637-75 and PKS 1510-089, with a total average level of 2.5% +/- 0.1% and 7.5% +/- 0.1%, respectively. There is no clear polarization in the broad Balmer emission lines and weak polarization in Mg II as the average level across all observations is 0.2% +/- 0.1% for Hbeta, 0.2% +/- 0.3% for Hgamma, and 0.6% +/- 0.2% for Mg II. We find that polarization measurements confirm the conclusions drawn from spectral energy distribution modeling of the disk-jet contributions to the emission as optical polarization and time variability for PKS 0637-75 are shown to be dominated by accretion disk emission while those of PKS 1510-089 are due to both disk and jet emission, with greater jet contribution during flaring states. △ Less

Submitted 25 June, 2024; originally announced June 2024.

Comments: 12 pages, 9 figures. Accepted for publication in The Astrophysical Journal

Journal ref: ApJ 968 (2024) 130 (12pp)

Exploring the AGN Fraction of a Sample of JWST's Little Red Dots at $5 < z < 8$: Overmassive Black Holes Are Strongly Favored

Authors: Emmanuel Durodola, Fabio Pacucci, Ryan C. Hickox

Abstract: JWST is revolutionizing our view of the early Universe by pushing the boundaries of detectable galaxies and black holes in redshift (upward) and mass (downward). The Little Red Dots (LRDs), detected by several surveys at $z > 4$, present a significant interpretational challenge, as their spectral energy distributions can mimic both AGN and stellar population templates. This study analyzes eight LR… ▽ More JWST is revolutionizing our view of the early Universe by pushing the boundaries of detectable galaxies and black holes in redshift (upward) and mass (downward). The Little Red Dots (LRDs), detected by several surveys at $z > 4$, present a significant interpretational challenge, as their spectral energy distributions can mimic both AGN and stellar population templates. This study analyzes eight LRDs from the CEERS survey, utilizing NIRCam and MIRI photometry. By performing SED fitting across a vast parameter space, we explore a broad range of AGN fractions, defined as the ratio of the monochromatic luminosities (AGN, galaxy, and dust) at 0.53 $μm$ rest-frame. We find that the SEDs of all the LRDs investigated are consistent with having significant AGN contributions, with the best-fitting AGN fractions ranging between 40% and 85%. Moreover, assuming these LRDs do indeed host AGN, we can place limits on their black hole masses using the inferred AGN bolometric luminosities and adopting the Eddington limit. We find that, independent of the specific AGN fraction adopted, the LRDs' black holes are significantly overmassive relative to their host galaxies -- by $\sim 2$ dex, and up to $\sim 4$ dex in the most extreme cases -- compared to the local $M_{\bullet} - M_{\star}$ relation. The presence of overmassive black holes in the high-$z$ Universe may provide the strongest evidence yet of heavy black hole seeding occurring during the cosmic dark ages. △ Less

Submitted 14 June, 2024; originally announced June 2024.

Comments: Submitted to The Astrophysical Journal. 9 pages, 3 figures

X-ray View of Little Red Dots: Do They Host Supermassive Black Holes?

Authors: Tonima Tasnim Ananna, Ákos Bogdán, Orsolya E. Kovács, Priyamvada Natarajan, Ryan C. Hickox

Abstract: The discovery of Little Red Dots (LRDs) -- a population of compact, high-redshift, dust-reddened galaxies -- is one of the most surprising results from JWST. However, the nature of LRDs is still debated: does the near-infrared emission originate from accreting supermassive black holes (SMBHs), or intense star formation? In this work, we utilize ultra-deep Chandra observations and study LRDs residi… ▽ More The discovery of Little Red Dots (LRDs) -- a population of compact, high-redshift, dust-reddened galaxies -- is one of the most surprising results from JWST. However, the nature of LRDs is still debated: does the near-infrared emission originate from accreting supermassive black holes (SMBHs), or intense star formation? In this work, we utilize ultra-deep Chandra observations and study LRDs residing behind the lensing galaxy cluster, Abell~2744. We probe the X-ray emission from individual galaxies but find that they remain undetected and provide SMBH mass upper limits of $\lesssim(1.5-16)\times10^{6}~\rm{M_{\odot}}$ assuming Eddington limited accretion. To increase the signal-to-noise ratios, we conduct a stacking analysis of the full sample with a total lensed exposure time of $\approx87$~Ms. We also bin the galaxies based on their stellar mass, lensing magnification, and detected broad-line H$α$ emission. For the LRDs exhibiting broad-line H$α$ emission, there is a hint of a stacked signal ($\sim2.6σ$), corresponding to a SMBH mass of $\sim3.2\times10^{6}~\rm{M_{\odot}}$. Assuming unobscured, Eddington-limited accretion, this BH mass is at least 1.5 orders of magnitude lower than that inferred from virial mass estimates using JWST spectra. Given galaxy-dominated stellar mass estimates, our results imply that LRDs do not host over-massive SMBHs and/or accrete at a few percent of their Eddington limit. However, alternative stellar mass estimates may still support that LRDs host over-massive BHs. The significant discrepancy between the JWST and Chandra data hints that the scaling relations used to infer the SMBH mass from the H$α$ line and virial relations may not be applicable for high-redshift LRDs. △ Less

Submitted 12 June, 2024; v1 submitted 29 April, 2024; originally announced April 2024.

Comments: Accepted by ApJL; DOI: 10.3847/2041-8213/ad5669

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

Recommendations for Early Definition Science with the Nancy Grace Roman Space Telescope

Authors: Robyn E. Sanderson, Ryan Hickox, Christopher M. Hirata, Matthew J. Holman, Jessica R. Lu, Ashley Villar

Abstract: The Nancy Grace Roman Space Telescope (Roman), NASA's next flagship observatory, has significant mission time to be spent on surveys for general astrophysics in addition to its three core community surveys. We considered what types of observations outside the core surveys would most benefit from early definition, given 700 hours of mission time in the first two years of Roman's operation. We recom… ▽ More The Nancy Grace Roman Space Telescope (Roman), NASA's next flagship observatory, has significant mission time to be spent on surveys for general astrophysics in addition to its three core community surveys. We considered what types of observations outside the core surveys would most benefit from early definition, given 700 hours of mission time in the first two years of Roman's operation. We recommend that a survey of the Galactic plane be defined early, based on the broad range of stakeholders for such a survey, the added scientific value of a first pass to obtain a baseline for proper motions complementary to Gaia's, and the significant potential synergies with ground-based surveys, notably the Legacy Survey of Space and Time (LSST) on Rubin. We also found strong motivation to follow a community definition process for ultra-deep observations with Roman. △ Less

Submitted 22 April, 2024; originally announced April 2024.

Comments: Final report of the committee convened by NASA to assess community proposals for early definition science with Roman

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

Binary black hole mergers and intermediate-mass black holes in dense star clusters with collisional runaways

Authors: Rujuta A. Purohit, Giacomo Fragione, Frederic A. Rasio, Grayson C. Petter, Ryan C. Hickox

Abstract: Intermediate-mass black holes (IMBHs) are believed to be the missing link between the supermassive black holes (BHs) found at the centers of massive galaxies and BHs formed through stellar core collapse. One of the proposed mechanisms for their formation is a collisional runaway process in high-density young star clusters, where an unusually massive object forms through repeated stellar collisions… ▽ More Intermediate-mass black holes (IMBHs) are believed to be the missing link between the supermassive black holes (BHs) found at the centers of massive galaxies and BHs formed through stellar core collapse. One of the proposed mechanisms for their formation is a collisional runaway process in high-density young star clusters, where an unusually massive object forms through repeated stellar collisions and mergers, eventually collapsing to form an IMBH. This seed IMBH could then grow further through binary mergers with other stellar-mass BHs. Here we investigate the gravitational-wave (GW) signals produced during these later IMBH--BH mergers. We use a state-of-the-art semi-analytic approach to study the stellar dynamics and to characterize the rates and properties of IMBH--BH mergers. We also study the prospects for detection of these mergers by current and future GW observatories, both space-based (LISA) and ground-based (LIGO Voyager, Einstein Telescope, and Cosmic Explorer). We find that most of the merger signals could be detected, with some of them being multi-band sources. Therefore, GWs represent a unique tool to test the collisional runaway scenario and to constrain the population of dynamically assembled IMBHs. △ Less

Submitted 13 March, 2024; v1 submitted 12 January, 2024; originally announced January 2024.

Comments: Accepted for publication in the Astronomical Journal

The High Energy X-ray Probe (HEX-P): The Future of Hard X-ray Dual AGN Science

Authors: Ryan W. Pfeifle, Peter G. Boorman, Kimberly A. Weaver, Johannes Buchner, Francesca Civano, Kristin Madsen, Daniel Stern, Núria Torres-Albà, Emanuele Nardini, Claudio Ricci, Stefano Marchesi, D. R. Ballantyne, Dominic Sicilian, Chien-Ting Chen, Elias Kammoun, Ryan C. Hickox, Javier A. García, Labani Mallick

Abstract: 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 lat… ▽ More 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. △ Less

Submitted 9 November, 2023; originally announced November 2023.

Comments: Under review for Frontiers in Astronomy and Space Sciences

The High Energy X-ray Probe (HEX-P): The Circum-nuclear Environment of Growing Supermassive Black Holes

Authors: P. G. Boorman, N. Torres-Albà, A. Annuar, S. Marchesi, R. Pfeifle, D. Stern, F. Civano, M. Baloković, J. Buchner, C. Ricci, D. M. Alexander, W. N. Brandt, M. Brightman, C. T. Chen, S. Creech, P. Gandhi, J. A. García, F. Harrison, R. Hickox, E. Kammoun, S. LaMassa, G. Lanzuisi, L. Marcotulli, K. Madsen, G. Matt , et al. (10 additional authors not shown)

Abstract: Ever since the discovery of the first Active Galactic Nuclei (AGN), substantial observational and theoretical effort has been invested into understanding how massive black holes have evolved across cosmic time. Circum-nuclear obscuration is now established as a crucial component, with almost every AGN observed known to display signatures of some level of obscuration in their X-ray spectra. But des… ▽ More Ever since the discovery of the first Active Galactic Nuclei (AGN), substantial observational and theoretical effort has been invested into understanding how massive black holes have evolved across cosmic time. Circum-nuclear obscuration is now established as a crucial component, with almost every AGN observed known to display signatures of some level of obscuration in their X-ray spectra. But despite more than six decades of effort, substantial open questions remain: How does the accretion power impact the structure of the circum-nuclear obscurer? What are the dynamical properties of the obscurer? Can dense circum-nuclear obscuration exist around intrinsically weak AGN? How many intermediate mass black holes occupy the centers of dwarf galaxies? In this paper, we showcase a number of next-generation prospects attainable with the High Energy X-ray Probe (https://hexp.org) to contribute towards solving these questions in the 2030s. The uniquely broad (0.2--80 keV) and strictly simultaneous X-ray passband of HEX-P makes it ideally suited for studying the temporal co-evolution between the central engine and circum-nuclear obscurer. Improved sensitivities and reduced background will enable the development of spectroscopic models complemented by current and future multi-wavelength observations. We show that the angular resolution of HEX-P both below and above 10 keV will enable the discovery and confirmation of accreting massive black holes at both low accretion power and low black hole masses even when concealed by thick obscuration. In combination with other next-generation observations of the dusty hearts of nearby galaxies, HEX-P will hence be pivotal in paving the way towards a complete picture of black hole growth and galaxy co-evolution. △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: To be submitted for peer review to a HEX-P Collection in Frontiers of Astronomy and Space Science. Constructive comments welcome!

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

Ionized Gas Extended Over 40 kpc in an Odd Radio Circle Host Galaxy

Authors: Alison L. Coil, Serena Perrotta, David S. N. Rupke, Cassandra Lochhaas, Christy A. Tremonti, Aleks Diamond-Stanic, Drummond Fielding, James Geach, Ryan C. Hickox, John Moustakas, Gregory H. Rudnick, Paul Sell, Kelly E. Whalen

Abstract: A new class of extragalactic astronomical sources discovered in 2021, named Odd Radio Circles (ORCs, Norris et al. 2021), are large rings of faint, diffuse radio continuum emission spanning ~1 arcminute on the sky. Galaxies at the centers of several ORCs have photometric redshifts of z~0.3-0.6, implying physical scales of several 100 kiloparsecs in diameter for the radio emission, the origin of wh… ▽ More A new class of extragalactic astronomical sources discovered in 2021, named Odd Radio Circles (ORCs, Norris et al. 2021), are large rings of faint, diffuse radio continuum emission spanning ~1 arcminute on the sky. Galaxies at the centers of several ORCs have photometric redshifts of z~0.3-0.6, implying physical scales of several 100 kiloparsecs in diameter for the radio emission, the origin of which is unknown. Here we report spectroscopic data on an ORC including strong [OII] emission tracing ionized gas in the central galaxy of ORC4 at z=0.4512. The physical extent of the [OII] emission is ~40 kpc in diameter, larger than expected for a typical early-type galaxy (Pandya et al, 2017) but an order of magnitude smaller than the large-scale radio continuum emission. We detect a ~200 km/s velocity gradient across the [OII] nebula, as well as a high velocity dispersion of ~180 km/s. The [OII] equivalent width (EW, ~50 Ang) is extremely high for a quiescent galaxy. The morphology, kinematics, and strength of the [OII] emission are consistent with the infall of shock ionized gas near the galaxy, following a larger-scale, outward moving shock driven by a galactic wind. Both the extended optical and radio emission, while observed on very different scales, may therefore result from the same dramatic event. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: 7 figures, accepted to Nature

Obscuration beyond the nucleus: infrared quasars can be buried in extreme compact starbursts

Authors: Carolina Andonie, David M. Alexander, Claire Greenwell, Annagrazia Puglisi, Brivael Laloux, Alba V. Alonso-Tetilla, Gabriela Calistro Rivera, Chris Harrison, Ryan C. Hickox, Melanie Kaasinen, Andrea Lapi, Iván E. López, Grayson Petter, Cristina Ramos Almeida, David J. Rosario, Francesco Shankar, Carolin Villforth

Abstract: In the standard quasar model, the accretion disk obscuration is due to the canonical dusty torus. Here, we argue that a substantial part of the quasar obscuration can come from the interstellar medium (ISM) when the quasars are embedded in compact starbursts. We use an obscuration-unbiased sample of 578 infrared (IR) quasars at $z\approx 1-3$ and archival ALMA submillimeter host galaxy sizes to in… ▽ More In the standard quasar model, the accretion disk obscuration is due to the canonical dusty torus. Here, we argue that a substantial part of the quasar obscuration can come from the interstellar medium (ISM) when the quasars are embedded in compact starbursts. We use an obscuration-unbiased sample of 578 infrared (IR) quasars at $z\approx 1-3$ and archival ALMA submillimeter host galaxy sizes to investigate the ISM contribution to the quasar obscuration. We calculate SFR and ISM column densities for the IR quasars and a control sample of submillimeter galaxies (SMGs) not hosting quasar activity and show that: (1) the quasar obscured fraction is constant up to $\rm SFR\approx 300 \: M_{\odot} \: yr^{-1}$, and then increases towards higher SFR, suggesting that the ISM obscuration plays a significant role in starburst host galaxies, and (2) at $\rm SFR\gtrsim 300 \: M_{\odot} \: yr^{-1}$, the SMGs and IR quasars have similarly compact submillimeter sizes ($R_{\rm e}\approx 0.5-3\rm \: kpc$) and, consequently, the ISM can heavily obscure the quasar, even reaching Compton-thick ($N_{\rm H}>10^{24} \rm \: cm^{-2}$) levels in extreme cases. Based on our results, we infer that $\approx 10-30\%$ of the IR quasars with $\rm SFR\gtrsim 300 \: M_{\odot} \: yr^{-1}$ are obscured solely by the ISM. △ Less

Submitted 10 October, 2023; v1 submitted 3 October, 2023; originally announced October 2023.

Comments: Accepted for publication in MNRAS Letters

The Kiloparsec Scale Influence of the AGN in NGC 1068 with SALT RSS Fabry-Pérot Spectroscopy

Authors: Raphael E. Hviding, Ryan Hickox, P. Väisänen, Rajin Ramphul, Kevin Hainline

Abstract: We present Fabry-Pérot (FP) imaging and longslit spectroscopy of the nearby Seyfert II galaxy NGC 1068 using the Robert Stobie Spectrograph (RSS) on the Southern African Large Telescope (SALT) to observe the impact of the central Active Galactic Nucleus (AGN) on the ionized gas in the galaxy on kiloparsec scales. With SALT RSS FP we are able to observe the H$α$+[N II] emission line complex over a… ▽ More We present Fabry-Pérot (FP) imaging and longslit spectroscopy of the nearby Seyfert II galaxy NGC 1068 using the Robert Stobie Spectrograph (RSS) on the Southern African Large Telescope (SALT) to observe the impact of the central Active Galactic Nucleus (AGN) on the ionized gas in the galaxy on kiloparsec scales. With SALT RSS FP we are able to observe the H$α$+[N II] emission line complex over a $\sim$2.6 arcmin$^2$ field of view. Combined with the longslit observation, we demonstrate the efficacy of FP spectroscopy for studying nearby Type II Seyfert galaxies and investigate the kiloparsec-scale ionized gas in NGC 1068. We confirm the results of previous work from the TYPHOON/Progressive Integral Step Method (PrISM) survey that the kiloparsec-scale ionized features in NGC 1068 are driven by AGN photoionization. We analyze the spatial variation of the AGN intensity to put forward an explanation for the shape and structure of the kiloparsec-scale ionization features. Using a toy model, we suggest the ionization features may be understood as a light-echo from a burst of enhanced AGN activity $\sim$2000 years in the past. △ Less

Submitted 12 July, 2023; originally announced July 2023.

Comments: 18 pages, 9 figures. Accepted for publication in The Astronomical Journal

Measuring the Cosmic X-ray Background in 3-20keV with Straylight from NuSTAR

Authors: Steven Rossland, Daniel Wik, Brian Grefenstette, Nico Cappelluti, Francesca Civano, Fabio Gastaldello, Roberto Gilli, Fiona Harrison, Ann Hornschemeier, Ryan Hickox, Roman Krivonos, Kristin Madsen, Silvano Molendi, Andrew Ptak, Daniel Stern, Andreas Zoglauer

Abstract: 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… ▽ More 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. △ Less

Submitted 19 April, 2023; v1 submitted 16 April, 2023; originally announced April 2023.

Comments: 23 pages, 15 figures, 4 tables

The Most Obscured AGNs in the XMM-SERVS Fields

Authors: Wei Yan, W. N. Brandt, Fan Zou, Shifu Zhu, Chien-Ting J. Chen, Ryan C. Hickox, Bin Luo, Qingling Ni, David M. Alexander, Franz E. Bauer, Cristian Vignali, Fabio Vito

Abstract: We perform X-ray spectral analyses to derive characteristics (e.g., column density, X-ray luminosity) of $\approx$10,200 active galactic nuclei (AGNs) in the XMM-Spitzer Extragalactic Representative Volume Survey (XMM-SERVS), which was designed to investigate the growth of supermassive black holes across a wide dynamic range of cosmic environments. Using physical torus models (e.g., Borus02) and a… ▽ More We perform X-ray spectral analyses to derive characteristics (e.g., column density, X-ray luminosity) of $\approx$10,200 active galactic nuclei (AGNs) in the XMM-Spitzer Extragalactic Representative Volume Survey (XMM-SERVS), which was designed to investigate the growth of supermassive black holes across a wide dynamic range of cosmic environments. Using physical torus models (e.g., Borus02) and a Bayesian approach, we uncover 22 representative Compton-thick (CT; $N_{\rm H} \;>\; 1.5\times10^{24}\; \rm cm^{-2}$) AGN candidates with good signal-to-noise ratios as well as a large sample of 136 heavily obscured AGNs. We also find an increasing CT fraction (\fct ) from low ($z<0.75$) to high ($z>0.75$) redshift. Our CT candidates tend to show hard X-ray spectral shapes and dust extinction in their SED fits, which may shed light on the connection between AGN obscuration and host-galaxy evolution. △ Less

Submitted 12 April, 2023; originally announced April 2023.

Comments: 12 pages, 9 figures, accepted for publication in ApJ

A High Fraction of Heavily X-ray Obscured Active Galactic Nuclei

Authors: Christopher M. Carroll, Tonima T. Ananna, Ryan C. Hickox, Alberto Masini, Roberto J. Assef, Daniel Stern, Chien-Ting J. Chen, Lauranne Lanz

Abstract: We present new estimates on the fraction of heavily X-ray obscured, Compton-thick (CT) active galactic nuclei (AGNs) out to a redshift of $z \leq$ 0.8. From a sample of 540 AGNs selected by mid-IR (MIR) properties in observed X-ray survey fields, we forward model the observed-to-intrinsic X-ray luminosity ratio ($R_{L_{\text{X}}}$) with a Markov chain Monte Carlo (MCMC) simulation to estimate the… ▽ More We present new estimates on the fraction of heavily X-ray obscured, Compton-thick (CT) active galactic nuclei (AGNs) out to a redshift of $z \leq$ 0.8. From a sample of 540 AGNs selected by mid-IR (MIR) properties in observed X-ray survey fields, we forward model the observed-to-intrinsic X-ray luminosity ratio ($R_{L_{\text{X}}}$) with a Markov chain Monte Carlo (MCMC) simulation to estimate the total fraction of CT AGNs ($f_{\text{CT}}$), many of which are missed in typical X-ray observations. We create model $N_{\text{H}}$ distributions and convert these to $R_{L_{\text{X}}}$ using a set of X-ray spectral models. We probe the posterior distribution of our models to infer the population of X-ray non-detected sources. From our simulation we estimate a CT fraction of $f_{\text{CT}}$ = $\text{0.555}^{+\text{0.037}}_{-\text{0.032}}$. We perform an X-ray stacking analysis for sources in Chandra X-ray Observatory fields and find that the expected soft (0.5-2 keV) and hard (2-7 keV) observed fluxes drawn from our model to be within 0.48 and 0.12 dex of our stacked fluxes, respectively. Our results suggests at least 50% of all MIR-selected AGNs, possibly more, are Compton-thick ($N_{\text{H}} \gtrsim$ 10$^{\text{24}}$ cm$^{-\text{2}}$), which is in excellent agreement with other recent work using independent methods. This work indicates that the total number of AGNs is higher than can be identified using X-ray observations alone, highlighting the importance of a multiwavelength approach. A high $f_{\text{CT}}$ also has implications for black hole (BH) accretion physics and supports models of BH and galaxy co-evolution that include periods of heavy obscuration. △ Less

Submitted 15 March, 2023; originally announced March 2023.

Comments: 14 pages, 6 figures, 1 table, plus appendix figures. Accepted for publication in ApJ

Kinematics, Structure, and Mass Outflow Rates of Extreme Starburst Galactic Outflows

Authors: Serena Perrotta, Alison L. Coil, David S. N. Rupke, Christy A. Tremonti, Julie D. Davis, Aleksandar M. Diamond-Stanic, James E. Geach, Ryan C. Hickox, John Moustakas, Gregory H. Rudnick, Paul H. Sell, Cameren N. Swiggum, Kelly E. Whalen

Abstract: We present results on the properties of extreme gas outflows in massive ($\rm M_* \sim$10$^{11} \ \rm M_{\odot}$), compact, starburst ($\rm SFR \sim$$200 \, \rm M_{\odot} \ yr^{-1}$) galaxies at z = $0.4-0.7$ with very high star formation surface densities ($\rm Σ_{SFR} \sim$$2000 \,\rm M_{\odot} \ yr^{-1} \ kpc^{-2}$). Using optical Keck/HIRES spectroscopy of 14 HizEA starburst galaxies we identi… ▽ More We present results on the properties of extreme gas outflows in massive ($\rm M_* \sim$10$^{11} \ \rm M_{\odot}$), compact, starburst ($\rm SFR \sim$$200 \, \rm M_{\odot} \ yr^{-1}$) galaxies at z = $0.4-0.7$ with very high star formation surface densities ($\rm Σ_{SFR} \sim$$2000 \,\rm M_{\odot} \ yr^{-1} \ kpc^{-2}$). Using optical Keck/HIRES spectroscopy of 14 HizEA starburst galaxies we identify outflows with maximum velocities of $820 - 2860$ \kmps. High-resolution spectroscopy allows us to measure precise column densities and covering fractions as a function of outflow velocity and characterize the kinematics and structure of the cool gas outflow phase (T $\sim$10$^4$ K). We find substantial variation in the absorption profiles, which likely reflects the complex morphology of inhomogeneously-distributed, clumpy gas and the intricacy of the turbulent mixing layers between the cold and hot outflow phases. There is not a straightforward correlation between the bursts in the galaxies' star formation histories and their wind absorption line profiles, as might naively be expected for starburst-driven winds. The lack of strong \mgii \ absorption at the systemic velocity is likely an orientation effect, where the observations are down the axis of a blowout. We infer high mass outflow rates of $\rm \sim$50 $-$ 2200 $\rm M_{\odot} \, yr^{-1}$, assuming a fiducial outflow size of 5 kpc, and mass loading factors of $η\sim$5 for most of the sample. %with $η\sim$20 for two galaxies. While these values have high uncertainties, they suggest that starburst galaxies are capable of ejecting very large amounts of cool gas that will substantially impact their future evolution. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Comments: Accepted for publication in The Astrophysical Journal

The Ionization and Dynamics of the Makani Galactic Wind

Authors: David S. N. Rupke, Alison L. Coil, Serena Perrotta, Julie D. Davis, Aleksandar M. Diamond-Stanic, James E. Geach, Ryan C. Hickox, John Moustakas, Grayson C. Petter, Gregory H. Rudnick, Paul H. Sell, Christy A. Tremonti, Kelly E. Whalen

Abstract: The Makani galaxy hosts the poster child of a galactic wind on scales of the circumgalactic medium. It consists of a two-episode wind in which the slow, outer wind originated 400 Myr ago (Episode I; R_I = 20-50 kpc) and the fast, inner wind is 7 Myr old (Episode II; R_II = 0-20 kpc). While this wind contains ionized, neutral, and molecular gas, the physical state and mass of the most extended phas… ▽ More The Makani galaxy hosts the poster child of a galactic wind on scales of the circumgalactic medium. It consists of a two-episode wind in which the slow, outer wind originated 400 Myr ago (Episode I; R_I = 20-50 kpc) and the fast, inner wind is 7 Myr old (Episode II; R_II = 0-20 kpc). While this wind contains ionized, neutral, and molecular gas, the physical state and mass of the most extended phase--the warm, ionized gas--is unknown. Here we present Keck optical spectra of the Makani outflow. These allow us to detect hydrogen lines out to r = 30-40 kpc and thus constrain the mass, momentum, and energy in the wind. Many collisionally-excited lines are detected throughout the wind, and their line ratios are consistent with 200-400 km/s shocks that power the ionized gas, with v_shock = $σ$_wind. Combining shock models, density-sensitive line ratios, and mass and velocity measurements, we estimate that the ionized mass and outflow rate in the Episode II wind could be as high as that of the molecular gas: M_II(HII) ~ M_II(H_2) = (1-2)x10^9 Msun and dM/dt_II(HII) ~ dM/dt_II(H_2) = 170-250 Msun/yr. The outer wind has slowed, so that dM/dt_I(HII) ~ 10 Msun/yr, but it contains more ionized gas: M_I(HII) = 5x10^9 Msun. The momentum and energy in the recent Episode II wind imply a momentum-driven flow (p ``boost" ~ 7) driven by the hot ejecta and radiation pressure from the Eddington-limited, compact starburst. Much of the energy and momentum in the older Episode I wind may reside in a hotter phase, or lie further into the CGM. △ Less

Submitted 28 February, 2023; originally announced March 2023.

Comments: ApJ, in press

Host Dark Matter Halos of WISE-selected Obscured & Unobscured Quasars: Evidence for Evolution

Authors: Grayson C. Petter, Ryan C. Hickox, David M. Alexander, Adam D. Myers, James E. Geach, Kelly E. Whalen, Carolina P. Andonie

Abstract: Obscuration in quasars may arise from steep viewing angles along the dusty torus, or instead may represent a distinct phase of supermassive black hole growth. We test these scenarios by probing the host dark matter halo environments of $\sim 1.4$ million WISE-selected obscured and unobscured quasars at $\langle z \rangle = 1.4$ using angular clustering measurements as well as cross-correlation mea… ▽ More Obscuration in quasars may arise from steep viewing angles along the dusty torus, or instead may represent a distinct phase of supermassive black hole growth. We test these scenarios by probing the host dark matter halo environments of $\sim 1.4$ million WISE-selected obscured and unobscured quasars at $\langle z \rangle = 1.4$ using angular clustering measurements as well as cross-correlation measurements of quasar positions with the gravitational lensing of the cosmic microwave background (CMB). We interpret these signals within a halo occupation distribution (HOD) framework to conclude that obscured systems reside in more massive effective halos ($ \sim 10^{12.9} h^{-1} M_{\odot}$) than their unobscured counterparts ($ \sim 10^{12.6} h^{-1} M_{\odot}$), though we do not detect a difference in the satellite fraction. We find excellent agreement between the clustering and lensing analyses and show that this implies the observed difference is robust to uncertainties in the obscured quasar redshift distribution, highlighting the power of combining angular clustering and weak lensing measurements. This finding appears in tension with models that ascribe obscuration exclusively to orientation of the dusty torus along the line-of-sight, and instead may be consistent with the notion that some obscured quasars are attenuated by galaxy-scale or circumnuclear material during an evolutionary phase. △ Less

Submitted 1 February, 2023; originally announced February 2023.

Comments: Accepted to The Astrophysical Journal

Monitoring observations of SMC X-1's excursions (MOOSE)-II: A new excursion accompanies spin-up acceleration

Authors: Chin-Ping Hu, Kristen C. Dage, William I. Clarkson, McKinley Brumback, Philip A. Charles, Daryl Haggard, Ryan C. Hickox, Tatehiro Mihara, Arash Bahramian, Rawan Karam, Wasundara Athukoralalage, Diego Altamirano, Joey Neilsen, Jamie Kennea

Abstract: SMC X-1 is a high-mass X-ray binary showing superorbital modulation with an unstable period. Previous monitoring shows three excursion events in 1996--1998, 2005--2007, and 2014--2016. The superorbital period drifts from >60 days to <40 days and then evolves back during an excursion. Here we report a new excursion event of SMC X-1 in 2020--2021, indicating that the superorbital modulation has an u… ▽ More SMC X-1 is a high-mass X-ray binary showing superorbital modulation with an unstable period. Previous monitoring shows three excursion events in 1996--1998, 2005--2007, and 2014--2016. The superorbital period drifts from >60 days to <40 days and then evolves back during an excursion. Here we report a new excursion event of SMC X-1 in 2020--2021, indicating that the superorbital modulation has an unpredictable, chaotic nature. We trace the spin-period evolution and find that the spin-up rate accelerated one year before the onset of this new excursion, which suggests a possible inside-out process connecting the spin-up acceleration and the superorbital excursion. This results in a deviation of the spin period residual, similar to the behaviour of the first excursion in 1996--1998. In further analysis of the pulse profile evolution, we find that the pulsed fraction shows a long-term evolution and may be connected to the superorbital excursion. These discoveries deepen the mystery of SMC X-1 because they cannot be solely interpreted by the warped disc model. Upcoming pointed observations and theoretical studies may improve our understanding of the detailed accretion mechanisms taking place. △ Less

Submitted 1 February, 2023; originally announced February 2023.

Comments: 7 pages, 3 figures, Accepted for publication in MNRAS

Black hole-galaxy co-evolution and the role of feedback

Authors: Pedro R. Capelo, Chiara Feruglio, Ryan C. Hickox, Francesco Tombesi

Abstract: Active galactic nuclei (AGN) are accreting supermassive black holes co-evolving with their host galaxies through a complex interplay of feeding and feedback. In this chapter, we first discuss AGN fuelling in galaxies, both in interacting and isolated systems, focusing on the role that instabilities have on the angular momentum budget of the gas. We then review observations and models of feedback t… ▽ More Active galactic nuclei (AGN) are accreting supermassive black holes co-evolving with their host galaxies through a complex interplay of feeding and feedback. In this chapter, we first discuss AGN fuelling in galaxies, both in interacting and isolated systems, focusing on the role that instabilities have on the angular momentum budget of the gas. We then review observations and models of feedback through AGN-driven winds from nuclear, sub-pc scales out to galactic and circumgalactic medium scales. We continue with an overview of surveys and statistical properties of the AGN population, before concluding with a discussion on the prospects for the future facilities, focusing in particular on Athena. △ Less

Submitted 1 November, 2022; originally announced November 2022.

Comments: Invited Chapter of Section "Active Galactic Nuclei in X- and Gamma-rays" (Section Editors: A. De Rosa and C. Vignali) of the "Handbook of X-ray and Gamma-ray Astrophysics" (Editors: C. Bambi and A. Santangelo), Springer Nature

Probing the Structure and Evolution of BASS AGN through Eddington Ratios

Authors: Tonima Tasnim Ananna, C. Megan Urry, Claudio Ricci, Priyamvada Natarajan, Ryan C. Hickox, Benny Trakhtenbrot, Ezequiel Treister, Anna K. Weigel, Yoshihiro Ueda, Michael J. Koss, F. E. Bauer, Matthew J. Temple, Mislav Balokovic, Richard Mushotzky, Connor Auge, David B. Sanders, Darshan Kakkad, Lia F. Sartori, Stefano Marchesi, Fiona Harrison, Daniel Stern, Kyuseok Oh, Turgay Caglar, Meredith C. Powell, Stephanie A. Podjed , et al. (1 additional authors not shown)

Abstract: We constrain the intrinsic Eddington ratio (\lamEdd ) distribution function for local AGN in bins of low and high obscuration (log NH <= 22 and 22 < log NH < 25), using the Swift-BAT 70-month/BASS DR2 survey. We interpret the fraction of obscured AGN in terms of circum-nuclear geometry and temporal evolution. Specifically, at low Eddington ratios (log lamEdd < -2), obscured AGN outnumber unobscure… ▽ More We constrain the intrinsic Eddington ratio (\lamEdd ) distribution function for local AGN in bins of low and high obscuration (log NH <= 22 and 22 < log NH < 25), using the Swift-BAT 70-month/BASS DR2 survey. We interpret the fraction of obscured AGN in terms of circum-nuclear geometry and temporal evolution. Specifically, at low Eddington ratios (log lamEdd < -2), obscured AGN outnumber unobscured ones by a factor of ~4, reflecting the covering factor of the circum-nuclear material (0.8, or a torus opening angle of ~ 34 degrees). At high Eddington ratios (\log lamEdd > -1), the trend is reversed, with < 30% of AGN having log NH > 22, which we suggest is mainly due to the small fraction of time spent in a highly obscured state. Considering the Eddington ratio distribution function of narrow-line and broad-line AGN from our prior work, we see a qualitatively similar picture. To disentangle temporal and geometric effects at high lamEdd, we explore plausible clearing scenarios such that the time-weighted covering factors agree with the observed population ratio. We find that the low fraction of obscured AGN at high lamEdd is primarily due to the fact that the covering factor drops very rapidly, with more than half the time is spent with < 10% covering factor. We also find that nearly all obscured AGN at high-lamEdd exhibit some broad-lines. We suggest that this is because the height of the depleted torus falls below the height of the broad-line region, making the latter visible from all lines of sight. △ Less

Submitted 15 October, 2022; originally announced October 2022.

Comments: Accepted by ApJL

The Space Density of Intermediate Redshift, Extremely Compact, Massive Starburst Galaxies

Authors: Kelly E. Whalen, Ryan C. Hickox, Alison L. Coil, Aleksandar M. Diamond-Stanic, James E. Geach, John Moustakas, Gregory H. Rudnick, David S. N. Rupke, Paul H. Sell, Christy A. Tremonti, Julie D. Davis, Serena Perrotta, Grayson C. Petter

Abstract: We present a measurement of the intrinsic space density of intermediate redshift ($z\sim0.5$), massive ($M_{*} \sim 10^{11} \ \text{M}_{\odot}$), compact ($R_{e} \sim 100$ pc) starburst ($Σ_{SFR} \sim 1000 \ \text{M}_{\odot} \ \text{yr}^{-1} \text{kpc}^{-1}$) galaxies with tidal features indicative of them having undergone recent major mergers. A subset of them host kiloparsec scale,… ▽ More We present a measurement of the intrinsic space density of intermediate redshift ($z\sim0.5$), massive ($M_{*} \sim 10^{11} \ \text{M}_{\odot}$), compact ($R_{e} \sim 100$ pc) starburst ($Σ_{SFR} \sim 1000 \ \text{M}_{\odot} \ \text{yr}^{-1} \text{kpc}^{-1}$) galaxies with tidal features indicative of them having undergone recent major mergers. A subset of them host kiloparsec scale, $>1000 \ \text{km}\ \text{s}^{-1}$ outflows and have little indication of AGN activity, suggesting that extreme star formation can be a primary driver of large-scale feedback. The aim for this paper is to calculate their space density so we can place them in a better cosmological context. We do this by empirically modeling the stellar populations of massive, compact starburst galaxies. We determine the average timescale for which galaxies that have recently undergone an extreme nuclear starburst would be targeted and included in our spectroscopically selected sample. We find that massive, compact starburst galaxies targeted by our criteria would be selectable for $\sim 148 ^{+27}_{-24}$ Myr and have an intrinsic space density $n_{\text{CS}} \sim (1.1^{+0.5}_{-0.3}) \times 10^{-6} \ \ \text{Mpc}^{-3}$. This space density is broadly consistent with our $z\sim0.5$ compact starbursts being the most extremely compact and star forming low redshift analogs of the compact star forming galaxies in the early Universe as well as them being the progenitors to a fraction of intermediate redshift post starburst and compact quiescent galaxies. △ Less

Submitted 27 September, 2022; originally announced September 2022.

Comments: 27 pages, 13 figures, accepted for publication in the Astrophysical Journal

NuSTAR observations of a heavily X-ray obscured AGN in the dwarf galaxy J144013+024744

Authors: Shrey. Ansh, Chien-Ting J. Chen, W. N. Brandt, Carol E. Hood, E. S. Kammoun, G. Lansbury, Stéphane Paltani, Amy E. Reines, C. Ricci, Douglas A. Swartz, Jonathan R. Trump, F. Vito, Ryan C. Hickox

Abstract: We present a multi-wavelength analysis of the dwarf Seyfert-2 galaxy J$144013+024744$, a candidate obscured active galactic nucleus (AGN) thought to be powered by an intermediate-mass black hole (IMBH, $M_\bullet \approx 10^{4-6} M_\odot$) of mass $M_{\bullet} \sim 10^{5.2}M_\odot$. To study its X-ray properties, we targeted J$144013+024744$ with NuSTAR for $\approx 100$ ks. The X-ray spectrum was… ▽ More We present a multi-wavelength analysis of the dwarf Seyfert-2 galaxy J$144013+024744$, a candidate obscured active galactic nucleus (AGN) thought to be powered by an intermediate-mass black hole (IMBH, $M_\bullet \approx 10^{4-6} M_\odot$) of mass $M_{\bullet} \sim 10^{5.2}M_\odot$. To study its X-ray properties, we targeted J$144013+024744$ with NuSTAR for $\approx 100$ ks. The X-ray spectrum was fitted with absorbed power law, Pexmon and a physical model (RXTorus). A Bayesian X-ray analysis was performed to estimate the posteriors. The phenomenological and the physical models suggest the AGN to be heavily obscured by a column density of $N_{\rm H} = (3.4-7.0)\times10^{23}$ cm$^{-2}$. In particular, the RXTorus model with a sub-solar metallicity suggests the obscuring column to be almost Compton-thick. We compared the $2-10$ keV intrinsic X-ray luminosity with the inferred X-ray luminosities based on empirical scaling relations for unobscured AGNs using $L_{\rm [OIV](25.89μ{\rm m})}$, $L_{[{\rm OIII}](5007 {\rm angstrom})}$, and $L_{6\rm μm}$ and found that the high-excitation $[{\rm OIV}]$ line provides a better estimate of the intrinsic $2-10$ keV X-ray luminosity ($L_{2-10}^{\rm int} \sim 10^{41.41}{\rm erg s}^{-1}$). Our results suggest that J$144013+024744$ is the first type-2 dwarf galaxy that shows X-ray spectroscopic evidence for obscuration. The column density that we estimated is among the highest measured to date for IMBH-powered AGNs, implying that a typical AGN torus geometry might extend to the low-mass end. This work has implications for constraining the black hole occupation fraction in dwarf galaxies using X-ray observations. △ Less

Submitted 20 September, 2022; originally announced September 2022.

Comments: Accepted for Publication in ApJ

Spectropolarimetric measurements of hidden broad lines in nearby megamaser galaxies: a lack of clear evidence for a correlation between black hole masses and virial products

Authors: Nora B. Linzer, Andy D. Goulding, Jenny E. Greene, Ryan C. Hickox

Abstract: High-accuracy black hole (BH) masses require excellent spatial resolution that is only achievable for galaxies within ~100 Mpc using present-day technology. At larger distances, BH masses are often estimated with single-epoch scaling relations for active galactic nuclei. This method requires only luminosity and the velocity dispersion of the broad line region (BLR) to calculate a virial product, a… ▽ More High-accuracy black hole (BH) masses require excellent spatial resolution that is only achievable for galaxies within ~100 Mpc using present-day technology. At larger distances, BH masses are often estimated with single-epoch scaling relations for active galactic nuclei. This method requires only luminosity and the velocity dispersion of the broad line region (BLR) to calculate a virial product, and an additional virial factor, $f$, to determine BH mass. The accuracy of these single-epoch masses, however, is unknown, and there are few empirical constraints on the variance of $f$ between objects. We attempt to calibrate single-epoch BH masses using spectropolarimetric measurements of nine megamaser galaxies from which we measure the velocity distribution of the BLR. We do not find strong evidence for a correlation between the virial products used for single-epoch masses and dynamical mass, neither for the megamaser sample alone or when combined with dynamical masses from reverberation mapping modeling. Furthermore, we find evidence that the virial parameter $f$ varies between objects, but we do not find strong evidence for a correlation with other observable parameters such as luminosity or broad line width. Although we cannot definitively rule out the existence of any correlation between dynamical mass and virial product, we find tension between allowed $f$ values for masers and those widely used in the literature. We conclude that the single-epoch method requires further investigation if it is to be used successfully to infer BH masses. △ Less

Submitted 18 October, 2022; v1 submitted 2 August, 2022; originally announced August 2022.

Comments: 27 pages, 17 figures, accepted by ApJ

A Probabilistic Method of Background Removal for High Energy Astrophysics Data

Authors: Steven Ehlert, Chien-Ting Chen, Doug A. Swartz, Ryan C. Hickox, Alexander A. Lutovinov, Andrey N. Semena, Roman Krivonos, Andrey E. Shtykovsky, Alexey Tkachenko

Abstract: We present a new statistical method for constructing background subtracted measurements from event list data gathered by X-ray and gamma ray observatories. This method was initially developed specifically to construct images that account for the high background fraction and low overall count rates observed in survey data from the Mikhail Pavlinsky ART-XC telescope aboard the Spektrum Röntgen Gamma… ▽ More We present a new statistical method for constructing background subtracted measurements from event list data gathered by X-ray and gamma ray observatories. This method was initially developed specifically to construct images that account for the high background fraction and low overall count rates observed in survey data from the Mikhail Pavlinsky ART-XC telescope aboard the Spektrum Röntgen Gamma (SRG) mission, although the mathematical underpinnings are valid for data taken with other imaging missions and analysis applications. This method fully accounts for the expected Poisson fluctuations in both the sky photon and non X-ray background count rates in a manner that does not result in unphysical negative counts. We derive the formulae for arbitrary confidence intervals for the source counts and show that our new measurement converges exactly to the standard background subtraction calculation in the high signal limit. Utilizing these results, we discuss several variants of images designed to optimize different science goals for both pointed and slewing telescopes. Using realistic simulated data of a galaxy cluster as observed by ART-XC we show that our method provides a more significant and robust detection of the cluster emission as compared to a standard background subtraction. We also demonstrate its advantages using real observations of a point source from the ART-XC telescope. These calculations may have widespread applications for a number of source classes observed with high energy telescopes. △ Less

Submitted 20 July, 2022; originally announced July 2022.

Comments: 15 pages, 9 figures, Accepted for Publication in MNRAS

Are Active Galactic Nuclei in Post-Starburst Galaxies Driving the Change or Along for the Ride?

Authors: Lauranne Lanz, Sofia Stepanoff, Ryan C. Hickox, Katherine Alatalo, K. Decker French, Kate Rowlands, Kristina Nyland, Phil Appleton, Mark Lacy, Anne Medling, John S. Mulchaey, Elizaveta Sazonova, Claudia Megan Urry

Abstract: We present an analysis of 10 ks snapshot Chandra observations of 12 shocked post-starburst galaxies, which provide a window into the unresolved question of active galactic nuclei (AGN) activity in post-starburst galaxies and its role in the transition of galaxies from actively star forming to quiescence. While 7/12 galaxies have statistically significant detections (with 2 more marginal detections… ▽ More We present an analysis of 10 ks snapshot Chandra observations of 12 shocked post-starburst galaxies, which provide a window into the unresolved question of active galactic nuclei (AGN) activity in post-starburst galaxies and its role in the transition of galaxies from actively star forming to quiescence. While 7/12 galaxies have statistically significant detections (with 2 more marginal detections), the brightest only obtained 10 photons. Given the wide variety of hardness ratios in this sample, we chose to pursue a forward modeling approach to constrain the intrinsic luminosity and obscuration of these galaxies rather than stacking. We constrain intrinsic luminosity of obscured power-laws based on the total number of counts and spectral shape, itself mostly set by the obscuration, with hardness ratios consistent with the data. We also tested thermal models. While all the galaxies have power-law models consistent with their observations, a third of the galaxies are better fit as an obscured power-law and another third are better fit as thermal emission. If these post-starburst galaxies, early in their transition, contain AGN, then these are mostly confined to a lower obscuration ($n_H \leq10^{23}$ cm$^{-2}$) and lower luminosity ($L_{2-10~ \rm keV}\leq10^{42}$erg s$^{-1}$). Two galaxies, however, are clearly best fit as significantly obscured AGN. At least half of this sample show evidence of at least low luminosity AGN activity, though none could radiatively drive out the remaining molecular gas reservoirs. Therefore, these AGN are more likely along for the ride, having been fed gas by the same processes driving the transition. △ Less

Submitted 6 July, 2022; v1 submitted 1 July, 2022; originally announced July 2022.

Comments: 16 pages, 6 figures, 3 tables, accepted for publication in ApJ

Monitoring Observations of SMC X-1's Excursions (MOOSE) I: Program Description and Initial High-State Spectral Results

Authors: Kristen C. Dage, McKinley Brumback, Joey Neilsen, Chin-Ping Hu, Diego Altamirano, Arash Bahramian, Philip A. Charles, William I. Clarkson, Daryl Haggard, Ryan C. Hickox, Jamie Kennea

Abstract: SMC X-1 has exhibited three super-orbital period excursions since the onset of X-ray monitoring beginning with RXTE's launch in 1995. NICER has recently probed a fourth observed excursion beginning in 2021 with our program Monitoring Observations of SMC X-1's Excursions (MOOSE). These sensitive new MOOSE data probe different super-orbital periods and phases within them. Spectral fits to the high-s… ▽ More SMC X-1 has exhibited three super-orbital period excursions since the onset of X-ray monitoring beginning with RXTE's launch in 1995. NICER has recently probed a fourth observed excursion beginning in 2021 with our program Monitoring Observations of SMC X-1's Excursions (MOOSE). These sensitive new MOOSE data probe different super-orbital periods and phases within them. Spectral fits to the high-state continuum during April 2021 to January 2022 show that the intrinsic spectral shapes are characterised by a soft (kT~0.19 keV) disc component and a hard (Gamma~0.7) power-law tail. When the 2021-2022 NICER observations, taken during an excursion, are compared to 2016 XMM-Newton observations (outside of an excursion), we find little evidence for intrinsic spectral variability across the high-states, but find evidence for a >3 sigma change in the absorption, although we caution that there may be calibration differences between the two instruments. Thus, over different lengths of super-orbital periods, we see little evidence for intrinsic spectral changes in the high-state. Upcoming studies of the pulse profiles may shed light on the mechanism behind the excursions. △ Less

Submitted 13 June, 2022; originally announced June 2022.

Comments: accepted to MNRAS

Understanding the Nature of an Unusual Post-Starburst Quasar with Exceptionally Strong Ne V Emission

Authors: Cameren Swiggum, Christy Tremonti, Serena Perrotta, Adam Schaefer, Ryan Hickox, Alison Coil, Paul Sell, Aleksandar Diamond-Stanic, Jalyn Krause, Gregory Mosby

Abstract: We present a $z = 0.94$ quasar, SDSS J004846.45-004611.9, discovered in the SDSS-III BOSS survey. A visual analysis of this spectrum reveals highly broadened and blueshifted narrow emission lines, in particular [Ne~V]$\lambda3426$ and [O~III]$\lambda5007$, with outflow velocities of 4000 km s$^{-1}$, along with unusually large [Ne V]$\lambda3426$/[Ne III]$\lambda3869$ ratios. The gas shows higher… ▽ More We present a $z = 0.94$ quasar, SDSS J004846.45-004611.9, discovered in the SDSS-III BOSS survey. A visual analysis of this spectrum reveals highly broadened and blueshifted narrow emission lines, in particular [Ne~V]$\lambda3426$ and [O~III]$\lambda5007$, with outflow velocities of 4000 km s$^{-1}$, along with unusually large [Ne V]$\lambda3426$/[Ne III]$\lambda3869$ ratios. The gas shows higher ionization at higher outflow velocities, indicating a connection between the powerful outflow and the unusual strength of the high ionization lines. The SED and the $i - \text{W3}$ color of the source reveal that it is likely a "core" Extremely Red Quasar (core ERQ); a candidate population of young AGN that are violently "blowing out" gas and dust from their centers. The dominance of host galaxy light in its spectrum and its fortuitous position in the SDSS S82 region allows us to measure its star formation history and investigate for variability for the first time in an ERQ. Our analysis indicates that SDSS J004846.45-004611.9 underwent a short-lived starburst phase 400 Myr ago and was subsequently quenched, possibly indicating a time-lag between star formation quenching and the onset of AGN activity. We also find that the strong extinction can be uniquely attributed to the AGN and does not persist in the host galaxy, contradicting a scenario where the source has recently transitioned from being a dusty sub-mm galaxy. In our relatively shallow photometric data, the source does not appear to be variable at $0.24-2.4~μ$m in the restframe, most likely due to the dominant contribution of host galaxy starlight at these wavelengths. △ Less

Submitted 10 March, 2022; originally announced March 2022.

Comments: 23 pages, 8 figures, accepted for publication in ApJ

Host Dark Matter Halos of SDSS Red and Blue Quasars: No Significant Difference in Large-scale Environment

Authors: Grayson C. Petter, Ryan C. Hickox, David M. Alexander, James E. Geach, Adam D. Myers, David J. Rosario, Victoria A. Fawcett, Lizelke Klindt, Kelly E. Whalen

Abstract: The observed optical colors of quasars are generally interpreted in one of two frameworks: unified models which attribute color to random orientation of the accretion disk along the line-of-sight, and evolutionary models which invoke connections between quasar systems and their environments. We test these schema by probing the dark matter halo environments of optically-selected quasars as a functi… ▽ More The observed optical colors of quasars are generally interpreted in one of two frameworks: unified models which attribute color to random orientation of the accretion disk along the line-of-sight, and evolutionary models which invoke connections between quasar systems and their environments. We test these schema by probing the dark matter halo environments of optically-selected quasars as a function of $g-i$ optical color by measuring the two-point correlation functions of $\sim$ 0.34 million eBOSS quasars as well as the gravitational deflection of cosmic microwave background photons around $\sim$ 0.66 million XDQSO photometric quasar candidates. We do not detect a trend of halo bias with optical color through either analysis, finding that optically-selected quasars at $0.8 < z < 2.2$ occupy halos of characteristic mass $M_{h}\sim 3\times 10^{12} \ h^{-1} M_{\odot}$ regardless of their color. This result implies that a quasar's large-scale halo environment is not strongly connected to its observed optical color. We also confirm findings of fundamental differences in the radio properties of red and blue quasars by stacking 1.4 GHz FIRST images at their positions, suggesting the observed differences cannot be attributed to orientation. Instead, the differences between red and blue quasars likely arise on nuclear-galactic scales, perhaps owing to reddening by a nuclear dusty wind. Finally, we show that optically-selected quasars' halo environments are also independent of their $r-W2$ optical-infrared colors, while previous work has suggested that mid-infrared-selected obscured quasars occupy more massive halos. We discuss implications of this result for models of quasar and galaxy co-evolution. △ Less

Submitted 19 January, 2022; originally announced January 2022.

Comments: Accepted for publication in The Astrophysical Journal

BASS XXX: Distribution Functions of DR2 Eddington-ratios, Black Hole Masses, and X-ray Luminosities

Authors: Tonima Tasnim Ananna, Anna K. Weigel, Benny Trakhtenbrot, Michael J. Koss, C. Megan Urry, Claudio Ricci, Ryan C. Hickox, Ezequiel Treister, Franz E. Bauer, Yoshihiro Ueda, Richard Mushotzky, Federica Ricci, Kyuseok Oh, Julian E. Mejia-Restrepo, Jakob Den Brok, Daniel Stern, Meredith C. Powell, Turgay Caglar, Kohei Ichikawa, O. Ivy Wong, Fiona A. Harrison, Kevin Schawinski

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

Submitted 6 February, 2022; v1 submitted 14 January, 2022; originally announced January 2022.

Comments: Accepted by APJS

The NuSTAR Extragalactic Survey of the James Webb Space Telescope North Ecliptic Pole Time-Domain Field

Authors: X. Zhao, F. Civano, F. M. Fornasini, D. M. Alexander, N. Cappelluti, C. T. Chen, S. H. Cohen, M. Elvis, P. Gandhi, N. A. Grogin, R. C. Hickox, R. A. Jansen, A. Koekemoer, G. Lanzuisi, W. P. Maksym, A. Masini, D. J. Rosario, M. J. Ward, C. N. A. Willmer, R. A. Windhorst

Abstract: We present the $NuSTAR$ extragalactic survey of the $James$ $Webb$ $Space$ $Telescope$ ($JWST$) North Ecliptic Pole (NEP) Time-Domain Field. The survey covers a $\sim$0.16 deg$^2$ area with a total exposure of 681 ks acquired in a total of nine observations from three epochs. The survey sensitivities at 20% of the area are 2.39, 1.14, 2.76, 1.52, and 5.20 $\times$ 10$^{-14}$ erg cm$^{-2}$ s… ▽ More We present the $NuSTAR$ extragalactic survey of the $James$ $Webb$ $Space$ $Telescope$ ($JWST$) North Ecliptic Pole (NEP) Time-Domain Field. The survey covers a $\sim$0.16 deg$^2$ area with a total exposure of 681 ks acquired in a total of nine observations from three epochs. The survey sensitivities at 20% of the area are 2.39, 1.14, 2.76, 1.52, and 5.20 $\times$ 10$^{-14}$ erg cm$^{-2}$ s$^{-1}$ in the 3-24, 3-8, 8-24, 8-16, and 16-24 keV bands, respectively. The NEP survey is one of the most sensitive extragalactic surveys with $NuSTAR$ so far. A total of 33 sources were detected above 95% reliability in at least one of the five bands. We present the number counts, log$N$-log$S$, measured in the hard X-ray 8-24 and 8-16 keV bands, uniquely accessible by $NuSTAR$ down to such faint fluxes. We performed source detection on the XMM-$Newton$ and $Chandra$ observations of the same field to search for soft X-ray counterparts of each $NuSTAR$ detection. The soft band positions were used to identify optical and infrared associations. We present the X-ray properties (hardness ratio and luminosity) and optical-to-X-ray properties of the detected sources. The measured fraction of candidate Compton-thick (N$\rm _H\ge10^{24} cm^{-2}$) active galactic nuclei, derived from the hardness ratio, is between 3% to 27%. As this survey was designed to have variability as its primary focus, we present preliminary results on multi-epoch flux variability in the 3-24 keV band. △ Less

Submitted 29 December, 2021; v1 submitted 28 September, 2021; originally announced September 2021.

Comments: 22 pages, 17 figures

Physical Properties of Massive Compact Starburst Galaxies with Extreme Outflows

Authors: Serena Perrotta, Erin R. George, Alison L. Coil, Christy A. Tremonti, David S. N. Rupke, Julie D. Davis, Aleksandar M. Diamond-Stanic, James E. Geach, Ryan C. Hickox, John Moustakas, Grayson C. Petter, Gregory H. Rudnick, Paul H. Sell, Cameren Swiggum, Kelly E. Whalen

Abstract: We present results on the nature of extreme ejective feedback episodes and the physical conditions of a population of massive ($\rm M_* \sim 10^{11} M_{\odot}$), compact starburst galaxies at z = 0.4-0.7. We use data from Keck/NIRSPEC, SDSS, Gemini/GMOS, MMT, and Magellan/MagE to measure rest-frame optical and near-IR spectra of 14 starburst galaxies with extremely high star formation rate surface… ▽ More We present results on the nature of extreme ejective feedback episodes and the physical conditions of a population of massive ($\rm M_* \sim 10^{11} M_{\odot}$), compact starburst galaxies at z = 0.4-0.7. We use data from Keck/NIRSPEC, SDSS, Gemini/GMOS, MMT, and Magellan/MagE to measure rest-frame optical and near-IR spectra of 14 starburst galaxies with extremely high star formation rate surface densities (mean $\rm Σ_{SFR} \sim 3000 \,M_{\odot} yr^{-1} kpc^{-2}$) and powerful galactic outflows (maximum speeds v$_{98} \sim$ 1000-3000 km s$^{-1}$). Our unique data set includes an ensemble of both emission [OII]$λλ$3726,3729, H$β$, [OIII]$λλ$4959,5007, H$α$, [NII]$λλ$6548,6583, and [SII]$λλ$6716,6731) and absorption MgII$λλ$2796,2803, and FeII$λ$2586) lines that allow us to investigate the kinematics of the cool gas phase (T$\sim$10$^4$ K) in the outflows. Employing a suite of line ratio diagnostic diagrams, we find that the central starbursts are characterized by high electron densities (median n$_e \sim$ 530 cm$^{-3}$), and high metallicity (solar or super-solar). We show that the outflows are most likely driven by stellar feedback emerging from the extreme central starburst, rather than by an AGN. We also present multiple intriguing observational signatures suggesting that these galaxies may have substantial Lyman continuum (LyC) photon leakage, including weak [SII] nebular emission lines. Our results imply that these galaxies may be captured in a short-lived phase of extreme star formation and feedback where much of their gas is violently blown out by powerful outflows that open up channels for LyC photons to escape. △ Less

Submitted 18 October, 2021; v1 submitted 4 June, 2021; originally announced June 2021.

Where Do Obscured AGN Fit in A Galaxy's Timeline?

Authors: Cassandra Hatcher, Allison Kirkpatrick, Francesca Fornasini, Francesca Civano, Erini Lambrides, Dale Kocesvski, Christopher M. Carroll, Mauro Giavalisco, Ryan Hickox, Zhiyuan Ji

Abstract: 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 \texti… ▽ More 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. △ Less

Submitted 28 May, 2021; originally announced May 2021.

Chandra Observations of Excess Fe K$α$ Line Emission in Galaxies with High Star Formation Rates: X-ray Reflection on Galaxy Scales?

Authors: Wei Yan, Ryan C. Hickox, Chien-Ting J. Chen, Claudio Ricci, Alberto Masini, Franz E. Bauer, David M. Alexander

Abstract: In active galactic nuclei (AGN), fluorescent Fe K$α$ (iron) line emission is generally interpreted as originating from obscuring material around a supermassive black hole (SMBH) on the scale of a few parsecs (pc). However, recent Chandra studies indicate the existence of iron line emission extending to kpc scales in the host galaxy. The connection between iron line emission and large-scale materia… ▽ More In active galactic nuclei (AGN), fluorescent Fe K$α$ (iron) line emission is generally interpreted as originating from obscuring material around a supermassive black hole (SMBH) on the scale of a few parsecs (pc). However, recent Chandra studies indicate the existence of iron line emission extending to kpc scales in the host galaxy. The connection between iron line emission and large-scale material can be spatially resolved directly only in nearby galaxies, but could be inferred in more distant AGNs by a connection between line emission and star-forming gas and dust that is more extended than the pc-scale torus. Here we present the results from a stacking analysis and X-ray spectral fitting performed on sources in the Chandra Deep Field South (CDFS) 7 Ms observations. From the deep stacked spectra, we select sources with stellar mass $\log(M_*/M_\odot)>10$ at $0.5<z<2$, obtaining 25 sources with high infrared luminosity ($ {\rm SFR}_{\rm FIR} \geq 17\;M_{\odot}\;{\rm yr}^{-1}$) and 32 sources below this threshold. We find that the equivalent width of the iron line EW(Fe) is a factor of three higher with 3$σ$ significance for high infrared luminosity measured from Herschel observations, indicating a connection between iron line emission and star-forming material on galaxy scales. We show that there is no significant dependence in EW(Fe) on $M_*$ or X-ray luminosity, suggesting the reflection of AGN X-ray emission over large scales in their host galaxies may be widespread. △ Less

Submitted 21 April, 2021; originally announced April 2021.

Comments: 8 pages, 5 figures, accepted by ApJ

Compact Starburst Galaxies with Fast Outflows: Central Escape Velocities and Stellar Mass Surface Densities from Multi-band Hubble Space Telescope Imaging

Authors: Aleksandar M. Diamond-Stanic, John Moustakas, Paul H. Sell, Christy A. Tremonti, Alison L. Coil, Julie D. Davis, James E. Geach, Sophia C. W. Gottlieb, Ryan C. Hickox, Amanda Kepley, Charles Lipscomb, Joshua Rines, Gregory H. Rudnick, Cristopher Thompson, Kingdell Valdez, Christian Bradna, Jordan Camarillo, Eve Cinquino, Senyo Ohene Serena Perrotta, Grayson C. Petter, David S. N. Rupke, Chidubem Umeh, Kelly E. Whalen

Abstract: We present multi-band Hubble Space Telescope imaging that spans rest-frame near-ultraviolet through near-infrared wavelengths (0.3-1.1 $μ$m) for 12 compact starburst galaxies at z=0.4-0.8. These massive galaxies (M_stellar ~ 10^11 M_Sun) are driving very fast outflows ($v_{max}$=1000-3000 km/s), and their light profiles are dominated by an extremely compact starburst component (half-light radius ~… ▽ More We present multi-band Hubble Space Telescope imaging that spans rest-frame near-ultraviolet through near-infrared wavelengths (0.3-1.1 $μ$m) for 12 compact starburst galaxies at z=0.4-0.8. These massive galaxies (M_stellar ~ 10^11 M_Sun) are driving very fast outflows ($v_{max}$=1000-3000 km/s), and their light profiles are dominated by an extremely compact starburst component (half-light radius ~ 100 pc). Our goal is to constrain the physical mechanisms responsible for launching these fast outflows by measuring the physical conditions within the central kiloparsec. Based on our stellar population analysis, the central component typically contributes $\approx$25% of the total stellar mass and the central escape velocities $v_{esc,central}\approx900$ km/s are a factor of two smaller than the observed outflow velocities. This requires physical mechanisms that can accelerate gas to speeds significantly beyond the central escape velocities, and it makes clear that these fast outflows are capable of traveling into the circumgalactic medium, and potentially beyond. We find central stellar densities comparable to theoretical estimates of the Eddington limit, and we estimate $Σ_1$ surface densities within the central kpc comparable to those of compact massive galaxies at $0.5<z<3.0$. Relative to "red nuggets" and "blue nuggets" at $z\sim2$, we find significantly smaller $r_e$ values at a given stellar mass, which we attribute to the dominance of a young stellar component in our sample and the better physical resolution for rest-frame optical observations at $z\sim0.6$ versus $z\sim2$. We compare to theoretical scenarios involving major mergers and violent disc instability, and we speculate that our galaxies are progenitors of power-law ellipticals in the local universe with prominent stellar cusps. △ Less

Submitted 22 February, 2021; originally announced February 2021.

Comments: accepted for publication in ApJ

A broad-band X-ray view of the precessing accretion disk and pre-eclipse dip in the pulsar Her X-1 with NuSTAR and XMM-Newton

Authors: McKinley C. Brumback, Ryan C. Hickox, Felix S. Fürst, Katja Pottschmidt, John A. Tomsick, Jörn Wilms, Rüdiger Staubert, Saeqa Vrtilek

Abstract: We present a broad-band X-ray timing study of the variations in pulse behavior with superorbital cycle in the low-mass X-ray binary Her X-1. This source shows a 35-day superorbital modulation in X-ray flux that is likely caused by occultation by a warped, precessing accretion disk. Our data set consists of four joint XMM-Newton and NuSTAR observations of Her X-1 which sample a complete superorbita… ▽ More We present a broad-band X-ray timing study of the variations in pulse behavior with superorbital cycle in the low-mass X-ray binary Her X-1. This source shows a 35-day superorbital modulation in X-ray flux that is likely caused by occultation by a warped, precessing accretion disk. Our data set consists of four joint XMM-Newton and NuSTAR observations of Her X-1 which sample a complete superorbital cycle. We focus our analysis on the first and fourth observations, which occur during the bright "main-on" phase, because these observations have strongly detected pulsations. We added an archival XMM-Newton observation during the "short-on" phase of the superorbital cycle since our observations at that phase are lower in signal to noise. We find that the energy-resolved pulse profiles show the same shape at similar superorbital phases and the profiles are consistent with expectations from a precessing disk. We demonstrate that a simple precessing accretion disk model is sufficient to reproduce the observed pulse profiles. The results of this model suggest that the similarities in the observed pulse profiles are due to reprocessing by a precessing disk that has returned to its original precession phase. We determine that the broad-band spectrum is well fit by an absorbed power law with a soft blackbody component, and show that the spectral continuum also exhibits dependence on the superorbital cycle. We also present a brief analysis of the energy resolved light curves of a pre-eclipse dip, which shows soft X-ray absorption and hard X-ray variability during the dip. △ Less

Submitted 9 February, 2021; originally announced February 2021.

Comments: 23 pages, 13 figures. Accepted to ApJ

Observational Evidence for Enhanced Black Hole Accretion in Giant Elliptical Galaxies

Authors: Michael McDonald, Brian R. McNamara, Michael S. Calzadilla, Chien-Ting Chen, Massimo Gaspari, Ryan C. Hickox, Erin Kara, Ilia Korchagin

Abstract: We present a study of the relationship between black hole accretion rate (BHAR) and star formation rate (SFR) in a sample of giant elliptical galaxies. These galaxies, which live at the centers of galaxy groups and clusters, have star formation and black hole activity that is primarily fueled by gas condensing out of the hot intracluster medium. For a sample of 46 galaxies spanning 5 orders of mag… ▽ More We present a study of the relationship between black hole accretion rate (BHAR) and star formation rate (SFR) in a sample of giant elliptical galaxies. These galaxies, which live at the centers of galaxy groups and clusters, have star formation and black hole activity that is primarily fueled by gas condensing out of the hot intracluster medium. For a sample of 46 galaxies spanning 5 orders of magnitude in BHAR and SFR, we find a mean ratio of log(BHAR/SFR) = -1.45 +/- 0.2, independent of the methodology used to constrain both SFR and BHAR. This ratio is significantly higher than most previously-published values for field galaxies. We investigate whether these high BHAR/SFR ratios are driven by high BHAR, low SFR, or a different accretion efficiency in radio galaxies. The data suggest that the high BHAR/SFR ratios are primarily driven by boosted black hole accretion in spheroidal galaxies compared to their disk counterparts. We propose that angular momentum of the cool gas is the primary driver in suppressing BHAR in lower mass galaxies, with massive galaxies accreting gas that has condensed out of the hot phase on nearly radial trajectories. Additionally, we demonstrate that the relationship between specific BHAR and SFR has much less scatter over 6 orders of magnitude in both parameters, due to competing dependence on morphology between the M_BH--M_* and BHAR--SFR relations. In general, active galaxies selected by typical techniques have sBHAR/sSFR ~ 10, while galactic nuclei with no clear AGN signatures have sBHAR/sSFR ~ 1, consistent with a universal M_BH--M_spheroid relation. △ Less

Submitted 16 December, 2020; originally announced December 2020.

Comments: 15 pages, 12 figure. Accepted for publication in ApJ

A Large Population of Luminous Active Galactic Nuclei Lacking X-ray Detections: Evidence for Heavy Obscuration?

Authors: Christopher M. Carroll, Ryan C. Hickox, Alberto Masini, Lauranne Lanz, Roberto J. Assef, Daniel Stern, Chien-Ting J. Chen, Tonima T. Ananna

Abstract: We present a large sample of infrared-luminous candidate active galactic nuclei (AGNs) that lack X-ray detections in Chandra, XMM-Newton, and NuSTAR fields. We selected all optically detected SDSS sources with redshift measurements, combined additional broadband photometry from WISE, UKIDSS, 2MASS, and GALEX, and modeled the spectral energy distributions (SEDs) of our sample sources. We parameteri… ▽ More We present a large sample of infrared-luminous candidate active galactic nuclei (AGNs) that lack X-ray detections in Chandra, XMM-Newton, and NuSTAR fields. We selected all optically detected SDSS sources with redshift measurements, combined additional broadband photometry from WISE, UKIDSS, 2MASS, and GALEX, and modeled the spectral energy distributions (SEDs) of our sample sources. We parameterize nuclear obscuration in our SEDs with $E(B\!-\!V)_{\text{AGN}}$ and uncover thousands of powerful obscured AGNs that lack X-ray counterparts, many of which are identified as AGN candidates based on straightforward WISE photometric criteria. Using the observed luminosity correlation between restframe 2-10 keV ($L_{\text{X}}$) and restframe AGN 6 $μ{\text{m}}$ ($L_{\text{MIR}}$), we estimate the intrinsic X-ray luminosities of our sample sources and combine these data with flux limits from X-ray catalogs to determine lower limits on nuclear obscuration. Using the ratio of intrinsic-to-observed X-ray luminosity ($R_{L_{\text{X}}}$), we find a significant fraction of sources with column densities approaching $N_{\text{H}}>$ 10$^{\text{24}}$ cm$^{-{\text{2}}}$, suggesting that multiwavelength observations are necessary to account for the population of heavily obscured AGNs. We simulate the underlying $N_{\text{H}}$ distribution for the X-ray non-detected sources in our sample through survival analysis, and confirm the presence of AGN activity via X-ray stacking. Our results point to a considerable population of extremely obscured AGNs undetected by current X-ray observatories. △ Less

Submitted 8 December, 2020; originally announced December 2020.

Comments: 15 pages, 6 tables, 10 figures. Accepted to ApJ. For the associated catalogs (tables 3 and 4), see http://dartgo.org/ccarroll_xagn_cat

The Dust-to-Gas Ratio and the Role of Radiation Pressure in Luminous, Obscured Quasars

Authors: Hyunsung D. Jun, Roberto J. Assef, Christopher M. Carroll, Ryan C. Hickox, Yonghwi Kim, Jaehyun Lee, Claudio Ricci, Daniel Stern

Abstract: The absence of high Eddington ratio, obscured Active Galactic Nuclei (AGN) in local ($z\lesssim0.1$) samples of moderate luminosity AGN has generally been explained to result from radiation pressure on the dusty gas governing the level of nuclear ($\lesssim10$pc) obscuration. However, very high accretion rates are routinely reported among obscured quasars at higher luminosities, and may require a… ▽ More The absence of high Eddington ratio, obscured Active Galactic Nuclei (AGN) in local ($z\lesssim0.1$) samples of moderate luminosity AGN has generally been explained to result from radiation pressure on the dusty gas governing the level of nuclear ($\lesssim10$pc) obscuration. However, very high accretion rates are routinely reported among obscured quasars at higher luminosities, and may require a different feedback mechanism. We compile constraints on obscuration and Eddington ratio for samples of X-ray, optical, infrared, and submm selected AGN at quasar luminosities. Whereas moderate luminosity, obscured AGN in the local universe have a range of lower Eddington ratios ($f_{\rm Edd} \sim 0.001-0.1$), the most luminous ($L_{\rm bol} \gtrsim 10^{46} $erg s$^{-1}$) IR/submm-bright, obscured quasars out to $z\sim3$ commonly have very high Eddington ratios ($f_{\rm Edd} \sim 0.1-1$). This apparent lack of radiation pressure feedback in luminous obscured quasars is likely coupled with AGN timescales, such that a higher fraction of luminous obscured quasars are seen due to the short timescale for which quasars are most luminous. Adopting quasar evolutionary scenarios, extended ($\sim10^{2-3}$pc) obscuration may work together with the shorter timescales to explain the observed fraction of obscured, luminous quasars, while outflows driven by radiation pressure will slowly clear this material over the AGN lifetime. △ Less

Submitted 3 January, 2021; v1 submitted 29 October, 2020; originally announced October 2020.

Comments: 15 pages, 6 figures, accepted to ApJ, typos corrected

The clustering of submillimeter galaxies detected with ALMA

Authors: Cristina Garcia-Vergara, Jacqueline Hodge, Joseph F. Hennawi, Axel Weiss, Julie Wardlow, Adam D. Myers, Ryan Hickox

Abstract: Previous studies measuring the clustering of submillimeter galaxies (SMGs) have based their measurements on single-dish detected sources, finding evidence for strong clustering. However, ALMA has revealed that, due to the coarse angular resolution of these instruments, single-dish sources can be comprised of multiple sources. This implies that the clustering inferred from single-dish surveys may b… ▽ More Previous studies measuring the clustering of submillimeter galaxies (SMGs) have based their measurements on single-dish detected sources, finding evidence for strong clustering. However, ALMA has revealed that, due to the coarse angular resolution of these instruments, single-dish sources can be comprised of multiple sources. This implies that the clustering inferred from single-dish surveys may be overestimated. Here, we measure the clustering of SMGs based on the ALESS survey, an ALMA follow-up of sources previously identified in the LABOCA ECDFS Submillimeter Survey (LESS). We present a method to measure the clustering of ALMA sources that have been previously identified using single-dish telescopes, based on forward modeling both the single-dish and the ALMA observations. We constrain upper limits for the median mass of halos hosting SMGs at $1<z<3$, finding $M_{\rm{halo}}\leq2.4\times10^{12}\,\rm\,M_{\odot}$ for SMGs with flux densities $S_{870}\geq4.0\,$mJy, which is at least $3.8^{+3.8}_{-2.6}$ times lower than the mass inferred based on the clustering of the LESS sources alone. This suggests that the strength of SMG clustering based on single-dish observations was overestimated and therefore SMGs might be hosted by dark matter halos less massive than has previously been estimated. By extrapolating our models down to flux densities of $S_{870}\geq1.2\,\rm\,$mJy, we find that such SMGs inhabit halos with median mass $M_{\rm{halo}}\leq3.2\times10^{11}\,\rm\,M_{\odot}$. We conclude that only the brightest ($S_{870}\gtrsim5-6\,\rm\,$mJy) SMGs would trace massive structures at $z\sim2$ and only SMGs with $S_{870}\gtrsim6\,\rm\,$mJy may be connected to massive local elliptical galaxies, quasars at intermediate redshifts and high-redshift star-forming galaxies, whereas fainter SMGs are unlikely linked to these populations. △ Less

Submitted 2 October, 2020; originally announced October 2020.

Comments: 20 pages, 14 figures, Accepted for publication in ApJ

Accretion History of AGN III: Radiative Efficiency and AGN Contribution to Reionization

Authors: Tonima Tasnim Ananna, C. Megan Urry, Ezequiel Treister, Ryan C. Hickox, Francesco Shankar, Claudio Ricci, Nico Cappelluti, Stefano Marchesi, Tracey Jane Turner

Abstract: The cosmic history of supermassive black hole (SMBH) growth is important for understanding galaxy evolution, reionization and the physics of accretion. Recent NuSTAR, Swift-BAT and \textit{Chandra} hard X-ray surveys have provided new constraints on the space density of heavily obscured Active Galactic Nuclei (AGN). Using the new X-ray luminosity function derived from these data, we here estimate… ▽ More The cosmic history of supermassive black hole (SMBH) growth is important for understanding galaxy evolution, reionization and the physics of accretion. Recent NuSTAR, Swift-BAT and \textit{Chandra} hard X-ray surveys have provided new constraints on the space density of heavily obscured Active Galactic Nuclei (AGN). Using the new X-ray luminosity function derived from these data, we here estimate the accretion efficiency of SMBHs and their contribution to reionization. We calculate the total ionizing radiation from active galactic nuclei (AGN) as a function of redshift, based on the X radiation and distribution of obscuring column density, converted to UV wavelengths. Limiting the luminosity function to unobscured AGN only, our results agree with current UV luminosity functions of unobscured AGN. For realistic assumptions about the escape fraction, the contribution of all AGN to cosmic reionization is $\sim4$ times lower than the galaxy contribution (23% at $z\sim6$). Our results also offer an observationally constrained prescription that can be used in simulations or models of galaxy evolution. To estimate the average efficiency with which supermassive black holes convert mass to light, we compare the total radiated energy, converted from X-ray light using a bolometric correction, to the most recent local black hole mass density. The most likely value, $η\sim 0.3-0.34$, approaches the theoretical limit for a maximally rotating Kerr black hole, $η=0.42$, implying that on average growing supermassive black holes are spinning rapidly. △ Less

Submitted 16 September, 2020; originally announced September 2020.

Comments: Accepted by ApJ, 13 pages, 5 figures

Deviations from the Infrared-Radio Correlation in Massive, Ultra-compact Starburst Galaxies

Authors: Grayson C. Petter, Amanda A. Kepley, Ryan C. Hickox, Gregory H. Rudnick, Christy A. Tremonti, Aleksandar M. Diamond-Stanic, James E. Geach, Alison L. Coil, Paul H. Sell, John Moustakas, David S. N. Rupke, Serena Perrotta, Kelly E. Whalen, Julie D. Davis

Abstract: Feedback through energetic outflows has emerged as a key physical process responsible for transforming star-forming galaxies into the quiescent systems observed in the local universe. To explore this process, this paper focuses on a sample of massive and compact merger remnant galaxies hosting high-velocity gaseous outflows ($|v| \gtrsim 10^{3}$ km s$^{-1}$), found at intermediate redshift (… ▽ More Feedback through energetic outflows has emerged as a key physical process responsible for transforming star-forming galaxies into the quiescent systems observed in the local universe. To explore this process, this paper focuses on a sample of massive and compact merger remnant galaxies hosting high-velocity gaseous outflows ($|v| \gtrsim 10^{3}$ km s$^{-1}$), found at intermediate redshift ($z \sim 0.6$). From their mid-infrared emission and compact morphologies, these galaxies are estimated to have exceptionally large star formation rate (SFR) surface densities ($Σ_{SFR} \sim 10^{3}$ $\mathrm{M_{\odot}}$ yr$^{-1}$ kpc$^{-2}$), approaching the Eddington limit for radiation pressure on dust grains. This suggests that star formation feedback may be driving the observed outflows. However, these SFR estimates suffer from significant uncertainties. We therefore sought an independent tracer of star formation to probe the compact starburst activity in these systems. In this paper, we present SFR estimates calculated using 1.5 GHz continuum Jansky Very Large Array observations for 19 of these galaxies. We also present updated infrared (IR) SFRs calculated from WISE survey data. We estimate SFRs from the IR to be larger than those from the radio for 16 out of 19 galaxies by a median factor of 2.5. We find that this deviation is maximized for the most compact galaxies hosting the youngest stellar populations, suggesting that compact starbursts deviate from the IR-radio correlation. We suggest that this deviation stems either from free-free absorption of synchrotron emission, a difference in the timescale over which each indicator traces star formation, or exceptionally hot IR-emitting dust in these ultra-dense galaxies. △ Less

Submitted 8 September, 2020; originally announced September 2020.

Comments: 15 pages, 9 figures, 2 tables; accepted for publication in ApJ

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

The Bolometric Quasar Luminosity Function at z = 0-7

Authors: Xuejian Shen, Philip F. Hopkins, Claude-André Faucher-Giguère, D. M. Alexander, Gordon T. Richards, Nicholas P. Ross, R. C. Hickox

Abstract: In this paper, we provide updated constraints on the bolometric quasar luminosity function (QLF) from $z=0$ to $z=7$. The constraints are based on an observational compilation that includes observations in the rest-frame IR, B band, UV, soft and hard X-ray in past decades. Our method follows Hopkins et al. 2007 with an updated quasar SED model and bolometric and extinction corrections. The new bes… ▽ More In this paper, we provide updated constraints on the bolometric quasar luminosity function (QLF) from $z=0$ to $z=7$. The constraints are based on an observational compilation that includes observations in the rest-frame IR, B band, UV, soft and hard X-ray in past decades. Our method follows Hopkins et al. 2007 with an updated quasar SED model and bolometric and extinction corrections. The new best-fit bolometric quasar luminosity function behaves qualitatively different from the Hopkins et al. 2007 model at high redshift. Compared with the old model, the number density normalization decreases towards higher redshift and the bright-end slope is steeper at $z\gtrsim 2$. Due to the paucity of measurements at the faint end, the faint end slope at $z\gtrsim 5$ is quite uncertain. We present two models, one featuring a progressively steeper faint-end slope at higher redshift and the other featuring a shallow faint-end slope at $z\gtrsim 5$. Further multi-band observations of the faint-end QLF are needed to distinguish between these models. The evolutionary pattern of the bolometric QLF can be interpreted as an early phase likely dominated by the hierarchical assembly of structures and a late phase likely dominated by the quenching of galaxies. We explore the implications of this model on the ionizing photon production by quasars, the CXB spectrum, the SMBH mass density and mass functions. The predicted hydrogen photoionization rate contributed by quasars is subdominant during the epoch of reionization and only becomes important at $z\lesssim 3$. The predicted CXB spectrum, cosmic SMBH mass density and SMBH mass function are generally consistent with existing observations. △ Less

Submitted 15 May, 2020; v1 submitted 8 January, 2020; originally announced January 2020.

Comments: Accepted by MNRAS. The code and data are available at https://bitbucket.org/ShenXuejian/quasarlf/src/master/

Physical Models for the Clustering of Obscured and Unobscured Quasars

Authors: Kelly E. Whalen, Ryan C. Hickox, Michael A. DiPompeo, Gordon T. Richards, Adam D. Myers

Abstract: Clustering measurements of obscured and unobscured quasars show that obscured quasars reside in more massive dark matter halos than their unobscured counterparts. These results are inconsistent with simple unified (torus) scenarios, but might be explained by models in which the distribution of obscuring material depends on Eddington ratio or galaxy stellar mass. We test these possibilities by cons… ▽ More Clustering measurements of obscured and unobscured quasars show that obscured quasars reside in more massive dark matter halos than their unobscured counterparts. These results are inconsistent with simple unified (torus) scenarios, but might be explained by models in which the distribution of obscuring material depends on Eddington ratio or galaxy stellar mass. We test these possibilities by constructing simple physical models to compare to observed AGN populations. We find that previously observed relationships between obscuration and Eddington ratio or stellar mass are not sufficient reproduce the observed quasar clustering results ($\langle \log M_{\text{halo}}/M_{\odot} \rangle = 12.94 ^{+ 0.10}_{- 0.11}$ and $\langle \log M_{\text{halo}}/M_{\odot} \rangle = 12.49 ^{+ 0.08}_{- 0.08}$ for obscured and unobscured populations, respectively) while maintaining the observed fraction of obscured quasars (30-65$\%$). This work suggests that evolutionary models, in which obscuration evolves on the typical timescale for black hole growth, are necessary to understand the observed clustering of mid-IR selected quasars. △ Less

Submitted 29 September, 2022; v1 submitted 3 December, 2019; originally announced December 2019.

Comments: 14 pages, 10 figures, accepted for publication in ApJ

Accretion History of AGN II: Constraints on AGN Spectral Parameters using the Cosmic X-ray Background

Authors: Tonima Tasnim Ananna, Ezequiel Treister, Claudia M. Urry, Claudio Ricci, Ryan Hickox, Nikhil Padmanabhan, Stefano Marchesi, Allison Kirkpatrick

Abstract: We constrain X-ray spectral shapes for the ensemble of AGN based on the shape of the Cosmic X-ray Background (CXB). Specifically, we rule out regions of X-ray spectral parameter space that do not reproduce the CXB in the energy range 1-100 keV. The key X-ray spectral parameters are the photon index, Γ; the cutoff energy, Ecutoff; and the reflection scaling factor, R. Assuming each parameter follow… ▽ More We constrain X-ray spectral shapes for the ensemble of AGN based on the shape of the Cosmic X-ray Background (CXB). Specifically, we rule out regions of X-ray spectral parameter space that do not reproduce the CXB in the energy range 1-100 keV. The key X-ray spectral parameters are the photon index, Γ; the cutoff energy, Ecutoff; and the reflection scaling factor, R. Assuming each parameter follows a Gaussian distribution, we first explore the parameter space using a Bayesian approach and a fixed X-ray luminosity function (XLF). For σ_E = 36 keV and σ_R = 0.14, fixed at the observed values from the Swift-BAT 70-month sample, we allow <R>, <Ecutoff > and <Γ> to vary subject to reproducing the CXB. We report results for σ_Γ = 0.1-0.5. In an alternative approach, we define the parameter distributions, then forward model to fit the CXB by perturbing the XLF using a neural network. This approach allows us to rule out parameter combinations that cannot reproduce the CXB for any XLF. The marginalized conditional probabilities for the four free parameters are: <R> = 0.99^{+0.11}_{-0.26}, <Ecutoff> = 118^{+24}_{-23}, σ_Γ = 0.101^{+0.097}_{-0.001} and <Γ> = 1.9^{+0.08}_{-0.09}. We provide an interactive online tool for users to explore any combination of <Ecutoff>, σ_E, <Γ>, σ_Γ, <R> and σ_R including different distributions for each absorption bin, subject to the integral CXB constraint. The distributions observed in many AGN samples can be ruled out by our analysis, meaning these samples can not be representative of the full AGN population. The few samples that fall within the acceptable parameter space are hard X-ray-selected, commensurate with their having fewer selection biases. △ Less

Submitted 25 November, 2019; originally announced November 2019.

Comments: Accepted by ApJ, 15 pages, 7 figures and 3 tables

A 100-kiloparsec wind feeding the circumgalactic medium of a massive compact galaxy

Authors: David S. N. Rupke, Alison Coil, James E. Geach, Christy Tremonti, Aleksandar M. Diamond-Stanic, Erin R. George, Ryan C. Hickox, Amanda A. Kepley, Gene Leung, John Moustakas, Gregory Rudnick, Paul H. Sell

Abstract: Ninety per cent of baryons are located outside galaxies, either in the circumgalactic or intergalactic medium. Theory points to galactic winds as the primary source of the enriched and massive circumgalactic medium. Winds from compact starbursts have been observed to flow to distances somewhat greater than ten kiloparsecs, but the circumgalactic medium typically extends beyond a hundred kiloparsec… ▽ More Ninety per cent of baryons are located outside galaxies, either in the circumgalactic or intergalactic medium. Theory points to galactic winds as the primary source of the enriched and massive circumgalactic medium. Winds from compact starbursts have been observed to flow to distances somewhat greater than ten kiloparsecs, but the circumgalactic medium typically extends beyond a hundred kiloparsecs. Here we report optical integral field observations of the massive but compact galaxy SDSS J211824.06+001729.4. The oxygen [O II] lines at wavelengths of 3726 and 3729 angstroms reveal an ionized outflow spanning 80 by 100 square kiloparsecs, depositing metal-enriched gas at 10,000 kelvin through an hourglass-shaped nebula that resembles an evacuated and limb-brightened bipolar bubble. We also observe neutral gas phases at temperatures of less than 10,000 kelvin reaching distances of 20 kiloparsecs and velocities of around 1,500 kilometres per second. This multi-phase outflow is probably driven by bursts of star formation, consistent with theory. △ Less

Submitted 29 October, 2019; originally announced October 2019.

Comments: Authors' version of a Letter published in Nature on 30 October 2019; journal version at https://www.nature.com/articles/s41586-019-1686-1

Modeling the precession of the warped inner accretion disk in the pulsars LMC X-4 and SMC X-1 with NuSTAR and XMM-Newton

Authors: McKinley C. Brumback, Ryan C. Hickox, Felix S. Fürst, Katja Pottschmidt, John A. Tomsick, Jörn Wilms

Abstract: We present a broad-band X-ray study of the effect of superorbital periods on X-ray spectra and pulse profiles in the neutron star X-ray binaries LMC X-4 and SMC X-1. These two sources display periodic or quasi-periodic variations in luminosity on the order of tens of days which are known to be superorbital, and are attributed to warped, precessing accretion disks. Using joint NuSTAR and XMM-Newton… ▽ More We present a broad-band X-ray study of the effect of superorbital periods on X-ray spectra and pulse profiles in the neutron star X-ray binaries LMC X-4 and SMC X-1. These two sources display periodic or quasi-periodic variations in luminosity on the order of tens of days which are known to be superorbital, and are attributed to warped, precessing accretion disks. Using joint NuSTAR and XMM-Newton observations that span a complete superorbital cycle, we examine the broad-band spectra of these sources and find the shape to be well described by an absorbed power law with a soft blackbody component. Changes in spectral shape and pulse profile shape are periodic with superorbital period, as expected from a precessing disk. We perform X-ray tomography using the changes in pulse profiles to model the geometry and kinematics of the inner accretion disk. Our simple beam and inner disk geometric model indicates that the long term changes in soft pulse shape and phase are consistent with reprocessed emission from a precessing inner disk. △ Less

Submitted 13 December, 2019; v1 submitted 23 September, 2019; originally announced September 2019.

Comments: Updated with version accepted to ApJ. 16 pages, 13 figures

Measuring the Obscuring Column of a Disk Megamaser AGN in a Nearby Merger

Authors: Alberto Masini, Andrea Comastri, Ryan C. Hickox, Michael Koss, Francesca Civano, Murray Brightman, Marcella Brusa, Giorgio Lanzuisi

Abstract: 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… ▽ More 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. △ Less

Submitted 11 July, 2019; originally announced July 2019.

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