Search

We present the host galaxies of four apparently nonrepeating fast radio bursts (FRBs), FRB 20181223C, FRB 20190418A, FRB 20191220A, and FRB 20190425A, reported in the first Canadian Hydrogen Intensity Mapping Experiment (CHIME/FRB) catalog. Our selection of these FRBs is based on a planned hypothesis testing framework where we search all CHIME/FRB Catalog-1 events that have low extragalactic dispersion measure ( 10°) and saved baseband data. We associate the selected FRBs with galaxies with moderate to high star formation rates located at redshifts between 0.027 and 0.071. We also search for possible multimessenger counterparts, including persistent compact radio and gravitational-wave sources, and find none. Utilizing the four FRB hosts from this study, along with the hosts of 14 published local Universe FRBs ( z < 0.1) with robust host association, we conduct an FRB host demographics analysis. We find all 18 local Universe FRB hosts in our sample to be spirals (or late-type galaxies), including the host of FRB 20220509G, which was previously reported to be elliptical. Using this observation, we scrutinize proposed FRB source formation channels and argue that core-collapse supernovae are likely the dominant channel to form FRB sources. Moreover, we infer no significant difference in the host properties of repeating and apparently nonrepeating FRBs in our local Universe FRB host sample. Finally, we find the burst rates of these four apparently nonrepeating FRBs to be consistent with those of the sample of localized repeating FRBs observed by CHIME/FRB. Therefore, we encourage further monitoring of these FRBs with more sensitive radio telescopes.

Recently we found compelling evidence for a gravitational-wave background with Hellings and Downs (HD) correlations in our 15 yr data set. These correlations describe gravitational waves as predicted by general relativity, which has two transverse polarization modes. However, more general metric theories of gravity can have additional polarization modes, which produce different interpulsar correlations. In this work, we search the NANOGrav 15 yr data set for evidence of a gravitational-wave background with quadrupolar HD and scalar-transverse (ST) correlations. We find that HD correlations are the best fit to the data and no significant evidence in favor of ST correlations. While Bayes factors show strong evidence for a correlated signal, the data does not strongly prefer either correlation signature, with Bayes factors ∼2 when comparing HD to ST correlations, and ∼1 for HD plus ST correlations to HD correlations alone. However, when modeled alongside HD correlations, the amplitude and spectral index posteriors for ST correlations are uninformative, with the HD process accounting for the vast majority of the total signal. Using the optimal statistic, a frequentist technique that focuses on the pulsar-pair cross-correlations, we find median signal-to-noise ratios of 5.0 for HD and 4.6 for ST correlations when fit for separately, and median signal-to-noise ratios of 3.5 for HD and 3.0 for ST correlations when fit for simultaneously. While the signal-to-noise ratios for each of the correlations are comparable, the estimated amplitude and spectral index for HD are a significantly better fit to the total signal, in agreement with our Bayesian analysis.

We report 10 fast radio bursts (FRBs) detected in the far sidelobe region (i.e., ≥5° off-meridian) of the Canadian Hydrogen Intensity Mapping Experiment (CHIME) from August 28 2018 to August 31 2021. We localize the bursts by fitting their spectra with a model of the CHIME/FRB synthesized beam response. We find that the far sidelobe events have on average ∼500 times greater fluxes than events detected in CHIME’s main lobe. We show that the sidelobe sample is therefore statistically ∼20 times closer than the main lobe sample. We find promising host galaxy candidates ( P _cc < 1%) for two of the FRBs, 20190112B and 20210310B, at distances of 38 and 16 Mpc, respectively. CHIME/FRB did not observe repetition of similar brightness from the uniform sample of 10 sidelobe FRBs in a total exposure time of 35,580 hr. Under the assumption of Poisson-distributed bursts, we infer that the mean repetition interval above the detection threshold of the far sidelobe events is longer than 11,880 hr, which is at least 2380 times larger than the interval from known CHIME/FRB detected repeating sources, with some caveats, notably that very narrowband events could have been missed. Our results from these far sidelobe events suggest one of two scenarios: either (1) all FRBs repeat and the repetition intervals span a wide range, with high-rate repeaters being a rare sub-population, or (2) non-repeating FRBs are a distinct population different from known repeaters.

Pulsar timing arrays (PTAs) are designed to detect low-frequency gravitational waves (GWs). GWs induce achromatic signals in PTA data, meaning that the timing delays do not depend on radio frequency. However, pulse arrival times are also affected by radio-frequency-dependent “chromatic” noise from sources such as dispersion measure (DM) and scattering delay variations. Furthermore, the characterization of GW signals may be influenced by the choice of chromatic noise model for each pulsar. To better understand this effect, we assess if and how different chromatic noise models affect the achromatic noise properties in each pulsar. The models we compare include existing DM models used by the North American Nanohertz Observatory for Gravitational waves (NANOGrav) and noise models used for the European PTA Data Release 2 (EPTA DR2). We perform this comparison using a subsample of six pulsars from the NANOGrav 15 yr data set, selecting the same six pulsars as from the EPTA DR2 six-pulsar data set. We find that the choice of chromatic noise model noticeably affects the achromatic noise properties of several pulsars. This is most dramatic for PSR J1713+0747, where the amplitude of its achromatic red noise lowers from ${\mathrm{log}}_{10}{A}_{\mathrm{RN}}=-{14.1}_{-0.1}^{+0.1}$ to $-{14.7}_{-0.5}^{+0.3}$ , and the spectral index broadens from ${\gamma }_{\mathrm{RN}}={2.6}_{-0.4}^{+0.5}$ to ${\gamma }_{\mathrm{RN}}={3.5}_{-0.9}^{+1.2}$ . We also compare each pulsar's noise properties with those inferred from the EPTA DR2, using the same models. From the discrepancies, we identify potential areas where the noise models could be improved. These results highlight the potential for custom chromatic noise models to improve PTA sensitivity to GWs.

Studying transient phenomena, such as individual pulses from pulsars, has garnered considerable attention in the era of astronomical big data. Of specific interest to this study are rotating radio transients (RRATs), nulling, and intermittent pulsars. This study introduces a new algorithm named LuNfit, tailored to correct the selection biases originating from the telescope and detection pipelines. Ultimately LuNfit estimates the intrinsic luminosity distribution and nulling fraction of the single pulses emitted by pulsars. LuNfit relies on Bayesian nested sampling so that the parameter space can be fully explored. Bayesian nested sampling also provides the additional benefit of simplifying model comparisons through the Bayes ratio. The robustness of LuNfit is shown through simulations and applying LuNfit onto pulsars with known nulling fractions. LuNfit is then applied to three RRATs, J0012+5431, J1538+1523, and J2355+1523, extracting their intrinsic luminosity distribution and burst rates. We find that their nulling fractions are 0.4(2), 0.749(5), and 0.995(2), respectively. We further find that a log-normal distribution likely describes the single pulse luminosity distribution of J0012+5431 and J1538+1523, while the Bayes ratio for J2355+1523 slightly favors an exponential distribution. We show the conventional method of correcting selection effects by “scaling up” the missed fraction of radio transients can be unreliable when the mean luminosity of the source is faint relative to the telescope sensitivity. Finally, we discuss the limitations of the current implementation of LuNfit while also delving into potential enhancements that would enable LuNfit to be applied to sources with complex pulse morphologies.

In 2021, a catalog of 536 fast radio bursts (FRBs) detected with the Canadian Hydrogen Intensity Mapping Experiment (CHIME) radio telescope was released by the CHIME/FRB Collaboration. This large collection of bursts, observed with a single instrument and uniform selection effects, has advanced our understanding of the FRB population. Here we update the results for 140 of these FRBs for which channelized raw voltage (“baseband”) data are available. With the voltages measured by the telescope’s antennas, it is possible to maximize the telescope sensitivity in any direction within the primary beam, an operation called “beamforming.” This allows us to increase the signal-to-noise ratios of the bursts and to localize them to subarcminute precision. The improved localizations are also used to correct the beam response of the instrument and to measure fluxes and fluences with an ∼10% uncertainty. Additionally, the time resolution is increased by 3 orders of magnitude relative to that in the first CHIME/FRB catalog, and, applying coherent dedispersion, burst morphologies can be studied in detail. Polarization information is also available for the full sample of 140 FRBs, providing an unprecedented data set to study the polarization properties of the population. We release the baseband data beamformed to the most probable position of each FRB. These data are analyzed in detail in a series of accompanying papers.

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.

We performed near-daily observations of the binary pulsars PSR J0218+4232, PSR J1518+4904, and PSR J2023+2853 with the Canadian Hydrogen Intensity Mapping Experiment (CHIME). For the first time, we detected the Shapiro time delay in all three pulsar binary systems, using only 2–4 yr of CHIME/Pulsar timing data. We measured the pulsar masses to be ${1.49}_{-0.20}^{+0.23}$ M _⊙ , ${1.470}_{-0.034}^{+0.030}$ M _⊙ , and ${1.50}_{-0.38}^{+0.49}$ M _⊙ , respectively. The companion mass to PSR J0218+4232 was found to be ${0.179}_{-0.016}^{+0.018}$ M _⊙ . We constrained the mass of the neutron star companion of PSR J1518+4904 to be ${1.248}_{-0.029}^{+0.035}$ M _⊙ , using the observed apsidal motion as a constraint on the mass estimation. The binary companion to PSR J2023+2853 was found to have a mass of ${0.93}_{-0.14}^{+0.17}$ M _⊙ ; in the context of the near-circular orbit, this mass estimate suggests that the companion to PSR J2023+2853 is likely a high-mass white dwarf. By comparing the timing model obtained for PSR J0218+4232 with previous observations, we found a significant change in the observed orbital period of the system of $\dot{{P}_{{\rm{b}}}}=0.14(2)\times {10}^{-12}$ ; we determined that this variation arises from “Shklovskii acceleration” due to the relative motion of the binary system, and used this measurement to estimate a distance of d = (6.7 ± 1.0) kpc to PSR J0218+4232. This work demonstrates the capability of high-cadence observations, enabled by the CHIME/Pulsar system, to detect and refine general-relativistic effects of binary pulsars over short observing timescales.

The millisecond pulsar J1713+0747 underwent a sudden and significant pulse shape change between 2021 April 16 and 17 (MJDs 59320 and 59321). Subsequently, the pulse shape gradually recovered over the course of several months. We report the results of continued multifrequency radio observations of the pulsar made using the Canadian Hydrogen Intensity Mapping Experiment and the 100 m Green Bank Telescope in a 3 yr period encompassing the shape change event, between 2020 February and 2023 February. As of 2023 February, the pulse shape had returned to a state similar to that seen before the event, but with measurable changes remaining. The amplitude of the shape change and the accompanying time-of-arrival residuals display a strong nonmonotonic dependence on radio frequency, demonstrating that the event is neither a glitch (the effects of which should be independent of radio frequency, ν ) nor a change in dispersion measure alone (which would produce a delay proportional to ν ^−2 ). However, it does bear some resemblance to the two previous “chromatic timing events” observed in J1713+0747, as well as to a similar event observed in PSR J1643−1224 in 2015.

We report on contemporaneous optical observations at ≈10 ms timescales from the fast radio burst (FRB) 20180916B of two repeat bursts (FRB 20201023 and FRB 20220908) taken with the ‘Alopeke camera on the Gemini-North telescope. These repeats have radio fluences of 2.8 and 3.5 Jy ms, respectively, approximately in the lower 50th percentile for fluence from this repeating burst. The ‘Alopeke data reveal no significant optical detections at the FRB position and we place 3 σ upper limits to the optical fluences of

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 for e _0 < 0.5 and η > 0.1. Similar 95% upper limits on the chirp mass are (1.98 ± 0.05) × 10 ^9 and (1.81 ± 0.01) × 10 ^9 M _☉ . 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.

We present the results of a Bayesian search for gravitational wave (GW) memory in the NANOGrav 12.5 yr data set. We find no convincing evidence for any gravitational wave memory signals in this data set. We find a Bayes factor of 2.8 in favor of a model that includes a memory signal and common spatially uncorrelated red noise (CURN) compared to a model including only a CURN. However, further investigation shows that a disproportionate amount of support for the memory signal comes from three dubious pulsars. Using a more flexible red-noise model in these pulsars reduces the Bayes factor to 1.3. Having found no compelling evidence, we go on to place upper limits on the strain amplitude of GW memory events as a function of sky location and event epoch. These upper limits are computed using a signal model that assumes the existence of a common, spatially uncorrelated red noise in addition to a GW memory signal. The median strain upper limit as a function of sky position is approximately 3.3 × 10 ^−14 . We also find that there are some differences in the upper limits as a function of sky position centered around PSR J0613−0200. This suggests that this pulsar has some excess noise that can be confounded with GW memory. Finally, the upper limits as a function of burst epoch continue to improve at later epochs. This improvement is attributable to the continued growth of the pulsar timing array.

We present results from angular cross correlations between select samples of CHIME/FRB repeaters and galaxies in three photometric galaxy surveys, which have shown correlations with the first CHIME/FRB catalog containing repeating and nonrepeating sources: Wide-field Infrared Survey Explorer (WISE) × SCOS, DESI-BGS, and DESI-LRG. We find a statistically significant correlation ( p -value 395 pc cm ^−3 and WISE × SCOS galaxies with redshift z > 0.275. We demonstrate that the correlation arises surprisingly because of a statistical association between FRB 20200320A (extragalactic DM ≈ 550 pc cm ^−3 ) and a galaxy group in the same dark matter halo at redshift z ≈ 0.32. We estimate that the host halo, along with an intervening halo at redshift z ≈ 0.12, accounts for at least ∼30% of the extragalactic DM. Our results strongly motivate incorporating galaxy group and cluster catalogs into direct host association pipelines for FRBs with $\lesssim 1^{\prime} $ localization precision, effectively utilizing the two-point information to constrain FRB properties such as their redshift and DM distributions. In addition, we find marginal evidence for a negative correlation at 99.4% confidence limit between a sample of repeating FRBs with baseband data (median extragalactic DM = 354 pc cm ^−3 ) and DESI-LRG galaxies with redshift 0.3 ≤ z < 0.45, suggesting that the repeaters might be more prone than apparent nonrepeaters to propagation effects in FRB–galaxy correlations due to intervening free electrons over angular scales ∼0.°5.

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.

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.

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.

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 ${ \mathcal M }\lt (1.41\pm 0.02)\times {10}^{9}\,{M}_{\odot }$ on the chirp mass of the supermassive black hole binary candidate 3C 66B.

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}\times {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.

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.

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.

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.

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.

FRB 20180916B is a repeating fast radio burst (FRB) with a 16.3 day periodicity in its activity. In this study, we present morphological properties of 60 FRB 20180916B bursts detected by CHIME/FRB between 2018 August and 2021 December. We recorded raw voltage data for 45 of these bursts, enabling microseconds time resolution in some cases. We studied variation of spectro-temporal properties with time and activity phase. We find that the variation in dispersion measure (DM) is ≲1 pc cm ^−3 and that there is burst-to-burst variation in scattering time estimates ranging from ∼0.16 to over 2 ms, with no discernible trend with activity phase for either property. Furthermore, we find no DM and scattering variability corresponding to the recent change in rotation measure from the source, which has implications for the immediate environment of the source. We find that FRB 20180916B has thus far shown no epochs of heightened activity as have been seen in other active repeaters by CHIME/FRB, with its burst count consistent with originating from a Poissonian process. We also observe no change in the value of the activity period over the duration of our observations and set a 1 σ upper limit of 1.5 × 10 ^−4 day day ^−1 on the absolute period derivative. Finally, we discuss constraints on progenitor models yielded by our results, noting that our upper limits on changes in scattering and DM as a function of phase do not support models invoking a massive binary companion star as the origin of the 16.3 day periodicity.

Fast radio bursts (FRBs) are a class of highly energetic, mostly extragalactic radio transients lasting for ∼milliseconds. While over 600 FRBs have been published so far, their origins are presently unclear, with some theories for extragalactic FRBs predicting accompanying high-energy emission. In this work, we use the Canadian Hydrogen Intensity Mapping Experiment (CHIME) Fast Radio Burst (CHIME/FRB) Project to explore whether any FRB-like radio emission coincides in space and time with 81 gamma-ray bursts (GRBs) detected between 2018 July 17 and 2019 July 8 by Swift/BAT and Fermi/GBM. We do not find any statistically significant coincident pairs within 3 σ of each other’s spatial localization regions and within a time difference of up to one week. In addition to searching for spatial matches between known FRBs and known GRBs, we use CHIME/FRB to constrain FRB-like (∼1–10 ms) radio emission before, at the time of, or after the reported high-energy emission at the position of 39 GRBs. For short gamma-ray bursts (SGRBs), we constrain the radio flux in the 400–800 MHz band to be under a few kJy for ∼5.5–12.5 hr post-high-energy burst. We use these limits to constrain models that predict FRB-like prompt radio emission after SGRBs. For long gamma-ray bursts, we constrain the radio flux to be under a few kJy from ∼6 hr pre-high-energy burst to ∼12 hr post-high-energy burst.

We report on improved sky localizations of 13 repeating fast radio bursts (FRBs) discovered by CHIME/FRB via the use of interferometric techniques on channelized voltages from the telescope. These so-called “baseband localizations” improve the localization uncertainty area presented in past studies by more than three orders of magnitude. The improved localization regions are provided for the full sample of FRBs to enable follow-up studies. The localization uncertainties, together with the limits on the source distances from their dispersion measures, allow us to identify likely host galaxies for two of the FRB sources. FRB 20180814A lives in a massive passive red spiral at z ∼ 0.068 with very little indication of star formation, while FRB 20190303A resides in a merging pair of spiral galaxies at z ∼ 0.064 undergoing significant star formation. These galaxies show very different characteristics, further confirming the presence of FRB progenitors in a variety of environments even among the repeating subclass.

We discover three new weak pulse components in two known pulsars, one in PSR J0304+1932 and two in PSR J1518+4904. These components are emitted about halfway between the main emission beam and the interpulse beam (beam from the opposite pole). They are separated from their main pulse peak by 99° ± 3° for J0304+1932 and 123.°6 ± 0.°7 (leading) and 93° ± 0.°4 (trailing) for J1518+4904. Their peak-intensity ratios to main pulses are ∼ 0.06% for J0304+1932 and ∼0.17% and ∼0.83% for J1518+4904. We also analyzed the flux fluctuations and profile variations of the emissions for the two pulsars. The results show correlations between the weak pulses and their main pulses, indicating that these emissions come from the same pole. We estimated the emission altitude of these weak pulses and derived a height of about half of the pulsar’s light-cylinder radius. These pulse components are a unique sample of high-altitude emissions from pulsars, and challenge the current pulsar emission models.

We present the discovery of 25 new repeating fast radio burst (FRB) sources found among CHIME/FRB events detected between 2019 September 30 and 2021 May 1. The sources were found using a new clustering algorithm that looks for multiple events colocated on the sky having similar dispersion measures (DMs). The new repeaters have DMs ranging from ∼220 to ∼1700 pc cm ^−3 , and include sources having exhibited as few as two bursts to as many as twelve. We report a statistically significant difference in both the DM and extragalactic DM (eDM) distributions between repeating and apparently nonrepeating sources, with repeaters having a lower mean DM and eDM, and we discuss the implications. We find no clear bimodality between the repetition rates of repeaters and upper limits on repetition from apparently nonrepeating sources after correcting for sensitivity and exposure effects, although some active repeating sources stand out as anomalous. We measure the repeater fraction over time and find that it tends to an equilibrium of ${2.6}_{-2.6}^{+2.9}$ % over our total time-on-sky thus far. We also report on 14 more sources, which are promising repeating FRB candidates and which merit follow-up observations for confirmation.

Of the more than 3000 radio pulsars currently known, only ∼300 are in binary systems, and only five of these consist of young pulsars with massive nondegenerate companions. We present the discovery and initial timing, accomplished using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope, of the sixth such binary pulsar, PSR J2108+4516, a 0.577 s radio pulsar in a 269 day orbit of eccentricity 0.09 with a companion of minimum mass 11 M _⊙ . Notably, the pulsar undergoes periods of substantial eclipse, disappearing from the CHIME 400–800 MHz observing band for a large fraction of its orbit, and displays significant dispersion measure and scattering variations throughout its orbit, pointing to the possibility of a circumstellar disk or very dense stellar wind associated with the companion star. Subarcsecond resolution imaging with the Karl G. Jansky Very Large Array unambiguously demonstrates that the companion is a bright, V ≃ 11 OBe star, EM* UHA 138, located at a distance of 3.26(14) kpc. Archival optical observations of EM* UHA 138 approximately suggest a companion mass ranging from 17.5 M _⊙ < M _c < 23 M _⊙ , in turn constraining the orbital inclination angle to 50.°3 ≲ i ≲ 58.°3. With further multiwavelength follow-up, PSR J2108+4516 promises to serve as another rare laboratory for the exploration of companion winds, circumstellar disks, and short-term evolution through extended-body orbital dynamics.

Fast radio bursts (FRBs) are brief, energetic, typically extragalactic flashes of radio emission whose progenitors are largely unknown. Although studying the FRB population is essential for understanding how these astrophysical phenomena occur, such studies have been difficult to conduct without large numbers of FRBs and characterizable observational biases. Using the recently released catalog of 536 FRBs published by the Canadian Hydrogen Intensity Mapping Experiment/Fast Radio Burst (CHIME/FRB) collaboration, we present a study of the FRB population that also calibrates for selection effects. Assuming a Schechter function, we infer a characteristic energy cut-off of ${E}_{\mathrm{char}}={2.38}_{-1.64}^{+5.35}\times {10}^{41}$ erg and a differential power-law index of γ = $-{1.3}_{-0.4}^{+0.7}$ . Simultaneously, we infer a volumetric rate of [ ${7.3}_{-3.8}^{+8.8}$ (stat.) ${}_{-1.8}^{+2.0}(\mathrm{sys}.)]\times {10}^{4}$ Gpc ^−3 yr ^−1 above a pivot energy of 10 ^39 erg and below a scattering timescale of 10 ms at 600 MHz, and find we cannot significantly constrain the cosmic evolution of the FRB population with star-formation rate. Modeling the host’s dispersion measure (DM) contribution as a log-normal distribution and assuming a total Galactic contribution of 80 pc cm ^−3 , we find a median value of ${\mathrm{DM}}_{\mathrm{host}}={84}_{-49}^{+69}$ pc cm ^−3 , comparable with values typically used in the literature. Proposed models for FRB progenitors should be consistent with the energetics and abundances of the full FRB population predicted by our results. Finally, we infer the redshift distribution of FRBs detected with CHIME, which will be tested with the localizations and redshifts enabled by the upcoming CHIME/FRB Outriggers project.

In this work, we present follow-up observations of two known repeating fast radio bursts (FRBs) and seven nonrepeating FRBs with complex morphology discovered with the CHIME/FRB project. These observations were conducted with the Arecibo Observatory 327 MHz receiver. We detected no additional bursts from these sources, nor did CHIME/FRB detect any additional bursts from these sources during our follow-up program. Based on these nondetections, we provide constraints on the repetition rate for all nine sources. We calculate repetition rates above 1 Jy using both a Poisson distribution of repetition and the Weibull distribution of repetition presented by Oppermann et al. For both distributions we find repetition upper limits of the order λ = 10 ^−2 –10 ^−1 hr ^−1 for all sources. These rates are much lower than those published for notable repeating FRBs like FRB 20121102A and FRB 20201124A, suggesting the possibility of a low-repetition subpopulation.

The CHIME/FRB project has detected hundreds of fast radio bursts (FRBs), providing an unparalleled population to statistically probe the foreground media that they illuminate. One such foreground medium is the ionized halo of the Milky Way (MW). We estimate the total Galactic electron column density from FRB dispersion measures (DMs) as a function of Galactic latitude using four different estimators, including ones that assume spherical symmetry of the ionized MW halo and ones that imply more latitudinal variation in density. Our observation-based constraints of the total Galactic DM contribution for ∣ b ∣ ≥ 30°, depending on the Galactic latitude and selected model, span 87.8–141 pc cm ^−3 . This constraint implies upper limits on the MW halo DM contribution that range over 52–111 pc cm ^−3 . We discuss the viability of various gas density profiles for the MW halo that have been used to estimate the halo’s contribution to DMs of extragalactic sources. Several models overestimate the DM contribution, especially when assuming higher halo gas masses (∼3.5 × 10 ^12 M _⊙ ). Some halo models predict a higher MW halo DM contribution than can be supported by our observations unless the effect of feedback is increased within them, highlighting the impact of feedback processes in galaxy formation.