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
A pulsar’s pulse profile gets broadened at low frequencies due to dispersion along the line of sight or due to multipath propagation. The dynamic nature of the interstellar medium makes both of these effects time-dependent and introduces slowly varying time delays in the measured times-of-arrival similar to those introduced by passing gravitational waves. In this article, we present an improved method to correct for such delays by obtaining unbiased dispersion measure (DM) measurements by using low-frequency estimates of the scattering parameters. We evaluate this method by comparing the obtained DM estimates with those, where scatter-broadening is ignored using simulated data. A bias is seen in the estimated DMs for simulated data with pulse-broadening with a larger variability for a data set with a variable frequency scaling index, $\alpha$, as compared to that assuming a Kolmogorov turbulence. Application of the proposed method removes this bias robustly for data with band averaged signal-to-noise ratio larger than 100. We report the measurements of the scatter-broadening time and $\alpha$ from analysis of PSR J1643$-$1224, observed with upgraded Giant Metrewave Radio Telescope as part of the Indian Pulsar Timing Array experiment. These scattering parameters were found to vary with epoch and $\alpha$ was different from that expected for Kolmogorov turbulence. Finally, we present the DM time-series after application of this technique to PSR J1643$-$1224.
more » « less- PAR ID:
- 10553970
- Author(s) / Creator(s):
- ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more »
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 535
- Issue:
- 1
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 1184-1192
- Size(s):
- p. 1184-1192
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract μ 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. -
more » « less
ABSTRACT With unparalleled rotational stability, millisecond pulsars (MSPs) serve as ideal laboratories for numerous astrophysical studies, many of which require precise knowledge of the distance and/or velocity of the MSP. Here, we present the astrometric results for 18 MSPs of the ‘MSPSR$\pi$’ project focusing exclusively on astrometry of MSPs, which includes the re-analysis of three previously published sources. On top of a standardized data reduction protocol, more complex strategies (i.e. normal and inverse-referenced 1D interpolation) were employed where possible to further improve astrometric precision. We derived astrometric parameters using sterne, a new Bayesian astrometry inference package that allows the incorporation of prior information based on pulsar timing where applicable. We measured significant (${>}3\, \sigma$) parallax-based distances for 15 MSPs, including 0.81 ± 0.02 kpc for PSR J1518+4904 – the most significant model-independent distance ever measured for a double neutron star system. For each MSP with a well-constrained distance, we estimated its transverse space velocity and radial acceleration. Among the estimated radial accelerations, the updated ones of PSR J1012+5307 and PSR J1738+0333 impose new constraints on dipole gravitational radiation and the time derivative of Newton’s gravitational constant. Additionally, significant angular broadening was detected for PSR J1643−1224, which offers an independent check of the postulated association between the HII region Sh 2-27 and the main scattering screen of PSR J1643−1224. Finally, the upper limit of the death line of γ-ray-emitting pulsars is refined with the new radial acceleration of the hitherto least energetic γ-ray pulsar PSR J1730−2304.
-
more » « less
Abstract The millisecond pulsar J1713+0747 underwent a sudden and significant pulse shape change between 2021 April 16 and 17 (MJDs 59320 and 59321). Subsequently, the pulse shape gradually recovered over the course of several months. We report the results of continued multifrequency radio observations of the pulsar made using the Canadian Hydrogen Intensity Mapping Experiment and the 100 m Green Bank Telescope in a 3 yr period encompassing the shape change event, between 2020 February and 2023 February. As of 2023 February, the pulse shape had returned to a state similar to that seen before the event, but with measurable changes remaining. The amplitude of the shape change and the accompanying time-of-arrival residuals display a strong nonmonotonic dependence on radio frequency, demonstrating that the event is neither a glitch (the effects of which should be independent of radio frequency,
ν ) nor a change in dispersion measure alone (which would produce a delay proportional toν −2). However, it does bear some resemblance to the two previous “chromatic timing events” observed in J1713+0747, as well as to a similar event observed in PSR J1643−1224 in 2015. -
more » « less
Abstract Broadband radio waves emitted from pulsars are distorted and delayed as they propagate toward the Earth due to interactions with the free electrons that compose the interstellar medium (ISM), with lower radio frequencies being more impacted than higher frequencies. Multipath propagation in the ISM results in both later times of arrival for the lower frequencies and causes the observed pulse to arrive with a broadened tail described via the pulse broadening function. We employ the CLEAN deconvolution technique to recover the pulse broadening timescale and by proxy the intrinsic pulse shape. This work expands upon previous descriptions of CLEAN deconvolution used in pulse broadening analyses by parameterizing the efficacy on simulated data and developing a suite of tests to establish which of a set of figures of merit leads to an automatic and consistent determination of the scattering timescale and its uncertainty. We compare our algorithm to the cyclic spectroscopy method of estimating the scattering timescale, specifically to the simulations performed in Dolch et al. (2021). We test our improved algorithm on the highly scattered millisecond pulsar J1903+0327, showing the scattering timescale to change over years, consistent with estimates of the refractive timescale of the pulsar.
-
more » « less
Abstract We investigate the use of bright single pulses from the Crab pulsar to determine separately the dispersion measure (DM) for the Main Pulse and Interpulse components. We develop two approaches using cross-correlation functions (CCFs). The first method computes the CCF of the total intensity of each of the 64 frequency channels with a reference channel and converts the time lag of maximum correlation into a DM. The second method separately computes the CCF between every pair of channels for each individual bright pulse and extracts an average DM from the distribution of all channel-pair DMs. Both methods allow the determination of the DM with a relative uncertainty of better than 10−5and provide robust estimates for the uncertainty of the best-fit value. We find differences in DM between the Main Pulse, the Low Frequency Interpulse, and the High Frequency Interpulse using both methods in a frequency range from 4 to 6 GHz. Earlier observations of the High Frequency Interpulse carried out by Hankins et al. (2016) resulted in DMHFIP–DMMPof 0.010 ± 0.016 pc cm−3. Our results indicate a DMHFIP–DMMPof 0.0127 ± 0.0011 pc cm−3(with DMcompbeing the DM value of the respective emission component), confirming earlier results with an independent method. During our studies we also find a relation between the brightness of single pulses in the High Frequency Interpulse and their DM. We also discuss the application of the developed methods on the identification of substructures in the case of Fast Radio Bursts.
