Timing Pulsars with CHIME and the NANOGrav 15yr Dataset

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.

Ensemble Pulsar Time Scale

Pulsar timing arrays (PTAs) have revealed galaxy-size gravitational waves (GWs) in the form of a stochastic gravitational wave background (SGWB), correlating the radio pulses emitted by millisecond pulsars. This discovery naturally leads to the question of the origin and the nature of the SGWB; the latter is synonymous to testing how quadrupolar the inter-pulsar spatial correlation is. In this paper, we investigate the nature of the SGWB by considering correlations beyond the Hellings–Downs (HD) curve of Einstein’s general relativity. We scrutinize the HD and non-Einsteinian GW correlations with the North American Nanohertz Observatory for Gravitational Waves and the Chinese PTA data, and find that both data sets allow a graviton mass of mg ≲ 1.04 × 10−22 eV/c2 and subluminal traveling waves. We discuss gravitational physics scenarios beyond general relativity that could host non-Einsteinian GW correlations in the SGWB and highlight the importance of the cosmic variance inherited from stochastic variations across realizations in interpreting PTA observations.

We search for an isotropic stochastic gravitational-wave background (GWB) in the newly released 11 year data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). While we find no evidence for a GWB, we place constraints on a population of inspiraling supermassive black hole (SMBH) binaries, a network of decaying cosmic strings, and a primordial GWB. For the first time, we find that the GWB constraints are sensitive to the solar system ephemeris (SSE) model used and that SSE errors can mimic a GWB signal. We developed an approach that bridges systematic SSE differences, producing the first pulsar-timing array (PTA) constraints that are robust against SSE errors. We thus place a 95% upper limit on the GW-strain amplitude of A GWB < 1.45 × 10−15 at a frequency of f = 1 yr−1 for a fiducial f −2/3 power-law spectrum and with interpulsar correlations modeled. This is a factor of ∼2 improvement over the NANOGrav nine-year limit calculated using the same procedure. Previous PTA upper limits on the GWB (as well as their astrophysical and cosmological interpretations) will need revision in light of SSE systematic errors. We use our constraints to characterize the combined influence on the GWB of the stellar mass density in galactic cores, the eccentricity of SMBH binaries, and SMBH–galactic-bulge scaling relationships. We constrain the cosmic-string tension using recent simulations, yielding an SSE-marginalized 95% upper limit of Gμ < 5.3 × 10−11—a factor of ∼2 better than the published NANOGrav nine-year constraints. Our SSE-marginalized 95% upper limit on the energy density of a primordial GWB (for a radiation-dominated post-inflation universe) is ΩGWB(f) h 2 < 3.4 × 10−10.

We report the detection of 21 cm emission at an average redshift z¯=2.3 in the cross-correlation of data from the Canadian Hydrogen Intensity Mapping Experiment (CHIME) with measurements of the Lyα forest from eBOSS. Data collected by CHIME over 88 days in the 400–500 MHz frequency band (1.8 < z < 2.5) are formed into maps of the sky and high-pass delay filtered to suppress the foreground power, corresponding to removing cosmological scales with k ∥ ≲ 0.13 Mpc−1 at the average redshift. Line-of-sight spectra to the eBOSS background quasar locations are extracted from the CHIME maps and combined with the Lyα forest flux transmission spectra to estimate the 21 cm–Lyα cross-correlation function. Fitting a simulation-derived template function to this measurement results in a 9σ detection significance. The coherent accumulation of the signal through cross-correlation is sufficient to enable a detection despite excess variance from foreground residuals ∼6–10 times brighter than the expected thermal noise level in the correlation function. These results are the highest-redshift measurement of 21 cm emission to date, and they set the stage for future 21 cm intensity mapping analyses at z > 1.8.

Pulsar timing array observations have found evidence for an isotropic gravitational-wave background with the Hellings–Downs angular correlations between pulsar pairs. This interpretation hinges on the measured shape of the angular correlations, which is predominantly quadrupolar under general relativity. Here we explore a more flexible parameterization: we expand the angular correlations into a sum of Legendre polynomials and use a Bayesian analysis to constrain their coefficients with the 15 yr pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). When including Legendre polynomials with multipoles ℓ ≥ 2, we only find a significant signal in the quadrupole with an amplitude consistent with general relativity and nonzero at the ∼95% confidence level and a Bayes factor of 200. When we include multipoles ℓ ≤ 1, the Bayes factor evidence for quadrupole correlations decreases by more than an order of magnitude due to evidence for a monopolar signal at approximately 4 nHz, which has also been noted in previous analyses of the NANOGrav 15 yr data. Further work needs to be done in order to better characterize the properties of this monopolar signal and its effect on the evidence for quadrupolar angular correlations.

Observations indicate that nearly all galaxies contain supermassive black holes at their centers. When galaxies merge, their component black holes form SMBH binaries (SMBHBs), which emit low-frequency gravitational waves (GWs) that can be detected by pulsar timing arrays. We have searched the North American Nanohertz Observatory for Gravitational Waves 11 yr data set for GWs from individual SMBHBs in circular orbits. As we did not find strong evidence for GWs in our data, we placed 95% upper limits on the strength of GWs from such sources. At = 8 nHz, we placed a sky-averaged upper limit of h 0 &lt; 7.3(3) × 10−15. We also developed a technique to determine the significance of a particular signal in each pulsar using “dropout” parameters as a way of identifying spurious signals. From these upper limits, we ruled out SMBHBs emitting GWs with = 8 nHz within 120 Mpc for , and within 5.5 Gpc for at our most sensitive sky location. We also determined that there are no SMBHBs with emitting GWs with = 2.8–317.8 nHz in the Virgo Cluster. Finally, we compared our strain upper limits to simulated populations of SMBHBs, based on galaxies in the Two Micron All-Sky Survey and merger rates from the Illustris cosmological simulation project, and found that only 34 out of 75,000 realizations of the local universe contained a detectable source.

The sensitivity and field of view of the Canadian Hydrogen Intensity Mapping Experiment (CHIME) has enabled its fast radio burst (FRB) backend to detect thousands of FRBs. However, the low angular resolution of CHIME prevents it from localizing most FRBs to their host galaxies. Very long baseline interferometry (VLBI) can readily provide the subarcsecond resolution needed to localize many FRBs to their hosts. Thus, we developed TONE: an interferometric array of eight 6-m dishes to serve as a pathfinder for the CHIME/FRB Outriggers project, which will use wide field-of-view cylinders to determine the sky positions for a large sample of FRBs, revealing their positions within their host galaxies to subarcsecond precision. In the meantime, TONE’s [Formula: see text] [Formula: see text]km baseline with CHIME proves to be an excellent testbed for the development and characterization of single-pulse VLBI techniques at the time of discovery. This work describes the TONE instrument, its sensitivity, and its astrometric precision in single-pulse VLBI. We believe that our astrometric errors are dominated by uncertainties in the clock measurements which build up between successive Crab pulsar calibrations which happen every [Formula: see text] [Formula: see text]h; the wider fields-of-view and higher sensitivity of the Outriggers will provide opportunities for higher-cadence calibration. At present, CHIME-TONE localizations of the Crab pulsar yield systematic one-dimensional localization errors of [Formula: see text] [Formula: see text]asec — considerably better than the resolution achievable with seeing limited followup observations in the optical/near-infrared from the ground.

Fast Radio Bursts (FRBs) are classified into repeaters and non-repeaters, with only a few per cent of the observed FRB population from the Canadian Hydrogen Intensity Mapping Experiment (CHIME) confirmed as repeaters. However, this figure represents only a lower limit due to the observational biases, and the true fraction of repeaters remains unknown. Correcting for these biases uncovers a notable decline in apparently non-repeating FRB detection rate as the CHIME operational time increases. This finding suggests that a significant portion of apparently non-repeating FRBs could in fact exhibit repetition when observed over more extended periods. A simple population model infers that the true repeater fraction likely exceeds 50 per cent with 99 per cent confidence, a figure substantially larger than the observed face value, even consistent with 100 per cent. This greater prevalence of repeaters had previously gone unnoticed due to their very low repetition rates (∼10−3.5 h−1 on average). Hence, theoretical FRB models must incorporate these low-rate repeaters. Furthermore, our results indicate a significantly higher repeater volume number density, potentially exceeding observed values by up to 104 times, which in turn impacts comparisons with potential FRB progenitors.

Gravitational wave was directly detected for the first time in 2015 after a century of its prediction by Einstein’s General Theory of Relativity in 1916. But, the continuous, low frequency Gravitational-wave Background (GWB) that is thought to be present throughout the universe has not been detected directly so far. The researchers at North American Nanohertz Observatory for Gravitational Waves (NANOGrav) have recently reported detection of a low-frequency signal that could be ‘Gravitational-wave Background (GWB)’.

Primordial black holes (PBHs) can constitute the predominant fraction of dark matter (DM) if PBHs reside in the currently unconstrained “sublunar” mass range. PBHs originating from scalar perturbations generated during inflation can naturally appear with a broad spectrum in a class of models. The resulting stochastic gravitational wave (GW) background generated from such PBH production can account for the recently reported North American Nanohertz Observatory for Gravitational Waves (NANOGrav) pulsar timing array data signal, and will be testable in future GW observations by interferometer-type experiments such as Laser Interferometer Space Antenna (LISA). We show that the broad mass function of such PBH DM is already being probed by Subaru Hyper Suprime-Cam (HSC) microlensing data and is consistent with a detected candidate event. Upcoming observations of HSC will be able to provide an independent definitive test of the stochastic GW signals originating from such PBH DM production scenarios.

The observations from pulsar timing arrays (PTAs), led by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), have provided opportunities to constrain primordial gravitational waves at low frequencies. In this paper, we analyze the best-fit parameter values for different Gauss-Bonnet Inflationary Gravitational Wave (GB-IGW) models with the PTArcade program, and we compare the results with the observations of NANOGrav, European Pulsar Timing Array (EPTA), Parkes Pulsar Timing Array (PPTA), and International Pulsar Timing Array (IPTA). We find the potential parameter n derived from the GB-IGW model is not necessarily positive. Instead, PTA data suggest the possibility of new parameter ranges. Meanwhile, the other parameters' results, such as the spectral indices of the scalar and tensor perturbation also have differences with traditional Cosmic Microwave Background analyses from Planck 2018, showing the reason for the different phenomena of Gravitational Wave power spectra expected. Additionally, the GB-IGW models we calculated are closer to the central values of multiple PTA datasets compared to the standard inflation model, making the GB-IGW model more likely to be detected by future space-based gravitational wave observatories.

The North American Nanohertz Observatory for Gravitational Waves has recently reported strong evidence for a stochastic common-spectrum process affecting the pulsar timing residuals in its 12.5-year data set. We demonstrate that this process admits an interpretation in terms of a stochastic gravitational-wave background emitted by a cosmic-string network in the early Universe. We study stable Nambu-Goto strings in dependence of their tension Gμ and loop size α and show that the entire viable parameter space will be probed by an array of future experiments.

We present an algorithm for the simultaneous measurement of a pulse\ntime-of-arrival (TOA) and dispersion measure (DM) from folded wideband pulsar\ndata. We extend the prescription from Taylor (1992) to accommodate a general\ntwo-dimensional template "portrait", the alignment of which can be used to\nmeasure a pulse phase and DM. We show that there is a dedispersion reference\nfrequency that removes the covariance between these two quantities, and note\nthat the recovered pulse profile scaling amplitudes can provide useful\ninformation. We experiment with pulse modeling by using a Gaussian-component\nscheme that allows for independent component evolution with frequency, a\n"fiducial component", and the inclusion of scattering. We showcase the\nalgorithm using our publicly available code on three years of wideband data\nfrom the bright millisecond pulsar J1824-2452A (M28A) from the Green Bank\nTelescope, and a suite of Monte Carlo analyses validates the algorithm. By\nusing a simple model portrait of M28A we obtain DM trends comparable to those\nmeasured by more standard methods, with improved TOA and DM precisions by\nfactors of a few. Measurements from our algorithm will yield precisions at\nleast as good as those from traditional techniques, but is prone to fewer\nsystematic effects and is without ad hoc parameters. A broad application of\nthis new method for dispersion measure tracking with modern large-bandwidth\nobserving systems should improve the timing residuals for pulsar timing array\nexperiments, like the North American Nanohertz Observatory for Gravitational\nWaves.\n

High precision calibration is essential for the new generation of radio interferometers looking for Baryon Acoustic Oscillation signatures in neutral hydrogen emissions which is a signal buried in foregrounds. Yet, differences in instrument design and non-redundancies in the subsystems consisting such instruments can pose great challenges to proper calibration. For instance, Canadian Hydrogen Intensity Mapping Experiment (CHIME) offers a large instantaneous field of view with fast mapping speed with a sensitivity equivalent to a single dish in terms of total collecting area. However, the calibration of its 1000+ individual receivers poses major obstacle to harness the full capability of the new technology. In future instruments, it is planned to use redundancy at the level required for science experiments (typically below 1%). Although, this idea is against the traditional knowledge of having the radio receivers accurate to 1/20th of the observed wavelength, we plan to meet this goal by solving the redundancy issue in the front-end subsystems using precise metrology and alignment methods. To achieve this goal, we work on a 2-element interferometer called Deep Dish Development Array (D3A), with a targeted precision of 1/1000th of a wavelength at 300 MHz, located at Dominion Radio Astrophysical Observatory, located in Penticton, BC, Canada. The D3A will serve as test bed for dish prototyping, composite dish repeatability, antennae back ends which can be scaled to larger arrays without sacrificing the precision. The two dishes are made out of fiber glass composite and measured in the shop condition and in the field after fabrication. The relative pointing accuracy was obtained as 0.02°. The elevation axes of the dishes were placed in East-West direction within 0.04° of error. The surface RMS error was obtained as 0.389 mm which meets the 1/1000th of observation wavelength. RMS error due to temperature variation and assembly error on the dishes were obtained at 1 mm. This article presents the metrology principles applied to obtain the results and challenges. The tests performed in D3A will be implemented in Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) and the Canadian Hydrogen Observatory and Radio Transient Detector (CHORD).