A publishing partnership

The following article is Open access

An X-Ray Census of Fast Radio Burst Host Galaxies: Constraints on Active Galactic Nuclei and X-Ray Counterparts

, , , , , , , , ,

Published 2023 November 14 © 2023. The Author(s). Published by the American Astronomical Society.
, , Citation T. Eftekhari et al 2023 ApJ 958 66DOI 10.3847/1538-4357/acf843

Download Article PDF
DownloadArticle ePub

You need an eReader or compatible software to experience the benefits of the ePub3 file format.

0004-637X/958/1/66

Abstract

We present the first X-ray census of fast radio burst (FRB) host galaxies to conduct the deepest search for active galactic nuclei (AGN) and X-ray counterparts to date. Our sample includes seven well-localized FRBs with unambiguous host associations and existing deep Chandra observations, including two events for which we present new observations. We find evidence for AGN in two FRB host galaxies based on the presence of X-ray emission coincident with their centers, including the detection of a luminous (LX ≈ 5 × 1042 erg s−1) X-ray source at the nucleus of FRB 20190608B's host, for which we infer an SMBH mass of MBH ∼ 108M and an Eddington ratio Lbol/LEdd ≈ 0.02, characteristic of geometrically thin disks in Seyfert galaxies. We also report nebular emission-line fluxes for 24 highly secure FRB hosts (including 10 hosts for the first time), and assess their placement on a BPT diagram, finding that FRB hosts trace the underlying galaxy population. We further find that the hosts of repeating FRBs are not confined to the star-forming locus, contrary to previous findings. Finally, we place constraints on associated X-ray counterparts to FRBs in the context of ultraluminous X-ray sources (ULXs), and find that existing X-ray limits for FRBs rule out ULXs brighter than LX ≳ 1040 erg s−1. Leveraging the CHIME/FRB catalog and existing ULX catalogs, we search for spatially coincident ULX–FRB pairs. We identify a total of 28 ULXs spatially coincident with the localization regions for 17 FRBs, but find that the DM-inferred redshifts for the FRBs are inconsistent with the ULX redshifts, disfavoring an association between these specific ULX–FRB pairs.

Export citation and abstractBibTeXRIS

1. Introduction

Fast radio bursts (FRBs) are bright, millisecond-duration flares of coherent radio emission (Lorimer et al. 2007; Thornton et al. 2013) whose large dispersion measures (DMs) point to an extragalactic origin. Despite over 500 FRBs published in the literature to date (e.g., Macquart et al. 2020; CHIME/FRB Collaboration et al. 2020, CHIME/FRB Collaboration et al. 2023; Law et al. 2023), the sources responsible for producing FRBs are highly debated, in large part due to the dearth of well-localized events (Eftekhari & Berger 2017) as well as the fact that some FRBs have been known to repeat (so-called "repeaters;" Spitler et al. 2016) while others appear as one-off events (Shannon et al. 2018; CHIME/FRB Collaboration et al. 2021). Although repeating and apparently nonrepeating FRBs exhibit a number of distinctions in their spectrotemporal properties (Pleunis et al. 2021)—and potentially their host galaxy properties (Gordon et al. 2023)—analyses of the FRB population cannot exclude the possibility that all FRBs intrinsically repeat (James 2023).

To reconcile this apparent diversity in burst properties and environments, a wide range of theoretical models have been put forth for FRB progenitors, ranging from compact object mergers (Margalit et al. 2019; Sridhar et al. 2021), and flares from highly magnetized stars (Margalit & Metzger 2018), to the super-Eddington accretion onto a neutron star or black hole (Sridhar et al. 2021, 2022), or even more exotic phenomena invoking cosmic strings (Costa et al. 2018) and primordial black holes (Abramowicz et al. 2018). At present, prevailing theoretical models invoke magnetars, whose extremely large magnetic fields are capable of storing enormous magnetospheric energies released in FRB flashes. Indeed, the discovery of a luminous, millisecond-duration burst from the known Galactic magnetar SGR 1935+2154 (CHIME/FRB Collaboration et al. 2020; Bochenek et al. 2020) supports a connection between extragalactic magnetars and FRBs (Margalit et al. 2020). Recent studies suggest that magnetars which formed via multiple formation channels—and not strictly through core-collapse supernovae—likely contribute to the known FRB population (e.g., Margalit et al. 2019; Kirsten et al. 2022; Pelliciari et al. 2023).

Technical upgrades to a number of FRB experiments in recent years have facilitated a small, but growing number of (sub)arcsecond localizations and robust host galaxy identifications for approximately three dozen events (e.g., Bhardwaj et al. 2021; Kirsten et al. 2022; Gordon et al. 2023; Law et al. 2023; Ravi et al. 2023; Ryder et al. 2023; Sharma et al. 2023). The first population studies demonstrate that FRB hosts span a wide range of stellar masses and star formation rates (SFRs) and trace the underlying population of field galaxies (Heintz et al. 2020; Bhandari et al. 2022; Gordon et al. 2023).

These early host galaxy compilations have also revealed that a nonnegligible fraction of FRB hosts occupy a distinct region of the Baldwin–Phillips–Terlevich (BPT; Baldwin et al. 1981) diagram, pointing to emission-line ratios in excess of typical star-forming galaxies (Heintz et al. 2020; Bhandari et al. 2022; Ibik et al. 2023). In particular, nonrepeating FRBs seem to indicate a preference for a subpopulation of low-ionization nuclear emission-line region (LINER) galaxies, while repeating FRBs lie exclusively along the star-forming branch. This distinction suggests that nonrepeating FRBs may be preferentially located in environments with enhanced levels of photoionization.

The nature of LINER galaxies, however, and the source of their ionizing photons have been under debate since their discovery (Heckman 1980). Existing theories associate LINERs with low-luminosity active galactic nuclei (AGN) or stellar ionization from evolved post–asymptotic giant branch stars (Singh et al. 2013). For FRB hosts, the latter scenario would implicate older stellar populations and hence FRB progenitors formed through channels with large average delay times between star formation and FRB source formation (e.g., short-duration gamma-ray bursts). Conversely, if the hard radiation fields observed in some FRB hosts are driven by central AGN, pointing to an elevated fraction of AGN in FRB hosts, this would lend support for theories that presuppose an enhanced rate of compact objects formed in AGN accretion disks (Cantiello et al. 2021; Perna et al. 2021a, 2021b; Jermyn et al. 2021) or FRB progenitor models that invoke AGN (Zhang 2017). Indeed, several of these theories have been put forth for relativistic transients, including both long- and short-duration gamma-ray bursts, but could be applied here given the higher stellar densities in AGN disks in general. Moreover, with only a small sample of FRB hosts that have been analyzed in this context (Heintz et al. 2020; Bhandari et al. 2022), a larger study is warranted to quantify the true occurrence rate of AGN (and LINER galaxies) in FRB hosts.

Concurrent to discerning the origin of the hard radiation fields prevalent in FRB hosts, X-ray observations can uniquely probe high-energy counterparts to FRBs. Constraints on periodicity in the burst rate from two repeating FRBs (Rajwade et al. 2020; CHIME/FRB Collaboration et al. 2020) have prompted theories that they are powered by relativistic flares along precessing jets associated with accretion onto compact objects (Sridhar et al. 2021). The observed (millisecond) durations and the extreme luminosities of typical FRBs, as well as the properties of their host galaxies, suggest that FRBs may be produced by stellar-mass ultraluminous X-ray sources (ULXs) powered by super-Eddington accretion onto a compact object (Sridhar et al. 2021). Within this framework, bright ULX-like X-ray emission should be associated with at least some FRBs. While deep limits on X-ray counterparts exist for a small number of well-localized FRB sources (e.g., Scholz et al. 2020; Kirsten et al. 2022; Pearlman et al. 2023) and targeted searches for FRBs from nearby galaxies hosting ULXs have been conducted (Pelliciari et al. 2023), there have been no concerted efforts thus far to search for ULXs coincident with the localization regions of the hundreds of FRBs detected across the sky.

Here we report the X-ray properties for all well-localized FRBs with existing Chandra X-ray observations (including newly obtained data for two events) to probe the presence of both AGN in FRB hosts and X-ray counterparts at the FRB locations. Our sample consists of seven FRBs, including five repeating sources. In Section 2, we present our observations and data reduction, which includes new Karl G. Jansky Very Large Array (VLA) observations for two FRBs. In Section 3, we combine radio and X-ray observations of our sample of hosts to constrain the multiwavelength properties of AGN in FRB host galaxies. We additionally present an updated BPT analysis for a uniformly modeled sample of 22 FRB host galaxies (in addition to two hosts from the literature). We leverage the X-ray observations to place limits on associated X-ray counterparts to FRBs in Section 4, and we further explore a possible connection between FRBs and ULXs by searching for FRBs beyond the original sample that are spatially coincident with cataloged ULXs. We summarize our results in Section 5. Throughout the paper, we use the latest Planck cosmological parameters for a flat Λ cold dark matter (ΛCDM) Universe, with H0 = 67.66 km s−1 Mpc−1, Ωm = 0.310, and Ωλ = 0.690 (Planck Collaboration et al. 2020).

2. Observations

2.1. Sample Selection

Our sample of FRB hosts includes all well-localized (∼a few to hundreds of milliarcseconds) events with unambiguous host associations (a probabilistic association of transients to their hosts (PATH); Aggarwal et al. 2021; with posterior probability ≳ 90%) and existing Chandra data. The resulting sample spans a redshift range of z ∼ 0.0–0.2418 30 and includes five FRBs with previously published X-ray data which we reanalyze here and two for which we present new Chandra observations (see Table 1). The FRB hosts in our sample span a wide range in both stellar masses () and SFRs (∼0.04–7.03 M yr−1), similar to the general population of FRB hosts (Gordon et al. 2023), although are biased toward repeating events (five out of seven FRBs).

Table 1. X-Ray and Radio Observations of the Host Galaxies of Localized FRBs

FRB z Average Observation EpochExposure Timelog(L2–10 keV)log(L5 GHz) NH References
  (MJD)(ks)(erg s−1)(erg s−1)(1022 cm−2) 
FRB 20121102A0.19272016 Nov 3099.9 a <41.239.10.631, 2, 3, 4
FRB 20180916B0.03372019 Dec 1032.7 b <40.0<36.10.615, 6
FRB 20190520B0.24182020 Sept 1414.9<41.939.10.267, 8
FRB 20190608B0.11782021 Dec 0619.8 38.10.06This work, 9
FRB 20200120E0.00082010 Jan 241205.3 c 37.00.0410, 11, 12
FRB 20200430A0.16082022 Apr 2434.8 d <41.4<37.70.04This work
FRB 20201124A0.09822021 Apr 2029.7 38.20.4513

Notes. X-ray luminosity limits and detections correspond to the host center. Limits at the FRB position for each source are given in Section 2. Limits correspond to 3σ.

a The total exposure time was spread across four epochs spanning 2015 Nov 23–2017 Nov 01. b The total exposure time was spread across two epochs spanning 2019 Dec 03–2019 Dec 18. c The total exposure time was spread across 58 epochs spanning 2000 May 07–2022 Jun 04. d The total exposure time was spread across two epochs spanning 2022 Apr 23–2022 Apr 24. References: (1) Scholz et al. (2016); (2) Chatterjee et al. (2017); (3) Scholz et al. (2017); (4) PI: S. Bogdanov; (5) Scholz et al. (2020); (6) Marcote et al. (2020); (7) J. Sydnor et al. (2023, in preparation); (8) Niu et al. (2022); (9) Bhandari et al. (2020); (10) Kirsten et al. (2022); (11) see Facilities statement; (12) Miller et al. (2010); and (13) Piro et al. (2021).

Download table as:  ASCIITypeset image

2.2. X-Ray Detection of the AGN in the Host Galaxy of FRB 20190608B

We obtained X-ray observations with the Chandra Advanced CCD Imaging Spectrometer (ACIS-S; Obs ID: 25251; PI: T. Eftekhari) of the field of FRB 20190608B (a nonrepeating FRB discovered and localized by the Australian Square Kilometre Array Pathfinder (ASKAP; Macquart et al. 2020) on UT 2021 December 6 with an exposure time of 19.8 ks. We analyzed the data using the CIAO software package (v4.13) and followed standard ACIS data filtering.

We detect a bright X-ray source offset from the FRB position by ≈1farcs6 with high significance (71σ) using a blind search with wavdetect (Figure 1). The source is detected in an elliptical region of size 2farcs8 × 2farcs5 (PA = 32°) centered at R.A. = 22h16m04fs858, decl. = −7d53m56fs297. To estimate the background, we create a 35'' radius region near the center of the chip and free from any obvious sources. We find a net source count rate of (108.3 ± 7.5) × 10−4 counts s−1 (0.5–8 keV). Given the spatial coincidence of the source with the host galaxy center and the unextended nature of the X-ray emission, we consider the source a putative AGN.

Figure 1. Refer to the following caption and surrounding text.

Figure 1. Left: Chandra X-ray detection (0.5–8 keV) of the host galaxy of FRB 20190608B. Contours from the Gaussian-smoothed HST (F300X) data (Chittidi et al. 2021) are overlaid in green. The red circle corresponds to the 1σ positional uncertainty on the FRB position (Day et al. 2021). Right: X-ray spectrum (gray points) of the AGN in the host galaxy of FRB 20190608B, modeled with an absorbed power law (green curve) with a best-fit power-law index Γ = 1.9 ± 0.3.

Standard image High-resolution image

We rebin the spectrum across the 2farcs8 × 2farcs5 source region using 2 counts bin–1 and fit the data using an absorbed power-law model (tbabs*ztbabs*pow in Xspec; Arnaud 1996) with a fixed Galactic absorption column density NH,MW = 5.9 ×1020 cm−2 (Willingale et al. 2013), solar abundances from Asplund et al. (2009), and use the maximum likelihood C-statistic for Poisson statistics (c-stat; Cash 1979). We find a photon index of Γ = 1.9 ± 0.3 (68% confidence) and an intrinsic column density NH = (4.7 ± 0.4) × 1021 cm−2 for the AGN in the host galaxy. The model fit is shown in the right panel of Figure 1. To measure the quality of the fit, we use the Cramer–von Mises (CvM) test statistic (statistic test cvm; Cramér 1928; Mises 1928) and XSPEC's "goodness" simulations to simulate realizations of the power law and determine what percentage of simulations has test statistics smaller than that of the data. A value of ∼50% indicates that the observed spectrum was produced by the model. We find that 43% of the realizations have a lower CvM statistic than the data, indicating the model is a good fit. Applying these spectral parameters to the net count rate and employing cflux, we find an unabsorbed X-ray flux FX = (1.4 ± 0.2) ×10−13 erg cm−2 s−1 (2–10 keV), corresponding to an X-ray luminosity LX = (5.4 ± 0.8) × 1042 erg s−1 at the redshift of the FRB host, z = 0.1178 (see Table 1).

To determine the location of the X-ray source relative to the host galaxy center, any substructure within the host galaxy, and the position of FRB 20190608B, we use the Hubble Space Telescope (HST) F300X image of the host galaxy (Chittidi et al. 2021), which enables a precise host centroid. Since there are not enough sources in common between the HST and Chandra images, we use a Dark Energy Camera Legacy Survey image (DECaLS; DR10; Dey et al. 2019) to obtain a common absolute astrometric frame. First, we perform relative astrometry between the HST image using Source Extractor (Bertin & Arnouts 1996) and point sources in common. We then register the Chandra field of view to the DECaLS astrometric system using sources in common. With both the HST and Chandra images tied to the same astrometric frame, we find that the X-ray source is coincident with the host galaxy center (Figure 1).

The observed properties of the X-ray source, namely its spatial coincidence with the host galaxy center, luminous X-ray emission at the level of LX = 5.4 × 1042 erg s−1, and an X-ray spectrum that is well characterized by an absorbed power law, indicate that the X-ray source is powered by an AGN. We note that we do not detect any counts at the FRB position, corresponding to a 3σ upper limit on the absorbed X-ray flux of FX < 3 × 10−15 erg cm−2 s−1 (2–10 keV).

2.3. Constraints on an AGN in the Host Galaxy of FRB20200430A

We obtained Chandra observations (PI: T. Eftekhari) of FRB 20200430A (a nonrepeating FRB detected and localized by ASKAP; Heintz et al. 2020) in two separate epochs with exposure times of 22.9 and 11.9 ks on UT 2022 April 23 (Obs ID: 25252) and 24 (Obs ID: 26397), respectively. We do not blindly detect X-ray emission at or near the position of FRB 20200430A in either epoch using wavdetect. We therefore generate a coadded, exposure-corrected image using merge_obs and perform a targeted extraction around the host center using an aperture radius of 1''. We perform background estimation in a similar manner to FRB 20190608B. The source is not detected in the merged event file corresponding to a 3σ upper limit on the count rate of 3.6 × 10−4 counts s−1 (assuming Poisson statistics; Gehrels 1986). For Γ = 2 (typical of AGN in nearby galaxies; Ho 2008) and NH,MW = 3.5 ×1020 cm−2, the 3σ limit on the 2−10 keV absorbed X-ray flux is FX < 3.2 × 10−15 erg cm−2 s−1, corresponding to an X-ray luminosity LX < 2.8 × 1041 erg s−1. We derive the same flux limit at the FRB position located ∼1'' from the host center.

2.4. Reanalysis of Archival Data

We also searched the positions of all localized FRBs in the Chandra archive to compile a complete list of X-ray observations for well-localized events (see Table 1). We find five FRBs with robust host galaxy associations with archival Chandra data (FRBs 20121102A, 20180916B, 20190520B, 20200120E, and 20201124A). With the exception of FRB20200120E, all of these were targeted for FRB science. For uniformity, we reanalyze the data using the above methodology. In all cases, hydrogen column densities are derived using the method of Willingale et al. (2013) and X-ray fluxes are estimated assuming a photoelectrically absorbed power-law spectrum with Γ = 2.

For FRB 20121102A, we combine four epochs of existing Chandra observations, including one ACIS-S exposure from 2015 November (Obs ID:18717; PI: P. Scholz) and three ACIS-I exposures from 2016 November and 2017 January and November (Obs ID: 19286, 19287, and 19288; PI: S. Bogdanov) corresponding to a total exposure time of 99.99 ks. Only two photons are detected within a 1'' radius aperture centered on the host position, consistent with a nondetection. The 3σ upper limit on the count rate from Poisson statistics is <1.09 × 10−4 counts s−1. For NH,MW = 6.3 × 1021 cm−2, the 3σ limit on the absorbed flux is FX < 1.2 × 10−15 erg cm−2 s−1 (2–10 keV), corresponding to an X-ray luminosity LX < 1.6 × 1041 erg s−1. We derive the same flux limit at the FRB position.

We independently analyze the two existing epochs of Chandra observations for FRB 20180916B following Scholz et al. (2020). We find four counts in a 1'' radius aperture centered on the FRB host position. Accounting for the background, the total source counts are 2.02 ± 1.01 (0.8σ significance), consistent with a nondetection. The 3σ Poisson upper limit on the count rate is 3.3 × 10−4 counts s−1. For NH,MW = 6.1 × 1021 cm−2, the 3σ limit on the absorbed flux is FX < 3.4 × 10−15 erg cm−2 s−1 (2–10 keV), corresponding to an X-ray luminosity LX < 9.9 × 1041 erg s−1. At the FRB position (offset from the host galaxy center by ≈7farcs8), we find a single count in a 1'' radius aperture, consistent with a nondetection. The 3σ Poisson upper limit on the count rate is therefore 2.7 × 10−4 counts s−1, corresponding to a 3σ limit on the absorbed flux of FX < 2.8 × 10−15 erg cm−2 s−1 (2–10 keV).

For FRB 20190520B, we analyze a single ACIS-S 14.9 ks exposure (Obs ID: 22370; PI: K. Aggarwal; Sydnor et al. 2023). We do not detect any counts in a 1'' radius aperture centered on the host center. The 3σ upper limit on the count rate is 4.4 × 10−4 counts s−1. For NH,MW = 2.56 × 1021 cm−2, the 3σ limit on the X-ray flux is FX < 3.7 × 10−15 erg cm−2 s−1 (2–10 keV), corresponding to an X-ray luminosity LX < 8.6 × 1041 erg s−1. We derive the same flux limit at the FRB position.

For FRB 20201124A, we detect five counts in a 1farcs5 radius aperture centered on the host position in a 29.7 ks exposure (Obs ID:25016; PI: L. Piro; Piro et al. 2021). We find a net source count rate of (1.6 ± 0.7) × 10−4 counts s−1 (2.4σ significance), consistent with a marginal detection, as reported in Piro et al. (2021). We convert this count rate to an X-ray flux with NH,MW = 4.5 × 1021 cm−2 and find FX = (1.6 ± 0.6) × 10−15 erg cm−2 s−1 (2–10 keV). At the distance of FRB 20201124A (z = 0.0982), the corresponding X-ray luminosity is LX = (4 ± 1.7) × 1040 erg s−1, consistent with the values reported in Piro et al. (2021) and a star formation origin, as demonstrated in their work. Given the small spatial offset between the FRB and the host center (0farcs7), we adopt the host galaxy X-ray flux as an upper limit for X-ray emission at the FRB position.

In the case of FRB 20200120E (located in M81, which has been observed extensively at X-ray wavelengths), we separately compile and merge all existing Chandra ACIS observations that cover the host and FRB position 31 using the python code superchandra.py (Kilpatrick 2023), which allows us to easily query and analyze all Chandra observations at a given set of input coordinates. We note that none of the ACIS observations comprising our merged event files were targeted for FRB science. We derive an absorbed X-ray flux at the host center (15'' aperture radius) of FX = (1.52 ± 0.01) × 10−11 erg cm−2 s−1 (2–10 keV) for NH,MW = 4 × 1020 cm−2, corresponding to an X-ray luminosity LX = (2.36 ± 0.02) × 1040 erg s−1. Here, NH,MW is derived with the superchandra.py script using the known linear relation between the line-of-sight optical extinction AV (Schlafly & Finkbeiner 2011) and the hydrogen column density (Güver & Özel 2009). We do not detect any counts in a 1'' radius aperture at the FRB position, corresponding to a 3σ limit on the absorbed X-ray flux of FX < 8.6 × 10−15 erg cm−2 s−1 (2–10 keV; see also Pearlman et al. 2023).

2.5. Radio Continuum Observations

We obtained radio observations of FRBs 20200430A (20A-157; PI: S. Bhandari) and 20190608B (22B-115; PI: T. Eftekhari) with the VLA at 6 GHz (C band) and 22 GHz (K band) on UT 2020 September 4 and UT 2022 April 24, respectively. The field of FRB 20200430A was observed in the B-configuration with a total on-source integration time of ∼1 hr. FRB 20190608B was observed in the A-configuration with a total on-source integration time of 56 minutes. Both observations utilized the 3-bit samplers providing the full 4 GHz of bandwidth at C band and 8 GHz of bandwidth at K band, not accounting for the excision of edge channels and RFI. For FRB 20200430A, we used 3C286 for the bandpass and flux density calibration and J1504+1029 for the complex gain calibration. For FRB 20190608B, we used 3C147 for the bandpass and flux density calibration and J2229-0832 for the complex gain calibration.

In both cases, we imaged the pipeline-calibrated measurement sets using 2048 × 2048 pixels at a scale of 0farcs2 pixel−1 and 0farcs02 pixel−1 at 6 and 22 GHz, respectively, using multifrequency synthesis (MFS; Sault & Wieringa 1994) and w-projection with 128 planes (Cornwell et al. 2008). We do not detect radio emission at the host center or the FRB location for either FRB, corresponding to rms values of 4.32 μJy (6 GHz) and 6.14 μJy (22 GHz) for FRB 20200430A and FRB 20190608B, respectively.

We note that FRB 20190608B was observed previously with the Australia Telescope Compact Array (ATCA) under project code C3211 (PI: R. M. Shannon) and with the VLA (19A-121; PI: S. Bhandari) at ∼6 GHz (Bhandari et al. 2020). A radio source is marginally detected in the ATCA observations with a flux density of ∼65 ± 15 μJy at 5.5 GHz. The VLA observations, taken in A-configuration, revealed resolved radio continuum emission with a peak brightness of ∼27 μJy beam−1 near the center of the host. The inferred SFR from the 5.5 GHz ATCA detection of (using the relation of Greiner et al. 2016) is below the optically inferred SFR of (Gordon et al. 2023), but we note that standard relations between radio flux and star formation result in differences of order a factor of two (e.g., Yun & Carilli 2002; Murphy et al. 2011). Moreover, given that the emission spans a range of spatial scales (as shown by comparing the peak brightness seen in the ATCA and VLA images), we emphasize that the radio flux densities should be regarded with caution. Indeed, for SFRradio = 3.9 M yr−1, the expected flux density at 22 GHz assuming a spectral index α = 0.7 is ∼24 μJy, which is above our 3σ limit of 18 μJy. This "missing flux" at higher frequencies is consistent with the higher angular resolution in our VLA 22 GHz observations (∼0farcs1) relative to the ATCA 5.5 GHz observations (20'' × 2''). Thus, we conclude that the most likely origin for the extended radio emission is star formation within the host galaxy.

We further supplement our sample of seven FRB hosts with radio continuum observations from the literature. This includes radio detections of the persistent radio sources (PRSs) associated with FRB 20121102A (F6 GHz = 203 μJy; Chatterjee et al. 2017) and FRB 20190520B (F3 GHz = 202 μJy; Niu et al. 2022); radio detections of the host galaxy of FRB 20201124A (F6 GHz = 221 μJy; Ravi et al. 2022); and a radio limit for FRB 20180916B (F1.6 GHz < 18 μJy (3σ); Marcote et al. 2020). For FRB 20200120E (M81), we adopt the average flux density (F8.4 GHz = 103.2 mJy) over seven epochs of observations as reported in Miller et al. (2010). We list the radio luminosity for each source (k-corrected to rest-frame 5 GHz assuming α = 0.7) in Table 1.

3. Multiwavelength Properties of the AGN in FRB Host Galaxies

3.1. The AGN in the Host of FRB20190608B

The host galaxy of FRB 20190608B (SDSS J221604.90-075355.9) is a grand design SB(r)c galaxy and is among the most massive (≈3.6 × 1010 M) FRB hosts discovered to date (Gordon et al. 2023). While the FRB location is coincident with a knot of UV emission that is likely due to star formation (Chittidi et al. 2021), model fits to the full spectral energy distribution (SED) indicate it is one of two FRB hosts with a poststarburst star formation history in which the galaxy has undergone a recent quenching event (Gordon et al. 2023). Such galaxies commonly host supermassive black holes (SMBHs) which are believed to regulate star formation via AGN-driven feedback (Martín-Navarro et al. 2018). Here, we report the first detection of luminous X-ray emission spatially coincident with the center of the host of FRB 20190608B. An AGN origin is consistent with the previous detection of broad Hα emission in the optical spectrum, indicative of a Type I Seyfert galaxy (Stern & Laor 2012; Chittidi et al. 2021). As such, the host of FRB 20190608B is the only one known to harbor a Type I AGN. Our X-ray detection provides further evidence for an AGN origin (Figure 1).

In Figure 2, we plot the rest-frame radio-to-X-ray luminosity ratio ) as a function of rest-frame X-ray luminosity LX = L2–10 keV for FRB 20190608B, where we conservatively adopt the 5.5 GHz flux density of the host previously ascribed wholly to star formation as an upper limit on any radio emission due to an AGN (see Section 2). For comparison, we also plot the values for several AGN samples from the literature, including low-luminosity AGN (LLAGN; d < 50 Mpc; She et al. 2017; Gross et al. 2019), X-ray-selected AGN (Ballo et al. 2012; LaMassa et al. 2016), quasars (Kellermann et al. 1989; Piconcelli et al. 2005; Gross et al. 2019), Seyferts (Ho & Peng 2001; Terashima & Wilson 2003; Panessa et al. 2007), and low-luminosity radio galaxies (LLRGs; 32 Chiaberge et al. 2005; Balmaverde & Capetti 2006). Vertical dashed lines at and roughly divide AGN into three luminosity regimes corresponding to LLAGN, Seyferts, and quasars (Brusa et al. 2007). We note however that Seyfert galaxies, particularly Type 2 Seyferts, can extend to much lower X-ray luminosities (Panessa et al. 2007), as shown in Figure 2.

Figure 2. Refer to the following caption and surrounding text.

Figure 2. Rest-frame radio-to-X-ray luminosity ratio (RX ) as a function of rest-frame X-ray luminosity (LX ) for FRB hosts (teal circles). Open circles correspond to FRBs with coincident PRSs, where we conservatively adopt the PRS luminosity as an upper limit on the radio luminosity of an associated AGN. Shown for comparison are values for field AGN, including LLAGN (She et al. 2017; Gross et al. 2019), X-ray selected AGN (Ballo et al. 2012; LaMassa et al. 2016), quasars (Kellermann et al. 1989; Piconcelli et al. 2005; Gross et al. 2019), Seyferts (Ho & Peng 2001; Terashima & Wilson 2003; Panessa et al. 2007), and LLRGs (Chiaberge et al. 2005; Balmaverde & Capetti 2006). Vertical dashed lines denote the nominal divisions between LLAGN and Seyferts at log(LX ∼ 41.76 erg s−1) and Seyferts and quasars at log(LX ∼ 43.76 erg s−1). The horizontal dashed line at RX = 4.5 roughly divides radio-loud and radio-quiet AGN.

Standard image High-resolution image

Our X-ray detection of the host of FRB 20190608B, coupled with a limit on the radio emission, places this source in the region corresponding to Seyfert galaxies, consistent with its optical classification. The value of RX < −4.6 is just below the nominal division between radio-loud and radio-quiet AGN (RX = −4.5; Terashima & Wilson 2003), although a diverse range of morphologies and properties for radio AGN challenge a simple classification (Panessa et al. 2019). The increase in RX for LLAGN with LX ≲ 1042 erg s−1 may indicate the presence of an advection-dominated accretion flow in these objects (Terashima & Wilson 2003), as was inferred for the LLAGN in M81 (Xu & Cao 2009) and consistent with its placement in this diagram.

We use the velocity dispersion at the center of the galaxy as derived from a Gaussian fit to the Hβ line (σ = 108 km s−1; Chittidi et al. 2021) and the Mσ relation of Zubovas & King (2012) for spiral galaxies to estimate the mass of the central SMBH. We infer an SMBH mass of MBH ∼ 108 M. For a derived stellar mass of 3.6 × 1010 M, the inferred black hole mass is consistent with scaling relations between central black hole masses and galaxy stellar masses (Kormendy & Ho 2013; Reines & Volonteri 2015). Finally, adopting a bolometric correction of Lbol = 16 LX (Ho 2008), we find an Eddington ratio of Lbol/LEdd ≈ 0.02, typical of Seyfert galaxies (Ho 2008). Moreover, this is consistent with the inferred value of RX < −4.6 (and placement within the radio-quiet regime) for the host of FRB 20190608B and reflects the fact that geometrically thin disks, which are prevalent in Seyfert galaxies, are less efficient at jet production when compared to LLAGN (Terashima & Wilson 2003).

3.2. X-Ray Constraints on AGN in the FRB Host Population

Based on the Chandra X-ray observations presented here and compiled from the literature, we find evidence for AGN in two FRB host galaxies: the host of FRB 20190608B (Section 3.1) and the host of FRB 20200120E (M81), which was already known to harbor an LLAGN (Miller et al. 2010). We note that X-ray emission is also detected from the host of FRB 20201124A but is consistent with a star formation origin (Piro et al. 2021). The detection of AGN in two out of seven FRB hosts corresponds to an occurrence rate of ∼29%. Alternatively, if we consider the hosts of repeating and nonrepeating events separately, we find occupation fractions of ∼20% and ∼50%, respectively. Conversely, the AGN fraction in the local Universe is low, ranging from 0.1% at z ∼ 0.2 to 1% at z ∼ 0.7 (Martini et al. 2009). At face value, it seems that host galaxies which harbor AGN are overrepresented compared to galaxies at similar redshifts. However, existing deep X-ray observations of FRB host galaxies are extremely sparse, with uncharacterized selection effects. Thus a robust assessment of the true AGN occupation fraction in FRB hosts is not yet tenable. Given that the two AGN in FRB hosts are found in both repeating and nonrepeating FRB host galaxies, we observe no trend in the prevalence of AGN with FRB type, albeit based on a small sample size.

To place constraints on the presence of AGN in the absence of X-ray detections, we combine the Chandra limits with radio observations (see Table 1) to plot the radio-to-X-ray luminosity RX against the X-ray luminosity LX for the remainder of sources in Figure 2. We find that the X-ray limits for the majority of FRBs in our sample allow us to rule out any AGN with LX ≲ 1041–1042 erg s−1, including some LLAGN. However, we cannot rule out the low-luminosity end of the LLAGN population centered around LX ∼ 1039 erg s−1, well below the sensitivity of existing X-ray observations for FRBs. For FRB 20200430A, our X-ray limit of LX < 2.5 × 1041 erg s−1 allows us to rule out luminous quasars and most X-ray selected AGN, but we note that we cannot exclude the possibility of a Seyfert host.

Existing radio limits for FRB hosts, coupled with X-ray nondetections, correspond to radio-to-X-ray luminosity ratios of RX > −4, and hence do not reach the regime of radio-quiet AGN. For FRBs 20121102A and 20190520B, the detection of coincident PRSs allows us to leverage the PRS luminosities as conservative upper limits on the radio luminosity of an associated AGN. While we cannot exclude the possibility of radio-loud LLAGN in these FRB hosts, we note that such sources are intrinsically rare in dwarf galaxies (e.g., Reines et al. 2020). We further note that for FRB 20200120E (M81), and indeed all FRBs in our sample (with the exception of FRB 20190608B), the radio and X-ray luminosities and limits are consistent with LLAGN.

3.3. Optical Emission-line Diagnostics for FRB Hosts

To determine the dominant source of ionization and distinguish between star-forming, LINER, and AGN galaxies in FRB hosts, we extract the nebular line fluxes for 22 highly secure FRB host galaxies (including 10 hosts for the first time; see Table 2) with high-quality spectra using the Prospector SED modeling code (Johnson et al. 2021) and following the prescription and sample selection of Gordon et al. (2023). We also include literature values for M81, the host of FRB 20200120E (Ho et al. 1996), and FRB 20220912A (Ravi et al. 2023). 33 We note that while SED fits exist for the host of FRB 20190711A, it is not included in our sample due to poor spectral data quality which precludes robust line emission diagnostics (Gordon et al. 2023). In Figure 3, we plot the emission-line ratios [O iii] λ5007/Hβ λ4861 and [N ii] λ6583/Hα λ6563 for FRB hosts on a BPT diagram (Baldwin et al. 1981) against the distribution of nearby (0.02 < z < 0.4) emission-line galaxies from the Sloan Digital Sky Survey (SDSS).

Figure 3. Refer to the following caption and surrounding text.

Figure 3. BPT classification diagram for FRB host galaxies, including repeating (orange) and nonrepeating (teal) events. Sources that are highlighted in the text are individually labeled and plotted as stars. Solid lines denote the typical demarcation between star-forming, AGN, and LINER galaxies (Kauffmann et al. 2003; Cid Fernandes et al. 2010). The gray-scale background corresponds to the density distribution of nearby (0.02 < z < 0.4) galaxies from SDSS. White contours show the KDE for each galaxy population used to model their probability density functions.

Standard image High-resolution image

Table 2. Nebular Emission-line Fluxes for FRB Hosts

FRBHα Hβ [N ii][O iii]
   6584 Å5007 Å
FRB20121102A <0.06
FRB20180301A 1.13
FRB20180916B <52.93 <74.20
FRB20180924B
FRB20181112A
FRB20190102C <0.88<0.42
FRB20190520B <0.02<0.06
FRB20190608B
FRB20190611B
FRB20190714A
FRB20191001A
FRB20200120E a 11,380198041802260
FRB20200430A
FRB20200906A
FRB20201124A 21.72
FRB20210117A
FRB20210320C
FRB20210410D <0.05
FRB20210807D <4.08
FRB20211127I
FRB20211203C
FRB20211212A <24.71 <9.51
FRB20220105A

Note. Fluxes are in units of 10−16 erg s−1 cm−2 and corrected for Galactic extinction. Not listed are the values for FRB 20220912A, for which we use the line ratios presented in Ravi et al. (2022). Limits correspond to 3σ.

a From Ho et al. (1996).

Download table as:  ASCIITypeset image

Using the standard demarcation lines between star-forming, AGN, and LINER galaxies for z ≲ 0.3 (Kauffmann et al. 2003; Cid Fernandes et al. 2010), we find that the majority of FRB hosts occupy the star-forming region, consistent with their SED classifications (Gordon et al. 2023), and are largely aligned with the locus of star-forming galaxies. Moreover, our analysis of a larger sample of hosts demonstrates that the repeating FRB hosts do not lie exclusively in the star-forming region, in contrast to previous studies (Heintz et al. 2020; Bhandari et al. 2020), and that the hosts of both repeating and nonrepeating FRBs populate all three regions of the diagram roughly according to the underlying galaxy population.

To quantify this, we compare the FRB host galaxy population, and the host populations for repeating and nonrepeating events, to the distributions of star-forming, AGN, and LINER galaxies using 2D Kolmogorov–Smirnov (K-S) tests to determine whether they emerge from the same underlying population. We first employ a kernel density estimation (KDE) to the SDSS sample to divide the parent population into individual galaxy classes. We next perform 2D K-S tests, comparing the host populations for repeating and nonrepeating events to each galaxy population, where we include all limits in the analysis. The results are summarized in Table 3.

Table 3. 2D K-S Tests for FRB Host Populations

Galaxy Type PKS (Repeaters) PKS (Nonrepeaters) PKS (All)
All Types0.011 0.168 0.058
SF0.0100.0110.003
AGN<0.001<0.001<0.001
LINER0.002<0.001<0.001
SF+LINER0.009 0.207 0.121
AGN+LINER0.002<0.001<0.001

Note. P-values derived from 2D K-S tests comparing the FRB host populations separately (repeaters and nonrepeaters) and collectively (all) to different galaxy populations according to the BPT diagram. Values in bold correspond to PKS > 0.05. SF: star forming.

Download table as:  ASCIITypeset image

Comparing the hosts of repeating and nonrepeating FRBs collectively to all galaxy types, we find PKS = 0.058, suggesting that we cannot reject the null hypothesis that FRB hosts are drawn from the same underlying galaxy population. We furthermore find that all FRB hosts are consistent with an underlying population of star-forming+LINER galaxies (PKS = 0.121), but statistically inconsistent with a population of AGN+LINER galaxies, indicating that star-forming galaxies are the dominant population.

With only three sources in the AGN locus (one of which is a limit), we reject the null hypothesis that FRB hosts are drawn from the same population as AGN. These objects include the hosts of FRB 20210807D, FRB 20210117A, and FRB 20190520B, although the latter does not exhibit detectable Hβ or [N ii] emission. Indeed, its location as an outlier on the BPT diagram, coupled with a dwarf host and the presence of a coincident PRS (Niu et al. 2022), bear a striking similarity to FRB 20121102A (Chatterjee et al. 2017). However, it is worth noting that the location of a galaxy on the BPT diagram is highly sensitive to the gas-phase metallicity (Kumari et al. 2021; Garg et al. 2022). In particular, low-metallicity galaxies exhibit a lower nitrogen abundance and a higher collisional excitation rate for [O iii], which leads to a decrease in the [N ii] flux, and an increase in the [O iii] flux, respectively (Garg et al. 2022). A third dwarf host, the host of FRB 20210117A, lies in a similar region of the diagram but hosts an apparently nonrepeating FRB and lacks a persistent radio counterpart (Bhandari et al. 2023).

Hosts for the X-ray-confirmed AGN, FRBs 20190608B and 20200120E (M81), lie within the LINER region. While the precise origin of LINER emission is not settled, studies in the last few decades have shown that AGN are present in roughly ∼ 50% of LINERs (Ho et al. 2001; Satyapal et al. 2004; Dudik et al. 2005). Here the presence of luminous X-ray emission coincident with the host nuclei for both sources points to a central AGN (as opposed to stellar ionization) as the mechanism driving the elevated levels of ionization relative to star-forming galaxies. While several other FRB hosts are located in the AGN/LINER regions, complementary X-ray constraints are lacking for these sources, with the exception of FRB 20180916B for which X-ray observations place a fairly deep limit on the presence of an AGN (see Section 2).

Our results are similar to previous findings in that we do find evidence for a subset of hosts residing in the LINER region and offset from the main locus of star-forming galaxies (Heintz et al. 2020; Bhandari et al. 2022). However, the larger sample presented here demonstrates that the hosts of repeating FRBs are not confined to the star-forming region, in contrast to these earlier studies. Indeed, our results show that the host galaxies for both repeating and nonrepeating FRBs are statistically consistent with the full underlying galaxy population, with star-forming galaxies comprising the dominant population.

4. FRBs from Accreting Compact Objects

4.1. Searching for X-Ray Counterparts

The discovery of an apparent modulation in the activity level of two luminous repeating FRB sources (FRB 20180916B, activity modulation period of 16 days; CHIME/FRB Collaboration et al. 2020; FRB 20121102A, activity modulation period of 160 days; Rajwade et al. 2020) has motivated progenitor theories invoking accreting stellar-mass compact objects (Waxman 2017; Deng et al. 2021; Sridhar et al. 2021)—most of which require super-Eddington mass transfer (in analogy with ULXs) to explain the luminous FRB population. In this scenario, a misalignment of the black hole spin axis with respect to the angular momentum axis of the accretion disk leads to Lens–Thirring precession of the disk and thus the jet axis (Middleton et al. 2019), driving periodicity in the FRB activity cycle on timescales of weeks to years. A prediction of this model is that FRBs—during their active cycle—should be accompanied by quasi-persistent X-ray counterparts.

While the discovery of periodicity is limited to only two FRBs thus far, it is useful to use our X-ray census to constrain the ULX-like progenitor scenario. In Figure 4 (left panel), we compare X-ray limits for the FRBs in our sample to the luminosity distribution for ULXs using the multimission catalog of ULXs from Walton et al. (2022). This catalog consists of 1843 ULX sources associated with 951 host galaxies compiled from XMM-Newton, Swift, and Chandra X-ray source catalogs. While the sample selection is highly nonuniform and incomplete, it represents the largest ULX catalog compiled to date.

Figure 4. Refer to the following caption and surrounding text.

Figure 4. Left: upper limits on the X-ray luminosities of the FRBs in our sample (vertical lines) compared to the luminosity distribution of ULXs from Walton et al. (2022). Not shown is FRB 20200120E, for which the X-ray luminosity limit is LX ≲ 1.3 × 1037 erg s−1. Right: X-ray limits (triangles) and detections (stars) for FRBs and FRB hosts respectively, as a function of their redshift. Shown for comparison are the fluxes of ULXs (circles), where gray lines illustrate the flux scaling with redshift. The horizontal dashed lines correspond to the sensitivity of Chandra and the planned NewAthena X-ray observatory. While existing X-ray facilities are sensitive to ULX emission below z ∼ 0.1, next-generation X-ray observatories such as NewAthena will enable searches for ULX-like counterparts from FRBs out to z ∼ 1.

Standard image High-resolution image

Current ULX surveys only identify them out to average distances of 75 Mpc (maximum distance of Mpc), while our X-ray limits on a few cosmological FRBs at z ≈ 0.0337–0.1608 are deep enough to probe the bright end of the ULX luminosity function (≈1040 erg s−1). The deepest limit comes from the closest repeating FRB 20200120E (at 3.6 Mpc), localized to a globular cluster in the nearby M81 galaxy. The X-ray luminosity at the location of the FRB is LX ≲ 1.3 × 1037 erg s−1 (Kirsten et al. 2022), corresponding to ∼0.01–0.1 LEdd for a 10 M black hole or neutron star accretor, respectively (see also Pearlman et al. 2023). This is consistent with the mass transfer rate required to power bursts from FRB 20200120E with a relatively fainter typical burst luminosity of LFRB ∼ 1037 erg s−1. We note that FRB 20200120E's localization to a globular cluster furthermore supports a connection with an old progenitor system such as accreting X-ray binaries (Kirsten et al. 2022). Conversely, the deep limit for FRB 20180916B (LX ≲ 8 × 1039 erg s−1; Scholz et al. 2020) corresponds to ∼5–50 LEdd (Sridhar et al. 2021). We note that the X-ray luminosity of the host of FRB 20190608B, which is due to a Type I AGN, is well above the ULX luminosity function as expected.

Future prospects for probing the peak of the ULX luminosity function at LX ≲ 1040 erg s−1 will be aided by the the wealth of low-DM FRBs expected to be discovered and localized by a number of existing and planned FRB facilities (including planned upgrades). In Figure 4 (right panel), we plot the X-ray fluxes of known ULXs (Walton et al. 2022) as a function of their redshift against X-ray fluxes and limits for well-localized FRBs and their hosts. Indeed, while the present-day sample of FRB hosts is dominated by sources at z ≳ 0.1, well above the population of ULXs, a number of FRB experiments promise to enhance the detection rate of precisely localized low-DM events, pushing into the low-redshift regime where ULXs are primarily detected. Such experiments (and upgrades) include ASKAP with its Commensal Real-time ASKAP Fast Transients Survey (CRAFT; Macquart et al. 2010) COherent upgrade (CRACO) system, the More TRAnsients and Pulsars (MeerTRAP; Sanidas et al. 2018) project on the MeerKAT telescope, and the Deep Synoptic Array (DSA-2000; Connor & Ravi 2023), which will also increase the detection horizon for FRBs out to z ∼ 5. The large field of view afforded by experiments like the Bustling Universe Radio Survey Telescope in Taiwan (BURSTT; Lin et al. 2022) and the Canadian Hydrogen Intensity Mapping FRB Experiment (CHIME/FRB; CHIME/FRB Collaboration et al. 2018), which will achieve subarcsecond localizations of FRBs beginning in early 2024 thanks to the addition of outrigger telescopes to the array, will significantly increase the rate of well-localized low-DM events. Above z ≳ 0.1, where the typical ULX flux scaling drops below the Chandra sensitivity threshold (see Figure 4), and as FRB experiments increase their horizon beyond the current z ∼ 1 (Ryder et al. 2023), it will be possible with next-generation X-ray observatories, such as NewAthena and the Advanced X-ray Imaging Satellite (AXIS), to perform productive searches for ULX-like counterparts.

4.2. Searching for Coincident ULX–FRB Pairs

Taking advantage of the latest ULX catalogs and the hundreds of known FRBs (albeit localized to tens of arcminutes), we search for ULXs coincident with FRBs using the multimission catalog of ULXs (Walton et al. 2022) described in Section 4.1 and the CHIME/FRB catalog (CHIME/FRB Collaboration et al. 2021). 34 The CHIME/FRB catalog includes 536 unique FRBs detected between 2018 July 25 to 2019 July 1. Among this sample, 62 events represent bursts from 18 previously reported repeating FRBs.

We first search for ULXs within 3° of a CHIME/FRB position using the FRB coordinates derived from the metadata headers (CHIME/FRB Collaboration et al. 2021). Given the non-Gaussian nature of the CHIME/FRB localization regions, we do not adopt the positional uncertainties provided in the metadata headers, and instead conservatively impose a 3σ = 3° localization uncertainty to account for the full extent of the CHIME localization regions. Our initial search returns 641 matches corresponding to 184 unique FRBs. We next manually check whether the matched ULXs actually lie within the FRB localization regions for each event using the HDF5 localization files. After applying this filter, we find a total of 28 ULXs spatially coincident with the localization regions for 17 FRBs. Interestingly, this includes one repeating FRB: FRB20190116A.

To assess whether the remaining FRB–ULX pairs are plausibly related, we plot the redshifts inferred from their DMs for the FRBs as a function of the ULX redshifts (Walton et al. 2022) in Figure 5. We estimate the FRB redshifts using the Macquart relation (Macquart et al. 2020) and two different values of DMhost: 0 pc cm−3 (shown in light blue on Figure 5) and 150 pc cm−3 (dark blue in Figure 5). 35 For the FRB DMs, we adopt the dm_exc_ne2001 values from the CHIME/FRB catalog, which account for the DMMW contribution assuming the NE2001 electron density model (Cordes & Lazio 2002). Given the low volumetric distance probed by the existing ULX catalogs, we find that the inferred redshifts for the FRBs in our cross-matched sample are inconsistent with the ULX redshifts, assuming that the true value of DMhost lies in the range 0–150 pc cm−3 and ignoring the scatter in the Macquart relation which may also impact our results (Baptista et al. 2023; James et al. 2022; James 2023; Simha et al. 2023). In particular, inhomogeneities in the baryon distribution of the CGM may bias a DM-inferred redshift lower or higher depending on the cosmic structure along an intervening sightline. We nevertheless conclude that an association between these particular ULX–FRB pairs is unlikely.

Figure 5. Refer to the following caption and surrounding text.

Figure 5. Inferred redshifts for the CHIME/FRB Catalog 1 FRBs (zFRB) with coincident ULXs in their localization regions as a function of the ULX redshifts (zULX). Dashed lines connect points with no DMhost contribution (light blue) to points with an assumed DMhost = 150 pc cm−3 (dark blue). The solid diagonal line indicates where zFRB = zULX. The horizontal dashed line at z ≈ 0.14 corresponds to the maximum observed redshift for a ULX from the complete ULX catalog.

Standard image High-resolution image

5. Conclusions

We have presented the first search for X-ray emission from a sample of FRB host galaxies. Our new Chandra observations, coupled with archival Chandra data, were used to conduct the deepest search for AGN and X-ray counterparts local to seven FRB sources. We also present an updated BPT analysis for 24 highly secure FRB hosts, presenting the nebular line properties for a large sample of FRB host galaxies for the first time. Our key results are summarized as follows.

  • 1.  
    Among the existing sample of seven well-localized FRBs with deep Chandra observations, we find that two of the host galaxies possess nuclear X-ray sources where the emission is powered by a central AGN, including the host of FRB 20200120E (M81), which hosts a known LLAGN. This corresponds to an AGN occupation fraction of ∼20% and ∼50% for repeating and nonrepeating FRBs, respectively, which may point to an elevated rate of AGN in FRB hosts (when compared to field galaxies at similar redshifts). However, these results are based on a small sample size—and uncharacterized selection effects—warranting X-ray observations of a larger sample of hosts.
  • 2.  
    Our detection of luminous (LX ≈ 5 × 1042 erg s−1) X-ray emission spatially coincident with the nucleus of the host of FRB 20190608B is consistent with its classification as a Type I Seyfert galaxy, and as such, is the only FRB host known to harbor a Type I AGN. We infer an SMBH mass of MBH ∼ 108 M and an Eddington ratio Lbol/LEdd ≈ 0.02, characteristic of geometrically thin disks in Seyfert galaxies which are less efficient at jet production. We do not detect X-ray emission at the FRB position (offset by ∼1farcs6 from the host center) and place a limit on the X-ray luminosity at the FRB location of LX ≲ 1041 erg s−1.
  • 3.  
    X-ray limits for the remainder of FRBs in our sample allow us to rule out AGN with LX ≲ 1041–1042 erg s−1, including some LLAGN. However, the majority of radio limits for these sources do not reach the regime of radio-quiet AGN.
  • 4.  
    Our updated BPT analysis of 24 FRB hosts demonstrates that FRB host galaxies trace the full underlying galaxy population but are dominated by star-forming galaxies which are largely aligned with the main star-forming locus. We find that the hosts of repeating FRBs are not located exclusively in the star-forming region, however, contrary to previous studies. Similar to earlier findings, we do find evidence for an elevated fraction of hosts occurring in the LINER region.
  • 5.  
    We demonstrate that existing X-ray limits for FRBs are sufficient for probing the bright end of the ULX luminosity function (LX ≳ 1040 erg s−1), but that the typical redshifts of FRB hosts (z ≳ 0.1) compared to the local Universe distances of ULXs renders detections at the peak of the luminosity function (LX ≲ 1040) challenging. Upcoming X-ray observatories like NewAthena and AXIS will extend the detection horizon for ULX-like counterparts to FRBs out to z ∼ 1. At the same time, upgrades to a number of FRB experiments will lend to an increased rate of low-DM FRBs which will be prime candidates for searching for ULXs at the peak of the ULX luminosity function.
  • 6.  
    Performing a cross-matching analysis between the CHIME/FRB catalog and the largest catalog of ULX sources compiled to date, we find a total of 28 ULXs spatially coincident with the localization regions for 17 FRBs. However, the DM-inferred redshifts for these FRBs imply distances ∼10–103 times larger than the distances to the ULXs, and hence we consider any associations unlikely.

Continued X-ray observations of FRB host galaxies will critically probe the occupation fraction of AGN. Understanding the AGN fraction in FRB host galaxies, as well as pinpointing the dominant source of ionizing radiation prevalent in hosts that are optically classified as LINERS, will provide new insight into the stellar populations that should ultimately drive FRB production. Searches for X-ray counterparts will benefit most from dedicated follow up of the most nearby events, where such observations can be most constraining.

Acknowledgments

We thank Kari A. Frank, Paul Green, Alice P. Curtin, Kristina Nyland, and Victoria Kaspi for helpful discussions. T.E. is supported by NASA through the NASA Hubble Fellowship grant HST-HF2-51504.001-A awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS5-26555. W.F. gratefully acknowledges support by the National Science Foundation under CAREER grant No. AST-2047919, the David and Lucile Packard Foundation, the Alfred P. Sloan Foundation, and the Research Corporation for Science Advancement through Cottrell Scholar Award #28284. N.S. acknowledges the support from NASA (grant No. 80NSSC22K0332), NASA FINESST (grant No. 80NSSC22K1597), and Columbia University Dean's fellowship. S.B. is supported by a Dutch Research Council (NWO) Veni Fellowship (VI.Veni.212.058). Y.D. is supported by the National Science Foundation Graduate Research Fellowship under grant No. DGE-1842165. A.T.D. and R.M.S. acknowledge support through Australian Research Council ARC DP DP220102305. R.M.S. acknowledges support through Australian Research Council Future Fellowship FT190100155. J.X.P., A.C.G., Y.D., W.F., T.E., N.T., and C.D.K. acknowledge support from NSF grants AST-1911140, AST-1910471, and AST-2206490 as members of the Fast and Fortunate for FRB Follow-up team. A.B.P. is a Banting Fellow, a McGill Space Institute (MSI) Fellow, and a Fonds de Recherche du Quebec—Nature et Technologies (FRQNT) postdoctoral fellow.

Facilities: the VLA observations presented here were obtained as part of programs VLA/20A-157 - , PI: Bhandari and VLA/22B-115 - , PI: Eftekhari. The VLA is operated by the National Radio Astronomy Observatory - , a facility of the National Science Foundation operated under cooperative agreement by Associated Universities - , Inc. The Chandra ACIS data for FRB20200120E used in this paper can be found in the Chandra Data Archive under the following Observation IDs for (1) the FRB position: 9540 - , PI: Jenkins; 4752 - , PI: Miller; 13728 and 14471 - , PI: Grise; 18051 - , 18052 - , 18817 - , 19685 - , 19992 and 19993 - , PI: Swartz; and (2) the host center: 735 - , 18047 - , 18048 - , 18049 - , 18050 - , 18051 - , 18052 - , 18053 - , 18054 - , 18817 - , 18875 - , 19685 - , 19688 - , 19981 - , 19982 - , 19991 - , 19992 and 19993 - , PI: Swartz; 5600 - , 5601 - , 6174 - , 6346 - , and 6347 - , PI: Canizares; 5935 - , 5936 - , 5937 - , 5938 - , 5939 - , 5940 - , 5941 - , 5942 - , 5943 - , 5944 - , 5945 - , 5946 - , 5947 - , 5948 and 5949 - , PI: Pooley; 4752 - , 6892 - , 6893 - , 6894 - , 6895 - , 6896 - , 6897 - , 6898 - , 6899 - , 6900 and 6901 - , PI: Miller; 9122 - , PI: Immler; 9540 - , 9805 - , PI: Jonker; 12301 - , PI: Dwarkadas; 13728 - , 14471 - , PI: Grise; 20624 - , PI: Barth; 21384 - , PI: Irwin; 26421 - , PI: Brightman - .

Software: astropy (Astropy Collaboration et al. 2022), CASA (McMullin et al. 2007), CIAO (v4.13) (Fruscione et al. 2006), matplotlib (Hunter 2007), numpy (Harris et al. 2020), pandas (Reback 2022), Prospector (Johnson et al. 2021), SAOImageDS9 (Joye & Mandel 2003), SExtractor (Bertin & Arnouts 1996), scipy (Virtanen et al. 2020), superchandra (Kilpatrick 2023), and xspec (Arnaud1996).

Footnotes

  • 30  

    The z = 0 source corresponds to FRB 20200120E located in the nearby M81 galaxy at d = 3.6 Mpc (Kirsten et al. 2022).

  • 31  

    Details of the observations included in each merged event file are given in the Facilities statement at the end of the paper.

  • 32  

    LLRGs are canonically unresolved, nonthermal radio cores in galaxies with optical classifications spanning Type 1 Seyferts, Type 2 Seyferts, and LINER galaxies.

  • 33  

    Line fluxes for FRBs 20200120E and 20220912A were obtained using a different method than the one we use here.

  • 34  
  • 35  

    This is comparable to the median value of estimated in James et al. (2022).

10.3847/1538-4357/acf843
undefined