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.
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.
We report the discovery of an eclipsing binary millisecond pulsar in the globular cluster M92 (NGC 6341) with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). PSR J1717+4308A, or ...M92A, has a pulse frequency of 316.5 Hz (3.16 ms) and a dispersion measure of 35.45 pc cm−3. The pulsar is a member of a binary system with an orbital period of 0.20 days around a low-mass companion that has a median mass of ∼0.18 M . From observations so far, at least two eclipsing events have been observed in each orbit. The longer one lasted for ∼5000 s in the orbital phase range 0.1-0.5. The other lasted for ∼500 s and occurred between 1000 and 2000 s before or after the longer eclipsing event. The lengths of these two eclipsing events also change. These properties suggest that J1717+4308A is a "red-back" system with a low-mass main-sequence or sub-giant companion. Timing observations of the pulsar and further searches of the data for additional pulsars are ongoing.
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
We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15 yr pulsar timing data set collected by the North American Nanohertz Observatory ...for Gravitational Waves. The correlations follow the Hellings–Downs pattern expected for a stochastic gravitational-wave background. The presence of such a gravitational-wave background with a power-law spectrum is favored over a model with only independent pulsar noises with a Bayes factor in excess of 10
14
, and this same model is favored over an uncorrelated common power-law spectrum model with Bayes factors of 200–1000, depending on spectral modeling choices. We have built a statistical background distribution for the latter Bayes factors using a method that removes interpulsar correlations from our data set, finding
p
= 10
−3
(≈3
σ
) for the observed Bayes factors in the null no-correlation scenario. A frequentist test statistic built directly as a weighted sum of interpulsar correlations yields
p
= 5 × 10
−5
to 1.9 × 10
−4
(≈3.5
σ
–4
σ
). Assuming a fiducial
f
−2/3
characteristic strain spectrum, as appropriate for an ensemble of binary supermassive black hole inspirals, the strain amplitude is
2.4
−
0.6
+
0.7
×
10
−
15
(median + 90% credible interval) at a reference frequency of 1 yr
−1
. The inferred gravitational-wave background amplitude and spectrum are consistent with astrophysical expectations for a signal from a population of supermassive black hole binaries, although more exotic cosmological and astrophysical sources cannot be excluded. The observation of Hellings–Downs correlations points to the gravitational-wave origin of this signal.
Abstract
The NANOGrav 15 yr data set shows evidence for the presence of a low-frequency gravitational-wave background (GWB). While many physical processes can source such low-frequency gravitational ...waves, here we analyze the signal as coming from a population of supermassive black hole (SMBH) binaries distributed throughout the Universe. We show that astrophysically motivated models of SMBH binary populations are able to reproduce both the amplitude and shape of the observed low-frequency gravitational-wave spectrum. While multiple model variations are able to reproduce the GWB spectrum at our current measurement precision, our results highlight the importance of accurately modeling binary evolution for producing realistic GWB spectra. Additionally, while reasonable parameters are able to reproduce the 15 yr observations, the implied GWB amplitude necessitates either a large number of parameters to be at the edges of expected values or a small number of parameters to be notably different from standard expectations. While we are not yet able to definitively establish the origin of the inferred GWB signal, the consistency of the signal with astrophysical expectations offers a tantalizing prospect for confirming that SMBH binaries are able to form, reach subparsec separations, and eventually coalesce. As the significance grows over time, higher-order features of the GWB spectrum will definitively determine the nature of the GWB and allow for novel constraints on SMBH populations.
Abstract
We present observations and timing analyses of 68 millisecond pulsars (MSPs) comprising the 15 yr data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). ...NANOGrav is a pulsar timing array (PTA) experiment that is sensitive to low-frequency gravitational waves (GWs). This is NANOGrav’s fifth public data release, including both “narrowband” and “wideband” time-of-arrival (TOA) measurements and corresponding pulsar timing models. We have added 21 MSPs and extended our timing baselines by 3 yr, now spanning nearly 16 yr for some of our sources. The data were collected using the Arecibo Observatory, the Green Bank Telescope, and the Very Large Array between frequencies of 327 MHz and 3 GHz, with most sources observed approximately monthly. A number of notable methodological and procedural changes were made compared to our previous data sets. These improve the overall quality of the TOA data set and are part of the transition to new pulsar timing and PTA analysis software packages. For the first time, our data products are accompanied by a full suite of software to reproduce data reduction, analysis, and results. Our timing models include a variety of newly detected astrometric and binary pulsar parameters, including several significant improvements to pulsar mass constraints. We find that the time series of 23 pulsars contain detectable levels of red noise, 10 of which are new measurements. In this data set, we find evidence for a stochastic GW background.
Abstract
We report on the 2020 reactivation of the energetic high magnetic field pulsar PSR J1846–0258 and its pulsar wind nebula (PWN) after 14 yr of quiescence with new Chandra and Green Bank ...Telescope observations. The emission of short-duration bursts from J1846–0258 was accompanied by an enhancement of X-ray persistent flux and significant spectral softening, similar to those observed during its first bursting episode in 2006. The 2020 pulsar spectrum is described by a power-law model with a photon index Γ = 1.7 ± 0.3 in comparison to a Γ = 1.2 ± 0.1 before outburst, and shows evidence of an emerging thermal component with blackbody temperature
kT
= 0.7 ± 0.1 keV. The 0.5–10 keV unabsorbed flux increased from
× 10
−12
erg cm
−2
s
−1
in quiescence to
× 10
−11
erg cm
−2
s
−1
following the outburst. We did not detect any radio pulsations from the pulsar at 2 GHz, and we place an upper limit of 7.1
μ
Jy and 55 mJy for the coherent pulsed emission and single pulses, respectively. The 2020 PWN spectrum, characterized by a photon index of 1.92 ± 0.04 and X-ray luminosity of (1.2 ± 0.1) × 10
35
erg s
−1
at a distance of 5.8 kpc, is consistent with those observed before the outburst. An analysis of regions closer to the pulsar shows small-scale time variabilities and brightness changes over the 20 yr period from 2000 to 2020, while the photon indices did not change. We conclude that the outburst in PSR J1846–0258 is a combination of crustal and magnetospheric effects, with no significant burst-induced variability in its PWN based on the current observations.
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 search for a first-order phase transition gravitational wave signal in 45 pulsars from the NANOGrav 12.5-year dataset. We find that the data can be modeled in terms of a strong first order phase ...transition taking place at temperatures below the electroweak scale. However, we do not observe any strong preference for a phase-transition interpretation of the signal over the standard astrophysical interpretation in terms of supermassive black hole mergers; but we expect to gain additional discriminating power with future datasets, improving the signal to noise ratio and extending the sensitivity window to lower frequencies. An interesting open question is how well gravitational wave observatories could separate such signals.