Some observable physical properties of the higher dimensional dS/AdS black holes in Einstein-bumblebee gravity theory

We study the greybody factors, quasinormal modes, and shadow of the higher dimensional de-Sitter (dS)/anti de-Sitter (AdS) black hole spacetimes derived from the Einstein-bumblebee gravity theory within the Lorentz symmetry breaking (LSB) framework. We specifically apply the semi-analytical WKB method and the time domain approach to study the scalar and Dirac perturbations of the black hole. In-depth researches are done on the effects of the LSB and dimensionality on the bosonic/fermionic greybody factors, quasinormal modes, and shadow of the higher dimensional bumblebee black hole. The results obtained are discussed, tabulated, and illustrated graphically.


I. INTRODUCTION
Since the Standard Model (SM) of particle physics [1] and older theories like general relativity (GR) [2], which describes how matter warps spacetime, cannot explain everything in the universe, including what occurs in the vicinity of a black hole, physicists are continually working to develop new and better ideas.Investigating that any retained concept such as Lorentz symmetry [3] may not be true in extreme cases is a very fruitful approach to explore for new physics.According to some gravitational wave models, the cosmos is not entirely symmetrical.Because of these ideas, the cosmos will always have extra elements that prevent it from perfectly adhering to the Lorentz symmetry.In other words, the cosmos would have a unique or favored orientation.These new models explain a theory known as "bumblebee gravity (BG)" [4][5][6].Its name comes from the alleged remark made by experts that bumblebees should not be allowed to fly since we did not know how their wings produced lift.Specifically, we do not fully comprehend how these gravity theories function and how they may be consistent with the universe that we currently observe.The possible use of bumblebee gravity models to explain dark energy [7], the phenomena that causes the cosmos to expand at an accelerated rate [8], is one of their most effective applications.It turns out that an effect that causes our universe to expand faster can be related to how much Lorentz symmetry our universe breaches.In addition, this notion seems very enticing because we do not know what is generating dark energy.In short, this is the direction that the bumblebee gravity theory is anticipated to contribute to GR and subsequently to the quantum gravity theory (QGT) [9,10].
Lorentz violations (LVs) affect both the predictions of the SM of particle physics and the basic predictions of special relativity, including the concept of relativity, the constancy of the speed of light in all inertial frames of reference, and time dilation.Test theories for special relativity and effective field theories, like the Standard Model Extension (SME) [11], have been developed to evaluate and forecast any violations [12].SME relates the SM to GR and includes additional features such as the LVs operating at the Planck scale [13].In other words, SME is an effective field theory integrating GR and SM on low energy scales.Other theories that propose the LVs besides SME include string theory [14], Einstein-aether theory [15,16], noncommutative field theory [17][18][19], loop QGT [20,21], brane-world scenarios [22,23], and massive gravity [24], also check [25].For the purpose of examining potential visible signs of a breach of particle Lorentz symmetry, the SME is experimentally accessible.In the model of SME, a spontaneous symmetry breaking potential caused by self-interacting tensor fields having vacuum expectation values (VEV), yields to the background tensor fields, which provides the local LV.An example of such a particular theory is the bumblebee field B µ , a self-interacting tensor field with a non-zero VEV with b µ that specifies a preferred direction in spacetime and spontaneously breaks the Lorentz symmetry.The potential of a bumblebee field can take many different forms.Among the other bumblebee's potentials, V (X) = kX 2 /2 is a smooth and functional potential, where k is a constant [26] and X = B µ B µ ± b 2 .When the bumblebee field matches its VEV, it has a minimum: V = 0 and V (X) = kX = 0 when X = 0.With this particular potential, the static black hole solutions in the Einstein-bumblebee gravity theory (EBGT) were derived by Bertolami et al. [27] and by Casana et al. [28].It was revealed by [37][38][39] that the solution obtained by Ding et al. [40], who claimed to have found a rotating black hole in the BG gravity model, was actually wrong.The first physically accepted slowly rotating black hole solution in the EBGT was wrong.The first physically accepted slowly rotating black hole solution in the EBGT was obtained by Ding and Chen [37].Subsequently, Jaha et al. [41] and Poulis and Soares [42] have recently managed to derive an arbitrarily spinning bumblebee black hole solution, which means crossing a difficult threshold in the EBGT.Another milestone in this regard was passed by [43], who achieved an exact higher dimensional anti-de Sitter (AdS) black hole solution in the EBGT.This AdS black hole can only exist with a special bumblebee potential having a linear functional form with a Lagrange-multiplier field λ.It is worth noting that, this additional field can be absorbed by the construction of an effective cosmological constant Λ e and is rigidly restricted by the equation of bumblebee motion.Furthermore, the obtained higher dimensional black hole of the EBGT is nothing but a Schwarzschild-AdS-like black hole solution since it cannot asymptotically approach anti-de Sitter spacetime, just as the Schwarzschild-like black hole [28].The bumblebee field has an impact on the black hole horizon location, in contrast to the Schwarzschild-like black hole [28].For the higher dimensional bumblebee metric, we compute very important observables in order to relate the bumblebee field to the spacetime geometry: the greybody factors (GbFs), quasinormal modes (QNMs), and shadow angular radius.Let us now briefly recall what the observables in question are: GbFs, which distinguish a black hole's thermal emission spectrum from a pure black-body spectrum, are functions of frequency, angular momentum, and black hole parameters.In other words, GbF is a quantity related to the quantum nature of a black hole, and there are different approches to computing the GbF [29][30][31][32][33][34][35][36].The modes of energy dissipation of a perturbed black hole or field are known as QNMs, and they characterize the perturbations of a field that dissipates with time.The solutions of the relevant perturbation equations that fulfill the boundary conditions necessary for purely incoming waves at the horizon and purely outgoing waves at infinity are represented by a black hole's QNMs.One must obtain the gravitational QNMs' spectra in order to examine the stability of higher-dimensional black hole solutions that might be present in nature.The stable and unstable black holes have a relationship with the damped and undamped states, respectively.Numerous techniques have been used to examine the QNMs frequencies, including the analytical method [44][45][46][47][48][49], WKB method [30,[50][51][52], Frobenius method [53], continuous fractions method [54], Mashhoon method [55], feedforward neural network method [56], and many more (for topical reviews, the reader is referred to [57][58][59]).In our work, we mainly focus on the WKB approximation method [50][51][52] to compute the QNMs.
The QNMs up to third order were first computed by Iyer and Will [51].Then, Konoplya [60] made it possible for us to compute the QNMs frequencies without using laborious numerical techniques, this resulted in a higher order contribution.Recently, the shadows of black holes have become one of the primary issues in physics.The reason for that is the Event Horizon Telescope Collaboration's debut photograph of a black hole, which was first released in 2019 [61,63,108].In fact, those pictures captured with the cutting-edge technology depict the shadow of M87 [64] and SgrA* [65,66], which are the supermassive black holes in the galaxy M87 and in the Milky Way Galaxy, respectively.But it was only during the last century that the first black hole's shadow was estimated.Synge [67] acquired what is known as the shadow of the Schwarzschild black hole today in the 1960s.Bardeen [68], thereafter extended Synge's work to the Kerr geometry.Recently, shadows have been considered for a number of black holes in a variety of scenarios.Recent researches suggest that there might be relationships between the shadow and black hole properties in general relativity or even in contexts outside of the Einsteinian paradigm [69][70][71][72].Shadow for the slowly rotating Kerr-like black hole in Einstein bumblebee gravity has been studied [73].We will be looking the effect of bumblee gravity on the shadow in higher dimension dS/AdS spacetime.
The main purpose of this paper is to study the perturbations of scalar and fermion fields in the higher dimensional dS/AdS black hole geometries of the EBGT and the shadows of those higher dimensional black holes.For the perturbations, we shall consider the Klein-Gordon and Dirac equations.The obtained wave equations allow for semi-analytical methods to be used for GbF and QNM analyses.Then, we consider the photon's orbit and radius of the shadow of the black hole.This article is structured as follows: In Sec.II, we briefly introduce the higher dimensional dS/AdS black hole solutions in the EBGT [43].Then, we discuss the GbFs of bosons in Sec.III.Section IV is devoted to Dirac equation of massless fermions on the higher dimensional dS/AdS black hole spacetime.We also compute the rigorous lower bounds on the fermionic GbFs.Then, the QNMs are studied in Sec.V. Null geodesics and shadow radius problems are discussed in Sec.VI.The purpose of Sec.VII is to investigate the connection between shadow radius and QNMs.In Sec.VIII, we check the effect of the bumblebee parameter on the shadow diameter via the real black holes.We draw our conclusions in Sec.IX.

II. HIGHER DIMENSIONAL dS/AdS BLACK HOLES IN EBGT
The bumblebee vector field B µ in the EBGT, has included a vacuum expectation value which is nonzero, in order to define a unconstrained Lorentz symmetry breaking via a given potential.The action of Einsteinbumblebee gravity in higher dimensions D ≥ 4 is given by [37,43], where, Λ is the cosmological constant.κ = 8πG D /c 4 , where G D = GΩ D−2 /4π [74] and The bumblebee field B µ takes a nonzero vacuum expectation value (VEV) < B µ >= b µ at these minima which tells us that vacuum of this model has a preferred direction in the spacetime.The vector b µ here is a constant function of spacetime which has a value b µ b µ = ∓b 2 , where the ± signs denote timelike or spacelike forms of the vector b µ .The bumblebee field strength is given by We get the following constraint on the B µν due to its antisymmetric nature [26], Varying the action (1) with respect to the metric yields the following field equation: where G µν = R µν −g µν R/2 and T B µν is known as the bumblebee energy momentum tensor, which is expressed by In the above expression, V represents the differentiation of V computed at x = B µ B µ ± b 2 .Then, varying the action (1) with respect to t, the bumblebee field gives the following field equation by assuming that there is no coupling between the bumblebee field and the Lagrangian of matter: We now suppose that there is no matter field and the bumblebee field is frozen at its VEV.Namely, we have (see, for example, [27,28]) Now, since we intend to include the cosmological constant in our theory, the non-zero cosmological constant requires the linear form of the potential as being stated in [43]: where, λ is a non-zero constant and considered as a Lagrange-multiplier field.The potential vanishes for condition (7) and the derivative of the potential V = λ 2 modifies the Einstein field equation.Hence, Eq. (4) recasts in [37,43], where, Now, we would like to construct a D-dimensional static and spherically symmetric metric in the EBGT.To this end, let us consider the following metric anstaz: Since, the spacetime considered has a strong radial variation compared to the temporal changes, we consider that the bumblebee field has a radial finite vacuum expectation value.Hence, the spacelike bumblebee field turns out to be b µ = (0, be ψ(r) , 0, 0...., 0), (12) where, b is a positive constant.The bumblebee field strength is defined by whose components and their divergences are now all zero.Therefore, From Eq. ( 6), we can see the projection of the Ricci tensor along the bumblebee field: Using Eq. ( 9), we will have three independent equations: where the prime symbol denotes differentiation of a function with respect to its argument and the Lorentzviolating parameter is given by L = b 2 ≥ 0. Now, Eq. ( 16) leads to the following metric function: where, where M is the mass of the black hole.To have the Schwarzchild-like solution [28] for Λ = 0, we set Hence, the bumblebee field reads and from Eqs. ( 15) and ( 18), one can see that the following expression for the cosmological constant should hold: which puts a constraint on the parameter λ from the potential (8).Therefore, λ is not a new degree of freedom in the theory.Moreover, one can define an effective cosmological constant Λ e as follows which means that Λ = (1 + L)Λ e .After all those computations, we get the final form of the metric as where the metric function reads It is clear from the metric function ( 26) that the event horizon is affected by the bumblebee field.The behavior of the metric function is illustrated in Fig. (1) for both dS and AdS spacetimes with different dimensions.One can observe that we have one horizon (event) for the AdS case however for the dS case double horizons appear: event horizon (inner) and cosmological (outer) horizon (see Fig. 1).It is worth noting that we have found the similar kinds of behaviors for the varying LSB parameter, which are depicted in Fig. 2, in which D = 4 is fixed for both for AdS and dS cases.

III. GbFs of BOSONS
Among the various methods for deriving the GbFs of the relevant black hole spacetimes, only a small number of cases renders possible to derive their precise analytical expressions.In this section, we will apply the rigorous bounds technique to procure the GbFs of higher dimensional dS/AdS black hole in EBGT.To this end, the excitation of uncharged and massless scalar fields is going to be determined by the Klein-Gordon equation:  where denotes the D'Alembert operator.So, Eq. ( 27) can be rewritten as in which for our D-dimensional metric (25) √ −g is given by In order to get separate radial and angular Klein-Gordon equations, let us apply the following ansatz [75]: where ω indicates frequency, l represents the azimuthal quantum number, and m(−l ≤ m ≤ l) is the spherical harmonic index.The angular equation yields the eigenvalue (λ) [76] as Using all of the above equations we can write the radial equation for the scalar field, where a prime mark denotes a derivative with respect to the radial coordinate, r.Applying the following transformation one can rewrite Eq. ( 32) as At this stage, by using the tortoise coordinate dr * = √ 1 + L dr f , we get a Schrödinger-like wave equation: where To understand the potential behaviors, we have plotted Eq. ( 36) for AdS and dS spacetimes in different dimensions (see Fig. 3).We have observed that for the AdS spacetime, the potential vanishes only once while for the dS spacetime, it vanishes twice, which is the outcome of the double horizons obtained in the dS spacetime.We have found similar behaviors for the LSB parameter (see Fig. 4 with fixed D = 4-dimension for AdS and dS backgrounds).The general semi-analytic bounds for the GbFs are given by where We have two conditions for the certain positive function h . After applying the conditions to Eq. ( 38), one may observe a direct proportionality between the GbFs and the effective potential, where the metric function plays a significant part in this process.Since there is no upper border in the integral of Eq. ( 37), without loss of generality one can set h = ω 2 − V ef f .Thus, Eq. (37) becomes which results in where h peak = ω 2 − V peak .Eq. ( 40) can also be rewritten as For evaluating V peak , first r peak should be determined for different sub-cases.The behaviours of the obtained GbFs for the scalar particles are depicted in Fig. (5) for the AdS black hole of the EBGT.What is interesting in Fig. (5) is that while the 4-dimensional black hole has the highest GbF values, the 5-dimensional black hole has the weakest GbF values.However, the other higher dimensions (D > 5) have GbF values between the 4 th and 5 th dimensions.In D > 5 dimensional black holes, the GbF values decrease as the dimension increases.Moreover, it is seen that the bosonic GbFs of the AdS bumblebee black hole are almost unaffected by the change in the LSB parameter.
For positive cosmological constant, by considering the second condition, Eq. ( 37) is expressed by whose integration is solvable.Thus, we have The behaviour of GbFs for Λ e > 0 (dS) is depicted in Fig. (6) to show the influences of the dimension and the LSB parameter on the GbFs of the higher dimensional dS black hole in the EBGT.The most important finding from Figs. (5) and ( 6) is that the GbF of the 4-dimensional black hole in the EBGT theory is higher than those of its higher dimensional versions.Namely, the GbF drastically reduces with the increasing dimensions, which means that the probability for detecting the thermal radiation of the higher dimensional black holes in the EBGT gets lower with D > 4.

IV. GbFs of FERMIONS
In this section, we shall investigate the GbFs of the Dirac particles i.e., fermions.To derive the 1-dimensional Schrödinger like wave equation, we apply a particular conformal transformation, which contains the Dirac Lagrangian invariant [77,78].Under the aforementioned conformal transformation [79,80], one has If we consider Ω = 1/r, the metric (25) becomes Thus, the t − r and (D − 2)-sphere parts of the metric are separated and whence the Dirac equation can be rewritten as γµ ∇µ ψ = 0, (49) which has the following expansion where (γ 5 ) 2 = 1: we can now omit the bar.Furthermore, let us consider χ (±) l as the eigenspinors for the (D − 2)-sphere [81].Then, we have where l = 0, 1, 2, 3... .We can also consider ψ as the orthogonal eigenspinors: Therefore, the Dirac equation ( 50) can be written as follows where γ 5 is the interaction term presented in the two dimensional Dirac equation.To tackle with the Dirac equation, we get help from the following auxiliary expressions: and where the σ i are the known Pauli matrices, . Therefore, we can write the spin connections as We will use the positive sign from now on, as both signs work similarly and we can use either one.Hence, we can explicitly rewrite the Dirac equations as follows After letting the following ansatz we get a simplified form of the Dirac equation: From the above expression, we get Recalling the tortoise coordinate , Eqs. ( 62) and ( 63) are simplified to We can now decouple the above equations as where These two potentials V 1,2 belong to the particle and anti-particles of Dirac fermions.We have shown the behaviour of the potentials V 1,2 for the different dimensions in dS and AdS space (see Figs. 7 and 8).Depending on the existence of the cosmological horizon in AdS/dS spacetimes, the potential vanishes at some radial distance as being observed in the scalar potential.Now, by considering the above potentials (68) in Eq. ( 42), the fermionic GbFs of the EBGT, for various dimensions, are obtained as the following: For D = 4: for D = 5:  and for D = 6: In the above results [Eqs.( 69)-( 71)], γ = Λ e (1 + L), the asymptotic series approach is utilized in order to facilitate the integration evaluations.That is why the GbFs are served in discrete forms for different dimensions.Since the integration is bounded between the outer and cosmological horizons, to define the GbFs of the dS black hole, we directly evaluate the integration (42) and get The behaviors of the Dirac GbFs for various dimensions and LSB parameters are illustrated in Fig.   that the most highest fermionic GbFs in both cases (dS/AdS) belong to the D = 4 case.The increase in the dimension decreases the GbFs.In order to have a wider perspective also on the impression of parameter L in the GFs of the bumblebee model refer to [29] as an example of constant dimension D = 4.

V. QNMs
The WKB (Wentzel, Hendrik Kramers and Léon Brillouin) approach is an approximate technique to solve the linear differential equations.The most significant utilization of the WKB approximation is to solve the time independent Schrödinger equation.In general, the equations conducting different types of non-rotating or static black hole QNM perturbations form in terms of the radial coordinate: see Eq. (35), in which ω stands for the complex QNMs.As is well-known, in quantum mechanics, the ω 2 parameter corresponds to 2m h2 E, where E indicates the particle energy.In this context, the effective potential is nothing but a barrier.In order to compute the QNMs, the appropriate boundary conditions must be considered at (r * → −∞) and (r * → ∞), which stand for the event horizon and spatial infinity, respectively.The ingoing modes represent waves moving away from the potential barrier.Namely, the ingoing waves (r * → −∞) correspond to the radiation which crosses the horizon into the black hole.On the other hand, since a QNM is occurred from a black hole's free oscillation, the ingoing modes at spatial infinity are ignored.In other words, only outgoing waves exist at spatial infinity, (r * → ∞).
To compute the QNMs with the WKB approach, we employ the following complex frequency expression [30] of fermionic QNMs, we have observed that the real terms for Λ < 0 correspond to the imaginary terms for Λ > 0. This implies that both positive and negative values of Λ exhibit similar behavior.

VI. NULL GEODESICS AND SHADOW RADIUS
In this section, we shall study the photon's orbit and radius of the shadow of the black hole [82].Let us first consider the Lagrangian L(x, ẋ) = (1/2)g µν ẋµ ẋν for the static spherically symmetric metric, which is given by ( 25) The spacetime has two conserved quantities which can be calculated by solving the Euler-Lagrange equations in the equatorial plane.Hence, the conserved quantities are obtained as where E and L are called the conserved specific energy and conserved specific angular momentum.For the photons, we can write the following equation: which gives the following equation by the aid of the conserved quantities where Now we define the impact parameter b = L/E.At the photon sphere radius r = r ph , the conditions dr/dφ| r ph = 0 (V ef f = 0) and V ef f = 0 should be satisfied.Hence, the impact parameter for the photon sphere is given by By using conditions dr/dφ| r ph = 0 and d 2 r/dφ 2 | r ph = 0 [83], we can find the radius of the photon sphere as where B(r) = r 2 f (r) .Therefore, Eq. ( 79) becomes Now, we define an angle α [84] between the null light ray and radial direction as follows which gives Using some trigonometric relations, one can get the following equation: We define b cr as the impact parameter at the critical impact parameter, thus we have Ultimately, the radius of the shadow for the static observer at r = r 0 [85] is found to be and by considering a mathematical constraint f (r 0 ) ≈ 1 for a static observer located at a special location (see for example [86,87] in which the appropriate normalization for the time-like Killing vector was applied), we can write the shadow radius [88] as follows Now, in terms of the celestial coordinate, the shadow radius is obtained as where (r 0 , θ 0 ) is the position of the observer at spatial infinity.Moreover, since we analyze the shadow of the black hole in the equatorial plane, the shadow of the radius is equivalent to the critical impact parameter of the photon sphere.Therefore, we have which is explicitly written as We have shown the variation of the shadow radius of the black hole with the space-time dimension in Fig. 12 (left) for Λ e = −0.01.After initially declining, the shadow radius then begins to rise.A similar effect is observed in the celestial coordinate, which is shown in Fig. 13 (left).Similarly, we have shown the variation of the shadow radius with the bumblebee parameter (L) in Fig. 12 (right) for D = 4 and Λ e = −0.01.The shadow radius in this case keeps on decreasing with L and the shadow of the black hole is depicted in celestial coordinate in Fig. 13 (right).We also looked at the shadow size for the dS space and took the observed cosmological constant value Λ e = 1.11 * 10 −52 m −2 .We first showed the variation of the shadow size with space-time dimension D in Fig. 14 (left) for M = L = 1.The shadow size variation is similar to the case of Ads but with a larger radius.A similar effect is observed in the celestial coordinate, which is shown in Fig. 15 (left).We also observed that the bumblebee gravity parameter does not affect the shadow size at all in any dimension and we particularly have plotted the shadow size for D = 4 case in Fig. 14 (right) for M = 1 and the shadow of the black hole is depicted in celestial coordinate in Fig. 15 (right).In this section, we will try to reveal the relation between the shadow radius and QNMs.As shown in [89][90][91][92][93], the real part of the QNMs at the eikonal limit corresponds to the angular velocity of the critical null circular orbit Ω c and the imaginary part of the QNMs is nothing but the Lyapunov exponent λ, which is used to determine the unstable timescale of the circular orbit [94].Namely, we have where the angular velocity is given by where r c is the radius of the circular null geodesics.Equation (93) allows us to write a relation between the QNMs and shadow radius at the eikonal limit [90] Re Expression ( 95) is only valid for the large orbital quantum numbers (l).But, Konoplya and Stuchik [95] demonstrated that this may not always be the case.The results seen in Table (III) are based on Eqs.(93) and (95).As can be deduced from Table (III), the results obtained are in agreement with Eq. ( 95), however they are inconsistent from the ones obtained for scalar QNMs which are tabulated in Table (I).Nevertheless, one can conclude that the relationship between the shadow radius and QNMs is reliable for l 1, which covers the majority of the cases in black hole physics.

VIII. CONSTRAINT ON BUMBLEBEE PARAMETER USING EHT DATA OBTAINED FOR M87* AND Sgr A*
A lot of effort has been paid to evaluating the effects of black holes on the astrophysical environment [101][102][103][104].In particular, we can point to supermassive black holes because they are typically found near the center of galaxies [105,106].The supermassive black hole is located at the center of the nearby gigantic elliptical galaxy Messier 87 (M87), also known as M87 * , according to astrophysical measurements cited in [107].The Event Horizon Telescope (EHT) was constructed in order to image the shadow of M87 * and the supermassive black holes in the heart of the Milky Way (Sgr A * ).Recent investigations have shown that M87 * has a shadow, as shown by [108][109][110].The results achieved, as we all know, were fantastic.In this section, we shall constrain the bumblebee parameter by using the data supplied by the EHT for M87 * and Sgr A * .As it was reported in [108], angular diameter of the M87 * black hole shadow is θ M87 * = 42 ± 3 µas, distance of M87 * from the Earth is measured as d M 87 * s = 16.8Mpc, and mass of the M87 * is M M87* = 6.5±0.90x10 9 M .Similarly, for Sgr A * the data for its shadow is given in the latest EHT paper [111].It is reported that the angular diameter of the Sgr A * shadow is θ Sgr A * = 48.7 ± 7 µas, distance of the Sgr A * from the Earth is d SgrA * s = 8277 ± 33 pc and mass of the Sgr A * black hole is M Sgr A * = 4.3 ± 0.013x10 6 M .Now by using those data and following the formula, one can calculate the diameter of the black hole shadow [112], From above, radial diameters of the shadow images for M87 * and Sgr A * can be obtained as d M87* sbh = (11 ± 1.5)M and d Sgr A * sbh = (9.5 ± 1.4)M , respectively.However, for our metric, the diameter of the black hole shadow can also be calculated from (92).Therefore, variation of the diameter of the black hole shadow with the bumblebee parameter for different dimensions can be easily studied as we depicted in Fig. 16 for Λ e = −0.01 and in Fig. 17 for Λ e = 1.11 * 10 −52 m −2 .We have considered 1σ and 2σ uncertainties to show the constraints properly.The pertinent figure clearly shows that the bumblebee parameter has a range for D = 4 that falls within a range of uncertainty, but for higher dimensions, this is not the case for both M87 * and Sgr A * for AdS case.However for the dS case, as we saw in the earlier plots L does not affect the shadow size and therefore for D = 4, we have all possibility to consider any L value and for higher dimensions, it again falls below the sigma regions.Consequently, in the future, if some observations are made, which suggest a smaller shadow size than what we have observed so far, the current study can thus reveal the importance of the existence of bumblebee gravity.

IX. CONCLUSIONS
In this work, we have performed a comprehensive discussion on GbFs in higher dimensional AdS/dS black hole spacetimes of the EBGT.The study has provided impressive results in higher dimensions when considering general relativity coupled to the bumblebee gravity.To compute the GbFs, we have considered the scalar and Dirac field perturbations.To analyze the obtained radial wave equations, we have employed the WKB approach up to sixth order to derive the GbFs.The effects of LSB or bumblebee parameter (L), cosmological constant (Λ e ), and dimension (D) on the GbFs have been thoroughly investigated, which provide significant impacts on the thermal radiation.It has been observed from Fig. ( 5) that the 4-dimensional AdS bumblebee black hole has the highest GbF values, whereas the 5-dimensional black hole has the weakest GbF values.Other higher dimensions (D > 5) have GbF values between the fourth and fifth dimensions, though.The GbF values drop as the dimension rises in D > 5 dimensional black holes.On the other hand, the bosonic GbFs of dS higher dimensional bumblebee black holes are shown in Fig. (6), which decrease regularly with increasing dimension.It is also understood that at higher dimensions, the LSB effect on the bosonic GbFs of the AdS higher dimensional bumblebee black holes is very weak compared to the dimension effect.The GbFs of the dS higher dimensional bumblebee black holes are effectively reduced by the increasing LSB, though.Fig. (9) shows the fermionic GbF behaviors of the higher (D ≥ 4) dimensional AdS/dS bumblebee black holes against the change in dimension and the LSB parameter.In general, irrespective of being a bosonic or fermionic perturbation, the most important findings obtained from Figs. 5, 6, and 9 are that the GbFs of the 4-dimensional dS/AdS bumblebee black hole in the EBGT theory are higher than those for D > 4 black holes.Namely, GbFs drastically reduce with the increasing dimensions, which means that the probability for detecting the thermal radiation of the higher dimensional ds/AdS black holes in the EBGT gets lower with D > 4. We have also noticed that the GbFs reach to 1 quicker in the scalar field perturbations compared to the fermionic field perturbations.This indicates that the bosonic thermal radiations can more likely reach to spatial infinity in comparison to the fermionic thermal emission.
Bosonic and fermionic QNMs with l = 0 and l = 1 cases of the higher dimensional dS/AdS bumblebee black holes are tabulated in Tables (I) an (II), respectively.We have inferred from those tables that both real and imaginary parts of the bosonic and fermionic QNMs increase with growing dimensionality.Almost the same behaviors obtained for the rising L parameter except for the damping mode are obtained at D = 6 and D = 8 dimensions for bosons and fermions, respectively.We have also extended our investigation to find a direct link between the QNMs frequencies and the shadows of the bumblebee black holes.The shadows of the higher dimensional ds/AdS bumblebee black holes have been studied in terms of null geodesics and spherical photon orbits.By changing the dimensionality, the black holes' shadow radii have been depicted and analyzed.It is clearly shown in Fig. 12 (left) that an increasing dimension (D ≥ 4) decreases the radius of the black hole's shadow.We have also shown that the latter result can also be obtained by means of the real part of the QNMs frequencies, which are valid in the eikonal limit.Finally, we have considered the null geodesics and obtained the black hole shadow radius with different dimension.Then, we have exhibited the outcomes in Fig. (13) (left).Similarly, we have shown the variation of the shadow radius with the bumblebee parameter (L) in Fig. 12 (right) for dimension D = 4.It is clear that shadow radius decreases with increasing the L. The same behaviour is depicted in the celestial coordinate in Fig. 13 (right).In the sequel, we have examined the variation of the black hole shadow diameter with the bumblebee parameter based on the data of real black holes ( M87 * and Sgr A * (see Fig. (16)).At the end of the day, we have shown that the bumblebee parameter has a reducing effect on the diameter of the black hole shadow.This result can be used as a tool to indirectly prove the existence of the bumblebee gravity theory.
In recent times, the search for proof of the existence of bumblebee gravity has has gained momentum.Especially, the research of Gu et al. [96] based on real black hole X-ray data is quite remarkable.In this context, this study will contribute to the examination of EBGT with possible optical and wave observations to be made in the future.Our research can be expanded to charged bumblebee black holes, which will require taking into account the bumblebee electrodynamics [97][98][99].A possible charged bumblebee black hole solution will allow us to examine the EBGT with linear/non-linear electrodynamics theories [100] and we will likely obtain more detailed results on the thermal radiation (GbFs), gravitational ringing (QNMs), and optical observations (lensing and shadow).This is the next stage of study that interests us.
is the area of a unit D − 2 sphere.From now on we will take G D = 1 and c = 1 for simplicity.b is a positive constant and L M represents the matter Lagrangian form.The strength of the non-minimal coupling of gravity with the bumblebee field B µ is determined by the coupling constant .The potential term V (B µ B ν ∓ b 2 ) represents Lorentz or CPT (charge, parity and time) violation.The potential has a minima at B µ B ν ± b 2 = 0 and V (b µ b µ ) = 0, which destroy the U (1) symmetry.

FIG. 6 :
FIG. 6: Graph of scalar GbFs for Λe > 0 (dS) having various D-dimensions and the bumblebee parameter L. As the solid lines are for L = 1, the dashed lines stand for L = 20.The physical parameter are chosen as M = l = 1 and Λe = 0.01

( 9 )
FIG. 9: Graph of Dirac GbFs for Λe = −0.1 (left) and Λe = 0.1 (right) in various dimensions D ≥ 4. As the solid lines represent GbFs for spin-up, the dashed and dotted lines stand for the GbFs of the spin-down particles.The physical parameters are chosen as M = l = 1.

FIG. 10 :
FIG. 10: Graph of QNMs for the scalar particle for various dimensions; Left figure stands for l = 0 in which the green line represents L = 1 and red line is for L = 1.5.Right figure stands for l = 1 in which the orange line is for = 1 and blue line exhibits L = 1.5.The physical parameters are chosen as M = 1 and Λe = −0.01.

FIG. 12 :
FIG. 12: In the left-hand side we have shown the shadow radius of the bumblebee black hole versus D-dimension.The physical parameters are chosen as M = L = 1 and Λe = −0.01.On the right-hand side, we have shown the shadow radius of the bumblebee black hole versus the bumblebee parameter L. The physical parameters are chosen as M = 1, D = 4, and Λe = −0.01.

FIG. 14 :
FIG. 14: In the left-hand side we have shown the shadow radius of the bumblebee black hole versus D-dimension.The physical parameters are chosen as M = L = 1 and Λe = 1.11 * 10 −52 m −2 .On the right-hand side, have shown the shadow radius of the bumblebee black hole versus the bumblebee parameter L. The physical parameters are chosen as M = 1, D = 4, and Λe = 1.11 * 10 −52 m −2 .

10 FIG. 16 :
FIG.16:The variation of the diameter of the shadow with respect to the bumblebee parameter for different dimensions.Here, we have considered 1σ and 2σ uncertainties for the M87 * (left) and Sgr A * (right) for Λe = −0.01.

10 FIG. 17 :
FIG.17:The variation of the diameter of the shadow with respect to the bumblebee parameter for different dimensions.Here, we have considered 1σ and 2σ uncertainties for the M87 * (left) and Sgr A * (right) for Λe = 1.11 * 10 −52 m −2 .

TABLE I :
Bosonic QNMs of various dimensional dS/AdS black holes ,

TABLE II :
Fermionic QNMs of various dimensional dS/AdS black holes ,