An efficient pipeline for joint gravitational wave searches from individual binaries and a gravitational wave background with Hamiltonian sampling

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Freedman, Gabriel E., Vigeland, Sarah J.
Format: Preprint
Published: 2024
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866910528709853184
author Freedman, Gabriel E.
Vigeland, Sarah J.
author_facet Freedman, Gabriel E.
Vigeland, Sarah J.
contents The pulsar timing array community has recently reported the first evidence of a low-frequency stochastic gravitational wave background. With longer observational timespans we expect to be able to resolve individual gravitational wave sources in our data alongside the background signal. The statistical modeling and Bayesian searches for such individual signals is a computationally taxing task that is the focus of many different avenues of methods development. We present a pipeline for performing efficient joint searches for gravitational waves originating from individual supermassive black hole binaries as well as a gravitational wave background using a Hamiltonian Monte Carlo sampling scheme. Hamiltonian sampling proposes samples based on the gradients of the model likelihood, and can both converge faster to more complicated and high-dimensional distributions as well as efficiently explore highly covariant parameter spaces such as the joint gravitational wave background and individual binary model. We show the effectiveness of our scheme by demonstrating accurate parameter estimation for simulated datasets containing low- (6 nHz) or high- (60 nHz) frequency binary sources. Additionally we show that our method is capable at more efficiently generating skymaps for individual binary sources -- maps displaying the upper limits on the gravitational wave strain of the source, $h_{0}$, as a function of sky location -- by sampling over larger portions of the full sky. Comparing against results for the NANOGrav 12.5-year dataset, we find similar reconstructed upper limits on the gravitational wave strain while simultaneously reducing the number of required analyses from 72 independent binned searches down to a single run.
format Preprint
id arxiv_https___arxiv_org_abs_2407_11135
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle An efficient pipeline for joint gravitational wave searches from individual binaries and a gravitational wave background with Hamiltonian sampling
Freedman, Gabriel E.
Vigeland, Sarah J.
Instrumentation and Methods for Astrophysics
General Relativity and Quantum Cosmology
The pulsar timing array community has recently reported the first evidence of a low-frequency stochastic gravitational wave background. With longer observational timespans we expect to be able to resolve individual gravitational wave sources in our data alongside the background signal. The statistical modeling and Bayesian searches for such individual signals is a computationally taxing task that is the focus of many different avenues of methods development. We present a pipeline for performing efficient joint searches for gravitational waves originating from individual supermassive black hole binaries as well as a gravitational wave background using a Hamiltonian Monte Carlo sampling scheme. Hamiltonian sampling proposes samples based on the gradients of the model likelihood, and can both converge faster to more complicated and high-dimensional distributions as well as efficiently explore highly covariant parameter spaces such as the joint gravitational wave background and individual binary model. We show the effectiveness of our scheme by demonstrating accurate parameter estimation for simulated datasets containing low- (6 nHz) or high- (60 nHz) frequency binary sources. Additionally we show that our method is capable at more efficiently generating skymaps for individual binary sources -- maps displaying the upper limits on the gravitational wave strain of the source, $h_{0}$, as a function of sky location -- by sampling over larger portions of the full sky. Comparing against results for the NANOGrav 12.5-year dataset, we find similar reconstructed upper limits on the gravitational wave strain while simultaneously reducing the number of required analyses from 72 independent binned searches down to a single run.
title An efficient pipeline for joint gravitational wave searches from individual binaries and a gravitational wave background with Hamiltonian sampling
topic Instrumentation and Methods for Astrophysics
General Relativity and Quantum Cosmology
url https://arxiv.org/abs/2407.11135