Detecting gravitational radiation from neutron stars using a six-parameter adaptive MCMC method

We present a Markov chain Monte Carlo technique for detecting gravitational radiation from a neutron star in laser interferometer data. The algorithm can estimate up to six unknown parameters of the target, including the rotation frequency and frequency derivative, using reparametrization, delayed rejection and simulated annealing. We highlight how a simple extension of the method, distributed over multiple computer processors, will allow for a search over a narrow frequency band. The ultimate goal of this research is to search for sources at a known locations, but uncertain spin parameters, such as may be found in SN1987A.Comment: Submitted to Classical and Quantum Gravity for GWDAW-8 proceeding

Inference on inspiral signals using LISA MLDC data

In this paper we describe a Bayesian inference framework for analysis of data obtained by LISA. We set up a model for binary inspiral signals as defined for the Mock LISA Data Challenge 1.2 (MLDC), and implemented a Markov chain Monte Carlo (MCMC) algorithm to facilitate exploration and integration of the posterior distribution over the 9-dimensional parameter space. Here we present intermediate results showing how, using this method, information about the 9 parameters can be extracted from the data.Comment: Accepted for publication in Classical and Quantum Gravity, GWDAW-11 special issu

A time-domain MCMC search and upper limit technique for gravitational waves of uncertain frequency from a targeted neutron star

It is computationally expensive to search the large parameter space associated with a gravitational wave signal of uncertain frequency, such as might be expected from the possible pulsar generated by SN1987A. To address this difficulty we have developed a Markov Chain Monte Carlo method that performs a time-domain Bayesian search for a signal over a 4 Hz frequency band and a spindown of magnitude of up to 1 × 10−9 Hz s−1. We use Monte Carlo simulations to set upper limits on signal amplitude with this technique, which we intend to apply to a gravitational wave search