Nested Sampling as a tool for LISA data analysis

Nested sampling is a technique for efficiently computing the probability of a data set under a particular hypothesis, also called the Bayesian Evidence or Marginal Likelihood, and for evaluating the posterior. MULTINEST is a multi-modal nested sampling algorithm which has been designed to efficiently explore and characterize posterior probability surfaces containing multiple secondary solutions. We have applied the MULTINEST algorithm to a number of problems in gravitational wave data analysis. In this article, we describe the algorithm and present results for several applications of the algorithm to analysis of mock LISA data. We summarise recently published results for a test case in which we searched for two non-spinning black hole binary merger signals in simulated LISA data. We also describe results obtained with MULTINEST in the most recent round of the Mock LISA Data Challenge (MLDC), in which the algorithm was used to search for and characterise both spinning supermassive black hole binary inspirals and bursts from cosmic string cusps. In all these applications, the algorithm found the correct number of signals and efficiently recovered the posterior probability distribution. Moreover, in most cases the waveform corresponding to the best a-posteriori parameters had an overlap in excess of 99% with the true signal.

[1]  Bernard F. Schutz,et al.  Data analysis of gravitational-wave signals from spinning neutron stars. I. The signal and its detection , 1998 .

[2]  N. Cornish,et al.  The LISA response function , 2002, gr-qc/0209011.

[3]  The search for massive black hole binaries with LISA , 2006, gr-qc/0612091.

[4]  Report on the second Mock LISA Data Challenge , 2007, 0711.2667.

[5]  J. Gair,et al.  The Mock LISA Data Challenges: from Challenge 1B to Challenge 3 , 2008, 0806.2110.

[6]  F. Feroz,et al.  Bayesian selection of sign μ within mSUGRA in global fits including WMAP5 results , 2008, 0807.4512.

[7]  A. Vecchio,et al.  Assigning confidence to inspiral gravitational wave candidates with Bayesian model selection , 2008, 0807.4483.

[8]  F. Feroz,et al.  The impact of priors and observables on parameter inferences in the constrained MSSM , 2008, 0809.3792.

[9]  F. Feroz,et al.  Bayesian modelling of clusters of galaxies from multifrequency‐pointed Sunyaev–Zel'dovich observations , 2008, 0811.1199.

[10]  J. Gair,et al.  Cosmic swarms: a search for supermassive black holes in the LISA data stream with a hybrid evolutionary algorithm , 2009, 0903.3733.

[11]  J. Gair,et al.  An algorithm for the detection of extreme mass ratio inspirals in LISA data , 2009, 0902.4133.

[12]  F. Feroz,et al.  MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physics , 2008, 0809.3437.

[13]  F. Feroz,et al.  Use of the MULTINEST algorithm for gravitational wave data analysis , 2009, 0904.1544.

[14]  Characterizing the Gravitational Wave Signature from Cosmic String Cusps , 2008, 0812.1590.