A Markov Chain Monte Carlo approach to the study of massive black hole binary systems with LISA

The Laser Interferometer Space Antenna (LISA) will produce a data stream containing a vast number of overlapping sources: from strong signals generated by the coalescence of massive black hole binary systems to much weaker radiation form sub-stellar mass compact binaries and extreme-mass ratio inspirals. It has been argued that the observation of weak signals could be hampered by the presence of loud ones and that they first need to be removed to allow such observations. Here we consider a different approach in which sources are studied simultaneously within the framework of Bayesian inference. We investigate the simplified case in which the LISA data stream contains radiation from a massive black hole binary system superimposed over a (weaker) quasi-monochromatic waveform generated by a white dwarf binary. We derive the posterior probability density function of the model parameters using an automatic Reversible Jump Markov Chain Monte Carlo algorithm (RJMCMC). We show that the information about the sources and noise are retrieved at the expected level of accuracy without the need of removing the stronger signal. Our analysis suggests that this approach is worth pursuing further and should be considered for the actual analysis of the LISA data.Comment: submitted to cqg as GWDAW-10 conference proceedings, 10 pages, 4 figures, some changes to plots and numerical detail

LISA Response Function and Parameter Estimation

We investigate the response function of LISA and consider the adequacy of its commonly used approximation in the high-frequency range of the observational band. We concentrate on monochromatic binary systems, such as white dwarf binaries. We find that above a few mHz the approxmation starts becoming increasingly inaccurate. The transfer function introduces additional amplitude and phase modulations in the measured signal that influence parameter estmation and, if not properly accounted for, lead to losses of signal-to-noise ratio.Comment: 4 pages, 2 figures, amaldi 5 conference proceeding

Detecting extreme mass ratio inspiral events in LISA data using the Hierarchical Algorithm for Clusters and Ridges (HACR)

One of the most exciting prospects for the Laser Interferometer Space Antenna (LISA) is the detection of gravitational waves from the inspirals of stellar-mass compact objects into supermassive black holes. Detection of these sources is an extremely challenging computational problem due to the large parameter space and low amplitude of the signals. However, recent work has suggested that the nearest extreme mass ratio inspiral (EMRI) events will be sufficiently loud that they might be detected using computationally cheap, template-free techniques, such as a time-frequency analysis. In this paper, we examine a particular time-frequency algorithm, the Hierarchical Algorithm for Clusters and Ridges (HACR). This algorithm searches for clusters in a power map and uses the properties of those clusters to identify signals in the data. We find that HACR applied to the raw spectrogram performs poorly, but when the data is binned during the construction of the spectrogram, the algorithm can detect typical EMRI events at distances of up to $\sim2.6$Gpc. This is a little further than the simple Excess Power method that has been considered previously. We discuss the HACR algorithm, including tuning for single and multiple sources, and illustrate its performance for detection of typical EMRI events, and other likely LISA sources, such as white dwarf binaries and supermassive black hole mergers. We also discuss how HACR cluster properties could be used for parameter extraction.Comment: 21 pages, 11 figures, submitted to Class. Quantum Gravity. Modified and shortened in light of referee's comments. Updated results consider tuning over all three HACR thresholds, and show 10-15% improvement in detection rat

Optimal statistic for detecting gravitational wave signals from binary inspirals with LISA

A binary compact object early in its inspiral phase will be picked up by its nearly monochromatic gravitational radiation by LISA. But even this innocuous appearing candidate poses interesting detection challenges. The data that will be scanned for such sources will be a set of three functions of LISA's twelve data streams obtained through time-delay interferometry, which is necessary to cancel the noise contributions from laser-frequency fluctuations and optical-bench motions to these data streams. We call these three functions pseudo-detectors. The sensitivity of any pseudo-detector to a given sky position is a function of LISA's orbital position. Moreover, at a given point in LISA's orbit, each pseudo-detector has a different sensitivity to the same sky position. In this work, we obtain the optimal statistic for detecting gravitational wave signals, such as from compact binaries early in their inspiral stage, in LISA data. We also present how the sensitivity of LISA, defined by this optimal statistic, varies as a function of sky position and LISA's orbital location. Finally, we show how a real-time search for inspiral signals can be implemented on the LISA data by constructing a bank of templates in the sky positions.Comment: 22 pages, 15 eps figures, Latex, uses iopart style/class files. Based on talk given at the 8th Gravitational Wave Data Analysis Workshop, Milwaukee, USA, December 17-20, 2003. Accepted for publication in Class. Quant. Gra

Extracting galactic binary signals from the first round of Mock LISA Data Challenges

We report on the performance of an end-to-end Bayesian analysis pipeline for detecting and characterizing galactic binary signals in simulated LISA data. Our principal analysis tool is the Blocked-Annealed Metropolis Hasting (BAM) algorithm, which has been optimized to search for tens of thousands of overlapping signals across the LISA band. The BAM algorithm employs Bayesian model selection to determine the number of resolvable sources, and provides posterior distribution functions for all the model parameters. The BAM algorithm performed almost flawlessly on all the Round 1 Mock LISA Data Challenge data sets, including those with many highly overlapping sources. The only misses were later traced to a coding error that affected high frequency sources. In addition to the BAM algorithm we also successfully tested a Genetic Algorithm (GA), but only on data sets with isolated signals as the GA has yet to be optimized to handle large numbers of overlapping signals.Comment: 13 pages, 4 figures, submitted to Proceedings of GWDAW-11 (Berlin, Dec. '06

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 Three-Stage Search for Supermassive Black Hole Binaries in LISA Data

Gravitational waves from the inspiral and coalescence of supermassive black-hole (SMBH) binaries with masses ~10^6 Msun are likely to be among the strongest sources for the Laser Interferometer Space Antenna (LISA). We describe a three-stage data-analysis pipeline designed to search for and measure the parameters of SMBH binaries in LISA data. The first stage uses a time-frequency track-search method to search for inspiral signals and provide a coarse estimate of the black-hole masses m_1, m_2 and of the coalescence time of the binary t_c. The second stage uses a sequence of matched-filter template banks, seeded by the first stage, to improve the measurement accuracy of the masses and coalescence time. Finally, a Markov Chain Monte Carlo search is used to estimate all nine physical parameters of the binary. Using results from the second stage substantially shortens the Markov Chain burn-in time and allows us to determine the number of SMBH-binary signals in the data before starting parameter estimation. We demonstrate our analysis pipeline using simulated data from the first LISA Mock Data Challenge. We discuss our plan for improving this pipeline and the challenges that will be faced in real LISA data analysis.Comment: 12 pages, 3 figures, submitted to Proceedings of GWDAW-11 (Berlin, Dec. '06

Intermediate and extreme mass-ratio inspirals — astrophysics, science applications and detection using LISA

Black hole binaries with extreme (gtrsim104:1) or intermediate (~102–104:1) mass ratios are among the most interesting gravitational wave sources that are expected to be detected by the proposed laser interferometer space antenna (LISA). These sources have the potential to tell us much about astrophysics, but are also of unique importance for testing aspects of the general theory of relativity in the strong field regime. Here we discuss these sources from the perspectives of astrophysics, data analysis and applications to testing general relativity, providing both a description of the current state of knowledge and an outline of some of the outstanding questions that still need to be addressed. This review grew out of discussions at a workshop in September 2006 hosted by the Albert Einstein Institute in Golm, Germany