A principal component approach to space-based gravitational wave astronomy

The current approach to data analysis for the Laser Interferometry Space Antenna (LISA) depends on the time delay interferometry observables (TDI) which have to be generated before any weak signal detection can be performed. These are linear combinations of the raw data with appropriate time shifts that lead to the cancellation of the laser frequency noises. This is possible because of the multiple occurrences of the same noises in the different raw data. Originally, these observables were manually generated starting with LISA as a simple stationary array and then adjusted to incorporate the antenna's motions. However, none of the observables survived the flexing of the arms in that they did not lead to cancellation with the same structure. The principal component approach is another way of handling these noises that was presented by Romano and Woan which simplified the data analysis by removing the need to create them before the analysis. This method also depends on the multiple occurrences of the same noises but, instead of using them for cancellation, it takes advantage of the correlations that they produce between the different readings. These correlations can be expressed in a noise (data) covariance matrix which occurs in the Bayesian likelihood function when the noises are assumed be Gaussian. Romano and Woan showed that performing an eigendecomposition of this matrix produced two distinct sets of eigenvalues that can be distinguished by the absence of laser frequency noise from one set. The transformation of the raw data using the corresponding eigenvectors also produced data that was free from the laser frequency noises. This result led to the idea that the principal components may actually be time delay interferometry observables since they produced the same outcome, that is, data that are free from laser frequency noise. The aims here were (i) to investigate the connection between the principal components and these observables, (ii) to prove that the data analysis using them is equivalent to that using the traditional observables and (ii) to determine how this method adapts to real LISA especially the flexing of the antenna. For testing the connection between the principal components and the TDI observables a 10x 10 covariance matrix containing integer values was used in order to obtain an algebraic solution for the eigendecomposition. The matrix was generated using fixed unequal arm lengths and stationary noises with equal variances for each noise type. Results confirm that all four Sagnac observables can be generated from the eigenvectors of the principal components. The observables obtained from this method however, are tied to the length of the data and are not general expressions like the traditional observables, for example, the Sagnac observables for two different time stamps were generated from different sets of eigenvectors. It was also possible to generate the frequency domain optimal AET observables from the principal components obtained from the power spectral density matrix. These results indicate that this method is another way of producing the observables therefore analysis using principal components should give the same results as that using the traditional observables. This was proven by fact that the same relative likelihoods (within 0.3%) were obtained from the Bayesian estimates of the signal amplitude of a simple sinusoidal gravitational wave using the principal components and the optimal AET observables. This method fails if the eigenvalues that are free from laser frequency noises are not generated. These are obtained from the covariance matrix and the properties of LISA that are required for its computation are the phase-locking, arm lengths and noise variances. Preliminary results of the effects of these properties on the principal components indicate that only the absence of phase-locking prevented their production. The flexing of the antenna results in time varying arm lengths which will appear in the covariance matrix and, from our toy model investigations, this did not prevent the occurrence of the principal components. The difficulty with flexing, and also non-stationary noises, is that the Toeplitz structure of the matrix will be destroyed which will affect any computation methods that take advantage of this structure. In terms of separating the two sets of data for the analysis, this was not necessary because the laser frequency noises are very large compared to the photodetector noises which resulted in a significant reduction in the data containing them after the matrix inversion. In the frequency domain the power spectral density matrices were block diagonals which simplified the computation of the eigenvalues by allowing them to be done separately for each block. The results in general showed a lack of principal components in the absence of phase-locking except for the zero bin. The major difference with the power spectral density matrix is that the time varying arm lengths and non-stationarity do not show up because of the summation in the Fourier transform.

Gravitational wave recoil in Robinson-Trautman spacetimes

We consider the gravitational recoil due to nonreflection-symmetric gravitational wave emission in the context of axisymmetric Robinson-Trautman spacetimes. We show that regular initial data evolve generically into a final configuration corresponding to a Schwarzschild black hole moving with constant speed. For the case of (reflection-)symmetric initial configurations, the mass of the remnant black hole and the total energy radiated away are completely determined by the initial data, allowing us to obtain analytical expressions for some recent numerical results that have appeared in the literature. Moreover, by using the Galerkin spectral method to analyze the nonlinear regime of the Robinson-Trautman equations, we show that the recoil velocity can be estimated with good accuracy from some asymmetry measures (namely the first odd moments) of the initial data. The extension for the nonaxisymmetric case and the implications of our results for realistic situations involving head-on collision of two black holes are also discussed.

The Schenberg spherical gravitational wave detector: the first commissioning runs

Here we present a status report of the first spherical antenna project equipped with a set of parametric transducers for gravitational detection. The Mario Schenberg, as it is called, started its commissioning phase at the Physics Institute of the University of Sao Paulo, in September 2006, under the full support of FAPESP. We have been testing the three preliminary parametric transducer systems in order to prepare the detector for the next cryogenic run, when it will be calibrated. We are also developing sapphire oscillators that will replace the current ones thereby providing better performance. We also plan to install eight transducers in the near future, six of which are of the two-mode type and arranged according to the truncated icosahedron configuration. The other two, which will be placed close to the sphere equator, will be mechanically non-resonant. In doing so, we want to verify that if the Schenberg antenna can become a wideband gravitational wave detector through the use of an ultra-high sensitivity non-resonant transducer constructed using the recent achievements of nanotechnology.

Gravitational wave pulse and soliton wave collision

We study the collision of a gravitational wave pulse and a soliton wave on a spatially homogeneous background. This collision is described by an exact solution of Einsteins equations in a vacuum which is generated from a nondiagonal seed by means of a soliton transformation. The effect produced by the soliton on the amplitude and polarization of the wave is considered.

Gravitational wave pulse and soliton wave collision

We study the collision of a gravitational wave pulse and a soliton wave on a spatially homogeneous background. This collision is described by an exact solution of Einsteins equations in a vacuum which is generated from a nondiagonal seed by means of a soliton transformation. The effect produced by the soliton on the amplitude and polarization of the wave is considered.

What can we learn about gravitational wave physics with an elastic spherial antenna?

A general formalism is set up to analyze the response of an arbitrary solid elastic body to an arbitrary metric gravitational wave (GW) perturbation, which fully displays the details of the interaction antenna wave. The formalism is applied to the spherical detector, whose sensitivity parameters are thereby scrutinized. A multimode transfer function is defined to study the amplitude sensitivity, and absorption cross sections are calculated for a general metric theory of GW physics. Their scaling properties are shown to be independent of the underlying theory, with interesting consequences for future detector design. The GW incidence direction deconvolution problem is also discussed, always within the context of a general metric theory of the gravitational field.

Proposed gravitational wave observatory based on solid elastic spheres

Spherical gravitational wave (GW) detectors offer a wealth of so far unexplored possibilities to detect gravitational radiation. We find that a sphere can be used as a powerful testbed for any metric theory of gravity, not only general relativity as considered so far, by making use of a deconvolution procedure for all the electric components of the Riemann tensor. We also find that the spheres cross section is large at two frequencies, and advantageous at higher frequencies in the sense that a single antenna constitutes a real xylophone in its own. Proposed GW networks will greatly benefit from this. The main features of a two large sphere observatory are reported.

Thermal and back-action noises in dual-sphere gravitational-wave detectors

We study the sensitivity limits of a broadband gravitational-wave detector based on dual resonators such as nested spheres. We determine both the thermal and back-action noises when the resonators displacements are read out with an optomechanical sensor. We analyze the contributions of all mechanical modes, using a new method to deal with the force-displacement transfer functions in the intermediate frequency domain between the two gravitational-wave sensitive modes associated with each resonator. This method gives an accurate estimate of the mechanical response, together with an evaluation of the estimate error. We show that very high sensitivities can be reached on a wide frequency band for realistic parameters in the case of a dual-sphere detector.

All-sky search algorithms for monochromatic signals in bar gravitational wave detector data

In this paper we design and develop several filtering strategies for the analysis of data generated by a resonant bar gravitational wave (GW) antenna, with the goal of assessing the presence (or absence) therein of long-duration monochromatic GW signals, as well as the eventual amplitude and frequency of the signals, within the sensitivity band of the detector. Such signals are most likely generated in the fast rotation of slightly asymmetric spinning stars. We develop practical procedures, together with a study of their statistical properties, which will provide us with useful information on the performance of each technique. The selection of candidate events will then be established according to threshold-crossing probabilities, based on the Neyman-Pearson criterion. In particular, it will be shown that our approach, based on phase estimation, presents a better signal-to-noise ratio than does pure spectral analysis, the most common approach.

On the energy transported by exact plane gravitational-wave solutions

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Search for gravitational wave ringdowns from perturbed intermediate mass black holes in LIGO-Virgo data from 2005-2010

We report results from a search for gravitational waves produced by perturbed intermediate mass black holes ( IMBH) in data collected by LIGO and Virgo between 2005 and 2010. The search was sensitive to astrophysical sources that produced damped sinusoid gravitational wave signals, also known as ringdowns, with frequency 50 <= f(0)/Hz <= 2000 and decay timescale 0.0001 less than or similar to tau/s less than or similar to 0.1 characteristic of those produced in mergers of IMBH pairs. No significant gravitational wave candidate was detected. We report upper limits on the astrophysical coalescence rates of IMBHs with total binary mass 50 <= M/ M circle dot <= 450 and component mass ratios of either 1: 1 or 4: 1. For systems with total mass 100 <= M/M circle dot <= 150, we report a 90% confidence upper limit on the rate of binary IMBH mergers with nonspinning and equal mass components of 6.9 x 10(-8) Mpc(-3) yr(-1). We also report a rate upper limit for ringdown waveforms from perturbed IMBHs, radiating 1% of their mass as gravitational waves in the fundamental, l = m = 2, oscillation mode, that is nearly three orders of magnitude more stringent than previous results.

The NINJA-2 project: detecting and characterizing gravitational waveforms modelled using numerical binary black hole simulations

The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave (GW) astrophysics communities. The purpose of NINJA is to study the ability to detect GWs emitted from merging binary black holes (BBH) and recover their parameters with next-generation GW observatories. We report here on the results of the second NINJA project, NINJA-2, which employs 60 complete BBH hybrid waveforms consisting of a numerical portion modelling the late inspiral, merger, and ringdown stitched to a post-Newtonian portion modelling the early inspiral. In a 'blind injection challenge' similar to that conducted in recent Laser Interferometer Gravitational Wave Observatory (LIGO) and Virgo science runs, we added seven hybrid waveforms to two months of data recoloured to predictions of Advanced LIGO (aLIGO) and Advanced Virgo (AdV) sensitivity curves during their first observing runs. The resulting data was analysed by GW detection algorithms and 6 of the waveforms were recovered with false alarm rates smaller than 1 in a thousand years. Parameter-estimation algorithms were run on each of these waveforms to explore the ability to constrain the masses, component angular momenta and sky position of these waveforms. We find that the strong degeneracy between the mass ratio and the BHs' angular momenta will make it difficult to precisely estimate these parameters with aLIGO and AdV. We also perform a large-scale Monte Carlo study to assess the ability to recover each of the 60 hybrid waveforms with early aLIGO and AdV sensitivity curves. Our results predict that early aLIGO and AdV will have a volume-weighted average sensitive distance of 300 Mpc (1 Gpc) for 10M circle dot + 10M circle dot (50M circle dot + 50M circle dot) BBH coalescences. We demonstrate that neglecting the component angular momenta in the waveform models used in matched-filtering will result in a reduction in sensitivity for systems with large component angular momenta. This reduction is estimated to be up to similar to 15% for 50M circle dot + 50M circle dot BBH coalescences with almost maximal angular momenta aligned with the orbit when using early aLIGO and AdV sensitivity curves.

Reconstruction of the gravitational wave signal h(t) during the Virgo science runs and independent validation with a photon calibrator

The Virgo detector is a kilometer-scale interferometer for gravitational wave detection located near Pisa (Italy). About 13 months of data were accumulated during four science runs (VSR1, VSR2, VSR3 and VSR4) between May 2007 and September 2011, with increasing sensitivity. In this paper, the method used to reconstruct, in the range 10 Hz-10 kHz, the gravitational wave strain time series h(t) from the detector signals is described. The standard consistency checks of the reconstruction are discussed and used to estimate the systematic uncertainties of the h(t) signal as a function of frequency. Finally, an independent setup, the photon calibrator, is described and used to validate the reconstructed h(t) signal and the associated uncertainties. The systematic uncertainties of the h(t) time series are estimated to be 8% in amplitude. The uncertainty of the phase of h(t) is 50 mrad at 10 Hz with a frequency dependence following a delay of 8 mu s at high frequency. A bias lower than 4 mu s and depending on the sky direction of the GW is also present.