Darwin-Riemann problems in general relativity

A review is given of recent results about the computation of irrotational Darwin-Riemann configurations in general relativity. Such configurations are expected to represent fairly well the late stages of inspiralling binary neutron stars.Comment: 20 pages, 11 PostScript figures, uses PTPTeX, to appear in the Proceedings of Yukawa International Seminar 99 "Black Holes and Gravitational Waves", edited by T. Nakamura & H. Kodama, Prog. Theor. Phys. Supp

Imaging a boson star at the Galactic center

Millimeter very long baseline interferometry will soon produce accurate images of the closest surroundings of the supermassive compact object at the center of the Galaxy, Sgr A*. These images may reveal the existence of a central faint region, the so-called shadow, which is often interpreted as the observable consequence of the event horizon of a black hole. In this paper, we compute images of an accretion torus around Sgr A* assuming this compact object is a boson star, i.e. an alternative to black holes within general relativity, with no event horizon and no hard surface. We show that very relativistic rotating boson stars produce images extremely similar to Kerr black holes, showing in particular shadow-like and photon-ring-like structures. This result highlights the extreme difficulty of unambiguously telling the existence of an event horizon from strong-field images.Comment: 21 pages, 9 figures, accepted in CQG; main difference wrt previous version is the last paragraph of the conclusio

Hyperboloidal evolution of test fields in three spatial dimensions

We present the numerical implementation of a clean solution to the outer boundary and radiation extraction problems within the 3+1 formalism for hyperbolic partial differential equations on a given background. Our approach is based on compactification at null infinity in hyperboloidal scri fixing coordinates. We report numerical tests for the particular example of a scalar wave equation on Minkowski and Schwarzschild backgrounds. We address issues related to the implementation of the hyperboloidal approach for the Einstein equations, such as nonlinear source functions, matching, and evaluation of formally singular terms at null infinity.Comment: 10 pages, 8 figure

An Efficient Pseudospectral Method for the Computation of the Self-force on a Charged Particle: Circular Geodesics around a Schwarzschild Black Hole

The description of the inspiral of a stellar-mass compact object into a massive black hole sitting at a galactic centre is a problem of major relevance for the future space-based gravitational-wave observatory LISA (Laser Interferometer Space Antenna), as the signals from these systems will be buried in the data stream and accurate gravitational-wave templates will be needed to extract them. The main difficulty in describing these systems lies in the estimation of the gravitational effects of the stellar-mass compact object on his own trajectory around the massive black hole, which can be modeled as the action of a local force, the self-force. In this paper, we present a new time-domain numerical method for the computation of the self-force in a simplified model consisting of a charged scalar particle orbiting a nonrotating black hole. We use a multi-domain framework in such a way that the particle is located at the interface between two domains so that the presence of the particle and its physical effects appear only through appropriate boundary conditions. In this way we eliminate completely the presence of a small length scale associated with the need of resolving the particle. This technique also avoids the problems associated with the impact of a low differentiability of the solution in the accuracy of the numerical computations. The spatial discretization of the field equations is done by using the pseudospectral collocation method and the time evolution, based on the method of lines, uses a Runge-Kutta solver. We show how this special framework can provide very efficient and accurate computations in the time domain, which makes the technique amenable for the intensive computations required in the astrophysically-relevant scenarios for LISA.Comment: 15 pages, 9 figures, Revtex 4. Minor changes to match published versio

Merger of black hole-neutron star binaries in full general relativity

We present our latest results for simulation for merger of black hole (BH)-neutron star (NS) binaries in full general relativity which is performed preparing a quasicircular state as initial condition. The BH is modeled by a moving puncture with no spin and the NS by the $\Gamma$-law equation of state with $\Gamma=2$ and corotating velocity field as a first step. The mass of the BH is chosen to be $\approx 3.2 M_{\odot}$ or $4.0M_{\odot}$, and the rest-mass of the NS $\approx 1.4 M_{\odot}$ with relatively large radius of the NS $\approx 13$--14 km. The NS is tidally disrupted near the innermost stable orbit but $\sim 80$--90% of the material is swallowed into the BH and resulting disk mass is not very large as $\sim 0.3M_{\odot}$ even for small BH mass $\sim 3.2M_{\odot}$. The result indicates that the system of a BH and a massive disk of $\sim M_{\odot}$ is not formed from nonspinning BH-NS binaries irrespective of BH mass, although a disk of mass $\sim 0.1M_{\odot}$ is a possible outcome for this relatively small BH mass range as $\sim 3$--4$M_{\odot}$. Our results indicate that the merger of low-mass BH and NS may form a central engine of short-gamma-ray bursts.Comment: 14 pages. To appear in a special issue of Classical and Quantum Gravity: New Frontiers in Numerical Relativit

Quasi-equilibrium binary black hole sequences for puncture data derived from helical Killing vector conditions

We construct a sequence of binary black hole puncture data derived under the assumptions (i) that the ADM mass of each puncture as measured in the asymptotically flat space at the puncture stays constant along the sequence, and (ii) that the orbits along the sequence are quasi-circular in the sense that several necessary conditions for the existence of a helical Killing vector are satisfied. These conditions are equality of ADM and Komar mass at infinity and equality of the ADM and a rescaled Komar mass at each puncture. In this paper we explicitly give results for the case of an equal mass black hole binary without spin, but our approach can also be applied in the general case. We find that up to numerical accuracy the apparent horizon mass also remains constant along the sequence and that the prediction for the innermost stable circular orbit is similar to what has been found with the effective potential method.Comment: 6 pages, 3 figures, 1 tabl

Binary black hole initial data for numerical general relativity based on post-Newtonian data

With the goal of taking a step toward the construction of astrophysically realistic initial data for numerical simulations of black holes, we for the first time derive a family of fully general relativistic initial data based on post-2-Newtonian expansions of the 3-metric and extrinsic curvature without spin. It is expected that such initial data provide a direct connection with the early inspiral phase of the binary system. We discuss a straightforward numerical implementation, which is based on a generalized puncture method. Furthermore, we suggest a method to address some of the inherent ambiguity in mapping post-Newtonian data onto a solution of the general relativistic constraints.Comment: 13 pages, 8 figures, RevTex

Gauge conditions for binary black hole puncture data based on an approximate helical Killing vector

We show that puncture data for quasicircular binary black hole orbits allow a special gauge choice that realizes some of the necessary conditions for the existence of an approximate helical Killing vector field. Introducing free parameters for the lapse at the punctures we can satisfy the condition that the Komar and ADM mass agree at spatial infinity. Several other conditions for an approximate Killing vector are then automatically satisfied, and the 3-metric evolves on a timescale smaller than the orbital timescale. The time derivative of the extrinsic curvature however remains significant. Nevertheless, quasicircular puncture data are not as far from possessing a helical Killing vector as one might have expected.Comment: 11 pages, 6 figures, 2 table

Analysis of LIGO data for gravitational waves from binary neutron stars

We report on a search for gravitational waves from coalescing compact binary systems in the Milky Way and the Magellanic Clouds. The analysis uses data taken by two of the three LIGO interferometers during the first LIGO science run and illustrates a method of setting upper limits on inspiral event rates using interferometer data. The analysis pipeline is described with particular attention to data selection and coincidence between the two interferometers. We establish an observational upper limit of $\mathcal{R}<$1.7 \times 10^{2}$per year per Milky Way Equivalent Galaxy (MWEG), with 90coalescence rate of binary systems in which each component has a mass in the range 1--3$M_\odot$.Comment: 17 pages, 9 figure