0000000000361707
AUTHOR
Pedro J. Montero
Gravitational Waves from the Papaloizou-Pringle Instability in Black-Hole-Torus Systems
Black hole (BH)--torus systems are promising candidates for the central engine of gamma-ray bursts (GRBs), and also possible outcomes of the collapse of supermassive stars to supermassive black holes (SMBHs). By three-dimensional general relativistic numerical simulations, we show that an $m=1$ nonaxisymmetric instability grows for a wide range of self-gravitating tori orbiting BHs. The resulting nonaxisymmetric structure persists for a timescale much longer than the dynamical one, becoming a strong emitter of large amplitude, quasiperiodic gravitational waves. Our results indicate that both, the central engine of GRBs and newly formed SMBHs, can be strong gravitational wave sources observa…
NADA: A new code for studying self-gravitating tori around black holes
We present a new two-dimensional numerical code called Nada designed to solve the full Einstein equations coupled to the general relativistic hydrodynamics equations. The code is mainly intended for studies of self-gravitating accretion disks (or tori) around black holes, although it is also suitable for regular spacetimes. Concerning technical aspects the Einstein equations are formulated and solved in the code using a formulation of the standard 3+1 (ADM) system, the so-called BSSN approach. A key feature of the code is that derivative terms in the spacetime evolution equations are computed using a fourth-order centered finite difference approximation in conjunction with the Cartoon metho…
3-D collapse of rotating stars to Kerr black holes
We study gravitational collapse of uniformly rotating neutron stars to Kerr black holes, using a new three-dimensional, fully general relativistic hydrodynamics code, which uses high-resolution shock-capturing techniques and a conformal traceless formulation of the Einstein equations. 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 su…
General Relativistic Simulations of Accretion Disks Around Tilted Kerr Black Holes
We simulate the dynamics of self-gravitating accretion disks around tilted Kerr black holes (BH) in full 3D general relativity. For this purpose we employ the EinsteinToolkit, using the thorn McLachlan for the evolution of the spacetime via the BSSN formalism of the Einstein equations and the thorn GRHydro for the evolution of the hydrodynamics, using a 3D Cartesian mesh with adaptive mesh refinement. We investigate the effects of the tilt angle between the disk angular momentum and BH spin vector on the dynamics of these systems as the disk evolves in the tilted spacetime. By evolving the spacetime and matter fields, we are able to observe how both BH and disk react and evolve in the tilte…
Numerical relativity simulations of thick accretion disks around tilted Kerr black holes
In this work we present 3D numerical relativity simulations of thick accretion disks around tilted Kerr BH. We investigate the evolution of three different initial disk models with a range of initial black hole spin magnitudes and tilt angles. For all the disk-to-black hole mass ratios considered (0.044-0.16) we observe significant black hole precession and nutation during the evolution. This indicates that for such mass ratios, neglecting the self-gravity of the disks by evolving them in a fixed background black hole spacetime is not justified. We find that the two more massive models are unstable against the Papaloizou-Pringle (PP) instability and that those PP-unstable models remain unst…
On the dynamics of tilted black hole-torus systems
We present results from three-dimensional, numerical relativity simulations of a {\it tilted} black hole-thick accretion disc system. The simulations are analysed using tracer particles in the disc which are advected with the flow. Such tracers, which we employ in these new simulations for the first time, provide a powerful means to analyse in detail the complex dynamics of tilted black hole-torus systems. We show how its use helps to gain insight in the overall dynamics of the system, discussing the origin of the observed black hole precession and the development of a global non-axisymmetric $m=1$ mode in the disc. Our three-dimensional simulations show the presence of quasi-periodic oscil…
Quasi-stationary solutions of self-gravitating scalar fields around black holes
Recent perturbative studies have shown the existence of long-lived, quasi-stationary configurations of scalar fields around black holes. In particular, such configurations have been found to survive for cosmological timescales, which is a requirement for viable dark matter halo models in galaxies based on such type of structures. In this paper we perform a series of numerical relativity simulations of dynamical non-rotating black holes surrounded by self-gravitating scalar fields. We solve numerically the coupled system of equations formed by the Einstein and the Klein-Gordon equations under the assumption of spherical symmetry using spherical coordinates. Our results confirm the existence …
Quasistationary solutions of self-gravitating scalar fields around black holes
Recent perturbative studies have shown the existence of long-lived, quasistationary configurations of scalar fields around black holes. In particular, such configurations have been found to survive for cosmological time scales, which is a requirement for viable dark matter halo models in galaxies based on such types of structures. In this paper we perform a series of numerical relativity simulations of dynamical nonrotating black holes surrounded by self-gravitating scalar fields. We solve numerically the coupled system of equations formed by the Einstein and the Klein-Gordon equations under the assumption of spherical symmetry using spherical coordinates. Our results confirm the existence …
On the runaway instability of self-gravitating torus around black holes
Black holes surrounded by self-gravitating tori are astrophysical systems which may naturally form following the core collapse of a massive star or the merger of two neutron stars. We present here results from fully general relativistic numerical simulations of such systems in order to assess the influence of the torus self-gravity on the onset of the so-called runaway instability. This instability, which might drive the rapid accretion of the disk on shorter timescales than those required to power a relativistic fireball, potentially challenges current models of gamma-ray bursts. Our simulations indicate that the self-gravity of the torus does not actually favour the onset of the instabili…
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 spacetime 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. Firstly, the induced quasi-periodic variation of the mass quadrupole moment makes relativistic tori of nuclear matter densities, as those formed duri…
Dynamical formation of a Reissner-Nordström black hole with scalar hair in a cavity
In a recent Letter [Sanchis-Gual et al., Phys. Rev. Lett. 116, 141101 (2016)], we presented numerical relativity simulations, solving the full Einstein--Maxwell--Klein-Gordon equations, of superradiantly unstable Reissner-Nordstr\"om black holes (BHs), enclosed in a cavity. Low frequency, spherical perturbations of a charged scalar field trigger this instability. The system's evolution was followed into the nonlinear regime, until it relaxed into an equilibrium configuration, found to be a hairy BH: a charged horizon in equilibrium with a scalar field condensate, whose phase is oscillating at the (final) critical frequency. Here, we investigate the impact of adding self-interactions to the …
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 modelled 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 dynam…
Measuring the black hole spin direction in 3D Cartesian numerical relativity simulations
We show that the so-called flat-space rotational Killing vector method for measuring the Cartesian components of a black hole spin can be derived from the surface integral of Weinberg's pseudotensor over the apparent horizon surface when using Gaussian normal coordinates in the integration. Moreover, the integration of the pseudotensor in this gauge yields the Komar angular momentum integral in a foliation adapted to the axisymmetry of the spacetime. As a result, the method does not explicitly depend on the evolved lapse $\ensuremath{\alpha}$ and shift ${\ensuremath{\beta}}^{i}$ on the respective time slice, as they are fixed to Gaussian normal coordinates while leaving the coordinate label…
Influence of self-gravity on the runaway instability of black-hole-torus systems.
Results from the first fully general relativistic numerical simulations in axisymmetry of a system formed by a black hole surrounded by a self-gravitating torus in equilibrium are presented, aiming to assess the influence of the torus self-gravity on the onset of the runaway instability. We consider several models with varying torus-to-black hole mass ratio and angular momentum distribution orbiting in equilibrium around a non-rotating black hole. The tori are perturbed to induce the mass transfer towards the black hole. Our numerical simulations show that all models exhibit a persistent phase of axisymmetric oscillations around their equilibria for several dynamical timescales without the …
Quasistationary solutions of scalar fields around accreting black holes
Massive scalar fields can form long-lived configurations around black holes. These configurations, dubbed quasi-bound states, have been studied both in the linear and nonlinear regimes. In this paper we show that quasi-bound states can form in a dynamical scenario in which the mass of the black hole grows significantly due to the capture of infalling matter. We solve the Klein-Gordon equation numerically in spherical symmetry, mimicking the evolution of the spacetime through a sequence of analytic Schwarzschild black hole solutions of increasing mass. It is found that the frequency of oscillation of the quasi-bound states decreases as the mass of the black hole increases. In addition, accre…
Dynamics of magnetized relativistic tori oscillating around black holes
We present a numerical study of the dynamics of magnetized, relativistic, non-self-gravitating, axisymmetric tori orbiting in the background spacetimes of Schwarzschild and Kerr black holes. The initial models have a constant specific angular momentum and are built with a non-zero toroidal magnetic field component, for which equilibrium configurations have recently been obtained. In this work we extend our previous investigations which dealt with purely hydrodynamical thick discs, and study the dynamics of magnetized tori subject to perturbations which, for the values of the magnetic field strength considered here, trigger quasi-periodic oscillations lasting for tens of orbital periods. Ove…
Dynamics of oscillating magnetized relativistic tori around a Schwarzschild black hole
We present a comprehensive numerical study of the dynamics of magnetized relativistic axisymmetric tori orbiting in the background spacetime of a Schwarzschild black hole. The tori are modeled as having a purely toroidal magnetic field and a constant distribution of the specific angular momentum. Following previous investigations of tori in a purely hydrodynamical context, the dynamics of these objects has been studied upon the introduction of a perturbation which, for the values of the magnetic field considered here, triggers quasi-periodic oscillations (QPOs) lasting tens of orbital periods. As in the hydrodynamical case, the spectral distribution of the eigenfrequencies shows the presenc…
Comparison between the fCCZ4 and BSSN formulations of Einstein equations in spherical polar coordinates
Recently, we generalized a covariant and conformal version of the Z4 system of the Einstein equations using a reference metric approach, that we denote as fCCZ4. We successfully implemented and tested this approach in a 1D code that uses spherical coordinates and assumes spherical symmetry, obtaining from one to three orders of magnitude reduction of the Hamiltonian constraint violations with respect to the BSSN formulation in tests involving neutron star spacetimes. In this work, we show preliminary results obtained with the 3D implementation of the fCCZ4 formulation in a fully 3D code using spherical polar coordinates.
Explosion and Final State of an Unstable Reissner-Nordström Black Hole
A Reissner-Nordstr\"om black hole (BH) is superradiantly unstable against spherical perturbations of a charged scalar field, enclosed in a cavity, with frequency lower than a critical value. We use numerical relativity techniques to follow the development of this unstable system -- dubbed a charged BH bomb -- into the non-linear regime, solving the full Einstein--Maxwell--Klein-Gordon equations, in spherical symmetry. We show that: $i)$ the process stops before all the charge is extracted from the BH; $ii)$ the system settles down into a hairy BH: a charged horizon in equilibrium with a scalar field condensate, whose phase is oscillating at the (final) critical frequency. For low scalar fie…
Fully Covariant and Conformal Formulation of the Z4 System Compared to the BSSN Formulation in Spherical Symmetry
We have generalized a covariant and conformal version of the Z4 system of the Einstein equations by adopting a reference metric approach, that we denote as fCCZ4, well suited for curvilinear as well as Cartesian coordinates. We implement this formalism in spherical polar coordinates under the assumption of spherical symmetry using a partially-implicit Runge-Kutta (PIRK) method, without using any regularization scheme, and show that our code can evolve both vacuum and non-vacuum spacetimes without encountering instabilities. We have performed several tests and compared the Hamiltonian constraint violations of the fCCZ4 system, for different choices of certain free parameters, with these of B…
Quasistationary solutions of self-gravitating scalar fields around collapsing stars
Recent work has shown that scalar fields around black holes can form long-lived, quasistationary configurations surviving for cosmological time scales. Scalar fields thus cannot be discarded as viable candidates for dark matter halo models in galaxies around central supermassive black holes (SMBHs). One hypothesized formation scenario of most SMBHs at high redshift is the gravitational collapse of supermassive stars (SMSs) with masses of $\ensuremath{\sim}{10}^{5}\text{ }\text{ }{\mathrm{M}}_{\ensuremath{\bigodot}}$. Any such scalar field configurations must survive the gravitational collapse of a SMS in order to be a viable model of physical reality. To check for the postcollapse survival …
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 dynami…
Fully covariant and conformal formulation of the Z4 system in a reference-metric approach: Comparison with the BSSN formulation in spherical symmetry
We adopt a reference-metric approach to generalize a covariant and conformal version of the Z4 system of the Einstein equations. We refer to the resulting system as ``fully covariant and conformal", or fCCZ4 for short, since it is well suited for curvilinear as well as Cartesian coordinates. We implement this fCCZ4 formalism in spherical polar coordinates under the assumption of spherical symmetry using a partially-implicit Runge-Kutta (PIRK) method and show that our code can evolve both vacuum and non-vacuum spacetimes without encountering instabilities. Our method does not require regularization of the equations to handle coordinate singularities, nor does it depend on constraint-preservi…