We discuss a procedure to build new sequences of magnetised, equilibrium tori around Kerr black holes which combines two approaches previously considered in the literature. For simplicity we have assumed that the test-fluid approximation holds, and hence we neglected the self-gravity of the fluid. The models were built assuming a particular form of the angular momentum distribution from which the location and morphology of equipotential surfaces can be computed. This ansatz includes, in particular, the constant angular momentum case originally employed in the construction of thick tori -or Polish doughnuts -and it has already been used to build equilibrium sequences of purely hydrodynamical models. We discuss the properties of the new models and their dependence on the initial parameters. These new sequences can be used as initial data for magnetohydrodynamical evolutions in general relativity.

Testing the true nature of black holes -the no-hair hypothesis -will become increasingly more precise in the next few years as new observational data is collected in both the gravitational wave channel and the electromagnetic channel. In this paper we consider numerically generated spacetimes of Kerr black holes with synchronised scalar hair and build stationary models of magnetized thick disks (or tori) around them. Our approach assumes that the disks are not self-gravitating, they obey a polytropic equation of state, the distribution of their specific angular momentum is constant, and they are marginally stable, i.e. the disks completely fill their Roche lobe. Moreover, contrary to existing approaches in the literature, our models are thermodinamically relativist, as the specific enthalpy of the fluid can adopt values significantly larger than unity. We study the dependence of the morphology and properties of the accretion tori on the type of black hole considered, from purely Kerr black holes with varying degrees of spin parameter, namely from a Schwarzschild black hole to a nearly extremal Kerr case, to Kerr black holes with scalar hair with different ADM mass and horizon angular velocity. Comparisons between the disk properties for both types of black holes are presented. The sequences of magnetized, equilibrium disks models discussed in this study can be used as initial data for numerical relativity codes to investigate their dynamical (non-linear) stability and used in tandem with ray-tracing codes to obtain synthetic images of black holes (i.e. shadows) in astrophysically relevant situations where the light source is provided by an emitting accretion disk.PACS numbers: 95.30. Sf, 04.70.Bw, 04.25.dg 95.30.Qd

We investigate stationary, self-gravitating, magnetised disks (or tori) around black holes. The models are obtained by numerically solving the coupled system of the Einstein equations and the equations of ideal general-relativistic magnetohydrodynamics. The mathematical formulation and numerical aspects of our approach are similar to those reported in previous works modeling stationary self-gravitating perfect-fluid tori, but the inclusion of magnetic fields represents a new ingredient. Following previous studies of purely hydrodynamical configurations, we construct our models assuming Keplerian rotation in the disks and both spinning and spinless black holes. We focus on the case of a toroidal distribution of the magnetic field and build a large set of models corresponding to a wide range of values of the magnetisation parameter, starting with weakly magnetised disks and ending at configurations in which the magnetic pressure dominates over the thermal one. In all our models, the magnetic field affects the equilibrium structure of the torus mainly due to the magnetic pressure. In particular, an increasing contribution of the magnetic field shifts the location of the maximum of the rest-mass density towards inner regions of the disk. The total mass of the system and the angular momentum are affected by the magnetic field in a complex way, that depends on the black hole spin and the location of the inner radius of the disk. The non-linear dynamical stability of the solutions presented in this paper will be reported elsewhere.

No abstract

We build equilibrium solutions of magnetised thick discs around a highly spinning Kerr black hole and evolve these initial data up to a final time of about 100 orbital periods. The numerical simulations reported in this paper solve the general relativistic magnetohydrodynamics equations using the BHAC code and are performed in axisymmetry. Our study assumes nonself-gravitating, polytropic, constant angular momentum discs endowed with a purely toroidal magnetic field. In order to build the initial data we consider three approaches, two of which incorporate the magnetic field in a self-consistent way and a third approach in which the magnetic field is included as a perturbation on to an otherwise purely hydrodynamical solution. To test the dependence of the evolution on the initial data, we explore four representative values of the magnetisation parameter spanning from almost hydrodynamical discs to very strongly magnetised tori. The initial data are perturbed to allow for mass and angular momentum accretion on to the black hole. Notable differences are found in the long-term evolutions of the initial data. In particular, our study reveals that highly magnetised discs are unstable, and hence prone to be fully accreted and expelled, unless the magnetic field is incorporated into the initial data in a self-consistent way.