José A. Font

Morphology and Dynamics of Relativistic Jets

We present a comprehensive analysis of the morphology and dynamics of relativistic pressure-matched axisymmetric jets. The numerical simulations have been carried out with a high-resolution shock-capturing hydrocode based on an approximate relativistic Riemann solver derived from the spectral decomposition of the Jacobian matrices of relativistic hydrodynamics. We discuss the dependence of the jet morphology on several parameters, paying special attention to the relativistic effects caused by high Lorentz factors and large internal energies of the beam flow. The parameter space of our analysis is spanned by the ratio of the beam and ambient medium rest mass density (η), the beam Mach number (Mb), the beam Lorentz factor (Wb), and the adiabatic index (γ) of the equation of state (assuming an ideal gas). Both the ultrarelativistic regime (Wb ≥ 20) and the hypersonic regime (relativistic Mach number greater than 100) have been studied. Our results show that the enhancement of the effective inertial mass of the beam due to relativistic effects (through the specific enthalpy and the Lorentz factor) makes relativistic jets significantly more stable than Newtonian jets. We find that relativistic jets propagate very efficiently through the ambient medium, at speeds that agree very well with those obtained from an estimate based on a one-dimensional momentum balance. The propagation efficiency of a relativistic jet is an increasing function of the beam flow velocity. Relativistic jets seem to give rise to two different morphologies, according to the relevance of relativistic effects. Hot beams (i.e., with internal energies comparable to the beam rest-mass energy) show little internal structure (as they are almost in pressure equilibrium with their surroundings) and relatively smooth cocoons forming lobes near the head of the jet. Highly supersonic models, in which the kinematic relativistic effects due to high beam flow Lorentz factors dominate, display extended cocoons that are overpressured with respect to the environment. The cocoon thickness decreases, and its mean pressure increases with increasing beam Lorentz factor.

Numerical 3+1 General Relativistic Magnetohydrodynamics: A Local Characteristic Approach

We present a general procedure to solve numerically the general relativistic magnetohydrodynamics (GRMHD) equations within the framework of the 3+1 formalism. The work reported here extends our previous investigation in general relativistic hydrodynamics (Banyuls et al. 1997) where magnetic fields were not considered. The GRMHD equations are written in conservative form to exploit their hyperbolic character in the solution procedure. All theoretical ingredients necessary to build up high-resolution shock-capturing schemes based on the solution of local Riemann problems (i.e., Godunov-type schemes) are described. In particular, we use a renormalized set of regular eigenvectors of the flux Jacobians of the relativistic MHD equations. In addition, the paper describes a procedure based on the equivalence principle of general relativity that allows the use of Riemann solvers designed for special relativistic MHD in GRMHD. Our formulation and numerical methodology are assessed by performing various test simulations recently considered by different authors. These include magnetized shock tubes, spherical accretion onto a Schwarzschild black hole, equatorial accretion onto a Kerr black hole, and magnetized thick disks accreting onto a black hole and subject to the magnetorotational instability.

Three-dimensional relativistic simulations of rotating neutron-star collapse to a Kerr black hole

We present a new three-dimensional fully general-relativistic hydrodynamics code using high-resolution shock-capturing techniques and a conformal traceless formulation of the Einstein equations. Besides presenting a thorough set of tests which the code passes with very high accuracy, we discuss its application to the study of the gravitational collapse of uniformly rotating neutron stars to Kerr black holes. The initial stellar models are modeled as relativistic polytropes which are either secularly or dynamically unstable and with angular velocities which range from slow rotation to the mass-shedding limit. We investigate the gravitational collapse by carefully studying not only the dynamics of the matter, but also that of the trapped surfaces, i.e., of both the apparent and event horizons formed during the collapse. The use of these surfaces, together with the dynamical horizon framework, allows for a precise measurement of the black-hole mass and spin. The ability to successfully perform these simulations for sufficiently long times relies on excising a region of the computational domain which includes the singularity and is within the apparent horizon. The dynamics of the collapsing matter is strongly influenced by the initial amount of angular momentum in the progenitor star and, for initial models with sufficiently high angular velocities, the collapse can lead to the formation of an unstable disc in differential rotation. All of the simulations performed with uniformly rotating initial data and a polytropic or ideal-fluid equation of state show no evidence of shocks or of the presence of matter on stable orbits outside the black hole.

Relativistic simulations of rotational core collapse. II. Collapse dynamics and gravitational radiation

We have performed hydrodynamic simulations of relativistic rotational supernova core collapse in axisymmetry and have computed the gravitational radiation emitted by such an event. The Einstein equations are formulated using the confor- mally flat metric approximation, and the corresponding hydrodynamic equations are written as a first-order flux-conservative hyperbolic system. Details of the methodology and of the numerical code have been given in an accompanying paper. We have simulated the evolution of 26 models in both Newtonian and relativistic gravity. The initial configurations are dierentially rotat- ing relativistic 4=3-polytropes in equilibrium which have a central density of 10 10 gc m 3 . Collapse is initiated by decreasing the adiabatic index to some prescribed fixed value. The equation of state consists of a polytropic and a thermal part for a more realis- tic treatment of shock waves. Any microphysics like electron capture and neutrino transport is neglected. Our simulations show that the three dierent types of rotational supernova core collapse and gravitational waveforms identified in previous Newtonian simulations (regular collapse, multiple bounce collapse, and rapid collapse) are also present in relativistic gravity. However, ro- tational core collapse with multiple bounces is only possible in a much narrower parameter range in relativistic gravity. The relativistic models cover almost the same range of gravitational wave amplitudes (4 10 21 h TT 3 10 20 for a source at a distance of 10 kpc) and frequencies (60 Hz 1000 Hz) as the corresponding Newtonian ones. Averaged over all models, the total energy radiated in the form of gravitational waves is 8:2 10 8 Mc 2 in the relativistic case, and 3:6 10 8 Mc 2 in the Newtonian case. For all collapse models that are of the same type in both Newtonian and relativistic gravity, the gravitational wave signal is of lower amplitude. If the collapse type changes, either weaker or stronger signals are found in the relativistic case. For a given model, relativistic gravity can cause a large increase of the characteristic signal frequency of up to a factor of five, which may have important consequences for the signal detection. Our study implies that the prospects for detection of gravita- tional wave signals from axisymmetric supernova rotational core collapse do not improve when taking into account relativistic gravity. The gravitational wave signals obtained in our study are within the sensitivity range of the first generation laser interfer- ometer detectors if the source is located within the Local Group. An online catalogue containing the gravitational wave signal amplitudes and spectra of all our models is available at the URL http://www.mpa-garching.mpg.de/Hydro/hydro.html.

Accurate evolutions of unequal-mass neutron-star binaries: properties of the torus and short GRB engines

We present new results from accurate and fully general-relativistic simulations of the coalescence of unmagnetized binary neutron stars with various mass ratios. The evolution of the stars is followed through the inspiral phase, the merger, and the prompt collapse to a black hole, up until the appearance of a thick accretion disc, which is studied as it enters and remains in a regime of quasi-steady accretion. Although a simple ideal-fluid equation of state with Γ = 2 is used, this work presents a systematic study within a fully general-relativistic framework of the properties of the resulting black-hole–torus system produced by the merger of unequal-mass binaries. More specifically, we show that (1) the mass of the torus increases considerably with the mass asymmetry, and equal-mass binaries do not produce significant tori if they have a total baryonic mass Mtot 3.7 M⊙; (2) tori with masses Mtor ~ 0.2 M⊙ are measured for binaries with Mtot ~ 3.4 M⊙ and mass ratios q ~ 0.75–0.85; (3) the mass of the torus can be estimated by the simple expression , involving the maximum mass for the binaries and coefficients constrained from the simulations, and suggesting that the tori can have masses as large as for Mtot ~ 2.8 M⊙ and q ~ 0.75–0.85; (4) using a novel technique to analyze the evolution of the tori, we find no evidence for the onset of non-axisymmetric instabilities and that very little, if any, of their mass is unbound; (5) finally, for all the binaries considered, we compute the complete gravitational waveforms and the recoils imparted to the black holes, discussing the prospects of the detection of these sources for a number of present and future detectors.

Numerical Hydrodynamics and Magnetohydrodynamics in General Relativity

This article presents a comprehensive overview of numerical hydrodynamics and magneto-hydrodynamics (MHD) in general relativity. Some significant additions have been incorporated with respect to the previous two versions of this review (2000, 2003), most notably the coverage of general-relativistic MHD, a field in which remarkable activity and progress has occurred in the last few years. Correspondingly, the discussion of astrophysical simulations in general-relativistic hydrodynamics is enlarged to account for recent relevant advances, while those dealing with general-relativistic MHD are amply covered in this review for the first time. The basic outline of this article is nevertheless similar to its earlier versions, save for the addition of MHD-related issues throughout. Hence, different formulations of both the hydrodynamics and MHD equations are presented, with special mention of conservative and hyperbolic formulations well adapted to advanced numerical methods. A large sample of numerical approaches for solving such hyperbolic systems of equations is discussed, paying particular attention to solution procedures based on schemes exploiting the characteristic structure of the equations through linearized Riemann solvers. As previously stated, a comprehensive summary of astrophysical simulations in strong gravitational fields is also presented. These are detailed in three basic sections, namely gravitational collapse, black-hole accretion, and neutron-star evolutions; despite the boundaries, these sections may (and in fact do) overlap throughout the discussion. The material contained in these sections highlights the numerical challenges of various representative simulations. It also follows, to some extent, the chronological development of the field, concerning advances in the formulation of the gravitational field, hydrodynamics and MHD equations and the numerical methodology designed to solve them. To keep the length of this article reasonable, an effort has been made to focus on multidimensional studies, directing the interested reader to earlier versions of the review for discussions on one-dimensional works.

Numerical Hydrodynamics in General Relativity

The current status of numerical solutions for the equations of ideal general relativistic hydrodynamics is reviewed. Different formulations of the equations are presented, with special mention of conservative and hyperbolic formulations well-adapted to advanced numerical methods. A representative sample of available numerical schemes is discussed and particular emphasis is paid to solution procedures based on schemes exploiting the characteristic structure of the equations through linearized Riemann solvers. A comprehensive summary of relevant astrophysical simulations in strong gravitational fields, including gravitational collapse, accretion onto black holes and evolution of neutron stars, is also presented.

Relativistic simulations of rotational core collapse - I. Methods, initial models, and code tests

We describe an axisymmetric general relativistic code for rotational core collapse. The code evolves the coupled system of metric and fluid equations using the ADM 3 + 1 formalism and a conformally flat metric approximation of the Einstein equations. Within this approximation the ADM 3 + 1 equations reduce to a set of ve coupled non-linear elliptic equations for the metric components. The equations are discretized on a 2D grid in spherical polar coordinates and are solved by means of a Newton-Raphson iteration using a block elimination scheme to solve the diagonally dominant, sparse linear system arising within each iteration step. The relativistic hydrodynamics equations are formulated as a rst-order flux-conservative hyperbolic system and are integrated using high-resolution shock-capturing schemes based on Riemann solvers. We assess the quality of the conformally flat metric approximation for relativistic core collapse and present a comprehensive set of tests that the code successfully passed. The tests include relativistic shock tubes, the preservation of the rotation prole and of the equilibrium of rapidly and dierentially rotating neutron stars (approximated as rotating polytropes), spherical relativistic core collapse, and the conservation of rest-mass and angular momentum in dynamic spacetimes. The application of the code to relativistic rotational core collapse, with emphasis on the gravitational waveform signature, is presented in an accompanying paper.

Three Dimensional Numerical General Relativistic Hydrodynamics ; 1, Formulations, Methods, and Code Tests

This is the first in a series of papers on the construction and validation of a three-dimensional code for general relativistic hydrodynamics, and its application to general relativistic astrophysics. This paper studies the consistency and convergence of our general relativistic hydrodynamic treatment and its coupling to the spacetime evolutions described by the full set of Einstein equations with a perfect fluid source. The numerical treatment of the general relativistic hydrodynamic equations is based on high resolution shock capturing schemes. These schemes rely on the characteristic information of the system. A spectral decomposition for general relativistic hydrodynamics suitable for a general spacetime metric is presented. Evolutions based on three different approximate Riemann solvers coupled to four different discretizations of the Einstein equations are studied and compared. The coupling between the hydrodynamics and the spacetime (the right and left hand side of the Einstein equations) is carried out in a treatment which is second order accurate in {\it both} space and time. Convergence tests for all twelve combinations with a variety of test beds are studied, showing consistency with the differential equations and correct convergence properties. The test-beds examined include shocktubes, Friedmann-Robertson-Walker cosmology tests, evolutions of self-gravitating compact (TOV) stars, and evolutions of relativistically boosted TOV stars. Special attention is paid to the numerical evolution of strongly gravitating objects, e.g., neutron stars, in the full theory of general relativity, including a simple, yet effective treatment for the surface region of the star (where the rest mass density is abruptly dropping to zero).

Quasi-periodic accretion and gravitational waves from oscillating ‘toroidal neutron stars’ around a Schwarzschild black hole

We present general relativistic hydrodynamics simulations of constant specific angular momentum tori orbiting a Schwarzschild black hole. These tori are expected to form as a result of stellar gravitational collapse, binary neutron star merger or disruption, can reach very high rest-mass densities and behave effectively as neutron stars but with a toroidal topology (i.e. ‘toroidal neutron stars’). Here our attention is focused on the dynamical response of these objects to axisymmetric perturbations. We show that upon the introduction of perturbations, these systems either become unstable to the runaway instability or exhibit a regular oscillatory behaviour, resulting in a quasi-periodic variation of the accretion rate as well as of the mass quadrupole. The latter, in particular, is responsible for the emission of intense gravitational radiation for which the signal-to-noise ratio at the detector is comparable to or larger than the typical one expected in stellar-core collapse, making these new sources of gravitational waves potentially detectable. We discuss a systematic investigation of the parameter space in both the linear and non-linear regimes, providing estimates of how the gravitational radiation emitted depends on the mass of the torus and on the strength of the perturbation.