More Like this

1. Quality over quantity: Optimizing pulsar timing array analysis for stochastic and continuous gravitational wave signals

   https://doi.org/10.1093/mnras/stac3237  

   Speri, Lorenzo ; Porayko, Nataliya K. ; Falxa, Mikel ; Chen, Siyuan ; Gair, Jonathan R. ; Sesana, Alberto ; Taylor, Stephen R. ( November 2022 , Monthly Notices of the Royal Astronomical Society)

   The search for gravitational waves using Pulsar Timing Arrays (PTAs) is a computationally expensive complex analysis that involves source-specific noise studies. As more pulsars are added to the arrays, this stage of PTA analysis will become increasingly challenging. Therefore, optimizing the number of included pulsars is crucial to reduce the computational burden of data analysis. Here, we present a suite of methods to rank pulsars for use within the scope of PTA analysis. First, we use the maximization of the signal-to-noise ratio as a proxy to select pulsars. With this method, we target the detection of stochastic and continuous gravitational wave signals. Next, we present a ranking that minimizes the coupling between spatial correlation signatures, namely monopolar, dipolar, and Hellings & Downs correlations. Finally, we also explore how to combine these two methods. We test these approaches against mock data using frequentist and Bayesian hypothesis testing. For equal-noise pulsars, we find that an optimal selection leads to an increase in the log-Bayes factor two times steeper than a random selection for the hypothesis test of a gravitational wave background versus a common uncorrelated red noise process. For the same test but for a realistic European PTA (EPTA) data set, a subset of 25 pulsars selected out of 40 can provide a log-likelihood ratio that is 89 % of the total, implying that an optimally selected subset of pulsars can yield results comparable to those obtained from the whole array. We expect these selection methods to play a crucial role in future PTA data combinations.

    

   more » « less
2. The NANOGrav 12.5 yr Data Set: A Computationally Efficient Eccentric Binary Search Pipeline and Constraints on an Eccentric Supermassive Binary Candidate in 3C 66B

   https://doi.org/10.3847/1538-4357/ad1f61  

   Agazie, Gabriella ; Arzoumanian, Zaven ; Baker, Paul T. ; Bécsy, Bence ; Blecha, Laura ; Blumer, Harsha ; Brazier, Adam ; Brook, Paul R. ; Burke-Spolaor, Sarah ; Casey-Clyde, J. Andrew ; et al ( March 2024 , The Astrophysical Journal)

   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.

    

   more » « less
3. Multi-messenger Approaches to Supermassive Black Hole Binary Detection and Parameter Estimation. II. Optimal Strategies for a Pulsar Timing Array

   https://doi.org/10.3847/1538-4357/acb492  

   Liu, Tingting ; Cohen, Tyler ; McGrath, Casey ; Demorest, Paul B. ; Vigeland, Sarah J. ( March 2023 , The Astrophysical Journal)

   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.

    

   more » « less
4. The NANOGrav 15 yr Data Set: Detector Characterization and Noise Budget

   https://doi.org/10.3847/2041-8213/acda88  

   Agazie, Gabriella ; Anumarlapudi, Akash ; Archibald, Anne M. ; Arzoumanian, Zaven ; Baker, Paul T. ; Bécsy, Bence ; Blecha, Laura ; Brazier, Adam ; Brook, Paul R. ; Burke-Spolaor, Sarah ; et al ( June 2023 , The Astrophysical Journal Letters)

   Pulsar timing arrays (PTAs) are galactic-scale gravitational wave (GW) detectors. Each individual arm, composed of a millisecond pulsar, a radio telescope, and a kiloparsecs-long path, differs in its properties but, in aggregate, can be used to extract low-frequency GW signals. We present a noise and sensitivity analysis to accompany the NANOGrav 15 yr data release and associated papers, along with an in-depth introduction to PTA noise models. As a first step in our analysis, we characterize each individual pulsar data set with three types of white-noise parameters and two red-noise parameters. These parameters, along with the timing model and, particularly, a piecewise-constant model for the time-variable dispersion measure, determine the sensitivity curve over the low-frequency GW band we are searching. We tabulate information for all of the pulsars in this data release and present some representative sensitivity curves. We then combine the individual pulsar sensitivities using a signal-to-noise ratio statistic to calculate the global sensitivity of the PTA to a stochastic background of GWs, obtaining a minimum noise characteristic strain of 7 × 10⁻¹⁵at 5 nHz. A power-law-integrated analysis shows rough agreement with the amplitudes recovered in NANOGrav’s 15 yr GW background analysis. While our phenomenological noise model does not model all known physical effects explicitly, it provides an accurate characterization of the noise in the data while preserving sensitivity to multiple classes of GW signals.

    

   more » « less
5. Precision Timing of PSR J0437–4715 with the IAR Observatory and Implications for Low-frequency Gravitational Wave Source Sensitivity

   https://doi.org/10.3847/1538-4357/abeb64  

   Lam, M. T. ; Hazboun, J. S. ( April 2021 , The Astrophysical Journal)

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

   more » « less