The information content of gravitational wave harmonics in compact binary inspiral

The nonlinear aspect of gravitational wave generation that produces power at harmonics of the orbital frequency, above the fundamental quadrupole frequency, is examined to see what information about the source is contained in these higher harmonics. We use an order (4/2) post-Newtonian expansion of the gravitational wave waveform of a binary system to model the signal seen in a spaceborne gravitational wave detector such as the proposed LISA detector. Covariance studies are then performed to determine the ultimate accuracy to be expected when the parameters of the source are fit to the received signal. We find three areas where the higher harmonics contribute crucial information that breaks degeneracies in the model and allows otherwise badly-correlated parameters to be separated and determined. First, we find that the position of a coalescing massive black hole binary in an ecliptic plane detector, such as OMEGA, is well-determined with the help of these harmonics. Second, we find that the individual masses of the stars in a chirping neutron star binary can be separated because of the mass dependence of the harmonic contributions to the wave. Finally, we note that supermassive black hole binaries, whose frequencies are too low to be seen in the detector sensitivity window for long, may still have their masses, distances, and positions determined since the information content of the higher harmonics compensates for the information lost when the orbit-induced modulation of the signal does not last long enough to be apparent in the data.Comment: 13 pages, 5 figure

Periodic harmonic functions on lattices and points count in positive characteristic

This survey addresses pluri-periodic harmonic functions on lattices with values in a positive characteristic field. We mention, as a motivation, the game "Lights Out" following the work of Sutner, Goldwasser-Klostermeyer-Ware, Barua-Ramakrishnan-Sarkar, Hunzikel-Machiavello-Park e.a.; see also 2 previous author's preprints for a more detailed account. Our approach explores harmonic analysis and algebraic geometry over a positive characteristic field. The Fourier transform allows us to interpret pluri-periods of harmonic functions on lattices as torsion multi-orders of points on the corresponding affine algebraic variety.Comment: These are notes on 13p. based on a talk presented during the meeting "Analysis on Graphs and Fractals", the Cardiff University, 29 May-2 June 2007 (a sattelite meeting of the programme "Analysis on Graphs and its Applications" at the Isaac Newton Institute from 8 January to 29 June 2007

An explicit harmonic code for black-hole evolution using excision

We describe an explicit in time, finite-difference code designed to simulate black holes by using the excision method. The code is based upon the harmonic formulation of the Einstein equations and incorporates several features regarding the well-posedness and numerical stability of the initial-boundary problem for the quasilinear wave equation. After a discussion of the equations solved and of the techniques employed, we present a series of testbeds carried out to validate the code. Such tests range from the evolution of isolated black holes to the head-on collision of two black holes and then to a binary black hole inspiral and merger. Besides assessing the accuracy of the code, the inspiral and merger test has revealed that individual apparent horizons can touch and even intersect. This novel feature in the dynamics of the marginally trapped surfaces is unexpected but consistent with theorems on the properties of apparent horizons

Solving Einstein's Equations With Dual Coordinate Frames

A method is introduced for solving Einstein's equations using two distinct coordinate systems. The coordinate basis vectors associated with one system are used to project out components of the metric and other fields, in analogy with the way fields are projected onto an orthonormal tetrad basis. These field components are then determined as functions of a second independent coordinate system. The transformation to the second coordinate system can be thought of as a mapping from the original ``inertial'' coordinate system to the computational domain. This dual-coordinate method is used to perform stable numerical evolutions of a black-hole spacetime using the generalized harmonic form of Einstein's equations in coordinates that rotate with respect to the inertial frame at infinity; such evolutions are found to be generically unstable using a single rotating coordinate frame. The dual-coordinate method is also used here to evolve binary black-hole spacetimes for several orbits. The great flexibility of this method allows comoving coordinates to be adjusted with a feedback control system that keeps the excision boundaries of the holes within their respective apparent horizons.Comment: Updated to agree with published versio

Implementation of higher-order absorbing boundary conditions for the Einstein equations

We present an implementation of absorbing boundary conditions for the Einstein equations based on the recent work of Buchman and Sarbach. In this paper, we assume that spacetime may be linearized about Minkowski space close to the outer boundary, which is taken to be a coordinate sphere. We reformulate the boundary conditions as conditions on the gauge-invariant Regge-Wheeler-Zerilli scalars. Higher-order radial derivatives are eliminated by rewriting the boundary conditions as a system of ODEs for a set of auxiliary variables intrinsic to the boundary. From these we construct boundary data for a set of well-posed constraint-preserving boundary conditions for the Einstein equations in a first-order generalized harmonic formulation. This construction has direct applications to outer boundary conditions in simulations of isolated systems (e.g., binary black holes) as well as to the problem of Cauchy-perturbative matching. As a test problem for our numerical implementation, we consider linearized multipolar gravitational waves in TT gauge, with angular momentum numbers l=2 (Teukolsky waves), 3 and 4. We demonstrate that the perfectly absorbing boundary condition B_L of order L=l yields no spurious reflections to linear order in perturbation theory. This is in contrast to the lower-order absorbing boundary conditions B_L with L<l, which include the widely used freezing-Psi_0 boundary condition that imposes the vanishing of the Newman-Penrose scalar Psi_0.Comment: 25 pages, 9 figures. Minor clarifications. Final version to appear in Class. Quantum Grav

Simulation of Binary Black Hole Spacetimes with a Harmonic Evolution Scheme

A numerical solution scheme for the Einstein field equations based on generalized harmonic coordinates is described, focusing on details not provided before in the literature and that are of particular relevance to the binary black hole problem. This includes demonstrations of the effectiveness of constraint damping, and how the time slicing can be controlled through the use of a source function evolution equation. In addition, some results from an ongoing study of binary black hole coalescence, where the black holes are formed via scalar field collapse, are shown. Scalar fields offer a convenient route to exploring certain aspects of black hole interactions, and one interesting, though tentative suggestion from this early study is that behavior reminiscent of "zoom-whirl" orbits in particle trajectories is also present in the merger of equal mass, non-spinning binaries, with appropriately fine-tuned initial conditions.Comment: 16 pages, 14 figures; replaced with published versio