A new code for parameter estimation in searches for gravitational waves from known pulsars

We describe the consistency testing of a new code for gravitational wave signal parameter estimation in known pulsar searches. The code uses an implementation of nested sampling to explore the likelihood volume. Using fake signals and simulated noise we compare this to a previous code that calculated the signal parameter posterior distributions on both a grid and using a crude Markov chain Monte Carlo (MCMC) method. We define a new parameterisation of two orientation angles of neutron stars used in the signal model (the initial phase and polarisation angle), which breaks a degeneracy between them and allows more efficient exploration of those parameters. Finally, we briefly describe potential areas for further study and the uses of this code in the future.Comment: Accepted for proceedings of Amaldi 9 meetin

Studying stellar binary systems with the Laser Interferometer Space Antenna using Delayed Rejection Markov chain Monte Carlo methods

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.Comment: 12 pages, 4 figures, accepted in CQG (GWDAW-13 proceedings

Compensation of Strong Thermal Lensing in High Optical Power Cavities

In an experiment to simulate the conditions in high optical power advanced gravitational wave detectors such as Advanced LIGO, we show that strong thermal lenses form in accordance with predictions and that they can be compensated using an intra-cavity compensation plate heated on its cylindrical surface. We show that high finesse ~1400 can be achieved in cavities with internal compensation plates, and that the cavity mode structure can be maintained by thermal compensation. It is also shown that the measurements allow a direct measurement of substrate optical absorption in the test mass and the compensation plate.Comment: 8 page

Supplement: Going the Distance: Mapping Host Galaxies of LIGO and Virgo Sources in Three Dimensions Using Local Cosmography and Targeted Follow-up

This is a supplement to the Letter of Singer et al. (https://arxiv.org/abs/1603.07333), in which we demonstrated a rapid algorithm for obtaining joint 3D estimates of sky location and luminosity distance from observations of binary neutron star mergers with Advanced LIGO and Virgo. We argued that combining the reconstructed volumes with positions and redshifts of possible host galaxies can provide large-aperture but small field of view instruments with a manageable list of targets to search for optical or infrared emission. In this Supplement, we document the new HEALPix-based file format for 3D localizations of gravitational-wave transients. We include Python sample code to show the reader how to perform simple manipulations of the 3D sky maps and extract ranked lists of likely host galaxies. Finally, we include mathematical details of the rapid volume reconstruction algorithm.Comment: For associated data release, see http://asd.gsfc.nasa.gov/Leo.Singer/going-the-distanc

A Bayesian parameter estimation approach to pulsar time-of-arrival analysis

The increasing sensitivities of pulsar timing arrays to ultra-low frequency (nHz) gravitational waves promises to achieve direct gravitational wave detection within the next 5-10 years. While there are many parallel efforts being made in the improvement of telescope sensitivity, the detection of stable millisecond pulsars and the improvement of the timing software, there are reasons to believe that the methods used to accurately determine the time-of-arrival (TOA) of pulses from radio pulsars can be improved upon. More specifically, the determination of the uncertainties on these TOAs, which strongly affect the ability to detect GWs through pulsar timing, may be unreliable. We propose two Bayesian methods for the generation of pulsar TOAs starting from pulsar "search-mode" data and pre-folded data. These methods are applied to simulated toy-model examples and in this initial work we focus on the issue of uncertainties in the folding period. The final results of our analysis are expressed in the form of posterior probability distributions on the signal parameters (including the TOA) from a single observation.Comment: 16 pages, 4 figure

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

Targeted search for continuous gravitational waves: Bayesian versus maximum-likelihood statistics

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.Comment: 12 pages, presented at GWDAW13, to appear in CQ

Report on the first round of the Mock LISA Data Challenges

The Mock LISA Data Challenges (MLDCs) have the dual purpose of fostering the development of LISA data analysis tools and capabilities, and demonstrating the technical readiness already achieved by the gravitational-wave community in distilling a rich science payoff from the LISA data output. The first round of MLDCs has just been completed: nine data sets containing simulated gravitational wave signals produced either by galactic binaries or massive black hole binaries embedded in simulated LISA instrumental noise were released in June 2006 with deadline for submission of results at the beginning of December 2006. Ten groups have participated in this first round of challenges. Here we describe the challenges, summarise the results, and provide a first critical assessment of the entries.Comment: Proceedings report from GWDAW 11. Added author, added reference, clarified some text, removed typos. Results unchanged; Removed author, minor edits, reflects submitted versio

Physical and biological variability in the Antarctic Polar Frontal Zone: report on research cruise 103 of the MV SA Agulhas

A detailed hydrographic and biological survey was carried out in the region of the South-west Indian Ridge during April 2002. Hydrographic data revealed that the Andrew Bain Fracture Zone, centred at 30oE, 50oS, functions as an important choke point to the flow of the Antarctic Circumpolar Current, resulting in the convergence of the Antarctic Polar Front (APF) and the southern branch of the Sub-Antarctic Front (SSAF). Total chlorophyll-a concentration and zooplankton biomass were highest at stations occupied in the vicinity of two frontal features represented by the APF and SSAF. These data suggest that the region of the South-west Indian Ridge is an area of elevated biological activity and probably acts as an important offshore feeding area for the top predators on the Prince Edward Islands