Draft version November 13, 2020 Typeset using LATEX twocolumn style in AASTeX63 Fundamental Reference AGN Monitoring Experiment (FRAMEx) I: Jumping Out of the Plane with the VLBA Travis C. Fischer ,1, 2 Nathan J. Secrest ,2 Megan C. Johnson ,2 Bryan N. Dorland ,2 Phillip J. Cigan ,2 Luis C. Fernandez ,3 Lucas R. Hunt ,2 Michael Koss ,4 Henrique R. Schmitt ,5 and Norbert Zacharias 2 1AURA for ESA, Space Telescope Science Institute, Baltimore, MD, USA, 3700 San Martin Drive, Baltimore, MD 21218, USA 2U.S. Naval Observatory, 3450 Massachusetts Ave NW, Washington, DC 20392-5420, USA 3Department of Physics and Astronomy, George Mason University, MS3F3, 4400 University Drive, Fairfax, VA 22030, USA 4Eureka Scientific, 2452 Delmer Street Suite 100, Oakland, CA 94602-3017, USA 5Naval Research Laboratory, Washington, DC 20375, USA ABSTRACT We present the first results from the Fundamental Reference AGN Monitoring Experiment (FRAMEx), an observational campaign dedicated to understanding the physical processes that affect the apparent positions and morphologies of AGNs. In this work, we obtained simultaneous Swift X-ray Telescope (XRT) and Very Long Baseline Array (VLBA) radio observations for a snapshot campaign of 25 local AGNs that form a volume-complete sample with hard X-ray (14{195 keV) luminosities above 1042 erg s−1, out to a distance of 40 Mpc. Despite achieving an observation depth of ∼ 20 µJy, we find that 16 of 25 AGNs in our sample are not detected with the VLBA on milli-arcsecond (sub-parsec) scales, and the corresponding core radio luminosity upper limits are systematically below predictions from the Fundamental Plane of black hole activity. Using archival Jansky Very Large Array (VLA) radio measurements, our sample jumps back onto the Fundamental Plane, suggesting that extended radio emission is responsible for the apparent correlation between radio emission, X-ray emission, and black hole mass. We suggest that this discrepancy is likely due to extra-nuclear radio emission pro- duced via interactions between the AGN and host environment. We compare VLBA observations of AGNs to VLA observations of nearby Galactic black holes (GBHs) and we find a mass-independent −6 correlation between radio and X-ray luminosities of black holes of L6 cm/L2−10 keV ∼ 10 , in line with predictions for coronal emission, but allowing for the possibility of truly radio silent AGNs. Keywords: Radio astrometry (1337), Active galaxies (17), Radio active galactic nuclei (2134), X-ray active galactic nuclei (2035) 1. INTRODUCTION back" from the AGN into the surrounding interstellar Over the past several decades, supermassive black medium via either the deposition of mechanical energy holes (SMBHs), which manifest as active galactic nu- from a collimated jet, or radiative energy in the form clei (AGNs), have gained profound importance to sev- of AGN winds (for an essential review of this topic, see Kormendy & Ho 2013). Similarly, the discovery of in- arXiv:2011.06570v1 [astro-ph.GA] 12 Nov 2020 eral fields of astrophysics. With the discovery of the M − σ relation (Ferrarese & Merritt 2000; Gebhardt creasingly high-redshift quasars, such as the recent an- BH 8 et al. 2000) at the turn of the millennium, it became nouncement of one powered by an 8×10 M SMBH at a clear that SMBHs and their host galaxies co-evolve in redshift of 7.54 (Ba~nadoset al. 2018) when the universe some fashion, likely involving the mutual feeding of star was only a few hundred million years old, has raised seri- formation and SMBH accretion, combined with \feed- ous questions about the formation mechanisms of black holes in the early universe (e.g., Volonteri 2010), and the distribution of quasars across the sky may even provide Corresponding author: Travis Fischer independent tests of our cosmological assumptions (e.g., tfi[email protected] Singal 2019, and references therein). 2 Fischer et al. The nearly uniform distribution of quasars across the redshifts of 1{2 correspond to a minimum distance sep- sky, combined with their negligible proper motions and aration of ∼ 8 parsec at the milli-arcsecond (mas) reso- significant luminosity at almost every wavelength, have lution of very long baseline interferometry (VLBI) net- additionally made them ideal objects with which to re- works such as the Very Long Baseline Array (VLBA). alize a quasi-inertial celestial reference frame, with the The accretion regions of AGNs, which are typically con- first quasar-based realization of the International Ce- fined to within ∼ 1 parsec, are therefore much smaller lestial Reference System (ICRS), the International Ce- in angular extent than what is typically observable, lim- lestial Reference Frame (ICRF; Ma et al. 1998), being iting research to AGN emission on large scales. It is adopted at the same time as the discovery of SMBH- therefore worthwhile to study more nearby AGNs, where galaxy scaling relations. Successive refinements of the mas-scale observations taken with facilities such as the ICRF have improved both the density of these refer- VLBA can resolve features subtending only a few tenths ence sources and their respective astrometric precision of a parsec, providing a direct probe into the accretion (ICRF2: Fey et al. 2015, and ICRF3: Charlot et al., in physics of AGNs. prep.). The ICRF has historically been a radio refer- In this work, we present results from an initial snap- ence frame with quasar positions defined at 8.6 GHz. shot survey of a volume-complete sample of 25 nearby However, for the first time, the ICRF3 is now multi- (< 40 Mpc) AGNs, in which we obtained simultaneous wavelength with positions defined at 24 GHz (1.2 cm) observations with the VLBA and the Neil Gehrels Swift and 32 GHz (0.9 cm) (Charlot et al., in prep. ). If the Observatory X-ray Telescope (XRT), as part of a long- ICRF is to be viewed as the fundamental reference frame term program to monitor the behavior of AGNs at sub- for all astronomical positions and proper motions, there parsec scales. The simultaneity of these data mitigates should be no significant disagreement between the ra- the uncertainty introduced by variability that is inher- dio positions of ICRF sources and their counterparts at ent to studying AGNs on small physical scales, allowing other wavelengths. This has unfortunately not been the a more robust assessment of the relationship between case, with evidence of the existence of significant optical- emission mechanisms at different wavelengths and ener- radio positional offsets being published soon after the gies. Our sample was selected from the Swift 105-month release of ICRF1 (da Silva Neto et al. 2002), ICRF2 Burst Alert Telescope (BAT) catalog (Oh et al. 2018), (Assafin et al. 2013; Orosz & Frey 2013; Zacharias & enabling obscuration-independent selection of a statisti- Zacharias 2014) and, along with the second data release cally representative sample of nearby AGNs above some of Gaia, ICRF3 (Makarov et al. 2017). intrinsic luminosity limit. Earlier, pre-Gaia work hinted at an astrophysical The primary focus of this work is the relationship be- nature to optical-radio offsets (Zacharias & Zacharias tween the radio luminosity, X-ray luminosity, and black 2014), but the unprecedented astrometric precision of hole mass of AGNs, i.e. the Fundamental Plane of black Gaia has put the astrophysical nature of these offsets hole activity (e.g., Merloni et al. 2003, hereafter the FP), on firm empirical ground (e.g., Harris et al. 2016; Ko- which potentially serves as a parameter space unifying valev et al. 2017; Petrov et al. 2019; Plavin et al. 2019). black hole accretion across all black hole masses, includ- It is therefore imperative to understand how the physi- ing Galactic black holes. We additionally explore the re- cal phenomena associated with AGN activity, such as lationship between the core radio and X-ray luminosities jets, intrinsic variability, variable obscuration, binary of these AGNs more broadly, considering the potential SMBHs, host galaxy morphology, etc. affect the appar- bimodality between radio-loud and radio-quiet AGNs, ent positions of AGNs and quasars across the electro- as well as the prevalence of potentially truly radio-silent magnetic spectrum. AGNs. Towards this goal, the U.S. Naval Observatory and its collaborators from partner institutions have begun 2. METHODOLOGY the Fundamental Reference AGN Monitoring Experi- 2.1. Sample Selection ment (FRAMEx), a research framework that includes Our AGN sample was selected from the Swift BAT the physical characterization and monitoring of AGN 105-month catalog (Oh et al. 2018). Because this cata- phenomena that may affect the apparent positions of log has a uniform flux limit (∼ 8 × 10−12 erg cm−2 s−1) AGNs at different wavelengths (see Dorland et al.(2020) across the sky, and hard X-rays (14{195 keV) reliably for a more detailed overview of the FRAMEx collab- pick out bona-fide AGNs, even in the presence of heavy oration and discussion of the relevant angular scales). absorption (N 1024 cm−2), constructing a volume- While direct studies of the quasars that comprise the H . complete sample of AGNs above some luminosity thresh- ICRF are an essential part of FRAMEx, their typical old is straightforward. Additionally, AGNs from the FRAMEx I: Volume-Complete Radio and X-ray Snapshot 3 BAT catalog have extensive multi-wavelength cover- age, and have the benefit of intense recent study (e.g., Trakhtenbrot et al. 2017; Shimizu et al. 2018; Powell et al. 2018; B¨aret al. 2019). We choose AGNs above a hard X-ray luminosity of 1042 erg s−1, which corresponds to a distance limit of ∼ 40 Mpc for a volume-complete sample using a flat −1 −1 ΛCDM cosmology with H0 = 70 km s Mpc and ΩM = 0:3. This luminosity threshold was chosen to sam- ple both moderate and high luminosity AGNs, and to provide a sample of AGNs close enough to give exquisite physical resolution on the AGN.