Search

The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational-wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTAs that constitute the International Pulsar Timing Array (IPTA). We show that despite making different modeling choices, there is no significant difference in the GWB parameters that are measured by the different PTAs, agreeing within 1 σ . The pulsar noise parameters are also consistent between different PTAs for the majority of the pulsars included in these analyses. We bridge the differences in modeling choices by adopting a standardized noise model for all pulsars and PTAs, finding that under this model there is a reduction in the tension in the pulsar noise parameters. As part of this reanalysis, we “extended” each PTA’s data set by adding extra pulsars that were not timed by that PTA. Under these extensions, we find better constraints on the GWB amplitude and a higher signal-to-noise ratio for the Hellings–Downs correlations. These extensions serve as a prelude to the benefits offered by a full combination of data across all pulsars in the IPTA, i.e., the IPTA’s Data Release 3, which will involve not just adding in additional pulsars but also including data from all three PTAs where any given pulsar is timed by more than a single PTA.

In this work, we present polarization profiles for 23 millisecond pulsars observed at 820 and 1500 MHz with the Green Bank Telescope as part of the NANOGrav pulsar timing array. We calibrate the data using Mueller matrix solutions calculated from observations of PSRs B1929+10 and J1022+1001. We discuss the polarization profiles, which can be used to constrain pulsar emission geometry, and present both the first published radio polarization profiles for nine pulsars and the discovery of very low-intensity average profile components (“microcomponents”) in four pulsars. We obtain the Faraday rotation measures for each pulsar and use them to calculate the Galactic magnetic field parallel to the line of sight for different lines of sight through the interstellar medium. We fit for linear and sinusoidal trends in time in the dispersion measure and Galactic magnetic field and detect magnetic field variations with a period of 1 yr in some pulsars, but overall find that the variations in these parameters are more consistent with a stochastic origin.

Low-frequency gravitational-wave experiments require the highest timing precision from an array of the most stable millisecond pulsars. Several known sources of noise on short timescales in single radio pulsar observations are well described by a simple model of three components: template fitting from a finite signal-to-noise ratio, pulse phase/amplitude jitter from single-pulse stochasticity, and scintillation errors from short-timescale interstellar scattering variations. Currently template-fitting errors dominate, but as radio telescopes push toward higher signal-to-noise ratios, jitter becomes the next dominant term for most millisecond pulsars. Understanding the statistics of jitter becomes crucial for properly characterizing arrival time uncertainties. We characterize the radio frequency dependence of jitter using data on 48 pulsars in the North American Nanohertz Observatory for Gravitational Waves timing program. We detect significant jitter in 43 of the pulsars and test several functional forms for its frequency dependence; we find significant frequency dependence for 30 pulsars. We find moderate correlations of rms jitter with pulse width (R = 0.62) and number of profile components (R = 0.40); the single-pulse rms jitter is typically ≈1% of pulse phase. The average frequency dependence for all pulsars using a power-law model has index −0.42. We investigate the jitter variations for the interpulse of PSR B1937+21 and find no significant deviations from the main pulse rms jitter. We also test the time variation of jitter in two pulsars and find that systematics likely bias the results for high-precision pulsars. Pulsar timing array analyses must properly model jitter as a significant component of the noise within the detector. [ABSTRACT FROM AUTHOR]

Using the Very Long Baseline Array, we have made astrometric observations of the binary pulsar B1913+16 spanning an 18-month period in 2014–2015. From these observations we make the first determination of the annual geometric parallax of B1913+16, measuring mas (68% confidence interval). The inferred parallax probability distribution differs significantly from a Gaussian. Using our parallax measurement and prior information on the spatial and luminosity distributions of the millisecond pulsar population, we infer a distance of kpc, which is significantly closer than the 9.8 ± 3.1 kpc suggested by the pulsar’s dispersion measure (DM) and analyses of the ionized interstellar medium. While the relatively low significance of the parallax detection (∼3σ) currently precludes an improved test of general relativity using the orbital decay of PSR B1913+16, ongoing observations with improved control of systematic astrometric errors could reach the 10% distance uncertainty required for this goal. The proper motion measured by our Very Long Baseline Interferometry astrometry differs substantially from that obtained by pulsar timing, a discrepancy that has also been found between the proper motion measurements made by interferometers and pulsar timing for some other pulsars, which we speculate is the result of timing noise or DM variations in the timing data set. Our parallax and proper motion measurements yield a transverse velocity of km s−1 in the solar reference frame. Analysis incorporating galactic rotation and solar motion finds that the space velocity of the pulsar relative to its standard of rest has a component km s−1 perpendicular to the galactic plane and components on the order of 100 km s−1 parallel to the galactic plane. [ABSTRACT FROM AUTHOR]

We analyze dispersion measure (DM) variations of 37 millisecond pulsars in the nine-year North American Nanohertz Observatory for Gravitational Waves (NANOGrav) data release and constrain the sources of these variations. DM variations can result from a changing distance between Earth and the pulsar, inhomogeneities in the interstellar medium, and solar effects. Variations are significant for nearly all pulsars, with characteristic timescales comparable to or even shorter than the average spacing between observations. Five pulsars have periodic annual variations, 14 pulsars have monotonically increasing or decreasing trends, and 14 pulsars show both effects. Of the four pulsars with linear trends that have line-of-sight velocity measurements, three are consistent with a changing distance and require an overdensity of free electrons local to the pulsar. Several pulsars show correlations between DM excesses and lines of sight that pass close to the Sun. Mapping of the DM variations as a function of the pulsar trajectory can identify localized interstellar medium features and, in one case, an upper limit to the size of the dispersing region of 4 au. Four pulsars show roughly Kolmogorov structure functions (SFs), and another four show SFs less steep than Kolmogorov. One pulsar has too large an uncertainty to allow comparisons. We discuss explanations for apparent departures from a Kolmogorov-like spectrum, and we show that the presence of other trends and localized features or gradients in the interstellar medium is the most likely cause. [ABSTRACT FROM AUTHOR]

Gravitational wave (GW) astronomy using a pulsar timing array requires high-quality millisecond pulsars (MSPs), correctable interstellar propagation delays, and high-precision measurements of pulse times of arrival. Here we identify noise in timing residuals that exceeds that predicted for arrival time estimation for MSPs observed by the North American Nanohertz Observatory for Gravitational Waves. We characterize the excess noise using variance and structure function analyses. We find that 26 out of 37 pulsars show inconsistencies with a white-noise-only model based on the short timescale analysis of each pulsar, and we demonstrate that the excess noise has a red power spectrum for 15 pulsars. We also decompose the excess noise into chromatic (radio-frequency-dependent) and achromatic components. Associating the achromatic red-noise component with spin noise and including additional power-spectrum-based estimates from the literature, we estimate a scaling law in terms of spin parameters (frequency and frequency derivative) and data-span length and compare it to the scaling law of Shannon & Cordes. We briefly discuss our results in terms of detection of GWs at nanohertz frequencies. [ABSTRACT FROM AUTHOR]

The use of pulsars as astrophysical clocks for gravitational wave (GW) experiments demands the highest possible timing precision. Pulse times of arrival (TOAs) are limited by stochastic processes that occur in the pulsar itself, along the line of sight through the interstellar medium, and in the measurement process. On timescales of seconds to hours, the TOA variance exceeds that from template-fitting errors due to additive noise. We assess contributions to the total variance from two additional effects: amplitude and phase jitter intrinsic to single pulses and changes in the interstellar impulse response from scattering. The three effects have different dependencies on time, frequency, and pulse signal-to-noise ratio. We use data on 37 pulsars from the North American Nanohertz Observatory for GWs to assess the individual contributions to the overall intraday noise budget for each pulsar. We detect jitter in 22 pulsars and estimate the average value of rms jitter in our pulsars to be of pulse phase. We examine how jitter evolves as a function of frequency and find evidence for evolution. Finally, we compare our measurements with previous noise parameter estimates and discuss methods to improve GW detection pipelines. [ABSTRACT FROM AUTHOR]

We report 21-year timing of one of the most precise pulsars: PSR J1713+0747. Its pulse times of arrival are well modeled by a comprehensive pulsar binary model including its three-dimensional orbit and a noise model that incorporates short- and long-timescale correlated noise such as jitter and red noise. Its timing residuals have weighted root mean square ∼92 ns. The new data set allows us to update and improve previous measurements of the system properties, including the masses of the neutron star (1.31 ± 0.11 M⊙) and the companion white dwarf (0.286 ± 0.012 M⊙) as well as their parallax distance 1.15 ± 0.03 kpc. We measured the intrinsic change in orbital period, , is −0.20 ± 0.17 ps s−1, which is not distinguishable from zero. This result, combined with the measured of other pulsars, can place a generic limit on potential changes in the gravitational constant G. We found that is consistent with zero [(−0.6 ± 1.1) × 10−12 yr−1, 95% confidence] and changes at least a factor of 31 (99.7% confidence) more slowly than the average expansion rate of the universe. This is the best limit from pulsar binary systems. The of pulsar binaries can also place limits on the putative coupling constant for dipole gravitational radiation (95% confidence). Finally, the nearly circular orbit of this pulsar binary allows us to constrain statistically the strong-field post-Newtonian parameters Δ, which describes the violation of strong equivalence principle, and , which describes a breaking of both Lorentz invariance in gravitation and conservation of momentum. We found, at 95% confidence, and based on PSR J1713+0747. [ABSTRACT FROM AUTHOR]