ABSTRACT We analyze 24 binary radio pulsars in the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) nine-year data set. We make 14 significant measurements of the Shapiro ...delay, including new detections in four pulsar-binary systems (PSRs J0613−0200, J2017+0603, J2302+4442, and J2317+1439), and derive estimates of the binary-component masses and orbital inclination for these MSP-binary systems. We find a wide range of binary pulsar masses, with values as low as for PSR J1918−0642 and as high as for PSR J1614−2230 (both 68.3% credibility). We make an improved measurement of the Shapiro timing delay in the PSR J1918−0642 and J2043+1711 systems, measuring the pulsar mass in the latter system to be (68.3% credibility) for the first time. We measure secular variations of one or more orbital elements in many systems, and use these measurements to further constrain our estimates of the pulsar and companion masses whenever possible. In particular, we used the observed Shapiro delay and periastron advance due to relativistic gravity in the PSR J1903+0327 system to derive a pulsar mass of (68.3% credibility). We discuss the implications that our mass measurements have on the overall neutron-star mass distribution, and on the "mass/orbital-period" correlation due to extended mass transfer.
ABSTRACT Using the nine-year radio-pulsar timing data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), collected at Arecibo Observatory and the Green Bank ...Telescope, we have measured the positions, proper motions, and parallaxes for 37 millisecond pulsars. We report twelve significant parallax measurements and distance measurements, and eighteen lower limits on distance. We compare these measurements to distances predicted by the NE2001 interstellar electron density model and find them to be in general agreement. We use measured orbital-decay rates and spin-down rates to confirm two of the parallax distances and to place distance upper limits on other sources; these distance limits agree with the parallax distances with one exception, PSR J1024-0719, which we discuss at length. Using the proper motions of the 37 NANOGrav pulsars in combination with other published measurements, we calculate the velocity dispersion of the millisecond pulsar population in Galactocentric coordinates. We find the radial, azimuthal, and perpendicular dispersions to be 46, 40, and 24 , respectively, in a model that allows for high-velocity outliers; or 81, 58, and 62 for the full population. These velocity dispersions are far smaller than those of the canonical pulsar population, and are similar to older Galactic disk populations. This suggests that millisecond pulsar velocities are largely attributable to their being an old population rather than being artifacts of their birth and evolution as neutron star binary systems. The components of these velocity dispersions follow similar proportions to other Galactic populations, suggesting that our results are not biased by selection effects.
We have made observations of 98 low Galactic latitude pulsars to measure pulse broadening caused by multipath propagation through the interstellar medium. Data were collected with the 305 m Arecibo ...telescope at four radio frequencies between 430 and 2380 MHz. We used a CLEAN-based algorithm to deconvolve interstellar pulse broadening from the measured pulse shapes. We employed two distinct pulse-broadening functions (PBFs): PBF sub(1) is appropriate for a thin screen of scattering material between the Earth and a pulsar, while PBF sub(2) is appropriate for scattering material uniformly distributed along the line of sight from the Earth to a pulsar. We found that some observations were better fitted by PBF sub(1) and some by PBF sub(2). Pulse-broadening times ( tau sub(d)) are derived from fits of PBFs to the data and are compared with the predictions of a smoothed model of the Galactic electron distribution. Several lines of sight show excess broadening, which we model as clumps of high- density scattering material. A global analysis of all available data finds that the pulse broadening scales with frequency, nu , as tau sub(d) proportional nu super(-alpha), where alpha approx 3.9 plus or minus 0.2. This is somewhat shallower than the value alpha = 4.4 expected from a Kolmogorov medium but could arise if the spectrum of turbulence has an inner cutoff at approx300-800 km. A few objects follow particularly shallow scaling laws (the mean scaling index ang10alphaang0 approx 3.1 plus or minus 0.1 and approx3.8 plus or minus 0.2, respectively, for PBF sub(1) and PBF sub(2)), which may arise from large-scale refraction or from the truncation of scattering screens transverse to the Earth-pulsar line of sight.
ABSTRACT PSR J1024-0719 is a millisecond pulsar that was long thought to be isolated. However, puzzling results concerning its velocity, distance, and low rotational period derivative have led to a ...reexamination of its properties. We present updated radio timing observations along with new and archival optical data which show that PSR J1024-0719 is most likely in a long-period (2-20 kyr) binary system with a low-mass ( ), low-metallicity ( dex) main-sequence star. Such a system can explain most of the anomalous properties of this pulsar. We suggest that this system formed through a dynamical exchange in a globular cluster that ejected it into a halo orbit, which is consistent with the low observed metallicity for the stellar companion. Further astrometric and radio timing observations such as measurement of the third period derivative could strongly constrain the range of orbital parameters.
ABSTRACT We report on an effort to extract and monitor interstellar scintillation parameters in regular timing observations collected for the North American Nanohertz Observatory for Gravitational ...Waves pulsar timing array. Scattering delays are measured by creating dynamic spectra for each pulsar and observing epoch of wide-band observations centered near 1500 MHz and carried out at the Green Bank Telescope and the Arecibo Observatory. The ∼800 MHz wide frequency bands imply dramatic changes in scintillation bandwidth across the bandpass, and a stretching routine has been included to account for this scaling. For most of the 10 pulsars for which the scaling has been measured, the bandwidths scale with frequency less steeply than expected for a Kolmogorov medium. We find estimated scattering delay values that vary with time by up to an order of magnitude. The mean measured scattering delays are similar to previously published values and are slightly higher than predicted by interstellar medium models. We investigate the possibility of increasing the timing precision by mitigating timing errors introduced by the scattering delays. For most of the pulsars, the uncertainty in the time of arrival of a single timing point is much larger than the maximum variation of the scattering delay, suggesting that diffractive scintillation remains as only a negligible part of their noise budget.
Abstract
Using Bayesian analyses we study the solar electron density with the NANOGrav 11 yr pulsar timing array (PTA) data set. Our model of the solar wind is incorporated into a global fit starting ...from pulse times of arrival. We introduce new tools developed for this global fit, including analytic expressions for solar electron column densities and open source models for the solar wind that port into existing PTA software. We perform an ab initio recovery of various solar wind model parameters. We then demonstrate the richness of information about the solar electron density,
n
E
, that can be gleaned from PTA data, including higher order corrections to the simple 1/
r
2
model associated with a free-streaming wind (which are informative probes of coronal acceleration physics), quarterly binned measurements of
n
E
and a continuous time-varying model for
n
E
spanning approximately one solar cycle period. Finally, we discuss the importance of our model for chromatic noise mitigation in gravitational-wave analyses of pulsar timing data and the potential of developing synergies between sophisticated PTA solar electron density models and those developed by the solar physics community.
Abstract Noise characterization for pulsar-timing applications accounts for interstellar dispersion by assuming a known frequency dependence of the delay it introduces in the times of arrival (TOAs). ...However, calculations of this delay suffer from misestimations due to other chromatic effects in the observations. The precision in modeling dispersion is dependent on the observed bandwidth. In this work, we calculate the offsets in infinite-frequency TOAs due to misestimations in the modeling of dispersion when using varying bandwidths at the Green Bank Telescope. We use a set of broadband observations of PSR J1643−1224, a pulsar with unusual chromatic timing behavior. We artificially restricted these observations to a narrowband frequency range, then used both the broad- and narrowband data sets to calculate residuals with a timing model that does not account for time variations in the dispersion. By fitting the resulting residuals to a dispersion model and comparing the fits, we quantify the error introduced in the timing parameters due to using a reduced frequency range. Moreover, by calculating the autocovariance function of the parameters, we obtained a characteristic timescale over which the dispersion misestimates are correlated. For PSR J1643−1224, which has one of the highest dispersion measures (DM) in the NANOGrav pulsar timing array, we find that the infinite-frequency TOAs suffer from a systematic offset of ∼22 μ s due to incomplete frequency sampling, with correlations over about one month. For lower-DM pulsars, the offset is ∼7 μ s. This error quantification can be used to provide more robust noise modeling in the NANOGrav data, thereby increasing the sensitivity and improving the parameter estimation in gravitational wave searches.
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 at1500 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 at820 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 are 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 nanoseconds or less.
PSR J0751+1807 is a millisecond pulsar in a circular 6 hr binary system with a helium white dwarf secondary. Through high-precision pulse timing measurements with the Arecibo and Effelsberg radio ...telescopes, we have detected the decay of its orbit due to emission of gravitational radiation. This is the first detection of the relativistic orbital decay of a low-mass, circular binary pulsar system. The measured rate of change in orbital period, corrected for acceleration biases, is P( )= (-6.4 c 0.9) x 10(-14). Interpreted in the context of general relativity, and combined with measurement of Shapiro delay, it implies a pulsar mass of 2.1 c 0.2 M(z), the most massive pulsar measured. This adds to the emerging trend toward relatively high neutron star masses in neutron star-white dwarf binaries. In addition, there is some evidence for an inverse correlation between pulsar mass and orbital period in these systems. We consider alternatives to the general relativistic analysis of the data, and we use the pulsar timing data to place limits on violations of the strong equivalence principle.
We present the polarization pulse profiles for 28 pulsars observed with the Arecibo Observatory by the North American Nanohertz Observatory for Gravitational Waves timing project at 2.1 GHz, 1.4 GHz, ...and 430 MHz. These profiles represent some of the most sensitive polarimetric millisecond pulsar profiles to date, revealing the existence of microcomponents (that is, pulse components with peak intensities much lower than the total pulse peak intensity). Although microcomponents have been detected in some pulsars previously, we present microcomponents for PSR B1937+21, PSR J1713+0747, and PSR J2234+0944 for the first time. These microcomponents can have an impact on pulsar timing, geometry, and flux density determination. We present rotation measures for all 28 pulsars, determined independently at different observation frequencies and epochs, and find the Galactic magnetic fields derived from these rotation measures to be consistent with current models. These polarization profiles were made using measurement equation template matching, which allows us to generate the polarimetric response of the Arecibo Observatory on an epoch-by-epoch basis. We use this method to describe its time variability and find that the polarimetric responses of the Arecibo Observatory's 1.4 and 2.1 GHz receivers vary significantly with time.