Pulsar timing experiments require high-fidelity template profiles in order to minimize the biases in pulse time-of-arrival (TOA) measurements and their uncertainties. Efforts to acquire more precise ...TOAs given the fixed effective area of telescopes, finite receiver noise, and limited integration time have led pulsar astronomers to the solution of implementing ultra-wideband receivers. This solution, however, has run up against the problem that pulse profile shapes evolve with frequency, which raises the question of how to properly measure and analyze TOAs obtained using template-matching methods. This paper proposes a new method for one facet of this problem, that of template profile generation, and demonstrates it on the well-timed millisecond pulsar J1713+0747. Specifically, we decompose pulse profile evolution into a linear combination of basis eigenvectors, the coefficients of which change slowly with frequency such that their evolution is modeled simply by a sum of low-degree piecewise polynomial spline functions. These noise-free, high-fidelity, frequency-dependent templates can be used to make measurements of so-called "wideband TOAs" simultaneously with an estimate of the instantaneous dispersion measure. The use of wideband TOAs is becoming important for pulsar timing array experiments, as the volume of data sets comprised of conventional, subbanded TOAs are quickly becoming unwieldy for the Bayesian analyses needed to uncover latent gravitational wave signals. Although motivated by high-precision timing experiments, our technique is applicable in more general pulsar observations.
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
We report on Bayesian estimation of the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, conditional on pulse-profile modeling of Neutron Star Interior ...Composition Explorer X-ray Timing Instrument event data. We condition on informative pulsar mass, distance, and orbital inclination priors derived from the joint North American Nanohertz Observatory for Gravitational Waves and Canadian Hydrogen Intensity Mapping Experiment/Pulsar wideband radio timing measurements of Fonseca et al. We use XMM-Newton European Photon Imaging Camera spectroscopic event data to inform our X-ray likelihood function. The prior support of the pulsar radius is truncated at 16 km to ensure coverage of current dense matter models. We assume conservative priors on instrument calibration uncertainty. We constrain the equatorial radius and mass of PSR J0740+6620 to be
12.39
−
0.98
+
1.30
km and
2.072
−
0.066
+
0.067
M
⊙
respectively, each reported as the posterior credible interval bounded by the 16% and 84% quantiles, conditional on surface hot regions that are non-overlapping spherical caps of fully ionized hydrogen atmosphere with uniform effective temperature; a posteriori, the temperature is
log
10
(
T
K
)
=
5.99
−
0.06
+
0.05
for each hot region. All software for the X-ray modeling framework is open-source and all data, model, and sample information is publicly available, including analysis notebooks and model modules in the Python language. Our marginal likelihood function of mass and equatorial radius is proportional to the marginal joint posterior density of those parameters (within the prior support) and can thus be computed from the posterior samples.
ELEMENTARY WIDEBAND TIMING OF RADIO PULSARS Pennucci, Timothy T; Demorest, Paul B; Ransom, Scott M
The Astrophysical journal,
08/2014, Volume:
790, Issue:
2
Journal Article
Peer reviewed
We present an algorithm for the simultaneous measurement of a pulse time-of-arrival (TOA) and dispersion measure (DM) from folded wideband pulsar data. We extend the prescription from Taylor's 1992 ...work to accommodate a general two-dimensional template "portrait," the alignment of which can be used to measure a pulse phase and DM. We show that there is a dedispersion reference frequency that removes the covariance between these two quantities and note that the recovered pulse profile scaling amplitudes can provide useful information. We experiment with pulse modeling by using a Gaussian-component scheme that allows for independent component evolution with frequency, a "fiducial component," and the inclusion of scattering. We showcase the algorithm using our publicly available code on three years of wideband data from the bright millisecond pulsar J1824-2452A (M28A) from the Green Bank Telescope, and a suite of Monte Carlo analyses validates the algorithm. By using a simple model portrait of M28A, we obtain DM trends comparable to those measured by more standard methods, with improved TOA and DM precisions by factors of a few. Measurements from our algorithm will yield precisions at least as good as those from traditional techniques, but is prone to fewer systematic effects and is without ad hoc parameters. A broad application of this new method for dispersion measure tracking with modern large-bandwidth observing systems should improve the timing residuals for pulsar timing array experiments, such as the North American Nanohertz Observatory for Gravitational Waves.
We present high-precision timing data over time spans of up to 11 years for 45 millisecond pulsars observed as part of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) ...project, aimed at detecting and characterizing low-frequency gravitational waves. The pulsars were observed with the Arecibo Observatory and/or the Green Bank Telescope at frequencies ranging from 327 MHz to 2.3 GHz. Most pulsars were observed with approximately monthly cadence, and six high-timing-precision pulsars were observed weekly. All were observed at widely separated frequencies at each observing epoch in order to fit for time-variable dispersion delays. We describe our methods for data processing, time-of-arrival (TOA) calculation, and the implementation of a new, automated method for removing outlier TOAs. We fit a timing model for each pulsar that includes spin, astrometric, and (for binary pulsars) orbital parameters; time-variable dispersion delays; and parameters that quantify pulse-profile evolution with frequency. The timing solutions provide three new parallax measurements, two new Shapiro delay measurements, and two new measurements of significant orbital-period variations. We fit models that characterize sources of noise for each pulsar. We find that 11 pulsars show significant red noise, with generally smaller spectral indices than typically measured for non-recycled pulsars, possibly suggesting a different origin. A companion paper uses these data to constrain the strength of the gravitational-wave background.
Many physically motivated extensions to general relativity (GR) predict substantial deviations in the properties of spacetime surrounding massive neutron stars. We report the measurement of a 2.01 ± ...0.04 solar mass (M⊙) pulsar in a 2.46-hour orbit with a 0.172 ± 0.003 M⊙ white dwarf. The high pulsar mass and the compact orbit make this system a sensitive laboratory of a previously untested strong-field gravity regime. Thus far, the observed orbital decay agrees with GR, supporting its validity even for the extreme conditions present in the system. The resulting constraints on deviations support the use of GR-based templates for ground-based gravitational wave detectors. Additionally, the system strengthens recent constraints on the properties of dense matter and provides insight to binary stellar astrophysics and pulsar recycling.
The NANOGrav Collaboration reported strong Bayesian evidence for a common-spectrum stochastic process in its12.5 yr pulsar timing array data set, with median characteristic strain amplitude at ...periods of a year of A(yr) = 1.92(+0.75,-0.55) x 10^(-15). However, evidence for the quadrupolar Hellings & Downs interpulsar correlations, which are characteristic of gravitational-wave signals, was not yet significant. We emulate and extend the NANOGrav data set, injecting a wide range of stochastic gravitational-wave background(GWB)signals that encompass a variety of amplitudes and spectral shapes, and quantify three key milestones.(I)Given the amplitude measured in the 12.5 yr analysis and assuming this signal is a GWB, we expect to accumulate robust evidence of an interpulsar-correlated GWB signal with 15–17 yr of data, i.e., an additional 2–5 yr from the 12.5 yr data set.(II)At the initial detection, we expect a fractional uncertainty of 40% on the power-law strain spectrum slope, which is sufficient to distinguish a GWB of supermassive black hole binary origin from some models predicting more exotic origins.(III)Similarly, the measured GWB amplitude will have an uncertainty of 44% upon initial detection, allowing us to arbitrate between some population models of supermassive black hole binaries. In addition, power-law models are distinguishable from those having low-frequency spectral turnovers once 20 yr of data are reached. Even though our study is based on the NANOGrav data, we also derive relations that allow for a generalization to other pulsar timing array data sets. Most notably, by combining the data of individual arrays into the International Pulsar Timing Array, all of these milestones can be reached significantly earlier.
Abstract
In 2018 an ultra–wide-bandwidth low-frequency (UWL) receiver was installed on the 64 m Parkes Radio Telescope, enabling observations with an instantaneous frequency coverage from 704 to 4032 ...MHz. Here we present the analysis of a 3 yr data set of 35 ms pulsars observed with the UWL by the Parkes Pulsar Timing Array, using wide-band timing methods. The two key differences compared to typical narrowband methods are (1) generation of two-dimensional templates accounting for pulse shape evolution with frequency and (2) simultaneous measurements of the pulse time of arrival (TOA) and dispersion measure (DM). This is the first time that wide-band timing has been applied to a uniform data set collected with a single large fractional bandwidth receiver, for which such techniques were originally developed. As a result of our study, we present a set of profile evolution models and new timing solutions, including initial noise analysis. Precision of our TOA and DM measurements is in the range of 0.005–2.08
μ
s and (0.043–14.24) × 10
−4
cm
−3
pc, respectively, with 94% of the pulsars achieving a median TOA uncertainty of less than 1
μ
s.
Abstract
Modeling of frequency-dependent effects, contributed by the turbulence in the free electron density of interstellar plasma, is required to enable the detection of the expected imprints from ...the stochastic gravitational-wave (GW) background in pulsar timing data. In this work, we present an investigation of temporal variations of interstellar medium for a set of millisecond pulsars (MSPs) with the upgraded Giant Metrewave Radio Telescope (GMRT) aided by large fractional bandwidth at lower observing frequencies. Contrary to the conventional narrowband analysis using a frequency-invariant template profile, we applied
PulsePortraiture
-based wide-band timing analysis while correcting for the evolution of the pulsar profile with frequency. Implementation of the
PulsePortraiture
-based wide-band timing method for the GMRT-discovered MSPs to probe the dispersion measure (DM) variations resulted in a DM precision of 10
−4
pc cm
−3
. In general, we achieve similar DM and timing precision from wide-band timing compared to the narrowband timing with matching temporal variations of DMs. This wide-band timing study of newly discovered MSPs over a wide frequency range highlights the effectiveness of profile modeling at low frequencies and probes the potential of using them in a pulsar timing array.
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.