Search

The regularity of pulsar emissions becomes apparent once we reference the pulses’ times of arrivals to the inertial rest frame of the solar system. It follows that errors in the determination of Earthʼs position with respect to the solar system barycenter can appear as a time-correlated bias in pulsar-timing residual time series, affecting the searches for low-frequency gravitational waves performed with pulsar-timing arrays. Indeed, recent array data sets yield different gravitational-wave background upper limits and detection statistics when analyzed with different solar system ephemerides. Crucially, the ephemerides do not generally provide usable error representations. In this article, we describe the motivation, construction, and application of a physical model of solar system ephemeris uncertainties, which focuses on the degrees of freedom (Jupiterʼs orbital elements) most relevant to gravitational-wave searches with pulsar-timing arrays. This model, BAYESEPHEM, was used to derive ephemeris-robust results in NANOGravʼs 11 yr stochastic-background search, and it provides a foundation for future searches by NANOGrav and other consortia. The analysis and simulations reported here suggest that ephemeris modeling reduces the gravitational-wave sensitivity of the 11 yr data set and that this degeneracy will vanish with improved ephemerides and with pulsar-timing data sets that extend well beyond a single Jovian orbital period.

An ensemble of inspiraling supermassive black hole binaries should produce a stochastic background of very low frequency gravitational waves. This stochastic background is predicted to be a power law, with a gravitational-wave strain spectral index of −2/3, and it should be detectable by a network of precisely timed millisecond pulsars, widely distributed on the sky. This paper reports a new “time slicing” analysis of the 11 yr data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) using 34 millisecond pulsars. Methods to flag potential “false-positive” signatures are developed, including techniques to identify responsible pulsars. Mitigation strategies are then presented. We demonstrate how an incorrect noise model can lead to spurious signals, and we show how independently modeling noise across 30 Fourier components, spanning NANOGrav’s frequency range, effectively diagnoses and absorbs the excess power in gravitational-wave searches. This results in a nominal, and expected, progression of our gravitational-wave statistics. Additionally, we show that the first interstellar medium event in PSR J1713+0747 pollutes the common red-noise process with low spectral index noise, and we use a tailored noise model to remove these effects.

The mergers of supermassive black hole binaries (SMBHBs) promise to be incredible sources of gravitational waves (GWs).While the oscillatory part of the merger gravitational waveform will be outside the frequency sensitivity range of pulsar timing arrays, the nonoscillatory GW memory effect is detectable. Further, any burst of GWs will produce GW memory, making memory a useful probe of unmodeled exotic sources and new physics. We searched the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11 yr data set for GW memory. This data set is sensitive to very low-frequency GWs of ∼3 to 400 nHz (periods of ∼11 yr–1 month). Finding no evidence for GWs, we placed limits on the strain amplitude of GW memory events during the observation period. We then used the strain upper limits to place limits on the rate of GW memory causing events. At a strain of 2.5 × 10−14, corresponding to the median upper limit as a function of source sky position, we set a limit on the rate of GW memory events at <0.4 yr−1. That strain corresponds to an SMBHB merger with reduced mass of ηM ~2 × 1010Mmoonand inclination of ι = π/3 at a distance of 1 Gpc. As a test of our analysis, we analyzed the NANOGrav 9 yr data set as well. This analysis found an anomolous signal, which does not appear in the 11 yr data set. This signal is not a GW, and its origin remains unknown.

The Neutron Star Interior Composition Explorer (NICER) is an X-ray astrophysics payload on the International Space Station. It enables unprecedented high-precision timing of millisecond pulsars (MSPs) without the pulse broadening and delays due to dispersion and scattering within the interstellar medium that plague radio timing. We present initial timing results from a year of data on the MSPs PSR B1937+21 and PSR J0218+4232, and nine months of data on PSR B1821−24. NICER time-of-arrival uncertainties for the three pulsars are consistent with theoretical lower bounds and simulations based on their pulse shape templates and average source and background photon count rates. To estimate timing stability, we use the σz measure, which is based on the average of the cubic coefficients of polynomial fits to subsets of timing residuals. So far we are achieving timing stabilities σz≈3×10−14 for PSR B1937+21 and on the order of 10−12 for PSRs B1821−24 and J0218+4232. Within the span of our NICER data we do not yet see the characteristic break point in the slope of σz; detection of such a break would indicate that further improvement in the cumulative root-mean-square timing residual is limited by timing noise. We see this break point in our comparison radio data sets for PSR B1821−24 and PSR B1937+21 on timescales of >2 yr.

We report the timing results for PSR J2234+0611, a 3.6 ms pulsar in a 32 day, eccentric (e = 0.13) orbit with a helium white dwarf. The precise timing and eccentric nature of the orbit allow measurements of an unusual number of parameters: (a) a precise proper motion of 27.10(3) mas yr−1 and a parallax of 1.05(4) mas resulting in a pulsar distance of 0.95(4) kpc; enabling an estimate of the transverse velocity, 123(5) km s−1. Together with previously published spectroscopic measurements of the systemic radial velocity, this allows a 3D determination of the system’s velocity; (b) precise measurements of the rate of advance of periastron yields a total system mass of 1.6518 +0.0033-0.0035 Msun; (c) a Shapiro delay measurement, h3 = 82 ± 14 ns, despite the orbital inclination not being near 90°; combined with the measurement of the total mass yields a pulsar mass of 1.353 +0.014-0.017 Msun and a companion mass of 0.298 +0.015-0.0120 Msun; (d) we measure precisely the secular variation of the projected semimajor axis and detect a significant annual orbital parallax; together these allow a determination of the 3D orbital geometry of the system, including an unambiguous orbital inclination (i = 138.7 +2.5-2.2 deg) and a position angle for the line of nodes ( omega = 44 +5-4 deg). We discuss the component masses to investigate the hypotheses previously advanced to explain the origin of eccentric MSPs. The unprecedented determination of the 3D position, motion, and orbital orientation of the system, plus the precise pulsar and WD masses and the latter’s optical detection make this system a unique test of our understanding of white dwarfs and their atmospheres.