Source author record

Sarah Burke-Spolaor

Sarah Burke-Spolaor appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

32works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

32 published item(s)

preprint2026arXiv

The NANOGrav 15 yr Data Set: Piecewise Power-Law Reconstruction of the Gravitational-Wave Background

The NANOGrav 15-year (NG15) data set provides evidence for a gravitational-wave background (GWB) signal at nanohertz frequencies, which is expected to originate either from a cosmic population of inspiraling supermassive black-hole binaries or new particle physics in the early Universe. A firm identification of the source of the NG15 signal requires an accurate reconstruction of its frequency spectrum. In this paper, we provide such a spectral characterization of the NG15 signal based on a piecewise power-law (PPL) ansatz that strikes a balance between existing alternatives in the literature. Our PPL reconstruction is more flexible than the standard constant-power-law model, which describes the GWB spectrum in terms of only two parameters: an amplitude A and a spectral index gamma. Concurrently, it better approximates physically realistic GWB spectra -- especially those of cosmological origin -- than the free spectral model, since the latter allows for arbitrary variations in the GWB amplitude from one frequency bin to the next. Our PPL reconstruction of the NG15 signal relies on individual PPL models with a fixed number of internal nodes (i.e., constant power law, broken power law, doubly broken power law, etc.) that are ultimately combined in a Bayesian model average. The data products resulting from our analysis provide the basis for fast refits of spectral GWB models.

preprint2022arXiv

Deep VLBI Observations Challenge Previous Evidence of a Binary Supermassive Black Hole Residing in the Seyfert Galaxy NGC 7674

Previous Ku-band (15 GHz) imaging with data obtained from the Very Long Baseline Array (VLBA) had shown two compact, sub-pc components at the location of a presumed kpc-scale radio core in the Seyfert galaxy NGC 7674. It was then presumed that these two unresolved and compact components were dual radio cores corresponding to two supermassive black holes (SMBHs) accreting surrounding gas and launching radio-bright relativistic jets. However, utilizing the original VLBA dataset used to claim the detection of a binary SMBH, in addition to later multi-epoch/multi-frequency datatsets obtained from both the VLBA and the European VLBI Network, we find no evidence to support the presence of a binary SMBH. We place stringent upper limits to the flux densities of any sub-pc-scale radio cores which are at least an order of magnitude lower than the original VLBI radio-core detections, directly challenging the original binary SMBH detection claim. With this in mind, we discuss the possible reasons for the non-detection of any VLBI radio cores in our imaging, the possibility of a binary SMBH still residing in NGC 7674, and the prospect of future observations shedding further light on the true nature of this active galactic nucleus.

preprint2022arXiv

Quasars with Periodic Variability: Capabilities and Limitations of Bayesian Searches for Supermassive Black Hole Binaries in Time-Domain Surveys

Supermassive black hole binaries (SMBHBs) are an inevitable consequence of galaxy mergers. At subparsec separations, they are practically impossible to resolve, and the most promising technique is to search for quasars with periodic variability. However, searches for quasar periodicity in time-domain data are challenging due to the stochastic variability of quasars. In this paper, we used Bayesian methods to disentangle periodic SMBHB signals from intrinsic damped random walk (DRW) variability in active galactic nuclei light curves. We simulated a wide variety of realistic DRW and DRW+sine light curves. Their observed properties are modeled after the Catalina Real-time Transient Survey (CRTS) and expected properties of the upcoming Legacy Survey of Space and Time (LSST) from the Vera C. Rubin Observatory. Through a careful analysis of parameter estimation and Bayesian model selection, we investigated the range of parameter space for which binary systems can be detected. We also examined which DRW signals can mimic periodicity and be falsely classified as binary candidates. We found that periodic signals are more easily detectable if the period is short or the amplitude of the signal is large compared to the contribution of the DRW noise. We saw similar detection rates both in the CRTS and LSST-like simulations, while the false-detection rate depends on the quality of the data and is minimal in LSST. Our idealized simulations provide an excellent way to uncover the intrinsic limitations in quasar periodicity searches and set the stage for future searches for SMBHBs.

preprint2022arXiv

Searching For Gravitational Waves From Cosmological Phase Transitions With The NANOGrav 12.5-year dataset

We search for a first-order phase transition gravitational wave signal in 45 pulsars from the NANOGrav 12.5 year dataset. We find that the data can be modeled in terms of a strong first order phase transition taking place at temperatures below the electroweak scale. However, we do not observe any strong preference for a phase-transition interpretation of the signal over the standard astrophysical interpretation in terms of supermassive black holes mergers; but we expect to gain additional discriminating power with future datasets, improving the signal to noise ratio and extending the sensitivity window to lower frequencies. An interesting open question is how well gravitational wave observatories could separate such signals.

preprint2021arXiv

Characterizing the FRB host galaxy population and its connection to transients in the local and extragalactic Universe

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.

preprint2021arXiv

The host galaxy and persistent radio counterpart of FRB 20201124A

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.

preprint2021arXiv

The NANOGrav 12.5-year Data Set: Search For An Isotropic Stochastic Gravitational-Wave Background

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.

preprint2020arXiv

A Demonstration of Extremely Low Latency $γ$-ray, X-Ray & UV Follow-Up of a Millisecond Radio Transient

We report results of a novel high-energy follow-up observation of a potential Fast Radio Burst. The radio burst was detected by VLA/realfast and followed-up by the Neil Gehrels Swift Observatory in very low latency utilizing new operational capabilities of Swift (arXiv:2005.01751), with pointed soft X-ray and UV observations beginning at T0+32 minutes, and hard X-ray/gamma-ray event data saved around T0. These observations are $>10$x faster than previous X-ray/UV follow-up of any radio transient to date. No emission is seen coincident with the FRB candidate at T0, with a 0.2s fluence $5σ$ upper limit of $1.35\times10^{-8}$ erg cm$^{-2}$ (14-195 keV) for a SGR 1935+2154-like flare, nor at T0+32 minutes down to $3σ$ upper limits of 22.18 AB mag in UVOT u band, and $3.33\times10^{-13}$ erg cm$^{-2}$ s$^{-1}$ from 0.3-10 keV for the 2 ks observation. The candidate FRB alone is not significant enough to be considered astrophysical, so this note serves as a technical demonstration. These new Swift operational capabilities will allow future FRB detections to be followed up with Swift at even lower latencies than demonstrated here: 15-20 minutes should be regularly achievable, and 5-10 minutes occasionally achievable. We encourage FRB detecting facilities to release alerts in low latency to enable this science.

preprint2020arXiv

Chandra Observations of Abell 2261 Brightest Cluster Galaxy, a Candidate Host to a Recoiling Black Hole

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.

preprint2020arXiv

FETCH: A deep-learning based classifier for fast transient classification

With the upcoming commensal surveys for Fast Radio Bursts (FRBs), and their high candidate rate, usage of machine learning algorithms for candidate classification is a necessity. Such algorithms will also play a pivotal role in sending real-time triggers for prompt follow-ups with other instruments. In this paper, we have used the technique of Transfer Learning to train the state-of-the-art deep neural networks for classification of FRB and Radio Frequency Interference (RFI) candidates. These are convolutional neural networks which work on radio frequency-time and dispersion measure-time images as the inputs. We trained these networks using simulated FRBs and real RFI candidates from telescopes at the Green Bank Observatory. We present 11 deep learning models, each with an accuracy and recall above 99.5% on our test dataset comprising of real RFI and pulsar candidates. As we demonstrate, these algorithms are telescope and frequency agnostic and are able to detect all FRBs with signal-to-noise ratios above 10 in ASKAP and Parkes data. We also provide an open-source python package FETCH (Fast Extragalactic Transient Candidate Hunter) for classification of candidates, using our models. Using FETCH, these models can be deployed along with any commensal search pipeline for real-time candidate classification.

preprint2020arXiv

Multi-Messenger Gravitational Wave Searches with Pulsar Timing Arrays: Application to 3C66B Using the NANOGrav 11-year Data Set

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.

preprint2020arXiv

VLA/Realfast Detection of a Burst from FRB180916.J0158+65 and Tests for Periodic Activity

We report on the detection of a burst from FRB180916 by realfast/VLA and present software for interpreting fast radio bursts (FRB) periodicity. We demonstrate a range of periodicity analyses with bursts from FRB180916, FRB121102 and FRB180814. Our results for FRB180916 and FRB121102 are consistent with published results. For FRB180814, we did not detect any significant periodic episodes. The realfast-detected and other high-frequency bursts for FRB180916 tend to lie at the beginning of the activity window, indicating a possible phase-frequency relation. The python package $\texttt{frbpa}$ can be used to reproduce and expand on this analysis to test models for repeating FRBs.

preprint2016arXiv

Constraints on Black Hole/Host Galaxy Co-evolution and Binary Stalling Using Pulsar Timing Arrays

Pulsar timing arrays are now setting increasingly tight limits on the gravitational wave background from binary supermassive black holes. But as upper limits grow more constraining, what can be implied about galaxy evolution? We investigate which astrophysical parameters have the largest impact on strain spectrum predictions and provide a simple framework to directly translate between measured values for the parameters of galaxy evolution and PTA limits on the gravitational wave background of binary supermassive black holes. We find that the most influential observable is the relation between a host galaxy's central bulge and its central black hole, $\mbox{$M_{\bullet}$-$M_{\rm bulge}$}$, which has the largest effect on the mean value of the characteristic strain amplitude. However, the variance of each prediction is dominated by uncertainties in the galaxy stellar mass function. Using this framework with the best published PTA limit, we can set limits on the shape and scatter of the $\mbox{$M_{\bullet}$-$M_{\rm bulge}$}$ relation. We find our limits to be in contention with strain predictions using two leading measurements of this relation. We investigate several possible reasons for this disagreement. If we take the $\mbox{$M_{\bullet}$-$M_{\rm bulge}$}$ relations to be correct within a simple power-law model for the gravitational wave background, then the inconsistency is reconcilable by allowing for an additional "stalling" time between a galaxy merger and evolution of a binary supermassive black hole to sub-parsec scales, with lower limits on this timescale of $\sim 1-2$ Gyr.

preprint2016arXiv

Gravitational-wave cosmology across 29 decades in frequency

Quantum fluctuations of the gravitational field in the early Universe, amplified by inflation, produce a primordial gravitational-wave background across a broad frequency band. We derive constraints on the spectrum of this gravitational radiation, and hence on theories of the early Universe, by combining experiments that cover 29 orders of magnitude in frequency. These include Planck observations of cosmic microwave background temperature and polarization power spectra and lensing, together with baryon acoustic oscillations and big bang nucleosynthesis measurements, as well as new pulsar timing array and ground-based interferometer limits. While individual experiments constrain the gravitational-wave energy density in specific frequency bands, the combination of experiments allows us to constrain cosmological parameters, including the inflationary spectral index, $n_t$, and the tensor-to-scalar ratio, $r$. Results from individual experiments include the most stringent nanohertz limit of the primordial background to date from the Parkes Pulsar Timing Array, $Ω_{\rm gw}(f)<2.3\times10^{-10}$. Observations of the cosmic microwave background alone limit the gravitational-wave spectral index at 95\% confidence to $n_t\lesssim5$ for a tensor-to-scalar ratio of $r = 0.11$. However, the combination of all the above experiments limits $n_t<0.36$. Future Advanced LIGO observations are expected to further constrain $n_t<0.34$ by 2020. When cosmic microwave background experiments detect a non-zero $r$, our results will imply even more stringent constraints on $n_t$ and hence theories of the early Universe.

preprint2016arXiv

Limits on Fast Radio Bursts from Four Years of the V-FASTR Experiment

The V-FASTR experiment on the Very Long Baseline Array was designed to detect dispersed pulses of milliseconds duration, such as fast radio bursts (FRBs). We use all V-FASTR data through February 2015 to report V-FASTR's upper limits on the rates of FRBs, and compare these with re-derived rates from Parkes FRB detection experiments. V-FASTR's operation at lambda=20 cm allows direct comparison with the 20 cm Parkes rate, and we derive a power-law limit of γ<-0.4 (95% confidence limit) on the index of FRB source counts, N(>S)\propto S^γ. Using the previously measured FRB rate and the unprecedented amount of survey time spent searching for FRBs at a large range of wavelengths (0.3 cm > λ> 90 cm), we also place frequency-dependent limits on the spectral distribution of FRBs. The most constraining frequencies place two-point spectral index limits of α_{20cm}^{4cm} < 5.8 and α_{90cm}^{20cm} > -7.6, where fluence F \propto f^αif we assume true the burst rate reported by Champion et al. (2016) of R(F~0.6 Jy ms) = 7 x 10^3 sky^{-1} day^{-1} (for bursts of ~3 ms duration). This upper limit on αsuggests that if FRBs are extragalactic but non-cosmological, that on average they are not experiencing excessive free-free absorption due to a medium with high optical depth (assuming temperature ~8,000 K), which excessively invert their low-frequency spectrum. This in turn implies that the dispersion of FRBs arises in either or both of the intergalactic medium or the host galaxy, rather than from the source itself.

preprint2016arXiv

Rare Event Statistics Applied to Fast Radio Bursts

Statistical interpretation of sparsely sampled event rates has become vital for new transient surveys, particularly those aimed at detecting fast radio bursts (FRBs). We provide an accessible reference for a number of simple, but critical, statistical questions relevant for current transient and FRB research and utilizing the negative binomial model for counts in which the count rate parameter is uncertain or randomly biased from one study to the next. We apply these methods to re-assess and update results from previous FRB surveys, finding as follows. 1) Thirteen FRBs detected across five high-Galactic-latitude (> 30$^\circ$) surveys are highly significant $(p = 5\times 10^{-5})$ evidence of a higher rate relative to the single FRB detected across four low-latitude (< 5$^\circ$) surveys, even after accounting for effects that dampen Galactic plane sensitivity. High- vs. mid-latitude (5 to 15$^\circ$) is marginally significant $(p = 0.03)$. 2) A meta analysis of twelve heterogeneous surveys gives an FRB rate of 2866 sky$^{-1}$ day$^{-1}$ above 1 Jy at high Galactic latitude (95% confidence 1121 to 7328) and 285 sky$^{-1}$ day$^{-1}$ at low/mid latitudes (95% from 48 to 1701). 3) Using the Parkes HTRU high-latitude setup requires 193 observing hours to achieve 50% probability of detecting an FRB and 937 hours to achieve 95% probability, based on the ten detections of (Champion et al. 2016) and appropriately accounting for uncertainty in the unknown Poisson rate. 4) Two quick detections at Parkes from a small number of high-latitude fields (Ravi et al. 2015; Petroff et al. 2015) tentatively favor a look long survey style relative to the scan wide HTRU survey, but only at $p = 0.07$ significance.

preprint2015arXiv

Gravitational-Wave Detection and Astrophysics with Pulsar Timing Arrays

We have begun an exciting era for gravitational wave detection, as several world-leading experiments are breaching the threshold of anticipated signal strengths. Pulsar timing arrays (PTAs) are pan-Galactic gravitational wave detectors that are already cutting into the expected strength of gravitational waves from cosmic strings and binary supermassive black holes in the nHz-$μ$Hz gravitational wave band. These limits are leading to constraints on the evolutionary state of the Universe. Here, we provide a broad review of this field, from how pulsars are used as tools for detection, to astrophysical sources of uncertainty in the signals PTAs aim to see, to the primary current challenge areas for PTA work. This review aims to provide an up-to-date reference point for new parties interested in the field of gravitational wave detection via pulsar timing.

preprint2015arXiv

Next Generation Very Large Array Memo No. 9 Science Working Group 4: Time Domain, Fundamental Physics, and Cosmology

We report here on key science topics for the Next Generation Very Large Array in the areas of time domain, fundamental physics, and cosmology. Key science cases considered are pulsars in orbit around the Galactic Center massive black hole, Sagittarius A*, electromagnetic counterparts to gravitational waves, and astrometric cosmology. These areas all have the potential for ground-breaking and transformative discovery. Numerous other topics were discussed during the preparation of this report and some of those discussions are summarized here, as well. There is no doubt that further investigation of the science case will reveal rich and compelling opportunities.

preprint2015arXiv

The NANOGrav Nine-year Data Set: Limits on the Isotropic Stochastic Gravitational Wave Background

We compute upper limits on the nanohertz-frequency isotropic stochastic gravitational wave background (GWB) using the 9-year data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) collaboration. We set upper limits for a GWB from supermassive black hole binaries under power law, broken power law, and free spectral coefficient GW spectrum models. We place a 95\% upper limit on the strain amplitude (at a frequency of yr$^{-1}$) in the power law model of $A_{\rm gw} < 1.5\times 10^{-15}$. For a broken power law model, we place priors on the strain amplitude derived from simulations of Sesana (2013) and McWilliams et al. (2014). We find that the data favor a broken power law to a pure power law with odds ratios of 22 and 2.2 to one for the McWilliams and Sesana prior models, respectively. The McWilliams model is essentially ruled out by the data, and the Sesana model is in tension with the data under the assumption of a pure power law. Using the broken power-law analysis we construct posterior distributions on environmental factors that drive the binary to the GW-driven regime including the stellar mass density for stellar-scattering, mass accretion rate for circumbinary disk interaction, and orbital eccentricity for eccentric binaries, marking the first time that the shape of the GWB spectrum has been used to make astrophysical inferences. We then place the most stringent limits so far on the energy density of relic GWs, $Ω_\mathrm{gw}(f)\,h^2 < 4.2 \times 10^{-10}$, yielding a limit on the Hubble parameter during inflation of $H_*=1.6\times10^{-2}~m_{Pl}$, where $m_{Pl}$ is the Planck mass. Our limit on the cosmic string GWB, $Ω_\mathrm{gw}(f)\, h^2 < 2.2 \times 10^{-10}$, translates to a conservative limit of $Gμ<3.3\times 10^{-8}$ - a factor of 4 better than the joint Planck and high-$l$ CMB data from other experiments.

preprint2014arXiv

A hunt for dual radio active galactic nuclei in the VLASS

This whitepaper describes how the VLASS could be designed in a manner to allow the identification of candidate dual active galactic nuclei (AGN) at separations <7 kpc. Dual AGN represent a clear marker of two supermassive black holes within an ongoing merger. A dual AGN survey will provide a wealth of studies in structure growth and gravitational-wave science. Radio wavelengths are ideal for identifying close pairs, as disturbed stellar and gaseous material can obscure their presence in optical and shorter wavelengths. With sufficiently high resolution and sensitivity, a large-scale radio imaging survey like the VLASS will uncover many of these systems and provide the means to broadly study the radio properties of candidate dual systems revealed at other wavelengths. We determine that the ideal survey for our purposes will be at as high a resolution as possible, with significantly more science return in A array at L-band or higher, or B array at C-band or higher. We describe a range of potential survey parameters within this document. Based on the analysis outlined in this whitepaper, our ideal survey would create a catalogue of $\gtrsim$100 dual AGN in either: 1) a medium-sensitivity (~1 mJy detection threshold), wide-field (few thousand square degree) survey, or 2) a high-sensitivity (~10 uJy threshold) survey of several hundred square degrees.

preprint2014arXiv

A Millisecond Interferometric Search for Fast Radio Bursts with the Very Large Array

We report on the first millisecond timescale radio interferometric search for the new class of transient known as fast radio bursts (FRBs). We used the Very Large Array (VLA) for a 166-hour, millisecond imaging campaign to detect and precisely localize an FRB. We observed at 1.4 GHz and produced visibilities with 5 ms time resolution over 256 MHz of bandwidth. Dedispersed images were searched for transients with dispersion measures from 0 to 3000 pc/cm3. No transients were detected in observations of high Galactic latitude fields taken from September 2013 though October 2014. Observations of a known pulsar show that images typically had a thermal-noise limited sensitivity of 120 mJy/beam (8 sigma; Stokes I) in 5 ms and could detect and localize transients over a wide field of view. Our nondetection limits the FRB rate to less than 7e4/sky/day (95% confidence) above a fluence limit of 1.2 Jy-ms. Assuming a Euclidean flux distribution, the VLA rate limit is inconsistent with the published rate of Thornton et al. We recalculate previously published rates with a homogeneous consideration of the effects of primary beam attenuation, dispersion, pulse width, and sky brightness. This revises the FRB rate downward and shows that the VLA observations had a roughly 60% chance of detecting a typical FRB and that a 95% confidence constraint would require roughly 500 hours of similar VLA observing. Our survey also limits the repetition rate of an FRB to 2 times less than any known repeating millisecond radio transient.

preprint2013arXiv

A framework for interpreting fast radio transients search experiments: application to the V-FASTR experiment

We define a framework for determining constraints on the detection rate of fast transient events from a population of underlying sources, with a view to incorporating beam shape, frequency effects, scattering effects, and detection efficiency into the metric. We then demonstrate a method for combining independent datasets into a single event rate constraint diagram, using a probabilistic approach to the limits on parameter space. We apply this new framework to present the latest results from the V-FASTR experiment, a commensal fast transients search using the Very Long Baseline Array (VLBA). In the 20 cm band, V-FASTR now has the ability to probe the regions of parameter space of importance for the observed Lorimer and Keane fast radio transient candidates, by combining the information from observations with differing bandwidths, and properly accounting for the source dispersion measure, VLBA antenna beam shape, experiment time sampling, and stochastic nature of events. We then apply the framework to combine the results of the V-FASTR and ATA Fly's Eye experiments, demonstrating their complementarity. Expectations for fast transients experiments for the SKA Phase I dish array are then computed, and the impact of large differential bandwidths is discussed.

preprint2013arXiv

Multi-messenger approaches to binary supermassive black holes in the "continuous-wave" regime

Pulsar timing arrays are sensitive to gravitational waves from supermassive black hole (SMBH) binaries at orbital separations of << 1pc. There is currently an observational paucity of such systems, although they are central figures in studies of galaxy evolution, merger dynamics, and active nucleus formation. We review the prospects of detecting SMBH binaries through electromagnetic radiative processes thought to be associated with galaxy mergers and late-stage binary evolution. We then discuss the scientific goals of joint pulsar timing and electromagnetic studies of these systems, including the facilitation of binary parameter estimation, identifying galactic hosts of gravitational wave emitters, and relevant studies of merger dynamics and cosmology. The use of upcoming high-precision timing arrays with the International Pulsar Timing Array and the Square Kilometre Array, combined with ongoing electromagnetic observing campaigns to identify active SMBH binaries, provide generous possibilities for multi-messenger astrophysics in the near future.

preprint2013arXiv

Real Time Event Detection in Astronomical Data Streams: Lessons from the VLBA

A new generation of observational science instruments is dramatically increasing collected data volumes in a range of fields. These instruments include the Square Kilometre Array (SKA), Large Synoptic Survey Telescope (LSST), terrestrial sensor networks, and NASA satellites participating in "decadal survey" missions. Their unprecedented coverage and sensitivity will likely reveal wholly new categories of unexpected and transient events. Commensal methods passively analyze these data streams, recognizing anomalous events of scientific interest and reacting in real time. We report on a case example: V-FASTR, an ongoing commensal experiment at the Very Long Baseline Array (VLBA) that uses online adaptive pattern recognition to search for anomalous fast radio transients. V-FASTR triages a millisecond-resolution stream of data and promotes candidate anomalies for further offline analysis. It tunes detection parameters in real time, injecting synthetic events to continually retrain itself for optimum performance. This self-tuning approach retains sensitivity to weak signals while adapting to changing instrument configurations and noise conditions. The system has operated since July 2011, making it the longest-running real time commensal radio transient experiment to date.

preprint2013arXiv

The Gravitational-Wave Discovery Space of Pulsar Timing Arrays

Recent years have seen a burgeoning interest in using pulsar timing arrays (PTAs) as gravitational-wave (GW) detectors. To date, that interest has focused mainly on three particularly promising source types: supermassive--black-hole binaries, cosmic strings, and the stochastic background from early-Universe phase transitions. In this paper, by contrast, our aim is to investigate the PTA potential for discovering unanticipated sources. We derive significant constraints on the available discovery space based solely on energetic and statistical considerations: we show that a PTA detection of GWs at frequencies above ~3.e-5 Hz would either be an extraordinary coincidence or violate "cherished beliefs;" we show that for PTAs GW memory can be more detectable than direct GWs, and that, as we consider events at ever higher redshift, the memory effect increasingly dominates an event's total signal-to-noise ratio. The paper includes also a simple analysis of the effects of pulsar red noise in PTA searches, and a demonstration that the effects of periodic GWs in the 10^-8 -- 10^-4.5 Hz band would not be degenerate with small errors in standard pulsar parameters (except in a few narrow bands).

preprint2012arXiv

Multi-wavelength Observations of the Radio Magnetar PSR J1622-4950 and Discovery of its Possibly Associated Supernova Remnant

We present multi-wavelength observations of the radio magnetar PSR J1622-4950 and its environment. Observations of PSR J1622-4950 with Chandra (in 2007 and 2009) and XMM (in 2011) show that the X-ray flux of PSR J1622-4950 has decreased by a factor of ~50 over 3.7 years, decaying exponentially with a characteristic time of 360 +/- 11 days. This behavior identifies PSR J1622-4950 as a possible addition to the small class of transient magnetars. The X-ray decay likely indicates that PSR J1622-4950 is recovering from an X-ray outburst that occurred earlier in 2007, before the 2007 Chandra observations. Observations with the Australia Telescope Compact Array show strong radio variability, including a possible radio flaring event at least one and a half years after the 2007 X-ray outburst that may be a direct result of this X-ray event. Radio observations with the Molonglo Observatory Synthesis Telescope reveal that PSR J1622-4950 is 8' southeast of a diffuse radio arc, G333.9+0.0, which appears non-thermal in nature and which could possibly be a previously undiscovered supernova remnant. If G333.9+0.0 is a supernova remnant then the estimates of its size and age, combined with the close proximity and reasonable implied velocity of PSR J1622-4950, suggests that these two objects could be physically associated.

preprint2011arXiv

Polarization observations of 20 millisecond pulsars

Polarization profiles are presented for 20 millisecond pulsars that are being observed as part of the Parkes Pulsar Timing Array project. The observations used the Parkes multibeam receiver with a central frequency of 1369 MHz and the Parkes digital filterbank pulsar signal-processing system PDFB2. Because of the large total observing time, the summed polarization profiles have very high signal/noise ratios and show many previously undetected profile features. Thirteen of the 20 pulsars show emission over more than half of the pulse period. Polarization variations across the profiles are complex and the observed position angle variations are generally not in accord with the rotating-vector model for pulsar polarization. Never-the-less, the polarization properties are broadly similar to those of normal (non-millisecond) pulsars, suggesting that the basic radio emission mechanism is the same in both classes of pulsar. The results support the idea that radio emission from millisecond pulsars originates high in the pulsar magnetosphere, probably close to the emission regions for high-energy X-ray and gamma-ray emission. Rotation measures were obtained for all 20 pulsars, eight of which had no previously published measurements.

preprint2011arXiv

Rotation measure variations for 20 millisecond pulsars

We report on variations in the mean position angle of the 20 millisecond pulsars being observed as part of the Parkes Pulsar Timing Array (PPTA) project. It is found that the observed variations are dominated by changes in the Faraday rotation occurring in the Earth's ionosphere. Two ionospheric models are used to correct for the ionospheric contribution and it is found that one based on the International Reference Ionosphere gave the best results. Little or no significant long-term variation in interstellar RM was found with limits typically about 0.1 rad m$^{-2}$ yr$^{-1}$ in absolute value. In a few cases, apparently significant RM variations over timescales of a few 100 days or more were seen. These are unlikely to be due to localised magnetised regions crossing the line of sight since the implied magnetic fields are too high. Most probably they are statistical fluctuations due to random spatial and temporal variations in the interstellar electron density and magnetic field along the line of sight.

preprint2010arXiv

A Radio Census of Binary Supermassive Black Holes

Using archival VLBI data for 3114 radio-luminous active galactic nuclei, we searched for binary supermassive black holes using a radio spectral index mapping technique which targets spatially resolved, double radio-emitting nuclei. Only one source was detected as a double nucleus. This result is compared with a cosmological merger rate model and interpreted in terms of (1) implications for post-merger timescales for centralisation of the two black holes, (2) implications for the possibility of "stalled" systems, and (3) the relationship of radio activity in nuclei to mergers. Our analysis suggests that the binary evolution of paired supermassive black holes (both of masses >= 1e8 Msun) spends less than 500 Myr in progression from the merging of galactic stellar cores to within the purported stalling radius for supermassive black hole pairs. The data show no evidence for an excess of stalled binary systems at small separations. We see circumstantial evidence that the relative state of radio emission between paired supermassive black holes is correlated within orbital separations of 2.5 kpc.

preprint2010arXiv

Radio Bursts with Extragalactic Spectral Characteristics Show Terrestrial Origins

Three years ago, the report of a solitary radio burst was thought to be the first discovery of a rare, impulsive event of unknown extragalactic origin (Lorimer et al. 2007). The extragalactic interpretation was based on the swept-frequency nature of the event, which followed the dispersive delay expected from an extragalactic pulse. We report here on the detection of 16 pulses, the bulk of which exhibit a frequency sweep with a shape and magnitude resembling the Lorimer Burst. These new events were detected in a sidelobe of the Parkes Telescope and are of clearly terrestrial origin, with properties unlike any known sources of terrestrial broad-band radio emission. The new detections cast doubt on the extragalactic interpretation of the original burst, and call for further sophistication in radio-pulse survey techniques to identify the origin of the anomalous terrestrial signals and definitively distinguish future extragalactic pulse detections from local signals. The ambiguous origin of these seemingly dispersed, swept-frequency signals suggest that radio-pulse searches using multiple detectors will be the only experiments able to provide definitive information about the origin of new swept-frequency radio burst detections.

preprint2010arXiv

The Australia Telescope 20GHz (AT20G) Survey: analysis of the extragalactic source sample

The Australia Telescope 20 GHz (AT20G) survey is a blind survey of the whole Southern sky at 20 GHz with follow-up observations at 4.8, 8.6, and 20 GHz carried out with the Australia Telescope Compact Array (ATCA). In this paper we present an analysis of radio spectral properties in total intensity and polarisation, sizes, optical identifications, and redshifts of the sample of the 5808 extragalactic sources in the survey catalogue of confirmed sources over the whole Southern sky excluding the strip at Galactic latitude |b|<1.5deg. The sample has a flux density limit of 40 mJy. Completeness has been measured as a function of scan region and flux density. Averaging over the whole survey area the follow-up survey is 78% complete above 50mJy and 93% complete above 100mJy. 3332 sources with declination <-15deg have good quality almost simultaneous observations at 4.8, 8.6, and 20GHz. The spectral analysis shows that the sample is dominated by flat-spectrum sources. The fraction of flat-spectrum sources decreases from 81% for 20GHz flux densities S>500mJy, to 60% for S<100mJy. There is also a clear spectral steepening at higher frequencies with the median spectral index decreasing from -0.16 between 4.8 and 8.6GHz to -0.28 between 8.6 and 20GHz. Simultaneous observations in polarisation are available for all the sources at all the frequencies. 768 sources have a good quality detection of polarised flux density at 20GHz; 467 of them were also detected in polarisation at 4.8 and/or at 8.6GHz so that it has been possible to compare the spectral behaviour in total intensity and polarisation. We have found that the polarised fraction increases slightly with frequency and decreases with flux density. Cross matches and comparisons have been made with other catalogues at lower radio frequencies, and in the optical, X-ray and gamma-ray bands. Redshift estimates are available for 825 sources.

preprint2009arXiv

The Australia Telescope 20 GHz Survey: The Source Catalogue

We present the full source catalogue from the Australia Telescope 20 GHz (AT20G) Survey. The AT20G is a blind radio survey carried out at 20 GHz with the Australia Telescope Compact Array (ATCA) from 2004 to 2008, and covers the whole sky south of declination 0 deg. The AT20G source catalogue presented here is an order of magnitude larger than any previous catalogue of high-frequency radio sources, and includes 5890 sources above a 20 GHz flux-density limit of 40 mJy. All AT20G sources have total intensity and polarisation measured at 20 GHz, and most sources south of declination -15 deg also have near-simultaneous flux-density measurements at 5 and 8 GHz. A total of 1559 sources were detected in polarised total intensity at one or more of the three frequencies. We detect a small but significant population of non-thermal sources that are either undetected or have only weak detections in low-frequency catalogues. We introduce the term Ultra-Inverted Spectrum (UIS) to describe these radio sources, which have a spectral index alpha(5, 20) > +0.7 and which constitute roughly 1.2 per cent of the AT20G sample. The 20 GHz flux densities measured for the strongest AT20G sources are in excellent agreement with the WMAP 5-year source catalogue of Wright et al. (2009), and we find that the WMAP source catalogue is close to complete for sources stronger than 1.5 Jy at 23 GHz.