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 1014, 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.

The 15 yr pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) shows positive evidence for the presence of a low-frequency gravitational-wave (GW) background. In this paper, we investigate potential cosmological interpretations of this signal, specifically cosmic inflation, scalar-induced GWs, first-order phase transitions, cosmic strings, and domain walls. We find that, with the exception of stable cosmic strings of field theory origin, all these models can reproduce the observed signal. When compared to the standard interpretation in terms of inspiraling supermassive black hole binaries (SMBHBs), many cosmological models seem to provide a better fit resulting in Bayes factors in the range from 10 to 100. However, these results strongly depend on modeling assumptions about the cosmic SMBHB population and, at this stage, should not be regarded as evidence for new physics. Furthermore, we identify excluded parameter regions where the predicted GW signal from cosmological sources significantly exceeds the NANOGrav signal. These parameter constraints are independent of the origin of the NANOGrav signal and illustrate how pulsar timing data provide a new way to constrain the parameter space of these models. Finally, we search for deterministic signals produced by models of ultralight dark matter (ULDM) and dark matter substructures in the Milky Way. We find no evidence for either of these signals and thus report updated constraints on these models. In the case of ULDM, these constraints outperform torsion balance and atomic clock constraints for ULDM coupled to electrons, muons, or gluons.

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.

Pulsar timing arrays (PTAs) are galactic-scale gravitational wave (GW) detectors. Each individual arm, composed of a millisecond pulsar, a radio telescope, and a kiloparsecs-long path, differs in its properties but, in aggregate, can be used to extract low-frequency GW signals. We present a noise and sensitivity analysis to accompany the NANOGrav 15 yr data release and associated papers, along with an in-depth introduction to PTA noise models. As a first step in our analysis, we characterize each individual pulsar data set with three types of white-noise parameters and two red-noise parameters. These parameters, along with the timing model and, particularly, a piecewise-constant model for the time-variable dispersion measure, determine the sensitivity curve over the low-frequency GW band we are searching. We tabulate information for all of the pulsars in this data release and present some representative sensitivity curves. We then combine the individual pulsar sensitivities using a signal-to-noise ratio statistic to calculate the global sensitivity of the PTA to a stochastic background of GWs, obtaining a minimum noise characteristic strain of 7 × 10−15 at 5 nHz. A power-law-integrated analysis shows rough agreement with the amplitudes recovered in NANOGrav’s 15 yr GW background analysis. While our phenomenological noise model does not model all known physical effects explicitly, it provides an accurate characterization of the noise in the data while preserving sensitivity to multiple classes of GW signals.

We report the first discovery of a fast radio burst (FRB), FRB 20200125A, by the Green Bank Northern Celestial Cap (GBNCC) Pulsar Survey conducted with the Green Bank Telescope at 350 MHz. FRB 20200125A was detected at a Galactic latitude of 58°. 43 with a dispersion measure of 179 pc cm 3 , while electron density models predict a maximum Galactic contribution of 25 pc cm 3 along this line of sight. Moreover, no apparent Galactic foreground sources of ionized gas that could account for the excess DM are visible in multiwavelength surveys of this region. This argues that the source is extragalactic. The maximum redshift for the host galaxy is z max =0.17, corresponding to a maximum comoving distance of approximately 750 Mpc. The measured peak flux density for FRB 20200125A is 0.37 Jy, and we measure a pulse width of 3.7 ms, consistent with the distribution of FRB widths observed at higher frequencies. Based on this detection and assuming a Euclidean flux density distribution of FRBs, we calculate an all-sky rate at 350 MHz of-+ 3.4 10 3.3 15.43 FRBs sky −1 day −1 above a peak flux density of 0.42 Jy for an unscattered pulse having an intrinsic width of 5 ms, consistent with rates reported at higher frequencies, albeit with large uncertainties. Given the recent improvements in our single-pulse search pipeline, we also revisit the GBNCC survey sensitivity to various burst properties. Finally, we find no evidence of strong interstellar scattering in FRB 20200125A, adding to the growing evidence that some FRBs have circumburst environments where free-free absorption and scattering are not significant.

The Green Bank North Celestial Cap (GBNCC) pulsar survey will cover the entire northern sky (δ > −40°) at 350 MHz, and is one of the most uniform and sensitive all-sky pulsar surveys to date. We have created a pipeline to reanalyze GBNCC survey data to take a 350 MHz census of all pulsars detected by the survey, regardless of their discovery survey. Of the 1413 pulsars in the survey region, we were able to recover 670. For these we present measured signal-to-noise ratios (S/N), flux densities, pulse widths, profiles, and where appropriate, refined measurements of dispersion measures (DMs) (656 out of 670) and new or improved spectral indices (339 out of 670 total, 47 new, 292 improved). We also measure the period-pulse width relation at 350 MHz to scale as . Detection scans for several hundred sources were reanalyzed in order to inspect pulsars’ single pulse behavior and 223 were found to exhibit evidence of nulling. With a detailed analysis of measured and expected S/N values and the evolving radio frequency interference environment at 350 MHz, we assess the GBNCC survey’s sensitivity as a function of spin period, DM, and sky position. We find the sky-averaged limiting flux density of the survey to be 0.74 mJy. Combining this analysis with PsrPopPy pulsar population simulations, we predict 60/5 nonrecycled/MSP discoveries in the survey’s remaining 21,000 pointings, and we begin to place constraints on population model parameters.

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 Green Bank North Celestial Cap survey is a 350 MHz all-sky survey for pulsars and fast radio transients using the Robert C. Byrd Green Bank Telescope. To date, the survey has discovered over 190 pulsars, including 33 millisecond pulsars and 24 rotating radio transients. Several exotic pulsars have been discovered in the survey, including PSR J1759+5036, a binary pulsar with a 176 ms spin period in an orbit with a period of 2.04 days, an eccentricity of 0.3, and a projected semi-major axis of 6.8 light seconds. Using seven years of timing data, we are able to measure one post–Keplerian parameter, advance of periastron, which has allowed us to constrain the total system mass to 2.62 ± 0.03 M ⊙. This constraint, along with the spin period and orbital parameters, suggests that this is a double neutron star system, although we cannot entirely rule out a pulsar-white dwarf binary. This pulsar is only detectable in roughly 45% of observations, most likely due to scintillation. However, additional observations are required to determine whether there may be other contributing effects.