by Cornish, Neil J.
4 pages, no figures

Theoretical studies in gravitational wave astronomy often require the calculation of Fisher Information Matrices and Likelihood functions, which in a direct approach entail the costly step of computing gravitational waveforms. Here I describe an alternative technique that sidesteps the need to compute full waveforms, resulting in significant computational savings. I describe how related techniques can be used to speed up Bayesian inference applied to real gravitational wave data.

by Key, Joey Shapiro and Cornish, Neil J.
11 pages, 19 figures

The Laser Interferometer Space Antenna (LISA) is designed to detect gravitational wave signals from astrophysical sources, including those from coalescing binary systems of compact objects such as black holes. Colliding galaxies have central black holes that sink to the center of the merged galaxy and begin to orbit one another and emit gravitational waves. Some galaxy evolution models predict that the binary black hole system will enter the LISA band with significant orbital eccentricity, while other models suggest that the orbits will already have circularized. Using a full seventeen parameter waveform model that includes the effects of orbital eccentricity, spin precession and higher harmonics, we investigate how well the source parameters can be inferred from simulated LISA data. Defining the reference eccentricity as the value one year before merger, we find that for typical LISA sources, it will be possible to measure the eccentricity to an accuracy of parts in a thousand. The accuracy with which the eccentricity can be measured depends only very weakly on the eccentricity, making it possible to distinguish circular orbits from those with very small eccentricities. LISA measurements of the orbital eccentricity can provide strong constraints on theories of galaxy mergers in the early universe.

by Finn, Lee Samuel and Lommen, Andrea N.
43 pages, 13 figures, submitted to ApJ.

Efforts to detect gravitational waves by timing an array of pulsars have focused traditionally on stationary gravitational waves: e.g., stochastic or periodic signals. Gravitational wave bursts — signals whose duration is much shorter than the observation period — will also arise in the pulsar timing array waveband. Sources that give rise to detectable bursts include the formation or coalescence of supermassive black holes (SMBHs), the periapsis passage of compact objects in highly elliptic or unbound orbits about a SMBH, or cusps on cosmic strings. Here we describe how pulsar timing array data may be analyzed to detect and characterize these bursts. Our analysis addresses, in a mutually consistent manner, a hierarchy of three questions: \emph{i}) What are the odds that a dataset includes the signal from a gravitational wave burst? \emph{ii}) Assuming the presence of a burst, what is the direction to its source? and \emph{iii}) Assuming the burst propagation direction, what is the burst waveform’s time dependence in each of its polarization states? Applying our analysis to synthetic data sets we find that we can \emph{detect} gravitational waves even when the radiation is too weak to either localize the source of infer the waveform, and \emph{detect} and \emph{localize} sources even when the radiation amplitude is too weak to permit the waveform to be determined. While the context of our discussion is gravitational wave detection via pulsar timing arrays, the analysis itself is directly applicable to gravitational wave detection using either ground or space-based detector data.

by Cornish, Neil J. and Key, Joey Shapiro
9 pages

Several scenarios have been proposed in which the orbits of binary black holes enter the band of a gravitational wave detector with significant eccentricity. To avoid missing these signals or biasing the parameter estimation it is important that we consider waveform models that account for eccentricity. The ingredients needed to compute post-Newtonian (PN) waveforms produced by spinning black holes inspiralling on quasi-eccentric orbits have been available for almost two decades at 2 PN order, and this work has recently been extended to 2.5 PN order. However, the computational cost of directly implementing these waveforms is high, requiring many steps per orbit to evolve the system of coupled differential equations. Here we employ a separation of timescales and a generalized Keplarian parameterization of the orbits to produce efficient waveforms describing spinning black hole binaries with arbitrary spin orientations on quasi-eccentric orbits to 1.5 PN order. Our solution includes the spin contributions to the decay of the semi-major axis and eccentricity. We outline a scheme for extending our approach to higher post-Newtonian order.

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 Petiteau, Antoine and Shang, Yu and Babak, Stanislav and Feroz, Farhan
25 pages, 9 figures

Coalescing massive Black Hole binaries are the strongest and probably the most important gravitational wave sources in the LISA band. The spin and orbital precessions bring complexity in the waveform and make the likelihood surface richer in structure as compared to the non-spinning case. We introduce an extended multimodal genetic algorithm which utilizes the properties of the signal and the detector response function to analyze the data from the third round of mock LISA data challenge (MLDC 3.2). The performance of this method is comparable, if not better, to already existing algorithms. We have found all five sources present in MLDC 3.2 and recovered the coalescence time, chirp mass, mass ratio and sky location with reasonable accuracy. As for the orbital angular momentum and two spins of the Black Holes, we have found a large number of widely separated modes in the parameter space with similar maximum likelihood values.

by Khanna, Gaurav and McKennon, Justin
14 pages, 4 figures

In this work, we make use of the OpenCL framework to accelerate an EMRI modeling application using the hardware accelerators — Cell BE and Tesla CUDA GPU. We describe these compute technologies and our parallelization approach in detail, present our performance results, and then compare them with those from our previous implementations based on the native CUDA and Cell SDKs. The OpenCL framework allows us to execute identical source-code on both architectures and yet obtain strong performance gains that are comparable to what can be derived from the native SDKs.

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 Feroz, Farhan and Gair, Jonathan R. and Graff, Philip and Hobson, Michael P and Lasenby, Anthony
21 pages, 11 figures, submitted to CQG

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.

by Whelan, John T. and Prix, Reinhard and Khurana, Deepak
26 pages, 11 figures

We describe the application of an F-statistic search for continuous gravitational waves to the search for galactic white-dwarf binaries in the Second Mock LISA Data Challenge. The search method employs a hierarchical template-grid based exploration of the parameter space, using a coincidence step to distinguish between primary (”true”) and secondary maxima, followed by a final (multi-TDI) “zoom” stage to provide an accurate parameter estimation of the final candidates. Suitably tuned, the pipeline is able to extract 1989 true signals with only 5 false alarms. The use of the rigid adiabatic approximation allows recovery of signal parameters comparable to statistical expectations, although there is still some systematic excess above expected statistical errors due to Gaussian noise. An experimental iterative pipeline with seven rounds of subtraction and re-analysis allows us to increase the number of signals recovered, up to a total of 3419 with 29 false alarms.

by Harry, Ian and Allen, Bruce and Sathyaprakash, B. S.
14 pages, 11 figures

This paper presents an algorithm for constructing matched-filter template banks in an arbitrary parameter space. The method places templates at random, then removes those which are “too close” together. The properties and optimality of stochastic template banks generated in this manner are investigated for some simple models. The effectiveness of these template banks for gravitational wave searches for binary inspiral waveforms is also examined. The properties of a stochastic template bank are then compared to the deterministically placed template banks that are currently used in gravitational wave data analysis.

by Klein, Antoine and Jetzer, Philippe and Sereno, Mauro
20 pages, 20 figures. Submitted to Phys. Rev. D

Gravitational waves emitted by binary systems in the inspiral phase carry a complicated structure, consisting in a superposition of different harmonics of the orbital frequency, the amplitude of each of them taking the form of a Post-Newtonian series. In addition to that, spinning binaries experience couplings which induce a precession of the orbital angular momentum and of the individual spins. So far, all studies of the measurement accuracy of gravitational wave experiments for comparable-mass binary systems have considered either spinless binaries, or spinning binaries without subdominant harmonics in the waveform as well as no amplitude modulations. In this paper, we consider supermassive black hole binaries as expected to be observed with the planned space-based interferometer LISA, and study the measurement accuracy for several astrophysically interesting parameters obtainable taking into account the full 2PN waveform for spinning bodies, as well as spin-precession effects. We find that for binaries with a total mass in the range 10^5 M_Sun < M 10 for M < 10^7 M_Sun, 1.5 – 5 times higher than with the restricted waveform. We computed that the full waveform allows to use supermassive black hole binaries as standard sirens up to a redshift of z = 1.6, about 0.4 larger than what previous studies allowed.

by Prix, Reinhard and Krishnan, Badri
12 pages, presented at GWDAW13, to appear in CQG

We investigate the Bayesian framework for detection of continuous gravitational waves (GWs) in the context of targeted searches, where the phase evolution of the GW signal is assumed to be known, while the four amplitude parameters are unknown. We show that the orthodox maximum-likelihood statistic (known as F-statistic) can be rediscovered as a Bayes factor with an unphysical prior in amplitude parameter space. We introduce an alternative detection statistic (”B-statistic”) using the Bayes factor with a more natural amplitude prior, namely an isotropic probability distribution for the orientation of GW sources. Monte-Carlo simulations of targeted searches show that the resulting Bayesian B-statistic is more powerful in the Neyman-Pearson sense (i.e. has a higher expected detection probability at equal false-alarm probability) than the frequentist F-statistic.

by Stroeer, Alexander and Veitch, John
18 pages, 8 figures, 3 tables, submitted to PRD

The Laser Interferometer Space Antenna (LISA) defines new demands on data analysis efforts in its all-sky gravitational wave survey, recording simultaneously thousands of galactic compact object binary foreground sources and tens to hundreds of background sources like binary black hole mergers and extreme mass ratio inspirals. We approach this problem with an adaptive and fully automatic Reversible Jump Markov Chain Monte Carlo sampler, able to sample from the joint posterior density function (as established by Bayes theorem) for a given mixture of signals “out of the box”, handling the total number of signals as an additional unknown parameter beside the unknown parameters of each individual source and the noise floor. We show in examples from the LISA Mock Data Challenge implementing the full response of LISA in its TDI description that this sampler is able to extract monochromatic Double White Dwarf signals out of colored instrumental noise and additional foreground and background noise successfully in a global fitting approach. We introduce 2 examples with fixed number of signals (MCMC sampling), and 1 example with unknown number of signals (RJ-MCMC), the latter further promoting the idea behind an experimental adaptation of the model indicator proposal densities in the main sampling stage. We note that the experienced runtimes and degeneracies in parameter extraction limit the shown examples to the extraction of a low but realistic number of signals.

by Reisswig, Christian and Husa, Sascha and Rezzolla, Luciano and Dorband, Ernst Nils and Pollney, Denis and Seiler, Jennifer
18 pages, 11 figures

Binary black-hole systems with spins aligned or anti-aligned to the orbital angular momentum provide the natural ground to start detailed studies of the influence of strong-field spin effects on gravitational wave observations of coalescing binaries. Furthermore, such systems may be the preferred end-state of the inspiral of generic supermassive binary black-hole systems. In view of this, we have computed the inspiral and merger of a large set of binary systems of equal-mass black holes with spins parallel to the orbital angular momentum but otherwise arbitrary. Our attention is particularly focused on the gravitational-wave emission so as to quantify how much spin effects contribute to the signal-to-noise ratio, to the horizon distances, and to the relative event rates for the representative ranges in masses and detectors. As expected, the signal-to-noise ratio increases with the projection of the total black hole spin in the direction of the orbital momentum. We find that equal-spin binaries with maximum spin aligned with the orbital angular momentum are more than “three times as loud” as the corresponding binaries with anti-aligned spins, thus corresponding to event rates up to 30 times larger. We also consider the waveform mismatch between the different spinning configurations and find that, within our numerical accuracy, binaries with opposite spins S_1=-S_2 cannot be distinguished whereas binaries with spin S_1=S_2 have clearly distinct gravitational-wave emissions. Finally, we derive a simple expression for the energy radiated in gravitational waves and find that the binaries always have efficiencies E_rad/M > 3.6%, which can become as large as E_rad/M = 10% for maximally spinning binaries with spins aligned with the orbital angular momentum.

by Lindblom, Lee
9 pages, 7 figures

Accuracy standards have been developed to ensure that the waveforms used for gravitational-wave data analysis are good enough to serve their intended purposes. These standards place constraints on certain norms of the frequency-domain representations of the waveform errors. Examples are given here of possible misinterpretations and misapplications of these standards, whose effect could be to vitiate the quality control they were intended to enforce. Suggestions are given for ways to avoid these problems.

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 Petiteau, Antoine and Yu, Shang and Babak, Stanislav
10 pages, 4 figures, proceeding for GWDAW13 (Puerto Rico)

We use a genetic algorithm to analyze the data from the third round of the mock LISA data challenge. These data consist of gaussian stationary instrumental noise, a Galactic background and four to six signals from the inspiralling spinning BHs in quasi-circular orbits. We present a particular implementation of the genetic algorithm which uses properties of the signal and the response function. We discuss the results of a preliminary search for a single signal in the instrumental noise.

by Feroz, Farhan and Gair, Jonathan R and Hobson, Michael P and Porter, Edward K
16 pages, 4 figures, submitted to Class. Quantum Grav

We describe an application of the MultiNest algorithm to gravitational wave data analysis. MultiNest is a multimodal nested sampling algorithm designed to efficiently evaluate the Bayesian evidence and return posterior probability densities for likelihood surfaces containing multiple secondary modes. The algorithm employs a set of live points which are updated by partitioning the set into multiple overlapping ellipsoids and sampling uniformly from within them. This set of live points climbs up the likelihood surface through nested iso-likelihood contours and the evidence and posterior distributions can be recovered from the point set evolution. The algorithm is model-independent in the sense that the specific problem being tackled enters only through the likelihood computation, and does not change how the live point set is updated. In this paper, we consider the use of the algorithm for gravitational wave data analysis by searching a simulated LISA data set containing two non-spinning supermassive black hole binary signals. The algorithm is able to rapidly identify all the modes of the solution and recover the true parameters of the sources to high precision.

by Sathyaprakash, B. S. and Schutz, B. F.
137 pages, 16 figures, Published version

Gravitational wave detectors are already operating at interesting sensitivity levels, and they have an upgrade path that should result in secure detections by 2014. We review the physics of gravitational waves, how they interact with detectors (bars and interferometers), and how these detectors operate. We study the most likely sources of gravitational waves and review the data analysis methods that are used to extract their signals from detector noise. Then we consider the consequences of gravitational wave detections and observations for physics, astrophysics, and cosmology.