Rapidly evaluating the compact-binary likelihood function via interpolation

R. J.E. Smith, Chad Richard Hanna, I. Mandel, A. Vecchio

Research output: Contribution to journalArticle

3 Citations (Scopus)

Abstract

Bayesian parameter estimation on gravitational waves from compact-binary coalescences (CBCs) typically requires millions of template waveform computations at different values of the parameters describing the binary. Sampling techniques such as Markov chain Monte Carlo and nested sampling evaluate likelihoods and, hence, compute template waveforms, serially; thus, the total computational time of the analysis scales linearly with that of template generation. Here we address the issue of rapidly computing the likelihood function of CBC sources with nonspinning components. We show how to efficiently compute the continuous likelihood function on the three-dimensional subspace of parameters on which it has a nontrivial dependence - the chirp mass, symmetric mass ratio and coalescence time - via interpolation. Subsequently, sampling this interpolated likelihood function is a significantly cheaper computational process than directly evaluating the likelihood; we report improvements in computational time of two to three orders of magnitude while keeping likelihoods accurate to 0.025%. Generating the interpolant of the likelihood function over a significant portion of the CBC mass space is computationally expensive but highly parallelizable, so the wall time can be very small relative to the time of a full parameter-estimation analysis.

Original languageEnglish (US)
Article number044074
JournalPhysical Review D - Particles, Fields, Gravitation and Cosmology
Volume90
Issue number4
DOIs
StatePublished - Aug 29 2014

Fingerprint

interpolation
coalescing
templates
sampling
waveforms
Markov chains
chirp
mass ratios
gravitational waves

All Science Journal Classification (ASJC) codes

  • Nuclear and High Energy Physics
  • Physics and Astronomy (miscellaneous)

Cite this

@article{eb1baf8bce934fa5a19e18cda71c77a7,
title = "Rapidly evaluating the compact-binary likelihood function via interpolation",
abstract = "Bayesian parameter estimation on gravitational waves from compact-binary coalescences (CBCs) typically requires millions of template waveform computations at different values of the parameters describing the binary. Sampling techniques such as Markov chain Monte Carlo and nested sampling evaluate likelihoods and, hence, compute template waveforms, serially; thus, the total computational time of the analysis scales linearly with that of template generation. Here we address the issue of rapidly computing the likelihood function of CBC sources with nonspinning components. We show how to efficiently compute the continuous likelihood function on the three-dimensional subspace of parameters on which it has a nontrivial dependence - the chirp mass, symmetric mass ratio and coalescence time - via interpolation. Subsequently, sampling this interpolated likelihood function is a significantly cheaper computational process than directly evaluating the likelihood; we report improvements in computational time of two to three orders of magnitude while keeping likelihoods accurate to 0.025{\%}. Generating the interpolant of the likelihood function over a significant portion of the CBC mass space is computationally expensive but highly parallelizable, so the wall time can be very small relative to the time of a full parameter-estimation analysis.",
author = "Smith, {R. J.E.} and Hanna, {Chad Richard} and I. Mandel and A. Vecchio",
year = "2014",
month = "8",
day = "29",
doi = "10.1103/PhysRevD.90.044074",
language = "English (US)",
volume = "90",
journal = "Physical Review D",
issn = "0556-2821",
number = "4",

}

Rapidly evaluating the compact-binary likelihood function via interpolation. / Smith, R. J.E.; Hanna, Chad Richard; Mandel, I.; Vecchio, A.

In: Physical Review D - Particles, Fields, Gravitation and Cosmology, Vol. 90, No. 4, 044074, 29.08.2014.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Rapidly evaluating the compact-binary likelihood function via interpolation

AU - Smith, R. J.E.

AU - Hanna, Chad Richard

AU - Mandel, I.

AU - Vecchio, A.

PY - 2014/8/29

Y1 - 2014/8/29

N2 - Bayesian parameter estimation on gravitational waves from compact-binary coalescences (CBCs) typically requires millions of template waveform computations at different values of the parameters describing the binary. Sampling techniques such as Markov chain Monte Carlo and nested sampling evaluate likelihoods and, hence, compute template waveforms, serially; thus, the total computational time of the analysis scales linearly with that of template generation. Here we address the issue of rapidly computing the likelihood function of CBC sources with nonspinning components. We show how to efficiently compute the continuous likelihood function on the three-dimensional subspace of parameters on which it has a nontrivial dependence - the chirp mass, symmetric mass ratio and coalescence time - via interpolation. Subsequently, sampling this interpolated likelihood function is a significantly cheaper computational process than directly evaluating the likelihood; we report improvements in computational time of two to three orders of magnitude while keeping likelihoods accurate to 0.025%. Generating the interpolant of the likelihood function over a significant portion of the CBC mass space is computationally expensive but highly parallelizable, so the wall time can be very small relative to the time of a full parameter-estimation analysis.

AB - Bayesian parameter estimation on gravitational waves from compact-binary coalescences (CBCs) typically requires millions of template waveform computations at different values of the parameters describing the binary. Sampling techniques such as Markov chain Monte Carlo and nested sampling evaluate likelihoods and, hence, compute template waveforms, serially; thus, the total computational time of the analysis scales linearly with that of template generation. Here we address the issue of rapidly computing the likelihood function of CBC sources with nonspinning components. We show how to efficiently compute the continuous likelihood function on the three-dimensional subspace of parameters on which it has a nontrivial dependence - the chirp mass, symmetric mass ratio and coalescence time - via interpolation. Subsequently, sampling this interpolated likelihood function is a significantly cheaper computational process than directly evaluating the likelihood; we report improvements in computational time of two to three orders of magnitude while keeping likelihoods accurate to 0.025%. Generating the interpolant of the likelihood function over a significant portion of the CBC mass space is computationally expensive but highly parallelizable, so the wall time can be very small relative to the time of a full parameter-estimation analysis.

UR - http://www.scopus.com/inward/record.url?scp=84908436170&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84908436170&partnerID=8YFLogxK

U2 - 10.1103/PhysRevD.90.044074

DO - 10.1103/PhysRevD.90.044074

M3 - Article

AN - SCOPUS:84908436170

VL - 90

JO - Physical Review D

JF - Physical Review D

SN - 0556-2821

IS - 4

M1 - 044074

ER -