Search

The cosmic merger history of supermassive black hole binaries (SMBHBs) is expected to produce a low-frequency gravitational wave background (GWB). Here we investigate how signs of the discrete nature of this GWB can manifest in pulsar timing arrays through excursions from, and breaks in, the expected $f_{\mathrm{GW}}^{-2/3}$ power-law of the GWB strain spectrum. To do this, we create a semi-analytic SMBHB population model, fit to NANOGrav's 15 yr GWB amplitude, and with 1,000 realizations we study the populations' characteristic strain and residual spectra. Comparing our models to the NANOGrav 15 yr spectrum, we find two interesting excursions from the power-law. The first, at $2 \; \mathrm{nHz}$, is below our GWB realizations with $p$-value significance $p = 0.05$ to $0.06$ ($\approx 1.8 \sigma - 1.9 \sigma$). The second, at $16 \; \mathrm{nHz}$, is above our GWB realizations with $p = 0.04$ to $0.15$ ($\approx 1.4 \sigma - 2.1 \sigma$). We explore the properties of a loud SMBHB which could cause such an excursion. Our simulations also show that the expected number of SMBHBs decreases by three orders of magnitude, from $\sim 10^6$ to $\sim 10^3$, between $2\; \mathrm{nHz}$ and $20 \; \mathrm{nHz}$. This causes a break in the strain spectrum as the stochasticity of the background breaks down at $26^{+28}_{-19} \; \mathrm{nHz}$, consistent with predictions pre-dating GWB measurements. The diminished GWB signal from SMBHBs at frequencies above the $26$~nHz break opens a window for PTAs to detect continuous GWs from individual SMBHBs or GWs from the early universe., Comment: 10 pages, 8 figures, 1 appendix, submitted to ApJ

Recently we found compelling evidence for a gravitational wave background with Hellings and Downs (HD) correlations in our 15-year data set. These correlations describe gravitational waves as predicted by general relativity, which has two transverse polarization modes. However, more general metric theories of gravity can have additional polarization modes which produce different interpulsar correlations. In this work we search the NANOGrav 15-year data set for evidence of a gravitational wave background with quadrupolar Hellings and Downs (HD) and Scalar Transverse (ST) correlations. We find that HD correlations are the best fit to the data, and no significant evidence in favor of ST correlations. While Bayes factors show strong evidence for a correlated signal, the data does not strongly prefer either correlation signature, with Bayes factors $\sim 2$ when comparing HD to ST correlations, and $\sim 1$ for HD plus ST correlations to HD correlations alone. However, when modeled alongside HD correlations, the amplitude and spectral index posteriors for ST correlations are uninformative, with the HD process accounting for the vast majority of the total signal. Using the optimal statistic, a frequentist technique that focuses on the pulsar-pair cross-correlations, we find median signal-to-noise-ratios of 5.0 for HD and 4.6 for ST correlations when fit for separately, and median signal-to-noise-ratios of 3.5 for HD and 3.0 for ST correlations when fit for simultaneously. While the signal-to-noise-ratios for each of the correlations are comparable, the estimated amplitude and spectral index for HD are a significantly better fit to the total signal, in agreement with our Bayesian analysis., Comment: 11 pages, 5 figures

Pulsar timing arrays (PTAs) are galactic-scale gravitational wave detectors. Each individual arm, composed of a millisecond pulsar, a radio telescope, and a kiloparsecs-long path, differs in its properties but, in aggregate, can be used to extract low-frequency gravitational wave (GW) signals. We present a noise and sensitivity analysis to accompany the NANOGrav 15-year data release and associated papers, along with an in-depth introduction to PTA noise models. As a first step in our analysis, we characterize each individual pulsar data set with three types of white noise parameters and two red noise parameters. These parameters, along with the timing model and, particularly, a piecewise-constant model for the time-variable dispersion measure, determine the sensitivity curve over the low-frequency GW band we are searching. We tabulate information for all of the pulsars in this data release and present some representative sensitivity curves. We then combine the individual pulsar sensitivities using a signal-to-noise-ratio statistic to calculate the global sensitivity of the PTA to a stochastic background of GWs, obtaining a minimum noise characteristic strain of $7\times 10^{-15}$ at 5 nHz. A power law-integrated analysis shows rough agreement with the amplitudes recovered in NANOGrav's 15-year GW background analysis. While our phenomenological noise model does not model all known physical effects explicitly, it provides an accurate characterization of the noise in the data while preserving sensitivity to multiple classes of GW signals., Comment: 67 pages, 73 figures, 3 tables; published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

Mrk 1018 is a nearby changing-look AGN that has oscillated between spectral Type 1.9 and Type 1 over a period of 40 years. Recently, a recoiling supermassive black hole (rSMBH) scenario has been proposed to explain the spectral and flux variability observed in this AGN. Detections of rSMBHs are important for understanding the processes by which SMBH binaries merge and how rSMBHs influence their galactic environment through feedback mechanisms. However, conclusive identification of any rSMBHs has remained elusive to date. In this paper, we present an analysis of 6.5 years of multi-frequency Very Long Baseline Array (VLBA) monitoring of Mrk 1018. We find that the radio emission is compact down to 2.4 pc, and displays flux density and spectral variability over the length of our campaign, typical of a flat spectrum radio core. We observe proper motion in RA of the radio core at -36.4 $\pm$ 8.6 $\mu$as yr$^{-1}$ (4.2$\sigma$), or $0.10c \pm 0.02c$ at the redshift of Mrk 1018. No significant proper motion is found in DEC (31.3 $\pm$ 25.1 $\mu$as yr$^{-1}$). We discuss possible physical mechanisms driving the proper motion, including a rSMBH. We conclude that the apparent velocity we observe of the VLBI radio core is too high to reconcile with theoretical predictions of rSMBH velocities and that the proper motion is most likely dominated by an unresolved, outflowing jet component. Future observations may yet reveal the true nature of Mrk 1018. However, our observations are not able to confirm it as a true rSMBH., Comment: 17 pages, 7 figures, 3 tables; accepted for publication in ApJ

Pulsar timing array collaborations, such as the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), are seeking to detect nanohertz gravitational waves emitted by supermassive black hole binaries formed in the aftermath of galaxy mergers. We have searched for continuous waves from individual circular supermassive black hole binaries using the NANOGrav's recent 12.5-year data set. We created new methods to accurately model the uncertainties on pulsar distances in our analysis, and we implemented new techniques to account for a common red noise process in pulsar timing array data sets while searching for deterministic gravitational wave signals, including continuous waves. As we found no evidence for continuous waves in our data, we placed 95\% upper limits on the strain amplitude of continuous waves emitted by these sources. At our most sensitive frequency of 7.65 nanohertz, we placed a sky-averaged limit of $h_0 < $ $(6.82 \pm 0.35) \times 10^{-15}$, and $h_0 <$ $(2.66 \pm 0.15) \times 10^{-15}$ in our most sensitive sky location. Finally, we placed a multi-messenger limit of $\mathcal{M} <$ $(1.41 \pm 0.02) \times 10^9 M_\odot$ on the chirp mass of the supermassive black hole binary candidate 3C~66B., Comment: 22 pages, 12 figures. Accepted by ApJL

We search NANOGrav's 12.5-year data set for evidence of a gravitational wave background (GWB) with all the spatial correlations allowed by general metric theories of gravity. We find no substantial evidence in favor of the existence of such correlations in our data. We find that scalar-transverse (ST) correlations yield signal-to-noise ratios and Bayes factors that are higher than quadrupolar (tensor transverse, TT) correlations. Specifically, we find ST correlations with a signal-to-noise ratio of 2.8 that are preferred over TT correlations (Hellings and Downs correlations) with Bayesian odds of about 20:1. However, the significance of ST correlations is reduced dramatically when we include modeling of the Solar System ephemeris systematics and/or remove pulsar J0030$+$0451 entirely from consideration. Even taking the nominal signal-to-noise ratios at face value, analyses of simulated data sets show that such values are not extremely unlikely to be observed in cases where only the usual TT modes are present in the GWB. In the absence of a detection of any polarization mode of gravity, we place upper limits on their amplitudes for a spectral index of $\gamma = 5$ and a reference frequency of $f_\text{yr} = 1 \text{yr}^{-1}$. Among the upper limits for eight general families of metric theories of gravity, we find the values of $A^{95\%}_{TT} = (9.7 \pm 0.4)\times 10^{-16}$ and $A^{95\%}_{ST} = (1.4 \pm 0.03)\times 10^{-15}$ for the family of metric spacetime theories that contain both TT and ST modes., Comment: 24 pages, 18 figures, 3 appendices. Please send any comments/questions to Nima Laal (laaln@oregonstate.edu)

We present the localization and host galaxies of one repeating and two apparently non-repeating Fast Radio Bursts. FRB20180301A was detected and localized with the Karl G. Jansky Very Large Array to a star-forming galaxy at $z=0.3304$. FRB20191228A, and FRB20200906A were detected and localized by the Australian Square Kilometre Array Pathfinder to host galaxies at $z=0.2430$ and $z=0.3688$, respectively. We combine these with 13 other well-localized FRBs in the literature, and analyze the host galaxy properties. We find no significant differences in the host properties of repeating and apparently non-repeating FRBs. FRB hosts are moderately star-forming, with masses slightly offset from the star-forming main-sequence. Star formation and low-ionization nuclear emission-line region (LINER) emission are major sources of ionization in FRB host galaxies, with the former dominant in repeating FRB hosts. FRB hosts do not track stellar mass and star formation as seen in field galaxies (more than 95% confidence). FRBs are rare in massive red galaxies, suggesting that progenitor formation channels are not solely dominated by delayed channels which lag star formation by Gigayears. The global properties of FRB hosts are indistinguishable from core-collapse supernovae (CCSNe) and short gamma-ray bursts (SGRBs) hosts, and the spatial offset (from galaxy centers) of FRBs is mostly inconsistent with that of the Galactic neutron star population (95% confidence). The spatial offsets of FRBs (normalized to the galaxy effective radius) also differ from those of globular clusters (GCs) in late- and early-type galaxies with 95% confidence., Comment: 27 pages, 12 figures, 8 tables. Accepted for publication in the Astrophysical Journal

The physical properties of fast radio burst (FRB) host galaxies provide important clues towards the nature of FRB sources. The 16 FRB hosts identified thus far span three orders of magnitude in mass and specific star-formation rate, implicating a ubiquitously occurring progenitor object. FRBs localised with ~arcsecond accuracy also enable effective searches for associated multi-wavelength and multi-timescale counterparts, such as the persistent radio source associated with FRB 20121102A. Here we present a localisation of the repeating source FRB 20201124A, and its association with a host galaxy (SDSS J050803.48+260338.0, z=0.098) and persistent radio source. The galaxy is massive ($\sim3\times10^{10} M_{\odot}$), star-forming (few solar masses per year), and dusty. Very Large Array and Very Long Baseline Array observations of the persistent radio source measure a luminosity of $1.2\times10^{29}$ erg s$^{-1}$ Hz$^{-1}$, and show that is extended on scales $\gtrsim50$ mas. We associate this radio emission with the ongoing star-formation activity in SDSS J050803.48+260338.0. Deeper, more detailed observations are required to better utilise the milliarcsecond-scale localisation of FRB 20201124A reported from the European VLBI Network, and determine the origin of the large dispersion measure ($150-220$ pc cm$^{-3}$) contributed by the host. SDSS J050803.48+260338.0 is an order of magnitude more massive than any galaxy or stellar system previously associated with a repeating FRB source, but is comparable to the hosts of so far non-repeating FRBs, further building the link between the two apparent populations., Comment: 9 pages, 6 figures, 1 table, submitted to MNRAS

In several previous studies, quasars exhibiting broad emission lines with >1000 km/s velocity offsets with respect to the host galaxy rest frame have been discovered. One leading hypothesis for the origin of these velocity-offset broad lines is the dynamics of a binary supermassive black hole (SMBH). We present high-resolution radio imaging of 34 quasars showing these velocity-offset broad lines with the Very Long Baseline Array (VLBA), aimed at finding evidence for the putative binary SMBHs (such as dual radio cores), and testing the competing physical models. We detect exactly half of the target sample from our VLBA imaging, after implementing a 5 detection limit. While we do not resolve double radio sources in any of the targets, we obtain limits on the instantaneous projected separations of a radio-emitting binary for all of the detected sources under the assumption that a binary still exists within our VLBA angular resolution limits. We also assess the likelihood that a radio-emitting companion SMBH exists outside of our angular resolution limits, but its radio luminosity is too weak to produce a detectable signal in the VLBA data. Additionally, we compare the precise sky positions afforded by these data to optical positions from both the SDSS and Gaia DR2 source catalogs. We find projected radio/optical separations on the order of 10 pc for three quasars. Finally, we explore how future multi-wavelength campaigns with optical, radio, and X-ray observatories can help discriminate further between the competing physical models., Comment: Accepted for publishing by ApJ

Supermassive black hole binaries (SMBHBs) should form frequently in galactic nuclei as a result of galaxy mergers. At sub-parsec separations, binaries become strong sources of low-frequency gravitational waves (GWs), targeted by Pulsar Timing Arrays (PTAs). We used recent upper limits on continuous GWs from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11yr dataset to place constraints on putative SMBHBs in nearby massive galaxies. We compiled a comprehensive catalog of ~44,000 galaxies in the local universe (up to redshift ~0.05) and populated them with hypothetical binaries, assuming that the total mass of the binary is equal to the SMBH mass derived from global scaling relations. Assuming circular equal-mass binaries emitting at NANOGrav's most sensitive frequency of 8nHz, we found that 216 galaxies are within NANOGrav's sensitivity volume. We ranked the potential SMBHBs based on GW detectability by calculating the total signal-to-noise ratio (S/N) such binaries would induce within the NANOGrav array. We placed constraints on the chirp mass and mass ratio of the 216 hypothetical binaries. For 19 galaxies, only very unequal-mass binaries are allowed, with the mass of the secondary less than 10 percent that of the primary, roughly comparable to constraints on a SMBHB in the Milky Way. Additionally, we were able to exclude binaries delivered by major mergers (mass ratio of at least 1/4) for several of these galaxies. We also derived the first limit on the density of binaries delivered by major mergers purely based on GW data., Comment: Submitted to ApJ. Send comments to Maria Charisi (maria.charisi@nanograv.org)

Fast Radio Burst FRB 20180916B in its host galaxy SDSS J015800.28+654253.0 at 149 Mpc is by far the closest-known FRB with a robust host galaxy association. The source also exhibits a 16.35-day period in its bursting. Here we present optical and infrared imaging as well as integral field spectroscopy observations of FRB 20180916B with the WFC3 camera on the Hubble Space Telescope and the MEGARA spectrograph on the 10.4-m Gran Telescopio Canarias. The 60-90 milliarcsecond (mas) resolution of the Hubble imaging, along with the previous 2.3-mas localization of FRB 20180916B, allow us to probe its environment with a 30-60 pc resolution. We constrain any point-like star-formation or HII region at the location of FRB 20180916B to have an H$\alpha$ luminosity $L_\mathrm{H\alpha} \lesssim 10^{37}\,\mathrm{erg\,s^{-1}}$ and, correspondingly, constrain the local star-formation rate to be $\lesssim10^{-4}\,\mathrm{M_\odot\,yr^{-1}}$. The constraint on H$\alpha$ suggests that possible stellar companions to FRB 20180916B should be of a cooler, less massive spectral type than O6V. FRB 20180916B is 250 pc away (in projected distance) from the brightest pixel of the nearest young stellar clump, which is $\sim380$\,pc in size (full-width at half maximum). With the typical projected velocities of pulsars, magnetars, or neutron stars in binaries (60-750 km s$^{-1}$), FRB 20180916B would need 800 kyr to 7 Myr to traverse the observed distance from its presumed birth site. This timescale is inconsistent with the active ages of magnetars ($\lesssim10$ kyr). Rather, the inferred age and observed separation are compatible with the ages of high-mass X-ray binaries and gamma-ray binaries, and their separations from the nearest OB associations., Comment: Updated version: Updated Figure 2. 16 pages, 4 figures, 1 table. Published in ApJ Letters. Comments welcome

We use Chandra X-ray observations to look for evidence of a recoiling black hole from the brightest cluster galaxy in Abell 2261 (A2261-BCG). A2261-BCG is a strong candidate for a recoiling black hole because of its large, flat stellar core, revealed by Hubble Space Telescope imaging observations. We took 100-ksec observations with Chandra and combined it with 35 ksec of archival observations to look for low-level accretion onto a black hole of expected mass $M\sim10^{10}\ M_{\scriptscriptstyle \odot}$ that could possibly be located in one of four off-center stellar knots near the galaxy's center or else in the optical center of the galaxy or in the location of radio emission. We found no X-ray emission arising from a point source in excess of the cluster gas and can place limits on the accretion of any black hole in the central region to a 2-7 keV flux below $4.3 \times 10^{-16}\ \mathrm{erg\ s^{-1}\ cm^{-2}}$, corresponding to a bolometric Eddington fraction of about $10^{-6}$. Thus there is either no $10^{10}\ M_{\scriptscriptstyle \odot}$ black hole in the core of A2261-BCG, or it is accreting at a low level. We also discuss the morphology of the X-ray emitting gas in the cluster and how its asymmetry is consistent with a large dynamic event., Comment: Accepted to AAS Journals. 12 pages of main text and figures, plus references and appendix tables and figures

The NANOGrav Collaboration reported strong Bayesian evidence for a common-spectrum stochastic process in its 12.5-yr pulsar timing array dataset, with median characteristic strain amplitude at periods of a year of $A_{\rm yr} = 1.92^{+0.75}_{-0.55} \times 10^{-15}$. However, evidence for the quadrupolar Hellings \& Downs interpulsar correlations, which are characteristic of gravitational wave signals, was not yet significant. We emulate and extend the NANOGrav dataset, injecting a wide range of stochastic gravitational wave background (GWB) signals that encompass a variety of amplitudes and spectral shapes, and quantify three key milestones: (I) Given the amplitude measured in the 12.5 yr analysis and assuming this signal is a GWB, we expect to accumulate robust evidence of an interpulsar-correlated GWB signal with 15--17 yrs of data, i.e., an additional 2--5 yrs from the 12.5 yr dataset; (II) At the initial detection, we expect a fractional uncertainty of $40\%$ on the power-law strain spectrum slope, which is sufficient to distinguish a GWB of supermassive black-hole binary origin from some models predicting more exotic origins;(III) Similarly, the measured GWB amplitude will have an uncertainty of $44\%$ upon initial detection, allowing us to arbitrate between some population models of supermassive black-hole binaries. In addition, power-law models are distinguishable from those having low-frequency spectral turnovers once 20~yrs of data are reached. Even though our study is based on the NANOGrav data, we also derive relations that allow for a generalization to other pulsar-timing array datasets. Most notably, by combining the data of individual arrays into the International Pulsar Timing Array, all of these milestones can be reached significantly earlier., Comment: 15 pages, 7 figures

We search for an isotropic stochastic gravitational-wave background (GWB) in the $12.5$-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. Our analysis finds strong evidence of a stochastic process, modeled as a power-law, with common amplitude and spectral slope across pulsars. The Bayesian posterior of the amplitude for an $f^{-2/3}$ power-law spectrum, expressed as the characteristic GW strain, has median $1.92 \times 10^{-15}$ and $5\%$--$95\%$ quantiles of $1.37$--$2.67 \times 10^{-15}$ at a reference frequency of $f_\mathrm{yr} = 1 ~\mathrm{yr}^{-1}$. The Bayes factor in favor of the common-spectrum process versus independent red-noise processes in each pulsar exceeds $10,000$. However, we find no statistically significant evidence that this process has quadrupolar spatial correlations, which we would consider necessary to claim a GWB detection consistent with general relativity. We find that the process has neither monopolar nor dipolar correlations, which may arise from, for example, reference clock or solar system ephemeris systematics, respectively. The amplitude posterior has significant support above previously reported upper limits; we explain this in terms of the Bayesian priors assumed for intrinsic pulsar red noise. We examine potential implications for the supermassive black hole binary population under the hypothesis that the signal is indeed astrophysical in nature., Comment: 25 pages, 14 figures, 5 tables, 3 appendices. Published in The Astrophysical Journal Letters. Please send any comments/questions to Joseph Simon (joe.simon@nanograv.org). Jupyter notebook tutorials and some MCMC chain files are available at https://github.com/nanograv/12p5yr_stochastic_analysis

The "Space VLBI 2020: Science and Technology Futures" meeting was the second in The Future of High-Resolution Radio Interferometry in Space series. The first meeting (2018 September 5--6; Noordwijk, the Netherlands) focused on the full range of science applications possible for very long baseline interferometry (VLBI) with space-based antennas. Accordingly, the observing frequencies (wavelengths) considered ranged from below 1~MHz (> 300 m) to above 300~GHz (< 1 mm). For this second meeting, the focus was narrowed to mission concepts and the supporting technologies to enable the highest angular resolution observations at frequencies of 30~GHz and higher (< 1 cm). This narrowing of focus was driven by both scientific and technical considerations. First, results from the RadioAstron mission and the Event Horizon Telescope (EHT) have generated considerable excitement for studying the inner portions of black hole (BH) accretion disks and jets and testing elements of the General Theory of Relativity (GR). Second, the technologies and requirements involved in space-based VLBI differ considerably between 100~MHz and 100~GHz; a related consideration is that there are a number of existing instruments or mission concepts for frequencies of approximately 100~MHz and below, while it has been some time since attention has been devoted to space VLBI at frequencies above 10~GHz. This conference summary attempts to capture elements of presentations and discussions that occurred., Comment: Nine pages + references and conference schedule appendix; conference Web site http://www.cvent.com/events/space-vlbi-2020/event-summary-c7b30de0068e454591a66e5639f86668.aspx

When galaxies merge, the supermassive black holes in their centers may form binaries and, during the process of merger, emit low-frequency gravitational radiation in the process. In this paper we consider the galaxy 3C66B, which was used as the target of the first multi-messenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of the source of the source, it has been theorized to contain a supermassive black hole binary. Its apparent 1.05-year orbital period would place the gravitational wave emission directly in the pulsar timing band. Since the first pulsar timing array study of 3C66B, revised models of the source have been published, and timing array sensitivities and techniques have improved dramatically. With these advances, we further constrain the chirp mass of the potential supermassive black hole binary in 3C66B to less than $(1.65\pm0.02) \times 10^9~{M_\odot}$ using data from the NANOGrav 11-year data set. This upper limit provides a factor of 1.6 improvement over previous limits, and a factor of 4.3 over the first search done. Nevertheless, the most recent orbital model for the source is still consistent with our limit from pulsar timing array data. In addition, we are able to quantify the improvement made by the inclusion of source properties gleaned from electromagnetic data to `blind' pulsar timing array searches. With these methods, it is apparent that it is not necessary to obtain exact a priori knowledge of the period of a binary to gain meaningful astrophysical inferences., Comment: 14 pages, 6 figures. Accepted by ApJ

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

We use a semi-analytic model to explore the potential impact of a brief and violent period of radio-loud accretion onto black-holes (The Radio Scream) during the Cosmic Dawn on the HI hyperfine 21 cm signal. We find that radio emission from super-massive black hole seeds can impact the global 21 cm signal at the level of tens to hundreds of percent provided that they were as radio loud as $z\approx1$ black holes and obscured by gas with column depths of $N_\text{H}\gtrsim 10^{23}$ cm$^{-2}$. We determine plausible sets of parameters that reproduce some of the striking features of the EDGES absorption feature including its depth, timing, and side steepness while producing radio/X-ray backgrounds and source counts that are consistent with published limits. Scenarios yielding a dramatic 21 cm signature also predict large populations of $\sim \mu$Jy point sources that will be detectable in future deep surveys from the Square Kilometer Array (SKA). Thus, 21 cm measurements, complemented by deep point source surveys, have the potential to constrain optimistic scenarios where super-massive black-hole progenitors were radio-loud., Comment: 21 pages. 15 figures. Accepted to MNRAS. Replacement matches accepted version

The discovery of Fast Radio Bursts (FRBs) at cosmological distances has opened a powerful window on otherwise unseen matter in the Universe. In the 2020s, observations of $>10^{4}$ FRBs will assess the baryon contents and physical conditions in the hot/diffuse circumgalactic, intracluster, and intergalactic medium, and test extant compact-object dark matter models., Comment: Science white paper submitted to the Astro2020 Decadal Survey. 15 pages, 3 color figures

Pulsar timing array (PTA) collaborations in North America, Australia, and Europe, have been exploiting the exquisite timing precision of millisecond pulsars over decades of observations to search for correlated timing deviations induced by gravitational waves (GWs). PTAs are sensitive to the frequency band ranging just below 1 nanohertz to a few tens of microhertz. The discovery space of this band is potentially rich with populations of inspiraling supermassive black-holes binaries, decaying cosmic string networks, relic post-inflation GWs, and even non-GW imprints of axionic dark matter. This article aims to provide an understanding of the exciting open science questions in cosmology, galaxy evolution, and fundamental physics that will be addressed by the detection and study of GWs through PTAs. The focus of the article is on providing an understanding of the mechanisms by which PTAs can address specific questions in these fields, and to outline some of the subtleties and difficulties in each case. The material included is weighted most heavily towards the questions which we expect will be answered in the near-term with PTAs; however, we have made efforts to include most currently anticipated applications of nanohertz GWs., Comment: Accepted to The Astronomy and Astrophysics Review (A&ARv)

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

The Next-Generation Very Large Array (ngVLA) has the potential to be a workhorse for the discovery and study of paired supermassive black holes either at large separations (dual) or in tightly bound systems (binary). In this chapter, we outline the science case for the study of these supermassive pairs, and summarize discovery methods that can be used at radio wavelengths to discover them: including morphological, spectral, and time-domain identifications. One critical aspect of this work is that multi-messenger binary black hole studies may be possible with the ngVLA when combined with gravitational-wave searches using pulsar timing array techniques. However, long-baseline interferometery (>>1000 km) will make this possibility more likely by expanding the redshift range at which radio emission arising from two separate black holes may be resolved and studied., Comment: 12 pages, 3 figures; Chapter in the ngVLA science book

We present VLA images and HST/STIS spectra of sources within the center of the brightest cluster galaxy (BCG) in Abell 2261. These observations were obtained to test the hypothesis that its extremely large, flat core reflects the ejection of its supermassive black hole. Spectra of three of the four most luminous "knots" embedded in the core were taken to test whether one may represent stars bound to a displaced massive black hole. The three knots have radial velocity offsets dV < ~150 km/s from the BCG. Knots 2 and 3 show kinematics, colors, and stellar masses consistent with infalling low-mass galaxies or larger stripped cluster members. Large errors in the stellar velocity dispersion of Knot 1, however, mean that we cannot rule out the hypothesis that it hosts a high-mass black hole. A2261-BCG has a compact, relic radio-source offset by 6.5 kpc (projected) from the optical core's center, but no active radio core that would pinpoint the galaxy's central black hole to a tight 10 GHz flux limit <3.6 uJy. Its spectrum and morphology are suggestive of an active galactic nucleus that switched off >48 Myr ago, with an equipartition condition magnetic field of 15 uG. These observations are still consistent with the hypothesis that the nuclear black hole has been ejected from its core, but the critical task of locating the supermassive black hole or demonstrating that A2261-BCG lacks one remains to be done., Comment: 12 pages, 9 figures, submitted to ApJ

This paper presents Very Long Baseline Array (VLBA) observations of 13 double-peaked [O III] emission-line type-2 Active Galactic Nuclei (AGNs) at redshifts 0.06 < z < 0.41 (with a median redshift of z~0.15) identified in the Sloan Digital Sky Survey. Such double-peaked emission-line objects may result from jets or outflows from the central engine or from a dual AGN. The VLBA provides an angular resolution of <~10 pc at the distance of many of these galaxies, sufficient to resolve the radio emission from extremely close dual AGNs and to contribute to understanding the origin of double-peaked [O III] emission lines. Of the 13 galaxies observed at 3.6 cm (8.4 GHz), we detect six at a 1\sigma\ sensitivity level of ~0.15 mJy/beam, two of which show clear jet structures on scales ranging from a few milliarcseconds to tens of milliarcseconds (corresponding to a few pc to tens of pc at a median redshift of 0.15). We suggest that radio-loud double-peaked emission-line type-2 AGNs may be indicative of jet produced structures, but a larger sample of double-peaked [O III] AGNs with high angular resolution radio observations will be required to confirm this suggestion., Comment: 14 pages, 7 figures; ApJ in press

Supermassive black hole binaries (SMBHBs) in the 10 million to 10 billion $M_\odot$ range form in galaxy mergers, and live in galactic nuclei with large and poorly constrained concentrations of gas and stars. There are currently no observations of merging SMBHBs--- it is in fact possible that they stall at their final parsec of separation and never merge. While LIGO has detected high frequency GWs, SMBHBs emit GWs in the nanohertz to millihertz band. This is inaccessible to ground-based interferometers, but possible with Pulsar Timing Arrays (PTAs). Using data from local galaxies in the 2 Micron All-Sky Survey, together with galaxy merger rates from Illustris, we find that there are on average $91\pm7$ sources emitting GWs in the PTA band, and $7\pm2$ binaries which will never merge. Local unresolved SMBHBs can contribute to GW background anisotropy at a level of $\sim20\%$, and if the GW background can be successfully isolated, GWs from at least one local SMBHB can be detected in 10 years., Comment: submitted to Nature Astronomy (reformatted for arXiv)

The precise localization of the repeating fast radio burst (FRB 121102) has provided the first unambiguous association (chance coincidence probability $p\lesssim3\times10^{-4}$) of an FRB with an optical and persistent radio counterpart. We report on optical imaging and spectroscopy of the counterpart and find that it is an extended ($0.6^{\prime\prime}-0.8^{\prime\prime}$) object displaying prominent Balmer and [OIII] emission lines. Based on the spectrum and emission line ratios, we classify the counterpart as a low-metallicity, star-forming, $m_{r^\prime} = 25.1$ AB mag dwarf galaxy at a redshift of $z=0.19273(8)$, corresponding to a luminosity distance of 972 Mpc. From the angular size, the redshift, and luminosity, we estimate the host galaxy to have a diameter $\lesssim4$ kpc and a stellar mass of $M_*\sim4-7\times 10^{7}\,M_\odot$, assuming a mass-to-light ratio between 2 to 3$\,M_\odot\,L_\odot^{-1}$. Based on the H$\alpha$ flux, we estimate the star formation rate of the host to be $0.4\,M_\odot\,\mathrm{yr^{-1}}$ and a substantial host dispersion measure depth $\lesssim 324\,\mathrm{pc\,cm^{-3}}$. The net dispersion measure contribution of the host galaxy to FRB 121102 is likely to be lower than this value depending on geometrical factors. We show that the persistent radio source at FRB 121102's location reported by Marcote et al (2017) is offset from the galaxy's center of light by $\sim$200 mas and the host galaxy does not show optical signatures for AGN activity. If FRB 121102 is typical of the wider FRB population and if future interferometric localizations preferentially find them in dwarf galaxies with low metallicities and prominent emission lines, they would share such a preference with long gamma ray bursts and superluminous supernovae., Comment: 12 pages, 3 figures, Published in ApJ Letters. V2: Corrected mistake in author list

A community meeting on the topic of "Radio Astronomy in the LSST Era" was hosted by the National Radio Astronomy Observatory in Charlottesville, VA (2013 May 6--8). The focus of the workshop was on time domain radio astronomy and sky surveys. For the time domain, the extent to which radio and visible wavelength observations are required to understand several classes of transients was stressed, but there are also classes of radio transients for which no visible wavelength counterpart is yet known, providing an opportunity for discovery. From the LSST perspective, the LSST is expected to generate as many as 1 million alerts nightly, which will require even more selective specification and identification of the classes and characteristics of transients that can warrant follow up, at radio or any wavelength. The LSST will also conduct a deep survey of the sky, producing a catalog expected to contain over 38 billion objects in it. Deep radio wavelength sky surveys will also be conducted on a comparable time scale, and radio and visible wavelength observations are part of the multi-wavelength approach needed to classify and understand these objects. Radio wavelengths are valuable because they are unaffected by dust obscuration and, for galaxies, contain contributions both from star formation and from active galactic nuclei. The workshop touched on several other topics, on which there was consensus including the placement of other LSST "Deep Drilling Fields," inter-operability of software tools, and the challenge of filtering and exploiting the LSST data stream. There were also topics for which there was insufficient time for full discussion or for which no consensus was reached, which included the procedures for following up on LSST observations and the nature for future support of researchers desiring to use LSST data products., Comment: Conference summary, 29 pages, 1 figure; to be published in the Publ. Astron. Soc. Pacific; full science program and presentations available at http://science.nrao.edu/science/event/RALSST2013

This paper considers the suitability of a number of emerging and future instruments for the study of radio recombination lines (RRLs) at frequencies below 200 MHz. These lines arise only in low-density regions of the ionized interstellar medium, and they may represent a frequency-dependent foreground for next-generation experiments trying to detect H I signals from the Epoch of Reionization and Dark Ages ("21-cm cosmology"). We summarize existing decametre-wavelength observations of RRLs, which have detected only carbon RRLs. We then show that, for an interferometric array, the primary instrumental factor limiting detection and study of the RRLs is the areal filling factor of the array. We consider the Long Wavelength Array (LWA-1), the LOw Frequency ARray (LOFAR), the low-frequency component of the Square Kilometre Array (SKA-lo), and a future Lunar Radio Array (LRA), all of which will operate at decametre wavelengths. These arrays offer digital signal processing, which should produce more stable and better defined spectral bandpasses; larger frequency tuning ranges; and better angular resolution than that of the previous generation of instruments that have been used in the past for RRL observations. Detecting Galactic carbon RRLs, with optical depths at the level of 10^-3, appears feasible for all of these arrays, with integration times of no more than 100 hr. The SKA-lo and LRA, and the LWA-1 and LOFAR at the lowest frequencies, should have a high enough filling factor to detect lines with much lower optical depths, of order 10^-4 in a few hundred hours. The amount of RRL-hosting gas present in the Galaxy at the high Galactic latitudes likely to be targeted in 21-cm cosmology studies is currently unknown. If present, however, the spectral fluctuations from RRLs could be comparable to or exceed the anticipated H I signals., Comment: 9 pages; Astron. & Astrophys., in press

We report detection of strong circularly polarized emission from the transient bursting source GCRT J1745-3009 based on new analysis of 325 MHz GMRT observations conducted on 28 September 2003. We place 8 Solar radius as the upper limit on the size of the emission region. The implied high brightness temperature required for an object beyond 1 pc and the high fraction of circular polarization firmly establish the emission as coherent. Electron cyclotron or plasma emission from a highly subsolar magnetically dominated dwarf located less than 4 kpc away could have given rise to the GCRT radio emission., Comment: 17 pages, 3 figures, accepted for publication in ApJ Letters