Geometrically thick equilibrium tori around a dyonic black hole with quasi-topological electromagnetism

We study geometrically thick and non-self gravitating equilibrium tori orbiting a static spherically symmetric dyonic black hole with quasi-topological electromagnetism. Our results show that the electric and magnetic charges together with the coupling parameter in the quasi-topological electromagnetism lead to a much richer class of equilibrium tori. There are a range of parameters which allow for the existence of double tori. The properties of double equilibrium tori become far richer. There exist transitions between single torus and double tori solutions as we change the specific angular momentum of the fluid. These properties of equilibrium tori could help to understand the dyonic black hole and its thick accretion disk.


I. INTRODUCTION
The presence of the bright zones in the images of the supermassive black holes M87* [1][2][3][4][5][6] and Sgr A* [7] means that a black hole at the centre of a galaxy must be surrounded by an accretion disk.This is because conversion of gravitational energy into heat and radiation in the matter accretion is the most efficient so that the strong electromagnetic radiation emitted by disk can illuminate regions near black holes.In the real astrophysical systems, the matter accretion is a highly complicated dynamic process and its complete description must resort to high precise numerical calculations, such as general relativistic magnetohydrodynamics (GRMHD) simulations [8][9][10][11][12][13].However, in the past few decades, a simple model of geometrically thick and stationary tori orbiting black holes has been attracted a lot of attention.In this model, matter is assumed to be stationarily rotating and without actually approaching the black hole [14][15][16][17][18][19][20][21][22].Moreover, the self-gravity of the fluid body is neglected.Interestingly, such stationarily rotating perfect fluid tori, known as Polish doughnuts, are exact analytical solutions of the relativistic Euler equation [23][24][25].Due to the fluid being in equilibrium, Polish doughnuts are often used as an initial condition for numerical simulations of accretion flows.Additionally, the features and configurations of these geometrically thick equilibrium tori carry a lot of important information on the spacetime in strong field regions because the motion of fluid is very close to the event horizon of black holes.Thus, study of geometrically thick equilibrium tori can offer a potential way of probing spacetime characteristic properties imprinted on the equilibrium tori.
The geometrically thick non-self gravitating equilibrium tori orbiting black holes have been studied for many spacetimes in general relativity and in other alternative theories of gravity [26][27][28][29][30][31][32][33].Recently, equilibrium tori have been investigated for a spherically symmetric black hole in Born-Infeld teleparallel gravity, which show that there is only a single torus as in the Schwarzschild spacetime, but the teleparallel gravity parameter leads to that the size of the torus becomes small [34].The parameterised Rezzolla-Zhidenko black hole is found to have a much richer class of equilibrium tori [35].There exist standard single-torus and non-standard double-tori solutions within the allowed range of parameters and the transitions between single torus and double tori solutions can occur by regulating the specific angular momentum of the fluid.Moreover, the magnetized equilibrium tori around Kerr black hole with scalar hair have been studied with the constant angular momentum model [36] and the nonconstant angular momentum model [37].The properties of the stationary thick tori differed from that in the Kerr case [27] could be used to further constrain the nohair hypothesis with future observations.The geometrically thick equilibrium tori with the constant angular momentum model are investigated for a static spherically symmetric black hole in f (R)-gravity with a Yukawalike modification [38].It is shown that configurations of the tori own the notable differences from those in the usual black hole in the general relativity.Moreover, the magnetised equilibrium tori around a Kerr black hole are also studied with the non-constant specific angular momentum distribution model [39].The configurations of the geometrically thick tori have been also studied in the background of compact object with a quadrupole moment [40] and of a binary black hole system [41].
Recently, a dyonic black hole solution is obtained in the frame of a quasi-topological electromagnetism [42], which is a higher-order extension with the bilinear norm of Maxwell theory.Although under some appropriate conditions, quasi-topological terms have no contribution to Maxwell equation and energy-momentum tensor, they can nontrivially modify the dyonic solutions.The dyonic black hole solution carries both electric and magnetic charges.Although magnetic charges have not been observed in nature, they can arise as topologically nontrivial classical solution in some spontaneously broken gauge theories.It looks very difficult to generate artificially dyonic black holes because it would involve producing or gathering lots of magnetic monopoles and then collapsing them into a black hole.One plausible mechanism is to produce first a large number of monopoles and anti-monopoles in the early universe, then at larger scales one needs large primordial fluctuations that produce black holes [43].The dyonic black hole solution [42] in quasi-topological electromagnetism exhibits some unusual properties.In certain parameter region, the black hole solution possess four horizons and three photon spheres.The energy condition is analyzed for the existence of three photon spheres in the dyonic black hole [44].Especially, this dyonic black hole solution can be used to construct Dyson-like spheres around the black hole [45], at which a massive particle remains at rest with respect to a static asymptotic static observer.These studies shed new light on understanding the dyonic black hole solution with quasi-topological electromagnetism.The main motivation of this paper is to study the geometrically thick equilibrium tori around the dyonic black hole with quasi-topological electromagnetism and to see what new properties of the equilibrium tori in this case.
The paper is organized as follows: In Sec.II, we briefly introduce the dyonic black hole solution and Polish doughnut model of thick accretion disks.In Sec.III, we present properties of equilibrium tori around the dyonic black hole.Finally, we present a summary.

II. DYONIC BLACK HOLE AND POLISH DOUGHNUT MODEL OF THICK ACCRETION DISKS
Let us now to briefly review the dyonic black hole [42], which is a static spherically symmetric black hole solution in the quasi-topological electromagnetism.The action has a form [42] S = 1 16π where the field strength is with where 2 F 1 (a, b, c; x) is the Hypergeometric function.This solution owns three integration constants M , q and p, which are respectively related to the mass, electric and magnetic charges.The corresponding ansatz of the Maxwell field is The electric and magnetic charges are given by [42] where It is obvious that the electric and magnetic charges enter the metric asymmetrically and then the electromagnetic duality breaks down by the factor α 2 .From Eq.( 4), we find that there exist only the radial electric field and the radial magnetic field.The corresponding electric intensity E and magnetic induction intensity B in the local orthonormal frame can be expressed as [46] Fig. (1) shows distributions of electric and magnetic fields for the dyonic black hole with an specific set of FIG.1: Distributions of electric and magnetic fields for the dyonic black hole with an specific set of parameters p = 0.14, q = 1.1, α1 = 1 and α2 = 9. parameters p = 0.14, q = 1.1, α 1 = 1 and α 2 = 9.The magnetic induction intensity monotonically decreases with the distance away from the central black hole, while the electric field intensity first increases and then decreases due to the effect of the magnetic charge parameter.Moreover, the electric and magnetic field energy can be expressed as [42] Since the energy H is nonnegative in the case p = 0, the parameter α 1 is nonnegative.However, when α 1 = 0, the electric potential ϕ(r) = ϕ 0 − qr 3 12α2 p2 is divergent at spatial infinite.Thus, the parameter α 1 is strictly positive.The metric function f (r) can be expanded as where F(r, p, q, α 1 , α 2 ) is a non-zero function.The function f (r) can be reduce to that in the Schwarzschild spacetime only if α 2 = 0 and α 1 = ±iq/p.Thus, for the real non-zero value of α 1 , the solution (2) does not degenerate to the Schwarzschild one.For the convenience, we rescale the parameters as q = q/ √ α 1 , p = p√ α 1 , and α = α 2 /α 2 1 , which makes the metric function f (r) more simple, i.e., f (r To visualize the spacetime (2), we present the embedding diagrams in Fig. 2 where the equatorial slice θ = π/2 at some fix moment in time t =constant is embed into three-dimensional Euclidean space ds 2 = dz 2 + dr 2 + r 2 dϕ 2 with dz = ± 1/f − 1dr.Moreover, the metric ( 2) is asymptotically flat and satisfies the dominant energy condition [42].Comparing with the usual Reissner-Nordström black hole, the black hole (2) could have some interesting spacetime structures, for example, there are four black hole horizons and three photon spheres in certain parameter regions.Moreover, in the case of pure electric charge (i.e., p = 0) or pure magnetic charge (i.e., q = 0), the metric of the dyonic black hole solution (2) can the Reissner-Nordström metric with pure electric charge or pure magnetic charge.However, when electric and magnetic charges emerge simultaneously, the dyonic black hole solution (2) differs from the Reissner-Nordström solution in general relativity.Thus, the further studies of the dyonic black hole solution (2) could deepen the understanding gravity and quasi-topological electromagnetism.In Fig. 3, we also present changes of Hawking temperature of the dyonic black hole (2) with different parameters.It is well known that the black hole with higher Hawking temperature has higher evaporation rate and then has shorter lifetime.Fig. 3 shows that the Hawking temperature of the dyonic black hole (2) increases with the parameter α and then its lifetime decreases for fixed p = 0.14 and q = 0.6.With the increase of p, we can obtain that the lifetime of the black hole first decreases and then increases for fixed q = 0.6 and α = 9.0.With the increase of q, for p = 0.14 and α = 9.0, the lifetime of the black hole first decreases and then increases, finally decreases again.
In order to study the equilibrium non-selfgravitating accretion tori, we must discuss the marginally stable orbit and the marginally bound orbit obtained from the metric (2), which are two essential quantities for determining the thick disk model in a given spacetime.Lets us now focus on the equatorial geodesics with θ = π 2 for a test timelike particle.Combining the two Killing vectors of the black hole background (2), one can obtain two conserved quantities of particle moving along the geodesics, namely, its energy and angular momentum, With these conserved quantities, the motion equation of the timelike particle moving in the equatorial plane can be simplified as where the effective potential is The circular orbit of a timelike particle is determined by V ef f = E 2 and V ′ ef f = 0.The stable circular orbits also meet the second-order derivative of the effective potential V ′′ ef f > 0. The marginally bound orbit is the innermost unstable circular orbit for a timelike particle, which is determined by V ef f = 1 and V ′ ef f = 0. orbit r mb and the outer event horizon of the black hole.Theoretically, the choice of these parameters for a black hole is arbitrary as long as its event horizon exists.The main motivation of choosing such set of parameters is that there are multiple horizons for the black hole and the corresponding complex dynamical properties of test particles could yield richer configurations of thick equilibrium tori.Moreover, for fixed p, q and α, the relationship 2Q e Q m /ℏ = n ( where n is an integer ) can also be ensured by adjusting the factor ), for fixed p = 0.14 and q = 1.1, we find that there exist two marginally stable orbits r + ms and r − ms in certain range of α.With the increase of α, the r + ms monotonically increases, but r − ms first decreases and then increases.With the further increase of α, the inner marginally stable orbit r − ms vanishes and the outer marginally stable orbit r + ms still exists.For the marginally bound orbit, as α < 8.830 for the fixed p = 0.14 and q = 1.1, there exist only a marginally bound orbit and its radius r − mb increases with α.As α > 8.830, there exist another marginally bound orbit with larger radius r + mb .As α increases up to 9.683, the inner marginally bound orbit disappears.For fixed q = 1.1 and α = 9, we find the changes of r − ms , r ± mb and r H with p are similar to their changes with α, but the r + ms decreases with p.For fixed p = 0.14 and α = 9, we find that r + ms , r ± mb and r H decrease with the electric charge parameter q, but r − ms first decreases and then increases.These complex properties of particle's motion orbits are helpful to probe possible double equilibrium tori around the dyonic black hole.Moreover, we find that the orbit radii r ± ms and r ± mb exist simultaneously only within a narrow range of parameters p, q and α.In the Kerr black hole spacetime, it is found that there also exist the outer and inner marginally stable orbits r + ms , r − ms , the outer and inner marginally bound orbits r + mb , r − mb .This is caused by the black hole spin and the test particle co-rotating or counter-rotating around the black hole.These co-rotating and counter-rotating rings yield some equilibrium configurations in ringed accretion disks [47,48].It is differen from those in the static dyonic black hole spacetime (2) where the existences of these characterized orbits and the corresponding double equilibrium tori are caused by gravitational effects from quasi-topological electromagnetism rather than the black hole rotation.Moreover, in the Kerr black hole, the orbit radii r ± ms and r ± mb always exist simultaneously for arbitrary non-zero spin value.This is another difference from that in the case of the dyonic black hole with quasi-topological electromagnetism (2).
Let us now to briefly review the Polish doughnuts model about the equilibrium torus of electrically neutral and stationary rotating fluids around a black hole [14][15][16][17][18][19][20][21][22].For the perfect fluid, its stress energy tensor can be expressed as where P and w are the pressure and the enthalpy density of the fluid, respectively.u µ is the four velocity of the fluid particle, which satisfies the relationship u µ u µ = −1.In the Polish doughnuts model, the fluid in the thick disk is assumed to be a barotropic perfect fluid with positive pressure and its self-gravity is negligible so that the influence of disk on background spacetime is negligible.Moreover, the fluid is assumed to be axisymmetric and stationary, and the rotation of perfect fluid is restricted to be in the azimuthal direction.
With these assumptions, the four velocity and the stress energy tensor of the perfect fluid can be expressed as From the conservation for the perfect fluid ∇ µ T µ ν = 0, one can obtain [14-22] In the background spacetime (2), the redshift factor of the fluid particle u t has a form The quantities Ω and l are respectively the angular velocity and specific angular momentum of the fluid particle For a barotropic fluid, the enthalpy is a function of p and the right side of Eq.( 16) is an exact differential.
According to the von Zeipel theorem [14][15][16][17][18][19][20][21][22], one can get Ω = Ω(l) and then can obtain a solution by integrating Eq.( 16) to obtain The subscript "in" denotes that the quantity is evaluated at at the inner edge of the disk.The potential W ef f is given by which determines the topologies of equipotential surfaces in the disk.As in [49], we assume that a barotropic equation of state takes a form P = Kw κ , where K and κ are constants.Inserting into Eq.( 20), one can easily obtain Thus, once the potenial W ef f and the specific enthalpy h = w/ρ are given, one can get the enthalpy density w and the rest-mass density ρ distribution in the disk.Here, we set h = 1 for simplicity [49].Moreover, we here focus on only the case in which the fluid has a constant specific angular momentum l because dl = 0 and the fluid angular velocity Ω in Eq.( 18) becomes a function only related to spacetime metric, which makes the calculation of equilibrium tori particularly simple.In this simple model, the potential W ef f can be further simplified as W ef f = ln |u t |.Such stationarily rotating perfect-fluid tori with this property are known as "Polish doughnuts".

III. EQUILIBRIUM TORI AROUND THE DYONIC BLACK HOLE
In this section, we will study the electrically neutral equilibrium tori around the dyonic black hole and probe effects of the electric, magnetic charge and coupling parameters on the equilibrium tori.As in the previous discussion, the specific angular momentum l plays an important role in the potential W ef f and has a great impact on the fluid equilibrium tori around black holes [14][15][16][17][18][19][20][21][22]. Figure .(5)shows the changes of the specific angular momentum l of the fluid particles with the rescaled electric, magnetic and coupling parameters (q, p, α) in the dyonic black hole spacetime (2).When the specific angular momentum is equal to that at the marginally bound orbit, l = l mb , one can find that a cusp is located at the marginally closed surface that just extents to infinity.As l > l mb , the outermost equipotential surfaces are connect with the event horizon and the tori are effectively accreting without crossing the cusp.Actually, the maximum of the potential W ef f in this case is larger than the corresponding value at spatial infinity and there is no stable torus in these regions.As the specific angular momentum is less than the value of the marginally stable orbit, l < l ms , there is no equilibrium torus.From Fig. 5, one can find that the curves l(r) for different parameters can be classified as three types: (i) The curve l(r) has only a minimum value l ms and the specific angular momentum at the marginally bound orbit l mb is larger than the minimum value, l mb > l ms , and then there may be only a single tori around the black hole, which is similar to that in the Schwarzschild black hole spacetime.(ii) The curve l(r) has two minimum values, l + ms and l − ms , and l mb is larger than both of the two minimum values, l mb > l ± ms , then there may be double tori around the black hole, which does not appear in the Schwarzschild spacetime.(iii) The curve l(r) has two minimum values, l mb is in between these two minimum values, i.e., l + ms > l mb > l − ms , then there may still be only a torus in this case.In Fig. 6, we show the parameter space for tori solutions around the dyonic black hole (2).The green-shaded region corresponds to the case where the single-torus solutions can be found, the blue-shaded region is the region in which double-tori solutions are possible.FIG.6: Multidimensional parameter space for tori solutions around the dyonic black hole (2).In the blue region there may exist double tori and two cusps.The green region indicates that there is only a single torus and one cusp, and the white region indicates that there is no equilibrium torus.Fig. 7 presents the changes of W ef f in the equatorial plane with the radial coordinate x = 1 − r H /r for different specific angular momentum l.The cusps and the tori centres are respectively located at the positions of the local minima and maxima of W ef f .From Fig. 7, one can find that the cusps move toward the black hole with the increase of l, while the tori centres move outward.In the dyonic black hole spacetime (2) with p = 0.14, q = 1.1 and α = 10, one can find that the potential W ef f has only one local maximum and minimum for certain specific angular momentum l.This means that there exist only a cusp and single equilibrium torus for the fluid with the constant specific angular momentum l around the black hole.In the cases with p = 0.14, q = 1.1 and α = 9 or α = 9.5, we find that there exist two local maxima at r = r ± cusp and two minima at r = r ± max for the potential W ef f in the some range of l.Thus, there are two cusps and two equilibrium torus around the black hole.Moreover, from Fig. 7, we also note that for the dyonic black hole with parameters p = 0.14, q = 1.1 and α = 9, there are double equilibrium tori if 2.838 < l < 2.914, and the centres of the inner and outer tori are located at r = r − max and r = r + max .However, as l decreases down to in the range 2.755 < l < 2.838, the outer torus with the centre r = r + max vanishes, but the inner torus is remained, and then there is a single equilibrium torus around the black hole.For the dyonic black hole with parameters p = 0.14, q = 1.1 and α = 9.5, we find that there also exist double equilibrium tori if 2.977 < l < 3.009.
Similarly, as l decreases down to in the range 2.838 < l < 2.977, there is a single equilibrium torus around the black hole.The difference from the former is that here the inner torus vanishes and the outer torus is remained.These are also shown in Figs. 8 and 9, in which the equipotential surfaces of the effective potential  I.For the cases with double equilibrium tori, if the value of the potential W ef f at the inner cusp is less than one at the outer cusp, i.e., W ef f (r − cusp ) < W ef f (r + cusp ), we find that the equipotential surface of the inner cusp is inside the equipotential surface of the outer one.As shown in the left panel (Model A) in Fig. 9, for the fluid particle with W ef f (r − cusp ) < W ef f < W ef f (r + cusp ), it can be accreted into the black hole if the particle is located inside the inner torus even if its equipotential is less than one at the outer cusp.However, if this fluid particle is located inside the outer torus, it can still keep equilibrium although its potential is more than one at the inner cusp.These are different from that in the single cusp case in which can move along the circular orbit and cannot be accreted into the black hole as shown in the right panel (Model E) in Fig. 9, which is also different from that in the case with single cusp.In Figs. 8 and 9, we also note that the fluid particles in the gray regions can not form equilibrium tori.
Thus, in the double equilibrium tori case, the occurrence of fluid particle accretion depends on the potential values at the two cusps and the position of fluid particle.In Fig. 10, we also present the distributions of mass density and pressure in equilibrium tori by solving Eq.( 21) with the equation of state P = Kw κ .Here, we set κ = 4/3 and w rmax = 1, and then the coefficient K can be calculated by K = W rmax − W in for different models, where W rmax is the value of W ef f at the centre r = r max of the equilibrium torus.Comparing Fig. 10 with Figs. 8 and 9, one can find that the distributions of mass density and pressure in equilibrium tori are similar to the distribution of the effective potential W ef f .Fig. 11 presents the whole set of possible tori solutions appeared in the dyonic black hole spacetime (2) as 0.14259 < p < 0.39124, there is only a single torus in this case.
Finally, we briefly discuss the possible connection between the parameters for double tori in the dyonic black hole spacetime (2) and in the parametric Rezzolla-Zhidenko spacetime studied by Cassing et.al [35].
Making use of the compactified radial coordinate x = 1 − r H /r, one can find that the metric ansatz for the dyonic black hole spacetime (2) can be expressed as where A(x) can be expanded as a similar form in [35,50], i.e., where Comparing with the Rezzolla-Zhidenko black hole case, it is easy to find that the parameter ϵ has the same form.The Rezzolla-Zhidenko expansion coefficients ϵ, a 0 and a 1 are the functions of q, p and α because the event horizon r H also depends on α for the dyonic black hole (2).Moreover, we find that coefficients ϵ, a 0 and a 1 are related by a 1 = a 0 + ϵ in this case.Thus, these results could help to understand equilibrium tori around the dyonic black hole and its thick accretion disk.

IV. SUMMARY
It is useful to study non-self gravitating equilibrium tori around black holes because they have been widely applied to numerical simulation of accretion flows into black holes as a kind of initial conditions.We study geometrically thick non-self gravitating equilibrium tori orbiting a static spherically symmetric dyonic black hole.Within the allowed space of parameters, we find that there exist standard single torus solutions and non-standard double tori solutions.For the single torus solutions, the properties of the equilibrium torus including the cusp and the torus centre are very similar to those in the Schwarzschild black hole spacetime.
For the double tori solutions, the properties of equilibrium tori and the accretion near black hole become far richer.When the value of the potential W ef f at the inner cusp is less than one at the outer cusp, the fluid particle with W ef f (r − cusp ) < W ef f < W ef f (r + cusp ) can be accreted into the black hole as it is inside the inner torus even if its potential is less than one at the outer cusp.However, if the fluid particle is located inside the outer torus, it can still keep equilibrium although its equipotential is more than one at the inner cusp.
When the value of the potential W ef f at the inner cusp is larger than one at the outer cusp, the fluid particles with W ef f (r + cusp ) < W ef f < W ef f (r − cusp ) can move along the circular orbit and cannot be accreted into the black hole.Thus, in the double equilibrium tori case, the occurrence of fluid particle accretion depends on the potential values at the two cusps and the position of fluid particle.More interestingly, the transitions between single torus and double tori solutions may occur by changing the specific angular momentum of the fluid.
Thus, the electric and magnetic charges and the coupling parameter in the quasitopological electromagnetic theory lead to a much richer class of equilibrium tori than in the usual static spherically symmetric cases.The possible observable effects originating from the equilibrium tori around the dyonic black hole could provide unique tests of general relativity.
V. ACKNOWLEDGMENTS

Fig.( 4 )
Fig.(4) shows effects of the parameters α, p and q on the marginally stable orbit r ms , the marginally bound

FIG. 5 :
FIG. 5: Changes of the specific angular momentum l with the circular orbit radius x = 1− r H r for different parameters, where rH is the outermost horizon of the black hole.The black dots indicate the values of l mb .

-0. 4 -0. 3 -0. 2 -FIG. 7 :
FIG.7: Changes of the effective potential W ef f in the equatorial plane with the radial coordinate x for different specific angular momentum l in the dyonic black hole spacetime (2).

FIG. 8 :
FIG. 8: Equipotential surfaces of the effective potential W ef f shown with Cartesian coordinates x = r sin θ cos ϕ and z = r cos θ for black hole parameters p = 0.14 and q = 1.1.The specific angular momentum is set to l = 2.88 in the left panel and is set to l = 3 in the middle and right panels.The parameter α is respectively set to α = 9, 9.5 and 10 from left to right.

FIG. 9 :
FIG. 9: Equipotential surfaces of the effective potential W ef f shown with Cartesian coordinates x = r sin θ cos ϕ and z = r cos θ for black hole parameters p = 0.14, q = 1.1 and α = 9.The specific angular momentum are set to l = 2.88, l = 2.885865 and l = 2.89 in the panels from the left to the right, which respectively correspond to the Models A, D and E. Model α l rH r + cusp vector potential.The coupling parameters α 1 and α 2 are two non-negative constants,