Highly accurate calculation of rotating neutron stars: Detailed description of the numerical methods

We give a detailed description of the recently developed multi-domain spectral method for constructing highly accurate general-relativistic models of rapidly rotating stars. For both "ordinary" and "critical" configurations, it is exhibited by means of representative examples, how the accuracy improves as the order of the approximation increases. Apart from homogeneous fluid bodies, we also discuss models of polytropic and strange stars.Comment: 22 pages, 4 figures, 9 tables, version accepted by A&

Numerical implementation of isolated horizon boundary conditions

We study the numerical implementation of a set of boundary conditions derived from the isolated horizon formalism, and which characterize a black hole whose horizon is in quasi-equilibrium. More precisely, we enforce these geometrical prescriptions as inner boundary conditions on an excised sphere, in the numerical resolution of the Conformal Thin Sandwich equations. As main results, we firstly establish the consistency of including in the set of boundary conditions a "constant surface gravity" prescription, interpretable as a lapse boundary condition, and secondly we assess how the prescriptions presented recently by Dain et al. for guaranteeing the well-posedness of the Conformal Transverse Traceless equations with quasi-equilibrium horizon conditions extend to the Conformal Thin Sandwich elliptic system. As a consequence of the latter analysis, we discuss the freedom of prescribing the expansion associated with the ingoing null normal at the horizon.Comment: 11 pages, 5 figures, references added and correcte

Equilibrium Configurations of Homogeneous Fluids in General Relativity

By means of a highly accurate, multi-domain, pseudo-spectral method, we investigate the solution space of uniformly rotating, homogeneous and axisymmetric relativistic fluid bodies. It turns out that this space can be divided up into classes of solutions. In this paper, we present two new classes including relativistic core-ring and two-ring solutions. Combining our knowledge of the first four classes with post-Newtonian results and the Newtonian portion of the first ten classes, we present the qualitative behaviour of the entire relativistic solution space. The Newtonian disc limit can only be reached by going through infinitely many of the aforementioned classes. Only once this limiting process has been consummated, can one proceed again into the relativistic regime and arrive at the analytically known relativistic disc of dust.Comment: 8 pages, colour figures, v3: minor additions including one reference, accepted by MNRA

Non-existence of stationary two-black-hole configurations

We resume former discussions of the question, whether the spin-spin repulsion and the gravitational attraction of two aligned black holes can balance each other. To answer the question we formulate a boundary value problem for two separate (Killing-) horizons and apply the inverse (scattering) method to solve it. Making use of results of Manko, Ruiz and Sanabria-G\'omez and a novel black hole criterion, we prove the non-existence of the equilibrium situation in question.Comment: 15 pages, 3 figures; Contribution to Juergen Ehlers Memorial Issue (GeRG journal

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

Relativistic Dyson Rings and Their Black Hole Limit

In this Letter we investigate uniformly rotating, homogeneous and axisymmetric relativistic fluid bodies with a toroidal shape. The corresponding field equations are solved by means of a multi-domain spectral method, which yields highly accurate numerical solutions. For a prescribed, sufficiently large ratio of inner to outer coordinate radius, the toroids exhibit a continuous transition to the extreme Kerr black hole. Otherwise, the most relativistic configuration rotates at the mass-shedding limit. For a given mass-density, there seems to be no bound to the gravitational mass as one approaches the black-hole limit and a radius ratio of unity.Comment: 13 pages, 1 table, 5 figures, v2: some discussion and two references added, accepted for publication in Astrophys. J. Let

Eccentric binary black-hole mergers: The transition from inspiral to plunge in general relativity

We study the transition from inspiral to plunge in general relativity by computing gravitational waveforms of non-spinning, equal-mass black-hole binaries. We consider three sequences of simulations, starting with a quasi-circular inspiral completing 1.5, 2.3 and 9.6 orbits, respectively, prior to coalescence of the holes. For each sequence, the binding energy of the system is kept constant and the orbital angular momentum is progressively reduced, producing orbits of increasing eccentricity and eventually a head-on collision. We analyze in detail the radiation of energy and angular momentum in gravitational waves, the contribution of different multipolar components and the final spin of the remnant. We find that the motion transitions from inspiral to plunge when the orbital angular momentum L=L_crit is about 0.8M^2. For L<L_crit the radiated energy drops very rapidly. Orbits with L of about L_crit produce our largest dimensionless Kerr parameter for the remnant, j=J/M^2=0.724. Generalizing a model recently proposed by Buonanno, Kidder and Lehner to eccentric binaries, we conjecture that (1) j=0.724 is the maximal Kerr parameter that can be obtained by any merger of non-spinning holes, and (2) no binary merger (even if the binary members are extremal Kerr black holes with spins aligned to the orbital angular momentum, and the inspiral is highly eccentric) can violate the cosmic censorship conjecture.Comment: Added sequence of long inspirals to the study. To match published versio

A universal constraint between charge and rotation rate for degenerate black holes surrounded by matter

We consider stationary, axially and equatorially symmetric systems consisting of a central rotating and charged degenerate black hole and surrounding matter. We show that $a^2+Q^2=M^2$ always holds provided that a continuous sequence of spacetimes can be identified, leading from the Kerr-Newman solution in electrovacuum to the solution in question. The quantity $a=J/M$ is the black hole's intrinsic angular momentum per unit mass, $Q$ its electric charge and $M$ the well known black hole mass parameter introduced by Christodoulou and Ruffini.Comment: 19 pages, 2 figures, replaced with published versio

Black Holes Surrounded by Uniformly Rotating Rings

Highly accurate numerical solutions to the problem of Black Holes surrounded by uniformly rotating rings in axially symmetric, stationary spacetimes are presented. The numerical methods developed to handle the problem are discussed in some detail. Related Newtonian problems are described and numerical results provided, which show that configurations can reach an inner mass-shedding limit as the mass of the central object increases. Exemplary results for the full relativistic problem for rings of constant density are given and the deformation of the event horizon due to the presence of the ring is demonstrated. Finally, we provide an example of a system for which the angular momentum of the central Black Hole divided by the square of its mass exceeds one.Comment: 12 pages, 14 figures, revtex, v4: minor changes, Eq. (17) corrected, corresponds to version in PR

Dirichlet Boundary Value Problems of the Ernst Equation

We demonstrate how the solution to an exterior Dirichlet boundary value problem of the axisymmetric, stationary Einstein equations can be found in terms of generalized solutions of the Backlund type. The proof that this generalization procedure is valid is given, which also proves conjectures about earlier representations of the gravitational field corresponding to rotating disks of dust in terms of Backlund type solutions.Comment: 22 pages, to appear in Phys. Rev. D, Correction of a misprint in equation (4