Building a stochastic template bank for detecting massive black hole binaries

Coalescence of two massive black holes is the strongest and most promising source for LISA. In fact, gravitational signal from the end of inspiral and merger will be detectable throughout the Universe. In this article we describe the first step in the two-step hierarchical search for gravitational wave signal from the inspiraling massive BH binaries. It is based on the routinely used in the ground base gravitational wave astronomy method of filtering the data through the bank of templates. However we use a novel Monte-Carlo based (stochastic) method to lay a grid in the parameter space, and we use the likelihood maximized analytically over some parameters, known as F-statistic, as a detection statistic. We build a coarse template bank to detect gravitational wave signals and to make preliminary parameter estimation. The best candidates will be followed up using Metropolis-Hasting stochastic search to refine the parameter estimation. We demonstrate the performance of the method by applying it to the Mock LISA data challenge 1B (training data set).Comment: revtex4, 8 figure

Evolutionary quantum cosmology in a gauge-fixed picture

We study the classical and quantum models of a flat Friedmann-Robertson-Walker (FRW) space-time, coupled to a perfect fluid, in the context of the consensus and a gauge-fixed Lagrangian frameworks. It is shown that, either in the usual or in the gauge-fixed actions, the evolution of the universe based on the classical cosmology represents a late time power law expansion, coming from a big-bang singularity in which the scale factor goes to zero for the standard matter, and tending towards a big-rip singularity in which the scale factor diverges for the phantom fluid. We then employ the familiar canonical quantization procedure in the given cosmological setting to find the cosmological wave functions in the corresponding minisuperspace. Using a gauge-fixed (reduced) Lagrangian, we show that, it may lead to a Schr\"{o}dinger equation for the quantum-mechanical description of the model under consideration, the eigenfunctions of which can be used to construct the time dependent wave function of the universe. We use the resulting wave function in order to investigate the possibility of the avoidance of classical singularities due to quantum effects by means of the many-worlds and ontological interpretation of quantum cosmology.Comment: 15 pages, 10 figures, typos corrected, Refs. adde

Detection Strategies for Extreme Mass Ratio Inspirals

The capture of compact stellar remnants by galactic black holes provides a unique laboratory for exploring the near horizon geometry of the Kerr spacetime, or possible departures from general relativity if the central cores prove not to be black holes. The gravitational radiation produced by these Extreme Mass Ratio Inspirals (EMRIs) encodes a detailed map of the black hole geometry, and the detection and characterization of these signals is a major scientific goal for the LISA mission. The waveforms produced are very complex, and the signals need to be coherently tracked for hundreds to thousands of cycles to produce a detection, making EMRI signals one of the most challenging data analysis problems in all of gravitational wave astronomy. Estimates for the number of templates required to perform an exhaustive grid-based matched-filter search for these signals are astronomically large, and far out of reach of current computational resources. Here I describe an alternative approach that employs a hybrid between Genetic Algorithms and Markov Chain Monte Carlo techniques, along with several time saving techniques for computing the likelihood function. This approach has proven effective at the blind extraction of relatively weak EMRI signals from simulated LISA data sets.Comment: 10 pages, 4 figures, Updated for LISA 8 Symposium Proceeding

A hierarchical search for gravitational waves from supermassive black hole binary mergers

We present a method to search for gravitational waves from coalescing supermassive binary black holes in LISA data. The search utilizes the $\mathcal{F}$-statistic to maximize over, and determine the values of, the extrinsic parameters of the binary system. The intrinsic parameters are searched over hierarchically using stochastically generated multi-dimensional template banks to recover the masses and sky locations of the binary. We present the results of this method applied to the mock LISA data Challenge 1B data set.Comment: 11 pages, 2 figures, for GWDAW-12 proceedings edition of CQ

Cosmic Swarms: A search for Supermassive Black Holes in the LISA data stream with a Hybrid Evolutionary Algorithm

We describe a hybrid evolutionary algorithm that can simultaneously search for multiple supermassive black hole binary (SMBHB) inspirals in LISA data. The algorithm mixes evolutionary computation, Metropolis-Hastings methods and Nested Sampling. The inspiral of SMBHBs presents an interesting problem for gravitational wave data analysis since, due to the LISA response function, the sources have a bi-modal sky solution. We show here that it is possible not only to detect multiple SMBHBs in the data stream, but also to investigate simultaneously all the various modes of the global solution. In all cases, the algorithm returns parameter determinations within $5\sigma$ (as estimated from the Fisher Matrix) of the true answer, for both the actual and antipodal sky solutions.Comment: submitted to Classical & Quantum Gravity. 19 pages, 4 figure

A How-To for the Mock LISA Data Challenges

The LISA International Science Team Working Group on Data Analysis (LIST-WG1B) is sponsoring several rounds of mock data challenges, with the purpose of fostering development of LISA data-analysis capabilities, and of demonstrating technical readiness for the maximum science exploitation of the LISA data. The first round of challenge data sets were released at this Symposium. We describe the models and conventions (for LISA and for gravitational-wave sources) used to prepare the data sets, the file format used to encode them, and the tools and resources available to support challenge participants.Comment: 10 pages, 1 figure, in Proceedings of the Sixth International LISA Symposium (AIP, 2006

Random template placement and prior information

In signal detection problems, one is usually faced with the task of searching a parameter space for peaks in the likelihood function which indicate the presence of a signal. Random searches have proven to be very efficient as well as easy to implement, compared e.g. to searches along regular grids in parameter space. Knowledge of the parameterised shape of the signal searched for adds structure to the parameter space, i.e., there are usually regions requiring to be densely searched while in other regions a coarser search is sufficient. On the other hand, prior information identifies the regions in which a search will actually be promising or may likely be in vain. Defining specific figures of merit allows one to combine both template metric and prior distribution and devise optimal sampling schemes over the parameter space. We show an example related to the gravitational wave signal from a binary inspiral event. Here the template metric and prior information are particularly contradictory, since signals from low-mass systems tolerate the least mismatch in parameter space while high-mass systems are far more likely, as they imply a greater signal-to-noise ratio (SNR) and hence are detectable to greater distances. The derived sampling strategy is implemented in a Markov chain Monte Carlo (MCMC) algorithm where it improves convergence.Comment: Proceedings of the 8th Edoardo Amaldi Conference on Gravitational Waves. 7 pages, 4 figure

The Mock LISA Data Challenges: from Challenge 3 to Challenge 4

The Mock LISA Data Challenges are a program to demonstrate LISA data-analysis capabilities and to encourage their development. Each round of challenges consists of one or more datasets containing simulated instrument noise and gravitational waves from sources of undisclosed parameters. Participants analyze the datasets and report best-fit solutions for the source parameters. Here we present the results of the third challenge, issued in Apr 2008, which demonstrated the positive recovery of signals from chirping Galactic binaries, from spinning supermassive--black-hole binaries (with optimal SNRs between ~ 10 and 2000), from simultaneous extreme-mass-ratio inspirals (SNRs of 10-50), from cosmic-string-cusp bursts (SNRs of 10-100), and from a relatively loud isotropic background with Omega_gw(f) ~ 10^-11, slightly below the LISA instrument noise.Comment: 12 pages, 2 figures, proceedings of the 8th Edoardo Amaldi Conference on Gravitational Waves, New York, June 21-26, 200

Classifying LISA gravitational wave burst signals using Bayesian evidence

We consider the problem of characterisation of burst sources detected with the Laser Interferometer Space Antenna (LISA) using the multi-modal nested sampling algorithm, MultiNest. We use MultiNest as a tool to search for modelled bursts from cosmic string cusps, and compute the Bayesian evidence associated with the cosmic string model. As an alternative burst model, we consider sine-Gaussian burst signals, and show how the evidence ratio can be used to choose between these two alternatives. We present results from an application of MultiNest to the last round of the Mock LISA Data Challenge, in which we were able to successfully detect and characterise all three of the cosmic string burst sources present in the release data set. We also present results of independent trials and show that MultiNest can detect cosmic string signals with signal-to-noise ratio (SNR) as low as ~7 and sine-Gaussian signals with SNR as low as ~8. In both cases, we show that the threshold at which the sources become detectable coincides with the SNR at which the evidence ratio begins to favour the correct model over the alternative.Comment: 21 pages, 11 figures, accepted by CQG; v2 has minor changes for consistency with accepted versio

LISA as a dark energy probe

Recently it was shown that the inclusion of higher signal harmonics in the inspiral signals of binary supermassive black holes (SMBH) leads to dramatic improvements in parameter estimation with the Laser Interferometer Space Antenna (LISA). In particular, the angular resolution becomes good enough to identify the host galaxy or galaxy cluster, in which case the redshift can be determined by electromagnetic means. The gravitational wave signal also provides the luminosity distance with high accuracy, and the relationship between this and the redshift depends sensitively on the cosmological parameters, such as the equation-of-state parameter $w=p_{\rm DE}/\rho_{\rm DE}$ of dark energy. With a single binary SMBH event at $z < 1$ having appropriate masses and orientation, one would be able to constrain $w$ to within a few percent. We show that, if the measured sky location is folded into the error analysis, the uncertainty on $w$ goes down by an additional factor of 2-3, leaving weak lensing as the only limiting factor in using LISA as a dark energy probe.Comment: 11pages, 1 Table, minor changes in text, accepted for publication in Classical and Quantum Gravity (special issue for proceedings of 7th LISA symposium