skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: The Host Galaxies of Pulsar Timing Array Sources: Converting Supermassive Black Hole Binary Parameters into Electromagnetic Observables
Abstract Pulsar timing arrays (PTAs) are approaching the sensitivity required to resolve gravitational waves (GWs) from individual supermassive black hole (SMBH) binaries. However, the large uncertainty in source localization will make the identification of its host environment challenging. We show how to convert the posterior probability function of binary parameters inferred by GW analyses into distributions of apparent magnitudes of the host galaxy. We do so for a scenario in which the host environment is a regular early-type galaxy, and one in which it is an active galactic nucleus. We estimate the reach of PTAs in the near and intermediate future, and estimate whether the binary hosts will be detectable in all-sky electromagnetic (EM) surveys. A PTA with a baseline of 20 yr and 116 pulsars, resembling the upcoming data release of the International Pulsar Timing Array, can detect binaries out to a luminosity distance of 2 Gpc (corresponding to a redshift ofz ∼ 0.36), while a PTA with a baseline of 30 yr and 200 pulsars can reach out to distances slightly greater than 3 Gpc (z ∼ 0.53). We find that the host galaxies of all binaries detectable with a baseline of 20 yr are expected to be present in the Wide-field Infrared Survey Explorer and SuperCOSMOS surveys, if they lie outside the plane of the Milky Way. The Two Micron All Sky Survey becomes incomplete for hosts of binaries more massive than 109.8Mat a luminosity distance greater than 1 Gpc. The EM surveys become slightly more incomplete when PTAs with longer baselines and therefore improved sensitivities are considered.  more » « less
Award ID(s):
2307719
PAR ID:
10661523
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
IOP
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
990
Issue:
1
ISSN:
0004-637X
Page Range / eLocation ID:
46
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Pulsar timing arrays (PTAs) are Galactic-scale gravitational wave (GW) detectors consisting of precisely timed pulsars distributed across the sky. Within the decade, PTAs are expected to detect nanohertz GWs emitted by close-separation supermassive black hole binaries (SMBHBs), thereby opening up the low-frequency end of the GW spectrum for science. Individual SMBHBs which power active galactic nuclei are also promising multi-messenger sources; they may be identified via theoretically predicted electromagnetic (EM) signatures and be followed up by PTAs for GW observations. In this work, we study the detection and parameter estimation prospects of a PTA which targets EM-selected SMBHBs. Adopting a simulated Galactic millisecond pulsar population, we envisage three different pulsar timing campaigns which observe three mock sources at different sky locations. We find that an all-sky PTA which times the best pulsars is an optimal and feasible approach to observe EM-selected SMBHBs and measure their source parameters to high precision (i.e., comparable to or better than conventional EM measurements). We discuss the implications of our findings in the context of future PTA experiments with the planned Deep Synoptic Array-2000 and the multi-messenger studies of SMBHBs such as the well-known binary candidate OJ 287. 
    more » « less
  2. Abstract Supermassive black hole binaries (SMBHBs) are thought to form in galaxy mergers, possessing the potential to produce electromagnetic (EM) radiation as well as gravitational waves (GWs) detectable with pulsar timing arrays (PTAs). Once GWs from individually resolved SMBHBs are detected, the identification of the host galaxy will be a major challenge due to the ambiguity in possible EM signatures and the poor localization capability of PTAs. To aid EM observations in choosing follow-up sources, we use NANOGrav’s galaxy catalog to quantify the number of plausible hosts in both realistic and idealistic scenarios. We outline a host identification pipeline that injects a single-source GW signal into a simulated PTA data set, recovers the signal using production-level techniques, quantifies the localization region and number of galaxies contained therein, and finally imposes cuts on the galaxies using parameter estimates from the GW search. In an ideal case, the 90% credible areas span 29–241 deg2, containing about 14–341 galaxies. After cuts, the number of galaxies remaining ranges from 22 at worst to one true host at best. In a realistic case, these areas range from 287 to 530 deg2and enclose about 285–1238 galaxies. After cuts, the number of galaxies is 397 at worst and 27 at best. While the signal-to-noise ratio is the primary determinant of the localization area of a given source, we find that the area is also influenced by the proximity to nearby pulsars on the sky and the binary chirp mass. 
    more » « less
  3. Abstract The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational-wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTAs that constitute the International Pulsar Timing Array (IPTA). We show that despite making different modeling choices, there is no significant difference in the GWB parameters that are measured by the different PTAs, agreeing within 1σ. The pulsar noise parameters are also consistent between different PTAs for the majority of the pulsars included in these analyses. We bridge the differences in modeling choices by adopting a standardized noise model for all pulsars and PTAs, finding that under this model there is a reduction in the tension in the pulsar noise parameters. As part of this reanalysis, we “extended” each PTA’s data set by adding extra pulsars that were not timed by that PTA. Under these extensions, we find better constraints on the GWB amplitude and a higher signal-to-noise ratio for the Hellings–Downs correlations. These extensions serve as a prelude to the benefits offered by a full combination of data across all pulsars in the IPTA, i.e., the IPTA’s Data Release 3, which will involve not just adding in additional pulsars but also including data from all three PTAs where any given pulsar is timed by more than a single PTA. 
    more » « less
  4. Abstract Statistical anisotropy in the nanohertz-frequency gravitational wave background (GWB) is expected to be detected by pulsar timing arrays (PTAs) in the near future. By developing a frequentist statistical framework that intrinsically restricts the GWB power to be positive, we establish scaling relations for multipole-dependent anisotropy decision thresholds that are a function of the noise properties, timing baselines, and cadences of the pulsars in a PTA. We verify that (i) a larger number of pulsars, and (ii) factors that lead to lower uncertainty on spatial cross-correlation measurements between pulsars, lead to a higher overall GWB signal-to-noise ratio, and lower anisotropy decision thresholds with which to reject the null hypothesis of isotropy. Using conservative simulations of realistic NANOGrav data sets, we predict that an anisotropic GWB with angular power C l =1 > 0.3 C l =0 may be sufficient to produce tension with isotropy at the p = 3 × 10 −3 (∼3 σ ) level in near-future NANOGrav data with a 20 yr baseline. We present ready-to-use scaling relationships that can map these thresholds to any number of pulsars, configuration of pulsar noise properties, or sky coverage. We discuss how PTAs can improve the detection prospects for anisotropy, as well as how our methods can be adapted for more versatile searches. 
    more » « less
  5. 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 with PTAs. 
    more » « less