Showing 1–50 of 95 results for author: Taylor, S R

Spectral Variance in a Stochastic Gravitational Wave Background From a Binary Population

Authors: William G. Lamb, Stephen R. Taylor

Abstract: A population of compact object binaries emitting gravitational waves that are not individually resolvable will form a stochastic gravitational wave signal. While the expected spectrum over population realizations is well known from Phinney (2001), its higher order moments have not been fully studied before or computed in the case of arbitrary binary evolution. We calculate analytic scaling relatio… ▽ More A population of compact object binaries emitting gravitational waves that are not individually resolvable will form a stochastic gravitational wave signal. While the expected spectrum over population realizations is well known from Phinney (2001), its higher order moments have not been fully studied before or computed in the case of arbitrary binary evolution. We calculate analytic scaling relationships as a function of gravitational-wave frequency for the statistical variance, skewness, and kurtosis of a stochastic gravitational-wave signal over population realizations. If the time derivative of the binary orbital frequency can be expressed as a power-law in frequency, we find that these moment quantities also take the form of power-law relationships. We also develop a numerical population synthesis framework against which we compare our analytic results, finding excellent agreement. These new scaling relationships provide physical context to understanding spectral fluctuations in a gravitational-wave background signal and may provide additional information that can aid in explaining the origin of the nanohertz-frequency signal observed by pulsar timing array campaigns. △ Less

Submitted 18 July, 2024; v1 submitted 8 July, 2024; originally announced July 2024.

Comments: Accepted to ApJL

Host Galaxy Demographics Of Individually Detectable Supermassive Black-hole Binaries with Pulsar Timing Arrays

Authors: Katharine Cella, Stephen R. Taylor, Luke Zoltan Kelley

Abstract: Massive black hole binaries (MBHBs) produce gravitational waves (GWs) that are detectable with pulsar timing arrays. We determine the properties of the host galaxies of simulated MBHBs at the time they are producing detectable GW signals. The population of MBHB systems we evaluate is from the Illustris cosmological simulations taken in tandem with post processing semi-analytic models of environmen… ▽ More Massive black hole binaries (MBHBs) produce gravitational waves (GWs) that are detectable with pulsar timing arrays. We determine the properties of the host galaxies of simulated MBHBs at the time they are producing detectable GW signals. The population of MBHB systems we evaluate is from the Illustris cosmological simulations taken in tandem with post processing semi-analytic models of environmental factors in the evolution of binaries. Upon evolving to the GW frequency regime accessible by pulsar timing arrays, we calculate the detection probability of each system using a variety of different values for pulsar noise characteristics in a plausible near-future International Pulsar Timing Array dataset. We find that detectable systems have host galaxies that are clearly distinct from the overall binary population and from most galaxies in general. With conservative noise factors, we find that host stellar metallicity, for example, peaks at twice solar metallicity as opposed to the total population of galaxies which peaks at ~0.6 solar metallicity. Additionally, the most detectable systems are much brighter in magnitude and more red in color than the overall population, indicating their likely identity as large ellipticals with diminished star formation. These results can be used to develop effective search strategies for identifying host galaxies and electromagnetic counterparts following GW detection by pulsar timing arrays. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: Submitted to CQG Focus Issue

Spatial and Spectral Characterization of the Gravitational-wave Background with the PTA Optimal Statistic

Authors: Kyle A. Gersbach, Stephen R. Taylor, Patrick M. Meyers, Joseph D. Romano

Abstract: Pulsar timing arrays (PTAs) have made tremendous progress and are now showing strong evidence for the gravitational-wave background (GWB). Further probing the origin and characteristics of the GWB will require more generalized analysis techniques. Bayesian methods are most often used but can be computationally expensive. On the other hand, frequentist methods, like the PTA Optimal Statistic (OS),… ▽ More Pulsar timing arrays (PTAs) have made tremendous progress and are now showing strong evidence for the gravitational-wave background (GWB). Further probing the origin and characteristics of the GWB will require more generalized analysis techniques. Bayesian methods are most often used but can be computationally expensive. On the other hand, frequentist methods, like the PTA Optimal Statistic (OS), are more computationally efficient and can produce results that are complementary to Bayesian methods, allowing for stronger statistical cases to be built from a confluence of different approaches. In this work we expand the capabilities of the OS through a technique we call the Per-Frequency Optimal Statistic (PFOS). The PFOS removes the underlying power-law assumption inherent in previous implementations of the OS, and allows one to estimate the GWB spectrum in a frequency-by-frequency manner. We have also adapted a recent generalization from the OS pipeline into the PFOS, making it capable of accurately characterizing the spectrum in the intermediate and strong GW signal regimes using only a small fraction of the necessary computational resources when compared with fully-correlated Bayesian methods, while also empowering many new types of analyses not possible before. We find that even in the strong GW signal regime, where the GWB dominates over noise in all frequencies, the injected value of the signal lies within the 50th-percentile of the PFOS uncertainty distribution in 41-45% of simulations, remaining 3$σ$-consistent with unbiased estimation. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: 23 pages, 9 figures,submitted to Physical Review D

Identifying Host Galaxies of Supermassive Black Hole Binaries Found by PTAs

Authors: Polina Petrov, Stephen R. Taylor, Maria Charisi, Chung-Pei Ma

Abstract: Supermassive black hole binaries (SMBHBs) present us with exciting opportunities for multi-messenger science. These systems are thought to form naturally in galaxy mergers and therefore have the potential to produce electromagnetic (EM) radiation as well as gravitational waves (GWs) detectable with pulsar timing arrays (PTAs). Once GWs from individually resolved SMBHBs are detected, the identifica… ▽ More Supermassive black hole binaries (SMBHBs) present us with exciting opportunities for multi-messenger science. These systems are thought to form naturally in galaxy mergers and therefore have the potential to produce electromagnetic (EM) radiation as well as gravitational waves (GWs) detectable with pulsar timing arrays (PTAs). Once GWs from individually resolved SMBHBs are detected, the identification of the host galaxy will be a major challenge due to the ambiguity in possible EM signatures and the poor localization capability of PTAs. In order to aid EM observations in choosing which sources to follow up, we attempt to quantify the number of plausible hosts in both realistic and idealistic scenarios. We outline a host galaxy identification pipeline that injects a single-source GW signal into a simulated PTA dataset, uses production-level techniques to recover the signal, quantifies the localization region and number of galaxies contained therein, and finally imposes cuts on the galaxies using the binary parameters estimated from the GW search. In an ideal case, we find that the 90% credible areas span 29 deg^2 to 241 deg^2, containing about 14 to 341 galaxies. After cuts, the number of galaxies remaining ranges from 22 at worst to 1 (the true host) at best. In a more realistic case, if the signal is sufficiently localized, the sky areas range from 287 deg^2 to 530 deg^2 and enclose about 285 to 1238 galaxies. After cuts, the number of galaxies is 397 at worst and 27 at best. While the signal-to-noise ratio is the primary determinant of the localization area of a given source, we find that the size of the area is also influenced by the proximity of nearby pulsars on the sky and the chirp mass of the source. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: Submitted to ApJ

Overlap reduction functions for pulsar timing arrays and astrometry

Authors: Keisuke Inomata, Marc Kamionkowski, Celia M. Toral, Stephen R. Taylor

Abstract: We present an efficient technique for calculating the angular two-point correlation functions (or ''overlap reduction functions'') induced by gravitational waves in both the pulse arrival times of pulsars and in the angular deflections of distant sources. In the most general case, there are six auto- and cross-correlations for the pulse arrival times and the two components of the angular deflectio… ▽ More We present an efficient technique for calculating the angular two-point correlation functions (or ''overlap reduction functions'') induced by gravitational waves in both the pulse arrival times of pulsars and in the angular deflections of distant sources. In the most general case, there are six auto- and cross-correlations for the pulse arrival times and the two components of the angular deflection. We provide results for spin-2 (i.e., general-relativistic) gravitational waves as well as the spin-1 modes that may arise in alternative-gravity theories. These calculations can be easily implemented for future analysis or study, and we provide code to do so. △ Less

Submitted 31 May, 2024; originally announced June 2024.

Comments: 20 pages, 6 figures, code available at https://github.com/KeisukeInomata0/pyORFs

The NANOGrav 15 yr Data Set: Looking for Signs of Discreteness in the Gravitational-wave Background

Authors: Gabriella Agazie, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Lucas Brown, Sarah Burke-Spolaor, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Heling Deng, Timothy Dolch, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca, Gabriel E. Freedman, Nate Garver-Daniels , et al. (58 additional authors not shown)

Abstract: The cosmic merger history of supermassive black hole binaries (SMBHBs) is expected to produce a low-frequency gravitational wave background (GWB). Here we investigate how signs of the discrete nature of this GWB can manifest in pulsar timing arrays through excursions from, and breaks in, the expected $f_{\mathrm{GW}}^{-2/3}$ power-law of the GWB strain spectrum. To do this, we create a semi-analyt… ▽ More The cosmic merger history of supermassive black hole binaries (SMBHBs) is expected to produce a low-frequency gravitational wave background (GWB). Here we investigate how signs of the discrete nature of this GWB can manifest in pulsar timing arrays through excursions from, and breaks in, the expected $f_{\mathrm{GW}}^{-2/3}$ power-law of the GWB strain spectrum. To do this, we create a semi-analytic SMBHB population model, fit to NANOGrav's 15 yr GWB amplitude, and with 1,000 realizations we study the populations' characteristic strain and residual spectra. Comparing our models to the NANOGrav 15 yr spectrum, we find two interesting excursions from the power-law. The first, at $2 \; \mathrm{nHz}$, is below our GWB realizations with $p$-value significance $p = 0.05$ to $0.06$ ($\approx 1.8 σ- 1.9 σ$). The second, at $16 \; \mathrm{nHz}$, is above our GWB realizations with $p = 0.04$ to $0.15$ ($\approx 1.4 σ- 2.1 σ$). We explore the properties of a loud SMBHB which could cause such an excursion. Our simulations also show that the expected number of SMBHBs decreases by three orders of magnitude, from $\sim 10^6$ to $\sim 10^3$, between $2\; \mathrm{nHz}$ and $20 \; \mathrm{nHz}$. This causes a break in the strain spectrum as the stochasticity of the background breaks down at $26^{+28}_{-19} \; \mathrm{nHz}$, consistent with predictions pre-dating GWB measurements. The diminished GWB signal from SMBHBs at frequencies above the $26$~nHz break opens a window for PTAs to detect continuous GWs from individual SMBHBs or GWs from the early universe. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: 10 pages, 8 figures, 1 appendix, submitted to ApJ

Fast Parameter Inference on Pulsar Timing Arrays with Normalizing Flows

Authors: David Shih, Marat Freytsis, Stephen R. Taylor, Jeff A. Dror, Nolan Smyth

Abstract: Pulsar timing arrays (PTAs) perform Bayesian posterior inference with expensive MCMC methods. Given a dataset of ~10-100 pulsars and O(10^3) timing residuals each, producing a posterior distribution for the stochastic gravitational wave background (SGWB) can take days to a week. The computational bottleneck arises because the likelihood evaluation required for MCMC is extremely costly when conside… ▽ More Pulsar timing arrays (PTAs) perform Bayesian posterior inference with expensive MCMC methods. Given a dataset of ~10-100 pulsars and O(10^3) timing residuals each, producing a posterior distribution for the stochastic gravitational wave background (SGWB) can take days to a week. The computational bottleneck arises because the likelihood evaluation required for MCMC is extremely costly when considering the dimensionality of the search space. Fortunately, generating simulated data is fast, so modern simulation-based inference techniques can be brought to bear on the problem. In this paper, we demonstrate how conditional normalizing flows trained on simulated data can be used for extremely fast and accurate estimation of the SGWB posteriors, reducing the sampling time from weeks to a matter of seconds. △ Less

Submitted 18 October, 2023; originally announced October 2023.

Comments: 8 pages, 3 figures

The NANOGrav 15-year data set: Search for Transverse Polarization Modes in the Gravitational-Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Jeremy Baier, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Dallas DeGan, Paul B. Demorest , et al. (74 additional authors not shown)

Abstract: Recently we found compelling evidence for a gravitational wave background with Hellings and Downs (HD) correlations in our 15-year data set. These correlations describe gravitational waves as predicted by general relativity, which has two transverse polarization modes. However, more general metric theories of gravity can have additional polarization modes which produce different interpulsar correl… ▽ More Recently we found compelling evidence for a gravitational wave background with Hellings and Downs (HD) correlations in our 15-year data set. These correlations describe gravitational waves as predicted by general relativity, which has two transverse polarization modes. However, more general metric theories of gravity can have additional polarization modes which produce different interpulsar correlations. In this work we search the NANOGrav 15-year data set for evidence of a gravitational wave background with quadrupolar Hellings and Downs (HD) and Scalar Transverse (ST) correlations. We find that HD correlations are the best fit to the data, and no significant evidence in favor of ST correlations. While Bayes factors show strong evidence for a correlated signal, the data does not strongly prefer either correlation signature, with Bayes factors $\sim 2$ when comparing HD to ST correlations, and $\sim 1$ for HD plus ST correlations to HD correlations alone. However, when modeled alongside HD correlations, the amplitude and spectral index posteriors for ST correlations are uninformative, with the HD process accounting for the vast majority of the total signal. Using the optimal statistic, a frequentist technique that focuses on the pulsar-pair cross-correlations, we find median signal-to-noise-ratios of 5.0 for HD and 4.6 for ST correlations when fit for separately, and median signal-to-noise-ratios of 3.5 for HD and 3.0 for ST correlations when fit for simultaneously. While the signal-to-noise-ratios for each of the correlations are comparable, the estimated amplitude and spectral index for HD are a significantly better fit to the total signal, in agreement with our Bayesian analysis. △ Less

Submitted 18 October, 2023; originally announced October 2023.

Comments: 11 pages, 5 figures

The NANOGrav 12.5-year data set: A computationally efficient eccentric binary search pipeline and constraints on an eccentric supermassive binary candidate in 3C 66B

Authors: Gabriella Agazie, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Harsha Blumer, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Belinda D. Cheeseboro, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Lankeswar Dey, Timothy Dolch, Justin A. Ellis, Robert D. Ferdman, Elizabeth C. Ferrara , et al. (63 additional authors not shown)

Abstract: 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-year period and low redshift ($\sim0.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 contin… ▽ More 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-year period and low redshift ($\sim0.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 on 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-year 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 eccentricity $e_0<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\pm3.7$ ns for the Earth term-only and $81.74\pm0.86$ ns for the Earth+Pulsar term searches for $e_0<0.5$ and $η>0.1$. Similar 95% upper limits on the chirp mass are $(1.98 \pm 0.05) \times 10^9\,M_{\odot}$ and $(1.81 \pm 0.01) \times 10^9\,M_{\odot}$. These upper limits, while less stringent than those calculated from a circular binary search in the NANOGrav 12.5-year data set, are consistent with the SMBHB model of 3C 66B developed from electromagnetic observations. △ Less

Submitted 15 January, 2024; v1 submitted 29 September, 2023; originally announced September 2023.

Comments: 27 Pages, 10 Figures, 1 Table, Accepted for publication in ApJ

How to Detect an Astrophysical Nanohertz Gravitational-Wave Background

Authors: Bence Bécsy, Neil J. Cornish, Patrick M. Meyers, Luke Zoltan Kelley, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Katerina Chatziioannou, Tyler Cohen, James M. Cordes, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch , et al. (71 additional authors not shown)

Abstract: Analysis of pulsar timing data have provided evidence for a stochastic gravitational wave background in the nHz frequency band. The most plausible source of such a background is the superposition of signals from millions of supermassive black hole binaries. The standard statistical techniques used to search for such a background and assess its significance make several simplifying assumptions, nam… ▽ More Analysis of pulsar timing data have provided evidence for a stochastic gravitational wave background in the nHz frequency band. The most plausible source of such a background is the superposition of signals from millions of supermassive black hole binaries. The standard statistical techniques used to search for such a background and assess its significance make several simplifying assumptions, namely: i) Gaussianity; ii) isotropy; and most often iii) a power-law spectrum. However, a stochastic background from a finite collection of binaries does not exactly satisfy any of these assumptions. To understand the effect of these assumptions, we test standard analysis techniques on a large collection of realistic simulated datasets. The dataset length, observing schedule, and noise levels were chosen to emulate the NANOGrav 15-year dataset. Simulated signals from millions of binaries drawn from models based on the Illustris cosmological hydrodynamical simulation were added to the data. We find that the standard statistical methods perform remarkably well on these simulated datasets, despite their fundamental assumptions not being strictly met. They are able to achieve a confident detection of the background. However, even for a fixed set of astrophysical parameters, different realizations of the universe result in a large variance in the significance and recovered parameters of the background. We also find that the presence of loud individual binaries can bias the spectral recovery of the background if we do not account for them. △ Less

Submitted 1 December, 2023; v1 submitted 8 September, 2023; originally announced September 2023.

Comments: 14 pages, 8 figures, version matching published paper

Journal ref: ApJ 959 9 (2023)

Comparing recent PTA results on the nanohertz stochastic gravitational wave background

Authors: The International Pulsar Timing Array Collaboration, G. Agazie, J. Antoniadis, A. Anumarlapudi, A. M. Archibald, P. Arumugam, S. Arumugam, Z. Arzoumanian, J. Askew, S. Babak, M. Bagchi, M. Bailes, A. -S. Bak Nielsen, P. T. Baker, C. G. Bassa, A. Bathula, B. Bécsy, A. Berthereau, N. D. R. Bhat, L. Blecha, M. Bonetti, E. Bortolas, A. Brazier, P. R. Brook, M. Burgay , et al. (220 additional authors not shown)

Abstract: The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTA… ▽ More The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTAs that constitute the International Pulsar Timing Array (IPTA). We show that despite making different modeling choices, there is no significant difference in the GWB parameters that are measured by the different PTAs, agreeing within $1σ$. The pulsar noise parameters are also consistent between different PTAs for the majority of the pulsars included in these analyses. We bridge the differences in modeling choices by adopting a standardized noise model for all pulsars and PTAs, finding that under this model there is a reduction in the tension in the pulsar noise parameters. As part of this reanalysis, we "extended" each PTA's data set by adding extra pulsars that were not timed by that PTA. Under these extensions, we find better constraints on the GWB amplitude and a higher signal-to-noise ratio for the Hellings and Downs correlations. These extensions serve as a prelude to the benefits offered by a full combination of data across all pulsars in the IPTA, i.e., the IPTA's Data Release 3, which will involve not just adding in additional pulsars, but also including data from all three PTAs where any given pulsar is timed by more than as single PTA. △ Less

Submitted 1 September, 2023; originally announced September 2023.

Comments: 21 pages, 9 figures, submitted to ApJ

The NANOGrav 12.5-year Data Set: Search for Gravitational Wave Memory

Authors: Gabriella Agazie, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Harsha Blumer, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Dallas DeGan, Paul B. Demorest, Timothy Dolch, Brendan Drachler, Justin A. Ellis , et al. (65 additional authors not shown)

Abstract: 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 (Bayes factor = 2.8). As such, 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 sign… ▽ More 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 (Bayes factor = 2.8). As such, 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 \times 10^{-14}$. 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 which 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. △ Less

Submitted 25 July, 2023; originally announced July 2023.

Comments: 29 pages, 5 figures

The NANOGrav 15-year Gravitational-Wave Background Analysis Pipeline

Authors: Aaron D. Johnson, Patrick M. Meyers, Paul T. Baker, Neil J. Cornish, Jeffrey S. Hazboun, Tyson B. Littenberg, Joseph D. Romano, Stephen R. Taylor, Michele Vallisneri, Sarah J. Vigeland, Ken D. Olum, Xavier Siemens, Justin A. Ellis, Rutger van Haasteren, Sophie Hourihane, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Bence Bécsy, J. Andrew Casey-Clyde , et al. (71 additional authors not shown)

Abstract: This paper presents rigorous tests of pulsar timing array methods and software, examining their consistency across a wide range of injected parameters and signal strength. We discuss updates to the 15-year isotropic gravitational-wave background analyses and their corresponding code representations. Descriptions of the internal structure of the flagship algorithms \texttt{Enterprise} and \texttt{P… ▽ More This paper presents rigorous tests of pulsar timing array methods and software, examining their consistency across a wide range of injected parameters and signal strength. We discuss updates to the 15-year isotropic gravitational-wave background analyses and their corresponding code representations. Descriptions of the internal structure of the flagship algorithms \texttt{Enterprise} and \texttt{PTMCMCSampler} are given to facilitate understanding of the PTA likelihood structure, how models are built, and what methods are currently used in sampling the high-dimensional PTA parameter space. We introduce a novel version of the PTA likelihood that uses a two-step marginalization procedure that performs much faster when the white noise parameters remain fixed. We perform stringent tests of consistency and correctness of the Bayesian and frequentist analysis software. For the Bayesian analysis, we test prior recovery, injection recovery, and Bayes factors. For the frequentist analysis, we test that the cross-correlation-based optimal statistic, when modified to account for a non-negligible gravitational-wave background, accurately recovers the amplitude of the background. We also summarize recent advances and tests performed on the optimal statistic in the literature from both GWB detection and parameter estimation perspectives. The tests presented here validate current and future analyses of PTA data. △ Less

Submitted 7 July, 2023; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 30 pages, 10 figures, 1 table; Companion paper to "The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background"; For questions or comments, please email [email protected]

The NANOGrav 15-year Data Set: Search for Anisotropy in the Gravitational-Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Brendan Drachler, Elizabeth C. Ferrara, William Fiore , et al. (68 additional authors not shown)

Abstract: 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 dataset. 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… ▽ More 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 dataset. 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, and place a Bayesian $95\%$ upper limit on the level of broadband anisotropy such that $(C_{l>0} / C_{l=0}) < 20\%$. 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 dataset, and show that this dataset 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. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 19 pages, 11 figures; submitted to Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email [email protected]

The NANOGrav 15-year Data Set: Constraints on Supermassive Black Hole Binaries from the Gravitational Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Paul T. Baker, Bence Bécsy, Laura Blecha, Alexander Bonilla, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Katerina Chatziioannou, Belinda D. Cheeseboro, Siyuan Chen, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Curt J. Cutler , et al. (89 additional authors not shown)

Abstract: The NANOGrav 15-year data set shows evidence for the presence of a low-frequency gravitational-wave background (GWB). While many physical processes can source such low-frequency gravitational waves, here we analyze the signal as coming from a population of supermassive black hole (SMBH) binaries distributed throughout the Universe. We show that astrophysically motivated models of SMBH binary popul… ▽ More The NANOGrav 15-year data set shows evidence for the presence of a low-frequency gravitational-wave background (GWB). While many physical processes can source such low-frequency gravitational waves, here we analyze the signal as coming from a population of supermassive black hole (SMBH) binaries distributed throughout the Universe. We show that astrophysically motivated models of SMBH binary populations are able to reproduce both the amplitude and shape of the observed low-frequency gravitational-wave spectrum. While multiple model variations are able to reproduce the GWB spectrum at our current measurement precision, our results highlight the importance of accurately modeling binary evolution for producing realistic GWB spectra. Additionally, while reasonable parameters are able to reproduce the 15-year observations, the implied GWB amplitude necessitates either a large number of parameters to be at the edges of expected values, or a small number of parameters to be notably different from standard expectations. While we are not yet able to definitively establish the origin of the inferred GWB signal, the consistency of the signal with astrophysical expectations offers a tantalizing prospect for confirming that SMBH binaries are able to form, reach sub-parsec separations, and eventually coalesce. As the significance grows over time, higher-order features of the GWB spectrum will definitively determine the nature of the GWB and allow for novel constraints on SMBH populations. △ Less

Submitted 18 July, 2023; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: Accepted by Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email [email protected]. Edited to fix two equation typos (Eq.13 & 21), and minor text typos

The NANOGrav 15-year Data Set: Search for Signals from New Physics

Authors: Adeela Afzal, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Jose Juan Blanco-Pillado, Laura Blecha, Kimberly K. Boddy, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, Maria Charisi, Shami Chatterjee, Katerina Chatziioannou, Belinda D. Cheeseboro, Siyuan Chen, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie , et al. (98 additional authors not shown)

Abstract: The 15-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) shows positive evidence for the presence of a low-frequency gravitational-wave (GW) background. In this paper, we investigate potential cosmological interpretations of this signal, specifically cosmic inflation, scalar-induced GWs, first-order phase transitions, cosmic string… ▽ More The 15-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) shows positive evidence for the presence of a low-frequency gravitational-wave (GW) background. In this paper, we investigate potential cosmological interpretations of this signal, specifically cosmic inflation, scalar-induced GWs, first-order phase transitions, cosmic strings, and domain walls. We find that, with the exception of stable cosmic strings of field theory origin, all these models can reproduce the observed signal. When compared to the standard interpretation in terms of inspiraling supermassive black hole binaries (SMBHBs), many cosmological models seem to provide a better fit resulting in Bayes factors in the range from 10 to 100. However, these results strongly depend on modeling assumptions about the cosmic SMBHB population and, at this stage, should not be regarded as evidence for new physics. Furthermore, we identify excluded parameter regions where the predicted GW signal from cosmological sources significantly exceeds the NANOGrav signal. These parameter constraints are independent of the origin of the NANOGrav signal and illustrate how pulsar timing data provide a new way to constrain the parameter space of these models. Finally, we search for deterministic signals produced by models of ultralight dark matter (ULDM) and dark matter substructures in the Milky Way. We find no evidence for either of these signals and thus report updated constraints on these models. In the case of ULDM, these constraints outperform torsion balance and atomic clock constraints for ULDM coupled to electrons, muons, or gluons. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 74 pages, 31 figures, 4 tables; published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email [email protected]

The NANOGrav 15-Year Data Set: Detector Characterization and Noise Budget

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. Decesar, Paul B. Demorest, Timothy Dolch, Brendan Drachler, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca , et al. (66 additional authors not shown)

Abstract: Pulsar timing arrays (PTAs) are galactic-scale gravitational wave 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 gravitational wave (GW) signals. We present a noise and sensitivity analysis to accompany the NANOGrav 15-year data release and associated… ▽ More Pulsar timing arrays (PTAs) are galactic-scale gravitational wave 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 gravitational wave (GW) signals. We present a noise and sensitivity analysis to accompany the NANOGrav 15-year 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\times 10^{-15}$ at 5 nHz. A power law-integrated analysis shows rough agreement with the amplitudes recovered in NANOGrav's 15-year 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. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 67 pages, 73 figures, 3 tables; published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email [email protected]

The NANOGrav 15-year Data Set: Observations and Timing of 68 Millisecond Pulsars

Authors: Gabriella Agazie, Md Faisal Alam, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Laura Blecha, Victoria Bonidie, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Bence Bécsy, Christopher Chapman, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Brendan Drachler , et al. (75 additional authors not shown)

Abstract: We present observations and timing analyses of 68 millisecond pulsars (MSPs) comprising the 15-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). NANOGrav is a pulsar timing array (PTA) experiment that is sensitive to low-frequency gravitational waves. This is NANOGrav's fifth public data release, including both "narrowband" and "wideband" time-of-arrival… ▽ More We present observations and timing analyses of 68 millisecond pulsars (MSPs) comprising the 15-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). NANOGrav is a pulsar timing array (PTA) experiment that is sensitive to low-frequency gravitational waves. This is NANOGrav's fifth public data release, including both "narrowband" and "wideband" time-of-arrival (TOA) measurements and corresponding pulsar timing models. We have added 21 MSPs and extended our timing baselines by three years, now spanning nearly 16 years for some of our sources. The data were collected using the Arecibo Observatory, the Green Bank Telescope, and the Very Large Array between frequencies of 327 MHz and 3 GHz, with most sources observed approximately monthly. A number of notable methodological and procedural changes were made compared to our previous data sets. These improve the overall quality of the TOA data set and are part of the transition to new pulsar timing and PTA analysis software packages. For the first time, our data products are accompanied by a full suite of software to reproduce data reduction, analysis, and results. Our timing models include a variety of newly detected astrometric and binary pulsar parameters, including several significant improvements to pulsar mass constraints. We find that the time series of 23 pulsars contain detectable levels of red noise, 10 of which are new measurements. In this data set, we find evidence for a stochastic gravitational-wave background. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 90 pages, 74 figures, 6 tables; published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email [email protected]

The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Becsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, Maria Charisi, Shami Chatterjee, Katerina Chatziioannou, Belinda D. Cheeseboro, Siyuan Chen, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Curt J. Cutler, Megan E. DeCesar , et al. (89 additional authors not shown)

Abstract: We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15-year pulsar-timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. The correlations follow the Hellings-Downs pattern expected for a stochastic gravitational-wave background. The presence of such a gravitational-wave background with a power-law-spectr… ▽ More We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15-year pulsar-timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. The correlations follow the Hellings-Downs pattern expected for a stochastic gravitational-wave background. The presence of such a gravitational-wave background with a power-law-spectrum is favored over a model with only independent pulsar noises with a Bayes factor in excess of $10^{14}$, and this same model is favored over an uncorrelated common power-law-spectrum model with Bayes factors of 200-1000, depending on spectral modeling choices. We have built a statistical background distribution for these latter Bayes factors using a method that removes inter-pulsar correlations from our data set, finding $p = 10^{-3}$ (approx. $3σ$) for the observed Bayes factors in the null no-correlation scenario. A frequentist test statistic built directly as a weighted sum of inter-pulsar correlations yields $p = 5 \times 10^{-5} - 1.9 \times 10^{-4}$ (approx. $3.5 - 4σ$). Assuming a fiducial $f^{-2/3}$ characteristic-strain spectrum, as appropriate for an ensemble of binary supermassive black-hole inspirals, the strain amplitude is $2.4^{+0.7}_{-0.6} \times 10^{-15}$ (median + 90% credible interval) at a reference frequency of 1/(1 yr). The inferred gravitational-wave background amplitude and spectrum are consistent with astrophysical expectations for a signal from a population of supermassive black-hole binaries, although more exotic cosmological and astrophysical sources cannot be excluded. The observation of Hellings-Downs correlations points to the gravitational-wave origin of this signal. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 30 pages, 18 figures. Published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email [email protected]

Exploring the Capabilities of Gibbs Sampling in Pulsar Timing Arrays

Authors: Nima Laal, William G Lamb, Joseph D. Romano, Xavier Siemens, Stephen R. Taylor, Rutger van Haasteren

Abstract: We explore the use of Gibbs sampling in estimating the noise properties of individual pulsars and illustrate its effectiveness using the NANOGrav 11-year data set. We find that Gibbs sampling noise modeling (GM) is more efficient than the current standard Bayesian techniques (SM) for single pulsar analyses by yielding model parameter posteriors with average effective-sample-size ratio (GM/SM) of 6… ▽ More We explore the use of Gibbs sampling in estimating the noise properties of individual pulsars and illustrate its effectiveness using the NANOGrav 11-year data set. We find that Gibbs sampling noise modeling (GM) is more efficient than the current standard Bayesian techniques (SM) for single pulsar analyses by yielding model parameter posteriors with average effective-sample-size ratio (GM/SM) of 6 across all parameters and pulsars. Furthermore, the output of GM contains posteriors for the Fourier coefficients that can be used to characterize the underlying red noise process of any pulsar's timing residuals, which are absent in current implementations of SM. Through simulations, we demonstrate the potential for such coefficients to measure the spatial cross-correlations between pulsar pairs produced by a gravitational wave background. △ Less

Submitted 7 September, 2023; v1 submitted 20 May, 2023; originally announced May 2023.

Comments: Published

Journal ref: Phys. Rev. D 108, 063008 (2023)

Efficient large-scale, targeted gravitational-wave probes of supermassive black-hole binaries

Authors: Maria Charisi, Stephen R. Taylor, Caitlin A. Witt, Jessie Runnoe

Abstract: Supermassive black hole binaries are promising sources of low-frequency gravitational waves (GWs) and bright electromagnetic emission. Pulsar timing array searches for resolved binaries are complex and computationally expensive and so far limited to only a few sources. We present an efficient approximation that empowers large-scale targeted multi-messenger searches by neglecting GW signal componen… ▽ More Supermassive black hole binaries are promising sources of low-frequency gravitational waves (GWs) and bright electromagnetic emission. Pulsar timing array searches for resolved binaries are complex and computationally expensive and so far limited to only a few sources. We present an efficient approximation that empowers large-scale targeted multi-messenger searches by neglecting GW signal components from the pulsar term. This Earth-term approximation provides similar constraints on the total mass and GW frequency of the binary, yet is $>100$ times more efficient. △ Less

Submitted 7 April, 2023; originally announced April 2023.

Comments: Comments welcome

The Need For Speed: Rapid Refitting Techniques for Bayesian Spectral Characterization of the Gravitational Wave Background Using PTAs

Authors: William G. Lamb, Stephen R. Taylor, Rutger van Haasteren

Abstract: Pulsar timing arrays (PTAs) have recently found evidence for a nanohertz-frequency stochastic gravitational-wave background (SGWB). Constraining its spectral characteristics will reveal its origins. To achieve this, we must understand how data and modeling conditions in each pulsar influence the precision and accuracy of SGWB spectral recovery, typically requiring many Bayesian analyses on real da… ▽ More Pulsar timing arrays (PTAs) have recently found evidence for a nanohertz-frequency stochastic gravitational-wave background (SGWB). Constraining its spectral characteristics will reveal its origins. To achieve this, we must understand how data and modeling conditions in each pulsar influence the precision and accuracy of SGWB spectral recovery, typically requiring many Bayesian analyses on real data sets and large-scale simulations that are slow and computationally taxing. To combat this, we have developed several new rapid approaches that operate on intermediate SGWB analysis products. These techniques refit SGWB spectral models to previously computed Bayesian posterior estimates of the timing power spectra. We test our new techniques on simulated PTA data sets and the NANOGrav 12.5-year data set, where in the latter our refit posterior achieves a Hellinger distance -- bounded between 0 for identical distributions and 1 for zero overlap -- from the current full production-level pipeline that is < 0.1. Our techniques are ~ $10^2$--$10^4$ times faster than the production-level likelihood and scale much more favorably (sub-linearly) as a PTA is expanded with new pulsars or observations. Our techniques also allow us to demonstrate conclusively that SGWB spectral characterization in PTA data sets is driven by the longest-timed pulsars and the best-measured power spectral densities. Indeed, the common-process spectral properties found in the NANOGrav 12.5-year data set are given by analyzing only the ~14 longest-timed pulsars out of the full 45 pulsar array, and we find that the 'shallowing' of the common-process power-law model occurs when gravitational-wave frequencies higher than ~50 nanohertz are included. The implementation of our techniques is openly available as a software suite to allow fast and flexible PTA SGWB spectral characterization and model selection. △ Less

Submitted 9 October, 2023; v1 submitted 27 March, 2023; originally announced March 2023.

Comments: 22 pages, 13 figures. Accepted to Physical Review D

Journal ref: Phys. Rev. D 108, 103019 (2023)

Searching for continuous Gravitational Waves in the second data release of the International Pulsar Timing Array

Authors: M. Falxa, S. Babak, P. T. Baker, B. Bécsy, A. Chalumeau, S. Chen, Z. Chen, N. J. Cornish, L. Guillemot, J. S. Hazboun, C. M. F. Mingarelli, A. Parthasarathy, A. Petiteau, N. S. Pol, A. Sesana, S. B. Spolaor, S. R. Taylor, G. Theureau, M. Vallisneri, S. J. Vigeland, C. A. Witt, X. Zhu, J. Antoniadis, Z. Arzoumanian, M. Bailes , et al. (102 additional authors not shown)

Abstract: The International Pulsar Timing Array 2nd data release is the combination of datasets from worldwide collaborations. In this study, we search for continuous waves: gravitational wave signals produced by individual supermassive black hole binaries in the local universe. We consider binaries on circular orbits and neglect the evolution of orbital frequency over the observational span. We find no evi… ▽ More The International Pulsar Timing Array 2nd data release is the combination of datasets from worldwide collaborations. In this study, we search for continuous waves: gravitational wave signals produced by individual supermassive black hole binaries in the local universe. We consider binaries on circular orbits and neglect the evolution of orbital frequency over the observational span. We find no evidence for such signals and set sky averaged 95% upper limits on their amplitude h 95 . The most sensitive frequency is 10nHz with h 95 = 9.1 10-15 . We achieved the best upper limit to date at low and high frequencies of the PTA band thanks to improved effective cadence of observations. In our analysis, we have taken into account the recently discovered common red noise process, which has an impact at low frequencies. We also find that the peculiar noise features present in some pulsars data must be taken into account to reduce the false alarm. We show that using custom noise models is essential in searching for continuous gravitational wave signals and setting the upper limit. △ Less

Submitted 19 March, 2023; originally announced March 2023.

Generalized optimal statistic for characterizing multiple correlated signals in pulsar timing arrays

Authors: Shashwat C. Sardesai, Sarah J. Vigeland, Kyle A. Gersbach, Stephen R. Taylor

Abstract: The optimal statistic (OS) is a frequentist estimator for the amplitude and significance of a spatially-correlated signal in pulsar timing array (PTA) data, and it is widely used to search for the gravitational wave background (GWB). However, the OS cannot perfectly distinguish between different spatial correlations. In this paper, we introduce the multiple component optimal statistic (MCOS): a ge… ▽ More The optimal statistic (OS) is a frequentist estimator for the amplitude and significance of a spatially-correlated signal in pulsar timing array (PTA) data, and it is widely used to search for the gravitational wave background (GWB). However, the OS cannot perfectly distinguish between different spatial correlations. In this paper, we introduce the multiple component optimal statistic (MCOS): a generalization of the OS that allows for multiple correlations to be simultaneously fit to the data. We use simulated data to show that this method more accurately recovers injected spatially correlated signals, and in particular reduces the risk of a false detection of a signal with the wrong spatial correlation. We also demonstrate that this method can be used to recover multiple correlated signals. △ Less

Submitted 17 January, 2024; v1 submitted 16 March, 2023; originally announced March 2023.

Comments: 11 pages, 9 figures, 4 tables

Journal ref: Phys. Rev. D 108 (2023) 124081

The NANOGrav 12.5-year Data Set: Bayesian Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries

Authors: Zaven Arzoumanian, Paul T. Baker, Laura Blecha, Harsha Blumer, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Bence Bécsy, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Siyuan Chen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Brendan Drachler, Justin A. Ellis, E. C. Ferrara, William Fiore, Emmanuel Fonseca, Gabriel E. Freedman , et al. (53 additional authors not shown)

Abstract: Pulsar timing array collaborations, such as the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), are seeking to detect nanohertz gravitational waves emitted by supermassive black hole binaries formed in the aftermath of galaxy mergers. We have searched for continuous waves from individual circular supermassive black hole binaries using the NANOGrav's recent 12.5-year data s… ▽ More Pulsar timing array collaborations, such as the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), are seeking to detect nanohertz gravitational waves emitted by supermassive black hole binaries formed in the aftermath of galaxy mergers. We have searched for continuous waves from individual circular supermassive black hole binaries using the NANOGrav's recent 12.5-year data set. We created new methods to accurately model the uncertainties on pulsar distances in our analysis, and we implemented new techniques to account for a common red noise process in pulsar timing array data sets while searching for deterministic gravitational wave signals, including continuous waves. As we found no evidence for continuous waves in our data, we placed 95\% upper limits on the strain amplitude of continuous waves emitted by these sources. At our most sensitive frequency of 7.65 nanohertz, we placed a sky-averaged limit of $h_0 < $ $(6.82 \pm 0.35) \times 10^{-15}$, and $h_0 <$ $(2.66 \pm 0.15) \times 10^{-15}$ in our most sensitive sky location. Finally, we placed a multi-messenger limit of $\mathcal{M} <$ $(1.41 \pm 0.02) \times 10^9 M_\odot$ on the chirp mass of the supermassive black hole binary candidate 3C~66B. △ Less

Submitted 6 June, 2023; v1 submitted 9 January, 2023; originally announced January 2023.

Comments: 22 pages, 12 figures. Accepted by ApJL

Quality over Quantity: Optimizing pulsar timing array analysis for stochastic and continuous gravitational wave signals

Authors: Lorenzo Speri, Nataliya K. Porayko, Mikel Falxa, Siyuan Chen, Jonathan R. Gair, Alberto Sesana, Stephen R. Taylor

Abstract: 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 presen… ▽ More 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 EPTA dataset, 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. △ Less

Submitted 6 November, 2022; originally announced November 2022.

An unusual pulse shape change event in PSR J1713+0747 observed with the Green Bank Telescope and CHIME

Authors: Ross J. Jennings, James M. Cordes, Shami Chatterjee, Maura A. McLaughlin, Paul B. Demorest, Zaven Arzoumanian, Paul T. Baker, Harsha Blumer, Paul R. Brook, Tyler Cohen, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Timothy Dolch, Elizabeth C. Ferrara, Emmanuel Fonseca, Deborah C. Good, Jeffrey S. Hazboun, Megan L. Jones, David L. Kaplan, Michael T. Lam, T. Joseph W. Lazio, Duncan R. Lorimer, Jing Luo, Ryan S. Lynch , et al. (19 additional authors not shown)

Abstract: The millisecond pulsar J1713+0747 underwent a sudden and significant pulse shape change between April 16 and 17, 2021 (MJDs 59320 and 59321). Subsequently, the pulse shape gradually recovered over the course of several months. We report the results of continued multi-frequency radio observations of the pulsar made using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) and the 100-meter G… ▽ More The millisecond pulsar J1713+0747 underwent a sudden and significant pulse shape change between April 16 and 17, 2021 (MJDs 59320 and 59321). Subsequently, the pulse shape gradually recovered over the course of several months. We report the results of continued multi-frequency radio observations of the pulsar made using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) and the 100-meter Green Bank Telescope (GBT) in a three-year period encompassing the shape change event, between February 2020 and February 2023. As of February 2023, 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 TOA residuals display a strong non-monotonic 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 (DM) alone (which would produce a delay proportional to $ν^{-2}$). However, it does bear some resemblance to the two previous "chromatic timing events" observed in J1713+0747 (Demorest et al. 2013; Lam et al. 2016), as well as to a similar event observed in PSR J1643-1224 in 2015 (Shannon et al. 2016). △ Less

Submitted 31 January, 2024; v1 submitted 21 October, 2022; originally announced October 2022.

Comments: 18 pages, 8 figures. Submitted to ApJ. Data available at https://doi.org/10.5281/zenodo.7236459

Disentangling Multiple Stochastic Gravitational Wave Background Sources in PTA Datasets

Authors: Andrew R. Kaiser, Nihan S. Pol, Maura A. McLaughlin, Siyuan Chen, Jeffrey S. Hazboun, Luke Zoltan Kelley, Joseph Simon, Stephen R. Taylor, Sarah J. Vigeland, Caitlin A. Witt

Abstract: With strong evidence of a common-spectrum stochastic process in the most recent datasets from the NANOGrav Collaboration, the European Pulsar Timing Array (PTA), Parkes PTA, and the International PTA, it is crucial to assess the effects of the several astrophysical and cosmological sources that could contribute to the stochastic gravitational wave background (GWB). Using the same dataset creation… ▽ More With strong evidence of a common-spectrum stochastic process in the most recent datasets from the NANOGrav Collaboration, the European Pulsar Timing Array (PTA), Parkes PTA, and the International PTA, it is crucial to assess the effects of the several astrophysical and cosmological sources that could contribute to the stochastic gravitational wave background (GWB). Using the same dataset creation and injection techniques as in Pol et al. (2021), we assess the separability of multiple GWBs by creating single and multiple GWB source datasets. We search for these injected sources using Bayesian PTA analysis techniques to assess recovery and separability of multiple astrophysical and cosmological backgrounds. For a GWB due to supermassive black hole binaries and an underlying weaker background due to primordial gravitational waves with a GW energy density ratio of $Ω_{\mathrm{PGW}}/Ω_{\mathrm{SMBHB}} = 0.5$, the Bayes' factor for a second process exceeds unity at 17 years, and increases with additional data. At 20 years of data, we are able to constrain the spectral index and amplitude of the weaker GWB at this density ratio to a fractional uncertainty of 64% and 110%, respectively, using current PTA methods and techniques. Using these methods and findings, we outline a basic protocol to search for multiple backgrounds in future PTA datasets. △ Less

Submitted 5 August, 2022; v1 submitted 3 August, 2022; originally announced August 2022.

Comments: 16 pages, 10 figures, 3 tables, accepted in ApJ

Forecasting pulsar timing array sensitivity to anisotropy in the stochastic gravitational wave background

Authors: Nihan Pol, Stephen R. Taylor, Joseph D. Romano

Abstract: Statistical anisotropy in the nanohertz-frequency gravitational-wave background (GWB) is expected to be detected by pulsar timing arrays (PTAs) in the near future. By developing a frequentist statistical framework that intrinsically restricts the GWB power to be positive, we establish scaling relations for multipole-dependent anisotropy decision thresholds that are a function of the noise properti… ▽ More Statistical anisotropy in the nanohertz-frequency gravitational-wave background (GWB) is expected to be detected by pulsar timing arrays (PTAs) in the near future. By developing a frequentist statistical framework that intrinsically restricts the GWB power to be positive, we establish scaling relations for multipole-dependent anisotropy decision thresholds that are a function of the noise properties, timing baselines, and cadences of the pulsars in a PTA. We verify that $(i)$ a larger number of pulsars, and $(ii)$ factors that lead to lower uncertainty on the cross-correlation measurements between pulsars, lead to a higher overall GWB signal-to-noise ratio, and lower anisotropy decision thresholds with which to reject the null hypothesis of isotropy. Using conservative simulations of realistic NANOGrav datasets, we predict that an anisotropic GWB with angular power $C_{l=1} > 0.3\,C_{l=0}$ may be sufficient to produce tension with isotropy at the $p = 3\times10^{-3}$ ($\sim3σ$) level in near-future NANOGrav data with a $20$~yr baseline. We present ready-to-use scaling relationships that can map these thresholds to any number of pulsars, configuration of pulsar noise properties, and sky coverage. We discuss how PTAs can improve the detection prospects for anisotropy, as well as how our methods can be adapted for more versatile searches. △ Less

Submitted 20 June, 2022; originally announced June 2022.

Comments: Submitted to ApJ. Comments welcome

Deep learning and Bayesian inference of gravitational-wave populations: Hierarchical black-hole mergers

Authors: Matthew Mould, Davide Gerosa, Stephen R. Taylor

Abstract: The catalog of gravitational-wave events is growing, and so are our hopes of constraining the underlying astrophysics of stellar-mass black-hole mergers by inferring the distributions of, e.g., masses and spins. While conventional analyses parametrize this population with simple phenomenological models, we propose an emulation-based approach that can compare astrophysical simulations against gravi… ▽ More The catalog of gravitational-wave events is growing, and so are our hopes of constraining the underlying astrophysics of stellar-mass black-hole mergers by inferring the distributions of, e.g., masses and spins. While conventional analyses parametrize this population with simple phenomenological models, we propose an emulation-based approach that can compare astrophysical simulations against gravitational-wave data. We combine state-of-the-art deep-learning techniques with hierarchical Bayesian inference and exploit our approach to constrain the properties of repeated black-hole mergers from the gravitational-wave events in the most recent LIGO/Virgo catalog. Deep neural networks allow us to (i) construct a flexible single-channel population model that accurately emulates simple parametrized numerical simulations of hierarchical mergers, (ii) estimate selection effects, and (iii) recover the branching ratios of repeated-merger generations. Among our results, we find the following: The distribution of host-environment escape speeds favors values less than $100~\mathrm{km\,s^{-1}}$ but is relatively flat, with around $37\%$ of first-generation mergers retained in their host environments; first-generation black holes are born with a maximum mass that is compatible with current estimates from pair-instability supernovae; there is multimodal substructure in both the mass and spin distributions, which, in our model, can be explained by repeated mergers; and binaries with a higher-generation component make up at least $14\%$ of the underlying population. Though these results are inferred through emulation of a simplified model, the deep-learning pipeline we present is readily applicable to realistic astrophysical simulations △ Less

Submitted 3 November, 2022; v1 submitted 7 March, 2022; originally announced March 2022.

Comments: 29 pages, 17 figures, 5 tables. Accepted in PRD

Journal ref: Phys. Rev. D 106, 103013 (2022)

A Parallelized Bayesian Approach To Accelerated Gravitational-Wave Background Characterization

Authors: Stephen R. Taylor, Joseph Simon, Levi Schult, Nihan Pol, William G. Lamb

Abstract: The characterization of nanohertz-frequency gravitational waves (GWs) with pulsar-timing arrays requires a continual expansion of datasets and monitored pulsars. Whereas detection of the stochastic GW background is predicated on measuring a distinctive pattern of inter-pulsar correlations, characterizing the background's spectrum is driven by information encoded in the power spectra of the individ… ▽ More The characterization of nanohertz-frequency gravitational waves (GWs) with pulsar-timing arrays requires a continual expansion of datasets and monitored pulsars. Whereas detection of the stochastic GW background is predicated on measuring a distinctive pattern of inter-pulsar correlations, characterizing the background's spectrum is driven by information encoded in the power spectra of the individual pulsars' time series. We propose a new technique for rapid Bayesian characterization of the stochastic GW background that is fully parallelized over pulsar datasets. This Factorized Likelihood (FL) technique empowers a modular approach to parameter estimation of the GW background, multi-stage model selection of a spectrally-common stochastic process and quadrupolar inter-pulsar correlations, and statistical cross-validation of measured signals between independent pulsar sub-arrays. We demonstrate the equivalence of this technique's efficacy with the full pulsar-timing array likelihood, yet at a fraction of the required time. Our technique is fast, easily implemented, and trivially allows for new data and pulsars to be combined with legacy datasets without re-analysis of the latter. △ Less

Submitted 30 March, 2022; v1 submitted 16 February, 2022; originally announced February 2022.

Comments: 14 pages, 6 figures. Matches version accepted by PRD

Gravitational Wave Statistics for Pulsar Timing Arrays: Examining Bias from Using a Finite Number of Pulsars

Authors: Aaron D. Johnson, Sarah J. Vigeland, Xavier Siemens, Stephen R. Taylor

Abstract: Recently, many different pulsar timing array (PTA) collaborations have reported strong evidence for a common stochastic process in their data sets. The reported amplitudes are in tension with previously computed upper limits. In this paper, we investigate how using a subset of a set of pulsars biases Bayesian upper limit recovery. We generate 500 simulated PTA data sets based on the NANOGrav 11-ye… ▽ More Recently, many different pulsar timing array (PTA) collaborations have reported strong evidence for a common stochastic process in their data sets. The reported amplitudes are in tension with previously computed upper limits. In this paper, we investigate how using a subset of a set of pulsars biases Bayesian upper limit recovery. We generate 500 simulated PTA data sets based on the NANOGrav 11-year data set with an injected stochastic gravitational wave background (GWB). We then compute upper limits by sampling individual pulsar likelihoods, and combine them through a factorized version of the PTA likelihood to obtain upper limits on the GWB amplitude using different numbers of pulsars. We find that it is possible to recover an upper limit (95\% credible interval) \textit{below} the injected value, and that it is significantly more likely for this to occur when using a subset of pulsars to compute the upper limit. When picking pulsars to induce the maximum possible bias, we find that the 95\% Bayesian upper limit recovered is below the injected value in 10.6\% (53 of 500) of realizations. Further, we find that if we choose a subset of pulsars in order to obtain a lower upper limit than when using the full set of pulsars, the distribution of upper limits obtained from these 500 realizations is shifted to lower amplitude values. △ Less

Submitted 27 June, 2022; v1 submitted 25 January, 2022; originally announced January 2022.

Comments: 10 pages, 6 figures; updated to journal version

The International Pulsar Timing Array second data release: Search for an isotropic Gravitational Wave Background

Authors: J. Antoniadis, Z. Arzoumanian, S. Babak, M. Bailes, A. -S. Bak Nielsen, P. T. Baker, C. G. Bassa, B. Becsy, A. Berthereau, M. Bonetti, A. Brazier, P. R. Brook, M. Burgay, S. Burke-Spolaor, R. N. Caballero, J. A. Casey-Clyde, A. Chalumeau, D. J. Champion, M. Charisi, S. Chatterjee, S. Chen, I. Cognard, J. M. Cordes, N. J. Cornish, F. Crawford , et al. (101 additional authors not shown)

Abstract: We searched for an isotropic stochastic gravitational wave background in the second data release of the International Pulsar Timing Array, a global collaboration synthesizing decadal-length pulsar-timing campaigns in North America, Europe, and Australia. In our reference search for a power law strain spectrum of the form $h_c = A(f/1\,\mathrm{yr}^{-1})^α$, we found strong evidence for a spectrally… ▽ More We searched for an isotropic stochastic gravitational wave background in the second data release of the International Pulsar Timing Array, a global collaboration synthesizing decadal-length pulsar-timing campaigns in North America, Europe, and Australia. In our reference search for a power law strain spectrum of the form $h_c = A(f/1\,\mathrm{yr}^{-1})^α$, we found strong evidence for a spectrally-similar low-frequency stochastic process of amplitude $A = 3.8^{+6.3}_{-2.5}\times10^{-15}$ and spectral index $α= -0.5 \pm 0.5$, where the uncertainties represent 95\% credible regions, using information from the auto- and cross-correlation terms between the pulsars in the array. For a spectral index of $α= -2/3$, as expected from a population of inspiralling supermassive black hole binaries, the recovered amplitude is $A = 2.8^{+1.2}_{-0.8}\times10^{-15}$. Nonetheless, no significant evidence of the Hellings-Downs correlations that would indicate a gravitational-wave origin was found. We also analyzed the constituent data from the individual pulsar timing arrays in a consistent way, and clearly demonstrate that the combined international data set is more sensitive. Furthermore, we demonstrate that this combined data set produces comparable constraints to recent single-array data sets which have more data than the constituent parts of the combination. Future international data releases will deliver increased sensitivity to gravitational wave radiation, and significantly increase the detection probability. △ Less

Submitted 11 January, 2022; originally announced January 2022.

Comments: 17 pages, 12 figures, accepted in MNRAS

Controlling Outlier Contamination In Multimessenger Time-domain Searches For Supermasssive Binary Black Holes

Authors: Qiaohong Wang, Stephen R. Taylor

Abstract: Time-domain datasets of many varieties can be prone to statistical outliers that result from instrumental or astrophysical anomalies. These can impair searches for signals within the time series and lead to biased parameter estimation. Versatile outlier mitigation methods tuned toward multimessenger time-domain searches for supermassive binary black holes have yet to be fully explored. In an effor… ▽ More Time-domain datasets of many varieties can be prone to statistical outliers that result from instrumental or astrophysical anomalies. These can impair searches for signals within the time series and lead to biased parameter estimation. Versatile outlier mitigation methods tuned toward multimessenger time-domain searches for supermassive binary black holes have yet to be fully explored. In an effort to perform robust outlier isolation with low computational costs, we propose a Gibbs sampling scheme. This provides structural simplicity to outlier modeling and isolation, as it requires minimal modifications to adapt to time-domain modeling scenarios with pulsar-timing array or photometric data. We robustly diagnose outliers present in simulated pulsar-timing datasets, and then further apply our methods to pulsar J$1909$$-$$3744$ from the NANOGrav 9-yr Dataset. We also explore the periodic binary-AGN candidate PG$1302$$-$$102$ using datasets from the Catalina Real-time Transient Survey, All-Sky Automated Survey for Supernovae, and the Lincoln Near-Earth Asteroid Research. We present our findings and outline future work that could improve outlier modeling and isolation for multimessenger time-domain searches. △ Less

Submitted 16 September, 2022; v1 submitted 10 December, 2021; originally announced December 2021.

Comments: 14 pages, 10 figures. Matches version accepted in MNRAS

Multi-messenger Time-domain Signatures Of Supermassive Black Hole Binaries

Authors: Maria Charisi, Stephen R. Taylor, Jessie Runnoe, Tamara Bogdanovic, Jonathan R. Trump

Abstract: Supermassive black hole binaries (SMBHBs) are a natural outcome of galaxy mergers and should form frequently in galactic nuclei. Sub-parsec binaries can be identified from their bright electromagnetic emission, e.g., Active Galactic Nuclei (AGN) with Doppler shifted broad emission lines or AGN with periodic variability, as well as from the emission of strong gravitational radiation. The most massi… ▽ More Supermassive black hole binaries (SMBHBs) are a natural outcome of galaxy mergers and should form frequently in galactic nuclei. Sub-parsec binaries can be identified from their bright electromagnetic emission, e.g., Active Galactic Nuclei (AGN) with Doppler shifted broad emission lines or AGN with periodic variability, as well as from the emission of strong gravitational radiation. The most massive binaries (with total mass >10^8 M_sol) emit in the nanohertz band and are targeted by Pulsar Timing Arrays (PTAs). Here we examine the synergy between electromagnetic and gravitational wave signatures of SMBHBs. We connect both signals to the orbital dynamics of the binary and examine the common link between them, laying the foundation for joint multi-messenger observations. We find that periodic variability arising from relativistic Doppler boost is the most promising electromagnetic signature to connect with GWs. We delineate the parameter space (binary total mass/chirp mass versus binary period/GW frequency) for which joint observations are feasible. Currently multi-messenger detections are possible only for the most massive and nearby galaxies, limited by the sensitivity of PTAs. However, we demonstrate that as PTAs collect more data in the upcoming years, the overlapping parameter space is expected to expand significantly. △ Less

Submitted 27 October, 2021; originally announced October 2021.

Comments: submitted to MNRAS

Quasars with Periodic Variability: Capabilities and Limitations of Bayesian Searches for Supermassive Black Hole Binaries in Time-Domain Surveys

Authors: Caitlin A. Witt, Maria Charisi, Stephen R. Taylor, Sarah Burke-Spolaor

Abstract: Supermassive black hole binaries (SMBHBs) are an inevitable consequence of galaxy mergers. At subparsec separations, they are practically impossible to resolve, and the most promising technique is to search for quasars with periodic variability. However, searches for quasar periodicity in time-domain data are challenging due to the stochastic variability of quasars. In this paper, we used Bayesian… ▽ More Supermassive black hole binaries (SMBHBs) are an inevitable consequence of galaxy mergers. At subparsec separations, they are practically impossible to resolve, and the most promising technique is to search for quasars with periodic variability. However, searches for quasar periodicity in time-domain data are challenging due to the stochastic variability of quasars. In this paper, we used Bayesian methods to disentangle periodic SMBHB signals from intrinsic damped random walk (DRW) variability in active galactic nuclei light curves. We simulated a wide variety of realistic DRW and DRW+sine light curves. Their observed properties are modeled after the Catalina Real-time Transient Survey (CRTS) and expected properties of the upcoming Legacy Survey of Space and Time (LSST) from the Vera C. Rubin Observatory. Through a careful analysis of parameter estimation and Bayesian model selection, we investigated the range of parameter space for which binary systems can be detected. We also examined which DRW signals can mimic periodicity and be falsely classified as binary candidates. We found that periodic signals are more easily detectable if the period is short or the amplitude of the signal is large compared to the contribution of the DRW noise. We saw similar detection rates both in the CRTS and LSST-like simulations, while the false-detection rate depends on the quality of the data and is minimal in LSST. Our idealized simulations provide an excellent way to uncover the intrinsic limitations in quasar periodicity searches and set the stage for future searches for SMBHBs. △ Less

Submitted 17 August, 2022; v1 submitted 14 October, 2021; originally announced October 2021.

Comments: 18 pages, 9 figures. Accepted by ApJ

The NANOGrav 12.5-year data set: Search for Non-Einsteinian Polarization Modes in theGravitational-Wave Background

Authors: Zaven Arzoumanian, Paul T. Baker, Harsha Blumer, Bence Becsy, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Maria Charisi, Shami Chatterjee, Siyuan Chen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Dallas M. DeGan, Paul B. Demorest, Timothy Dolch, Brendan Drachler, Justin A. Ellis, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca, Nathan Garver-Daniels, Peter A. Gentile , et al. (46 additional authors not shown)

Abstract: We search NANOGrav's 12.5-year data set for evidence of a gravitational wave background (GWB) with all the spatial correlations allowed by general metric theories of gravity. We find no substantial evidence in favor of the existence of such correlations in our data. We find that scalar-transverse (ST) correlations yield signal-to-noise ratios and Bayes factors that are higher than quadrupolar (ten… ▽ More We search NANOGrav's 12.5-year data set for evidence of a gravitational wave background (GWB) with all the spatial correlations allowed by general metric theories of gravity. We find no substantial evidence in favor of the existence of such correlations in our data. We find that scalar-transverse (ST) correlations yield signal-to-noise ratios and Bayes factors that are higher than quadrupolar (tensor transverse, TT) correlations. Specifically, we find ST correlations with a signal-to-noise ratio of 2.8 that are preferred over TT correlations (Hellings and Downs correlations) with Bayesian odds of about 20:1. However, the significance of ST correlations is reduced dramatically when we include modeling of the Solar System ephemeris systematics and/or remove pulsar J0030$+$0451 entirely from consideration. Even taking the nominal signal-to-noise ratios at face value, analyses of simulated data sets show that such values are not extremely unlikely to be observed in cases where only the usual TT modes are present in the GWB. In the absence of a detection of any polarization mode of gravity, we place upper limits on their amplitudes for a spectral index of $γ= 5$ and a reference frequency of $f_\text{yr} = 1 \text{yr}^{-1}$. Among the upper limits for eight general families of metric theories of gravity, we find the values of $A^{95\%}_{TT} = (9.7 \pm 0.4)\times 10^{-16}$ and $A^{95\%}_{ST} = (1.4 \pm 0.03)\times 10^{-15}$ for the family of metric spacetime theories that contain both TT and ST modes. △ Less

Submitted 29 September, 2021; originally announced September 2021.

Comments: 24 pages, 18 figures, 3 appendices. Please send any comments/questions to Nima Laal ([email protected])

Journal ref: The Astrophysical Journal Letters, vol. 923, no. 2, p. L22, Dec. 2021

Discriminating between different scenarios for the formation and evolution of massive black holes with LISA

Authors: Alexandre Toubiana, Kaze W. K. Wong, Stanislav Babak, Enrico Barausse, Emanuele Berti, Jonathan R. Gair, Sylvain Marsat, Stephen R. Taylor

Abstract: Electromagnetic observations have provided strong evidence for the existence of massive black holes in the center of galaxies, but their origin is still poorly known. Different scenarios for the formation and evolution of massive black holes lead to different predictions for their properties and merger rates. LISA observations of coalescing massive black hole binaries could be used to reverse engi… ▽ More Electromagnetic observations have provided strong evidence for the existence of massive black holes in the center of galaxies, but their origin is still poorly known. Different scenarios for the formation and evolution of massive black holes lead to different predictions for their properties and merger rates. LISA observations of coalescing massive black hole binaries could be used to reverse engineer the problem and shed light on these mechanisms. In this paper, we introduce a pipeline based on hierarchical Bayesian inference to infer the mixing fraction between different theoretical models by comparing them to LISA observations of massive black hole mergers. By testing this pipeline against simulated LISA data, we show that it allows us to accurately infer the properties of the massive black hole population as long as our theoretical models provide a reliable description of the Universe. We also show that measurement errors, including both instrumental noise and weak lensing errors, have little impact on the inference. △ Less

Submitted 28 October, 2021; v1 submitted 25 June, 2021; originally announced June 2021.

Comments: Matches PRD version, minor changes

The Nanohertz Gravitational Wave Astronomer

Authors: Stephen R. Taylor

Abstract: Gravitational waves are a radically new way to peer into the darkest depths of the cosmos. Pulsars can be used to make direct detections of gravitational waves through precision timing. When a gravitational wave passes between a pulsar and the Earth, it stretches and squeezes the intermediate space-time, leading to deviations of the measured pulse arrival times away from model expectations. Combin… ▽ More Gravitational waves are a radically new way to peer into the darkest depths of the cosmos. Pulsars can be used to make direct detections of gravitational waves through precision timing. When a gravitational wave passes between a pulsar and the Earth, it stretches and squeezes the intermediate space-time, leading to deviations of the measured pulse arrival times away from model expectations. Combining the data from many Galactic pulsars can corroborate such a signal, and enhance its detection significance. This technique is known as a Pulsar Timing Array (PTA). Here I provide an overview of PTAs as a precision gravitational-wave detection instrument, then review the types of signal and noise processes that we encounter in typical pulsar data analysis. I take a pragmatic approach, illustrating how searches are performed in real life, and where possible directing the reader to codes or techniques that they can explore for themselves. The goal is to provide theoretical background and practical recipes for data exploration that allow the reader to join in the exciting hunt for very low frequency gravitational waves. △ Less

Submitted 27 May, 2021; originally announced May 2021.

Comments: Draft of a short technical book to be published later this year by Taylor & Francis. 156 pages. Comments and errata are welcome

Searching For Gravitational Waves From Cosmological Phase Transitions With The NANOGrav 12.5-year dataset

Authors: Zaven Arzoumanian, Paul T. Baker, Harsha Blumer, Bence Bécsy, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Maria Charisi, Shami Chatterjee, Siyuan Chen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Justin A. Ellis, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca, Nathan Garver-Daniels, Peter A. Gentile, Deborah C. Good, Jeffrey S. Hazboun , et al. (40 additional authors not shown)

Abstract: We search for a first-order phase transition gravitational wave signal in 45 pulsars from the NANOGrav 12.5 year dataset. We find that the data can be modeled in terms of a strong first order phase transition taking place at temperatures below the electroweak scale. However, we do not observe any strong preference for a phase-transition interpretation of the signal over the standard astrophysical… ▽ More We search for a first-order phase transition gravitational wave signal in 45 pulsars from the NANOGrav 12.5 year dataset. We find that the data can be modeled in terms of a strong first order phase transition taking place at temperatures below the electroweak scale. However, we do not observe any strong preference for a phase-transition interpretation of the signal over the standard astrophysical interpretation in terms of supermassive black holes mergers; but we expect to gain additional discriminating power with future datasets, improving the signal to noise ratio and extending the sensitivity window to lower frequencies. An interesting open question is how well gravitational wave observatories could separate such signals. △ Less

Submitted 11 January, 2022; v1 submitted 28 April, 2021; originally announced April 2021.

Comments: 13 pages, 4 figures. v2: updated to match published version

Journal ref: Phys.Rev.Lett. 127 (2021) 25, 251302

Bayesian Forecasts for Dark Matter Substructure Searches with Mock Pulsar Timing Data

Authors: Vincent S. H. Lee, Stephen R. Taylor, Tanner Trickle, Kathryn M. Zurek

Abstract: Dark matter substructure, such as primordial black holes (PBHs) and axion miniclusters, can induce phase shifts in pulsar timing arrays (PTAs) measurements due to gravitational effects. In order to gain a more realistic forecast for the detectability of such models of dark matter with PTAs, we propose a Bayesian inference framework to search for phase shifts generated by PBHs and perform the analy… ▽ More Dark matter substructure, such as primordial black holes (PBHs) and axion miniclusters, can induce phase shifts in pulsar timing arrays (PTAs) measurements due to gravitational effects. In order to gain a more realistic forecast for the detectability of such models of dark matter with PTAs, we propose a Bayesian inference framework to search for phase shifts generated by PBHs and perform the analysis on mock PTA data. For most PBH masses the constraints on the dark matter abundance agree with previous (frequentist) analyses (without mock data) to $\mathcal{O}(1)$ factors. This further motivates a dedicated search for PBHs (and dense small scale structures) in the mass range from $10^{-8}\,M_{\odot}$ to well above $10^2\,M_{\odot}$ with the Square Kilometer Array. Moreover, with a more optimistic set of timing parameters, future PTAs are predicted to constrain PBHs down to $10^{-11}\,M_{\odot}$. Lastly, we discuss the impact of backgrounds, such as Supermassive Black Hole Mergers, on detection prospects, suggesting a future program to separate a dark matter signal from other astrophysical sources. △ Less

Submitted 1 October, 2021; v1 submitted 12 April, 2021; originally announced April 2021.

Comments: 22 pages, 7 figures; v2: arguments revised, results and figures unchanged, matches journal version

Report number: CALT-TH-2021-016

Journal ref: Journal of Cosmology and Astroparticle Physics 08 (2021) 025

Mapping the Gravitational-wave Sky with LISA: A Bayesian Spherical Harmonic Approach

Authors: Sharan Banagiri, Alexander Criswell, Tommy Kuan, Vuk Mandic, Joseph D. Romano, Stephen R. Taylor

Abstract: The millihertz gravitational-wave frequency band is expected to contain a rich symphony of signals with sources ranging from galactic white dwarf binaries to extreme mass ratio inspirals. Many of these gravitational-wave signals will not be individually resolvable. Instead, they will incoherently add to produce stochastic gravitational-wave confusion noise whose frequency content will be governed… ▽ More The millihertz gravitational-wave frequency band is expected to contain a rich symphony of signals with sources ranging from galactic white dwarf binaries to extreme mass ratio inspirals. Many of these gravitational-wave signals will not be individually resolvable. Instead, they will incoherently add to produce stochastic gravitational-wave confusion noise whose frequency content will be governed by the dynamics of the sources. The angular structure of the power of the confusion noise will be modulated by the distribution of the sources across the sky. Measurement of this structure can yield important information about the distribution of sources on galactic and extra-galactic scales, their astrophysics and their evolution over cosmic timescales. Moreover, since the confusion noise is part of the noise budget of LISA, mapping it will also be essential for studying resolvable signals. In this paper, we present a Bayesian algorithm to probe the angular distribution of the stochastic gravitational-wave confusion noise with LISA using a spherical harmonic basis. We develop a technique based on Clebsch-Gordan coefficients to mathematically constrain the spherical harmonics to yield a non-negative distribution, making them optimal for expanding the gravitational-wave power and amenable to Bayesian inference. We demonstrate these techniques using a series of simulations and analyses, including recovery of simulated distributed and localized sources of gravitational-wave power. We also apply this method to map the gravitational-wave foreground from galactic white-dwarfs using a simplified model of the galactic white dwarf distribution. △ Less

Submitted 12 September, 2021; v1 submitted 1 March, 2021; originally announced March 2021.

Comments: 13 pages, 6 figures

Journal ref: Mon. Notices Royal Astron. Soc. 2021, 09, 0035-8711

The NANOGrav 11yr Data Set: Limits on Supermassive Black Hole Binaries in Galaxies within 500Mpc

Authors: Zaven Arzoumanian, Paul T. Baker, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Bence Becsy, Maria Charisi, Shami Chatterjee, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Rodney D. Elliott, Justin A. Ellis, Elizabeth C. Ferrara, Emmanuel Fonseca, Nathan Garver-Daniels, Peter A. Gentile, Deborah C. Good, Jeffrey S. Hazboun, Kristina Islo, Ross J. Jennings , et al. (32 additional authors not shown)

Abstract: Supermassive black hole binaries (SMBHBs) should form frequently in galactic nuclei as a result of galaxy mergers. At sub-parsec separations, binaries become strong sources of low-frequency gravitational waves (GWs), targeted by Pulsar Timing Arrays (PTAs). We used recent upper limits on continuous GWs from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11yr dataset to… ▽ More Supermassive black hole binaries (SMBHBs) should form frequently in galactic nuclei as a result of galaxy mergers. At sub-parsec separations, binaries become strong sources of low-frequency gravitational waves (GWs), targeted by Pulsar Timing Arrays (PTAs). We used recent upper limits on continuous GWs from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11yr dataset to place constraints on putative SMBHBs in nearby massive galaxies. We compiled a comprehensive catalog of ~44,000 galaxies in the local universe (up to redshift ~0.05) and populated them with hypothetical binaries, assuming that the total mass of the binary is equal to the SMBH mass derived from global scaling relations. Assuming circular equal-mass binaries emitting at NANOGrav's most sensitive frequency of 8nHz, we found that 216 galaxies are within NANOGrav's sensitivity volume. We ranked the potential SMBHBs based on GW detectability by calculating the total signal-to-noise ratio (S/N) such binaries would induce within the NANOGrav array. We placed constraints on the chirp mass and mass ratio of the 216 hypothetical binaries. For 19 galaxies, only very unequal-mass binaries are allowed, with the mass of the secondary less than 10 percent that of the primary, roughly comparable to constraints on a SMBHB in the Milky Way. Additionally, we were able to exclude binaries delivered by major mergers (mass ratio of at least 1/4) for several of these galaxies. We also derived the first limit on the density of binaries delivered by major mergers purely based on GW data. △ Less

Submitted 7 January, 2021; originally announced January 2021.

Comments: Submitted to ApJ. Send comments to Maria Charisi ([email protected])

Astrophysics Milestones For Pulsar Timing Array Gravitational Wave Detection

Authors: Nihan S. Pol, Stephen R. Taylor, Luke Zoltan Kelley, Sarah J. Vigeland, Joseph Simon, Siyuan Chen, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Shami Chatterjee, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca, Nathan Garver-Daniels, Deborah C. Good , et al. (27 additional authors not shown)

Abstract: The NANOGrav Collaboration reported strong Bayesian evidence for a common-spectrum stochastic process in its 12.5-yr pulsar timing array dataset, with median characteristic strain amplitude at periods of a year of $A_{\rm yr} = 1.92^{+0.75}_{-0.55} \times 10^{-15}$. However, evidence for the quadrupolar Hellings \& Downs interpulsar correlations, which are characteristic of gravitational wave sign… ▽ More The NANOGrav Collaboration reported strong Bayesian evidence for a common-spectrum stochastic process in its 12.5-yr pulsar timing array dataset, with median characteristic strain amplitude at periods of a year of $A_{\rm yr} = 1.92^{+0.75}_{-0.55} \times 10^{-15}$. However, evidence for the quadrupolar Hellings \& Downs interpulsar correlations, which are characteristic of gravitational wave signals, was not yet significant. We emulate and extend the NANOGrav dataset, injecting a wide range of stochastic gravitational wave background (GWB) signals that encompass a variety of amplitudes and spectral shapes, and quantify three key milestones: (I) Given the amplitude measured in the 12.5 yr analysis and assuming this signal is a GWB, we expect to accumulate robust evidence of an interpulsar-correlated GWB signal with 15--17 yrs of data, i.e., an additional 2--5 yrs from the 12.5 yr dataset; (II) At the initial detection, we expect a fractional uncertainty of $40\%$ on the power-law strain spectrum slope, which is sufficient to distinguish a GWB of supermassive black-hole binary origin from some models predicting more exotic origins;(III) Similarly, the measured GWB amplitude will have an uncertainty of $44\%$ upon initial detection, allowing us to arbitrate between some population models of supermassive black-hole binaries. In addition, power-law models are distinguishable from those having low-frequency spectral turnovers once 20~yrs of data are reached. Even though our study is based on the NANOGrav data, we also derive relations that allow for a generalization to other pulsar-timing array datasets. Most notably, by combining the data of individual arrays into the International Pulsar Timing Array, all of these milestones can be reached significantly earlier. △ Less

Submitted 24 March, 2021; v1 submitted 22 October, 2020; originally announced October 2020.

Comments: 15 pages, 7 figures

The NANOGrav 12.5-year Data Set: Search For An Isotropic Stochastic Gravitational-Wave Background

Authors: Zaven Arzoumanian, Paul T. Baker, Harsha Blumer, Bence Becsy, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Shami Chatterjee, Siyuan Chen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Justin A. Ellis, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca, Nathan Garver-Daniels, Peter A. Gentile, Deborah C. Good, Jeffrey S. Hazboun, A. Miguel Holgado , et al. (36 additional authors not shown)

Abstract: We search for an isotropic stochastic gravitational-wave background (GWB) in the $12.5$-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. Our analysis finds strong evidence of a stochastic process, modeled as a power-law, with common amplitude and spectral slope across pulsars. The Bayesian posterior of the amplitude for an $f^{-2/3}$ power-… ▽ More We search for an isotropic stochastic gravitational-wave background (GWB) in the $12.5$-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. Our analysis finds strong evidence of a stochastic process, modeled as a power-law, with common amplitude and spectral slope across pulsars. The Bayesian posterior of the amplitude for an $f^{-2/3}$ power-law spectrum, expressed as the characteristic GW strain, has median $1.92 \times 10^{-15}$ and $5\%$--$95\%$ quantiles of $1.37$--$2.67 \times 10^{-15}$ at a reference frequency of $f_\mathrm{yr} = 1 ~\mathrm{yr}^{-1}$. The Bayes factor in favor of the common-spectrum process versus independent red-noise processes in each pulsar exceeds $10,000$. However, we find no statistically significant evidence that this process has quadrupolar spatial correlations, which we would consider necessary to claim a GWB detection consistent with general relativity. We find that the process has neither monopolar nor dipolar correlations, which may arise from, for example, reference clock or solar system ephemeris systematics, respectively. The amplitude posterior has significant support above previously reported upper limits; we explain this in terms of the Bayesian priors assumed for intrinsic pulsar red noise. We examine potential implications for the supermassive black hole binary population under the hypothesis that the signal is indeed astrophysical in nature. △ Less

Submitted 8 January, 2021; v1 submitted 9 September, 2020; originally announced September 2020.

Comments: 25 pages, 14 figures, 5 tables, 3 appendices. Published in The Astrophysical Journal Letters. Please send any comments/questions to Joseph Simon ([email protected]). Jupyter notebook tutorials and some MCMC chain files are available at https://github.com/nanograv/12p5yr_stochastic_analysis

Journal ref: The Astrophysical Journal Letters, Volume 905, Number 2 (2020)

Inferring the properties of a population of compact binaries in presence of selection effects

Authors: Salvatore Vitale, Davide Gerosa, Will M. Farr, Stephen R. Taylor

Abstract: Shortly after a new class of objects is discovered, the attention shifts from the properties of the individual sources to the question of their origin: do all sources come from the same underlying population, or several populations are required? What are the properties of these populations? As the detection of gravitational waves is becoming routine and the size of the event catalog increases, fin… ▽ More Shortly after a new class of objects is discovered, the attention shifts from the properties of the individual sources to the question of their origin: do all sources come from the same underlying population, or several populations are required? What are the properties of these populations? As the detection of gravitational waves is becoming routine and the size of the event catalog increases, finer and finer details of the astrophysical distribution of compact binaries are now within our grasp. This Chapter presents a pedagogical introduction to the main statistical tool required for these analyses: hierarchical Bayesian inference in the presence of selection effects. All key equations are obtained from first principles, followed by two examples of increasing complexity. Although many remarks made in this Chapter refer to gravitational-wave astronomy, the write-up is generic enough to be useful to researchers and graduate students from other fields. △ Less

Submitted 18 March, 2022; v1 submitted 10 July, 2020; originally announced July 2020.

Comments: 57 pages. Updated and revised w.r.t. v1. Includes new section (5.2). v2. adds back the glossary, that was lost in previous version

Pulsar Timing Array Constraints on the Merger Timescale of Subparsec Supermassive Black Hole Binary Candidates

Authors: Khai Nguyen, Tamara Bogdanovic, Jessie C. Runnoe, Stephen R. Taylor, Alberto Sesana, Michael Eracleous, Steinn Sigurdsson

Abstract: We estimate the merger timescale of spectroscopically-selected, subparsec supermassive black hole binary (SMBHB) candidates by comparing their expected contribution to the gravitational wave background (GWB) with the sensitivity of current pulsar timing array (PTA) experiments and in particular, with the latest upper limit placed by the North American Nanohertz Observatory for Gravitational Waves… ▽ More We estimate the merger timescale of spectroscopically-selected, subparsec supermassive black hole binary (SMBHB) candidates by comparing their expected contribution to the gravitational wave background (GWB) with the sensitivity of current pulsar timing array (PTA) experiments and in particular, with the latest upper limit placed by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). We find that the average timescale to coalescence of such SMBHBs is $\langle t_{\rm evol} \rangle > 6\times 10^4\,$yr, assuming that their orbital evolution in the PTA frequency band is driven by emission of gravitational waves. If some fraction of SMBHBs do not reside in spectroscopically detected active galaxies, and their incidence in active and inactive galaxies is similar, then the merger timescale could be $\sim 10$ times longer, $\langle t_{\rm evol} \rangle > 6\times 10^5\,$yr. These limits are consistent with the range of timescales predicted by theoretical models and imply that all the SMBHB candidates in our spectroscopic sample could be binaries without violating the observational constraints on the GWB. This result illustrates the power of the multi-messenger approach, facilitated by the PTAs, in providing an independent statistical test of the nature of SMBHB candidates discovered in electromagnetic searches. △ Less

Submitted 9 October, 2020; v1 submitted 22 June, 2020; originally announced June 2020.

Comments: Published in ApJL (9 pages, 3 figures)

Journal ref: ApJL, 900, L42 (2020)

From Bright Binaries To Bumpy Backgrounds: Mapping Realistic Gravitational Wave Skies With Pulsar-Timing Arrays

Authors: Stephen R. Taylor, Rutger van Haasteren, Alberto Sesana

Abstract: Within the next several years, pulsar-timing array programs will likely usher in the next era of gravitational-wave astronomy through the detection of a stochastic background of nanohertz-frequency gravitational waves, originating from a cosmological population of inspiraling supermassive binary black holes. While the source positions will likely be isotropic to a good approximation, the gravitati… ▽ More Within the next several years, pulsar-timing array programs will likely usher in the next era of gravitational-wave astronomy through the detection of a stochastic background of nanohertz-frequency gravitational waves, originating from a cosmological population of inspiraling supermassive binary black holes. While the source positions will likely be isotropic to a good approximation, the gravitational-wave angular power distribution will be anisotropic, with the most massive and/or nearby binaries producing signals that may resound above the background. We study such a realistic angular power distribution, developing fast and accurate sky-mapping strategies to localize pixels and extended regions of excess power while simultaneously modeling the background signal from the less massive and more distant ensemble. We find that power anisotropy will be challenging to discriminate from isotropy for realistic gravitational-wave skies, requiring SNR $>10$ in order to favor anisotropy with $10:1$ posterior odds in our case study. Amongst our techniques, modeling the population signal with multiple point sources in addition to an isotropic background provides the most physically-motivated and easily interpreted maps, while spherical-harmonic modeling of the square-root power distribution, $P(\hatΩ)^{1/2}$, performs best in discriminating from overall isotropy. Our techniques are modular and easily incorporated into existing pulsar-timing array analysis pipelines. △ Less

Submitted 26 August, 2020; v1 submitted 8 June, 2020; originally announced June 2020.

Comments: 19 pages, 7 figures. Code available at https://github.com/stevertaylor/gworf. Matches version accepted by PRD

Journal ref: Phys. Rev. D 102, 084039 (2020)

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

Authors: Zaven Arzoumanian, Paul T. Baker, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Bence Becsy, Maria Charisi, Shami Chatterjee, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Rodney D. Elliott, Justin A. Ellis, Robert D. Ferdman, Elizabeth C. Ferrara, Emmanuel Fonseca, Nathan Garver-Daniels, Peter A. Gentile, Deborah C. Good, Jeffrey S. Hazboun , et al. (34 additional authors not shown)

Abstract: 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 th… ▽ More 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. △ Less

Submitted 12 August, 2020; v1 submitted 14 May, 2020; originally announced May 2020.

Comments: 14 pages, 6 figures. Accepted by ApJ

The NANOGrav 12.5-year Data Set: Wideband Timing of 47 Millisecond Pulsars

Authors: Md F. Alam, Zaven Arzoumanian, Paul T. Baker, Harsha Blumer, Keith E. Bohler, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Keeisi Caballero, Richard S. Camuccio, Rachel L. Chamberlain, Shami Chatterjee, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Justin A. Ellis, Robert D. Ferdman, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca, Yhamil Garcia , et al. (45 additional authors not shown)

Abstract: We present a new analysis of the profile data from the 47 millisecond pulsars comprising the 12.5-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), which is presented in a parallel paper (Alam et al. 2021a; NG12.5). Our reprocessing is performed using "wideband" timing methods, which use frequency-dependent template profiles, simultaneous time-of-arrival… ▽ More We present a new analysis of the profile data from the 47 millisecond pulsars comprising the 12.5-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), which is presented in a parallel paper (Alam et al. 2021a; NG12.5). Our reprocessing is performed using "wideband" timing methods, which use frequency-dependent template profiles, simultaneous time-of-arrival (TOA) and dispersion measure (DM) measurements from broadband observations, and novel analysis techniques. In particular, the wideband DM measurements are used to constrain the DM portion of the timing model. We compare the ensemble timing results to NG12.5 by examining the timing residuals, timing models, and noise model components. There is a remarkable level of agreement across all metrics considered. Our best-timed pulsars produce encouragingly similar results to those from NG12.5. In certain cases, such as high-DM pulsars with profile broadening, or sources that are weak and scintillating, wideband timing techniques prove to be beneficial, leading to more precise timing model parameters by 10-15%. The high-precision, multi-band measurements of several pulsars indicate frequency-dependent DMs. Compared to the narrowband analysis in NG12.5, the TOA volume is reduced by a factor of 33, which may ultimately facilitate computational speed-ups for complex pulsar timing array analyses. This first wideband pulsar timing data set is a stepping stone, and its consistent results with NG12.5 assure us that such data sets are appropriate for gravitational wave analyses. △ Less

Submitted 18 December, 2020; v1 submitted 13 May, 2020; originally announced May 2020.

Comments: 62 pages, 55 figures, 5 tables, 3 appendices. Data available at http://nanograv.org/data/ and via DOI 10.5281/zenodo.4312887

Journal ref: The Astrophysical Journal Supplement Series, 252, 5 (2021)