Upper bound on the GUP parameter using the black hole shadow

An upper bound on the parameter that provides a generalized uncertainty principle (GUP) is obtained from the black hole shadow. With the aid of a recent constraint between regular black holes and the GUP parameter, it is indicated a relation between this parameter and the deviation from circularity of the black hole shadow. In the case of the recent announcement of the M87* results from the Event Horizon Telescope collaboration, a deviation from circularity ≲10%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lesssim 10\%$$\end{document} imposes a GUP parameter β0<1090\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta _0 <10^{90}$$\end{document}.


Introduction
The Event Horizon Telescope collaboration announced the first black hole image ever captured [1,2]. By using interferometry, the collaboration built the shadow of the supermassive black hole (M87*) at the center of the Messier 87 galaxy. As pointed out by the collaboration, the black hole shadow is well-described by general relativity by adopting the Kerr metric in order to interpret the phenomenon. Moreover, the supermassive M87* -whose mass is M = (6.5 ± 0.7) × 10 9 M and is r obs = (16.8 ± 0.8)Mpc distant from us -presents an almost circular shadow. Accordingly, the shadow of M87* indicates a deviation 10% from circularity. Even with the good agreement with the Kerr metric, options are not totally ruled out. Due to the uncertainty on both the rotation parameter and observation angle, options to Kerr metric are still candidates for M87* shadow [3]. This article explores that avenue and applies the M87* parameters to Kerr-like objects, like rotating regular black holes.
Shadows of black holes have been studied since the pioneer work of Synge [4], in which the shadow of the Schwarzschild black hole was obtained. Bardeen [5] built the first shadow for a rotating black hole, namely, for the Kerr metric. Since then, whether in the general relativity context or beyond, shadows of different black holes have been puba e-mail: nevesjcs@if.usp.br (corresponding author) lished, like shadows for the Reissner-Nordström black hole [6], the Kerr-Newman black hole [7], black holes with a cosmological constant [8], Kerr-Newman-NUT black hole [9], Kerr-Newman-Kasuya black hole [10], braneworld black holes [11,12], rotating wormholes [13] and for regular black holes [14][15][16], which are focused on this article. It is worth emphasizing that the results exposed here do not consider an accretion disk around the black hole. Accretion disks modify the shadow's silhouette, and models that consider such influences are being developed [17][18][19]. The results presented in this article can be used in future researches that include accretion disks, whether made of plasma or dark matter.
As an alternative to singular black holes, regular black holes are solutions of the gravitational field equations without a singularity inside the event horizon. Once again, Bardeen was pionner when constructed the very first regular metric that describes a regular black hole [20]. The Bardeen black hole is spherical, i.e., a nonrotating black hole, and is described by a mass function that depends on the radial coordinate. Inside the event horizon, the Bardeen metric hides a de Sitter core instead of a Schwarzschild singularity. A de Sitter core inside the event horizon avoids either the point singularity of spherical black holes [21][22][23][24][25][26][27][28][29] or the ring singularity in geometries with axial symmetry [30][31][32][33][34][35]. As is well-known, the de Sitter core provides energy conditions violations in order to avoid consequences of the singularities theorems. 1 Then, the Bardeen regular black hole violates the strong energy condition, and, according to our work [35], those with rotation ignore the weak energy condition.
Besides a mass function that depends on the radial coordinate, regular black holes possess a mass function that depends on the parameter related directly to the regular geometry or the absence of a singularity. In Ayón-Beato and Garcia's work [37], such a parameter is a charge, and the Bardeen regular black hole is conceived of as a solution of general relativity coupled to a nonlinear electrodynamics for those authors. However, recently, our work [38] suggested another interpretation to the Bardeen metric. By using a generalized uncertainty principle (GUP), we computed quantum corrections to the Hawking temperature for the Schwarzschild black hole. Thus, we pointed out that, at second order, Bardeen's regular black hole may be view as a quantum-corrected Schwarzschild black hole. This new interpretation led to a constraint between a metric parameter, namely, the parameter in the mass function that produces regularity, and the GUP parameter. With the aid of such a relation, the GUP parameter will be estimated using the black hole shadow. The regular black holes used in this shadow calculation come from a class of rotating regular black holes with a cosmological constant published in our article [35]. GUPs appear, for example, in candidates for quantum gravity theories and in the deformed quantum mechanics, where the canonical commutation relations are modified in order to provide a minimal length scale and, consequently, its contribution to empirical results, like in the modified hydrogen-atom spectrum [39], Lamb shift, Landau Levels and in the scanning tunneling microscope [40,41]. 2 The dimensionless GUP parameter (also called quantum gravity parameter) deforms the Heisenberg relation, and with β 0 = 0, where β 0 is the mentioned parameter, Heisenberg's uncertainty relation is recovered. There is a debate on the value of the GUP parameter [40,43,44]. Assuming that the GUP parameter is β 0 ∼ 1 implies that effects of β 0 are hard or too small to be detected. But without such an a priori assumption, it is possible to obtain upper bounds on the GUP parameter by using recent experiments. Like Ref. [43], where the upper bounds were built with the aid of the light deflection and perihelion precession, and Ref. [44], where gravitational waves were adopted, in this work an option in the strong gravitational field regime (the shadow of M87*) is used in order to obtain an upper bound on β 0 . In general, as we will see, gravitational systems offer worse upper bounds than quantum options, like the Lamb shift, Landau Levels or the scanning tunneling microscope used in Ref. [40].
As I said, the deviation from circularity, reported by the collaboration, will be used in order to estimate the GUP parameter. It is worth emphasizing that the class of rotating regular black holes used here and presented in Ref. [35] generalizes earlier Kerr-like regular solutions because it possesses a general mass function m(r ) and a cosmological constant. The shadow of that class has its silhouette presented here for the first time in the literature. With the M87* parameters, by assuming that the angle between the black hole rotation axis and the observer is θ obs = 17 • (in agreement with observed jets supposedly aligned with the rotation axis [1,45]), it is indicated the β 0 < 10 90 as upper bound on the GUP parameter. 2 See Ref. [42] for a review on GUPs.
The structure of this paper is as follows: in Sect. 2, the geodesic equations for the class of rotating regular black holes and equations that provide the shadow's silhouette were derived. In Sect. 3, two shadow's observables are indicated, oblateness and root-mean-square distance from the average radius of the shadow, which gives the deviation from circularity, such that the latter was computed in order to provide an upper bound on the GUP parameter in Sect. 4. In Sect. 5, the final remarks are made. 3

The shadow of rotating regular black holes
Let us obtain the shadow of rotating regular black holes with a cosmological constant in this section. Firstly, the spacetime metric and its geodesic equations are shown, then the silhouette equations are computed.

Spacetime metric and geodesic equations
In this article, shadows are obtained from a class of rotating regular black holes developed in Ref. [35]. The spacetime metric of that class -using the Boyer-Lindquist coordinates -is given by where Δ θ = 1 + Λ 3 a 2 cos 2 θ, Σ = r 2 + a 2 cos 2 θ, and The constant Λ is the cosmological constant, a is the rotation parameter and, for the class of regular black holes studied here, the black hole mass depends on the radial coordinate r , namely, The mass function (4) provides black holes without a singularity. 4 For different values of the integer q, we have wellknown regular black holes. For example, for the spherical case, q = 2 provides the Bardeen black hole [20], and q = 3 produces the Hayward regular metric [27], but the mass function also gives axisymmetric regular black holes or Kerr-like black holes. The parameter M is the mass parameter (for large values of r , m(r ) ∼ M), and r 0 is a length parameter that provides regular metrics, parameter conceived of as a microscopical constant related to both the GUP parameter and the Planck length according to our work Ref. [38] (in Sect. 4 such a parameter will be briefly discussed). Due to cosmological observations [49], I will focus on the positive cosmological constant in this article. In this case, regarding r 0 M, the function Δ r provides three roots: the inner horizon r − , the event horizon r + , and the cosmological horizon r ++ . Then, the spacetime structure reads In order to construct the black hole shadows of the class of rotating regular black holes given by Eq. (1), the geodesic equations are needed. Geodesics for the Kerr metric were obtained by Carter [50], who showed the separability of the geodesic equations. Carter argued that a test particle in the Kerr spacetime possesses four constants of motion along geodesics. Accordingly, one has the two Killing vector fields ξ t and ξ φ with their respective constants, the mass of the test particle and the fourth constant, which is called Carter constant (indicated by K ). Both the Kerr-anti-de Sitter and Kerr-de Sitter spacetimes present those same four constants for test particles along geodesics, whether in the general relativity realm [51] or in the brane world, as we indicated in Ref. [52]. The class of rotating regular black holes given by (1) also presents these constants.
The geodesic equations for the spacetime (1) are obtained from the Hamilton-Jacobi equation: where σ is a parameter related to the affine parameter τ (i.e., τ = δσ , with δ playing the role of mass of the test particle along the geodesic). The function S is the Jacobi action, which is related to the generalized momentum through From the two Killing vector fields ξ t and ξ φ given by the geometry with axial symmetry (1), we have the particle's constants of motion E and L, namely, energy and angular momentum, respectively, Following Carter, it is assumed that S may be written as where plus and minus mean the de Sitter and anti-de Sitter cases, respectively, 5 and expressions for S θ and S r will be omitted here. As is known in general relativity, from the Lagrangian L , the generalized momentum is defined as p μ ≡ ∂L ∂ẋ μ = g μνẋ ν , where dot means ordinary derivative with respect to the parameter σ . Therefore, from that definition and Eq. (9), substituted into Eq. (7), one has the geodesic equations for the metric (1), which in the coordinate basis are where Σ is given by Eq.
(2). The functions P, R and Θ are written as As I said, the parameter δ represents the mass of the particle along the geodesic, so δ = 0 in the case studied here. The constant Q is related to Carter's constant K , that is to say, which vanishes for equatorial orbits.

The shadow's silhouette
In particular, for the shadow phenomenon, only photons orbits, or null geodesics, will be adopted (δ = 0). And as the cosmological context is considered in this article, the cosmological constant will be assumed positive, i.e., Λ > 0. Following [9], two new parameters are defined as which are constants in the shadow's silhouette. Such a silhouette is given by unstable photon orbits with r = r p constant outside the event horizon, i.e., r p > r + . In such orbits, photons may either fall into the black hole or escape to the observer position. Thus, according to geodesic equationṙ , we have R(r p ) = R (r p ) = 0 in order to provide the unstable orbits (the symbol means derivative with respect to r ). Typically, for a black hole with rotation r p = r p− and r p = r p+ (with r p− ≤ r p+ ), i.e., there are both a minimum and a maximum value for r p , and the edge of the black hole shadow should be built for those values of r p . That is, as we will see, the left and the right sides of the shadow can be different for rotating black holes due to the spacetime dragging. On the other hand, for the Schwarzschild black hole r p− = r p+ = 3M (in the Schwarzschild case, M is the black hole mass or the Arnowitt-Deser-Misner mass), and the shadow is perfectly symmetrical. The equations that involve R(r p ) and its derivative lead to and As we will see, these conserved quantities are part of the equations that "draw" the shadow. It is worth emphasizing that the metric (1) is not necessarily asymptotically flat, it is either asymptotically de Sitter or anti-de Sitter for Λ = 0. Thus, a given observer is not at infinity describing the black hole shadow, indeed he/she is in the domain of outer communication, which is the region defined in between the event horizon r + and the cosmological horizon r ++ in the de Sitter case. Therefore, the observer will be located at the point with coordinates (r obs , θ o ), according to Fig. 1. In terms of the observer position, because of the axial symmetry, the black hole shadow depends only on the radial and polar coordinates. The coordinate θ o stands for the observer angle, which is the position in which the shadow is observed in relation to the rotation axis, and the parameter r obs is the distance from the black hole. Following [9], it is adopted the orthonormal tetrad e in order to describe the shadow silhouette. That is, the null congruences coming from r p− ≤ r ≤ r p+ reaches the observer and are projected onto the tetrad, then the shadow phenomenon is described by using e μ a . Such a tetrad is written with the aid of the coordinate basis vectors (∂ t , ∂ r , ∂ θ , ∂ φ ), According to Grenzebach et al. [9], the direction of e 3 points toward the black hole (see Fig. 1), and e 0 is the observer's four-velocity, then in this case the observer is not necessarily at rest. The tetrad is chosen such that e 0 ± e 3 is tangential to the principal null congruence direction for a metric like (1). The null congruence, that is to say, the light rays that come from the region defined by r p , are conceived of as curves γ such that their tangent vectors are given bẏ in the coordinate basis. According to Ref. [9], at the observer position,γ is described by the tetrad e μ a , i.e., γ = ζ (−e 0 + sin α cos βe 1 + sin α sin βe 2 + cos αe 3 ) , (23) where the new angles α and β are celestial coordinates as indicated in Fig. 1. The description of the black hole shadow is made by using these coordinates. The factor ζ is obtained from Eqs. (22)-(23) and the tetrad equations, such that From the factor ζ , by comparing terms in Eqs. (22) and (23) with the aid of the equations for the tetrad, the celestial coordinates are straightforwardly written as Then, using the geodesic equations forṙ andφ, substituting them into Eqs. (25)-(26), we obtain simple expressions for α and β, namely, As we can see, the equations that describe the celestial coordinates depend on the observer position (r obs , θ o ) and parameters of the photon orbits, also known as photon "sphere". In order to describe the shadow's silhouette it is appropriate to define the Cartesian coordinates y(r p ) = −2 tan α(r p ) 2 cos β(r p ) .
Therefore, the parametric equations (29) and (30) draw the shadow's silhouette. According to Fig. 2, the shadow is symmetrical to the x-axis. For a Kerr black hole with mass M, the silhouette depends on the rotation parameter a and the observer polar angle θ o , considering an observer at infinity. On the other hand, for the de Sitter case, the observer is in the domain of outer communication (r = r obs ), but still far away from the black hole, and the silhouette depends on the cosmological constant as well. In particular, according to Fig. 3, the shadow will depend on the parameters of the mass function (4) for regular black holes as well. As we will see, the parameter r 0 , which generates regular metrics, can increase the distortion or the deviation from circularity and decrease the size of the shadow. As is known from the Kerr metric, the black hole shadow (its form or silhouette) is strongly sensible to the parameters a and θ o . However, for the class of regular black holes given by (1), we have two new parameters that modify the shadow: q and r 0 . Here it is shown just cases in which the rotation parameter is smaller than the black hole mass, i.e., a 2 < M 2 , and, in order to produce shadows compatible with M87*, it is adopted r 0 M, which, alongside a 2 < M 2 , provides the spacetime structure with three horizons indicated by (5). As we can see from Fig. 2 and Fig. 3, for θ o = 0, the more rotation, the more moved in the positive x-direction the shadow is. Such a motion points toward the rotation direction. In particular, the silhouette is deformed on the left (assuming that the rotation is from left to right) for large values of a, θ o and for r 0 . This difference between the left and the right sides is due to the photon orbits and the spacetime dragging. On the left, photons travel in the same direction of the black hole rotation, on the right they travel in the opposite direction. Points on the left and on the right in the shadow's silhouette are given by r p− and r p+ , respectively. In Fig. 3, we see the parameter r 0 can increase the shadow deformation for large values of θ o . On the other hand, small values of θ o can produce highly symmetrical shadows even for large values of a. Assuming that the M87* is observed at 17 • , we note that

Oblateness and deviation from circularity
In this section, two observables are built, oblateness and deviation from circularity. But following the Event Horizon Telescope collaboration, I will use the latter in order to compute an upper bound on the GUP parameter.

Oblateness
Following Refs. [3,[53][54][55], let us construct two observables for the black hole shadow. The first one indicates the difference between the x and y axes. Such a difference D is the oblateness. The second one, which was used for the Event Horizon Telescope collaboration, provides the shadow's deviation from circularity using the root-meansquare of the shadow's radius.
In order to evaluate D, one adopts an approximation for the Cartesian celestial coordinates x and y due to the large distance of the observer, r obs = (16.8 ± 0.8)Mpc, which turns the value of α into a small quantity. Thus, Eqs. (29) and (30) are rewritten as y(r p ) − sin α(r p ) cos β(r p ) .
The oblateness is simply defined as in which where r p− is solution of sin β = 1, and r p+ is given by sin β = −1. According to Fig. 4, the values x L and x R (left and right) indicate maximum and minimum values assumed by the shadow silhouette on the x-axis. On the other hand, y m is obtained from (32) and its maximum, i.e., dy(r p ) And the roots of Eq. (35) provide the value of r at which y m = y(r p ).
The parameter r = r p (r p− < r p < r p+ ) indicates the silhouette's maximum value on the y-axis. And, as we can see in all shadows shown here, ±y m are maximum and minimum values of y due to axial symmetry. Moreover, as we will see, the parameters x L , x R , and y m will be useful in order to construct the deviation from circularity using the root-meansquare of the shadow's radius. By using Eq. (33) for members of the class of rotating regular black holes studied here, Fig. 5 indicates a strong dependence between oblateness and the rotation parameter a. The more rotation, the more deformed the shadow is. Moreover, the parameter r 0 , which is related to the GUP parameter, modifies D, decreasing the oblateness.

Deviation from circularity
According to the Event Horizon Telescope collaboration [1], a measure for the shadow deformation is provided by the deviation from circularity, which is conceived of as deviation from the root-mean-square of the radius l(r p ). Following Bambi et al. [3] (but here another parameterization is adopted 6 ), the shadow radius l(r p ) is defined as with x c = x(r p ) and y c = 0 (see Fig. 4). The average radius (root-mean-square) is given by According to the cited authors, who follow the report from the collaboration, the deviation from circularity is conceived of as the root-mean-square distance from the average radius l R M S , that is to say, However, as I said, another parametrization is adopted and suggested here, demanding that (x R − x L ) is in the denominator of (39). Like the oblateness, ΔC increases with the black hole rotation and the observation angle (see Fig. 5). The Event Horizon Telescope reported ΔC < 0.1 assuming a Kerr metric as the geometry of M87*. Here Kerr-like objects are adopted, and using the M87* parameters the deviation from circularity given by Eq. (39) will provide an upper bound on the GUP parameter. The reported deviation from circularity < 10%, assuming the interval a * 0.99, imposes an upper bound on the GUP parameter, i.e., β 0 < 10 90 . In this graphic, it is adopted q = 1 in the mass function (4), which renders the best upper bound on β 0 , and Λ 1.1 × 10 −52 m

Estimating the quantum gravity parameter
In order to constrain the GUP parameter, it is necessary to relate it to the spacetime metric. In Ref. [38], it is applied a GUP in order to calculate quantum corrections to the black hole temperature. From the GUP where l p = G c 3 ≈ 10 −35 m is the Planck length, and β 0 is the so-called GUP parameter or the dimensionless quantum gravity parameter (as we can see, it is straightforward that for β 0 → 0 we have the standard uncertainty relation), we derived the quantum-corrected temperature for the Schwarzschild black hole and, interestingly, such a temperature is the Hawking temperature for the Bardeen regular black hole (up to second order approximation in l p /r + ). Thus, the Bardeen metric was interpreted as a quantum-corrected Schwarzschild black hole under the assumption that the GUP parameter could be related to a metric parameter (r 0 , the parameter that generates regular black holes introduced in Sect. 2), namely, Assuming that the relation (41) is valid for the entire class of regular black holes, for any q in Eq. (4), whether spherical or axisymmetric black holes, then the shadow of M87* can provide an upper bound on r 0 and, consequently, on β 0 . That is to say, the deviation from circularity (39) depends on the spacetime metric and its parameters like r 0 , thus, by using the relation (41), ΔC will depend on β 0 as well. Consequently, the upper bound on the GUP parameter can be Table 1 Upper bounds on the GUP parameter β 0 according to some approaches. Scanning tunneling microscope, Lamb shift and Landau Levels results are found in Ref. [40]. Gravitational waves result is found in Ref. [44], and the upper bounds from the light deflection, pulsar PRS B 1913+16 and perihelion precession were obtained in Ref. [43]. It is worth mentioning that the work on gravitational waves [44]  obtained by using ΔC(β 0 ) and the constraint provided by the Event Horizon Telescope Collaboration for the deviation from circularity. Therefore, by using ΔC(β 0 ), given by Eq. (39), and assuming the M87* parameters, i.e., M = (6.5 ± 0.7) × 10 9 M and r obs = (16.8 ± 0.8)Mpc, the deviation from circularity reported by the Collaboration, ΔC(β 0 ) < 0.1, imposes for a * 0.99 (with a * = a/M) and θ o = 17 • , for q = 1, 2 or q = 3 in the mass function (4). As we can see in Fig. 6, the constraint on the deviation from circularity gives an upper bound on the GUP parameter, i.e., curves that satisfy ΔC(β 0 ) < 0.1 with a * 0.99 are possible only for values given by Eq. (42). On the other hand, for a * > 0.99, the deviation from circularity is always ΔC(β 0 ) > 0.1. Thus, for θ o = 17 • , the parameter a * > 0.99 is ruled out according to M87* shadow. As I pointed out in Introduction, the value θ o = 17 • is conceived of as the angle between the jet and the observer, assuming that the jet direction is orthogonal to the M87* equatorial plane. The range for the rotation parameter presented here is in agreement with studies like [3], where the Kerr metric is adopted in order to describe the M87* shadow.
As we can see in Table 1, β 0 < 10 90 is the worst value compared to other upper bounds provided by different approaches, like the Lamb shift, Landau levels and scanning tunneling microscope in Ref. [40], or using light deflection, pulsar PRS B 1913+16 and perihelion precession in Ref. [43], or from the gravitational waves in Ref. [44]. According to Das and Vagenas [40], the scanning tunneling microscope delivers the best one, β 0 < 10 21 . On the other hand, like the present work, the authors of Refs. [43] and [44] present upper bounds based on gravitational phenomena, which rendered worse values than those provided by quantum systems. Quantum approaches have provided more stringent upper bounds on the GUP parameter.

Final remarks
The Event Horizon Telescope announced the first image of a black hole. According to the collaboration, the black hole shadow is well-described by the Kerr metric. However, it is argued that other options are still possible within the M87* parameters. Bambi et al. [3], for example, did not rule out a superspinar as a candidate for the geometry that produces the M87* shadow. In this article, a class of rotating regular black holes (thought of as a slight deviation from the Kerr metric) is adopted in order to describe the reported shadow.
The class of rotating regular black holes studied here presents a parameter according to which regular geometries or black holes without a singularity are generated. In our work [38], such a parameter from the spacetime metric was linked to the GUP parameter. GUPs appear in theories of quantum gravity and their parameters, which deform or generalize the Heisenberg principle, have been estimated. From a reported deviation from circularity of the M87* shadow (< 10%), an upper bound on the GUP parameter was computed. Alongside gravitational approaches adopted in order to obtain an upper bound on the GUP parameter, the value obtained here, β 0 < 10 90 , indicates that quantum options to estimate such a parameter are better alternatives than the gravitational options like gravitational waves, light deflection, pulsar PRS B 1913+16, perihelion precession, and the black hole shadow studied here. In general, quantum approaches provide more stringent upper bounds on the GUP parameter to date.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: This manuscript has no associated data.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .