skip to main content

Title: Multimessenger Gravitational-wave Searches with Pulsar Timing Arrays: Application to 3C 66B Using the NANOGrav 11-year Data Set
Abstract When galaxies merge, the supermassive black holes in their centers may form binaries and emit low-frequency gravitational radiation in the process. In this paper, we consider the galaxy 3C 66B, which was used as the target of the first multimessenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data 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 3C 66B, 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 3C 66B to less than (1.65 ± 0.02) × 10 9   M ⊙ 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 more » made by the inclusion of source properties gleaned from electromagnetic data over “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. « less
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; « less
Award ID(s):
Publication Date:
Journal Name:
The Astrophysical Journal
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Pulsar timing array (PTA) experiments are becoming increasingly sensitive to gravitational waves (GWs) in the nanohertz frequency range, where the main astrophysical sources are supermassive black hole binaries (SMBHBs), which are expected to form following galaxy mergers. Some of these individual SMBHBs may power active galactic nuclei, and thus their binary parameters could be obtained electromagnetically, which makes it possible to apply electromagnetic (EM) information to aid the search for a GW signal in PTA data. In this work, we investigate the effects of such an EM-informed search on binary detection and parameter estimation by performing mock data analyses on simulated PTA data sets. We find that by applying EM priors, the Bayes factor of some injected signals with originally marginal or sub-threshold detectability (i.e., Bayes factor ∼1) can increase by a factor of a few to an order of magnitude, and thus an EM-informed targeted search is able to find hints of a signal when an uninformed search fails to find any. Additionally, by combining EM and GW data, one can achieve an overall improvement in parameter estimation, regardless of the source’s sky location or GW frequency. We discuss the implications for the multi-messenger studies of SMBHBs withmore »PTAs.« less

    Supermassive black hole binaries (SMBHBs) are a natural outcome of galaxy mergers and should form frequently in galactic nuclei. Sub-parsec binaries can be identified from their bright electromagnetic emission, e.g. Active Galactic Nuclei (AGNs) with Doppler shifted broad emission lines or AGN with periodic variability, as well as from the emission of strong gravitational radiation. The most massive binaries (with total mass >108M⊙) emit in the nanohertz band and are targeted by Pulsar Timing Arrays (PTAs). Here we examine the synergy between electromagnetic and gravitational wave signatures of SMBHBs. We connect both signals to the orbital dynamics of the binary and examine the common link between them, laying the foundation for joint multimessenger observations. We find that periodic variability arising from relativistic Doppler boost is the most promising electromagnetic signature to connect with GWs. We delineate the parameter space (binary total mass/chirp mass versus binary period/GW frequency) for which joint observations are feasible. Currently multimessenger detections are possible only for the most massive and nearby galaxies, limited by the sensitivity of PTAs. However, we demonstrate that as PTAs collect more data in the upcoming years, the overlapping parameter space is expected to expand significantly.


    We searched for an isotropic stochastic gravitational wave background in the second data release of the International Pulsar Timing Array, a global collaboration synthesizing decadal-length pulsar-timing campaigns in North America, Europe, and Australia. In our reference search for a power-law strain spectrum of the form $h_c = A(f/1\, \mathrm{yr}^{-1})^{\alpha }$, we found strong evidence for a spectrally similar low-frequency stochastic process of amplitude $A = 3.8^{+6.3}_{-2.5}\times 10^{-15}$ and spectral index α = −0.5 ± 0.5, where the uncertainties represent 95 per cent credible regions, using information from the auto- and cross-correlation terms between the pulsars in the array. For a spectral index of α = −2/3, as expected from a population of inspiralling supermassive black hole binaries, the recovered amplitude is $A = 2.8^{+1.2}_{-0.8}\times 10^{-15}$. None the less, no significant evidence of the Hellings–Downs correlations that would indicate a gravitational-wave origin was found. We also analysed the constituent data from the individual pulsar timing arrays in a consistent way, and clearly demonstrate that the combined international data set is more sensitive. Furthermore, we demonstrate that this combined data set produces comparable constraints to recent single-array data sets which have more data than the constituent parts of the combination. Future international data releases will deliver increasedmore »sensitivity to gravitational wave radiation, and significantly increase the detection probability.

    « less
  4. Abstract We examine the light curves of two quasars, motivated by recent suggestions that a supermassive black hole binary (SMBHB) can exhibit sharp lensing spikes. We model the variability of each light curve as due to a combination of two relativistic effects: the orbital relativistic Doppler boost and gravitational binary self-lensing. In order to model each system we extend previous Doppler plus self-lensing models to include eccentricity. The first quasar is identified in optical data as a binary candidate with a 20-yr period (Ark 120), and shows a prominent spike. For this source, we rule out the lensing hypothesis and disfavor the Doppler-boost hypothesis due to discrepancies in the measured vs. recovered values of the binary mass and optical spectral slope. The second source, which we nickname Spikey, is the rare case of an active galactic nucleus (AGN) identified in Kepler’s high-quality, high-cadence photometric data. For this source, we find a model, consisting of a combination of Doppler modulation and a narrow symmetric lensing spike, consistent with an eccentric SMBHB with mass Mtot = 3 × 107M⊙, rest-frame orbital period T = 418 days, eccentricity e = 0.5, and seen at an inclination 8○ from edge-on. This interpretation can be testedmore »by monitoring Spikey for periodic behavior and recurring flares in the next few years. In preparation for such monitoring we present the first X-ray observations of this object taken by the Neil Gehrels Swift observatory.« less
  5. Abstract The nanohertz frequency band explored by pulsar timing arrays provides a unique discovery space for gravitational wave (GW) signals. In addition to signals from anticipated sources, such as those from supermassive black hole binaries, some previously unimagined sources may emit transient GWs (a.k.a. bursts) with unknown morphology. Unmodeled transients are not currently searched for in this frequency band, and they require different techniques from those currently employed. Possible sources of such GW bursts in the nanohertz regime are parabolic encounters of supermassive black holes, cosmic string cusps and kinks, or other, as-yet-unknown phenomena. In this paper we present BayesHopperBurst , a Bayesian search algorithm capable of identifying generic GW bursts by modeling both coherent and incoherent transients as a sum of Morlet–Gabor wavelets. A trans-dimensional reversible jump Markov chain Monte Carlo sampler is used to select the number of wavelets best describing the data. We test BayesHopperBurst on various simulated datasets including different combinations of signals and noise transients. Its capability to run on real data is demonstrated by analyzing data of the pulsar B1855 + 09 from the NANOGrav 9 year dataset. Based on a simulated dataset resembling the NANOGrav 12.5 year data release, we predict that atmore »our most sensitive time–frequency location we will be able to probe GW bursts with a root-sum-squared amplitude higher than ∼5 × 10 −11  Hz −1/2 , which corresponds to ∼40 M ⊙ c 2 emitted in GWs at a fiducial distance of 100 Mpc.« less