Search for: All records

Pulsar timing arrays (PTAs) are Galactic-scale gravitational wave (GW) detectors consisting of precisely timed pulsars distributed across the sky. Within the decade, PTAs are expected to detect nanohertz GWs emitted by close-separation supermassive black hole binaries (SMBHBs), thereby opening up the low-frequency end of the GW spectrum for science. Individual SMBHBs which power active galactic nuclei are also promising multi-messenger sources; they may be identified via theoretically predicted electromagnetic (EM) signatures and be followed up by PTAs for GW observations. In this work, we study the detection and parameter estimation prospects of a PTA which targets EM-selected SMBHBs. Adopting a simulated Galactic millisecond pulsar population, we envisage three different pulsar timing campaigns which observe three mock sources at different sky locations. We find that an all-sky PTA which times the best pulsars is an optimal and feasible approach to observe EM-selected SMBHBs and measure their source parameters to high precision (i.e., comparable to or better than conventional EM measurements). We discuss the implications of our findings in the context of future PTA experiments with the planned Deep Synoptic Array-2000 and the multi-messenger studies of SMBHBs such as the well-known binary candidate OJ 287.

 

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ν⁻²). 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.

 

Abstract Of the more than 3000 radio pulsars currently known, only ∼300 are in binary systems, and only five of these consist of young pulsars with massive nondegenerate companions. We present the discovery and initial timing, accomplished using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope, of the sixth such binary pulsar, PSR J2108+4516, a 0.577 s radio pulsar in a 269 day orbit of eccentricity 0.09 with a companion of minimum mass 11 M ⊙ . Notably, the pulsar undergoes periods of substantial eclipse, disappearing from the CHIME 400–800 MHz observing band for a large fraction of its orbit, and displays significant dispersion measure and scattering variations throughout its orbit, pointing to the possibility of a circumstellar disk or very dense stellar wind associated with the companion star. Subarcsecond resolution imaging with the Karl G. Jansky Very Large Array unambiguously demonstrates that the companion is a bright, V ≃ 11 OBe star, EM* UHA 138, located at a distance of 3.26(14) kpc. Archival optical observations of EM* UHA 138 approximately suggest a companion mass ranging from 17.5 M ⊙ < M c < 23 M ⊙ , in turn constraining the orbital inclination angle to 50.°3 ≲ i ≲ 58.°3. With further multiwavelength follow-up, PSR J2108+4516 promises to serve as another rare laboratory for the exploration of companion winds, circumstellar disks, and short-term evolution through extended-body orbital dynamics. 

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. 

The radio galaxy 3C 66B has been hypothesized to host a supermassive black hole binary (SMBHB) at its center based on electromagnetic observations. Its apparent 1.05 yr period and low redshift (∼0.02) make it an interesting testbed to search for low-frequency gravitational waves (GWs) using pulsar timing array (PTA) experiments. This source has been subjected to multiple searches for continuous GWs from a circular SMBHB, resulting in progressively more stringent constraints on its GW amplitude and chirp mass. In this paper, we develop a pipeline for performing Bayesian targeted searches for eccentric SMBHBs in PTA data sets, and test its efficacy by applying it to simulated data sets with varying injected signal strengths. We also search for a realistic eccentric SMBHB source in 3C 66B using the NANOGrav 12.5 yr data set employing PTA signal models containing Earth term-only as well as Earth+pulsar term contributions using this pipeline. Due to limitations in our PTA signal model, we get meaningful results only when the initial eccentricitye₀< 0.5 and the symmetric mass ratioη> 0.1. We find no evidence for an eccentric SMBHB signal in our data, and therefore place 95% upper limits on the PTA signal amplitude of 88.1 ± 3.7 ns for the Earth term-only and 81.74 ± 0.86 ns for the Earth+pulsar term searches fore₀< 0.5 andη> 0.1. Similar 95% upper limits on the chirp mass are (1.98 ± 0.05) × 10⁹and (1.81 ± 0.01) × 10⁹M_☉. These upper limits, while less stringent than those calculated from a circular binary search in the NANOGrav 12.5 yr data set, are consistent with the SMBHB model of 3C 66B developed from electromagnetic observations.

 

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has reported evidence for the presence of an isotropic nanohertz gravitational-wave background (GWB) in its 15 yr data set. However, if the GWB is produced by a population of inspiraling supermassive black hole binary (SMBHB) systems, then the background is predicted to be anisotropic, depending on the distribution of these systems in the local Universe and the statistical properties of the SMBHB population. In this work, we search for anisotropy in the GWB using multiple methods and bases to describe the distribution of the GWB power on the sky. We do not find significant evidence of anisotropy. By modeling the angular power distribution as a sum over spherical harmonics (where the coefficients are not bound to always generate positive power everywhere), we find that the Bayesian 95% upper limit on the level of dipole anisotropy is (C_l=1/C_l=0) < 27%. This is similar to the upper limit derived under the constraint of positive power everywhere, indicating that the dipole may be close to the data-informed regime. By contrast, the constraints on anisotropy at higher spherical-harmonic multipoles are strongly prior dominated. We also derive conservative estimates on the anisotropy expected from a random distribution of SMBHB systems using astrophysical simulations conditioned on the isotropic GWB inferred in the 15 yr data set and show that this data set has sufficient sensitivity to probe a large fraction of the predicted level of anisotropy. We end by highlighting the opportunities and challenges in searching for anisotropy in pulsar timing array data.

 

We present the results of a Bayesian search for gravitational wave (GW) memory in the NANOGrav 12.5 yr data set. We find no convincing evidence for any gravitational wave memory signals in this data set. We find a Bayes factor of 2.8 in favor of a model that includes a memory signal and common spatially uncorrelated red noise (CURN) compared to a model including only a CURN. However, further investigation shows that a disproportionate amount of support for the memory signal comes from three dubious pulsars. Using a more flexible red-noise model in these pulsars reduces the Bayes factor to 1.3. Having found no compelling evidence, we go on to place upper limits on the strain amplitude of GW memory events as a function of sky location and event epoch. These upper limits are computed using a signal model that assumes the existence of a common, spatially uncorrelated red noise in addition to a GW memory signal. The median strain upper limit as a function of sky position is approximately 3.3 × 10⁻¹⁴. We also find that there are some differences in the upper limits as a function of sky position centered around PSR J0613−0200. This suggests that this pulsar has some excess noise that can be confounded with GW memory. Finally, the upper limits as a function of burst epoch continue to improve at later epochs. This improvement is attributable to the continued growth of the pulsar timing array.

 

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 an 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 .