Olindo Zanotti

Dynamics of oscillating relativistic tori around Kerr black holes

We present a comprehensive numerical study of the dynamics of relativistic axisymmetric accretion tori with a power-law distribution of specific angular momentum orbiting in the background space‐time of a Kerr black hole. By combining general relativistic hydrodynamics simulations with a linear perturbative approach we investigate the main dynamical properties of these objects over a large parameter space. The astrophysical implications of our results extend and improve two interesting results that have been recently reported in the literature. First, the induced quasi-periodic variation of the mass quadrupole moment makes relativistic tori of nuclear matter densities, as those formed during the last stages of binary neutron star mergers, promising sources of gravitational radiation, potentially detectable by interferometric instruments. Secondly, p-mode oscillations in relativistic tori of low rest-mass densities could be used to explain high-frequency quasi-periodic oscillations observed in X-ray binaries containing a black hole candidate under conditions more generic than those considered so far. Ke yw ords: accretion, accretion discs ‐ gravitational waves ‐ hydrodynamics ‐ relativity.

An improved exact Riemann solver for relativistic hydrodynamics

A Riemann problem with prescribed initial conditions will produce one of three possible wave patterns corresponding to the propagation of the different discontinuities that will be produced once the system is allowed to relax. In general, when solving the Riemann problem numerically, the determination of the specific wave pattern produced is obtained through some initial guess which can be successively discarded or improved. We here discuss a new procedure, suitable for implementation in an exact Riemann solver in one dimension, which removes the initial ambiguity in the wave pattern. In particular we focus our attention on the relativistic velocity jump between the two initial states and use this to determine, through some analytic conditions, the wave pattern produced by the decay of the initial discontinuity. The exact Riemann problem is then solved by means of calculating the root of a nonlinear equation. Interestingly, in the case of two rarefaction waves, this root can even be found analytically. Our procedure is straightforward to implement numerically and improves the efficiency of numerical codes based on exact Riemann solvers

Dynamics of thick discs around Schwarzschild-de Sitter black holes

We consider the effects of a cosmological constant on the dynamics of constant angular momentum discs orbiting Schwarzschild-de Sitter black holes. The motivation behind this study is to investigate whether the presence of a radial force contrasting the black holes gravitational attraction can influence the occurrence of the runaway instability, a robust feature of the dynamics of constant angular momentum tori in Schwarzschild and Kerr spacetimes. In addition to the inner cusp near the black hole horizon through which matter can accrete onto the black hole, in fact, a positive cosmological constant introduces also an outer cusp through which matter can leave the torus without accreting onto the black hole. To assess the impact of this outflow on the development of the instability we have performed time-dependent and axisymmetric hydrodynamical simulations of equilibrium initial configurations in a sequence of background spacetimes of Schwarzschild-de Sitter black holes with increasing masses. The simulations have been performed with an unrealistic value for the cosmological constant which, however, yields sufficiently small discs to be resolved accurately on numerical grids and thus provides a first qualitative picture of the dynamics. The calculations, carried out for a wide range of initial conditions, show that the mass-loss from the outer cusp can have a considerable impact on the instability, with the latter being rapidly suppressed if the outflow is large enough.

Gravitational waves from oscillating accretion tori : Comparison between different approaches

Quasiperiodic oscillations of high density thick accretion disks orbiting a Schwarzschild black hole have been recently addressed as interesting sources of gravitational waves. The aim of this paper is to compare the gravitational waveforms emitted from these sources when computed using (variations of) the standard quadrupole formula and gauge-invariant metric perturbation theory. To this goal we evolve representative disk models using an existing general relativistic hydrodynamics code which has been previously employed in investigations of such astrophysical systems. Two are the main results of this work: First, for stable and marginally stable disks, no excitation of the black hole quasinormal modes is found. Second, we provide a simple, relativistic modification of the Newtonian quadrupole formula which, in certain regimes, yields excellent agreement with the perturbative approach. This holds true as long as back-scattering of GWs is negligible. Otherwise, any functional form of the quadrupole formula yields systematic errors {approx}10%.

Conformal and covariant Z4 formulation of the Einstein equations: strongly hyperbolic first-order reduction and solution with discontinuous Galerkin schemes

We present a strongly hyperbolic first-order formulation of the Einstein equations based on the conformal and covariant Z4 system (CCZ4) with constraint-violation damping, which we refer to as FO-CCZ4. As CCZ4, this formulation combines the advantages of a conformal and traceless formulation, with the suppression of constraint violations given by the damping terms, but being first order in time and space, it is particularly suited for a discontinuous Galerkin (DG) implementation. The strongly hyperbolic first-order formulation has been obtained by making careful use of first and second-order ordering constraints. A proof of strong hyperbolicity is given for a selected choice of standard gauges via an analytical computation of the entire eigenstructure of the FO-CCZ4 system. The resulting governing partial differential equations system is written in non-conservative form and requires the evolution of 58 unknowns. A key feature of our formulation is that the first-order CCZ4 system decouples into a set of pure ordinary differential equations and a reduced hyperbolic system of partial differential equations that contains only linearly degenerate fields. We implement FO-CCZ4 in a high-order path-conservative arbitrary-high-order-method-using-derivatives (ADER)-DG scheme with adaptive mesh refinement and local time-stepping, supplemented with a third-order ADER-WENO subcell finite-volume limiter in order to deal with singularities arising with black holes. We validate the correctness of the formulation through a series of standard tests in vacuum, performed in one, two and three spatial dimensions, and also present preliminary results on the evolution of binary black-hole systems. To the best of our knowledge, these are the first successful three-dimensional simulations of moving punctures carried out with high-order DG schemes using a first-order formulation of the Einstein equations.

Model for an optically thick torus in local thermodynamic equilibrium around a black hole

We propose a simple model for an optically thick radiative torus in local thermodynamic equilibrium around a Kerr black hole. The hydrodynamics structure, which is not affected by the radiation field, is the same as for the so--called polish doughnuts. Under the assumption of isentropic fluid and polytropic equation of state, a simple stationary and axisymmetric solution to the relativistic radiation hydrodynamics equations is possible, for which the temperature of the torus scales like the specific enthalpy. The astrophysical relevance of the model is briefly discussed.