R. Vaulin

Testing gravitational-wave searches with numerical relativity waveforms: results from the first Numerical INJection Analysis (NINJA) project

The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave data analysis communities. The purpose of NINJA is to study the sensitivity of existing gravitational-wave search algorithms using numerically generated waveforms and to foster closer collaboration between the numerical relativity and data analysis communities. We describe the results of the first NINJA analysis which focused on gravitational waveforms from binary black hole coalescence. Ten numerical relativity groups contributed numerical data which were used to generate a set of gravitational-wave signals. These signals were injected into a simulated data set, designed to mimic the response of the initial LIGO and Virgo gravitational-wave detectors. Nine groups analysed this data using search and parameter-estimation pipelines. Matched filter algorithms, un-modelled-burst searches and Bayesian parameter estimation and model-selection algorithms were applied to the data. We report the efficiency of these search methods in detecting the numerical waveforms and measuring their parameters. We describe preliminary comparisons between the different search methods and suggest improvements for future NINJA analyses.

Macroscopic Effects of the Quantum Trace Anomaly

The low energy effective action of gravity in any even dimension generally acquires nonlocal terms associated with the trace anomaly, generated by the quantum fluctuations of massless fields. The local auxiliary field description of this effective action in four dimensions requires two additional scalar fields, not contained in classical general relativity, which remain relevant at macroscopic distance scales. The auxiliary scalar fields depend upon boundary conditions for their complete specification, and therefore carry global information about the geometry and macroscopic quantum state of the gravitational field. The scalar potentials also provide coordinate invariant order parameters describing the conformal behavior and divergences of the stress tensor on event horizons. We compute the stress tensor due to the anomaly in terms of its auxiliary scalar potentials in a number of concrete examples, including the Rindler wedge, the Schwarzschild geometry, and de Sitter spacetime. In all of these cases, a small number of classical order parameters completely determine the divergent behaviors allowed on the horizon, and yield qualitatively correct global approximations to the renormalized expectation value of the quantum stress tensor.

Status of NINJA: the Numerical INJection Analysis project

The 2008 NRDA conference introduced the Numerical INJection Analysis project (NINJA), a new collaborative effort between the numerical relativity community and the data analysis community. NINJA focuses on modeling and searching for gravitational wave signatures from the coalescence of binary system of compact objects. We review the scope of this collaboration and the components of the first NINJA project, where numerical relativity groups, shared waveforms and data analysis teams applied various techniques to detect them when embedded in colored Gaussian noise.

Sensitivity Comparison of Searches for Binary Black Hole Coalescences with Ground-based Gravitational-Wave Detectors

Searches for gravitational-wave transients from binary black hole coalescences typically rely on one of two approaches: matched filtering with templates and morphology-independent excess power searches. Multiple algorithmic implementations in the analysis of data from the first generation of ground-based gravitational wave interferometers have used different strategies for the suppression of non-Gaussian noise transients, and targeted different regions of the binary black hole parameter space. In this paper we compare the sensitivity of three such algorithms: matched filtering with full coalescence templates, matched filtering with ringdown templates and a morphology-independent excess power search. The comparison is performed at a fixed false alarm rate and relies on Monte-carlo simulations of binary black hole coalescences for spinning, non-precessing systems with total mass 25-350 solar mass, which covers the parameter space of stellar mass and intermediate mass black hole binaries. We find that in the mass range of 25 -100 solar mass the sensitive distance of the search, marginalized over source parameters, is best with matched filtering to full waveform templates, to within 10 percent at a false alarm rate of 3 events per year. In the mass range of 100-350 solar mass, the same comparison favors the morphology-independent excess power search to within 20 percent. The dependence on mass and spin is also explored.

Likelihood-ratio ranking of gravitational-wave candidates in a non-Gaussian background.

We describe a general approach to detection of transient gravitational-wave signals in the presence of non-Gaussian background noise. We prove that under quite general conditions, the ratio of the likelihood of observed data to contain a signal to the likelihood of it being a noise fluctuation provides optimal ranking for the candidate events found in an experiment. The likelihood-ratio ranking allows us to combine different kinds of data into a single analysis. We apply the general framework to the problem of unifying the results of independent experiments and the problem of accounting for non-Gaussian artifacts in the searches for gravitational waves from compact binary coalescence in LIGO data. We show analytically and confirm through simulations that in both cases applying the likelihood-ratio ranking results in an improved analysis.

Quantum geometrodynamics of the Bianchi IX cosmological model

The canonical quantum theory of gravity?quantum geometrodynamics (QG)?is applied to the homogeneous Bianchi type IX cosmological model. As a result, a framework for the quantum theory of homogeneous cosmologies is developed. We show that the theory is internally consistent and prove that it possesses the correct classical limit (the theory of general relativity). To emphasize the special role that the constraints play in this new theory, we compare it to the traditional ADM square-root and Wheeler?DeWitt quantization schemes. We show that, unlike traditional approaches, QG leads to a well-defined Schr?dinger equation for the wavefunction of the universe that is inherently coupled to the expectation value of the constraint equations. This coupling to the constraints is responsible for the appearance of a coherent spacetime picture. Thus, the physical meaning of the constraints of the theory is quite different from Diracs interpretation. In light of this distinctive feature of the theory, we re-address the question of the dark energy effects in the Bianchi IX cosmological model for highly non-classical quantum states. We show that, at least for this model, for any choice of the initial wavefunction, quantum corrections will not produce accelerated expansion of the universe.

Detecting transient gravitational waves in non-Gaussian noise with partially redundant analysis methods

There is a broad class of astrophysical sources that produce detectable, transient, gravitational waves. Some searches for transient gravitational waves are tailored to known features of these sources. Other searches make few assumptions about the sources. Typically events are observable with multiple search techniques. This work describes how to combine the results of searches that are not independent, treating each search as a classifier for a given event. This will be shown to improve the overall sensitivity to gravitational-wave events while directly addressing the problem of consistent interpretation of multiple trials.

EFFECTS OF QUANTIZED FIELDS ON THE SPACETIME GEOMETRIES OF STATIC SPHERICALLY SYMMETRIC BLACK HOLES

Analytic approximations for the stress-energy of quantized fields in the Hartle-Hawking state in static black hole spacetimes predict divergences on the event horizon of the black hole for a number of important cases. Such divergences, if real, could substantially alter the spacetime geometry near the event horizon, possibly preventing the black hole from existing. The results of three investigations of these types of effects are presented. The first involves a new analytic approximation for conformally invariant fields in Reissner-Nordstrom spacetimes which is finite on the horizon. The second focuses on the stress-energy of massless scalar fields in Schwarzschild-de Sitter black holes. The third focuses on the stress-energy of massless scalar fields in zero temperature black hole geometries that could be solutions to the semiclassical backreaction equations near the event horizon of the black hole.

A Matched Filter Technique for Slow Radio Transient Detection and First Demonstration with the Murchison Widefield Array

Many astronomical sources produce transient phenomena at radio frequencies, but the transient sky at low frequencies (<300 MHz) remains relatively unexplored. Blind surveys with new widefield radio instruments are setting increasingly stringent limits on the transient surface density on various timescales. Although many of these instruments are limited by classical confusion noise from an ensemble of faint, unresolved sources, one can in principle detect transients below the classical confusion limit to the extent that the classical confusion noise is independent of time. We develop a technique for detecting radio transients that is based on temporal matched filters applied directly to time series of images rather than relying on source-finding algorithms applied to individual images. This technique has well-defined statistical properties and is applicable to variable and transient searches for both confusion-limited and non-confusion-limited instruments. Using the Murchison Widefield Array as an example, we demonstrate that the technique works well on real data despite the presence of classical confusion noise, sidelobe confusion noise, and other systematic errors. We searched for transients lasting between 2 minutes and 3 months. We found no transients and set improved upper limits on the transient surface density at 182 MHz for flux densities between ~20--200 mJy, providing the best limits to date for hour- and month-long transients.