Geodesic structure of a rotating regular black hole

We examine the dynamics of particles around a rotating regular black hole. In particular we focus on the effects of the characteristic length parameter of the spinning black hole on the motion of the particles by solving the equation of orbital motion. We have found that there is a fourth constant of motion that determines the dynamics of orbits out the equatorial plane similar as in the Kerr black hole. Through detailed analyses of the corresponding effective potentials for massive particles the possible orbits are numerically simulated. A comparison with the trajectories in a Kerr spacetime shows that the differences appear when the black holes rotate slowly for large values of the characteristic length parameter.


I. INTRODUCTION
Black holes are perhaps the most fascinating and mysterious objects in the Universe.Although classical black holes are well described by General Relativity they present a singularity where the theory is unable to give a complete description.Singularities thus, are regarded as indicating the breakdown of the theory.It has been argued that quantum modifications of general relativity may correct the singular behavior of black hole solutions and a complete description of quantum gravity will circumvent these problems, however, such formulation is still missing.
Although singularities are hidden by event horizons thus being harmless for the description of physics outside black holes, studying regular black holes is a problem of interest.Several attempts have been made to eliminate the singularity on black holes by replacing it at the Planck scale curvature by the de Sitter geometry.For instance, effective theories introduce a minimal cut-off length, close to the Planck length, yielding a singularity-free space time.There are also some classical approaches that cure the singularity of the spacetime by coupling gravity to an external form of matter, sometimes modeled by some form of nonlinear electrodynamics.The first regular black hole solution was reported by Bardeen in [1] and since then, many metrics that are spherically symmetric, static, asymptotically flat and have regular centers have been found [2][3][4][5][6][7][8].Ayon-Beato and García [9] and later Burinskii and Hilderbrandt [10] showed that the Bardeen solution represents a black hole with a magnetic monopole in a nonlinear electrodynamics.
Among the most noticeable models of regular black holes was the spherically symmetric spacetime proposed by Hayward in 2006 [11].Fan and Wang showed in [12,13] that this spacetime results from Einstein's equations in the presence of a magnetic charge in a nonlinear electromagnetic field as a source.Later, rotating versions of the Hayward metric were introduced by Bambi and Modesto [14].The metrics derived by Bambi and Modesto were constructed using the Newman-Janis algorithm commonly used to produce axisymmetric spacetimes from a spherically symmetric solution, see [15] for a detailed derivation.The resulting metric still represents a regular black hole, nevertheless it was shown in [16][17][18] that the matter content needed to produce the geometry is only an approximation to its corresponding Hayward nonrotating magnetic monopole source.
The interior of a black hole is a hidden region by definition.The surface that connects the interior to the exterior is the horizon.Thus, a way to peer into the interior is to study the near horizon region.While the study of the physical properties of black holes represents a wide field, the most basic study to be performed is a study of the geodesics in these spacetimes.The motion of test particles around black holes has been studied extensively over the years.For most of the common black holes, e.g.Schwarzschild, Reissner-Nordström, or Kerr, analytic solutions to the equations of motion can be given in terms of elliptic functions [19,20].For a non rotating Hayward black hole a geodesic study in the equatorial plane was made by Abbas and Sabiullah in [21].The geodesic equations for more general black holes become complicated and solutions to the orbits in a closed form are unavailable.Still, the symmetries of the spacetime may be used to simplify the geodesic equations and provide relevant information about the motion of particles around the black hole.The purpose of the present work is to analyze the motion of massive test particles in the vicinity of a rotating regular black hole.We use effective potential methods to characterize the motion and present numerical solution for the equations of motion out of the equatorial plane.These results are of particular interest if one considers for instance accretion processes.Usually, studies of accretion focus only on particles moving in the equatorial plane, because this is the easiest case, but particles initially moving around some particular orbit may be perturbed by a kick along the θ direction and give rise to a much richer orbits.Furthermore, studies of geodesics out of the equatorial plane are important from the gravitational wave point of view.In a Kerr spacetime for instance, the evolution of the trajectories is entirely determined by the radiated energy and angular momentum.
One can expect that astrophysical black holes are rotating.A progenitor massive star has non vanishing angular momentum.Even if some part of angular momentum is lost during the the black hole formation process, the resulting black hole would be rotating.Up to now the Kerr solution of the Einstein's equation has played a major role in the description of astrophysical black holes [22].The recent observations with the Event Horizon Telescope have provided us with an image of the shadow of the super-massive black hole M87 which is consistent with the shadow of a Kerr black hole.The study of alternatives to the Kerr solution as the one presented here and the study of geodesics around regular rotating black holes, is a useful test bed for exploring minimal departures from classical black hole geometries, additionally such studies may contribute to a better understanding of the data coming in the future [23].
The paper is organized as follows: In section II we present some basic properties of the rotating Hayward metric proposed by Bambi and Modesto [14], and present a characterization of the horizons and the ergospheres.In section III we derive the equations of motion using the Hamilton Jacobi approach paying particular attention on the conserved quantities and in the separability of the equations.We also present some properties of radial and polar motions using the effective potentials and highlight some features of the radial and angular motion.In section IV we show the numerical solutions of the equations of motion and discuss the trajectories of the particles in the configuration space.Finally, in section V we discuss the results and present some concluding remarks.Through the paper we use geometrical units with G = c = 1.

II. SPACETIME PROPERTIES
The spherically symmetric metric proposed by Hayward [11] is given by with where m is the Arnowitt-Desser-Misner mass and the parameter , is of the order of the Planck length, for = 0 the metric (2.1) reduces to Schwarzschild.The Hayward spacetime is asymptotically flat and in the limit r → 0 (2.1) behaves as a de Sitter metric with a cosmological constant Λ = 3 2 [11].Fan and Wang have shown that the Hayward metric can be obtained as a solution of Einstein's equations coupled with a nonlinear electrodynamics with a magnetic monopole as a source [13].Depending on the relative values of m and the metric (2.1) can have one, two or zero horizons.
By using the Newman-Janis algorithm Bambi and Modesto [14] obtained a family of possible generalizations of Hayward metric to include rotation.In this work we will focus on a metric that has the same form as the Kerr metric with a varying mass.
where the functions ∆, Σ and M are defined as: (2.4) The length parameter g, is of the order of the Planck length and is a measure of the deviation of the Kerr spacetime.The components of the metric inverse, which will become useful later, are ) The spin parameter a is related with the total angular momentum J, by a = J/m.For g = 0 equations (2.3) and (2.5) reduce to the Kerr solution in Boyer-Lindquist coordinates.As described in [14] the spacetime metric (2.3) is regular at r = 0 for g = 0.A rigorous analysis about the regularity of this spacetime can be found in Torres and Fayos [25].In particular, it has been argued that an extension of the spacetime with values r < 0 does present a singularity [26].To avoid such singularity we focus in the region covered by r ≥ 0.

Horizons
The horizons are defined by the relation g rr → ∞ which is equivalent to set ∆ = r 2 − 2M (r)r + a 2 = 0 . (2.6) This equation allows for real solutions for the radius depending on the values of the parameters a and g.The outermost radius determines the event horizon location.In the limit g → 0 there are two horizons for values 0 < a < m.In the extreme case, a = m the two horizons coincide.For the nonrotating case, the equation defining the location of the horizons becomes which implies a critical mass m c = 3 √ 3 /4 and a critical radius r c = √ 3 , such that the horizons are degenerate at r = r c when m = m c .When m < m c there are two horizons whereas there are no horizons when m > m c [11].When a = 0 Eq.(2.6) becomes a polynomial of fifth order in r and it becomes necessary to find the roots numerically.Fig. 1 shows the behavior of ∆ as a function of r for several values of the spin parameter a.The intersection with the horizontal axis determine the position of the horizons, and depending on the relative values of a and g the spacetime posses, two, one or none horizons.The left plot in the first row of Fig. 1 displays the behaviour of ∆ for the non rotating case a = 0.For values g < 1.058 the spacetime has two horizons, for the threshold value g = 1.058 the spacetime has one horizon and no horizons are present for g > 1.058.For positive values of the spin a the threshold value of g changes as shown in the rest of the plots in Fig. 1.For instance, in the right plot of the third row for a spacetime with a = 0.9, the threshold value is g = 0.492.It can be seen that, as the value of a approaches to 1, the threshold value of g decreases.Further information can be obtained by keeping the value of g fixed.Fig. 2 shows a surface plot of ∆ as a function of a and r for g = 0.9.From the figure one can see that, as the value of a increases from zero to a critical value a c < 1, two horizons are present.For values of a between a c and one there are no horizons.Remarkably and unlike the Kerr case, the degeneration of the horizon (the extreme black hole) happens for values of a less than one.

Ergosphere
There is another important surface of rotating black holes.The static limit.When a particle crosses the static limit the nature of the particle changes; a timelike geodesic becomes spacelike and a spacelike geodesic becomes timelike.The static limit surface is defined by the equation The ergosphere is a region located between the event horizon and the static limit surface.In this region particles can extract energy from the black hole via the Penrose process [27,28].A throughout study of the ergoregion of a rotating Hayward BH was performed in [29].The authors conclude that the ergoregion is enlarged compared to Kerr as the value of the parameter g increases.Fig. 3 illustrates the behavior of the shape and extension of the ergosphere of the Hayward rotating black hole for various values of the rotation parameter a and the parameter g.In order to plot these surfaces, we have used the spheroidal-like coordinates (2.11) In Fig. 3 we show a sequence of ergoregions and horizons of the rotating Hayward black hole varying the parameter g.The labels r + and r − at the surfaces in Fig. 3 denote the outer and inner horizons, while E + and E − represent the outer and inner ergospheres.The plots correspond to values of a = 0.9 and g running from g = 0 to g = 0.49.Note that, as the value of g increases the two horizons merge into one.
In order to give a quantitative measure of the effect of g on the ergospheres and horizons we compare the ratio between the length of the geodesic between the north and south pole and the equatorial circle ε(g; R) := lp le , where and dφ .
(2.13) FIG.2: Surface ∆ as a function of r and a for g = 0.9.The intersection of the surface with the plane z = 0 for a given value of a determines the position of the horizons.

FIG. 3:
The ergoregion is delimited by the limit surface defined by equation (2.8).The Kerr black hole corresponds to g = 0.The rotation parameter used in the plots is a = 0.9.The effect of the parameter g is more dramatic in the inner horizon r − and in the inner ergosphere E − .
The surface R will take the values r + , r − , E + E − respectively.
Fig. 4 shows the ratio ε as a function of g for the horizons r ± and ergospheres E ± .One can see that ε(g; r + ) and ε(g; E + ) decreases as g grows, while ε(g; r − ) and ε(g; E − ) are increasing functions.The Kerr case in this plot would correspond to constant horizontal lines at the value ε(0; R).What we can infer from the ratio ε(g; R) is that the oblate surfaces r + and E + get extended at the equator while the surfaces r − and E − become narrower and extended at the poles in accordance with the behavior shown in Fig. 3.

FIG. 4:
The ratio between the length of the pole to pole geodesic and the equatorial circle ε(g; R) for a rotation parameter a = 0.9 for the outer(inner) horizons and outer(inner) ergospheres.

III. GEODESIC MOTION OF TIMELIKE PARTICLES
In order to give a general outlook of the geodesic motion of massive particles in the rotating Hayward spacetime, we use three different but otherwise equivalent formalisms each one suitable to show a specific aspect of the motion.First, we appeal to the Hamilton-Jacobi formalism to show that a fourth integral of motion associated to the Carter constant exist and the problem is completely integrable.Then, once we have proven the motion is separable into radial and polar parts we describe the motion by means of effective potentials.Finally in order to integrate numerically the equations, we use the hamiltonian formalism.

Hamilton Jacobi approach
In the stationary and axially symmetric Kerr spacetime the geodesic equations are completely integrable.There are two obvious conserved quantities given by the symmetries of the spacetime; the energy and the azimuthal angular momentum.Furthermore, the condition p µ p µ = −µ 2 gives another integral of motion and a fourth integral was discovered by Carter [30].This fourth integral can be obtained by separating variables in the Hamilton Jacobi equation of motion.
Let us consider the geodesic motion of a massive test particle in the background of the metric (2.3).The Hamilton-Jacobi equation for this particle has the form: where τ is the proper time.Because the metric does not depend on t nor φ, we can introduce two integrals of motion along the trajectories, namely the energy E and the axial angular momentum L z as: As in the Kerr black hole, the Hamilton Jacobi function S for a geodesic in the Hayward rotating black hole can be written in the following separated form where we have set the mass of the particle µ = 1.Substituting the components of the inverse metric (2.5), the ansatz (3.4) and the corresponding momenta into Eq.(3.1) we obtain: Simplifying the equation (3.5), introducing a separation constant k representing an additional constant of motion, and the Carter's constant through the relation 2 , we arrive at the following differential equations for S r (r) and S θ (θ): where The resulting expressions for the radial and angular momentum are Formally, the solution of the Hamilton Jacobi equation is the principal function S of the form: While this solves the problem theoretically, it is difficult to understand the orbit structure analytically and numerical solutions become necessary.
In order to get a relation between coordinates and momenta we use the lagrangian: and thus where dot denote derivation with respect to the proper time τ .Explicitly, by the spacetime under study One can thus, by using Eqs.(3.9) and (3.13), write the equations of the geodesic motion in the first order form in the lagrangian formalism as However, equations (3.14) and (3.15) present terms with square roots and it is known that these terms causes difficulties in numerical solutions due to the change of signs in the turning points, see for instance [31].For this reason becomes preferable to reformulate the equations of motion using other equivalent approaches.
Before going deeper into the solutions of the equations of motion one can split the analysis considering the motion in r and θ separately and get some insight about the general motion.
In the following we show that equations (3.14), (3.15) can be used to determine the properties of the motion by using effective potentials.The effective potential method is an advantageous tool used to understand and characterize the motion of particles.Most importantly, this method allows us to infer the properties of the motion without any integration avoiding the numerical difficulties.

A. General features of radial motion
The radial motion for timelike particles moving along geodesics in the equatorial plane θ = π/2, on a non-rotating Hayward black hole was described by [32][33][34].The rotating case in the equatorial plane was also considered in [29].Here, we summarize some of the results derived by [29] in the equatorial plane by means of the effective potential method.
Taking the square of Eq. (3.14) and setting Q = 0, the right hand side becomes a polynomial of second order in E. Hence, after some algebra one can rewrite the equation of motion for r as where the potential functions V ± in terms of the metric coefficients are Substituting the expressions (2.5) one gets: This expression reduces to Kerr in the limit for M = cte.[35,36].We can discuss the qualitative features of the motion of massive particles by plotting V ±.Physically acceptable orbits correspond to particles with energy E, greater than V + .
In Fig. 5 are shown three different potential functions V + , for spacetimes with no horizon, one horizon and two horizons obtained with three different values of the scale length g = 0.7, 0.49, 0.2 respectively, with fixed angular momentum L z = 3.0 and a = 0.9.The vertical dotted line at r = 1.13 indicates the single event horizon for g = 0.49 while the line at r = 1.42 indicates the location of the external event horizon for g = 0.2.For the last case there is a relative maximum in the potential indicating a (unstable) circular orbit.Notice that for a fixed value of a, the external event horizon moves to smaller values of r as the value of g increases.For a = 0 we recover the Hayward non-rotating space time [11].The potential functions for spacetimes with g = 0.49 and g = 0.7 are quite different from the case with g = 0.2 In the former cases a particle with energy E 1 will reach a minimum radius and then escape towards infinity, while a particle with the same energy moving in a black hole with g = 0.2 will fall.A particle with energy E 2 cannot traverse the potential barrier for black holes with g = 0.49 and g = 0.7 but for a black hole with g = 0.2 the particle will follow a circular orbit.This different behaviour in the dynamics of the particles may be used to constraint the values of the parameter g.In Fig. 6 it is shown the potential for different values of L z with FIG.5: Potential function V + , for different values of g with a = 0.9 and L z = 2.5.For the potential with g = 0.7 there is no horizon and as the values of g decreases, the external horizon tends to the value given by the Kerr limit r g = 0.2 and a = 0.9.For this case, the spacetime posses two horizons and the values of the angular momentum are representative to show the behavior of the potential in the near horizon region.The vertical dotted line indicates the position of the external horizon.As the angular momentum of the particle decreases, the potential develops a local maximum and a local minimum allowing unstable and stable circular orbits respectively.The particles with energy between the local maximum and minimum of V + have a bounded orbit whose radius lies between the turning points.
In the limiting case of L z = 0 (radial motion) the potential V + in Eq. (3.20) becomes monotone without a maximum or minimum.From its derivative, the equation V + = 0 has no real solutions for any value of the parameters, thus when L z = 0 the particles will always fall into the black hole.In Fig. 7  the rotating parameter a.A general feature is that for large r the curves are asymptotic to zero and, like in the Kerr black hole, the effect of rotation is only noticeable for small r.Furthermore, the effect of g is also relevant for small values of r.
A salient feature of the potential function V + , is that for aL z < 0, it takes negative values as we approach the horizon.In the region where the potential is negative, the energy of the particles E may be also negative and extraction FIG.7: Potential function V + for different values of a and g = 0.2.The angular momentum of the particle is L z = 2.5.
of energy from the black hole is possible via the Penrose process.In Fig. 8 we plot V + for a negative value of the angular momentum L z = −2.5 and for the rotation parameter a = 0.6, a = 0.9 and the limit a = 0.The vertical dashed lines are the positions of the event horizons.In the non rotating limit, the potential function is always positive, whereas for a = 0.6 and a = 0.9 there is a region close the horizon where motion with negative energy is allowed.

B. General features of polar motion
As compared to the orbits in the equatorial plane, polar motion has a richer variety of orbits.A qualitative description of the motion can be performed analyzing the equation (3.15).Let us introduce a new independent variable u := cos θ.The square of the right hand side of equation (3.15) in terms of u becomes a quartic polynomial where From (3.21) we conclude that the motion is only possible for f (u) ≥ 0. Furthermore, evaluating f (u) and its derivatives f (u) and f (u), at the extreme values one infers that the particle can reach the axis u 2 = 1 (θ = 0 or θ = π) if and only if L z = 0. Following the analysis performed by Carter for the Kerr black hole [30] we found the types of θ motion can be classified according to the sign of Q.
• If Q < 0 then f (0) < 0 and f (1) ≤ 0, then the case of interest occurs when In general the motion is oscillatory between 0 < u 1 ≤ u ≤ u 2 where u 1 and u 2 are the two positive zeroes of f (u).In this case however, the particle never crosses the equatorial plane.
• If Q > 0 the particle moves in an oscillatory way crossing the equatorial plane with the angle θ lying in the range θ 0 ≤ θ ≤ π − θ 0 where cos θ 0 = u 0 and u 0 is the unique zero of f (u) in the range 0 < u ≤ 1. Aditionally, if L z = 0 and Q − a 2 (E 2 − 1) > 0, the particle will remain on the axis θ = 0.
In Fig. 9 we plot f (u) with Q = −2.0,Q = 0 and Q = 2.0.In this plot, we can identify the possible trajectories in the parameter space according to the behavior of f (u).The motion is only possible in the non shaded region.
FIG. 9: Motion of particles is permitted in the region where f (u) ≥ 0. For Q > 0 the values that u can take are between u = 0 and u 0 , the movement is bounded and the particle can cross the equatorial plane.For Q = 0 the particle remains in the equatorial plane u = 0.For Q < 0 the trajectory is bound between u 1 and u 2 but the particle cannot cross the equatorial plane.

IV. TRAJECTORIES OF PARTICLES IN THE CONFIGURATION SPACE
In this section, we consider the trajectories of timelike particles in a Hayward rotating black hole in more detail.The dynamics of test particles is governed by the geodesic equations (3.14), (3.15), (3.16), (3.17).However, as described in [37] for the Kerr spacetime, the square terms in ṙ and θ in the equations of motion cause problems in a numerical implementation when determining the turning points because there is a change of sign in ṙ and θ in those points.To deal with this problem we employ the Hamilton formulation to find the trajectories in the configuration space.The hamiltonian and the equations of motion are As before, the energy and the projection along the z axis of the angular momentum are constant of motion p t = −E and p φ = L z .The equation for ṗr is Using the normalization of the four momentum p µ p µ = −1 (recalling we have set the mass of the particles µ = 1) we can solve for Substituting this expression in (4.4) we get ṗr = 1 2 g θθ .( This equation involves only radial derivatives of the metric coefficients and constants of motion.The equation for ṗθ becomes Finally, in order to close the system, the equations for the coordinates in terms of the momenta and constants of motion are: θ = g θθ p θ , (4.10) φ = g φφ L z − g tφ E . (4.11) Note that since these equations do not contain square roots, they constitute a smoothly differentiable system of equations even at turning points and they can be integrated directly.We solve equations (4.4), (4.7), (4.8), (4.9), (4.10) and (4.11) numerically directly in the proper time using the Mathematica numerical integrator NDSolve with a variable step-size Runge-Kutta integrator which, at each level takes a sequence of steps in the independent variable and uses an adaptive procedure to determine the size of these steps.NDSolve reduces the size step to track the solution accurately [38].

Motion in the equatorial plane
To start, let us focus on the radial motion of particles moving in the potential described in Fig. 5.We solve the equations of motion numerically setting up the initial conditions in such a way that the particles have the energies E 1 and E 2 that correspond to the horizontal lines in Fig. 5.In Fig. 10, we show the trajectory of a particle with energy E 1 = 1.05.As one can infer from the potential, the trajectory is unbounded for a spacetime with g = 0.49 and a spacetime with g = 0.7.Particles with energies E > E 1 in these spacetimes will reach a minimum radius and then escape to infinity.However, for a spacetime with g = 0.2 the particle with the same energy will fall into the black hole.In Fig. 11, we plot the trajectories of test particles for a rotating Hayward black hole with angular momentum a = 0.9 and for g = 0.2, 0.49 and 0.7.The angular momentum of the particles is fixed at L z = 2.5 in all cases.The particles have energy E 2 = 0.961, that correspond to the horizontal line E 2 in Fig. 5.In Fig. 11 (left panel) for g = 0.2 the test particle moves on a unstable circular orbit with radius r = 1.98, any small disturbance will make the test particle out of their original orbit.For the cases with g = 0.49 and g = 0.7 (middle and right panels) the particle with the same energy and angular momentum moves on non circular bound orbits.For these orbits the precession direction is clockwise and the precession speed is faster than on the Kerr black hole with the same angular momentum.The orbits presented serve to exemplify the effect of the parameter g on the trajectories of particles; for small values of g a particle will fall into the horizon whereas a particle for larger values of g, with the same initial FIG.10: Motion of a particle with energy E 1 = 1.04 for a spacetime with g = 0.49 and g = 0.7.In the former case the spacetime has one horizon and in the later there is no horizon.
FIG. 11: Orbits for the radial potential with different values of g.The left panel corresponds to a spacetime with a single horizon, and the particle follows a (unstable) circular orbit.The central panel correspond to a spacetime with a single horizon and the trajectory of the particle is bound.The right panel shows the orbit of a particle moving on a spacetime with no horizons.conditions, will escape to infinity after encounter the potential barrier.In Fig. 12 we plot the trajectories of particles on a rotating Hayward black hole with parameters a = 0.9 and g = 0.2.The figures illustrate the trajectories for particles with angular momentum L z = 2.5, 1, 5, 0.5.Like in the Kerr black hole, as the angular momentum of the particles decreases, the potential barrier vanishes yielding to trajectories that eventually fall into the black hole.The trajectories correspond to particles with constant energy represented by the horizontal lines in Fig. 6.The dots on the intersection of the potential function and the energy in Fig. 6 denote the initial position of the particles r = 2.5 in all cases.The trajectory for the particle with L z = 2.5 is the only bound orbit with radius between r = 2.5 and r = 10.2.For particles with L z = 1.5 and L z = 0.5 the trajectories are unbounded and the particles will fall into the black hole as can be infer from the potential function in Fig. 6.
FIG. 12: Orbits for particles with different values of Lz.For the parameters a = 0.9 and g = 0.2 the metric corresponds to a black hole with two horizons.
In Fig. 13 we present the orbits of particles with angular momentum L z = 2.5, initial position r 0 = 2.5 which corresponds to the dots in the figure.The constant energy is indicated for each particle with an horizontal line in Fig. 7.The parameters of the black holes are g = 0.2 and a = 0.9, 0.6, 0. As illustrated in Fig. 7, for the value L z = 2.5 the potential barrier vanishes as a tends to zero and the local maximum and minimum of the potential tend to merge and eventually fade out for a = 0.In this process, bound orbits disappear leaving only trajectories that fall into the black hole.For a = 0.6 and a = 0 the particles fall into the black hole while the particle moving around the black hole with a = 0.9 remains orbiting between two radii r = 2.5 and r = 10.2.FIG.13: Trajectories of particles orbiting a black hole with different values of a.In each case the black hole have two event horizons.For a = 0 and a = 0.6 the particle falls into the black hole and for a = 0.9 the trajectory is bounded between two radii.
In Fig. 14 we show three plots for bound trajectories of particles moving around a black hole g = 0, 0.6 and g = 0.96 the last case corresponds to a limit value in which the spacetime has one single horizon.The spin of the black hole is a = 0.4.The initial position is r 0 = 7.0, the magnitude of the initial velocity is v 0 = 0.45 and the angular momentum of the particles is L z = 3.5.The dots in Fig. 14 corresponds to the initial position r 0 .For these cases, the values of a and L z were chosen to show noticeable changes in the trajectories as g varies.The time of evolution is τ = 2780.As g increases the difference in the trajectories becomes more noticeable.For these cases, the change in the orbits with respect to Kerr black holes may become relevant for instance in the gravitational wave emission [39,40].For instance when in a binary system where binaries are assumed to consist of a Kerr massive black hole and a small compact star which is taken to be a point particle.The differences in the trajectories of massive particles in the equatorial plane, between the Kerr black hole and the rotating Hayward black hole, are the change in the angle of precesion.As g increases, the precesion increases.This difference is more noticeable for small values of the spin parameter a, for nearly extreme black holes (a ∼ 1) the differences become negligible.Also we have found that the parameter g may play an important role in the dynamics of massive particles, for a specific value of a and the angular momentum, a particle would orbit the black hole or fall into it depending exclusively on the parameter g.E.g. if a particle falls into a Kerr black hole, we can choose a value of g in the rotating Hayward black hole such that the particle, with the same conditions, will not fall into the black hole.FIG.14: Orbits of particles moving in the equatorial plane around black holes with spin a = 0.4.The particles have angular momentum L z = 3.5.The motion in the equatorial plane corresponds to a Carter's constant Q = 0.

Motion out the equatorial plane
Some studies of geodesics in the equatorial plane have been done for regular black holes [21].Nevertheless, there are not such studies for the motion out the equatorial plane.It is important to make a fully study of the orbits out of the plane in order to make accurate comparisons with the Kerr black hole and the data obtained from astronomical observations.Here we calculate numerically the trajectories of particles moving on the rotating Hayward spacetime out the equatorial plane.As described above the orbits can be classified according the sign of the Carter constant Q.
In Fig. 15 we plot the trajectories for a particle with a negative Carter's constant Q = −0.2.The two plots of the left panel corresponds to the Kerr spacetime with g = 0, the plots in middle panel is a spacetime with two horizons g = 0.2, and the right panel is a spacetime with a single horizon determined by g = 0.49.In all cases the black hole has spin a = 0.9.For all particles, the initial velocity is 0.96 and the initial position r 0 = 5.0 and is represented with a dot in the plots.The time evolution in units of M was τ = 160.The first row corresponds to a projection of the orbit in the plane x-y.The black hole spin axis is z and is rotating counterclockwise.The second row corresponds to a projection in the z-x plane.The motion of the particle occurs between two angles θ 1 and θ 2 , which correspond to u 1 and u 2 of the function potential f (u) in section III.These angles change according to the value of g, as g increases the difference between the two angles becomes slightly small.For g = 0 the value of θ lies between θ 1 = 0.316π and θ 2 = 0.233π.In Fig. 16, the particles initially satisfy the condition Q > 0, the initial velocity is v 0 = 0.52 and the initial position is r 0 = 6.For the black hole with g = 0, the particle is restricted to move between the angles θ 1 = 0.280π and θ 2 = 0.719π.For the cases illustrated, g = 0 is the Kerr black hole and g = 0.6, g = 0.8 are rotating Hayward black holes with two and one horizon.For rapidly rotating black holes, the effect of g on the trajectories is almost negligible but for black holes with low spin, the trajectories change with respect to Kerr.As the value of g increases the effect is more noticeable.However, there is a limit on g in which the horizons of the black hole vanish yielding a spacetime with no horizons.

V. DISCUSSION AND CONCLUDING REMARKS
In this work we have described the properties of the horizons and ergospheres for a rotating black hole geometry proposed by Bambi and Modesto in Ref. [14].The proposal was made as a generalization of the Hayward geometry [11] and is characterized by a length scale g and the spin a.We have studied the timelike geodesics around such rotating black holes.We discussed possible types of orbits in these spacetimes using effective potentials.We examined its properties paying particular attention to the changes due to the parameter g the spin of the black hole a and the angular momentum of the particles L z .We have described the properties of the geodesics using the Hamilton-Jacobi approach showing that for the rotating Hayward black hole proposed by Bambi and Modesto, there is a fourth constant of motion similar to the Carter constant in the Kerr black hole.We showed that the motion of particles out of the equatorial plane is well characterized in terms of this constant.Furthermore, we found the trajectories of the particles integrating numerically the equations of motion considering particles lying in and out the equatorial plane.We compare these trajectories with those of the Kerr black hole with the same mass and spin.We found slight differences in the trajectories in both spacetimes.In particular, for slowly rotating black holes the differences become noticeable for large values of the parameter g.These differences may become relevant in studies of gravitational wave FIG.16: Two views of the trajectories of particles moving around rotating black holes.The spin of the black holes is a = 0.6.The spacetime with g = 0.836 corresponds to a black hole with a single horizon.The initial data of the particles were set such that the Carter's constant Q > 0.
emission of compact objects orbiting massive black holes.It is the general belief that astrophysical black holes are of the Kerr type and the astrophysical importance of studies of regular black holes may have limited applications.However, the study of orbits around rotating black holes is important from conceptual and theoretical points of view since it is a medium to seek for differences among the black holes specially in the near horizon region where the parameters of the theory are important.Little is known about the orbits in compact rotating objects other than Kerr.If such objects were to represent serious candidates for compact astrophysical objects, the presence of representative signatures would be highly valuable for their identification.

FIG. 1 :
FIG. 1: Possible spacetime structures for rotating regular black holes depend on the values of g and a.The zeroes of ∆ determines the presence of two, one or none horizons.

FIG. 6 :
FIG. 6: Potential function V + for different values of L z .The black hole parameters are g = 0.2 and a = 0.9.The vertical line refers to the external event horizon.

FIG. 15 :
FIG. 15: Different views of the trajectories of particles moving in spacetimes with different g.The spin of the black hole is a = 0.9 and the particles satisfy the condition Q < 0. The figures in the right column correspond to a spacetime with a single horizon