New Numerical Methods to Evaluate Homogeneous Solutions of the Teukolsky Equation

We discuss a numerical method to compute the homogeneous solutions of the Teukolsky equation which is the basic equation of the black hole perturbation method. We use the formalism developed by Mano, Suzuki and Takasugi, in which the homogeneous solutions of the radial Teukolsky equation are expressed in terms of two kinds of series of special functions, and the formulas for the asymptotic amplitudes are derived explicitly.Although the application of this method was previously limited to the analytical evaluation of the homogeneous solutions, we find that it is also useful for numerical computation. We also find that so-called "renormalized angular momentum parameter", $\nu$, can be found only in the limited region of $\omega$ for each $l,m$ if we assume $\nu$ is real (here, $\omega$ is the angular frequency, and $l$ and $m$ are degree and order of the spin-weighted spheroidal harmonics respectively). We also compute the flux of the gravitational waves induced by a compact star in a circular orbit on the equatorial plane around a rotating black hole. We find that the relative error of the energy flux is about $10^{-14}$ which is much smaller than the one obtained by usual numerical integration methods.Comment: 36 pages,7 figure

An Improved Search Method for Gravitational Ringing of Black Holes

A black hole has characteristic quasi-normal modes that will be excited when it is formed or when the geometry is perturbed. The state of a black hole when the quasi-normal modes are excited is called the gravitational ringing, and detections of it will be a direct confirmation of the existence of black holes. To detect it, a method based on matched filtering needs to be developed. Generically, matched filtering requires a large number of templates, because one has to ensure a proper match of a real gravitational wave with one of template waveforms to keep the detection efficiency as high as possible. On the other hand, the number of templates must be kept as small as possible under limited computational costs. In our previous paper, assuming that the gravitational ringing is dominated by the least-damped (fundamental) mode with the least imaginary part of frequency, we constructed an efficient method for tiling the template space. However, the dependence of the template space metric on the initial phase of a wave was not taken into account. This dependence arises because of an unavoidable mismatch between the parameters of a signal waveform and those given discretely in the template space. In this paper, we properly take this dependence into account and present an improved, efficient search method for gravitational ringing of black holes.Comment: 19 pages, 9 figure

Coherent versus coincidence detection of gravitational wave signals from compact inspiraling binaries

We compare two multi-detector detection strategies, namely, the coincidence and the coherent, for the detection of spinless inspiraling compact binary gravitational wave signals. The coincident strategy treats the detectors as if they are isolated - compares individual detector statistics with their respective thresholds while the coherent strategy combines the detector network data {\it phase coherently} to obtain a single detection statistic which is then compared with a single threshold. In the case of geographically separated detectors, we also consider an {\it enhanced} coincidence strategy because the usual (naive) coincidence strategy yields poor results for misaligned detectors. For simplicity, we consider detector pairs having the same power spectral density of noise, as that of initial LIGO and also assume the noise to be stationary and Gaussian. We compare the performances of the methods by plotting the \emph{receiver operating characteristic} (ROC) for the two strategies. A single astrophysical source as well as a distribution of sources is considered. We find that the coherent strategy performs better than the two coincident strategies under the assumptions of stationary Gaussian detector noise.Comment: Based on the presentation at the 1st Galileo Xu Guangqi conference, Shanghai

Detecting gravitational waves from precessing binaries of spinning compact objects. II. Search implementation for low-mass binaries

Detection template families (DTFs) are built to capture the essential features of true gravitational waveforms using a small set of phenomenological waveform parameters. Buonanno, Chen, and Vallisneri [Phys. Rev. D 67, 104025 (2003)] proposed the ``BCV2'' DTF to perform computationally efficient searches for signals from precessing binaries of compact stellar objects. Here we test the signal-matching performance of the BCV2 DTF for asymmetric--mass-ratio binaries, and specifically for double--black-hole binaries with component masses (m1,m2): (6~12Msun, 1~3Msun), and for black-hole--neutron-star binaries with component masses (m1,m2) = (10Msun, 1.4Msun); we take all black holes to be maximally spinning. We find a satisfactory signal-matching performance, with fitting factors averaging between 0.94 and 0.98. We also scope out the region of BCV2 parameters needed for a template-based search, we evaluate the template match metric, we discuss a template-placement strategy, and we estimate the number of templates needed for searches at the LIGO design sensitivity. In addition, after gaining more insight in the dynamics of spin--orbit precession, we propose a modification of the BCV2 DTF that is parametrized by physical (rather than phenomenological) parameters. We test this modified ``BCV2P'' DTF for the (10Msun, 1.4Msun) black-hole--neutron-star system, finding a signal-matching performance comparable to the BCV2 DTF, and a reliable parameter-estimation capability for target-binary quantities such as the chirp mass and the opening angle (the angle between the black-hole spin and the orbital angular momentum).Comment: 18 pages, 15 figure

Multipole particle in relativity

We discuss the motion of extended objects in a spacetime by considering a gravitational field created by these objects. We define multipole moments of the objects as a classification by Lie group SO(3). Then, we construct an energy-momentum tensor for the objects and derive equations of motion from it. As a result, we reproduce the Papapetrou equations for a spinning particle. Furthermore, we will show that we can obtain more simple equations than the Papapetrou equations by changing the center-of-mass.Comment: 22 pages, 2 figures. Accepted for publication in Phys. Rev.

Gravitational radiation from infall into a black hole: Regularization of the Teukolsky equation

The Teukolsky equation has long been known to lead to divergent integrals when it is used to calculate the gravitational radiation emitted when a test mass falls into a black hole from infinity. Two methods have been used in the past to remove those divergent integrals. In the first, integrations by parts are carried out, and the infinite boundary terms are simply discarded. In the second, the Teukolsky equation is transformed into another equation which does not lead to divergent integrals. The purpose of this paper is to show that there is nothing intrinsically wrong with the Teukolsky equation when dealing with non-compact source terms, and that the divergent integrals result simply from an incorrect choice of Green's function. In this paper, regularization of the Teukolsky equation is carried out in an entirely natural way which does not involve modifying the equation.Comment: ReVTeX, 23 page

Post-Newtonian expansion of gravitational waves from a particle in circular orbits around a rotating black hole: Up to O(v8)O(v^8) beyond the quadrupole formula

Extending a method developed by Sasaki in the Schwarzschild case and by Shibata, Sasaki, Tagoshi, and Tanaka in the Kerr case, we calculate the post-Newtonian expansion of the gravitational wave luminosities from a test particle in circular orbit around a rotating black hole up to $O(v^8)$ beyond the quadrupole formula. The orbit of a test particle is restricted on the equatorial plane. We find that spin dependent terms appear in each post-Newtonian order, and that at $O(v^6)$ they have a significant effect on the orbital phase evolution of coalescing compact binaries. By comparing the post-Newtonian formula of the luminosity with numerical results we find that, for $30M\lesssim r \lesssim 100M$, the spin dependent terms at $O(v^6)$ and $O(v^7)$ improve the accuracy of the post-Newtonian formula significantly, but those at $O(v^8)$ do not improve.Comment: 27 pages, revtex, 6 figures, submitted to Physical Review

Chi-square test on candidate events from CW signal coherent searches

In a blind search for continuous gravitational wave signals scanning a wide frequency band one looks for candidate events with significantly large values of the detection statistic. Unfortunately, a noise line in the data may also produce a moderately large detection statistic. In this paper, we describe how we can distinguish between noise line events and actual continuous wave (CW) signals, based on the shape of the detection statistic as a function of the signal's frequency. We will analyze the case of a particular detection statistic, the F statistic, proposed by Jaranowski, Krolak, and Schutz. We will show that for a broad-band 10 hour search, with a false dismissal rate smaller than 1e-6, our method rejects about 70 % of the large candidate events found in a typical data set from the second science run of the Hanford LIGO interferometer.Comment: proceedings of GWDAW8, 2003 conference, 12pages, 6 figure

Gravitational waves from inspiralling compact binaries with magnetic dipole moments

We investigate the effects of the magnetic dipole-dipole coupling and the electromagnetic radiation on the frequency evolution of gravitational waves from inspiralling binary neutron stars with magnetic dipole moments. This study is motivated by the discovery of the superstrongly magnetized neutron stars, i.e., magnetar. We derive the contributions of the magnetic fields to the accumulated cycles in gravitational waves as $N_{mag} \sim 6 \times 10^{-3} (H/10^{16}{\rm G})^{2}$, where $H$ denotes the strength of the polar magnetic fields of each neutron star in the binary system. It is found that the effects of the magnetic fields will be negligible for the detection and the parameter estimation of gravitational waves, if the upper limit for magnetic fields of neutron stars are less than $\sim 10^{16}$G, which is the maximum magnetic field observed in the soft gamma repeaters and the anomalous X-ray pulsars up to date. We also discuss the implications of electromagnetic radiation from the inspiralling binary neutron stars for the precursory X-ray emission prior to the gamma ray burst observed by the Ginga satellite.Comment: 15 pages, no figures, accepted for publication in Ap