Neutrinos play an important role in compact star astrophysics: neutrino-heating is one of the main ingredients in core-collapse supernovae, neutrino-matter interactions determine the composition of matter in binary neutron star mergers and have among others a strong impact on conditions for heavy element nucleosynthesis and neutron star cooling is dominated by neutrino emission except for very old stars. Many works in the last decades have shown that in dense matter medium effects considerably change the neutrino-matter interaction rates, whereas many astrophysical simulations use analytic approximations which are often far from reproducing more complete calculations. In this work we present a scheme which allows to incorporate improved rates, for charged current interactions, into simulations and show as an example some results for core-collapse supernovae, where a noticeable difference is found in the location of the neutrinospheres of the low-energy neutrinos in the early post-bounce phase.

The impact of electron-capture (EC) cross sections on neutron-rich nuclei on the dynamics of core-collapse during infall and early post-bounce is studied performing spherically symmetric simulations in general relativity using a multigroup scheme for neutrino transport and full nuclear distributions in extended nuclear statistical equilibrium models. We thereby vary the prescription for EC rates on individual nuclei, the nuclear interaction for the EoS, the mass model for the nuclear statistical equilibrium distribution and the progenitor model. In agreement with previous works, we show that the individual EC rates are the most important source of uncertainty in the simulations, while the other inputs only marginally influence the results. A recently proposed analytic formula to extrapolate microscopic results on stable nuclei for EC rates to the neutron rich region, with a functional form motivated by nuclear-structure data and parameters fitted from large scale shell model calculations, is shown to lead to a sizable (16%) reduction of the electron fraction at bounce compared to more primitive prescriptions for the rates, leading to smaller inner core masses and slower shock propagation. We show that the EC process involves $\approx$ 130 different nuclear species around 86 Kr mainly in the N = 50 shell closure region, and establish a list of the most important nuclei to be studied in order to constrain the global rates.

Using an axisymmetric numerical code, we perform an extensive study of the magnetic field configurations in non-rotating neutron stars, varying the mass, magnetic field strength and the equation of state. We find that the monopolar (spherically symmetric) part of the norm of the magnetic field can be described by a single profile, that we fit by a simple eighth-order polynomial, as a function of the star's radius. This new generic profile applies remarkably well to all magnetized neutron star configurations built on hadronic equations of state. We then apply this profile to build magnetized neutron stars in spherical symmetry, using a modified Tolman-Oppenheimer-Volkov (TOV) system of equations. This new formalism produces slightly better results in terms of mass-radius diagrams than previous attempts to add magnetic terms to these equations. However, we show that such approaches are less accurate than usual, non-magnetized TOV models, and that consistent models must depart from spherical symmetry. Thus, our `universal' magnetic field profile is intended to serve as a tool for nuclear physicists to obtain estimates of magnetic field inside neutron stars, as a function of radial depth, in order to deduce its influence on composition and related properties. It possesses the advantage of being based on magnetic field distributions from realistic self-consistent computations, which are solutions of Maxwell's equations.

The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration.

In this work we present a newly constructed equation of state (EoS) --applicable to stellar core collapse and neutron star mergers--, including the entire baryon octet. Our EoS is compatible with the main constraints from nuclear physics and, in particular, with a maximum mass for cold beta-equilibrated neutron stars of 2 solar masses in agreement with recent observations. As an application of our new EoS, we compute numerical stationary models for rapidly (rigidly) rotating hot neutron stars. We consider maximum masses of hot stars, such as proto-neutron stars or hypermassive neutron stars in the post-merger phase of binary neutron star coalescence. The universality of I-Q-relations at nonzero temperature for fast rotating models, comparing a purely nuclear EoS with its counterparts containing Lambda-hyperons or the entire baryon octet, respectively, is discussed, too. We find that the I-Q universality is broken when thermal effects become important, whatever the value of entropy gradients in our models. Thus, the use of I-Q relations for the analysis of proto-neutron stars or merger remnant data, including gravitational wave signals from the last stages of binary neutron star mergers, should be regarded with care.

We present a numerical model for uniformly rotating superfluid neutron stars, for the first time with realistic microphysics including entrainment, in a fully general relativistic framework. We compute stationary and axisymmetric configurations of neutron stars composed of two fluids, namely superfluid neutrons and charged particles (protons and electrons), rotating with different rates around a common axis. Both fluids are coupled by entrainment, a non-dissipative interaction which in case of a non-vanishing relative velocity between the fluids, causes the fluid momenta being not aligned with the respective fluid velocities. We extend the formalism by Comer and Joynt (2003) in order to calculate the equation of state (EoS) and entrainment parameters for an arbitrary relative velocity. The resulting entrainment matrix fulfills all necessary sum rules and in the limit of small relative velocity our results agree with Fermi liquid theory ones, derived to lowest order in the velocity. This formalism is applied to two new nuclear equations of state which are implemented in the numerical model. We are able to obtain precise equilibrium configurations. Resulting density profiles and moments of inertia are discussed employing both EoSs, showing the impact of entrainment and the dependence on the EoS.

In this paper, we study in detail the role of general relativity on the global dynamics of giant pulsar glitches as exemplified by Vela. For this purpose, we carry out numerical simulations of the spin up triggered by the sudden unpinning of superfluid vortices. In particular, we compute the exchange of angular momentum between the core neutron superfluid and the rest of the star within a two-fluid model including both (non-dissipative) entrainment effects and (dissipative) mutual friction forces. Our simulations are based on a quasi-stationary approach using realistic equations of state (EoSs). We show that the evolution of the angular velocities of both fluids can be accurately described by an exponential law. The associated characteristic rise time $\tau_{\text{r}}$, which can be precisely computed from stationary configurations only, has a form similar to that obtained in the Newtonian limit. However, general relativity changes the structure of the star and leads to additional couplings between the fluids due to frame-dragging effects. As a consequence, general relativity can have a large impact on the actual value of $\tau_{\text{r}}$: the errors incurred by using Newtonian gravity are thus found to be as large as $\sim 40 \%$ for the models considered. Values of the rise time are calculated for Vela and compared with current observational limits. Finally, we study the amount of gravitational waves emitted during a glitch. Simple expressions are obtained for the corresponding characteristic amplitudes and frequencies. The detectability of glitches through gravitational wave observatories is briefly discussed.

We develop a detailed and self-consistent numerical model for extremely-magnetised white dwarfs, which have been proposed as progenitors of overluminous Type Ia supernovae. This model can describe fully-consistent equilibria of magnetic stars in axial symmetry, with rotation, general-relativistic effects and realistic equations of state (including electron-ion interactions and taking into account Landau quantisation of electrons due to the magnetic field). We study the influence of each of these ingredients onto the white dwarf structure and, in particular, on their maximum mass. We perform an extensive stability analysis of such objects, with their highest surface magnetic fields reaching $\sim 10^{13}~G$ (at which point the star adopts a torus-like shape). We confirm previous speculations that although very massive strongly magnetised white dwarfs could potentially exist, the onset of electron captures and pycnonuclear reactions may severely limit their stability. Finally, the emission of gravitational waves by these objects is addressed, showing no possibility of detection by the currently planned space-based detector eLISA.

The existence of gravitational radiation is a natural prediction of any relativistic description of the gravitational interaction. In this chapter, we focus on gravitational waves, as predicted by Einstein's general theory of relativity. First, we introduce those mathematical concepts that are necessary to properly formulate the physical theory, such as the notions of manifold, vector, tensor, metric, connection and curvature. Second, we motivate, formulate and then discuss Einstein's equation, which relates the geometry of spacetime to its matter content. Gravitational waves are later introduced as solutions of the linearized Einstein equation around flat spacetime. These waves are shown to propagate at the speed of light and to possess two polarization states. Gravitational waves can interact with matter, allowing for their direct detection by means of laser interferometers. Finally, Einstein's quadrupole formulas are derived and used to show that nonspherical compact objects moving at relativistic speeds are powerful gravitational wave sources.

The explosion of core-collapse supernova depends on a sequence of events taking place in less than a second in a region of a few hundred kilometers at the center of a supergiant star, after the stellar core approaches the Chandrasekhar mass and collapses into a proto-neutron star, and before a shock wave is launched across the stellar envelope. Theoretical efforts to understand stellar death focus on the mechanism which transforms the collapse into an explosion. Progress in understanding this mechanism is reviewed with particular attention to its asymmetric character. We highlight a series of successful studies connecting observations of supernova remnants and pulsars properties to the theory of core-collapse using numerical simulations. The encouraging results from first principles models in axisymmetric simulations is tempered by new puzzles in 3D. The diversity of explosion paths and the dependence on the pre-collapse stellar structure is stressed, as well as the need to gain a better understanding of hydrodynamical and MHD instabilities such as SASI and neutrino-driven convection. The shallow water analogy of shock dynamics is presented as a comparative system where buoyancy effects are absent. This dynamical system can be studied numerically and also experimentally with a water fountain. The potential of this complementary research tool for supernova theory is analyzed. We also review its potential for public outreach in science museums.

We present a self-consistent model for the study of the structure of a neutron star in strong magnetic fields. Starting from a microscopic Lagrangian, this model includes the effect of the magnetic field on the equation of state, the interaction of the electromagnetic field with matter (magnetisation), and anisotropies in the energy-momentum tensor, as well as general relativistic aspects. We build numerical axisymmetric stationary models and show the applicability of the approach with one example quark matter equation of state (EoS) often employed in the recent literature for studies of strongly magnetised neutron stars. For this EoS, the effect of inclusion of magnetic field dependence or the magnetisation do not increase the maximum mass significantly in contrast to what has been claimed by previous studies.

We present a new code, Lorene's Ghost (for Lorene's gravitational handling of spectral transport) developed to treat the problem of neutrino transport in supernovae with the use of spectral methods. First, we derive the expression for the nonrelativistic Liouville operator in doubly spherical coordinates (r, theta, phi, epsilon, Theta, Phi)$, and further its general relativistic counterpart. We use the 3 + 1 formalism with the conformally flat approximation for the spatial metric, to express the Liouville operator in the Eulerian frame. Our formulation does not use any approximations when dealing with the angular arguments (theta, phi, Theta, Phi), and is fully energy-dependent. This approach is implemented in a spherical shell, using either Chebyshev polynomials or Fourier series as decomposition bases. It is here restricted to simplified collision terms (isoenergetic scattering) and to the case of a static fluid. We finish this paper by presenting test results using basic configurations, including general relativistic ones in the Schwarzschild metric, in order to demonstrate the convergence properties, the conservation of particle number and correct treatment of some general-relativistic effects of our code. The use of spectral methods enables to run our test cases in a six-dimensional setting on a single processor.

We present numerical simulations of stellar core-collapse with spherically symmetric, general relativistic hydrodynamics up to black hole formation. Using the CoCoNuT code, with a newly developed grey leakage scheme for the neutrino treatment, we investigate the effects of including pions and \Lambda-hyperons into the equation of state at high densities and temperatures on the black hole formation process. Results show non-negligible differences between the models with reference equation of state without any additional particles and models with the extended ones. For the latter, the maximum masses supported by the proto-neutron star are smaller and the collapse to a black hole occurs earlier. A phase transition to hyperonic matter is observed when the progenitor allows for a high enough accretion rate onto the proto-neutron star. Rough estimates of neutrino luminosity from these collapses are given, too.

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.

In stellar core-collapse events matter is heated and compressed to densities above nuclear matter saturation density. For progenitors stars with masses above about 25 solar masses, which eventually form a black hole, the temperatures and densities reached during the collapse are so high that a traditional description in terms of electrons, nuclei, and nucleons is no longer adequate. We present here an improved equation of state which contains in addition pions and hyperons. They become abundant in the high temperature and density regime. We study the different constraints on such an equation of state, coming from both hyperonic data and observations of neutron star properties. In order to test the zero-temperature versions, we perform numerical simulations of the collapse of a neutron star with such additional particles to a black hole. We discuss the influence of the additional particles on the thermodynamic properties within the hot versions of the equation of state and we show that in regimes relevant to core-collapse and black hole formation, the effects of pions and hyperons on pressure, internal energy and sound speed are not negligible.

The Modified Newtonian Dynamics (MOND) generically predicts a violation of the strong version of the equivalence principle. As a result the gravitational dynamics of a system depends on the external gravitational field in which the system is embedded. This so-called external field effect is shown to imply the existence of an anomalous quadrupolar correction, along the direction of the external galactic field, in the gravitational potential felt by planets in the Solar System. We compute this effect by a numerical integration of the MOND equation in the presence of an external field, and deduce the secular precession of the perihelion of planets induced by this effect. We find that the precession effect is rather large for outer gaseous planets, and in the case of Saturn is comparable to, and in some cases marginally excluded by published residuals of precession permitted by the best planetary ephemerides.

The Modified Newtonian Dynamics (MOND) has been formulated as a modification of the Poisson equation for the Newtonian gravitational field. This theory generically predicts a violation of the strong version of the equivalence principle, and as a result the gravitational dynamics of a system depends on the external gravitational field in which the system is embedded. This so-called external field effect has been recently shown to imply the existence of an anomalous quadrupolar correction, along the direction of the external Galactic field, in the gravitational potential felt by planets in the Solar System. In this paper we confirm the existence of this effect by a numerical integration of the MOND equation in the presence of an external field, and compute the secular precession of the perihelion of planets induced by this effect. We find that the precession effect is rather large for outer gaseous planets, and in the case of Saturn is comparable to published residuals of precession obtained by Saturn range tracking data. The effect is much smaller for inner planets, but in the case of the Earth it appears to be in conflict for most of the MOND functions $\mu(y)$ with the very good constraint on the perihelion precession obtained from Jupiter VLBI data. The MOND functions that are compatible with this constraint appear to have a very rapid transition from the MONDian regime to the Newtonian one.

We present a new spectral-method-based algorithm for finding apparent horizons in three-dimensional space-like hypersurfaces without symmetries. While there are already a wide variety of algorithms for finding apparent horizons, our new algorithm does not suffer from the same weakness as previous spectral apparent horizon finders: namely the monopolar coefficient ($\ell=0$ in terms of the spherical harmonics decomposition) needed to be determined by a root-finding procedure. Hence, this leads to a much faster and more robust spectral apparent horizon finder. The finder is tested with the Kerr-Schild and Brill-Lindquist data. Our finder is accurate and is as efficient as the currently fastest methods developed recently by Schnetter (2003 Class. Quantum Grav. {\bf 20}, 4719) and Thornburg (2004 Class. Quantum Grav. {\bf 21}, 743). At typical resolutions it takes only 0.5 second to find the apparent horizon of a Kerr-Schild black hole with $a=0.9M$ to the accuracy $\sim 10^{-5}$ for the fractional error in the horizon's location on a 2 GHz processor.

A new numerical code for computing stationary axisymmetric rapidly rotating stars in general relativity is presented. The formulation is based on a fully constrained-evolution scheme for 3+1 numerical relativity using the Dirac gauge and maximal slicing. We use both the polytropic and MIT bag model equations of state to demonstrate that the code can construct rapidly rotating neutron star and strange star models. We compare numerical models obtained by our code and a well-established code, which uses a different gauge condition, and show that the two codes agree to high accuracy.

We present a new 3D general relativistic hydrodynamics code for simulations of stellar core collapse to a neutron star, as well as pulsations and instabilities of rotating relativistic stars. It uses spectral methods for solving the metric equations, assuming the conformal flatness approximation for the three-metric. The matter equations are solved by high-resolution shock-capturing schemes. We demonstrate that the combination of a finite difference grid and a spectral grid can be successfully accomplished. This "Mariage des Maillages" (French for grid wedding) approach results in high accuracy of the metric solver and allows for fully 3D applications using computationally affordable resources, and ensures long term numerical stability of the evolution. We compare our new approach to two other, finite difference based, methods to solve the metric equations. A variety of tests in 2D and 3D is presented, involving highly perturbed neutron star spacetimes and (axisymmetric) stellar core collapse, demonstrating the ability to handle spacetimes with and without symmetries in strong gravity. These tests are also employed to assess gravitational waveform extraction, which is based on the quadrupole formula.

We have developed a theoretical model and a numerical code for stationary rotating superfluid neutron stars in full general relativity. The underlying two-fluid model is based on Carter's covariant multi-fluid hydrodynamic formalism. The two fluids, representing the superfluid neutrons on one hand, and the protons and electrons on the other, are restricted to uniform rotation around a common axis, but are allowed to have different rotation rates. We have performed extensive tests of the numerical code, including quantitative comparisons to previous approximative results for these models. The results presented here are the first ``exact'' calculations of such models in the sense that no approximations (other than that inherent in a discretized numerical treatment) are used. Using this code we reconfirm the existence of prolate-oblate shaped configurations. We studied the dependency of the Kepler rotation limit and of the mass-density relation on the relative rotation rate. We further demonstrate how one can simulate a (albeit fluid) neutron-star ``crust'' by letting one fluid extend further outwards than the other, which results in interesting cases where the Kepler limit is actually determined by the outermost but slower fluid.

We describe recent progress in the numerical study of the structure of rapidly rotating superfluid neutron star models in full general relativity. The superfluid neutron star is described by a model of two interpenetrating and interacting fluids, one representing the superfluid neutrons and the second consisting of the remaining charged particles (protons, electrons, muons). We consider general stationary configurations where the two fluids can have different rotation rates around a common rotation axis. The previously discovered existence of configurations with one fluid in a prolate shape is confirmed.

The collapse of spherical neutron stars is studied in General Relativity. The initial state is a stable neutron star to which an inward radial kinetic energy has been added through some velocity profile. For two different equations of state and two different shapes of velocity profiles, it is found that neutron stars can collapse to black holes for high enough inward velocities, provided that their masses are higher than some minimal value, depending on the equation of state. For a polytropic equation of state of the form $p=K\rho^\gamma $, with $\gamma = 2$ it is found to be $1.16 (\frac{K}{0.1})^{0.5} \msol$, whereas for a more realistic one, it reads $0.36 \msol $. In some cases of collapse forming a black hole, part of the matter composing the initial neutron star can be ejected through a shock, leaving only a fraction of the initial mass to form a black hole. Therefore, black holes of very small masses can be formed and, in particular, the mass scaling relation, as a function of initial velocity, takes the form discovered by Choptuik for critical collapses.

Tensor-scalar theory of gravity allows the generation of gravitational waves from astrophysical sources, like Supernov\ae{}, even in the spherical case. That motivated us to study the collapse of a degenerate stellar core, within tensor-scalar gravity, leading to the formation of a neutron star through a bounce and the formation of a shock. We discuss in this paper the effects of the scalar field on the evolution of the system, as well as the appearance of strong non-perturbative effects of this scalar field (the so-called ``spontaneous scalarization''). As a main result, we describe the resulting gravitational monopolar radiation (form and amplitude) and discuss the possibility of its detection by the gravitational detectors currently under construction, taking into account the existing constraints on the scalar field. From the numerical point of view it is worthy to point out that we have developed a combined code which uses pseudo-spectral methods, for the evolution of the scalar field, and High Resolution Shock-Capturing schemes, for the evolution of the hydrodynamical system. Although this code has been used to integrate the field equations of that theory of gravity, in the spherically symmetric case, a by-product of the present work is to gain experience for an ulterior extension to multidimensional problems in Numerical Relativity of such numerical strategy.

Spherical neutron star models are studied within tensor-scalar theories of gravity. Particularly, it is shown that, under some conditions on the second derivative of the coupling function and on star's mass, for a given star there exist two strong-scalar-field solutions as well as the usual weak-field one. This last solution happens to be unstable and a star, becoming massive enough to allow for all three solutions, evolves to reach one of the strong field configurations. This transition is dynamically computed and it appears that the star radiates away the difference in energy between both states (a few $10^{-3} M_\odot c^2$) as gravitational radiation. Since part of the energy ($\sim 10^{-5} M_\odot c^2$) is injected into the star as kinetic energy, the velocity of star's surface can reach up to $10^{-2} c$. The waveform of this monopolar radiation is shown as well as the oscillations undergone by the star. These oscillations are also studied within the slowly-rotating approximation, in order to estimate an order of magnitude of the resulting quadrupolar radiation.

Complete tensor-scalar and hydrodynamic equations are presented and integrated, for a self-gravitating perfect fluid. The initial conditions describe unstable-equilibrium neutron star configuration, with a polytropic equation of state. They are necessary in order to follow the gravitational collapse (including full hydrodynamics) of this star toward a black hole and to study the resulting scalar gravitational wave. The amplitude of this wave, as well as the radiated energy dramatically increase above some critical value of the parameter of the coupling function, due to the spontaneous scalarization, an effect not present in Brans-Dicke theory. In most cases, the pressure of the collapsing fluid does not have a significant impact on the resulting signal. These kind of sources are not likely to be observed by future laser interferometric detectors (such as VIRGO or LIGO) of gravitational waves, if they are located at more than a few 100 kpc. However, spontaneous scalarization could be constrained if such a gravitational collapse is detected by its quadrupolar gravitational signal, since this latter is quite lower than the monopolar one.

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.