skip to main content

Title: Scattering Delay Mitigation in High-accuracy Pulsar Timing: Cyclic Spectroscopy Techniques

We simulate scattering delays from the interstellar medium to examine the effectiveness of three estimators in recovering these delays in pulsar timing data. Two of these estimators use the more traditional process of fitting autocorrelation functions to pulsar dynamic spectra to extract scintillation bandwidths, while the third estimator uses the newer technique of cyclic spectroscopy on baseband pulsar data to recover the interstellar medium’s impulse response function. We find that either fitting a Lorentzian or Gaussian distribution to an autocorrelation function or recovering the impulse response function from the cyclic spectrum are, on average, accurate in recovering scattering delays, although autocorrelation function estimators have a large variance, even at high signal-to-noise ratio (S/N). We find that, given sufficient S/N, cyclic spectroscopy is more accurate than both Gaussian and Lorentzian fitting for recovering scattering delays at specific epochs, suggesting that cyclic spectroscopy is a superior method for scattering estimation in high-quality data.

more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Medium: X Size: Article No. 191
Article No. 191
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Radio pulsar signals are significantly perturbed by their propagation through the ionized interstellar medium. In addition to the frequency-dependent pulse times of arrival due to dispersion, pulse shapes are also distorted and shifted, having been scattered by the inhomogeneous interstellar plasma, affecting pulse arrival times. Understanding the degree to which scattering affects pulsar timing is important for gravitational-wave detection with pulsar timing arrays (PTAs), which depend on the reliability of pulsars as stable clocks with an uncertainty of ∼100 ns or less over ∼10 yr or more. Scattering can be described as a convolution of the intrinsic pulse shape with an impulse response function representing the effects of multipath propagation. In previous studies, the technique of cyclic spectroscopy has been applied to pulsar signals to deconvolve the effects of scattering from the original emitted signals, increasing the overall timing precision. We present an analysis of simulated data to test the quality of deconvolution using cyclic spectroscopy over a range of parameters characterizing interstellar scattering and pulsar signal-to-noise ratio (S/N). We show that cyclic spectroscopy is most effective for high S/N and/or highly scattered pulsars. We conclude that cyclic spectroscopy could play an important role in scattering correction to distant populations of highly scattered pulsars not currently included in PTAs. For future telescopes and for current instruments such as the Green Bank Telescope upgraded with the ultrawide bandwidth receiver, cyclic spectroscopy could potentially double the number of PTA-quality pulsars. 
    more » « less
  2. 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

    The parabolic structure of the secondary or conjugate spectra of pulsars is often the result of isolated one-dimensional (or at least highly anisotropic) lenses in the interstellar medium (ISM). The curvature of these features contains information about the velocities of the Earth, ISM, and pulsar along the primary axis of the lens. As a result, measuring variations in the curvature over the course of a year, or the orbital period for pulsars in binaries, can constrain properties of the screen and pulsar. In particular, the pulsar distance and orbital inclination for binary systems can be found for multiple screens or systems with prior information on sin(i). By mapping the conjugate spectra into a space where the main arc and inverted arclets are straight lines, we are able to make use of the full information content from the inverted arclet curvatures, amplitudes, and phases using eigenvectors to uniquely and optimally retrieve phase information. This allows for a higher precision measurement than the standard Hough transform for systems where these features are available. Our technique also directly yields the best-fitting one-dimensional impulse response function for the interstellar lens given in terms of the Doppler shift, time delay, and magnification of images on the sky as seen from a single observatory. This can be extended for use in holographic imaging of the lens by combining multiple telescopes. We present examples of this new method for both simulated data and actual observations of PSR B0834+06.

    more » « less
  4. Mariño, Inés P. (Ed.)
    In many physiological systems, real-time endogeneous and exogenous signals in living organisms provide critical information and interpretations of physiological functions; however, these signals or variables of interest are not directly accessible and must be estimated from noisy, measured signals. In this paper, we study an inverse problem of recovering gas exchange signals of animals placed in a flow-through respirometry chamber from measured gas concentrations. For large-scale experiments (e.g., long scans with high sampling rate) that have many uncertainties (e.g., noise in the observations or an unknown impulse response function), this is a computationally challenging inverse problem. We first describe various computational tools that can be used for respirometry reconstruction and uncertainty quantification when the impulse response function is known. Then, we address the more challenging problem where the impulse response function is not known or only partially known. We describe nonlinear optimization methods for reconstruction, where both the unknown model parameters and the unknown signal are reconstructed simultaneously. Numerical experiments show the benefits and potential impacts of these methods in respirometry. 
    more » « less

    Observations of pulsar scintillation are among the few astrophysical probes of very small-scale (≲ 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 overpressurized plasma in local ISM inhomogeneities, such as supernova remnants, H ii regions, and bow shocks. Here we present a survey of eight pulsars conducted at the Five-hundred-metre 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~{{\ \rm per\ cent}}$ 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 ∼105 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.

    more » « less