skip to main content

Title: Deconvolving Pulsar Signals with Cyclic Spectroscopy: A Systematic Evaluation
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 more » 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. « less
; ; ; ; ; ; ; ; ; ;
Award ID(s):
2020265 2009759
Publication Date:
Journal Name:
The Astrophysical Journal
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We have used X-ray data from the Neutron Star Interior Composition Explorer (NICER) to search for long-timescale temporal correlations (“red noise”) in the pulse times of arrival (TOAs) from the millisecond pulsars PSR J1824−2452A and PSR B1937+21. These data more closely track intrinsic noise because X-rays are unaffected by the radio-frequency-dependent propagation effects of the interstellar medium. Our search yields strong evidence (natural log Bayes factor of 9.634 ± 0.016) for red noise in PSR J1824−2452A, but the search is inconclusive for PSR B1937+21. In the interest of future X-ray missions, we devise and implement a method to simulatemore »longer and higher-precision X-ray data sets to determine the timing baseline necessary to detect red noise. We find that the red noise in PSR B1937+21 can be reliably detected in a 5 yr mission with a TOA error of 2 μ s and an observing cadence of 20 observations per month compared to the 5 μ s TOA error and 11 observations per month that NICER currently achieves in PSR B1937+21. We investigate detecting red noise in PSR B1937+21 with other combinations of observing cadences and TOA errors. We also find that time-correlated red noise commensurate with an injected stochastic gravitational-wave background having an amplitude of A GWB = 2 × 10 −15 and spectral index of timing residuals of γ GWB = 13/3 can be detected in a pulsar with similar TOA precision to PSR B1937+21. This is with no additional red noise in a 10 yr mission that observes the pulsar 15 times per month and has an average TOA error of 1 μ s.« less
  2. Abstract Over the past few decades, the measurement precision of some pulsar timing experiments has advanced from ∼10 μ s to ∼10 ns, revealing many subtle phenomena. Such high precision demands both careful data handling and sophisticated timing models to avoid systematic error. To achieve these goals, we present PINT ( P INT I s N ot T empo3 ), a high-precision Python pulsar timing data analysis package, which is hosted on GitHub and available on the Python Package Index (PyPI) as pint-pulsar . PINT is well tested, validated, object oriented, and modular, enabling interactive data analysis and providing anmore »extensible and flexible development platform for timing applications. It utilizes well-debugged public Python packages (e.g., the N um P y and A stropy libraries) and modern software development schemes (e.g., version control and efficient development with git and GitHub) and a continually expanding test suite for improved reliability, accuracy, and reproducibility. PINT is developed and implemented without referring to, copying, or transcribing the code from other traditional pulsar timing software packages (e.g., Tempo / Tempo2 ) and therefore provides a robust tool for cross-checking timing analyses and simulating pulse arrival times. In this paper, we describe the design, use, and validation of PINT , and we compare timing results between it and Tempo and Tempo2 .« less
  3. Abstract Pulsar timing is a process of iteratively fitting pulse arrival times to constrain the spindown, astrometric, and possibly binary parameters of a pulsar, by enforcing integer numbers of pulsar rotations between the arrival times. Phase connection is the process of unambiguously determining those rotation numbers between the times of arrival while determining a pulsar timing solution. Pulsar timing currently requires a manual process of step-by-step phase connection performed by individuals. In an effort to quantify and streamline this process, we created the Algorithmic Pulsar Timer (APT), an algorithm that can accurately phase connect and time isolate pulsars. Using themore »statistical F-test and knowledge of parameter uncertainties and covariances, the algorithm decides what new data to include in a fit, when to add additional timing parameters, and which model to attempt in subsequent iterations. Using these tools, the algorithm can phase-connect timing data that previously required substantial manual effort. We tested the algorithm on 100 simulated systems, with a 99% success rate. APT combines statistical tests and techniques with a logical decision-making process, very similar to the manual one used by pulsar astronomers for decades, and some computational brute force, to automate the often tricky process of isolated pulsar phase connection, setting the foundation for automated fitting of binary pulsar systems.« less
  4. Abstract Detecting continuous nanohertz gravitational waves (GWs) generated by individual close binaries of supermassive black holes (CB-SMBHs) is one of the primary objectives of pulsar timing arrays (PTAs). The detection sensitivity is slated to increase significantly as the number of well-timed millisecond pulsars will increase by more than an order of magnitude with the advent of next-generation radio telescopes. Currently, the Bayesian analysis pipeline using parallel tempering Markov Chain Monte Carlo has been applied in multiple studies for CB-SMBH searches, but it may be challenged by the high dimensionality of the parameter space for future large-scale PTAs. One solution ismore »to reduce the dimensionality by maximizing or marginalizing over uninformative parameters semianalytically, but it is not clear whether this approach can be extended to more complex signal models without making overly simplified assumptions. Recently, the method of diffusive nested (DNest) sampling has shown capability in coping with high dimensionality and multimodality effectively in Bayesian analysis. In this paper, we apply DNest to search for continuous GWs in simulated pulsar timing residuals and find that it performs well in terms of accuracy, robustness, and efficiency for a PTA including  ( 10 2 ) pulsars. DNest also allows a simultaneous search of multiple sources elegantly, which demonstrates its scalability and general applicability. Our results show that it is convenient and also highly beneficial to include DNest in current toolboxes of PTA analysis.« less
  5. Abstract While observations of many high-precision radio pulsars of order ≲1 μ s across the sky are needed for the detection and characterization of a stochastic background of low-frequency gravitational waves (GWs), sensitivity to single sources of GWs requires even higher timing precision. The Argentine Institute of Radio Astronomy (IAR; Instituto Argentino de Radioastronomía) has begun observations of the brightest known millisecond pulsar, J0437−4715. Even though the two antennas are smaller than other single-dish telescopes previously used for pulsar timing array (PTA) science, the IAR’s capability to monitor this pulsar daily, coupled with the pulsar’s brightness, allows for high-precision measurementsmore »of pulse-arrival time. While upgrades of the facility are currently underway, we show that modest improvements beyond current plans will provide IAR with unparalleled sensitivity to this pulsar. The most stringent upper limits on single GW sources come from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). Observations of PSR J0437−4715 will provide a significant sensitivity increase in NANOGrav’s “blind spot” in the sky where fewer pulsars are currently being observed. With state-of-the-art instrumentation installed, we estimate the array’s sensitivity will improve by a factor of ≈2–4 over 10 yr for 20% of the sky with the inclusion of this pulsar, as compared to a static version of the PTA used in NANOGrav’s most recent limits. More modest instrumentation results in factors of ≈1.4–3. We identify four other candidate pulsars as suitable for inclusion in PTA efforts. International PTA efforts will also benefit from inclusion of these data, given the potential achievable sensitivity.« less