No Arabic abstract
We extract interstellar scintillation parameters for pulsars observed by the NANOGrav radio pulsar timing program. Dynamic spectra for the observing epochs of each pulsar were used to obtain estimates of scintillation timescales, scintillation bandwidths, and the corresponding scattering delays using a stretching algorithm to account for frequency-dependent scaling. We were able to measure scintillation bandwidths for 28 pulsars at 1500 MHz and 15 pulsars at 820 MHz. We examine scaling behavior for 17 pulsars and find power-law indices ranging from $-0.7$ to $-3.6$, though these may be biased shallow due to insufficient frequency resolution at lower frequencies. We were also able to measure scintillation timescales for six pulsars at 1500 MHz and seven pulsars at 820 MHz. There is fair agreement between our scattering delay measurements and electron-density model predictions for most pulsars. We derive interstellar scattering-based transverse velocities assuming isotropic scattering and a scattering screen halfway between the pulsar and earth. We also estimate the location of the scattering screens assuming proper motion and interstellar scattering-derived transverse velocities are equal. We find no correlations between variations in scattering delay and either variations in dispersion measure or flux density. For most pulsars for which scattering delays were measurable, we find that time of arrival uncertainties for a given epoch are larger than our scattering delay measurements, indicating that variable scattering delays are currently subdominant in our overall noise budget but are important for achieving precisions of tens of ns or less.
We present a new analysis of the profile data from the 47 millisecond pulsars comprising the 12.5-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), which is presented in a parallel paper (Alam et al. 2021a; NG12.5). Our reprocessing is performed using wideband timing methods, which use frequency-dependent template profiles, simultaneous time-of-arrival (TOA) and dispersion measure (DM) measurements from broadband observations, and novel analysis techniques. In particular, the wideband DM measurements are used to constrain the DM portion of the timing model. We compare the ensemble timing results to NG12.5 by examining the timing residuals, timing models, and noise model components. There is a remarkable level of agreement across all metrics considered. Our best-timed pulsars produce encouragingly similar results to those from NG12.5. In certain cases, such as high-DM pulsars with profile broadening, or sources that are weak and scintillating, wideband timing techniques prove to be beneficial, leading to more precise timing model parameters by 10-15%. The high-precision, multi-band measurements of several pulsars indicate frequency-dependent DMs. Compared to the narrowband analysis in NG12.5, the TOA volume is reduced by a factor of 33, which may ultimately facilitate computational speed-ups for complex pulsar timing array analyses. This first wideband pulsar timing data set is a stepping stone, and its consistent results with NG12.5 assure us that such data sets are appropriate for gravitational wave analyses.
We search for an isotropic stochastic gravitational-wave background (GWB) in the $12.5$-year 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. The Bayesian posterior of the amplitude for an $f^{-2/3}$ power-law spectrum, expressed as the characteristic GW strain, has median $1.92 times 10^{-15}$ and $5%$--$95%$ quantiles of $1.37$--$2.67 times 10^{-15}$ at a reference frequency of $f_mathrm{yr} = 1 ~mathrm{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.
We perform the first search for an isotropic non-tensorial gravitational-wave background (GWB) allowed in general metric theories of gravity in the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 12.5-year data set. By modeling the GWB as a power-law spectrum, we find strong Bayesian evidence for a spatially correlated process with scalar transverse (ST) correlations whose Bayes factor versus the spatially uncorrelated common-spectrum process is $99pm 7$, but no statistically significant evidence for the tensor transverse, vector longitudinal and scalar longitudinal polarization modes. The median and the $90%$ equal-tail amplitudes of ST mode are $mathcal{A}_{mathrm{ST}}= 1.06^{+0.35}_{-0.28} times 10^{-15}$, or equivalently the energy density parameter per logarithm frequency is $Omega_{mathrm{GW}}^{mathrm{ST}} = 1.54^{+1.20}_{-0.71} times 10^{-9}$, at frequency of 1/year.
We present time-of-arrival (TOA) measurements and timing models of 47 millisecond pulsars (MSPs) 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 semi-major 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 report flux density measurements for all pulsars observed. Searches for stochastic and continuous gravitational waves using these data will be subjects of forthcoming publications.
The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has recently reported evidence for the presence of a common stochastic signal across their array of pulsars. The origin of this signal is still unclear. One of the possibilities is that it is due to a stochastic gravitational wave background (SGWB) in the $sim 1-10,{rm nHz}$ frequency region. Taking the NANOGrav observational result at face value, we show that this signal would be fully consistent with a SGWB produced by an unresolved population of in-spiralling massive black hole binaries (MBHBs) predicted by current theoretical models. Considering an astrophysically agnostic model we find that the MBHB merger rate is loosely constrained to the range $10^{-11} - 2$ $mathrm{Mpc}^{-3},mathrm{Gyr}^{-1}$. Including additional constraints from galaxy pairing fractions and MBH-bulge scaling relations, we find that the MBHB merger rate is $10^{-5} - 5times10^{-4}$ $mathrm{Mpc}^{-3},mathrm{Gyr}^{-1}$, the MBHB merger time-scale is $le 3,mathrm{Gyr}$ and the norm of the $M_mathrm{BH}-M_mathrm{bulge}$ relation $ge 1.2times 10^{8},M_odot$ (all intervals quoted at 90% confidence). Regardless of the astrophysical details of MBHB assembly, this result would imply that a sufficiently large population of massive black holes pair up, form binaries and merge within a Hubble time.