Pulsar timing arrays (PTAs) are searching for gravitational waves from supermassive black hole binaries (SMBHBs). Here we show how future PTAs could use a detection of gravitational waves from individually resolved SMBHB sources to produce a purely gravitational wave-based measurement of the Hubble constant. This is achieved by measuring two separate distances to the same source from the gravitational wave signal in the timing residual: the luminosity distance DL through frequency evolution effects, and the parallax distance Dpar through wavefront curvature (Fresnel) effects. We present a generalized timing residual model including these effects in an expanding universe. Of these two distances, Dpar is challenging to measure due to the pulsar distance wrapping problem, a degeneracy in the Earth-pulsar distance and gravitational wave source parameters that requires highly precise, sub-parsec level, pulsar distance measurements to overcome. However, in this paper we demonstrate that combining the knowledge of two SMBHB sources in the timing residual largely removes the wrapping cycle degeneracy. Two sources simultaneously calibrate the PTA by identifying the distances to the pulsars, which is useful in its own right, and allow recovery of the source luminosity and parallax distances which results in a measurement of the Hubble constant. We find that, with optimistic PTAs in the era of the Square Kilometre Array, two fortuitous SMBHB sources within a few hundred Mpc could be used to measure the Hubble constant with a relative uncertainty on the order of 10 per cent.
- Award ID(s):
- 1912649
- NSF-PAR ID:
- 10253440
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 505
- Issue:
- 3
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- 4531 to 4554
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
ABSTRACT -
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 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.more » « less
-
Abstract The radio galaxy 3C 66B has been hypothesized to host a supermassive black hole binary (SMBHB) at its center based on electromagnetic observations. Its apparent 1.05 yr period and low redshift (∼0.02) make it an interesting testbed to search for low-frequency gravitational waves (GWs) using pulsar timing array (PTA) experiments. This source has been subjected to multiple searches for continuous GWs from a circular SMBHB, resulting in progressively more stringent constraints on its GW amplitude and chirp mass. In this paper, we develop a pipeline for performing Bayesian targeted searches for eccentric SMBHBs in PTA data sets, and test its efficacy by applying it to simulated data sets with varying injected signal strengths. We also search for a realistic eccentric SMBHB source in 3C 66B using the NANOGrav 12.5 yr data set employing PTA signal models containing Earth term-only as well as Earth+pulsar term contributions using this pipeline. Due to limitations in our PTA signal model, we get meaningful results only when the initial eccentricity
e 0< 0.5 and the symmetric mass ratioη > 0.1. We find no evidence for an eccentric SMBHB signal in our data, and therefore place 95% upper limits on the PTA signal amplitude of 88.1 ± 3.7 ns for the Earth term-only and 81.74 ± 0.86 ns for the Earth+pulsar term searches fore 0< 0.5 andη > 0.1. Similar 95% upper limits on the chirp mass are (1.98 ± 0.05) × 109and (1.81 ± 0.01) × 109M ☉. These upper limits, while less stringent than those calculated from a circular binary search in the NANOGrav 12.5 yr data set, are consistent with the SMBHB model of 3C 66B developed from electromagnetic observations. -
Abstract Future searches for gravitational waves from space will be sensitive to double compact objects in our Milky Way. We present new simulations of the populations of double black holes (BHBHs), BH neutron stars (BHNSs), and double neutron stars (NSNSs) that will be detectable by the planned space-based gravitational-wave detector called Laser Interferometer Space Antenna (LISA). For our estimates, we use an empirically informed model of the metallicity-dependent star formation history of the Milky Way. We populate it using an extensive suite of binary population-synthesis predictions for varying assumptions relating to mass transfer, common-envelope, supernova kicks, remnant masses, and wind mass-loss physics. For a 4(10) yr LISA mission, we predict between 30–370(50–550) detections over these variations, out of which 6–154 (9–238) are BHBHs, 2–198 (3–289) are BHNSs, and 3–35 (4–57) are NSNSs. We expect that about 50% (60%) can be distinguished from double white dwarf sources based on their mass or eccentricity and localization. Specifically, for about 10% (15%), we expect to be able to determine chirp masses better than 10%. For 13% (13%), we expect sky-localizations better than 1°. We discuss how the variations in the physics assumptions alter the distribution of properties of the detectable systems, even when the detection rates are unchanged. We further discuss the possibility of multimessenger observations of pulsar populations with the Square Kilometre Array and assess the benefits of extending the LISA mission.
-
Abstract 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 NANOGrav’s recent 12.5 yr 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 nHz, we placed a sky-averaged limit of
h 0< (6.82 ± 0.35) × 10−15, andh 0< (2.66 ± 0.15) × 10−15in our most sensitive sky location. Finally, we placed a multimessenger limit of on the chirp mass of the supermassive black hole binary candidate 3C 66B.