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

A Foundation for Gravitational-Wave Population Inference within the LISA Global Fit

Authors: Alexander W. Criswell, Sharan Banagiri, Vera Delfavero, Maria Jose Bustamante-Rosell, Stephen R. Taylor, Robert Rosati

Abstract: Population inference in gravitational-wave astronomy allows us to connect individual detections to the astrophysics of compact objects and their environments. Current approaches employed for population inference with LIGO-Virgo-KAGRA data approximate evaluation of the hierarchical population likelihood via post-processing of individual-event posteriors. However, the case of the Laser Interferomete… ▽ More Population inference in gravitational-wave astronomy allows us to connect individual detections to the astrophysics of compact objects and their environments. Current approaches employed for population inference with LIGO-Virgo-KAGRA data approximate evaluation of the hierarchical population likelihood via post-processing of individual-event posteriors. However, the case of the Laser Interferometer Space Antenna (LISA) will be more complex for two main reasons: the transdimensional "global fit" approach to LISA data analysis which models all signals and noise simultaneously, and the presence of both individually-resolved signals and the unresolved stochastic ``Galactic foreground" arising from the Galactic binary population, which induces a circular dependence between the resolved and unresolved systems and our ability to detect the former. These challenges are not without opportunity; LISA's data will contain every mHz compact binary in the Milky Way -- either individually or within the Galactic foreground -- with great potential for Galactic and stellar astrophysics. We therefore propose an alternative approach: direct evaluation of the full hierarchical population likelihood within the LISA global fit. We develop a statistical formalism for joint inference of individually-resolved gravitational-wave sources, an unresolved stochastic foreground, and a shared, underlying astrophysical population, present PELARGIR, a prototype GPU-accelerated population inference module for the LISA global fit, demonstrate the formalism and PELARGIR via a toy model analysis, and lay out a roadmap towards an astrophysically-motivated LISA global fit with embedded population inference. While we apply the formalism here to the population of LISA Galactic binaries, it is applicable across the gravitational-wave spectrum with use cases in pulsar timing and next-generation terrestrial observatories. △ Less

Submitted 3 April, 2026; originally announced April 2026.

Comments: 26 pages, 7 figures; to be submitted to ApJ

Gravitational Wave Measurement of the Mbh-Mbulge Intrinsic Scatter at High Redshift

Authors: Cayenne Matt, Kayhan Gültekin, Gabriella Agazie, Nikita Agarwal, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Jeremy G. 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 , et al. (87 additional authors not shown)

Abstract: The observed GWB spectrum is higher in amplitude than model predictions by a factor of 2-3. Using a semi-analytic model, we evaluate the effect of a high-scatter supermassive black hole (SMBH) scaling relation (Mbh-Mbulge) on models of the nanoHertz gravitational wave background (GWB). By implementing an intrinsic scatter of the Mbh-Mbulge relation, which is larger at higher redshift, but matches… ▽ More The observed GWB spectrum is higher in amplitude than model predictions by a factor of 2-3. Using a semi-analytic model, we evaluate the effect of a high-scatter supermassive black hole (SMBH) scaling relation (Mbh-Mbulge) on models of the nanoHertz gravitational wave background (GWB). By implementing an intrinsic scatter of the Mbh-Mbulge relation, which is larger at higher redshift, but matches local observations, we find that the amplitude of GWB models increases to be consistent with the low-frequency end of the GWB spectrum. This amplitude increase is not uniform across frequencies, a strongly evolving scatter preferentially increases the number density of the most massive SMBHs which, in the GWB spectrum, minimizes the strength of the low-frequency turnover. Our models with positively evolving intrinsic scatter can reproduce the electromagnetically observed overmassive SMBHs at 4 < z < 6 without changing the Mbh-Mbulge normalization though we find that including moderate normalization evolution improves fits to the GWB data. We conclude that the Mbh-Mbulge relation which best describes the available GWB and electromagnetic data sets has intrinsic scatter which evolves as epsilon(z) = epsilon_0 + (0.56 +/- 0.4) log10(1 + z) and normalization which evolves as alpha(z) = alpha_0 (1 + z)^(0.84 +/- 0.35). The results of this work imply that the Mbh-Mbulge relation we see today is not universal throughout cosmic time and that a diversity of seeding models and growth mechanisms may be at play in the early stages of SMBH-galaxy evolution. △ Less

Submitted 11 March, 2026; originally announced March 2026.

Comments: Submitted to AAS journals

On the angular localization of gravitational-wave signals by pulsar timing arrays

Authors: Stephen R. Taylor

Abstract: We provide a complete study of the factors influencing gravitational-wave signal localization using pulsar timing arrays. We derive analytical expressions for the Cramér-Rao sky localization precision that delineate the impact of the angular proximity, $ξ$, between the pulsar and the gravitational wave source, and the precision, $σ_L$, with which pulsar distances are known. Interference between th… ▽ More We provide a complete study of the factors influencing gravitational-wave signal localization using pulsar timing arrays. We derive analytical expressions for the Cramér-Rao sky localization precision that delineate the impact of the angular proximity, $ξ$, between the pulsar and the gravitational wave source, and the precision, $σ_L$, with which pulsar distances are known. Interference between the Earth and pulsar terms creates rapid angular oscillations for sky-coordinate Fisher matrix elements that aids localization, which is complemented by more broadly varying antenna response gradient information. The relative importance of these factors depends on whether pulsar distances are known precisely [i.e., $σ_L\leqλ_\mathrm{GW}/(1-\cosξ)$] or imprecisely, respectively. If the former, tightening pulsar distance precisions improves signal localization according to $ΔΩ_\mathrm{sky}\proptoσ_L^2$ until the Earth-pulsar system reaches its diffraction limit. If the latter, localization precision is degraded, but more pulsars in close proximity to the source is the best means of improving. With $α$ indexing pulsars, this scales as $ΔΩ_\mathrm{sky}~\propto~(\sum_α\mathrm{SNR}_α^2/ξ_α^2)^{-1}$ in the small-angle limit of the unmarginalized Fisher matrix, and we derive the analytic generalization to any angle between a pulsar and the source. Finally, we study a scenario where pulsar-term phases are treated as nuisance variables that are unconnected to binary or PTA properties. This phase-decoupled scenario, which is how all PTA continuous wave searches are currently conducted, delivers localization performance similar to the antenna-response--driven case, and does not exhibit significant improvement as pulsar distance precisions are tightened. △ Less

Submitted 10 March, 2026; originally announced March 2026.

Comments: 14 pages, 9 figures

The NANOGrav 15 yr Data Set: Piecewise Power-Law Reconstruction of the Gravitational-Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Jeremy G. Baier, Paul T. Baker, Bence Bécsy, Amit Bhoonah, 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, Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest , et al. (86 additional authors not shown)

Abstract: The NANOGrav 15-year (NG15) data set provides evidence for a gravitational-wave background (GWB) signal at nanohertz frequencies, which is expected to originate either from a cosmic population of inspiraling supermassive black-hole binaries or new particle physics in the early Universe. A firm identification of the source of the NG15 signal requires an accurate reconstruction of its frequency spec… ▽ More The NANOGrav 15-year (NG15) data set provides evidence for a gravitational-wave background (GWB) signal at nanohertz frequencies, which is expected to originate either from a cosmic population of inspiraling supermassive black-hole binaries or new particle physics in the early Universe. A firm identification of the source of the NG15 signal requires an accurate reconstruction of its frequency spectrum. In this paper, we provide such a spectral characterization of the NG15 signal based on a piecewise power-law (PPL) ansatz that strikes a balance between existing alternatives in the literature. Our PPL reconstruction is more flexible than the standard constant-power-law model, which describes the GWB spectrum in terms of only two parameters: an amplitude A and a spectral index gamma. Concurrently, it better approximates physically realistic GWB spectra -- especially those of cosmological origin -- than the free spectral model, since the latter allows for arbitrary variations in the GWB amplitude from one frequency bin to the next. Our PPL reconstruction of the NG15 signal relies on individual PPL models with a fixed number of internal nodes (i.e., constant power law, broken power law, doubly broken power law, etc.) that are ultimately combined in a Bayesian model average. The data products resulting from our analysis provide the basis for fast refits of spectral GWB models. △ Less

Submitted 14 January, 2026; originally announced January 2026.

Comments: 21 pages, 7 figures, 1 table

FrankenStat I: a New Approach to Pulsar Timing Array Data Combination

Authors: David Wright, Kalista Wayt, Jeffrey S. Hazboun, Xavier Siemens, Rutger van Haasteren, Levi Schult, Stephen R. Taylor

Abstract: In 2023, after more than two decades of searching, pulsar timing array (PTA) collaborations around the world announced evidence for a stochastic gravitational wave background. It was quickly followed by work from the International Pulsar Timing Array (IPTA), demonstrating that the results of regional collaborations were consistent with each other. The combination of these datasets is still ongoing… ▽ More In 2023, after more than two decades of searching, pulsar timing array (PTA) collaborations around the world announced evidence for a stochastic gravitational wave background. It was quickly followed by work from the International Pulsar Timing Array (IPTA), demonstrating that the results of regional collaborations were consistent with each other. The combination of these datasets is still ongoing and represents a significant investment of time and expertise. In that IPTA comparison, authors of this letter combined the separate datasets in the standard PTA optimal detection statistic for cross-correlations incoherently, that is, the data was combined without fitting a merged timing model across all PTA datasets, treating datasets of the same pulsar as independent, and neglecting the "same pulsar, different datasets" cross-correlations. This work refines that method by extending its core ideas beyond detection statistics and into a full, general data-combination method. We have demonstrated its efficacy and extreme efficiency on simulated data. This new method, \textit{FrankenStat}, is very similar in sensitivity and parameter-constraining power to traditional data combination methods while completing the full data combination in just a few minutes. △ Less

Submitted 16 December, 2025; originally announced December 2025.

Comments: 8 pages, 4 figures

Multimessenger Probes of the Supermassive Black Hole Binary Population: The Role of Pulsar Timing Arrays

Authors: Nima Laal, Stephen R. Taylor, Cayenne Matt, Kayhan Gultekin

Abstract: By inferring the gravitational wave background (GWB) from a population of supermassive black hole binaries (SMBHBs), pulsar timing arrays (PTAs) enable the study of massive black holes. In many ways, PTAs manifest the promise of a multimessenger approach to astronomy: they can constrain SMBHB population characteristics that are otherwise difficult to constrain using electromagnetic observations, s… ▽ More By inferring the gravitational wave background (GWB) from a population of supermassive black hole binaries (SMBHBs), pulsar timing arrays (PTAs) enable the study of massive black holes. In many ways, PTAs manifest the promise of a multimessenger approach to astronomy: they can constrain SMBHB population characteristics that are otherwise difficult to constrain using electromagnetic observations, such as hardening mechanisms at sub-parsec separations. In this work, we quantify this multimessenger promise using Bayesian inference of many realizations of simulated PTA data, while adopting a model for the SMBHBs that has been successfully applied to the 15-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). Our analyses of 200 realistic, simulated NANOGrav data sets show that there is a greater than 50\% chance of reducing the prior uncertainty in the SMBHB hardening rate by more than 50\%, and in the SMBHB evolutionary lifetime by 25--75\%. Additionally, there is an 88\% chance that PTA data can reduce the prior uncertainty in the characteristic mass variable of the galaxy stellar mass function (GSMF) by 25--50\%. For $M_{\text{BH}}$--$M_{\text{Bulge}}$ parameters (in a model without redshift evolution) and the overall normalization parameter of the GSMF, PTA data can provide only marginal information gain beyond the constraints from electromagnetic observations. Our work delineates the domains over which electromagnetic and gravitational-wave data constrain the demographics and dynamics of the supermassive black-hole binary population, offering a clearer picture of the impact of population multi-messenger astrophysics probes with PTAs. △ Less

Submitted 12 December, 2025; originally announced December 2025.

Comments: Under review for publication

The NANOGrav 12.5-year Data Set: Chromatic Noise Characterization & Mitigation with Time-Domain Kernels

Authors: Jeffrey S. Hazboun, Joseph Simon, Jeremy Baier, Bjorn Larsen, Daniel J. Oliver, Paul T. Baker, Bence Bécsy, Siyuan Chen, Alberto Diaz Hernandez, Justin A. Ellis, A. Miguel Holgado, Kristina Islo, Aaron Johnson, Andrew R. Kaiser, Nima Laal, Alexander McEwen, Nihan S. Pol, Joey Shapiro Key, Min Young Kim, Matthew Samson, Brent J. Shapiro-Albert, Jerry P. Sun, Stephen R. Taylor, Caitlin A. Witt, Jeremy Volpe , et al. (42 additional authors not shown)

Abstract: Pulsar timing arrays (PTAs) have recently entered the detection era, quickly moving beyond the goal of simply improving sensitivity at the lowest frequencies for the sake of observing the stochastic gravitational wave background (GWB), and focusing on its accurate spectral characterization. While all PTA collaborations around the world use Fourier-domain Gaussian processes to model the GWB and int… ▽ More Pulsar timing arrays (PTAs) have recently entered the detection era, quickly moving beyond the goal of simply improving sensitivity at the lowest frequencies for the sake of observing the stochastic gravitational wave background (GWB), and focusing on its accurate spectral characterization. While all PTA collaborations around the world use Fourier-domain Gaussian processes to model the GWB and intrinsic long time-correlated (red) noise, techniques to model the time-correlated radio frequency-dependent (chromatic) processes have varied from collaboration to collaboration. Here we test a new class of models for PTA data, Gaussian processes based on time-domain kernels that model the statistics of the chromatic processes starting from the covariance matrix. As we will show, these models can be effectively equivalent to Fourier-domain models in mitigating chromatic noise. This work presents a method for Bayesian model selection across the various choices of kernel as well as deterministic chromatic models for non-stationary chromatic events and the solar wind. As PTAs turn towards high frequency (>1/yr) sensitivity, the size of the basis used to model these processes will need to increase, and these time-domain models present some computational efficiencies compared to Fourier-domain models. △ Less

Submitted 27 November, 2025; originally announced November 2025.

Comments: 30 pages, 12 Figures, 6 Tables

Finite Populations & Finite Time: The Non-Gaussianity of a Gravitational Wave Background

Authors: William G. Lamb, Jeremy M. Wachter, Andrea Mitridate, Shashwat C. Sardesai, Bence Bécsy, Emily L. Hagen, Stephen R. Taylor, Luke Zoltan Kelley

Abstract: Strong evidence for an isotropic, Gaussian gravitational wave background (GWB) has been found by multiple pulsar timing arrays (PTAs). The GWB is expected to be sourced by a finite population of supermassive black hole binaries (SMBHBs) emitting in the PTA sensitivity band, and astrophysical inference of PTA data sets suggests a GWB signal that is at the higher end of GWB spectral amplitude estima… ▽ More Strong evidence for an isotropic, Gaussian gravitational wave background (GWB) has been found by multiple pulsar timing arrays (PTAs). The GWB is expected to be sourced by a finite population of supermassive black hole binaries (SMBHBs) emitting in the PTA sensitivity band, and astrophysical inference of PTA data sets suggests a GWB signal that is at the higher end of GWB spectral amplitude estimates. However, current inference analyses make simplifying assumptions, such as modeling the GWB as Gaussian, assuming that all SMBHBs only emit at frequencies that are integer multiples of the total observing time, and ignoring the interference between the signals of different SMBHBs. In this paper, we build analytical and numerical models of an astrophysical GWB without the above approximations, and compare the statistical properties of its induced PTA signal to those of a signal produced by a Gaussian GWB. We show that finite population and windowing effects introduce non-Gaussianities in the PTA signal, which are currently unmodeled in PTA analyses. △ Less

Submitted 4 February, 2026; v1 submitted 12 November, 2025; originally announced November 2025.

Comments: 19 pages, 8 figures. Submitted to Phys. Rev. D. Version 2 - corrected an equation and added comparison between analytical and numerical results

The Dawn of Gravitational Wave Astronomy at Light-year Wavelengths: Insights from Pulsar Timing Arrays

Authors: Stephen R. Taylor

Abstract: Arrays of precisely-timed millisecond pulsars are used to search for gravitational waves with periods of months to decades. Gravitational waves affect the path of radio pulses propagating from a pulsar to Earth, causing the arrival times of those pulses to deviate from expectations based on the physical characteristics of the pulsar system. By correlating these timing residuals in a pulsar timing… ▽ More Arrays of precisely-timed millisecond pulsars are used to search for gravitational waves with periods of months to decades. Gravitational waves affect the path of radio pulses propagating from a pulsar to Earth, causing the arrival times of those pulses to deviate from expectations based on the physical characteristics of the pulsar system. By correlating these timing residuals in a pulsar timing array (PTA), one can search for a statistically isotropic background of gravitational waves by revealing evidence for a distinctive pattern predicted by General Relativity, known as the Hellings \& Downs curve. On June 29 2023, five regional PTA collaborations announced the first evidence for GWs at light-year wavelengths, predicated on support for this correlation pattern with statistical significances ranging from $\sim\!2-4σ$. The amplitude and shape of the recovered GW spectrum has also allowed many investigations of the expected source characteristics, ranging from a cosmic population of supermassive binary black holes to numerous processes in the early Universe. In the future, we expect to resolve signals from individual binary systems of supermassive black holes, and probe fundamental assumptions about the background, including its polarization, anisotropy, Gaussianity, and stationarity, all of which will aid efforts to discriminate its origin. In tandem with new facilities like DSA-2000 and the SKA, fueling further observations by regional PTAs and the International Pulsar Timing Array, PTAs have extraordinary potential to be engines of nanohertz GW discovery. △ Less

Submitted 11 November, 2025; originally announced November 2025.

Comments: Invited review for Springer Nature Astronomy Awardee Prize Collection. Matches version for publication

Expectations for the first supermassive black-hole binary resolved by PTAs I: Model efficacy

Authors: Levi Schult, Polina Petrov, Stephen R. Taylor, Nihan Pol, Nima Laal, Maria Charisi, Chung-Pei Ma

Abstract: One of the most promising targets for Pulsar Timing Arrays (PTAs) is identifying an individual supermassive black hole binary (SMBHB) out of the population of binaries theorized to produce a gravitational wave background (GWB). In this work, we emulate an evolving PTA dataset, complete with an increasing number of pulsars and timing baseline, into which we inject a single binary on top of a Gaussi… ▽ More One of the most promising targets for Pulsar Timing Arrays (PTAs) is identifying an individual supermassive black hole binary (SMBHB) out of the population of binaries theorized to produce a gravitational wave background (GWB). In this work, we emulate an evolving PTA dataset, complete with an increasing number of pulsars and timing baseline, into which we inject a single binary on top of a Gaussian GWB signal. We vary the binary's source parameters, including sky position and frequency, and create an ensemble of simulated datasets with which we assess current Bayesian binary search techniques. We apply two waveform-based template models and a frequency-resolved anisotropy search to these simulations to understand how they compare in their detection and characterization abilities. We find that a template-based search including the full gravitational-wave signal structure (i.e., both Earth and pulsar effects of an incident GW) returns the highest Bayes Factors (BF), exceeding our estimator's capabilities by (S/N)~9-19, and has the most robust parameter estimation. Our anisotropy model attains a realization-median BF>10 at 7<(S/N)<15. Interestingly, despite being a deterministic model, the Earth-term template struggles to detect and characterize low-frequency binaries (5 nHz). These binaries require higher (S/N)~16-19 to reach the same BF threshold. This is likely due to neglected confusion effects between the pulsar and Earth terms. By contrast, the frequency-resolved anisotropy model shows promise for parameter estimation despite treating a binary's GW signal as excess directional GW power without phase modeling. Sky location and frequency parameter constraints returned by the anisotropy model are only surpassed by the Earth term template model at (S/N)~12-13. Milestones for a first detection using the full-signal GW model are included in a companion paper (Petrov et al. 2025). △ Less

Submitted 7 October, 2025; v1 submitted 1 October, 2025; originally announced October 2025.

Comments: 25 pages, 6 figures, 3 tables. Submitted to Physical Review D

Expectations for the first supermassive black-hole binary resolved by PTAs II: Milestones for binary characterization

Authors: Polina Petrov, Levi Schult, Stephen R. Taylor, Nihan Pol, Nima Laal, Maria Charisi, Chung-Pei Ma

Abstract: Following the recent evidence for a gravitational wave (GW) background found by pulsar timing array (PTA) experiments, the next major science milestone is resolving individual supermassive black hole binaries (SMBHBs). The detection of these systems could arise via searches using a power-based GW anisotropy model or a deterministic template model. In Schult et al. 2025, we compared the efficacy of… ▽ More Following the recent evidence for a gravitational wave (GW) background found by pulsar timing array (PTA) experiments, the next major science milestone is resolving individual supermassive black hole binaries (SMBHBs). The detection of these systems could arise via searches using a power-based GW anisotropy model or a deterministic template model. In Schult et al. 2025, we compared the efficacy of these models in constraining the GW signal from a single SMBHB using realistic, near-future PTA datasets, and found that the full-signal deterministic continuous wave (CW) search may achieve detection and characterization first. Here, we continue our analyses using only the CW model given its better performance, focusing now on characterization milestones. We examine the order in which CW parameters are constrained as PTA data are accumulated and the signal-to-noise ratio (S/N) grows. We also study how these parameter constraints vary across sources of different sky locations and GW frequencies. We find that the GW frequency and strain are generally constrained at the same time (or S/N), closely followed by the sky location, and later the chirp mass (if the source is highly evolving) and inclination angle. At fixed S/N, sources at higher frequencies generally achieve better precision on the GW frequency, chirp mass, and sky location. The time (and S/N) at which the signal becomes constrained is dependent on the sky location and frequency of the source, with the effects of pulsar terms and PTA geometry playing crucial roles in source detection and localization. △ Less

Submitted 7 October, 2025; v1 submitted 1 October, 2025; originally announced October 2025.

Comments: Submitted to Physical Review D

Mapping the Gravitational-wave Background Across the Spectrum with a Next-Generation Anisotropic Per-frequency Optimal Statistic

Authors: Kyle A. Gersbach, Stephen R. Taylor, Bence Bécsy, Anna-Malin Lemke, Andrea Mitridate, Nihan Pol

Abstract: With pulsar timing arrays (PTAs) having observed a gravitational wave background (GWB) at nanohertz frequencies, the focus of the field is shifting towards determining and characterizing its origin. While the primary candidate is a population of GW-emitting supermassive black hole binaries (SMBHBs), many other cosmological processes could produce a GWB with similar spectral properties as have been… ▽ More With pulsar timing arrays (PTAs) having observed a gravitational wave background (GWB) at nanohertz frequencies, the focus of the field is shifting towards determining and characterizing its origin. While the primary candidate is a population of GW-emitting supermassive black hole binaries (SMBHBs), many other cosmological processes could produce a GWB with similar spectral properties as have been measured. One key argument to help differentiate an SMBHB GWB from a cosmologically sourced one is its level of anisotropy; a GWB sourced by a finite population will likely exhibit greater anisotropy than a cosmological GWB through finite source effects (``shot noise'') and potentially large-scale structure. Current PTA GWB anisotropy detection methods often use the frequentist PTA optimal statistic for its fast estimation of pulsar pair correlations and relatively low computational overhead compared to spatially-correlated Bayesian analyses. However, there are critical limitations with the status quo approach. In this paper, we improve this technique by incorporating three recent advancements: accounting for covariance between pulsar pairwise estimates of correlated GWB power; the per-frequency optimal statistic to dissect the GWB across the spectrum; and constructing null-hypothesis statistical distributions that include cosmic variance. By combining these methods, our new pipeline can localize GWB anisotropies to specific frequencies, through which anisotropy detection prospects -- while impacted by cosmic variance -- are shown to improve in our simulations from a $p$-value of $\sim0.2$ in a broadband search to $\sim0.01$ in the per-frequency search. Our methods are already incorporated in community-available code and ready to deploy on forthcoming PTA datasets. △ Less

Submitted 8 September, 2025; originally announced September 2025.

Flexible Spectral Separation of Multiple Isotropic and Anisotropic Stochastic Gravitational Wave Backgrounds in LISA

Authors: Alexander W. Criswell, Sharan Banagiri, Jessica Lawrence, Levi Schult, Steven Rieck, Stephen R. Taylor, Vuk Mandic

Abstract: The Laser Interferometer Space Antenna (LISA) will observe mHz gravitational waves from a wide variety of astrophysical sources. Of these, some will be characterizable as individual deterministic signals; the remainder will overlap to create astrophysical confusion noise. These sources of confusion noise are known as stochastic gravitational wave backgrounds (SGWBs). LISA data is expected to inclu… ▽ More The Laser Interferometer Space Antenna (LISA) will observe mHz gravitational waves from a wide variety of astrophysical sources. Of these, some will be characterizable as individual deterministic signals; the remainder will overlap to create astrophysical confusion noise. These sources of confusion noise are known as stochastic gravitational wave backgrounds (SGWBs). LISA data is expected to include several such astrophysical SGWBs, including the notable Galactic binary foreground, SGWBs from white dwarf binary populations in satellite galaxies of the Milky Way, and the SGWB from extragalactic stellar-origin binary black holes far from merger. To characterize these astrophysical signals and attempt to seek out possible underlying backgrounds of cosmological origin, it will be necessary to separate the contribution of each SGWB from that of the others. Crucially, several of these SGWBs are expected to be highly anisotropic on the sky, providing a powerful tool for spectral separation. To this end, we present BLIP 2.0: a flexible, GPU-accelerated framework for simulation and Bayesian analysis of arbitrary combinations of isotropic and anisotropic SGWBs. We leverage these capabilities to demonstrate for the first time spectral separation of the Galactic foreground, the Large Magellanic Cloud SGWB, and the SGWB from extragalactic stellar-origin binaries, and show a proof-of-concept for placing upper limits on the detection of an underlying isotropic cosmological SGWB in the presence of multiple astrophysical foregrounds. △ Less

Submitted 27 August, 2025; originally announced August 2025.

Comments: 21 pages, 5 figures; to be submitted to Phys. Rev. D

The NANOGrav 15 yr Data Set: Targeted Searches for Supermassive Black Hole Binaries

Authors: Nikita Agarwal, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Jeremy G. Baier, Paul T. Baker, Bence Becsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, J. Andrew Casey-Clyde, Yu-Ting Chang, Maria Charisi, Shami Chatterjee, Tyler Cohen, Paolo Coppi, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter , et al. (94 additional authors not shown)

Abstract: We present the first targeted searches for continuous gravitational waves (CWs) from 114 active galactic nuclei (AGN) that may host supermassive black hole binaries, using the NANOGrav 15 yr data set. By incorporating electromagnetic priors on sky location, distance, redshift, and CW frequency, our strain and chirp mass upper limits are typically improved by a factor of $\sim 2$ (median 2.2) relat… ▽ More We present the first targeted searches for continuous gravitational waves (CWs) from 114 active galactic nuclei (AGN) that may host supermassive black hole binaries, using the NANOGrav 15 yr data set. By incorporating electromagnetic priors on sky location, distance, redshift, and CW frequency, our strain and chirp mass upper limits are typically improved by a factor of $\sim 2$ (median 2.2) relative to all-sky limits at the same frequency. Bayesian comparisons against a model including only a Hellings-Downs correlated background disfavors a CW signal for all targets, with a mean Bayes factor of $0.73 \pm 0.32$. Two targets have Bayes factors slightly above unity, but coherence tests, random targeting experiments, and a conservative accounting of the 114-target trials factor all indicate that they are consistent with noise. We use these two candidates as worked examples to illustrate an end-to-end targeted CW search analysis and a suite of follow up tests that future promising candidates would need to pass. We find that the electromagnetic interpretations of both candidates are ambiguous, and we update the constraints on a putative binary in 3C 66B, ruling out part of its previously allowed parameter space. Ultimately, our results demonstrate the current sensitivity of targeted pulsar timing array searches for CWs and define a roadmap for future multimessenger CW detections. △ Less

Submitted 13 January, 2026; v1 submitted 22 August, 2025; originally announced August 2025.

Comments: Accepted for publication in ApJ Letters

Journal ref: The Astrophysical Journal Letters, Volume 998, Issue 1, id.L11, 28 pp. (2026)

The Host Galaxies of PTA Sources: Converting Supermassive BH Binary Parameters into EM Observables

Authors: Niccolo Veronesi, Maria Charisi, Stephen R Taylor, Jessie Runnoe, Daniel J D'Orazio

Abstract: Pulsar timing arrays (PTAs) are approaching the sensitivity required to resolve gravitational waves (GWs) from individual supermassive black hole (SMBH) binaries. However, the large uncertainty in source localization will make the identification of its host environment challenging. We show how to convert the posterior probability function of binary parameters inferred by GW analyses into distribut… ▽ More Pulsar timing arrays (PTAs) are approaching the sensitivity required to resolve gravitational waves (GWs) from individual supermassive black hole (SMBH) binaries. However, the large uncertainty in source localization will make the identification of its host environment challenging. We show how to convert the posterior probability function of binary parameters inferred by GW analyses into distributions of apparent magnitudes of the host galaxy. We do so for a scenario in which the host environment is a regular early-type galaxy, and one in which it is an active galactic nucleus. We estimate the reach of PTAs in the near and intermediate future, and estimate whether the binary hosts will be detectable in all-sky electromagnetic (EM) surveys. A PTA with a baseline of 20 yr and 116 pulsars, resembling the upcoming data release of the International Pulsar Timing Array, can detect binaries out to a luminosity distance of 2 Gpc (corresponding to a redshift of $z\sim0.36$), while a PTA with a baseline of 30 yr and 200 pulsars can reach out to distances slightly greater than 3 Gpc ($z\sim0.53$). We find that the host galaxies of all binaries detectable with a baseline of 20 yr are expected to be present in the Wide-field Infrared Survey Explorer and SuperCOSMOS surveys, if they lie outside the plane of the Milky Way. The Two Micron All Sky Survey becomes incomplete for hosts of binaries more massive than $10^{9.8}{\rm M}_\odot$ at a luminosity distance greater than 1 Gpc. The EM surveys become slightly more incomplete when PTAs with longer baselines and therefore improved sensitivities are considered. △ Less

Submitted 25 August, 2025; v1 submitted 16 May, 2025; originally announced May 2025.

Comments: 14 pages, 7 figures

Journal ref: The Astrophysical Journal, 990:46, 2025 September 1

Rapid Construction of Joint Pulsar Timing Array Datasets: The Lite Method

Authors: Bjorn Larsen, Chiara M. F. Mingarelli, Paul T. Baker, Jeffrey S. Hazboun, Siyuan Chen, Levi Schult, Stephen R. Taylor, Joseph Simon, John Antoniadis, Jeremy Baier, R. Nicolaos Caballero, Aurélien Chalumeau, Zu-Cheng Chen, Ismael Cognard, Debabrata Deb, Valentina Di Marco, Timothy Dolch, Innocent O. Eya, Elizabeth C. Ferrara, Kyle A. Gersbach, Deborah C. Good, Huanchen Hu, Agastya Kapur, Shubham Kala, Michael Kramer , et al. (19 additional authors not shown)

Abstract: The International Pulsar Timing Array (IPTA)'s second data release (IPTA DR2) combines decades of observations of 65 millisecond pulsars from 7 radio telescopes. IPTA datasets should be the most sensitive datasets to nanohertz gravitational waves (GWs), but take years to assemble, often excluding valuable recent data. To address this, we introduce the IPTA "Lite" analysis, where a Figure of Merit… ▽ More The International Pulsar Timing Array (IPTA)'s second data release (IPTA DR2) combines decades of observations of 65 millisecond pulsars from 7 radio telescopes. IPTA datasets should be the most sensitive datasets to nanohertz gravitational waves (GWs), but take years to assemble, often excluding valuable recent data. To address this, we introduce the IPTA "Lite" analysis, where a Figure of Merit is used to select an optimal PTA dataset to analyze for each pulsar, enabling immediate access to new data and preliminary results prior to full combination. We test the capabilities of the Lite analysis using IPTA DR2, finding that "DR2 Lite" can be used to detect the common red noise process with an amplitude of $A = 4.8^{+1.8}_{-1.8} \times 10^{-15}$ at $γ= 13/3$. This amplitude is slightly large in comparison to the combined analysis, and likely biased high as DR2 Lite is more sensitive to systematic errors from individual pulsars than the full dataset. Furthermore, although there is no strong evidence for Hellings-Downs correlations in IPTA DR2, we still find the full dataset is better at resolving Hellings-Downs correlations than DR2 Lite. Alongside the Lite analysis, we also find that analyzing a subset of pulsars from IPTA DR2, available at a hypothetical "early" stage of combination (EDR2), yields equally competitive results as the full dataset. Looking ahead, the Lite method will enable rapid synthesis of the latest PTA data, offering preliminary GW constraints before the superior full dataset combinations are available. △ Less

Submitted 11 September, 2025; v1 submitted 26 March, 2025; originally announced March 2025.

Comments: MNRAS accepted version

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 542, Issue 4, October 2025, Pages 3028-3048

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

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Jeremy G. Baier, Paul T. Baker, Bence Becsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, 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, Heling Deng, Lankeswar Dey , et al. (80 additional authors not shown)

Abstract: We present the results of a search for nonlinear gravitational wave memory in the NANOGrav 15-year data set. We find no significant evidence for memory signals in the dataset, with a maximum Bayes factor of 3.1 in favor of a model including memory. We therefore place upper limits on the strain of potential gravitational wave memory events as a function of sky location and observing epoch. We find… ▽ More We present the results of a search for nonlinear gravitational wave memory in the NANOGrav 15-year data set. We find no significant evidence for memory signals in the dataset, with a maximum Bayes factor of 3.1 in favor of a model including memory. We therefore place upper limits on the strain of potential gravitational wave memory events as a function of sky location and observing epoch. We find upper limits that are not always more constraining than previous NANOGrav results. We show that it is likely due to the increase in common red noise between the 12.5-year and 15-year NANOGrav datasets. △ Less

Submitted 28 February, 2025; v1 submitted 25 February, 2025; originally announced February 2025.

Towards robust gravitational wave detections from individual supermassive black hole binaries

Authors: Bence Bécsy, Neil J. Cornish, Polina Petrov, Xavier Siemens, Stephen R. Taylor, Sarah J. Vigeland, Caitlin A. Witt

Abstract: The recent discovery of the stochastic gravitational-wave background via pulsar timing arrays will likely be followed by the detection of individual black hole binaries that stand out above the background. However, to confidently claim the detection of an individual binary, we need not only more and better data, but also more sophisticated analysis techniques. In this paper, we develop two new app… ▽ More The recent discovery of the stochastic gravitational-wave background via pulsar timing arrays will likely be followed by the detection of individual black hole binaries that stand out above the background. However, to confidently claim the detection of an individual binary, we need not only more and better data, but also more sophisticated analysis techniques. In this paper, we develop two new approaches that can help us more robustly ascertain if a candidate found by a search algorithm is indeed an individual supermassive black hole binary. One of these is a coherence test that directly compares the full signal model to an incoherent version of that. The other is a model scrambling approach that builds null distributions of our detection statistic and compares that with the measured value to quantify our confidence in signal coherence. Both of these rely on finding the coherence between pulsars characteristic to gravitational waves from a binary system. We test these methods on simple simulated datasets and find that they work well in correctly identifying both true gravitational waves and false positives. However, as expected for such a flexible and simple signal model, confidently identifying signal coherence is significantly harder than simply finding a candidate in most scenarios. Our analyses also indicate that the confidence with which we can identify a true signal depends not only on the signal-to-noise ratio, but also on the number of contributing pulsars and the amount of frequency evolution shown by the signal. △ Less

Submitted 25 February, 2025; originally announced February 2025.

Comments: 21 pages, 7 figures, 2 tables, submitted to CQG

The NANOGrav 15 Yr Data Set: Removing Pulsars One by One from the Pulsar Timing Array

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Jeremy G. Baier, Paul T. Baker, Bence Becsy, 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, Heling Deng, Lankeswar Dey, Timothy Dolch , et al. (80 additional authors not shown)

Abstract: Evidence has emerged for a stochastic signal correlated among 67 pulsars within the 15-year pulsar-timing data set compiled by the NANOGrav collaboration. Similar signals have been found in data from the European, Indian, Parkes, and Chinese PTAs. This signal has been interpreted as indicative of the presence of a nanohertz stochastic gravitational wave background. To explore the internal consiste… ▽ More Evidence has emerged for a stochastic signal correlated among 67 pulsars within the 15-year pulsar-timing data set compiled by the NANOGrav collaboration. Similar signals have been found in data from the European, Indian, Parkes, and Chinese PTAs. This signal has been interpreted as indicative of the presence of a nanohertz stochastic gravitational wave background. To explore the internal consistency of this result we investigate how the recovered signal strength changes as we remove the pulsars one by one from the data set. We calculate the signal strength using the (noise-marginalized) optimal statistic, a frequentist metric designed to measure correlated excess power in the residuals of the arrival times of the radio pulses. We identify several features emerging from this analysis that were initially unexpected. The significance of these features, however, can only be assessed by comparing the real data to synthetic data sets. After conducting identical analyses on simulated data sets, we do not find anything inconsistent with the presence of a stochastic gravitational wave background in the NANOGrav 15-year data. The methodologies developed here can offer additional tools for application to future, more sensitive data sets. While this analysis provides an internal consistency check of the NANOGrav results, it does not eliminate the necessity for additional investigations that could identify potential systematics or uncover unmodeled physical phenomena in the data. △ Less

Submitted 23 May, 2025; v1 submitted 22 November, 2024; originally announced November 2024.

Comments: 21 pages, 11 figures, 2 tables

The NANOGrav 15 yr Data Set: Harmonic Analysis of the Pulsar Angular Correlations

Authors: Gabriella Agazie, Jeremy G. Baier, Paul T. Baker, Bence Becsy, Laura Blecha, Kimberly K. Boddy, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, 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, Lankeswar Dey, Timothy Dolch, Elizabeth C. Ferrara, William Fiore , et al. (64 additional authors not shown)

Abstract: Pulsar timing array observations have found evidence for an isotropic gravitational wave background with the Hellings-Downs angular correlations, expected from general relativity. This interpretation hinges on the measured shape of the angular correlations, which is predominately quadrupolar under general relativity. Here we explore a more flexible parameterization: we expand the angular correlati… ▽ More Pulsar timing array observations have found evidence for an isotropic gravitational wave background with the Hellings-Downs angular correlations, expected from general relativity. This interpretation hinges on the measured shape of the angular correlations, which is predominately quadrupolar under general relativity. Here we explore a more flexible parameterization: we expand the angular correlations into a sum of Legendre polynomials and use a Bayesian analysis to constrain their coefficients with the 15-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). When including Legendre polynomials with multipoles $\ell \geq 2$, we only find a significant signal in the quadrupole with an amplitude consistent with general relativity and non-zero at the $\sim 95\%$ confidence level and a Bayes factor of 200. When we include multipoles $\ell \leq 1$, the Bayes factor evidence for quadrupole correlations decreases by more than an order of magnitude due to evidence for a monopolar signal at approximately 4 nHz which has also been noted in previous analyses of the NANOGrav 15-year data. Further work needs to be done in order to better characterize the properties of this monopolar signal and its effect on the evidence for quadrupolar angular correlations. △ Less

Submitted 20 November, 2024; originally announced November 2024.

Comments: 15 pages, 6 figures

Deep Neural Emulation of the Supermassive Black-hole Binary Population

Authors: Nima Laal, Stephen R. Taylor, Luke Zoltan Kelley, Joseph Simon, Kayhan Gultekin, David Wright, Bence Becsy, J. Andrew Casey-Clyde, Siyuan Chen, Alexander Cingoranelli, Daniel J. D'Orazio, Emiko C. Gardiner, William G. Lamb, Cayenne Matt, Magdalena S. Siwek, Jeremy M. Wachter

Abstract: While supermassive black-hole (SMBH)-binaries are not the only viable source for the low-frequency gravitational wave background (GWB) signal evidenced by the most recent pulsar timing array (PTA) data sets, they are expected to be the most likely. Thus, connecting the measured PTA GWB spectrum and the underlying physics governing the demographics and dynamics of SMBH-binaries is extremely importa… ▽ More While supermassive black-hole (SMBH)-binaries are not the only viable source for the low-frequency gravitational wave background (GWB) signal evidenced by the most recent pulsar timing array (PTA) data sets, they are expected to be the most likely. Thus, connecting the measured PTA GWB spectrum and the underlying physics governing the demographics and dynamics of SMBH-binaries is extremely important. Previously, Gaussian processes (GPs) and dense neural networks have been used to make such a connection by being built as conditional emulators; their input is some selected evolution or environmental SMBH-binary parameters and their output is the emulated mean and standard deviation of the GWB strain ensemble distribution over many Universes. In this paper, we use a normalizing flow (NF) emulator that is trained on the entirety of the GWB strain ensemble distribution, rather than only mean and standard deviation. As a result, we can predict strain distributions that mirror underlying simulations very closely while also capturing frequency covariances in the strain distributions as well as statistical complexities such as tails, non-Gaussianities, and multimodalities that are otherwise not learnable by existing techniques. In particular, we feature various comparisons between the NF-based emulator and the GP approach used extensively in past efforts. Our analyses conclude that the NF-based emulator not only outperforms GPs in the ease and computational cost of training but also outperforms in the fidelity of the emulated GWB strain ensemble distributions. △ Less

Submitted 15 November, 2024; originally announced November 2024.

Inference on inner galaxy structure via gravitational waves from supermassive binaries

Authors: Yifan Chen, Matthias Daniel, Daniel J. D'Orazio, Xuanye Fan, Andrea Mitridate, Laura Sagunski, Xiao Xue, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Jeremy G. Baier, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish , et al. (85 additional authors not shown)

Abstract: The detection of a stochastic gravitational wave background by pulsar-timing arrays indicates the presence of a population of supermassive black hole binaries. Although the observed spectrum generally matches predictions for orbital evolution driven by gravitational-wave emission in circular orbits, there is a preference for a spectral turnover at the lowest observed frequencies, which may point t… ▽ More The detection of a stochastic gravitational wave background by pulsar-timing arrays indicates the presence of a population of supermassive black hole binaries. Although the observed spectrum generally matches predictions for orbital evolution driven by gravitational-wave emission in circular orbits, there is a preference for a spectral turnover at the lowest observed frequencies, which may point to substantial hardening during a transition from early environmental influences to later stages dominated by emission. In the vicinity of these binaries, the ejection of stars or dark matter particles through gravitational three-body slingshots efficiently extracts orbital energy, leading to a low-frequency turnover in the spectrum. Here we model how the gravitational-wave spectrum depends on the initial inner galactic profile before scouring by binary ejections while accounting for a range of initial binary eccentricities. By analysing the NANOGrav 15-year data, we find that a parsec-scale galactic-centre density of around $10^6 M_{\odot} \mathrm{pc}^{-3}$ is favoured across most of the parameter space, thus shedding light on the environmental effects that shape black hole evolution and the combined matter density near galaxy centres. △ Less

Submitted 4 February, 2026; v1 submitted 8 November, 2024; originally announced November 2024.

Comments: 19 pages, 6 figures, accepted version in Nature Astronomy

Solving the PTA Data Analysis Problem with a Global Gibbs Scheme

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

Abstract: The announcement in the summer of 2023 about the discovery of evidence for a gravitational wave background (GWB) using pulsar timing arrays (PTAs) has ignited both the PTA and the larger scientific community's interest in the experiment and the scientific implications of its findings. As a result, numerous scientific works have been published analyzing and further developing various aspects of the… ▽ More The announcement in the summer of 2023 about the discovery of evidence for a gravitational wave background (GWB) using pulsar timing arrays (PTAs) has ignited both the PTA and the larger scientific community's interest in the experiment and the scientific implications of its findings. As a result, numerous scientific works have been published analyzing and further developing various aspects of the experiment, from performing tests of gravity to improving the efficiency of the current data analysis techniques. In this regard, we contribute to the recent advancements in the field of PTAs by presenting the most general, agnostic, per-frequency Bayesian search for a low-frequency (red) noise process in these data. Our new method involves the use of a conjugate Jeffrey's-like multivariate prior which allows one to model all unique parameters of the global PTA-level red noise covariance matrix as a separate model parameter for which a marginalized posterior-probability distribution can be found using Gibbs sampling. Even though perfecting the implementation of the Gibbs sampling and mitigating the numerical stability challenges require further development, we show the power of this new method by analyzing realistic and theoretical PTA simulated data sets. We show how our technique is consistent with the more restricted standard techniques in recovering both the auto and the cross-spectrum of pulsars' low-frequency (red) noise. Furthermore, we highlight ways to approximately characterize a GWB (both its auto- and cross-spectrum) using Fourier coefficient estimates from single-pulsar and so-called CURN (common uncorrelated red noise) analyses via analytic draws from a specific Inverse-Wishart distribution. △ Less

Submitted 21 April, 2025; v1 submitted 15 October, 2024; originally announced October 2024.

Comments: Published at https://link.aps.org/doi/10.1103/PhysRevD.111.063067

The NANOGrav 15 yr Data Set: Running of the Spectral Index

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Jeremy George Baier, 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, Heling Deng, Lankeswar Dey, Timothy Dolch , et al. (80 additional authors not shown)

Abstract: The NANOGrav 15-year data provides compelling evidence for a stochastic gravitational-wave (GW) background at nanohertz frequencies. The simplest model-independent approach to characterizing the frequency spectrum of this signal consists in a simple power-law fit involving two parameters: an amplitude A and a spectral index γ. In this paper, we consider the next logical step beyond this minimal sp… ▽ More The NANOGrav 15-year data provides compelling evidence for a stochastic gravitational-wave (GW) background at nanohertz frequencies. The simplest model-independent approach to characterizing the frequency spectrum of this signal consists in a simple power-law fit involving two parameters: an amplitude A and a spectral index γ. In this paper, we consider the next logical step beyond this minimal spectral model, allowing for a running (i.e., logarithmic frequency dependence) of the spectral index, γ_run(f) = γ+ β\ln(f/f_ref). We fit this running-power-law (RPL) model to the NANOGrav 15-year data and perform a Bayesian model comparison with the minimal constant-power-law (CPL) model, which results in a 95% credible interval for the parameter βconsistent with no running, β\in [-0.80,2.96], and an inconclusive Bayes factor, B(RPL vs. CPL) = 0.69 +- 0.01. We thus conclude that, at present, the minimal CPL model still suffices to adequately describe the NANOGrav signal; however, future data sets may well lead to a measurement of nonzero β. Finally, we interpret the RPL model as a description of primordial GWs generated during cosmic inflation, which allows us to combine our results with upper limits from big-bang nucleosynthesis, the cosmic microwave background, and LIGO-Virgo-KAGRA. △ Less

Submitted 30 January, 2025; v1 submitted 19 August, 2024; originally announced August 2024.

Comments: 18 pages, 4 figures, 2 tables. v2: matches version published in ApJ Letters

Journal ref: Astrophys. J. Lett. 978 (2025) 2, L29

The NANOGrav 15 yr data set: Posterior predictive checks for gravitational-wave detection with pulsar timing arrays

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Jeremy George Baier, Paul T. Baker, Bence Bécsy, 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, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Heling Deng, Lankeswar Dey , et al. (77 additional authors not shown)

Abstract: Pulsar-timing-array experiments have reported evidence for a stochastic background of nanohertz gravitational waves consistent with the signal expected from a population of supermassive--black-hole binaries. Their analyses assume power-law spectra for intrinsic pulsar noise and for the background, as well as a Hellings--Downs cross-correlation pattern among the gravitational-wave--induced residual… ▽ More Pulsar-timing-array experiments have reported evidence for a stochastic background of nanohertz gravitational waves consistent with the signal expected from a population of supermassive--black-hole binaries. Their analyses assume power-law spectra for intrinsic pulsar noise and for the background, as well as a Hellings--Downs cross-correlation pattern among the gravitational-wave--induced residuals across pulsars. These assumptions may not be realized in actuality. We test them in the NANOGrav 15 yr data set using Bayesian posterior predictive checks. After fitting our fiducial model to real data, we generate a population of simulated data-set replications. We use the replications to assess whether the optimal-statistic significance, inter-pulsar correlations, and spectral coefficients are extreme. We recover Hellings--Downs correlations in simulated data sets at significance levels consistent with the correlations measured in the NANOGrav 15 yr data set. A similar test on spectral coefficients shows that their values in real data are not extreme compared to their distributions across replications. We also evaluate the evidence for the stochastic background using posterior-predictive versions of the frequentist optimal statistic and of Bayesian model comparison, and find comparable significance (3.2 $σ$ and 3 $σ$ respectively) to what was previously reported for the standard statistics. We conclude with novel visualizations of the reconstructed gravitational waveforms that enter the residuals for each pulsar. Our analysis strengthens confidence in the identification and characterization of the gravitational-wave background. △ Less

Submitted 13 March, 2025; v1 submitted 29 July, 2024; originally announced July 2024.

Comments: 20 pages, 14 Figures, accepted to PRD

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 due to finite source effects. 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 5 August, 2024; v1 submitted 8 July, 2024; originally announced July 2024.

Comments: 8 pages, 3 figures. Published in ApJL

Journal ref: ApJL 971 L10 (2024)

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 \textit{Illustris} cosmological simulations taken in tandem with post processing semi-analytic models of e… ▽ 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 \textit{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 $\sim 2 Z_\odot$ as opposed to the total population of galaxies which peaks at $\sim 0.6 Z_\odot$. 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 25 January, 2025; v1 submitted 1 July, 2024; originally announced July 2024.

Comments: Included in CQG Focus Issue on Pulsar Timing Arrays

Journal ref: Classical and Quantum Gravity, Volume 42, Issue 2, 2025

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 5 February, 2025; v1 submitted 17 June, 2024; originally announced June 2024.

Comments: 23 pages, 9 figures, Accepted by 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 22 November, 2024; v1 submitted 6 June, 2024; originally announced June 2024.

Comments: Published in ApJ

Journal ref: ApJ 976 129 (2024)

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 26 September, 2024; v1 submitted 31 May, 2024; originally announced June 2024.

Comments: 21 pages, 8 figures, v2: minor changes, published version, 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, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Jeremy George Baier, 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, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Heling Deng, Timothy Dolch , et al. (75 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~\mathrm{nHz}$ break opens a window for PTAs to detect continuous GWs from individual SMBHBs or GWs from the early universe. △ Less

Submitted 19 November, 2024; v1 submitted 10 April, 2024; originally announced April 2024.

Comments: 13 pages, 8 figures, 2 appendices, in press at ApJ

Journal ref: The Astrophysical Journal, Volume 978, Number 1 (2024)

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

Journal ref: Astrophysical Journal, 2024

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 Methods

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: Pulsar timing arrays (PTAs) use an array of millisecond pulsars to search for gravitational waves in the nanohertz regime in pulse time of arrival data. This paper presents rigorous tests of PTA methods, examining their consistency across the relevant parameter space. We discuss updates to the 15-year isotropic gravitational-wave background analyses and their corresponding code representations. De… ▽ More Pulsar timing arrays (PTAs) use an array of millisecond pulsars to search for gravitational waves in the nanohertz regime in pulse time of arrival data. This paper presents rigorous tests of PTA methods, examining their consistency across the relevant parameter space. 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 Enterprise and 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 in gravitational wave searches, reducing the required resources facilitating the computation of Bayes factors via thermodynamic integration and sampling a large number of realizations for computing Bayesian false-alarm probabilities. We perform stringent tests of consistency and correctness of the Bayesian and frequentist analysis methods. For the Bayesian analysis, we test prior recovery, simulation recovery, and Bayes factors. For the frequentist analysis, we test that the 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 analyses of PTA data. △ Less

Submitted 12 May, 2025; 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 comments@nanograv.org; updated to journal version

Journal ref: Phys. Rev. D 109, 103012 (2024)

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 comments@nanograv.org

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 comments@nanograv.org. 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 comments@nanograv.org

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 comments@nanograv.org

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 comments@nanograv.org

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 comments@nanograv.org

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