Constraining the tidal charge of brane black holes using their shadows

A constraint on the tidal charge generated within a brane world is shown. Using the shadow of a rotating black hole in a brane context in order to describe the M87* parameters recently announced by the Event Horizon Telescope Collaboration, the deviation from circularity of the reported shadow produces an upper bound on the bulk's nonlocal effect, which is conceived of as a tidal charge in the four-dimensional brane induced by the five-dimensional bulk. Therefore, a deviation from circularity $\lesssim 10\%$ leads to an upper bound on the tidal charge $\lesssim 0.004M^2$.


I. INTRODUCTION
The recent shadow of the supermassive black hole M87* at the center of the Messier 87 galaxy built by the Event Horizon Telescope Collaboration (EHT) [1,2] is a new window for studies in strong gravitational field. The shadow cast by black holes depends on the metric parameters, and the main candidate for the M87* shadow is the Kerr geometry, despite the existence of alternatives to the M87* shadow like superspinars [3] and regular black holes [4,5]. The shadow reported by the collaboration presents a deviation from circularity less than ten percent, i.e., ∆C 10%. Using constraints coming from the EHT, researchers have obtained upper bounds on the rotation parameter of M87* [3], extra dimension length [6], cosmological parameters [7] and on the parameter of the generalized uncertainty principle [4]. The importance of the shadow phenomenon is due to use of an appropriate geometry for studies on the strong gravitational field, a geometry that describes a rotating black hole.
Shadows of objects have been studied since the Synge work [8] on the shadow of the nonrotating Schwarzschild black hole. In the same direction, Bardeen [9] built the first shadow for the Kerr geometry, which is a rotating black hole. Since those pioneer works, shadows of different black holes have been published, like shadows for the Reissner-Nordström black hole [10], the Kerr-Newman black hole [11], black holes with a cosmological constant [12], Kerr-Newman-NUT black hole [13], Kerr-Newman-Kasuya black hole [14], regular black holes [4,[15][16][17] and for braneworld black holes [18,19]. Shadows of rotating black holes with a cosmological constant in a brane context will be focused on this work in order to describe the M87* data and produce an upper bound on the tidal charge induced by the extra dimension.
The brane world adopted here is the second Randall-Sundrum model. The Randall-Sundrum model one (RS-I) [20] tries to solve the hierarchy problem, namely, the huge difference between the electroweak scale (∼ TeV) * Electronic address: nevesjcs@if.usp.br and the four-dimensional Planck scale (∼ 10 16 TeV). According to the RS-I, that problem is solved with introduction of a warp factor in the bulk metric, which is conceived of as a five-dimensional asymptotically antide Sitter spacetime with two four-dimensional branes as boundaries. One of those branes is our four-dimensional universe, and the other one is a hidden universe. In the RS-I, the extra dimension-from the bulk-is compact with a finite radius. On the other hand, the second Randall-Sundrum model (RS-II) [21], focused on this article, assumes only a brane, and the extra dimension radius is infinite. 1 In order to construct solutions of black holes in a Randall-Sundrum scenario, it is necessary to consider new four-dimensional gravitational field equations. Following Shiromizu et al. [24], induced equations were obtained from the Gauss-Codazzi equations in five dimensions. According to the authors, by assuming the RS-II-type model and the Z 2 symmetry in the bulk, the effective gravitational field equations are obtained on the brane. With those field equations on the brane, spherically and axially symmetric solutions (black holes with and without rotation, wormholes) were widely studied [25][26][27][28][29][30][31][32][33][34][35][36].
Shadows of rotating black holes in a RS-II-type scenario used here were built by Eiroa and Sendra [19] from our work [29], where such solutions or geometries were obtained for the first time. In such a context, the brane field equations carry the bulk's influence on the four-dimensional spacetime. In regard to rotating black holes in the RS-II-type model, the geometry that describes those objects presents a tidal charge, commonly conceived of as a nonlocal bulk's influence on the brane [22,29,35]. Therefore, I will constrain that tidal charge using the shadow silhouette and its deformation or deviation from circularity. By adopting the M87* parameters with the angle between the black hole rotation axis and the observer position θ obs = 17 • , which is in agreement with the observed jets supposedly aligned with the black hole rotation axis [1,37], the deviation from circularity reported by the collaboration ( 10%) provides an upper bound on the tidal charge 0.004M 2 . This result, as I pointed out, is obtained from a rotating black hole. On the other hand, constraints on the tidal charge have been reported using, in general, metrics with spherical geometry, especially a Reissner-Nordström-type geometry (with a tidal charge instead of a Coulomb charge [38][39][40][41][42]) that appears in brane contexts as it was reported in Ref. [34]. An exception is Banerjee et al. [43] in which rotating black holes (without a cosmological constant) are used in the brane context. According to those authors, negative values of the tidal charge are favored with the M87* data. However, positive values are not ruled out. In this sense, this work presents an upper bound on the tidal charge (the positive one), and, as we will see, the normalized constraint evaluated here using a rotating black hole with a cosmological constant is the second best found in the literature.
The structure of this paper is as follows: in Section II, the gravitational field equations and the spacetime metric that describes a rotating black hole with a cosmological constant in a brane world are presented. The equations to draw the shadow are shown. In Section III, it is indicated the procedure in order to evaluate the shadow's deviation from circularity. An upper bound on the tidal charge is obtained from that procedure. The final remarks are made in Sec. IV. 2

II. SHADOWS OF ROTATING BLACK HOLES WITH A COSMOLOGICAL CONSTANT IN A BRANE WORLD
This section deals with a solution of rotating black holes in a brane world with a cosmological constant published by us in Ref. [29], whose shadows were constructed in Ref. [19]. I intent to point out some features of that shadows to apply them in the next section for the deviation from circularity analysis. More details of the geodesics and shadows calculations are found in the two references mentioned above.

A. Brane world context: the gravitational field equations
In order to build rotating black holes with a cosmological constant in a different realm from general relativity, new field equations are necessary. In the brane context, the Einstein field equations are replaced by field equations projected onto the brane. Then, the brane is conceived of as a four-dimensional world immersed in a five-dimensional bulk, which is an asymptotically anti-de Sitter spacetime. The field equations on the brane were elegantly obtained by Shiromizu, Maeda, and Sasaki [24]. The induced vacuum field equations on the brane are written as in which G µν is the four-dimensional Einstein tensor related to the brane metric g µν , and Λ 4D is the fourdimensional brane cosmological constant. The tensor E µν is proportional to the (traceless) projection of the five-dimensional Weyl tensor onto the brane. It is worth noting that the vacuum field equations given by Eq. (1) reduce to the vacuum Einstein equations in the low energy limit.
By using the effective field equations, written without specifying the traceless tensor E µν , the trace of Eq. (1) delivers a constraint on the Ricci scalar given by In order to solve Eq. (2) and produce rotating black holes, it is assumed the axial symmetry in the brane indicated by the Kerr-Schild-(anti)-de Sitter Ansatz [44]. Accordingly, where the symbol ds 2 Λ indicates the AdS (anti-de-Sitter) or dS (de Sitter) metric. It is worth emphasizing that the main interest here, due to cosmological observations, is within the asymptotically dS rotating black holes in the brane. Lastly, the remaining parameters in Eq. (3) are: H(r, θ) is a function of the radial and polar coordinates (when spherical coordinates are adopted), and l µ stands for a null vector field. Substituting the Ansatz (3) into the constraint (2), one has H(r, θ) and the complete metric in the Kerr-Schild form (3). That is to say, give us a spacetime with axial symmetry that will be more evident with the aid of the Boyer-Lindquist coordinates. In such coordinates, the parameters M , a, and q will be properly interpreted.

B. Spacetime metric and geodesic equations
The spacetime metric (3) in the Boyer-Lindquist coordinates (t, r, θ, φ) is given by (details in Ref. [29]) where and From the Boyer-Lindquist coordinates, the parameters M , a, and q are interpreted as mass parameter, rotation parameter, and charge parameter, respectively. Indeed, according to Refs. [45,46], the parameter M of the Kerr-Newman-(A)-dS is related to the mass in asymptotically dS and AdS solutions, in such spacetimes the black hole mass is ±M/Ξ 2 (plus for AdS and minus for dS). As I pointed out, the brane is assumed empty, a vacuum brane. Then, it is common to conceive of q as a tidal charge (see [22,29,35]), i.e., it is a nonlocal bulk's influence on the brane. The charge feature of q is evident once one has the components of E µν : 3 In the brane world, the traceless tensor E µν is thought of as a geometric component of Eq. (1). As I said, such a tensor is a part of the five-dimensional Weyl tensor projected onto the brane and gives a geometric influence on the vacuum brane. However, Eq. (1) may be interpreted as the Einstein field equations with a traceless energy-momentum tensor as well (T µ µ = E µ µ = 0). In the general relativity realm, Eq. (9) is the Kerr-Newman-(A)-dS energy-momentum tensor since one replaces q by with Q 2 playing the role of the Coulomb charge. This is the main argument in favor of a tidal charge interpretation of q. But, contrary to Kerr-Newman-(A)-dS geometries, negative and positive values of the tidal charge present different spacetime structures (at least their scales). On the other hand, in the Kerr-Newman-(A)-dS spacetimes, positive and negative values of Q provide the same localization for horizons and the ergosphere. Because of the tidal charge is given by q instead of q 2 , negative values of q increase the ergosphere radius compared to positive ones (see details in Refs. [29,35]). In this sense, for q < 0, brane black holes are more appropriate to the Penrose mechanism of extracting energy from the black hole rotation.
As I am concerned with the cosmological context, the positive cosmological constant will be considered in this article, in agreement with the observations of an accelerating universe [47]. In this case, following Ref. [29], regarding a 2 < M 2 and q min < q < q max , the function ∆ r provides three positive roots, which are null surfaces: the inner horizon r − , the event horizon r + , and the cosmological horizon r ++ . Therefore, the spacetime structure is indicated by On the other hand, zeros of the metric component g tt provide us Killing surfaces and the limit of the ergosphere. With the M87* mass parameter M = (6.5±0.7)×10 9 M ⊙ for a 2 < M 2 , i.e., avoiding an overspinning black hole, the metric term g tt has one real root in between the event and the cosmological horizon (r = S). Thus, with this new element, the spacetime structure reads The region r + < r < S indicates the ergosphere. It is commonly claimed [34,35] that a negative value of the tidal charge is more "natural" since it confines the gravitational field near the brane. As we will see, this statement presents a new effect when it is considered the black hole shadow. For negative values of q, the brane black hole presents a less deformed shadow when it is compared to the shadow generated by a positive value of the tidal charge. A clear deviation from circularity is possible only for positive values of the tidal charge. According to previous works, positive values of the tidal charge generate more modifications on the brane like turning a spacelike singularity into a timelike singularity [34]. Due to this ability, positive tidal charges modify strongly the shadow's shape.
It is worth noting that the upper bound on the tidal charge obtained in Sec. III, by using the black hole shadow, will be better than q max evaluated in order to allow the spacetime structure given by (12).

C. Shadow and its silhouette
The black hole shadow is drawn from the geodesic equations. Geodesics for a rotating black hole, namely the Kerr metric, were obtained by Carter [48], who indicated the separability of the geodesic equations and argued that a test particle in the Kerr spacetime presents four constants of motion along the geodesics. Such constants are the two Killing vector fields ξ t and ξ φ with their respective constants, the mass of the test particle along the geodesics, and the last constant is called the Carter constant (in general, it is indicated by K). The Kerranti-de Sitter and Kerr-de Sitter spacetimes [49] and the solution with rotation in a brane world adopted in this article present those same four constants for test particles along geodesics [29].
From the two Killing vector fields ξ t and ξ φ indicated directly by the geometry with axial symmetry given by Eq. (6), one has the particle's constants of motion E and L along the geodesic, energy and angular momentum, respectively, All the geodesic equations for the spacetime given by Eq.
(6) were obtained in Ref. [29] using the Hamilton-Jacobi equation. For our purpose, only the radial and azimuthal components are useful: with and dot indicates derivative with respect to an affine parameter. The parameter δ is the particle's mass (one of the four constants along the geodesics), and for δ = 0, the case studied here, one has null geodesics. Following [4,13,19], two new parameters, which are constants in the shadow's silhouette, are defined as The shadow silhouette or its form is provided by unstable photon orbits with r = r p > r + constant outside the event horizon. For those orbits, photons may eventually either fall into the black hole or escape to the observer position r obs ≫ r p . Therefore, according to the geodesic equation for the radial coordinate,ṙ, one has R(r p ) = R ′ (r p ) = 0 ( in which the symbol ′ means derivative with respect to r) in order to deliver the unstable orbits that draw the shadow silhouette. In general, for black holes with rotation, r = r p− and r = r p+ , where r p− ≤ r p+ , and the black hole shadow is built from an interval of unstable orbits in which r p− is a minimum and r p+ is a maximum value. It is worth noting that for the Schwarzschild black hole r p− = r p+ = 3M and the shadow is circular, with M playing the role of the Arnowitt-Deser-Misner mass for such a black hole. But due to rotation, the Kerr black hole and Kerr-like black holes do not present circular shadows for any set of parameters. As we will see, those shadows are deformed by large black hole rotation and, in the brane context, the tidal charge q contributes substantially to the shadow's shape.
The main equations for the shadow calculation that involve the radial coordinate or R(r p ) and its derivative lead to As we will see, these conserved quantities are part of the shadow's equations, are constant in the shadow silhouette. The constants η(r p ) and ξ(r p ) are terms of the celestial coordinates that are used by the observer in order to draw the shadow. As we saw, the geometry studied here is not asymptotically flat, it is asymptotically dS. Therefore, an observer is not at infinity describing the shadow. For Λ 4D > 0, there is a cosmological horizon in the brane. In this sense, the observer is located at the so-called domain of outer communication, which is the region defined between the event horizon r + and the cosmological horizon r ++ in asymptotically dS spacetimes. Then, our observer is located at the point with coordinates (r obs , θ o ). Because of the axial symmetry, the black hole shadow (its silhouette) depends only on both the radial and polar coordinates. The coordinate θ o , which is the position in which the shadow is observed in relation to the rotation axis, stands for the observer angle, and the parameter r obs is the distance from the black hole. In this article, one assumes the M87* parameters to obtain a constraint on the tidal charge, that is to say, θ o = 17 • and r obs = (16.8±0.8)Mpc. As we are dealing with an asymptotically dS spacetime, it is also adopted the orthonormal tetrad e µ a = (e µ 0 , e µ 1 , e µ 2 , e µ 3 ) in order to describe the shadow silhouette (see Refs. [13,19]) in the domain of outer communication, according to the observer position. Then, the null congruences coming from r p− ≤ r ≤ r p+ reaches the distant observer and are projected onto the tetrad. The shadow phenomenon is described by using e µ a , which is written with the aid of the coordinate basis vectors (∂ t , ∂ r , ∂ θ , ∂ φ ), that is to say, . (19) According to the authors of Ref. [13], the direction of the vector e 3 points toward the black hole, and the component e 0 is the observer's four-velocity. Therefore, in this case, the observer is not necessarily at rest. It is worth mentioning that the orthonormal tetrad is chosen such that e 0 ± e 3 is tangential to the principal null congruence direction for a metric with axial symmetry like the metric studied here. Following Refs. [13,19], the light rays that come from the region defined by r p (the shadow's silhouette) are conceived of as curves γ or a null congruence such that their tangent vectors in the coordinate basis arė γ = (ṫ∂ t ,ṙ∂ r ,θ∂ θ ,φ∂ φ ). When projected onto the orthonormal tetrad,γ readṡ γ = ζ (−e 0 + sin α cos βe 1 + sin α sin βe 2 + cos αe 3 ) , (20) and that vector is given at the observer's position by using the chosen tetrad (19). The new angles α and β are the so-called celestial coordinates, and the description of the black hole shadow is entirely made by using these coordinates (see Fig. 1 from Ref. [4] for details). Indeed, as we will see, the celestial coordinates will be replaced by Cartesian coordinates in order to draw the shadow.
The factor ζ and the celestial coordinates α and β are straightforwardly obtained from Eq. (20),γ = (ṫ∂ t ,ṙ∂ r ,θ∂ θ ,φ∂ φ ) and the orthonormal tetrad equations, in such a way that and, at the same time, where the geodesic equations for the componentsṙ andφ, given by Eq. (14), were used (see details in Ref. [19]). It is worth emphasizing that the equations for the celestial coordinates depend on the observer's position (r obs , θ o ) and parameters of orbits r = r p that characterize the so-called photon "spheres". It is common in articles [4,13,19] that deal with asymptotically dS and AdS black holes to replace the celestial coordinates by Cartesian coordinates in order to describe the shadow's shape. In the Cartesian coordi- nates, the shadow is drawn by the parametric equations Before constructing the deviation from circularity, it is appropriate to adopt an approximation for the Cartesian coordinates x and y, for the value of α is a small quantity due to the large distance of the observer. Then, Eqs. (24) and (25) are rewritten as x(r p ) ≃ − sin (α(r p )) sin (β(r p )) , (26) y(r p ) ≃ − sin (α(r p )) cos (β(r p )) .
According to Fig. 1, the shadow is symmetrical to the x-axis. As is well known, for a Kerr black hole, the shadow's shape will depend strongly on the rotation parameter a and the observer angle θ o , if we consider an observer at infinity. On the other hand, as we saw, for an asymptotically dS black hole, the observer is in the domain of outer communication, and the silhouette depends on the cosmological constant as well. Moreover, as we can see in Fig. 1, due to the brane context, the tidal charge is another important parameter that influences the shadow.
Shadows for the metric (6) were widely built in Ref. [19]. Since the main point in this article is the tidal charge, let us focus on the influence of q on the shadow by emphasizing the angle θ o = 17 • reported by the EHT, according to which this is the polar angle between the observed jets along the rotation axis of M87* and our position. In the cited article of Eiroa and Sendra [19] such an information was not known yet.
According to shadows published in the literature, the more rotation, the more moved in the positive x-direction the shadow is. That is, that motion points toward the rotation direction. In particular, by assuming that the rotation is from left to right, the shadow's silhouette is deformed on the left for large values of a, θ o and, interestingly, for positive values of q. The difference between the left and the right sides is due to the photon orbits and the spacetime dragging of rotating black holes. On the left of the shadow, photons travel in the same direction of the black hole rotation, and on the right they travel in the opposite direction. From the geodesic equations, maximal points on the left and on the right in the shadow's silhouette are given by r p− and r p+ , respectively.
In Fig 1, it is shown just cases in which the rotation parameter is smaller than the black hole mass, i.e., a 2 < M 2 , in order to produce shadows compatible with M87* and a spacetime structure like (12). Notice that Fig. 1 indicates that negative values of q produces larger shadows than positive ones considering the same parameters (M, a, θ o , Λ 4D , |q|, r obs ). Indeed, positive values of the tidal charge induces a deformed silhouette, effect emphasized for large angles of observation ( Fig. 1 on the  bottom). Large values of negative tidal charges do not change the shadow circularity even for large values of θ o . Indeed, the deviation from circularity approaches zero as −q increases. Therefore, a clear deviation from circularity is just due to positive tidal charges, justifying an upper bound from the deviation only for positive values of q. This conclusion on the effect of positive and negative tidal charges is in agreement with previous works [22,35] in which the effect of positive tidal charges is considered more radical. And the strong effect of a large positive tidal charge motivated the authors of [43] to claim that a negative value of the tidal charge is more acceptable for the M87* parameters. However, following those authors, it is not possible to rule out positive values for the tidal charge to date.

III. AN UPPER BOUND ON THE TIDAL CHARGE
In this section, the observable defined as deviation from circularity is built. Following the EHT [1,2], I will use the deviation from circularity in order to compute an upper bound on the positive tidal charge q.

A. Calculating the deviation from circularity
Following Refs. [3,4,50,51], the shadow's deviation from circularity is provided by using the root-meansquare of the shadow's radius l(r p ), which is defined as The point (x c , y c ) is at the "center" of the shadow and not at the center of the Cartesian coordinates. Therefore, y c = 0 and x c = x(r p ′ ), where r = r ′ p (r p− < r p ′ < r p+ ) is the value of orbits in which the y-coordinate is maximum, i.e., in such orbits the Cartesian coordinate assumes y m = y(r p ′ ). The value y m is obtained from Eq. From the shadow's radius l(r p ), the average radius or root-mean-square reads According to Refs. [3,4,50,51], the desired deviation from circularity is given by the root-mean-square distance from the average radius l RMS , i.e., it is written as 4 (30) in which two new points of the shadow's silhouette are adopted: x L and x R (the most negative and the most positive values assumed by the shadow's silhouette on the x-axis). According to Eq. (26), these points are given by x L = − sin α(r p− ) and x R = sin α(r p+ ), where r p− is solution of sin β = 1, and r p+ is solution of sin β = −1 (see Refs. [4,51] for details). Some comments on the deviation from circularity ∆C: it increases with the black hole rotation and the observation angle (see Fig. 1), besides a positive tidal charge also increases ∆C. As mentioned before, the EHT reported ∆C 0.1 assuming a Kerr metric as the geometry of the M87* black hole. In this work, instead, a rotating black hole in a brane world with a cosmological constant is adopted. With an appropriate value for the tidal charge, it is possible to reproduce the M87* shadow.

B. Constraining the tidal charge
According to Eq. (30), the deviation from circularity depends on the metric parameters. Therefore, by using the M87* parameters (namely, M = (6.5 ± 0.7) × 10 9 M ⊙ and r obs = (16.8 ± 0.8)Mpc) alongside the value of the four-dimensional cosmological constant (Λ 4D ≃ 1.1 × 10 −52 m), the reported deviation from circularity imposes the following upper bound on the positive tidal charge: The above upper bound is obtained for a * 0.99 (with a * = a/M ) and θ o = 17 • . It is important to note that with a * > 0.99 and the M87* parameters, the spacetime structure with horizons given by Eq. (12) is not achieved. Thus, for θ o = 17 • and the reported deviation, the parameter a * > 0.99 is ruled out. According to EHT collaboration [1], the rotation parameter of M87* is still in dispute.
As mentioned before, there exists a maximum value of the tidal charge in order to produce the spacetime structure (12). Using the M87* parameters such a maximum value is q max 0.019M 2 for a * 0.99. As we can see, the constraint from the shadow improves that value. a Using light deflection around the Sun [38]. b Using trajectories of bright stars around the black hole [41]. c Using gravitational lensing of stars around Sagittarius A* [42]. d Using quasiperiodic oscillations in binary systems [39]. e Using gravitational lensing [40].
According to Table I, other constraints on the tidal charge have been reported in articles. Contrary to this work, the major part uses a nonrotating geometry in order to describe astrophysical black holes and constrain the tidal charge. The best normalized value q/M 2 comes from light deflection of the Sun [38]. A supermassive black hole, Sagittarius A*, was also used in [40][41][42], but the normalized upper bound reported here is even better.

IV. FINAL REMARKS
The most commented image in science in 2019 was the first image of a black hole announced by the EHT collaboration [1,2]. According to collaboration, the black hole shadow is well-described by the Kerr metric. However, as we saw, it is argued that options to the Kerr geometry are still possible within the M87* parameters. For example, Bambi et al. [3] do not rule out a superspinar or an overspinning black hole as a candidate for the object that produces the M87* shadow, and Ref. [4] suggests a rotating regular black hole as a candidate for M87*. Here, another option is still suggested: a brane black hole with rotation and cosmological constant.
The brane world context adopted here is the RS-IItype model, in which our four-dimensional universe is immersed in a five-dimensional bulk. As is well known in the large literature, the bulk acts nonlocally on the brane, generating a tidal charge even in the empty brane. Using a rotating black hole with a cosmological in the vacuum brane, the tidal charge was constrained with the aid of the shadow deviation from circularity from the central supermassive black hole M87* reported by the EHT. A deviation from circularity 10% indicated a normalized upper bound on the tidal charge 0.004, which is the second best result for the normalized tidal charge to date (the best one is 0.003 for a solar system result).