A Note on Geodesics in Hayward Metric

We study timelike and null geodesics in a non-singular black hole metric proposed by Hayward. The metric contains an additional length-scale parameter $\ell$ and approaches the Schwarzschild metric at large radii while approaches a constant at small radii so that the singularity is resolved. We tabulate the various critical values of $\ell$ for timelike and null geodesics: the critical values for the existence of horizon, marginally stable circular orbit and photon sphere. We find the photon sphere exists even if the horizon is absent and two marginally stable circular orbits appear if the photon sphere is absent and a stable circular orbit for photons exists for a certain range of $\ell$. We visualize the image of a black hole and find that blight rings appear even if the photon sphere is absent.


I. INTRODUCTION
A non-singular model of a black hole was proposed by Hayward [1]. As such, it may be a solution of an ultraviolet complete gravity such as [2]. The geodesics in such a non-singular spacetime would illuminate some basics aspects of the spacetime, as the geodesics in the Schwarzschild spacetime do [3].
In this paper, we summarize the results of our study of the properties of geodesics in the geometry described by the Hayward metric in a self-contained manner. Some of the results may be already known (e.g., [4,5]), but the detailed study of both the properties of marginally stable circular orbits and the behaviour of null geodesics are new as far as we know.
In Sec. II, after giving the radii of the horizons of the spacetime, we study the properties of timelike geodesics and null geodesics. Sec. III is devoted to summary. In Appendix A, we summarize the properties of the geodesics in the Reissner-Nordström metric.
We use the units of G = c = 1.

II. GEODESICS IN HAYWARD METRIC
The Hayward metric is given by [1] ds 2 = −F (r)dt 2 + dr 2 F (r) + r 2 dΩ 2 (1) where M is a mass parameter and ℓ is a length-scale parameter. The metric approaches 1 − 2M/r as r → ∞ and approaches unity smoothly as r → 0 and hence is non-singular.
The metric is "minimal" in the sense that it contains the least number of free parameters (ℓ only) with the desired properties (regularity at the center such that F (r) → 1 + O(r 2 ) and Schwarzschild asymptotic behaviour at large radii. In fact, Frolov has recently shown that if F (r) is a rational function of r, the order of the polynomials must be larger than 2 and an F (r) constructed out of polynomials of order 3 contains two free parameters in general [6].
We may compute the effective energy momentum for the metric Eq. (2) via T µν = 1 8π G µν as done in [1]. 1 Then, the energy density ρ = −T t t , the radial pressure p r = T r r and the tangential pressure p T = T θ θ = T φ φ are given by [1] and we find that the weak energy condition is satisfied: ρ > 0, ρ + p r = 0, ρ + p T ≥ 0 and that the strong energy condition is violated for r < ( Henceforth, we normalize the length scale by M and introduce the dimensionless param- and use r = r/M as a dimensionless variable.

A. Horizon
The location of a horizon is determined by F (r) = 0, and horizons (outer horizon r + and inner horizon r − ) exist for 0 ≤ a ≤ a H [1], where The upper bound on a implies the lower bound on M, M ≥ a H ℓ. The implications of this lower bound on the formation and the evaporation of black holes are discussed in [1].

B. Timelike Geodesics
From the spherical symmetry, we may restrict ourselves to the equatorial plane without loss of generality. In terms of two conserved quantities, the energy E = F (r)ṫ and the angular momentum L = r 2φ , the timelike geodesics satisfy the energy equation whereṫ = dt/dτ with τ being the proper time.
A marginally stable circular orbit (MSCO) is determined by the condition V ′ = V ′′ = 0 which reduces to the following equation for r 2 and the innermost such an orbit is called the innermost stable circular orbit (ISCO). Then, 2 If F (r) = 1 − Q n−1 (r)/P n (r), where P n (r) and Q n−1 (r) are polynomials of order n and n − 1, respectively, the degree of the equation becomes 3n − 3.
We find that Eq. (9) has one positive real solution for a < a H and three positive real solutions for a H < a < a I and one positive real solution for a > a I , where a I is given by However, both two of the solutions for a H < a < a P and one solution for a > a P are turned out to have imaginary L and are unphysical, where a P = 25 24 so there is no ISCO for a > a I . After all, there is one ISCO for a < a P and two MSCOs for a P < a < a I and no ISCO for a > a I . In Fig. 1 and Fig.2, we show the radii of ISCO or MSCO (blue curve) and the angular momentum. The angular momentum of the inner MSCO for a H < a < a I is arbitrary large as a → a P .
In Fig. 3, we show the allowed region for the stable circular orbit (gray). The light gray region in the middle is the unstable orbit d 2 V /dr 2 < 0. The lower white region is forbidden because the angular momentum is imaginary or the radius is inside the horizon.
In order to study the properties of circular orbits, we calculate the energy of a particle in circular orbit as a function the radius for a = 0.5, 0.9, 1 as shown in Fig. 4 . The extrema of the energy correspond to MSCOs: for a = 1 the minimum is the outer MSCO and the maximum is the inner MSCO. The region in between two MSCOs (dashed curve) is unstable because d 2 V /dr 2 < 0 there. The outer MSCO can be reached from a particle in circular orbit outside by its emitting energy and angular momentum. On the other hand, since the energy of the inner MSCO is larger than 1 (the energy at large r), the inner MSCO is unbound and may be regarded as an "excited state" that cannot reached from a particle in circular orbit outside. In this sense, only the outer MSCO may be physically important. For a = 0.9, there is also a stable circular orbit with large L and E > 1.

C. Null Geodesics
As in the case of the timelike geodesics, in terms of E and L, the null geodesics satisfy whereṙ = dr/dλ with λ being the affine parameter.
Since L/E is the impact parameter at large r, the (local) maximum of the effective potential V null (r) determines the capture cross section for photons and hence the size of shadow of a black hole: photons with smaller impact parameter will be captured by the black hole. The location of the (local) maximum of V null (r), or the radius of unstable circular orbits of photons (photon sphere [7]), r P , is determined by Similarly to the case of ISCO, we find that the maximum of V null (r) exists if a < a P , where a P is defined in Eq. (12). In Fig. 1 we show the radius of photon sphere (red solid). There also appears a stable circular orbit of photon (red dashed) for a < a P . The appearance of a stable circular photon orbit is inevitable in the presence of the photon sphere (local maximum of V null )) since V null positively diverges as r → 0. 3 In Fig. 5, we show the radius of the shadow b P defined by In order to study the effect of the geometry on the propagation of light rays, we first compute the deflection of light in the Hayward metric. We consider the deflection of light in the equatorial plane. Then from Eq. (13) and L = r 2φ , in terms of the impact parameter b = L/E, up to the turning point of the deflection, φ satisfies F (r 0 )/r 2 0 , the deflection angle ∆φ is then given by We show the result only up to O((M/r 0 ) 6 ), although we can calculate arbitrarily higher orders. The coefficients for a = 0 fully agree with [9], but the sign of the coefficient of (M/r 0 ) 4 including a is different from [4] in which only the numerical values of the coefficient are given and higher-order terms are not given. In terms of the impact parameter b using b −2 = F (r 0 )/r 2 0 , the deflection angle up to O((M/b) 6 ) is given by Again, we find complete agreement with [9] up to this order if a = 0. Since the effect of nonzero a appears only for O((M/b) 4 ) and beyond, it is difficult to detect the effect by the However, the presence of a does affect the existence or nonexistence of the photon sphere and hence affects the behaviour of light rays traveling around the photon sphere. In Fig.   6, we show the trajectories of light rays coming in from upper right for a = 0, 0.5, 0.9, 1. A black circle is the (outer) horizon, red circles are the radii of photon spheres. In Fig. 7, we show the deflection angle as a function of the impact parameter b for a = 0, 0.5, 0.9, 1. If the impact parameter is slightly larger than the radius of a shadow b P , the deflection angle ∆φ is logarithmically divergent as ∆φ ∼ 1 If a = a P , i.e., r P = 25/12, the divergent behavior is changed into Quite similar behaviour is also found for the Reissner-Nordström metric where critical values corresponding to a H , a P appear depending on the value of the charge (AppendixA).
The images of a "black hole" are shown in Fig 8. Photons in black regions never reach the observer due to the presence of a horizon or the deflection angle being ±π/2(modulo 2π).
Photons in white regions travel around the central region more than once, while photons in gray regions reach the observer straightforwardly. The image for a = 0.5 is little different from that of Schwarzschild black hole (a = 0): A blight ring surrounding a black disk (shadow of a black hole) appears. Interestingly, for a = 0.9(a H < a < a P ), the photon sphere exists even though the horizon is absent; a black disk disappears and a black doughnut appears instead, and a ring image persists. More interestingly, even for a = 1(> a P ), the ring persists although the photon sphere is absent as shown in Fig. 8. Therefore, the existence of a blight ring image does not necessarily imply the existence of a photon sphere. Of course, for astrophysical black holes a ≃ 10 −38 (ℓ/ℓ p )(M ⊙ /M) ≪ 1, where ℓ p is the Planck length, and these phenomena would be expected only for Planck-scale primordial black holes.

III. SUMMARY
In this paper, we have studied the timelike and null geodesics in a non-singular black hole geometry proposed by Hayward that involves a parameter ℓ and found several interesting features of the geometry: the existence or non-existence of the horizon, the photon sphere and the ISCO. We also have found that two marginally stable circular orbits appear for a H < ℓ/M < a P although the inner orbit is unbound. The existence of a horizon and/or a photon sphere significantly affects the behaviour of the null geodesics. We have found that a black doughnut appears if the horizon is absent and that blight rings appear even if the photon sphere is absent.
One may think that bright regions in a shadow image are due to the existence of the photon sphere. However, as shown in Fig.8, bright regions also can appear in the spacetime with no photon sphere if the parameter a is slightly larger than a P . For such parameters, while there are no photon sphere, ∆φ can still take a value larger than 2π (see Fig.7). Similar behavior can be found in the Reissner-Nordström metric. These results suggest that such the shadow images are universal behavior for parameters slightly larger than the critical value where photon sphere marginally exists. The parameter ℓ is currently only loosely constrained by the solar-system experiment since the deflection angle at the post-Newtonian order is the same as that of Schwarzschild. The Event Horizon Telescope, a long baseline interferometer experiment, will be able to resolve black holes at horizon scales with the angular resolution of 20µas which corresponds to a size of 3.6M for SgA* at the galactic center [11]. If the shadow of a black hole is observed by such a telescope, we may put an O(1) constraint on ℓ/M. The Reissner-Nordström metric is given by where M, Q are mass and electric charge parameters, respectively. We study the properties of geodesics on this metric and compare the result with the case of the Hayward metric.
Since the electric charge parameter appears in the metric as Q 2 , we need to consider Q ≥ 0 only. If 1 ≥ Q/M, there are two horizons at r = r ± = M 2 ± M 2 − Q 2 . In Fig.9, we have plotted the regions of stable and unstable circular orbits. The qualitative feature is same as the case of the Hayward metric Fig.3. When Q/M takes Q H /M = 1, Q P /M = 3/(2 √ 2) ≃ 1.061, Q I /M = √ 5/2 ≃ 1.118, the horizon becomes extremal, the photon sphere marginally exists, the ISCO marginally exists, respectively. The locations of the photon sphere and the shadow image for the Reissner-Nordström metric are also discussed in [12]. Compared to the case of the Hayward metric (see Fig.7), the qualitative behavior is almost the same. The only difference can be found near the region for small impact parameter b ≃ 0 in overcharged cases Q/M > 1, where the deflection angle ∆φ is almost −π. This is because the null geodesic is reflected by the repulsive force near the origin r = 0. If we consider Q > Q P , there is no photon sphere. However, for a parameter value slightly larger than Q P (see the third figure in Fig.10), ∆φ can still take a value larger than 2π.
This means that their shadow image is similar to the fourth figure in Fig.8.