Our existence in the Universe resulted from a rare combination of circumstances. The same must hold for any highly developed extraterrestrial civilisation, and if they have ever existed in the Milky Way, they would likely be scattered over large distances in space and time. However, all technologically advanced species must be aware of the unique property of the galactic centre: it hosts Sagittarius A* (Sgr A*), the closest supermassive black hole to anyone in the Galaxy. A civilisation with sufficient technical know-how may have placed material in orbit around Sgr A* for research, energy extraction, and communication purposes. In either case, its orbital motion will necessarily be a source of gravitational waves. We show that a Jupiter-mass probe on the retrograde innermost stable circular orbit around Sgr A* emits, depending on the black hole spin, at a frequency of $f_{GW} = 0.63 - 1.07$ mHz and with a power of $P_{GW}=2.7 \times\, 10^{36} - 2.0 \times\, 10^{37}$ erg/s. We discuss that the energy output of a single star is sufficient to stabilise the location of an orbiting probe for a billion years against gravitational wave induced orbital decay. Placing and sustaining a device near Sgr A* is therefore astrophysically possible. Such a probe will emit an unambiguously artificial continuous gravitational wave signal that is observable with LISA-type detectors.

Context. The Event Horizon Telescope (EHT) collaboration recently obtained first images of the surroundings of the supermassive compact object M87* at the center of the galaxy M87. Aims. We want to develop a simple analytic disk model for the accretion flow of M87*. Compared to general-relativistic magnetohydrodynamic (GRMHD) models, it has the advantage of being independent of the turbulent character of the flow, and controlled by only few easy-to-interpret, physically meaningful parameters. We want to use this model to predict the image of M87* assuming that it is either a Kerr black hole, or an alternative compact object. Methods. We compute the synchrotron emission from the disk model and propagate the resulting light rays to the far-away observer by means of relativistic ray tracing. Such computations are performed assuming different spacetimes (Kerr, Minkowski, non-rotating ultracompact star, rotating boson star or Lamy spinning wormhole). We perform numerical fits of these models to the EHT data. Results. We discuss the highly-lensed features of Kerr images and show that they are intrinsically linked to the accretion-flow properties, and not only to gravitation. This fact is illustrated by the notion of secondary ring that we introduce. Our model of spinning Kerr black hole predicts mass and orientation consistent with the EHT interpretation. The non-Kerr images result in similar quality of the numerical fits and may appear very similar to Kerr images, once blurred to the EHT resolution. This implies that a strong test of the Kerr spacetime may be out of reach with the current data. We notice that future developments of the EHT could alter this situation. Conclusions. Our results show the importance of studying alternatives to the Kerr spacetime in order to be able to test the Kerr paradigm unambiguously.

A new code for computing fully general relativistic solutions of strongly magnetized rapidly rotating compact stars is developed as a part of the COCAL (Compact Object CALculator) code. The full set of Einstein's equations, Maxwell's equations and magnetohydrodynamic equations are consistently solved assuming perfect conductivity, stationarity, and axisymmetry, and strongly magnetized solutions associated with mixed poloidal and toroidal components of magnetic fields are successfully obtained in generic (non-circular) spacetimes. We introduce the formulation of the problem and the numerical method in detail, then present examples of extremely magnetized compact star solutions and their convergence tests. It is found that, in extremely magnetized stars, the stellar matter can be expelled from the region of strongest toroidal fields. Hence we conjecture that a toroidal electro-vacuum region may appear inside of the extremely magnetized compact stars, which may seem like the neutron star becoming the strongest toroidal solenoid coil in the universe.

We present the first fully relativistic study of gravitational radiation from bodies in circular equatorial orbits around the massive black hole at the Galactic Center, Sgr A* and we assess the detectability of various kinds of objects by the gravitational wave detector LISA. Our computations are based on the theory of perturbations of the Kerr spacetime and take into account the Roche limit induced by tidal forces in the Kerr metric. The signal-to-noise ratio in the LISA detector, as well as the time spent in LISA band, are evaluated. We have implemented all the computational tools in an open-source SageMath package, within the Black Hole Perturbation Toolkit framework. We find that white dwarfs, neutrons stars, stellar black holes, primordial black holes of mass larger than $10^{-4} M_\odot$, main-sequence stars of mass lower than $\sim 2.5\, M_\odot$ and brown dwarfs orbiting Sgr A* are all detectable in one year of LISA data with a signal-to-noise ratio above 10 during at least $10^5$ years in the slow inspiral towards either the innermost stable circular orbit (compact objects) or the Roche limit (main-sequence stars and brown dwarfs). The longest times in-band, of the order of $10^6$ years, are achieved for primordial black holes of mass $\sim 10^{-3} M_\odot$ down to $10^{-5} M_\odot$, depending on the spin of Sgr A*, as well as for brown dwarfs, just followed by white dwarfs and low mass main-sequence stars. The long time in-band of these objects makes Sgr A* a valuable target for LISA. We also consider bodies on close circular orbits around the massive black hole in the nucleus of the nearby galaxy M32 and find that, among them, compact objects and brown dwarfs stay for $10^3$ to $10^4$ years in LISA band with a 1-year signal-to-noise ratio above 10.

We show that the rotating generalization of Hayward's non-singular black hole previously studied in the literature is geodesically incomplete, and that its straightforward extension leads to a singular spacetime. We present another extension, which is devoid of any curvature singularity. The obtained metric depends on three parameters and, depending on their values, yields an event horizon or not. These two regimes, named respectively regular rotating Hayward black hole and naked rotating wormhole, are studied both numerically and analytically. In preparation for the upcoming results of the Event Horizon Telescope, the images of an accretion torus around Sgr A*, the supermassive object at the center of the Galaxy, are computed. These images contain, even in the absence of a horizon, a central faint region which bears a resemblance to the shadow of Kerr black holes and emphasizes the difficulty of claiming the existence of an event horizon from the analysis of strong-field images. The frequencies of the co- and contra-rotating orbits at the innermost stable circular orbit (ISCO) in this geometry are also computed, in the hope that quasi-periodic oscillations may permit to compare this model with Kerr's black hole on observational grounds.

Thermal dominated X-ray spectra of neutron stars in quiescent transient X-ray binaries and neutron stars that undergo thermonuclear bursts are sensitive to mass and radius. The mass-radius relation of neutron stars depends on the equation of state that governs their interior. Constraining this relation accurately is thus of fundamental importance to understand the nature of dense matter. In this context we introduce a pipeline to calculate realistic model spectra of rotating neutron stars with hydrogen and helium atmospheres. An arbitrarily fast rotating neutron star with a given equation of state generates the spacetime in which the atmosphere emits radiation. We use the Lorene/nrotstar code to compute the spacetime numerically and the ATM24 code to solve the radiative transfer equations self-consistently. Emerging specific intensity spectra are then ray-traced through the neutron star's spacetime from the atmosphere to a distant observer with the Gyoto code. Here, we present and test our fully relativistic numerical pipeline. To discuss and illustrate the importance of realistic atmosphere models we compare our model spectra to simpler models like the commonly used isotropic color-corrected blackbody emission. We highlight the importance of considering realistic model-atmosphere spectra together with relativistic ray tracing to obtain accurate predictions. We also insist on the crucial impact of the star's rotation on the observables. Finally, we close a controversy that has been appearing in the literature in the recent years regarding the validity of the ATM24 code.

The second-generation beam combiner at the Very Large Telescope (VLT), GRAVITY, observes the stars orbiting the compact object located at the center of our galaxy, with an unprecedented astrometric accuracy of 10 $\mu$as. The nature of this compact source is still unknown since black holes are not the only candidates explaining the four million solar masses at the Galactic center. Boson stars are such an alternative model to black holes. This paper focuses on the study of trajectories of stars orbiting a boson star and a Kerr black hole. We put in light strong differences between orbits obtained in both metrics when considering stars with sufficiently close pericenters to the compact object, typically $\lesssim 30~M$. Discovery of closer stars to the Galactic center than the S2 star by the GRAVITY instrument would thus be a powerful tool to possibly constrain the nature of the central source.

Carter and Lichnerowicz have established that barotropic fluid flows are conformally geodesic and obey Hamilton's principle. This variational approach can accommodate neutral, or charged and poorly conducting, fluids. We show that, unlike what has been previously thought, this approach can also accommodate perfectly conducting magnetofluids, via the Bekenstein-Oron description of ideal magnetohydrodynamics. When Noether symmetries associated with Killing vectors or tensors are present in geodesic flows, they lead to constants of motion polynomial in the momenta. We generalize these concepts to hydrodynamic flows. Moreover, the Hamiltonian descriptions of ideal magnetohydrodynamics allow one to cast the evolution equations into a hyperbolic form useful for evolving rotating or binary compact objects with magnetic fields in numerical general relativity. Conserved circulation laws, such as those of Kelvin, Alfv\'en and Bekenstein-Oron, emerge simply as special cases of the Poincar\'e-Cartan integral invariant of Hamiltonian systems. We use this approach to obtain an extension of Kelvin's theorem to baroclinic (non-isentropic) fluids, based on a temperature-dependent time parameter. We further extend this result to perfectly or poorly conducting baroclinic magnetoflows. Finally, in the barotropic case, such magnetoflows are shown to also be geodesic, albeit in a Finsler (rather than Riemann) space.

We address the astrophysical imaging of a family of deformed Kerr black holes (BHs). These are stationary, asymptotically flat black hole (BH) spacetimes, that are solutions of General Relativity minimally coupled to a massive, complex scalar field: Kerr BHs with scalar hair (KBHsSH). Such BHs bifurcate from the vacuum Kerr solution and can be regarded as a horizon within a rotating boson star. In a recent letter, it was shown that KBHsSH can exhibit very distinct shadows from the ones of their vacuum counterparts. The setup therein, however, considered the light source to be a celestial sphere sufficiently far away from the BH. Here, we analyse KBHsSH surrounded by an emitting torus of matter, simulating a more realistic astrophysical environment, and study the corresponding lensing of light as seen by a very far away observer, to appropriately model ground-based observations of Sgr A*. We find that the differences in imaging between KBHsSH and comparable vacuum Kerr BHs remain, albeit less dramatic than those observed for the corresponding shadows in the previous setup. In particular, we highlight two observables that might allow differentiating KBHsSH and Kerr BHs. The first is the angular size of the photon ring (in a Kerr spacetime) or lensing ring (in a KBHSH spacetime), the latter being significantly smaller for sufficiently non-Kerr-like spacetimes. The second is the existence of an edge in the intensity distribution (the photon ring in Kerr spacetime). This edge can disappear for very non-Kerr-like KBHsSH. It is plausible, therefore, that sufficiently precise Very Long Baseline Interferometric observations of BH candidates can constrain this model.

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.

Accretion disks play an important role in the evolution of their relativistic inner compact objects. The emergence of a new generation of interferometers will allow to resolve these accretion disks and provide more information about the properties of the central gravitating object. Due to this instrumental leap forward it is crucial to investigate the accretion disk physics near various types of inner compact objects now to deduce later constraints on the central objects from observations. A possible candidate for the inner object is the boson star. Here, we will try to analyze the differences between accretion structures surrounding boson stars and black holes. We aim at analysing the physics of circular geodesics around boson stars and study simple thick accretion tori (so-called Polish doughnuts) in the vicinity of these stars. We realize a detailed study of the properties of circular geodesics around boson stars. We then perform a parameter study of thick tori with constant angular momentum surrounding boson stars. This is done using the boson star models computed by a code constructed with the spectral solver library KADATH. We demonstrate that all the circular stable orbits are bound. In the case of a constant angular momentum torus, a cusp in the torus surface exists only for boson stars with a strong gravitational scalar field. Moreover, for each inner radius of the disk, the allowed specific angular momentum values lie within a constrained range which depends on the boson star considered. We show that the accretion tori around boson stars have different characteristics than in the vicinity of a black hole. With future instruments it could be possible to use these differences to constrain the nature of compact objects.

Stationary and axisymmetric solutions of relativistic rotating stars with strong mixed poloidal and toroidal magnetic fields are obtained numerically. Because of the mixed components of the magnetic field, the underlying stationary and axisymmetric spacetimes are no longer circular. These configurations are computed from the full set of the Einstein-Maxwell equations, Maxwell's equations and from first integrals and integrability conditions of the magnetohydrodynamic equilibrium equations. After a brief introduction of the formulation of the problem, we present the first results for highly deformed magnetized rotating compact stars.

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

In the case involving particles the necessary and sufficient condition for the Penrose process to extract energy from a rotating black hole is absorption of particles with negative energies and angular momenta. No torque at the black-hole horizon occurs. In this article we consider the case of arbitrary fields or matter described by an unspecified, general energy-momentum tensor $T_{\mu \nu}$ and show that the necessary and sufficient condition for extraction of a black hole's rotational energy is analogous to that in the mechanical Penrose process: absorption of negative energy and negative angular momentum. We also show that a necessary condition for the Penrose process to occur is for the Noether current (the conserved energy-momentum density vector) to be spacelike or past directed (timelike or null) on some part of the horizon. In the particle case, our general criterion for the occurrence of a Penrose process reproduces the standard result. In the case of relativistic jet-producing "magnetically arrested disks" we show that the negative energy and angular-momentum absorption condition is obeyed when the Blandford-Znajek mechanism is at work, and hence the high energy extraction efficiency up to $\sim 300\%$ found in recent numerical simulations of such accretion flows results from tapping the black hole's rotational energy through the Penrose process. We show how black-hole rotational energy extraction works in this case by describing the Penrose process in terms of the Noether current.

The equations governing null and timelike geodesics are derived within the 3+1 formalism of general relativity. In addition to the particle's position, they encompass an evolution equation for the particle's energy leading to a 3+1 expression of the redshift factor for photons. An important application is the computation of images and spectra in spacetimes arising from numerical relativity, via the ray-tracing technique. This is illustrated here by images of numerically computed stationary neutron stars and dynamical neutron stars collapsing to a black hole.

The Galactic center is an ideal laboratory to study strong-field general relativistic phenomena, as the supermassive black hole Sgr A* has the biggest angular Schwarzschild radius among all black holes. This article presents three different ways of using the immediate surroundings of Sgr A* as a probe in order to either constrain its spin, or even test the very nature of this compact object.

We calculate the "observed at infinity" image and spectrum of the accretion structure in Sgr A*, by modelling it as an optically thin, constant angular momentum ion torus in hydrodynamic equilibrium. The physics we consider includes a two-temperature plasma, a toroidal magnetic field, as well as radiative cooling by bremsstrahlung, synchrotron and inverse Compton processes. Our relativistic model has the virtue of being fully analytic and very simple, depending only on eight tunable parameters: the black hole spin and the inclination of the spin axis to our line of sight, the torus angular momentum, the polytropic index, the magnetic to total pressure ratio, the central values of density and electron temperature and the ratio of electron to ion temperatures. The observed image and spectrum are calculated numerically using the ray-tracing code GYOTO. Our results demonstrate that the ion torus model is able to account for the main features of the accretion structure surrounding Sgr A*.

We report work towards a relativistic formulation for modeling strongly magnetized neutron stars, rotating or in a close circular orbit around another neutron star or black hole, under the approximations of helical symmetry and ideal MHD. The quasi-stationary evolution is governed by the first law of thermodynamics for helically symmetric systems, which is generalized to include magnetic fields. The formulation involves an iterative scheme for solving the Einstein-Maxwell and relativistic MHD-Euler equations numerically. The resulting configurations for binary systems could be used as self-consistent initial data for studying their inspiral and merger.

The ability of the near future second generation VLTI instrument GRAVITY to constrain the properties of the Galactic center black hole is investigated. The Galactic center infrared flares are used as probes of strong-field gravity, within the framework of the hot spot model according to which the flares are the signature of a blob of gas orbiting close to the black hole's innermost stable circular orbit. Full general relativistic computations are performed, together with realistic observed data simulations, that lead to conclude that GRAVITY could be able to constrain the black hole's inclination parameter.

GYOTO, a general relativistic ray-tracing code, is presented. It aims at computing images of astronomical bodies in the vicinity of compact objects, as well as trajectories of massive bodies in relativistic environments. This code is capable of integrating the null and timelike geodesic equations not only in the Kerr metric, but also in any metric computed numerically within the 3+1 formalism of general relativity. Simulated images and spectra have been computed for a variety of astronomical targets, such as a moving star or a toroidal accretion structure. The underlying code is open source and freely available. It is user-friendly, quickly handled and very modular so that extensions are easy to integrate. Custom analytical metrics and astronomical targets can be implemented in C++ plug-in extensions independent from the main code.

These lecture notes are intended to introduce the theory of rotating stars in general relativity. The focus is put on the theoretical foundations, with a detailed discussion of the spacetime symmetries, the choice of coordinates and the derivation of the equations of structure from the Einstein equation. The global properties of rotating stars (mass, angular momentum, redshifts, orbits, etc.) are also introduced.

A fully geometrical treatment of general relativistic magnetohydrodynamics (GRMHD) is developed under the hypotheses of perfect conductivity, stationarity and axisymmetry. The spacetime is not assumed to be circular, which allows for greater generality than the Kerr-type spacetimes usually considered in GRMHD. Expressing the electromagnetic field tensor solely in terms of three scalar fields related to the spacetime symmetries, we generalize previously obtained results in various directions. In particular, we present the first relativistic version of the Soloviev transfield equation, subcases of which lead to fully covariant versions of the Grad-Shafranov equation and of the Stokes equation in the hydrodynamical limit. We have also derived, as another subcase of the relativistic Soloviev equation, the equation governing magnetohydrodynamical equilibria with purely toroidal magnetic fields in stationary and axisymmetric spacetimes.

Binary systems of compact objects with electromagnetic field are modeled by helically symmetric Einstein-Maxwell spacetimes with charged and magnetized perfect fluids. Previously derived thermodynamic laws for helically-symmetric perfect-fluid spacetimes are extended to include the electromagnetic fields, and electric currents and charges; the first law is written as a relation between the change in the asymptotic Noether charge $\dl Q$ and the changes in the area and electric charge of black holes, and in the vorticity, baryon rest mass, entropy, charge and magnetic flux of the magnetized fluid. Using the conservation laws of the circulation of magnetized flow found by Bekenstein and Oron for the ideal magnetohydrodynamic (MHD) fluid, and also for the flow with zero conducting current, we show that, for nearby equilibria that conserve the quantities mentioned above, the relation $\dl Q=0$ is satisfied. We also discuss a formulation for computing numerical solutions of magnetized binary compact objects in equilibrium with emphasis on a first integral of the ideal MHD-Euler equation.

A new method is described for constructing initial data for a binary neutron-star (BNS) system in quasi-equilibrium circular orbit. Two formulations for non-conformally flat data, waveless (WL) and near-zone helically symmetric (NHS), are introduced; in each formulation, the Einstein-Euler system, written in 3+1 form on an asymptotically flat spacelike hypersurface, is exactly solved for all metric components, including the spatially non-conformally flat potentials, and for irrotational flow. A numerical method applicable to both formulations is explained with an emphasis on the imposition of a spatial gauge condition. Results are shown for solution sequences of irrotational BNS with matter approximated by parametrized equations of state that use a few segments of polytropic equations of state. The binding energy and total angular momentum of solution sequences computed within the conformally flat -- Isenberg-Wilson-Mathews (IWM) -- formulation are closer to those of the third post-Newtonian (3PN) two point particles up to the closest orbits, for the more compact stars, whereas sequences resulting from the WL/NHS formulations deviate from the 3PN curve even more for the sequences with larger compactness. We think it likely that this correction reflects an overestimation in the IWM formulation as well as in the 3PN formula, by $\sim 1$ cycle in the gravitational wave phase during the last several orbits. The work suggests that imposing spatial conformal flatness results in an underestimate of the quadrupole deformation of the components of binary neutron-star systems in the last few orbits prior to merger.

Quite recently, some new mathematical approaches to black holes have appeared in the literature. They do not rely on the classical concept of event horizon -- which is very global, but on the local concept of hypersurfaces foliated by trapped surfaces. After a brief introduction to these new horizons, we focus on a viscous fluid analogy that can be developed to describe their dynamics, in a fashion similar to the membrane paradigm introduced for event horizons in the seventies, but with a significant change of sign of the bulk viscosity.

Recent observations of XTE J1739-285 suggest that it contains a neutron star rotating at 1122 Hz. Such rotational frequency would be the first for which the effects of rotation are significant. We study the consequences of very fast rotating neutron stars for the potentially observable quantities as stellar mass and pulsar period.

We present the general relativistic calculation of the energy release associated with a first order phase transition (PT) at the center of a rotating neutron star (NS). The energy release, E_rel, is equal to the difference in mass-energies between the initial (normal) phase configuration and the final configuration containing a superdense matter core, assuming constant total baryon number and the angular momentum. The calculations are performed with the use of precise pseudo-spectral 2-D numerical code; the polytropic equations of state (EOS) as well as realistic EOSs (Skyrme interactions, Mean Field Theory kaon condensate) are used. The results are obtained for a broad range of metastability of initial configuration and size of the new superdense phase core in the final configuration. For a fixed ``overpressure'', dP, defined as the relative excess of central pressure of a collapsing metastable star over the pressure of the equilibrium first-order PT, the energy release up to numerical accuracy does not depend on the stellar angular momentum and coincides with that for nonrotating stars with the same dP. When the equatorial radius of the superdense phase core is much smaller than the equatorial radius of the star, analytical expressions for the E_rel can be obtained: E_rel is proportional to dP^2.5 for small dP. At higher dP, the results of 1-D calculations of E_rel(dP) for non-rotating stars reproduce with very high precision exact 2-D results for fast-rotating stars. The energy release-angular momentum independence for a given overpressure holds also for the so-called ``strong'' PTs (that destabilise the star against the axi-symmetric perturbations), as well as for PTs with ``jumping'' over the energy barrier.

We calculate the energy release associated with a strong first-order phase transition, from normal phase N to an "exotic" superdense phase S, in a rotating neutron star. Such a phase transition, accompanied by a density jump rho_N --> rho_S, is characterized by rho_S/rho_N > 3/2(1+P_0/rho_N c^2), where P_0 is the pressure, at which phase transition occurs. Configurations with small S-phase cores are then unstable and collapse into stars with large S-phase cores. The energy release is equal to the difference in mass-energies between the initial (normal) configuration and the final configuration containing an S-phase core, total stellar baryon mass and angular momentum being kept constant. The calculations of the energy release are based on precise numerical 2-D calculations. Polytropic equations of state (EOSs) as well as realistic EOS with strong first-order phase transition due to kaon condensation are used. For polytropic EOSs, a large parameter space is studied. For a fixed "overpressure", dP, defined as the relative excess of central pressure of collapsing metastable star over the pressure of equilibrium first-order phase transition, the energy release E_rel does not depend on the stellar angular momentum. It coincides with that for nonrotating stars with the same dP. Therefore, results of 1-D calculations of E_rel(dP) for non-rotating stars can be used to predict, with very high precision, the outcome of much harder to perform 2-D calculations for rotating stars with the same dP. This result holds also for dP_min < dP < 0, corresponding to phase transitions with climbing over the energy barrier separating metastable N-phase configurations from those with an S-phase core. Such phase transitions could be realized in the cores of newly born, hot, pulsating neutron stars.

These lecture notes provide some introduction to the 3+1 formalism of general relativity, which is the foundation of most modern numerical relativity. The text is rather self-contained, with detailed calculations and numerous examples. Contents: 1. Introduction, 2. Geometry of hypersurfaces, 3. Geometry of foliations, 4. 3+1 decomposition of Einstein equation, 5. 3+1 equations for matter and electromagnetic field, 6. Conformal decomposition, 7. Asymptotic flatness and global quantities, 8. The initial data problem, 9. Choice of foliation and spatial coordinates, 10. Evolution schemes.

Equilibria of binary neutron stars in close circular orbits are computed numerically in a waveless formulation: The full Einstein-relativistic-Euler system is solved on an initial hypersurface to obtain an asymptotically flat form of the 4-metric and an extrinsic curvature whose time derivative vanishes in a comoving frame. Two independent numerical codes are developed, and solution sequences that model inspiraling binary neutron stars during the final several orbits are successfully computed. The binding energy of the system near its final orbit deviates from earlier results of third post-Newtonian and of spatially conformally flat calculations. The new solutions may serve as initial data for merger simulations and as members of quasiequilibrium sequences to generate gravitational wave templates, and may improve estimates of the gravitational-wave cutoff frequency set by the last inspiral orbit.

We investigate the viscosity driven instability in rotating relativistic stars by means of an iterative approach. We focus on polytropic rotating equilibrium stars and impose an m=2 perturbationin the lapse. We vary both the stiffness of the equation of state and the compactness of the star to study those effects on the value of the threshold. For a uniformly rotating star, the criterion T/W, where T is the rotational kinetic energy and W is the gravitational binding energy, mainly depends on the compactness of the star and takes values around 0.13 ~ 0.16, which differ slightly from that of Newtonian incompressible stars (~ 0.14). For differentially rotating stars, the critical value of T/W is found to span the range 0.17 - 0.25. This is significantly larger than the uniformly rotating case with the same compactness of the star. Finally we discuss a possibility of detecting gravitational waves from viscosity driven instability with ground-based interferometers.

We calculate energy release associated with a first order phase transition at the center of a rotating neutron star. The results are based on precise numerical 2-D calculations, in which both the polytropic equations of state (EOS) as well as realistic EOS of the normal phase are used. Presented results are obtained for a broad range of metastability of initial configuration and size of the new superdense phase core in the final configuration. For small radii of the superdense phase core analytical expressions for the energy release are obtained. For a fixed "overpressure" dP (the relative excess of central pressure of collapsing metastable star over the pressure of equilibrium first-order phase transition) the energy release remarkably does not depend on the stellar angular momentum and coincides with that for nonrotating stars with the same dP. The energy release is proportional to dP^2.5 for small dPs, when sufficiently precise brute force 2-D numerical calculations are out of question. For higher dPs, results of 1-D calculations of energy release for non-rotating stars are shown to reproduce, with very high precision, the exact 2-D results for rotating stars.

We derive from Einstein equation an evolution law for the area of a trapping or dynamical horizon. The solutions to this differential equation show a causal behavior. Moreover, in a viscous fluid analogy, the equation can be interpreted as an energy balance law, yielding to a positive bulk viscosity. These two features contrast with the event horizon case, where the non-causal evolution of the area and the negative bulk viscosity require teleological boundary conditions. This reflects the local character of trapping horizons as opposed to event horizons. Interpreting the area as the entropy, we propose to use an area/entropy evolution principle to select a unique dynamical horizon and time slicing in the Cauchy evolution of an initial marginally trapped surface.

This lecture provides some introduction to perfect fluid dynamics within the framework of general relativity. The presentation is based on the Carter-Lichnerowicz approach. It has the advantage over the more traditional approach of leading very straightforwardly to important conservation laws, such as the relativistic generalizations of Bernoulli's theorem or Kelvin's circulation theorem. It also permits to get easily first integrals of motion which are particularly useful for computing equilibrium configurations of relativistic stars in rotation or in binary systems. The presentation is relatively self-contained and does not require any a priori knowledge of general relativity. In particular, the three types of derivatives involved in relativistic hydrodynamics are introduced in detail: this concerns the Lie, exterior and covariant derivatives.

The back-bending phenomenon for compact stars is studied by means of analytical equations of state, for both constant-pressure phase transitions and the transitions through the mixed-phase region. We restrict ourselves to the case of normal rotating configurations, with baryon mass below the maximum allowable baryon mass for non-rotating stars. We use high-precision 2-D multi-domain spectral code LORENE to search the parameter space for possible instability regions, and possible changes in the stability character of rotating stars with phase transitions in their cores. Conditions on the density jump in constant-pressure phase transitions, leading to the existence of the unstable segments in the evolutionary sequences of spinning down isolated normal neutron stars, are derived. Conjectures concerning the existence of two disjoint families of non-rotating and rotating stationary configurations of neutron stars are formulated. Particular case of EOSs leading to marginal instability of static and rotating configurations is also studied: marginal instability point in non-rotating configurations continues to exist in all evolutionary spin-down tracks. The fate of rotating stars entering the region of instability is discussed. The change in radius, energy release, and spin-up associated with the corequake in rotating neutron star, triggered by the instability, are calculated. The energy release is found to be very weakly dependent on the angular momentum of collapsing star.

The isolated horizon formalism recently introduced by Ashtekar et al. aims at providing a quasi-local concept of a black hole in equilibrium in an otherwise possibly dynamical spacetime. In this formalism, a hierarchy of geometrical structures is constructed on a null hypersurface. On the other side, the 3+1 formulation of general relativity provides a powerful setting for studying the spacetime dynamics, in particular gravitational radiation from black hole systems. Here we revisit the kinematics and dynamics of null hypersurfaces by making use of some 3+1 slicing of spacetime. In particular, the additional structures induced on null hypersurfaces by the 3+1 slicing permit a natural extension to the full spacetime of geometrical quantities defined on the null hypersurface. This 4-dimensional point of view facilitates the link between the null and spatial geometries. We proceed by reformulating the isolated horizon structure in this framework. We also reformulate previous works, such as Damour's black hole mechanics, and make the link with a previous 3+1 approach of black hole horizon, namely the membrane paradigm. We explicit all geometrical objects in terms of 3+1 quantities, putting a special emphasis on the conformal 3+1 formulation. This is in particular relevant for the initial data problem of black hole spacetimes for numerical relativity. Illustrative examples are provided by considering various slicings of Schwarzschild and Kerr spacetimes.

An identity is derived from Einstein equation for any hypersurface H which can be foliated by spacelike two-dimensional surfaces. In the case where the hypersurface is null, this identity coincides with the two-dimensional Navier-Stokes-like equation obtained by Damour in the membrane approach to a black hole event horizon. In the case where H is spacelike or null and the 2-surfaces are marginally trapped, this identity applies to Hayward's trapping horizons and to the related dynamical horizons recently introduced by Ashtekar and Krishnan. The identity involves a normal fundamental form (normal connection 1-form) of the 2-surface, which can be viewed as a generalization to non-null hypersurfaces of the Hajicek 1-form used by Damour. This 1-form is also used to define the angular momentum of the horizon. The generalized Damour-Navier-Stokes equation leads then to a simple evolution equation for the angular momentum.

We present the first relativistic calculations of the final phase of inspiral of a binary system consisting of two stars built predominantely of strange quark matter (strange quark stars). We study the precoalescing stage within the Isenberg-Wilson-Mathews approximation of general relativity using a multidomain spectral method. A hydrodynamical treatment is performed under the assumption that the flow is either rigidly rotating or irrotational, taking into account the finite density at the stellar surface -- a distinctive feature with respect to the neutron star case. The gravitational-radiation driven evolution of the binary system is approximated by a sequence of quasi-equilibrium configurations at fixed baryon number and decreasing separation. We find that the innermost stable circular orbit (ISCO) is given by an orbital instability both for synchronized and irrotational systems. This constrasts with neutron stars for which the ISCO is given by the mass-shedding limit in the irrotational case. The gravitational wave frequency at the ISCO, which marks the end of the inspiral phase, is found to be 1400 Hz for two irrotational 1.35 Msol strange stars and for the MIT bag model of strange matter with massless quarks and a bag constant B=60 MeV/fm^3. Detailed comparisons with binary neutrons star models, as well as with third order Post-Newtonian point-mass binaries are given.

We present calculations of quasiequilibrium sequences of irrotational binary neutron stars based on realistic equations of state (EOS) for the whole neutron star interior. Three realistic nuclear EOSs of various softness and based on different microscopic models have been joined with a recent realistic EOS of the crust, giving in this way three different EOSs of neutron-star interior. Computations of quasiequilibrium sequences are performed within the Isenberg-Wilson-Mathews approximation to general relativity. For all evolutionary sequences, the innermost stable circular orbit (ISCO) is found to be given by mass-shedding limit (Roche lobe overflow). The EOS dependence on the last orbits is found to be quite important: for two 1.35 M_sol neutron stars, the gravitational wave frequency at the ISCO (which marks the end of the inspiral phase) ranges from 800 Hz to 1200 Hz, depending upon the EOS. Detailed comparisons with 3rd order post-Newtonian results for point-mass binaries reveals a very good agreement until hydrodynamical effects (dominated by high-order functions of frequency) become important, which occurs at a frequency ranging from 500 Hz to 1050 Hz, depending upon the EOS.

We investigate the evolution of rigidly and differentially rotating protoneutron stars (PNSs) during the first twenty seconds of their life. We solve the equations describing stationary axisymmetric configurations in general relativity coupled to a finite temperature, relativistic equation of state, to obtain a sequence of quasi-equilibrium configurations describing the evolution of newly born neutron stars. Our estimates show that the scale of variation of the angular velocity in a PNSs is of the order of 7-10 km. We obtain the maximum rotation frequency that can be reached as the protoneutron stars deleptonizes and cools down, as well as other relevant parameters such as total angular momentum or |T/W|. Our study shows that imposing physical constraints (conservation of baryonic mass and angular momentum) and choosing reasonable thermodynamical profiles as the star evolves gives results consistent with the energetics of more complex simulations of non-rotating PNSs. It appears to be unlikely that PNSs formed in nearly axisymmetric core collapse reach the critical angular velocity to undergo the bar mode instability. They could, however, undergo secular or low |T/W| rotational instabilities a few seconds after birth, resulting in a strong emission of gravitational waves retarded with respect to the neutrino luminosity peak. We also found that the geometry of strongly differentially rotating protoneutron stars can become toroidal-like for large values of the angular velocity, before reaching the mass shedding limit.

The effect of the hyperon softening of the equation of state (EOS) of dense matter on the spin evolution of isolated neutron stars is studied for a broad set of hyperonic EOSs. We use a multidomain 2-D code based on a spectral method, and show how important the precision of solving the equations of stationary motion is for the stability analysis. For some EOSs, the hyperon softening leads to spin-up by the angular momentum loss, in the form of the back bending phenomenon, for a rather broad range of stellar baryon mass. We show that large segments of the evolutionary tracks exhibiting the back bending behaviour in the moment-of-inertia - rotation-frequency plane, are actually unstable and therefore not astrophysically relevant. We show also that during the spin-up - angular-momentum-loss epoch, an isolated neutron star (e.g., a radio pulsar) can lose a sizable part of its initial angular momentum without significantly changing its rotation period. We propose also simple arguments and criteria allowing to connect the presence of a back bending epoch with the mass-radius relations and the stiffness and/or softness of the nucleon and hyperon EOSs of the neutron star core.

We study the viscosity driven (Jacobi-like) bar mode instability of rapidly rotating strange stars in general relativity. A triaxial, "bar shaped" compact star could be an efficient source of continuous wave gravitational radiation in the frequency range of the forthcoming interferometric detectors. We locate the secular instability point along several constant baryon mass sequences of uniformly rotating strange stars described by the MIT bag model. Contrary to neutron stars, strange stars with T/|W| (the ratio of the rotational kinetic energy to the absolute value of the gravitational potential energy) much lower than the corresponding value for the mass-shed limit can be secularly unstable to bar mode formation if shear viscosity is high enough to damp out any deviation from uniform rotation. The instability develops for a broad range of gravitational masses and rotational frequencies of strange quark stars. It imposes strong constraints on the lower limit of the frequency at the innermost stable circular orbit around rapidly rotating strange stars. The above results are robust for all linear self-bound equations of state assuming the growth time of the instability is faster than the damping timescale. We discuss astrophysical scenarios where triaxial instabilities (r-mode and viscosity driven instability) could be relevant in strange stars described by the standard MIT bag model of normal quark matter. Taking into account actual values of viscosities in strange quark matter and neglecting the magnetic field we show that Jacobi-like instability cannot develop in any astrophysicaly interesting temperature windows. The main result is that strange quark stars described by the MIT bag model can be accelerated to very high frequency in Low Mass X-ray binaries if the strange quark mass is ~ 200 MeV or higher.

Quasiequilibrium sequences of binary neutron stars are numerically calculated in the framework of the Isenberg-Wilson-Mathews (IWM) approximation of general relativity. The results are presented for both rotation states of synchronized spins and irrotational motion, the latter being considered as the realistic one for binary neutron stars just prior to the merger. We assume a polytropic equation of state and compute several evolutionary sequences of binary systems composed of different-mass stars as well as identical-mass stars with adiabatic indices gamma=2.5, 2.25, 2, and 1.8. From our results, we propose as a conjecture that if the turning point of binding energy (and total angular momentum) locating the innermost stable circular orbit (ISCO) is found in Newtonian gravity for some value of the adiabatic index gamma_0, that of the ADM mass (and total angular momentum) should exist in the IWM approximation of general relativity for the same value of the adiabatic index.

We present the first computations of quasiequilibrium binary neutron stars with different mass components in general relativity, within the Isenberg-Wilson-Mathews approximation. We consider both cases of synchronized rotation and irrotational motion. A polytropic equation of state is used with the adiabatic index gamma=2. The computations have been performed for the following combinations of stars: (M/R)_{star 1} vs. (M/R)_{star 2} = 0.12 vs. (0.12, 0.13, 0.14), 0.14 vs. (0.14, 0.15, 0.16), 0.16 vs. (0.16, 0.17, 0.18), and 0.18 vs. 0.18, where (M/R) denotes the compactness parameter of infinitely separated stars of the same baryon number. It is found that for identical mass binary systems there is no turning point of the binding energy (ADM mass) before the end point of the sequence (mass shedding point) in the irrotational case, while there is one before the end point of the sequence (contact point) in the synchronized case. On the other hand, in the different mass case, the sequence ends by the tidal disruption of the less massive star (mass shedding point). It is then more difficult to find a turning point in the ADM mass. Furthermore, we find that the deformation of each star depends mainly on the orbital separation and the mass ratio and very weakly on its compactness. On the other side, the decrease of the central energy density depends on the compactness of the star and not on that of the companion.

We perform some numerical study of the secular triaxial instability of rigidly rotating homogeneous fluid bodies in general relativity. In the Newtonian limit, this instability arises at the bifurcation point between the Maclaurin and Jacobi sequences. It can be driven in astrophysical systems by viscous dissipation. We locate the onset of instability along several constant baryon mass sequences of uniformly rotating axisymmetric bodies for compaction parameter $M/R = 0-0.275$. We find that general relativity weakens the Jacobi like bar mode instability, but the stabilizing effect is not very strong. According to our analysis the critical value of the ratio of the kinetic energy to the absolute value of the gravitational potential energy $(T/|W|)_{\rm crit}$ for compaction parameter as high as 0.275 is only 30% higher than the Newtonian value. The critical value of the eccentricity depends very weakly on the degree of relativity and for $M/R=0.275$ is only 2% larger than the Newtonian value at the onset for the secular bar mode instability. We compare our numerical results with recent analytical investigations based on the post-Newtonian expansion.

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.

Recycling strange stars to millisecond periods is studied within the framework of general relativity. We employ equations of state of strange quark matter based on the MIT Bag Model, with massive strange quarks and lowest order QCD interactions. The presence of the crust of normal matter is taken into account, with a bottom density assumed to be equal to the neutron-drip one. The calculations are performed by solving the exact 2-D equations for rigidly rotating stationary configurations in general relativity. Evolutionary tracks of accreting strange stars are computed, and their dependence on the initial strange star mass and on the fraction of the angular momentum transferred to the star from the infalling matter, is studied. The differences between recycling strange stars and neutron stars are pointed out.

We present the first results from a new method for computing spacetimes representing corotating binary black holes in circular orbits. The method is based on the assumption of exact equilibrium. It uses the standard 3+1 decomposition of Einstein equations and conformal flatness approximation for the 3-metric. Contrary to previous numerical approaches to this problem, we do not solve only the constraint equations but rather a set of five equations for the lapse function, the conformal factor and the shift vector. The orbital velocity is unambiguously determined by imposing that, at infinity, the metric behaves like the Schwarzschild one, a requirement which is equivalent to the virial theorem. The numerical scheme has been implemented using multi-domain spectral methods and passed numerous tests. A sequence of corotating black holes of equal mass is calculated. Defining the sequence by requiring that the ADM mass decrease is equal to the angular momentum decrease multiplied by the orbital angular velocity, it is found that the area of the apparent horizons is constant along the sequence. We also find a turning point in the ADM mass and angular momentum curves, which may be interpreted as an innermost stable circular orbit (ISCO). The values of the global quantities at the ISCO, especially the orbital velocity, are in much better agreement with those from third post-Newtonian calculations than with those resulting from previous numerical approaches.

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 study equilibrium sequences of close binary systems on circular orbits and composed of different mass stars with polytropic equation of state in Newtonian gravity. The solving method is a multidomain spectral method which we have recently developed. The computations are performed for both cases of synchronized and irrotational binary systems with adiabatic indices $\gamma = 3, 2.5, 2.25, 2$ and 1.8, and for three mass ratios: $M_1/M_2 =0.5, 0.2$ and 0.1. It is found that the equilibrium sequences always terminate at a mass shedding limit (appearance of a cusp on the surface of the less massive star). For synchronized binaries, this contrasts with the equal mass case, where the sequences terminate instead by a contact configuration. Regarding the turning point of the total angular momentum (or total energy) along a sequence, we find that it is difficult to get it for small mass ratios. Indeed, we do not find any turning points for $M_1/M_2 \le 0.5$ in the irrotational case. However, in the synchronized case, it becomes easier again to find one for mass ratios much smaller than $M_1/M_2 \sim 0.2$.

Observations of kHz quasi-periodic oscillations (QPOs) in the X-ray fluxes of low-mass X-ray binaries (LMXBs) have been used in attempts to constrain the external metric of the compact members of these binaries, as well as their masses and the equation of state of matter at supranuclear denisties. We compute the maximum orbital frequency of stable circular motion around uniformly rotating strange stars described by the MIT bag model. The calculations are performed for both normal and supramassive constant baryon mass sequences of strange stars rotating at all possible rates. We find the lower limits on the maximum orbital frequency and discuss them for a range of masses and for all rotational frequencies allowed in the model considered. We show that for slowly and moderately rotating strange stars the maximum value of orbital frequency can be a good indicator of the mass of the compact object. However, for rapidly rotating strange stars the same value of orbital frequency in the innermost stable circular orbit is obtained for stars with masses ranging from that of a planetoid to about three solar masses. At sufficiently high rotation rates of the strange star, the rotational period alone constrains the stellar mass to a surprisingly narrow range.

We study equilibrium sequences of close binary systems composed of identical polytropic stars in Newtonian gravity. The solving method is a multi-domain spectral method which we have recently developed. An improvement is introduced here for accurate computations of binary systems with stiff equation of state ($\gamma > 2$). The computations are performed for both cases of synchronized and irrotational binary systems with adiabatic indices $\gamma=3,~2.5,~2.25,~2$ and 1.8. It is found that the turning points of total energy along a constant-mass sequence appear only for $\gamma \ge 1.8$ for synchronized binary systems and $\gamma \ge 2.3$ for irrotational ones. In the synchronized case, the equilibrium sequences terminate by the contact between the two stars. On the other hand, for irrotational binaries, it is found that the sequences terminate at a mass shedding limit which corresponds to a detached configuration.

Calculations of the properties of rotating strange stars with crusts are performed within the framework of general relativity. We employ an equation of state (EOS) of strange quark matter based on the MIT Bag Model with massive strange quarks and lowest order QCD interactions. The crust is described by the BPS equation of state. A significant increase of the stellar radius is found close to the Keplerian (mass-shedding limit) configuration. This leads to the disappearance of the gap between the stellar surface and the innermost stable circular orbit (ISCO) at very high rotation rates, for a rather broad range of stellar masses. The Keplerian configuration for the strange star with crust corresponds to values of J, T/W which are about 20% smaller than in the case of bare strange stars. Because the Keplerian configuration is achieved due to the increase of the stellar oblateness, the Keplerian frequency (of the rotation) remains almost unaltered. The lack of the gap close to the Keplerian rotation could imply a more stringent limit on the frequency of the ISCO, if the existence of such a gap is supported by observations, as in the source 4U 1820-30 with the upper QPO frequency 1.07 kHz. During the spin-down of a magnetized strange quark star with crust, the crust matter is absorbed in the equatorial region by the strange matter core. The deconfinement of absorbed crust matter is a strongly exothermic process, which would influence the cooling curve of this compact object.

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).

It is shown that for very rapidly rotating low mass strange stars the marginally stable orbit is located above the stellar surface. This effect is explained by the very important role of the oblateness of the rotating strange star. The comparison with some ``academic'' examples is presented. This feature is purely Newtonian in its nature and has nothing to do with relativistic marginally stable orbit. The effect is very large and cannot be treated in a perturbative way. It seems that strange stars as a very dense self-bound objects are the only possibility in Nature to represent these toy models.

It has been suggested that observations of quasi periodic oscillations (QPOs) in the X-ray fluxes from low mass X-ray binaries can be used to constrain the mass of the compact object and the equation of state of its matter. A specific suggestion that the kHz QPO frequency saturates at the maximum orbital frequency has been widely considered. We examine rotating compact stars described by a new model of strange quark matter (Dey et al. 1998). We calculate the maximum orbital frequencies for both normal and supramassive strange stars described by the Dey model, and present these frequencies for sequences of equilibrium models with constant baryon mass. The maximum orbital frequencies for these compact objects are always higher than the kHz QPO frequencies observed to date.

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.

We study effects of the strange quark mass and of the QCD coupling constant in the MIT bag model of quark matter with lowest order quark-gluon interactions on the rapid rotation of strange stars. The influence of rotation on global parameters of strange stars is greater than in the case of the neutron star models. We show that independently of the strange quark mass and the coupling constans the ratio of the rotational kinetic energy to the absolute value of the gravitational potential energy for rotating strange stars is significantly higher than for ordinary neutron stars. This might indicate that rapidly rotating strange stars could be important sources of gravitational waves.

We compute numerical models of uniformly rotating strange stars (SS) in general relativity for the recently proposed QCD-based equation of state (EOS) of strange quark matter (Dey et al. 1998). Static models based on this EOS are characterised by a larger surface redshift than strange stars within the MIT bag model. The frequencies of the fastest rotating configurations described by Dey model are much higher than these for neutron stars (NS) and for the simplest SS MIT bag model. We determine a number of physical parameters for such stars and compare them with those obtained for NS. We construct constant baryon mass equilibrium sequences both normal and supramassive. Similarly to the NS a supramassive SS, prior to collapse to a black hole, spins up as it loses angular momentum. We find the upper limits on maximal masses and maximal frequencies of the rotating configurations. We show that the maximal rotating frequency for each of considered evolutionary sequences is never the Keplerian one. A normal and low mass supramassive strange stars gaining angular momentum always slows down just before reaching the Keplerian limit. For a high mass supramassive SS sequence the Keplerian configuration is the one with the lowest rotational frequency in the sequence. The value of $T/W$ for rapidly rotating SS of any mass is significantly higher than those for ordinary NS. For Keplerian configurations it increases as mass decreases. The results are robust for all linear self-bound equations of state.

Exact calculations of innermost stable circular orbit (ISCO) around rotating strange stars are performed within the framework of general relativity. Equations of state (EOS) of strange quark matter based on the MIT Bag Model with massive strange quarks and lowest order QCD interactions, are used. The presence of a solid crust of normal matter on rotating, mass accreting strange stars in LMXBs is taken into account. It is found that, contrary to neutron stars, above some minimum mass a gap always separates the ISCO and stellar surface, independently of the strange star rotation rate. For a given baryon mass of strange star, we calculate the ISCO frequency as function of stellar rotation frequency, from static to Keplerian configuration. For masses close to the maximum mass of static configurations the ISCO frequencies for static and Keplerian configurations are similar. However, for masses significantly lower than the maximum mass of static configurations, the minimum value of the ISCO frequency is reached in the Keplerian limit. Presence of a solid crust increases the ISCO frequency for the Keplerian configuration by about ten percent compared to that for a bare strange star of the same mass.

Exact models of uniformly rotating strange stars, built of self bound quark matter, are calculated within the framework of general relativity. This is made possible thanks to a new numerical technique capable to handle the strong density discontinuity at the surface of these stars. Numerical calculations are done for a simple MIT bag model equation of state of strange quark matter. Evolutionary sequences of models of rotating strange stars at constant baryon mass are calculated. Maximally rotating configurations of strange stars are determined, assuming that the rotation frequency is limited by the mass shedding and the secular instability with respect to axisymmetric perturbations. Exact formulae which give the dependence of the maximum rotation frequency, and of the maximum mass and corresponding radius of rotating configurations, on the value of the bag constant, are obtained. The values of T/W for rapidly rotating massive strange stars are significantly higher than those for ordinary neutron stars. This might indicate particular susceptibility of rapidly rotating strange stars to triaxial instabilities.

We present results of numerical computations of quasiequilibrium sequences of binary neutron stars with zero vorticity, in the general relativistic framework. The Einstein equations are solved under the assumption of a conformally flat spatial 3-metric (Wilson-Mathews approximation). The evolution of the central density of each star is monitored as the orbit shrinks in response to gravitational wave emission. For a compactification ratio M/R=0.14, the central density remains rather constant (with a slight increase, below 0.1%) before decreasing. For a higher compactification ratio M/R=0.17 (i.e. stars closer to the maximum mass configuration), a very small density increase (at most 0.3%) is observed before the decrease. This effect remains within the error induced by the conformally flat approximation. It can be thus concluded that no substantial compression of the stars is found, which would have indicated a tendency to individually collapse to black hole prior to merger. Moreover, no turning point has been found in the binding energy or angular momentum along evolutionary sequences, which may indicate that these systems do not have any innermost stable circular orbit (ISCO).

We report on general relativistic calculations of quasiequilibrium configurations of binary neutron stars in circular orbits with zero vorticity. These configurations are expected to represent realistic situations as opposed to corotating configurations. The Einstein equations are solved under the assumption of a conformally flat spatial 3-metric (Wilson-Mathews approximation). The velocity field inside the stars is computed by solving an elliptical equation for the velocity scalar potential. Results are presented for sequences of constant baryon number (evolutionary sequences). Although the central density decreases much less with the binary separation than in the corotating case, it still decreases. Thus, no tendency is found for the stars to individually collapse to black hole prior to merger.

We present spectral methods developed in our group to solve three-dimensional partial differential equations. The emphasis is put on equations arising from astrophysical problems in the framework of general relativity.

A multi-domain spectral method for computing very high precision 3-D stellar models is presented. The boundary of each domain is chosen in order to coincide with a physical discontinuity (e.g. the star's surface). In addition, a regularization procedure is introduced to deal with the infinite derivatives on the boundary that may appear in the density field when stiff equations of state are used. Consequently all the physical fields are smooth functions on each domain and the spectral method is absolutely free of any Gibbs phenomenon, which yields to a very high precision. The power of this method is demonstrated by direct comparison with analytical solutions such as MacLaurin spheroids and Roche ellipsoids. The relative numerical error reveals to be of the order of $10^{-10}$. This approach has been developed for the study of relativistic inspiralling binaries. It may be applied to a wider class of astrophysical problems such as the study of relativistic rotating stars too.

Various formalisms proposed recently for irrotational binary systems in general relativity are compared and explicit relations between them are exhibited. It is notably shown that the formalisms of (i) Teukolsky, (ii) Shibata and (iii) Bonazzola et al. (as corrected by Asada) are equivalent, i.e. yield exactly the same solution, although the former two are simpler than the latter one.

We conduct a direct comparison of three different representative numerical codes for constructing models of rapidly rotating neutron stars in general relativity. Our aim is to evaluate the accuracy of the codes and to investigate how the accuracy is affected by the choice of interpolation, domain of integration and equation of state. In all three codes, the same physical parameters, equations of state and interpolation method are used. We construct 25 selected models for polytropic equations of state and 22 models with realistic neutron star matter equations of state. The three codes agree well with each other (typical agreement is better than 0.1% to 0.01%) for most models, except for the extreme assumption of uniform density stars. We conclude that the codes can be used for the construction of highly accurate initial data configurations for polytropes of index N > 0.5 (which typically correspond to realistic neutron stars), when the domain of integration includes all space and for realistic equations with no phase transitions. With the exception of the uniform density case, the obtained values of physical parameters for the models considered in this paper can be regarded as ``standard'' and we display them in detail for all models.

Both observational and theoretical rates of binary neutron star coalescence give low prospects for detection of a single event by the initial LIGO/VIRGO interferometers. However, by utilizing at the best all the a priori information on the expected signal, a positive detection can be achieved. This relies on the hypothesis that $\gamma$-ray bursts are the electromagnetic signature of neutron star coalescences. The information about the direction of the source can then be used to add in phase the signals from different detectors in order (i) to increase the signal-to-noise ratio and (ii) to make the noise more Gaussian. Besides, the information about the time of arrival can be used to drastically decrease the observation time and thereby the false alarm rate. Moreover the fluence of the $\gamma$-ray emission gives some information about the amplitude of the gravitational signal. One can then add the signals from $\sim 10^4$ observation boxes ($\sim$ number of $\gamma$-ray bursts during 10 years) to yield a positive detection. Such a detection, based on the Maximum a Posteriori Probability Criterium, is a minimal one, in the sense that no information on the position and time of the events, nor on any parameter of the model, is collected. The advantage is that this detection requires an improvement of the detector sensitivity by a factor of only $\sim 1.5$ with respect to the initial LIGO/VIRGO interferometers, and that, if positive, it will confirm the $\gamma$-ray burst model.

An analytical scheme and a numerical method in order to study the effects of general relativity on the viscosity driven secular bar mode instability of rapidly rotating stars are presented. The approach consists in perturbing an axisymmetric and stationary configuration and studying its evolution by constructing a series of triaxial quasi-equilibrium configurations. These are obtained by solution of an approximate set of field equations where only the dominant non-axisymmetric terms are taken into account. The progress with respect to our former investigation consists in a higher relativistic order of the non-axisymmetric terms included into the computation, namely the fully three-dimensional treatment of the vector part of the space-time metric tensor as opposed to the scalar part, solely, in the former case. The scheme is applied to rotating stars built on a polytropic equation of state and compared to our previous results. The 3D-vector part turns out to inhibit the symmetry breaking efficiently. Nevertheless, the bar mode instability is still possible for an astrophysically relevant mass of M_ns=1.4 M_sun when a stiff polytropic equation of state with an adiabatic index of gamma=2.5 is employed. Triaxial neutron stars may be efficient emitters of gravitational waves and are thus potentially interesting sources for the forthcoming laser interferometric gravitational wave detectors such as LIGO, VIRGO and GEO600. From a numerical point of view, the solution of the three-dimensional minimal-distortion shift vector equation in spherical coordinates is an important achievement of our code.

A general relativistic version of the Euler equation for perfect fluid hydrodynamics is applied to a system of two neutron stars orbiting each other. In the quasi-equilibrium phase of the evolution of this system, a first integral of motion can be derived for certain velocity fields of the neutron star fluid including the (academic) case of co-rotation with respect to the orbital motion (synchronized binaries) and the realistic case of counter-rotation with respect to the orbital motion. The velocity field leading to this latter configuration can be computed by solving three-dimensional vector and scalar Poisson equations.

We study the possibility to detect the gravitational wave background generated by all the neutron stars in the Galaxy with only one gravitational wave interferometric detector. The proposed strategy consists in squaring the detector's output and searching for a sidereal modulation. The shape of the squared signal is computed for a disk and a halo distribution of neutron stars. The required noise stability of the interferometric detector is discussed. We argue that a possible population of old neutron stars, originating from a high stellar formation rate at the birth of the Galaxy and not emitting as radio pulsars, could be detected by the proposed technique in the low frequency range of interferometric experiments.

Triaxial neutron stars may be important sources of gravitational radiation for the forthcoming generation of interferometric gravitational wave detectors such as LIGO, VIRGO, and GEO600. We investigate the viscosity triggered bar mode secular instability of rapidly rotating neutron stars by means of a perturbation analysis of numerically constructed ``exact'' general relativistic axisymmetric star models. In the theoretical approach, only the dominant parts of the nonaxisymmetric terms of the 3D-Einstein equations are taken into account. A comparison of our results with previous studies of Newtonian polytropic stars confirms James' classical result gamma_crit=2.238 for the critical polytropic index. Beyond the Newtonian regime, gamma_crit reveals a slight increase toward highly relativistic configurations. Six out of twelve employed realistic dense matter equations of state admit the spontaneous symmetry breaking for masses above 1.6 M_sol.

We show that accreting neutron stars in binary systems or in Landau-Thorne-Zytkow objects are good candidates for continuous gravitational wave emission. Their gravitational radiation is strong enough to be detected by the next generation of detectors having a typical noise of 10^{-23} Hz^{-1/2}.

Contents: 1. Introduction 2. Spontaneous symmetry breaking 2.1 Review of classical results about Maclaurin/Jacobi ellipsoids 2.2 Spontaneous breaking of symmetry: a general phenomenom 2.3 Previous results for compressible Newtonian stars 2.4 Generation of gravitational waves 2.5 Finding the equilibrium configurations of a rotating star in the Newtonian regime 2.6 Extension to general relativity 2.7 First integral of fluid motion in general relativity 2.8 Gravitational field equations 2.9 Numerical results 3. CW emission from pulsars 4. Discussion and conclusion Appendix: Relativistic hydrodynamics in an accelerated frame

Continuous wave gravitational radiation from isolated rotating neutron stars is discussed. The general waveform and orders of magnitude for the amplitude are presented for various known pulsars. The specific case of gravitational radiation resulting from the distortion induced by the stellar magnetic field is presented. Finally some preliminary results about the signal from the whole population of neutron stars in the Galaxy are discussed.

Several applications of spectral methods to problems related to the relativistic astrophysics of compact objects are presented. Based on a proper definition of the analytical properties of regular tensorial functions we have developed a spectral method in a general sphericallike coordinate system. The applications include the investigation of spherically symmetric neutron star collapse as well as the solution of the coupled 2D-Einstein-Maxwell equations for magnetized, rapidly rotating neutron stars. In both cases the resulting codes are efficient and give results typically several orders of magnitude more accurate than equivalent codes based on finite difference schemes. We further report the current status of a 3D-code aiming at the simulation of non-axisymmetric neutron star collapse where we have chosen a tensor based numerical scheme.

The gravitational wave emission by a distorted rotating fluid star is computed. The distortion is supposed to be symmetric around some axis inclined with respect to the rotation axis. In the general case, the gravitational radiation is emitted at two frequencies: $\Omega$ and $2\Omega$, where $\Omega$ is the rotation frequency. The obtained formul\ae\ are applied to the specific case of a neutron star distorted by its own magnetic field. Assuming that the period derivative $\dot P$ of pulsars is a measure of their magnetic dipole moment, the gravitational wave amplitude can be related to the observable parameters $P$ and $\dot P$ and to a factor $\beta$ which measures the efficiency of a given magnetic dipole moment in distorting the star. $\beta$ depends on the nuclear matter equation of state and on the magnetic field distribution. The amplitude at the frequency $2\Omega$, expressed in terms of $P$, $\dot P$ and $\beta$, is independent of the angle $\alpha$ between the magnetic axis and the rotation axis, whereas at the frequency $\Omega$, the amplitude increases as $\alpha$ decreases. The value of $\beta$ for specific models of magnetic field distributions has been computed by means of a numerical code giving self-consistent models of magnetized neutron stars within general relativity. It is found that the distortion at fixed magnetic dipole moment is very dependent of the magnetic field distribution; a stochastic magnetic field or a superconductor stellar interior greatly increases $\beta$ with respect to the uniformly magnetized perfect conductor case and might lead to gravitational waves detectable by the VIRGO or LIGO interferometers. The amplitude modulation of the signal induced by the daily rotation of the Earth has been computed and specified to the case of the Crab pulsar and VIRGO

We investigate the effects of general relativity upon the non-axisymmetric ``bar'' mode secular instability of rapidly rotating stars, i.e. the relativistic and compressible analog of the transition from Maclaurin spheroids to Jacobi ellipsoids. Our method consists in perturbing a stationary axisymmetric configuration, constructed by a 2-D general relativistic numerical code, and taking into account only the dominant terms in the non-axisymmetric part of the 3-D relativistic equations. For a polytropic equation of state, we have determined, as a function of the degree of relativity, the critical adiabatic index $\gamma_{\rm crit}$ above which rapidly rotating stars can break their axial symmetry. A by-product of the present study is the confirmation of the Newtonian value $\gamma_{\rm crit} = 2.238$ obtained by James (1964). We have also considered neutron star models contructed upon twelve nuclear matter equations of state taken from the literature. We found that five equations of state from this sample allow the symmetry breaking for sufficiently high rotation velocities. For the others, the Keplerian velocity (mass-shedding from the equator) is reached before the axisymmetry is broken. Rotating neutron stars that break their axial symmetry can be an important source of gravitational waves for the LIGO/VIRGO interferometric detectors.

We present the first numerical solutions of the coupled Einstein-Maxwell equations describing rapidly rotating neutron stars endowed with a magnetic field. These solutions are fully relativistic and self-consistent, all the effects of the electromagnetic field on the star's equilibrium (Lorentz force, spacetime curvature generated by the electromagnetic stress-energy) being taken into account. The magnetic field is axisymmetric and poloidal. Five dense matter equations of state are employed. The partial differential equation system is integrated by means of a pseudo-spectral method. Various tests passed by the numerical code are presented. The effects of the magnetic field on neutron stars structure are then investigated, especially by comparing magnetized and non-magnetized configurations with the same baryon number. The deformation of the star induced by the magnetic field is important only for huge values of B (B>10^{10} T). The maximum mass as well as the maximum rotational velocity are found to increase with the magnetic field. The maximum allowable poloidal magnetic field is of the order of 10^{14} T (10^{18} G) and is reached when the magnetic pressure is comparable to the fluid pressure at the centre of the star. For such values, the maximum mass of neutron stars is found to increase by 13 to 29 % (depending upon the EOS) with respect to the maximum mass of non-magnetized stars.