GEMS embeddings of Hayward regular black holes in massless and massive gravities

After finding a solution for the Hayward regular black hole (HRBH) in massive gravity, we embed the (3+1)-dimensional HRBHs both in massless and in massive gravities into (5+2)- and (6+3)-dimensional Minkowski spacetimes, respectively. Here, massive gravity denotes that a graviton acquires a mass holographically by broken momentum conservation in the HRBH. The original HRBH has no holographically added gravitons, which we call massless. Making use of newly found embedding coordinates, we obtain desired Unruh temperatures and compare them with the Hawking and local fiducial temperatures, showing that the Unruh effect for a uniformly accelerated observer in a higher dimensional flat spacetime is equal to the Hawking effect for a fiducial observer in a black hole spacetime. We also obtain freely falling temperatures of the HRBHs in massless and massive gravities seen by freely falling observers, which remain finite even at the event horizons while becoming the Hawking temperatures in asymptotic infinity.


I. INTRODUCTION
Black holes are among the most mysterious and fascinating objects in our universe, observationally as well as theoretically.Theoretically, it dates back to a century ago when Schwarzschild found a spherically symmetric vacuum solution of Einstein's gravitational field equations, and in the same year, Einstein suggested the existence of gravitational waves as a natural outcome of his general relativity.Observationally, precisely a century later, LIGO with Virgo collaboration [1] detected the gravitational waves, representing the merger of two stellar mass black holes.In 2019, the Event Horizon Telescope (EHT) [2] finally detected the direct image of a black hole and its event horizon.These remarkable discoveries provide some constraints on the modified theories of gravity, such as the existence of a tight bound on the graviton mass [1], the speed of the gravitational wave [3], black hole parameters and properties [4][5][6], and so on, in the light of observational data.
The generalization of the Schwarzschild black hole solution to the electrically charged one [7,8] was done immediately after the work, and the more general and realistic rotating black hole solutions [9,10] were found in the 1960s.However, all of these black hole solutions have a curvature singularity at r = 0 in which spacetime is geodesically incomplete.According to the singularity theorem proved by Hawking and Penrose [11], singularity is inevitable.Moreover, the usual laws of physics break down at the singularity and many physicists believe that quantum gravity effects would work near the singularity.Nevertheless, since we do not have a complete quantum gravity theory yet, another line of work to avoid singularity has been pursued on regular black holes.Such works began with Gliner [12] and Sakharov [13] who proposed a way to avoid the singularity in terms of the matter source, which has a de Sitter core with an equation of state ρ = −p at the center of the spacetime.To avoid the singularity problem, Bardeen [14] also proposed a model of a regular black hole obeying the weak energy condition and thus the model does not obey at least one condition of Hawking-Penrose's singularity theorem.However, the Bardeen's regular black hole solution is not an exact solution to Einstein's equations.After more than three decades, Ayon-Beato and Garcia [15] showed that the Bardeen's metric could be interpreted as a magnetic solution to Einstein's equations coupled to nonlinear electrodynamics.Since then, there has been a lot of work on regular black holes, including Dymnikova [16,17], Bronnikov [18], Hayward [19], Ayon-Beato and Garcia [20], and more [21][22][23][24][25][26][27][28][29][30].These studies have inspired further investigations related to such black holes, for example, regarding particle geodesics [31][32][33][34][35][36][37][38][39], the shadows of regular black holes [40][41][42][43], and the quasi-normal modes [44][45][46][47][48].The thermodynamics and phase transitions for regular black holes have also been studied widely in Refs.[49][50][51][52][53][54][55][56][57][58][59][60][61][62][63].
On the other hand, as is well known, Einstein's general relativity (GR) is a theory of a massless graviton.However, quantum gravity phenomenology [64] at extreme limits has pushed forward to search for alternatives to GR, one of which is to introduce a massive graviton to GR. Historically, it started with Fierz and Pauli [65] who developed a massive theory by extending GR with a quadratic mass term.However, the theory suffers from the Boulware-Deser ghost problem [66] and the van Dam, Veltman, and Zakharov (vDVZ) discontinuity [67,68] in the massless graviton limit.The vDVZ discontinuity was cured by the Vainshtein mechanism [69].After half a century, the notorious Boulware-Deser ghost problem was at last solved by de Rham, Gabadadze, and Trolley (dRGT) [70,71] to have a ghost free massive gravity, which has nonlinearly interacting mass terms constructed from the metric coupled with a symmetric reference metric tensor.These new terms with properly tuned coefficients make it avoid the ghosts order by order.To all orders, the complete absence of the Boulware-Deser ghost was subsequently proven by Hassan and Rosen by a Hamiltonian analysis of the untruncated theory [72,73] and other works [74][75][76][77][78][79].Since then, the dRGT massive gravity has led to new astronomical and cosmological applications for modified gravity [80][81][82][83][84][85][86][87][88][89] including the black hole shadow [90,91].In particular, Hendi et al. [92] have obtained a fascinating result of having the allowed regions of the massive parameters by comparing the black hole shadow in the dRGT massive gravity with the EHT data of M87 * .On the other hand, Vegh [93] further elaborated on another nonlinear massive gravity with a special singular reference metric to apply it to gauge/gravity duality.The modification in the reference metric in the dRGT massive gravity keeps the diffeomorphism symmetry for coordinates (t, r) intact but breaks it in angular directions so that gravitons acquire the mass because of broken momentum conservation [94][95][96].As a result, momentum dissipates as the graviton may behave like a lattice, and it can avoid divergent conductivity.Since then, this Vegh's type of massive gravity, called holographic massive gravity, has been extensively exploited to investigate many interesting models in gravity [97][98][99][100][101][102][103][104][105][106][107][108][109] Very recently, we have studied the tidal effects [110] and statistical entropy [111] in the Schwarzschild black hole in holographic massive gravity.
We have also studied the charged BTZ [112] and Schwarzschild black holes [113] in holographic massive gravity in the global embedding Minkowski spacetime (GEMS) scheme.According to the GEMS scheme, any low dimensional Riemannian manifold can be locally isometrically embedded in a higher dimensional flat one [114][115][116].This can make us have a complete analytic extension of manifolds, or we can use it for visualizing or deriving physical properties of the embedded spacetimes, such as a unified description of Hawking [117] and Unruh effects [118].In this line of work, Deser and Levin [119][120][121] firstly showed that the Hawking temperature for a fiducial observer in a curved spacetime can be considered as the Unruh one for a uniformly accelerated observer in a higher-dimensional GEMS embedded flat spacetime.Since then, there has been much work on the GEMS approach to confirm these ideas in various other spacetimes [122][123][124][125][126][127][128][129][130][131][132][133][134] and an interesting extension to embedding gravity [135][136][137][138]. Later, Brynjolfsson and Thorlacius [139] introduced a local temperature measured by a freely falling observer in the GEMS method.We have also studied various interesting curved spacetimes [140,141] to investigate local temperatures and their equivalence to Hawking ones.
The main goal of this paper is to construct and analyze the GEMS embeddings of spacetimes having regular black holes in massless and massive gravity, whose embeddings are neither found nor even tried at all, as far as we know.Moreover, a recent study on the geodesic completion of a regular black hole [142] supports the needs of the GEMS embeddings of regular black holes.In this respect, it would be interesting to embed a regular black hole with massive gravitons into a higher dimensional flat spacetime.In this paper, we will consider the Hayward regular black hole (HRBH) in massless gravity as a representative of regular black holes and extend it to massive gravity to embed it in higher dimensional flat spacetimes.Here, the HRBH in massless gravity means the original Hayward black hole, while the HRBH in massive gravity is the one having massive gravitons obtained from the consideration of Vegh's type of massive gravity.We note that when the Hayward parameter vanishes, the HRBH in massless gravity becomes the Schwarzschild black hole, and when massive gravitons are turned off, the HRBH in massive gravity is reduced to the HRBH in massless gravity.
The remainder of the paper is organized as follows.In Sec.II, we will newly find solutions to the HRBH in holographic massive gravity.We will first briefly summarize the known solution of the HRBH in massless gravity and generalize it to one in holographic massive gravity.Then, we will show that the Kretschmann scalar for the HRBH in massive gravity is not regular near r = 0.In Sec.III, we will construct the GEMS embeddings of the HRBH both in massless and in massive gravity.As a result, making use of embedding coordinates, we will obtain Unruh, Hawking, and freely falling temperatures seen by different observers.Conclusions are drawn in Sec.IV.Lastly, since embedding coordinates of these regular black holes are very complicated, we have separately listed them and then shown their limits explicitly from the massive to the massless cases in Appendix A.

II. HRBH IN MASSIVE GRAVITY
The (3+1)-dimensional HRBH in holographic massive gravity is described by the action where R is the scalar curvature of the metric g µν , m is the black hole mass, l is a length-scale Hayward parameter present in the Hayward solution, m is a graviton mass1 , c a are the coupling constants, f µν is a fixed symmetric tensor usually called the reference metric, and U a are symmetric polynomial potentials of the eigenvalue of the matrix K µ ν ≡ √ g µα f αν given as Here, the square root in K means ( √ A) µ α ( √ A) α ν = A µ ν and square brackets denote the trace [K] = K µ µ .Indices are raised and lowered with the dynamical metric g µν , while the reference metric f µν is a non-dynamical, fixed symmetric tensor that is introduced to construct nontrivial interaction terms in holographic massive gravity.
Variation of the action (2.1) with respect to the metric g µν leads to the equations of motion given by with U −a = 0 and U 0 = 1.When one considers the spherically symmetric black hole solution ansatz as with the following degenerate reference metric one can find Note that the choice of the reference metric in Eq. (2.5) preserves general covariance in (t, r) but not in the angular directions.This gives the symmetric potentials as Therefore, we see that there are no contributions from c 3 and c 4 terms which would appear in (4+1) and (5+1)dimensional spacetimes, respectively.Then, we finally obtain the new solution of the HRBH in massive gravity as with graviton mass parameters R = c 0 c 1 m2 /4 and C = c 2 0 c 2 m2 .Note that m is an integration constant related to the mass of the black hole and c 0 is a positive constant.In the limit of R → 0 and C → 0, it becomes the HRBH in massless gravity, and when l → 0, it is further reduced to the Schwarzschild metric as expected.On the other hand, when R = 0 and C = 0 with l = 0, it becomes the Schwarzschild black hole in massive gravity [113].
It is appropriate to comment on the two terms due to massive gravitons.Firstly, the C term in Eq. (2.8) reminds us of a monopole solution introduced by Barriola and Vilenkin [145], which comes from a topological defect in the early Universe as a result of a gauge symmetry breaking.On the other hand, the R term is not uncommon in gravity theories, which also arises in, such as the dRGT massive gravity [70,71], conformal gravity [146], and f (R) gravity [147].Physically, the linear term in Eq. (2.8) stands for a deviation between the solution and the HRBH spacetime in massless gravity, as in the modified Newtonian dynamics (MOND) studies a deviation of massive bodies in the solar system from the Newtonian mechanics [148].The linear term also affects the radius of the photon sphere in a black hole, the size of the shadow, which is smaller than the one in a Schwarzschild spacetime [149].In addition, Hendi et al. [92] have shown that the black hole shadow in the dRGT massive gravity consistent with EHT data is given for small R < 0.072 and negative −0.3 < C < −0.03 by comparing the black hole shadow with the data of EHT collaboration.We note that the additional C and R terms in Eq. (2.8) are obtained from the consideration of massive gravitons, while the other theories have different causes.
It is well known that the HRBH solution in massless gravity (R = 0 and C = 0) has two horizons for m > m * ≡ 3 √ 3l/4, one for m = m * , and none for m < m * [19].However, for the case of the HRBH in massive gravity, these are modified by R and C. While we will find their exact solutions of f (r H ) = 0 in the following subsections, we summarize here the numbers of event horizons of the HRBH in holographic massive gravity according to various values of R and C in Tables I and II.In the Tables, max(2) or max(3) mean that there exist maximum two or three horizons, respectively, satisfying f (r H ) = 0 in the given range of R and C. Otherwise, a given number of event horizons are allowed in the range.In order to describe these clearly, as an example, in Fig. 1, we have plotted f (r)-graphs for −1 < C < 0 with R for the case of m > m * where one can see the changes in the number of event horizons according to R and C, respectively.TABLE I: Numbers of event horizons of the HRBH in holographic massive gravity according to R and C when m > m * .Note that in the table, max(2) or max(3) denote the cases that maximum two or three horizons are formed satisfying f (rH ) = 0 in the given range, respectively.
One can see from the Tables I, II, and Fig. 1 that when m > m * , it would be physically interesting in both all ranges of C with R > 0 and all ranges of R with C > −1.In those ranges, there exists not only an outer, but also at least one inner event horizon.When m = m * , the extremal case for the HRBH in massless gravity, it remains extremal for C ≤ −1 and R > 0 in the HRBH in massive gravity.However, for all ranges of R with −1 < C < 0 TABLE II: Numbers of event horizons of the HRBH in holographic massive gravity according to R and C when m = m * .and R < 0 with C ≥ 0, it changes to have two and more event horizons.As a result, according to the graviton mass parameters R and C, it is expected that the thermodynamics of the HRBH in massive gravity is differently described from the HRBH in massless gravity.
On the other hand, from the metric solution (2.8), one can find that the event horizon r H determines the mass as . (2.9) In Fig. 2, we have plotted the mass function by comparing the HRBH in massless gravity with ones in holographic massive gravity.As a result, one can see that for a fixed R (or C), as C (or R) decreases, the outer event horizon r H increases in the HRBH in massive gravity.Finally, it is appropriate to comment on the Kretschmann scalar, which is known to be finite for regular black holes, for the HRBH in massive gravity.First of all, for the HRBH in massless gravity [27], the Kretschmann scalar is given by At the center of the spacetime as r → 0, due to the Hayward parameter l, it is finite as which shows that the HRBH in massless gravity is regular everywhere with no curvature singularity.On the other hand, as l → 0, it becomes the Kretschmann scalar of the Schwarzschild spacetime in massless gravity, as expected.Now, when massive gravitons are introduced as the HRBH in massive gravity having the lapse function (2.8), we have new Kretschmann scalar modified by R and C as (2.13) In the massless limit of R = 0 and C = 0, it becomes the Kretschmann scalar (2.10) of the HRBH in massless gravity.However, near the center as r → 0, the last term in Eq. (2.13) becomes dominant so it diverges as Thus, the regular structure of the HRBH in massless gravity has been changed to singular due to massive gravitons, less-severely divergent than the Schwarzschild case though.Note that in the massive Schwarzschild limit as l → 0, the Kretschmann scalar is reduced to As a result, we have found that introducing massive gravitons to the HRBH in massless gravity affects the spacetime structure from regular to singular.This singular structure of spacetime can be relieved by enlarging the dimensions of the spacetime according to the well established GEMS scheme, which we will consider in Section III.

A. Solutions of the HRBH in massless gravity
Let us first briefly summarize the solutions of the HRBH in massless gravity satisfying f (r H ) = 0 in Eq. (2.8) where rH denotes the horizon radius of the black hole in massless gravity.In the massless case with the graviton mass parameters R = C = 0, event horizons can be found from the following cubic equation rewritten as By following the general procedure from Eq. (2.32) to Eq. (2.37) presented in the next subsection, one finally gets the solutions of the HRBH in massless gravity as whose behaviors are depicted in Fig. 3. Here, cos ψ is defined as Note that the solutions (2.17  FIG.3: Solutions for the HRBH in massless gravity: the blue curve is for rH2, the black curve for rH3 and the red curve for rH1 with l = 1.Here, we see that rH2 corresponds to an outer horizon (r+) and rH3 to an inner horizon (r−).rH1 is negative so discarded.

B. Solution of the HRBH in massive gravity
In this subsection, we will newly find general solutions of the HRBH in holographic massive gravity (2.8) satisfying f (r H ) = 0 where r H denotes the horizon radius of the black hole in massive gravity.First of all, one can rewrite f (r H ) = 0 as which is a quartic equation with R = 0, compared with the cubic one in the previous massless case.After dividing by 2R and by introducing a new variable y as one can find the standard quartic equation written as where Note that at this stage one cannot simply take the limit of R → 0 since the coefficients are all divergent.Thus, if one wants to find the limit, it should consider the equation (2.20) again from the start.The quartic equation (2.22) can be solved by adding and subtracting xy 2 + x 2 /4 to Eq. (2.22) as which can be rewritten as Thus, if we demand the last two terms to vanish as one can have These are the product of two quadratic equations whose roots can be easily obtained separately.As a result, we have the following four roots for the quartic equation (2.22) with In Eq. (2.29), note that x 0 is a root of the cubic equation (2.26) rewritten as Then, we have where b 1 , b 2 , and b 3 are the coefficients of the standard cubic equation of These give us to yield Now, we define cos ψ [143,144] as (2.35) Making use of q 1 and q 2 in Eq. (2.34), we find Then, the solutions for the cubic equation (2.30) are given by (2.37) Thus, making use of these solutions, we finally have the following four roots of the quartic equation (2.22) as where Here, we have chosen x 1 as x 0 in p 1 .In Fig. 4, we depict a set of solutions for the HRBH in massive gravity for R = 0.01 and C = −0.1.As explained in Table I, for the chosen R and C, we expect that there are two event horizons, and Fig. 4 shows the same behavior that there are two, i.e., one is an outer r H1 and the other is an inner horizon r H2 , respectively.The remaining two r H3 , r H4 are of no physical meanings for event horizons since they are negative.In Fig. 5, we also draw the solution for the HRBH in massive gravity compared with the HRBH in massless gravity, where one can see how massive gravitons change r H (m).
Note that r Hi (i = 1, 2, 3, 4) in Eq. (2.38) denote the event horizons of the HRBH in holographic massive gravity, while rHi (i = 1, 2, 3) in Eqs.(2.17) stand for those of the HRBH in massless gravity.It is also appropriate to comment that it is not possible directly to get the solutions of the HRBH in massless gravity by taking R → 0 and C → 0 from the solutions (2.38) of the quartic equation.Note that this can be understood due to the fact that they are obtained from the implicit condition of R = 0 as in Eqs.(2.20) and (2.21).

III. GEMS EMBEDDING OF HRBH A. HRBH in massless gravity
The (3+1)-dimensional HRBH in massless gravity can be embedded in a (5+2)-dimensional Minkowski spacetime given by where embedding coordinates are obtained as with In the above embedding functions (3.2), kH is the surface gravity defined as where ξ µ is a Killing vector and rH is the event horizon of the HRBH in massless gravity.Note that in the limit of l → 0, z 6 vanishes, and z 5 becomes Therefore, with z 0 • • • z 4 expressed in the same limit, the embedding coordinates (3.2) are correctly reduced to the well-known (5+1)-dimensional GEMS embeddings of the Schwarzschild black hole [114].
Here, we note that from the original spacetime metric, the Hawking temperature TH seen by an asymptotic observer and a local fiducial temperature measured by an observer who rests at a distance from the black hole are simply found as respectively.Now, let us consider the Unruh effect in the embedded flat spacetime, which originally states that accelerated observers or detectors with an acceleration a along the x direction by following the trajectory a −2 = x 2 − t 2 measure the Unruh temperature given by 2πT = a [118].To apply this for in a higher dimensional flat spacetime, we notice that the static detectors (r, θ, φ = constant) in the original curved spacetime are described by a fixed point in the (z 2 , z 3 , z 4 , z 5 , z 6 ) plane on the GEMS embedded spacetime.Then, an observer who is uniformly accelerated in the (5+2)-dimensional flat spacetime, follows a hyperbolic trajectory described by Thus, one can find the Unruh temperature for the uniformly accelerated observer in the (5+2)-dimensional flat spacetime as (3.9) This corresponds to the fiducial temperature (3.7) for the observer located at a distance from the HRBH in massless gravity.The Hawking temperature TH seen by an asymptotic observer can be obtained as As a result, one can see that the Hawking effect for a fiducial observer in a black hole spacetime is equal to the Unruh effect for a uniformly accelerated observer in a higher-dimensional flat spacetime.Now, let us find a freely falling acceleration and corresponding temperature in the (5+2)-dimensional embedded flat spacetime.For an observer who is freely falling from rest r = r 0 at τ = 0, the equations of motion are where (−) sign is for inward motion.Then, making use of the embedding coordinates in Eq. (3.2) and the geodesic equations in Eq. (3.11), one can explicitly find a freely falling acceleration ā7 in the GEMS embedded (5+2)dimensional spacetime as where with Note here that r 0 is replaced with r in Eq. (3.12).
According to the Unruh's prescription, the freely falling acceleration gives us the freely falling temperature measured by the freely falling observer as Making use of the dimensionless parameters x = rH /r and b = l/r H , the squared freely falling temperature can be written as where As r → ∞, the freely falling temperature TFF is reduced to the Hawking temperature (3.6).Moreover, as l → 0, it becomes the freely falling temperature of the Schwarzschild black hole as where rH = 2m is the radius of the event horizon of the Schwarzschild black hole.Note also that as r → ∞, T Sch FF becomes the Hawking temperature of the Schwarzschild black hole, while as r → rH , the freely falling temperature becomes not diverge, but remains finite at the event horizon.In Fig. 6, we have depicted the ratio of the squared freely falling temperatures to the squared Hawking temperature, T 2 FF / T 2 H .For comparison purposes, it is shown in Fig. 6(a) that the freely falling temperature of the Schwarzschild black hole, which is a prototype of a singular black hole, is finite at the event horizon of x = 1 (r = rH ), while it becomes the Hawking temperature at asymptotic infinity of x → 0 (r → ∞).On the other hand, as in Fig. 6(b), for the HRBH in massless gravity, the freely falling temperatures are finite only when 0 < b < b c (≡ 0.3568) at the event horizons, and when b = b c , it vanishes at the event horizon.Meanwhile, the freely falling temperatures become the Hawking temperature at asymptotic infinity.The parameter b also has an upper bound coming from the Hawking temperature (3.6), which is defined when rH ≥ √ 3l (or b ≤ 1/ √ 3 = 0.5774) where T 2 FF / T 2 H diverges.In between b c < b < 0.5774, one can see in Fig. 6(b) that the freely falling temperatures of the HRBH behave quite differently from the Schwarzschild singular black hole.They are going up and down, and then become negative.It is well known that the negativity of squared freely falling temperatures is not entirely prohibited, which means that there is no thermal radiation.This is allowed for a geodesic observer who follows a spacelike motion similar to the case of the Schwarzschild-AdS black hole in massless gravity [113,119,139].

B. HRBH in massive gravity
Now, after a lengthy calculation, we newly find that the (3+1)-dimensional HRBH in massive gravity can be embedded in a (6+3)-dimensional Minkowski spacetime as whose embedding coordinates are explicitly written as Here, the surface gravity of the HRBH in massive gravity is given by Also, H 4 , H 5 and H 6 are defined as ) ) are given in Appendix A. Note that the above GEMS embedding is explicitly carried out under the assumption of R ≥ 0 and C ≥ −1, which is physically much more interesting than others as shown in Figs. 1 and 2, and Tables I and II.It seems appropriate to comment that in the massless gravity limit of R → 0 and C → 0, H 4 and H 5 become unity, and H 6 becomes H 0 in Eq. (3.3).Also, when subtracting the embedding coordinates z 6 from z 5 , the coefficients of u i − ūi (i = 1, 2, • • •, 7) in the numerator of z 5 − z 6 are reduced to Eq. (A.9) in the Appendix.In the same way, by subtracting the embedding coordinates z 8 from z 7 , the coefficients of v i − vi (i = 1, 2, • • •, 6) in the numerator of z 7 − z 8 becomes Eq.(A.10).Since H 6 is reduced to H 0 in the massless limit, one can find that they are exactly the same coefficients in z 5 and z 6 , respectively, of the HRBH in massless gravity in Eqs.(3.2).As a result, one can see that when R → 0 and C → 0, the (6+3)-dimensional embedding coordinates (3.21) of the HRBH in massive gravity have proper limits of the (5+2)-dimensional ones (3.2) in massless gravity.Moreover, when one additionally takes the limit of l → 0, z 7 and z 8 identically vanish, and the combination of z 5 and z 6 by subtraction is finally reduced to where rH = 2m is the event horizon of the Schwarzschild black hole.Thus, one can see that the embedding coordinate of z 5 − z 6 is finally reduced to z 5 , which is one of the spacelike embedding coordinates of the Schwarzschild black hole in the (5+1)-dimensional GEMS scheme.
Here, we note again that from the original spacetime metric, the Hawking temperature T H seen by an asymptotic observer can be found as and a local fiducial temperature measured by an observer who rests at a distance from the black hole is given by In the massless limit of R → 0 and C → 0, one can easily find that it becomes the fiducial temperature (3.7) since H 5 → 1 and H 6 → H 0 .On the other hand, in order to investigate the Unruh effect in the GEMS embedded flat spacetime, we note that the static detectors (r, θ, φ = constant) in the original curved spacetime are described by a fixed point in the (z 2 , z 3 , z 4 , z 5 , z 6 , z 7 , z 8 ) plane on the GEMS embedded spacetime.Then, an observer who is uniformly accelerated in the (6+3)-dimensional flat spacetime, follows a hyperbolic trajectory in (z 0 , z 1 ) described by Thus, as before, one can arrive at the Unruh temperature for the uniformly accelerated observer in the (6+3)dimensional flat spacetime as This corresponds to the fiducial temperature for the observer located at a distance from the HRBH in massive gravity.The Hawking temperature T H seen by an asymptotic observer can be obtained as As a result, one can see that the Hawking effect for a fiducial observer in a black hole spacetime is equal to the Unruh effect for a uniformly accelerated observer in a higher-dimensional flat spacetime.Now, let us find a freely falling acceleration and corresponding temperature in the (6+3)-dimensional embedded flat spacetime.For an observer who is freely falling from rest r = r 0 at τ = 0, the equations of motion are Then, making use of the embedding coordinates in Eq. (3.21) and the geodesic equations in Eq. (3.30), one can explicitly find a freely falling acceleration a 9 in the GEMS embedded (6+3)-dimensional spacetime as Here, the dashed black curve is for the HRBH in massless gravity as in Fig. 6 with b = 0.3, and the dotted black curve is for the Schwarzschild black hole in massless gravity.Also, the vertical dotted lines are drawn at event horizons. where One can easily check that in the massless limit of R → 0 and C → 0, D 3 and D 4 are reduced to D 1 and D 2 , respectively in Eq. (3.13) and thus the freely falling acceleration ā2 9 becomes ā2 7 in Eq. (3.12) of the HRBH in massless gravity.According to the Unruh's prescription, one can find T FF measured by the freely falling observer from the freely falling acceleration as (3.34) In the massless limit of R → 0 and C → 0, this is exactly the same with the previous one, the freely falling temperature of the HRBH in massless gravity in Eq. (3.15).Moreover, as r → ∞, the freely falling temperature T FF is reduced to the Hawking temperature.Now, making use of the dimensionless parameters x = r H /r, b = l/r H and d = Rr H , the squared freely falling temperature can be written as where with In Fig. 7, we have depicted the ratio of the squared freely falling temperatures to the squared Hawking temperature, T 2 FF /T 2 H for the HRBH in massive gravity.One can see that at the event horizon the freely falling temperatures are all finite, while the fiducial temperature diverges [139].In the limit of b → 0 (or l → 0), which corresponds to the case of the Schwarzschild black hole in massive gravity, the freely falling temperature is reduced to This is exactly the same with T 2 FF in Ref. [113].Furthermore, in the massless limit of d → 0 (or R → 0) and C → 0, it becomes the freely falling temperature of the Schwarzschild black hole in massless gravity [139].

IV. DISCUSSION
In this paper, we have newly studied the Hayward regular black hole (HRBH) in massive gravity, which is a modification of the HRBH in massless gravity to have nonzero mass of gravitons as proposed by Vegh in the framework of holography.By solving Einstein's equations, we have found a solution of the HRBH in massive gravity and analyzed the novel structures of event horizons classified by the graviton mass parameters R and C, qualitatively.Concretely, when m > m * , physically interesting ranges of R and C lie in both all C with R > 0 and all R with C > −1.In those ranges, there exist an outer event horizon and at least one inner event horizon.When m = m * , which is the extremal case for the HRBH in massless gravity, it remains extremal for C ≤ −1 and R > 0. However, for all R with −1 < C < 0 and R < 0 with C ≥ 0 in the HRBH in massive gravity, it changes to have two and more event horizons.Therefore, according to the specific values in the graviton mass parameters R and C, we expect that the thermodynamics of the HRBH in massive gravity would be different from the one of the HRBH in massless gravity.Furthermore, to find full information on event horizons, we have also analytically solved the metric function f (r H ) = 0 of the HRBH in massive gravity, which is a non-trivial quartic equation due to massive gravitons, while the HRBH in massless gravity is the cubic equation.
After exploiting the geometric property of the spacetime, we have proceeded to study the GEMS embeddings of the HRBHs in massless and massive gravity where the former is geodesically incomplete and the latter is singular due to massive gravitons.As a result, we have globally embedded the (3+1)-dimensional HRBH in massless gravity into a (5+2)-dimensional flat Minkowski spacetime, and the (3+1)-dimensional HRBH in massive gravity into a much higher (6+3)-dimensional flat Minkowski spacetime, respectively, where the difference in embedding dimensions comes from whether or not holographically introduced massive gravitons exist.Furthermore, making use of newly obtained embedding coordinates, we have directly obtained Unruh temperatures and compared them with the Hawking and local fiducial temperatures, showing that the Unruh effect for a uniformly accelerated observer in a higher-dimensionally embedded flat spacetime is equal to the Hawking effect for a fiducial observer in the corresponding curved spacetime.We have also obtained freely falling temperatures of the HRBH in massless and massive gravities seen by freely falling observers following their geodesic trajectories, which remain finite even at the event horizons.These are different from the Hawking temperatures divergent at the event horizons.
Finally, it seems appropriate to comment that the solutions of the HRBH in massless gravity is not possibly obtained simply by taking R → 0 and C → 0 due to their implicit condition of R = 0 as seen in Eqs.(2.20) and (2.21).On the other hand, the GEMS embedding of the HRBH in massive gravity can be reduced to the corresponding GEMS embedding of the HRBH in massless gravity in the limit of R → 0 and C → 0 by redefining some embedding coordinates.As a result, when R → 0 and C → 0, the (6+3)-dimensional embedding coordinates (3.21) of the HRBH in massive gravity have proper limits of the (5+2)-dimensional ones (3.2) in massless gravity, which may be a characteristic of the Hayward nonsingular black hole in massive gravity.

FIG. 4 :
FIG.4: Solutions for the HRBH in massive gravity: in (a) the blue curve is for rH1, the black curve for rH2, which correspond to an outer horizon (r+) and to an inner horizon (r−), respectively.In (b), the curves are for rH3 (upper) and rH4 (lower), which are negative so discarded.Here, we set R = 0.01, C = −0.1 with l = 1.

FIG. 5 :
FIG.5: Solutions for the HRBH in massless and massive gravity: the blue and black solid curves are for massive gravity, while the dotted curves for massless gravity.Here, we set R = 0.01, C = −0.1 with l = 1.

FIG. 6 : 2 FF / T 2 H
FIG. 6: Squared freely falling temperatures T 2 FF / T 2 H for the HRBH in massless gravity drawn by a dimensionless parameter x (= rH/r).(a) the freely falling temperature of the Schwarzschild black hole in massless gravity which corresponds to b = 0 (or, l = 0).(b) the freely falling temperatures of the HRBH in massless gravity for b = 0.1, 0.2, 0.3, bc, 0.4 from top to bottom.Here, the vertical dotted lines are drawn at event horizons.

FIG. 7 : 2 FF /T 2 H
FIG. 7: Squared freely falling temperatures T 2 FF /T 2 H for the HRBH in massive gravity drawn by a dimensionless parameter x (= rH/r).(a) the freely falling temperature of the HRBH in massive gravity for d = 0.01, 0.03, 0.05, 0.07 from top to bottom with a fixed C = −0.1 and b = 0.3.(b) the freely falling temperature of the HRBH in massive gravity for C = −0.4,−0.3, −0.2, −0.1 from top to bottom with a fixed d = 0.07 and b = 0.3.Here, the dashed black curve is for the HRBH in massless gravity as in Fig.6with b = 0.3, and the dotted black curve is for the Schwarzschild black hole in massless gravity.Also, the vertical dotted lines are drawn at event horizons.