skip to main content


Title: Redeveloping a CLEAN Deconvolution Algorithm for Scatter-broadened Radio Pulsar Signals
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
Award ID(s):
2020265 2330164
PAR ID:
10490825
Author(s) / Creator(s):
;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
962
Issue:
2
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 131
Size(s):
Article No. 131
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We report the discovery of a young, highly scattered pulsar in a search for highly circularly polarized radio sources as part of the Australian Square Kilometre Array Pathfinder Variables and Slow Transients survey. In follow-up observations with the Parkes radio telescope, Murriyang, we identified PSR J1032−5804 and measured a period of 78.7 ms, a dispersion measure of 819 ± 4 pc cm−3, a rotation measure of −2000 ± 1 rad m−2, and a characteristic age of 34.6 kyr. We found a pulse scattering timescale at 3 GHz of ∼22 ms, implying a timescale at 1 GHz of ∼3845 ms, which is the third most scattered pulsar known and explains its nondetection in previous pulsar surveys. We discuss the identification of a possible pulsar wind nebula and supernova remnant in the pulsar’s local environment by analyzing the pulsar spectral energy distribution and the surrounding extended emission from multiwavelength images. Our result highlights the possibility of identifying extremely scattered pulsars from radio continuum images. Ongoing and future large-scale radio continuum surveys will offer us an unprecedented opportunity to find more extreme pulsars (e.g., highly scattered, highly intermittent, and highly accelerated), which will enhance our understanding of the characteristics of pulsars and the interstellar medium.

     
    more » « less
  2. ABSTRACT

    We present radio pulsar emission beam analyses and models with the primary intent of examining pulsar beam geometry and physics over the broadest band of radio frequencies reasonably obtainable. We consider a set of well-studied pulsars that lie within the Arecibo sky. These pulsars stand out for the broad frequency range over which emission is detectable, and have been extensively observed at frequencies up to 4.5 GHz and down to below 100 MHz. We utilize published profiles to quantify a more complete picture of the frequency evolution of these pulsars using the core/double-cone emission beam model as our classification framework. For the low-frequency observations, we take into account measured scattering time-scales to infer intrinsic versus scatter broadening of the pulse profile. Lastly, we discuss the populational trends of the core/conal class profiles with respect to intrinsic parameters. We demonstrate that for this subpopulation of pulsars, core and conal dominated profiles cluster together into two roughly segregated $P{\!-\!}\dot{P}$ populations, lending credence to the proposal that an evolution in the pair-formation geometries is responsible for core/conal emission and other emission effects such as nulling and mode changing.

     
    more » « less
  3. ABSTRACT

    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
  4. 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
  5. Abstract

    In 2017 the Event Horizon Telescope (EHT) observed the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz (λ= 1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT data sets. The high observing frequency means that pulsars—which typically exhibit steep emission spectra—are expected to be very faint. However, it also negates pulse scattering, an effect that could hinder pulsar detections in the Galactic center. Additionally, magnetars or a secondary inverse Compton emission could be stronger at millimeter wavelengths than at lower frequencies. We present a search for pulsars close to Sgr A* using the data from the three most sensitive stations in the EHT 2017 campaign: the Atacama Large Millimeter/submillimeter Array, the Large Millimeter Telescope, and the IRAM 30 m Telescope. We apply three detection methods based on Fourier-domain analysis, the fast folding algorithm, and single-pulse searches targeting both pulsars and burst-like transient emission. We use the simultaneity of the observations to confirm potential candidates. No new pulsars or significant bursts were found. Being the first pulsar search ever carried out at such high radio frequencies, we detail our analysis methods and give a detailed estimation of the sensitivity of the search. We conclude that the EHT 2017 observations are only sensitive to a small fraction (≲2.2%) of the pulsars that may exist close to Sgr A*, motivating further searches for fainter pulsars in the region.

     
    more » « less