We present the first fully relativistic numerical calculations of differentially rotating Strange Quark Stars models for broad ranges of the maximum density and of the degree of differential rotation. Our simulations are performed with the very accurate and stable multi-domain spectral code FlatStar and use the MIT Bag model for describing strange quark matter. Our calculations based on a thorough exploration of the solution space show that the maximum mass of strange stars depends on both the degree of differential rotation and a type of solution, similarly to neutron stars described by a polytropic equations of state. The highest increase of the maximum mass (compared to the value for a non-rotating star) is obtained for models with a low degree of differential rotation. This highest mass is over four times larger than that of the equivalent non rotating configuration. Comparing our results with calculations done for realistic models of neutron stars, we conclude that with the help of differential rotation, strange stars can sustain masses much larger than stars made from nuclear matter for low degree of differential rotation which reinforces the hope of demonstrating, or of ruling out, the existence of strange matter through the observation by gravitational waves, by gamma rays or by neutrinos of the massive material object born from the merger of a compact binary system or during some supernova events.

We study the main astrophysical properties of differentially rotating neutron stars described as stationary and axisymmetric configurations of a moderately stiff $\Gamma=2$ polytropic fluid. The high level of accuracy and of stability of our relativistic multidomain pseudo-spectral code enables us to explore the whole solution space for broad ranges of the degree of differential rotation, but also of the stellar density and oblateness. Staying within an astrophysicaly motivated range of rotation profiles, we investigate the characteristics of neutron stars with maximal mass for all types of families of differentially rotating relativistic objects identified in a previous article Ansorg, Gondek-Rosinsla, Villain (2009). We find that the maximum mass depends on both the degree of differential rotation and on the type of solution. It turns out that the maximum allowed mass can be up to 4 times higher than what it is for non-rotating stars with the same equation of state. Such values are obtained for a modest degree of differential rotation but for one of the newly discovered type of solutions. Since such configurations of stars are not that extreme, this result may have important consequences for the gravitational wave signal to expect from coalescing neutron star binaries or from some supernovae events.

We prove that in Einstein-Maxwell theory the inequality $(8\pi J)^2+(4\pi Q^2)^2 < A^2$ holds for any sub-extremal axisymmetric and stationary black hole with arbitrary surrounding matter. Here $J, Q$, and $A$ are angular momentum, electric charge, and horizon area of the black hole, respectively.

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.

Binary black-hole systems with spins aligned with the orbital angular momentum are of special interest, as studies indicate that this configuration is preferred in nature. If the spins of the two bodies differ, there can be a prominent beaming of the gravitational radiation during the late plunge, causing a recoil of the final merged black hole. We perform an accurate and systematic study of recoil velocities from a sequence of equal-mass black holes whose spins are aligned with the orbital angular momentum, and whose individual spins range from a = +0.584 to -0.584. In this way we extend and refine the results of a previous study and arrive at a consistent maximum recoil of 448 +- 5 km/s for anti-aligned models as well as to a phenomenological expression for the recoil velocity as a function of spin ratio. This relation highlights a nonlinear behavior, not predicted by the PN estimates, and can be readily employed in astrophysical studies on the evolution of binary black holes in massive galaxies. An essential result of our analysis is the identification of different stages in the waveform, including a transient due to lack of an initial linear momentum in the initial data. Furthermore we are able to identify a pair of terms which are largely responsible for the kick, indicating that an accurate computation can be obtained from modes up to l=3. Finally, we provide accurate measures of the radiated energy and angular momentum, finding these to increase linearly with the spin ratio, and derive simple expressions for the final spin and the radiated angular momentum which can be easily implemented in N-body simulations of compact stellar systems. Our code is calibrated with strict convergence tests and we verify the correctness of our measurements by using multiple independent methods whenever possible.

To investigate the imprint on the gravitational-wave emission from extreme mass-ratio inspirals in non-pure Kerr spacetimes, we have studied the ``kludge'' waveforms generated in highly-accurate, numerically-generated spacetimes containing a black hole and a self-gravitating, homogeneous torus with comparable mass and spin. In order to maximize their impact on the produced waveforms, we have considered tori that are compact, massive and close to the central black hole, investigating under what conditions the LISA experiment could detect their presence. Our results show that for a large portion of the space of parameters the waveforms produced by EMRIs in these black hole-torus systems are indistinguishable from pure-Kerr waveforms. Hence, a ``confusion problem'' will be present for observations carried out over a timescale below or comparable to the dephasing time.

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.

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.

Using a multi domain spectral method, we investigate systematically the general-relativistic model for axisymmetric uniformly rotating, homogeneous fluid bodies generalizing the analytically known Maclaurin and Schwarzschild solutions. Apart from the curves associated with these solutions and a further curve of configurations that rotate at the mass shedding limit, two more curves are found to border the corresponding two parameter set of solutions. One of them is a Newtonian lens shaped sequence bifurcating from the Maclaurin spheroid sequence, while the other one corresponds to highly relativistic bodies with an infinite central pressure. The properties of the configuration for which both the gravitational and the baryonic masses, moreover angular velocity, angular momentum as well as polar red shift obtain their maximal values are discussed in detail. In particular, by comparison with the static Schwarzschild solution, we obtain an increase of 34.25% in the gravitational mass. Moreover, we provide exemplarily a discussion of angular velocity and gravitational mass on the entire solution class.

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.

We give a thorough investigation of sequences of uniformly rotating, homogeneous axisymmetric Newtonian equilibrium configurations that bifurcate from highly flattened Maclaurin spheroids. Each one of these sequences possesses a mass-shedding limit. Starting at this point, the sequences proceed towards the Maclaurin sequence and beyond. The first sequence leads to the well known Dyson rings, whereas the end points of the higher sequences are characterized by the formation of a two-body system, either a core-ring system (for the second, the fourth etc. sequence) or a two-ring system (for the third, the fifth etc. sequence). Although the general qualitative picture drawn by Eriguchi and Hachisu in the eighties has been confirmed, slight differences turned out in the interpretation of the origin of the first two-ring sequence and in the general appearance of fluid bodies belonging to higher sequences.

A new spectral code for constructing general-relativistic models of rapidly rotating stars with an unprecedented accuracy is presented. As a first application, we reexamine uniformly rotating homogeneous stars and compare our results with those obtained by several previous codes. Moreover, representative relativistic examples corresponding to highly flattened rotating bodies are given.