Abstract
Evidence for a low-frequency stochastic gravitational-wave background has recently been reported based on analyses of pulsar timing array data. The most likely source of such a background is ...a population of supermassive black hole binaries, the loudest of which may be individually detected in these data sets. Here we present the search for individual supermassive black hole binaries in the NANOGrav 15 yr data set. We introduce several new techniques, which enhance the efficiency and modeling accuracy of the analysis. The search uncovered weak evidence for two candidate signals, one with a gravitational-wave frequency of ∼4 nHz, and another at ∼170 nHz. The significance of the low-frequency candidate was greatly diminished when Hellings–Downs correlations were included in the background model. The high-frequency candidate was discounted due to the lack of a plausible host galaxy, the unlikely astrophysical prior odds of finding such a source, and since most of its support comes from a single pulsar with a commensurate binary period. Finding no compelling evidence for signals from individual binary systems, we place upper limits on the strain amplitude of gravitational waves emitted by such systems. At our most sensitive frequency of 6 nHz, we place a sky-averaged 95% upper limit of 8 × 10
−15
on the strain amplitude. We also calculate an exclusion volume and a corresponding effective radius, within which we can rule out the presence of black hole binaries emitting at a given frequency.
The NANOGrav Collaboration reported strong Bayesian evidence for a common-spectrum stochastic process in its12.5 yr pulsar timing array data set, with median characteristic strain amplitude at ...periods of a year of A(yr) = 1.92(+0.75,-0.55) x 10^(-15). However, evidence for the quadrupolar Hellings & Downs interpulsar correlations, which are characteristic of gravitational-wave signals, was not yet significant. We emulate and extend the NANOGrav data set, injecting a wide range of stochastic gravitational-wave background(GWB)signals that encompass a variety of amplitudes and spectral shapes, and quantify three key milestones.(I)Given the amplitude measured in the 12.5 yr analysis and assuming this signal is a GWB, we expect to accumulate robust evidence of an interpulsar-correlated GWB signal with 15–17 yr of data, i.e., an additional 2–5 yr from the 12.5 yr data set.(II)At the initial detection, we expect a fractional uncertainty of 40% on the power-law strain spectrum slope, which is sufficient to distinguish a GWB of supermassive black hole binary origin from some models predicting more exotic origins.(III)Similarly, the measured GWB amplitude will have an uncertainty of 44% upon initial detection, allowing us to arbitrate between some population models of supermassive black hole binaries. In addition, power-law models are distinguishable from those having low-frequency spectral turnovers once 20 yr of data are reached. Even though our study is based on the NANOGrav data, we also derive relations that allow for a generalization to other pulsar timing array data sets. Most notably, by combining the data of individual arrays into the International Pulsar Timing Array, all of these milestones can be reached significantly earlier.
Abstract
The CHIME/FRB Project has recently released its first catalog of fast radio bursts (FRBs), containing 492 unique sources. We present results from angular cross-correlations of CHIME/FRB ...sources with galaxy catalogs. We find a statistically significant (
p
-value ∼ 10
−4
, accounting for look-elsewhere factors) cross-correlation between CHIME FRBs and galaxies in the redshift range 0.3 ≲
z
≲ 0.5, in three photometric galaxy surveys: WISE × SCOS, DESI-BGS, and DESI-LRG. The level of cross-correlation is consistent with an order-one fraction of the CHIME FRBs being in the same dark matter halos as survey galaxies in this redshift range. We find statistical evidence for a population of FRBs with large host dispersion measure (∼400 pc cm
−3
) and show that this can plausibly arise from gas in large halos (
M
∼ 10
14
M
⊙
), for FRBs near the halo center (
r
≲ 100 kpc). These results will improve in future CHIME/FRB catalogs, with more FRBs and better angular resolution.
Abstract Noise characterization for pulsar-timing applications accounts for interstellar dispersion by assuming a known frequency dependence of the delay it introduces in the times of arrival (TOAs). ...However, calculations of this delay suffer from misestimations due to other chromatic effects in the observations. The precision in modeling dispersion is dependent on the observed bandwidth. In this work, we calculate the offsets in infinite-frequency TOAs due to misestimations in the modeling of dispersion when using varying bandwidths at the Green Bank Telescope. We use a set of broadband observations of PSR J1643−1224, a pulsar with unusual chromatic timing behavior. We artificially restricted these observations to a narrowband frequency range, then used both the broad- and narrowband data sets to calculate residuals with a timing model that does not account for time variations in the dispersion. By fitting the resulting residuals to a dispersion model and comparing the fits, we quantify the error introduced in the timing parameters due to using a reduced frequency range. Moreover, by calculating the autocovariance function of the parameters, we obtained a characteristic timescale over which the dispersion misestimates are correlated. For PSR J1643−1224, which has one of the highest dispersion measures (DM) in the NANOGrav pulsar timing array, we find that the infinite-frequency TOAs suffer from a systematic offset of ∼22 μ s due to incomplete frequency sampling, with correlations over about one month. For lower-DM pulsars, the offset is ∼7 μ s. This error quantification can be used to provide more robust noise modeling in the NANOGrav data, thereby increasing the sensitivity and improving the parameter estimation in gravitational wave searches.
ABSTRACT
The Canadian Hydrogen Mapping Experiment (CHIME) is a radio telescope located in British Columbia, Canada. The large field of view allows CHIME/FRB to be an exceptional pulsar and rotating ...radio transient (RRAT) finding machine, despite saving only the metadata of incoming Galactic events. We have developed a pipeline to search for pulsar/RRAT candidates using density-based spatial clustering of applications with noise (dbscan), a clustering algorithm. Follow-up observations are then scheduled with the more sensitive CHIME/Pulsar instrument capable of near-daily high-time resolution spectra observations. We have developed the CHIME/Pulsar Single Pulse Pipeline to automate the processing of CHIME/Pulsar search-mode data. We report the discovery of 21 new Galactic sources, with 14 RRATs, 6 isolated long-period pulsars, and 1 binary system. Owing to CHIME/Pulsar’s observations, we have obtained timing solutions for 8 of the 14 RRATs along with all the regular pulsars and the binary system. Notably, we report that the binary system is in a long orbit of 412 d with a minimum companion mass of 0.1303 solar masses and no evidence of an optical companion within 10″ of the pulsar position. This highlights that working synergistically with CHIME/FRB’s large survey volume CHIME/Pulsar can obtain arc second localizations for low-burst rate RRATs through pulsar timing. We find that the properties of our newly discovered RRATs are consistent with those of the presently known population. They tend to have lower burst rates than those found in previous surveys, which is likely due to survey bias rather than the underlying population.
Abstract
The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has reported evidence for the presence of an isotropic nanohertz gravitational-wave background (GWB) in its 15 yr ...data set. However, if the GWB is produced by a population of inspiraling supermassive black hole binary (SMBHB) systems, then the background is predicted to be anisotropic, depending on the distribution of these systems in the local Universe and the statistical properties of the SMBHB population. In this work, we search for anisotropy in the GWB using multiple methods and bases to describe the distribution of the GWB power on the sky. We do not find significant evidence of anisotropy. By modeling the angular power distribution as a sum over spherical harmonics (where the coefficients are not bound to always generate positive power everywhere), we find that the Bayesian 95% upper limit on the level of dipole anisotropy is (
C
l
=1
/
C
l
=0
) < 27%. This is similar to the upper limit derived under the constraint of positive power everywhere, indicating that the dipole may be close to the data-informed regime. By contrast, the constraints on anisotropy at higher spherical-harmonic multipoles are strongly prior dominated. We also derive conservative estimates on the anisotropy expected from a random distribution of SMBHB systems using astrophysical simulations conditioned on the isotropic GWB inferred in the 15 yr data set and show that this data set has sufficient sensitivity to probe a large fraction of the predicted level of anisotropy. We end by highlighting the opportunities and challenges in searching for anisotropy in pulsar timing array data.
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
, and
h
0
< (2.66 ± 0.15) × 10
−15
in our most sensitive sky location. Finally, we placed a multimessenger limit of
<
(
1.41
±
0.02
)
×
10
9
M
⊙
on the chirp mass of the supermassive black hole binary candidate 3C 66B.
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.
Abstract
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.
Abstract
Analyses of pulsar timing data have provided evidence for a stochastic gravitational wave background in the nanohertz frequency band. The most plausible source of this background is the ...superposition of signals from millions of supermassive black hole binaries. The standard statistical techniques used to search for this background and assess its significance make several simplifying assumptions, namely (i) Gaussianity, (ii) isotropy, and most often, (iii) a power-law spectrum. However, a stochastic background from a finite collection of binaries does not exactly satisfy any of these assumptions. To understand the effect of these assumptions, we test standard analysis techniques on a large collection of realistic simulated data sets. The data-set length, observing schedule, and noise levels were chosen to emulate the NANOGrav 15 yr data set. Simulated signals from millions of binaries drawn from models based on the Illustris cosmological hydrodynamical simulation were added to the data. We find that the standard statistical methods perform remarkably well on these simulated data sets, even though their fundamental assumptions are not strictly met. They are able to achieve a confident detection of the background. However, even for a fixed set of astrophysical parameters, different realizations of the universe result in a large variance in the significance and recovered parameters of the background. We also find that the presence of loud individual binaries can bias the spectral recovery of the background if we do not account for them.