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 non-self-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.

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.

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 assume that the test-fluid approximation holds, and hence we neglect the self-gravity of the fluid. The models are 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.