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.

This paper is devoted to the computation of compact binaries composed of one black hole and one neutron star. The objects are assumed to be on exact circular orbits. Standard 3+1 decomposition of Einstein equations is performed and the conformal flatness approximation is used. The obtained system of elliptic equations is solved by means of multi-domain spectral methods. Results are compared with previous work both in the high mass ratio limit and for one neutron star with very low compactness parameter. The accuracy of the present code is shown to be greater than with previous codes. Moreover, for the first time, some sequences containing one neutron star of realistic compactness are presented and discussed.

Relativistic spin-orbit and spin-spin couplings has been shown to modify the gravitational waveforms expected from inspiraling binaries with a black hole and a neutron star. As a result inspiral signals may be missed due to significant losses in signal-to-noise ratio, if precession effects are ignored in gravitational-wave searches. We examine the sensitivity of the anticipated loss of signal-to-noise ratio on two factors: the accuracy of the precessing waveforms adopted as the true signals and the expected distributions of spin-orbit tilt angles, given the current understanding of their physical origin. We find that the results obtained using signals generated by approximate techniques are in good agreement with the ones obtained by integrating the 2PN equations. This shows that a complete account of all high-order post-Newtonian effects is usually not necessary for the determination of detection efficiencies. Based on our current astrophysical expectations, large tilt angles are not favored and as a result the decrease in detection rate varies rather slowly with respect to the black hole spin magnitude and is within 20--30% of the maximum possible values.

In a recent investigation of the effects of precession on the anticipated detection of gravitational-wave inspiral signals from compact object binaries with moderate total masses, we found that (i) if precession is ignored, the inspiral detection rate can decrease by almost a factor of 10, and (ii) previously proposed ``mimic'' templates cannot improve the detection rate significantly (by more than a factor of 2). In this paper we propose a new family of templates that can improve the detection rate by factors of 5--6 in cases where precession is most important. Our proposed method for these new ``mimic'' templates involves a hierarchical scheme of efficient, two-parameter template searches that can account for a sequence of spikes that appear in the residual inspiral phase, after one corrects for the any oscillatory modification in the phase. We present our results for two cases of compact object masses (10 and 1.4 solar masses and 7 and 3 solar masses) as a function of spin properties. Although further work is needed to fully assess the computational efficiency of this newly proposed template family, we conclude that these ``spiky templates'' are good candidates for a family of precession templates used in realistic searches, that can improve detection rates of inspiral events.

We consider the problem of searching for gravitational waves emitted during the inspiral phase of binary systems when the orbital plane precesses due to relativistic spin-orbit coupling. Such effect takes place when the spins of the binary members are misaligned with respect to the orbital angular momentum. As a first step we assess the importance of precession specifically for the first-generation of LIGO detectors. We investigate the extent of the signal-to-noise ratio reduction and, hence, detection rate that occurs when precession effects are not accounted for in the template waveforms. We restrict our analysis to binary systems that undergo the so-called simple precession and have a total mass close to 10 solar mass. We find that for binary systems with rather high mass ratios (e.g., a 1.4 solar mass neutron star and a 10 solar mass black hole) the detection rate can decrease by almost an order of magnitude. Current astrophysical estimates of the rate of binary inspiral events suggest that LIGO could detect at most a few events per year, and therefore the reduction of the detection rate even by a factor of a few is critical. In the second part of our analysis, we examine whether the effect of precession could be included in the templates by capturing the main features of the phase modulation through a small number of extra parameters. Specifically we examine and tested for the first time the 3-parameter family suggested by Apostolatos. We find that, even though these ``mimic'' templates improve the detection rate, they are still inadequate in recovering the signal-to-noise ratio at the desired level. We conclude that a more complex template family is needed in the near future, still maintaining the number of additional parameters as small as possible in order to reduce the computational costs.

Coalescing binary neutron stars (NS) are expected to be an important source of gravitational waves (GW) detectable by laser interferometers. We present here a simple method for determining the compactness ratio M/R of NS based on the observed deviation of the GW energy spectrum from point-mass behavior at the end of an inspiral event. Our method is based on the properties of quasi-equilibrium binary NS sequences and does not require the computation of the full GW signal h(t). Combined with the measurement of the NS masses from the GW signal during inspiral, the determination of M/R will allow very strong constraints to be placed on the equation of state of nuclear matter at high densities.

We compare recent numerical results, obtained within a ``helical Killing vector'' (HKV) approach, on circular orbits of corotating binary black holes to the analytical predictions made by the effective one body (EOB) method (which has been recently extended to the case of spinning bodies). On the scale of the differences between the results obtained by different numerical methods, we find good agreement between numerical data and analytical predictions for several invariant functions describing the dynamical properties of circular orbits. This agreement is robust against the post-Newtonian accuracy used for the analytical estimates, as well as under choices of resummation method for the EOB ``effective potential'', and gets better as one uses a higher post-Newtonian accuracy. These findings open the way to a significant ``merging'' of analytical and numerical methods, i.e. to matching an EOB-based analytical description of the (early and late) inspiral, up to the beginning of the plunge, to a numerical description of the plunge and merger. We illustrate also the ``flexibility'' of the EOB approach, i.e. the possibility of determining some ``best fit'' values for the analytical parameters by comparison with numerical data.

Binary black hole systems in the pre-coalescence stage are numerically constructed by demanding that the associated spacetime admits a helical Killing vector. Comparison with third order post-Newtonian calculations indicates a rather good agreement until the innermost stable circular orbit.

We present a new approach to the problem of binary black holes in the pre-coalescence stage, i.e. when the notion of orbit has still some meaning. Contrary to previous numerical treatments which are based on the initial value formulation of general relativity on a (3-dimensional) spacelike hypersurface, our approach deals with the full (4-dimensional) spacetime. This permits a rigorous definition of the orbital angular velocity. Neglecting the gravitational radiation reaction, we assume that the black holes move on closed circular orbits, which amounts to endowing the spacetime with a helical Killing vector. We discuss the choice of the spacetime manifold, the desired properties of the spacetime metric, as well as the choice of the rotation state for the black holes. As a simplifying assumption, the space 3-metric is approximated by a conformally flat one. The problem is then reduced to solving five of the ten Einstein equations, which are derived here, as well as the boundary conditions on the black hole surfaces and at spatial infinity. We exhibit the remaining five Einstein equations and propose to use them to evaluate the error induced by the conformal flatness approximation. The orbital angular velocity of the system is computed through a requirement which reduces to the classical virial theorem at the Newtonian limit.

We present a spectral method for solving elliptic equations which arise in general relativity, namely three-dimensional scalar Poisson equations, as well as generalized vectorial Poisson equations of the type $\Delta \vec{N} + \lambda \vec{\nabla}(\vec{\nabla}\cdot \vec{N}) = \vec{S}$ with $\lambda \not= -1$. The source can extend in all the Euclidean space ${\bf R}^3$, provided it decays at least as $r^{-3}$. A multi-domain approach is used, along with spherical coordinates $(r,\theta,\phi)$. In each domain, Chebyshev polynomials (in $r$ or $1/r$) and spherical harmonics (in $\theta$ and $\phi$) expansions are used. If the source decays as $r^{-k}$ the error of the numerical solution is shown to decrease at least as $N^{-2(k-2)}$, where $N$ is the number of Chebyshev coefficients. The error is even evanescent, i.e. decreases as $\exp(-N)$, if the source does not contain any spherical harmonics of index $l\geq k -3$ (scalar case) or $l\geq k-5$ (vectorial case).

We present a numerical method to compute quasiequilibrium configurations of close binary neutron stars in the pre-coalescing stage. A hydrodynamical treatment is performed under the assumption that the flow is either rigidly rotating or irrotational. The latter state is technically more complicated to treat than the former one (synchronized binary), but is expected to represent fairly well the late evolutionary stages of a binary neutron star system. As regards the gravitational field, an approximation of general relativity is used, which amounts to solving five of the ten Einstein equations (conformally flat spatial metric). The obtained system of partial differential equations is solved by means of a multi-domain spectral method. Two spherical coordinate systems are introduced, one centered on each star; this results in a precise description of the stellar interiors. Thanks to the multi-domain approach, this high precision is extended to the strong field regions. The computational domain covers the whole space so that exact boundary conditions are set to infinity. Extensive tests of the numerical code are performed, including comparisons with recent analytical solutions. Finally a constant baryon number sequence (evolutionary sequence) is presented in details for a polytropic equation of state with gamma=2.