by Cohen, Michael I. and Cutler, Curt and Vallisneri, Michele
Submitted to CQG; 28 pages, 10 figures; higher-resolution plots available at http://www.vallis.org/publications/cosmicstrings

A network of observable, macroscopic cosmic (super-)strings may have formed in the early universe. If so, the cusps that generically develop on cosmic-string loops emit bursts of gravitational radiation that could be detectable by both ground- and space-based gravitational-wave interferometers. Here we report on two versions of a LISA-oriented string-burst search pipeline that we have developed and tested within the context of the Mock LISA Data Challenges. The two versions rely on the publicly available MultiNest and PyMC software packages, respectively. To reduce the effective dimensionality of the search space, our implementations use the F-statistic to analytically maximize over the signal’s amplitude and polarization, A and psi, and use the FFT to search quickly over burst arrival times t_C. The standard F-statistic is essentially a frequentist statistic that maximizes the likelihood; we also demonstrate an approximate, Bayesian version of the F-statistic that incorporates realistic priors on A and psi. We calculate how accurately LISA can expect to measure the physical parameters of string-burst sources. To understand LISA’s angular resolution for string-burst sources, we draw maps of the waveform fitting factor [maximized over (A psi, t_C)] as a function of sky position; these maps dramatically illustrate why (for LISA) inferring the correct sky location of the emitting string loop will often be practically impossible. We also identify and elucidate several symmetries that are imbedded in this search problem, and we derive the distribution of cut-off frequencies f_max for observable bursts.

by Manca, Gian Mario and Vallisneri, Michele
RevTeX4, 21 pages, 9 PDF figures

The efficient placement of signal templates in source-parameter space is a crucial requisite for exhaustive matched-filtering searches of modeled gravitational-wave sources. Unfortunately, the current placement algorithms based on regular parameter-space meshes are difficult to generalize beyond simple signal models with few parameters. Various authors have suggested that a general, flexible, yet efficient alternative can be found in randomized placement strategies such as random placement and stochastic placement, which enhances random placement by selectively rejecting templates that are too close to others. In this article we explore several theoretical and practical issues in randomized placement: the size and performance of the resulting template banks; the effects of parameter-space boundaries; the use of quasi-random (self avoiding) number sequences; most important, the implementation of these algorithms in curved signal manifolds with and without the use of a Riemannian signal metric, which may be difficult to obtain. Specifically, we show how the metric can be replaced with a discrete triangulation-based representation of local geometry. We argue that the broad class of randomized placement algorithms offers a promising answer to many search problems, but that the specific choice of a scheme and its implementation details will still need to be fine-tuned separately for each problem.

by McWilliams, Sean T. and Thorpe, James Ira and Baker, John G. and Kelly, Bernard J.
16 pages, 9 figures, submitted to Phys. Rev. D

We investigate the precision with which the parameters describing the characteristics and location of nonspinning black hole binaries can be measured with the Laser Interferometer Space Antenna (LISA). By using complete waveforms including the inspiral, merger and ringdown portions of the signals, we find that LISA will have far greater precision than previous estimates for nonspinning mergers that ignored the merger and ringdown. Our analysis covers nonspinning waveforms with moderate mass ratios, q >= 1/10, and total masses 10^5 < M/M_{Sun} < 10^7. We compare the parameter uncertainties using the Fisher matrix formalism, and establish the significance of mass asymmetry and higher-order content to the predicted parameter uncertainties resulting from inclusion of the merger. In real-time observations, the later parts of the signal lead to significant improvements in sky-position precision in the last hours and even the final minutes of observation. For comparable mass systems with total mass M/M_{Sun} = ~10^6, we find that the increased precision resulting from including the merger is comparable to the increase in signal-to-noise ratio. For the most precise systems under investigation, half can be localized to within O(10 arcmin), and 10% can be localized to within O(1 arcmin).

by Trias, Miquel and Vecchio, Alberto and Veitch, John
11 pages, 4 figures, submitted to GWDAW-13 proceedings

Bayesian analysis of LISA data sets based on Markov chain Monte Carlo methods has been shown to be a challenging problem, in part due to the complicated structure of the likelihood function consisting of several isolated local maxima that dramatically reduces the efficiency of the sampling techniques. Here we introduce a new fully Markovian algorithm, a Delayed Rejection Metropolis-Hastings Markov chain Monte Carlo method, to efficiently explore these kind of structures and we demonstrate its performance on selected LISA data sets containing a known number of stellar-mass binary signals embedded in Gaussian stationary noise.

by Gair, Jonathan R. and Porter, Edward K.
submitted to Classical & Quantum Gravity. 19 pages, 4 figures

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 $latex 5\sigma$ (as estimated from the Fisher Matrix) of the true answer, for both the actual and antipodal sky solutions.

by Babak, Stanislav and Gair, Jonathan R. and Porter, Edward K.
14 pages, 4 figures

The gravitational wave signal from a compact object spiralling toward a massive black hole (MBH) is thought to be one of the most difficult sources to detect in the LISA data stream. Due to the large parameter space of possible signals and many orbital cycles spent in the sensitivity band of LISA, it has been estimated previously that of the order of 10^{35} templates would be required for a fully coherent search with a template grid, which is computationally impossible. Here we describe an algorithm based on a constrained Metropolis-Hastings stochastic search which allows us to find and accurately estimate parameters of isolated EMRI signals buried in Gaussian instrumental noise. We illustrate the effectiveness of the algorithm with results from searches of the Mock LISA Data Challenge round 1B data sets.