Abstract
We search for an isotropic stochastic gravitational-wave background (GWB) in the 12.5 yr 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. Under our fiducial model, the Bayesian posterior of the amplitude for an
f
−2/3
power-law spectrum, expressed as the characteristic GW strain, has median 1.92 × 10
−15
and 5%–95% quantiles of 1.37–2.67 × 10
−15
at a reference frequency of
f
yr
=
1
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.
We present time-of-arrival (TOA) measurements and timing models of 47 millisecond pulsars observed from 2004 to 2017 at the Arecibo Observatory and the Green Bank Telescope by the North American ...Nanohertz Observatory for Gravitational Waves (NANOGrav). The observing cadence was three to four weeks for most pulsars over most of this time span, with weekly observations of six sources. These data were collected for use in low-frequency gravitational wave searches and for other astrophysical purposes. We detail our observational methods and present a set of TOA measurements, based on "narrowband" analysis, in which many TOAs are calculated within narrow radio-frequency bands for data collected simultaneously across a wide bandwidth. A separate set of "wideband" TOAs will be presented in a companion paper. We detail a number of methodological changes, compared to our previous work, which yield a cleaner and more uniformly processed data set. Our timing models include several new astrometric and binary pulsar measurements, including previously unpublished values for the parallaxes of PSRs J1832−0836 and J2322+2057, the secular derivatives of the projected semimajor orbital axes of PSRs J0613−0200 and J2229+2643, and the first detection of the Shapiro delay in PSR J2145−0750. We report detectable levels of red noise in the time series for 14 pulsars. As a check on timing model reliability, we investigate the stability of astrometric parameters across data sets of different lengths. We also report flux density measurements for all pulsars observed. Searches for stochastic and continuous gravitational waves using these data will be subjects of forthcoming publications.
We present a new analysis of the profile data from the 47 millisecond pulsars comprising the 12.5 yr data set of the North American Nanohertz Observatory for Gravitational Waves, which is presented ...in a parallel paper (Alam et al., hereafter NG12.5). Our reprocessing is performed using "wideband" timing methods, which use frequency-dependent template profiles, simultaneous time-of-arrival (TOA) and dispersion measure (DM) measurements from broadband observations, and novel analysis techniques. In particular, the wideband DM measurements are used to constrain the DM portion of the timing model. We compare the ensemble timing results to those in NG12.5 by examining the timing residuals, timing models, and noise-model components. There is a remarkable level of agreement across all metrics considered. Our best-timed pulsars produce encouragingly similar results to those from NG12.5. In certain cases, such as high-DM pulsars with profile broadening or sources that are weak and scintillating, wideband timing techniques prove to be beneficial, leading to more precise timing model parameters by 10%-15%. The high-precision, multiband measurements of several pulsars indicate frequency-dependent DMs. Compared to the narrowband analysis in NG12.5, the TOA volume is reduced by a factor of 33, which may ultimately facilitate computational speed-ups for complex pulsar timing array analyses. This first wideband pulsar timing data set is a stepping stone, and its consistent results with NG12.5 assure us that such data sets are appropriate for gravitational wave analyses.
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) × 109 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.
Abstract
We search NANOGrav’s 12.5 yr 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
γ
= 5 and a reference frequency of
f
yr
= 1 yr
−1
. Among the upper limits for eight general families of metric theories of gravity, we find the values of
A
TT
95
%
=
(
9.7
±
0.4
)
×
10
−
16
and
A
ST
95
%
=
(
1.4
±
0.03
)
×
10
−
15
for the family of metric spacetime theories that contain both TT and ST modes.
Pulsar timing arrays (PTAs) are uniquely poised to detect the nanohertz-frequency gravitational waves from supermassive black hole binaries (SMBHBs) formed during galaxy merger. Efforts are underway ...to observe three species of gravitational signal from these systems: the stochastic ensemble, individual, adiabatic binary inspirals, and bursts with memory. This dissertation discusses all three. A typical Bayesian search for evidence of a stochastic gravitational wave background from the superposition of many unresolvable SMBHB inspirals requires weeks to months to deliver results. This is due in part to the inclusion of inter-pulsar spatial and temporal correlations induced in PTA data by such a signal. By integrating a simplified Bayesian search into an existing frequentist statistic, we are able to create a robust background amplitude estimator that requires minimal CPU time and does not compromise the key information gleaned from a full Bayesian analysis. As PTA sensitivity increases, individual binary inspirals will rise above the stochastic background, promising information about local SMBHBs. PTAs, like the North American Nanohertz Observatory for Gravitational Waves, regularly conduct searches for these single-frequency sources, with the latest results coming from the 9- and 11-year data sets. Although detection is still in the future, results are already informing binary candidate properties. Finally, an SMBHB coalescence is theoretically accessible to PTAs through bursts with memory, a purely General Relativistic phenomenon which imparts a permanent spacetime deformation and affects the coalescence signal amplitude at leading quadrupole order. Simulations parameterized by astrophysical observables from galaxy mergers out to z=3 predict the rates and signal-to-noise ratios for bursts occurring in the PTA-band. Extending the synthesized population to include less massive SMBHBs shows space-based interferometers may also observe this atypical signature.
Supermassive black hole binaries (SMBHBs) should form frequently in galactic nuclei as a result of galaxy mergers. At subparsec separations, binaries become strong sources of low-frequency ...gravitational waves (GWs), targeted by Pulsar Timing Arrays. We used recent upper limits on continuous GWs from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11 yr data set 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 8 nHz, 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 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% that of the primary, roughly comparable to constraints on an SMBHB in the Milky Way. However, we demonstrated that the (typically large) uncertainties in the mass measurements can weaken the upper limits on the chirp mass. 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.
Gravitational wave memory is theorized to arise from the integrated history of gravitational wave emission, and manifests as a spacetime deformation in the wake of a propagating gravitational wave. ...We explore the detectability of the memory signals from a population of coalescencing supermassive black hole binaries with pulsar timing arrays and the Laser Interferometer Space Antenna (LISA). We find that current pulsar timing arrays have poor prospects, but it is likely that between 1 and 10 memory events with signal-to-noise ratio in excess of 5 will occur within LISA's planned 4-year mission.