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.

Continued timing observations of the double pulsar PSR J0737–3039A/B, which consists of two active radio pulsars (A and B) that orbit each other with a period of 2.45 h in a mildly eccentric (e=0.088) binary system, have led to large improvements in the measurement of relativistic effects in this system. With a 16-yr data span, the results enable precision tests of theories of gravity for strongly self-gravitating bodies and also reveal new relativistic effects that have been expected but are now observed for the first time. These include effects of light propagation in strong gravitational fields which are currently not testable by any other method. In particular, we observe the effects of retardation and aberrational light bending that allow determination of the spin direction of the pulsar. In total, we detect seven post-Keplerian parameters in this system, more than for any other known binary pulsar. For some of these effects, the measurement precision is now so high that for the first time we have to take higher-order contributions into account. These include the contribution of the A pulsar’s effective mass loss (due to spin-down) to the observed orbital period decay, a relativistic deformation of the orbit, and the effects of the equation of state of superdense matter on the observed post-Keplerian parameters via relativistic spin-orbit coupling. We discuss the implications of our findings, including those for the moment of inertia of neutron stars, and present the currently most precise test of general relativity’s quadrupolar description of gravitational waves, validating the prediction of general relativity at a level of 1.3×10^{-4} with 95% confidence. We demonstrate the utility of the double pulsar for tests of alternative theories of gravity by focusing on two specific examples and also discuss some implications of the observations for studies of the interstellar medium and models for the formation of the double pulsar system. Finally, we provide context to other types of related experiments and prospects for the future.

We report on eight millisecond pulsars (MSPs) in binary systems discovered with the Arecibo L-Band Feed Array (PALFA) pulsar survey. Phase-coherent timing solutions derived from 2.5–5 yr of observations carried out at the Arecibo and Jodrell Bank observatories are provided. PSR J1921+1929 is a 2.65 ms pulsar in a 39.6 day orbit for which we detect γ-ray pulsations in archival Fermi data. PSR J1928+1245 is a very low-mass-function system with an orbital period of 3.3 hr that belongs to the non-eclipsing black widow population. We also present PSR J1932+1756, the longest-orbital-period (41.5 days) intermediate-mass binary pulsar known to date. In light of the numerous discoveries of binary MSPs over the past years, we characterize the Galactic distribution of known MSP binaries in terms of binary class. Our results support and strengthen previous claims that the scatter in the Galactic scale height distribution correlates inversely with the binary mass function. We provide evidence of observational biases against detecting the most recycled pulsars near the Galactic plane, which overestimates the scale height of lighter systems. A possible bimodality in the mass function of MSPs with massive white dwarfs is also reported. [ABSTRACT FROM AUTHOR]

In this paper, we present the results of timing observations of PSRs J1949+3106 and J1950+2414, two binary millisecond pulsars (MSPs) discovered in data from the Arecibo ALFA pulsar survey (PALFA). The timing parameters include precise measurements of the proper motions of both pulsars, which show that PSR J1949+3106 has a transversal motion very similar to that of an object in the local standard of rest. The timing also includes measurements of the Shapiro delay and the rate of advance of periastron for both systems. Assuming general relativity, these allow estimates of the masses of the components of the two systems; for PSR J1949+3106, the pulsar mass is and the companion mass for PSR J1950+2414 and (all values have 68.3% confidence limits). We use these masses and proper motions to investigate the evolutionary history of both systems: PSR J1949+3106 is likely the product of a low-kick supernova; PSR J1950+2414 is a member of a new class of eccentric MSP binaries with an unknown formation mechanism. We discuss the proposed hypotheses for the formations of these systems in light of our new mass measurements. [ABSTRACT FROM AUTHOR]

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 h0 < 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. [ABSTRACT FROM AUTHOR]

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has observed dozens of millisecond pulsars for over a decade. We have accrued a large collection of dispersion measure (DM) measurements sensitive to the total electron content between Earth and the pulsars at each observation. All lines of sight cross through the solar wind (SW), which produces correlated DM fluctuations in all pulsars. We develop and apply techniques for extracting the imprint of the SW from the full collection of DM measurements in the recently released NANOGrav 11 yr data set. We filter out long-timescale DM fluctuations attributable to structure in the interstellar medium and carry out a simultaneous analysis of all pulsars in our sample that can differentiate the correlated signature of the wind from signals unique to individual lines of sight. When treating the SW as spherically symmetric and constant in time, we find the electron number density at 1 au to be 7.9 ± 0.2 cm−3. We find our data to be insensitive to long-term variation in the density of the wind. We argue that our techniques paired with a high-cadence, low-radio-frequency observing campaign of near-ecliptic pulsars would be capable of mapping out large-scale latitudinal structure in the wind. [ABSTRACT FROM AUTHOR]

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]

Access to 50 years of data has led to the discovery of pulsar emission and rotation variability on timescales of months and years. Most of this long-term variability has been seen in long-period pulsars, with relatively little focus on recycled millisecond pulsars. We have analyzed a 38-pulsar subset of the 45 millisecond pulsars in the NANOGrav 11-year data set, in order to review their pulse profile stability. The most variability, on any timescale, is seen in PSRs J1713+0747, B1937+21, and J2145−0750. The strongest evidence for long-timescale pulse profile changes is seen in PSRs B1937+21 and J1643−1224. We have focused our analyses on these four pulsars in an attempt to elucidate the causes of their profile variability. Effects of scintillation seem to be responsible for the profile modifications of PSR J2145−0750. We see evidence that imperfect polarization calibration contributes to the profile variability of PSRs J1713+0747 and B1937+21, along with radio frequency interference around 2 GHz, but find that propagation effects also have an influence. The changes seen in PSR J1643−1224 have been reported previously, yet elude explanation beyond their astrophysical nature. Regardless of cause, unmodeled pulse profile changes are detrimental to the accuracy of pulsar timing and must be incorporated into the timing models where possible. [ABSTRACT FROM AUTHOR]

The frequency dependence of radio pulse arrival times provides a probe of structures in the intervening media. Demorest et al. was the first to show a short-term (∼100–200 days) reduction in the electron content along the line of sight to PSR J1713+0747 in data from 2008 (approximately MJD 54750) based on an apparent dip in the dispersion measure of the pulsar. We report on a similar event in 2016 (approximately MJD 57510), with average residual pulse-arrival times ≈−3.0, −1.3, and −0.7 μs at 820, 1400, and 2300 MHz, respectively. Timing analyses indicate possible departures from the standard ν−2 dispersive-delay dependence. We discuss and rule out a wide variety of potential interpretations. We find the likeliest scenario to be lensing of the radio emission by some structure in the interstellar medium, which causes multiple frequency-dependent pulse arrival-time delays. [ABSTRACT FROM AUTHOR]

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 AGWB < 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) h2 < 3.4 × 10−10. [ABSTRACT FROM AUTHOR]

We present a timing and glitch analysis of the young X-ray pulsar PSR J0537−6910, located within the Large Magellanic Cloud, using 13 yr of data from the now-decommissioned Rossi X-ray Timing Explorer. Rotating with a spin period of 16 ms, PSR J0537−6910 is the fastest-spinning and most energetic young pulsar known. It also displays the highest glitch activity of any known pulsar. We have found 42 glitches over the data span, corresponding to a glitch rate of 3.2 yr−1, with an overall glitch activity rate of . The high glitch frequency has allowed us to study the glitch behavior in ways that are inaccessible in other pulsars. We observe a strong linear correlation between spin frequency glitch magnitude and wait time to the following glitch. We also find that the post-glitch spin-down recovery is well described by a single two-component model fit to all glitches for which we have adequate input data. This consists of an exponential amplitude , decay timescale s, and linear slope . The latter slope corresponds to a second frequency derivative , from which we find an implied braking index . We also present a maximum likelihood technique for searching for periods in event-time data, which we used to both confirm previously published values and determine rotation frequencies in later observations. We discuss the implied constraints on glitch models from the observed behavior of this system, which we argue cannot be fully explained in the context of existing theories. [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]

We report on the discovery and timing observations of 29 distant long-period pulsars found in the ongoing Arecibo L-band Feed Array pulsar survey. Following discovery with the Arecibo Telescope, confirmation and timing observations of these pulsars over several years at Jodrell Bank Observatory have yielded high-precision positions and measurements of rotation and radiation properties. We have used multi-frequency data to measure the interstellar scattering properties of some of these pulsars. Most of the pulsars have properties that mirror those of the previously known pulsar population, although four show some notable characteristics. PSRs J1907+0631 and J1925+1720 are young and are associated with supernova remnants or plerionic nebulae: J1907+0631 lies close to the center of SNR G40.5−0.5, while J1925+1720 is coincident with a high-energy Fermi γ-ray source. One pulsar, J1932+1500, is in a surprisingly eccentric, 199 day binary orbit with a companion having a minimum mass of 0.33 M⊙. Several of the sources exhibit timing noise, and two, PSRs J0611+1436 and J1907+0631, have both suffered large glitches, but with very different post-glitch rotation properties. In particular, the rotational period of PSR J0611+1436 will not recover to its pre-glitch value for about 12 years, a far greater recovery timescale than seen following any other large glitches. [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]

We compute upper limits on the nanohertz-frequency isotropic stochastic gravitational wave background (GWB) using the 9 year data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) collaboration. Well-tested Bayesian techniques are used to set upper limits on the dimensionless strain amplitude (at a frequency of 1 yr−1) for a GWB from supermassive black hole binaries of . We also parameterize the GWB spectrum with a broken power-law model by placing priors on the strain amplitude derived from simulations of Sesana and McWilliams et al. Using Bayesian model selection we find that the data favor a broken power law to a pure power law with odds ratios of 2.2 and 22 to one for the Sesana and McWilliams prior models, respectively. Using the broken power-law analysis we construct posterior distributions on environmental factors that drive the binary to the GW-driven regime including the stellar mass density for stellar-scattering, mass accretion rate for circumbinary disk interaction, and orbital eccentricity for eccentric binaries, marking the first time that the shape of the GWB spectrum has been used to make astrophysical inferences. Returning to a power-law model, we place stringent limits on the energy density of relic GWs, . Our limit on the cosmic string GWB, , translates to a conservative limit on the cosmic string tension with , a factor of four better than the joint Planck and high-l cosmic microwave background data from other experiments. [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 present timing solutions and spin properties of the young pulsar PSR B0540–69 from analysis of 15.8 years of data from the Rossi X-Ray Timing Explorer. We perform a partially phase-coherent timing analysis in order to mitigate the pronounced effects of timing noise in this pulsar. We also perform fully coherent timing over large subsets of the data set in order to arrive at a more precise solution. In addition to the previously reported first glitch undergone by this pulsar, we find a second glitch, which occurred at MJD 52927 ± 4, with fractional changes in spin frequency Δν/ν = (1.64 ± 0.05) × 10−9 and spin-down rate (taken from our fully coherent analysis). We measure a braking index that is consistent over the entire data span, with a mean value n = 2.129 ± 0.012, from our partially coherent timing analysis. We also investigated the emission behavior of this pulsar, and have found no evidence for significant flux changes, flares, burst-type activity, or pulse profile shape variations. While there is strong evidence for the much-touted similarity of PSR B0540–69 to the Crab pulsar, they nevertheless differ in several aspects, including glitch activity, where PSR B0540–69 can be said to resemble certain other very young pulsars. It seems clear that the specific processes governing the formation, evolution, and interiors of this population of recently born neutron stars can vary significantly, as reflected in their observed properties. [ABSTRACT FROM AUTHOR]

Among efforts to detect gravitational radiation, pulsar timing arrays are uniquely poised to detect “memory” signatures, permanent perturbations in spacetime from highly energetic astrophysical events such as mergers of supermassive black hole binaries. The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) observes dozens of the most stable millisecond pulsars using the Arecibo and Green Bank radio telescopes in an effort to study, among other things, gravitational wave memory. We herein present the results of a search for gravitational wave bursts with memory (BWMs) using the first five years of NANOGrav observations. We develop original methods for dramatically speeding up searches for BWM signals. In the directions of the sky where our sensitivity to BWMs is best, we would detect mergers of binaries with reduced masses of out to distances of 30 Mpc; such massive mergers in the Virgo cluster would be marginally detectable. We find no evidence for BWMs. However, with our non-detection, we set upper limits on the rate at which BWMs of various amplitudes could have occurred during the time spanned by our data—e.g., BWMs with amplitudes greater than 10−13 must encounter the Earth at a rate less than 1.5 yr−1. [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]

A correction to the article "Timing of Five Millisecond Pulsars Discovered in the Palfa Survey" that was published in the April 2015 issue is presented.