Pulsar Scintillation through Thick and Thin: Bow Shocks, Bubbles, and the Broader Interstellar Medium

Observations of pulsar scintillation are among the few astrophysical probes of very small-scale ($\lesssim$ au) phenomena in the interstellar medium (ISM). In particular, characterization of scintillation arcs, including their curvature and intensity distributions, can be related to interstellar turbulence and potentially over-pressurized plasma in local ISM inhomogeneities, such as supernova remnants, HII regions, and bow shocks. Here we present a survey of eight pulsars conducted at the Five-hundred-meter Aperture Spherical Telescope (FAST), revealing a diverse range of scintillation arc characteristics at high sensitivity. These observations reveal more arcs than measured previously for our sample. At least nine arcs are observed toward B1929$+$10 at screen distances spanning $\sim 90\%$ of the pulsar's $361$ pc path-length to the observer. Four arcs are observed toward B0355$+$54, with one arc yielding a screen distance as close as $\sim10^5$ au ($<1$ pc) from either the pulsar or the observer. Several pulsars show highly truncated, low-curvature arcs that may be attributable to scattering near the pulsar. The scattering screen constraints are synthesized with continuum maps of the local ISM and other well-characterized pulsar scintillation arcs, yielding a three-dimensional view of the scattering media in context.Comment: 20 pages, 14 figures. Submitted to MNRAS and comments welcome. Interactive version of Figure 12 available at https://stella-ocker.github.io/scattering_ism3d_ocker202

Testing Theories of Gravitation Using 21-Year Timing of Pulsar Binary J1713+0747

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 similar to 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-circle dot) and the companion white dwarf (0.286 +/- 0.012 M-circle dot) as well as their parallax distance 1.15 +/- 0.03 kpc. We measured the intrinsic change in orbital period, (P) over dot(b)(Int), is -0.20 +/- 0.17 ps s(-1), which is not distinguishable from zero. This result, combined with the measured (P) over dot(b)(Int) of other pulsars, can place a generic limit on potential changes in the gravitational constant G. We found that (G) over dot/G is consistent with zero [(-0.6 +/- 1.1) x 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 (G) over dot/G limit from pulsar binary systems. The (P) over dot(b)(Int) of pulsar binaries can also place limits on the putative coupling constant for dipole gravitational radiation kappa(D) = (-0.9 +/- 3.3) 10(-4) (95% confidence). Finally, the nearly circular orbit of this pulsar binary allows us to constrain statistically the strong-field post-Newtonian parameters Delta, which describes the violation of strong equivalence principle, and (alpha) over cap (3), which describes a breaking of both Lorentz invariance in gravitation and conservation of momentum. We found, at 95% confidence, Delta <0.01 and (3) <2 x 10(-20) based on PSR J1713+0747

The NANOGrav Nine-year Data Set:Noise Budget for Pulsar Arrival Times on Intraday Timescales

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 $\sim 1\%$ 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

The NANOGrav Nine-year Data Set: Measurement and Analysis of Variations in Dispersion Measures

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

The NANOGrav 11-Year Data Set: Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries

Observations indicate that nearly all galaxies contain supermassive black holes (SMBHs) 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 (PTAs). We have searched the recently-released North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year 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 as a function of GW frequency and sky location. We placed a sky-averaged upper limit on the GW strain of $h_0 < 7.3(3) \times 10^{-15}$ at $f_\mathrm{gw}= 8$ nHz. 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 in measurements from individual pulsars. We used our upper limits on the GW strain to place lower limits on the distances to individual SMBHBs. At the most-sensitive sky location, we ruled out SMBHBs emitting GWs with $f_\mathrm{gw}= 8$ nHz within 120 Mpc for $\mathcal{M} = 10^9 \, M_\odot$, and within 5.5 Gpc for $\mathcal{M} = 10^{10} \, M_\odot$. We also determined that there are no SMBHBs with $\mathcal{M} > 1.6 \times 10^9 \, M_\odot$ emitting GWs in the Virgo Cluster. Finally, we estimated the number of potentially detectable sources given our current strain upper limits based on galaxies in Two Micron All-Sky Survey (2MASS) and merger rates from the Illustris cosmological simulation project. Only 34 out of 75,000 realizations of the local Universe contained a detectable source, from which we concluded it was unsurprising that we did not detect any individual sources given our current sensitivity to GWs.Comment: 10 pages, 11 figures. Accepted by Astrophysical Journal. Please send any comments/questions to S. J. Vigeland ([email protected]

A 24 HR global campaign to assess precision timing of the millisecond pulsar J1713+0747

The radiomillisecond pulsar J1713+0747 is regarded as one of the highest-precision clocks in the sky and is regularly timed for the purpose of detecting gravitational waves. The International Pulsar Timing Array Collaboration undertook a 24 hr global observation of PSR J1713+0747 in an effort to better quantify sources of timing noise in this pulsar, particularly on intermediate (1.24 hr) timescales. We observed the pulsar continuously over 24 hr with the Arecibo, Effelsberg, GMRT, Green Bank, LOFAR, Lovell, Nancay, Parkes, and WSRT radio telescopes. The combined pulse times-of-arrival presented here provide an estimate of what sources of timing noise, excluding DM variations, would be present as compared to an idealized N√ improvement in timing precision, where N is the number of pulses analyzed. In the case of this particular pulsar, we find that intrinsic pulse phase jitter dominates arrival time precision when the signal-to-noise ratio of single pulses exceeds unity, as measured using the eight telescopes that observed at L band/1.4 GHz.We present first results of specific phenomena probed on the unusually long timescale (for a single continuous observing session) of tens of hours, in particular interstellar scintillation, and discuss the degree to which scintillation and profile evolution affect precision timing. This paper presents the data set as a basis for future, deeper studies

The Green Bank North Celestial Cap Survey. IX. Timing Follow-up for 128 Pulsars

The Green Bank North Celestial Cap survey is one of the largest and most sensitive searches for pulsars and transient radio objects. Observations for the survey have finished; priorities have shifted toward long-term monitoring of its discoveries. In this study, we have developed a pipeline to handle large data sets of archival observations and connect them to recent, high-cadence observations taken using the Canadian Hydrogen Intensity Mapping Experiment telescope. This pipeline handles data for 128 pulsars and has produced measurements of spin, positional, and orbital parameters that connect data over observation gaps as large as 2000 days. We have also measured glitches in the timing residuals for five of the pulsars included and proper motion for 19 sources (13 new). We include updates to orbital parameters for 19 pulsars, including nine previously unpublished binaries. For two of these binaries, we provide updated measurements of post-Keplerian binary parameters, which result in much more precise estimates of the total masses of both systems. For PSR J0509+3801, the much improved measurement of the Einstein delay yields much improved mass measurements for the pulsar and its companion, 1.399(6) M⊙ and 1.412(6) M⊙, respectively. For this system, we have also obtained a measurement of the orbital decay due to the emission of gravitational waves

The NANOGrav 11 yr data set: Constraints on planetary masses around 45 millisecond pulsars

We search for extrasolar planets around millisecond pulsars using pulsar timing data and seek to determine the minimum detectable planetary masses as a function of orbital period. Using the 11 yr data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), we look for variations from our models of pulse arrival times due to the presence of exoplanets. No planets are detected around the millisecond pulsars in the NANOGrav 11 yr data set, but taking into consideration the noise levels of each pulsar and the sampling rate of our observations, we develop limits that show we are sensitive to planetary masses as low as that of the moon. We analyzed potential planet periods, P, in the range 7 days < P < 2000 days, with somewhat smaller ranges for some binary pulsars. The planetary-mass limit for our median-sensitivity pulsar within this period range is