Time-domain Implementation of the Optimal Cross-Correlation Statistic for Stochastic Gravitational-Wave Background Searches in Pulsar Timing Data

Supermassive black hole binaries, cosmic strings, relic gravitational waves from inflation, and first order phase transitions in the early universe are expected to contribute to a stochastic background of gravitational waves in the 10^(-9) Hz-10^(-7) Hz frequency band. Pulsar timing arrays (PTAs) exploit the high precision timing of radio pulsars to detect signals at such frequencies. Here we present a time-domain implementation of the optimal cross-correlation statistic for stochastic background searches in PTA data. Due to the irregular sampling typical of PTA data as well as the use of a timing model to predict the times-of-arrival of radio pulses, time-domain methods are better suited for gravitational wave data analysis of such data. We present a derivation of the optimal cross-correlation statistic starting from the likelihood function, a method to produce simulated stochastic background signals, and a rigorous derivation of the scaling laws for the signal-to-noise ratio of the cross-correlation statistic in the two relevant PTA regimes: the weak signal limit where instrumental noise dominates over the gravitational wave signal at all frequencies, and a second regime where the gravitational wave signal dominates at the lowest frequencies.Comment: 12 pages, 3 figures, submitted to Physical Review

PINT: A Modern Software Package for Pulsar Timing

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 (PINT Is Not Tempo3), 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 NumPy and Astropy 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

Your: Your Unified Reader

The advancement in signal processing and GPU based systems has enabled new transient detectors at various telescopes to perform much more sensitive searches than their predecessors. Typically the data output from the telescopes is in one of the two commonly used formats: psrfits and Sigproc filterbank. Software developed for transient searches often only works with one of these two formats, limiting their general applicability. Therefore, researchers have to write custom scripts to read/write the data in their format of choice before they can begin any data analysis relevant for their research. \textsc{Your} (Your Unified Reader) is a python-based library that unifies the data processing across multiple commonly used formats. \textsc{Your} implements a user-friendly interface to read and write in the data format of choice. It also generates unified metadata corresponding to the input data file for a quick understanding of observation parameters and provides utilities to perform common data analysis operations. \textsc{Your} also provides several state-of-the-art radio frequency interference mitigation (RFI) algorithms, which can now be used during any stage of data processing (reading, writing, etc.) to filter out artificial signals.Comment: 3 pages, Published in JOSS, Github: https://github.com/thepetabyteproject/you

SEXTANT X-Ray Pulsar Navigation Demonstration: Additional On-Orbit Results

The Station Explorer for X-ray Timing and Navigation Technology (SEXTANT) is a technology demonstration enhancement to the Neutron-star Interior Composition Explorer (NICER) mission, a NASA Astrophysics Explorer Mission of Opportunity to the International Space Station, launched in June of 2017. In late 2017, SEXTANT successfully completed a first demonstration of in-space and autonomous X-ray pulsar navigation (XNAV). This form of navigation relies on processing faint signals from millisecond pulsars-rapidly rotating neutron stars that appear to pulsate in the X-ray band-and could potentially provide a GPS-like navigation capability applicable throughout the solar-system and beyond. In this work, we briefly review prior SEXTANT results and then present new results focusing on: making use of the high- flux but rotationally unstable Crab pulsar, and using XNAV to estimate position, velocity, and time in the presence of an imperfect local clock

Multi-Messenger Gravitational Wave Searches with Pulsar Timing Arrays: Application to 3C66B Using the NANOGrav 11-year Data Set

When galaxies merge, the supermassive black holes in their centers may form binaries and, during the process of merger, emit low-frequency gravitational radiation in the process. In this paper we consider the galaxy 3C66B, which was used as the target of the first multi-messenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of the source of the source, it has been theorized to contain a supermassive black hole binary. Its apparent 1.05-year orbital period would place the gravitational wave emission directly in the pulsar timing band. Since the first pulsar timing array study of 3C66B, revised models of the source have been published, and timing array sensitivities and techniques have improved dramatically. With these advances, we further constrain the chirp mass of the potential supermassive black hole binary in 3C66B to less than $(1.65\pm0.02) \times 10^9~{M_\odot}$ using data from the NANOGrav 11-year data set. This upper limit provides a factor of 1.6 improvement over previous limits, and a factor of 4.3 over the first search done. Nevertheless, the most recent orbital model for the source is still consistent with our limit from pulsar timing array data. In addition, we are able to quantify the improvement made by the inclusion of source properties gleaned from electromagnetic data to `blind' pulsar timing array searches. With these methods, it is apparent that it is not necessary to obtain exact a priori knowledge of the period of a binary to gain meaningful astrophysical inferences.Comment: 14 pages, 6 figures. Accepted by Ap

A Distant Fast Radio Burst Associated with Its Host Galaxy by the Very Large Array

We present the discovery and subarcsecond localization of a new fast radio burst (FRB) by the Karl G. Jansky Very Large Array (VLA) and realfast search system. The FRB was discovered on 2019 June 14 with a dispersion measure of 959 pc cm⁻³. This is the highest DM of any localized FRB and its measured burst fluence of 0.6 Jy ms is less than nearly all other FRBs. The source is not detected to repeat in 15 hr of VLA observing and 153 hr of CHIME/FRB observing. We describe a suite of statistical and data quality tests we used to verify the significance of the event and its localization precision. Follow-up optical/infrared photometry with Keck and Gemini associate the FRB with a pair of galaxies with r ∼ 23 mag. The false-alarm rate for radio transients of this significance that are associated with a host galaxy is roughly 3×10⁻⁴ hr⁻¹. The two putative host galaxies have similar photometric redshifts of z_(phot) ∼ 0.6, but different colors and stellar masses. Comparing the host distance to that implied by the dispersion measure suggests a modest (~ 50 pc/cm⁻³) electron column density associated with the FRB environment or host galaxy/galaxies

The NANOGrav 11-Year Data Set: Arecibo Observatory Polarimetry And Pulse Microcomponents

We present the polarization pulse profiles for 28 pulsars observed with the Arecibo Observatory by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) timing project at 2.1 GHz, 1.4 GHz, and 430 MHz. These profiles represent some of the most sensitive polarimetric millisecond pulsar profiles to date, revealing the existence of microcomponents (that is, pulse components with peak intensities much lower than the total pulse peak intensity). Although microcomponents have been detected in some pulsars previously, we present microcomponents for PSRs B1937+21, J1713+0747, and J2234+0944 for the first time. These microcomponents can have an impact on pulsar timing, geometry, and flux density determination. We present rotation measures for all 28 pulsars, determined independently at different observation frequencies and epochs, and find the Galactic magnetic fields derived from these rotation measures to be consistent with current models. These polarization profiles were made using measurement equation template matching, which allows us to generate the polarimetric response of the Arecibo Observatory on an epoch-by-epoch basis. We use this method to describe its time variability, and find that the polarimetric responses of the Arecibo Observatory's 1.4 and 2.1 GHz receivers vary significantly with time.Comment: 41 pages, 20 figure