The response of interferometric gravitational wave detectors

The derivation of the response function of an interferometric gravitational wave detector is a paradigmatic calculation in the field of gravitational wave detection. Surprisingly, the standard derivation of the response wave detectors makes several unjustifiable assumptions, both conceptual and quantitative, regarding the coordinate trajectory and coordinate velocity of the null geodesic the light travels along. These errors, which appear to have remained unrecognized for at least 35 years, render the "standard" derivation inadequate and misleading as an archetype calculation. Here we identify the flaws in the existing derivation and provide, in full detail, a correct derivation of the response of a single-bounce Michelson interferometer to gravitational waves, following a procedure that will always yield correct results; compare it to the "standard", but incorrect, derivation; show where the earlier mistakes were made; and identify the general conditions under which the "standard" derivation will yield correct results. By a fortuitous set of circumstances, not generally so, the final result is the same in the case of Minkowski background spacetime, synchronous coordinates, transverse-traceless gauge metric perturbations, and arm mirrors at coordinate rest.Comment: 10 pages, one figure, as accepted to PR

Removing non-stationary, non-harmonic external interference from gravitational wave interferometer data

We describe a procedure to identify and remove a class of non-stationary and non-harmonic interference lines from gravitational wave interferometer data. These lines appear to be associated with the external electricity main supply, but their amplitudes are non-stationary and they do not appear at harmonics of the fundamental supply frequency. We find an empirical model able to represent coherently all the non-harmonic lines we have found in the power spectrum, in terms of an assumed reference signal of the primary supply input signal. If this signal is not available then it can be reconstructed from the same data by making use of the coherent line removal algorithm that we have described elsewhere. All these lines are broadened by frequency changes of the supply signal, and they corrupt significant frequency ranges of the power spectrum. The physical process that generates this interference is so far unknown, but it is highly non-linear and non-stationary. Using our model, we cancel the interference in the time domain by an adaptive procedure that should work regardless of the source of the primary interference. We have applied the method to laser interferometer data from the Glasgow prototype detector, where all the features we describe in this paper were observed. The algorithm has been tuned in such a way that the entire series of wide lines corresponding to the electrical interference are removed, leaving the spectrum clean enough to detect signals previously masked by them. Single-line signals buried in the interference can be recovered with at least 75 % of their original signal amplitude.Comment: 14 pages, 5 figures, Revtex, psfi

A new numerical method to construct binary neutron star initial data

We present a new numerical method for the generation of binary neutron star initial data using a method along the lines of the the Wilson-Mathews or the closely related conformal thin sandwich approach. Our method uses six different computational domains, which include spatial infinity. Each domain has its own coordinates which are chosen such that the star surfaces always coincide with domain boundaries. These properties facilitate the imposition of boundary conditions. Since all our fields are smooth inside each domain, we are able to use an efficient pseudospectral method to solve the elliptic equations associated with the conformal thin sandwich approach. Currently we have implemented corotating configurations with arbitrary mass ratios, but an extension to arbitrary spins is possible. The main purpose of this paper is to introduce our new method and to test our code for several different configurations.Comment: 18 pages, 8 figures, 1 tabl

Lagrangian perfect fluids and black hole mechanics

The first law of black hole mechanics (in the form derived by Wald), is expressed in terms of integrals over surfaces, at the horizon and spatial infinity, of a stationary, axisymmetric black hole, in a diffeomorphism invariant Lagrangian theory of gravity. The original statement of the first law given by Bardeen, Carter and Hawking for an Einstein-perfect fluid system contained, in addition, volume integrals of the fluid fields, over a spacelike slice stretching between these two surfaces. When applied to the Einstein-perfect fluid system, however, Wald's methods yield restricted results. The reason is that the fluid fields in the Lagrangian of a gravitating perfect fluid are typically nonstationary. We therefore first derive a first law-like relation for an arbitrary Lagrangian metric theory of gravity coupled to arbitrary Lagrangian matter fields, requiring only that the metric field be stationary. This relation includes a volume integral of matter fields over a spacelike slice between the black hole horizon and spatial infinity, and reduces to the first law originally derived by Bardeen, Carter and Hawking when the theory is general relativity coupled to a perfect fluid. We also consider a specific Lagrangian formulation for an isentropic perfect fluid given by Carter, and directly apply Wald's analysis. The resulting first law contains only surface integrals at the black hole horizon and spatial infinity, but this relation is much more restrictive in its allowed fluid configurations and perturbations than that given by Bardeen, Carter and Hawking. In the Appendix, we use the symplectic structure of the Einstein-perfect fluid system to derive a conserved current for perturbations of this system: this current reduces to one derived ab initio for this system by Chandrasekhar and Ferrari.Comment: 26 pages LaTeX-2

On the motion of spinning test particles in plane gravitational waves

The Mathisson-Papapetrou-Dixon equations for a massive spinning test particle in plane gravitational waves are analysed and explicit solutions constructed in terms of solutions of certain linear ordinary differential equations. For harmonic waves this system reduces to a single equation of Mathieu-Hill type. In this case spinning particles may exhibit parametric excitation by gravitational fields. For a spinning test particle scattered by a gravitational wave pulse, the final energy-momentum of the particle may be related to the width, height, polarisation of the wave and spin orientation of the particle.Comment: 11 page

Scattering of Spinning Test Particles by Plane Gravitational and Electromagnetic Waves

The Mathisson-Papapetrou-Dixon (MPD) equations for the motion of electrically neutral massive spinning particles are analysed, in the pole-dipole approximation, in an Einstein-Maxwell plane-wave background spacetime. By exploiting the high symmetry of such spacetimes these equations are reduced to a system of tractable ordinary differential equations. Classes of exact solutions are given, corresponding to particular initial conditions for the directions of the particle spin relative to the direction of the propagating background fields. For Einstein-Maxwell pulses a scattering cross section is defined that reduces in certain limits to those associated with the scattering of scalar and Dirac particles based on classical and quantum field theoretic techniques. The relative simplicity of the MPD approach and its use of macroscopic spin distributions suggests that it may have advantages in those astrophysical situations that involve strong classical gravitational and electromagnetic environments.Comment: Submitted to Classical and Quantum Gravity. 12 page

Computational Resources to Filter Gravitational Wave Data with P-approximant Templates

The prior knowledge of the gravitational waveform from compact binary systems makes matched filtering an attractive detection strategy. This detection method involves the filtering of the detector output with a set of theoretical waveforms or templates. One of the most important factors in this strategy is knowing how many templates are needed in order to reduce the loss of possible signals. In this study we calculate the number of templates and computational power needed for a one-step search for gravitational waves from inspiralling binary systems. We build on previous works by firstly expanding the post-Newtonian waveforms to 2.5-PN order and secondly, for the first time, calculating the number of templates needed when using P-approximant waveforms. The analysis is carried out for the four main first-generation interferometers, LIGO, GEO600, VIRGO and TAMA. As well as template number, we also calculate the computational cost of generating banks of templates for filtering GW data. We carry out the calculations for two initial conditions. In the first case we assume a minimum individual mass of $1 M_{\odot}$ and in the second, we assume a minimum individual mass of $5 M_{\odot}$. We find that, in general, we need more P-approximant templates to carry out a search than if we use standard PN templates. This increase varies according to the order of PN-approximation, but can be as high as a factor of 3 and is explained by the smaller span of the P-approximant templates as we go to higher masses. The promising outcome is that for 2-PN templates the increase is small and is outweighed by the known robustness of the 2-PN P-approximant templates.Comment: 17 pages, 8 figures, Submitted to Class.Quant.Gra

Detecting gravitational waves from test-mass bodies orbiting a Kerr black hole with P-approximant templates

In this study we apply post-Newtonian (T-approximants) and resummed post-Newtonian (P-approximants) to the case of a test-particle in equatorial orbit around a Kerr black hole. We compare the two approximants by measuring their effectualness (i.e. larger overlaps with the exact signal), and faithfulness (i.e. smaller biases while measuring the parameters of the signal) with the exact (numerical) waveforms. We find that in the case of prograde orbits, T-approximant templates obtain an effectualness of ~0.99 for spins q < 0.75. For 0.75 < q < 0.95, the effectualness drops to about 0.82. The P-approximants achieve effectualness of > 0.99 for all spins up to q = 0.95. The bias in the estimation of parameters is much lower in the case of P-approximants than T-approximants. We find that P-approximants are both effectual and faithful and should be more effective than T-approximants as a detection template family when q > 0. For q < 0 both T- and P-approximants perform equally well so that either of them could be used as a detection template family. However, for parameter estimation, the P-approximant templates still outperforms the T-approximants.Comment: 11 Pages - 9 figures. Accepted for publication. Proceedings of GWDAW 9. Special edition of Classical and Quantum Gravit

On the Solution to the "Frozen Star" Paradox, Nature of Astrophysical Black Holes, non-Existence of Gravitational Singularity in the Physical Universe and Applicability of the Birkhoff's Theorem

Oppenheimer and Snyder found in 1939 that gravitational collapse in vacuum produces a "frozen star", i.e., the collapsing matter only asymptotically approaches the gravitational radius (event horizon) of the mass, but never crosses it within a finite time for an external observer. Based upon our recent publication on the problem of gravitational collapse in the physical universe for an external observer, the following results are reported here: (1) Matter can indeed fall across the event horizon within a finite time and thus BHs, rather than "frozen stars", are formed in gravitational collapse in the physical universe. (2) Matter fallen into an astrophysical black hole can never arrive at the exact center; the exact interior distribution of matter depends upon the history of the collapse process. Therefore gravitational singularity does not exist in the physical universe. (3) The metric at any radius is determined by the global distribution of matter, i.e., not only by the matter inside the given radius, even in a spherically symmetric and pressureless gravitational system. This is qualitatively different from the Newtonian gravity and the common (mis)understanding of the Birkhoff's Theorem. This result does not contract the "Lemaitre-Tolman-Bondi" solution for an external observer.Comment: 8 pages, 4 figures, invited plenary talk at "The first Galileo-Xu Guangqi conference", Shanghai, China, 2009. To appear in International Journal of Modern Physics D (2010

Gravitational wave astronomy

The first decade of the new millenium should see the first direct detections of gravitational waves. This will be a milestone for fundamental physics and it will open the new observational science of gravitational wave astronomy. But gravitational waves already play an important role in the modeling of astrophysical systems. I review here the present state of gravitational radiation theory in relativity and astrophysics, and I then look at the development of detector sensitivity over the next decade, both on the ground (such as LIGO) and in space (LISA). I review the sources of gravitational waves that are likely to play an important role in observations by first- and second-generation interferometers, including the astrophysical information that will come from these observations. The review covers some 10 decades of gravitational wave frequency, from the high-frequency normal modes of neutron stars down to the lowest frequencies observable from space. The discussion of sources includes recent developments regarding binary black holes, spinning neutron stars, and the stochastic background.Comment: 29 pages, 2 figures, as submitted for special millenium issue of Classical and Quantum Gravit