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

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

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

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 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

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.

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

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

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

Evolution of black hole and galaxy growth in a semi-numerical galaxy formation model

Authors: Mackenzie L. Jones, R. C. Hickox, S. J. Mutch, D. J. Croton, A. F. Ptak, M. A. DiPompeo

Abstract: We present a simple semi-numerical model designed to explore black hole growth and galaxy evolution. This method builds on a previous model for black hole accretion that uses a semi-numerical galaxy formation model and universal Eddington ratio distribution to describe the full AGN population by independently connecting galaxy and AGN growth to the evolution of the host dark matter halos. We fit o… ▽ More We present a simple semi-numerical model designed to explore black hole growth and galaxy evolution. This method builds on a previous model for black hole accretion that uses a semi-numerical galaxy formation model and universal Eddington ratio distribution to describe the full AGN population by independently connecting galaxy and AGN growth to the evolution of the host dark matter halos. We fit observed X-ray luminosity functions up to a redshift of z ~ 4, as well as investigate the evolution of the Eddington ratio distributions. We find that the Eddington ratio distribution evolves with redshift such that the slope of the low-Eddington accretion rate distribution increases with cosmic time, consistent with the behavior predicted in hydrodynamical simulations for galaxies with different gas fractions. We also find that the evolution of our average Eddington ratio is correlated with observed star formation histories, supporting a picture in which black holes and galaxies evolve together in a global sense. We further confirm the impact of luminosity limits on observed galaxy and halo properties by applying selection criteria to our fiducial model and comparing to surveys across a wide range of redshifts. △ Less

Submitted 1 July, 2019; originally announced July 2019.

Comments: 9 pages, 8 figures, Accepted to ApJ

Resolving the cosmic X-ray background with a next-generation high-energy X-ray observatory

Authors: R. C. Hickox, F. Civano, D. R. Ballantyne, M. Balokovic, P. G. Boorman, W. N. Brandt, R. E. A. Canning, F. Fornasini, P. Gandhi, M. L. Jones, G. B. Lansbury, L. Lanz, G. Lanzuisi, K. K. Madsen, S. Marchesi, A. Masini, T. Ananna, D. Stern, C. Ricci

Abstract: The cosmic X-ray background (CXB), which peaks at an energy of ~30 keV, is produced primarily by emission from accreting supermassive black holes (SMBHs). The CXB therefore serves as a constraint on the integrated SMBH growth in the Universe and the accretion physics and obscuration in active galactic nuclei (AGNs). This paper gives an overview of recent progress in understanding the high-energy (… ▽ More The cosmic X-ray background (CXB), which peaks at an energy of ~30 keV, is produced primarily by emission from accreting supermassive black holes (SMBHs). The CXB therefore serves as a constraint on the integrated SMBH growth in the Universe and the accretion physics and obscuration in active galactic nuclei (AGNs). This paper gives an overview of recent progress in understanding the high-energy (>~10 keV) X-ray emission from AGNs and the synthesis of the CXB, with an emphasis on results from NASA's NuSTAR hard X-ray mission. We then discuss remaining challenges and open questions regarding the nature of AGN obscuration and AGN physics. Finally, we highlight the exciting opportunities for a next-generation, high-resolution hard X-ray mission to achieve the long-standing goal of resolving and characterizing the vast majority of the accreting SMBHs that produce the CXB. △ Less

Submitted 27 May, 2019; originally announced May 2019.

Comments: Science White paper submitted to Astro2020 Decadal Survey; 5 pages, 3 figures, plus references and cover page

Hot Dust-Obscured Galaxies with Excess Blue Light

Authors: R. J. Assef, M. Brightman, D. J. Walton, D. Stern, F. E. Bauer, A. W. Blain, T. Diaz-Santos, P. R. M. Eisenhardt, R. C. Hickox, H. D. Jun, A. Psychogyios, C. -W. Tsai, J. W. Wu

Abstract: Hot Dust-Obscured Galaxies (Hot DOGs) are among the most luminous galaxies in the Universe. Powered by highly obscured, possibly Compton-thick, active galactic nuclei (AGNs), Hot DOGs are characterized by SEDs that are very red in the mid-IR yet dominated by the host galaxy stellar emission in the UV and optical. An earlier study identified a sub-sample of Hot DOGs with significantly enhanced UV e… ▽ More Hot Dust-Obscured Galaxies (Hot DOGs) are among the most luminous galaxies in the Universe. Powered by highly obscured, possibly Compton-thick, active galactic nuclei (AGNs), Hot DOGs are characterized by SEDs that are very red in the mid-IR yet dominated by the host galaxy stellar emission in the UV and optical. An earlier study identified a sub-sample of Hot DOGs with significantly enhanced UV emission. One target, W0204-0506, was studied in detail and, based on Chandra observations, it was concluded that the enhanced emission was most likely due to either extreme unobscured star-formation (${\rm SFR}>1000~M_{\odot}~\rm yr^{-1}$) or to light from the highly obscured AGN scattered by gas or dust into our line of sight. Here, we present a follow-up study of W0204-0506 as well as two more Hot DOGs with excess UV emission. For the two new objects we obtained Chandra/ACIS-S observations, and for all three targets we obtained HST/WFC3 F555W and F160W imaging. We conclude that the excess UV emission is primarily dominated by light from the central highly obscured, hyper-luminous AGN that has been scattered into our line of sight. We cannot rule out, however, that star-formation may significantly contribute to the UV excess of W0204-0506. △ Less

Submitted 4 June, 2020; v1 submitted 10 May, 2019; originally announced May 2019.

Comments: 17 pages, 12 figures, 3 tables. Accepted for publication in The Astrophysical Journal

The halo mass of optically-luminous quasars at z=1-2 measured via gravitational deflection of the cosmic microwave background

Authors: J. E. Geach, J. A. Peacock, A. D. Myers, R. C. Hickox, M. C. Burchard, M. L. Jones

Abstract: We measure the average deflection of cosmic microwave background photons by quasars at $\langle z \rangle =1.7$. Our sample is selected from the Sloan Digital Sky Survey to cover the redshift range $0.9\leq z\leq2.2$ with absolute i-band magnitudes of $M_i\leq-24$ (K-corrected to z=2). A stack of nearly 200,000 targets reveals an 8$σ$ detection of Planck's estimate of the lensing convergence towar… ▽ More We measure the average deflection of cosmic microwave background photons by quasars at $\langle z \rangle =1.7$. Our sample is selected from the Sloan Digital Sky Survey to cover the redshift range $0.9\leq z\leq2.2$ with absolute i-band magnitudes of $M_i\leq-24$ (K-corrected to z=2). A stack of nearly 200,000 targets reveals an 8$σ$ detection of Planck's estimate of the lensing convergence towards the quasars. We fit the signal with a model comprising a Navarro-Frenk-White density profile and a 2-halo term accounting for correlated large scale structure, which dominates the observed signal. The best-fitting model is described by an average halo mass $\log_{10}(M_{\rm h}/h^{-1}M_\odot)=12.6\pm0.2$ and linear bias $b=2.7\pm0.3$ at $z=1.7$, in excellent agreement with clustering studies. We also report of a hint, at a 90% confidence level, of a correlation between the convergence amplitude and luminosity, indicating that quasars brighter than $M_i\lesssim -26$ reside in halos of typical mass ${M_{\rm h}\approx 10^{13}\,h^{-1}M_\odot}$, scaling roughly as ${M_{\rm h}\propto L_{\rm opt}^{3/4}}$ at ${M_i\lesssim-24}$, in good agreement with physically-motivated quasar demography models. Although we acknowledge this luminosity dependence is a marginal result, the observed $M_{\rm h}$-$L_{\rm opt}$ relationship could be interpreted as a reflection of the cutoff in the distribution of black hole accretion rates towards high Eddington ratios: the weak trend of $M_{\rm h}$ with $L_{\rm opt}$ observed at low luminosity becomes stronger for the most powerful quasars, which tend to be accreting close to the Eddington limit. △ Less

Submitted 19 February, 2019; originally announced February 2019.

Comments: Accepted for publication in ApJ

NuSTAR and Keck Observations of Heavily Obscured Quasars Selected by WISE

Authors: Wei Yan, Ryan C. Hickox, Kevin N. Hainline, Daniel Stern, George Lansbury, David M. Alexander, Raphael E. Hviding, Roberto J. Assef, David R. Ballantyne, Michael A. DiPompeo, Lauranne Lanz, Christopher M. Carroll, Michael Koss, Isabella Lamperti, Francesca Civano, Agnese Del Moro, Poshak Gandhi, Adam D. Myers

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

Submitted 6 November, 2018; originally announced November 2018.

Comments: Accepted for publication in the Astrophysical Journal, 11 pages, 10 figures

Investigating the Covering Fraction Distribution of Swift/BAT AGN with X-ray and IR Observations

Authors: Lauranne Lanz, Ryan C. Hickox, Mislav Baloković, Taro Shimizu, Claudio Ricci, Andy D. Goulding, David R. Ballantyne, Franz E. Bauer, Chien-Ting J. Chen, Agnese Del Moro, Duncan Farrah, Michael J. Koss, Stephanie LaMassa, Alberto Masini, Luca Zappacosta

Abstract: We present an analysis of a sample of 69 local obscured Swift/Burst Alert Telescope active galactic nuclei (AGNs) with X-ray spectra from NuSTAR and infrared (IR) spectral energy distributions from Herschel and WISE. We combine this X-ray and IR phenomenological modeling and find a significant correlation between reflected hard X-ray emission and IR AGN emission, with suggestive indications that t… ▽ More We present an analysis of a sample of 69 local obscured Swift/Burst Alert Telescope active galactic nuclei (AGNs) with X-ray spectra from NuSTAR and infrared (IR) spectral energy distributions from Herschel and WISE. We combine this X-ray and IR phenomenological modeling and find a significant correlation between reflected hard X-ray emission and IR AGN emission, with suggestive indications that this correlation may be stronger than the one between intrinsic hard X-ray and IR emissions. This relation between the IR and reflected X-ray emission suggests that both are the result of the processing of intrinsic emission from the corona and accretion disk by the same structure. We explore the resulting implications on the underlying distribution of covering fraction for all AGNs, by generating mock observables for the reflection parameter and IR luminosity ratio using empirical relations found for the covering fraction with each quantity. We find that the observed distributions of the reflection parameter and IR-to-X-ray ratio are reproduced with broad distributions centered around covering fractions of at least ~40%-50%, whereas narrower distributions match our observations only when centered around covering fractions of ~70%-80%. Our results are consistent with both independent estimates of the covering fractions of individual objects and the typical covering fraction obtained on the basis of obscured fractions for samples of AGNs. These results suggest that the level of reprocessing in AGNs, including X-ray reflection, is related in a relatively straightforward way to the geometry of the obscuring material. △ Less

Submitted 8 April, 2019; v1 submitted 6 November, 2018; originally announced November 2018.

Comments: Published in the Astrophysical Journal: Lanz, L. et al. 2019, ApJ, 870, 26; 19 pages; 9 figures

The NuSTAR Extragalactic Surveys: unveiling rare, buried AGNs and detecting the contributors to the peak of the Cosmic X-ray Background

Authors: Alberto Masini, Andrea Comastri, Francesca Civano, Ryan C. Hickox, Christopher M. Carroll, Hyewon Suh, William N. Brandt, Michael A. DiPompeo, Fiona A. Harrison, Daniel Stern

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

Submitted 28 September, 2018; originally announced October 2018.

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

Violent quenching: molecular gas blown to 1000 km/s during a major merger

Authors: J. E. Geach, C. Tremonti, A. M. Diamond-Stanic, P. H. Sell, A. A. Kepley, A. L. Coil, G. Rudnick, R. C. Hickox, J. Moustakas, Yujin Yang

Abstract: We present Atacama Large Millimeter/submillimeter Array observations of a massive (M_stars~10^11 M_Sun) compact (r_e,UV~100 pc) merger remnant at z=0.66 that is driving a 1000 km/s outflow of cool gas, with no observational trace of an active galactic nucleus (AGN). We resolve molecular gas on scales of approximately 1-2 kpc, and our main finding is the discovery of a wing of blueshifted CO(2-1) e… ▽ More We present Atacama Large Millimeter/submillimeter Array observations of a massive (M_stars~10^11 M_Sun) compact (r_e,UV~100 pc) merger remnant at z=0.66 that is driving a 1000 km/s outflow of cool gas, with no observational trace of an active galactic nucleus (AGN). We resolve molecular gas on scales of approximately 1-2 kpc, and our main finding is the discovery of a wing of blueshifted CO(2-1) emission out to -1000 km/s relative to the stars. We argue that this is the molecular component of a multiphase outflow, expelled from the central starburst within the past 5 Myr through stellar feedback, although we cannot rule out previous AGN activity as a launching mechanism. If the latter is true, then this is an example of a relic multiphase AGN outflow. We estimate a molecular mass outflow rate of approximately 300 M_Sun/yr, or about one third of the 10 Myr-averaged star formation rate. This system epitomizes the multiphase 'blowout' episode following a dissipational major merger - a process that has violently quenched central star formation and supermassive black hole growth. △ Less

Submitted 23 October, 2018; v1 submitted 25 July, 2018; originally announced July 2018.

Comments: Accepted for publication in ApJ Letters: The Astrophysical Journal Letters, Volume 864, Number 1

Discovery of pulsation dropout and turn-on during the high state of the accreting X-ray pulsar LMC X-4

Authors: McKinley C. Brumback, Ryan C. Hickox, Matteo Bachetti, Ralf Ballhausen, Felix S. Fürst, Sean Pike, Katja Pottschmidt, John A. Tomsick, Jörn Wilms

Abstract: Two \textit{NuSTAR} observations of the luminous X-ray pulsar LMC X-4 in October and November 2015 captured several bright accretion flares from this source, which has a long history of stable pulse and superorbital behavior. We present a timing analysis of these data in which we detect a rapid pulse turn-on in association with the accretion flares, during which the source reaches super-Eddington… ▽ More Two \textit{NuSTAR} observations of the luminous X-ray pulsar LMC X-4 in October and November 2015 captured several bright accretion flares from this source, which has a long history of stable pulse and superorbital behavior. We present a timing analysis of these data in which we detect a rapid pulse turn-on in association with the accretion flares, during which the source reaches super-Eddington luminosities. Pulsations, which are normally seen from this source, are found to only occur for approximately one hour before and during the bright flares. Beyond one hour before and after the flares, we find pulsations to be weak or nonexistent, with fractional RMS amplitudes of less than 0.05. At the onset of the flare, the pulse profiles exhibit a phase shift of 0.25 cycles that could be associated with a change in the emission geometry. This increase in pulse strength occurring well before the flare cannot be explained by the propeller effect, and potentially offers a connection between the magnetic properties of pulsars that accrete close to their Eddington limits and ULX pulsars. △ Less

Submitted 28 June, 2018; v1 submitted 25 June, 2018; originally announced June 2018.

Comments: 7 pages, 4 figures, accepted to ApJL, updated references

Obscured Active Galactic Nuclei

Authors: Ryan C. Hickox, David M. Alexander

Abstract: Active Galactic Nuclei (AGN) are powered by the accretion of material onto a supermassive black hole (SMBH), and are among the most luminous objects in the Universe. However, the huge radiative power of most AGN cannot be seen directly, as the accretion is "hidden" behind gas and dust that absorbs many of the characteristic observational signatures. This obscuration presents an important challenge… ▽ More Active Galactic Nuclei (AGN) are powered by the accretion of material onto a supermassive black hole (SMBH), and are among the most luminous objects in the Universe. However, the huge radiative power of most AGN cannot be seen directly, as the accretion is "hidden" behind gas and dust that absorbs many of the characteristic observational signatures. This obscuration presents an important challenge for uncovering the complete AGN population and understanding the cosmic evolution of SMBHs. In this review we describe a broad range of multi-wavelength techniques that are currently employed to identify obscured AGN, and assess the reliability and completeness of each technique. We follow with a discussion of the demographics of obscured AGN activity, explore the nature and physical scales of the obscuring material, and assess the implications of obscured AGN for observational cosmology. We conclude with an outline of the prospects for future progress from both observations and theoretical models, and highlight some of the key outstanding questions. △ Less

Submitted 12 June, 2018; originally announced June 2018.

Comments: 50 pages, 16 figures. Invited review to appear in Annual Reviews of Astronomy and Astrophysics, Volume 56. An animated version of the schematic of AGN spectral energy distributions shown in Figure 4 is available at http://www.dartmouth.edu/~hickox/hickox_alexander_AGN.php