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

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

The gravitational wave signal from a binary neutron star inspiral contains information on the nuclear equation of state. This information is contained in a combination of the tidal polarizability parameters of the two neutron stars and is clearest in the late inspiral, just before merger. We use the recently defined tidal extension of the effective one-body formalism to construct a controlled analytical description of the frequency-domain phasing of neutron star inspirals up to merger. Exploiting this analytical description we find that the tidal polarizability parameters of neutron stars can be measured by the advanced LIGO-Virgo detector network from gravitational wave signals having a reasonable signal-to-noise ratio of $\rho=16$. This measurability result seems to hold for all the nuclear equations of state leading to a maximum mass larger than $1.97M_\odot$. We also propose a promising new way of extracting information on the nuclear equation of state from a coherent analysis of an ensemble of gravitational wave observations of separate binary merger events.

A flexible spectral code for the study of general relativistic magnetohydrodynamics is presented. Aiming at investigating the physics of slowly rotating magnetized compact stars, this new code makes use of various physically motivated approximations. Among them, the relativistic anelastic approximation is a key ingredient of the current version of the code. In this article, we mainly outline the method, putting emphasis on algorithmic techniques that enable to benefit as much as possible of the non-dissipative character of spectral methods, showing also a potential astrophysical application and providing a few illustrative tests.

Direct observations of gravitational waves will open in the near future new windows on the Universe. Among the expected sources, instabilities of rotating compact astrophysical objects are waited to be detected with some impatience as this will sign the birth of ``gravitational waves asteroseismology'', a crucial way to improve our knowledge of matter equation of state in conditions that cannot be reproduced in a lab. However, the theoretical work needed to really get informations from to-be-detected signals is still quite large, numerical simulations having become a necessary key ingredient. This article tries to provide a short overview of the main physical topics involved in this field (general relativity, gravitational waves, instabilities of rotating fluids, {\it etc.}), concluding with a brief description of the work that was done in Paris-Meudon Observatory by Silvano Bonazzola and collaborators.

The influence of nucleons superfluidity on the beta relaxation time of degenerate neutron star cores, composed of neutrons, protons and electrons, is investigated. We numerically calculate the implied reduction factors for both direct and modified Urca reactions, with isotropic pairing of protons or anisotropic pairing of neutrons. We find that due to the non-zero value of the temperature and/or to the vanishing of anisotropic gaps in some directions of the phase-space, superfluidity does not always completely inhibit beta relaxation, allowing for some reactions if the superfluid gap amplitude is not too large in respect to both the typical thermal energy and the chemical potential mismatch. We even observe that if the ratio between the critical temperature and the actual temperature is very small, a suprathermal regime is reached for which superfluidity is almost irrelevant. On the contrary, if the gap is large enough, the composition of the nuclear matter can stay frozen for very long durations, unless the departure from beta equilibrium is at least as important as the gap amplitude. These results are crucial for precise estimation of the superfluidity effect on the cooling/slowing-down of pulsars and we provide online subroutines to be implemented in codes for simulating such evolutions.

With (non-barotropic) equations of state valid even when the neutron, proton and electron content of neutron star cores is not in beta equilibrium, we study inertial and composition gravity modes of relativistic rotating neutron stars. We solve the relativistic Euler equations in the time domain with a three dimensional numerical code based on spectral methods, in the slow rotation, relativistic Cowling and anelastic approximations. Principally, after a short description of the gravity modes due to smooth composition gradients, we focus our analysis on the question of how the inertial modes are affected by non-barotropicity of the nuclear matter. In our study, the deviation with respect to barotropicity results from the frozen composition of non-superfluid matter composed of neutrons, protons and electrons, when beta equilibrium is broken by millisecond oscillations. We show that already for moderatly fast rotating stars the increasing coupling between polar and axial modes makes those two cases less different than for very slowly rotating stars. In addition, as we directly solve the Euler equations, without coupling only a few number of spherical harmonics, we always found, for the models that we use, a discrete spectrum for the $l = m = 2$ inertial mode. Finally, we find that, for non-barotropic stars, the frequency of this mode, which is our main focus, decreases in a non-negligible way, whereas the time dependence of the energy transfer between polar and axial modes is substantially different due to the existence of low-frequencies gravity modes.

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.

We present a new hydro code based on spectral methods using spherical coordinates. The first version of this code aims at studying time evolution of inertial modes in slowly rotating neutron stars. In this article, we introduce the anelastic approximation, developed in atmospheric physics, using the mass conservation equation to discard acoustic waves. We describe our algorithms and some tests of the linear version of the code, and also some preliminary linear results. We show, in the Newtonian framework with differentially rotating background, as in the relativistic case with the strong Cowling approximation, that the main part of the velocity quickly concentrates near the equator of the star. Thus, our time evolution approach gives results analogous to those obtained by Karino {\it et al.} \cite{karino01} within a calculation of eigenvectors. Furthermore, in agreement with the work of Lockitch {\it et al.} \cite{lockandf01}, we found that the velocity seems to always get a non-vanishing polar part.